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} \]
    の初期値を0.1に指定する
  • (3) npdosao=1を指定してブラベー格子中の1つのNi原子の各軌道に射影した状態密度を計算する

全エネルギー

$ grep ETOT output_file | 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

状態密度

dos_Ni.png

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

$ grep AO_LDOS output_file
AO_LDOS: atm spn  energy   s     px     py     pz    dzz dxx-yy    dxy    dyz    dzx ...
AO_LDOS:   1 1 -15.00  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000 ...
AO_LDOS:   1 1 -14.96  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000 ...
AO_LDOS:   1 1 -14.92  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000 ...
AO_LDOS:   1 1 -14.88  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000 ...
AO_LDOS:   1 1 -14.84  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000 ...
AO_LDOS:   1 1 -14.80  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000 ...
AO_LDOS:   1 1 -14.76  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000 ...
AO_LDOS:   1 1 -14.72  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000 ...
AO_LDOS:   1 1 -14.68  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000 ...
AO_LDOS:   1 1 -14.64  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000 ...
AO_LDOS:   1 1 -14.60  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000 ...
AO_LDOS:   1 1 -14.56  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000 ...
AO_LDOS:   1 1 -14.52  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000 ...
AO_LDOS:   1 1 -14.48  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000 ...
AO_LDOS:   1 1 -14.44  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000 ...
AO_LDOS:   1 1 -14.40  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000 ...
AO_LDOS:   1 1 -14.36  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000 ...
AO_LDOS:   1 1 -14.32  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000 ...
AO_LDOS:   1 1 -14.28  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000 ...
AO_LDOS:   1 1 -14.24  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000 ...
AO_LDOS:   1 1 -14.20  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000 ...
AO_LDOS:   1 1 -14.16  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000 ...
AO_LDOS:   1 1 -14.12  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000 ...
AO_LDOS:   1 1 -14.08  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000 ...
AO_LDOS:   1 1 -14.04  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000 ...
AO_LDOS:   1 1 -14.00  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000 ...
...
aoldos_Ni.png
トップ   新規 一覧 単語検索 最終更新   ヘルプ   最終更新のRSS