第一原理分子動力学プログラム STATE Senri Wiki
開始行:
* STATE クイックリファレンス [#e91b1b9c]
** 入力ファイルの例 (CO分子) [#s044357d]
0 0 0 0 0 0 : dummy line (6 i...
5.50 20.00 2 2 2 : GMAX, GMAXP, NT...
1 0 : space group num...
6.00 4.00 4.00 90.00 90.00 90.00 : a, b, c, alpha,...
1 1 1 1 1 1 : knx, kny, knz, ...
1 0 : NCORD, NINV
0.0000 0.0000 0.0000 1 1 1 : cps, iwei, imdt...
2.2000 0.0000 0.0000 1 1 2 : cps, iwei, imdt...
6 0.1500 51577.50 3 1 0.d0 : IATOMN, ALFA, A...
8 0.1500 51577.50 3 1 0.d0 : IATOMN, ALFA, A...
0 0 0 0 0 : ICOND, INIPOS, ...
0 1 : IPRE, IPRI
200 200 0 57200.00 0 : NMD1, NMD2, ite...
3 1 : way_mix, mix_what
0 8 0.8 : starting mixing...
0.60 0.50 0.60 0.70 1.00 : DTIM1, DTIM2, D...
30.00 2 1 0.10D-08 1.d-06 : DTIO, IMDALG, I...
0.0010 0.10D+02 0 : WIDTH, FORCCR, ...
ggapbe 1 : XCTYPE, nspin
1.00 3 : destm, n_stm
102 : NBZTYP
0 0 0 : NKX, NKY, NKZ...
0 0 0 : NKX2, NKY2, NKZ...
8 : NEG (# of bands)
1 : NEXTST (1: G-sp...
0 : 0; random numbe...
2 : imsd (2: Davids...
0 : eval. eko diff....
0 : npdosao
0 0.0 : SM_N doping
*** ダミー行 [#b28085e8]
0 0 0 0 0 0 : dummy line (6 i...
6個の整数を指定。歴史的経緯で存在するダミー行。~
※過去のバージョンでは1カラム目の変数が1の場合は基本並進ベ...
※現状STATEでは無視されているが、ユーティリティプログラムr...
*** カットオフと原子数 [#lad659e0]
5.50 20.00 2 2 2 : GMAX, GMAXP, NT...
-GMAX: 波動関数を展開する平面波のカットオフエネルギーを指...
-GMAXP: 補強電荷を展開する平面波のカットオフエネルギー。E...
GMAX=5.00 GMAXP=15.0 の条件で計算を行う場合は、一般的にpb...
GMAX=5.50-6.0 GMAXP=20.0 の条件で計算を行う場合はpbe1のps...
※これは一般的なルールなのでエネルギーと力のカットオフエネ...
-NTYP: 原子の種類~
-NATM: 空間反転対称性(NINV=1)がある場合の非等価な原子の...
-NATM2:全原子数、空間反転対称性が無い場合はNATM=NATM2であ...
*** 格子の対称性と種類 [#u75baf4e]
1 0 : space group num...
-NUM_SPACE_GROUP: 空間群の番号~
-TYPE: ブラベー格子の種類を番号で指定:
//-- 0 : simple, 1: body center, 2: face center, 3: a-fac...
|#|bravais lattice type|h
|0|simple|
|1|body-centered|
|2|face-centered|
|3|a-face-centered|
|4|b-face-centered|
|5|c-face-centered|
|6|rhombohedral|
*** [#wf17c7ac]
6.00 4.00 4.00 90.00 90.00 90.00 : a, b, c, alpha,...
-a, b, c: 格子ベクトルの長さ~
-alpha, beta, gamma: 格子ベクトルのなす角度
これらの代わりに以下のように"Cartesian"と記述するとそれ以...
Cartesian
6.00 0.00 0.00
0.00 4.00 0.00
0.00 0.00 4.00
*** [#db665620]
1 1 1 1 1 1 : knx, kny, knz, ...
Brillouin zone積分で用いるk点のメッシュ(knx, kny, knz)...
4-6カラムでのk点のオフセットの指定方法は、ガンマ点を含む...
六方晶の場合には対称性の破れが起こるためk点はオフセットし...
*** 原子位置の指定方法と空間反転対称性 [#h3df27cd]
1 0 : NCORD, NINV
//-原子位置の指定方法と空間反転対称性の有無を指定~
-NCORD: 1: カーテシアン座標で原子座標を入力(単位はBohr)...
-NINV: 0: 系に空間反転対称性がない、1: 空間反転対称性があ...
このオプションを使用すると、空間反転対称性が存在する場合...
*** 原子位置と原子種 [#o40afb1e]
0.0000 0.0000 0.0000 1 1 1 : cps, iwei, imdt...
2.2000 0.0000 0.0000 1 1 2 : cps, iwei, imdt...
//- 原子位置と原子種の指定
- 1-3カラム: cps(katm,3) (pos(katm,3)): 原子座標。NCORDが...
- 4カラム: iwei(katm): 空間反転対称性について等価な原子数
- 5カラム: imdtyp(katm): Hellmann-Feynman力に従って原子を...
- 6カラム: ityp(katm): 原子の種類を番号で指定。原子番号等...
*** 原子(擬イオン)に関する情報 [#w5f29409]
6 0.1500 51577.50 3 1 0.d0 : IATOMN, ALFA, A...
8 0.1500 51577.50 3 1 0.d0 : IATOMN, ALFA, A...
//-原子(擬イオン)に関する情報を入力する。意味のある値は
-1カラム: 原子番号
-2カラム: 初期電荷(ダミー)
-3カラム: 原子の質量(Da(旧のa.m.u.)(=水素原子が約1))。後...
有限温度の分子動力学を実行するときには現実の質量に設定す...
- 4, 5カラム: Local potentialの角運動量とウルトラソフト擬...
- 6カラム: zeta=(rho(up)-rho(down))/(rho(up)+rho(down))の...
-6カラム: 磁化率の初期値(一般的に分極してそうなものは0.5...
-その他: ilocはローカルポテンシャルの軌道量子数l+1、ivan...
*** 計算の継続(初期化)に関するオプション [#e0c8091e]
0 0 0 0 0 : ICOND, INIPOS, ...
|ICOND||h
|0|初期波動関数は乱数から発生させる。scratchからSCF計算を...
|1|初期波動関数と電荷密度は前の計算のものを使用する。SCF...
|2|初期波動関数は乱数から、電荷密度は前の計算のものを使用...
|3|初期波動関数とポテンシャルは古いものから出発するが、電...
|4|電荷密度をfixしたままSCF計算を行う。ICOND=2と同じ。|
|9, 11|実空間での電荷密度を出力. (9は全電荷密度、11はsof...
|10|STM|
|12|DOS_CAL 電子状態を読み込んで修正計算せずにDOSを計算す...
|14, 24|[[Partial density of state (PDOS)>計算コードの使...
|15|WFN_FFT 実空間での波動関数の出力|
|115|WFN_FFT_BAND 実空間での波動関数の出力 (バンド計算 (I...
|17, 117|[[Crystal Orbital Overlap Population (COOP Analy...
|18|CHG_DIFF_ANA 電荷密度差の表示(全系と、それを2つに分...
|20|INITIAL_MG 詳細不明|
|22, 23|[[BAND STRUCTURE>/計算コードの使い方/STATE/BandSt...
|33, 133 (13)|ALDOS XY面に平行な面でSuperCellをスライスし...
|40|規約ブリルアンゾーン内で波動関数(zag.data)と電荷密...
|41|41: nfkpt.dataで指定した経路上で波動関数と電荷密度をG...
~
|INIPOS|原子位置|h
|0|原子位置を入力ファイルから読み込む|
|1|restart.dataから原子座標を読み込んでMDまたは構造最適化...
|2|GEOMETRYから原子座標を読み込んで計算を継続する(restar...
- INIVEL: 原子の速度
--INIPOSと同様、1の場合はrestart.dataから、2の場合はGEOME...
- ININOS: 能勢の温度制御
- INIACC: accumulator:物理量の統計平均の計算の継続に関す...
計算の継続する例 (1): SCF計算の継続
-SCF計算のみを継続したいとき(geometryは入力ファイルから...
1 0 0 0 0 : ICOND,INIPOS,IN...
-構造最適化を継続 (restart.dataを参照)
1 1 0 0 0 : ICOND,INIPOS,IN...
-構造最適化を継続 (GEOMETRYを参照 (restart.dataは構造最適...
1 2 0 0 0 : ICOND,INIPOS,INI...
-構造最適化を継続 (restart.dataを参照、波動関数を初期化)
0 1 0 0 0 : ICOND,INIPOS,IN...
-MDの継続(IMDALG <=0)
1 1 1 0 0 : ICOND,INIPOS,IN...
*** ストレス計算、出力レベル[#q99d8327]
0 1 : IPRE, IPRI
- 圧力(ストレス)の計算(現在は実装されていないのでダミ...
- 標準出力のレベル(通常は1、1よりも大きいと多くの情報が...
*** 繰り返しの回数と計算時間に関するオプション [#wd00e91f]
200 200 0 57200.00 0 : NMD1, NMD2, ite...
- NMD1: 電子系のSCF計算の回数の上限を指定する。
- NMD2: 分子動力学 (構造最適化) のステップ数の上限を指定...
- iter_last:
- CPUMAX: CPU時間の上限を秒単位で指定
- ifstop:
*** ミキシング [#p971f063]
3 1 : way_mix, mix_what
|way_mix|ミキシングの方法|h
|1|simple mixing|
|2|Broyden|
|3|Broyden2|
|4|DFP|
|5|Pulay|
|6|Blugel|
~
|mix_what|ミキシングの対象|h
|1|電荷密度|
|2|ポテンシャル|
*** ミキシングに関するパラメータ [#na67c744]
0 8 0.8 : starting mixing...
//-ミキシングに関するパラメータ
- start mixing: dummy
- kbmix: ミキシングに用いる履歴の回数
- alpha=mix_alpha: ミキシングパラメータ
*** RMMのミキシングに関するパラメータ [#m243cb27]
0.60 0.50 0.60 0.70 1.00 : DTIM1, DTIM2, D...
- DTIM2: RMMで使用、その他はダミー変数.
- dtim_last: ダミー変数
*** 構造最適化または分子動力学(MD)に関するオプション [#c9...
30.00 2 1 0.10D-08 1.d-06 : DTIO,IMDALG,IEX...
- DTIO: イオン系のMDで用いる時間刻み(タイムステップ)(a....
- IMDALG: MDまたは構造最適化の方法を指定。正の値は構造最...
- [[Molecular dynamics>計算コードの使い方/STATE/Molecular...
|IMDALG|構造最適化の方法|h
|1|Newtonian dynamics|
|2|Quenched molecular dynamics|
|3|Vibrational mode analysis(nfvibrate.dataというファイ...
|4|GDIIS|
|5|TS search by GDIIS|
|6|NEB [[Nudged Elastic Band method>計算コードの使い方/ST...
|7|CINEB [[Climbing Image NEB method>計算コードの使い方/S...
|0|Newtonian dynamics|
|-1|有限温度のNewtonian dynamics(MVELSC=0:Microcanonical...
|-2|Langevin MD|
- IEXPL: dummy
- EDELTA: 電子系の収束の閾値を規定。収束は一原子あたりの...
*** [#y79be363]
0.0010 0.10D+02 0 : WIDTH, FORCCR, ...
- WIDTH: Fermi面のぼかしの幅. 正の場合はデルタ関数を多項...
- FORCCR: Forceの閾値。Forceがこの値よりも小さくなると構...
- ISTRESS: ストレスの計算用オプション。現在は使用されてい...
*** [#e2571185]
ggapbe 1 : XCTYPE, nspin
-XCTYPE: 交換相関汎関数を指定.
-- ggapbe: Perdew, Burke, Ernzerhof GGA (1996)
-- ggapw91: Perdew-Wang GGA (1991)
-- ldapw91: Perdew-Wang L(S)DA (1991)
-- RPBE: revised PBE of Hammer et al.
-- revPBE: revised PBE of Zhang and Yang
-- WC: Wu-Cohen GGA
-- vdW-DF1/vdW-DF2/optB88-vdW/optB86b-vdW/rev-vdW-DF2
-NSPIN: 1ならば非磁性、2でcollinearな磁性体の計算を行う。
*** STM用オプション [#k57e8267]
1.00 3 : destm, n_stm
*** [#f318170e]
102 : NBZTYP
- WIDTH < -10.0の時のみ有効。101が推奨値。
*** ダミー行 [#b8e34609]
0 0 0 : NKX, NKY, NKZ...
0 0 0 : NKX2, NKY2, NKZ...
*** 計算に用いるバンドの数 [#z175352f]
8 : NEG (# of bands)
STATEの計算では、常に非占有状態を含めて計算を行う必要があ...
*** nonlocal pseudopotentialの計算方法 [#m2ffee89]
1 : NEXTST (1: G-sp...
- 0: R-space (Davidsonの時は使用しない)
- 1: G-space
*** ダミー行 [#p6cd59ce]
0 : 0; random numbe...
*** 対角化法の選択 [#z0450187]
2 : imsd (2: Davids...
- 1: RMM-DIIS
- 2: Davidson
大規模な計算をするときにはRMM-DIISとReal space projection...
*** [#oa32d26e]
0 : eval. eko diff....
前回のiterationにおける擬固有値との比較を行うかどうかを指...
*** PDOSの計算オプション [#jb83af81]
0 ...
0でない場合にPDOSの計算を行う。但しnpdosaoは原子の数より...
*** f電子用の半経験的パラメータ用のオプション [#ya0ef43c]
0 0.0 ...
1コラム目は整数を、2コラム目は実数を入力。現在使用してい...
* &aname(AOLDOSsetting){PDOS/AOLDOSの設定}; [#vd7c4b40]
AOLDOS。Layer-resolved DOSの計算の指定は以下を参照。
*** 原子軌道に射影したPDOS(Partial Density Of States、も...
(以下、pdos用のオプション。仮にnpdosao = 5とし、通し番号...
5 : NPDOSAO
1
3
4
9
10
-15.00 5.00 0.20 501 :EPDOS(1), EPDOS(...
2.2 0.3 :RAD, WIDTH FOR T...
1.0 0.3 :RAD, WIDTH FOR T...
0.2 12 :DR, NR
&OTHERS
GAUSSDOS
&END
- npdosao: DOSを投影させる原子数。その次の行にnpdosao個の...
- EPDOS(1,2): pdosを表示するエネルギー領域(Fermi準位との...
- EPDOS(3): gaussianを用いる場合の標準偏差、tetrahedronを...
-- GaussianBroadeningを使うためには次のオプション指定が必...
&OTHERS
GAUSSDOS
&END
- NPDOSE: EPDOS(1)~EPDOS(2)の間のエネルギーメッシュ数
- rad: DOSを投影させる原子中心の波動関数(擬波動関数)の...
(bohr単位。通常、原子の価電子半径程度を指定。e.g. H:0.6 ...
- width: 上記のrad前後で、原子の擬波動関数をFermi分布関...
- DR,NR: (おそらく)未使用~
(元々は、擬波動関数同士やプロジェクターbetaとの重なり積分...
* Layer-resolved PDOSの計算 [#accbb8c2]
xy面に平行な面でスライスしてできた領域内に含まれる電子のD...
いったん電子状態を求めてzaj.data,potential.dataを出力した...
&ALDOS
PRINT_WEIGHT
FERMI_LEVEL -0.14395389
NUMBER_OF_LAYER 7
LAYER
-7.7938
-5.3974
-2.9894
-0.5968
1.41385
3.832933
6.64160692845
10.60270297025
EMIN -15.D0
EMAX 5.D0
EWIDTH 0.2D0
NPDOSE 401
&END
終了行:
* STATE クイックリファレンス [#e91b1b9c]
** 入力ファイルの例 (CO分子) [#s044357d]
0 0 0 0 0 0 : dummy line (6 i...
5.50 20.00 2 2 2 : GMAX, GMAXP, NT...
1 0 : space group num...
6.00 4.00 4.00 90.00 90.00 90.00 : a, b, c, alpha,...
1 1 1 1 1 1 : knx, kny, knz, ...
1 0 : NCORD, NINV
0.0000 0.0000 0.0000 1 1 1 : cps, iwei, imdt...
2.2000 0.0000 0.0000 1 1 2 : cps, iwei, imdt...
6 0.1500 51577.50 3 1 0.d0 : IATOMN, ALFA, A...
8 0.1500 51577.50 3 1 0.d0 : IATOMN, ALFA, A...
0 0 0 0 0 : ICOND, INIPOS, ...
0 1 : IPRE, IPRI
200 200 0 57200.00 0 : NMD1, NMD2, ite...
3 1 : way_mix, mix_what
0 8 0.8 : starting mixing...
0.60 0.50 0.60 0.70 1.00 : DTIM1, DTIM2, D...
30.00 2 1 0.10D-08 1.d-06 : DTIO, IMDALG, I...
0.0010 0.10D+02 0 : WIDTH, FORCCR, ...
ggapbe 1 : XCTYPE, nspin
1.00 3 : destm, n_stm
102 : NBZTYP
0 0 0 : NKX, NKY, NKZ...
0 0 0 : NKX2, NKY2, NKZ...
8 : NEG (# of bands)
1 : NEXTST (1: G-sp...
0 : 0; random numbe...
2 : imsd (2: Davids...
0 : eval. eko diff....
0 : npdosao
0 0.0 : SM_N doping
*** ダミー行 [#b28085e8]
0 0 0 0 0 0 : dummy line (6 i...
6個の整数を指定。歴史的経緯で存在するダミー行。~
※過去のバージョンでは1カラム目の変数が1の場合は基本並進ベ...
※現状STATEでは無視されているが、ユーティリティプログラムr...
*** カットオフと原子数 [#lad659e0]
5.50 20.00 2 2 2 : GMAX, GMAXP, NT...
-GMAX: 波動関数を展開する平面波のカットオフエネルギーを指...
-GMAXP: 補強電荷を展開する平面波のカットオフエネルギー。E...
GMAX=5.00 GMAXP=15.0 の条件で計算を行う場合は、一般的にpb...
GMAX=5.50-6.0 GMAXP=20.0 の条件で計算を行う場合はpbe1のps...
※これは一般的なルールなのでエネルギーと力のカットオフエネ...
-NTYP: 原子の種類~
-NATM: 空間反転対称性(NINV=1)がある場合の非等価な原子の...
-NATM2:全原子数、空間反転対称性が無い場合はNATM=NATM2であ...
*** 格子の対称性と種類 [#u75baf4e]
1 0 : space group num...
-NUM_SPACE_GROUP: 空間群の番号~
-TYPE: ブラベー格子の種類を番号で指定:
//-- 0 : simple, 1: body center, 2: face center, 3: a-fac...
|#|bravais lattice type|h
|0|simple|
|1|body-centered|
|2|face-centered|
|3|a-face-centered|
|4|b-face-centered|
|5|c-face-centered|
|6|rhombohedral|
*** [#wf17c7ac]
6.00 4.00 4.00 90.00 90.00 90.00 : a, b, c, alpha,...
-a, b, c: 格子ベクトルの長さ~
-alpha, beta, gamma: 格子ベクトルのなす角度
これらの代わりに以下のように"Cartesian"と記述するとそれ以...
Cartesian
6.00 0.00 0.00
0.00 4.00 0.00
0.00 0.00 4.00
*** [#db665620]
1 1 1 1 1 1 : knx, kny, knz, ...
Brillouin zone積分で用いるk点のメッシュ(knx, kny, knz)...
4-6カラムでのk点のオフセットの指定方法は、ガンマ点を含む...
六方晶の場合には対称性の破れが起こるためk点はオフセットし...
*** 原子位置の指定方法と空間反転対称性 [#h3df27cd]
1 0 : NCORD, NINV
//-原子位置の指定方法と空間反転対称性の有無を指定~
-NCORD: 1: カーテシアン座標で原子座標を入力(単位はBohr)...
-NINV: 0: 系に空間反転対称性がない、1: 空間反転対称性があ...
このオプションを使用すると、空間反転対称性が存在する場合...
*** 原子位置と原子種 [#o40afb1e]
0.0000 0.0000 0.0000 1 1 1 : cps, iwei, imdt...
2.2000 0.0000 0.0000 1 1 2 : cps, iwei, imdt...
//- 原子位置と原子種の指定
- 1-3カラム: cps(katm,3) (pos(katm,3)): 原子座標。NCORDが...
- 4カラム: iwei(katm): 空間反転対称性について等価な原子数
- 5カラム: imdtyp(katm): Hellmann-Feynman力に従って原子を...
- 6カラム: ityp(katm): 原子の種類を番号で指定。原子番号等...
*** 原子(擬イオン)に関する情報 [#w5f29409]
6 0.1500 51577.50 3 1 0.d0 : IATOMN, ALFA, A...
8 0.1500 51577.50 3 1 0.d0 : IATOMN, ALFA, A...
//-原子(擬イオン)に関する情報を入力する。意味のある値は
-1カラム: 原子番号
-2カラム: 初期電荷(ダミー)
-3カラム: 原子の質量(Da(旧のa.m.u.)(=水素原子が約1))。後...
有限温度の分子動力学を実行するときには現実の質量に設定す...
- 4, 5カラム: Local potentialの角運動量とウルトラソフト擬...
- 6カラム: zeta=(rho(up)-rho(down))/(rho(up)+rho(down))の...
-6カラム: 磁化率の初期値(一般的に分極してそうなものは0.5...
-その他: ilocはローカルポテンシャルの軌道量子数l+1、ivan...
*** 計算の継続(初期化)に関するオプション [#e0c8091e]
0 0 0 0 0 : ICOND, INIPOS, ...
|ICOND||h
|0|初期波動関数は乱数から発生させる。scratchからSCF計算を...
|1|初期波動関数と電荷密度は前の計算のものを使用する。SCF...
|2|初期波動関数は乱数から、電荷密度は前の計算のものを使用...
|3|初期波動関数とポテンシャルは古いものから出発するが、電...
|4|電荷密度をfixしたままSCF計算を行う。ICOND=2と同じ。|
|9, 11|実空間での電荷密度を出力. (9は全電荷密度、11はsof...
|10|STM|
|12|DOS_CAL 電子状態を読み込んで修正計算せずにDOSを計算す...
|14, 24|[[Partial density of state (PDOS)>計算コードの使...
|15|WFN_FFT 実空間での波動関数の出力|
|115|WFN_FFT_BAND 実空間での波動関数の出力 (バンド計算 (I...
|17, 117|[[Crystal Orbital Overlap Population (COOP Analy...
|18|CHG_DIFF_ANA 電荷密度差の表示(全系と、それを2つに分...
|20|INITIAL_MG 詳細不明|
|22, 23|[[BAND STRUCTURE>/計算コードの使い方/STATE/BandSt...
|33, 133 (13)|ALDOS XY面に平行な面でSuperCellをスライスし...
|40|規約ブリルアンゾーン内で波動関数(zag.data)と電荷密...
|41|41: nfkpt.dataで指定した経路上で波動関数と電荷密度をG...
~
|INIPOS|原子位置|h
|0|原子位置を入力ファイルから読み込む|
|1|restart.dataから原子座標を読み込んでMDまたは構造最適化...
|2|GEOMETRYから原子座標を読み込んで計算を継続する(restar...
- INIVEL: 原子の速度
--INIPOSと同様、1の場合はrestart.dataから、2の場合はGEOME...
- ININOS: 能勢の温度制御
- INIACC: accumulator:物理量の統計平均の計算の継続に関す...
計算の継続する例 (1): SCF計算の継続
-SCF計算のみを継続したいとき(geometryは入力ファイルから...
1 0 0 0 0 : ICOND,INIPOS,IN...
-構造最適化を継続 (restart.dataを参照)
1 1 0 0 0 : ICOND,INIPOS,IN...
-構造最適化を継続 (GEOMETRYを参照 (restart.dataは構造最適...
1 2 0 0 0 : ICOND,INIPOS,INI...
-構造最適化を継続 (restart.dataを参照、波動関数を初期化)
0 1 0 0 0 : ICOND,INIPOS,IN...
-MDの継続(IMDALG <=0)
1 1 1 0 0 : ICOND,INIPOS,IN...
*** ストレス計算、出力レベル[#q99d8327]
0 1 : IPRE, IPRI
- 圧力(ストレス)の計算(現在は実装されていないのでダミ...
- 標準出力のレベル(通常は1、1よりも大きいと多くの情報が...
*** 繰り返しの回数と計算時間に関するオプション [#wd00e91f]
200 200 0 57200.00 0 : NMD1, NMD2, ite...
- NMD1: 電子系のSCF計算の回数の上限を指定する。
- NMD2: 分子動力学 (構造最適化) のステップ数の上限を指定...
- iter_last:
- CPUMAX: CPU時間の上限を秒単位で指定
- ifstop:
*** ミキシング [#p971f063]
3 1 : way_mix, mix_what
|way_mix|ミキシングの方法|h
|1|simple mixing|
|2|Broyden|
|3|Broyden2|
|4|DFP|
|5|Pulay|
|6|Blugel|
~
|mix_what|ミキシングの対象|h
|1|電荷密度|
|2|ポテンシャル|
*** ミキシングに関するパラメータ [#na67c744]
0 8 0.8 : starting mixing...
//-ミキシングに関するパラメータ
- start mixing: dummy
- kbmix: ミキシングに用いる履歴の回数
- alpha=mix_alpha: ミキシングパラメータ
*** RMMのミキシングに関するパラメータ [#m243cb27]
0.60 0.50 0.60 0.70 1.00 : DTIM1, DTIM2, D...
- DTIM2: RMMで使用、その他はダミー変数.
- dtim_last: ダミー変数
*** 構造最適化または分子動力学(MD)に関するオプション [#c9...
30.00 2 1 0.10D-08 1.d-06 : DTIO,IMDALG,IEX...
- DTIO: イオン系のMDで用いる時間刻み(タイムステップ)(a....
- IMDALG: MDまたは構造最適化の方法を指定。正の値は構造最...
- [[Molecular dynamics>計算コードの使い方/STATE/Molecular...
|IMDALG|構造最適化の方法|h
|1|Newtonian dynamics|
|2|Quenched molecular dynamics|
|3|Vibrational mode analysis(nfvibrate.dataというファイ...
|4|GDIIS|
|5|TS search by GDIIS|
|6|NEB [[Nudged Elastic Band method>計算コードの使い方/ST...
|7|CINEB [[Climbing Image NEB method>計算コードの使い方/S...
|0|Newtonian dynamics|
|-1|有限温度のNewtonian dynamics(MVELSC=0:Microcanonical...
|-2|Langevin MD|
- IEXPL: dummy
- EDELTA: 電子系の収束の閾値を規定。収束は一原子あたりの...
*** [#y79be363]
0.0010 0.10D+02 0 : WIDTH, FORCCR, ...
- WIDTH: Fermi面のぼかしの幅. 正の場合はデルタ関数を多項...
- FORCCR: Forceの閾値。Forceがこの値よりも小さくなると構...
- ISTRESS: ストレスの計算用オプション。現在は使用されてい...
*** [#e2571185]
ggapbe 1 : XCTYPE, nspin
-XCTYPE: 交換相関汎関数を指定.
-- ggapbe: Perdew, Burke, Ernzerhof GGA (1996)
-- ggapw91: Perdew-Wang GGA (1991)
-- ldapw91: Perdew-Wang L(S)DA (1991)
-- RPBE: revised PBE of Hammer et al.
-- revPBE: revised PBE of Zhang and Yang
-- WC: Wu-Cohen GGA
-- vdW-DF1/vdW-DF2/optB88-vdW/optB86b-vdW/rev-vdW-DF2
-NSPIN: 1ならば非磁性、2でcollinearな磁性体の計算を行う。
*** STM用オプション [#k57e8267]
1.00 3 : destm, n_stm
*** [#f318170e]
102 : NBZTYP
- WIDTH < -10.0の時のみ有効。101が推奨値。
*** ダミー行 [#b8e34609]
0 0 0 : NKX, NKY, NKZ...
0 0 0 : NKX2, NKY2, NKZ...
*** 計算に用いるバンドの数 [#z175352f]
8 : NEG (# of bands)
STATEの計算では、常に非占有状態を含めて計算を行う必要があ...
*** nonlocal pseudopotentialの計算方法 [#m2ffee89]
1 : NEXTST (1: G-sp...
- 0: R-space (Davidsonの時は使用しない)
- 1: G-space
*** ダミー行 [#p6cd59ce]
0 : 0; random numbe...
*** 対角化法の選択 [#z0450187]
2 : imsd (2: Davids...
- 1: RMM-DIIS
- 2: Davidson
大規模な計算をするときにはRMM-DIISとReal space projection...
*** [#oa32d26e]
0 : eval. eko diff....
前回のiterationにおける擬固有値との比較を行うかどうかを指...
*** PDOSの計算オプション [#jb83af81]
0 ...
0でない場合にPDOSの計算を行う。但しnpdosaoは原子の数より...
*** f電子用の半経験的パラメータ用のオプション [#ya0ef43c]
0 0.0 ...
1コラム目は整数を、2コラム目は実数を入力。現在使用してい...
* &aname(AOLDOSsetting){PDOS/AOLDOSの設定}; [#vd7c4b40]
AOLDOS。Layer-resolved DOSの計算の指定は以下を参照。
*** 原子軌道に射影したPDOS(Partial Density Of States、も...
(以下、pdos用のオプション。仮にnpdosao = 5とし、通し番号...
5 : NPDOSAO
1
3
4
9
10
-15.00 5.00 0.20 501 :EPDOS(1), EPDOS(...
2.2 0.3 :RAD, WIDTH FOR T...
1.0 0.3 :RAD, WIDTH FOR T...
0.2 12 :DR, NR
&OTHERS
GAUSSDOS
&END
- npdosao: DOSを投影させる原子数。その次の行にnpdosao個の...
- EPDOS(1,2): pdosを表示するエネルギー領域(Fermi準位との...
- EPDOS(3): gaussianを用いる場合の標準偏差、tetrahedronを...
-- GaussianBroadeningを使うためには次のオプション指定が必...
&OTHERS
GAUSSDOS
&END
- NPDOSE: EPDOS(1)~EPDOS(2)の間のエネルギーメッシュ数
- rad: DOSを投影させる原子中心の波動関数(擬波動関数)の...
(bohr単位。通常、原子の価電子半径程度を指定。e.g. H:0.6 ...
- width: 上記のrad前後で、原子の擬波動関数をFermi分布関...
- DR,NR: (おそらく)未使用~
(元々は、擬波動関数同士やプロジェクターbetaとの重なり積分...
* Layer-resolved PDOSの計算 [#accbb8c2]
xy面に平行な面でスライスしてできた領域内に含まれる電子のD...
いったん電子状態を求めてzaj.data,potential.dataを出力した...
&ALDOS
PRINT_WEIGHT
FERMI_LEVEL -0.14395389
NUMBER_OF_LAYER 7
LAYER
-7.7938
-5.3974
-2.9894
-0.5968
1.41385
3.832933
6.64160692845
10.60270297025
EMIN -15.D0
EMAX 5.D0
EWIDTH 0.2D0
NPDOSE 401
&END
ページ名: