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

• バックアップ一覧
• 差分 を表示
• 現在との差分 を表示
• ソース を表示
• Ni へ行く。
  • 1 (2016-03-02 (水) 17:30:24)
  • 2 (2016-03-02 (水) 23:56:59)
  • 3 (2016-03-05 (土) 19:17:49)
  • 4 (2016-03-06 (日) 17:51:28)

bulk Ni †

入力ファイル †

0 0 0 0 0 0
5.00000 15.0000 1 1 1                           : gmax gmaxp ntyp natm natm2
221 2                                           : num_space_group, type
6.700000 6.700000 6.700000 90.000 90.000 90.000 : a b c alpha beta gamma
12 12 12 1 1 1                                  : k_mesh
0 0                                             : ncord ninv
0.000000d0 0.000000d0 0.000000d0 1 0 1          : iwei imdtyp ityp
28 0.5000 22104.64 6 1 0.1                      : type 1atomn alfa amion iloc ivan ... (1)
0 0 0 0 0                                       : icond
0 1                                             : ipre ipri
30 30 0 84200.00 0                              : nmd1 nmd2 last_iter cpumax ifstop
6 1                                             : SIMPLE=1,BROYD1=2,BROYD2=3,DFP=4,PULAY=5
0 20 0.30                                       : iter_start  kbxmix mix_alpha
0.20 0.30 0.20 0.20 0.20                        : dtim1 dtim2 dtim3 dtim4 dtim
300.00 4 1 0.50D-09                             : dtio imdalg iexpl edelta
-0.0020 0.50D+03 0                              : width forccr istress
ggapbe 2                                        : xctype,kspin ....................... (2)
2.00                                            : destm
101                                             : nbztyp
4 4 4
4 4 4
10                                              : keg
1                                               : nextst(mb)
0                                               : 0 random numbers 1: matrix
2                                               : imsd  1=MSD, 2=Preconditioned
0                                               : evaluate_eko_diff
1                                               : npdosao ............................ (3)
1
-15.00 5.00 0.20 501
2.40000 0.200000
0.20000 14
0 0.000                                         : sm_n, dopping

説明 †

• (2) ニッケルは磁性体であるためkspin=2を指定してそれぞれのスピンに対して収束計算を行う
• (1) 6番目のパラメータでNi原子のスピン分極密度
  \[ \zeta=\frac{\rho_\uparrow-\rho_\downarrow}{\rho_\uparrow+\rho_\downarrow} \]
  の初期値を指定する
• (3) npdosao=1を指定してブラベー格子中の1つの原子に対して各軌道に射影した状態密度を計算する

全エネルギー †

$ grep ETOT nfout_1|grep "[0-9]"
ETOT:   1    -48.03422182  0.4803E+02  0.1740E-01
ETOT:   2    -48.09018411  0.5596E-01  0.2913E-01
ETOT:   3    -48.25418459  0.1640E+00  0.1789E-01
ETOT:   4    -48.33744359  0.8326E-01  0.7657E-02
ETOT:   5    -48.35584074  0.1840E-01  0.1601E-02
ETOT:   6    -48.35584552  0.4775E-05  0.8964E-03
ETOT:   7    -48.35603490  0.1894E-03  0.4295E-03
ETOT:   8    -48.35609065  0.5575E-04  0.7521E-04
ETOT:   9    -48.35608964  0.1011E-05  0.1006E-03
ETOT:  10    -48.35609257  0.2925E-05  0.2678E-04
ETOT:  11    -48.35609278  0.2145E-06  0.4263E-05
ETOT:  12    -48.35609277  0.9573E-08  0.1797E-05
ETOT:  13    -48.35609278  0.5008E-08  0.2724E-06
ETOT:  14    -48.35609278  0.1381E-09  0.1338E-06
ETOT:  15    -48.35609278  0.9063E-09  0.4832E-07
ETOT:  16    -48.35609278  0.1112E-09  0.3447E-07
  ETOT(Q) + SM_energy =  -48.3560927768421       0.698797641137880

状態密度 †

各軌道に射影した状態密度 †