<span style>To follow up on Nicola's point, the bulk modulus of sapphire is 240 GPa or 2400 kbar (</span><a href="http://www.mt-berlin.com/frames_cryst/descriptions/sapphire.htm" target="_blank" style>http://www.mt-berlin.com/frames_cryst/descriptions/sapphire.htm</a><span style>).  So a 1 kbar error corresponds to a very small change in volume.</span><div style>
<br></div><div style>If we use the quick and dirty Birch equation of state:</div><div style><br></div><div style>P(V) = 3/2 K0 [ (V0/V)^(7/3) - (V0/V)^(5/3)]</div><div style><br></div><div style>with K0 = 2400 kbar and ask what volume will produce a 1 kbar change in pressure we get</div>
<div style><br></div><div style>delta V/V0 = +/- 0.0004 </div><div style><br></div><div style>Considering the normal errors in DFT, it's not worth trying to converge the stress to the 1 kbar accuracy you're trying to achieve.<div class="yj6qo ajU" style="outline-style:none;outline-width:initial;outline-color:initial;padding-top:10px;padding-right:0px;padding-bottom:10px;padding-left:0px;width:22px;margin-top:2px;margin-right:0px;margin-bottom:0px;margin-left:0px">
<div id=":t1" class="ajR" tabindex="0" style="background-color:rgb(241,241,241);border-top-width:1px;border-right-width:1px;border-bottom-width:1px;border-left-width:1px;border-top-style:solid;border-right-style:solid;border-bottom-style:solid;border-left-style:solid;border-top-color:rgb(221,221,221);border-right-color:rgb(221,221,221);border-bottom-color:rgb(221,221,221);border-left-color:rgb(221,221,221);clear:both;line-height:6px;outline-style:none;outline-width:initial;outline-color:initial;width:20px">
<img class="ajT" src="https://mail.google.com/mail/u/0/images/cleardot.gif" style="background-image: url(https://ssl.gstatic.com/ui/v1/icons/mail/ellipsis.png); background-attachment: initial; background-origin: initial; background-clip: initial; background-color: initial; height: 8px; opacity: 0.3; width: 20px; background-position: initial initial; background-repeat: no-repeat no-repeat; "></div>
</div><div class="im adL" style="color:rgb(80,0,80)"><br><br><div class="gmail_quote">On Mon, May 7, 2012 at 10:07 AM, Nicola Marzari<span dir="ltr"><<a href="mailto:nicola.marzari@epfl.ch" target="_blank" style="color:rgb(17,85,204)">nicola.marzari@epfl.ch</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin-top:0px;margin-right:0px;margin-bottom:0px;margin-left:0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><br>How much does 1 kbar error translates into an error in lattice<br>
parameter? (keep atoms fixed, using relative coordinates, cutoff fixed,<br>and expand celldm(1) by 0.3% - what's the change in stress? that change<br>should be very well converged)<br><br></blockquote></div></div><span class="HOEnZb adL"><font color="#888888">-- <div>
Michael Mehl</div><div>US Naval Research Laboratory</div><div>Washington DC</div><div>(Home email)<br></div></font></span></div><br><div class="gmail_quote"><br></div>