PWscf(atomic PP)

  ここではatomic(PWscf)を用いた PAW&USPP 入力ファイルとその結果を掲載していく。
--------------------------------------------------------------------------------
■ atomic
  atomic で作られたPAWポテンシャルとその入力ファイルは下記に掲載されている。  必要となる原子のボタンを押してファイルを開き 間の部分を入力ファイルにすればよい。GIPAW形式にもなっている。GIPAW形式の擬ポテンシャルとしては下記がある。 ただし、こちらは入力ファイルがどのように作られているかは分からない。
--------------------------------------------------------------------------------
■ relationship output with input(出力と入力の関係)

◆ generation configuration (config is OK in occupied state)
-----
/
3
-----

◆ number_of_wfc="5"
-----
$test
/
5
-----

◆ valence configureation (other config is OK in unoccupied state)
-----
$test
/
5
(here)
-----

◆ has_wfc="T" -> lsave_wfc=.true.,

◆ core_correction="T" -> nlcc=.true.,

◆ l_local="2" -> lloc=2,

Note
◆ iswitch=3 -> PP generation
◆ rel=1 -> scalar relativistic, rel=2 -> full relativistic, rel=0 -> non relativistic
◆ pseudotype=1 -> norm-conserving (single-projector), pseudotype=2 -> norm-conserving (multi-projector) , pseudotype=3 -> ultrasoft
◆ tm=.true. -> Troullier-Martins, tm=.false. -> rrkj
※ in tm=.false. and error case, set rho(0) = 0.001

[1] ftp://ftp.ictp.it/shared/qetutorial/espresso-4.2.1/Doc/INPUT_LD1.html 
[2] http://indico.ictp.it/event/a11191/session/17/contribution/9/material/0/0.pdf
[3] https://stuff.mit.edu/afs/athena/software/espresso_v5.0.1/i386_ubuntu1104/espresso-5.0.1/atomic/Doc/INPUT_LD1.html

pseudotype = 3 -> USPP
lpaw -> PAW
weite_coulomb -> 1/r
--------------------------------------------------------------------------------
■ Co hole (NC)
$input
  title='Coh',
  zed=27.0,
  config='1s1 2s2 2p6 3s2 3p6 3d7 4s0 4p0 4d0',
  iswitch=3,
  dft='PBE'
  rel=1,
/
$inputp
  pseudotype=2,
  file_pseudopw='Co.star1s-pbe-mt_gipaw.UPF',
  lloc=2,
  tm=.true.,
  lgipaw_reconstruction=.true.,
/
3
3S 1 0 2.00 0.00 1.100 1.100
3P 2 1 6.00 0.00 1.100 1.100
3D 3 2 7.00 0.00 1.100 1.100
--------------------------------------------------------------------------------
■ Co hole (US)
Co.pbe-n-van_gipaw.UPF (Method: Rappe Rabe Kaxiras Joannopoulos): 1s2 -> 1s1 (rrkj -> tm)
 &input
   title='Co h',
   zed=27.,
   rel=1,
   config='1s1 2s2 2p6 3s2 3p6 3d8 4s1 4p0 4d0 5s0 5p0',
   iswitch=3,
   dft='PBE'
 /
 &inputp
   pseudotype=3,
   file_pseudopw='Co.star1s-pbe-n-tm_gipaw.UPF',
   rcore=1.2,
   lloc=-1,
   rcloc=2.2000,
   rmatch_augfun= 1.3,
   tm=.true.
   nlcc=.true.,
   new_core_ps=.true.,
   lgipaw_reconstruction=.true.
   use_paw_as_gipaw=.true.
 /
6
4S  1 0 1.00  0.00  2.20  2.50  0.0
4S  1 0 0.00  0.10  2.20  2.50  0.0
4P  2 1 0.00  0.00  2.40  2.60  0.0
4P  2 1 0.00  0.05  2.40  2.60  0.0
3D  3 2 8.00  0.00  1.50  2.10  0.0
3D  3 2 0.00  0.05  1.20  2.10  0.0
$test
/
6
4S  1 0 1.00  0.00  2.20  2.50  0.0
5S  1 0 0.00  0.00  2.40  2.60  0.0
4P  2 1 0.00  0.00  2.40  2.60  0.0
5P  2 1 0.00  0.00  2.40  2.60  0.0
3D  3 2 8.00  0.00  2.40  2.60  0.0
4D  3 2 0.00  0.00  2.40  2.60  0.0
--------------------------------------------------------------------------------
■ output
□ All-electron run
Etot = Ekin + Encl + Eh + Exc

□ Testing the pseudopotential
Etotps = Ekin + Encl + Ehrt + Ecxc - Ecc + Epseu

□ PP_HEADER generated
total_psenergy="-9.033924183255684E+001"
total_psenergy=etotps
etosps: total pseudo-valence energy of PP

References
http://www.quantum-espresso.org/pseudopotentials/unified-pseudopotential-format/
--------------------------------------------------------------------------------
■ Fe2VAl


Al.in (PAW)


Al.ld1.in (USPP)


■ Fe (PAW)
Fe.in (atomic/examples/paw_examples/input/Fe.in)
----
 &input
    title='Fe'
      zed=26,
    config='[Ar] 4s2 3d6 4p0'
    iswitch=3,
    dft='PBE'
    nld=3,
        eminld=-10,
        emaxld=10
        deld=0.01d0,
        rlderiv=2.3,
 /
 &inputp
   pseudotype=3,
   nlcc=.true.,
     new_core_ps = .true.
     rcore=1.2000
   lloc=-2,
     rcloc=2.2000
   file_pseudopw='Fe.pbe-paw_kj.UPF'
   zval=8
   lpaw=.true.
     lnc2paw=.false.
     which_augfun ='BESSEL'
     rmatch_augfun= 2.3000
   author='"Lorenzo Paulatto" '
   tm = .true.
!nX n  l   occ   nrg   rmin   rcut
 /
6
3D  3 2 6.00    0.00    1.50    2.2000
3D  3 2 0.00   -0.40    1.50    2.2000
4S  1 0 2.00    0.00    1.50    2.2000
4S  1 0 0.00    0.30    1.50    2.2000
4P  2 1 0.00    0.00    1.50    2.2000
4P  2 1 0.00    0.70    1.50    2.2000
----
Command: $HOME/espresso-5.1/bin/ld1.x < Fe.in


■ V (PAW)
V.in (base: atomic/examples/paw_examples/input/Fe.in)
----
 &input
    title='V'
      zed=23,
    config='[Ar] 4s2 3d3 4p0'
    iswitch=3,
    dft='PBE'
    nld=3,
        eminld=-10,
        emaxld=10
        deld=0.01d0,
        rlderiv=2.3,
 /
 &inputp
   pseudotype=3,
   nlcc=.true.,
     new_core_ps = .true.
     rcore=1.2000
   lloc=-2,
     rcloc=2.2000
   file_pseudopw='V.pbe-paw_kj.UPF'
   zval=5
   lpaw=.true.
     lnc2paw=.false.
     which_augfun ='BESSEL'
     rmatch_augfun= 2.3000
   tm = .true.
!nX n  l   occ   nrg   rmin   rcut
 /
6
3D  3 2 3.00    0.00    1.50    2.2000
3D  3 2 0.00   -0.40    1.50    2.2000
4S  1 0 2.00    0.00    1.50    2.5000
4S  1 0 0.00    0.30    1.50    2.5000
4P  2 1 0.00    0.00    1.50    2.5000
4P  2 1 0.00    0.70    1.50    2.5000
----
Command: $HOME/espresso-5.1/bin/ld1.x < V.in


■ Al (PAW)
Al.in (base: atomic/examples/paw_examples/input/si.in)
-----
 &input
    title='Al'
      zed=13,
    config='[Ne] 3s2 3p1 3d-1'
    iswitch=3,
    dft='PBE'
    nld=3,
        eminld=-4,
        emaxld=4,
        deld=0.01d0,
        rlderiv=1.85,
 /
 &inputp
   pseudotype=3,
   nlcc=.true.,
     new_core_ps = .true.
     rcore=1.2
   lloc=2,
   file_pseudopw='Al.pbe-paw_kj.UPF'
   zval=3
   lpaw=.true.
     lnc2paw=.false.
     which_augfun ='BESSEL'
     rmatch_augfun= 1.80
     tm=.true.
   !bessel_vloc=.true.
!   lgipaw_reconstruction=.true.
!nX n  l   occ   nrg   rmin   rcut
 /
5
3S  1  0  2.00  0.00  1.50  1.80
3S  1  0  0.00  2.00  1.50  1.80
3P  2  1  1.00  0.00  1.50  1.80
3P  2  1  0.00  2.00  1.50  1.80
3D  3  2 -1.00  0.20  1.65  1.65
-----
Command: $HOME/espresso-5.1/bin/ld1.x < Al.in


■ Al (USPP)
Al.ld1.in
-----
 &input
   atom='Al',  dft='PBE',  zed=13.0, config='1s2 2s2 2p6 3s2 3p1',
   iswitch=3
 /
 &inputp
   pseudotype=3, nlcc=.true.,rcore=1.4, lloc=1,
   file_pseudopw='Al.pbe-nsp-rrkj-us.UPF',
 /
2
3S 3 0 2.00  0.00 1.7 1.75
3P 3 1 1.00  0.00 1.7 1.75
&test
   configts(1)='3s2 3p1',
/
-----
Command: $HOME/espresso-5.1/bin/ld1.x < Al.ld1.in
--------------------------------------------------------------------------------
■ NCPPの作成(調査中)
1. カットオフ半径の推定
  Ti_r.ld1.in
  -----
 &input
   atom='Ti', dft='PBE', config='[Ar] 3d2 4s2 4p0'
/
  -----
2. $HOME/espresso-5.1/bin/ld1.x < Ti_r.ld1.in
3. 下記の結果を得る。
  s(3S/3S) =  1.000000  =   1.0069  =    1.1699  r(max) =   0.8702
  s(3P/3P) =  1.000000  =   1.0860  =    1.3907  r(max) =   0.8985
  s(3D/3D) =  1.000000  =   1.6171  =    3.5729  r(max) =   0.9811
  s(4S/4S) =  1.000000  =   3.5138  =   14.2491  r(max) =   2.9123
  s(4P/4P) =  1.000000  =   4.8653  =   27.9369  r(max) =   3.8227
  3S、3P、及び3Dはほぼ同じ r (a.u. 単位)の値である。擬ポテンシャルにするには4S や 4P は余りにも大きい。3Dより内側の軌道を内殻にすれば良さそうである。
※ カットオフ半径 rc は典型的に ρcore(rc) >= 2*ρvalence(rc) となるものが選択される。 ρcore(r) = ρvalence(r)近傍で非線形の効果が重要になってくる。 内殻補正のカットオフ半径は小さい方が精度が高くなるが、平面波の大きなカットオフ値が必要になる。
4. 入力ファイルの作成
  Ti.ld1.in
  -----
&input
   atom='Ti',  dft='PBE',  config='[Ar] 3d2 4s2 4p0',
   iswitch=3
/
 &inputp
   pseudotype=1, nlcc=.true., rcore=2.0, lloc=1,
   file_pseudopw='Ti.pbe-n-rrkj.UPF',
   rho0=0.001,
/
3
3D 3 2 2.00  0.00 1.3 1.3
4S 4 0 2.00  0.00 2.9 2.9
4P 4 1 0.00  0.00 2.9 2.9
 &test
   configts(1)='3d2 4s2 4p0'
/
  -----
  ここで各変数の意味は下記のようになる。
  • iswitch=3 : 擬ポテンシャルを生成する
  • pseudotype=1 : ノルム保存型擬ポテンシャルを生成する
  • nlcc = . true.  : 非線形内殻補正を入れる
  • rcore  : 内殻での電荷の平滑化の範囲(半径でa.u.単位)。設定している全ての軌道に対するカットオフ半径よりも小さな値が望ましいようである(要確認)。
  • lloc  : l軌道の局所参照ポテンシャルを用いる。私の経験では、3Dや4Dと記述している部分で最外殻のlの軌道が lloc の値に対応していないとエラーがよく発生した。
  • 下の方にある右側の数値 : カットオフ半径。カットオフ半径は上記で得られた r を参考する。カットオフ半径が小さいと多くのカットオフを必要とする。
※ エラー(phi has nodes before r_c)が出て上手くいかない場合は、TM型の擬ポテンシャルとする tm=.ture. か 小さなρ(r = 0)  となるように(例えば、rho0=0.001) を &inputp 内に記述する。これはRRKJでは3つか4つのベッセル関数を用いているためである。指数関数を用いた多項式によるTM型であれば問題ないが非常に頑丈だ が粗野である。rho0=0.001 は不運にもs軌道に対してしか働かない。
5. $HOME/espresso-5.1/bin/ld1.x < Ti.ld1.in
6. 下記の結果を得る。
  Wfc   3D  rcut= 1.296  Estimated cut-off energy=      137.82 Ry
  This function has    0 nodes for 0 < r <    1.296
 これは、3D でのカットオフが137.82 Ry 必要となることを意味する。これは計算できないカットオフの大きさではないが、比較的大きなカットオフ値であり計算が大変になる。
  3 2     3D   1( 2.00)       -0.40319       -0.40457        0.00138
  1 0     4S   1( 1.00)       -0.38394       -0.38420        0.00026
  2 1     4P   1( 1.00)       -0.15248       -0.15237       -0.00011
  dEtot_ae =       0.226061 Ry
  dEtot_ps =       0.226250 Ry,   Delta E=      -0.000189 Ry
  上記の記述も得られる。これらは全電子計算(All Electron: AE)と擬ポテンシャル(PP)でのエネルギー差である。右側にある数値が小さいことを確認する。1 mRy 程度である。Delta E は 0.2 mRy 程度である。Delta E は 0.001 Ryのオーダーが非常に良いが、0.01 Ry のオーダーでも受け入れられる精度である。あなたが欲しいエネルギーの精度とカットオフを rcore や カットオフ半径を変えて探すと良い。
7. 準内殻を得たい場合には、例えば下記のようにすれば良い。
  -----
&input
   atom='Ti',  dft='PBE',  config='[Ar] 3d2 4s2 4p0',
   iswitch=3
/
&inputp
   pseudotype=1, nlcc=.true., lloc=-1, rcloc=1.3
   file_pseudopw='Ti.pbe-sp-rrkj.UPF'
/
3
3S 3 0 2.00  0.00 1.1 1.1
3P 3 1 6.00  0.00 1.2 1.2
3D 3 2 2.00  0.00 1.3 1.3
&test
   configts(1)='3s2 3p6 3d2 4s2 4p0',
/
  -----
 3D, 4S, 4P の部分が 3S, 3P, 3D になっていることが分かる。エラーが出ずに欲しいエネルギー精度で小さいカットオフが得られるように「rcore、lloc 、各軌道のカットオフ半径」 や「lloc=-1 と rcloc 」を調整する。
※ lloc=-1 とすると各軌道でのカットオフ半径の記述が無視されて、rcloc の値が採用される。
※ rcore は自動で計算されるが、s軌道が最外殻の系では rcore を指定しないと計算が走らないことが多い。
8. check
  -----
  gnuplot
plot [-2:2][-20:20] 'ld1.dlog' u 1:2 w l lt 1, 'ld1.dlog' u 1:3 w l lt 2,\
                    'ld1.dlog' u 1:4 w l lt 3, 'ld1ps.dlog' u 1:2 lt 1, \
                  'ld1ps.dlog' u 1:3     lt 2, 'ld1ps.dlog' u 1:4 lt 3
  -----
  対数微分が内殻の軌道が存在するエネルギー近傍以外でほとんど一致していることを確認する。上手く線に乗っていれば、散乱(位相シフト)が再現できる。
 ※ 上手くフィットしない場合は、3S などでの 左から5つめの 0.00 を適切な値にすると上手くいくかもしれない。
  -----
  gnuplot 
plot [0:5] 'ld1ps.wfc' u 1:2 lt 1    , 'ld1ps.wfc' u 1:3 lt 3    , \
           'ld1ps.wfc' u 1:4 lt 2    , 'ld1ps.wfc' u 1:5 lt 1 w l, \
           'ld1ps.wfc' u 1:6 lt 3 w l, 'ld1ps.wfc' u 1:7 lt 2 w l
  -----
 設定したカットオフ半径以上で、線と一致していれば、上手く動径分布関数を再現できていることが分かる。


■ PAWの作成(調査中)
PAW (atomic input form): http://theossrv1.epfl.ch/Main/Pseudopotentials

1. Fe.in
  -----
 &input
    title='Fe'
      zed=26,
    config='[Ar] 4s2 3d6 4p0'
    iswitch=3,
    dft='PBE'
    nld=3,
        eminld=-10,
        emaxld=10
        deld=0.01d0,
        rlderiv=2.3,
 /
 &inputp
   pseudotype=3,
   nlcc=.true.,
     new_core_ps = .true.
     rcore=1.2000
   lloc=-2,
     rcloc=2.2000
   file_pseudopw='Fe.pbe-paw_kj.UPF'
   zval=8
   lpaw=.true.
     lnc2paw=.false.
     which_augfun ='BESSEL'
     rmatch_augfun= 2.3000
   author='"Lorenzo Paulatto" '
   tm = .true.
!nX n  l   occ   nrg   rmin   rcut
 /
6
3D  3 2 6.00    0.00    1.50    2.2000
3D  3 2 0.00   -0.40    1.50    2.2000
4S  1 0 2.00    0.00    1.50    2.2000
4S  1 0 0.00    0.30    1.50    2.2000
4P  2 1 0.00    0.00    1.50    2.2000
4P  2 1 0.00    0.70    1.50    2.2000
  -----
  • nld=3 : the number of logarithmic derivatives to be calculated
  • eminld=-10 : Energy range (min energy, in Ry) at which logarithmic derivatives are calculated
  • emaxld=10 : Energy range (max energy, in Ry) at which logarithmic derivatives are calculated
  • deld=0.01d0 : Delta e (Ry) of energy for logarithmic derivatives
  • rlderiv=2.3 : radius (a.u.) at which logarithmic derivatives are calucalted
  • new_core_ps = .true. : pseudized the core charege with bessel functions
  • zval=8 : zval is automatically calculated from available data. if the value of zval is provided in input, it will be checked versus the calculated value. The only case in which you need to explicitly provide the value of zval for noninteger zval (i.e. half core-hole psudo-potentials)
  • lpaw=.true. : produce a PAW dataset, experimental feature only for psudotype=3
  • lnc2paw=.false. :
  • which_augfun ='BESSEL' : Use Bessel functions to pseudize the Q
  • rmatch_augfun= 2.3000 : pseudization radius for the augmentation functions. Presently it has the same value for all L.
2. $HOME/espresso-5.1/bin/ld1.x < Fe.in

■ 重要な記述[1, 6]
------
rcloc = matching radius (a.u.) for local pseudo-potential (no default) 
  Must be specified only if lloc=-1, otherwise the corresponding value of rcut is used. 
rcore= matching radius (a.u.) for the smoothing of the core charge
  If not specified, the matching radius is determined
  by the condition rho_core(rcore) = 2*rho_valence(rcore)
------
-
in transition metals, whose typical outer electronic configuration is something like (n = main quantum number) ndi(n + 1)sj(n + 1)pk, it is not always evident that the ns and np states (``semicore states'') can be safely put into the core. The problem is that nd states are localized in the same spatial region as ns and np states, deeper than (n + 1)s and (n + 1)p states. This may lead to poor transferability. Typically, PP's with semicore states in the core work well in solids with weak or metallic bonding, but perform poorly in compounds with a stronger (chemical) type of bonding.
-
Heavy alkali metals (Rb, Cs, maybe also K) have a large polarizable core. PP's with just one electron may not always give satisfactory results.
-
In some II-VI and III-V semiconductors, such as ZnSe and GaN, the contribution of the d states of the cation to the bonding is not negligible and may require explicit inclusion of those d states into the valence.
------
NC-PP's
even less for 4f electrons : rc = 0.7-0.8 a.u. : >=100 Ry
second-row(2p) elements: rc =1.1-1.2 a.u. : 50-70 Ry
3d transiton metals : rc =1.1-1.2 a.u. : > 80 Ry
highest l have slower convergence for the same rc

US-PP's
second-row (2p) elements : NC rc / US rc = 1.3 / 1.5 a.u.
3d transiton metals : NC rc / US rc = 1.7 / 2.2 a.u.
do not forget that the sum of the rc of two atoms should not exceed the typical bond length of those atoms.
------
  A big problem in NC-PP's is how strike a compromise between softness and transferability, especially for difficult elements. The basic question: ``how much should I push rc outwards in order to have reasonable results with a reasonable PW cutoff''. has no clear-cut answer. The choice of rc at the outermost maximum for ``difficult'' elements (those described in Sec.2.2.1): typically 0.7-0.8 a.u, even less for 4f electrons, yields very hard PP's (more than 100 Ry needed in practical calculations). With a little bit of experience one can say that for second-row (2p) elements, rc = 1.1 - 1.2 will yield reasonably good results for 50-70 Ry PW kinetic energy cutoff; for 3d transition metals, the same rc will require > 80 Ry cutoff (highest l have slower convergence for the same rc). The above estimates are for TM pseudization. RRKJ pseudization will yield an estimate of the required cutoff. For multiple-projectors PP's, the rc of unbound states may be chosen in the same range as for bound states. Use small rc and don't try to push them outwards: the US pseudization will take care of softness. US pseudization radii can be chosen much larger than NC ones (e.g. 1.3÷ 1.5 a.u. for second-row 2p elements, 1.7÷ 2.2 a.u. for 3d transition metals), but do not forget that the sum of the rc of two atoms should not exceed the typical bond length of those atoms.
  Note that it is the hardest atom that determines the PW cuto in a solid or molecule. Do not waste time trying to nd optimally soft PP's for element X if element Y is harder then element X.
------
  the pseudopotential web-page at www.quantum-espresso.org/  with the label "*star1s*_gipaw*" for K-edges, "*star2s*_gipaw*" for L1-edges and so on.
[1] http://www.quantum-espresso.org/wp-content/uploads/Doc/pseudo-gen/node1.html
[2] https://sites.google.com/site/dceresoli/
[3] http://qe-forge.org/pipermail/pw_forum/2007-October/081886.html
[4] http://www.abinit.org/documentation/helpfiles/for-v7.8/tutorial/lesson_paw2.html
[5] http://www.abinit.org/documentation/helpfiles/for-v7.8/tutorial/lesson_paw3.html
[6] http://www.nist.gov/pml/data/dftdata/index.cfm
[7] http://qe-forge.org/pipermail/pw_forum/2011-September/096300.html 
[8] http://www.fisica.uniud.it/~giannozz/QE-Tutorial/tutorial_pseudo.html
[9] ftp://ftp.ictp.it/shared/qetutorial/espresso-4.2.1/Doc/INPUT_LD1.html
[10] http://qe-forge.org/pipermail/pw_forum/2009-October/089260.html
--------------------------------------------------------------------------------
■ under testing ( 1s core-hole type)
USPP: http://www.quantum-espresso.org/pseudopotentials/#

■ Fe_1sh.in
-----
 &input
    title='Fe'
      zed=26,
    config='1s1 2s2 2p6 3s2 3p6 4s2 3d6 4p0'
    iswitch=3,
    dft='PBE'
    nld=3,
        eminld=-10,
        emaxld=10
        deld=0.01d0,
        rlderiv=2.3,
 /
 &inputp
   pseudotype=3,
   nlcc=.true.,
     new_core_ps = .true.
     rcore=1.2000
   lloc=-2,
     rcloc=2.2000
   file_pseudopw='Fe.pbe-star1s-gipaw_kj.UPF'
   zval=9
   lpaw=.true.
     lnc2paw=.false.
     which_augfun ='BESSEL'
     rmatch_augfun= 2.3000
   tm = .true.
   !bessel_vloc=.true.
   lgipaw_reconstruction=.true.
!nX n  l   occ   nrg   rmin   rcut
 /
6
3D  3 2 6.00    0.00    1.50    2.2000
3D  3 2 0.00   -0.40    1.50    2.2000
4S  1 0 2.00    0.00    1.50    2.2000
4S  1 0 0.00    0.30    1.50    2.2000
4P  2 1 0.00    0.00    1.50    2.2000
4P  2 1 0.00    0.70    1.50    2.2000
-----


■ V_1sh.in
-----
 &input
    title='V'
      zed=23,
    config='1s1 2s2 2p6 3s2 3p6 4s2 3d3 4p0'
    iswitch=3,
    dft='PBE'
    nld=3,
        eminld=-10,
        emaxld=10
        deld=0.01d0,
        rlderiv=2.3,
 /
 &inputp
   pseudotype=3,
   nlcc=.true.,
     new_core_ps = .true.
     rcore=1.2000
   lloc=-2,
     rcloc=2.2000
   file_pseudopw='V.pbe-star1s-gipaw_kj.UPF'
   zval=6
   lpaw=.true.
     lnc2paw=.false.
     which_augfun ='BESSEL'
     rmatch_augfun= 2.3000
   tm = .true.
   lgipaw_reconstruction=.true.
!nX n  l   occ   nrg   rmin   rcut
 /
6
3D  3 2 3.00    0.00    1.50    2.2000
3D  3 2 0.00   -0.40    1.50    2.2000
4S  1 0 2.00    0.00    1.50    2.5000
4S  1 0 0.00    0.30    1.50    2.5000
4P  2 1 0.00    0.00    1.50    2.5000
4P  2 1 0.00    0.70    1.50    2.5000
-----


■ Al_1sh.in
-----
 &input
    title='Al'
      zed=13,
    config='1s1 2s2 2p6 3s2 3p1 3d-1'
    iswitch=3,
    dft='PBE'
    nld=3,
        eminld=-4,
        emaxld=4,
        deld=0.01d0,
        rlderiv=1.85,
 /
 &inputp
   pseudotype=3,
   nlcc=.true.,
     new_core_ps = .true.
     rcore=1.2
   lloc=2,
   file_pseudopw='Al.pbe-star1s-gipaw_kj.UPF'
   zval=4
   lpaw=.true.
     lnc2paw=.false.
     which_augfun ='BESSEL'
     rmatch_augfun= 1.80
     tm=.true.
   !bessel_vloc=.true.
   lgipaw_reconstruction=.true.
!nX n  l   occ   nrg   rmin   rcut
 /
5
3S  1  0  2.00  0.00  1.50  1.80
3S  1  0  0.00  2.00  1.50  1.80
3P  2  1  1.00  0.00  1.50  1.80
3P  2  1  0.00  2.00  1.50  1.80
3D  3  2 -1.00  0.20  1.65  1.65
-----


■ Mg_1sh.in (Reference: http://theossrv1.epfl.ch/Main/Pseudopotentials)
-----
  &input
    title='Mg'
      zed=12
    config='1s1 2s2 2p6 3s2 3p0 3d0'
    iswitch=3,
    dft='PBE'
    nld=3,
        eminld=-4,
        emaxld=4,
        deld=0.01d0,
        rlderiv=1.85,
 /
 &inputp
   pseudotype=3,
   nlcc=.true.,
     new_core_ps = .true.
     rcore=1.2
   lloc=2,
   file_pseudopw='Mg.pbe-star1s-gipaw_kj.UPF'
   zval=3
   lpaw=.true.
     lnc2paw=.false.
     which_augfun ='BESSEL'
     rmatch_augfun= 1.80
     tm=.true.
   !bessel_vloc=.true.
   lgipaw_reconstruction=.true.
!nX n  l   occ   nrg   rmin   rcut
 /
5
3S  1  0  2.00  0.00  1.50  1.80
3S  1  0  0.00  2.00  1.50  1.80
3P  2  1  0.00  0.00  1.50  1.80
3P  2  1  0.00  2.00  1.50  1.80
3D  3  2  0.00  0.20  1.65  1.65
-----
--------------------------------------------------------------------------------

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