<div style="line-height:1.7;color:#000000;font-size:14px;font-family:arial"> Dear QE Developer and Users<br>I am using DFPT+QHA to predict model's properties at high temperature. According to the paper, <i>Ab initio lattice dynamics and phase transformations of ZrO2, Phy Rev B 71, 064301</i>, free energy <i>v.s.</i> volume figure can be plotted at certain temperature for EOS-fitting and so on.<br>Having determined the equlibrium lattice parameter, 8.18 a.u. , of TiC(FCC) by both EOS-fitting and vc-relax, I chose 20 different lattice paramters (7.98~8.38a.u.) and calculated phonon dispersion and thermodynamic properties. However, the free energy kept decreasing as a monotone decreasing function of lattice parameter, even when imaginary frequency began to occur at 8.88a.u..<br>My DFPT calculation script shows below:<br><i></i><span style="color: rgb(0, 0, 255);"><i><span style="color: rgb(0, 0, 255);">cat >TiC_${INT}.scf.in <<EOF</span><br><span style="color: rgb(0, 0, 255);">&control</span><br><span style="color: rgb(0, 0, 255);">     calculation = 'scf' ,</span><br><span style="color: rgb(0, 0, 255);">     prefix = 'TiC_${INT}' ,</span><br><span style="color: rgb(0, 0, 255);">     outdir = '$TempDIR/' ,</span><br><span style="color: rgb(0, 0, 255);">     pseudo_dir = '$PseudoDIR/' ,</span><br><span style="color: rgb(0, 0, 255);">     disk_io = 'low' ,</span><br><span style="color: rgb(0, 0, 255);">     etot_conv_thr = 1.0d-9 ,</span><br><span style="color: rgb(0, 0, 255);">     tstress = .true. ,</span><br><span style="color: rgb(0, 0, 255);">/</span><br><span style="color: rgb(0, 0, 255);">&system</span><br><span style="color: rgb(0, 0, 255);">     ibrav = 2 ,</span><br><span style="color: rgb(0, 0, 255);">     celldm(1) = $OP ,</span><br><span style="color: rgb(0, 0, 255);">     ecutwfc = 42.5 ,</span><br><span style="color: rgb(0, 0, 255);">     ecutrho = 400 ,</span><br><span style="color: rgb(0, 0, 255);">     nat = 2 ,</span><br><span style="color: rgb(0, 0, 255);">     ntyp = 2 ,</span><br><span style="color: rgb(0, 0, 255);">     occupations = 'smearing' ,</span><br><span style="color: rgb(0, 0, 255);">     smearing = 'mp' ,</span><br><span style="color: rgb(0, 0, 255);">     degauss = 0.02 ,</span><br><span style="color: rgb(0, 0, 255);">/</span><br><span style="color: rgb(0, 0, 255);">&electrons</span><br><span style="color: rgb(0, 0, 255);">     conv_thr = 1.0d-9 ,</span><br><span style="color: rgb(0, 0, 255);">     diagonalization = 'cg' ,</span><br><span style="color: rgb(0, 0, 255);">     startingpot = 'atomic' ,</span><br><span style="color: rgb(0, 0, 255);">     startingwfc = 'atomic+random' ,</span><br><span style="color: rgb(0, 0, 255);">     mixing_mode = 'plain' ,</span><br><span style="color: rgb(0, 0, 255);">/</span><br><span style="color: rgb(0, 0, 255);">ATOMIC_SPECIES</span><br><span style="color: rgb(0, 0, 255);">  Ti  47.900  Ti.USPP_PW91.UPF</span><br><span style="color: rgb(0, 0, 255);">   C  12.011  C.USPP_PW91n.UPF</span><br><span style="color: rgb(0, 0, 255);">ATOMIC_POSITIONS crystal</span><br><span style="color: rgb(0, 0, 255);">  Ti  0.00  0.00  0.00</span><br><span style="color: rgb(0, 0, 255);">   C  0.50  0.50  0.50</span><br><span style="color: rgb(0, 0, 255);">K_POINTS automatic</span><br><span style="color: rgb(0, 0, 255);">  11  11  11   0  0  0</span><br><span style="color: rgb(0, 0, 255);">EOF</span><br><span style="color: rgb(0, 0, 255);">mpirun -n 6 pw.x -npool 2 <TiC_${INT}.scf.in >TiC_${INT}.scf.out</span><br><br><span style="color: rgb(0, 0, 255);">cat >TiC_${INT}.ph.in <<EOF</span><br><span style="color: rgb(0, 0, 255);">Phonon of TiC at $OP Bohr</span><br><span style="color: rgb(0, 0, 255);">&inputph</span><br><span style="color: rgb(0, 0, 255);">     outdir = '$TempDIR/' ,</span><br><span style="color: rgb(0, 0, 255);">     prefix = 'TiC_${INT}' ,</span><br><span style="color: rgb(0, 0, 255);">     nq1 = 3 ,</span><br><span style="color: rgb(0, 0, 255);">     nq2 = 3 ,</span><br><span style="color: rgb(0, 0, 255);">     nq3 = 3 ,</span><span style="color: rgb(0, 0, 255);"></span><br><span style="color: rgb(0, 0, 255);">     reduce_io = .true. ,</span><br><span style="color: rgb(0, 0, 255);">     ldisp = .true. ,</span><br><span style="color: rgb(0, 0, 255);">     tr2_ph = 1.0d-13 ,</span><br><span style="color: rgb(0, 0, 255);">     fildyn = 'TiC_${INT}.dyn' ,</span><br><span style="color: rgb(0, 0, 255);">/</span><br><span style="color: rgb(0, 0, 255);">EOF</span><br><span style="color: rgb(0, 0, 255);">mpirun -n 6 ph.x -npool 2 <TiC_${INT}.ph.in >TiC_${INT}.ph.out</span></i></span><br>QHA is calculated by Prof. Isaev's code which released with QE package.<br>Any suggestion will be welcome. <br>Thanks.<br><br><br><div>--<br>GAO Zhe<br>CMC Lab, MSE, SNU, Seoul, S.Korea<br>        
</div></div><br><br><span title="neteasefooter"><span id="netease_mail_footer"></span></span>