<div dir="ltr">Finding the bulk modulus wth ev.x is essentially calculating<div><br></div><div>B(V) = V E''(V).</div><div><br></div><div>What the heck is E(V)?  That's actually rather simple:  E(V) is the minimum energy of the the system at fixed volume.  For your system,</div><div><br></div><div>E(V) = min eng(V,c/a)     for all values of c/a at fixed volume V,</div><div><br></div><div>where eng(V,c/a) is the energy you get from the relax mode for given values of V and c/a.</div><div><br></div><div>For your tetragonal system, you want E(V) for several different volumes.  Here, V = a*a*c, or a*a*c/2 if you've got a body-centered system.</div><div><br></div><div>So:</div><div><br></div><div>1) Pick a volume V and initial value c/a. </div><div><br></div><div>2) Use pw.x's relax mode to find the equilibrium positions of the atoms for this cell.  Save the energy.</div><div><br></div><div>3) Adjust c/a -- keeping the volume fixed -- until you find the minimum.  That is, the lowest energy eng(V,c/a) for all values of c/a at this particular volume.  Plot eng(V,c/a) versus c/a at fixed V if that helps you see what's going on.  There will be at least one minimum, somewhere, for any reasonable choice of volume.</div><div><br></div><div>4) Repeat for enough volumes to get a good number from ev.x.  Ideally your volumes should bracket the equilibrium volume.  Again, plotting helps.</div><div><br></div><div>5) Use the pair (V,E(V)) as input to ev.x.  Out will come the equilibrium volume and bulk modulus.</div><div><br></div><div>Note that your procedure only gives you the minimum energy structure at one volume -- equilibrium.</div><div><br></div><div>This is ''easily'' generalized to other types of lattices, though you might want to look at the cell_dofree='shape" option for that.</div><div><br></div><div>Alternatively, you can compute the bulk modulus from the elastic constants, but for a tetragonal system this is not trivial -- read up on it in the literature.  I have the formula in several of my papers from the early 1990s, but unfortunately I don't have the references on this computer.</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Sat, Jun 10, 2017 at 2:50 PM, Nadire Nayir <span dir="ltr"><<a href="mailto:nadire.nayr42@gmail.com" target="_blank">nadire.nayr42@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div class="gmail_default" style="font-family:'times new roman',serif;font-size:small;color:rgb(76,17,48)">Dear all,</div><div class="gmail_default" style="font-family:'times new roman',serif;font-size:small;color:rgb(76,17,48)">I deal with the cell optimization and the equation of states of the trigonal structures (especially, GeO2). To get the optimized structure, I started from the experimental structure and performed a series of the relax calculations with "relax" keyword to it. During the relax calculations, first, I changed celldm(1)=a and celldm(3)=c/a parameters by keeping c parameter constant and I found the optimized celldm(1), after that, I kept celldm(1) fixed at its optimized value and I started to change only celldm(3) parameter until getting the optimized c parameter. </div><div class="gmail_default" style="font-family:'times new roman',serif;font-size:small;color:rgb(76,17,48)"><br></div><div class="gmail_default" style="font-family:'times new roman',serif;font-size:small;color:rgb(76,17,48)">After getting the optimized celldm(1) and celldm(3), I calculated the equation of states of the structure by applying SCF calculation and finally, I used ev.x to calculate the bulk modulus with the non-cubic option, which is ~140 GPa.</div><div class="gmail_default" style="font-family:'times new roman',serif;font-size:small;color:rgb(76,17,48)"><br></div><div class="gmail_default" style="font-family:'times new roman',serif;font-size:small;color:rgb(76,17,48)">However, the bulk modulus value I found is so different from the experimental value which is 42 GPa.  </div><div class="gmail_default" style="font-family:'times new roman',serif;font-size:small;color:rgb(76,17,48)"><br></div><div class="gmail_default" style="font-family:'times new roman',serif;font-size:small;color:rgb(76,17,48)">When I applied the same optimization process to the cubic structures (Si-diamond and Ge-diamond), I achieved the almost same results with their experimental values. </div><div class="gmail_default" style="font-family:'times new roman',serif;font-size:small;color:rgb(76,17,48)"><br></div><div class="gmail_default" style="font-family:'times new roman',serif;font-size:small;color:rgb(76,17,48)">I did not understand why this procedure does not work for the trigonal structures. I wonder, is there anything that I missed? </div><div class="gmail_default" style="font-family:'times new roman',serif;font-size:small;color:rgb(76,17,48)"><br></div><div class="gmail_default" style="font-family:'times new roman',serif;font-size:small;color:rgb(76,17,48)">By the way, I also applied vc-relax calculation to the trigonal structure, again, I got the larger bulk modulus than the experimental one!</div><div class="gmail_default" style="font-family:'times new roman',serif;font-size:small;color:rgb(76,17,48)"><br></div><div class="gmail_default" style="font-family:'times new roman',serif;font-size:small;color:rgb(76,17,48)">I would greatly appreciate if you could kindly help me with this problem.</div><div class="gmail_default" style="font-family:'times new roman',serif;font-size:small;color:rgb(76,17,48)"><br></div><div class="gmail_default" style="font-family:'times new roman',serif;font-size:small;color:rgb(76,17,48)">Best regards,</div><div class="gmail_default" style="font-family:'times new roman',serif;font-size:small;color:rgb(76,17,48)">Nadire</div><span class="HOEnZb"><font color="#888888"><div class="gmail_default" style="font-family:'times new roman',serif;font-size:small;color:rgb(76,17,48)"><pre style="font-family:courier,'courier new',monospace;font-size:14px;white-space:pre-wrap;word-wrap:break-word;margin-top:0px;margin-bottom:0px;color:rgb(0,0,0);line-height:19.6px"><br></pre></div>-- <br><div class="m_3198701642783865246gmail_signature"><div dir="ltr">Nadire Nayir<div>Research Assistant</div><div><div>Middle East Technical University</div><div>Physics Department</div></div></div></div>
</font></span></div>
<br>______________________________<wbr>_________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
<a href="http://pwscf.org/mailman/listinfo/pw_forum" rel="noreferrer" target="_blank">http://pwscf.org/mailman/<wbr>listinfo/pw_forum</a><br></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><span><div dir="ltr" style="margin-left:0pt"><span><div dir="ltr" style="margin-left:0pt"><table style="border:none;border-collapse:collapse"><colgroup><col width="73"><col width="240"></colgroup><tbody><tr style="height:0px"><td style="border-left:solid #ffffff 1px;border-right:solid #ffffff 1px;border-bottom:solid #ffffff 1px;border-top:solid #ffffff 1px;vertical-align:top;padding:7px 7px 7px 7px"><br><img src="https://lh4.googleusercontent.com/cWahY61df5ZQz3h0L--aoIVBK4j_u5sWfktSfeRgcbi3bi_PKu-Vtt8OPXlJEEX7_J8UVfbP6EYLXroCVHzqegMVgjl4KD3yBpnJZp0zwht66TTfYleiGnAbqPTj7TEwVxNRrdA6" width="57" height="80" style="border:none" alt="USNA_Gold_Seal.png"></td><td style="border-left:solid #ffffff 1px;border-right:solid #ffffff 1px;border-bottom:solid #ffffff 1px;border-top:solid #ffffff 1px;vertical-align:top;padding:7px 7px 7px 7px"><p dir="ltr" style="line-height:1.2;margin-top:0pt;margin-bottom:0pt"><span style="font-size:14.6667px;font-family:Arial;color:rgb(0,0,0);vertical-align:baseline;white-space:pre-wrap;background-color:transparent">Michael J. Mehl, Ph.D.     </span></p><p dir="ltr" style="line-height:1.2;margin-top:0pt;margin-bottom:0pt"><span style="font-size:14.6667px;font-family:Arial;color:rgb(0,0,0);vertical-align:baseline;white-space:pre-wrap;background-color:transparent">Kinnear Chair in Physics</span></p><p dir="ltr" style="line-height:1.2;margin-top:0pt;margin-bottom:0pt"><span style="font-size:14.6667px;font-family:Arial;color:rgb(0,0,0);vertical-align:baseline;white-space:pre-wrap;background-color:transparent">The United States Naval Academy</span></p><p dir="ltr" style="line-height:1.2;margin-top:0pt;margin-bottom:0pt"><span style="font-size:14.6667px;font-family:Arial;color:rgb(0,0,0);vertical-align:baseline;white-space:pre-wrap;background-color:transparent">Mail Stop 9C</span></p><p dir="ltr" style="line-height:1.2;margin-top:0pt;margin-bottom:0pt"><span style="font-size:14.6667px;font-family:Arial;color:rgb(0,0,0);vertical-align:baseline;white-space:pre-wrap;background-color:transparent">572 Holloway Road</span></p><p dir="ltr" style="line-height:1.2;margin-top:0pt;margin-bottom:0pt"><span style="font-size:14.6667px;font-family:Arial;color:rgb(0,0,0);vertical-align:baseline;white-space:pre-wrap;background-color:transparent">Chauvenet Hall 257</span></p><p dir="ltr" style="line-height:1.2;margin-top:0pt;margin-bottom:0pt"><span style="font-size:14.6667px;font-family:Arial;color:rgb(0,0,0);vertical-align:baseline;white-space:pre-wrap;background-color:transparent">Annapolis MD 21402</span></p><p dir="ltr" style="line-height:1.2;margin-top:0pt;margin-bottom:0pt"><span style="font-size:14.6667px;font-family:Arial;color:rgb(0,0,0);vertical-align:baseline;white-space:pre-wrap;background-color:transparent">(410)293-6685</span></p><p dir="ltr" style="line-height:1.2;margin-top:0pt;margin-bottom:0pt"><span style="font-size:14.6667px;font-family:Arial;color:rgb(0,0,0);vertical-align:baseline;white-space:pre-wrap;background-color:transparent"><a href="mailto:mmehl@usna.edu" target="_blank">mmehl@usna.edu</a></span></p><p dir="ltr" style="line-height:1.2;margin-top:0pt;margin-bottom:0pt"><a href="http://aflow.org/CrystalDatabase/" target="_blank">Library of Crystallographic Prototypes</a></p></td></tr></tbody></table></div></span></div></span></div><span><br></span></div></div></div></div></div></div></div>
</div>