<div dir="ltr"><div>Greetings!</div><div>I am using MOLPRO 2012 and would like to have ROHF for Li2 with the following:</div><div><br></div><div>alpha occupancy:  3   0   0   0   1   0   0   0</div><div>beta  occupancy:   1   0   0   0   1   0   0   0</div><div><br></div><div>But even if I have:</div><div>occ,3,0,0,0,1,0,0,0<br></div><div><br></div><div>The ROHF program prints:</div><div>"Orbital guess generated from atomic densities. Full valence occupancy:</div><div> 2   0   0   0   2   0   0   0<br></div><div><br></div><div><div>Initial alpha occupancy:  2   0   0   0   2   0   0   0</div><div>Initial beta  occupancy:   1   0   0   0   1   0   0   0",</div></div><div><br></div><div>and the SCF converges to this (incorrect) solution.</div><div><br></div><div>I was not able to fix this by using "sym 1" in the wf card, and/or</div><div>by adding  "closed,1,0,0,0,1,0,0,0".</div><div><br></div><div>Said another way, the input file:</div><div>========</div><div><div>geometry={angstrom;Li1;Li2,Li1,3.065}</div><div>basis = cc-pVDZ</div><div>wf,charge=0,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>rhf</div></div><div>========<br></div><div>gives an ROHF occupancy with sym=5 and:</div><div><div>alpha occupancy:  2   0   0   0   2   0   0   0</div><div>beta  occupancy:   1   0   0   0   1   0   0   0</div></div><div><br></div><div>regardless of whether or not I have "sym=1" and/or the "closed" card specified as above.</div><div><br></div><div>What can I do to force the correct initial guess to have the occupancies I want? </div><div><br></div><div>Respectfully yours,</div><div>Nike Dattani</div>
</div>