第一原理分子動力学プログラム STATE Senri Wiki

• バックアップ一覧
• 差分 を表示
• 現在との差分 を表示
• ソース を表示
• 計算コードの使い方/STATE/Cu111 へ行く。
  • 1 (2019-01-07 (月) 06:47:04)
  • 2 (2019-01-07 (月) 07:33:34)
  • 3 (2019-01-07 (月) 15:44:16)
  • 4 (2019-01-07 (月) 22:41:40)
  • 5 (2019-01-28 (月) 15:47:16)

銅(111)表面 †

ここではPBEで最適化したバルクの格子定数を用いて作成した6原子層から成る銅(111)表面のSCF計算とバンド計算の例を示す。

SCF計算 †

• 入力ファイル (nfinp_scf)

    0  0  0  0  0   0
     6.00000 20.00000  1   6   6         : GMAX GMAXP NTYP NATM NATM2
     1    0                              : num_space_group braivis_lattice_type
   Cartesian
        4.865948314378      0.000000000000      0.000000000000
       -2.432974157189      4.214034853753      0.000000000000
        0.000000000000      0.000000000000     79.460603233215
     12  12  1  1  1  1                  : K_mesh
      1   0                              : NCORD NINV : IWEI  IMDTYP  ITYP
        0.000000000000      0.000000000000      9.932575404152    1    1    1
        0.000000000000      2.809356569169      5.959545242491    1    1    1
        0.000000000000     -2.809356569169      1.986515080830    1    1    1
        0.000000000000      0.000000000000     -1.986515080830    1    1    1
        0.000000000000      2.809356569169     -5.959545242491    1    1    1
        0.000000000000     -2.809356569169     -9.932575404152    1    1    1
     29  0.5000  22104.64 6 1 0.2        : IATOMN ALFA AMION ILOC IVAN
      0 0 0 0 0                          : ICOND INIPOS INIVEL ININOS INIACC
       0    1                            : ipre  ipri
      200   200    0   3200.00    0      : nmd1 nmd2 last_iter cpumax ifstop
       6    1                            : WAY_MIX MIX_WHAT
       0   20 0.30                       : iter_start  KBXMIX   MIX_ALPHA
       0.20  0.30  0.20  0.20  0.20      : dtim1 dtim2 dtim3 dtim4 dtim
     200.00     4     1    5.00D-12      : dtio  imdalg  iexpl  edelta
     -0.0010  1.00D+03    0              : width forccr istress
  ggapbe      1                          : xctype kspin
       2.00                              : destm
     101                                 : nbztyp
       4    4    4                       : DUMMY
       4    4    4                       : DUMMY
      48                                 : keg
       1                                 : nextst
       0                                 : 0 random numbers 1: matrix
       2                                 : imsd
       0                                 : evaluate_eko_diff
       0                                 : npdosao
       0   0.000                         : sm_n  dopping

  STATEの実行

  mpirun -np 8 ./STATE < nfinp_scf > nfout_scf

バンド計算 †