## SOC matrix element between high multiplicities

(It is not necessary to upload input files, they can be found in the output file.)

manishrana
Posts: 1
Joined: 14 Jun 2018, 12:06
First name(s): Manish
Last name(s): Kumar
Affiliation: IISER Bhopal
Country: India

### SOC matrix element between high multiplicities

Dear users and developers,

I want to calculate the SOC matrix elements between quintet and septet spin multiplicities. But on this forum, I have only seen discussion on SOC matrix elements between singlet and triplet. Can I calculate the SOC matrix elements between higher spin multiplicities in Dalton?

Thanks

taylor
Posts: 600
Joined: 15 Oct 2013, 05:37
First name(s): Peter
Middle name(s): Robert
Last name(s): Taylor
Affiliation: Tianjin University
Country: China

### Re: SOC matrix element between high multiplicities

The spin-orbit operator can be thought of as having the action DeltaS = +/-1 when operating on a given wave function. If that wave function has S=0 then there is only one result: DeltaS=+1 (so singlet/triplet coupling). If you operate on a wave function with S other than zero (or 1/2) you can get results from DeltaS=+1 and DeltaS=-1. So if you operate on quintets you can get triplets and septets. As long as you are comfortable sorting through the results of the response calculation and identifying which final states are which it should be possible to do what you want.

Be aware that the only computational level that can be used for this is MCSCF. There are no methods in Dalton for including both spin-orbit and dynamical correlation for open-shell states that I am aware of. One could perhaps use open-shell DFT but I would be very sceptical of the quality of the results because the target quantities would not have been used in parametrizing the functional.

Best regards
Pete

olav
Posts: 145
Joined: 28 Aug 2013, 06:20
First name(s): Olav
Last name(s): Vahtras

### Re: SOC matrix element between high multiplicities

Theoretically you can only do this with Dalton in the context of a CI response calculation, assuming e.g. that you have a converged quintet reference state and run a CI linear response with (frozen orbital) calculation on top of that, plus, you would have to manually analyse the multiplicity of the excited states based on their determinant expansions. If the question is if there is an input option to accomplish this in a black-box manner, the answer is no.

Regards,
Olav

### Who is online

Users browsing this forum: No registered users and 9 guests