<div dir="ltr"><div style="font-family:arial,sans-serif;font-size:13px">I am wondering if there is an easy way to retain unrestricted type symmetry-breaking during the SCF procedure using Molpro. I have included below an example for H2, where DIIS takes me back to the RHF solution even though there a lower energy UHF solution.  In this case it would be easy to generate the correct wavefunction but I need this for computing UHF solutions of transition metal systems where it isn't possible to write down the correct wavefunction.    </div><div style="font-family:arial,sans-serif;font-size:13px"><br></div><div style="font-family:arial,sans-serif;font-size:13px">Thank you,</div><div style="font-family:arial,sans-serif;font-size:13px">EJS</div><div style="font-family:arial,sans-serif;font-size:13px"><br></div><div style="font-family:arial,sans-serif;font-size:13px">***, H2-Example</div><div dir="ltr" style="font-family:arial,sans-serif;font-size:13px">memory,1500,M</div><div dir="ltr" style="font-family:arial,sans-serif;font-size:13px">NOSYM</div><div dir="ltr" style="font-family:arial,sans-serif;font-size:13px">geometry={</div><div dir="ltr" style="font-family:arial,sans-serif;font-size:13px">2</div><div dir="ltr" style="font-family:arial,sans-serif;font-size:13px"><br></div><div dir="ltr" style="font-family:arial,sans-serif;font-size:13px">H 0 0 0</div><div dir="ltr" style="font-family:arial,sans-serif;font-size:13px">H 0 0 1.5</div><div dir="ltr" style="font-family:arial,sans-serif;font-size:13px">}</div><div dir="ltr" style="font-family:arial,sans-serif;font-size:13px">basis=sto-3g</div><div dir="ltr" style="font-family:arial,sans-serif;font-size:13px">NOENEST=yes</div><div dir="ltr" style="font-family:arial,sans-serif;font-size:13px">{rhf,maxdis=50,iptyp='DIIS',maxit=500,orbital=2100.2,accu=10,energy=1e-12}</div><div dir="ltr" style="font-family:arial,sans-serif;font-size:13px"><br></div><div dir="ltr" style="font-family:arial,sans-serif;font-size:13px">{MERGE</div><div dir="ltr" style="font-family:arial,sans-serif;font-size:13px">ORBITAL,2100.2</div><div dir="ltr" style="font-family:arial,sans-serif;font-size:13px">MOVE,1.1,,1.1</div><div dir="ltr" style="font-family:arial,sans-serif;font-size:13px">ROT,1.1,2.1,30</div><div dir="ltr" style="font-family:arial,sans-serif;font-size:13px">SAVE,2100.2</div><div dir="ltr" style="font-family:arial,sans-serif;font-size:13px">}</div><div dir="ltr" style="font-family:arial,sans-serif;font-size:13px">NOENEST=yes</div><div dir="ltr" style="font-family:arial,sans-serif;font-size:13px">{uhf,maxdis=50,iptyp='DIIS',maxit=500,orbital=2100.2,accu=10,energy=1e-12}</div></div><div class="gmail_extra"><br clear="all"><div><div dir="ltr">-------------------------------------------------------------------<br>Dr. Eric Jon Sundstrom                      <br><p dir="ltr" style="line-height:1;margin-top:0pt;margin-bottom:0pt"><span style="font-family:Arial;vertical-align:baseline;white-space:pre-wrap">Postdoctoral Researcher</span></p><p dir="ltr" style="line-height:1;margin-top:0pt;margin-bottom:0pt"><span style="font-family:Arial;vertical-align:baseline;white-space:pre-wrap">Thomas F. Miller III Group </span></p>Division of Chemistry and Chemical Engineering<br>California Institute of Technology<div>ph: 765.491.5765  <br>email: <a href="mailto:esundstr@gmail.com" target="_blank">esundstr@gmail.com</a></div></div></div>
<br><div class="gmail_quote">On Tue, Sep 23, 2014 at 1:01 PM, Eric Jon Sundstrom <span dir="ltr"><<a href="mailto:esundstr@gmail.com" target="_blank">esundstr@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div dir="ltr"><div>I am wondering if there is an easy way to retain unrestricted type symmetry-breaking during the SCF procedure using Molpro. I have included below an example for H2, where DIIS takes me back to the RHF solution even though there a lower energy UHF solution.  In this case it would be easy to generate the correct wavefunction but I need this for computing UHF solutions of transition metal systems where it isn't possible to write down the correct wavefunction.    </div><div><br></div><div>Thank you,</div><div>EJS</div><div><br></div><div>***, H2-Example</div><div dir="ltr">memory,1500,M</div><div dir="ltr">NOSYM</div><div dir="ltr">geometry={</div><div dir="ltr">2</div><div dir="ltr"><br></div><div dir="ltr">H 0 0 0</div><div dir="ltr">H 0 0 1.5</div><div dir="ltr">}</div><div dir="ltr">basis=sto-3g</div><div dir="ltr">NOENEST=yes</div><div dir="ltr">{rhf,maxdis=50,iptyp='DIIS',maxit=500,orbital=2100.2,accu=10,energy=1e-12}</div><div dir="ltr"><br></div><div dir="ltr">{MERGE</div><div dir="ltr">ORBITAL,2100.2</div><div dir="ltr">MOVE,1.1,,1.1</div><div dir="ltr">ROT,1.1,2.1,30</div><div dir="ltr">SAVE,2100.2</div><div dir="ltr">}</div><div dir="ltr">NOENEST=yes</div><div dir="ltr">{uhf,maxdis=50,iptyp='DIIS',maxit=500,orbital=2100.2,accu=10,energy=1e-12}</div><div dir="ltr"><br></div></div></div>
</div>
</blockquote></div><br></div>