
<html><head><meta http-equiv="Content-Type" content="text/html charset=us-ascii"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class=""><div class="">Dear QE developers,</div><div class=""><br class=""></div><div class="">this is a feature request, which I am sending to the mailing list in lieu of access to the QE bug tracker on <a href="http://qe-forge.org" class="">qe-forge.org</a></div><div class=""><br class=""></div><div class=""><br class=""></div><div class="">I am sure many of you are familiar with the AiiDA project (<a href="http://www.aiida.net/" class="">http://www.aiida.net/</a>), which provides not only a very convenient python infrastructure for running QE calculations, but also a database for a consistent storage of results.</div><div class=""> </div><div class="">We have recently come across a "conflict of philosophy" between QE (more specifically, pp.x) and AiiDA that causes considerable computational/data overhead.</div><div class=""><br class=""></div><div class="">Basically, we need to plot lots of KS orbitals for different nanostructures. </div><div class="">pp.x is designed to produce 1 cube file per execution and AiiDA is designed to produce 1 working directory per calculation. </div><div class="">I.e. when computing 100 KS orbitals from a band structure calculation, using AiiDA the information contained in the band structure directory is replicated 100 times. </div><div class=""> </div><div class="">Our solution so far has been to modify the source of pp.x in such a way that one run can produce N KS orbitals. This has the added benefit of saving a lot of "startup time" as compared to running pp.x N times (in a test run generating 21 KS orbitals, the speedup was ~3x but this does vary, of course).</div><div class="">We have already contacted AiiDA developers Nicolas Mounet and Giovanni Pizzi and they are in favor of this solution as well.</div><div class=""><br class=""></div><div class="">I attach a patch against QE 6.0 (git SHA 93ced686f67fe1f3250824358015e2ffde364960) that adds plot_num 22 and 23 in pp.x.</div><div class="">plot_num 23 allows to plot N KS orbitals and plot_num 22 allows to plot the local density of states on a grid of energies between emin and emax (useful for scanning tunneling spectroscopy simulations). The patch starts with a commit message that contains sample input files.</div><div class=""><br class=""></div><div class="">I think the main question is, whether you would consider deviating from the "1 run = 1 output file" policy. </div><div class="">If yes, then one can think of the details regarding the implementation (the patch, including the input format for the KS orbitals, can certainly be improved).</div><div class=""><br class=""></div><div class="">Best wishes,</div><div class="">Leopold Talirz, Carlo Pignedoli and Sasha Yakutovich</div><div class=""><br class=""></div><div class=""></div></body></html>