<div dir="ltr">Dear Jacky, and Kirk,<div>Both of your methods worked!</div><div><br></div><div>In summary, using MOLPRO 2012:</div><div><br></div><div>==========================</div><div><span style="font-size:12.8px">geometry={angstrom;Li1;Li2,</span><wbr style="font-size:12.8px"><span style="font-size:12.8px">Li1,3.065}</span><br style="font-size:12.8px"><span style="font-size:12.8px">basis = cc-pVDZ</span><br style="font-size:12.8px"><span style="font-size:12.8px">wf,charge=0,spin=2,sym=1</span><br style="font-size:12.8px"><span style="font-size:12.8px">occ,   3,0,0,0,1,0,0,0</span><br style="font-size:12.8px"><span style="font-size:12.8px">closed,1,0,0,0,1,0,0,0</span><br style="font-size:12.8px"><span style="font-size:12.8px">core</span><br style="font-size:12.8px"><span style="font-size:12.8px">rhf</span><br></div><div>==========================<span style="font-size:12.8px"><br></span></div><div><br></div><div>Gives: RHF STATE 1.1 Energy =  -12.635662909597<br></div><div><br></div><div>If I wrap up the RHF part in brackets:</div><div><br></div><div><div>==========================</div><div><div> {rhf</div><div> wf,<span style="font-size:12.8px">charge=0</span>,spin=2,sym=1</div><div> occ,3,0,0,0,1,0,0,0</div><div> closed,1,0,0,0,1,0,0,0</div><div> core</div><div> }</div></div><div>==========================<span style="font-size:12.8px"><br></span></div></div><div><br></div><div>I get: RHF STATE 1.1 Energy =  -14.821433808867</div><div><br></div><div>Even if I wrap up everything but the "occ" part:</div><div><div><br></div><div>==========================<br></div><div>wf,<span style="font-size:12.8px">charge=0</span>,spin=2,sym=1</div><div>closed,1,0,0,0,1,0,0,0</div><div>core</div><div><br></div><div>{rhf</div><div>occ,3,0,0,0,1,0,0,0</div><div>}</div></div><div>==========================<br></div><div><br></div><div>I get the correct energy and occupations.</div><div><br></div><div>It appears that the SCF program simply does not remember the "occ" cards defined before the RHF command, even though it DOES remember charge, spin, sym, and closed.</div><div><br></div><div>Is this enough to be considered a bug? <br>Or was "occ" never designed to work the same way as "closed" ?</div><div><br></div><div>Does my first input work in MOLPRO 2015 ?</div><div><br></div><div>==========================</div><div><br></div><div>Another puzzle is that even with Kirk's input format:</div><div><br></div><div><span style="font-size:12.8px">angstrom</span><br style="font-size:12.8px"><span style="font-size:12.8px">geometry={Li1;Li2,Li1,3.065}</span><br style="font-size:12.8px"><span style="font-size:12.8px">basis = aug-cc-pCV5Z</span><br style="font-size:12.8px"><span style="font-size:12.8px">{rhf;occ,3,0,0,0,1,0,0,0</span><br style="font-size:12.8px"><span style="font-size:12.8px">closed,1,0,0,0,1,0,0,0</span><br style="font-size:12.8px"><span style="font-size:12.8px">wf,6,1,2</span><br style="font-size:12.8px"><span style="font-size:12.8px">}</span><br></div><div><br></div><div>The occupations are correct, but the RHF energy is few cm-1 higher than CFOUR, </div><div>and the CCSD(T) energy is 0.1 Hartree higher than CFOUR at aug-cc-pCV5Z level:</div><div><br></div><div>RHF__CFOUR = -14.824,356</div><div>RHF_MOLPRO = -14.824,341</div><div><br></div><div>CCSD(T)__CFOUR = -14.91920357</div><div>CCSD(T)_MOLPRO = -14.82979310<br></div><div><br></div><div>Respectfully yours,</div><div>Nike Dattani</div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div></div>