<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
<meta name="Generator" content="Microsoft Word 15 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0in;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri","sans-serif";}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:#0563C1;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:#954F72;
text-decoration:underline;}
span.EmailStyle17
{mso-style-type:personal-compose;
font-family:"Calibri","sans-serif";
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
font-family:"Calibri","sans-serif";}
@page WordSection1
{size:8.5in 11.0in;
margin:1.0in 1.0in 1.0in 1.0in;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang="EN-US" link="#0563C1" vlink="#954F72">
<div class="WordSection1">
<p class="MsoNormal">Hi, <o:p></o:p></p>
<p class="MsoNormal">I am trying to use RASPT2 for direct dynamics calculations, therefore analytical gradients would be very useful. I did a gradient test calculation, resulting in the following error message (see input below):
<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">nci (cpmcscf)= 490 nci(rs2)= 187<o:p></o:p></p>
<p class="MsoNormal">? Error<o:p></o:p></p>
<p class="MsoNormal">? wrong nci<o:p></o:p></p>
<p class="MsoNormal">? The problem occurs in rs2grad_load<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">I should also mention that I’m using a level shift. This error happens after the RS2 energy has been calculated and the CP-CASPT2 equations have been solved. It appears to me that another call to the MULTI program is made, where the configuration
restrictions are not passed correctly. 187 is the no. of CSFs in the RASSCF reference, while 490 is the corresponding no. for the full CASSCF(8e7o) wavefunction. Are analytical gradients only available for CASSCF references, not for general MCSCF?<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">I tried both versions 2012.1 and 2015.1 without success. Any thoughts or suggestions on this issue?<o:p></o:p></p>
<p class="MsoNormal">Thanks, <o:p></o:p></p>
<p class="MsoNormal">Mark<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">The following is my input file: <o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">***,gradient_test_raspt2<o:p></o:p></p>
<p class="MsoNormal">file,2,record2,new<o:p></o:p></p>
<p class="MsoNormal">memory,800,M<o:p></o:p></p>
<p class="MsoNormal">angstrom<o:p></o:p></p>
<p class="MsoNormal">geomtyp=xyz<o:p></o:p></p>
<p class="MsoNormal">geometry={<o:p></o:p></p>
<p class="MsoNormal"> 9<o:p></o:p></p>
<p class="MsoNormal"> m05<o:p></o:p></p>
<p class="MsoNormal"> C 0.09401600 0.77190600 -0.18413000<o:p></o:p></p>
<p class="MsoNormal"> C -1.14074500 -0.45052700 0.27666600<o:p></o:p></p>
<p class="MsoNormal"> H 0.11432900 0.74310700 -1.26932600<o:p></o:p></p>
<p class="MsoNormal"> H -0.20211800 1.70842000 0.27790800<o:p></o:p></p>
<p class="MsoNormal"> H -0.65384200 -1.29464400 -0.24692500<o:p></o:p></p>
<p class="MsoNormal"> H -1.04027900 -0.49614400 1.37312500<o:p></o:p></p>
<p class="MsoNormal"> O 1.27204600 0.35159800 0.39429100<o:p></o:p></p>
<p class="MsoNormal"> O 1.87822000 -0.61969200 -0.26472400<o:p></o:p></p>
<p class="MsoNormal"> O -2.25269500 -0.01465300 -0.20748500<o:p></o:p></p>
<p class="MsoNormal"> }<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">basis=sto-3g<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">{mcscf<o:p></o:p></p>
<p class="MsoNormal">closed,16<o:p></o:p></p>
<p class="MsoNormal">occ, 23<o:p></o:p></p>
<p class="MsoNormal">wf,40,1,0<o:p></o:p></p>
<p class="MsoNormal">}<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">put,molden,sto3g.molden<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">basis=cc-pVDZ<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">{mcscf<o:p></o:p></p>
<p class="MsoNormal">rotate,14.1,18.1,0.<o:p></o:p></p>
<p class="MsoNormal">rotate,23.1,24.1,0.<o:p></o:p></p>
<p class="MsoNormal">closed,16<o:p></o:p></p>
<p class="MsoNormal">occ, 23<o:p></o:p></p>
<p class="MsoNormal">wf,40,1,0<o:p></o:p></p>
<p class="MsoNormal">restrict,4,6,17.1,18.1,19.1<o:p></o:p></p>
<p class="MsoNormal">restrict,0,2,22.1,23.1<o:p></o:p></p>
<p class="MsoNormal">}<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">put,molden,ccpvdz.molden<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">{rs2,shift=0.20<o:p></o:p></p>
<p class="MsoNormal">closed,16<o:p></o:p></p>
<p class="MsoNormal">occ, 23<o:p></o:p></p>
<p class="MsoNormal">wf,40,1,0<o:p></o:p></p>
<p class="MsoNormal">restrict,4,6,17.1,18.1,19.1<o:p></o:p></p>
<p class="MsoNormal">restrict,0,2,22.1,23.1<o:p></o:p></p>
<p class="MsoNormal">maxiter,200,200}<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">force<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">---<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">----------------------------------------------------<o:p></o:p></p>
<p class="MsoNormal">Dr. Mark Pfeifle<o:p></o:p></p>
<p class="MsoNormal" style="margin-bottom:12.0pt">Chemical Sciences & Engineering Division – Fundamental Interactions
<br>
Argonne National Laboratory <br>
9700 S Cass Ave <br>
Lemont, IL, 60439 <u><span style="color:blue"><br>
</span></u><a href="mailto:mpfeifle@anl.gov"><span style="color:#0082BF">mpfeifle@anl.gov</span></a>
<br>
+1 (630) 252-3717 <o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
</body>
</html>