<div dir="auto">Dear Nicolas, <div dir="auto">What do you mean by direct numerical derivation? How many points did you use for the fit, and for the derivation? My guess is that the chisq (chi squared?) factor is purely accidental.</div><div dir="auto"><br></div><div dir="auto">Regards <br><br><div data-smartmail="gmail_signature" dir="auto">-- <br>Lorenzo Paulatto CNRS /SU <br>Written on a virtual keyboard with real fingers</div></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sun, 31 Mar 2019, 05:23 Nicolas Leconte, <<a href="mailto:lecontenicolas0@uos.ac.kr">lecontenicolas0@uos.ac.kr</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Dear QE users,<div><br></div><div>I have recently used the ev.x executable to obtain the equation of state from my energy versus volume data.</div><div><br></div><div>When I compare the Birch Murnaghan (B-M) 3d order equation of state results with a direct numerical derivation, I obtain the best matching between both methods. Yet, there is still some mismatch.</div><div><br></div><div>I tried reducing my data to only include smaller pressures and reduce the impact of higher order terms, but it doesn't really change much, while I thought it would.</div><div><br></div><div>Then, while playing with my output data, I noticed that simply multiplying my pressure curve from B-M with the value of chisq, I get perfect agreement between both methods.</div><div><br></div><div>Am I missing something as on how to use the output from this script. It feels like too good a coincidence.</div><div><br></div><div>Now, I wonder what precludes me from using this chisq as a prefactor to the bulk modulus B0 (<a href="https://en.wikipedia.org/wiki/Birch%E2%80%93Murnaghan_equation_of_state" target="_blank" rel="noreferrer">https://en.wikipedia.org/wiki/Birch%E2%80%93Murnaghan_equation_of_state</a>) to obtain a proper agreement with the numerically derived curve.</div><div><br></div><div>Regards,</div><div>Nicolas</div></div>
_______________________________________________<br>
users mailing list<br>
<a href="mailto:users@lists.quantum-espresso.org" target="_blank" rel="noreferrer">users@lists.quantum-espresso.org</a><br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a></blockquote></div>