<html><head></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">Dear Anirban,<div><br></div><div>it sounds to me like in the end Gaussian is optimizing to a symmetry-broken solution.  In Molpro you could try reducing the symmetry from D2h to C2v and giving it the occupations that Gaussian ends up with.  Probably though you really don't want this since it's presumably just an artifact of using a UHF wavefunction in this case.</div><div><br></div><div>best regards,</div><div><br></div><div>Kirk</div><div><br></div><div><div><div>On Sep 17, 2011, at 11:47 AM, anirban mandal wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite">

<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<div><div style="color:#000; background-color:#fff; font-family:times new roman, new york, times, serif;font-size:12pt"><div style="RIGHT: auto"><span style="RIGHT: auto">Hi all,</span></div>
<div style="RIGHT: auto"><span style="RIGHT: auto"></span> </div>
<div style="RIGHT: auto"><span style="RIGHT: auto">I have a quick question about the calculations in Gaussian (98 and 03) using the condition "stable=opt". I did an UHF calculation for a Ru2 triplet using both Molpro 2010 and Gaussian (same basis set). First the energies didn't match, because they didn't converge to the same state. So I did dome orbital reordering in Gaussian to get the same state as obtained by Molpro, and the energies matched.</span></div>
<div style="RIGHT: auto"><span style="RIGHT: auto"></span> </div>
<div style="RIGHT: auto"><span style="RIGHT: auto">Now the problem is, there is an option in Gaussian called "stable=opt". If specified, it checks for internal and external instabilities of the wavefunction and if instabilities found, it again minimizes the energy. If I use that option, I do get lower energy than Molpro or Gamess everytime, but the orbitals now look really really strange. It looks like the orbitals get localized. For example, if I have two Ru atoms, the mo's have contributions from the s, p , and d orbitals of both Ru atoms. One the calculations are with stable=opt, some orbitals have contribution from one Ru atom, others have contribution from the other Ru atom. I don't know how to reproduce those orbitals in Molpro. I was wondering if someone had the same problem before, and how to actually handle it. Thank you guys.</span></div>
<div></div>
<div style="RIGHT: auto"> </div>
<div style="RIGHT: auto">Best regards,</div>
<div style="RIGHT: auto"> </div>
<div style="RIGHT: auto">Anirban.<var id="yui-ie-cursor"></var></div></div></div>_______________________________________________<br>Molpro-user mailing list<br><a href="mailto:Molpro-user@molpro.net">Molpro-user@molpro.net</a><br>http://www.molpro.net/mailman/listinfo/molpro-user<br></blockquote></div><br></div></body></html>