Differences

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

--- quantumespresso:部分状態密度 [2025/02/09 22:50] – koudai
+++ quantumespresso:部分状態密度 [2025/02/26 03:04] (current) – [注意] koudai
@@ Line 22: / Line 22: @@
   nspin = 2,  starting_magnetization(1)=0.7,
   ecutwfc = 24.0, ecutrho = 288.0,
-  occupations = 'tetrahedra_opt'
+  occupations = 'smearing'
+  smearing = 'mp'
+  degauss = 0.02
 /
@@ Line 47: / Line 49: @@
 k点を細かく取って、それぞれの点でのエネルギーを求めます。
 Ni.scf.inをコピーして作成します。
-calculations='nscf'にして、k点の数を多めに取ります。
+変更箇所は以下の2ヶ所です。
+  * <nowiki>calculations='nscf'</nowiki>にする
+  * k点の数を多めに取る
 <file - Ni.nscf.in>
@@ Line 62: / Line 66: @@
   starting_magnetization(1) = 0.7,
   ecutwfc = 24.0, ecutrho = 288.0, nbnd = 8,
-  occupations = 'tetrahedra_opt'
+  occupations = 'smearing'
+  smearing = 'mp'
+  degauss = 0.02
 /
@@ Line 110: / Line 116: @@
 最後に、状態密度を軌道ごとに分解します。
+Ni.dos.inをコピーしてNi.pdos.inというファイル名にしてください。
+変更箇所はありません。
 <file - Ni.pdos.in>
@@ Line 190: / Line 198: @@
 すると 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'とすると、状態密度は自動的にテトラヘドロン法を使ったものになります。