(Hopefully the test input will come through this time...)<br><br>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>memory,8,m<br>gthresh,energy=1d-12,twoint=1d-18,prefac=1d-20<br>geomtyp=xyz<br>geom={<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>
!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></div>