<div dir="ltr">On Tue, Jul 28, 2015 at 5:26 PM, anyy zsj <span dir="ltr"><<a href="mailto:zsjanyy@gmail.com" target="_blank">zsjanyy@gmail.com</a>></span> wrote:<div class="gmail_extra"><div class="gmail_quote"><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">   1.4989396432   -6.0466776328<div><div>   1.5430233198   -6.0795705718</div><div>   <a href="tel:1.5879629631" value="+15879629631" target="_blank">1.5879629631</a>   -6.0993455319</div><div>   1.6804430727   -5.<a href="tel:9859764266" value="+19859764266" target="_blank">9859764266</a></div></div><div>The first column is volume and the second energy with Ry. After ev.x I input Ang, noncubic, 4, it gives me a results of 0.40 for minimal volume. It is certainly wrong.</div></div></blockquote><div><br></div><div>of course it is: ev.x assumes that the input data corresponds to a physical EOS. Your data isn't physical: the volume is way too small<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><br></div><div>Another example for the below data:</div><div><div>10.1      -48.570860572794118</div><div>10.3      -49.368013045588235</div><div>10.5      -49.714550792647060</div><div>10.7      -49.667482497058820</div><div>10.9      -49.305953287500003</div></div><div>The first column is lattice constant with angstrom and the second energy with Ry.</div></div></blockquote><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div> After ev.x I input Ang, sc, 4, it gives me 10.50369, while I get a different value 10.57517 by fitting to the formula in PRB 28,5480 (1983) of Murnaghan's equation.: <a href="http://gilgamesh.cheme.cmu.edu/doc/software/jacapo/appendices/appendix-eos.html" target="_blank">http://gilgamesh.cheme.cmu.edu/doc/software/jacapo/appendices/appendix-eos.html</a></div></div></blockquote><div><br> <div>sure that the first column is in angstrom and not a.u.? This is what I get assuming it is in a.u.:<br><br> # equation of state: murnaghan.        chisq = 0.5402D-05<br># a0 = 10.5740 a.u., k0 =13776 kbar, dk0 =  4.17 d2k0 =  0.000 emin =  -49.73649<br># a0 =  5.59554 Ang, k0 =1377.7 GPa,  V0 = ******* (a.u.)^3,  V0 =175.197 A^3 <br></div><br></div><div>note the small value of the |chi|^2 (chisq) of the fit<br><br></div><div>Paolo<br><br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div>I have tested many groups of data but the discrepancy is there for those large lattice parameters. Does ev.x use another equation or something else I have missing? Thank you very much.<br></div></div>
<br>_______________________________________________<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/listinfo/pw_forum</a><br></blockquote></div><br clear="all"><br>-- <br><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><span><span><font color="#888888">Paolo Giannozzi, Dept. Chemistry&Physics&Environment,<br>
Univ. Udine, via delle Scienze 208, 33100 Udine, Italy<br>
Phone <a href="tel:%2B39-0432-558216" value="+390432558216" target="_blank">+39-0432-558216</a>, fax <a href="tel:%2B39-0432-558222" value="+390432558222" target="_blank">+39-0432-558222</a></font></span></span></div></div></div></div>
</div></div>