Quantum ESPRESSO入力ファイル作成手順７．局所状態密度の計算

クロスアビリティ Winmostarサポートチームです。
WinmostarはQuautum ESPRESSOを簡単かつ高度に利用するための統合GUI環境となっています。学生は無料で利用でき、学生以外も無料トライアルを入手することができますので、是非 [WinmostarのWebサイト](https://winmostar.com/jp/) をご覧ください。

#0. はじめに
第一原理計算ソフトQuantunm ESPRESSO <https://www.quantum-espresso.org> の入力ファイルの書き方について、複数回にわたって解説していきます。今回は第７回目です。局所状態密度の計算方法について解説したいと思います。

> この記事では、初歩的な内容しか紹介しません。
> 入力ファイルの詳細について知りたい方は、マニュアル http://www.quantum-espresso.org/Doc/INPUT_PW.html をご参照下さい。


#1. 計算手順
局所状態密度の計算手順は、基本的には全状態密度の計算手順 <https://qiita.com/xa_member/items/7aa75304b2ce7f90c8d4> と同じです。「SCF計算」および「Non-SCF計算」を完了した後、「状態密度データの生成」を行います。全の状態密度の計算では、状態密度データの生成に実行モジュール"dos.x"を用います。一方、局所状態密度の場合には、dos.xの代わりに"projwfc.x"を使用します。


#2. Non-SCF計算における注意事項

##2-1. Tetrahedron法の使用不可
全状態密度の計算では、Tetrahedron法が利用可能でした。Non-SCF計算の入力ファイルにおける occupations = "tetrahedra" という設定が、これに該当します。残念ながら、局所状態密度の計算においては、Tetrahedron法は利用できません。Non-SCF計算の入力ファイルでは、

```
&SYSTEM
    occupations = "smearing"
/
```

と設定する必要がるので、留意してください。


##2-2. 波動関数データの集約

また、全状態密度データの生成(dos.x)には、バンドエネルギーのみが入力情報となります。一方、局所状態密度データの生成(projwfc.x)では、バンドエネルギーに加えて波動関数の情報も必要となります。Quantum ESPRESSO(Ver.6.2以前)のデフォルト設定では、波動関数の情報は MPI並列のプロセス毎に入出力されます。したがって、デフォルト設定では、Non-SCF計算におけるMPI並列数をprojwfc.xにおける並列数に合致させる必要があります。これでは少々不便ですよね。そこで、全プロセスに分散された波動関数データを一か所に集約して入出力することで、この問題を回避します。具体的には、Non-
SCF計算における入力ファイルにて、

```
&CONTROL
    wf_collect = .TRUE.
/
```

と設定することで、波動関数データを集約できます。Quantum ESPRESSOのVer.6.3以降では、デフォルトで wf_collect = .TRUE. なので当該の設定は不要です。


#3. 状態密度データの生成
Non-SCF計算が完了したら、実行モジュール"projwfc.x"を用いて局所状態密度のデータを生成します。projwfc.xの実行には、以下のような入力ファイルを用います。

```
&PROJWFC
    emax   =  50.0
    emin   = -50.0
    deltae =  0.01
/
```

dos.x(全状態密度のデータ生成)の場合と同様の設定内容になります。emaxおよびeminは、それぞれ、状態密度データのエネルギー範囲の上限および下限値です。deltaeは状態密度データのエネルギー刻み幅です。

projwfc.xを実行すると、各原子・各軌道に対して局所状態密度のデータが生成されます。当該データは、原子および軌道の組み合わせ毎に個別のファイルに出力されます。ファイル名は、pwscf.pdos_atm#<i>原子番号</i>(<i>原子名</i>)_wfc#<i>軌道番号</i>(<i>軌道名</i>) です。例えば、１番目のSi原子のp軌道の場合、ファイル名は pwscf.pdos_atm#1(Si)_wfc2(p) となります。ここで、原子番号は、SCFおよびNon-SCF計算における入力ファイルのATOMIC_POSITIONSで指定した順番を意味します。当該ファイルには、以下のようなデータが記されます。

```
# E (eV)   ldos(E)   pdos(E)    pdos(E)    pdos(E)   
 -6.331  0.348E-10  0.293E-10  0.273E-11  0.273E-11
 -6.321  0.580E-10  0.474E-10  0.529E-11  0.529E-11
 -6.311  0.967E-10  0.764E-10  0.102E-10  0.102E-10
 -6.301  0.161E-09  0.123E-09  0.193E-10  0.193E-10
 -6.291  0.270E-09  0.197E-09  0.363E-10  0.363E-10
 -6.281  0.451E-09  0.316E-09  0.676E-10  0.676E-10
```

一列目が、エネルギーです。三列目以降が、軌道毎の状態密度です。p軌道の場合には、pz, px, py の３つデータがそれぞれ、三列目、四列目、五列目に出力されます。二列目は、全軌道の状態密度の総和です。この例では、三列目と四列目と五列目の値の和が、二列目に出力されます。p, d, f軌道について、球面調和関数の出力順番は下表の通りです。

| 軌道 | 出力順番 |
|-----|---------|
|  p  | z, x, y |
|  d  | z2, xz, yz, x2-y2, xy |
|  f  | z3, xz2, yz2, zx2-zy2, xyz, x3-3xy2, 3yx2-y3 |

また、スピン分極DFTの場合には、同ファイル内に上向きスピンと下向きスピンのデータが別々に出力されます。

---------------------------------------
次回は、結晶のバンド構造計算について解説したいと思います。
なお、Winmostarを利用するとQuautum ESPRESSOをGUI上から簡単かつ高度に活用することができます。学生は無料で利用でき、学生以外も無料トライアルを入手することができますので、是非 [WinmostarのWebサイト](https://winmostar.com/jp/) をご覧ください。

Quantum ESPRESSO入力ファイル作成手順シリーズ
[１．結晶構造の作成](https://qiita.com/xa_member/items/727c1a62930611babaf7) 
[２．SCF計算の設定](https://qiita.com/xa_member/items/a098e94d633a1d29f1c0) 
[３．擬ポテンシャルファイルの選択方法](https://qiita.com/xa_member/items/006d9dd44c662a17903d) 
[４．原子座標のみの最適化](https://qiita.com/xa_member/items/dc860b1656f7dfe93c8a)
[５．格子ベクトル及び格子内部の原子座標の最適化](https://qiita.com/xa_member/items/4c667817e50facd8a770)
[６．状態密度の計算](https://qiita.com/xa_member/items/7aa75304b2ce7f90c8d4)
[７．局所状態密度の計算](https://qiita.com/xa_member/items/95c7bd64981f7329b94b)
[８．バンド構造の計算](https://qiita.com/xa_member/items/a519c53b63df75a68894)
[９．バンド数の設定](https://qiita.com/xa_member/items/d923ea537328cdba8b10)
[１０．van der Waals相互作用](https://qiita.com/xa_member/items/bd64362329c7235bf248)
[１１．LDA+U法](https://qiita.com/xa_member/items/ab3e1333e92bc424f5a7)
[１２．NEB法](https://qiita.com/xa_member/items/e3b63a1fdca12f95e8cb)
[１３．Phonon計算（特定q点）](https://qiita.com/xa_member/items/77e464c01c5db1e847ea)
[１４．Phonon計算（バンド構造）](https://qiita.com/xa_member/items/5c02fcae08431588659d)
[１５．擬ポテンシャルの作成](https://qiita.com/xa_member/items/5feb8f432170141c89cd)
[１６．擬ポテンシャルのテスト](https://qiita.com/xa_member/items/3902c1660fb7064a4c34)
[１７．SCF計算の入力ファイルの実例](https://qiita.com/xa_member/items/febde40c8a17b7b83cc5)