第一原理分子動力学プログラム STATE Senri Wiki
開始行:
* はじめに [#r9af3044]
従来のSTATEではGGA等で収束させた電荷密度と原子配置を用い...
* 使い方 [#rfb133e2]
** 準備 [#s4e48aee]
プログラムからRomán-Pérez-Solerアルゴリズム用のvdW kernel...
~/STATE/src/
等にvdwdphi.datをおいて、計算を行うディレクトリにvdwdphi....
ln -fs ~/STATE/src/vdwdphi.dat
ln -fs ~/STATE/src/kernel.dat
ない場合は、smithの/home/hamamoto/STATE/srcからコピーし...
** Post-processing vdW-DF（非推奨）[#mb073ffa]
GGAで電子密度と原子構造を求めたあとicond=26に変更して計算...
従来通りnon-self-consistentなvdW-DFの計算が行えます。
この場合、実空間で非局所相関エネルギーを計算するので計算...
入力ファイルの下に
&VDW-DF
FFTVDW
&END
を追加することで、Román-Pérez-Solerアルゴリズムを用いるこ...
vdW-DF1以外のvdW密度汎関数を用いる場合には交換汎関数と非...
例えばvdW-DF1に対して交換汎関数を改良したoptB86b-vdWを用...
&VDW-DF
FFTVDW
x optb86b
&END
のように変更します。xの次の文字列は交換汎関数を表します。
さらに交換汎関数に加え非局所相関汎関数も変更したrev-vdW-D...
&VDW-DF
FFTVDW
VDW-DF2
x b86r
&END
とします。
VDW-DF2を省略するとVDW-DF1の非局所相関汎関数が指定されま...
交換汎関数と非局所相関汎関数の組み合わせは下表を参考にし...
*** 注意 [#xcfedb7e]
STATE 5.5.4以降では交換汎関数と非局所相関汎関数の組み合わ...
下表のvdW密度汎関数名をxctypeに指定すればPost-processing ...
この方法では&color(red){非局所相関エネルギーが常にvdW-DF1...
vdW-DF2およびrev-vdW-DF2を用いる際には注意して下さい。
最新版のSTATE 5.6.0でもこのバグは修正されていないので、と...
** Self-consistent vdW-DF（推奨）[#v5a8c6f4]
一方、icond=0(初めから) or 1(続きから)のとき、xctypeとし...
Self-consistentの場合は&VDW-DFのブロックでFFTVDWを指定す...
表中で非局所相関汎関数のvdW-DF1とvdW-DF2の違いは、パラメ...
|xctype|非局所相関|交換|文献|h
|vdw-df または drsll|vdW-DF1|revPBE|[[M. Dion et al., Phy...
|vdw-df2 または lmkll|vdW-DF2|PW86R|[[K. Lee et al., Phys...
|c09 または c09-vdw または drsllc|vdW-DF1|C09x|[[V. R. Co...
|c09-vdw2 または lmkllc|vdW-DF2|C09x||
|optpbe または optpbe-vdw|vdW-DF1|optPBE|[[J. Klimeš et a...
|optb88 または optb88-vdw または kbm|vdW-DF1|optB88|[[J. ...
|optb86b または optb86b-vdw|vdW-DF1|optB86b|[[J. Klimeš e...
|rev-vdw-df2 または lmkllh|vdW-DF2|B86R|[[I. Hamada, Phys...
//Román-Pérez-Solerアルゴリズムを適用するには、入力ファイ...
// &VDW-DF
// FFTVDW
// &END
*** 注意 [#da6c43b0]
STATE 5.5.4および5.6.0ではB86交換ポテンシャルの計算のバグ...
optB86b-vdWおよびrev-vdW-DF2を用いる際には修正版を入手し...
** vdW kernelの分割 [#l5f17611]
Román-Pérez-SolerアルゴリズムではvdW kernelをq0に関して分...
例えばbulkの軽元素の場合、q0の最大値は3程度なので
&VDW-DF
qcut 5.d0
nqs 10
&END
とします。一方、bulkの遷移金属元素ではq0の最大値が7程度に...
&VDW-DF
qcut 10.d0
nqs 20
&END
とします。軽元素であっても孤立系や表面系ではq0の分布が長...
分割数nqsは大きければ大きいほど良いというわけでもなく、ke...
特に大きな系では計算機のメモリに応じてnqsの大きさを制限す...
経験的には原子数50個程度の系でnqsを20より大きくしても
全エネルギーの変化は数meV程度に収まるようです。
何も指定しない場合はqcut=3.d0、nqs=10が設定されますが、次...
終了行:
* はじめに [#r9af3044]
従来のSTATEではGGA等で収束させた電荷密度と原子配置を用い...
* 使い方 [#rfb133e2]
** 準備 [#s4e48aee]
プログラムからRomán-Pérez-Solerアルゴリズム用のvdW kernel...
~/STATE/src/
等にvdwdphi.datをおいて、計算を行うディレクトリにvdwdphi....
ln -fs ~/STATE/src/vdwdphi.dat
ln -fs ~/STATE/src/kernel.dat
ない場合は、smithの/home/hamamoto/STATE/srcからコピーし...
** Post-processing vdW-DF（非推奨）[#mb073ffa]
GGAで電子密度と原子構造を求めたあとicond=26に変更して計算...
従来通りnon-self-consistentなvdW-DFの計算が行えます。
この場合、実空間で非局所相関エネルギーを計算するので計算...
入力ファイルの下に
&VDW-DF
FFTVDW
&END
を追加することで、Román-Pérez-Solerアルゴリズムを用いるこ...
vdW-DF1以外のvdW密度汎関数を用いる場合には交換汎関数と非...
例えばvdW-DF1に対して交換汎関数を改良したoptB86b-vdWを用...
&VDW-DF
FFTVDW
x optb86b
&END
のように変更します。xの次の文字列は交換汎関数を表します。
さらに交換汎関数に加え非局所相関汎関数も変更したrev-vdW-D...
&VDW-DF
FFTVDW
VDW-DF2
x b86r
&END
とします。
VDW-DF2を省略するとVDW-DF1の非局所相関汎関数が指定されま...
交換汎関数と非局所相関汎関数の組み合わせは下表を参考にし...
*** 注意 [#xcfedb7e]
STATE 5.5.4以降では交換汎関数と非局所相関汎関数の組み合わ...
下表のvdW密度汎関数名をxctypeに指定すればPost-processing ...
この方法では&color(red){非局所相関エネルギーが常にvdW-DF1...
vdW-DF2およびrev-vdW-DF2を用いる際には注意して下さい。
最新版のSTATE 5.6.0でもこのバグは修正されていないので、と...
** Self-consistent vdW-DF（推奨）[#v5a8c6f4]
一方、icond=0(初めから) or 1(続きから)のとき、xctypeとし...
Self-consistentの場合は&VDW-DFのブロックでFFTVDWを指定す...
表中で非局所相関汎関数のvdW-DF1とvdW-DF2の違いは、パラメ...
|xctype|非局所相関|交換|文献|h
|vdw-df または drsll|vdW-DF1|revPBE|[[M. Dion et al., Phy...
|vdw-df2 または lmkll|vdW-DF2|PW86R|[[K. Lee et al., Phys...
|c09 または c09-vdw または drsllc|vdW-DF1|C09x|[[V. R. Co...
|c09-vdw2 または lmkllc|vdW-DF2|C09x||
|optpbe または optpbe-vdw|vdW-DF1|optPBE|[[J. Klimeš et a...
|optb88 または optb88-vdw または kbm|vdW-DF1|optB88|[[J. ...
|optb86b または optb86b-vdw|vdW-DF1|optB86b|[[J. Klimeš e...
|rev-vdw-df2 または lmkllh|vdW-DF2|B86R|[[I. Hamada, Phys...
//Román-Pérez-Solerアルゴリズムを適用するには、入力ファイ...
// &VDW-DF
// FFTVDW
// &END
*** 注意 [#da6c43b0]
STATE 5.5.4および5.6.0ではB86交換ポテンシャルの計算のバグ...
optB86b-vdWおよびrev-vdW-DF2を用いる際には修正版を入手し...
** vdW kernelの分割 [#l5f17611]
Román-Pérez-SolerアルゴリズムではvdW kernelをq0に関して分...
例えばbulkの軽元素の場合、q0の最大値は3程度なので
&VDW-DF
qcut 5.d0
nqs 10
&END
とします。一方、bulkの遷移金属元素ではq0の最大値が7程度に...
&VDW-DF
qcut 10.d0
nqs 20
&END
とします。軽元素であっても孤立系や表面系ではq0の分布が長...
分割数nqsは大きければ大きいほど良いというわけでもなく、ke...
特に大きな系では計算機のメモリに応じてnqsの大きさを制限す...
経験的には原子数50個程度の系でnqsを20より大きくしても
全エネルギーの変化は数meV程度に収まるようです。
何も指定しない場合はqcut=3.d0、nqs=10が設定されますが、次...
ページ名:
![[PukiWiki]](image/logo_state-senri3_80x80.png "[PukiWiki]"){kind=logo}
