Dear Prof. Juerg Hutter,<div>Thank you very much for reply,</div><div>I have increased
EPS_SCHWARZ_FORCES to 1.0e-4, the overall optimization efficiency is much faster.</div><div>However I am not sure if it is safe option, is there any rule of thumb for how to set up EPS_SCHWARZ_FORCES ?</div><div>Thank you in advance.</div><div><br>Best,<br>Geng</div><div><div><br><div><br></div><div><div> &SCREENING</div><div> EPS_SCHWARZ 1.0E-6</div><div> SCREEN_ON_INITIAL_P T</div><div> EPS_SCHWARZ_FORCES 1.0E-4</div><div> SCREEN_P_FORCES T</div><div> &END SCREENING</div><div><br></div><br></div></div></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">在2021年1月5日星期二 UTC-8 上午1:53:20<jgh> 写道：<br/></div><blockquote class="gmail_quote" style="margin: 0 0 0 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">Hi
<br>
<br>Derivative ERI's are only used once (in the gradient calculation)
<br>are therefore not stored. Typical cost is 4x for the normal ERIs.
<br>The advantage for the derivative ERI calculation is, that it
<br>can be screened much more efficiently.
<br>Have a look at the section
<br>
<br>CP2K_INPUT / FORCE_EVAL / DFT / XC / HF / SCREENING
<br>
<br>EPS_SCHWARZ_FORCES
<br>SCREEN_P_FORCES
<br>
<br>Optimal settings of these Keywords can make a huge difference.
<br>
<br>In addition, timings are very dependent on the ADMM basis and
<br>the radial cutoff if you are using a truncated operator.
<br>
<br>regards
<br>
<br>Juerg Hutter
<br>--------------------------------------------------------------
<br>Juerg Hutter Phone : <a href="tel:+41%2044%20635%2044%2091" value="+41446354491" target="_blank" rel="nofollow">++41 44 635 4491</a>
<br>Institut für Chemie C FAX : <a href="tel:+41%2044%20635%2068%2038" value="+41446356838" target="_blank" rel="nofollow">++41 44 635 6838</a>
<br>Universität Zürich E-mail: <a href data-email-masked rel="nofollow">h...@chem.uzh.ch</a>
<br>Winterthurerstrasse 190
<br>CH-8057 Zürich, Switzerland
<br>---------------------------------------------------------------
<br>
<br>-----<a href data-email-masked rel="nofollow">c...@googlegroups.com</a> wrote: -----
<br>To: "cp2k" <<a href data-email-masked rel="nofollow">c...@googlegroups.com</a>>
<br>From: "Sun Geng"
<br>Sent by: <a href data-email-masked rel="nofollow">c...@googlegroups.com</a>
<br>Date: 01/04/2021 09:31PM
<br>Subject: [CP2K:14472] Number of sph. DERIV's calculated on the fly
<br>
<br>Dear CP2K users,
<br>
<br>I am using CP2K with the ADMM method for structure optimization.
<br>The system is an inorganic bulk with around 100 atoms in a cubic box 12x12x12 angstrom^3.
<br>I am able to set up the calculation resulting smooth SCF convergence.
<br>The "Number of sph. ERI's calculated on the fly" is 0 and hence the SCF converges quickly. The first SCF takes about 1800 seconds and the following SCF steps take only an impressive 20 seconds. In the end, CP2K can converge the SCF in approximate 2500 seconds.
<br>However, I notice that in between the optimization steps, the derivatives of integrals are still calculated on the fly (see the output below), and this step appears to be very slow. CP2K did not print how long it took to compute the derivatives, but I can tell by checking the updating time in the output. My estimation is that it takes more than 2 hours to calculate the derivatives.
<br>
<br> HFX_MEM_INFO| Number of cart. primitive DERIV's calculated: 11208181891200
<br> HFX_MEM_INFO| Number of sph. DERIV's calculated: 140954718240
<br> HFX_MEM_INFO| Number of sph. DERIV's stored in-core: 0
<br> HFX_MEM_INFO| Number of sph. DERIV's calculated on the fly: 140954718240
<br> HFX_MEM_INFO| Total memory consumption DERIV's RAM [MiB]: 0
<br> HFX_MEM_INFO| Whereof max-vals [MiB]: 1
<br> HFX_MEM_INFO| Total compression factor DERIV's RAM: 0.00
<br>
<br>
<br>Finally, the overall efficiency of CP2K for optimization is much slower than that of VASP code for this system (both of them used HSE06 functional). In VASP, although an SCF step takes 200 seconds and SCF converged in 15 steps (resulting in 3000 seconds per optimization step), the forces are calculated with negligible cost in VASP.
<br>
<br>I wonder are there any tricks in the cp2k inputs that can speed up the derivatives of the integrals? Or any pieces of advice that I can use in using CP2K/ADMM methods?
<br>My personal impression is that CP2K/ADMM is not the optimal choice (in terms of optimization efficiency) for a small system (maybe within 150~200 atoms, I did not test though) compared with VASP code. Maybe, for a larger system, the overhead of calculating forces in CP2K/ADMM is relatively small, and in that case, CP2K/ADMM will be more efficient than VASP? Please let me know if this is a realistic point for choosing VASP and CP2K.
<br>
<br>Thank you very much in advance!
<br>
<br>Best
<br>Geng
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br> --
<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 data-email-masked rel="nofollow">cp...@googlegroups.com</a>.
<br> To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/e8bc30d7-297c-4217-839d-29a86a8f38f3n%40googlegroups.com" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=zh-CN&q=https://groups.google.com/d/msgid/cp2k/e8bc30d7-297c-4217-839d-29a86a8f38f3n%2540googlegroups.com&source=gmail&ust=1609948976939000&usg=AFQjCNFYVib5K67GctkRrQPikamormim2g">https://groups.google.com/d/msgid/cp2k/e8bc30d7-297c-4217-839d-29a86a8f38f3n%40googlegroups.com</a>.
<br>
<br></blockquote></div>