*************************************************** Input file *************************************************** # # File Name # System.CurrrentDir ./ # default=./ System.Name Si_PBE13 Log.print OFF # ON|OFF System.UseRestartfile yes # NO|YES, default=NO System.Restartfile Si_PBE13 # default=null # # Calculation type # eq.type dirac # sch|sdirac|dirac calc.type vps # ALL|VPS|PAO xc.type GGA # LDA|GGA # # Atom # AtomSpecies 14 max.ocupied.N 3 total.electron 14.0 valence.electron 4.0 # # parameters for solving 1D-differential equations # grid.xmin -8.0 # default=-7.0 rmin(a.u.)=exp(grid.xmin) grid.xmax 3.5 # default= 2.5 rmax(a.u.)=exp(grid.xmax) grid.num 10000 # default=4000 grid.num.output 500 # default=2000 # # SCF # scf.maxIter 40 # default=40 scf.Mixing.Type Simple # Simple|GR-Pulay scf.Init.Mixing.Weight 0.10 # default=0.300 scf.Min.Mixing.Weight 0.001 # default=0.001 scf.Max.Mixing.Weight 0.500 # default=0.800 scf.Mixing.History 7 # default=5 scf.Mixing.StartPulay 4 # default=6 scf.criterion 1.0e-12 # default=1.0e-9 # # Pseudo potetial, cutoff (A.U.) # vps.type MBK # BHS|TM number.vps 6 Blochl.projector.num 4 # default=1 which means KB-form local.type polynomial # Simple|Polynomial local.part.vps 1 # default=0 local.cutoff 1.50 # default=smallest_cutoff_vps local.origin.ratio 1.50 # default=3.0 log.deri.RadF.calc on # ON|OFF log.deri.MinE -2.0 # default=-3.0 (Hartree) log.deri.MaxE 2.0 # default= 2.0 (Hartree) log.deri.num 70 # default=50 ghost.check off # ON|OFF # # Core electron density for partial core correction # pcc.ratio=rho_core/rho_V, # pcc.ratio.origin = rho_core(orgin)/rho_core(ip) # charge.pcc.calc on # ON|OFF pcc.ratio 0.1 # default=1.0 pcc.ratio.origin 9.0 # default=6.0 # # Pseudo atomic orbitals # maxL.pao 3 # default=2 num.pao 5 # default=7 radial.cutoff.pao 7.5 # default=5.0 (Bohr) height.of.wall 20000.0 # default=4000.0 (Hartree) rising.edge 0.2 # default=0.5(Bohr),r1=rc-rising.edge search.LowerE -3.000 # default=-3.000 (Hartree) search.UpperE 20.000 # default=20.000 (Hartree) num.of.partition 1200 # default=300 matching.point.ratio 0.67 # default=0.67 ***************************************************** SCF history in all electron calculations ***************************************************** SCF= 1 Eeigen=-165.7065558976287 (Hartree) NormRD= 0.0875178485593 SCF= 2 Eeigen=-166.2697923111511 (Hartree) NormRD= 0.1312482839868 SCF= 3 Eeigen=-166.2658053362534 (Hartree) NormRD= 0.1310122145882 SCF= 4 Eeigen=-164.1457929276084 (Hartree) NormRD= 0.0267726641962 SCF= 5 Eeigen=-163.7221321568353 (Hartree) NormRD= 0.0063016424936 SCF= 6 Eeigen=-163.6400762169141 (Hartree) NormRD= 0.0015174895908 SCF= 7 Eeigen=-163.6179149555693 (Hartree) NormRD= 0.0003665407741 SCF= 8 Eeigen=-163.6088878751362 (Hartree) NormRD= 0.0000883711150 SCF= 9 Eeigen=-163.6043785207718 (Hartree) NormRD= 0.0000212640691 SCF= 10 Eeigen=-163.6019768430669 (Hartree) NormRD= 0.0000051136757 SCF= 11 Eeigen=-163.6006823633554 (Hartree) NormRD= 0.0000012307788 SCF= 12 Eeigen=-163.5999879884237 (Hartree) NormRD= 0.0000002967803 SCF= 13 Eeigen=-163.5996189293184 (Hartree) NormRD= 0.0000000717445 SCF= 14 Eeigen=-163.5994245993319 (Hartree) NormRD= 0.0000000173947 SCF= 15 Eeigen=-163.5993230938176 (Hartree) NormRD= 0.0000000042308 SCF= 16 Eeigen=-163.5992704087638 (Hartree) NormRD= 0.0000000010325 SCF= 17 Eeigen=-163.5992431905402 (Hartree) NormRD= 0.0000000002528 SCF= 18 Eeigen=-163.5992291740722 (Hartree) NormRD= 0.0000000000621 SCF= 19 Eeigen=-163.5992273527629 (Hartree) NormRD= 0.0000000000473 SCF= 20 Eeigen=-163.5992272535524 (Hartree) NormRD= 0.0000000000466 SCF= 21 Eeigen=-163.5992272410163 (Hartree) NormRD= 0.0000000000465 SCF= 22 Eeigen=-163.5992272284920 (Hartree) NormRD= 0.0000000000464 SCF= 23 Eeigen=-163.5992272159801 (Hartree) NormRD= 0.0000000000463 SCF= 24 Eeigen=-163.5992272034794 (Hartree) NormRD= 0.0000000000462 SCF= 25 Eeigen=-163.5992271909911 (Hartree) NormRD= 0.0000000000461 SCF= 26 Eeigen=-163.5992271785144 (Hartree) NormRD= 0.0000000000460 SCF= 27 Eeigen=-163.5992271660504 (Hartree) NormRD= 0.0000000000459 SCF= 28 Eeigen=-163.5992271535988 (Hartree) NormRD= 0.0000000000458 SCF= 29 Eeigen=-163.5992271411600 (Hartree) NormRD= 0.0000000000457 SCF= 30 Eeigen=-163.5992271287333 (Hartree) NormRD= 0.0000000000456 SCF= 31 Eeigen=-163.5992271163173 (Hartree) NormRD= 0.0000000000456 SCF= 32 Eeigen=-163.5992271039145 (Hartree) NormRD= 0.0000000000455 SCF= 33 Eeigen=-163.5992270915240 (Hartree) NormRD= 0.0000000000454 SCF= 34 Eeigen=-163.5992270791444 (Hartree) NormRD= 0.0000000000453 SCF= 35 Eeigen=-163.5992270667774 (Hartree) NormRD= 0.0000000000452 SCF= 36 Eeigen=-163.5992270544222 (Hartree) NormRD= 0.0000000000451 SCF= 37 Eeigen=-163.5992270420788 (Hartree) NormRD= 0.0000000000450 SCF= 38 Eeigen=-163.5992270297484 (Hartree) NormRD= 0.0000000000449 SCF= 39 Eeigen=-163.5992270174290 (Hartree) NormRD= 0.0000000000448 SCF= 40 Eeigen=-163.5992270051217 (Hartree) NormRD= 0.0000000000000 ***************************************************** Eigenvalues (Hartree) in the all electron calculation ***************************************************** j=l+1/2 j=l-1/2 n= 1 l= 0 -65.5902236996133 -65.5902236996133 n= 2 l= 0 -5.1223958802110 -5.1223958802110 n= 2 l= 1 -3.5014036466474 -3.5251688983298 n= 3 l= 0 -0.3970611127118 -0.3970611127118 n= 3 l= 1 -0.1494948780044 -0.1506531071133 ***************************************************** Energies (Hartree) in the all electron calculation ***************************************************** Eeigen = -163.5992270051217 Ekin = 290.5954270049695 EHart = 132.1957321546315 Exc = -20.6350025306832 Eec = -691.9886813845599 Etot = Ekin + EHart + Exc + Eec Etot = -289.8325247556421 *********************************************************** ** DATA for factorized norm conserving pseudo potentials ** *********************************************************** j.dependent.pseudo.potentials on L=0, dif of log deris for all electrons = 1.0210121651 1.0210121651 L=0, dif of log deris for semi local = 0.3046228060 0.3046228060 L=1, dif of log deris for all electrons = 1.7159598155 0.0797214682 L=1, dif of log deris for semi local = 49.4170347403 10.8866643511 L=2, dif of log deris for all electrons = 46.9527901883 43.8475258369 L=2, dif of log deris for semi local = 3.3381077195 3.1326933355 *********************************************************** ** Core electron densities for PCC ** ***********************************************************