<div dir="ltr">Dear Molpro Users/Developers,<div><br></div><div>I have been trying to compute potential energy curves for the dissociation of C2 molecule for the three low-lying singlet states. I run SA-CASSCF for 3 states and specify LQUANT,0,0,2, because I would like to compute the two lowest Sigma^+_g states and one Delta_g state. The SA-CASSCF computation runs fine and produces qualitatively correct potential energy curves. However, when computing potential energy curves near conical intersection in the region of 2.9-3.2 Angstrom, state-specific NEVPT2 computations for the two lowest-energy states converge to the same root:</div><div><br></div><div>             1.1                 2.1                 3.1<br>2.8 -75.6722945 -75.6489181 -75.6359653</div>2.9 <b>-75.6433987 -75.6433986</b> -75.6299956<br>3.0 <b>-75.6360085 -75.6360085</b> -75.6217944<div>3.1 <b>-75.6274115 -75.6274115</b> -75.6114991<br>3.2 <b>-75.6181054 -75.6181054</b> -75.5990394</div><div><br></div><div><div>Similar situation is observed with state-specific CASPT2. As far as I understand, the problem arises at the initial CI step, which produces the same reference state in both state-specific computations, despite the fact that the reference SA-CASSCF energies are different.</div><div><br></div><div>Is there a way to avoid this problem to generate correct NEVPT2 potential energy curves near conical intersection? </div><div><br></div><div>I run MOLPRO 2010.1. Here's my input file:</div><div><br></div> memory,500,m<br> gprint,orbitals=5,basis,civector<br> gthresh,energy=1.0d-10,orbital=1.0d-8<br><br> do i=1,32<br><br> BOHR<br><br> Geometry = {<br> C<br> C 1 R<br> }<br><br> R = 1.8 + 0.1 * (i - 1)<br><br> basis=cc-pvdz<br><br> {rhf;orbital,2100.2}<br> {multi,maxit=50;start,2100.2;closed,1,0,0,0,1,0,0,0;occ,3,1,1,0,3,1,1,0;wf,12,1,0;state,3;lquant,0,0,2}<br> {nevpt2,maxiti=500;core,0;wf,12,1,0;state,1,1}<br> {nevpt2,maxiti=500;core,0;wf,12,1,0;state,1,2}<br> {nevpt2,maxiti=500;core,0;wf,12,1,0;state,1,3}</div><div><br> e(i) = energy<br> rvec(i) = R<br><br> enddo<br><br> table,rvec,e</div><div><br><div><br clear="all"><div><div class="gmail_signature"><div dir="ltr"><div>Best wishes,<br></div><div><div><br></div><div>Alexander Sokolov</div><font color="#888888"><div>Postdoctoral Research Associate</div><div>Princeton University</div><div>Princeton, NJ</div></font></div></div></div></div>
</div></div></div>