forum@abinit.org
Subject: The ABINIT Users Mailing List ( CLOSED )
List archive
- From: "matthieu verstraete" <matthieu.jean.verstraete@gmail.com>
- To: "Petr Šesták" <sestak@kn.vutbr.cz>
- Cc: pganesh@ciw.edu, forum@abinit.org
- Subject: Re: Re: [abinit-forum] constrained optimization
- Date: Wed, 7 Jan 2009 18:00:43 +0100
- Domainkey-signature: a=rsa-sha1; c=nofws; d=gmail.com; s=gamma; h=message-id:date:from:to:subject:cc:in-reply-to:mime-version :content-type:references; b=KOjeWMRODHHTCA5ysAFOY/4LCNr/29z+/ejCAlEJ3cisZh3Zf2SIdbXxgY5AGsRuDN HhkYQGxDPjBN9HqhnlWS+6+/V+EEA+F43OX8JCaMpeImigGSC/v6wUGalZRvpIWkTcgf FgVNqmD+7B9plhyaa7gdF2ryVF/17RjSXwtRk=
Thanks Petr, this is great. I hope someone will be able to include it as a new form of optcell in abinit (so we don't have to deal with scripting). Is this ok with you? And have you considered doing it yourself?
I have mainly looked at the pdf summary so far. It appears quite similar to the standard way of updating the rprim, apart perhaps from the compliance tensor bit. Is this what you fix in order to keep certain axes constant and not others in the monoclinic case? Otherwise, how do you determine which vectors you want to keep fixed and which will relax? I didn't find a line which zeros certain elements of the delta stress tensor (probably all but the 33 one).
Your i,j indices are in cartesian coordinates?
cheers, and Happy New Year
Matthieu
On Wed, Jan 7, 2009 at 3:47 PM, Petr Šesták <sestak@kn.vutbr.cz> wrote:
Hi Matthieu and P. Ganesh,
In the attachement I am sending you my relaxation script (Perl) along with my input files (for better understanding) and the schematic diagram of the entire relaxation cycle. I apologize for the delay in my response (caused by the Chrismas). If you need any further clarification do not hesitate to contact me.
If I clearly understood the PGanesh question why did I not make a relaxation of ionic position, so the anwser is as you have written - the ionic positions were given by the crystal symetry. I had only two atoms in the basis. Now I am doing some computations on a more complicated structure and I am planing to extend this script by a posibillity of full relaxation (by alternating the translating vectors relaxations and the relaxations of ionic positions).
To the second question "about changing the RPRIM according to stress tensor from previous computation": This procedure can be seen in the script and the increments of RPRIM values are also affected by the convergence parameter (see the array @du).
Best regards
Ing. Peter Sestak
Institute of Physical Engineering
Faculty of Mechanical Engineering
Brno University of Technology
abinit_config.in
abinit_config.out
abinit_configi
abinit_configo
abinit_config
22-Ti.GGA.fhi
28-Ni.GGA.fhi
--
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
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
- Re: Re: [abinit-forum] constrained optimization, Petr Šesták, 01/07/2009
- Re: Re: [abinit-forum] constrained optimization, matthieu verstraete, 01/07/2009
- Re: Re: Re: [abinit-forum] constrained optimization, Petr Sestak, 01/08/2009
- Re: [abinit-forum] constrained optimization, PGanesh, 01/08/2009
- Re: Re: [abinit-forum] constrained optimization, Petr Sestak, 01/12/2009
- Re: [abinit-forum] constrained optimization, PGanesh, 01/08/2009
- Re: Re: Re: [abinit-forum] constrained optimization, Petr Sestak, 01/08/2009
- Re: Re: [abinit-forum] constrained optimization, matthieu verstraete, 01/07/2009
Archive powered by MHonArc 2.6.15.