<div dir="ltr">Dear Molpro users,<div><br></div><div>Let say I have HF orbitals of a molecule and I want to to CCSD(T) of this molecule.</div><div>Now I would like to correlate just some orbitals. including some valance (but not all valance) and some low-lying core orbitals. So I have to rotate the core orbitals with the valance orbitals that I do not want to correlate and do only one iteration. For example:</div><div><br></div><div><pre style="white-space:pre-wrap;color:rgb(0,0,0)">{hf
wf,108,1,0
shift,-1.0,-0.5
maxit,1
start,2101,2
closed,18,9,9,3,9,3,3,0
occ,   18,9,9,3,9,3,3,0
rotate,9.1,16.1,90
rotate,4.2,9.2,90
rotate,4.3,9.3,90
rotate,4.5,9.5,90</pre><pre style="white-space:pre-wrap;color:rgb(0,0,0)"><span style="font-family:arial">}</span></pre><pre style="white-space:pre-wrap;color:rgb(0,0,0)"><span style="font-family:arial">where 9.1, 4.2, 4.3, 4.5 are core orbitals and 16.1, 9.2, 9.3 and 9.5 are valance orbitals that I do not want to correlate.</span></pre><pre style="white-space:pre-wrap;color:rgb(0,0,0)"><span style="font-family:arial">After that I do a CCSD(T) calculation based on this HF calculation.</span></pre><pre style="white-space:pre-wrap;color:rgb(0,0,0)"><pre style="white-space:pre-wrap">{ccsd(t)
core,16,9,9,2,9,2,2,0
maxit,50
shift,1
}</pre></pre><pre style="white-space:pre-wrap"><font face="arial"><font color="#000000">I would like to ask if it is make sense to do this? </font></font></pre><pre style="white-space:pre-wrap"><font face="arial"><font color="#000000">I also notice that after rotation and do the HF calculation, the core orbitals (</font></font><span style="font-family:arial;color:rgb(0,0,0)">9.1, 4.2, 4.3, 4.5</span><span style="color:rgb(0,0,0);font-family:arial">)  are now high-lying and have the energy equal the the energy of the valence orbitals that were rotated. I am not sure if this is correct. T</span><span style="color:rgb(0,0,0);font-family:arial">o be clear:</span></pre><pre style="white-space:pre-wrap"><span style="color:rgb(0,0,0);font-family:arial">Before rotation (full HF calculation):</span></pre><pre style="white-space:pre-wrap"><pre style="white-space:pre-wrap;color:rgb(0,0,0)">ORBITAL  9.1 (core)    eig=   -4.4491
ORBITAL 16.1 (valance) eig=   -0.9945
<br></pre></pre><pre style="white-space:pre-wrap"><font face="arial"><font color="#000000">After rotation, HF one iteration:</font></font></pre><pre style="white-space:pre-wrap"><pre style="white-space:pre-wrap;color:rgb(0,0,0)">ORBITAL  9.1 (valance)  eig=   -4.4491
ORBITAL 16.1 (core)     eig=   -0.9945</pre></pre><pre style="white-space:pre-wrap"><font face="arial"><font color="#000000"><br></font></font></pre><pre style="white-space:pre-wrap"><font face="arial"><font color="#000000">Best regards,</font></font></pre><pre style="white-space:pre-wrap"><font face="arial"><font color="#000000">Quan </font></font></pre></div><div><font face="arial"><font color="#000000"><br></font></font></div></div>