<div dir="ltr">So is it still possible to use OT with ADDED_MOS or must multiplicity be fixed. I am modeling a molecule interacting with a surface, I wanted to let the spin state be allowed to change during the interaction if needed. I don't want to give up on OT since it is so fast compared to diagonalization for my large system.<br><br>On Wednesday, August 12, 2015 at 10:45:28 AM UTC-6, Matt W wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr">I can confirm what you say, but I'm pretty sure it is incorrect behaviour of cp2k to run with OT and relax_multiplicity.  <div><br></div><div>The relax_multiplicity flag must be getting ignored when the code does an atomic guess.<div><br></div><div>Matt<br><br>On Wednesday, August 12, 2015 at 4:51:07 PM UTC+1, Christopher O'Brien wrote:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Matt,<div>Thanks for you response.  Although I don't exactly understand how the OT minimizer works, there is no error generated by CP2K or any mention of a conflict between changing multiplicity in the documentation for the OT minimizer.  </div><div><br></div><div>However, I don't think there is a conflict with using OT and varying multiplicity because using the same input, without the WFN file used as an initial guess, runs without error.   <br><br>On Wednesday, August 12, 2015 at 12:58:06 AM UTC-6, Matt W wrote:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi,<div><br></div><div>I think the problem is coming from the relax multiplicity keyword in your input - that only seems to make sense to me using extra MOs and not the OT optimizer as it requires having extra MOs it can put electrons into to vary the multiplicity.</div><div><br></div><div>Matt<br><br>On Tuesday, August 11, 2015 at 8:05:28 PM UTC+1, Christopher O'Brien wrote:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">OT calculations cannot be restarted in CP2K.<div><br></div><div>The error listed below results from using an existing WFN file and the appropriate coordinates from another CP2K MD run to start a new METADYNAMICS run. In my DFT section I have:</div><div><div style="background-color:rgb(250,250,250);border:1px solid rgb(187,187,187);word-wrap:break-word"><code><div><span style="color:#660">&</span><span style="color:#000">FORCE_EVAL<br><br> METHOD QS <br></span><span style="color:#660"> &</span><span style="color:#000">DFT <br><br>    BASIS_SET_FILE_NAME BASIS_MOLOPT <br>    POTENTIAL_FILE_NAME GTH_POTENTIALS <br>    WFN_RESTART_FILE_NAME XXX</span><span style="color:#660">-</span><span style="color:#000">RESTART</span><span style="color:#660">.</span><span style="color:#000">wfn <br>    UKS<br>    RELAX_MULTIPLICITY </span><span style="color:#066">0.15</span><span style="color:#000"><br><br></span><span style="color:#000"><span style="color:rgb(102,102,0)">  &</span><span>SCF</span><br><br>    </span><span style="color:#660">&</span><span style="color:#000">OT <br><br>      PRECONDITIONER FULL_SINGLE_INVERSE<br>      MINIMIZER CG</span></div><div><span style="color:#000"><br>    </span><span style="color:#660">&</span><span style="color:#008">END</span><span style="color:#000"> OT<br></span><span style="color:#000"><br>    EPS_SCF </span><span style="color:#066">1.0E-6</span><span style="color:#000"><br>    SCF_GUESS RESTART<br>    MAX_SCF </span><span style="color:#066">100</span><span style="color:#000"><br><br>  </span><span style="color:#660">&</span><span style="color:#008">END</span><span style="color:#000"> SCF<br><br><span> </span><span style="color:rgb(102,102,0)">&</span><span style="color:rgb(0,0,136)">END</span><span> DFT<br></span><br> &SUBSYS ...<br><br></span><font color="#666600"><span style="color:#660">&</span><span style="color:#008">END</span><span style="color:#000"> FORCE_EVAL</span></font><span style="color:#000"><br><br></span></div></code></div><div>However, after searching the manual I can't find an option to add MOs to OT calculations. I also don't understand why I need to considering that, in the output of the calculation I am starting from, all spins are already paired. <br></div><div><br></div><div style="background-color:rgb(250,250,250);border:1px solid rgb(187,187,187);word-wrap:break-word"><code><div><span style="color:#000"> </span><span style="color:#606">Spin</span><span style="color:#000"> </span><span style="color:#066">1</span><span style="color:#000"> <br><br><br> </span><span style="color:#606">Number</span><span style="color:#000"> of electrons</span><span style="color:#660">:</span><span style="color:#000"> </span><span style="color:#066">1993</span><span style="color:#000"> <br><br> </span><span style="color:#606">Number</span><span style="color:#000"> of occupied orbitals</span><span style="color:#660">:</span><span style="color:#000"> </span><span style="color:#066">1993</span><span style="color:#000"> <br><br> </span><span style="color:#606">Number</span><span style="color:#000"> of molecular orbitals</span><span style="color:#660">:</span><span style="color:#000"> </span><span style="color:#066">1993</span><span style="color:#000"> <br><br><br> </span><span style="color:#606">Spin</span><span style="color:#000"> </span><span style="color:#066">2</span><span style="color:#000"> <br><br><br> </span><span style="color:#606">Number</span><span style="color:#000"> of electrons</span><span style="color:#660">:</span><span style="color:#000"> </span><span style="color:#066">1993</span><span style="color:#000"> <br><br> </span><span style="color:#606">Number</span><span style="color:#000"> of occupied orbitals</span><span style="color:#660">:</span><span style="color:#000"> </span><span style="color:#066">1993</span><span style="color:#000"> <br><br> </span><span style="color:#606">Number</span><span style="color:#000"> of molecular orbitals</span><span style="color:#660">:</span><span style="color:#000"> </span><span style="color:#066">1993</span><span style="color:#000"> <br> <br><br> </span><span style="color:#606">Number</span><span style="color:#000"> of orbital functions</span><span style="color:#660">:</span><span style="color:#000"> </span><span style="color:#066">15770</span><span style="color:#000"> <br><br> </span><span style="color:#606">Number</span><span style="color:#000"> of independent orbital functions</span><span style="color:#660">:</span><span style="color:#000"> </span><span style="color:#066">15770</span><span style="color:#000"> <br><br><br> </span><span style="color:#606">Extrapolation</span><span style="color:#000"> method</span><span style="color:#660">:</span><span style="color:#000"> initial_guess <br><br><br> </span><span style="color:#660">***</span><span style="color:#000"> </span><span style="color:#066">04</span><span style="color:#660">:</span><span style="color:#066">08</span><span style="color:#660">:</span><span style="color:#066">41</span><span style="color:#000"> WARNING </span><span style="color:#008">in</span><span style="color:#000"> qs_mo_types</span><span style="color:#660">:</span><span style="color:#000">set_mo_occupation_<wbr>2 </span><span style="color:#660">::</span><span style="color:#000"> </span><span style="color:#606">All</span><span style="color:#000"> alpha </span><span style="color:#606">MOs</span><span style="color:#000"> </span><span style="color:#660">***</span><span style="color:#000"> <br><br> </span><span style="color:#660">***</span><span style="color:#000"> are occupied</span><span style="color:#660">.</span><span style="color:#000"> </span><span style="color:#606">Add</span><span style="color:#000"> more alpha </span><span style="color:#606">MOs</span><span style="color:#000"> to allow </span><span style="color:#008">for</span><span style="color:#000"> a higher multiplicity </span><span style="color:#660">***</span><span style="color:#000"><br><br><br></span></div></code></div></div></div></blockquote></div></div></blockquote></div></div></blockquote></div></div></div></blockquote></div>