SPR-KKR (ARPES)

xbandを用いてARPESの入力を見てみましたが、基本的にディフォルトの値で良い様です。ただし、range of polar angleとPhoton Ene.は入力が必要です。
THETA MIN 0.0
THETA MAX 0.0
NT 1
EWORK_EV 5.0
except above setting, glibc (double free or corruption) error is occured.

パラメーターの意味を理解していても、値が発散せずに計算が成功するパラーメーターを見つけるのが難しい。
幾つか計算して成功体験を積み重ねることが必要。

success
------------------
THETA MIN -6
THETA MAX 6
NT 13
------------------


IQ:表面の最上層である特殊な原子位置
MILLER HKL 3D: 表面のミラー指数
※ RSLAB_FMT: slab モデルのことか? 私は使ったことが無い。


IMV_FIN_EV: ポテンシャルの虚数部分(eV, 終状態)
IMV_INI_EV: ポテンシャルの虚数部分(eV, 始状態)
EWORK_EV: バルクのインナーポテンシャル(eV)
NE: nuber of energy points
EMIN_EV: EFに関連したエネルギーウィンドウの最小値(そのまま「9999.0」だとディフォルト値になる)
EMAX_EV: EFに関連したエネルギーウィンドウの最大値(そのまま「9999.0」だとディフォルト値になる)
※ 本来は、 EMIN_EVとEMAX_EVを -2 と 1 にするなどにしたいが、私が計算したときにはエラーが出た。

実験値や計算で想定される値 (EWORK_EV の参考値)
----------------------------
Fe2-xV1+xAl: inner potential 13.4 eV, work function 9.2 eV
FeGa3: inner potential 10.0 eV, work function 4.5 eV
Multilayer Graphene: 16.5 eV
Ga1-xMnxAs: inner potential 7.7 eV, work function 5 eV
TaAs: inner potential 15 eV
SnTe: inner potential 8.5 eV
CoSi2/Si(111): inner potential 14.8 eV
CeIrIn5: inner potential 15 eV
CeRu2Si2: inner potential 12 eV
CeNi2Ge2: inner potential 13.8 eV
SrVO3: inner potential 10 eV
SrTiO3: inner potential 12 eV
BaFe2As2: inner potential 13 eV
CaFe2As2: inner potential 20 eV
Bi2Se3:  inner potential 10 eV
URu2Si2: inner potential 12 eV
EuB6: inner potential 11.2 eV
Ba(Fe1-xCox)2As2: inner potential 15 eV
LiFeAs: inner potential 15.4 eV
OuGoGa5: inner potential 10 eV
----------------------------


NT と NP は分割数


入射光
Photon Ene. = 入射光のエネルギー (eV)
Polarisation: 偏光

この「Structure」についてはマニュアルと記述が異なるので完全に把握できず。

NLAT_G_VEC: 逆格子ベクトルの数
N_LAYER: 光電子の計算におけるバルク層の数
※ 推測
LAYDBL FIN: bulk layer, final state ?
LAYDBL INI: bulk layer, initial state ?
------------------------------------------------------------------------------
■ 計算が走った入力条件
(Si/Al組成比を変えた計算も行なったが計算途中でクラッシュした)
いまのところ成功しているのは、Fe2VAlにおいて、入射光のエネルギー、検出器の角度、層数、インナーポテンシャルかえての計算が可能であるということ。私の場合は検出器の角度の分割数は20を超えると計算の途中でクラッシュしました)

After error occurred, remove case_ARPES.tau and case_ARPES.spc files.
-------------------------------------------------------------------
success condition
case_SCF.inp
--------------
CONTROL  DATASET     = AlVFe2
         ADSI        = SCF
         POTFIL      = AlVFe2.pot
         PRINT = 0   
 
MODE     OP = LDA+U-MF  
         LOPT  = { n d d d }
         UEFF1 = 0.00 eV   JEFF1 = 0.00 eV
         UEFF2 = 1.34 eV   JEFF2 = 0.00 eV
         UEFF3 = 1.80 eV   JEFF3 = 0.00 eV
         UEFF4 = 1.80 eV   JEFF4 = 0.00 eV
         LLOYD
 
TAU      BZINT= POINTS  NKTAB= 30
 
ENERGY   GRID={5}  NE={50}
         ImE=0.001 Ry   EMIN=-0.2 Ry
 
SCF      NITER=200 MIX=0.20 VXC=PBE
         TOL=0.00010  MIXOP=0.20  ISTBRY=1
         QIONSCL=0.80
         NOSSITER
--------------

case_ARPES.inp
--------------
CONTROL  DATASET     = AlVFe2
         ADSI        = ARPES
         POTFIL      = AlVFe2.pot
         PRINT = 0   
 
MODE     OP = LDA+U-MF  
         LOPT  = { n d d d }
         UEFF1 = 0.00 eV   JEFF1 = 0.00 eV
         UEFF2 = 1.34 eV   JEFF2 = 0.00 eV
         UEFF3 = 1.80 eV   JEFF3 = 0.00 eV
         UEFF4 = 1.80 eV   JEFF4 = 0.00 eV
         LLOYD
 
SITES    NL = { 4 }
 
TAU      BZINT= POINTS  NKTAB= 50
 
ENERGY   GRID={1}  NE={60}
         ImE=0.0 Ry
      EMINEV=-5
      EMAXEV=0.5
      EWORK_EV=13.4
     IMV_INI_EV=0.1 IMV_FIN_EV=1.0
 
TASK     ARPES
         STRVER=1
         IQ_AT_SURF=1
         MILLER_HKL={0,0,1}

SPEC_PH THETA=45.0
         PHI=90.0
         POL_P=P
         EPHOT=21.0

SPEC_STR N_LAYDBL={10,10}
         NLAGT_G_VEC=30 N_LAYER=30
         SURF_BAR={0.25,0.25}

SPEC_EL THETA={-20,20}
         PHI={0.0,0.0}
         NT=13
         NP=1
         POL_E=PZ
--------------

case_ARPES.inp
--------------
CONTROL  DATASET     = AlVFe2
         ADSI        = ARPES
         POTFIL      = AlVFe2.pot
         PRINT = 0   
 
MODE     OP = LDA+U-MF  
         LOPT  = { n d d d }
         UEFF1 = 0.00 eV   JEFF1 = 0.00 eV
         UEFF2 = 1.34 eV   JEFF2 = 0.00 eV
         UEFF3 = 1.80 eV   JEFF3 = 0.00 eV
         UEFF4 = 1.80 eV   JEFF4 = 0.00 eV
         LLOYD
 
SITES    NL = { 4 }
 
TAU      BZINT= POINTS  NKTAB= 50
 
ENERGY   GRID={1}  NE={60}
         ImE=0.0 Ry
      EMINEV=-5
      EMAXEV=0.5
      EWORK_EV=13.4
     IMV_INI_EV=0.1 IMV_FIN_EV=1.0
 
TASK     ARPES
         STRVER=1
         IQ_AT_SURF=1
         MILLER_HKL={0,0,1}

SPEC_PH THETA=45.0
         PHI=90.0
         POL_P=P
         EPHOT=21.0

SPEC_STR N_LAYDBL={1,1}
         NLAGT_G_VEC=10 N_LAYER=10
         SURF_BAR={0.25,0.25}

SPEC_EL THETA={-20,20}
         PHI={0.0,0.0}
         NT=17
         NP=1
         POL_E=PZ
--------------

case_ARPES.inp
--------------
CONTROL  DATASET     = AlVFe2
         ADSI        = ARPES
         POTFIL      = AlVFe2.pot
         PRINT = 0   
 
MODE     OP = LDA+U-MF  
         LOPT  = { n d d d }
         UEFF1 = 0.00 eV   JEFF1 = 0.00 eV
         UEFF2 = 1.34 eV   JEFF2 = 0.00 eV
         UEFF3 = 1.80 eV   JEFF3 = 0.00 eV
         UEFF4 = 1.80 eV   JEFF4 = 0.00 eV
         LLOYD
 
SITES    NL = { 4 }
 
TAU      BZINT= POINTS  NKTAB= 50
 
ENERGY   GRID={1}  NE={60}
         ImE=0.0 Ry
      EMINEV=-5
      EMAXEV=2
      EWORK_EV=13.4
     IMV_INI_EV=0.1 IMV_FIN_EV=1.0
 
TASK     ARPES
         STRVER=1
         IQ_AT_SURF=1
         MILLER_HKL={0,0,1}

SPEC_PH THETA=45.0
         PHI=90.0
         POL_P=P
         EPHOT=21.0

SPEC_STR N_LAYDBL={1,1}
         NLAGT_G_VEC=10 N_LAYER=10
         SURF_BAR={0.25,0.25}

SPEC_EL THETA={-20,20}
         PHI={0.0,0.0}
         NT=17
         NP=1
         POL_E=PZ
--------------
---------------------------------------------------------------------

success condition
----------------------------------------------------------
case_SCF.inp
--------------
CONTROL  DATASET     = AlVFe2
         ADSI        = SCF
         POTFIL      = AlVFe2.pot
         PRINT = 0   
 
MODE     OP = LDA+U-MF  
         LOPT  = { n d d d }
         UEFF1 = 0.00 eV   JEFF1 = 0.00 eV
         UEFF2 = 1.34 eV   JEFF2 = 0.00 eV
         UEFF3 = 1.80 eV   JEFF3 = 0.00 eV
         UEFF4 = 1.80 eV   JEFF4 = 0.00 eV
         LLOYD
 
TAU      BZINT= POINTS  NKTAB= 50
 
ENERGY   GRID={5}  NE={32}
         ImE=0.0 Ry   EMIN=-0.2 Ry
 
SCF      NITER=200 MIX=0.20 VXC=VWN
         TOL=0.00001  MIXOP=0.20  ISTBRY=1
         QIONSCL=1.0
         NOSSITER
--------------

case_dos.inp
--------------
CONTROL  DATASET     = AlVFe2
         ADSI        = DOS
         POTFIL      = AlVFe2.pot
         PRINT = 0   
 
MODE     OP = LDA+U-MF  
         LOPT  = { n d d d }
         UEFF1 = 0.00 eV   JEFF1 = 0.00 eV
         UEFF2 = 1.34 eV   JEFF2 = 0.00 eV
         UEFF3 = 1.80 eV   JEFF3 = 0.00 eV
         UEFF4 = 1.80 eV   JEFF4 = 0.00 eV
         LLOYD
 
TAU      BZINT= POINTS  NKTAB= 100
 
ENERGY   GRID={3}  NE={300}
         ImE=0.001 Ry   EMIN=-0.2   EMAX=1.0 Ry
 
TASK     DOS
--------------

case_ARPES.inp
--------------
CONTROL  DATASET     = AlVFe2
         ADSI        = ARPES
         POTFIL      = AlVFe2.pot
         PRINT = 0   
 
MODE     OP = LDA+U-MF  
         LOPT  = { n d d d }
         UEFF1 = 0.00 eV   JEFF1 = 0.00 eV
         UEFF2 = 1.34 eV   JEFF2 = 0.00 eV
         UEFF3 = 1.80 eV   JEFF3 = 0.00 eV
         UEFF4 = 1.80 eV   JEFF4 = 0.00 eV
         LLOYD
 
SITES    NL = { 4 }
 
TAU      BZINT= POINTS  NKTAB= 100
 
ENERGY   GRID={1}  NE={100}
         ImE=0.001 Ry   EMIN=-0.2   EMAX=1.0 Ry
      EWORK_EV=5.0
     IMV_INI_EV=0.1 IMV_FIN_EV=1.0
 
TASK     ARPES
         STRVER=1
         IQ_AT_SURF=1
         MILLER_HKL={0,0,1}

SPEC_PH THETA=45.0
         PHI=90.0
         POL_P=P
         EPHOT=21.0

SPEC_STR N_LAYDBL={10,10}
         NLAGT_G_VEC=30 N_LAYER=30
         SURF_BAR={0.25,0.25}

SPEC_EL THETA={-6,6}
         PHI={0.0,0.0}
         NT=18
         NP=1
         POL_E=PZ
--------------

old condition
----------------------------------------------------------
case_SCF.inp
---------------------
CONTROL  DATASET     = AlVFe2
         ADSI        = SCF
         POTFIL      = AlVFe2.pot
         PRINT = 0   
 
MODE     OP = LDA+U-MF  
         LOPT  = { n d d d }
         UEFF1 = 0.00 eV   JEFF1 = 0.00 eV
         UEFF2 = 1.34 eV   JEFF2 = 0.00 eV
         UEFF3 = 1.80 eV   JEFF3 = 0.00 eV
         UEFF4 = 1.80 eV   JEFF4 = 0.00 eV
         LLOYD
 
TAU      BZINT= POINTS  NKTAB= 50
 
ENERGY   GRID={5}  NE={32}
         ImE=0.0 Ry   EMIN=-0.2 Ry
 
SCF      NITER=200 MIX=0.20 VXC=VWN
         TOL=0.00001  MIXOP=0.20  ISTBRY=1
         QIONSCL=1.0
         NOSSITER
---------------------

case_DOS.inp
---------------------
CONTROL  DATASET     = AlVFe2
         ADSI        = DOS
         POTFIL      = AlVFe2.pot
         PRINT = 0   
 
MODE     OP = LDA+U-MF  
         LOPT  = { n d d d }
         UEFF1 = 0.00 eV   JEFF1 = 0.00 eV
         UEFF2 = 1.34 eV   JEFF2 = 0.00 eV
         UEFF3 = 1.80 eV   JEFF3 = 0.00 eV
         UEFF4 = 1.80 eV   JEFF4 = 0.00 eV
         LLOYD
 
TAU      BZINT= POINTS  NKTAB= 100
 
ENERGY   GRID={3}  NE={300}
         ImE=0.001 Ry   EMIN=-0.2   EMAX=1.0 Ry
 
TASK     DOS
---------------------


case_DOS.inp
---------------------
CONTROL  DATASET     = AlVFe2
         ADSI        = ARPES
         POTFIL      = AlVFe2.pot
         PRINT = 0   
 
MODE     OP = LDA+U-MF  
         LOPT  = { n d d d }
         UEFF1 = 0.00 eV   JEFF1 = 0.00 eV
         UEFF2 = 1.34 eV   JEFF2 = 0.00 eV
         UEFF3 = 1.80 eV   JEFF3 = 0.00 eV
         UEFF4 = 1.80 eV   JEFF4 = 0.00 eV
         LLOYD
 
SITES    NL = { 4 }
 
TAU      BZINT= POINTS  NKTAB= 100
 
ENERGY   GRID={1}  NE={300}
         ImE=0.001 Ry   EMIN=-0.2   EMAX=1.0 Ry
      EWORK_EV=5.0
     IMV_INI_EV=0.1 IMV_FIN_EV=1.0
 
TASK     ARPES
         STRVER=1
         IQ_AT_SURF=1
         MILLER_HKL={0,0,1}

SPEC_PH THETA=45.0
         PHI=90.0
         POL_P=P
         EPHOT=21.0

SPEC_STR N_LAYDBL={10,10}
         NLAGT_G_VEC=30 N_LAYER=30
         SURF_BAR={0.25,0.25}

SPEC_EL THETA={-6,6}
         PHI={0.0,0.0}
         NT=13
         NP=1
         POL_E=PZ
---------------------

case_ARPES.inp
---------------------
CONTROL  DATASET     = AlVFe2
         ADSI        = ARPES
         POTFIL      = AlVFe2.pot
         PRINT = 0   
 
MODE     OP = LDA+U-MF  
         LOPT  = { n d d d }
         UEFF1 = 0.00 eV   JEFF1 = 0.00 eV
         UEFF2 = 1.34 eV   JEFF2 = 0.00 eV
         UEFF3 = 1.80 eV   JEFF3 = 0.00 eV
         UEFF4 = 1.80 eV   JEFF4 = 0.00 eV
         LLOYD
 
SITES    NL = { 4 }
 
TAU      BZINT= POINTS  NKTAB= 100
 
ENERGY   GRID={1}  NE={300}
         ImE=0.001 Ry   EMIN=-0.2   EMAX=1.0 Ry
      EWORK_EV=5.0
     IMV_INI_EV=0.1 IMV_FIN_EV=1.0
 
TASK     ARPES
         STRVER=1
         IQ_AT_SURF=1
         MILLER_HKL={0,0,1}

SPEC_PH THETA=45.0
         PHI=90.0
         POL_P=P
         EPHOT=21.0

SPEC_STR N_LAYDBL={10,10}
         NLAGT_G_VEC=30 N_LAYER=30
         SURF_BAR={0.25,0.25}

SPEC_EL THETA={-6,6}
         PHI={0.0,0.0}
         NT=13
         NP=1
         POL_E=PZ
---------------------

output data
------------------------------
case_ARPES_data.spc
case_ARPES_data.spc_tot_spc.ps
case_ARPES_data.spc_tot_spc.gp
------------------------------

error
------------------------------
SPEC_EL THETA={-6,6}
         PHI={0.0,0.0}
         NT=26
------------------------------
------------------------------------------------------------------------------


QRコード
携帯用QRコード
アクセス数
ページビュー数
[無料でホームページを作成] [通報・削除依頼]