<div style="line-height:1.7;color:#000000;font-size:14px;font-family:arial"> Dear QE developers and users:<br>Hi.<br>I am using Quantum-ESPRESSO v4.2.1 to deal with DFPT and QHA. Because of the computational limitation of my PC, I chose VCA method to compare the phonon dispersion and thermo-properties change as a function of compositions.<br>The system I calculated was ZrC-ZrN. Since based on our experiments, we have known that the Zr(CN) solid-solution is possible to produce. PW91 xc-functional has been selected at 13x13x13 k-points mesh with ecutwfc=40Ry and ecutrho=400Ry. The q-mesh was 4x4x4. Structure of ZrC, ZrN and Zr(CN) is Fm-3m (NaCl Structure). The Vanderbilt ultrasoft pseudo-potential of C and N were combined by virtual.x in certain ratio.<br>However, I found the formation energy:<br><i></i><span style="color: rgb(0, 0, 255);"><i style="color: rgb(0, 0, 255);">Ff[T, Zr(CNx)] = F[T, Zr(CNx)] - (1-x) * F[T, ZrC] - x * F[T, ZrN]</i></span><br>was always positive number.<br>Then, I checked ZrC, ZrN and Zr(CN0.5) carefully by both of 'vc-relax' and E-V curve for total energy calculation. Finally, I found that <i style="color: rgb(0, 0, 255);">E[Zr(CN0.5)] - 0.5 * E[ZrC] - 0.5 * E[ZrN]</i> was about<i style="color: rgb(255, 0, 0);"> +0.4Ry/cell</i>, that means <i style="color: rgb(255, 0, 0);">+525.088kJ/mol</i>. According to thermodynamics principle, Zr(CN0.5) cannot be prepared. However, when I tried supercell-doping method, the formation energy was negative ( about <i><span style="color: rgb(255, 0, 0);">-10kJ/mol</span></i> ).<br>Therefore, I want to know why VCA method gives positive formation energy. Is it reasonable? If it should be like this, then can we believe the phonon result from VCA as reliable one?<br>Any suggestion will be welcomed. Thanks.<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>