*************************************************** Input file *************************************************** # # File Name # System.CurrrentDir ./ # default=./ System.Name K_CA13 Log.print off # ON|OFF System.UseRestartfile yes # NO|YES, default=NO System.Restartfile K_CA13 # default=null # # Calculation type # eq.type dirac # sch|sdirac|dirac calc.type VPS # ALL|VPS|PAO xc.type LDA # LDA|GGA # # Atom # AtomSpecies 19 max.occupied.N 4 total.electron 19.0 valence.electron 9.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 5 Blochl.projector.num 4 # default=1 which means KB-form local.type Polynomial # Simple|Polynomial local.part.vps 1 # default=0 local.cutoff 1.30 # 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 8.0 # default=6.0 # # Pseudo atomic orbitals # maxL.pao 3 # default=2 num.pao 15 # default=7 radial.cutoff.pao 9.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=-443.2064744025636 (Hartree) NormRD=475.4199218391971 SCF= 2 Eeigen=-442.9419593613625 (Hartree) NormRD=474.9224373145460 SCF= 3 Eeigen=-433.0447145823588 (Hartree) NormRD= 28.8633161839837 SCF= 4 Eeigen=-377.1941141408989 (Hartree) NormRD= 20.7767730816322 SCF= 5 Eeigen=-343.6432537516608 (Hartree) NormRD= 12.4355486415934 SCF= 6 Eeigen=-364.4109524439998 (Hartree) NormRD= 0.9130673053335 SCF= 7 Eeigen=-347.7781949944694 (Hartree) NormRD= 0.1335752140557 SCF= 8 Eeigen=-354.7715660421555 (Hartree) NormRD= 0.0350490803299 SCF= 9 Eeigen=-351.3786411539654 (Hartree) NormRD= 0.0072863964091 SCF= 10 Eeigen=-352.9800144708491 (Hartree) NormRD= 0.0016773931388 SCF= 11 Eeigen=-352.2218636809779 (Hartree) NormRD= 0.0003739837342 SCF= 12 Eeigen=-352.5819808791696 (Hartree) NormRD= 0.0000841835576 SCF= 13 Eeigen=-352.4115767845223 (Hartree) NormRD= 0.0000189059024 SCF= 14 Eeigen=-352.4924259921460 (Hartree) NormRD= 0.0000042460134 SCF= 15 Eeigen=-352.4541316059518 (Hartree) NormRD= 0.0000009539915 SCF= 16 Eeigen=-352.4722875850600 (Hartree) NormRD= 0.0000002142559 SCF= 17 Eeigen=-352.4636842150960 (Hartree) NormRD= 0.0000000481330 SCF= 18 Eeigen=-352.4677621691179 (Hartree) NormRD= 0.0000000108113 SCF= 19 Eeigen=-352.4658295235915 (Hartree) NormRD= 0.0000000024286 SCF= 20 Eeigen=-352.4667455196602 (Hartree) NormRD= 0.0000000005455 SCF= 21 Eeigen=-352.4663113892317 (Hartree) NormRD= 0.0000000001225 SCF= 22 Eeigen=-352.4664236331104 (Hartree) NormRD= 0.0000000000047 SCF= 23 Eeigen=-352.4664263280364 (Hartree) NormRD= 0.0000000000038 SCF= 24 Eeigen=-352.4664263734854 (Hartree) NormRD= 0.0000000000038 SCF= 25 Eeigen=-352.4664264188319 (Hartree) NormRD= 0.0000000000038 SCF= 26 Eeigen=-352.4664264640920 (Hartree) NormRD= 0.0000000000038 SCF= 27 Eeigen=-352.4664265092691 (Hartree) NormRD= 0.0000000000038 SCF= 28 Eeigen=-352.4664265543637 (Hartree) NormRD= 0.0000000000038 SCF= 29 Eeigen=-352.4664265993758 (Hartree) NormRD= 0.0000000000037 SCF= 30 Eeigen=-352.4664266443032 (Hartree) NormRD= 0.0000000000037 SCF= 31 Eeigen=-352.4664266891489 (Hartree) NormRD= 0.0000000000037 SCF= 32 Eeigen=-352.4664267339127 (Hartree) NormRD= 0.0000000000037 SCF= 33 Eeigen=-352.4664267785926 (Hartree) NormRD= 0.0000000000037 SCF= 34 Eeigen=-352.4664268231905 (Hartree) NormRD= 0.0000000000037 SCF= 35 Eeigen=-352.4664268677064 (Hartree) NormRD= 0.0000000000037 SCF= 36 Eeigen=-352.4664269121408 (Hartree) NormRD= 0.0000000000036 SCF= 37 Eeigen=-352.4664269564932 (Hartree) NormRD= 0.0000000000036 SCF= 38 Eeigen=-352.4664270007634 (Hartree) NormRD= 0.0000000000036 SCF= 39 Eeigen=-352.4664270449519 (Hartree) NormRD= 0.0000000000036 SCF= 40 Eeigen=-352.4664270890603 (Hartree) NormRD= 0.0000000000036 SCF= 41 Eeigen=-352.4664271330859 (Hartree) NormRD= 0.0000000000036 SCF= 42 Eeigen=-352.4664271770314 (Hartree) NormRD= 0.0000000000036 SCF= 43 Eeigen=-352.4664272208958 (Hartree) NormRD= 0.0000000000036 SCF= 44 Eeigen=-352.4664272646792 (Hartree) NormRD= 0.0000000000035 SCF= 45 Eeigen=-352.4664273083812 (Hartree) NormRD= 0.0000000000035 SCF= 46 Eeigen=-352.4664273520034 (Hartree) NormRD= 0.0000000000035 SCF= 47 Eeigen=-352.4664273955455 (Hartree) NormRD= 0.0000000000035 SCF= 48 Eeigen=-352.4664274390064 (Hartree) NormRD= 0.0000000000035 SCF= 49 Eeigen=-352.4664274823885 (Hartree) NormRD= 0.0000000000035 SCF= 50 Eeigen=-352.4664275256893 (Hartree) NormRD= 0.0000000000035 SCF= 51 Eeigen=-352.4664275689119 (Hartree) NormRD= 0.0000000000034 SCF= 52 Eeigen=-352.4664276120533 (Hartree) NormRD= 0.0000000000034 SCF= 53 Eeigen=-352.4664276551155 (Hartree) NormRD= 0.0000000000034 SCF= 54 Eeigen=-352.4664276980997 (Hartree) NormRD= 0.0000000000034 SCF= 55 Eeigen=-352.4664277410027 (Hartree) NormRD= 0.0000000000034 SCF= 56 Eeigen=-352.4664277838278 (Hartree) NormRD= 0.0000000000034 SCF= 57 Eeigen=-352.4664278265752 (Hartree) NormRD= 0.0000000000034 SCF= 58 Eeigen=-352.4664278692417 (Hartree) NormRD= 0.0000000000034 SCF= 59 Eeigen=-352.4664279118306 (Hartree) NormRD= 0.0000000000033 SCF= 60 Eeigen=-352.4664279543401 (Hartree) NormRD= 0.0000000000000 ***************************************************** Eigenvalues (Hartree) in the all electron calculation ***************************************************** j=l+1/2 j=l-1/2 n= 1 l= 0 -128.9885636565145 -128.9885636565145 n= 2 l= 0 -12.9350647098204 -12.9350647098204 n= 2 l= 1 -10.2455869663476 -10.3515018536315 n= 3 l= 0 -1.2911939953272 -1.2911939953272 n= 3 l= 1 -0.6877322949852 -0.6976843561944 n= 4 l= 0 -0.0892668175401 -0.0892668175401 ***************************************************** Energies (Hartree) in the all electron calculation ***************************************************** Eeigen = -352.4664279543401 Ekin = 604.5104918974652 EHart = 258.4295197497485 Exc = -31.8062096557754 Eec = -1431.7688710625475 Etot = Ekin + EHart + Exc + Eec Etot = -600.6350690711093 *********************************************************** ** DATA for factorized norm conserving pseudo potentials ** *********************************************************** j.dependent.pseudo.potentials on L=0, dif of log deris for all electrons = 249913.2321147912 249913.2321147912 L=0, dif of log deris for semi local = 266977.5256753487 266977.5256753487 L=1, dif of log deris for all electrons = 41.4070496126 307161.7783857044 L=1, dif of log deris for semi local = 2064.5870817490 1434077.5914837022 L=2, dif of log deris for all electrons = 0.1185640807 0.1217525669 L=2, dif of log deris for semi local = 0.0014858973 0.0015772501 *********************************************************** ** Core electron densities for PCC ** ***********************************************************