User Tools

Site Tools


quantumespresso:dft_u

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revision Previous revision
Next revision
Previous revision
quantumespresso:dft_u [2020/03/09 20:22]
koudai [DFT計算]
quantumespresso:dft_u [2021/06/27 22:04] (current)
Line 25: Line 25:
  
 まずはUを入れずに計算を行う。 まずはUを入れずに計算を行う。
 +
 擬ポテンシャルはウルトラソフト型のO.pbe-rrkjus.UPFとNi.pbe-nd-rrkjus.UPFを使う。 擬ポテンシャルはウルトラソフト型のO.pbe-rrkjus.UPFとNi.pbe-nd-rrkjus.UPFを使う。
 単位胞は反強磁性秩序を示すので2倍となる。 単位胞は反強磁性秩序を示すので2倍となる。
Line 30: Line 31:
 ここではそれぞれをNi1とNi2で分けて指定している。 ここではそれぞれをNi1とNi2で分けて指定している。
  
 +nspin=2のときは、絶縁体であってもoccupationsの指定が必要である。
  
 <file - NiO.scf.in> <file - NiO.scf.in>
Line 50: Line 52:
    starting_magnetization(2)= 0.5,     starting_magnetization(2)= 0.5, 
    starting_magnetization(3)=-0.5,    starting_magnetization(3)=-0.5,
 +   tot_magnetization = 0.0
 / /
 &electrons &electrons
Line 69: Line 72:
  4 4 4 1 1 1  4 4 4 1 1 1
 </file> </file>
- +tot_magnetizationは単位胞内の磁化の合を表していて、反強磁性の場合ゼロである。 
-スピン分極させる算ではoccupations, smearing, degaussの指定が必要である。+もしこれを指定しないのであれば、半導体・絶縁体であってもoccupations, smearing, degaussの指定が必要である。
 実行は次の通り。 実行は次の通り。
  
Line 106: Line 109:
    ntyp= 3,    ntyp= 3,
    ecutwfc = 30.0, ecutrho = 240.0,    ecutwfc = 30.0, ecutrho = 240.0,
-   occupations = 'tetrahedra', 
    nspin=2,     nspin=2, 
    starting_magnetization(2)= 0.5,     starting_magnetization(2)= 0.5, 
    starting_magnetization(3)=-0.5,    starting_magnetization(3)=-0.5,
 +   tot_magnetization = 0.0
 / /
 &electrons &electrons
Line 143: Line 146:
   $ dos.x < NiO.dos.in > NiO.dos.out   $ dos.x < NiO.dos.in > NiO.dos.out
  
-状態密度はnio.dosに出力される。+状態密度はNiO.dosに出力される。
  
  
Line 169: Line 172:
    ntyp= 3,    ntyp= 3,
    ecutwfc = 30.0, ecutrho = 240.0,    ecutwfc = 30.0, ecutrho = 240.0,
-   occupations='smearing',  
-   smearing='mp',  
-   degauss=0.02, 
    nspin=2,     nspin=2, 
    starting_magnetization(2)= 0.5,     starting_magnetization(2)= 0.5, 
    starting_magnetization(3)=-0.5,    starting_magnetization(3)=-0.5,
 +   tot_magnetization = 0.0,
    lda_plus_u = .true.,    lda_plus_u = .true.,
    Hubbard_U(2) = 8.0,    Hubbard_U(2) = 8.0,
quantumespresso/dft_u.1583752971.txt.gz · Last modified: 2021/06/27 22:00 (external edit)