PMTbulk結果

From Ecal

Jump to: navigation, search

Contents

bulk

  • 注意事項

以下は、 core spillout により計算が破綻しているものもどう破綻するのか見せるために載せてあります。core spilloutはRが小さい時に起きます。以下R*=0.7はcore spilloutにより破綻している場合が多いのがわかります。

ここにあった「 基底の決めかたに関して(小谷による)」はdiscussionへ移動

SrTiO3

↓いくつか線があるのはkmxaを変えている。pwemaxを大きくしたときにはkmxaを大きくしないと計算が破綻することがある。ただし、kmxaはtotal energyにはほとんど影響はない。

  • Sr Ti O  : R= 3.000 2.090 1.595

↓計算できる範囲。kmxaが大きい方は計算可能。 pwemaxは2,4,6だが、表示のために+0.1ずつずらしてある。

C diamond

  • C  : R= 1.457

(小谷:がたつきは収束の甘さが原因?エネルギースケールが小さいから問題ないと思うし、これ以上にやる必要も ないようにも思う。)


R mol

  • C  : R= 0.906

Co

  • Co  : R= 2.352

change R

  • 0.7, Co  : R= 1.646
  • 0.8, Co  : R= 1.881
  • 0.9, Co  : R= 2.116
  • 1.0, Co  : R= 2.352

EH,EH=-0.1,-2.0

R mol

R=2.352では大きすぎて失敗する。

EH shift 0.1

  • 0.7, Co  : R= 1.646
  • 0.8, Co  : R= 1.881
  • 0.9, Co  : R= 2.116
  • 1.0, Co  : R= 2.352

EH shift 0.1, EH[2:4]

  • 0.7, Co  : R= 1.646
  • 0.8, Co  : R= 1.881
  • 0.9, Co  : R= 2.116
  • 1.0, Co  : R= 2.352


Cu

  • Cu  : R= 2.403

↓pwemax=2

change R

  • 0.7, Cu  : R= 1.682
  • 0.8, Cu  : R= 1.923
  • 0.9, Cu  : R= 2.163


EH,EH2=-0.1,-2.0

Fe/bcc ferro

  • Fe  : R= 2.258

EH[2:4]

change R

  • 0.7, Fe  : R= 1.581
  • 0.8, Fe  : R= 1.807
  • 0.9, Fe  : R= 2.033


Rが小さい程安定するが、total energyが上がる。

↓収束が甘いらしい。

↓Rが小さいとmin.が存在しにくくなる。

R mol

  • Fe  : R= 1.806

EH shift 0.1

  • 0.7, Fe  : R= 1.581
  • 0.8, Fe  : R= 1.807
  • 0.9, Fe  : R= 2.033
  • 1.0, Fe  : R= 2.258

EH shift 0.1, EH[2:4]

  • 0.7, Fe  : R= 1.581
  • 0.8, Fe  : R= 1.807
  • 0.9, Fe  : R= 2.033
  • 1.0, Fe  : R= 2.258

a dependence

ref. Wang et al. PRL54, 1852 (1985).

GaAs

  • Ga As  : R= 2.254 2.362


change R

core spilloutにより失敗する場合の処置

以下の計算例で、As:R=1.653と小さくした例ではギャップがなくなったまま収束しており失敗してます。Asを小さくした計算は4dのfrozen core近似がまずくなります。 Asで、sm core charge = 3.360443 (sphere) + 0.190795 (spillout) = 3.551239 で、(spillout) が大きいです。この0.190795はfrozen coreのMT外側の電荷数です。それで

ATOM= As Z= 33 R= 1.653
    RSMH=   0.827 0.827 0.827 0.827 EH=  -0.1 -0.1 -0.1 -0.1 
    RSMH2=  0.827 0.827 0.827 0.827 EH2= -2.0 -2.0 -2.0 -2.0 
    KMXA={kmxa}  LMX=3 LMXA=4
    MMOM= 0 0 0 0 
    PZ=0,0,13.9 P=0,0,4.2

と PZ=0,0,13.9 P=0,0,4.2を付け加えて,3dを対角化に入れてやると、30回まわって、

   it 30  of 30    ehf=   -8410.075100   ehk=   -8410.075122
  From last iter    ehf=   -8410.075097   ehk=   -8410.075121
  diffe(q)= -0.000003 (0.000140)    tol= 0.000001 (0.000001)   more=F

でした。収束はがたつきながらで、あと30+22回で、

c ehf=-8410.0751296 ehk=-8410.0751294 

を得ました(収束がおそすぎますが)。これはMTが大きいときの値とmRyまで一致してます。バンドギャップ0.5589eVです。 今はctrlgen.pyはGaはlocal orbtalがありAsにはlocal orbitalがない設定になっています。 *そもそもなぜPZがGaにあってAsにないか?ですがこれは、恣意的なものです。GaAsのGa4dセミコアはvalenceの直下にあり、 コアとして扱うのが不適当な場合が多いからです。--TakaoKotani 01:06, 28 December 2010 (JST)


  • 0.7, Ga As  : R= 1.578 1.653
  • 0.8, Ga As  : R= 1.803 1.889
  • 0.9, Ga As  : R= 2.029 2.126


  • (小谷)エネルギーは低い程よい、という意味では、大きめの方がよい。しかし収束安定性(何回で収束するのかにも反映する)も重要。エネルギーは差分に意味があるので。

↓8,9のみ。

Na

  • Na  : R= 3.000

change R

  • 0.7, Na  : R= 2.414
  • 0.8, Na  : R= 2.759
  • 0.9, Na  : R= 3.000

R mol

  • Na  : R= 2.825

EH shift 0.1

  • 0.6, Na  : R= 2.069
  • 0.7, Na  : R= 2.414
  • 0.8, Na  : R= 2.759
  • 1.0, Na  : R= 3.000


EH shift 0.1, EH[2:4]

  • 0.6, Na  : R= 2.069
  • 0.7, Na  : R= 2.414
  • 0.8, Na  : R= 2.759
  • 1.0, Na  : R= 3.000


NaCl

  • Na Cl  : R= 2.639 2.690

Naのs channelを除かないと計算ができない。

change R

  • 0.7, Na Cl  : R= 1.847 1.883
  • 0.8, Na Cl  : R= 2.111 2.152
  • 0.9, Na Cl  : R= 2.375 2.421
  • 1.0, Na Cl  : R= 2.639 2.690
  • 1.1, Na Cl  : R= 2.903 2.959

↓Na*=1.1にしてCl*=0.7-0.9とCl*1.1にしてNa*=0.7-0.8まで計算。 pwemax=2のみ計算可能。点をみやすくするためpwemaxをすこしづつずらしている。

R mol

  • Na Cl  : R= 2.825 1.640

Si

  • Si  : R= 2.222

change R

  • 0.7, Si  : R= 1.555
  • 0.8, Si  : R= 1.777
  • 0.9, Si  : R= 1.999
  • 1.0, Si  : R= 2.222


EH shift 0.1

  • 0.6, Si  : R= 1.555
  • 0.7, Si  : R= 1.555
  • 0.8, Si  : R= 1.777
  • 0.9, Si  : R= 1.999
  • 1.0, Si  : R= 2.222

SiO2

  • Si O  : R= 1.649 1.343

change R

  • 0.7, Si O  : R= 1.154 0.940
  • 0.8, Si O  : R= 1.319 1.075
  • 0.9, Si O  : R= 1.484 1.209


TiO2

  • T1 O1  : R= 2.063 1.579

R mol

  • T1 O1  : R= 2.063 1.579

↓上と同じだ。

ZrO2

  • Zr O  : R= 2.331 1.623

GdN

Gd g-channelを入れると収束する。

  • Gd N  : R= 2.772 1.930

Ge(diamond)

  • Ge  : R= 2.312

Au

s channelを全て除く。

  • Au  : R= 2.726


alpha-Al2O3

  • Al O  : R= 1.937 1.571

C graphite

  • C  : R= 1.344

↑pwemaxの表示単位を間違えた。pwemax=[2:12]

R mol

  • C  : R= 0.906

CaCl2

  • Ca Cl  : R= 2.777 2.402

↓cacl2の間違い

EH shift 0.1

(小谷:R*=0.7でcore spiloutの心配はないか?)

  • R*=0.7, Ca,Cl: R= 1.944,R= 1.682
  • R*=0.8, Ca,Cl: R= 2.221,R= 1.922
  • R*=0.9, Ca,Cl: R= 2.499,R= 2.162
  • R*=1.0, Ca,Cl: R= 2.777,R= 2.402


EH shift 0.1, EH[2:4]

  • R*=0.7, Ca,Cl: R= 1.944,R= 1.682
  • R*=0.8, Ca,Cl: R= 2.221,R= 1.922
  • R*=0.9, Ca,Cl: R= 2.499,R= 2.162
  • R*=1.0, Ca,Cl: R= 2.777,R= 2.402

CaVO3(cubic)

  • V Ca O  : R= 2.106 3.000 1.609

↓pwemax=2

SrVO3(cubic)

EH shift 0.1

R*=0.7, V,Sr,O: R=1.436, 2.489, 1.107

R*=0.8, V,Sr,O: R=1.641, 2.845, 1.265

R*=0.9, V,Sr,O: R=1.846, 3.000, 1.423

R*=1.0, V,Sr,O: R=2.051, 3.000, 1.581

LaCuO4 HTT

EH shift 0.1

R*=0.7 La,Cu,O1,O2: R= 1.908, 1.370, 1.172, 1.215

R*=0.8 La,Cu,O1,O2: R= 2.180, 1.566. 1.340, 1.389

R*=0.9 La,Cu,O1,O2: R= 2.453, 1.762, 1.507, 1.562

R*=1.0 La,Cu,O1,O2: R= 2.725, 1.958, 1.675, 1.736

MgO

  • Mg O  : R= 2.176 1.803

change R

  • 0.7, Mg O  : R= 1.523 1.262
  • 0.8, Mg O  : R= 1.741 1.443
  • 0.9, Mg O  : R= 1.958 1.623
  • 1.0, Mg O  : R= 2.176 1.803

EH shift 0.1

  • 0.7, Mg O  : R= 1.523 1.262
  • 0.8, Mg O  : R= 1.741 1.443
  • 0.9, Mg O  : R= 1.958 1.623
  • 1.0, Mg O  : R= 2.176 1.803

NH4ClO4

  • Cl O1 O2 O3 N2 H1 H2 H3  : R= 1.507 1.269 1.201 1.219 1.190 0.662 0.731 0.682

Na2Co2O4

  • Co Na2 O  : R= 2.004 2.683 1.612

a toy model of Na0.74CoO2.

MODE=12 計算できず。

change R

  • 0.7, Co Na O  : R= 1.403 1.878 1.129
  • 0.8, Co Na O  : R= 1.603 2.147 1.290
  • 0.9, Co Na O  : R= 1.804 2.415 1.451


R mol

  • Co Na O  : R= 2.004 2.825 1.058

"R mol"はEH,EH2=-0.1,-2.0

BiVO4

  • Bi V O  : R= 3.000 1.828 1.435

Change R

以下のR依存性は、(おそらく)、VのMTが小さいために、3s,3pコアのfrozen charge近似が破綻しているため。 R=0.9でもspiloutが0.561と相当に大きい。 lmfのsiploutにも表示されるがlmfaのcoretail: q=0.658, rho(rmt)=0.561にも表示される。 まずは、SPEC_ATOMnのVのセクションにPZ=0,13.9 P=4.2,4.2を挿入し、 local orbital(というよりもMTO)を入れてMT依存性がないことを確認する必要がある。--TakaoKotani 14:07, 28 December 2010 (JST)

  • 0.7, Bi V O  : R= 2.236 1.280 1.004
  • 0.8, Bi V O  : R= 2.556 1.462 1.148
  • 0.9, Bi V O  : R= 2.875 1.645 1.291


FeS2,ortho

  • Fe S  : R= 1.613 1.665

kmxa=5で計算可能

FeS2,ortho,AF

FeS2,cubic

  • Fe S  : R= 2.241 2.037 (小谷:なぜcubicとRがかなり違うのか?)

pwemax=2のみ収束。

HfO2

  • Hf O1 O2  : R= 2.264 1.576 1.821

LaGaO3

  • Ga La O  : R= 2.036 2.997 1.626

Li/C2

  • Li  : R= 1.470

LiFePO4

change R, paramagnetic

  • R*=0.7, Li Fe P O1 O2 O3  : R= 1.860 1.791 1.121 0.908 0.899 0.937
  • R*=0.8, Li Fe P O1 O2 O3  : R= 2.125 2.046 1.281 1.037 1.027 1.070
  • R*=0.9, Li Fe P O1 O2 O3  : R= 2.391 2.302 1.441 1.167 1.156 1.204
  • R*=1.0, Li Fe P O1 O2 O3  : R= 2.657 2.558 1.601 1.297 1.284 1.338

change R, ferromagnetic

  • R*=0.7, Li Fe P O1 O2 O3  : R= 1.860 1.791 1.121 0.908 0.899 0.937
  • R*=0.8, Li Fe P O1 O2 O3  : R= 2.125 2.046 1.281 1.037 1.027 1.070
  • R*=0.9, Li Fe P O1 O2 O3  : R= 2.391 2.302 1.441 1.167 1.156 1.204
  • R*=1.0, Li Fe P O1 O2 O3  : R= 2.657 2.558 1.601 1.297 1.284 1.338


LiのRだけを標準値から変えても計算できない。


Li,Fe両方を変えると計算できる。

YBCO

  • Y Ba Cu1 Cu2 O2 O4  : R= 2.843 3.000 1.829 1.981 1.692 1.575

kmxa=5

beta-Na2V6O15

  • V1 V2 V3 Na O1 O2 O3 O4 O5 O6 O7 O8  : R= 1.751 1.699 1.795 1.847 1.861 1.785 1.883 1.385 1.590 1.349 1.789 1.415

InVO4

  • In V O O2  : R= 2.508 1.758 1.578 1.390

kmxa=5

CaB6

EH shift 0.1

R*=0.7, Ca,B:R=2.742, 1.102

R*=0.8, Ca,B:R=3.0, 1.260

R*=0.9, Ca,B:R=3.0, 1.417

R*=1.0, Ca,B:R=3.0, 1.575

EH shift 0.1, Ca:EH[2:4]


NiO,cubic,paramagnetic

  • Ni O  : R= 2.193 1.755

R mol

  • Ni O  : R= 2.193 1.058



NiO, AF II

change R

  • 0.7, Mnup Mndn O  : R= 2.188 2.188 1.752
  • 0.8, Mnup Mndn O  : R= 1.751 1.751 1.401
  • 0.9, Mnup Mndn O  : R= 1.969 1.969 1.577
  • 1.0, Mnup Mndn O  : R= 2.188 2.188 1.752

EH shift 0.1

  • 0.7, Mnup Mndn O  : R= 2.188 2.188 1.752
  • 0.8, Mnup Mndn O  : R= 1.751 1.751 1.401
  • 0.9, Mnup Mndn O  : R= 1.969 1.969 1.577
  • 1.0, Mnup Mndn O  : R= 2.188 2.188 1.752

EH shift 0.1, Ni:EH[2:4]

  • 0.7, Mnup Mndn O  : R= 2.188 2.188 1.752
  • 0.8, Mnup Mndn O  : R= 1.751 1.751 1.401
  • 0.9, Mnup Mndn O  : R= 1.969 1.969 1.577
  • 1.0, Mnup Mndn O  : R= 2.188 2.188 1.752

TiS2

  • Ti S  : R= 1.357 1.335

C60(sc)

  • C  : R= 1.292

現状ではサイズが大きすぎる。

Ag

  • Ag  : R= 2.733

YH3

  • Y HO HT  : R= 2.887 2.187 1.507

CN(9.0)

  • C  : R= 1.355

重い

CN(6,6)

  • C  : R= 1.354

重い。

pwemax=2は収束。

LSMO

EH=-0.1,-2.0でpwemax=2は収束

R*=0.7, Mn,La,O: 1.437, 2.499, 1.126

R*=0.8, Mn,La,O: 1.642, 2.856, 1.287

R*=0.9, Mn,La,O: 1.848, 3.000, 1.448

R*=1.0, Mn,La,O: 2.053,3.000,1.609

R*=1.0が赤

change R

change R, La-Sr,Mn:EH[2:4]

EH shift 0.1

EH shift 0.1, La-Sr,Mn[2:4]

MnO AF

(R*=0.7が小さすぎないか?違うところへ行ってるように見える)。 Rが小さい時に妙な依存性に見える場合はcore spilloutの性であることが多い。

R*=0.7, Mn,O: R= 1.669, 1.271

R*=0.8, Mn,O: R= 1.908, 1.452

R*=0.9, Mn,O: R= 2.146, 1.634

R*=1.0, Mn,O: R= 2.385, 1.815

EH[2:4]

Personal tools