<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40"> <head> <meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1"> <meta name="Generator" content="Microsoft Word 15 (filtered medium)"> <style></style> </head> <body lang="en-CH" link="blue" vlink="purple" style="word-wrap:break-word"> <div class="WordSection1"> <p class="MsoNormal"><span lang="DE-CH" style="font-size:11.0pt;mso-fareast-language:EN-US">Hi<o:p></o:p></span></p> <p class="MsoNormal"><span lang="DE-CH" style="font-size:11.0pt;mso-fareast-language:EN-US"><o:p> </o:p></span></p> <p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;mso-fareast-language:EN-US">The GMAX parameter for Ewald is a bit sloppy in these inputs. Suggested values are 1 to 2 points per Angstrom based on the cell edge lengths, i.e. GMAX 20 is numerically more appropriate. In the example I gave you the GMAX value is automatically adapted to the cell size.<o:p></o:p></span></p> <p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;mso-fareast-language:EN-US"><o:p> </o:p></span></p> <p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;mso-fareast-language:EN-US">Best<o:p></o:p></span></p> <p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;mso-fareast-language:EN-US"><o:p> </o:p></span></p> <p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;mso-fareast-language:EN-US">Matthias<o:p></o:p></span></p> <p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;mso-fareast-language:EN-US"><o:p> </o:p></span></p> <div id="mail-editor-reference-message-container"> <div> <div style="border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0cm 0cm 0cm"> <p class="MsoNormal" style="mso-margin-top-alt:0cm;margin-right:0cm;margin-bottom:12.0pt;margin-left:36.0pt"> <b><span style="font-size:12.0pt;font-family:"Aptos",sans-serif;color:black">From: </span></b><span style="font-size:12.0pt;font-family:"Aptos",sans-serif;color:black">cp2k@googlegroups.com <cp2k@googlegroups.com> on behalf of Léon Luntadila Lufungula <Leon.luntadilalufungula@uantwerpen.be><br> <b>Date: </b>Monday, 22 January 2024 at 14:13<br> <b>To: </b>cp2k <cp2k@googlegroups.com><br> <b>Subject: </b>Re: [CP2K:19788] Cons. Qty for different ensembles + Pressure fluctuation<o:p></o:p></span></p> </div> <p class="MsoNormal" style="margin-left:36.0pt"><span style="font-size:11.0pt">P.S. These calculations were done based of the input files in https://www.cp2k.org/exercises:2018_uzh_acpc2:mol_sol, because the Fist regtest-water folder was not yet available at the time, but the FF params are correct, the coordinates are just different than the one molecule box in the tests directory.<o:p></o:p></span></p> <div> <p class="MsoNormal" style="margin-left:36.0pt"><span style="font-size:11.0pt"><o:p> </o:p></span></p> </div> <div> <div> <p class="MsoNormal" style="margin-left:36.0pt"><span style="font-size:11.0pt">On Monday 22 January 2024 at 14:03:06 UTC+1 Léon Luntadila Lufungula wrote:<o:p></o:p></span></p> </div> <blockquote style="border:none;border-left:solid #CCCCCC 1.0pt;padding:0cm 0cm 0cm 6.0pt;margin-left:4.8pt;margin-right:0cm"> <p class="MsoNormal" style="margin-left:36.0pt"><span style="font-size:11.0pt">Dear Matthias,<o:p></o:p></span></p> <div> <p class="MsoNormal" style="margin-left:36.0pt"><span style="font-size:11.0pt"><o:p> </o:p></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36.0pt"><span style="font-size:11.0pt">I have similar results for a box of 256 water molecules:<o:p></o:p></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36.0pt"><span style="font-size:11.0pt">NVT:<o:p></o:p></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36.0pt"><span style="font-size:11.0pt"><img id="_x0000_i1031" src="https://groups.google.com/group/cp2k/attach/721f9288cda1e/NVT-256.png?part=0.1&view=1"><br> NPT: (fluctuations are a bit smaller but still in the range of several kbar...)<o:p></o:p></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36.0pt"><span style="font-size:11.0pt"><img id="_x0000_i1030" src="https://groups.google.com/group/cp2k/attach/721f9288cda1e/NPT-256.png?part=0.3&view=1"><img id="_x0000_i1029" src="https://groups.google.com/group/cp2k/attach/721f9288cda1e/Pressure-256.png?part=0.2&view=1"><o:p></o:p></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36.0pt"><span style="font-size:11.0pt"><o:p> </o:p></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36.0pt"><span style="font-size:11.0pt">You say that the sample should be equilibrated well in advance before starting the sampling run, so do you mean I should run an NVT simulation for longer until it equilibrates before running an NPT from the thermalized end structure? Or do you mean both NVT and NPT should be run for longer until the Cons. Qty. is also equilibrated? Because in both cases the temperature seems to be equilibrated quite well and the pressure in the NPT ensemble is also constant albeit with very large oscillations...<o:p></o:p></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36.0pt"><span style="font-size:11.0pt"><o:p> </o:p></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36.0pt"><span style="font-size:11.0pt">Kind regards,<o:p></o:p></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36.0pt"><span style="font-size:11.0pt">Léon<o:p></o:p></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36.0pt"><span style="font-size:11.0pt"><o:p> </o:p></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36.0pt"><span style="font-size:11.0pt"><o:p> </o:p></span></p> </div> <div> <div> <p class="MsoNormal" style="margin-left:36.0pt"><span style="font-size:11.0pt">On Monday 22 January 2024 at 13:40:05 UTC+1 Krack Matthias wrote:<o:p></o:p></span></p> </div> <blockquote style="border:none;border-left:solid #CCCCCC 1.0pt;padding:0cm 0cm 0cm 6.0pt;margin-left:4.8pt;margin-right:0cm"> <div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:36.0pt"> <span lang="DE-CH" style="font-size:11.0pt">Dear Léon</span><span style="font-size:11.0pt"><o:p></o:p></span></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:36.0pt"> <span lang="DE-CH" style="font-size:11.0pt"> </span><span style="font-size:11.0pt"><o:p></o:p></span></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:36.0pt"> <span lang="EN-US" style="font-size:11.0pt">a 32 H2O cell is still a relatively small system which can cause large fluctuations. I suggest to try at least the 6x6x6 = 216 H2O <a href="https://github.com/cp2k/cp2k/blob/master/tests/Fist/regtest-water/H2O-TIP3P-6x6x6.inp" target="_blank"> input</a>. The sample should be equilibrated well in advance before starting the sampling run.</span><span style="font-size:11.0pt"><o:p></o:p></span></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:36.0pt"> <span lang="EN-US" style="font-size:11.0pt"> </span><span style="font-size:11.0pt"><o:p></o:p></span></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:36.0pt"> <span lang="EN-US" style="font-size:11.0pt">HTH</span><span style="font-size:11.0pt"><o:p></o:p></span></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:36.0pt"> <span lang="EN-US" style="font-size:11.0pt"> </span><span style="font-size:11.0pt"><o:p></o:p></span></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:36.0pt"> <span lang="EN-US" style="font-size:11.0pt">Matthias</span><span style="font-size:11.0pt"><o:p></o:p></span></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:36.0pt"> <span lang="EN-US" style="font-size:11.0pt"> </span><span style="font-size:11.0pt"><o:p></o:p></span></p> </div> </div> <div> <div> <div> <div> <div style="border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0cm 0cm 0cm"> <p class="MsoNormal" style="mso-margin-top-alt:auto;margin-bottom:12.0pt;margin-left:72.0pt"> <b><span style="font-size:12.0pt;font-family:"Aptos",sans-serif;color:black">From: </span></b><span style="font-size:12.0pt;font-family:"Aptos",sans-serif;color:black">cp...@googlegroups.com <cp...@googlegroups.com> on behalf of Léon Luntadila Lufungula <Leon.luntad...@uantwerpen.be><br> <b>Date: </b>Monday, 22 January 2024 at 12:56<br> <b>To: </b>cp2k <cp...@googlegroups.com><br> <b>Subject: </b>[CP2K:19784] Cons. Qty for different ensembles + Pressure fluctuation</span><span style="font-size:11.0pt"><o:p></o:p></span></p> </div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:72.0pt"> <span style="font-size:11.0pt">Dear all,<o:p></o:p></span></p> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:72.0pt"> <span style="font-size:11.0pt"> <o:p></o:p></span></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:72.0pt"> <span style="font-size:11.0pt">I ran several simulations on a box of 32 water molecules with the TIP3P force field:<o:p></o:p></span></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:72.0pt"> <span style="font-size:11.0pt">1) 1 ns NVE trajectory with timestep 0.5 fs<o:p></o:p></span></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:72.0pt"> <span style="font-size:11.0pt">2) 1 ns NVT trajectory at 298.15K with timestep 0.5 fs and thermostat time constant of 100fs<o:p></o:p></span></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:72.0pt"> <span style="font-size:11.0pt"> 3) 1 ns NPT_I trajectory restarted from the final structure of the NVT trajectory at 298.15K and 1 bar with timestep 0.5 and thermostat and barostat time constants of 100 fs<o:p></o:p></span></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:72.0pt"> <span style="font-size:11.0pt"> <o:p></o:p></span></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:72.0pt"> <span style="font-size:11.0pt">I plotted the Cons. Qty [a.u.] vs. Time [fs] for the three simulations and saw that the constant of motion was only more or less constant for the NVE trajectory. The NVT and NPT trajectories both had strong oscillations and significant drift of the constant of motion (see below) and I was wondering if this is normal behaviour or that this means that the system is not yet equilibrated? I would think that the Cons. Qty. should remain constant but I can also understand that the thermostat and barostat might influence its behaviour in a way.<o:p></o:p></span></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:72.0pt"> <span style="font-size:11.0pt">NVE:<o:p></o:p></span></p> </div> </div> </div> </div> </div> <div> <div> <div> <div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:72.0pt"> <span style="font-size:11.0pt"><img border="0" width="605" height="432" style="width:6.302in;height:4.5in" id="_x0000_i1028" src="https://groups.google.com/group/cp2k/attach/720b7b929b503/NVE.png?part=0.1&view=1"><o:p></o:p></span></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:72.0pt"> <span style="font-size:11.0pt">NVT:<o:p></o:p></span></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:72.0pt"> <span style="font-size:11.0pt"><img border="0" width="613" height="432" style="width:6.3854in;height:4.5in" id="_x0000_i1027" src="https://groups.google.com/group/cp2k/attach/720b7b929b503/NVT.png?part=0.4&view=1"><o:p></o:p></span></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:72.0pt"> <span style="font-size:11.0pt">NPT:<o:p></o:p></span></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:72.0pt"> <span style="font-size:11.0pt"><img border="0" width="613" height="432" style="width:6.3854in;height:4.5in" id="_x0000_i1026" src="https://groups.google.com/group/cp2k/attach/720b7b929b503/NPT.png?part=0.3&view=1"><o:p></o:p></span></p> </div> </div> </div> </div> </div> <div> <div> <div> <div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:72.0pt"> <span style="font-size:11.0pt">The simulations of the water box were small test cases, but I also need to know what behaviour to expect when going to my production runs with my systems of interest...<o:p></o:p></span></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:72.0pt"> <span style="font-size:11.0pt"> <o:p></o:p></span></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:72.0pt"> <span style="font-size:11.0pt">Additionally, in the NPT ensemble, the pressure seems to fluctuate very strongly from -20000 to +20000 bar which seems to be quite excessive to me. Is this normal? I've also seen this asked in <a href="https://groups.google.com/g/cp2k/c/gxA8r6wmyF4/m/k8QT7iglAAAJ" target="_blank"> another post</a>, but that person did not get a reply unfortunately....<o:p></o:p></span></p> </div> </div> </div> </div> </div> <div> <div> <div> <div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:72.0pt"> <span style="font-size:11.0pt"><img border="0" width="600" height="432" style="width:6.25in;height:4.5in" id="_x0000_i1025" src="https://groups.google.com/group/cp2k/attach/720b7b929b503/Pressure.png?part=0.2&view=1"><o:p></o:p></span></p> </div> </div> </div> </div> </div> <div> <div> <div> <div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:72.0pt"> <span style="font-size:11.0pt"> <o:p></o:p></span></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:72.0pt"> <span style="font-size:11.0pt">Thanks in advance for the help!<o:p></o:p></span></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:72.0pt"> <span style="font-size:11.0pt"> <o:p></o:p></span></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:72.0pt"> <span style="font-size:11.0pt">Kind regards,<o:p></o:p></span></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:72.0pt"> <span style="font-size:11.0pt">Léon<o:p></o:p></span></p> </div> </div> </div> </div> </div> <div> <div> <div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:72.0pt"> <span style="font-size:11.0pt">-- <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 cp2k+uns...@googlegroups.com.<br> To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/76d3e94b-4a67-4fcc-8861-78a566ebb107n%40googlegroups.com?utm_medium=email&utm_source=footer" target="_blank"> https://groups.google.com/d/msgid/cp2k/76d3e94b-4a67-4fcc-8861-78a566ebb107n%40googlegroups.com</a>.<o:p></o:p></span></p> </div> </div> </div> </div> </blockquote> </div> </blockquote> </div> <p class="MsoNormal" style="margin-left:36.0pt"><span style="font-size:11.0pt">-- <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/e6d4002e-91a7-45cd-817c-adc66616d1d6n%40googlegroups.com?utm_medium=email&utm_source=footer"> https://groups.google.com/d/msgid/cp2k/e6d4002e-91a7-45cd-817c-adc66616d1d6n%40googlegroups.com</a>.<o:p></o:p></span></p> </div> </div> </div> </body> </html> <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/ZRAP278MB082736BE0B4ECAEF2E055D0EF4752%40ZRAP278MB0827.CHEP278.PROD.OUTLOOK.COM?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/ZRAP278MB082736BE0B4ECAEF2E055D0EF4752%40ZRAP278MB0827.CHEP278.PROD.OUTLOOK.COM</a>.<br />