<div dir="ltr">Dear Paolo,<br><br>Thank you very much for your inputs. Before I try some modifications with the code, I wanted to ask:Does pp.x even print out the numerical contributions to the total energy? Whenever I have used it, I've only generated potentials in the form of FFT grids. This could tell me the variation of potential along different surfaces but this is not exactly what I am aiming for. I am looking to deconstruct the total energy from my DFT calculation into its different components, including exchange and correlation.<br><br>Thanks again for your help.<br><br>Sincerely,<br>Shivesh Sivakumar<br>University of Washington-Seattle<br>WA-98105  <br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, Jan 19, 2021 at 1:52 AM Paolo Giannozzi <<a href="mailto:p.giannozzi@gmail.com">p.giannozzi@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div dir="ltr">On Tue, Jan 19, 2021 at 10:15 AM Shivesh Sivakumar <<a href="mailto:shiveshsivakumar@gmail.com" target="_blank">shiveshsivakumar@gmail.com</a>> wrote:<br></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><br>I tried something like 'gfortran v_of_rho.f90 <input> > <out>' but it didn't work out </div></blockquote><div><br></div><div>it's not that simple: that subroutine requires a bunch of stuff to be initialized and must be called with the needed input and output variables. You may have better luck by modifying one of the many post processing codes in PP/src that read final SCF data: they call 'read_file_new' that calls 'v_of_rho' to rebuild the potential from the charge. IIRC one of those codes (ppacf.f90) actually performs some splitting of exchange and correlation terms, in addition to many other operations.<br></div><div><br></div><div>Paolo<br></div><div><br></div></div>-- <br><div dir="ltr"><div dir="ltr"><div><div dir="ltr"><div>Paolo Giannozzi, Dip. Scienze Matematiche Informatiche e Fisiche,<br>Univ. Udine, via delle Scienze 206, 33100 Udine, Italy<br>Phone +39-0432-558216, fax +39-0432-558222<br><br></div></div></div></div></div></div>
_______________________________________________<br>
Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu" rel="noreferrer" target="_blank">www.max-centre.eu</a>)<br>
users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a><br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a></blockquote></div>