<html><head><style type="text/css"><!-- DIV {margin:0px;} --></style></head><body><div style="font-family:times new roman,new york,times,serif;font-size:10pt">Dear Chen,<br><br>Can you please send me your FC files, so that I can play around?<br><br>Bests,<br>Eyvaz.<br>-------------------------------------------------------------------<br>Prof. Eyvaz Isaev, <br>Department of Physics, Chemistry, and Biology (IFM), Linkoping University, Sweden <br>Theoretical Physics Department, Moscow State Institute of Steel & Alloys, Russia, <br>isaev@ifm.liu.se, eyvaz_isaev@yahoo.com<div><br></div><div style="font-family: times new roman,new york,times,serif; font-size: 10pt;"><br><div style="font-family: times new roman,new york,times,serif; font-size: 12pt;"><font face="Tahoma" size="2"><hr size="1"><b><span style="font-weight: bold;">From:</span></b> Chen Li <lichen@caltech.edu><br><b><span style="font-weight: bold;">To:</span></b>
pw_forum@pwscf.org<br><b><span style="font-weight: bold;">Sent:</span></b> Fri, January 7, 2011 8:34:31 PM<br><b><span style="font-weight: bold;">Subject:</span></b> [Pw_forum] Conservation of projected phonon DOS<br></font><br>
Dear All,<div><br></div><div>I'm wondering if anyone else had encountered similar problem before. Any suggestions will be greatly appreciated.</div><div><br></div><div>I have been using the QHA code and our own code to calculated the atomic projected phonon DOS using the phonon eigenvectors from QE phonon calculation. Both methods gave very similar results, so I have no doubt about that part. </div>
<div><br></div><div>For two different structures, ScF3 (cubic) and ZrO2 (monoclinic), the projected phonon DOS didn't have the correct ratio between species. For example, one will expect that oxygen will have twice as much contribution as from Zr in ZrO2, but the actual ratio is about 1.4 instead. Each atom in unit cell is supposed to have 3 degrees of freedom but Zr ended up with 3.7 while oxygen is about 2.6. (It is strange that the total degrees of freedom still add up to 9, as expected. </div>
<div><br></div><div>I has used the our code on the eigenvectors from VASP calculation and it gives the correct ratio. (I even adapted the VASP eigenvector format to be used with QHA code. The results are good too.) It seems the eigenvectors from QE may have some problem or maybe there is some factors to them? I'm pretty confident about the convergence of the calculations. It is also unlikely the convergence of calculation will give that much difference because it should not change the phonon DOS conservation anyway.</div>
<div><br></div><div>Regards,<br clear="all">Chen Li<br>Materials Science<br>California Institute of Technology<br>MC 138-78 1200 E California Blvd<br>Pasadena, CA 91125<br>
</div>
</div></div>
</div><br>
</body></html>