Dear Prof. Blaha,
thanks for the fast reply. I will try that later. Currently calculations
are running. I wanted to calculate the 19F chemical shift for TlF3 just
as a model compound for experimental / computational shift correlations.
So it is not that important for my work.
Thanks again and 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.fechtelk...@ruhr-uni-bochum.de
Web Page:
https://www.ruhr-uni-bochum.de/kristallographie/kc/mitarbeiter/fechtelkord/
_______________________________________________
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