<div dir="ltr"><div>Hello,</div><div> </div><div>I'm running a MRAQCC single point energy calculation with the following input:</div><div> </div><div>{mcscf,maxit=100<br>occ,26<br>closed,20<br>wf,nelec=46,spin=0,charge=0,sym=1<br>

 print,orbitals,virtuals<br>cpmcscf,hess<br>}</div><div> </div><div>{aqcc,maxit=250<br>core,6<br>wf,nelec=46,spin=0,charge=0,sym=1<br>natorb,print=20<br>}</div><div> </div><div>The natural orbitals from the AQCC step suggest that the six core orbitals have double occupancy but the remaining 14 closed orbitals do not. Also the virtual orbitals from number 27 onwards have non-zero occupancies. Is there an error in my input that is resulting in the closed orbitals not maintaining an occupancy of two? Why are the virtual orbitals outside the active space occupied?</div>

<div> </div><div>I obtain the same solution whether I include the occ and closed commands in the AQCC block or not.</div><div> </div><div>Thank you for your help.</div><div>David.<br clear="all"></div><div>Postgraduate research student</div>

<div><a href="mailto:david.l.scarborough@gmail.com">david.l.scarborough@gmail.com</a></div></div>