Hubbard-U-Rechnungen
Bei Verwendung eines U-Parameters für einige Atome bei Start der Rechnung LDA+U auswählen.
Dann werden die U-Parameter für die einzelnen Atome in die Datei *.inorb
eingegeben:
1 2 0 (nmod, natorb, ipr)
PRATT 1.0 (BROYD/PRATT, mixing)
1 1 2 (iatom nlorb, lorb)
2 1 2 (iatom nlorb, lorb)
1 (nsic 0..AFM, 1..SIC, 2..HFM)
0.52 0.00 (U J (Ry) Note: we recommend to use U_eff = U-J and J=0)
0.52 0.00 (U J)
Wie bei der DOS wird also für jedes Atom das Orbital und der U-Parameter angegeben.
- nmod: 1
- natorg: Zahl der Atome, denen man ein U verpassen will
- ipr: 0
- iatom, nlorg, lorb: Atomnummer, Zahl der Orbitale mit U, Nebenquantenzahl (2=d, 3=f)
der mit U behafteten Orbitale
- nsic, AFM, HFM: 1=AFM
- U, J: für jedes Atom eine Zeile mit U und J (Werte in [Ry])
Zusätzlich benötigt man eine Datei *.indm. Diese sollte wie folgt aussehen:
-9. (Emin cutoff energy)
1 (number of atoms for which density matrix is calculated)
1 1 3 (index of 1st atom, number of L's, L1)
0 0 (r-index, (l,s)index)
1 2 0 (nmod, natorb, ipr)
Dabei bedeutet: