<div dir="ltr">Dear QE developers,<div>     I have been using Quantum Espresso for more than a decade and recently a paper (PRB 91, 024107 (2015)) caught my attention, which shows that the total energy, pressure and cell relaxation are ill-defined in charged systems (even if a compensating charge background is present). The paper shows that the total energy and pressure of a charged system depend on the average electrostatic potential convention, and thus is code-dependent (QE, ABINIT and VASP use different conventions).</div><div>    I think the paper is theoretically correct. However in practice, we have recently tested doped ferroelectric materials using both QE and VASP. In this charged system, we find both QE and VASP yield qualitatively consistent structural properties.</div><div>    Since relaxation of a charged system is very common in literature using QE or VASP (e.g. PRL 109, 247601 (2012) and many others), I am wondering whether the relaxed structure of a charged system from QE is reliable. According to PRB 91, 024107, it is not legitimate to do cell relaxation of a charged system. However both QE and VASP do allow to calculate total energy and perform cell relaxation in a charged system. And from literature, it seems that the relaxed structure of a charge system from QE or VASP calculations is "not too wrong".</div><div>    So what is your insight on this issue? Shall we prohibit ourselves from doing cell relaxation of a charged system in QE, or can we still use the standard procedure but with the caution that the optimized structure of a charged system from QE might not be accurate (due to the QE convention of average electrostatic potential)?</div><div>     Thank you very much.</div><div><br></div><div>Dr. Hanghui Chen</div><div>Department of Physics </div><div>NYU Shanghai and New York University</div><div>     </div></div>