Skip to Content.
Sympa Menu

forum - Re: [abinit-forum] cell/ion relaxation results

forum@abinit.org

Subject: The ABINIT Users Mailing List ( CLOSED )

List archive

Re: [abinit-forum] cell/ion relaxation results


Chronological Thread 
  • From: "matthieu verstraete" <matthieu.jean.verstraete@gmail.com>
  • To: forum@abinit.org
  • Subject: Re: [abinit-forum] cell/ion relaxation results
  • Date: Mon, 19 Jan 2009 11:19:10 +0100
  • Domainkey-signature: a=rsa-sha1; c=nofws; d=gmail.com; s=gamma; h=message-id:date:from:to:subject:in-reply-to:mime-version :content-type:references; b=rmVphmQc5eKHOnAJ+GSeqYzIRJi7a6cxch0FfCVGvFwFsSUOrF84X4KJWRWywr9/d+ D0HfWzQ7g2rK91iTLbBCylED8S4Ma+q60Sn/Lb9X5ybn5Ubzzw3Od1crF/1dPJMjyUpL 6QQAxU50wOBPV28nfFRWspHZ92wjO2AjHuVjE=

As Pierre-Matthieu said, we need more information to be accurate, but a few months ago something like this was discussed on the forum: it is probable that you have changed the ecutsm or dilatmx parameters between your 2 runs. This changes
1) the actual spheres of G vectors used
2) the contribution of the last G vectors (with ecutsm their contribution goes smoothly to 0 for the last shell)

both of these will affect the energy and stress calculated. If you are well converged then the contribution from the last shells of G should be negligible, but usually this is not the case. Try to repeat your 2nd calculation with the exact same ecutsm and dilatmx (and ecut of course) to see if it improves agreement. Also, you should try to exhaust the convergence, to see if you are close to the correct value.

Matthieu

On Mon, Jan 19, 2009 at 8:27 AM, Anglade Pierre-Matthieu <anglade@gmail.com> wrote:
Hi,

It is quite likely possible to help you. However, to proceed, one would need : both of your input file and relevant part of the relaxation output file.
Do you consider possible that your calculation parameters were unconverged ? Or that your scf/non-scf loops were unconverged ?

regards

PMA


On Mon, Jan 19, 2009 at 6:38 AM, kong <kongl@physics.rutgers.edu> wrote:
Hi,

I relaxed a system using optcell =2 and ionmov=3. After the desired convergence
is achieved, I copied the relaxed lattice parameters( acell/angles) and
atomic positions(xangst)  from the output file and start a new calculation
without
relaxation. All the other input variables are same as before.  I am very
surprised
that the new total energy is about 0.2Ha lower than  the relaxed results,
but the stress tensor is three order of magnitude larger!

I note that ngfft/boxcut(ratio) in the relaxation is different. But using the
same ngfft(which also leads to the same boxcut(ratio)) in the input file does
not improve the consistency.
Would anyone help me out on it? Thanks!

Lingzhu



--
Pierre-Matthieu Anglade




--
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
Dr. Matthieu Verstraete

European Theoretical Spectroscopy Facility (ETSF) 
Dpto. Fisica de Materiales, 
U. del Pais Vasco,            
Centro Joxe Mari Korta, Av. de Tolosa, 72,   Phone: +34-943018393
E-20018 Donostia-San Sebastian, Spain        Fax  : +34-943018390

Mail : matthieu.jean.verstraete@gmail.com
http://www-users.york.ac.uk/~mjv500



Archive powered by MHonArc 2.6.15.

Top of Page