Dear all,<div><br></div><div>We are trying to implement a workflow that performs cell optimisation. An important part of the workflows is the Pulay stress [1] analysis.- estimation of the error in the stress tensor introduced by the incompleteness of the basis set. </div><div><br></div><div>Essentially, to compute the Pulay stress one needs to re-run the SCF (re-initialising the plane-wave basis set) with the final geometry.</div><div><br></div><div>At the same time, I believe, re-initialisation of the basis set happens before the very final SCF run (the one that starts with "Reevaluating energy at the minimum"). Can someone confirm that this is correct?</div><div><br></div><div>If that is correct, is there a possibility to print the stress tensor at this very final SCF step? We didn't find an option to enable that. Having this information printed would allow us to avoid re-running the SCF just to compute the final stress (the Pulay stress).</div><div><br></div><div>Thanks,</div><div>Best regards,</div><div>Sasha</div><div><br></div><div>[1] https://iopscience.iop.org/article/10.1088/0953-8984/2/19/007</div><div><br></div>