[Wien] [WIEN2k] forrtl IO error in x_nmr_lapw for Heavy metal structures (TlF3, HgF2)

[Michael Fechtelkord]

Dear Prof. Blaha,


I recompiled lapw1 with LOMAX = 4 but then the scf cycle fails. The real 
problem was the used cut off energy of -11 Ry in init_lapw. That 
introduces also too much orbitals in x_nmr -mode in1 and a qtl-b errors 
in the first two loops. After using the default value of -6 Ry for core 
/ valence separation, scf cycles converge without QTL-B errors and the 
x_nmr initializations starts with less atomic orbitals and  the 
calculations (x_nmr -p) no longer produces forrtl I/O errors.


Thanks again for the help!


Best regards,

Michael Fechtelkord


Am 12.11.2023 um 23:28 schrieb Peter Blaha:
> Once I've seen your in1 file, the solution is probably very simple:
>
> I did not know that you included the 4f states of Tl (near -8 Ry) as 
> valence.
> The nmr code constructs by default NMR-local orbitals up to 
> "l-exception" + 1, i.e. up to l=4 when you have l=3 states listed in 
> the regular case.in1
>
> While this is possible, it requires to recompile lapw1/2,nmr with a 
> modified parameter  LOMAX =4 (param.inc in lapw1, in other codes in 
> modules - do a search).
>
> This is necessary if you handle 4f elements or early 5d metals, 
> however, I very much doubt that it is a good idea to include the 4f 
> states for Tl (with RMT=2.5) as valence.  I would not use -ecut -11.
> All it produces is noise as the 4f convergence can be quite 
> problematic and SO effects might be of importance.
>
> Best regards
> Peter Blaha
>
> Am 12.11.2023 um 22:12 schrieb Michael Fechtelkord:
>> Lieber Herr Blaha,
>>
>> schon mal vorab herzlichen Dank für die schnelle Hilfe auch am 
>> Wochenende. Anbei die gewünschten Daten und folgendermaßen bin ich 
>> vorgegangen:
>>
>> im Verzeichnis TlF3
>>
>> 1) cif2struct TlF3.cif
>>
>> 2) Kontrolle und Nachbearbeitung mit struct generator in w2web
>>
>> 3) rmt gesetzt mit 0% Reduktion in w2web struct Generator (set 
>> automatically RMT and continue editing)
>>
>> 4) Structfile abgeschlossen (save file and cleanup)
>>
>> Weiter im Terminalfenster:
>>
>> 5) init_lapw -b -rkmax 7 -numk 1000 -ecut -11 (endete mit ok)
>>
>> 6) run_lapw -p -ec 0.0001 -cc 0.0001 (konvergierte nach ca. 13 Zyklen)
>>
>> 7) save_lapw TlF3_pbe_rkmax_7_numk_1000_ecut_11_cc_0001
>>
>> 8) x kgen auf 10000 k points (habe es auch mit weniger probiert, 
>> daran liegt es wohl nicht)
>>
>> 9)  x_nmr_lapw -mode in1
>>
>> 10)  x_nmr_lapw -p
>>
>> Ich hänge auch das cif File und das machines File der Vollständigkeit 
>> halber an. NTMATMAX ist 40000, NUME 6000, OMP_NUM_THREADS 2
>>
>>
>> Sollten Sie zusätzliche Daten benötigen, schreiben Sie mich kurz an.
>>
>>
>> Viele Dank schon mal und einen guten Wochenstart
>>
>> wünscht
>>
>> Michael Fechtelkord
>>
>>
>
-- 
Dr. Michael Fechtelkord

Institut für Geologie, Mineralogie und Geophysik
Ruhr-Universität Bochum
Universitätsstr. 150
D-44780 Bochum

Phone: +49 (234) 32-24380
Fax:  +49 (234) 32-04380
Email: Michael.Fechtelkord at ruhr-uni-bochum.de
Web Page: https://www.ruhr-uni-bochum.de/kristallographie/kc/mitarbeiter/fechtelkord/