Dear developers and users,<div><br></div><div>We are looking at the equilibrium density for an AlCl3 molten salt melt by measuring pressure values over 10 ps in NVT at a variety of volumes. We are comparing the pressures for a variety of functionals with various vdW corrections: PBE-D3, optB86b-vdW, optB88-vdW, SCAN-D3, and SCAN-RVV10. The data are shown in the attached plot with the mean pressure (± 1 standard deviation) for each volume (the plot for PBE-D3-smooth can be ignored).</div><div><br></div><div>It seems optB86b-vdW, optB88-vdW, and SCAN-RVV10 are reporting unexpected behavior in the stress tensors. Both optB86b and optB88 seem to be the most "wrong" with negative pressures across the board and an unusual <i>increase</i> in the pressure as the system's volume is expanded (all systems have the same number of particles). SCAN-RVV10 at least follows a believable trend and isn't too far off, so it may just not be appropriate for this system.</div><div><br></div><div>We are wondering if the stress tensor calculation is properly implemented for these non-local vdW functionals, since both PBE-D3 and SCAN-D3 are producing believable trends.</div><div><br></div><div>For reference, I have included the input files for the 3 problematic functionals I mentioned for the volumes corresponding to V_exp in the attached plots. The versions used for each of those runs are listed below:</div><div><ul><li>optB86b-vdW: CP2K-8.2 </li><li>optB88-vdW: CP2K-6.1</li><li>SCAN-RVV10: CP2K-2022.2</li></ul>Any help or insight is appreciated!<br></div><div><br></div><div>Thank you,</div><div><br></div><div>Luke Gibson</div>

<p></p>

-- <br />
You received this message because you are subscribed to the Google Groups "cp2k" group.<br />
To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp2k+unsubscribe@googlegroups.com">cp2k+unsubscribe@googlegroups.com</a>.<br />
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/afc955e0-daad-40f4-8bce-87b08f87cb8dn%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/afc955e0-daad-40f4-8bce-87b08f87cb8dn%40googlegroups.com</a>.<br />