<div>Hi,<br></div><div><br></div><div>I think ROKS should work with diagonalization as well as OT. There's a short thread "ROKS error" where Matthias implies that at least high spin should be working. Whether the eigenvalues of the orbitals are really properly defined I don't know.</div><div><br></div><div>Removing &OT and adding an &DIAGONALIZATION section I get energies and occupations printed from &MO</div><div><br></div><div>Cheers,</div><div><br></div><div>Matt </div><br>On Friday, May 25, 2012 9:24:01 AM UTC+1, Peter Tentscher wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;">Dear Matt,
<br>
<br>thanks for that information.
<br>Things are a bit more complicated than I thought.
<br>Ultimately, what I want is the SOMO energy from
<br>a ROKS calculation. That is, I do need OT in order
<br>to converge to the right solution.
<br>I tried using MO_CUBES in conjunction with OT/ROKS,
<br>but what I get is
<br>"  Unclear how we define MOs in the restricted case ... skipping"
<br>so this doesn't solve my problem.
<br>I also tried to use the *-RESTART.wfn from an OT calculation
<br>as a guess for the corresponding ROKS calculation without
<br>OT. But the subsequent SCF procedure then leads away from
<br>the former OT solution, i.e. raising the energy.
<br>
<br>Is there any way to get orbital energies from a ROKS calculation
<br>using OT?
<br>
<br>Best,
<br>
<br>Peter
<br>
<br>On 05/23/2012 10:39 PM, Matt W wrote:
<br>> Hi,
<br>>
<br>> seems you're a victim of some not fully synchronized algorithms in the 
<br>> code.
<br>>
<br>> "Controls the printing of the molecular orbitals.Note that this is 
<br>> only functional with diagonalization based methods, in particular not 
<br>> with OT (see MO_CUBES)"
<br>>
<br>> Turn OT off. Or try pdos or MO_CUBES in the print section of DFT.
<br>>
<br>> Matt
<br>>
<br>
<br></blockquote>