Hi users,
I am running Al16 N16 scf (pw.x) and phG (ph.x) calculations for 32 atoms. Both
scf and phG calculations run smoothly. I need to get d33 (piezoelectric
response), e33 (stress constant) and c33 (elastic constant). I know these are
related by (d33=e33/c33). Can you please suggest a simple way to get to d33
using any input parameters in phG perhaps?
Please find both scf and phG files below.
Best,
Moh
Research Student
University of New Mexico
pw.x scf.in:
&CONTROL
calculation = 'scf'
etot_conv_thr = 3.2000000000d-04
forc_conv_thr = 1.0000000000d-04
outdir = './out/'
prefix = 'aiida'
pseudo_dir = '.'
tprnfor = .true.
tstress = .true.
verbosity = 'high'
/
&SYSTEM
degauss = 1.4699723600d-02
ecutrho = 3.5279352913d+01
ecutwfc = 4.4099191141d+00
ibrav = 0
nat = 32
nosym = .false.
ntyp = 2
occupations = 'smearing'
smearing = 'cold'
/
&ELECTRONS
conv_thr = 6.4000000000d-09
electron_maxstep = 80
mixing_beta = 4.0000000000d-01
/
ATOMIC_SPECIES
Al 26.981538 Al.pbesol-n-kjpaw_psl.1.0.0.UPF
N 14.0067 N.pbesol-theos.UPF
ATOMIC_POSITIONS crystal
Al 0.1666666720 0.3333333430 0.4996435050
Al 0.1666666720 0.3333333430 0.9996435050
Al 0.1666666720 0.8333333730 0.4996435050
Al 0.1666666720 0.8333333730 0.9996435050
Al 0.6666666870 0.3333333430 0.4996435050
Al 0.6666666870 0.3333333430 0.9996435050
Al 0.6666666870 0.8333333730 0.4996435050
Al 0.6666666870 0.8333333730 0.9996435050
Al 0.3333333430 0.1666666720 0.2496435050
Al 0.3333333430 0.1666666720 0.7496435050
Al 0.3333333430 0.6666666870 0.2496435050
Al 0.3333333430 0.6666666870 0.7496435050
Al 0.8333333730 0.1666666720 0.2496435050
Al 0.8333333730 0.1666666720 0.7496435050
Al 0.8333333730 0.6666666870 0.2496435050
Al 0.8333333730 0.6666666870 0.7496435050
N 0.1666666720 0.3333333430 0.1903564930
N 0.1666666720 0.3333333430 0.6903564930
N 0.1666666720 0.8333333730 0.1903564930
N 0.1666666720 0.8333333730 0.6903564930
N 0.6666666870 0.3333333430 0.1903564930
N 0.6666666870 0.3333333430 0.6903564930
N 0.6666666870 0.8333333730 0.1903564930
N 0.6666666870 0.8333333730 0.6903564930
N 0.3333333430 0.1666666720 0.4403564930
N 0.3333333430 0.1666666720 0.9403564930
N 0.3333333430 0.6666666870 0.4403564930
N 0.3333333430 0.6666666870 0.9403564930
N 0.8333333730 0.1666666720 0.4403564930
N 0.8333333730 0.1666666720 0.9403564930
N 0.8333333730 0.6666666870 0.4403564930
N 0.8333333730 0.6666666870 0.9403564930
K_POINTS automatic
6 6 4 0 0 0
CELL_PARAMETERS angstrom
6.2571763992 0.0000000000 0.0000000000
-3.1285881996 5.4188737177 0.0000000000
0.0000000000 0.0000000000 10.0339097977
ph.x phG.in:
phonons of AlN at Gamma
&inputph
tr2_ph=1.0d-14, ! threshold for self consistency
prefix='aiida', ! same as in pw.x
epsil=.false., ! set true for dielectric constant
amass(1)=26.98,
amass(2)= 14.00,
search_sym = .false.,
fildyn='aln.dynG', ! is an output file not an input
outdir='./out/', ! same as in pw.x
/
0.0 0.0 0.0
_______________________________________________
Quantum ESPRESSO is supported by MaX (www.max-centre.eu)
users mailing list [email protected]
https://lists.quantum-espresso.org/mailman/listinfo/users
