*************************************************** Input file *************************************************** # # File Name # System.CurrrentDir ./ # default=./ System.Name Dy_CA13_OC Log.print OFF # ON|OFF System.UseRestartfile yes # NO|YES, default=NO System.Restartfile Dy_CA13_OC # default=null # # Calculation type # eq.type dirac # sch|sdirac|dirac calc.type VPS # ALL|VPS|PAO xc.type LDA # LDA|GGA # # Atom # AtomSpecies 66 max.occupied.N 6 total.electron 66.0 valence.electron 11.0 # # parameters for solving 1D-differential equations # grid.xmin -10.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 1000 # 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.900 # default=0.800 scf.Mixing.History 7 # default=5 scf.Mixing.StartPulay 4 # default=6 scf.criterion 1.0e-13 # default=1.0e-9 # # Pseudo potetial, cutoff (A.U.) # vps.type MBK # BHS|TM number.vps 5 Blochl.projector.num 3 # default=1 which means KB-form local.type polynomial # Simple|Polynomial local.part.vps 1 # default=0 local.cutoff 1.80 # default=smallest_cutoff_vps local.origin.ratio 1.70 # 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.8 # default=1.0 pcc.ratio.origin 8.0 # default=6.0 # # Pseudo atomic orbitals # maxL.pao 4 # default=2 num.pao 15 # default=7 radial.cutoff.pao 8.0 # default=5.0 (Bohr) height.of.wall 30000.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 3000 # default=300 matching.point.ratio 0.67 # default=0.67 PAO.potential.charge 1.0 # default=0.0 ***************************************************** SCF history in all electron calculations ***************************************************** SCF= 1 Eeigen=-8458.5059226693156 (Hartree) NormRD=101382.6802424750786 SCF= 2 Eeigen=-8187.8390468485177 (Hartree) NormRD=97235.8138515305473 SCF= 3 Eeigen=-8180.6499701936555 (Hartree) NormRD=92234.0652225289959 SCF= 4 Eeigen=-8296.3597859784095 (Hartree) NormRD=23204.7034832649042 SCF= 5 Eeigen=-7643.4233364766360 (Hartree) NormRD=2454.6201210020922 SCF= 6 Eeigen=-7483.8077362074855 (Hartree) NormRD=192.5540511434957 SCF= 7 Eeigen=-7599.4524893415810 (Hartree) NormRD= 16.1523170808528 SCF= 8 Eeigen=-7399.4779771827707 (Hartree) NormRD= 38.6144311134100 SCF= 9 Eeigen=-7633.1821034386739 (Hartree) NormRD= 2.2714685650882 SCF= 10 Eeigen=-7564.2568580975212 (Hartree) NormRD= 0.3921140187102 SCF= 11 Eeigen=-7527.7291075469984 (Hartree) NormRD= 0.1877255288754 SCF= 12 Eeigen=-7506.4046633033013 (Hartree) NormRD= 0.0174310550542 SCF= 13 Eeigen=-7509.8384426502826 (Hartree) NormRD= 0.0152597079398 SCF= 14 Eeigen=-7499.2966055735114 (Hartree) NormRD= 0.0486248045460 SCF= 15 Eeigen=-7503.8018394545361 (Hartree) NormRD= 0.0047035744342 SCF= 16 Eeigen=-7505.2986193575416 (Hartree) NormRD= 0.0003264786294 SCF= 17 Eeigen=-7505.6878039319372 (Hartree) NormRD= 0.0000214456746 SCF= 18 Eeigen=-7505.7681739688096 (Hartree) NormRD= 0.0000027426721 SCF= 19 Eeigen=-7505.7820636702336 (Hartree) NormRD= 0.0000007474436 SCF= 20 Eeigen=-7505.7842538338418 (Hartree) NormRD= 0.0000001915727 SCF= 21 Eeigen=-7505.7860703169690 (Hartree) NormRD= 0.0000000306033 SCF= 22 Eeigen=-7505.7871471863018 (Hartree) NormRD= 0.0000000031804 SCF= 23 Eeigen=-7505.7874288132507 (Hartree) NormRD= 0.0000000002139 SCF= 24 Eeigen=-7505.7876915362949 (Hartree) NormRD= 0.0000000000152 SCF= 25 Eeigen=-7505.7876581334067 (Hartree) NormRD= 0.0000000000075 SCF= 26 Eeigen=-7505.7876579654776 (Hartree) NormRD= 0.0000000000074 SCF= 27 Eeigen=-7505.7876579035501 (Hartree) NormRD= 0.0000000000074 SCF= 28 Eeigen=-7505.7876578419382 (Hartree) NormRD= 0.0000000000074 SCF= 29 Eeigen=-7505.7876577805901 (Hartree) NormRD= 0.0000000000074 SCF= 30 Eeigen=-7505.7876577195611 (Hartree) NormRD= 0.0000000000074 SCF= 31 Eeigen=-7505.7876576587832 (Hartree) NormRD= 0.0000000000074 SCF= 32 Eeigen=-7505.7876575983155 (Hartree) NormRD= 0.0000000000073 SCF= 33 Eeigen=-7505.7876575381479 (Hartree) NormRD= 0.0000000000073 SCF= 34 Eeigen=-7505.7876574782540 (Hartree) NormRD= 0.0000000000073 SCF= 35 Eeigen=-7505.7876574186266 (Hartree) NormRD= 0.0000000000073 SCF= 36 Eeigen=-7505.7876573593257 (Hartree) NormRD= 0.0000000000073 SCF= 37 Eeigen=-7505.7876573002868 (Hartree) NormRD= 0.0000000000072 SCF= 38 Eeigen=-7505.7876572415444 (Hartree) NormRD= 0.0000000000072 SCF= 39 Eeigen=-7505.7876571830639 (Hartree) NormRD= 0.0000000000072 SCF= 40 Eeigen=-7505.7876571248844 (Hartree) NormRD= 0.0000000000072 SCF= 41 Eeigen=-7505.7876570669687 (Hartree) NormRD= 0.0000000000072 SCF= 42 Eeigen=-7505.7876570093404 (Hartree) NormRD= 0.0000000000072 SCF= 43 Eeigen=-7505.7876569519585 (Hartree) NormRD= 0.0000000000071 SCF= 44 Eeigen=-7505.7876568949005 (Hartree) NormRD= 0.0000000000071 SCF= 45 Eeigen=-7505.7876568381253 (Hartree) NormRD= 0.0000000000071 SCF= 46 Eeigen=-7505.7876567815883 (Hartree) NormRD= 0.0000000000071 SCF= 47 Eeigen=-7505.7876567253470 (Hartree) NormRD= 0.0000000000071 SCF= 48 Eeigen=-7505.7876566693803 (Hartree) NormRD= 0.0000000000070 SCF= 49 Eeigen=-7505.7876566136538 (Hartree) NormRD= 0.0000000000070 SCF= 50 Eeigen=-7505.7876565582274 (Hartree) NormRD= 0.0000000000070 SCF= 51 Eeigen=-7505.7876565030674 (Hartree) NormRD= 0.0000000000070 SCF= 52 Eeigen=-7505.7876564481712 (Hartree) NormRD= 0.0000000000070 SCF= 53 Eeigen=-7505.7876563935606 (Hartree) NormRD= 0.0000000000070 SCF= 54 Eeigen=-7505.7876563392301 (Hartree) NormRD= 0.0000000000069 SCF= 55 Eeigen=-7505.7876562851161 (Hartree) NormRD= 0.0000000000069 SCF= 56 Eeigen=-7505.7876562312804 (Hartree) NormRD= 0.0000000000069 SCF= 57 Eeigen=-7505.7876561777230 (Hartree) NormRD= 0.0000000000069 SCF= 58 Eeigen=-7505.7876561244666 (Hartree) NormRD= 0.0000000000069 SCF= 59 Eeigen=-7505.7876560714285 (Hartree) NormRD= 0.0000000000069 SCF= 60 Eeigen=-7505.7876560186778 (Hartree) NormRD= 0.0000000000000 ***************************************************** Eigenvalues (Hartree) in the all electron calculation ***************************************************** j=l+1/2 j=l-1/2 n= 1 l= 0 -1967.9850717663357 -1967.9850717663357 n= 2 l= 0 -327.8081283295002 -327.8081283295002 n= 2 l= 1 -282.5401172981612 -311.9106454965491 n= 3 l= 0 -73.1737820916426 -73.1737820916426 n= 3 l= 1 -60.1226676851699 -66.2187440086825 n= 3 l= 2 -46.7834960829244 -48.1156721984706 n= 4 l= 0 -14.6674249925999 -14.6674249925999 n= 4 l= 1 -10.5482742348444 -11.9225216856505 n= 4 l= 2 -5.7658683807044 -6.0155945858569 n= 4 l= 3 -0.3289636231643 -0.3636339771745 n= 5 l= 0 -1.9421818898850 -1.9421818898850 n= 5 l= 1 -0.9973803592232 -1.1827900017571 n= 5 l= 2 -0.0870100143399 -0.0976323007093 n= 6 l= 0 -0.1622016095516 -0.1622016095516 ***************************************************** Energies (Hartree) in the all electron calculation ***************************************************** Eeigen = -7505.7876560186778 Ekin = 13333.8939934904192 EHart = 4737.3954527459055 Exc = -256.1009917180393 Eec = -29974.6085475162108 Etot = Ekin + EHart + Exc + Eec Etot = -12159.4200929979233 *********************************************************** ** DATA for factorized norm conserving pseudo potentials ** *********************************************************** j.dependent.pseudo.potentials on L=0, dif of log deris for all electrons = 106.4500125352 106.4500125352 L=0, dif of log deris for semi local = 3548.8737866708 3548.8737866708 L=1, dif of log deris for all electrons = 20021.8726440627 11236.9452561142 L=1, dif of log deris for semi local = 25536.6949787733 157845.7363320642 L=2, dif of log deris for all electrons = 35.0799949198 3.1259167218 L=2, dif of log deris for semi local = 92.1848674416 8.3146293238 L=3, dif of log deris for all electrons = 16.9917073154 26.1900222317 L=3, dif of log deris for semi local = 3223.0565477913 1623.2845948410 *********************************************************** ** Core electron densities for PCC ** ***********************************************************