Well, at the end it is exactly as I said:
Your manual RMT settings are very bad. One of the rules is, that
smallest and largest RMTs must not be too different.
With your spheres you get "effective" R(Si)Kmax of more than 12 and this
gives numerical linear dependency.
With your struct file: RMT(N)=1.2 (Si)=2.1
init
set RKmax=7
run
----> gives the semicore error you describe due to some ghostbands.
init
run
save rkm6
set RKmax=7
run
---> runs through and converges.
However, increasing RKmax to 8 gives ghost bands again.
---------------------------------
setrmt
cp case.struct_setrmt case.struct ( similar RMTs for Si and N, around 1.6)
init -ecut -8 (to avoid core leakage)
set rkmax=7
run
---> converges without problems. One can increase Rkmax further to 8 and
even 9.
----------------------------------
PS: These ghostbands are located in the interstital, give no qtl-b
errors. Once such a state is taken into the density, you get these
"select"-errors.
Am 14.11.2023 um 19:19 schrieb hajar.nejatipoor--- via Wien:
Dr. Blaha
sometimes, *semicore* error appears in iteration3, sometime in 6, and
... (with changing rmts).
I tried with the struct attached here and the default init_lapw. After
finalized initialization, I changed RKm=7 (and Emax=3) in case.in1c
attached here, and run dstart. This time, semicore error was appeared in
the *first* iteration!
I have been confused why?
(I see that default init_lapw for above setting (RKm=7 and Emax=3))
contains 24 kpoints. and error is appeared.
(In another scf with above struct (RKm=7 and Emax=1.5), default number
of kpoints was 7! ). and scf ended successfully.
May this error be dependent on the number of k-point and the number of
cores in .machines file?
.machines file in two calculations contained 7 cores.
On Tuesday, November 14, 2023 at 03:08:50 PM GMT+3:30, Peter Blaha
<peter.bl...@tuwien.ac.at> wrote:
Again, your message gets too big. You must delete the older content.
-------------------------------------------------------
grep :DIS in case.scf:
:DIS : CHARGE DISTANCE ( 0.0122755 for atom 3 spin 1)
0.0083335
:DIS : CHARGE DISTANCE ( 0.0117894 for atom 3 spin 1)
0.0077543
:DIS : CHARGE DISTANCE ( 0.0405700 for atom 1 spin 1)
0.0200036
:DIS : CHARGE DISTANCE ( 0.2010006 for atom 1 spin 1)
0.0741310
:DIS : CHARGE DISTANCE ( 0.0164221 for atom 1 spin 1)
0.0107305
:DIS : CHARGE DISTANCE ( 0.1052329 for atom 1 spin 1)
0.0370176
:DIS : CHARGE DISTANCE ( 0.0075476 for atom 1 spin 1)
0.0021153
:DIS : CHARGE DISTANCE ( 0.0848258 for atom 1 spin 1)
0.0300654
:DIS : CHARGE DISTANCE ( 0.0018758 for atom 1 spin 1)
0.0007564
:DIS : CHARGE DISTANCE ( 0.0008796 for atom 3 spin 1)
0.0006306
:DIS : CHARGE DISTANCE ( 0.0013281 for atom 3 spin 1)
0.0008331
after iteration 3, the semicore error is appeared for rkm=7.
------------------------------------------------------------------
Nobody knows what you were sending. Is this from the RKM=6 calculation ?
You have done a scf with RKmax=6. This should be saved.
Then you have an empty scffile, and you should increase RKmax to 7 and
run the scf.
If the error occurs after iteration 3, we expect to see exactly 3 lines.
???
Did you ever save the rkm6 results ?
restore them in a new directory, increase rkmax and then run_lapw.
When it crashes, show us :dis.
--
--------------------------------------------------------------------------
Peter BLAHA, Inst.f. Materials Chemistry, TU Vienna, A-1060 Vienna
Phone: +43-1-58801-165300
Email: peter.bl...@tuwien.ac.at <mailto:peter.bl...@tuwien.ac.at>
WIEN2k: http://www.wien2k.at <http://www.wien2k.at>
WWW: http://www.imc.tuwien.ac.at <http://www.imc.tuwien.ac.at>
-------------------------------------------------------------------------
_______________________________________________
Wien mailing list
Wien@zeus.theochem.tuwien.ac.at <mailto:Wien@zeus.theochem.tuwien.ac.at>
http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien
<http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien>
SEARCH the MAILING-LIST at:
http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
<http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html>
_______________________________________________
Wien mailing list
Wien@zeus.theochem.tuwien.ac.at
http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien
SEARCH the MAILING-LIST at:
http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
--
--------------------------------------------------------------------------
Peter BLAHA, Inst.f. Materials Chemistry, TU Vienna, A-1060 Vienna
Phone: +43-1-58801-165300
Email: peter.bl...@tuwien.ac.at WIEN2k: http://www.wien2k.at
WWW: http://www.imc.tuwien.ac.at
-------------------------------------------------------------------------
_______________________________________________
Wien mailing list
Wien@zeus.theochem.tuwien.ac.at
http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien
SEARCH the MAILING-LIST at:
http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
