Thank you for your suggestion. After increasing EPS_PGF_ORB to 1E-6, the calculation proceeded smoothly.<br /><br /><div class="gmail_quote"><div dir="auto" class="gmail_attr">在2026年1月22日星期四 UTC+8 17:23:28<Augustin Bussy> 写道：<br/></div><blockquote class="gmail_quote" style="margin: 0 0 0 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;"><div>Hi Yisichi,</div><div>Your calculation is most likely not stuck, but proceeding. Hybrid DFT in the condensed phase can be very expensive, even when using ADMM and non-diffuse basis sets like pob-DZVP-rev2. Given an optimized input file, the only solution is to used more resources.<br><br>That being said, your input in not perfect. In particular, you are using and <b>extremely tight</b> value for EPS_PGF_ORB. The default value if 1E-5, and particularly accurate calculations can used 1E-6. I recommend you to change it. The MIN_PAIR_LIST_RADIUS option that you commented out is a much more efficient solution to the "Kohn Sham matrix occupancy warning": use this instead. For performance, make sure to also optimize the MAX_MEMORY keyword to the specifications of your computer. See the "<span style="font-family:sans-serif" role="presentation" dir="ltr">Hartree</span><span style="font-family:sans-serif" role="presentation" dir="ltr">−</span><span style="font-family:sans-serif" role="presentation" dir="ltr">Fock and Hybrid Density Functional Theory</span>" section of our latest paper (<a href="https://pubs.acs.org/doi/full/10.1021/acs.jpcb.5c05851" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=zh-CN&q=https://pubs.acs.org/doi/full/10.1021/acs.jpcb.5c05851&source=gmail&ust=1769176280220000&usg=AOvVaw3NlfIKRjvIJxGpdkhtmy2M">https://pubs.acs.org/doi/full/10.1021/acs.jpcb.5c05851</a>) to understand the meaning of these keywords.</div><div><br></div><div>Best,</div><div>Augustin</div><br><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Wednesday, 21 January 2026 at 17:41:23 UTC+1 <a href data-email-masked rel="nofollow">yis...@163.com</a> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div style="box-sizing:border-box;padding:0px;color:rgb(128,0,128);font-family:quote-cjk-patch,Inter,system-ui,-apple-system,BlinkMacSystemFont,"Segoe UI",Roboto,Oxygen,Ubuntu,Cantarell,"Open Sans","Helvetica Neue",sans-serif;font-variant-ligatures:no-contextual"><div style="color:var(--dsw-alias-label-primary)"><p style="margin:0px 0px 16px"><span style="color:var(--dsw-alias-label-primary)">To the CP2K Development Team,</span></p></div></div><div style="box-sizing:border-box;padding:0px;color:rgb(128,0,128);font-family:quote-cjk-patch,Inter,system-ui,-apple-system,BlinkMacSystemFont,"Segoe UI",Roboto,Oxygen,Ubuntu,Cantarell,"Open Sans","Helvetica Neue",sans-serif;font-variant-ligatures:no-contextual"><div style="color:var(--dsw-alias-label-primary)"><p style="margin:16px 0px">Hello,</p><p style="margin:16px 0px">I am a beginner in DFT calculations and am currently learning to use CP2K for calculations with the PBE0 hybrid functional. My research involves a Mo-doped NiFe-LDH system. I have constructed a 6×8×3 supercell containing a total of 720 atoms. Given the significant presence of Ni, Fe, and O in the system, when using the GTH pseudopotential with the DZVP-MOLOPT-SR basis set, convergence was only achieved with very high values of CUTOFF and REL_CUTOFF. Therefore, I intend to use the GAPW method to reduce computational cost and improve accuracy.</p><p style="margin:16px 0px">Currently, I am employing the pob-TZVP-rev2 basis set for all-electron calculations, along with pob-DZVP-rev2 as the ADMM basis set to accelerate the hybrid functional calculation. During the calculation, I first performed a pure functional calculation using the pob-TZVP-rev2 basis set. Then, I used the wavefunction generated from this calculation as the initial guess for the hybrid functional calculation. Unfortunately, the hybrid functional calculation does not seem to proceed—after running continuously for 4 hours, the calculation remains stuck at the interface shown in the attached image and does not appear to enter the self-consistent iteration process.</p><p style="margin:16px 0px">The supercomputer I am using is equipped with Intel Xeon Platinum Processors (3rd Gen, 64 cores) and 256 GB of memory. I would like to ask the developers for their insight into the following:</p><ol start="1" style="margin:16px 0px;padding-left:18px"><li><p style="margin:0px">What might be causing this issue?</p></li><li style="margin-top:6px"><p style="margin:0px">How should hybrid functional calculations be conducted for such a system?</p></li><li style="margin-top:6px"><p style="margin:0px">Could you recommend a more suitable basis set for this system?</p></li></ol><p style="margin:16px 0px">For your reference, I have attached the input and output files from my CP2K calculation at the end of this email.</p><p style="margin:16px 0px">Thank you for your attention and assistance.</p><p style="margin:16px 0px 0px">Best regards,</p><p style="margin:16px 0px 0px">yisichi</p><p style="margin:16px 0px 0px"><img alt="微信图片_20260122004005_251_55.png" width="534px" height="314px" src="https://groups.google.com/group/cp2k/attach/1da344bc278f5/%E5%BE%AE%E4%BF%A1%E5%9B%BE%E7%89%87_20260122004005_251_55.png?part=0.1&view=1"><br></p></div></div></blockquote></div></blockquote></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 visit <a href="https://groups.google.com/d/msgid/cp2k/36dc2999-ecf5-4c85-9343-f2a8dec6df4cn%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/36dc2999-ecf5-4c85-9343-f2a8dec6df4cn%40googlegroups.com</a>.<br />