[molpro-user] Spin-orbit coupling between Triplet and Singlet states

[Kirk Peterson]

Kelvin,

try:

{ci;hlsmat,ls,4500.1,4600.1;print,vls=1;option,matel=1}

best,

-Kirk

PS - I believe the problem in your previous input was not specifying enough arguments of the tranls option.


> On May 6, 2015, at 10:52 PM, Kelvin Lee <kin.lee at unsw.edu.au> wrote:
> 
> Dear Molpro mailing list,
> 
> 
> 
> I am trying to learn how to calculate the spin-orbit matrix element between a singlet and triplet state. As a starting point, I've been trying to calculate CH2 using CASSCF and MRCI, as there is a lot of published data on this molecule.
> 
> I've tried following the example given in the manual on the S-atom, but the results I'm getting on CH2 make very little sense to me, and I'd appreciate it if anyone would extend their help.
> 
> My input is as follows:
> 
> ***,Calculating spin-orbit coupling of CH2
> memory,148,M
> angstrom
> gthresh,orbital=1.d-8,energy=1.d-8,step=1.d-5,gradient=1.d-5
> nosym
> geom={
> C
> H1 C rCH
> H2 C rCH H1 aHCH
> }
>  RCH=                 1.12276332 ANGSTROM
>  AHCH=              102.15789308 DEGREE
> basis=dzp
> {uhf;occ,5;wf,8,0,2         ! UHF of triplet ground state
> orbital,4400.1
> }
> {multi;occ,10;wf,nelec=8,spin=2;wf,nelec=8,spin=0}
> 
> {ci;wf,nelec=8,spin=2;save,4500.1}        ! triplet wavefunction
> {ci;wf,nelec=8,spin=0;save,4600.1}        ! singlet wavefunction
> 
> {ci;tranls,4500.1,4600.1}                 ! calculate matrix element
> 
> 
> Which, at the end of the output gives me:
> 
>  Spin-orbit matrix elements for mean field operator:
> 
>  !MRCI expec          <1.1|LSZ|1.1>     0.000000000000i au =     0.000000006576i cm-1
> 
> If I'm interpreting this correctly, this is supposed to be my matrix element for these two states?
> The previously calculated values lie around ~12 cm-1 or so, so what unbelievably stupid mistake am I making?
> 
> Any form of help would be greatly appreciated.
> 
> 
> 
> 
> Thank you,
> 
> 
> Kelvin
> 
> 
> -- 
> KIN LONG KELVIN LEE | PhD Student
> Molecular Photonics Group | School of Chemistry | Faculty of Science
> University of New South Wales
> Room 417, Chemical Sciences Building F10 | NSW | AUSTRALIA | 2052
> T 9385-0457 |
> E kin.lee at unsw.edu.au <mailto:k.lee at chem.usyd.edu.au> |
> 
>  
> _______________________________________________
> Molpro-user mailing list
> Molpro-user at molpro.net
> https://urldefense.proofpoint.com/v1/url?u=http://www.molpro.net/mailman/listinfo/molpro-user&k=EWEYHnIvm0nsSxnW5y9VIw%3D%3D%0A&r=ndj2Vta%2FisUMMOFOMobgvg%3D%3D%0A&m=9UiJjnccwnmhaMZfbsX0%2F9dQDSJgRct1KmcWUVBvLFU%3D%0A&s=24c8332ab40af1e5499c50e976caa6ccff81d92a0035386188b8170b306bb918

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.molpro.net/pipermail/molpro-user/attachments/20150507/ddb9c853/attachment.html>