forum@abinit.org
Subject: The ABINIT Users Mailing List ( CLOSED )
List archive
- From: "matthieu verstraete" <matthieu.jean.verstraete@gmail.com>
- To: forum@abinit.org
- Subject: Re: [abinit-forum] natfix, iatfix
- Date: Sat, 10 Jan 2009 12:26:12 +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=o3Q1KECXZjAFSMlqv4dsoAy5KyU7Joyux9li/mhLoWEW07k8IsIBHrIX8GvZOpPbeJ u/zEnWr/6rCW9CCor4Vmow+8Nf2T2DpT+2KZbjfMwQF5yfC416+Pci+DOauBYKUE4o70 BKpw8DmFUjmS9qGaEOnV6/iuvFcLlhOg8kjVQ=
On Sat, Jan 10, 2009 at 7:02 AM, Julie Smart <smart.julie@rocketmail.com> wrote:
Hi all,
Hello
Working on bulk GaAs optimization, I have imposed natfix and iatfix inorder to fix "Ga" position. But after the optimization I have the output xred like:
0.3333333333333 0.6666666666667 1.0658143578E-15
0.6666666666667 0.3333333333333 0.500000000000
0.3333333333333 0.6666666666667 0.374956338791
0.6666666666667 0.3333333333333 0.874956338791
you mean the 1.e-15? That is machine precision: simple numerical roundoff errors will always give you this kind of thing. You can reset it to 0 if you like but it won't change anything. You can also see that the last 2 atoms are precisely 0.5 apart in z, down to 12 digits - this is amply enough precision. Finally, the details of the numerical precision are down to which platform and compiler you use, so if it's not the same as the reference machine for tnlo you can get these tiny differences. Abinit is one of the most portable and precise codes around, for these matters.
For symmetry purposes, abinit examines down to 1.e-10 in reduced coordinates.
Further, by symmetry, you do not need to fix the Ga here. However if you didn't you might end up with a shift of both atoms, e.g. z_Ga=0.2 and z_As=0.2+0.374956338791, which is physically the same but the numbers are not as aesthetic.
For symmetry purposes, abinit examines down to 1.e-10 in reduced coordinates.
Further, by symmetry, you do not need to fix the Ga here. However if you didn't you might end up with a shift of both atoms, e.g. z_Ga=0.2 and z_As=0.2+0.374956338791, which is physically the same but the numbers are not as aesthetic.
Why do I have the slight displace in the first atom z parameter (following exactly tnlo-1.in) ? Also, please let me know how we can set "strfact" in the optimization in order to get to the level of zero stress in optimization?
strfact is 0 by default, which relaxes to 0 stress (see the definition of strfact for units and components).
Matthieu
Matthieu
--
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
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
- [abinit-forum] natfix, iatfix, Julie Smart, 01/10/2009
- Re: [abinit-forum] natfix, iatfix, matthieu verstraete, 01/10/2009
Archive powered by MHonArc 2.6.15.