<table cellspacing="0" cellpadding="0" border="0" ><tr><td valign="top" style="font: inherit;"><DIV>Hi，</DIV>
<DIV> </DIV>
<DIV>Attached please find the example of geometry optimization with MRCI+SOC. The lowest Omega=3/2 state (2Pi + 2Delta) of the linear HePtF complex was optimized, in which the energies of 2Pi and 2Delta were calculated at the MRCISD+Q level whereas SOC was at the CASSCF level. Since the 2Sigma+ energy was calculated at the CASSCF level (because 2Sigma+ does not contribute to 3/2), the first two spinor states are 3/2 instead of the ground spinor state 1/2.</DIV>
<DIV> </DIV>
<DIV>At some points, the 5/2 state may be a little lower than 3/2. In this case you have to stop the calculation and optimize energy(3) or energy(4) instead. You can also modify the source codes of molpro to identify 3/2 and 5/2 and then swich between energy(1)(2) and energy(3)(4).</DIV>
<DIV> </DIV>
<DIV>The details can be found in Dalton Trans., 39, 2023 (2010)</DIV>
<DIV> </DIV>
<DIV>Best regards,</DIV>
<DIV>Wenli</DIV>
<DIV> </DIV>
<DIV><BR>--- On <B>Thu, 5/12/11, Jeong Sik Lim <I><lim.jeongsik@gmail.com></I></B> wrote:<BR></DIV>
<BLOCKQUOTE style="BORDER-LEFT: rgb(16,16,255) 2px solid; PADDING-LEFT: 5px; MARGIN-LEFT: 5px"><BR>From: Jeong Sik Lim <lim.jeongsik@gmail.com><BR>Subject: [molpro-user] energy gradient of MR-SO-CISD wavefunction<BR>To: molpro-user@molpro.net<BR>Date: Thursday, May 12, 2011, 12:19 AM<BR><BR>
<DIV id=yiv1359644115>
<STYLE><!--
#yiv1359644115  
 _filtered #yiv1359644115 {font-family:"Malgun Gothic";panose-1:2 11 5 3 2 0 0 2 0 4;}
 _filtered #yiv1359644115 {font-family:"Malgun Gothic";panose-1:2 11 5 3 2 0 0 2 0 4;}
#yiv1359644115  
#yiv1359644115 p.yiv1359644115MsoNormal, #yiv1359644115 li.yiv1359644115MsoNormal, #yiv1359644115 div.yiv1359644115MsoNormal
        {margin:0cm;margin-bottom:.0001pt;text-align:justify;text-justify:inter-ideograph;font-size:10.0pt;font-family:"Malgun Gothic";}
#yiv1359644115 a:link, #yiv1359644115 span.yiv1359644115MsoHyperlink
        {color:blue;text-decoration:underline;}
#yiv1359644115 a:visited, #yiv1359644115 span.yiv1359644115MsoHyperlinkFollowed
        {color:purple;text-decoration:underline;}
#yiv1359644115 span.yiv1359644115EmailStyle17
        {font-family:"Malgun Gothic";color:windowtext;}
#yiv1359644115 .yiv1359644115MsoChpDefault
        {font-family:"Malgun Gothic";}
#yiv1359644115  
 _filtered #yiv1359644115 {margin:3.0cm 72.0pt 72.0pt 72.0pt;}
#yiv1359644115 div.yiv1359644115WordSection1
        {}
--></STYLE>

<DIV class=yiv1359644115WordSection1>
<P class=yiv1359644115MsoNormal><SPAN lang=EN-US>Dear all,</SPAN></DIV>
<P class=yiv1359644115MsoNormal><SPAN lang=EN-US> </SPAN></DIV>
<P class=yiv1359644115MsoNormal><SPAN lang=EN-US>So far, I am doing spin-orbit calculation using ECP. But I am not sure about that energy gradient of MR-SO-CISD wavefunction of specific interested state can be abstracted for geometry optimization? If does, please inform me how to do.</SPAN></DIV>
<P class=yiv1359644115MsoNormal><SPAN lang=EN-US> </SPAN></DIV>
<P class=yiv1359644115MsoNormal><SPAN lang=EN-US>Thank you so much.</SPAN></DIV></DIV></DIV><BR>-----Inline Attachment Follows-----<BR><BR>
<DIV class=plainMail>_______________________________________________<BR>Molpro-user mailing list<BR><A href="http://us.mc151.mail.yahoo.com/mc/compose?to=Molpro-user@molpro.net" ymailto="mailto:Molpro-user@molpro.net">Molpro-user@molpro.net</A><BR><A href="http://www.molpro.net/mailman/listinfo/molpro-user" target=_blank>http://www.molpro.net/mailman/listinfo/molpro-user</A><BR></DIV></BLOCKQUOTE></td></tr></table>