<div dir="ltr">Dear cp2k community,<br><br>I am calculating the analytical stress tensor of a cubic diamond 3x3x3 supercell using the popt/psmp version of cp2k.<br>I observe that the diagonal elements of the analytical stress tensor depend on the number of mpi processes. Only when running cp2k with one single mpi task, the analytical stress tensor roughly agrees with the numerical one.<br>I encountered this issue with cp2k 8.1 + openmpi 4.0 (compiled with gcc 9.2), and also with intel mpi 2018.3 (intel compiler suite 18.0).<br><br>Any suggestions on how to fix/avoid this inconsistency would be appreciated. Thanks!<br><br>In the attachment, you can find the input file, an example of an output file, and the computed diagonal elements of the stress tensor for different number of mpi processes. The problem does not affect the off-diagonal entries.<br><br>Kind Regards,<br>Severin</div>