*************************************************** Input file *************************************************** # # File Name # System.CurrrentDir ./ # default=./ System.Name Nd_CA11 Log.print OFF # ON|OFF System.UseRestartfile yes # NO|YES, default=NO System.Restartfile Nd_CA11 # default=null # # Calculation type # eq.type dirac # sch|sdirac|dirac calc.type vps # ALL|VPS|PAO xc.type LDA # LDA|GGA # # Atom # AtomSpecies 60 max.occupied.N 6 total.electron 60.0 valence.electron 14.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.700 # 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 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.80 # default=smallest_cutoff_vps local.origin.ratio 3.00 # 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 2.0 # default=1.0 pcc.ratio.origin 10.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=-6522.0659579020412 (Hartree) NormRD=62087.7195934754054 SCF= 2 Eeigen=-6500.5982731849626 (Hartree) NormRD=62410.7530586206340 SCF= 3 Eeigen=-6501.3425399846119 (Hartree) NormRD=62278.8448638592308 SCF= 4 Eeigen=-6770.5014057689696 (Hartree) NormRD=16025.0252066721114 SCF= 5 Eeigen=-6820.6939383432500 (Hartree) NormRD=2063.0475326192673 SCF= 6 Eeigen=-6193.6861890734253 (Hartree) NormRD=282.7395470321869 SCF= 7 Eeigen=-5574.0912638784985 (Hartree) NormRD=268.1948398234395 SCF= 8 Eeigen=-6297.6804868591362 (Hartree) NormRD= 65.4408029558476 SCF= 9 Eeigen=-5965.1858577342782 (Hartree) NormRD= 5.0915400946162 SCF= 10 Eeigen=-5980.2191222481679 (Hartree) NormRD= 0.4270715294586 SCF= 11 Eeigen=-5962.8247426872422 (Hartree) NormRD= 1.8555835922612 SCF= 12 Eeigen=-5995.0038846153093 (Hartree) NormRD= 0.1508681797584 SCF= 13 Eeigen=-5986.6896594624850 (Hartree) NormRD= 0.0551071486477 SCF= 14 Eeigen=-5978.0270477203630 (Hartree) NormRD= 0.0088756830606 SCF= 15 Eeigen=-5973.4819223076920 (Hartree) NormRD= 0.0031713015570 SCF= 16 Eeigen=-5971.2667464801270 (Hartree) NormRD= 0.0012694919866 SCF= 17 Eeigen=-5971.1687334163107 (Hartree) NormRD= 0.0003137175694 SCF= 18 Eeigen=-5970.2415498219634 (Hartree) NormRD= 0.0001694146241 SCF= 19 Eeigen=-5971.2223069303909 (Hartree) NormRD= 0.0002568895009 SCF= 20 Eeigen=-5970.7345076003476 (Hartree) NormRD= 0.0000035031011 SCF= 21 Eeigen=-5970.6803255046889 (Hartree) NormRD= 0.0000002193467 SCF= 22 Eeigen=-5970.6759391402411 (Hartree) NormRD= 0.0000000972637 SCF= 23 Eeigen=-5970.6771396392032 (Hartree) NormRD= 0.0000000220677 SCF= 24 Eeigen=-5970.6785554688067 (Hartree) NormRD= 0.0000000023518 SCF= 25 Eeigen=-5970.6786764958788 (Hartree) NormRD= 0.0000000002820 SCF= 26 Eeigen=-5970.6793175139055 (Hartree) NormRD= 0.0000000000804 SCF= 27 Eeigen=-5970.6791765647749 (Hartree) NormRD= 0.0000000000245 SCF= 28 Eeigen=-5970.6791731509293 (Hartree) NormRD= 0.0000000000233 SCF= 29 Eeigen=-5970.6791729497891 (Hartree) NormRD= 0.0000000000232 SCF= 30 Eeigen=-5970.6791727495201 (Hartree) NormRD= 0.0000000000232 SCF= 31 Eeigen=-5970.6791725500407 (Hartree) NormRD= 0.0000000000231 SCF= 32 Eeigen=-5970.6791723514307 (Hartree) NormRD= 0.0000000000230 SCF= 33 Eeigen=-5970.6791721536019 (Hartree) NormRD= 0.0000000000229 SCF= 34 Eeigen=-5970.6791719566281 (Hartree) NormRD= 0.0000000000229 SCF= 35 Eeigen=-5970.6791717604683 (Hartree) NormRD= 0.0000000000228 SCF= 36 Eeigen=-5970.6791715651170 (Hartree) NormRD= 0.0000000000227 SCF= 37 Eeigen=-5970.6791713705679 (Hartree) NormRD= 0.0000000000226 SCF= 38 Eeigen=-5970.6791711768365 (Hartree) NormRD= 0.0000000000226 SCF= 39 Eeigen=-5970.6791709839117 (Hartree) NormRD= 0.0000000000225 SCF= 40 Eeigen=-5970.6791707917937 (Hartree) NormRD= 0.0000000000224 SCF= 41 Eeigen=-5970.6791706004797 (Hartree) NormRD= 0.0000000000224 SCF= 42 Eeigen=-5970.6791704099469 (Hartree) NormRD= 0.0000000000223 SCF= 43 Eeigen=-5970.6791702202390 (Hartree) NormRD= 0.0000000000222 SCF= 44 Eeigen=-5970.6791700313006 (Hartree) NormRD= 0.0000000000222 SCF= 45 Eeigen=-5970.6791698431507 (Hartree) NormRD= 0.0000000000221 SCF= 46 Eeigen=-5970.6791696558093 (Hartree) NormRD= 0.0000000000220 SCF= 47 Eeigen=-5970.6791694692365 (Hartree) NormRD= 0.0000000000220 SCF= 48 Eeigen=-5970.6791692834413 (Hartree) NormRD= 0.0000000000219 SCF= 49 Eeigen=-5970.6791690984346 (Hartree) NormRD= 0.0000000000218 SCF= 50 Eeigen=-5970.6791689142074 (Hartree) NormRD= 0.0000000000218 SCF= 51 Eeigen=-5970.6791687307523 (Hartree) NormRD= 0.0000000000217 SCF= 52 Eeigen=-5970.6791685480666 (Hartree) NormRD= 0.0000000000216 SCF= 53 Eeigen=-5970.6791683661559 (Hartree) NormRD= 0.0000000000216 SCF= 54 Eeigen=-5970.6791681850118 (Hartree) NormRD= 0.0000000000215 SCF= 55 Eeigen=-5970.6791680046181 (Hartree) NormRD= 0.0000000000214 SCF= 56 Eeigen=-5970.6791678250083 (Hartree) NormRD= 0.0000000000214 SCF= 57 Eeigen=-5970.6791676461307 (Hartree) NormRD= 0.0000000000213 SCF= 58 Eeigen=-5970.6791674680353 (Hartree) NormRD= 0.0000000000212 SCF= 59 Eeigen=-5970.6791672906720 (Hartree) NormRD= 0.0000000000212 SCF= 60 Eeigen=-5970.6791671140654 (Hartree) NormRD= 0.0000000000000 ***************************************************** Eigenvalues (Hartree) in the all electron calculation ***************************************************** j=l+1/2 j=l-1/2 n= 1 l= 0 -1590.3800735300401 -1590.3800735300401 n= 2 l= 0 -256.9249267508145 -256.9249267508145 n= 2 l= 1 -224.0046733481479 -243.0387194586771 n= 3 l= 0 -55.6863249507433 -55.6863249507433 n= 3 l= 1 -45.9777925340518 -49.7735787805127 n= 3 l= 2 -34.8549821347975 -35.6747912984140 n= 4 l= 0 -11.0485462907315 -11.0485462907315 n= 4 l= 1 -7.9944171141169 -8.8224780353282 n= 4 l= 2 -4.1921906471310 -4.3380568744829 n= 4 l= 3 -0.0547923992500 -0.0693823881213 n= 5 l= 0 -1.4814092897397 -1.4814092897397 n= 5 l= 1 -0.7749836986615 -0.8869963570898 n= 6 l= 0 -0.1336165806504 -0.1336165806504 ***************************************************** Energies (Hartree) in the all electron calculation ***************************************************** Eeigen = -5970.6791671140654 Ekin = 10454.6707728737420 EHart = 3747.6219303455809 Exc = -217.0451210533288 Eec = -23632.6162052652362 Etot = Ekin + EHart + Exc + Eec Etot = -9647.3686230992425 *********************************************************** ** DATA for factorized norm conserving pseudo potentials ** *********************************************************** j.dependent.pseudo.potentials on L=0, dif of log deris for all electrons = 138343.0449900448 138343.0449900448 L=0, dif of log deris for semi local = 133352.8436615304 133352.8436615304 L=1, dif of log deris for all electrons = 1167.4367537505 2606.3047391964 L=1, dif of log deris for semi local = 56.3853966756 350.8467724917 L=2, dif of log deris for all electrons = 14.6249332740 1.0042695592 L=2, dif of log deris for semi local = 0.0683480391 0.0043880300 L=3, dif of log deris for all electrons = 2107.2244038669 9892.5366120269 L=3, dif of log deris for semi local = 1314.6558271580 4888.9512052067 *********************************************************** ** Core electron densities for PCC ** ***********************************************************