[Wien] question about mixing scf calculations with different NUME in param.inc at compilation
Tomas Kana
kana at seznam.cz
Tue Dec 10 14:21:25 CET 2024
Dear Wien2k community,
I am calculating a Mg23Al30 structure as in input for Phonopy code.
I have build a supercell with 636 atoms in the unit cell.
case.scf gives a total number of 6720 electrons.
For the 14 scf iterations I used proper version
wien2k/21.1.matmax500k (the name
given by our system admins) that was compiled using
NMATMAX= 500000
NUME= 50000
in param.inc.
After running out of computer time, I continued the scf calculation
but I did not pay attention and used the version wien2k/21.1 with
NMATMAX= 50000
NUME= 6000
in param.inc
I know from previous posts that
NUME > "number of electrons"/2 + 20-30%
When analyzing the file case.scf
it has two warnings:
:WARN : WARNING: RKmax reduced due to NMATMAX
:WARN : WARNING: NE limited by NUME in param.inc
I do not mind the first warning, it only reduced RMT times Kmax
from 7.0 to 5.94.
I tried to analyze the second warning and found that
after 14 "right" cycles the interstitial (grep :NPC case.scf)
charge decreased from 742 to 711 (I give only integers, not the
numbers after decimal points).
I now plan tostop the calculations and use the density
calculated using wien2k/21.1 as starting density for calculation
using wien2k/21.1.matmax500k.
Is this correct?
Thank you for your answer.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://zeus.theochem.tuwien.ac.at/pipermail/wien/attachments/20241210/72b78d9a/attachment.htm>
More information about the Wien
mailing list