PMTdimer距離依存性結果

From Ecal

Jump to: navigation, search

Contents

計算条件

  • 15Ang^3 cube,PBE,しかし、大きすぎて計算が大変で失敗だった。(kinoの環境では)pwemax=6がメモリーが足りなくて計算ができない。途中で10Ang^3へ変更。
  • 目的:R,EHの設定、local orbital,valence orbitalの設定。

訂正&注意:kinoのenergy単位はauではなくてRy。混在すると間違うのでその単位はauのままにする。読みかえること。

以下<EH,EH2>,<EH> sectionと区別がないものは<EH,EH2>と2つ使った結果。


R*  R        / core(spillout)/ valence
10 R=  2.299 / 1s(0.000038)  /2s 2p 3d 4f 5g 

という意味。local orbialを用いている場合はvalenceで2つの同じLの軌道が表示される。coreがspilloutしていない場合はその軌道は表示されない。 coreが2つ表示される場合はspin splittingしているため。

収束性

MIX...のパラメタ依存性。

Ca2/+lo/EH,EH2

File:Ca2-mol.lo.mix_dep.png

Ni2/EH,EH2

File:Ni2-mol.mix_dep.png

Ni2/EH

File:Ni2-mol.s3.eh1.mix_dep.png

Cr2

File:Cr2-mol.mix_dep.png

Mn2

File:Mn2-mol.fsmom2.mix_dep.png

V2

File:Mn2-mol.lo.fsmom2.mix_dep.png

Ti2

File:Ti2-mol.lo.s3.mix_dep.png

Ge2

File:Ge2-mol.s3.mix_dep.png

Fe2

File:Fe2-mol.s7.mix_dep.png

  • spin momentは~6.8で7(integer)になりません。

H2

EH,EH2

File:H2-mol.a15.a_dep.png

EH

He2

EH,EH2

EH

Li2

05 R=  1.292 / /1s 2s 2p 3d 4f 5g 
06 R=  1.550 / /1s 2s 2p 3d 4f 5g 
07 R=  1.809 / /1s 2s 2p 3d 4f 5g 
08 R=  2.067 / /1s 2s 2p 3d 4f 5g 
09 R=  2.325 / /1s 2s 2p 3d 4f 5g 
10 R=  2.584 / /1s 2s 2p 3d 4f 5g 

EH,EH2

File:Li2-mol.a15.a_dep.pwe.png

EH

Be2

05 R=  1.150 / 1s(0.015184) /2s 2p 3d 4f 5g 
06 R=  1.380 / 1s(0.004762) /2s 2p 3d 4f 5g 
07 R=  1.609 / 1s(0.001466) /2s 2p 3d 4f 5g 
08 R=  1.839 / 1s(0.000441) /2s 2p 3d 4f 5g 
09 R=  2.069 / 1s(0.000131) /2s 2p 3d 4f 5g 
10 R=  2.299 / 1s(0.000038) /2s 2p 3d 4f 5g 
11 R=  2.529 / 1s(0.000011) /2s 2p 3d 4f 5g 

EH,EH2

EH

B2

05 R=  0.724 / 1s(0.044798) 1s(0.044798) /2s 2p 3d 4f 5g 
06 R=  0.869 / 1s(0.017842) 1s(0.017842) /2s 2p 3d 4f 5g 
07 R=  1.013 / 1s(0.006974) 1s(0.006974) /2s 2p 3d 4f 5g 
08 R=  1.158 / 1s(0.002656) 1s(0.002656) /2s 2p 3d 4f 5g 
09 R=  1.303 / 1s(0.000996) 1s(0.000996) /2s 2p 3d 4f 5g 
10 R=  1.448 / 1s(0.000369) 1s(0.000369) /2s 2p 3d 4f 5g 
11 R=  1.593 / 1s(0.000135) 1s(0.000135) /2s 2p 3d 4f 5g 
  • ↓S=4

File:B2-mol.a_dep.pwe.png

not fully converged.

B2/+lo

05 R=  0.724 / /1s 2s 2p 3d 4f 5g 
06 R=  0.869 / /1s 2s 2p 3d 4f 5g 
07 R=  1.013 / /1s 2s 2p 3d 4f 5g 
08 R=  1.158 / /1s 2s 2p 3d 4f 5g 
09 R=  1.303 / /1s 2s 2p 3d 4f 5g 
10 R=  1.448 / /1s 2s 2p 3d 4f 5g 
11 R=  1.593 / /1s 2s 2p 3d 4f 5g 

EH,EH2

  • S3

File:B2.lo.s3.a_dep.pwe.png

  • S3, 15Ang^3

File:B2.lo.s3.a15.a_dep.pwe.png

  • S5

File:B2.lo.s5.a_dep.pwe.png

s3とs5を重ねて書くとGaussian09と同じくs5の方が基底状態になる。図の色はGaussianの図の色と同じものを使っています。

File:B2.lo.s3s5.a_dep.pwe.png

この結果はMMOM(lmfaで初期値の指定に用いる。)に依らない。(PZ=1.9を使っておりcoreはないので当然。)


C2

05 R=  0.622 / 1s(0.036052) 1s(0.035742) /2s 2p 3d 4f 5g 
06 R=  0.747 / 1s(0.013511) 1s(0.013340) /2s 2p 3d 4f 5g 
07 R=  0.871 / 1s(0.004960) 1s(0.004882) /2s 2p 3d 4f 5g 
08 R=  0.996 / 1s(0.001768) 1s(0.001736) /2s 2p 3d 4f 5g 
09 R=  1.120 / 1s(0.000625) 1s(0.000613) /2s 2p 3d 4f 5g 
10 R=  1.245 / 1s(0.000217) 1s(0.000212) /2s 2p 3d 4f 5g 
11 R=  1.369 / 1s(0.000075) 1s(0.000073) /2s 2p 3d 4f 5g 

EH,EH2

File:C2-mol.a15.m2.a_dep.png

MMOM

MMOM(atomのpolarizatin)による依存はmarginalだが0ではない。

  • S1

  • S3

EH

N2

05 R=  0.523 / 1s(0.037145) /2s 2p 3d 4f 5g 
06 R=  0.628 / 1s(0.013958) /2s 2p 3d 4f 5g 
07 R=  0.732 / 1s(0.005141) /2s 2p 3d 4f 5g 
08 R=  0.837 / 1s(0.001835) /2s 2p 3d 4f 5g 
09 R=  0.942 / 1s(0.000644) /2s 2p 3d 4f 5g 
10 R=  1.046 / 1s(0.000226) /2s 2p 3d 4f 5g 
11 R=  1.151 / 1s(0.000078) /2s 2p 3d 4f 5g 

EH,EH2

EH

O2

05 R=  0.576 / 1s(0.009890) 1s(0.009788) /2s 2p 3d 4f 5g 
06 R=  0.691 / 1s(0.002716) 1s(0.002682) /2s 2p 3d 4f 5g 
07 R=  0.807 / 1s(0.000715) 1s(0.000705) /2s 2p 3d 4f 5g 
08 R=  0.922 / 1s(0.000187) 1s(0.000184) /2s 2p 3d 4f 5g 
09 R=  1.037 / 1s(0.000048) 1s(0.000047) /2s 2p 3d 4f 5g 
10 R=  1.152 / 1s(0.000012) 1s(0.000012) /2s 2p 3d 4f 5g 
11 R=  1.268 / 1s(0.000003) 1s(0.000003) /2s 2p 3d 4f 5g 

(確認のためKinoの条件を同じにした結果を載せる。)

EH,EH2

File:O2-mol.a15.a_dep.pwe.png

EH

F2

05 R=  0.676 / 1s(0.001131) 1s(0.001131) /2s 2p 3d 4f 5g 
06 R=  0.811 / 1s(0.000189) 1s(0.000189) /2s 2p 3d 4f 5g 
07 R=  0.946 / 1s(0.000031) 1s(0.000031) /2s 2p 3d 4f 5g 
08 R=  1.081 / 1s(0.000005) 1s(0.000005) /2s 2p 3d 4f 5g 
09 R=  1.216 / 1s(0.000001) 1s(0.000001) /2s 2p 3d 4f 5g 
10 R=  1.351 / /2s 2p 3d 4f 5g 
11 R=  1.486 / /2s 2p 3d 4f 5g 

EH,EH2

File:F2-mol.15Ang.a_dep.png

EH

Ne2/VWN

05 R=  1.205 / /2s 2p 3d 4f 5g 
06 R=  1.446 / /2s 2p 3d 4f 5g 
07 R=  1.687 / /2s 2p 3d 4f 5g 
08 R=  1.928 / /2s 2p 3d 4f 5g 
09 R=  2.169 / /2s 2p 3d 4f 5g 
10 R=  2.410 / /2s 2p 3d 4f 5g 
11 R=  2.651 / /2s 2p 3d 4f 5g 

Ne2/+lo

SPEC ATOM=Ne  P=3.5 PZ=12.8
05 R=  1.205 / /2s 3s 2p 3d 4f 5g 
06 R=  1.446 / /2s 3s 2p 3d 4f 5g 
07 R=  1.687 / /2s 3s 2p 3d 4f 5g 
08 R=  1.928 / /2s 3s 2p 3d 4f 5g 
09 R=  2.169 / /2s 3s 2p 3d 4f 5g 
10 R=  2.410 / /2s 3s 2p 3d 4f 5g 
11 R=  2.651 / /2s 3s 2p 3d 4f 5g 

EH,EH2

File:Ne2-mol.lo.lsda.a_dep.pwe.png

  • 15Ang^3

File:Ne2-mol.lo.a15.a_dep.pwe.png

EH

File:Ne2-mol.lo.lsda.eh1.a_dep.pwe.png

Na2

05 R=  1.460 / 2s(0.049518) /3s 2p 3p 3d 4f 5g 
06 R=  1.752 / 2s(0.017252) /3s 2p 3p 3d 4f 5g 
07 R=  2.044 / 2s(0.005817) /3s 2p 3p 3d 4f 5g 
08 R=  2.336 / 2s(0.001919) /3s 2p 3p 3d 4f 5g 
09 R=  2.628 / 2s(0.000626) /3s 2p 3p 3d 4f 5g 
10 R=  2.920 / 2s(0.000202) /3s 2p 3p 3d 4f 5g 
11 R=  3.000 / 2s(0.000148) /3s 2p 3p 3d 4f 5g 

EH,EH2

EH

Mg2

05 R=  1.655 / 2s(0.010145) /3s 2p 3p 3d 4f 5g 
06 R=  1.986 / 2s(0.002350) /3s 2p 3p 3d 4f 5g 
07 R=  2.317 / 2s(0.000529) /3s 2p 3p 3d 4f 5g 
08 R=  2.648 / 2s(0.000117) /3s 2p 3p 3d 4f 5g 
09 R=  2.979 / 2s(0.000025) /3s 2p 3p 3d 4f 5g 
10 R=  3.000 / 2s(0.000023) /3s 2p 3p 3d 4f 5g 
11 R=  3.000 / 2s(0.000023) /3s 2p 3p 3d 4f 5g 

EH,EH2

EH

Al2

05 R=  1.178 / 2s(0.041916) 2p(0.054765) 2s(0.041916) 2p(0.054724) /3s 3p 3d 4f 5g 
06 R=  1.413 / 2s(0.013310) 2p(0.021515) 2s(0.013296) 2p(0.021451) /3s 3p 3d 4f 5g 
07 R=  1.649 / 2s(0.004054) 2p(0.008231) 2s(0.004033) 2p(0.008160) /3s 3p 3d 4f 5g 
08 R=  1.884 / 2s(0.001210) 2p(0.003109) 2s(0.001196) 2p(0.003060) /3s 3p 3d 4f 5g 
09 R=  2.120 / 2s(0.000352) 2p(0.001153) 2s(0.000346) 2p(0.001125) /3s 3p 3d 4f 5g 
10 R=  2.355 / 2s(0.000101) 2p(0.000424) 2s(0.000099) 2p(0.000411) /3s 3p 3d 4f 5g 
11 R=  2.591 / 2s(0.000029) 2p(0.000154) 2s(0.000028) 2p(0.000148) /3s 3p 3d 4f 5g 

EH,EH2

EH

Si2

05 R=  1.029 / 2s(0.049941) 2p(0.057413) 2s(0.049911) 2p(0.057338) /3s 3p 3d 4f 5g 
06 R=  1.234 / 2s(0.016327) 2p(0.022368) 2s(0.016281) 2p(0.022269) /3s 3p 3d 4f 5g 
07 R=  1.440 / 2s(0.005090) 2p(0.008434) 2s(0.005052) 2p(0.008350) /3s 3p 3d 4f 5g 
08 R=  1.646 / 2s(0.001539) 2p(0.003111) 2s(0.001520) 2p(0.003062) /3s 3p 3d 4f 5g 
09 R=  1.851 / 2s(0.000457) 2p(0.001133) 2s(0.000449) 2p(0.001109) /3s 3p 3d 4f 5g 
10 R=  2.057 / 2s(0.000133) 2p(0.000405) 2s(0.000130) 2p(0.000394) /3s 3p 3d 4f 5g 
11 R=  2.263 / 2s(0.000038) 2p(0.000143) 2s(0.000037) 2p(0.000139) /3s 3p 3d 4f 5g

EH,EH2

File:Si2-mol.15Ang.a_dep.pwe.png

EH

P2

05 R=  0.904 / 2s(0.060626) 2p(0.062992) /3s 3p 3d 4f 5g 
06 R=  1.084 / 2s(0.020567) 2p(0.024698) /3s 3p 3d 4f 5g 
07 R=  1.265 / 2s(0.006613) 2p(0.009317) /3s 3p 3d 4f 5g 
08 R=  1.446 / 2s(0.002055) 2p(0.003428) /3s 3p 3d 4f 5g 
09 R=  1.627 / 2s(0.000622) 2p(0.001237) /3s 3p 3d 4f 5g 
10 R=  1.807 / 2s(0.000186) 2p(0.000442) /3s 3p 3d 4f 5g 
11 R=  1.988 / 2s(0.000055) 2p(0.000155) /3s 3p 3d 4f 5g 

EH,EH2

EH

S2

05 R=  0.915 / 2s(0.034919) 2p(0.036496) 2s(0.034919) 2p(0.036496) /3s 3p 3d 4f 5g 
06 R=  1.098 / 2s(0.010083) 2p(0.012265) 2s(0.010083) 2p(0.012265) /3s 3p 3d 4f 5g 
07 R=  1.281 / 2s(0.002768) 2p(0.003971) 2s(0.002768) 2p(0.003971) /3s 3p 3d 4f 5g 
08 R=  1.463 / 2s(0.000739) 2p(0.001259) 2s(0.000739) 2p(0.001259) /3s 3p 3d 4f 5g 
09 R=  1.646 / 2s(0.000191) 2p(0.000389) 2s(0.000191) 2p(0.000389) /3s 3p 3d 4f 5g 
10 R=  1.829 / 2s(0.000048) 2p(0.000118) 2s(0.000048) 2p(0.000118) /3s 3p 3d 4f 5g 
11 R=  2.012 / 2s(0.000012) 2p(0.000036) 2s(0.000012) 2p(0.000036) /3s 3p 3d 4f 5g 

EH,EH2

EH

Cl2

05 R=  0.970 / 2s(0.013968) 2p(0.015258) /3s 3p 3d 4f 5g 
06 R=  1.164 / 2s(0.003153) 2p(0.004081) /3s 3p 3d 4f 5g 
07 R=  1.358 / 2s(0.000676) 2p(0.001049) /3s 3p 3d 4f 5g 
08 R=  1.552 / 2s(0.000140) 2p(0.000263) /3s 3p 3d 4f 5g 
09 R=  1.746 / 2s(0.000028) 2p(0.000065) /3s 3p 3d 4f 5g 
10 R=  1.940 / 2s(0.000006) 2p(0.000016) /3s 3p 3d 4f 5g 
11 R=  2.134 / 2s(0.000001) 2p(0.000004) /3s 3p 3d 4f 5g 

EH,EH2

File:Cl2-mol.15Ang.a_dep.pwe.png

EH

Ar2

05 R=  1.625 / 2s(0.000027) 2p(0.000054) /3s 3p 3d 4f 5g 
06 R=  1.950 / 2s(0.000001) 2p(0.000004) /3s 3p 3d 4f 5g 
07 R=  2.275 / /3s 3p 3d 4f 5g 
08 R=  2.600 / /3s 3p 3d 4f 5g 
09 R=  2.925 / /3s 3p 3d 4f 5g 
10 R=  3.000 / /3s 3p 3d 4f 5g 
11 R=  3.000 / /3s 3p 3d 4f 5g 

EH,EH2

15Ang^3

EH

K2/4s 3p 4p

05 R=  1.890 / 2s(0.000001) 2p(0.000002) 3s(0.106241) /4s 3p 4p 3d 4f 5g 
06 R=  2.268 / 3s(0.038395) /4s 3p 4p 3d 4f 5g 
07 R=  2.646 / 3s(0.013163) /4s 3p 4p 3d 4f 5g 
08 R=  3.000 / 3s(0.004686) /4s 3p 4p 3d 4f 5g 
09 R=  3.000 / 3s(0.004686) /4s 3p 4p 3d 4f 5g 
10 R=  3.000 / 3s(0.004686) /4s 3p 4p 3d 4f 5g 
11 R=  3.000 / 3s(0.004686) /4s 3p 4p 3d 4f 5g 

EH,EH2

File:K2-mol.a15.a_dep.pwe.png

EH

K2/4s 4p

05 R=  1.890 / 2s(0.000001) 2p(0.000002) 3s(0.106241) 3p(0.197579) /4s 4p 3d 4f 5g 
06 R=  2.268 / 3s(0.038395) 3p(0.094714) /4s 4p 3d 4f 5g 
07 R=  2.646 / 3s(0.013163) 3p(0.043690) /4s 4p 3d 4f 5g 
08 R=  3.000 / 3s(0.004686) 3p(0.020720) /4s 4p 3d 4f 5g 
09 R=  3.000 / 3s(0.004686) 3p(0.020720) /4s 4p 3d 4f 5g 
10 R=  3.000 / 3s(0.004686) 3p(0.020720) /4s 4p 3d 4f 5g 
11 R=  3.000 / 3s(0.004686) 3p(0.020720) /4s 4p 3d 4f 5g 

EH,EH2

K2,+lo

SPEC ATOM=K  PZ=13.8 13.8 P=4.2 4.2
05 R=  1.890 / 2s(0.000001) 2p(0.000002) /3s 4s 3p 4p 3d 4f 5g 
06 R=  2.268 / /3s 4s 3p 4p 3d 4f 5g 
07 R=  2.646 / /3s 4s 3p 4p 3d 4f 5g 
08 R=  3.000 / /3s 4s 3p 4p 3d 4f 5g 
09 R=  3.000 / /3s 4s 3p 4p 3d 4f 5g 
10 R=  3.000 / /3s 4s 3p 4p 3d 4f 5g 

EH,EH2

EH

K2,lo(cMTO),kotani

PZを13.8などで入れる(PZの13.9などの少数点以下の数は初期条件にのみ関係。アルカリ金属はセミコアが浅く小数点以下を.8(Kのとき13.8)にしないと計算が落ちる時がある)。計算はすべてpwe2

  • 以下の15ang cellの場合と比較すると、12ang cellは不十分なことがわかる。なにしろ柔らかい結合である上に、原子間距離が4angと大きい。

  • PZとして3sの影響をみた。constant shiftがおこるのみ。

  • R*依存性は十分に小さい。とくにR*=0.8とR*=0.9はほとんど重なっている(確認したが正しいようだ)。R*=0.8はR=3.024 a.u.になる。

SrなどでもそうだがRを3.0a.u.以下程度にとるほうがよいようだ。

  • 15angのときと20angのときを比較;セルサイズを大きくすると原子間距離が、わずかに小さくなるように見える。正確な原子間距離決定には、20angでもまだ不十分かもしれないがコード整備が先だろう。

Ca2

05 R=  1.964 / 3s(0.049685) /4s 3p 4p 3d 4f 5g 
06 R=  2.357 / 3s(0.014120) /4s 3p 4p 3d 4f 5g 
07 R=  2.749 / 3s(0.003836) /4s 3p 4p 3d 4f 5g 
08 R=  3.000 / 3s(0.001640) /4s 3p 4p 3d 4f 5g 
09 R=  3.000 / 3s(0.001640) /4s 3p 4p 3d 4f 5g 
10 R=  3.000 / 3s(0.001640) /4s 3p 4p 3d 4f 5g 
11 R=  3.000 / 3s(0.001640) /4s 3p 4p 3d 4f 5g 

エネルギースケールが厳しい小さい。Chemical Physics Letters 416 (2005) 370–375の値だと、8.09a.u.×0.96=4.11ang.

File:Ca2-mol.a_dep.png

Ca2/nrel

  • REL=F

Ca2/+lo

SPEC ATOM=Ca Z=20  PZ=13.8,13.8 P=4.2,4.2
05 R=  1.964 / /3s 4s 3p 4p 3d 4f 5g 
06 R=  2.357 / /3s 4s 3p 4p 3d 4f 5g 
07 R=  2.749 / /3s 4s 3p 4p 3d 4f 5g 
08 R=  3.000 / /3s 4s 3p 4p 3d 4f 5g 
09 R=  3.000 / /3s 4s 3p 4p 3d 4f 5g 
10 R=  3.000 / /3s 4s 3p 4p 3d 4f 5g 
11 R=  3.000 / /3s 4s 3p 4p 3d 4f 5g 

File:Ca2-mol.lo.a_dep.png

File:Ca2-mol.lo.a15.a_dep.new.png

Sc2

05 R=  1.240 / 2s(0.000066) 2p(0.000094) 3s(0.300422) 3p(0.384089) 2s(0.000066) 2p(0.000093) 3s(0.299905) 3p(0.383767) /4s 4p 3d 4f 5g 
06 R=  1.489 / 2s(0.000004) 2p(0.000007) 3s(0.146347) 3p(0.219575) 2s(0.000004) 2p(0.000007) 3s(0.146090) 3p(0.219573) /4s 4p 3d 4f 5g 
07 R=  1.737 / 2p(0.000001) 3s(0.066315) 3p(0.118608) 2p(0.000001) 3s(0.066278) 3p(0.118873) /4s 4p 3d 4f 5g 
08 R=  1.985 / 3s(0.028626) 3p(0.061689) 3s(0.028674) 3p(0.062039) /4s 4p 3d 4f 5g 
09 R=  2.233 / 3s(0.011967) 3p(0.031305) 3s(0.012020) 3p(0.031612) /4s 4p 3d 4f 5g 
10 R=  2.481 / 3s(0.004897) 3p(0.015629) 3s(0.004934) 3p(0.015860) /4s 4p 3d 4f 5g 
11 R=  2.729 / 3s(0.001973) 3p(0.007710) 3s(0.001997) 3p(0.007873) /4s 4p 3d 4f 5g 

EH,EH2

EH

kotani's result

pwemax=2,R=2,SYMGRP e。W=0.0001にWを小さくする必要あり(W=0.001だと占有数SPLWTSの欄が0と1に分離しない場合がある)。 local orbitalを入れないといけない(PZ=0,3.9をつかう)。

FSMOM=2のとき奇妙なラインを描く。実際、格子定数の小さいところから大きくしていくと、どんどんLUMOーHOMOギャップが小さくなっていき、 Distance=2.36より大きいところでは、逆転してしまうようだ(いずれLUMOーHOMOギャップの詳細もだす)。 FSMOM=2の格子定数が大きいところでは収束に時間がかかる(初期条件にrst,mixで小さい格子定数での収束解からスタートするようなスクリプトである のは収束に回数がかかる一因。初期条件がよくない)。ただ、それだけでなく、gapが小さいときの不安定性もある。

  • 以下のラインは、FSMOM=0(cubicで格子の対称性を仮定したもの)、FSMOM=0,2,4(対称性をこわしたもの。初期条件に対称性の壊れたものを使わないと対称のまま収束する場合がある)。PBE-GGA(Gaussian)の結果と比較すると、FSMOM=0(Sz=0)の格子定数は少しずれてるように見える。

File:Sc2_PZ_W1e-4_SYMe.png

Ti2

05 R=  0.896 / 2s(0.001571) 2p(0.001625) 3s(0.586717) 3p(0.639307) 2s(0.001571) 2p(0.001625) 3s(0.586717) 3p(0.639307) /4s 4p 3d 4f 5g 
06 R=  1.076 / 2s(0.000196) 2p(0.000239) 3s(0.380720) 3p(0.453903) 2s(0.000196) 2p(0.000239) 3s(0.380720) 3p(0.453903) /4s 4p 3d 4f 5g 
07 R=  1.255 / 2s(0.000023) 2p(0.000034) 3s(0.227533) 3p(0.302420) 2s(0.000023) 2p(0.000034) 3s(0.227533) 3p(0.302420) /4s 4p 3d 4f 5g 
08 R=  1.434 / 2s(0.000003) 2p(0.000005) 3s(0.128297) 3p(0.192555) 2s(0.000003) 2p(0.000005) 3s(0.128297) 3p(0.192555) /4s 4p 3d 4f 5g 
09 R=  1.614 / 2p(0.000001) 3s(0.069210) 3p(0.118437) 2p(0.000001) 3s(0.069210) 3p(0.118437) /4s 4p 3d 4f 5g 
10 R=  1.793 / 3s(0.036388) 3p(0.071410) 3s(0.036388) 3p(0.071410) /4s 4p 3d 4f 5g 
11 R=  1.972 / 3s(0.018731) 3p(0.042362) 3s(0.018731) 3p(0.042362) /4s 4p 3d 4f 5g 

File:Ti2-mol.a_dep.pwe.png

+lo

SPEC ATOM=Ti Z=22 PZ=0,13.8,P=0,4.2 
05 R=  0.896 / 2s(0.001571) 2p(0.001625) 3s(0.586717) 2s(0.001571) 2p(0.001625) 3s(0.586717) /4s 3p 4p 3d 4f 5g 
06 R=  1.076 / 2s(0.000196) 2p(0.000239) 3s(0.380720) 2s(0.000196) 2p(0.000239) 3s(0.380720) /4s 3p 4p 3d 4f 5g 
07 R=  1.255 / 2s(0.000023) 2p(0.000034) 3s(0.227533) 2s(0.000023) 2p(0.000034) 3s(0.227533) /4s 3p 4p 3d 4f 5g 
08 R=  1.434 / 2s(0.000003) 2p(0.000005) 3s(0.128297) 2s(0.000003) 2p(0.000005) 3s(0.128297) /4s 3p 4p 3d 4f 5g 
09 R=  1.614 / 2p(0.000001) 3s(0.069210) 2p(0.000001) 3s(0.069210) /4s 3p 4p 3d 4f 5g 
10 R=  1.793 / 3s(0.036388) 3s(0.036388) /4s 3p 4p 3d 4f 5g 
11 R=  1.972 / 3s(0.018731) 3s(0.018731) /4s 3p 4p 3d 4f 5g 

kotani's resutls

MTO半径R=1.6とした。PZ=0,3.9は必須(13.9にすること)。FSMOM=8は格子定数が大きい領域で振動して収束しなかった。 異方性のある解に自動的に収束させるため、わざと格子を歪ませる(PLAT)。それでないと縮退がきれいにとけないということが起こり得る。W=0.0001(smearing)は小さい。この方が一般的にはよいようだが、local minimumから抜け出せないというようなことも考えうる。

File:Ti2kotani.png

V2

05 R=  0.822 / 2s(0.002140) 2p(0.002080) 3s(0.615996) 3p(0.658998) 2s(0.002140) 2p(0.002080) 3s(0.615996) 3p(0.658998) /4s 4p 3d 4f 5g 
06 R=  0.987 / 2s(0.000285) 2p(0.000323) 3s(0.410181) 3p(0.475197) 2s(0.000285) 2p(0.000323) 3s(0.410181) 3p(0.475197) /4s 4p 3d 4f 5g 
07 R=  1.151 / 2s(0.000037) 2p(0.000048) 3s(0.251477) 3p(0.321515) 2s(0.000037) 2p(0.000048) 3s(0.251477) 3p(0.321515) /4s 4p 3d 4f 5g 
08 R=  1.315 / 2s(0.000005) 2p(0.000007) 3s(0.145377) 3p(0.207808) 2s(0.000005) 2p(0.000007) 3s(0.145377) 3p(0.207808) /4s 4p 3d 4f 5g 
09 R=  1.480 / 2s(0.000001) 2p(0.000001) 3s(0.080340) 3p(0.129673) 2s(0.000001) 2p(0.000001) 3s(0.080340) 3p(0.129673) /4s 4p 3d 4f 5g 
10 R=  1.644 / 3s(0.043254) 3p(0.079296) 3s(0.043254) 3p(0.079296) /4s 4p 3d 4f 5g 

File:V2-mol.a_dep.pwe.png

+lo

SPEC ATOM=V Z=23  PZ=13.8 13.8 P=4.2 4.2
05 R=  0.822 / 2s(0.002140) 2p(0.002080) 2s(0.002140) 2p(0.002080) /3s 4s 3p 4p 3d 4f 5g 
06 R=  0.987 / 2s(0.000285) 2p(0.000323) 2s(0.000285) 2p(0.000323) /3s 4s 3p 4p 3d 4f 5g 
07 R=  1.151 / 2s(0.000037) 2p(0.000048) 2s(0.000037) 2p(0.000048) /3s 4s 3p 4p 3d 4f 5g 
08 R=  1.315 / 2s(0.000005) 2p(0.000007) 2s(0.000005) 2p(0.000007) /3s 4s 3p 4p 3d 4f 5g 
09 R=  1.480 / 2s(0.000001) 2p(0.000001) 2s(0.000001) 2p(0.000001) /3s 4s 3p 4p 3d 4f 5g 
10 R=  1.644 / /3s 4s 3p 4p 3d 4f 5g 

Cr2

07 R=  1.055 / 2s(0.000060) 2p(0.000072) 3s(0.286390) 3p(0.355106) /4s 4p 3d 4f 5g 
08 R=  1.206 / 2s(0.000008) 2p(0.000011) 3s(0.171437) 3p(0.236053) /4s 4p 3d 4f 5g 
09 R=  1.357 / 2s(0.000001) 2p(0.000002) 3s(0.098510) 3p(0.152081) /4s 4p 3d 4f 5g 
10 R=  1.507 / 3s(0.055189) 3p(0.096121) /4s 4p 3d 4f 5g 
11 R=  1.658 / 3s(0.030147) 3p(0.059596) /4s 4p 3d 4f 5g 

File:Cr2-mol.a_dep.pwe.png

Cr2/+lo

SPEC ATOM=Cr  PZ=0,13.8 P=0,4.2 MMOM=0 0 0 0  
05 R=  0.754 / 2s(0.002972) 2p(0.002731) 3s(0.652438) /4s 3p 4p 3d 4f 5g 
06 R=  0.904 / 2s(0.000435) 2p(0.000457) 3s(0.451316) /4s 3p 4p 3d 4f 5g 
07 R=  1.055 / 2s(0.000060) 2p(0.000072) 3s(0.286390) /4s 3p 4p 3d 4f 5g 
08 R=  1.206 / 2s(0.000008) 2p(0.000011) 3s(0.171437) /4s 3p 4p 3d 4f 5g 
09 R=  1.357 / 2s(0.000001) 2p(0.000002) 3s(0.098510) /4s 3p 4p 3d 4f 5g 
10 R=  1.507 / 3s(0.055189) /4s 3p 4p 3d 4f 5g 
11 R=  1.658 / 3s(0.030147) /4s 3p 4p 3d 4f 5g 

EH,EH2

EH

Mn2

07 R=  1.090 / 2s(0.000018) 2p(0.000022) 3s(0.204970) 3p(0.265148) 2s(0.000018) 2p(0.000022) 3s(0.204970) 3p(0.265148) /4s 4p 3d 4f 5g 
08 R=  1.246 / 2s(0.000002) 2p(0.000003) 3s(0.112508) 3p(0.163343) 2s(0.000002) 2p(0.000003) 3s(0.112508) 3p(0.163343) /4s 4p 3d 4f 5g 
09 R=  1.402 / 3s(0.059362) 3p(0.097594) 3s(0.059362) 3p(0.097594) /4s 4p 3d 4f 5g 
10 R=  1.557 / 3s(0.030586) 3p(0.057265) 3s(0.030586) 3p(0.057265) /4s 4p 3d 4f 5g 

EH,EH2

Mn2/lo+sp

07 R=  1.090 / 2s(0.000018) 2p(0.000023) 2s(0.000017) 2p(0.000022) /3s 4s 3p 4p 3d 4f 5g 
08 R=  1.246 / 2s(0.000002) 2p(0.000003) 2s(0.000002) 2p(0.000003) /3s 4s 3p 4p 3d 4f 5g 
09 R=  1.402 / /3s 4s 3p 4p 3d 4f 5g 
10 R=  1.557 / /3s 4s 3p 4p 3d 4f 5g 
11 R=  1.713 / /3s 4s 3p 4p 3d 4f 5g 

EH,EH2

comment:EHの値によってカーブが違う。

EH,EH2/S11

やりなおし。

Mn2/lo+p

SPEC ATOM=Mn  PZ=0 13.8 

File:result.Mn2-mol.lo+p.s11.a_dep.pwe.png

Mn2/S3/lo+sp

SPEC ATOM=Mn MMOM=0 0 3 0  PZ=3.9 3.9

MMOM

このスケールではmarginal。

Mn2/S11/lo+sp

SPEC ATOM=Mn MMOM=0 0 0 0  PZ=3.9 3.9

MMOM

このスケールでは差が見える。

この場合E(MMOM=0)>E(MMOM=3)>E(MMOM=5)になる。

Fe2

05 R=  0.950 / 2s(0.000065) 2p(0.000073) 3s(0.284777) 3p(0.342595) 2s(0.000065) 2p(0.000073) 3s(0.284777) 3p(0.342595) /4s 4p 3d 4f 5g 
06 R=  1.140 / 2s(0.000004) 2p(0.000005) 3s(0.136334) 3p(0.188144) 2s(0.000004) 2p(0.000005) 3s(0.136334) 3p(0.188144) /4s 4p 3d 4f 5g 
07 R=  1.330 / 3s(0.060707) 3p(0.097702) 3s(0.060707) 3p(0.097702) /4s 4p 3d 4f 5g 
08 R=  1.521 / 3s(0.025687) 3p(0.048782) 3s(0.025687) 3p(0.048782) /4s 4p 3d 4f 5g 
09 R=  1.711 / 3s(0.010585) 3p(0.023879) 3s(0.010585) 3p(0.023879) /4s 4p 3d 4f 5g 
10 R=  1.901 / 3s(0.004270) 3p(0.011504) 3s(0.004270) 3p(0.011504) /4s 4p 3d 4f 5g 
11 R=  2.091 / 3s(0.001696) 3p(0.005480) 3s(0.001696) 3p(0.005480) /4s 4p 3d 4f 5g 

EH,EH2

Kotani's Fe2

R*=0.9 , fsmom=2Sz,初期条件はCubicでないセルで作っている。fsmom=4(Sz=5)でdistance=2.1を越えると別の解へ移行するように見える。

File:Fe2Cubickotani1.txt

Fe/lo+3p

ATOM=Fe PZ=0 13.8 P=0 4.3

Fe/lo+3sp

ATOM=Fe PZ=13.8 13.8 P=4.3 4.3

Co2

06 R=  1.125 / 2s(0.000002) 2p(0.000003) 3s(0.116279) 3p(0.163331) 2s(0.000002) 2p(0.000003) 3s(0.116279) 3p(0.163331) /4s 4p 3d 4f 5g 
07 R=  1.312 / 3s(0.049937) 3p(0.082173) 3s(0.049937) 3p(0.082173) /4s 4p 3d 4f 5g 
08 R=  1.500 / 3s(0.020414) 3p(0.039812) 3s(0.020414) 3p(0.039812) /4s 4p 3d 4f 5g 
09 R=  1.687 / 3s(0.008141) 3p(0.018938) 3s(0.008141) 3p(0.018938) /4s 4p 3d 4f 5g 
10 R=  1.875 / 3s(0.003165) 3p(0.008838) 3s(0.003165) 3p(0.008838) /4s 4p 3d 4f 5g 

EH,EH2

EH

Co2/lo+3p

ATOM=Co PZ=0 13.8 P=0 4.3

Co2/lo+3sp

ATOM=Co PZ=13.8 13.8 P=4.3 4.3

Ni2

05 R=  1.005 / 2s(0.000007) 2p(0.000009) 3s(0.160012) 3p(0.209830) 2s(0.000007) 2p(0.000009) 3s(0.160012) 3p(0.209830) /4s 4p 3d 4f 5g 
06 R=  1.206 / 3s(0.063406) 3p(0.098507) 3s(0.063406) 3p(0.098507) /4s 4p 3d 4f 5g 
07 R=  1.407 / 3s(0.023590) 3p(0.044082) 3s(0.023590) 3p(0.044082) /4s 4p 3d 4f 5g 
08 R=  1.607 / 3s(0.008481) 3p(0.019234) 3s(0.008481) 3p(0.019234) /4s 4p 3d 4f 5g 
09 R=  1.808 / 3s(0.002956) 3p(0.008200) 3s(0.002956) 3p(0.008200) /4s 4p 3d 4f 5g 
10 R=  2.009 / 3s(0.001012) 3p(0.003451) 3s(0.001012) 3p(0.003451) /4s 4p 3d 4f 5g 

EH,EH2

EH,EH2, NREL

  • REL=F

REL=Fの効果はr_eが大きくなるだけらしい。

EH

Ni2,kotani

収束は、

  1. 初期条件をCubicでないセルで作る。W=0.01だとより安定して収束する。
  2. Cubicなセルでsymops iでそれを読み込んで収束計算。

の手順で行いうる。

fsmom=0,nspin=2では解が得られるがすこし強引なconstranitをかけたことになっていると思われる。

木野氏によるとGaussianでparaで収束させるとE(LUMO)<E(HOMO)となるとのこと(占有軌道に穴があく)。


以下は記録のためのこす(要修正)

収束させるのは厄介な点がある。W=0.01で異方的PLATで収束させれた。PZなし。おのおのatmファイルから出発。 考えられる危険性は、初期条件(これはまるい原子の重ね合わせ)を保持したまま収束してしまうこと。

注:W=0.01の計算では、オンサイトモーメントはなし。


次にこれで得たrstを読み込んでcubicなPLATで収束させた。W=0.001

fsmom=0でpwe4のとき二本のラインがある。青の□は、para(upとdownが一致している)。 この計算はたまたままちがってmixmを消さない計算で見出された。

それ以外のfsmom=0の計算(pwe2,pwe3,pwe4の*)は「paraではないがfsmom=0」という状態に収束している。 オンサイトでわずかにモーメント(0.01程度)をもつ。すなわち、ferroになりたいんだが、fixed moment法のため無理にモーメントを固定されているというような状態のようだ.

Cu2

05 R=  1.063 / 2s(0.000001) 2p(0.000002) 3s(0.102802) 3p(0.146556) /4s 4p 3d 4f 5g 
06 R=  1.276 / 3s(0.036004) 3p(0.062620) /4s 4p 3d 4f 5g 
07 R=  1.488 / 3s(0.012005) 3p(0.025822) /4s 4p 3d 4f 5g 
08 R=  1.701 / 3s(0.003851) 3p(0.010355) /4s 4p 3d 4f 5g 
09 R=  1.913 / 3s(0.001215) 3p(0.004108) /4s 4p 3d 4f 5g 
10 R=  2.126 / 3s(0.000376) 3p(0.001606) /4s 4p 3d 4f 5g 
11 R=  2.339 / 3s(0.000115) 3p(0.000623) /4s 4p 3d 4f 5g 

EH,EH2

File:Cu2-mol.15Ang.a_dep.png

relativisity off

EH

Cu2/+lo

SPEC ATOM=Cu Z=29 PZ=13.8, 13.8 P=4.2,4.2
05 R=  1.063 / 2s(0.000001) 2p(0.000002) /3s 4s 3p 4p 3d 4f 5g 
06 R=  1.276 / /3s 4s 3p 4p 3d 4f 5g 
07 R=  1.488 / /3s 4s 3p 4p 3d 4f 5g 
08 R=  1.701 / /3s 4s 3p 4p 3d 4f 5g 
09 R=  1.913 / /3s 4s 3p 4p 3d 4f 5g 
10 R=  2.126 / /3s 4s 3p 4p 3d 4f 5g 
11 R=  2.339 / /3s 4s 3p 4p 3d 4f 5g 

File:Cu2-mol.lo.a_dep.png

  • aは格子定数。a10=10Ang^3, a15=15Ang^3という意味。これからすると10Ang^3で収束している。
  • 核間距離依存性はEH,EH2には寄らない。
  • R*=08はpwemaxの依存性が大きい。(収束していない。)

Cu2 some cases with lo, kotani

以下の結果はPZは3.9などになっており、cMTOではない。参考として残す。

lmfaの結果, R*=10

core:
   1s    -650.39700         0.000       0.034       0.069     0.000000
   2s     -78.06623         0.070       0.197       0.307     0.000000
   2p     -67.36627         0.000       0.158       0.334     0.000000
   3s      -8.42846         0.288       0.613       0.893     0.000376
   3p      -5.30026         0.260       0.620       1.075     0.001606
 valence:      eval       node at      max at       c.t.p.   rho(r>rmt)
   4s      -0.33879         0.888       2.272       3.587     0.723146
   4p      -0.05123         0.975       3.659       8.223     0.939685
   3d      -0.37034         0.000       0.603       3.440     0.073097
   4f       0.01741         0.000      35.996      49.616*    1.000000
   5g       0.02449         0.000      37.475      49.616*    1.000000

(1)no cMTO, 4s,4p,3dがvMTO, pwemax=2 (cMTOはloと同じ意味)

ctrl.cu2

(2)3p3dをcMTO, 4s,4p,3dがvMTO, PZ=0,3.9,3.9 P=0,4.3,4.3 を追加

(3)3pをcMTO, 4s,4p,3dがvMTO, PZ=0,3.9 を追加

(4)3dをcMTO, 4s,4p,4dがvMTO, PZ=0,0,3.9 P=0,0,4.3 を追加

結論 プロットは、Gaussianの格子定数(2.250ang)のレベルenergyを基準にした。 基底数が最大なのは(2)のケース。これは(1)とそんなには離れてない。(<-意味不明)ほぼ同じに見える。一方で、(3)ではR依存性が大きくなる;予測としては、(<-?)局在性を持たせるためのAPWが不足している可能性がある。それで、3dをcMTOとしたケース(4)をやってみた;結果R依存性は小さくなった。(Kino:R依存性が小さいと収束している、というわけではない。)

(1)

(2)

(3)

(4)

Zn2

ATOM=Zn  PZ= 0 0 13.9      P= 0 0 4.2 
05 R=  1.516 / 3s(0.006967) 3p(0.015933) /4s 4p 3d 4d 4f 5g 
06 R=  1.819 / 3s(0.001230) 3p(0.003900) /4s 4p 3d 4d 4f 5g 
07 R=  2.122 / 3s(0.000209) 3p(0.000930) /4s 4p 3d 4d 4f 5g 
08 R=  2.425 / 3s(0.000035) 3p(0.000218) /4s 4p 3d 4d 4f 5g 
09 R=  2.728 / 3s(0.000006) 3p(0.000051) /4s 4p 3d 4d 4f 5g 
10 R=  3.000 / 3s(0.000001) 3p(0.000014) /4s 4p 3d 4d 4f 5g 
11 R=  3.000 / 3s(0.000001) 3p(0.000014) /4s 4p 3d 4d 4f 5g

EH,EH2

EH

Ga2

05 R=  1.162 / 3s(0.036127) 3p(0.058931) 3s(0.036127) 3p(0.058931) /4s 4p 3d 4d 4f 5g 
06 R=  1.394 / 3s(0.009394) 3p(0.019381) 3s(0.009394) 3p(0.019381) /4s 4p 3d 4d 4f 5g 
07 R=  1.626 / 3s(0.002322) 3p(0.006135) 3s(0.002322) 3p(0.006135) /4s 4p 3d 4d 4f 5g 
08 R=  1.859 / 3s(0.000553) 3p(0.001887) 3s(0.000553) 3p(0.001887) /4s 4p 3d 4d 4f 5g 
09 R=  2.091 / 3s(0.000130) 3p(0.000574) 3s(0.000130) 3p(0.000574) /4s 4p 3d 4d 4f 5g 
10 R=  2.323 / 3s(0.000030) 3p(0.000173) 3s(0.000030) 3p(0.000173) /4s 4p 3d 4d 4f 5g 
11 R=  2.556 / 3s(0.000007) 3p(0.000051) 3s(0.000007) 3p(0.000051) /4s 4p 3d 4d 4f 5g 

EH,EH2

Ge2

05 R=  1.144 / 3s(0.029378) 3p(0.047983) 3s(0.029318) 3p(0.047830) /4s 4p 3d 4d 4f 5g 
06 R=  1.373 / 3s(0.007116) 3p(0.014690) 3s(0.007084) 3p(0.014597) /4s 4p 3d 4d 4f 5g 
07 R=  1.602 / 3s(0.001637) 3p(0.004321) 3s(0.001623) 3p(0.004273) /4s 4p 3d 4d 4f 5g 
08 R=  1.831 / 3s(0.000364) 3p(0.001239) 3s(0.000360) 3p(0.001218) /4s 4p 3d 4d 4f 5g 
09 R=  2.060 / 3s(0.000079) 3p(0.000349) 3s(0.000078) 3p(0.000341) /4s 4p 3d 4d 4f 5g 
10 R=  2.289 / 3s(0.000017) 3p(0.000097) 3s(0.000017) 3p(0.000094) /4s 4p 3d 4d 4f 5g 
11 R=  2.518 / 3s(0.000004) 3p(0.000027) 3s(0.000003) 3p(0.000026) /4s 4p 3d 4d 4f 5g 

EH,EH2

EH

As2

05 R=  1.004 / 3s(0.051294) 3p(0.074330) 3d(0.154071) /4s 4p 4d 4f 5g 
06 R=  1.205 / 3s(0.014199) 3p(0.025089) 3d(0.080122) /4s 4p 4d 4f 5g 
07 R=  1.405 / 3s(0.003727) 3p(0.008119) 3d(0.041154) /4s 4p 4d 4f 5g 
08 R=  1.606 / 3s(0.000936) 3p(0.002534) 3d(0.020871) /4s 4p 4d 4f 5g 
09 R=  1.807 / 3s(0.000229) 3p(0.000775) 3d(0.010514) /4s 4p 4d 4f 5g 
10 R=  2.008 / 3s(0.000055) 3p(0.000233) 3d(0.005271) /4s 4p 4d 4f 5g 
11 R=  2.208 / 3s(0.000013) 3p(0.000070) 3d(0.002642) /4s 4p 4d 4f 5g 

EH,EH2

File:As2-mol.a15.a_dep.pwe.png

EH

Se2

05 R=  1.042 / 3s(0.030002) 3p(0.045788) 3d(0.102028) /4s 4p 4d 4f 5g 
06 R=  1.251 / 3s(0.007081) 3p(0.013365) 3d(0.046852) /4s 4p 4d 4f 5g 
07 R=  1.459 / 3s(0.001587) 3p(0.003743) 3d(0.021199) /4s 4p 4d 4f 5g 
08 R=  1.668 / 3s(0.000340) 3p(0.001010) 3d(0.009438) /4s 4p 4d 4f 5g 
09 R=  1.876 / 3s(0.000071) 3p(0.000268) 3d(0.004182) /4s 4p 4d 4f 5g 
10 R=  2.085 / 3s(0.000015) 3p(0.000070) 3d(0.001835) /4s 4p 4d 4f 5g 
11 R=  2.293 / 3s(0.000003) 3p(0.000018) 3d(0.000805) /4s 4p 4d 4f 5g 

EH,EH2

File:Se2-mol.a15.a_dep.pwe.png

EH

Br2

05 R=  1.100 / 3s(0.014425) 3p(0.023708) 3d(0.059456) 4s(0.853132) /5s 4p 4d 4f 5g 
06 R=  1.320 / 3s(0.002794) 3p(0.005784) 3d(0.023528) 4s(0.707070) /5s 4p 4d 4f 5g 
07 R=  1.540 / 3s(0.000511) 3p(0.001345) 3d(0.009113) 4s(0.541751) /5s 4p 4d 4f 5g 
08 R=  1.760 / 3s(0.000090) 3p(0.000303) 3d(0.003480) 4s(0.390969) /5s 4p 4d 4f 5g 
09 R=  1.980 / 3s(0.000015) 3p(0.000067) 3d(0.001317) 4s(0.269895) /5s 4p 4d 4f 5g 
10 R=  2.200 / 3s(0.000003) 3p(0.000015) 3d(0.000496) 4s(0.180273) /5s 4p 4d 4f 5g 

File:Br2.mol.a_dep.png

Kotani:4sをvalenceに使わないといけない(subs/defpq.Fのデフォルト設定に問題あり。書きなおし中)。 あるいは4sをMTOに、5sにMTOcoreを使ってどちらもvalenceにする。

Br2,4s valence

SPEC ATOM=Br Z=35 P=4.5
05 R=  1.100 / 3s(0.014425) 3p(0.023708) 3d(0.059456) /4s 4p 4d 4f 5g 
06 R=  1.320 / 3s(0.002794) 3p(0.005784) 3d(0.023528) /4s 4p 4d 4f 5g 
07 R=  1.540 / 3s(0.000511) 3p(0.001345) 3d(0.009113) /4s 4p 4d 4f 5g 
08 R=  1.760 / 3s(0.000090) 3p(0.000303) 3d(0.003480) /4s 4p 4d 4f 5g 
09 R=  1.980 / 3s(0.000015) 3p(0.000067) 3d(0.001317) /4s 4p 4d 4f 5g 
10 R=  2.200 / 3s(0.000003) 3p(0.000015) 3d(0.000496) /4s 4p 4d 4f 5g 
11 R=  2.420 / 3p(0.000003) 3d(0.000186) /4s 4p 4d 4f 5g 

EH,EH2

relativisity off

EH

Kr2

05 R=  1.774 / 3s(0.000043) 3p(0.000148) 3d(0.001698) /4s 4p 4d 4f 5g 
06 R=  2.129 / 3s(0.000002) 3p(0.000011) 3d(0.000297) /4s 4p 4d 4f 5g 
07 R=  2.484 / 3p(0.000001) 3d(0.000052) /4s 4p 4d 4f 5g 
08 R=  2.839 / 3d(0.000009) /4s 4p 4d 4f 5g 
09 R=  3.000 / 3d(0.000004) /4s 4p 4d 4f 5g 
10 R=  3.000 / 3d(0.000004) /4s 4p 4d 4f 5g 
11 R=  3.000 / 3d(0.000004) /4s 4p 4d 4f 5g 

EH,EH2

relativisity off

REL=F

平衡核間距離が大きくなる。

Personal tools