* bulk Ni [#o4499a7f] ** 入力ファイル [#b46caa2a] 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) 28 0.5000 22104.64 6 1 0.1 : iatomn alfa amion iloc ivan zeta1 ... (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) 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 : npdosao ............................. (3) 1 -15.00 5.00 0.20 501 2.40000 0.200000 0.20000 14 0 0.000 : sm_n, dopping *** 説明 [#e445114b] - (2) ニッケルは磁性体であるためkspin=2を指定してそれぞれのスピンに対して収束計算を行う - (1) 6番目のパラメータでNi原子のスピン分極密度 \[ \zeta=\frac{\rho_\uparrow-\rho_\downarrow}{\rho_\uparrow+\rho_\downarrow} \] の初期値を0.1に指定する - (3) npdosao=1を指定してブラベー格子中の1つのNi原子の各軌道に射影した状態密度を計算する ** 全エネルギー [#s9a81df8] $ 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 ** 状態密度 [#a4beeae0] #ref(http://www-cp.prec.eng.osaka-u.ac.jp/puki_state/graph/dos_Ni.png,center,nolink) ** 各軌道に射影した状態密度 [#ed2d8b16] $ 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 ... ... #ref(http://www-cp.prec.eng.osaka-u.ac.jp/puki_state/graph/aoldos_Ni.png,center,nolink)