Dear MOLPRO gurus:<br><br>Very recently I ran into a peculiar example of a transition state for which MOLPRO disagrees with every other code on the ROHF reference energy. See the simple input below. I tried every SCF convergence acceleration/coaxing trick in the book and a few that aren't, to no avail.<br>
<br>The other structures on the same PES do not exhibit the problem. Does anybody have a clue what I am overlooking?<br><br>(And yes, I know this is a highly multiref problem. Still, ROHF should give the same answer as Gaussian and ACES?)memory,8,m<br>
<br>Many thanks in advance!<br>Jan Martin<br><br><br>    4<br>E(UHF)=-262.355671204 au|Dpl=(0.5152,-0.4848,0) D||E(ROHF)=-262.338447756 au|Dpl=(0.9471,-1.9519,0) Debye<br>O    -0.832632    1.113558    0.000000<br>C     0.000000    0.281073    0.000000<br>
O    -0.367117   -1.293021    0.000000<br>O     1.199748   -0.031341    0.000000<br>}<br>basis=vdz<br>int<br>{uhf;wf,30,2,2;occ,13,3;closed,12,2}<br>! yields the correct energy<br>!UHF STATE 1.2 ENERGY               -262.355671266903<br>

!UHF STATE 1.2 DIPOLE MOMENT           0.20268671    -0.19074510     0.00000000<br>{hf;wf,30,2,2;occ,13,3;closed,12,2}<br>! yields a different energy, and dipole vector looks very different too <br>!RHF STATE 1.2 ENERGY               -262.331395031327<br>
!RHF STATE 1.2 DIPOLE MOMENT           0.10470623    -0.21621532     0.00000000<br><br><br>