<div style="line-height:1.7;color:#000000;font-size:14px;font-family:arial">Dear Furman:<br>I think, the easier and better way to make such a "compress" is using isostatic pressure, which can be specified by parameter "press" in "cell" namelist, especially for the un-cubic system. For example, in HCP structure, the cell is determined by "a" and "a/c", while the volume changing, not only "a", but also "a/c" ratio changes, too.<br>Therefore, you can do a batch of "vc-relax" tasks with different press and gather the final volume and total energy to fit EOS. The aim compress, as 10% or 20% shrinkage, can be found from EOS fitting result, where you can obtain the P-V curve.<br><br><div>--<br>GAO Zhe<br>CMC Lab, Materials Science & Engineering Department,<br>Seoul National University, South Korea<br>        
 </div><div id="divNeteaseMailCard"></div><br>At 2012-06-08 15:17:36,"David Furman" <sirok4@gmail.com> wrote:<br> <blockquote id="isReplyContent" style="PADDING-LEFT: 1ex; MARGIN: 0px 0px 0px 0.8ex; BORDER-LEFT: #ccc 1px solid"><style></style><div class="WordSection1"><p class="MsoNormal" style="text-align:left;direction:ltr;unicode-bidi:embed">Hello dear users and developers,<o:p></o:p></p><p class="MsoNormal" style="text-align:left;direction:ltr;unicode-bidi:embed"><o:p> </o:p></p><p class="MsoNormal" style="text-align:left;direction:ltr;unicode-bidi:embed">I'm trying to calculate the equation of state of Europium oxide, but I have some doubts regarding my procedure.<o:p></o:p></p><p class="MsoNormal" style="text-align:left;direction:ltr;unicode-bidi:embed">I have the xyz coordinates of the unit cell in angstroms which I put into atom_positions and also define the needed celldm's().<o:p></o:p></p><p class="MsoNormal" style="text-align:left;direction:ltr;unicode-bidi:embed"><o:p> </o:p></p><p class="MsoNormal" style="text-align:left;direction:ltr;unicode-bidi:embed">I'm interested to compress the material to different values (say, 0%, 10% and 20% volumetric compression) and calculate the total energy each time.<o:p></o:p></p><p class="MsoNormal" style="text-align:left;direction:ltr;unicode-bidi:embed"><o:p> </o:p></p><p class="MsoNormal" style="text-align:left;direction:ltr;unicode-bidi:embed">Do I have to change the given initial xyz coordinates in atom_positions AND the celldm() , or just change the celldm() accordingly to the compression and pw.x will do the proper changes to the positions of the atoms ?<o:p></o:p></p><p class="MsoNormal" style="text-align:left;direction:ltr;unicode-bidi:embed"><o:p> </o:p></p><p class="MsoNormal" style="text-align:left;direction:ltr;unicode-bidi:embed">I have done some calculations with uncompressed unit cell and with 20% compressed unit cell (by changing to the proper compressed celldm() values) but the difference in energy was about 0.5 Ry which i believe to be quite small. (Of course I checked previously for convergence).<o:p></o:p></p><p class="MsoNormal" style="text-align:left;direction:ltr;unicode-bidi:embed"><o:p> </o:p></p><p class="MsoNormal" style="text-align:left;direction:ltr;unicode-bidi:embed">Am I doing it wrong ?<o:p></o:p></p><p class="MsoNormal" style="text-align:left;direction:ltr;unicode-bidi:embed"><o:p> </o:p></p><p class="MsoNormal" style="text-align:left;direction:ltr;unicode-bidi:embed">Any help will be greatly appreciated <o:p></o:p></p><p class="MsoNormal" style="text-align:left;direction:ltr;unicode-bidi:embed"><o:p> </o:p></p><p class="MsoNormal" style="text-align:left;direction:ltr;unicode-bidi:embed"><b><span style="font-family:"Candara","sans-serif";color:#00B0F0">David Furman, Graduate Student</span></b><span style="font-family:"Candara","sans-serif";color:#00B0F0"> </span><span style="font-size:10.0pt;font-family:"Candara","sans-serif";color:#00B0F0">| Fritz Haber Research Center for <o:p></o:p></span></p><p class="MsoNormal" style="text-align:left;direction:ltr;unicode-bidi:embed"><span style="font-size:10.0pt;font-family:"Candara","sans-serif";color:#00B0F0">Molecular Dynamics</span><span style="font-family:"Candara","sans-serif";color:#00B0F0">, </span><span style="font-size:10.0pt;font-family:"Candara","sans-serif";color:#00B0F0">The Hebrew University of Jerusalem | </span><b><span style="font-family:"Candara","sans-serif";color:#00B0F0">+972-8656-8909</span></b><b><o:p></o:p></b></p><p class="MsoNormal" dir="RTL"><span dir="LTR"><o:p> </o:p></span></p></div></blockquote></div><br><br><span title="neteasefooter"><span id="netease_mail_footer"></span></span>