Page 1 of 1

WF convergence issue

Posted: Fri Mar 28, 2025 1:19 pm
by aldo_romero1

I have been trying to converge the wavefunction of a system I am dealing with. The energy goes down for few iterations as:
N E dE d eps ncg rms rms(c)
DAV: 1 -0.217186945372E+04 -0.21719E+04 0.97968E+13 18496 0.158E+12
DAV: 2 -0.979649664299E+02 0.20739E+04 -0.40086E+03 20968 0.447E+02
DAV: 3 -0.166678657388E+03 -0.68714E+02 -0.63448E+02 26664 0.969E+01
DAV: 4 -0.171380884559E+03 -0.47022E+01 -0.46862E+01 27940 0.209E+01
DAV: 5 -0.171568694410E+03 -0.18781E+00 -0.18775E+00 26256 0.411E+00 0.348E+03
RMM: 6 0.700459931069E+04 0.71762E+04 -0.35924E+08 29412 0.569E+06 0.425E+05

But then I get the error:

| ERROR FEXCF: supplied exchange-correlation table |
| is too small, maximal index : 6025 |
| |
| ----> I REFUSE TO CONTINUE WITH THIS SICK JOB ... BYE!!! <---- |

Interatomic distances are fine, cutoff is fine, K-mesh grid seems OK, though I have tried different values.

Any suggestions on how to tackle this problem?

Thanks!

PS Let me know if you need any further information.


Re: WF convergence issue

Posted: Fri Mar 28, 2025 2:32 pm
by marie-therese.huebsch

Hi, thanks for raising the issue.

Can you please share the files for a minimal reproducible example in accordance with the forum guidelines?

Cheers,
Marie-Therese


Re: WF convergence issue

Posted: Fri Mar 28, 2025 3:18 pm
by aldo_romero1

I am attaching a tgz file with all necessary information to reproduce my problem.

Thanks

-aldo.


Re: WF convergence issue

Posted: Mon Mar 31, 2025 9:46 am
by marie-therese.huebsch

Hi, thank you for sharing the files. I see that you have chosen NCORE=32. Though likely inefficient, that should be fine, since you are running on 32 ranks. Still some issue occurs due to the large value of NCORE. You can see already the first electronic step complains about a nonhermitian subspace matrix. I will look into the reason.

As a fix for your calculation, I suggest running a few steps (NELM=10) with NCORE= 2, 4, 8, 16 and select the one most efficient. For this test, particularly the number of ranks, PREC, ENCUT should remain the one you will use in the final calculation.

Could you please let me know if this helps?

Best regards,
Marie-Therese


Re: WF convergence issue

Posted: Tue Apr 01, 2025 2:05 pm
by aldo_romero1

Marie

it works for this case. I optimized the NCORE value, it run, it converged and it runs faster. Now I will translate that learning into the thousand of materials I am dealing with and I will get back to you to report. Thanks for your help!. Really appreciated.


Re: WF convergence issue

Posted: Wed Apr 02, 2025 7:00 am
by marie-therese.huebsch

Great! Glad I could help.