[QE-users] How to get piezoelectric response?

[Ashas Quazi]

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