<div style="line-height:1.7;color:#000000;font-size:14px;font-family:arial">&nbsp;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 &gt;TiC_${INT}.scf.in &lt;&lt;EOF</span><br><span style="color: rgb(0, 0, 255);">&amp;control</span><br><span style="color: rgb(0, 0, 255);">&nbsp;&nbsp;&nbsp;&nbsp; calculation = 'scf' ,</span><br><span style="color: rgb(0, 0, 255);">&nbsp;&nbsp;&nbsp;&nbsp; prefix = 'TiC_${INT}' ,</span><br><span style="color: rgb(0, 0, 255);">&nbsp;&nbsp;&nbsp;&nbsp; outdir = '$TempDIR/' ,</span><br><span style="color: rgb(0, 0, 255);">&nbsp;&nbsp;&nbsp;&nbsp; pseudo_dir = '$PseudoDIR/' ,</span><br><span style="color: rgb(0, 0, 255);">&nbsp;&nbsp;&nbsp;&nbsp; disk_io = 'low' ,</span><br><span style="color: rgb(0, 0, 255);">&nbsp;&nbsp;&nbsp;&nbsp; etot_conv_thr = 1.0d-9 ,</span><br><span style="color: rgb(0, 0, 255);">&nbsp;&nbsp;&nbsp;&nbsp; tstress = .true. ,</span><br><span style="color: rgb(0, 0, 255);">/</span><br><span style="color: rgb(0, 0, 255);">&amp;system</span><br><span style="color: rgb(0, 0, 255);">&nbsp;&nbsp;&nbsp;&nbsp; ibrav = 2 ,</span><br><span style="color: rgb(0, 0, 255);">&nbsp;&nbsp;&nbsp;&nbsp; celldm(1) = $OP ,</span><br><span style="color: rgb(0, 0, 255);">&nbsp;&nbsp;&nbsp;&nbsp; ecutwfc = 42.5 ,</span><br><span style="color: rgb(0, 0, 255);">&nbsp;&nbsp;&nbsp;&nbsp; ecutrho = 400 ,</span><br><span style="color: rgb(0, 0, 255);">&nbsp;&nbsp;&nbsp;&nbsp; nat = 2 ,</span><br><span style="color: rgb(0, 0, 255);">&nbsp;&nbsp;&nbsp;&nbsp; ntyp = 2 ,</span><br><span style="color: rgb(0, 0, 255);">&nbsp;&nbsp;&nbsp;&nbsp; occupations = 'smearing' ,</span><br><span style="color: rgb(0, 0, 255);">&nbsp;&nbsp;&nbsp;&nbsp; smearing = 'mp' ,</span><br><span style="color: rgb(0, 0, 255);">&nbsp;&nbsp;&nbsp;&nbsp; degauss = 0.02 ,</span><br><span style="color: rgb(0, 0, 255);">/</span><br><span style="color: rgb(0, 0, 255);">&amp;electrons</span><br><span style="color: rgb(0, 0, 255);">&nbsp;&nbsp;&nbsp;&nbsp; conv_thr = 1.0d-9 ,</span><br><span style="color: rgb(0, 0, 255);">&nbsp;&nbsp;&nbsp;&nbsp; diagonalization = 'cg' ,</span><br><span style="color: rgb(0, 0, 255);">&nbsp;&nbsp;&nbsp;&nbsp; startingpot = 'atomic' ,</span><br><span style="color: rgb(0, 0, 255);">&nbsp;&nbsp;&nbsp;&nbsp; startingwfc = 'atomic+random' ,</span><br><span style="color: rgb(0, 0, 255);">&nbsp;&nbsp;&nbsp;&nbsp; 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);">&nbsp; Ti&nbsp; 47.900&nbsp; Ti.USPP_PW91.UPF</span><br><span style="color: rgb(0, 0, 255);">&nbsp;&nbsp; C&nbsp; 12.011&nbsp; 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);">&nbsp; Ti&nbsp; 0.00&nbsp; 0.00&nbsp; 0.00</span><br><span style="color: rgb(0, 0, 255);">&nbsp;&nbsp; C&nbsp; 0.50&nbsp; 0.50&nbsp; 0.50</span><br><span style="color: rgb(0, 0, 255);">K_POINTS automatic</span><br><span style="color: rgb(0, 0, 255);">&nbsp; 11&nbsp; 11&nbsp; 11&nbsp;&nbsp; 0&nbsp; 0&nbsp; 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 &lt;TiC_${INT}.scf.in &gt;TiC_${INT}.scf.out</span><br><br><span style="color: rgb(0, 0, 255);">cat &gt;TiC_${INT}.ph.in &lt;&lt;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);">&amp;inputph</span><br><span style="color: rgb(0, 0, 255);">&nbsp;&nbsp;&nbsp;&nbsp; outdir = '$TempDIR/' ,</span><br><span style="color: rgb(0, 0, 255);">&nbsp;&nbsp;&nbsp;&nbsp; prefix = 'TiC_${INT}' ,</span><br><span style="color: rgb(0, 0, 255);">&nbsp;&nbsp;&nbsp;&nbsp; nq1 = 3 ,</span><br><span style="color: rgb(0, 0, 255);">&nbsp;&nbsp;&nbsp;&nbsp; nq2 = 3 ,</span><br><span style="color: rgb(0, 0, 255);">&nbsp;&nbsp;&nbsp;&nbsp; nq3 = 3 ,</span><span style="color: rgb(0, 0, 255);"></span><br><span style="color: rgb(0, 0, 255);">&nbsp;&nbsp;&nbsp;&nbsp; reduce_io = .true. ,</span><br><span style="color: rgb(0, 0, 255);">&nbsp;&nbsp;&nbsp;&nbsp; ldisp = .true. ,</span><br><span style="color: rgb(0, 0, 255);">&nbsp;&nbsp;&nbsp;&nbsp; tr2_ph = 1.0d-13 ,</span><br><span style="color: rgb(0, 0, 255);">&nbsp;&nbsp;&nbsp;&nbsp; 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 &lt;TiC_${INT}.ph.in &gt;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>