公開されていないものには、元々公開が不可能なものもあります、ベッセル 関数や、対角化、時間表示などのサブルーチンは、非公開の汎用ライブラリー 等のものを利用しているので、これらに関しては各自で確保してもらうしかな いと思います。
文字色が青色っぽいのが既に公開されているサブルーチンです。 +-main-----TIME------CLOCKM 時間計測用サブルーチン(なくてもよい) |--SIMP シンプソン積分用係数計算 |--INPUT1 座標データ入力用 |--INFOUT 計算情報出力(なくてもよい) |--LATTIC----HPSORT 基本逆格子ベクトル計算(LATTIC) |--GSTEP1 逆格子座標計算(2Gmax以下) |--HPSORT 逆格子座標を小さい順に並べる |--SYMM------OPGR------OPMTRX 対称性計算(SYMM) |--GSTEPF FFT用マッピング |--KSTEP-----KPMSF 対称性なし(表面用)サンプリングk点計算 | |--KPMWBZ 対称性なし(バルク用) | |--SCCM 単純立方用 | |--BCCM 体心立方用 | |--FCCM 面心立方用 | |--HEXM 六方最密充填用 | |--TETRAH 正方晶(ルチル構造)用 | |--APBO2 斜方晶(αPbO2構造)用 |--BASNUM 平面波用逆格子計算(Gmax以下) |--LATSCA LATTICストレス用(定圧) |--GSTSCA GSTEP1ストレス用(定圧) |--SYMSCA SYMMストレス用(定圧) |--GSFSCA GSTEPFストレス用(定圧) |--FORM フォームファクター計算用 |--PCC-------DSJNV 部分内殻補正用(PCC) |------DSJNV 局所擬ポテンシャル用(入力、フーリエ変換)(PSEUDO) |--PSELMD PSEUDO(非MDループ用) |--EWVEC エバルト計算用 |--EWVMD EWVEC(非MDループ用) |--INTCHG 初期電荷密度設定用 |--EVIN パラメータ入力(なくてもよい) |--XCFFT-----C3FFT <---- MFFT 交換相関項計算用(XCFTT) |-------DSJNV 非局所擬ポテンシャル部分(入力、フーリエ変換) |--DIAGON----CHOBSD 初期波動関数(固有ベクトル)設定〔対角化〕(DIAGON) |--KBINT-----DSJNV 非局所擬ポテンシャル(非初期用、フーリエ変換) |--MSD-------KBINT-----DSJNV カー・パリネロ本体(MSD) | |--FORZFB 非局所部分計算用(非MDループ用) | |--C3FFT <---- MFFT 高速フーリエ変換 | |--TIME------CLOCKM 時間計測(なくてもよい) |--FERMI-----WIDTH2 フェルミエネルギー計算用(FERMI) |--FORCE 力の計算(非局所部分)用 |--FORZFB 非局所部分計算用(非MDループ用) |--CHAVER----C3FFT <---- MFFT 電荷密度計算用(CHAVER) | |--CHGAVR 電荷密度対称性計算用 |--STRESS----XCFFT-----C3FFT <---- MFFT ストレス計算用(STRESS) | |--XSTPC-----C3FFT <---- MFFT (局所部分用) | |--STRNL ストレス計算(非局所部分用) |--FORLOC 力の計算用(局所部分) |--FORCES 力の対称性計算用 |--MD 古典的分子動力学部分用 |------XCFFT-----C3FFT <---- MFFT 全エネルギー計算用(ENERGY) |--CONV2 収束判定用 |--EVOUT 計算結果、パラメーター出力用(その1) |--EVOU2 計算結果、パラメーター出力用(その2)