Differences

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

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