<div dir="ltr">Hi Nic, <br><br>1)  I have tested the pseudo-potentials (psps) with small atomic separations. The pseudo-potentials give reasonable till two atoms are as close as 0.8 A. When the total energy starts to diverge, which is about Ionic iteration 17, the smallest inter-atomic distance is about 1.4 A. So I do not think the breakdown of the psps is the trigger of this divergence of the total energy, although it should be the reason for the enormous large negative energies. <br>
<br>2) I have tried to restart the calculation from where the divergence starts. Once I did that, the first scf calculation would give a reasonable energy around -41.7 Ry, and the calculation could run smoothly till the total energy converges most of the times. So what I deal with this problem right now is to restart calculations from where they start to have diverged total energies and it's been working well. However, I hope this is not the optimal way to solve this problem. :)<br>
<br>Thank you very much,<br>Hui<br><br><div class="gmail_quote">On Thu, Oct 16, 2008 at 3:09 PM, Nicola Marzari <span dir="ltr"><<a href="mailto:marzari@mit.edu" target="_blank">marzari@mit.edu</a>></span> wrote:<br>

<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;"><div>Hui Tang wrote:<br>
> Dear all,<br>
><br>
> I'm trying to do cell relaxation on my Mg-Boron systems. For some of the<br>
> structures, I got extremely strange results. The total energy diverges<br>
> starting from some ionic step. For example, here are the total energy in<br>
> each ionic step for one of the problematic run:<br>
><br>
<br>
</div>Two suggestions<br>
<br>
1) are atoms overlapping much more than they should ? I.e. are these<br>
new configurations reasonable, or are core overlaps very significant ?<br>
The pseudopotential approximation breaks down if cores overlap (a lot).<br>
<br>
2) is it a bug/memory overrun ? Take the first configuration with<br>
much lower energy, and rerun that from scratch. Do you get an<br>
energy that looks like -41.7 Ry, or much lower ?<br>
<br>
                        nic<br>
<br>
<br>
<br>
---------------------------------------------------------------------<br>
Prof Nicola Marzari   Department of Materials Science and Engineering<br>
13-5066   MIT   77 Massachusetts Avenue   Cambridge MA 02139-4307 USA<br>
tel 617.4522758 fax 2586534 <a href="mailto:marzari@mit.edu" target="_blank">marzari@mit.edu</a> <a href="http://quasiamore.mit.edu" target="_blank">http://quasiamore.mit.edu</a><br>
_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org" target="_blank">Pw_forum@pwscf.org</a><br>
<a href="http://www.democritos.it/mailman/listinfo/pw_forum" target="_blank">http://www.democritos.it/mailman/listinfo/pw_forum</a><br>
</blockquote></div><br></div>