<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=gb2312">
<style type="text/css" style="display:none;"> P {margin-top:0;margin-bottom:0;} </style>
</head>
<body dir="ltr">
<div class="elementToProof" style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
Hi,</div>
<div class="elementToProof" style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div class="elementToProof" style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
what you observe is somewhat expected. The limiting factor in your B3LYP calculation is the Hartree-Fock exchange, which computational performance depends strongly on the choice of EPS_PGF_ORB, EPS_SCHWARZ, CUTOFF_RADIUS and MAX_MEMORY.</div>
<div class="elementToProof" style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div class="elementToProof" style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
In an ideal CP2K run, when performing HFX calculations, all of the ERIs are initially calculated and stored in memory (during the first SCF step). This allows for much faster timings for the next SCF steps, because ERIs are simply fetched from memory. The keyword
MAX_MEMORY from the &HF section controls how much memory is dedicated to storing the ERIs (per MPI rank). If there are more ERIs than available memory, the rest is computed on the fly, at each SCF step. This is linked to your first question: when the memory
stabilizes, it means that you have reached that MAX_MEMORY upper limit, and the new ERIs are not stored. You should adapt the value of MAX_MEMORY to your hardware.</div>
<div class="elementToProof" style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div class="elementToProof" style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
For your second and third question: HFX performance in CP2K depends a lot on the choice of EPS_SCHWARZ, EPS_PGF_ORB and CUTOFF_RADIUS. These parameters control how many ERIs need to be calculated, and going from 10E-8 to 10E-10 EPS_SCHWARZ while increasing
the CUTOFF_RADIUS will lead to a large number of additional ERIs to be calculated (and CP2K might seems stuck because there is so much more to do). In most cases, I would say that EPS_SCHWARZ = 1.0E-8, EPS_PGF_ORB =1.0E-5 - 1.0E-6 and CUTOFF_RADIUS = 6 Angstrom
is accurate enough and efficient.</div>
<div class="elementToProof" style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div class="elementToProof" style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
You can also consider the ADMM method to further accelerate your calculation (<a href="https://pubs.acs.org/doi/full/10.1021/ct1002225" id="LPlnk354856">https://pubs.acs.org/doi/full/10.1021/ct1002225</a>).</div>
<div class="elementToProof" style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div class="elementToProof" style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
I hope that helps.</div>
<div class="elementToProof" style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
Best,</div>
<div class="elementToProof" style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
Augustin<br>
</div>
<div id="appendonsend"></div>
<hr style="display:inline-block;width:98%" tabindex="-1">
<div id="divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" style="font-size:11pt" color="#000000"><b>From:</b> cp2k@googlegroups.com <cp2k@googlegroups.com> on behalf of Roosevelt Gildner <roosevelt.gildner12@gmail.com><br>
<b>Sent:</b> Saturday, September 23, 2023 4:29 PM<br>
<b>To:</b> cp2k <cp2k@googlegroups.com><br>
<b>Subject:</b> [CP2K:19271] On the time&memory consumption of ERI calculation in CP2K hybrid functional single point</font>
<div> </div>
</div>
<div>
<div>Dear CP2K team,</div>
<div><br>
</div>
<div>I am attempting to calculate the single point energy of a huge periodic cluster (near 5k atoms, just contains C H O N, cell size 30+ angstrom, only Gamma point). Due to the demand of following analysis, I used B3LYP hybrid functional (with D3(BJ)) and
6-31G(d) basis for GAPW calculation.</div>
<div>I have prepared enough memory space to do this (TB scale). Meanwhile I want as accurate result (especially for DOS analysis, need unoccupied orbitals) as possible, but not too excessive time comsuption.</div>
<div>I have tried several parameter sets and found when the EPS_SCHWARZ and CUTOFF_RADIUS is respectively over 8 (exponent) and 10 angstrom, the time comsuption significant increases.</div>
<div>Below are the two typical examples:</div>
<div>A <br>
<p lang="en-US" style="margin:0in; font-family:Calibri; font-size:11pt">EPS_SCHWARZ=1E-8.</p>
<p lang="en-US" style="margin:0in; font-family:Calibri; font-size:11pt">SCREEN_ON_INITIAL_P=T</p>
<p lang="en-US" style="margin:0in; font-family:Calibri; font-size:11pt">CUTOFF_RADIUS=8.0.</p>
<p lang="en-US" style="margin:0in; font-family:Calibri; font-size:11pt">The total time comsuption is about 25 hours, memory peak 4.2TB(not on-the-fly). The HPC monitor tells me that the memory linear rised to the stable level and soon the first SCF step completed.(total
SCF=24 steps)</p>
<p lang="en-US" style="margin:0in; font-family:Calibri; font-size:11pt"><img alt="微信图片_20230923213258.png" width="532px" height="249px" data-outlook-trace="F:1|T:1" src="cid:0b2686b7-490f-4ebc-bc81-3a12d5e2354b"><br>
</p>
<p lang="en-US" style="margin:0in; font-family:Calibri; font-size:11pt">B</p>
<p lang="en-US" style="margin:0in; font-family:Calibri; font-size:11pt"></p>
<p lang="en-US" style="margin:0in; font-family:Calibri; font-size:11pt">EPS_SCHWARZ=1E-10.</p>
<p lang="en-US" style="margin:0in; font-family:Calibri; font-size:11pt">SCREEN_ON_INITIAL_P=T</p>
<p lang="en-US" style="margin:0in; font-family:Calibri; font-size:11pt">CUTOFF_RADIUS=10.0</p>
</div>
<div>The memory comsuption stablized at about 7TB. But it has been 1 day (now already 5 days) since the memory rised to the stable level, the first SCF step is not completed yet.</div>
<div><img alt="微信图片_20230922132420-2.png" width="532px" height="283px" data-outlook-trace="F:1|T:1" src="cid:f9fd47dd-12b0-4dcb-99f9-5f6843f721de"></div>
<div><br>
</div>
<div>Here are my confusions:</div>
<div>1 Whether the memory rising to the stable level means the ERI calculation is completed? if yes, why the example B does not complete the first SCF step even after 1-5 days? (the memory used is not 100%, so i believe it is not on-the-fly mode)And why do
example A and B take such different time comsuption?<br>
</div>
<div>2 I am afraid the CP2K mission B got stuck. but I am not sure. Is there any method to judge whether the mission is normal? And need I continue waiting?</div>
<div>3 Which parameter(s) is the key factor(s) affecting the time comsuption of hybrid functional calculation? EPS_DEFAULT? EPS_PGF_ORB? EPS_SCHWARZ? CUTOFF_RADIUS?<br>
</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/57a56625-e467-4a50-ab54-665e7df6d445n%40googlegroups.com?utm_medium=email&utm_source=footer">
https://groups.google.com/d/msgid/cp2k/57a56625-e467-4a50-ab54-665e7df6d445n%40googlegroups.com</a>.<br>
</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/GV0P278MB0370A4DAC28DDFF84DB2F29ED3FCA%40GV0P278MB0370.CHEP278.PROD.OUTLOOK.COM?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/GV0P278MB0370A4DAC28DDFF84DB2F29ED3FCA%40GV0P278MB0370.CHEP278.PROD.OUTLOOK.COM</a>.<br />