DBCSR| Multiplication driver SMM DBCSR| Multiplication stack size 1000 DBCSR| Multiplication stack buffers 0 DBCSR| Multiplication size stacks 3 **** **** ****** ** PROGRAM STARTED AT 2013-12-03 19:16:22.468 ***** ** *** *** ** PROGRAM STARTED ON spring.kcl.ac.uk ** **** ****** PROGRAM STARTED BY tong ***** ** ** ** ** PROGRAM PROCESS ID 24417 **** ** ******* ** PROGRAM STARTED IN /home/tong/tutorials/geometry_optimis ation/example_1 CP2K| version string: CP2K version 2.5 (Development Version) CP2K| source code revision number: git:071de34 CP2K| is freely available from http://www.cp2k.org/ CP2K| Program compiled at Thu Oct 24 15:34:30 BST 2013 CP2K| Program compiled on spring.kcl.ac.uk CP2K| Program compiled for Linux-spring-gfortran CP2K| Input file name H2O.inp GLOBAL| Force Environment number 1 GLOBAL| Basis set file name ./BASIS_SET GLOBAL| Geminal file name BASIS_GEMINAL GLOBAL| Potential file name ./POTENTIAL GLOBAL| MM Potential file name MM_POTENTIAL GLOBAL| Coordinate file name __STD_INPUT__ GLOBAL| Method name CP2K GLOBAL| Project name H2O GLOBAL| Preferred FFT library FFTW3 GLOBAL| Run type GEO_OPT GLOBAL| All-to-all communication in single precision F GLOBAL| FFTs using library dependent lengths F GLOBAL| Global print level LOW GLOBAL| Total number of message passing processes 1 GLOBAL| Number of threads for this process 1 GLOBAL| This output is from process 0 MEMORY| system memory details [Kb] MEMORY| rank 0 min max average MEMORY| MemTotal 65362240 65362240 65362240 65362240 MEMORY| MemFree 63578732 63578732 63578732 63578732 MEMORY| Buffers 430200 430200 430200 430200 MEMORY| Cached 866244 866244 866244 866244 MEMORY| Slab 303724 303724 303724 303724 MEMORY| SReclaimable 0 0 0 0 MEMORY| MemLikelyFree 64875176 64875176 64875176 64875176 GENERATE| Preliminary Number of Bonds generated: 0 GENERATE| Achieved consistency in connectivity generation. ******************************************************************************* ******************************************************************************* ** ** ** ##### ## ## ** ** ## ## ## ## ## ** ** ## ## ## ###### ** ** ## ## ## ## ## ##### ## ## #### ## ##### ##### ** ** ## ## ## ## ## ## ## ## ## ## ## ## ## ## ** ** ## ## ## ## ## ## ## #### ### ## ###### ###### ** ** ## ### ## ## ## ## ## ## ## ## ## ## ** ** ####### ##### ## ##### ## ## #### ## ##### ## ** ** ## ## ** ** ** ** ... make the atoms dance ** ** ** ** Copyright (C) by CP2K Developers Group (2000 - 2013) ** ** ** ******************************************************************************* SCF PARAMETERS Density guess: ATOMIC -------------------------------------------------------- max_scf: 200 max_scf_history: 0 max_diis: 4 -------------------------------------------------------- eps_scf: 1.00E-05 eps_scf_history: 0.00E+00 eps_diis: 1.00E-01 eps_eigval: 1.00E-05 -------------------------------------------------------- level_shift [a.u.]: 0.00 -------------------------------------------------------- Mixing method: PULAY_MIXING charge density mixing in g-space -------------------------------------------------------- No outer SCF ******************************************************************************* *** STARTING GEOMETRY OPTIMIZATION *** *** CONJUGATE GRADIENTS *** ******************************************************************************* Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Extrapolation method: initial_guess SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 NoMix/Diag. 0.50E+00 0.6 1.14126074 -17.0091284306 -1.70E+01 2 Pulay/Diag. 0.50E+00 1.1 0.37692421 -15.9319529123 1.08E+00 3 Pulay/Diag. 0.50E+00 1.1 0.08104039 -17.4315722690 -1.50E+00 4 Pulay/Diag. 0.50E+00 1.1 0.01022197 -17.1557819144 2.76E-01 5 Pulay/Diag. 0.50E+00 1.1 0.00392885 -17.1916748347 -3.59E-02 6 Pulay/Diag. 0.50E+00 1.1 0.00130980 -17.1693054239 2.24E-02 7 Pulay/Diag. 0.50E+00 1.1 0.00059051 -17.1649261163 4.38E-03 8 Pulay/Diag. 0.50E+00 1.1 0.00056162 -17.1645402028 3.86E-04 9 Pulay/Diag. 0.50E+00 1.1 0.00020889 -17.1637669166 7.73E-04 10 Pulay/Diag. 0.50E+00 1.1 0.00008745 -17.1638478825 -8.10E-05 11 Pulay/Diag. 0.50E+00 1.1 0.00003405 -17.1636663285 1.82E-04 12 Pulay/Diag. 0.50E+00 1.1 0.00001742 -17.1637269776 -6.06E-05 13 Pulay/Diag. 0.50E+00 1.1 0.00000474 -17.1636724752 5.45E-05 *** SCF run converged in 13 steps *** Electronic density on regular grids: -8.0000016694 -0.0000016694 Core density on regular grids: 7.9999999260 -0.0000000740 Total charge density on r-space grids: -0.0000017433 Total charge density g-space grids: -0.0000017433 Overlap energy of the core charge distribution: 0.00000008790619 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.88082901797401 Hartree energy: 17.93243606773645 Exchange-correlation energy: -4.14404710293828 Total energy: -17.16367247523647 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.163667082526779 -------- Informations at step = 0 ------------ Optimization Method = SD Total Energy = -17.1636670825 Used time = 16.051 --------------------------------------------------- -------------------------- OPTIMIZATION STEP: 1 -------------------------- Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 0 B(1) = 1.000000 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.00000319 -17.1636670825 -1.72E+01 *** SCF run converged in 1 steps *** Electronic density on regular grids: -8.0000016694 -0.0000016694 Core density on regular grids: 7.9999999260 -0.0000000740 Total charge density on r-space grids: -0.0000017433 Total charge density g-space grids: -0.0000017433 Overlap energy of the core charge distribution: 0.00000008790619 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.88083641142988 Hartree energy: 17.93243528069200 Exchange-correlation energy: -4.14404831664001 Total energy: -17.16366708252677 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.163667082526775 Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 0 B(1) = 2.000000 B(2) = -1.000000 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.09314893 -17.1596411102 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.0 0.05242231 -17.3364934962 -1.77E-01 3 Pulay/Diag. 0.50E+00 1.1 0.01721066 -17.1071065770 2.29E-01 4 Pulay/Diag. 0.50E+00 1.1 0.00188574 -17.1639028607 -5.68E-02 5 Pulay/Diag. 0.50E+00 1.1 0.00038991 -17.1591802039 4.72E-03 6 Pulay/Diag. 0.50E+00 1.1 0.00015544 -17.1608426971 -1.66E-03 7 Pulay/Diag. 0.50E+00 1.1 0.00004638 -17.1613983340 -5.56E-04 8 Pulay/Diag. 0.50E+00 1.2 0.00003099 -17.1614614029 -6.31E-05 9 Pulay/Diag. 0.50E+00 1.2 0.00001176 -17.1614750194 -1.36E-05 10 Pulay/Diag. 0.50E+00 1.2 0.00000549 -17.1614738322 1.19E-06 *** SCF run converged in 10 steps *** Electronic density on regular grids: -8.0000016489 -0.0000016489 Core density on regular grids: 7.9999997971 -0.0000002029 Total charge density on r-space grids: -0.0000018518 Total charge density g-space grids: -0.0000018518 Overlap energy of the core charge distribution: 0.00000000989725 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.69954411811351 Hartree energy: 18.06385524344190 Exchange-correlation energy: -4.09198265770862 Total energy: -17.16147383217079 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.161473832170792 Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 1 B(1) = 2.500000 B(2) = -2.000000 B(3) = 0.500000 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.35670575 -17.0997779604 -1.71E+01 2 Pulay/Diag. 0.50E+00 1.1 0.17256417 -16.4457222936 6.54E-01 3 Pulay/Diag. 0.50E+00 1.1 0.05554014 -17.2035105556 -7.58E-01 4 Pulay/Diag. 0.50E+00 1.1 0.01055538 -17.1068669901 9.66E-02 5 Pulay/Diag. 0.50E+00 1.2 0.00166654 -17.1409680444 -3.41E-02 6 Pulay/Diag. 0.50E+00 1.2 0.00065723 -17.1348781226 6.09E-03 7 Pulay/Diag. 0.50E+00 1.2 0.00017923 -17.1331525275 1.73E-03 8 Pulay/Diag. 0.50E+00 1.2 0.00016873 -17.1329573629 1.95E-04 9 Pulay/Diag. 0.50E+00 1.2 0.00004880 -17.1328642651 9.31E-05 10 Pulay/Diag. 0.50E+00 1.2 0.00002857 -17.1328976616 -3.34E-05 11 Pulay/Diag. 0.50E+00 1.2 0.00000781 -17.1328740965 2.36E-05 *** SCF run converged in 11 steps *** Electronic density on regular grids: -8.0000016501 -0.0000016501 Core density on regular grids: 7.9999997218 -0.0000002782 Total charge density on r-space grids: -0.0000019283 Total charge density g-space grids: -0.0000019283 Overlap energy of the core charge distribution: 0.00000234899953 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 13.24448189418320 Hartree energy: 17.69629833702023 Exchange-correlation energy: -4.24076613080539 Total energy: -17.13287409651727 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.132874096517270 ******************************************************************************* *** MNBRACK - NUMBER OF ENERGY EVALUATIONS : 3 *** ******************************************************************************* Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 2 B(1) = 2.800000 B(2) = -2.800000 B(3) = 1.200000 B(4) = -0.200000 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.81656101 -17.0612665052 -1.71E+01 2 Pulay/Diag. 0.50E+00 1.1 0.32206409 -18.3378181460 -1.28E+00 3 Pulay/Diag. 0.50E+00 1.1 0.10971166 -16.8855067378 1.45E+00 4 Pulay/Diag. 0.50E+00 1.1 0.01798009 -17.2034938506 -3.18E-01 5 Pulay/Diag. 0.50E+00 1.1 0.00242339 -17.1508777852 5.26E-02 6 Pulay/Diag. 0.50E+00 1.2 0.00081669 -17.1607007064 -9.82E-03 7 Pulay/Diag. 0.50E+00 1.2 0.00027040 -17.1633320812 -2.63E-03 8 Pulay/Diag. 0.50E+00 1.2 0.00019397 -17.1636753049 -3.43E-04 9 Pulay/Diag. 0.50E+00 1.2 0.00006461 -17.1637164178 -4.11E-05 10 Pulay/Diag. 0.50E+00 1.2 0.00002902 -17.1636878620 2.86E-05 11 Pulay/Diag. 0.50E+00 1.2 0.00001741 -17.1636771203 1.07E-05 12 Pulay/Diag. 0.50E+00 1.2 0.00000536 -17.1636746497 2.47E-06 *** SCF run converged in 12 steps *** Electronic density on regular grids: -8.0000016694 -0.0000016694 Core density on regular grids: 7.9999999260 -0.0000000740 Total charge density on r-space grids: -0.0000017433 Total charge density g-space grids: -0.0000017433 Overlap energy of the core charge distribution: 0.00000008790619 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.88081888982946 Hartree energy: 17.93244598082056 Exchange-correlation energy: -4.14404906235695 Total energy: -17.16367464971557 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.163667082526242 Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.65602547 -17.0592440578 -1.71E+01 2 Pulay/Diag. 0.50E+00 1.1 0.42094204 -15.6568085897 1.40E+00 3 Pulay/Diag. 0.50E+00 1.1 0.12091188 -17.5105859460 -1.85E+00 4 Pulay/Diag. 0.50E+00 1.1 0.01071531 -17.1551034320 3.55E-01 5 Pulay/Diag. 0.50E+00 1.2 0.00331085 -17.1821879153 -2.71E-02 6 Pulay/Diag. 0.50E+00 1.2 0.00133765 -17.1693084599 1.29E-02 7 Pulay/Diag. 0.50E+00 1.2 0.00044932 -17.1649104554 4.40E-03 8 Pulay/Diag. 0.50E+00 1.2 0.00031926 -17.1643558972 5.55E-04 9 Pulay/Diag. 0.50E+00 1.2 0.00010877 -17.1640915693 2.64E-04 10 Pulay/Diag. 0.50E+00 1.2 0.00005875 -17.1641580185 -6.64E-05 11 Pulay/Diag. 0.50E+00 1.2 0.00001449 -17.1640970467 6.10E-05 12 Pulay/Diag. 0.50E+00 1.2 0.00000563 -17.1641253297 -2.83E-05 *** SCF run converged in 12 steps *** Electronic density on regular grids: -8.0000016667 -0.0000016667 Core density on regular grids: 7.9999997526 -0.0000002474 Total charge density on r-space grids: -0.0000019141 Total charge density g-space grids: -0.0000019141 Overlap energy of the core charge distribution: 0.00000002993139 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.78640769380675 Hartree energy: 17.99968668426500 Exchange-correlation energy: -4.11732919182507 Total energy: -17.16412532973676 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164119921466458 Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.25607537 -17.1496587234 -1.71E+01 2 Pulay/Diag. 0.50E+00 1.1 0.14812337 -17.6889529824 -5.39E-01 3 Pulay/Diag. 0.50E+00 1.1 0.04632243 -17.0213739009 6.68E-01 4 Pulay/Diag. 0.50E+00 1.1 0.00488297 -17.1704574872 -1.49E-01 5 Pulay/Diag. 0.50E+00 1.1 0.00098556 -17.1590405561 1.14E-02 6 Pulay/Diag. 0.50E+00 1.2 0.00037165 -17.1629843709 -3.94E-03 7 Pulay/Diag. 0.50E+00 1.2 0.00012790 -17.1642104659 -1.23E-03 8 Pulay/Diag. 0.50E+00 1.2 0.00008048 -17.1643615861 -1.51E-04 9 Pulay/Diag. 0.50E+00 1.2 0.00002814 -17.1643652869 -3.70E-06 10 Pulay/Diag. 0.50E+00 1.2 0.00001260 -17.1643516040 1.37E-05 11 Pulay/Diag. 0.50E+00 1.2 0.00000609 -17.1643466946 4.91E-06 *** SCF run converged in 11 steps *** Electronic density on regular grids: -8.0000016641 -0.0000016641 Core density on regular grids: 7.9999997913 -0.0000002087 Total charge density on r-space grids: -0.0000018728 Total charge density g-space grids: -0.0000018728 Overlap energy of the core charge distribution: 0.00000004390668 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.81880460581665 Hartree energy: 17.97632913086666 Exchange-correlation energy: -4.12658992930718 Total energy: -17.16434669463202 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164344232801547 Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.01620034 -17.1643144601 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.00671754 -17.1421863988 2.21E-02 3 Pulay/Diag. 0.50E+00 1.1 0.00153561 -17.1709307686 -2.87E-02 4 Pulay/Diag. 0.50E+00 1.1 0.00006308 -17.1640374947 6.89E-03 5 Pulay/Diag. 0.50E+00 1.2 0.00001385 -17.1642837612 -2.46E-04 6 Pulay/Diag. 0.50E+00 1.2 0.00000890 -17.1643184929 -3.47E-05 *** SCF run converged in 6 steps *** Electronic density on regular grids: -8.0000016642 -0.0000016642 Core density on regular grids: 7.9999997958 -0.0000002042 Total charge density on r-space grids: -0.0000018684 Total charge density g-space grids: -0.0000018684 Overlap energy of the core charge distribution: 0.00000004471872 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.82041060723368 Hartree energy: 17.97519871941514 Exchange-correlation energy: -4.12703731838550 Total energy: -17.16431849293280 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164344750827542 Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.03406702 -17.1640999970 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.01669989 -17.1042489171 5.99E-02 3 Pulay/Diag. 0.50E+00 1.1 0.00562704 -17.1789823745 -7.47E-02 4 Pulay/Diag. 0.50E+00 1.1 0.00080420 -17.1628295901 1.62E-02 5 Pulay/Diag. 0.50E+00 1.2 0.00014566 -17.1651708651 -2.34E-03 6 Pulay/Diag. 0.50E+00 1.2 0.00005877 -17.1645696322 6.01E-04 7 Pulay/Diag. 0.50E+00 1.2 0.00001630 -17.1643757592 1.94E-04 8 Pulay/Diag. 0.50E+00 1.2 0.00001289 -17.1643537465 2.20E-05 9 Pulay/Diag. 0.50E+00 1.2 0.00000429 -17.1643432927 1.05E-05 *** SCF run converged in 9 steps *** Electronic density on regular grids: -8.0000016643 -0.0000016643 Core density on regular grids: 7.9999997977 -0.0000002023 Total charge density on r-space grids: -0.0000018665 Total charge density g-space grids: -0.0000018665 Overlap energy of the core charge distribution: 0.00000004502459 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.82098050186069 Hartree energy: 17.97477184026135 Exchange-correlation energy: -4.12720513395545 Total energy: -17.16434329272366 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164344686958589 ******************************************************************************* *** BRENT - NUMBER OF ENERGY EVALUATIONS : 5 *** ******************************************************************************* -------- Informations at step = 1 ------------ Optimization Method = SD Total Energy = -17.1643447508 Real energy change = -0.0006776683 Decrease in energy = YES Used time = 90.837 Convergence check : Max. step size = 0.0336570168 Conv. limit for step size = 0.0010000000 Convergence in step size = NO RMS step size = 0.0168136889 Conv. limit for RMS step = 0.0010000000 Convergence in RMS step = NO Max. gradient = 0.0182785685 Conv. limit for gradients = 0.0010000000 Conv. for gradients = NO RMS gradient = 0.0091312361 Conv. limit for RMS grad. = 0.0010000000 Conv. for gradients = NO --------------------------------------------------- Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.01125367 -17.1643176223 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.00576001 -17.1849798946 -2.07E-02 3 Pulay/Diag. 0.50E+00 1.1 0.00190123 -17.1591895975 2.58E-02 4 Pulay/Diag. 0.50E+00 1.1 0.00025408 -17.1647958475 -5.61E-03 5 Pulay/Diag. 0.50E+00 1.1 0.00004687 -17.1640832862 7.13E-04 6 Pulay/Diag. 0.50E+00 1.2 0.00001866 -17.1642752777 -1.92E-04 7 Pulay/Diag. 0.50E+00 1.2 0.00000508 -17.1643364606 -6.12E-05 *** SCF run converged in 7 steps *** Electronic density on regular grids: -8.0000016642 -0.0000016642 Core density on regular grids: 7.9999997958 -0.0000002042 Total charge density on r-space grids: -0.0000018684 Total charge density g-space grids: -0.0000018684 Overlap energy of the core charge distribution: 0.00000004471872 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.82039054547135 Hartree energy: 17.97520433100699 Exchange-correlation energy: -4.12704083590046 Total energy: -17.16433646061824 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164344750825521 -------------------------- OPTIMIZATION STEP: 2 -------------------------- Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.00178634 -17.1643440672 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.00091688 -17.1610761917 3.27E-03 3 Pulay/Diag. 0.50E+00 1.1 0.00029574 -17.1651676270 -4.09E-03 4 Pulay/Diag. 0.50E+00 1.1 0.00003862 -17.1642796677 8.88E-04 5 Pulay/Diag. 0.50E+00 1.2 0.00000742 -17.1643874695 -1.08E-04 *** SCF run converged in 5 steps *** Electronic density on regular grids: -8.0000016642 -0.0000016642 Core density on regular grids: 7.9999997958 -0.0000002042 Total charge density on r-space grids: -0.0000018684 Total charge density g-space grids: -0.0000018684 Overlap energy of the core charge distribution: 0.00000004471872 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.82035479376651 Hartree energy: 17.97518460709966 Exchange-correlation energy: -4.12703636918320 Total energy: -17.16438746951314 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164387469513144 Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.01883382 -17.1642290770 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.00144760 -17.1634885418 7.41E-04 3 Pulay/Diag. 0.50E+00 1.1 0.00076880 -17.1633552375 1.33E-04 4 Pulay/Diag. 0.50E+00 1.1 0.00065799 -17.1618509856 1.50E-03 5 Pulay/Diag. 0.50E+00 1.2 0.00010833 -17.1647392637 -2.89E-03 6 Pulay/Diag. 0.50E+00 1.2 0.00001867 -17.1643040430 4.35E-04 7 Pulay/Diag. 0.50E+00 1.2 0.00000510 -17.1642924469 1.16E-05 *** SCF run converged in 7 steps *** Electronic density on regular grids: -8.0000016383 -0.0000016383 Core density on regular grids: 7.9999996330 -0.0000003670 Total charge density on r-space grids: -0.0000020054 Total charge density g-space grids: -0.0000020054 Overlap energy of the core charge distribution: 0.00000004717883 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.82261556545921 Hartree energy: 17.97352730351612 Exchange-correlation energy: -4.12754481712457 Total energy: -17.16429244688525 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164292446885248 Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.08748390 -17.1626675359 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.00869882 -17.1744207651 -1.18E-02 3 Pulay/Diag. 0.50E+00 1.1 0.00492822 -17.1688112015 5.61E-03 4 Pulay/Diag. 0.50E+00 1.1 0.00406450 -17.1801885582 -1.14E-02 5 Pulay/Diag. 0.50E+00 1.2 0.00040090 -17.1624285631 1.78E-02 6 Pulay/Diag. 0.50E+00 1.2 0.00009738 -17.1637709580 -1.34E-03 7 Pulay/Diag. 0.50E+00 1.2 0.00002626 -17.1638620215 -9.11E-05 8 Pulay/Diag. 0.50E+00 1.2 0.00001440 -17.1638886655 -2.66E-05 9 Pulay/Diag. 0.50E+00 1.2 0.00000574 -17.1638927611 -4.10E-06 *** SCF run converged in 9 steps *** Electronic density on regular grids: -8.0000017083 -0.0000017083 Core density on regular grids: 8.0000001303 0.0000001303 Total charge density on r-space grids: -0.0000015780 Total charge density g-space grids: -0.0000015780 Overlap energy of the core charge distribution: 0.00000004628670 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.81660493871823 Hartree energy: 17.97844033432607 Exchange-correlation energy: -4.12604753450010 Total energy: -17.16389276108394 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.163892761083936 ******************************************************************************* *** MNBRACK - NUMBER OF ENERGY EVALUATIONS : 3 *** ******************************************************************************* Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.15723552 -17.1604095236 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.01562662 -17.1432337173 1.72E-02 3 Pulay/Diag. 0.50E+00 1.1 0.00898724 -17.1560328507 -1.28E-02 4 Pulay/Diag. 0.50E+00 1.1 0.00721788 -17.1350673719 2.10E-02 5 Pulay/Diag. 0.50E+00 1.2 0.00066790 -17.1668046428 -3.17E-02 6 Pulay/Diag. 0.50E+00 1.2 0.00011713 -17.1646401959 2.16E-03 7 Pulay/Diag. 0.50E+00 1.2 0.00003686 -17.1644159587 2.24E-04 8 Pulay/Diag. 0.50E+00 1.2 0.00002621 -17.1643583423 5.76E-05 9 Pulay/Diag. 0.50E+00 1.2 0.00000848 -17.1643494315 8.91E-06 *** SCF run converged in 9 steps *** Electronic density on regular grids: -8.0000016642 -0.0000016642 Core density on regular grids: 7.9999997958 -0.0000002042 Total charge density on r-space grids: -0.0000018684 Total charge density g-space grids: -0.0000018684 Overlap energy of the core charge distribution: 0.00000004471872 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.82038529157481 Hartree energy: 17.97519571173878 Exchange-correlation energy: -4.12703993359810 Total energy: -17.16434943148063 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164344750795269 Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.13386955 -17.1615670512 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.00966117 -17.1737748910 -1.22E-02 3 Pulay/Diag. 0.50E+00 1.1 0.00594804 -17.1716119666 2.16E-03 4 Pulay/Diag. 0.50E+00 1.1 0.00501011 -17.1839350126 -1.23E-02 5 Pulay/Diag. 0.50E+00 1.2 0.00065728 -17.1617524052 2.22E-02 6 Pulay/Diag. 0.50E+00 1.2 0.00015268 -17.1642248403 -2.47E-03 7 Pulay/Diag. 0.50E+00 1.2 0.00004198 -17.1643147907 -9.00E-05 8 Pulay/Diag. 0.50E+00 1.2 0.00002149 -17.1643473827 -3.26E-05 9 Pulay/Diag. 0.50E+00 1.2 0.00000818 -17.1643506375 -3.25E-06 *** SCF run converged in 9 steps *** Electronic density on regular grids: -8.0000016504 -0.0000016504 Core density on regular grids: 7.9999997078 -0.0000002922 Total charge density on r-space grids: -0.0000019427 Total charge density g-space grids: -0.0000019427 Overlap energy of the core charge distribution: 0.00000004558986 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.82151176777913 Hartree energy: 17.97433028642683 Exchange-correlation energy: -4.12730219139964 Total energy: -17.16435063751867 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164352299703623 Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.04824996 -17.1639997701 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.00397166 -17.1591825776 4.82E-03 3 Pulay/Diag. 0.50E+00 1.1 0.00250056 -17.1616821905 -2.50E-03 4 Pulay/Diag. 0.50E+00 1.1 0.00204121 -17.1560800792 5.60E-03 5 Pulay/Diag. 0.50E+00 1.2 0.00020231 -17.1651227524 -9.04E-03 6 Pulay/Diag. 0.50E+00 1.2 0.00004162 -17.1644354739 6.87E-04 7 Pulay/Diag. 0.50E+00 1.2 0.00001091 -17.1643772083 5.83E-05 8 Pulay/Diag. 0.50E+00 1.2 0.00000793 -17.1643612339 1.60E-05 *** SCF run converged in 8 steps *** Electronic density on regular grids: -8.0000016556 -0.0000016556 Core density on regular grids: 7.9999997426 -0.0000002574 Total charge density on r-space grids: -0.0000019130 Total charge density g-space grids: -0.0000019130 Overlap energy of the core charge distribution: 0.00000004515717 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.82106224146502 Hartree energy: 17.97466564791540 Exchange-correlation energy: -4.12719862255560 Total energy: -17.16436123393285 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164357819318919 Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.00294651 -17.1643568951 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.00035415 -17.1641481903 2.09E-04 3 Pulay/Diag. 0.50E+00 1.1 0.00015403 -17.1642082488 -6.01E-05 4 Pulay/Diag. 0.50E+00 1.1 0.00012071 -17.1639031901 3.05E-04 5 Pulay/Diag. 0.50E+00 1.2 0.00001802 -17.1644264066 -5.23E-04 6 Pulay/Diag. 0.50E+00 1.2 0.00000192 -17.1643591881 6.72E-05 *** SCF run converged in 6 steps *** Electronic density on regular grids: -8.0000016555 -0.0000016555 Core density on regular grids: 7.9999997417 -0.0000002583 Total charge density on r-space grids: -0.0000019138 Total charge density g-space grids: -0.0000019138 Overlap energy of the core charge distribution: 0.00000004516714 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.82107323008925 Hartree energy: 17.97465868239386 Exchange-correlation energy: -4.12720059978790 Total energy: -17.16435918805249 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164357823356635 Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.00802053 -17.1643473607 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.00075089 -17.1651422053 -7.95E-04 3 Pulay/Diag. 0.50E+00 1.1 0.00042236 -17.1648297687 3.12E-04 4 Pulay/Diag. 0.50E+00 1.1 0.00035232 -17.1657594648 -9.30E-04 5 Pulay/Diag. 0.50E+00 1.1 0.00004055 -17.1642053519 1.55E-03 6 Pulay/Diag. 0.50E+00 1.2 0.00000852 -17.1643479116 -1.43E-04 *** SCF run converged in 6 steps *** Electronic density on regular grids: -8.0000016554 -0.0000016554 Core density on regular grids: 7.9999997412 -0.0000002588 Total charge density on r-space grids: -0.0000019143 Total charge density g-space grids: -0.0000019143 Overlap energy of the core charge distribution: 0.00000004517292 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.82108857879664 Hartree energy: 17.97465659363417 Exchange-correlation energy: -4.12720258330051 Total energy: -17.16434791161162 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164357822016584 ******************************************************************************* *** BRENT - NUMBER OF ENERGY EVALUATIONS : 5 *** ******************************************************************************* -------- Informations at step = 2 ------------ Optimization Method = CG Total Energy = -17.1643578234 Real energy change = -0.0000130725 Decrease in energy = YES Used time = 84.833 Convergence check : Max. step size = 0.0084583772 Conv. limit for step size = 0.0010000000 Convergence in step size = NO RMS step size = 0.0041001011 Conv. limit for RMS step = 0.0010000000 Convergence in RMS step = NO Max. gradient = 0.0014597313 Conv. limit for gradients = 0.0010000000 Conv. for gradients = NO RMS gradient = 0.0007075880 Conv. limit for RMS grad. = 0.0010000000 Conv. in RMS gradients = YES --------------------------------------------------- Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.00252682 -17.1643567924 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.00022895 -17.1641115745 2.45E-04 3 Pulay/Diag. 0.50E+00 1.1 0.00013301 -17.1642100229 -9.84E-05 4 Pulay/Diag. 0.50E+00 1.1 0.00011032 -17.1639166656 2.93E-04 5 Pulay/Diag. 0.50E+00 1.1 0.00001229 -17.1644042656 -4.88E-04 6 Pulay/Diag. 0.50E+00 1.2 0.00000255 -17.1643612620 4.30E-05 *** SCF run converged in 6 steps *** Electronic density on regular grids: -8.0000016555 -0.0000016555 Core density on regular grids: 7.9999997417 -0.0000002583 Total charge density on r-space grids: -0.0000019138 Total charge density g-space grids: -0.0000019138 Overlap energy of the core charge distribution: 0.00000004516714 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.82107156098274 Hartree energy: 17.97465818333018 Exchange-correlation energy: -4.12720050558456 Total energy: -17.16436126201933 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164357823356244 -------------------------- OPTIMIZATION STEP: 3 -------------------------- Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.00011612 -17.1643578208 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.0 0.00000902 -17.1643628545 -5.03E-06 *** SCF run converged in 2 steps *** Electronic density on regular grids: -8.0000016554 -0.0000016554 Core density on regular grids: 7.9999997417 -0.0000002583 Total charge density on r-space grids: -0.0000019137 Total charge density g-space grids: -0.0000019137 Overlap energy of the core charge distribution: 0.00000004516714 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.82105997632089 Hartree energy: 17.97466712860896 Exchange-correlation energy: -4.12719945867212 Total energy: -17.16436285448997 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164362854489966 Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.00461243 -17.1643678820 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.00040361 -17.1635722282 7.96E-04 3 Pulay/Diag. 0.50E+00 1.1 0.00019446 -17.1643687561 -7.97E-04 4 Pulay/Diag. 0.50E+00 1.1 0.00011950 -17.1638664436 5.02E-04 5 Pulay/Diag. 0.50E+00 1.1 0.00000776 -17.1643963231 -5.30E-04 *** SCF run converged in 5 steps *** Electronic density on regular grids: -8.0000016539 -0.0000016539 Core density on regular grids: 7.9999997123 -0.0000002877 Total charge density on r-space grids: -0.0000019416 Total charge density g-space grids: -0.0000019416 Overlap energy of the core charge distribution: 0.00000004552175 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.82200131109398 Hartree energy: 17.97394011312336 Exchange-correlation energy: -4.12744724697194 Total energy: -17.16439632314768 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164396323147677 Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.00173076 -17.1643777394 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.0 0.00015096 -17.1646606755 -2.83E-04 3 Pulay/Diag. 0.50E+00 1.1 0.00008404 -17.1644101813 2.50E-04 4 Pulay/Diag. 0.50E+00 1.1 0.00005342 -17.1646043347 -1.94E-04 5 Pulay/Diag. 0.50E+00 1.1 0.00000379 -17.1643664682 2.38E-04 *** SCF run converged in 5 steps *** Electronic density on regular grids: -8.0000016515 -0.0000016515 Core density on regular grids: 7.9999996759 -0.0000003241 Total charge density on r-space grids: -0.0000019756 Total charge density g-space grids: -0.0000019756 Overlap energy of the core charge distribution: 0.00000004610784 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.82348828553160 Hartree energy: 17.97285991155271 Exchange-correlation energy: -4.12782416552370 Total energy: -17.16436646824639 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164366468246385 ******************************************************************************* *** MNBRACK - NUMBER OF ENERGY EVALUATIONS : 3 *** ******************************************************************************* Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.01629198 -17.1643381474 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.0 0.00134289 -17.1671420505 -2.80E-03 3 Pulay/Diag. 0.50E+00 1.1 0.00065009 -17.1643133656 2.83E-03 4 Pulay/Diag. 0.50E+00 1.1 0.00038072 -17.1660153805 -1.70E-03 5 Pulay/Diag. 0.50E+00 1.1 0.00002188 -17.1642923505 1.72E-03 6 Pulay/Diag. 0.50E+00 1.1 0.00000626 -17.1643482650 -5.59E-05 *** SCF run converged in 6 steps *** Electronic density on regular grids: -8.0000016539 -0.0000016539 Core density on regular grids: 7.9999997123 -0.0000002877 Total charge density on r-space grids: -0.0000019416 Total charge density g-space grids: -0.0000019416 Overlap energy of the core charge distribution: 0.00000004552175 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.82203800046740 Hartree energy: 17.97395300977982 Exchange-correlation energy: -4.12744877482764 Total energy: -17.16434826497350 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164370582764310 Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.02780149 -17.1642809191 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.0 0.00230834 -17.1596403680 4.64E-03 3 Pulay/Diag. 0.50E+00 1.1 0.00111811 -17.1644751494 -4.83E-03 4 Pulay/Diag. 0.50E+00 1.1 0.00065338 -17.1615545701 2.92E-03 5 Pulay/Diag. 0.50E+00 1.1 0.00003859 -17.1645166665 -2.96E-03 6 Pulay/Diag. 0.50E+00 1.1 0.00001113 -17.1644175312 9.91E-05 7 Pulay/Diag. 0.50E+00 1.1 0.00000330 -17.1643877160 2.98E-05 *** SCF run converged in 7 steps *** Electronic density on regular grids: -8.0000016528 -0.0000016528 Core density on regular grids: 7.9999996925 -0.0000003075 Total charge density on r-space grids: -0.0000019603 Total charge density g-space grids: -0.0000019603 Overlap energy of the core charge distribution: 0.00000004581312 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.82275441182310 Hartree energy: 17.97338766253627 Exchange-correlation energy: -4.12763929025080 Total energy: -17.16438771599315 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164375831954104 Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.01978728 -17.1643284892 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.0 0.00160598 -17.1675890749 -3.26E-03 3 Pulay/Diag. 0.50E+00 1.1 0.00080388 -17.1644326190 3.16E-03 4 Pulay/Diag. 0.50E+00 1.1 0.00048815 -17.1664926165 -2.06E-03 5 Pulay/Diag. 0.50E+00 1.1 0.00002913 -17.1642774325 2.22E-03 6 Pulay/Diag. 0.50E+00 1.1 0.00000752 -17.1643489054 -7.15E-05 *** SCF run converged in 6 steps *** Electronic density on regular grids: -8.0000016522 -0.0000016522 Core density on regular grids: 7.9999996800 -0.0000003200 Total charge density on r-space grids: -0.0000019722 Total charge density g-space grids: -0.0000019722 Overlap energy of the core charge distribution: 0.00000004599331 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.82322430615833 Hartree energy: 17.97307172966863 Exchange-correlation energy: -4.12775444129877 Total energy: -17.16434890539334 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164376681917446 Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.00504734 -17.1643735218 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.0 0.00041563 -17.1635406532 8.33E-04 3 Pulay/Diag. 0.50E+00 1.1 0.00020829 -17.1643774819 -8.37E-04 4 Pulay/Diag. 0.50E+00 1.1 0.00012226 -17.1638409337 5.37E-04 5 Pulay/Diag. 0.50E+00 1.1 0.00000678 -17.1644032478 -5.62E-04 *** SCF run converged in 5 steps *** Electronic density on regular grids: -8.0000016522 -0.0000016522 Core density on regular grids: 7.9999996804 -0.0000003196 Total charge density on r-space grids: -0.0000019718 Total charge density g-space grids: -0.0000019718 Overlap energy of the core charge distribution: 0.00000004598490 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.82316244295635 Hartree energy: 17.97307217501266 Exchange-correlation energy: -4.12774736584280 Total energy: -17.16440324780371 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164376683335345 Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.00100721 -17.1643765587 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.00008840 -17.1641886021 1.88E-04 3 Pulay/Diag. 0.50E+00 1.1 0.00005032 -17.1643788289 -1.90E-04 4 Pulay/Diag. 0.50E+00 1.1 0.00002576 -17.1642775891 1.01E-04 5 Pulay/Diag. 0.50E+00 1.2 0.00000235 -17.1643816868 -1.04E-04 *** SCF run converged in 5 steps *** Electronic density on regular grids: -8.0000016522 -0.0000016522 Core density on regular grids: 7.9999996808 -0.0000003192 Total charge density on r-space grids: -0.0000019714 Total charge density g-space grids: -0.0000019714 Overlap energy of the core charge distribution: 0.00000004597658 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.82315869207691 Hartree energy: 17.97309245368851 Exchange-correlation energy: -4.12774233258422 Total energy: -17.16438168675705 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164376680781903 ******************************************************************************* *** BRENT - NUMBER OF ENERGY EVALUATIONS : 5 *** ******************************************************************************* -------- Informations at step = 3 ------------ Optimization Method = CG Total Energy = -17.1643766833 Real energy change = -0.0000188600 Decrease in energy = YES Used time = 61.345 Convergence check : Max. step size = 0.0190351391 Conv. limit for step size = 0.0010000000 Convergence in step size = NO RMS step size = 0.0074837796 Conv. limit for RMS step = 0.0010000000 Convergence in RMS step = NO Max. gradient = 0.0033936444 Conv. limit for gradients = 0.0010000000 Conv. for gradients = NO RMS gradient = 0.0013342318 Conv. limit for RMS grad. = 0.0010000000 Conv. for gradients = NO --------------------------------------------------- Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.00073233 -17.1643766177 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.00006716 -17.1645204872 -1.44E-04 3 Pulay/Diag. 0.50E+00 1.1 0.00003656 -17.1643722406 1.48E-04 4 Pulay/Diag. 0.50E+00 1.1 0.00001889 -17.1644484783 -7.62E-05 5 Pulay/Diag. 0.50E+00 1.1 0.00000172 -17.1643727528 7.57E-05 *** SCF run converged in 5 steps *** Electronic density on regular grids: -8.0000016522 -0.0000016522 Core density on regular grids: 7.9999996804 -0.0000003196 Total charge density on r-space grids: -0.0000019718 Total charge density g-space grids: -0.0000019718 Overlap energy of the core charge distribution: 0.00000004598490 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.82318618339517 Hartree energy: 17.97307930185148 Exchange-correlation energy: -4.12774773815796 Total energy: -17.16437275284123 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164376683342120 -------------------------- OPTIMIZATION STEP: 4 -------------------------- Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.00015035 -17.1643766806 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.00000948 -17.1643926433 -1.60E-05 *** SCF run converged in 2 steps *** Electronic density on regular grids: -8.0000016522 -0.0000016522 Core density on regular grids: 7.9999996804 -0.0000003196 Total charge density on r-space grids: -0.0000019718 Total charge density g-space grids: -0.0000019718 Overlap energy of the core charge distribution: 0.00000004598490 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.82315711332945 Hartree energy: 17.97308738875923 Exchange-correlation energy: -4.12774664546194 Total energy: -17.16439264330320 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164392643303199 Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.00979991 -17.1643911150 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.00049919 -17.1654866309 -1.10E-03 3 Pulay/Diag. 0.50E+00 1.1 0.00034571 -17.1646724077 8.14E-04 4 Pulay/Diag. 0.50E+00 1.1 0.00024874 -17.1655127411 -8.40E-04 5 Pulay/Diag. 0.50E+00 1.2 0.00001016 -17.1643728152 1.14E-03 6 Pulay/Diag. 0.50E+00 1.2 0.00000259 -17.1643949037 -2.21E-05 *** SCF run converged in 6 steps *** Electronic density on regular grids: -8.0000016550 -0.0000016550 Core density on regular grids: 7.9999996667 -0.0000003333 Total charge density on r-space grids: -0.0000019883 Total charge density g-space grids: -0.0000019883 Overlap energy of the core charge distribution: 0.00000004539466 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.82201079961986 Hartree energy: 17.97392358413728 Exchange-correlation energy: -4.12743878692004 Total energy: -17.16439490368307 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164394903683075 Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.00371027 -17.1644326905 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.00020933 -17.1639758883 4.57E-04 3 Pulay/Diag. 0.50E+00 1.1 0.00011941 -17.1644147929 -4.39E-04 4 Pulay/Diag. 0.50E+00 1.1 0.00007827 -17.1640739371 3.41E-04 5 Pulay/Diag. 0.50E+00 1.2 0.00000191 -17.1644436306 -3.70E-04 *** SCF run converged in 5 steps *** Electronic density on regular grids: -8.0000016599 -0.0000016599 Core density on regular grids: 7.9999996190 -0.0000003810 Total charge density on r-space grids: -0.0000020409 Total charge density g-space grids: -0.0000020409 Overlap energy of the core charge distribution: 0.00000004440458 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.81981447354642 Hartree energy: 17.97548633372173 Exchange-correlation energy: -4.12685393631916 Total energy: -17.16444363056126 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164443630561259 Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.00803057 -17.1644481588 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.00049761 -17.1658014958 -1.35E-03 3 Pulay/Diag. 0.50E+00 1.1 0.00033142 -17.1645240943 1.28E-03 4 Pulay/Diag. 0.50E+00 1.1 0.00019375 -17.1652916455 -7.68E-04 5 Pulay/Diag. 0.50E+00 1.2 0.00001080 -17.1644244520 8.67E-04 6 Pulay/Diag. 0.50E+00 1.2 0.00000248 -17.1644467347 -2.23E-05 *** SCF run converged in 6 steps *** Electronic density on regular grids: -8.0000016674 -0.0000016674 Core density on regular grids: 7.9999995296 -0.0000004704 Total charge density on r-space grids: -0.0000021377 Total charge density g-space grids: -0.0000021377 Overlap energy of the core charge distribution: 0.00000004276685 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.81558645785949 Hartree energy: 17.97855487905384 Exchange-correlation energy: -4.12569756846444 Total energy: -17.16444673469909 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164446734699091 Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.04610639 -17.1642292369 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.00278904 -17.1572236993 7.01E-03 3 Pulay/Diag. 0.50E+00 1.1 0.00177151 -17.1640335377 -6.81E-03 4 Pulay/Diag. 0.50E+00 1.1 0.00108280 -17.1597500906 4.28E-03 5 Pulay/Diag. 0.50E+00 1.2 0.00005700 -17.1646210357 -4.87E-03 6 Pulay/Diag. 0.50E+00 1.2 0.00001423 -17.1645005942 1.20E-04 7 Pulay/Diag. 0.50E+00 1.2 0.00000338 -17.1644669537 3.36E-05 *** SCF run converged in 7 steps *** Electronic density on regular grids: -8.0000016641 -0.0000016641 Core density on regular grids: 7.9999995670 -0.0000004330 Total charge density on r-space grids: -0.0000020971 Total charge density g-space grids: -0.0000020971 Overlap energy of the core charge distribution: 0.00000004353221 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.81766053739344 Hartree energy: 17.97703236120221 Exchange-correlation energy: -4.12626934990524 Total energy: -17.16446695369222 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164466953692223 ******************************************************************************* *** MNBRACK - NUMBER OF ENERGY EVALUATIONS : 5 *** ******************************************************************************* Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.05345275 -17.1641457289 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.00308584 -17.1723254042 -8.18E-03 3 Pulay/Diag. 0.50E+00 1.1 0.00216252 -17.1653812031 6.94E-03 4 Pulay/Diag. 0.50E+00 1.1 0.00133505 -17.1702842690 -4.90E-03 5 Pulay/Diag. 0.50E+00 1.2 0.00007001 -17.1642626046 6.02E-03 6 Pulay/Diag. 0.50E+00 1.2 0.00001507 -17.1644048349 -1.42E-04 7 Pulay/Diag. 0.50E+00 1.2 0.00000399 -17.1644400604 -3.52E-05 *** SCF run converged in 7 steps *** Electronic density on regular grids: -8.0000016641 -0.0000016641 Core density on regular grids: 7.9999995670 -0.0000004330 Total charge density on r-space grids: -0.0000020971 Total charge density g-space grids: -0.0000020971 Overlap energy of the core charge distribution: 0.00000004353221 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.81767708380043 Hartree energy: 17.97704466135965 Exchange-correlation energy: -4.12627130314291 Total energy: -17.16444006036546 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164453496303871 Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.02306686 -17.1644001934 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.00141461 -17.1607402829 3.66E-03 3 Pulay/Diag. 0.50E+00 1.1 0.00085000 -17.1643344714 -3.59E-03 4 Pulay/Diag. 0.50E+00 1.1 0.00052483 -17.1621701082 2.16E-03 5 Pulay/Diag. 0.50E+00 1.2 0.00002522 -17.1645368184 -2.37E-03 6 Pulay/Diag. 0.50E+00 1.2 0.00000725 -17.1644794955 5.73E-05 *** SCF run converged in 6 steps *** Electronic density on regular grids: -8.0000016656 -0.0000016656 Core density on regular grids: 7.9999995460 -0.0000004540 Total charge density on r-space grids: -0.0000021196 Total charge density g-space grids: -0.0000021196 Overlap energy of the core charge distribution: 0.00000004314905 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.81663148804134 Hartree energy: 17.97776933753500 Exchange-correlation energy: -4.12598981831059 Total energy: -17.16447949550004 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164456475671834 Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.00288281 -17.1644556888 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.00018033 -17.1639952584 4.60E-04 3 Pulay/Diag. 0.50E+00 1.1 0.00012228 -17.1644308590 -4.36E-04 4 Pulay/Diag. 0.50E+00 1.1 0.00006691 -17.1641678286 2.63E-04 5 Pulay/Diag. 0.50E+00 1.2 0.00000436 -17.1644658094 -2.98E-04 *** SCF run converged in 5 steps *** Electronic density on regular grids: -8.0000016660 -0.0000016660 Core density on regular grids: 7.9999995447 -0.0000004553 Total charge density on r-space grids: -0.0000021213 Total charge density g-space grids: -0.0000021213 Overlap energy of the core charge distribution: 0.00000004306565 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.81641155982802 Hartree energy: 17.97794049683231 Exchange-correlation energy: -4.12592736318280 Total energy: -17.16446580937166 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164456561476300 Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.00466504 -17.1644542533 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.00028515 -17.1651964166 -7.42E-04 3 Pulay/Diag. 0.50E+00 1.1 0.00018387 -17.1644948068 7.02E-04 4 Pulay/Diag. 0.50E+00 1.1 0.00010781 -17.1649256383 -4.31E-04 5 Pulay/Diag. 0.50E+00 1.1 0.00000588 -17.1644413939 4.84E-04 *** SCF run converged in 5 steps *** Electronic density on regular grids: -8.0000016658 -0.0000016658 Core density on regular grids: 7.9999995461 -0.0000004539 Total charge density on r-space grids: -0.0000021197 Total charge density g-space grids: -0.0000021197 Overlap energy of the core charge distribution: 0.00000004309544 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.81651301668068 Hartree energy: 17.97788662396465 Exchange-correlation energy: -4.12595053173839 Total energy: -17.16444139391247 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164456555343371 ******************************************************************************* *** BRENT - NUMBER OF ENERGY EVALUATIONS : 4 *** ******************************************************************************* -------- Informations at step = 4 ------------ Optimization Method = CG Total Energy = -17.1644565615 Real energy change = -0.0000798781 Decrease in energy = YES Used time = 70.757 Convergence check : Max. step size = 0.0613507578 Conv. limit for step size = 0.0010000000 Convergence in step size = NO RMS step size = 0.0284870415 Conv. limit for RMS step = 0.0010000000 Convergence in RMS step = NO Max. gradient = 0.0013037662 Conv. limit for gradients = 0.0010000000 Conv. for gradients = NO RMS gradient = 0.0006053787 Conv. limit for RMS grad. = 0.0010000000 Conv. in RMS gradients = YES --------------------------------------------------- Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.00027241 -17.1644565534 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.00002168 -17.1645151917 -5.86E-05 3 Pulay/Diag. 0.50E+00 1.1 0.00001219 -17.1644538234 6.14E-05 4 Pulay/Diag. 0.50E+00 1.1 0.00000672 -17.1644831596 -2.93E-05 *** SCF run converged in 4 steps *** Electronic density on regular grids: -8.0000016660 -0.0000016660 Core density on regular grids: 7.9999995447 -0.0000004553 Total charge density on r-space grids: -0.0000021213 Total charge density g-space grids: -0.0000021213 Overlap energy of the core charge distribution: 0.00000004306565 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.81639147636391 Hartree energy: 17.97794347971079 Exchange-correlation energy: -4.12592761285237 Total energy: -17.16448315962686 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164456561476719 -------------------------- OPTIMIZATION STEP: 5 -------------------------- Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.00184245 -17.1644562020 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.00011743 -17.1641510686 3.05E-04 3 Pulay/Diag. 0.50E+00 1.1 0.00007371 -17.1644463359 -2.95E-04 4 Pulay/Diag. 0.50E+00 1.1 0.00004269 -17.1642727544 1.74E-04 5 Pulay/Diag. 0.50E+00 1.2 0.00000251 -17.1644633740 -1.91E-04 *** SCF run converged in 5 steps *** Electronic density on regular grids: -8.0000016660 -0.0000016660 Core density on regular grids: 7.9999995447 -0.0000004553 Total charge density on r-space grids: -0.0000021213 Total charge density g-space grids: -0.0000021213 Overlap energy of the core charge distribution: 0.00000004306565 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.81641366763618 Hartree energy: 17.97794091200880 Exchange-correlation energy: -4.12592745081641 Total energy: -17.16446337402062 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164463374020617 Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.04880213 -17.1635018277 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.00802888 -17.1373809831 2.61E-02 3 Pulay/Diag. 0.50E+00 1.1 0.00542529 -17.1653738113 -2.80E-02 4 Pulay/Diag. 0.50E+00 1.1 0.00251064 -17.1550957317 1.03E-02 5 Pulay/Diag. 0.50E+00 1.2 0.00022327 -17.1643901817 -9.29E-03 6 Pulay/Diag. 0.50E+00 1.2 0.00006044 -17.1639119840 4.78E-04 7 Pulay/Diag. 0.50E+00 1.2 0.00001955 -17.1638175802 9.44E-05 8 Pulay/Diag. 0.50E+00 1.2 0.00000905 -17.1638051718 1.24E-05 *** SCF run converged in 8 steps *** Electronic density on regular grids: -8.0000016133 -0.0000016133 Core density on regular grids: 7.9999992922 -0.0000007078 Total charge density on r-space grids: -0.0000023211 Total charge density g-space grids: -0.0000023211 Overlap energy of the core charge distribution: 0.00000007327724 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.84634899716681 Hartree energy: 17.95696837501079 Exchange-correlation energy: -4.13423207132387 Total energy: -17.16380517178386 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.163805171783864 Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.22922568 -17.1548345851 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.03626234 -17.2850130793 -1.30E-01 3 Pulay/Diag. 0.50E+00 1.1 0.02428203 -17.1462558936 1.39E-01 4 Pulay/Diag. 0.50E+00 1.1 0.01205951 -17.1987267248 -5.25E-02 5 Pulay/Diag. 0.50E+00 1.2 0.00102089 -17.1585288965 4.02E-02 6 Pulay/Diag. 0.50E+00 1.2 0.00025706 -17.1605824499 -2.05E-03 7 Pulay/Diag. 0.50E+00 1.2 0.00006785 -17.1610839560 -5.02E-04 8 Pulay/Diag. 0.50E+00 1.2 0.00003575 -17.1611412672 -5.73E-05 9 Pulay/Diag. 0.50E+00 1.2 0.00002070 -17.1611549828 -1.37E-05 10 Pulay/Diag. 0.50E+00 1.2 0.00001127 -17.1611510241 3.96E-06 11 Pulay/Diag. 0.50E+00 1.2 0.00000760 -17.1611270506 2.40E-05 *** SCF run converged in 11 steps *** Electronic density on regular grids: -8.0000017330 -0.0000017330 Core density on regular grids: 7.9999998891 -0.0000001109 Total charge density on r-space grids: -0.0000018439 Total charge density g-space grids: -0.0000018439 Overlap energy of the core charge distribution: 0.00000005637380 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.77869341649782 Hartree energy: 18.00745684865444 Exchange-correlation energy: -4.11438682620108 Total energy: -17.16112705058985 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.161127050589855 ******************************************************************************* *** MNBRACK - NUMBER OF ENERGY EVALUATIONS : 3 *** ******************************************************************************* Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.35339945 -17.1444941079 -1.71E+01 2 Pulay/Diag. 0.50E+00 1.1 0.06585056 -16.9400012412 2.04E-01 3 Pulay/Diag. 0.50E+00 1.1 0.04480844 -17.1832150890 -2.43E-01 4 Pulay/Diag. 0.50E+00 1.1 0.02055971 -17.0967547342 8.65E-02 5 Pulay/Diag. 0.50E+00 1.2 0.00184210 -17.1698930836 -7.31E-02 6 Pulay/Diag. 0.50E+00 1.2 0.00059353 -17.1657469071 4.15E-03 7 Pulay/Diag. 0.50E+00 1.2 0.00017722 -17.1645916324 1.16E-03 8 Pulay/Diag. 0.50E+00 1.2 0.00007730 -17.1644712603 1.20E-04 9 Pulay/Diag. 0.50E+00 1.2 0.00004234 -17.1644344774 3.68E-05 10 Pulay/Diag. 0.50E+00 1.2 0.00002259 -17.1644336304 8.47E-07 11 Pulay/Diag. 0.50E+00 1.2 0.00001311 -17.1644756971 -4.21E-05 12 Pulay/Diag. 0.50E+00 1.2 0.00000748 -17.1644480265 2.77E-05 *** SCF run converged in 12 steps *** Electronic density on regular grids: -8.0000016660 -0.0000016660 Core density on regular grids: 7.9999995447 -0.0000004553 Total charge density on r-space grids: -0.0000021213 Total charge density g-space grids: -0.0000021213 Overlap energy of the core charge distribution: 0.00000004306565 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.81643204996297 Hartree energy: 17.97793828310673 Exchange-correlation energy: -4.12592785667717 Total energy: -17.16444802645666 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164456561472882 Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.36993739 -17.1506340541 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.04418118 -17.2963310096 -1.46E-01 3 Pulay/Diag. 0.50E+00 1.1 0.03502972 -17.1728051826 1.24E-01 4 Pulay/Diag. 0.50E+00 1.1 0.01877120 -17.2351880274 -6.24E-02 5 Pulay/Diag. 0.50E+00 1.2 0.00165326 -17.1601551985 7.50E-02 6 Pulay/Diag. 0.50E+00 1.2 0.00030792 -17.1639770117 -3.82E-03 7 Pulay/Diag. 0.50E+00 1.2 0.00007868 -17.1643990707 -4.22E-04 8 Pulay/Diag. 0.50E+00 1.2 0.00005291 -17.1644539779 -5.49E-05 9 Pulay/Diag. 0.50E+00 1.2 0.00002962 -17.1644621617 -8.18E-06 10 Pulay/Diag. 0.50E+00 1.2 0.00001452 -17.1644465539 1.56E-05 11 Pulay/Diag. 0.50E+00 1.2 0.00001079 -17.1644240145 2.25E-05 12 Pulay/Diag. 0.50E+00 1.2 0.00000612 -17.1644495013 -2.55E-05 *** SCF run converged in 12 steps *** Electronic density on regular grids: -8.0000016377 -0.0000016377 Core density on regular grids: 7.9999993838 -0.0000006162 Total charge density on r-space grids: -0.0000022540 Total charge density g-space grids: -0.0000022540 Overlap energy of the core charge distribution: 0.00000005226004 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.83071767514271 Hartree energy: 17.96768926717895 Exchange-correlation energy: -4.12996594992828 Total energy: -17.16444950126142 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164438492138181 Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.11385974 -17.1628789790 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.01604405 -17.1140682334 4.88E-02 3 Pulay/Diag. 0.50E+00 1.1 0.01257779 -17.1635768100 -4.95E-02 4 Pulay/Diag. 0.50E+00 1.1 0.00646939 -17.1410573416 2.25E-02 5 Pulay/Diag. 0.50E+00 1.1 0.00055978 -17.1661169134 -2.51E-02 6 Pulay/Diag. 0.50E+00 1.2 0.00014433 -17.1648052030 1.31E-03 7 Pulay/Diag. 0.50E+00 1.2 0.00004555 -17.1645534645 2.52E-04 8 Pulay/Diag. 0.50E+00 1.2 0.00002064 -17.1645218909 3.16E-05 9 Pulay/Diag. 0.50E+00 1.2 0.00001152 -17.1645126986 9.19E-06 10 Pulay/Diag. 0.50E+00 1.2 0.00000553 -17.1645155009 -2.80E-06 *** SCF run converged in 10 steps *** Electronic density on regular grids: -8.0000016528 -0.0000016528 Core density on regular grids: 7.9999994648 -0.0000005352 Total charge density on r-space grids: -0.0000021880 Total charge density g-space grids: -0.0000021880 Overlap energy of the core charge distribution: 0.00000004613078 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.82283013639058 Hartree energy: 17.97329616166822 Exchange-correlation energy: -4.12775129917329 Total energy: -17.16451550089855 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164520255195153 Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.01348138 -17.1644930707 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.00343047 -17.1524130021 1.21E-02 3 Pulay/Diag. 0.50E+00 1.1 0.00162718 -17.1669321542 -1.45E-02 4 Pulay/Diag. 0.50E+00 1.1 0.00057456 -17.1629541245 3.98E-03 5 Pulay/Diag. 0.50E+00 1.1 0.00006024 -17.1646708787 -1.72E-03 6 Pulay/Diag. 0.50E+00 1.2 0.00001973 -17.1645592630 1.12E-04 7 Pulay/Diag. 0.50E+00 1.2 0.00000659 -17.1645245804 3.47E-05 *** SCF run converged in 7 steps *** Electronic density on regular grids: -8.0000016523 -0.0000016523 Core density on regular grids: 7.9999994631 -0.0000005369 Total charge density on r-space grids: -0.0000021892 Total charge density g-space grids: -0.0000021892 Overlap energy of the core charge distribution: 0.00000004627194 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.82304806049634 Hartree energy: 17.97313265238689 Exchange-correlation energy: -4.12781479359998 Total energy: -17.16452458035964 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164520335011918 Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.02624365 -17.1644328282 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.00441880 -17.1793173358 -1.49E-02 3 Pulay/Diag. 0.50E+00 1.1 0.00289721 -17.1632128539 1.61E-02 4 Pulay/Diag. 0.50E+00 1.1 0.00134755 -17.1689948440 -5.78E-03 5 Pulay/Diag. 0.50E+00 1.1 0.00011930 -17.1642171873 4.78E-03 6 Pulay/Diag. 0.50E+00 1.2 0.00003163 -17.1644636011 -2.46E-04 7 Pulay/Diag. 0.50E+00 1.2 0.00000934 -17.1645154722 -5.19E-05 *** SCF run converged in 7 steps *** Electronic density on regular grids: -8.0000016522 -0.0000016522 Core density on regular grids: 7.9999994624 -0.0000005376 Total charge density on r-space grids: -0.0000021898 Total charge density g-space grids: -0.0000021898 Overlap energy of the core charge distribution: 0.00000004631378 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.82312122727550 Hartree energy: 17.97309130540068 Exchange-correlation energy: -4.12783750524087 Total energy: -17.16451547216575 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164520329800375 ******************************************************************************* *** BRENT - NUMBER OF ENERGY EVALUATIONS : 5 *** ******************************************************************************* -------- Informations at step = 5 ------------ Optimization Method = SD Total Energy = -17.1645203350 Real energy change = -0.0000637735 Decrease in energy = YES Used time = 96.410 Convergence check : Max. step size = 0.0129135089 Conv. limit for step size = 0.0010000000 Convergence in step size = NO RMS step size = 0.0054077748 Conv. limit for RMS step = 0.0010000000 Convergence in RMS step = NO Max. gradient = 0.0061936998 Conv. limit for gradients = 0.0010000000 Conv. for gradients = NO RMS gradient = 0.0025937283 Conv. limit for RMS grad. = 0.0010000000 Conv. for gradients = NO --------------------------------------------------- Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.00850269 -17.1645110119 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.00140920 -17.1597753529 4.74E-03 3 Pulay/Diag. 0.50E+00 1.1 0.00093286 -17.1649097603 -5.13E-03 4 Pulay/Diag. 0.50E+00 1.1 0.00044059 -17.1630403614 1.87E-03 5 Pulay/Diag. 0.50E+00 1.1 0.00003873 -17.1646222353 -1.58E-03 6 Pulay/Diag. 0.50E+00 1.2 0.00001087 -17.1645408668 8.14E-05 7 Pulay/Diag. 0.50E+00 1.2 0.00000355 -17.1645225474 1.83E-05 *** SCF run converged in 7 steps *** Electronic density on regular grids: -8.0000016523 -0.0000016523 Core density on regular grids: 7.9999994631 -0.0000005369 Total charge density on r-space grids: -0.0000021892 Total charge density g-space grids: -0.0000021892 Overlap energy of the core charge distribution: 0.00000004627194 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.82304906897861 Hartree energy: 17.97313408682487 Exchange-correlation energy: -4.12781520358119 Total energy: -17.16452254742060 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164520335018693 -------------------------- OPTIMIZATION STEP: 6 -------------------------- Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.00073396 -17.1645202614 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.00012615 -17.1649655735 -4.45E-04 3 Pulay/Diag. 0.50E+00 1.1 0.00006467 -17.1644684046 4.97E-04 4 Pulay/Diag. 0.50E+00 1.1 0.00003298 -17.1646394972 -1.71E-04 5 Pulay/Diag. 0.50E+00 1.2 0.00000246 -17.1645124660 1.27E-04 *** SCF run converged in 5 steps *** Electronic density on regular grids: -8.0000016523 -0.0000016523 Core density on regular grids: 7.9999994631 -0.0000005369 Total charge density on r-space grids: -0.0000021892 Total charge density g-space grids: -0.0000021892 Overlap energy of the core charge distribution: 0.00000004627194 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.82305638473245 Hartree energy: 17.97313749892093 Exchange-correlation energy: -4.12781585004606 Total energy: -17.16451246603558 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164512466035578 Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.00881078 -17.1645448502 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.00013230 -17.1650395681 -4.95E-04 3 Pulay/Diag. 0.50E+00 1.1 0.00018442 -17.1648066269 2.33E-04 4 Pulay/Diag. 0.50E+00 1.1 0.00015167 -17.1651798746 -3.73E-04 5 Pulay/Diag. 0.50E+00 1.2 0.00001196 -17.1644956790 6.84E-04 6 Pulay/Diag. 0.50E+00 1.2 0.00000064 -17.1645497653 -5.41E-05 *** SCF run converged in 6 steps *** Electronic density on regular grids: -8.0000016469 -0.0000016469 Core density on regular grids: 7.9999994033 -0.0000005967 Total charge density on r-space grids: -0.0000022436 Total charge density g-space grids: -0.0000022436 Overlap energy of the core charge distribution: 0.00000004605680 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.82265481394913 Hartree energy: 17.97338758072965 Exchange-correlation energy: -4.12770166007671 Total energy: -17.16454976525596 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164549765255963 Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.00340909 -17.1645891705 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.00008617 -17.1643457869 2.43E-04 3 Pulay/Diag. 0.50E+00 1.1 0.00008119 -17.1645023662 -1.57E-04 4 Pulay/Diag. 0.50E+00 1.1 0.00006542 -17.1642999602 2.02E-04 5 Pulay/Diag. 0.50E+00 1.2 0.00000338 -17.1646033970 -3.03E-04 *** SCF run converged in 5 steps *** Electronic density on regular grids: -8.0000016393 -0.0000016393 Core density on regular grids: 7.9999993282 -0.0000006718 Total charge density on r-space grids: -0.0000023111 Total charge density g-space grids: -0.0000023111 Overlap energy of the core charge distribution: 0.00000004562188 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.82183721810896 Hartree energy: 17.97392077824134 Exchange-correlation energy: -4.12747089303896 Total energy: -17.16460339698163 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164603396981626 Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.06805402 -17.1641562975 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.00266348 -17.1735917616 -9.44E-03 3 Pulay/Diag. 0.50E+00 1.1 0.00190441 -17.1652533477 8.34E-03 4 Pulay/Diag. 0.50E+00 1.1 0.00131803 -17.1702548869 -5.00E-03 5 Pulay/Diag. 0.50E+00 1.2 0.00005063 -17.1642259702 6.03E-03 6 Pulay/Diag. 0.50E+00 1.2 0.00000537 -17.1644572823 -2.31E-04 *** SCF run converged in 6 steps *** Electronic density on regular grids: -8.0000016247 -0.0000016247 Core density on regular grids: 7.9999995024 -0.0000004976 Total charge density on r-space grids: -0.0000021223 Total charge density g-space grids: -0.0000021223 Overlap energy of the core charge distribution: 0.00000004104366 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.81291427822038 Hartree energy: 17.98044907083172 Exchange-correlation energy: -4.12493012652888 Total energy: -17.16445728234796 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164457282347961 ******************************************************************************* *** MNBRACK - NUMBER OF ENERGY EVALUATIONS : 4 *** ******************************************************************************* Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.23572683 -17.1606718714 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.00528072 -17.1460569769 1.46E-02 3 Pulay/Diag. 0.50E+00 1.1 0.00592455 -17.1573485726 -1.13E-02 4 Pulay/Diag. 0.50E+00 1.1 0.00476389 -17.1442096281 1.31E-02 5 Pulay/Diag. 0.50E+00 1.2 0.00026637 -17.1658766648 -2.17E-02 6 Pulay/Diag. 0.50E+00 1.2 0.00002493 -17.1646780430 1.20E-03 7 Pulay/Diag. 0.50E+00 1.2 0.00000686 -17.1646112575 6.68E-05 *** SCF run converged in 7 steps *** Electronic density on regular grids: -8.0000016393 -0.0000016393 Core density on regular grids: 7.9999993282 -0.0000006718 Total charge density on r-space grids: -0.0000023111 Total charge density g-space grids: -0.0000023111 Overlap energy of the core charge distribution: 0.00000004562188 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.82183899648170 Hartree energy: 17.97390871373897 Exchange-correlation energy: -4.12746846739720 Total energy: -17.16461125746950 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164589943296946 Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.32354018 -17.1572971743 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.00502140 -17.1869817081 -2.97E-02 3 Pulay/Diag. 0.50E+00 1.1 0.00850605 -17.1747191488 1.23E-02 4 Pulay/Diag. 0.50E+00 1.1 0.00571918 -17.1879217519 -1.32E-02 5 Pulay/Diag. 0.50E+00 1.2 0.00051990 -17.1623042300 2.56E-02 6 Pulay/Diag. 0.50E+00 1.2 0.00002421 -17.1645753245 -2.27E-03 7 Pulay/Diag. 0.50E+00 1.2 0.00000828 -17.1645926573 -1.73E-05 *** SCF run converged in 7 steps *** Electronic density on regular grids: -8.0000016244 -0.0000016244 Core density on regular grids: 7.9999992933 -0.0000007067 Total charge density on r-space grids: -0.0000023311 Total charge density g-space grids: -0.0000023311 Overlap energy of the core charge distribution: 0.00000004372544 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.81826465008933 Hartree energy: 17.97648692888402 Exchange-correlation energy: -4.12645373408194 Total energy: -17.16459265729798 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164603243182803 Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.26246941 -17.1597768449 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.00848287 -17.1369545381 2.28E-02 3 Pulay/Diag. 0.50E+00 1.1 0.00692765 -17.1595640756 -2.26E-02 4 Pulay/Diag. 0.50E+00 1.1 0.00518949 -17.1418557951 1.77E-02 5 Pulay/Diag. 0.50E+00 1.2 0.00021141 -17.1657183054 -2.39E-02 6 Pulay/Diag. 0.50E+00 1.2 0.00002417 -17.1647341660 9.84E-04 7 Pulay/Diag. 0.50E+00 1.2 0.00001017 -17.1646487824 8.54E-05 8 Pulay/Diag. 0.50E+00 1.2 0.00000687 -17.1646263578 2.24E-05 *** SCF run converged in 8 steps *** Electronic density on regular grids: -8.0000016292 -0.0000016292 Core density on regular grids: 7.9999992537 -0.0000007463 Total charge density on r-space grids: -0.0000023755 Total charge density g-space grids: -0.0000023755 Overlap energy of the core charge distribution: 0.00000004463412 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.81999541960618 Hartree energy: 17.97521086107663 Exchange-correlation energy: -4.12694213716007 Total energy: -17.16462635775798 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164619031240079 Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.14393262 -17.1631824163 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.00242056 -17.1743939278 -1.12E-02 3 Pulay/Diag. 0.50E+00 1.1 0.00356601 -17.1691280488 5.27E-03 4 Pulay/Diag. 0.50E+00 1.1 0.00263624 -17.1755872235 -6.46E-03 5 Pulay/Diag. 0.50E+00 1.1 0.00020247 -17.1636994195 1.19E-02 6 Pulay/Diag. 0.50E+00 1.2 0.00001045 -17.1645987653 -8.99E-04 7 Pulay/Diag. 0.50E+00 1.2 0.00000324 -17.1646128739 -1.41E-05 *** SCF run converged in 7 steps *** Electronic density on regular grids: -8.0000016293 -0.0000016293 Core density on regular grids: 7.9999992546 -0.0000007454 Total charge density on r-space grids: -0.0000023747 Total charge density g-space grids: -0.0000023747 Overlap energy of the core charge distribution: 0.00000004465610 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.82004323827729 Hartree energy: 17.97519120603899 Exchange-correlation energy: -4.12695681694754 Total energy: -17.16461287389000 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164619024126480 ******************************************************************************* *** BRENT - NUMBER OF ENERGY EVALUATIONS : 4 *** ******************************************************************************* -------- Informations at step = 6 ------------ Optimization Method = CG Total Energy = -17.1646190312 Real energy change = -0.0000986962 Decrease in energy = YES Used time = 74.884 Convergence check : Max. step size = 0.0581805997 Conv. limit for step size = 0.0010000000 Convergence in step size = NO RMS step size = 0.0233431802 Conv. limit for RMS step = 0.0010000000 Convergence in RMS step = NO Max. gradient = 0.0021850428 Conv. limit for gradients = 0.0010000000 Conv. for gradients = NO RMS gradient = 0.0008766814 Conv. limit for RMS grad. = 0.0010000000 Conv. in RMS gradients = YES --------------------------------------------------- Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.05129509 -17.1644362307 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.00126231 -17.1601905220 4.25E-03 3 Pulay/Diag. 0.50E+00 1.1 0.00129221 -17.1631740381 -2.98E-03 4 Pulay/Diag. 0.50E+00 1.1 0.00103380 -17.1601400071 3.03E-03 5 Pulay/Diag. 0.50E+00 1.1 0.00004971 -17.1648646832 -4.72E-03 6 Pulay/Diag. 0.50E+00 1.2 0.00000353 -17.1646343788 2.30E-04 *** SCF run converged in 6 steps *** Electronic density on regular grids: -8.0000016292 -0.0000016292 Core density on regular grids: 7.9999992537 -0.0000007463 Total charge density on r-space grids: -0.0000023755 Total charge density g-space grids: -0.0000023755 Overlap energy of the core charge distribution: 0.00000004463412 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.81998766075189 Hartree energy: 17.97521025648709 Exchange-correlation energy: -4.12694179478723 Total energy: -17.16463437882897 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164619031250190 -------------------------- OPTIMIZATION STEP: 7 -------------------------- Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.01035380 -17.1646115925 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.0 0.00023132 -17.1654560425 -8.44E-04 3 Pulay/Diag. 0.50E+00 1.1 0.00025709 -17.1649293265 5.27E-04 4 Pulay/Diag. 0.50E+00 1.1 0.00020683 -17.1655075916 -5.78E-04 5 Pulay/Diag. 0.50E+00 1.1 0.00001115 -17.1645650594 9.43E-04 6 Pulay/Diag. 0.50E+00 1.2 0.00000072 -17.1646161119 -5.11E-05 *** SCF run converged in 6 steps *** Electronic density on regular grids: -8.0000016292 -0.0000016292 Core density on regular grids: 7.9999992537 -0.0000007463 Total charge density on r-space grids: -0.0000023755 Total charge density g-space grids: -0.0000023755 Overlap energy of the core charge distribution: 0.00000004463412 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.81999992631328 Hartree energy: 17.97521753183770 Exchange-correlation energy: -4.12694306880797 Total energy: -17.16461611193772 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164616111937722 Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.03291282 -17.1638434922 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.01323909 -17.1151682992 4.87E-02 3 Pulay/Diag. 0.50E+00 1.1 0.00491466 -17.1747370274 -5.96E-02 4 Pulay/Diag. 0.50E+00 1.1 0.00088101 -17.1623904993 1.23E-02 5 Pulay/Diag. 0.50E+00 1.1 0.00013946 -17.1646417887 -2.25E-03 6 Pulay/Diag. 0.50E+00 1.2 0.00005533 -17.1641981713 4.44E-04 7 Pulay/Diag. 0.50E+00 1.2 0.00001456 -17.1640520868 1.46E-04 8 Pulay/Diag. 0.50E+00 1.2 0.00000749 -17.1640326754 1.94E-05 *** SCF run converged in 8 steps *** Electronic density on regular grids: -8.0000016396 -0.0000016396 Core density on regular grids: 7.9999993870 -0.0000006130 Total charge density on r-space grids: -0.0000022526 Total charge density g-space grids: -0.0000022526 Overlap energy of the core charge distribution: 0.00000008238221 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.87418130695149 Hartree energy: 17.93691191704169 Exchange-correlation energy: -4.14223543584034 Total energy: -17.16403267537979 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164032675379790 Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.15554495 -17.1590276552 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.06723486 -17.3955048816 -2.36E-01 3 Pulay/Diag. 0.50E+00 1.1 0.02454275 -17.0964608087 2.99E-01 4 Pulay/Diag. 0.50E+00 1.1 0.00372758 -17.1685087348 -7.20E-02 5 Pulay/Diag. 0.50E+00 1.1 0.00064673 -17.1600710662 8.44E-03 6 Pulay/Diag. 0.50E+00 1.2 0.00025750 -17.1621613007 -2.09E-03 7 Pulay/Diag. 0.50E+00 1.2 0.00006447 -17.1628674529 -7.06E-04 8 Pulay/Diag. 0.50E+00 1.2 0.00002989 -17.1629791424 -1.12E-04 9 Pulay/Diag. 0.50E+00 1.2 0.00003480 -17.1629602248 1.89E-05 10 Pulay/Diag. 0.50E+00 1.2 0.00000996 -17.1629932301 -3.30E-05 *** SCF run converged in 10 steps *** Electronic density on regular grids: -8.0000016155 -0.0000016155 Core density on regular grids: 7.9999994228 -0.0000005772 Total charge density on r-space grids: -0.0000021927 Total charge density g-space grids: -0.0000021927 Overlap energy of the core charge distribution: 0.00000001698392 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.73720386427723 Hartree energy: 18.03568838615832 Exchange-correlation energy: -4.10299495161532 Total energy: -17.16299323011069 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.162993230110693 ******************************************************************************* *** MNBRACK - NUMBER OF ENERGY EVALUATIONS : 3 *** ******************************************************************************* Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.24672054 -17.1520733535 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.11586557 -16.7476760020 4.04E-01 3 Pulay/Diag. 0.50E+00 1.1 0.04166585 -17.2625082962 -5.15E-01 4 Pulay/Diag. 0.50E+00 1.1 0.00674385 -17.1539164762 1.09E-01 5 Pulay/Diag. 0.50E+00 1.2 0.00122195 -17.1701988698 -1.63E-02 6 Pulay/Diag. 0.50E+00 1.2 0.00049940 -17.1662190090 3.98E-03 7 Pulay/Diag. 0.50E+00 1.2 0.00013057 -17.1648473193 1.37E-03 8 Pulay/Diag. 0.50E+00 1.2 0.00007046 -17.1646693316 1.78E-04 9 Pulay/Diag. 0.50E+00 1.2 0.00006269 -17.1646470441 2.23E-05 10 Pulay/Diag. 0.50E+00 1.2 0.00001888 -17.1645906569 5.64E-05 11 Pulay/Diag. 0.50E+00 1.2 0.00000912 -17.1646318420 -4.12E-05 *** SCF run converged in 11 steps *** Electronic density on regular grids: -8.0000016292 -0.0000016292 Core density on regular grids: 7.9999992537 -0.0000007463 Total charge density on r-space grids: -0.0000023755 Total charge density g-space grids: -0.0000023755 Overlap energy of the core charge distribution: 0.00000004463412 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.81998901870742 Hartree energy: 17.97521234569993 Exchange-correlation energy: -4.12694270511912 Total energy: -17.16463184199250 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164619031245365 Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.24313246 -17.1559216825 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.08484187 -17.4773142788 -3.21E-01 3 Pulay/Diag. 0.50E+00 1.1 0.03303578 -17.0932556610 3.84E-01 4 Pulay/Diag. 0.50E+00 1.1 0.00640358 -17.1780758500 -8.48E-02 5 Pulay/Diag. 0.50E+00 1.1 0.00089136 -17.1605215013 1.76E-02 6 Pulay/Diag. 0.50E+00 1.2 0.00034962 -17.1634352599 -2.91E-03 7 Pulay/Diag. 0.50E+00 1.2 0.00008861 -17.1643609013 -9.26E-04 8 Pulay/Diag. 0.50E+00 1.2 0.00004162 -17.1644900910 -1.29E-04 9 Pulay/Diag. 0.50E+00 1.2 0.00004763 -17.1644787454 1.13E-05 10 Pulay/Diag. 0.50E+00 1.2 0.00001375 -17.1645128048 -3.41E-05 11 Pulay/Diag. 0.50E+00 1.2 0.00000641 -17.1644816139 3.12E-05 *** SCF run converged in 11 steps *** Electronic density on regular grids: -8.0000016336 -0.0000016336 Core density on regular grids: 7.9999993150 -0.0000006850 Total charge density on r-space grids: -0.0000023186 Total charge density g-space grids: -0.0000023186 Overlap energy of the core charge distribution: 0.00000006054819 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.84679611377784 Hartree energy: 17.95615800897550 Exchange-correlation energy: -4.13454525127842 Total energy: -17.16448161389172 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164486006852030 Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.07062537 -17.1637532139 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.02871434 -17.0606679250 1.03E-01 3 Pulay/Diag. 0.50E+00 1.1 0.01078719 -17.1889852782 -1.28E-01 4 Pulay/Diag. 0.50E+00 1.1 0.00191652 -17.1610329933 2.80E-02 5 Pulay/Diag. 0.50E+00 1.2 0.00030590 -17.1660438404 -5.01E-03 6 Pulay/Diag. 0.50E+00 1.2 0.00012515 -17.1650188875 1.02E-03 7 Pulay/Diag. 0.50E+00 1.2 0.00003194 -17.1646754338 3.43E-04 8 Pulay/Diag. 0.50E+00 1.2 0.00001694 -17.1646301101 4.53E-05 9 Pulay/Diag. 0.50E+00 1.2 0.00001609 -17.1646258522 4.26E-06 10 Pulay/Diag. 0.50E+00 1.2 0.00000478 -17.1646121808 1.37E-05 *** SCF run converged in 10 steps *** Electronic density on regular grids: -8.0000016295 -0.0000016295 Core density on regular grids: 7.9999992563 -0.0000007437 Total charge density on r-space grids: -0.0000023732 Total charge density g-space grids: -0.0000023732 Overlap energy of the core charge distribution: 0.00000004560945 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.82189911881427 Hartree energy: 17.97385980647447 Exchange-correlation energy: -4.12748060578139 Total energy: -17.16461218079803 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164619844525824 Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.02136322 -17.1645362267 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.00960997 -17.1301572711 3.44E-02 3 Pulay/Diag. 0.50E+00 1.1 0.00342790 -17.1731320987 -4.30E-02 4 Pulay/Diag. 0.50E+00 1.1 0.00053858 -17.1638117961 9.32E-03 5 Pulay/Diag. 0.50E+00 1.2 0.00009558 -17.1650284272 -1.22E-03 6 Pulay/Diag. 0.50E+00 1.2 0.00003751 -17.1647350148 2.93E-04 7 Pulay/Diag. 0.50E+00 1.2 0.00001011 -17.1646354657 9.95E-05 8 Pulay/Diag. 0.50E+00 1.2 0.00000481 -17.1646208484 1.46E-05 *** SCF run converged in 8 steps *** Electronic density on regular grids: -8.0000016295 -0.0000016295 Core density on regular grids: 7.9999992563 -0.0000007437 Total charge density on r-space grids: -0.0000023732 Total charge density g-space grids: -0.0000023732 Overlap energy of the core charge distribution: 0.00000004561947 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.82191268380943 Hartree energy: 17.97384175350525 Exchange-correlation energy: -4.12748478540577 Total energy: -17.16462084838645 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164619844622784 Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.02645495 -17.1645008073 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.01099479 -17.2042401236 -3.97E-02 3 Pulay/Diag. 0.50E+00 1.1 0.00403937 -17.1549896936 4.93E-02 4 Pulay/Diag. 0.50E+00 1.1 0.00068028 -17.1657955664 -1.08E-02 5 Pulay/Diag. 0.50E+00 1.1 0.00011138 -17.1641311678 1.66E-03 6 Pulay/Diag. 0.50E+00 1.2 0.00004418 -17.1644861367 -3.55E-04 7 Pulay/Diag. 0.50E+00 1.2 0.00001123 -17.1646034607 -1.17E-04 8 Pulay/Diag. 0.50E+00 1.2 0.00000553 -17.1646195619 -1.61E-05 *** SCF run converged in 8 steps *** Electronic density on regular grids: -8.0000016295 -0.0000016295 Core density on regular grids: 7.9999992564 -0.0000007436 Total charge density on r-space grids: -0.0000023731 Total charge density g-space grids: -0.0000023731 Overlap energy of the core charge distribution: 0.00000004562944 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.82191790622931 Hartree energy: 17.97384812588453 Exchange-correlation energy: -4.12749509375853 Total energy: -17.16461956193008 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164619844558686 ******************************************************************************* *** BRENT - NUMBER OF ENERGY EVALUATIONS : 5 *** ******************************************************************************* -------- Informations at step = 7 ------------ Optimization Method = CG Total Energy = -17.1646198446 Real energy change = -0.0000008134 Decrease in energy = YES Used time = 98.869 Convergence check : Max. step size = 0.0014228446 Conv. limit for step size = 0.0010000000 Convergence in step size = NO RMS step size = 0.0006717304 Conv. limit for RMS step = 0.0010000000 Convergence in RMS step = YES Max. gradient = 0.0006671286 Conv. limit for gradients = 0.0010000000 Conv. in gradients = YES RMS gradient = 0.0003149540 Conv. limit for RMS grad. = 0.0010000000 Conv. in RMS gradients = YES --------------------------------------------------- Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.00878424 -17.1646065386 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.00367811 -17.1513741181 1.32E-02 3 Pulay/Diag. 0.50E+00 1.1 0.00134365 -17.1678281950 -1.65E-02 4 Pulay/Diag. 0.50E+00 1.1 0.00022566 -17.1642356067 3.59E-03 5 Pulay/Diag. 0.50E+00 1.2 0.00003763 -17.1647881120 -5.53E-04 6 Pulay/Diag. 0.50E+00 1.2 0.00001511 -17.1646668549 1.21E-04 7 Pulay/Diag. 0.50E+00 1.2 0.00000390 -17.1646261219 4.07E-05 *** SCF run converged in 7 steps *** Electronic density on regular grids: -8.0000016295 -0.0000016295 Core density on regular grids: 7.9999992563 -0.0000007437 Total charge density on r-space grids: -0.0000023732 Total charge density g-space grids: -0.0000023732 Overlap energy of the core charge distribution: 0.00000004561947 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.82190239533363 Hartree energy: 17.97384692821715 Exchange-correlation energy: -4.12748494519372 Total energy: -17.16462612193830 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164619844627957 -------------------------- OPTIMIZATION STEP: 8 -------------------------- Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.00109478 -17.1646196473 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.00044852 -17.1662206872 -1.60E-03 3 Pulay/Diag. 0.50E+00 1.1 0.00015469 -17.1642230626 2.00E-03 4 Pulay/Diag. 0.50E+00 1.1 0.00002506 -17.1646601342 -4.37E-04 5 Pulay/Diag. 0.50E+00 1.1 0.00000442 -17.1645986306 6.15E-05 *** SCF run converged in 5 steps *** Electronic density on regular grids: -8.0000016295 -0.0000016295 Core density on regular grids: 7.9999992563 -0.0000007437 Total charge density on r-space grids: -0.0000023731 Total charge density g-space grids: -0.0000023731 Overlap energy of the core charge distribution: 0.00000004561947 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.82192241648771 Hartree energy: 17.97385678692573 Exchange-correlation energy: -4.12748733376806 Total energy: -17.16459863064998 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164598630649976 Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.00082029 -17.1646200864 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.00006429 -17.1648239627 -2.04E-04 3 Pulay/Diag. 0.50E+00 1.1 0.00002977 -17.1645963676 2.28E-04 4 Pulay/Diag. 0.50E+00 1.1 0.00001465 -17.1646859271 -8.96E-05 5 Pulay/Diag. 0.50E+00 1.1 0.00000060 -17.1646174831 6.84E-05 *** SCF run converged in 5 steps *** Electronic density on regular grids: -8.0000016296 -0.0000016296 Core density on regular grids: 7.9999992566 -0.0000007434 Total charge density on r-space grids: -0.0000023730 Total charge density g-space grids: -0.0000023730 Overlap energy of the core charge distribution: 0.00000004553109 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.82172517161986 Hartree energy: 17.97398554767842 Exchange-correlation energy: -4.12743770198658 Total energy: -17.16461748307206 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164617483072057 Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.00031364 -17.1646203116 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.00002708 -17.1645351049 8.52E-05 3 Pulay/Diag. 0.50E+00 1.1 0.00001165 -17.1646328469 -9.77E-05 4 Pulay/Diag. 0.50E+00 1.1 0.00000512 -17.1645969945 3.59E-05 *** SCF run converged in 4 steps *** Electronic density on regular grids: -8.0000016299 -0.0000016299 Core density on regular grids: 7.9999992562 -0.0000007438 Total charge density on r-space grids: -0.0000023737 Total charge density g-space grids: -0.0000023737 Overlap energy of the core charge distribution: 0.00000004538719 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.82144843762548 Hartree energy: 17.97420304744701 Exchange-correlation energy: -4.12735797904042 Total energy: -17.16459699449558 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164596994495575 ******************************************************************************* *** MNBRACK - NUMBER OF ENERGY EVALUATIONS : 3 *** ******************************************************************************* Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.00283350 -17.1646191929 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.00019374 -17.1640493075 5.70E-04 3 Pulay/Diag. 0.50E+00 1.1 0.00009749 -17.1646675808 -6.18E-04 4 Pulay/Diag. 0.50E+00 1.1 0.00005653 -17.1643727971 2.95E-04 5 Pulay/Diag. 0.50E+00 1.2 0.00000194 -17.1646325824 -2.60E-04 *** SCF run converged in 5 steps *** Electronic density on regular grids: -8.0000016296 -0.0000016296 Core density on regular grids: 7.9999992566 -0.0000007434 Total charge density on r-space grids: -0.0000023730 Total charge density g-space grids: -0.0000023730 Overlap energy of the core charge distribution: 0.00000004553109 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.82171314483836 Hartree energy: 17.97398215526428 Exchange-correlation energy: -4.12743738215267 Total energy: -17.16463258243378 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164620168163019 Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.00484643 -17.1646174536 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.00033211 -17.1655979567 -9.81E-04 3 Pulay/Diag. 0.50E+00 1.1 0.00016712 -17.1645387255 1.06E-03 4 Pulay/Diag. 0.50E+00 1.1 0.00009627 -17.1650430935 -5.04E-04 5 Pulay/Diag. 0.50E+00 1.1 0.00000323 -17.1645997049 4.43E-04 *** SCF run converged in 5 steps *** Electronic density on regular grids: -8.0000016298 -0.0000016298 Core density on regular grids: 7.9999992564 -0.0000007436 Total charge density on r-space grids: -0.0000023734 Total charge density g-space grids: -0.0000023734 Overlap energy of the core charge distribution: 0.00000004545927 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.82159073913537 Hartree energy: 17.97409856348845 Exchange-correlation energy: -4.12739850706988 Total energy: -17.16459970490162 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164620302933983 Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.00340116 -17.1646189220 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.00023344 -17.1639412276 6.78E-04 3 Pulay/Diag. 0.50E+00 1.1 0.00011756 -17.1646773641 -7.36E-04 4 Pulay/Diag. 0.50E+00 1.1 0.00006771 -17.1643226782 3.55E-04 5 Pulay/Diag. 0.50E+00 1.2 0.00000220 -17.1646350243 -3.12E-04 *** SCF run converged in 5 steps *** Electronic density on regular grids: -8.0000016299 -0.0000016299 Core density on regular grids: 7.9999992563 -0.0000007437 Total charge density on r-space grids: -0.0000023736 Total charge density g-space grids: -0.0000023736 Overlap energy of the core charge distribution: 0.00000004541009 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.82146034561394 Hartree energy: 17.97416559209969 Exchange-correlation energy: -4.12737046151202 Total energy: -17.16463502430314 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164620329384700 Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.00078583 -17.1646202537 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.00005083 -17.1647546653 -1.34E-04 3 Pulay/Diag. 0.50E+00 1.1 0.00002769 -17.1646153945 1.39E-04 4 Pulay/Diag. 0.50E+00 1.1 0.00001689 -17.1646958045 -8.04E-05 5 Pulay/Diag. 0.50E+00 1.2 0.00000051 -17.1646168710 7.89E-05 *** SCF run converged in 5 steps *** Electronic density on regular grids: -8.0000016299 -0.0000016299 Core density on regular grids: 7.9999992563 -0.0000007437 Total charge density on r-space grids: -0.0000023736 Total charge density g-space grids: -0.0000023736 Overlap energy of the core charge distribution: 0.00000004541219 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.82147929293876 Hartree energy: 17.97416627924517 Exchange-correlation energy: -4.12737194271368 Total energy: -17.16461687103239 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164620329344935 ******************************************************************************* *** BRENT - NUMBER OF ENERGY EVALUATIONS : 4 *** ******************************************************************************* -------- Informations at step = 8 ------------ Optimization Method = CG Total Energy = -17.1646203294 Real energy change = -0.0000004848 Decrease in energy = YES Used time = 54.413 Convergence check : Max. step size = 0.0022367246 Conv. limit for step size = 0.0010000000 Convergence in step size = NO RMS step size = 0.0012689684 Conv. limit for RMS step = 0.0010000000 Convergence in RMS step = NO Max. gradient = 0.0001648410 Conv. limit for gradients = 0.0010000000 Conv. in gradients = YES RMS gradient = 0.0000935198 Conv. limit for RMS grad. = 0.0010000000 Conv. in RMS gradients = YES --------------------------------------------------- Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.00028237 -17.1646203198 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.00002437 -17.1646984689 -7.81E-05 3 Pulay/Diag. 0.50E+00 1.1 0.00000987 -17.1646094938 8.90E-05 *** SCF run converged in 3 steps *** Electronic density on regular grids: -8.0000016299 -0.0000016299 Core density on regular grids: 7.9999992563 -0.0000007437 Total charge density on r-space grids: -0.0000023736 Total charge density g-space grids: -0.0000023736 Overlap energy of the core charge distribution: 0.00000004541009 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.82149476897617 Hartree energy: 17.97415909577485 Exchange-correlation energy: -4.12737285809365 Total energy: -17.16460949384737 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164620329364844 -------------------------- OPTIMIZATION STEP: 9 -------------------------- Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.00024522 -17.1646203216 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.00005205 -17.1644572130 1.63E-04 3 Pulay/Diag. 0.50E+00 1.1 0.00001229 -17.1646683798 -2.11E-04 4 Pulay/Diag. 0.50E+00 1.1 0.00000157 -17.1646150567 5.33E-05 *** SCF run converged in 4 steps *** Electronic density on regular grids: -8.0000016299 -0.0000016299 Core density on regular grids: 7.9999992563 -0.0000007437 Total charge density on r-space grids: -0.0000023736 Total charge density g-space grids: -0.0000023736 Overlap energy of the core charge distribution: 0.00000004541009 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.82147825430897 Hartree energy: 17.97416779806135 Exchange-correlation energy: -4.12737060858916 Total energy: -17.16461505672358 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164615056723584 Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.00206200 -17.1646200622 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.00015080 -17.1641162836 5.04E-04 3 Pulay/Diag. 0.50E+00 1.1 0.00007781 -17.1646721278 -5.56E-04 4 Pulay/Diag. 0.50E+00 1.1 0.00003766 -17.1644651228 2.07E-04 5 Pulay/Diag. 0.50E+00 1.2 0.00000269 -17.1646304622 -1.65E-04 *** SCF run converged in 5 steps *** Electronic density on regular grids: -8.0000016291 -0.0000016291 Core density on regular grids: 7.9999992543 -0.0000007457 Total charge density on r-space grids: -0.0000023748 Total charge density g-space grids: -0.0000023748 Overlap energy of the core charge distribution: 0.00000004577491 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.82215636889406 Hartree energy: 17.97366606158680 Exchange-correlation energy: -4.12756239251214 Total energy: -17.16463046217120 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164630462171203 Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.00077281 -17.1646190278 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.00004909 -17.1647861983 -1.67E-04 3 Pulay/Diag. 0.50E+00 1.1 0.00002855 -17.1646124523 1.74E-04 4 Pulay/Diag. 0.50E+00 1.1 0.00001589 -17.1646861835 -7.37E-05 5 Pulay/Diag. 0.50E+00 1.2 0.00000092 -17.1646152060 7.10E-05 *** SCF run converged in 5 steps *** Electronic density on regular grids: -8.0000016279 -0.0000016279 Core density on regular grids: 7.9999992593 -0.0000007407 Total charge density on r-space grids: -0.0000023686 Total charge density g-space grids: -0.0000023686 Overlap energy of the core charge distribution: 0.00000004636955 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.82328010000049 Hartree energy: 17.97286638809598 Exchange-correlation energy: -4.12787119456545 Total energy: -17.16461520601426 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164615206014261 ******************************************************************************* *** MNBRACK - NUMBER OF ENERGY EVALUATIONS : 3 *** ******************************************************************************* Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.00682266 -17.1646165838 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.00065926 -17.1668187458 -2.20E-03 3 Pulay/Diag. 0.50E+00 1.1 0.00027022 -17.1642134312 2.61E-03 4 Pulay/Diag. 0.50E+00 1.1 0.00010540 -17.1650132577 -8.00E-04 5 Pulay/Diag. 0.50E+00 1.2 0.00000925 -17.1645850818 4.28E-04 *** SCF run converged in 5 steps *** Electronic density on regular grids: -8.0000016291 -0.0000016291 Core density on regular grids: 7.9999992543 -0.0000007457 Total charge density on r-space grids: -0.0000023748 Total charge density g-space grids: -0.0000023748 Overlap energy of the core charge distribution: 0.00000004577491 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.82219129564452 Hartree energy: 17.97367813019795 Exchange-correlation energy: -4.12756400753041 Total energy: -17.16458508182787 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164620403247092 Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.00906784 -17.1646137174 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.00086758 -17.1617363025 2.88E-03 3 Pulay/Diag. 0.50E+00 1.1 0.00035377 -17.1651523825 -3.42E-03 4 Pulay/Diag. 0.50E+00 1.1 0.00013963 -17.1640966425 1.06E-03 5 Pulay/Diag. 0.50E+00 1.2 0.00001216 -17.1646681349 -5.71E-04 6 Pulay/Diag. 0.50E+00 1.2 0.00000353 -17.1646347855 3.33E-05 *** SCF run converged in 6 steps *** Electronic density on regular grids: -8.0000016295 -0.0000016295 Core density on regular grids: 7.9999992556 -0.0000007444 Total charge density on r-space grids: -0.0000023738 Total charge density g-space grids: -0.0000023738 Overlap energy of the core charge distribution: 0.00000004559224 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.82180799625579 Hartree energy: 17.97391391045046 Exchange-correlation energy: -4.12746619190965 Total energy: -17.16463478552600 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164620450993205 Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.00386576 -17.1646192278 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.00036155 -17.1658187371 -1.20E-03 3 Pulay/Diag. 0.50E+00 1.1 0.00014592 -17.1644044265 1.41E-03 4 Pulay/Diag. 0.50E+00 1.1 0.00005953 -17.1648493779 -4.45E-04 5 Pulay/Diag. 0.50E+00 1.1 0.00000489 -17.1646003203 2.49E-04 *** SCF run converged in 5 steps *** Electronic density on regular grids: -8.0000016294 -0.0000016294 Core density on regular grids: 7.9999992564 -0.0000007436 Total charge density on r-space grids: -0.0000023730 Total charge density g-space grids: -0.0000023730 Overlap energy of the core charge distribution: 0.00000004563188 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.82190911331656 Hartree energy: 17.97386943358767 Exchange-correlation energy: -4.12748836691598 Total energy: -17.16460032029470 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164620454671539 Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.00002402 -17.1646204545 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.00001003 -17.1646603579 -3.99E-05 3 Pulay/Diag. 0.50E+00 1.1 0.00000973 -17.1646174466 4.29E-05 *** SCF run converged in 3 steps *** Electronic density on regular grids: -8.0000016294 -0.0000016294 Core density on regular grids: 7.9999992564 -0.0000007436 Total charge density on r-space grids: -0.0000023730 Total charge density g-space grids: -0.0000023730 Overlap energy of the core charge distribution: 0.00000004563410 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.82190558865917 Hartree energy: 17.97385702201710 Exchange-correlation energy: -4.12748955697458 Total energy: -17.16461744657905 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164620454658039 ******************************************************************************* *** BRENT - NUMBER OF ENERGY EVALUATIONS : 4 *** ******************************************************************************* -------- Informations at step = 9 ------------ Optimization Method = CG Total Energy = -17.1646204547 Real energy change = -0.0000001253 Decrease in energy = YES Used time = 48.710 Convergence check : Max. step size = 0.0014330696 Conv. limit for step size = 0.0010000000 Convergence in step size = NO RMS step size = 0.0006178345 Conv. limit for RMS step = 0.0010000000 Convergence in RMS step = YES Max. gradient = 0.0001481613 Conv. limit for gradients = 0.0010000000 Conv. in gradients = YES RMS gradient = 0.0000638763 Conv. limit for RMS grad. = 0.0010000000 Conv. in RMS gradients = YES --------------------------------------------------- Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.00054234 -17.1646204304 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.00007162 -17.1643766835 2.44E-04 3 Pulay/Diag. 0.50E+00 1.1 0.00002722 -17.1646773891 -3.01E-04 4 Pulay/Diag. 0.50E+00 1.1 0.00000748 -17.1645984250 7.90E-05 *** SCF run converged in 4 steps *** Electronic density on regular grids: -8.0000016294 -0.0000016294 Core density on regular grids: 7.9999992564 -0.0000007436 Total charge density on r-space grids: -0.0000023730 Total charge density g-space grids: -0.0000023730 Overlap energy of the core charge distribution: 0.00000004563188 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.82191737963337 Hartree energy: 17.97386217250768 Exchange-correlation energy: -4.12748747684677 Total energy: -17.16459842498868 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164620454673436 -------------------------- OPTIMIZATION STEP: 10 -------------------------- Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.00013929 -17.1646204530 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.00002874 -17.1647172673 -9.68E-05 3 Pulay/Diag. 0.50E+00 1.1 0.00000909 -17.1645938272 1.23E-04 *** SCF run converged in 3 steps *** Electronic density on regular grids: -8.0000016294 -0.0000016294 Core density on regular grids: 7.9999992564 -0.0000007436 Total charge density on r-space grids: -0.0000023730 Total charge density g-space grids: -0.0000023730 Overlap energy of the core charge distribution: 0.00000004563188 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.82192161568590 Hartree energy: 17.97386323981737 Exchange-correlation energy: -4.12748818238799 Total energy: -17.16459382716768 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164593827167678 Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.00085427 -17.1646203492 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.0 0.00009924 -17.1649914362 -3.71E-04 3 Pulay/Diag. 0.50E+00 1.1 0.00007034 -17.1645935115 3.98E-04 4 Pulay/Diag. 0.50E+00 1.1 0.00002894 -17.1647227586 -1.29E-04 5 Pulay/Diag. 0.50E+00 1.1 0.00000257 -17.1646145621 1.08E-04 *** SCF run converged in 5 steps *** Electronic density on regular grids: -8.0000016296 -0.0000016296 Core density on regular grids: 7.9999992565 -0.0000007435 Total charge density on r-space grids: -0.0000023731 Total charge density g-space grids: -0.0000023731 Overlap energy of the core charge distribution: 0.00000004537695 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.82145846438636 Hartree energy: 17.97417818119536 Exchange-correlation energy: -4.12736070714993 Total energy: -17.16461456210610 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164614562106099 Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.00032160 -17.1646195967 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.00003500 -17.1644868112 1.33E-04 3 Pulay/Diag. 0.50E+00 1.1 0.00002670 -17.1646251298 -1.38E-04 4 Pulay/Diag. 0.50E+00 1.1 0.00001148 -17.1645776147 4.75E-05 5 Pulay/Diag. 0.50E+00 1.1 0.00000100 -17.1646219066 -4.43E-05 *** SCF run converged in 5 steps *** Electronic density on regular grids: -8.0000016299 -0.0000016299 Core density on regular grids: 7.9999992606 -0.0000007394 Total charge density on r-space grids: -0.0000023693 Total charge density g-space grids: -0.0000023693 Overlap energy of the core charge distribution: 0.00000004497109 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.82073827385480 Hartree energy: 17.97468456096975 Exchange-correlation energy: -4.12715424051979 Total energy: -17.16462190663900 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164621906638999 Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.00193966 -17.1646195247 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.00025312 -17.1636900577 9.29E-04 3 Pulay/Diag. 0.50E+00 1.1 0.00016049 -17.1647252167 -1.04E-03 4 Pulay/Diag. 0.50E+00 1.1 0.00006154 -17.1644135534 3.12E-04 5 Pulay/Diag. 0.50E+00 1.1 0.00000572 -17.1646328706 -2.19E-04 *** SCF run converged in 5 steps *** Electronic density on regular grids: -8.0000016298 -0.0000016298 Core density on regular grids: 7.9999992597 -0.0000007403 Total charge density on r-space grids: -0.0000023701 Total charge density g-space grids: -0.0000023701 Overlap energy of the core charge distribution: 0.00000004508119 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.82092414926722 Hartree energy: 17.97454359247412 Exchange-correlation energy: -4.12721011154017 Total energy: -17.16463287063248 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164632870632484 ******************************************************************************* *** MNBRACK - NUMBER OF ENERGY EVALUATIONS : 4 *** ******************************************************************************* Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.00232221 -17.1646193498 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.0 0.00029662 -17.1657134389 -1.09E-03 3 Pulay/Diag. 0.50E+00 1.1 0.00019028 -17.1645035874 1.21E-03 4 Pulay/Diag. 0.50E+00 1.1 0.00007406 -17.1648722289 -3.69E-04 5 Pulay/Diag. 0.50E+00 1.1 0.00000678 -17.1646045187 2.68E-04 *** SCF run converged in 5 steps *** Electronic density on regular grids: -8.0000016298 -0.0000016298 Core density on regular grids: 7.9999992597 -0.0000007403 Total charge density on r-space grids: -0.0000023701 Total charge density g-space grids: -0.0000023701 Overlap energy of the core charge distribution: 0.00000004508119 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.82094846999847 Hartree energy: 17.97454895504206 Exchange-correlation energy: -4.12721144292346 Total energy: -17.16460451871658 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164619926985708 Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.00168382 -17.1646199364 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.00021015 -17.1638425078 7.77E-04 3 Pulay/Diag. 0.50E+00 1.1 0.00013537 -17.1646970760 -8.55E-04 4 Pulay/Diag. 0.50E+00 1.1 0.00005367 -17.1644333907 2.64E-04 5 Pulay/Diag. 0.50E+00 1.1 0.00000481 -17.1646314301 -1.98E-04 *** SCF run converged in 5 steps *** Electronic density on regular grids: -8.0000016297 -0.0000016297 Core density on regular grids: 7.9999992576 -0.0000007424 Total charge density on r-space grids: -0.0000023721 Total charge density g-space grids: -0.0000023721 Overlap energy of the core charge distribution: 0.00000004522853 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.82118475019351 Hartree energy: 17.97435941233305 Exchange-correlation energy: -4.12728509194123 Total energy: -17.16463143010098 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164620240586157 Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.00104687 -17.1646202320 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.00012654 -17.1650913432 -4.71E-04 3 Pulay/Diag. 0.50E+00 1.1 0.00008181 -17.1645788652 5.12E-04 4 Pulay/Diag. 0.50E+00 1.1 0.00003329 -17.1647399075 -1.61E-04 5 Pulay/Diag. 0.50E+00 1.2 0.00000288 -17.1646134741 1.26E-04 *** SCF run converged in 5 steps *** Electronic density on regular grids: -8.0000016296 -0.0000016296 Core density on regular grids: 7.9999992570 -0.0000007430 Total charge density on r-space grids: -0.0000023726 Total charge density g-space grids: -0.0000023726 Overlap energy of the core charge distribution: 0.00000004530260 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.82132973268008 Hartree energy: 17.97427065209444 Exchange-correlation energy: -4.12732335829929 Total energy: -17.16461347413699 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164620349906546 Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.00036815 -17.1646203780 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.00004074 -17.1644662228 1.54E-04 3 Pulay/Diag. 0.50E+00 1.1 0.00002650 -17.1646289049 -1.63E-04 4 Pulay/Diag. 0.50E+00 1.1 0.00001168 -17.1645746803 5.42E-05 5 Pulay/Diag. 0.50E+00 1.2 0.00000089 -17.1646227830 -4.81E-05 *** SCF run converged in 5 steps *** Electronic density on regular grids: -8.0000016296 -0.0000016296 Core density on regular grids: 7.9999992567 -0.0000007433 Total charge density on r-space grids: -0.0000023729 Total charge density g-space grids: -0.0000023729 Overlap energy of the core charge distribution: 0.00000004533974 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.82138655344844 Hartree energy: 17.97422279032439 Exchange-correlation energy: -4.12734162622891 Total energy: -17.16462278303118 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164620392715179 Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.00007951 -17.1646204104 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.00000609 -17.1646454420 -2.50E-05 *** SCF run converged in 2 steps *** Electronic density on regular grids: -8.0000016296 -0.0000016296 Core density on regular grids: 7.9999992566 -0.0000007434 Total charge density on r-space grids: -0.0000023730 Total charge density g-space grids: -0.0000023730 Overlap energy of the core charge distribution: 0.00000004535834 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.82139717628210 Hartree energy: 17.97419970553043 Exchange-correlation energy: -4.12735182326466 Total energy: -17.16464544200864 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164620411153422 Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.00003528 -17.1646204196 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.00001346 -17.1646502673 -2.98E-05 3 Pulay/Diag. 0.50E+00 1.1 0.00000442 -17.1646113309 3.89E-05 *** SCF run converged in 3 steps *** Electronic density on regular grids: -8.0000016296 -0.0000016296 Core density on regular grids: 7.9999992565 -0.0000007435 Total charge density on r-space grids: -0.0000023731 Total charge density g-space grids: -0.0000023731 Overlap energy of the core charge distribution: 0.00000004536764 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.82144768272752 Hartree energy: 17.97418750925669 Exchange-correlation energy: -4.12735602234184 Total energy: -17.16461133090482 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164620419635462 Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.00003086 -17.1646204236 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.00001371 -17.1645845494 3.59E-05 3 Pulay/Diag. 0.50E+00 1.1 0.00000311 -17.1646323844 -4.78E-05 *** SCF run converged in 3 steps *** Electronic density on regular grids: -8.0000016296 -0.0000016296 Core density on regular grids: 7.9999992565 -0.0000007435 Total charge density on r-space grids: -0.0000023731 Total charge density g-space grids: -0.0000023731 Overlap energy of the core charge distribution: 0.00000004537229 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.82143257374900 Hartree energy: 17.97418348409939 Exchange-correlation energy: -4.12735794172987 Total energy: -17.16463238442401 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164620423690280 ******************************************************************************* *** BRENT - NUMBER OF ENERGY EVALUATIONS : 7 *** ******************************************************************************* -------- Informations at step = 10 ------------ Optimization Method = CG Total Energy = -17.1646204237 Real energy change = 0.0000000310 Decrease in energy = NO Used time = 68.108 Convergence check : Max. step size = 0.0010438900 Conv. limit for step size = 0.0010000000 Convergence in step size = NO RMS step size = 0.0005033489 Conv. limit for RMS step = 0.0010000000 Convergence in RMS step = YES Max. gradient = 0.0001794506 Conv. limit for gradients = 0.0010000000 Conv. in gradients = YES RMS gradient = 0.0000865285 Conv. limit for RMS grad. = 0.0010000000 Conv. in RMS gradients = YES --------------------------------------------------- Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.00001691 -17.1646204237 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.00001101 -17.1646549086 -3.45E-05 3 Pulay/Diag. 0.50E+00 1.1 0.00000233 -17.1646092303 4.57E-05 *** SCF run converged in 3 steps *** Electronic density on regular grids: -8.0000016296 -0.0000016296 Core density on regular grids: 7.9999992565 -0.0000007435 Total charge density on r-space grids: -0.0000023731 Total charge density g-space grids: -0.0000023731 Overlap energy of the core charge distribution: 0.00000004537229 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.82145624222421 Hartree energy: 17.97418318301711 Exchange-correlation energy: -4.12735815502053 Total energy: -17.16460923032174 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164620423690334 -------------------------- OPTIMIZATION STEP: 11 -------------------------- Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.00001091 -17.1646204237 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.00000606 -17.1646006807 1.97E-05 *** SCF run converged in 2 steps *** Electronic density on regular grids: -8.0000016296 -0.0000016296 Core density on regular grids: 7.9999992565 -0.0000007435 Total charge density on r-space grids: -0.0000023731 Total charge density g-space grids: -0.0000023731 Overlap energy of the core charge distribution: 0.00000004537229 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.82146815898643 Hartree energy: 17.97418005835638 Exchange-correlation energy: -4.12735839746318 Total energy: -17.16460068066291 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164600680662907 Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.00092481 -17.1646202317 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.00022723 -17.1637927816 8.27E-04 3 Pulay/Diag. 0.50E+00 1.1 0.00011314 -17.1647869515 -9.94E-04 4 Pulay/Diag. 0.50E+00 1.1 0.00003418 -17.1645380556 2.49E-04 5 Pulay/Diag. 0.50E+00 1.2 0.00000384 -17.1646319301 -9.39E-05 *** SCF run converged in 5 steps *** Electronic density on regular grids: -8.0000016290 -0.0000016290 Core density on regular grids: 7.9999992545 -0.0000007455 Total charge density on r-space grids: -0.0000023745 Total charge density g-space grids: -0.0000023745 Overlap energy of the core charge distribution: 0.00000004586647 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.82231342072494 Hartree energy: 17.97355238113673 Exchange-correlation energy: -4.12760723193098 Total energy: -17.16463193011768 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164631930117682 Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.00035204 -17.1646185155 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.00008938 -17.1649425405 -3.24E-04 3 Pulay/Diag. 0.50E+00 1.1 0.00004240 -17.1645512785 3.91E-04 4 Pulay/Diag. 0.50E+00 1.1 0.00001252 -17.1646482229 -9.69E-05 5 Pulay/Diag. 0.50E+00 1.1 0.00000143 -17.1646140823 3.41E-05 *** SCF run converged in 5 steps *** Electronic density on regular grids: -8.0000016281 -0.0000016281 Core density on regular grids: 7.9999992499 -0.0000007501 Total charge density on r-space grids: -0.0000023782 Total charge density g-space grids: -0.0000023782 Overlap energy of the core charge distribution: 0.00000004669072 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.82374894837757 Hartree energy: 17.97253991673188 Exchange-correlation energy: -4.12801244815162 Total energy: -17.16461408226629 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164614082266290 ******************************************************************************* *** MNBRACK - NUMBER OF ENERGY EVALUATIONS : 3 *** ******************************************************************************* Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.00319148 -17.1646192320 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.00075778 -17.1673913679 -2.77E-03 3 Pulay/Diag. 0.50E+00 1.1 0.00038755 -17.1640794008 3.31E-03 4 Pulay/Diag. 0.50E+00 1.1 0.00011983 -17.1649168128 -8.37E-04 5 Pulay/Diag. 0.50E+00 1.1 0.00001319 -17.1645808666 3.36E-04 6 Pulay/Diag. 0.50E+00 1.2 0.00000468 -17.1646095132 -2.86E-05 *** SCF run converged in 6 steps *** Electronic density on regular grids: -8.0000016290 -0.0000016290 Core density on regular grids: 7.9999992545 -0.0000007455 Total charge density on r-space grids: -0.0000023745 Total charge density g-space grids: -0.0000023745 Overlap energy of the core charge distribution: 0.00000004586647 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.82233091077761 Hartree energy: 17.97355870319785 Exchange-correlation energy: -4.12760862714417 Total energy: -17.16460951321707 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164620324279021 Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.00422657 -17.1646185435 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.00100909 -17.1609327136 3.69E-03 3 Pulay/Diag. 0.50E+00 1.1 0.00051174 -17.1653437131 -4.41E-03 4 Pulay/Diag. 0.50E+00 1.1 0.00015800 -17.1642301197 1.11E-03 5 Pulay/Diag. 0.50E+00 1.2 0.00001749 -17.1646734577 -4.43E-04 6 Pulay/Diag. 0.50E+00 1.2 0.00000629 -17.1646353170 3.81E-05 *** SCF run converged in 6 steps *** Electronic density on regular grids: -8.0000016293 -0.0000016293 Core density on regular grids: 7.9999992550 -0.0000007450 Total charge density on r-space grids: -0.0000023743 Total charge density g-space grids: -0.0000023743 Overlap energy of the core charge distribution: 0.00000004561795 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.82187240192352 Hartree energy: 17.97386472673639 Exchange-correlation energy: -4.12748194532071 Total energy: -17.16463531695768 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164620470222580 Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.00202918 -17.1646200334 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.00048217 -17.1663811986 -1.76E-03 3 Pulay/Diag. 0.50E+00 1.1 0.00024088 -17.1642738588 2.11E-03 4 Pulay/Diag. 0.50E+00 1.1 0.00007470 -17.1648064916 -5.33E-04 5 Pulay/Diag. 0.50E+00 1.1 0.00000826 -17.1645951758 2.11E-04 *** SCF run converged in 5 steps *** Electronic density on regular grids: -8.0000016294 -0.0000016294 Core density on regular grids: 7.9999992554 -0.0000007446 Total charge density on r-space grids: -0.0000023740 Total charge density g-space grids: -0.0000023740 Overlap energy of the core charge distribution: 0.00000004555422 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.82179036365332 Hartree energy: 17.97395678767471 Exchange-correlation energy: -4.12745182674534 Total energy: -17.16459517577794 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164620476604895 Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.00043987 -17.1646204556 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.1 0.00009659 -17.1642662312 3.54E-04 3 Pulay/Diag. 0.50E+00 1.1 0.00004548 -17.1646866931 -4.20E-04 4 Pulay/Diag. 0.50E+00 1.1 0.00001528 -17.1645778057 1.09E-04 5 Pulay/Diag. 0.50E+00 1.1 0.00000160 -17.1646260706 -4.83E-05 *** SCF run converged in 5 steps *** Electronic density on regular grids: -8.0000016294 -0.0000016294 Core density on regular grids: 7.9999992554 -0.0000007446 Total charge density on r-space grids: -0.0000023740 Total charge density g-space grids: -0.0000023740 Overlap energy of the core charge distribution: 0.00000004555605 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.82176933408484 Hartree energy: 17.97394644419256 Exchange-correlation energy: -4.12745134851572 Total energy: -17.16462607059711 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164620476603666 ******************************************************************************* *** BRENT - NUMBER OF ENERGY EVALUATIONS : 4 *** ******************************************************************************* -------- Informations at step = 11 ------------ Optimization Method = SD Total Energy = -17.1646204766 Real energy change = -0.0000000529 Decrease in energy = YES Used time = 49.893 Convergence check : Max. step size = 0.0003393150 Conv. limit for step size = 0.0010000000 Convergence in step size = YES RMS step size = 0.0001493298 Conv. limit for RMS step = 0.0010000000 Convergence in RMS step = YES Max. gradient = 0.0001787448 Conv. limit for gradients = 0.0010000000 Conv. in gradients = YES RMS gradient = 0.0000786642 Conv. limit for RMS grad. = 0.0010000000 Conv. in RMS gradients = YES --------------------------------------------------- ******************************************************************************* *** GEOMETRY OPTIMIZATION COMPLETED *** ******************************************************************************* Reevaluating energy at the minimum Number of electrons: 8 Number of occupied orbitals: 4 Number of molecular orbitals: 4 Number of orbital functions: 23 Number of independent orbital functions: 23 Parameters for the always stable predictor-corrector (ASPC) method: ASPC order: 3 B(1) = 3.000000 B(2) = -3.428571 B(3) = 1.928571 B(4) = -0.571429 B(5) = 0.071429 Extrapolation method: ASPC SCF WAVEFUNCTION OPTIMIZATION Step Update method Time Convergence Total energy Change ------------------------------------------------------------------------------ 1 Pulay/Diag. 0.50E+00 0.5 0.00005615 -17.1646204762 -1.72E+01 2 Pulay/Diag. 0.50E+00 1.0 0.00000563 -17.1646347711 -1.43E-05 *** SCF run converged in 2 steps *** Electronic density on regular grids: -8.0000016293 -0.0000016293 Core density on regular grids: 7.9999992554 -0.0000007446 Total charge density on r-space grids: -0.0000023739 Total charge density g-space grids: -0.0000023739 Overlap energy of the core charge distribution: 0.00000004555422 Self energy of the core charge distribution: -43.83289054591484 Core Hamiltonian energy: 12.82175605770555 Hartree energy: 17.97395116120845 Exchange-correlation energy: -4.12745148966141 Total energy: -17.16463477110803 ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -17.164634771108034 ------------------------------------------------------------------------------- - - - DBCSR STATISTICS - - - ------------------------------------------------------------------------------- COUNTER CPU GPU GPU% number of processed stacks 4732 0 0.0 mutmuls inhomo. stacks 0 0 0.0 matmuls 5 x 5 x 4 2379 0 0.0 matmuls 5 x 4 x 5 1560 0 0.0 matmuls 13 x 5 x 4 793 0 0.0 matmuls 5 x 13 x 4 793 0 0.0 matmuls 13 x 4 x 5 780 0 0.0 matmuls 5 x 4 x 13 780 0 0.0 matmuls 13 x 13 x 4 793 0 0.0 matmuls 13 x 4 x 13 390 0 0.0 matmuls total 8268 0 0.0 flops total 4023136 0 0.0 ------------------------------------------------------------------------------- ------------------------------------------------------------------------------- - - - R E F E R E N C E S - - - ------------------------------------------------------------------------------- CP2K version 2.5 (Development Version), the CP2K developers group (2013). CP2K is freely available from http://www.cp2k.org/ . Krack, M. THEORETICAL CHEMISTRY ACCOUNTS, 114 (1-3), 145-152 (2005). Pseudopotentials for H to Kr optimized for gradient-corrected exchange-correlation functionals. http://dx.doi.org/10.1007/s00214-005-0655-y VandeVondele, J; Krack, M; Mohamed, F; Parrinello, M; Chassaing, T; Hutter, J. COMPUTER PHYSICS COMMUNICATIONS, 167 (2), 103-128 (2005). QUICKSTEP: Fast and accurate density functional calculations using a mixed Gaussian and plane waves approach. http://dx.doi.org/10.1016/j.cpc.2004.12.014 Frigo, M; Johnson, SG. PROCEEDINGS OF THE IEEE, 93 (2), 216-231 (2005). The design and implementation of FFTW3. http://dx.doi.org/10.1109/JPROC.2004.840301 Kolafa, J. JOURNAL OF COMPUTATIONAL CHEMISTRY, 25 (3), 335-342 (2004). Time-reversible always stable predictor-corrector method for molecular dynamics of polarizable molecules. http://dx.doi.org/10.1002/jcc.10385 Hartwigsen, C; Goedecker, S; Hutter, J. PHYSICAL REVIEW B, 58 (7), 3641-3662 (1998). Relativistic separable dual-space Gaussian pseudopotentials from H to Rn. http://dx.doi.org/10.1103/PhysRevB.58.3641 Lippert, G; Hutter, J; Parrinello, M. MOLECULAR PHYSICS, 92 (3), 477-487 (1997). A hybrid Gaussian and plane wave density functional scheme. http://dx.doi.org/10.1080/002689797170220 Goedecker, S; Teter, M; Hutter, J. PHYSICAL REVIEW B, 54 (3), 1703-1710 (1996). Separable dual-space Gaussian pseudopotentials. http://dx.doi.org/10.1103/PhysRevB.54.1703 ------------------------------------------------------------------------------- - - - T I M I N G - - - ------------------------------------------------------------------------------- SUBROUTINE CALLS ASD SELF TIME TOTAL TIME AVERAGE MAXIMUM AVERAGE MAXIMUM CP2K 1 1.0 0.024 0.024 819.011 819.011 cp_cg_main 1 4.0 0.002 0.002 818.002 818.002 geoopt_cg 1 3.0 0.000 0.000 818.002 818.002 cp_geo_opt 1 2.0 0.000 0.000 818.002 818.002 cp_eval_at 101 8.5 0.013 0.013 817.975 817.975 qs_energies_scf 101 10.1 0.009 0.009 783.829 783.829 linmin_gold 11 6.0 0.000 0.000 728.224 728.224 cg_linmin 11 5.0 0.000 0.000 728.224 728.224 scf_env_do_scf 101 11.1 0.008 0.008 682.818 682.818 scf_env_do_scf_inner_loop 630 12.2 0.090 0.090 682.808 682.808 qs_forces 63 9.3 0.002 0.002 553.288 553.288 fft_wrap_pw1pw2 7750 16.6 0.185 0.185 482.816 482.816 cg_dbrent 11 7.0 0.000 0.000 466.442 466.442 cg_deval1d 52 8.0 0.001 0.001 466.441 466.441 fft_wrap_pw1pw2_100 3478 17.0 20.045 20.045 438.244 438.244 qs_ks_build_kohn_sham_matrix 693 14.0 0.109 0.109 375.842 375.842 qs_ks_update_qs_env 630 13.2 0.009 0.009 341.883 341.883 fft3d_s 7751 18.6 303.609 303.609 303.670 303.670 cg_eval1d 37 8.0 0.000 0.000 261.782 261.782 cg_mnbrak 11 7.0 0.000 0.000 261.782 261.782 qs_rho_update_rho 731 13.4 0.008 0.008 235.683 235.683 calculate_rho_elec 731 14.4 11.532 11.532 235.675 235.675 density_rs2pw 731 15.4 0.030 0.030 221.952 221.952 gspace_mixing 529 13.3 5.758 5.758 136.355 136.355 sum_up_and_integrate 693 15.0 6.735 6.735 124.648 124.648 integrate_v_rspace 693 16.0 3.302 3.302 117.909 117.909 potential_pw2rs 693 17.0 0.476 0.476 114.482 114.482 qs_vxc_create 693 15.0 0.015 0.015 114.257 114.257 xc_vxc_pw_create 693 16.0 7.562 7.562 114.242 114.242 xc_rho_set_and_dset_create 693 17.0 0.034 0.034 106.677 106.677 xc_functional_eval 693 18.0 103.681 103.681 103.681 103.681 init_scf_run 101 11.1 0.012 0.012 86.429 86.429 scf_env_initial_rho_setup 101 12.1 0.004 0.004 86.386 86.386 pw_scatter_s 4032 18.7 85.025 85.025 85.025 85.025 pw_gather_s 3718 18.6 72.024 72.024 72.024 72.024 pulay_mixing 529 14.3 62.725 62.725 62.770 62.770 mixing_init 101 13.1 53.207 53.207 53.207 53.207 fft_wrap_pw1pw2_40 1424 18.1 1.544 1.544 38.573 38.573 pw_poisson_solve 693 15.0 23.163 23.163 34.507 34.507 qs_ks_update_qs_env_forces 63 10.3 0.000 0.000 33.992 33.992 wfi_extrapolate 101 13.1 0.011 0.011 32.792 32.792 pw_copy 5544 17.1 24.670 24.670 24.670 24.670 ------------------------------------------------------------------------------- **** **** ****** ** PROGRAM ENDED AT 2013-12-03 19:30:01.571 ***** ** *** *** ** PROGRAM RAN ON spring.kcl.ac.uk ** **** ****** PROGRAM RAN BY tong ***** ** ** ** ** PROGRAM PROCESS ID 24417 **** ** ******* ** PROGRAM STOPPED IN /home/tong/tutorials/geometry_optimis ation/example_1