
<html><head><meta http-equiv="Content-Type" content="text/html; charset=us-ascii"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">Dear Jakob,<div class=""><br class=""></div><div class="">Do you mean the Helmoltz free energy (F)? </div><div class=""><br class=""></div><div class="">As by definition, the Gibbs free energy (G) does not depend upon volume, but on pressure and temperature.</div><div class="">At p=0, then they are the same, but notice always that:</div><div class=""><br class=""></div><div class=""><ul class="MailOutline"><li class="">G(p,T) = min_{V}[F(V, T)+pV] </li></ul></div><div class=""><br class=""></div><div class="">Thus, at p=0, you still need to take the minimum of F at the fixed temperature.</div><div class=""><br class=""></div><div class="">Moreover, I guess you are using the quasi-harmonic approximation. There exist several flavours. Unless your system is isotropic, scaling the volume isotropically can be a source of error (due to a physical assumption, not due to the numerics).</div><div class=""><br class=""></div><div class="">Best,</div><div class="">Lorenzo Bastonero</div><div class=""><br class=""><div class="">

<meta charset="UTF-8" class=""><div dir="auto" style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0); letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class=""><div dir="auto" style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0); letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class=""><div dir="auto" style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0); letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class=""><div dir="auto" style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0); letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class=""><div dir="auto" style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0); letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class=""><div><div class="">********************************</div><div class="">Lorenzo Bastonero</div><div class=""><br class=""></div><div class="">PhD Student</div><div class="">U Bremen Excellence Chair,<br class="">Bremen Center for Computational Materials Science,<br class="">and MAPEX Center for Materials and Processes</div><div class=""><br class=""></div><div class="">University of Bremen</div><div class="">Faculty of Production Engineering</div><div class="">TAB-Building, Room 3.32</div><div class="">Am Fallturm 1</div><div class="">28359 Bremen, DE</div><br class="Apple-interchange-newline"><a href="http://www.hmi.uni-bremen.de/" class="">http://www.hmi.uni-bremen.de/</a></div></div></div></div></div></div>

</div>

<div><br class=""><blockquote type="cite" class=""><div class="">Il giorno 4 nov 2022, alle ore 09:20, Jakob Kraus <<a href="mailto:jakob.kraus@physik.tu-freiberg.de" class="">jakob.kraus@physik.tu-freiberg.de</a>> ha scritto:</div><br class="Apple-interchange-newline"><div class=""><div class="">Dear QE users,<br class=""><br class=""><br class="">I am trying to calculate heat capacities, and to this end, I have calculated the Gibbs energy as a function of volume<br class=""><br class="">at several temperatures in a range of 0-700 K.<br class=""><br class="">However, the functions that result from this do not show a single minimum, but rather a whole lot of oscillations, which<br class=""><br class="">seem to originate from the phonons, since the E(V) diagram at 0 K shows a single minimum as expected.<br class=""><br class="">I am thankful for any help you can offer with getting rid of these oscillations.<br class=""><br class="">I have attached input files and some plots to showcase the problem.<br class=""><br class=""><br class="">Best regards,<br class=""><br class=""><br class="">Jakob Kraus, TU Freiberg<br class=""><span id="cid:637D156A-0221-4117-96E7-92B52C7875FD"><G_Zr_BH4_0.0.png></span><span id="cid:9FAF9B70-5E90-4863-8238-B99F8B87FB8C"><E_Zr_BH4.png></span><span id="cid:E953F3A7-A318-4BE4-ACA2-BEB7DE8D57C9"><espresso.qi></span><span id="cid:136CDA8E-FB90-4F21-A113-D7AA9AD35CD0"><espresso.pwi></span><span id="cid:E8103B0F-6182-4463-9C17-B48A0B311D79"><espresso.phi></span><span id="cid:61409018-8555-4F01-AAA8-1C71B5C2FABA"><espresso.mdi></span>_______________________________________________<br class="">The Quantum ESPRESSO community stands by the Ukrainian<br class="">people and expresses its concerns about the devastating<br class="">effects that the Russian military offensive has on their<br class="">country and on the free and peaceful scientific, cultural,<br class="">and economic cooperation amongst peoples<br class="">_______________________________________________<br class="">Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu" class="">www.max-centre.eu</a>)<br class="">users mailing list <a href="mailto:users@lists.quantum-espresso.org" class="">users@lists.quantum-espresso.org</a><br class=""><a href="https://lists.quantum-espresso.org/mailman/listinfo/users" class="">https://lists.quantum-espresso.org/mailman/listinfo/users</a></div></div></blockquote></div><br class=""></div></body></html>