Skip to Content.
Sympa Menu

forum - Re: [abinit-forum] Why the Etotal has the parabola-function-relationship in respect to the increasment of the parameter *acell*?

forum@abinit.org

Subject: The ABINIT Users Mailing List ( CLOSED )

List archive

Re: [abinit-forum] Why the Etotal has the parabola-function-relationship in respect to the increasment of the parameter *acell*?


Chronological Thread 
  • From: "Anglade Pierre-Matthieu" <anglade@gmail.com>
  • To: forum@abinit.org
  • Subject: Re: [abinit-forum] Why the Etotal has the parabola-function-relationship in respect to the increasment of the parameter *acell*?
  • Date: Mon, 16 Apr 2007 12:07:16 +0200
  • Domainkey-signature: a=rsa-sha1; c=nofws; d=gmail.com; s=beta; h=received:message-id:date:from:to:subject:in-reply-to:mime-version:content-type:content-transfer-encoding:content-disposition:references; b=m5XYocgtQxsUU2uuRpzn1RYkJdnylwiEdM4cYqOMESvPAJ2R2yuioFahxhrBrStm3RGhKi+v2YAfqfQuG7UufA2AWgnY4upv8EDSljYQQ64Bz9MbgpCUj2H9rGJkx7q7Ao5Yib1VUlpScIq8Ol4RZKipfIv1CoQf9n8UFZGRo2E=

Hi,
You've just dicovered the so called equilibrium lattice parameter.
Would the function be monotonic it would means that the structure
under consideration either collapse spontaneously or explode. Which is
quite not an appropriate behavior for neutral atoms. Right ?
regards
PMA


On 4/16/07, Hongyi Zhao <hongyi.zhao@gmail.com> wrote:
Hi all,

I am wondering Why the Etotal has the parabola-function-relationship in
respect to the increasment of the parameter *acell*? For detail of my
problem, I describe it as follows:

The input file for abinit is like this:

-----------------
ndtset 15
acell: 9.6 9.6 9.6
acell+ 0.1 0.1 0.1
kptopt 1
ngkpt 6 6 6
toldfe 1.0d-6
iscf 5
nstep 15
diemac 12.0
ecut 9.0
enunit 2
nband 5
ixc 1
typat 1 1
znucl 14
natom 2
rprim 0.0 0.5 0.5
0.5 0.0 0.5
0.5 0.5 0.0
xred 0.0 0.0 0.0
0.25 0.25 0.25
-----------------

Then, I find that the the Etotal has the parabola-function-relationship
in respect to the increasment of the parameter *acell*, in detail, the
Etotal will take its minimum vaule when the *acell* is about 10.2 bohr.
Dut to my poor theoretical background, I can't understand this
phenomena. Who can give me some hints on this thing?

Thanks in advance.

---
Hongyi Zhao
GnuPG DSA: 0xD108493B




--
Pierre-Matthieu Anglade



Archive powered by MHonArc 2.6.16.

Top of Page