<div dir="ltr"><div class="gmail_default" style="font-family:arial,sans-serif">Dear Nanchen,</div><div class="gmail_default" style="font-family:arial,sans-serif"><br></div><div class="gmail_default" style="font-family:arial,sans-serif">Thank you very much for your sample input and result. I used a similar setting as yours to run GAPW for 1 O atom with the def2-SVP basis and 4 O atoms with the aug-cc-pVTZ basis. The def2-SVP calculation converged, but it took almost 500 SCF steps, and there's a weird oscillation of Convergence between 1e-4 and 1e-3. The aug-cc-pVTZ calculation converges more slowly and can't converge within 500 steps. Do you have any idea what causes the slow convergence and how to improve it? Thank you! Below is my CP2K version:</div><div class="gmail_default" style="font-family:arial,sans-serif"> CP2K version 2026.1 (Development Version)<br> Source code revision 7b94146043<br> cp2kflags: omp libxc elpa parallel scalapack<br> compiler: GCC version 14.3.0<br> compiler options:<br> -I /global/homes/h/hengyuan/hys_project/cp2k/build/src/start -I /glo<br> bal/homes/h/hengyuan/hys_project/cp2k/src -I /global/homes/h/hengyua<br> n/hys_project/cp2k/build/src -I /global/homes/h/hengyuan/hys_project<br> /cp2k/src/base -I /global/homes/h/hengyuan/hys_project/cp2k/src/comm<br> on -I /global/homes/h/hengyuan/hys_project/cp2k/src/motion -I /globa<br> l/homes/h/hengyuan/hys_project/cp2k/src/dbm -I /global/homes/h/hengy<br> uan/hys_project/cp2k/build/src/mod_files -I /global/homes/h/hengyuan<br> /hys_project/cp2k/spack/spack/opt/spack/view/include -I /global/home<br> s/h/hengyuan/hys_project/cp2k/spack/spack/opt/spack/__spack_path_pla<br> ceholder__/__spack_path_placeholder__/__spack_pa/linux-zen3/elpa-202<br> 6.02.001-s5vz6kbd3dbethregansxz4e5tlqsafv/include/elpa_openmp-2026.0<br> 2.001 -I /global/homes/h/hengyuan/hys_project/cp2k/spack/spack/opt/s<br> pack/__spack_path_placeholder__/__spack_path_placeholder__/__spack_p<br> a/linux-zen3/elpa-2026.02.001-s5vz6kbd3dbethregansxz4e5tlqsafv/inclu<br> de/elpa_openmp-2026.02.001/modules -I /global/cfs/cdirs/mp54/hys/cp2<br> k/spack/spack/opt/spack/__spack_path_placeholder__/__spack_path_plac<br> eholder__/__spack_pa/linux-zen3/mpich-4.3.2-6zwkkcll2y5mdmgz3at7bbns<br> liyxgfbz/include -march=znver3 -mmmx -mpopcnt -msse -msse2 -msse3 -m<br> ssse3 -msse4.1 -msse4.2 -mavx -mavx2 -msse4a -mno-fma4 -mno-xop -mfm<br> a -mno-avx512f -mbmi -mbmi2 -maes -mpclmul -mno-avx512vl -mno-avx512<br> bw -mno-avx512dq -mno-avx512cd -mno-avx512vbmi -mno-avx512ifma -mno-<br> avx512vpopcntdq -mno-avx512vbmi2 -mno-gfni -mvpclmulqdq -mno-avx512v<br> nni -mno-avx512bitalg -mno-avx512bf16 -mno-avx512vp2intersect -mno-3<br> dnow -madx -mabm -mno-cldemote -mclflushopt -mclwb -mclzero -mcx16 -<br> mno-enqcmd -mf16c -mfsgsbase -mfxsr -mno-hle -msahf -mno-lwp -mlzcnt<br> -mmovbe -mno-movdir64b -mno-movdiri -mmwaitx -mno-pconfig -mpku -mp<br> rfchw -mno-ptwrite -mrdpid -mrdrnd -mrdseed -mno-rtm -mno-serialize<br> -mno-sgx -msha -mshstk -mno-tbm -mno-tsxldtrk -mvaes -mno-waitpkg -m<br> wbnoinvd -mxsave -mxsavec -mxsaveopt -mxsaves -mno-amx-tile -mno-amx<br> -int8 -mno-amx-bf16 -mno-uintr -mno-hreset -mno-kl -mno-widekl -mno-<br> avxvnni -mno-avx512fp16 -mno-avxifma -mno-avxvnniint8 -mno-avxneconv<br> ert -mno-cmpccxadd -mno-amx-fp16 -mno-prefetchi -mno-raoint -mno-amx<br> -complex -mno-avxvnniint16 -mno-sm3 -mno-sha512 -mno-sm4 -mno-apxf -<br> mno-usermsr --param=l1-cache-size=32 --param=l1-cache-line-size=64 -<br> -param=l2-cache-size=512 -mtune=znver3 -g -O3 -O3 -Wno-deprecated-de<br> clarations -Wno-maybe-uninitialized -Wuninitialized -Wuse-without-on<br> ly -std=f2008 -fPIE -ffree-form -fimplicit-none -fno-omit-frame-poin<br> ter -fbacktrace -fallow-argument-mismatch -funroll-loops -fopenmp -f<br> preprocessed -J src -fpre-include=/usr/include/finclude/math-vector-<br> fortran.h</div></div><img src="https://email-signature-image.com/signature.gif?u=1605853&e=1110047184&v=a6c233c70bd0f8d5b2c5dbc0601211e9979662c174fbd9fe19c7374b90b5a5b6" alt="" width="0" height="0" style="width:2px;max-height:0;overflow:hidden"><br><div class="gmail_quote gmail_quote_container"><div dir="ltr" class="gmail_attr">On Sat, May 2, 2026 at 6:34 AM Nan Dong <<a href="mailto:chen.east.7@gmail.com">chen.east.7@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Hi Hengyuan, <div><br></div><div>I ran the GAPW calculations with def2-SVP for O, and it converged (within 300 SCF step) </div><div>Based on your def2-SVP.in, I changed something in SCF section, and attached input and output here. </div><div>Hope this will help. </div><div><br></div><div>Best,</div><div>Nanchen </div><div><br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">在2026年4月30日星期四 UTC+2 21:12:25<Hengyuan Shen> 写道：<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Dear Marcella,<div><br></div><div>Thank you very much for your suggestions. I switched to Fermi Dirac smearing, and the GPW calculation finishes normally now, with the correct total charge density. However, once I switch to GAPW calculation, the SCF still failed to converge, for either aug-cc-pVTZ or def2-SVP basis. I also tried to use all electron basis on only one O atom, but the convergence issue still persist. Do you have any suggestions on converging SCF for the GAPW calculations? Thank you.</div><div><br></div><div>Best regards,</div><div>Hengyuan<br><br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">在2026年4月30日星期四 UTC-7 04:44:14<Marcella Iannuzzi> 写道：<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Dear Hengyuan,<div><br></div><div>The smearing method that you are using is obviously not working properly, also in the case of GPW.</div><div>It is not normal that the occupation is smeared down to the first MO. You can add as many added MOs as you want, this does not change.</div><div><p style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo;font-size-adjust:none;font-kerning:auto;font-variant-alternates:normal;font-variant-ligatures:normal;font-variant-numeric:normal;font-variant-east-asian:normal;font-feature-settings:normal;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"> *** WARNING in qs_mo_occupation.F:121 :: Fermi-Dirac smearing includes ***</span></p><p style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo;font-size-adjust:none;font-kerning:auto;font-variant-alternates:normal;font-variant-ligatures:normal;font-variant-numeric:normal;font-variant-east-asian:normal;font-feature-settings:normal;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"> *** the first MO</span></p><p style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo;font-size-adjust:none;font-kerning:auto;font-variant-alternates:normal;font-variant-ligatures:normal;font-variant-numeric:normal;font-variant-east-asian:normal;font-feature-settings:normal;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><br></span></p></div><div>Indications of a wrong behaviour are also the oscillating convergence and the non-zero total charge density </div><div><p style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo;font-size-adjust:none;font-kerning:auto;font-variant-alternates:normal;font-variant-ligatures:normal;font-variant-numeric:normal;font-variant-east-asian:normal;font-feature-settings:normal;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"> Total charge density on r-space grids: -0.0006221769</span></p>
<p style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo;font-size-adjust:none;font-kerning:auto;font-variant-alternates:normal;font-variant-ligatures:normal;font-variant-numeric:normal;font-variant-east-asian:normal;font-feature-settings:normal;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"> Total charge density g-space grids: -0.0006221769</span></p><p style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo;font-size-adjust:none;font-kerning:auto;font-variant-alternates:normal;font-variant-ligatures:normal;font-variant-numeric:normal;font-variant-east-asian:normal;font-feature-settings:normal;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><br></span></p><p style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo;font-size-adjust:none;font-kerning:auto;font-variant-alternates:normal;font-variant-ligatures:normal;font-variant-numeric:normal;font-variant-east-asian:normal;font-feature-settings:normal;color:rgb(0,0,0)">My suggestion is to run the GPW calculation with the Fermi-Dirac smearing to confirm that one problem is the choice of the smearing method. </p><p style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo;font-size-adjust:none;font-kerning:auto;font-variant-alternates:normal;font-variant-ligatures:normal;font-variant-numeric:normal;font-variant-east-asian:normal;font-feature-settings:normal;color:rgb(0,0,0)"><br></p>Regards</div><div>Marcella</div><div><br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Thursday, April 30, 2026 at 5:26:31 AM UTC+2 <a rel="nofollow">hengyu...@gmail.com</a> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Dear CP2K community,<div><br></div><div>I tried to calculate the XPS binding energy of some O atoms on a Pt surface. I used a big all electron basis for the interested O atoms, a smaller all electron basis for other O atoms, and GTH-MOLOPT basis for Pt with pseudopotential. However, the SCF calculation could not converge. I then tried the ground state calculation with a similar GPAW method and GPW method. For GPAW method, the SCF convergence also failed, and gave me an absurd energy. The GPW method converged after about 100 steps. Strangely, both calculations showed warnings about Fermi-Dirac smearing, although I used METHFESSEL_PAXTON smearing and 1000 ADDED_MO. Any idea on this issue? Any suggestions would be appreciated. Thank you.</div><div><br></div><div>Best regards,</div><div>Hengyuan</div></blockquote></div></blockquote></div></blockquote></div>
<p></p>
-- <br>
You received this message because you are subscribed to a topic in the Google Groups "cp2k" group.<br>
To unsubscribe from this topic, visit <a href="https://groups.google.com/d/topic/cp2k/6--LsBoEEZw/unsubscribe" target="_blank">https://groups.google.com/d/topic/cp2k/6--LsBoEEZw/unsubscribe</a>.<br>
To unsubscribe from this group and all its topics, send an email to <a href="mailto:cp2k+unsubscribe@googlegroups.com" target="_blank">cp2k+unsubscribe@googlegroups.com</a>.<br>
To view this discussion visit <a href="https://groups.google.com/d/msgid/cp2k/84d452dd-bb63-4a23-b9b6-24bf2dcdede7n%40googlegroups.com?utm_medium=email&utm_source=footer" target="_blank">https://groups.google.com/d/msgid/cp2k/84d452dd-bb63-4a23-b9b6-24bf2dcdede7n%40googlegroups.com</a>.<br>
</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/CAB5xK-v_Z9VmV1ZdTjsDmQCWn2piQ_3HwtyHmR8GjQOX0U0pJg%40mail.gmail.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/CAB5xK-v_Z9VmV1ZdTjsDmQCWn2piQ_3HwtyHmR8GjQOX0U0pJg%40mail.gmail.com</a>.<br />