*************************************************** Input file *************************************************** # # File Name # System.CurrrentDir ./ # default=./ System.Name V_CA13 Log.print off # ON|OFF System.UseRestartfile yes # NO|YES, default=NO System.Restartfile V_CA13 # default=null # # Calculation type # eq.type dirac # sch|sdirac|dirac calc.type VPS # ALL|VPS|PAO xc.type LDA # LDA|GGA # # Atom # AtomSpecies 23 max.occupied.N 4 total.electron 23.0 valence.electron 13.0 # # parameters for solving 1D-differential equations # grid.xmin -8.0 # default=-7.0 rmin(a.u.)=exp(grid.xmin) grid.xmax 3.2 # default= 2.5 rmax(a.u.)=exp(grid.xmax) grid.num 12000 # default=4000 grid.num.output 500 # default=2000 # # SCF # scf.maxIter 60 # 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.800 # default=0.800 scf.Mixing.History 7 # default=5 scf.Mixing.StartPulay 7 # default=6 scf.criterion 1.0e-14 # default=1.0e-9 # # Pseudo potetial, cutoff (A.U.) # vps.type MBK # BHS|TM number.vps 7 Blochl.projector.num 4 # default=1 which means KB-form local.type Polynomial # Simple|Polynomial local.part.vps 1 # default=0 local.cutoff 1.40 # default=smallest_cutoff_vps local.origin.ratio 2.20 # 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 100 # 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.06 # default=1.0 pcc.ratio.origin 12.0 # default=6.0 # # Pseudo atomic orbitals # maxL.pao 3 # default=2 num.pao 15 # default=7 radial.cutoff.pao 6.0 # 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 60.000 # default=20.000 (Hartree) num.of.partition 2200 # default=300 matching.point.ratio 0.67 # default=0.67 ***************************************************** SCF history in all electron calculations ***************************************************** SCF= 1 Eeigen=-568.7357073073828 (Hartree) NormRD= 2.1787759276045 SCF= 2 Eeigen=-568.0479527770622 (Hartree) NormRD= 2.0242413389908 SCF= 3 Eeigen=-567.9304873886174 (Hartree) NormRD= 1.9635133660665 SCF= 4 Eeigen=-565.8560639874627 (Hartree) NormRD= 0.4310282019933 SCF= 5 Eeigen=-565.0706520267515 (Hartree) NormRD= 0.1327226367309 SCF= 6 Eeigen=-563.5690511986670 (Hartree) NormRD= 0.0032481268698 SCF= 7 Eeigen=-564.4542504319146 (Hartree) NormRD= 0.0010291486698 SCF= 8 Eeigen=-563.5073589629927 (Hartree) NormRD= 0.0009890002145 SCF= 9 Eeigen=-564.5019398658436 (Hartree) NormRD= 0.0010927136622 SCF= 10 Eeigen=-564.0059959091973 (Hartree) NormRD= 0.0000011505321 SCF= 11 Eeigen=-563.9884550145641 (Hartree) NormRD= 0.0000000207162 SCF= 12 Eeigen=-563.9880281249161 (Hartree) NormRD= 0.0000000076271 SCF= 13 Eeigen=-563.9880480356258 (Hartree) NormRD= 0.0000000009556 SCF= 14 Eeigen=-563.9880606111024 (Hartree) NormRD= 0.0000000000787 SCF= 15 Eeigen=-563.9880616197044 (Hartree) NormRD= 0.0000000000080 SCF= 16 Eeigen=-563.9880618976824 (Hartree) NormRD= 0.0000000000071 SCF= 17 Eeigen=-563.9880619024606 (Hartree) NormRD= 0.0000000000071 SCF= 18 Eeigen=-563.9880619072285 (Hartree) NormRD= 0.0000000000071 SCF= 19 Eeigen=-563.9880619119908 (Hartree) NormRD= 0.0000000000071 SCF= 20 Eeigen=-563.9880619167419 (Hartree) NormRD= 0.0000000000071 SCF= 21 Eeigen=-563.9880619214844 (Hartree) NormRD= 0.0000000000071 SCF= 22 Eeigen=-563.9880619262196 (Hartree) NormRD= 0.0000000000070 SCF= 23 Eeigen=-563.9880619309473 (Hartree) NormRD= 0.0000000000070 SCF= 24 Eeigen=-563.9880619356665 (Hartree) NormRD= 0.0000000000070 SCF= 25 Eeigen=-563.9880619403773 (Hartree) NormRD= 0.0000000000070 SCF= 26 Eeigen=-563.9880619450788 (Hartree) NormRD= 0.0000000000070 SCF= 27 Eeigen=-563.9880619497733 (Hartree) NormRD= 0.0000000000070 SCF= 28 Eeigen=-563.9880619544589 (Hartree) NormRD= 0.0000000000070 SCF= 29 Eeigen=-563.9880619591390 (Hartree) NormRD= 0.0000000000069 SCF= 30 Eeigen=-563.9880619638061 (Hartree) NormRD= 0.0000000000069 SCF= 31 Eeigen=-563.9880619684677 (Hartree) NormRD= 0.0000000000069 SCF= 32 Eeigen=-563.9880619731207 (Hartree) NormRD= 0.0000000000069 SCF= 33 Eeigen=-563.9880619777656 (Hartree) NormRD= 0.0000000000069 SCF= 34 Eeigen=-563.9880619824052 (Hartree) NormRD= 0.0000000000069 SCF= 35 Eeigen=-563.9880619870341 (Hartree) NormRD= 0.0000000000069 SCF= 36 Eeigen=-563.9880619916545 (Hartree) NormRD= 0.0000000000068 SCF= 37 Eeigen=-563.9880619962665 (Hartree) NormRD= 0.0000000000068 SCF= 38 Eeigen=-563.9880620008732 (Hartree) NormRD= 0.0000000000068 SCF= 39 Eeigen=-563.9880620054709 (Hartree) NormRD= 0.0000000000068 SCF= 40 Eeigen=-563.9880620100588 (Hartree) NormRD= 0.0000000000068 SCF= 41 Eeigen=-563.9880620146412 (Hartree) NormRD= 0.0000000000068 SCF= 42 Eeigen=-563.9880620192156 (Hartree) NormRD= 0.0000000000068 SCF= 43 Eeigen=-563.9880620237803 (Hartree) NormRD= 0.0000000000067 SCF= 44 Eeigen=-563.9880620283374 (Hartree) NormRD= 0.0000000000067 SCF= 45 Eeigen=-563.9880620328869 (Hartree) NormRD= 0.0000000000067 SCF= 46 Eeigen=-563.9880620374291 (Hartree) NormRD= 0.0000000000067 SCF= 47 Eeigen=-563.9880620419618 (Hartree) NormRD= 0.0000000000067 SCF= 48 Eeigen=-563.9880620464910 (Hartree) NormRD= 0.0000000000067 SCF= 49 Eeigen=-563.9880620510080 (Hartree) NormRD= 0.0000000000067 SCF= 50 Eeigen=-563.9880620555190 (Hartree) NormRD= 0.0000000000066 SCF= 51 Eeigen=-563.9880620600211 (Hartree) NormRD= 0.0000000000066 SCF= 52 Eeigen=-563.9880620645168 (Hartree) NormRD= 0.0000000000066 SCF= 53 Eeigen=-563.9880620690027 (Hartree) NormRD= 0.0000000000066 SCF= 54 Eeigen=-563.9880620734818 (Hartree) NormRD= 0.0000000000066 SCF= 55 Eeigen=-563.9880620779546 (Hartree) NormRD= 0.0000000000066 SCF= 56 Eeigen=-563.9880620824176 (Hartree) NormRD= 0.0000000000066 SCF= 57 Eeigen=-563.9880620868747 (Hartree) NormRD= 0.0000000000066 SCF= 58 Eeigen=-563.9880620913241 (Hartree) NormRD= 0.0000000000065 SCF= 59 Eeigen=-563.9880620957659 (Hartree) NormRD= 0.0000000000065 SCF= 60 Eeigen=-563.9880621001979 (Hartree) NormRD= 0.0000000000000 ***************************************************** Eigenvalues (Hartree) in the all electron calculation ***************************************************** j=l+1/2 j=l-1/2 n= 1 l= 0 -196.5261997848718 -196.5261997848718 n= 2 l= 0 -22.0680475181952 -22.0680475181952 n= 2 l= 1 -18.3723467150902 -18.6322937514404 n= 3 l= 0 -2.5632549996754 -2.5632549996754 n= 3 l= 1 -1.6029874106754 -1.6338310096939 n= 3 l= 2 -0.1961406370672 -0.1985676097077 n= 4 l= 0 -0.1783092319255 -0.1783092319255 ***************************************************** Energies (Hartree) in the all electron calculation ***************************************************** Eeigen = -563.9880621001979 Ekin = 957.2398638909700 EHart = 397.2336138872300 Exc = -43.1723046688482 Eec = -2258.5628848862707 Etot = Ekin + EHart + Exc + Eec Etot = -947.2617117769188 *********************************************************** ** DATA for factorized norm conserving pseudo potentials ** *********************************************************** j.dependent.pseudo.potentials on L=0, dif of log deris for all electrons = 4.9490085516 4.9490085516 L=0, dif of log deris for semi local = 983.9119790048 983.9119790048 L=1, dif of log deris for all electrons = 0.1260817005 0.2957205990 L=1, dif of log deris for semi local = 66.3163842702 204828.4446661940 L=2, dif of log deris for all electrons = 3.4933957250 2.3095953751 L=2, dif of log deris for semi local = 1452.3359018523 1523.1376649284 L=3, dif of log deris for all electrons = 0.0000505810 0.0130730468 L=3, dif of log deris for semi local = 0.3130436590 0.3287923216 *********************************************************** ** Core electron densities for PCC ** ***********************************************************