<br><br>On Thursday, March 22, 2012 1:17:07 PM UTC, Samuel Lamphier wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;">Thank you for your response, I used crystal maker to derive a 2x2x2
<br> supercell, I exported the orthogonal coordinates directly. However, I
<br> thought that if I incorporated a vacancy that I would see a much
<br> larger relaxation as I do using classical simulations. The change was
<br> on the order of 0.000001 for cp2k and during classical simulations, I
<br> observed changes of 0.1 to 1.5A during relaxation.
<br>
</blockquote><div> </div><div>Sorry, didn't notice that you had included an oxygen vacancy. In that case without seeing more of your output it is impossible to say what occurs.</div><div><br></div><div>Does the wavefunction converge?  Does CP2K report geometry optimization steps in the output? Perhaps you should quickly check you can relax a water molecule or similar just to be sure the problem is with this input.</div><div><br></div><div>BTW, I suspect that the Ce basis set and pseudo potential have never been used in the wild, so be very careful. </div><div><br></div><div>PS Does anyone fancy optimizing a GGA pseudo for Ce?</div>