
<html><body><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></body></html>