molecule追加された行はこの色です。削除された行はこの色です。integers)歴史的経緯でそのままになっている。gammaで与えていたようである。~※バージョン539tでは活きているようです。※現状STATEではおそらく無視されている。repeat.fではxyzのファイルを作る際に参照されている。STATEプログラムの改修者へ:この項目を入力から削除すると面倒なことになるので当面このパラメータを削除しないように。波動関数を展開する平面波のカットオフエネルギーを指定。補強電荷を展開する平面波のカットオフエネルギー。potentialセットを用いるとよい~※これは一般的なルールなのでエネルギーと力のカットオフエネルギーについての収束性は必ずチェックすべきである。無い場合はNATM=NATM2である(現在は常にNATM=NATM2とすること)。格子ベクトルのなす角度これらの代わりに以下のように"Cartesian"と記述するとそれ以下の行で格子ベクトルを直接カーテシアン座標で指定することができる。knz)を最初の3カラムで、k点のオフセットを次の3カラムで指定する。4-6カラムでのk点のオフセットの指定方法は、ガンマ点を含む場合には1を、ガンマ点を含まない場合(Monkhorst-Pack)は2を指定する。六方晶の場合にはk点はオフセットしない方が良いようである。空間反転対称性がある。~このオプションを使用すると、空間反転対称性が存在する場合には明示的に考慮する原子数をおよそ半分にすることが可能で、計算時間の短縮をすることが可能である。(FTMD)の場合、1000以上の数を指定するとその原子は熱浴に接することを指定。原子の種類を番号で指定。原子番号等は後に続く列に記述する。ZETA1//-原子(擬イオン)に関する情報を入力する。原子の質量(Da(旧のa.m.u.)(=水素原子が約1))。後述するが、上の値は現実のものとは異なる。MD(IMDALG=2)を行うときには原子質量は同じ値、例えばCの12.0107やSiの28.0855などにしておくと計算を安定して行うことができる。有限温度の分子動力学を実行するときには現実の質量に設定するべきである。zeta=(rho(up)-rho(down))/(rho(up)+rho(down))の初期値を与える。スピン分極の計算を行う際はnspin=2としてzetaとしては0でない有限の値を設定する。zeta=0とすると非磁性の解に収束するだろう。pseudopotentialを使用するかどうか、といったパラメータであるが、入力ファイルからではなく擬ポテンシャルのファイルから読むので実際にはダミーパラメータである。ちなみにAlFAという変数は現在は使用されていない。電荷密度を乱数から発生させるときに使用していたようである。INIACC|ICOND||h|0|初期波動関数は乱数から発生させる。scratchからSCF計算を開始する時に使用する。||1|初期波動関数と電荷密度は前の計算のものを使用する。SCF計算を継続する際に使用する。zaj.dataとpotential.dataが必要||2|初期波動関数は乱数から、電荷密度は前の計算のものを使用するが、電荷密度(ポテンシャル)はfixしたまま計算を行う。ポテンシャルをfixした計算(ICOND=2)を継続する際に用いる。potential.dataとzaj.dataが必要。||4|電荷密度をfixしたままSCF計算を行う。ICOND=2と同じ。バンド計算用(23は継続計算)。XY面に平行な面でSuperCellをスライスした断片各々のDOS(133はK点分解で出力、13は旧バージョン(使用しない))範囲の指定にはnfaldos.dataという名前のファイルを用意する。下記AOLDOSの後ろに記述する。XY面に平行な面でSuperCellをスライスした断片各々のDOS(133はK点分解で出力、13は旧バージョン(使用しない))範囲の指定にはnfaldos.dataという名前のファイルを用意する。下記AOLDOSの後ろに記述する。5.3.8b用)|~|INIPOS|原子位置|h|0|原子位置を入力ファイルから読み込む||1|restart.dataから原子座標を読み込んでMDまたは構造最適化を継続する||2|GEOMETRYから原子座標を読み込んで計算を継続する(restart.dataが必要)。restart.dataはバイナリーのデータで扱いにくいが、GEOMETRYはascii形式であるので、比較的扱いやすい。計算を現在のgeometryと異なるものから始めるときに使用すると便利である。原子の速度--INIPOSと同様、1の場合はrestart.dataから、2の場合はGEOMETRYから速度を読み込む。INIPOSの時と同様にrestart.dataが必要。標準出力のレベル(通常は1、1よりも大きいと多くの情報が出力される。電子系のSCF計算の回数の上限を指定する。分子動力学のステップ数の上限を指定する。NMD2+1が実際のMD(構造最適化)の回数の上限である。イオン系のMDで用いる時間刻み(タイムステップ)(a.u.)。(IMDALG=2)またはMDで有効。MDまたは構造最適化の方法を指定。dynamics]]も参考にしてください。電子系の収束の閾値を規定。収束は一原子あたりの全エネルギーで評価している。正の場合はデルタ関数を多項式で、負の場合にはHermite-Gaussianの一次式で近似する。この時はNBZTYPが有効になる。Forceがこの値よりも小さくなると構造最適化を終了する。IMDALG>0のみ有効である。ストレスの計算用オプション。現在は使用されていない。1ならば非磁性、2でcollinearな磁性体の計算を行う。-10.0の時のみ有効。bands)STATEの計算では、常に非占有状態を含めて計算を行う必要がある。(電子数/2)+4などとするのが良いだろう。projectionの組み合わせ(NEXTST=0、IMSD=1)が推奨されますが、初期波動関数は乱数から発生させるものではなく、Davidsonである程度収束させたものを使うのが良いでしょう。しなければ最安定な電子状態に収束しない可能性がある。yes前回のiterationにおける擬固有値との比較を行うかどうかを指定。npdosao0でない場合にPDOSの計算を行う。npdosaoは原子の数より小さくなければならない。molecule#lfe4ea7f]]を参照のこと。SM_dopping現在使用している人はいないため十分テストされていません。[#vd7c4b40]AOLDOS。DOSの計算の指定は、ここ野で説明ではなく次項を参照。[#eea82346](以下、pdos用のオプション。5とし、通し番号が1,3,4,9,10の原子の軌道に投影したDOSを計算、原子の種類は2種類。DOSを投影させる原子数。GaussianBroadeningを使うためには次のオプション指定が必要なので注意。DOSを投影させる原子中心の波動関数(擬波動関数)の、動径方向のカットオフ値を指定~(bohr単位。通常、原子の価電子半径程度を指定。[#accbb8c2]xy面に平行な面でスライスしてできた領域内に含まれる電子のDOSを算出する。電子状態を求めてzaj.data,potential.dataを出力した後に、ICOND=33として実行する。下記のようなnfaldos.dataというファイルを用意する。フェルミエネルギーはHartree単位で与える(直前のSCF計算の最後に得られるFermiEnergyをnfoutから読み取って書く)、Layerの境界はBohr単位のz座標である。Layerの数より一つ大きい数のz座標を与える。factor)を与える。