how to relax cell under hydrostatic pressure by abinit?

Total energy, geometry optimization, DFT+U, spin....

Moderator: bguster

Locked
nicheal
Posts: 3
Joined: Wed Apr 07, 2010 5:32 pm

how to relax cell under hydrostatic pressure by abinit?

Post by nicheal » Tue Nov 16, 2010 9:12 am

hello,
I want to add hydrostatic pressure to a system, but i don't know how to solve it?
I mean I have no conception between "strtarget" and hydrostatic, also the "optcell" and "ionmov"

sincerely,
Lichun

User avatar
jzwanzig
Posts: 504
Joined: Mon Aug 17, 2009 9:25 am

Re: how to relax cell under hydrostatic pressure by abinit?

Post by jzwanzig » Tue Nov 16, 2010 11:44 am

Well, strtarget is the desired final stress generated by the cell, in Voigt notation. So the longitudinal components are the first three entries and the shear components are the second three. The units are atomic, that is, Hartree/Bohr^3. So 1 GPa = 3.40E-5 in atomic units. So if you want the final structure to be under an external hydrostatic pressure of 1 GPa, set strtarget to -3.4E-5 -3.4E-5 -3.4E-5 0.0 0.0 0.0. If you want a negative external pressure (tension) change the signs to positive (remember that strtarget refers to internal stress). In order for this variable to do anything you will have to be relaxing the cell with ionmov 2 and optcell 1 or optcell 2.
Josef W. Zwanziger
Professor, Department of Chemistry
Canada Research Chair in NMR Studies of Materials
Dalhousie University
Halifax, NS B3H 4J3 Canada
jzwanzig@gmail.com

Locked