Differences

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

--- quantumespresso:部分状態密度 [2020/03/21 11:49] – [状態密度] koudai
+++ quantumespresso:部分状態密度 [2025/02/26 03:04] (current) – [注意] koudai
@@ Line 1: / Line 1: @@
 ====== 概要 ======
-強磁性体のNiをもとに、状態密度を軌道の重みで分けた部分状態密度を計算します
+強磁性体のNiをもとに、状態密度を軌道の重みで分けた部分状態密度を計算します。
 ====== 計算手順 ======
-擬ポテンシャルとして、ウルトラソフト型の Ni.pz-nd-rrkjus.UPF をダウンロードしました
+擬ポテンシャルとして、ウルトラソフト型の Ni.pz-nd-rrkjus.UPF をダウンロードしました。
 ===== SCF計算 =====
@@ Line 12: / Line 12: @@
 <file - Ni.scf.in>
 &control
-    calculation='scf'
+  calculation='scf'
-    restart_mode='from_scratch',
+  prefix='Ni'
-    prefix='Ni',
+  pseudo_dir = './pseudo/'
-    pseudo_dir = './pseudo/',
+  outdir='./tmp/'
-    outdir='./tmp/'
 /
 &system
-    ibrav=2, celldm(1) =6.48, nat=1, ntyp=1,
+  ibrav=2, celldm(1) =6.48, nat=1, ntyp=1,
-    nspin = 2,  starting_magnetization(1)=0.7,
+  nspin = 2,  starting_magnetization(1)=0.7,
-    ecutwfc = 24.0, ecutrho = 288.0,
+  ecutwfc = 24.0, ecutrho = 288.0,
-    occupations='smearing', smearing='mv', degauss=0.02
+  occupations = 'smearing'
+  smearing = 'mp'
+  degauss = 0.02
 /
 &electrons
-    conv_thr = 1.0e-10
+  conv_thr = 1.0e-10
-    mixing_beta = 0.7
+  mixing_beta = 0.7
 /
@@ Line 46: / Line 47: @@
 ===== 状態密度 =====
-k点を細かく取って、それぞれの点でのエネルギーを求める。
+k点を細かく取って、それぞれの点でのエネルギーを求めます。
-Ni.scf.inをコピーして作成する。
+Ni.scf.inをコピーして作成します。
+変更箇所は以下の2ヶ所です。
+  * <nowiki>calculations='nscf'</nowiki>にする
+  * k点の数を多めに取る
-<file - Ni.dos.in>
+<file - Ni.nscf.in>
 &control
-  calculation = ’nscf’
+  calculation = 'nscf'
-  prefix = ’ni’,
+  prefix = 'Ni',
-  pseudo_dir = ’./pseudo/’,
+  pseudo_dir = './pseudo/'
-  outdir = ’./tmp/’
+  outdir = './tmp/'
 /
@@ Line 62: / Line 66: @@
   starting_magnetization(1) = 0.7,
   ecutwfc = 24.0, ecutrho = 288.0, nbnd = 8,
-  occupations = ’tetrahedra’
+  occupations = 'smearing'
+  smearing = 'mp'
+  degauss = 0.02
 /
@@ Line 81: / Line 87: @@
 <code>
-$ pw.x < ni.dos.in > ni.dos.out
+$ pw.x < ni.nscf.in > ni.nscf.out
 </code>
-次に、状態密度を計算する。
+次に、状態密度を計算します。
-<file - Ni.dos2.in>
+<file - Ni.dos.in>
 &dos
-  outdir = ’./tmp/’
+  outdir = './tmp/'
-  prefix = ’Ni’
+  prefix = 'Ni'
-  fildos = ’Ni.dos’,
+  fildos = 'Ni.dos'
   Emin = 5.0, Emax = 25.0, DeltaE = 0.1
   ngauss = 1, degauss = 0.02
@@ Line 104: / Line 110: @@
 <code>
-$ dos.x < Ni.dos2.in > Ni.dos2.out
+$ dos.x < Ni.dos.in > Ni.dos.out
 </code>
-ni.dos というファイルに、アップスピンとダウンスピン、それを合計したものの状態密度が出力される。
+Ni.dos というファイルに、アップスピンとダウンスピン、それを合計したものの状態密度が出力されます。
-最後に、状態密度を軌道ごとに分解する
+最後に、状態密度を軌道ごとに分解します。
+Ni.dos.inをコピーしてNi.pdos.inというファイル名にしてください。
+変更箇所はありません。
 <file - Ni.pdos.in>
 &projwfc
-  outdir = ’./tmp/’
+  outdir = './tmp/'
   prefix = 'Ni'
   Emin = 5.0, Emax = 25.0, DeltaE = 0.1
@@ Line 124: / Line 132: @@
 </code>
-Ni.pdos_atm#1(Ni)_wfc#1(s) にs軌道の部分状態密度が、Ni.pdos_atm#1(Ni)_wfc#2(d) にd電子の部分状態密度が、Ni.pdos_totに全状態密度が出力される。
+Ni.pdos_atm#1(Ni)_wfc#1(s) にs軌道の部分状態密度が、Ni.pdos_atm#1(Ni)_wfc#2(d) にd電子の部分状態密度が、Ni.pdos_totに全状態密度が出力されます。
-出力されたファイルでのd軌道の部分状態密度であるが、左から順にd3z2−r2, dzx, dzy, dx2−y2, dxyである。
+出力されたファイルでのd軌道の部分状態密度は、左から順にd3z2−r2, dzx, dzy, dx2−y2, dxyです。
-詳しくはマニュアル /PP/Doc/INPUT PROJWFC.html を参照すること。
+詳しくは[[https://www.quantum-espresso.org/Doc/INPUT_PROJWFC.html|マニュアル]]を参照してください。
@@ Line 134: / Line 142: @@
 &control にある calculation を scf から bands に変更する。
 K_POINTSをバンドの経路に変更する
-<file - Ni.scf.in>
+<file - Ni.bands.in>
 &control
-    calculation='scf'
+    calculation='bands'
-    restart_mode='from_scratch',
+    prefix='Ni'
-    prefix='Ni',
+    pseudo_dir = './pseudo/'
-    pseudo_dir = './pseudo/',
     outdir='./tmp/'
 /
@@ Line 174: / Line 181: @@
 projwfc.xを使って部分状態密度を計算する
-<file - >
+<file - Ni.kpdos.in>
 &projwfc
-  outdir = ’./tmp/'
+  outdir = './tmp/'
-  prefix = ’Ni’
+  prefix = 'Ni'
   ngauss = 0, degauss = 0.036748
   DeltaE = 0.01
   kresolveddos = .true.
-  filpdos = ’Ni.k’
+  filpdos = 'Ni.k'
 /
 </file>
@@ Line 189: / Line 196: @@
 </code>
-すると ni.kpdos.out というファイルの他に、ni.k.pdos_tot と ni.k.pdos_atm#1(Ni)_wfc#1(s) と ni.k.pdos atm#1(Ni)_wfc#2(d) というファイルができるはずである。
+すると Ni.kpdos.out というファイルの他に、Ni.k.pdos_tot, Ni.k.pdos_atm#1(Ni)_wfc#1(s), Ni.k.pdos atm#1(Ni)_wfc#2(d) というファイルができるはずである。
 これはそれぞれ両方、 s電子のみ、d電子のみのk分解状態密度のデータである。
+===== 注意 =====
+nscf計算でoccupations='tetrahedra'とすると、状態密度は自動的にテトラヘドロン法を使ったものになります。