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 users@lists.quantum-espresso.org https://lists.quantum-espresso.org/mailman/listinfo/users