雑多な記録

User Tools

Site Tools

Trace: • dft_u • フォノンの状態密度と分散 • scf計算 • tb-mbj交換ポテンシャル
quantumespresso:dft_u

Differences

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

Link to this comparison view

Next revision 		Previous revision
quantumespresso:dft_u [2018/11/03 16:21]
koudai 作成 		quantumespresso:dft_u [2021/06/27 22:04] (current)
Line 11: Line 11:
  
  
-====== 計算手順 ====== +====== NiO ======
- +
-===== NiO =====+
  
 酸化ニッケル(NiO)は典型的な反強磁性体である。 ただし、電子相関が強いため、通常のDFTによるバンド計算ではバンドギャップの大きさを説明できない。 酸化ニッケル(NiO)は典型的な反強磁性体である。 ただし、電子相関が強いため、通常のDFTによるバンド計算ではバンドギャップの大きさを説明できない。
Line 19: Line 17:
  
   * W. L. Roth, Phys. Rev. 110, 1333 (1958).   * W. L. Roth, Phys. Rev. 110, 1333 (1958).
 +====== 計算手順 ======
 +
 +
  
  
Line 24: 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 29: Line 31:
 ここではそれぞれをNi1とNi2で分けて指定している。 ここではそれぞれをNi1とNi2で分けて指定している。
  
 +nspin=2のときは、絶縁体であってもoccupationsの指定が必要である。
  
 <file - NiO.scf.in> <file - NiO.scf.in>
 &control &control
    calculation='scf',    calculation='scf',
-   prefix='nio',+   prefix='NiO',
    pseudo_dir = './',    pseudo_dir = './',
    outdir='./'    outdir='./'
Line 49: 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 68: Line 72:
  4 4 4 1 1 1  4 4 4 1 1 1
 </file> </file>
 +tot_magnetizationは単位胞内の磁化の合計を表していて、反強磁性の場合はゼロである。 
 +もしこれを指定しないのであれば、半導体・絶縁体であってもoccupations, smearing, degaussの指定が必要である。
 実行は次の通り。 実行は次の通り。
  
Line 75: Line 80:
 出力ファイルを見れば全磁化がゼロで絶対磁化が有限の値を持つので、反強磁性秩序を起こしていることがわかる。 出力ファイルを見れば全磁化がゼロで絶対磁化が有限の値を持つので、反強磁性秩序を起こしていることがわかる。
  
-<file - nio.scf.out>+<file - NiO.scf.out>
 （略） （略）
      total magnetization       =     0.00 Bohr mag/cell      total magnetization       =     0.00 Bohr mag/cell
Line 104: 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 141: Line 146:
   $ dos.x < NiO.dos.in > NiO.dos.out   $ dos.x < NiO.dos.in > NiO.dos.out
  
-状態密度はnio.dosに出力される。+状態密度はNiO.dosに出力される。
  
  
Line 167: 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,
Line 218: Line 221:
 状態密度の計算はさきほどと同様である。 Uを入れることでバンドギャップが4eV程度となった。 状態密度の計算はさきほどと同様である。 Uを入れることでバンドギャップが4eV程度となった。
  
-dos.png +{{:quantumespresso:nio_u_dos.png?400|}}
quantumespresso/dft_u.1541229704.txt.gz · Last modified: 2021/06/27 22:00 (external edit)

Page Tools

Except where otherwise noted, content on this wiki is licensed under the following license: CC Attribution-Share Alike 4.0 International
CC Attribution-Share Alike 4.0 International Donate Powered by PHP Valid HTML5 Valid CSS Driven by DokuWiki