This shows you the differences between two versions of the page.
Both sides previous revision Previous revision Next revision | Previous revision | ||
wien2k:状態密度 [2020/06/22 18:42] koudai [手順] |
wien2k:状態密度 [2023/04/24 15:20] (current) koudai |
||
---|---|---|---|
Line 19: | Line 19: | ||
$ x lapw1 -p | $ x lapw1 -p | ||
</ | </ | ||
- | * スピン分極している場合は、アップスピンとダウンスピンの両方で計算します。以下、xによる計算が出てきた場合は同様です。出力ファイル名も最後にupあるいはdnが追加されます< | + | * スピン分極している場合は、アップスピンとダウンスピンの両方を計算する必要があります。以下、xによる計算が出てきた場合は同様です。出力ファイル名も最後にupあるいはdnが追加されます< |
$ x lapw1 -up -p | $ x lapw1 -up -p | ||
$ x lapw1 -dn -p | $ x lapw1 -dn -p | ||
</ | </ | ||
- (必要ならば)case.inqを編集します。局所座標を回すことができます(気が向いたらその書き方も説明しようと思います) | - (必要ならば)case.inqを編集します。局所座標を回すことができます(気が向いたらその書き方も説明しようと思います) | ||
- | - x qtlにより軌道ごとの電荷密度を計算します< | + | - QTLというプログラムを使って軌道ごとの電荷密度を計算します< |
- | $ x qtl | + | $ x qtl -p |
</ | </ | ||
- | * 今のところ qtl は並列計算には対応していないので -p をつけないでください。 | + | * 以下のようにLAPW2でも軌道ごとの電荷密度が計算できますが、case.inqは読み込みません< |
- | - case.qtl のヘッダーを見て、使用できる軌道を確認します。あるいは次のようにします< | + | - case.qtl のヘッダーを見て、使用できる軌道を確認します。あるいは次のようにして確認します< |
- | $ grep "tot" TiC.qtl | + | $ grep "JATOM" TiC.qtl |
- | | + | |
- | | + | |
</ | </ | ||
- DOSの計算の入力ファイル case.int を生成します。TiCの例です< | - DOSの計算の入力ファイル case.int を生成します。TiCの例です< | ||
Line 39: | Line 39: | ||
* 最後の end を忘れがちなので気をつけてください。 | * 最後の end を忘れがちなので気をつけてください。 | ||
* \$ configure_int_lapw だけで実行すれば対話形式で作成できます。 | * \$ configure_int_lapw だけで実行すれば対話形式で作成できます。 | ||
- | - 必要ならば case.int の2行目を編集します< | + | - 必要ならば case.int の2行目を編集します< |
TiC # | TiC # | ||
0.25767 0.002 1.2492922854 0.003 # | 0.25767 0.002 1.2492922854 0.003 # | ||
Line 59: | Line 59: | ||
</ | </ | ||
- 結果は case.dos1 と case.dos1ev に出力される(前者がエネルギー単位Ryで、後者がeV)ので、GNUPLOTなどでプロットします。 | - 結果は case.dos1 と case.dos1ev に出力される(前者がエネルギー単位Ryで、後者がeV)ので、GNUPLOTなどでプロットします。 | ||
+ | - 部分状態密度の軌道の種類が多いとcase.dos2, | ||
図はTiCの状態密度の計算例です。 | 図はTiCの状態密度の計算例です。 | ||
{{ : | {{ : |