unittest.out 27.8 KB
Newer Older
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 217 218 219 220 221 222 223 224 225 226 227 228 229 230 231 232 233 234 235 236 237 238 239 240 241 242 243 244 245 246 247 248 249 250 251 252 253 254 255 256 257 258 259 260 261 262 263 264 265 266 267 268 269 270 271 272 273 274 275 276 277 278 279 280 281 282 283 284 285 286 287 288 289 290 291 292 293 294 295 296 297 298 299 300 301 302 303 304 305 306 307 308 309 310 311 312 313 314 315 316 317 318 319 320 321 322 323 324 325 326 327 328 329 330 331 332 333 334 335 336 337 338 339 340 341 342 343 344 345 346 347 348 349 350 351 352 353 354 355 356 357 358 359 360 361 362 363 364 365 366 367 368 369 370 371 372 373 374 375 376 377 378 379 380 381 382 383 384 385 386 387 388 389 390 391 392 393 394 395 396 397 398 399 400 401 402 403 404 405 406 407 408 409 410 411 412 413 414 415 416 417 418 419 420 421 422 423 424 425 426 427 428 429 430 431 432 433 434 435 436 437 438 439 440 441 442 443 444 445 446 447 448 449 450 451 452 453 454 455 456 457 458 459 460 461
 DBCSR| Multiplication driver                                                SMM
 DBCSR| Multrec recursion limit                                              512
 DBCSR| Multiplication stack size                                           1000
 DBCSR| Multiplication size stacks                                             3
 DBCSR| Use subcommunicators                                                   T
 DBCSR| Use MPI combined types                                                 F
 DBCSR| Use MPI memory allocation                                              T
 DBCSR| Use Communication thread                                               T
 DBCSR| Communication thread load                                             87


  **** **** ******  **  PROGRAM STARTED AT               2016-04-20 16:17:33.006
 ***** ** ***  *** **   PROGRAM STARTED ON                   lauri-Lenovo-Z50-70
 **    ****   ******    PROGRAM STARTED BY                                 lauri
 ***** **    ** ** **   PROGRAM PROCESS ID                                 22465
  **** **  *******  **  PROGRAM STARTED IN /home/lauri/Dropbox/nomad-dev/nomad-l
                                           ab-base/parsers/cp2k/test/unittests/c
                                           p2k_2.6.2/sic/mauri_us

 CP2K| version string:                                        CP2K version 2.6.2
 CP2K| source code revision number:                                    svn:15893
 CP2K| is freely available from                             http://www.cp2k.org/
 CP2K| Program compiled at                           ke 4.11.2015 08.48.42 +0200
 CP2K| Program compiled on                                   lauri-Lenovo-Z50-70
 CP2K| Program compiled for                          Linux-x86-64-gfortran_basic
 CP2K| Input file name                                              mauri_us.inp

 GLOBAL| Force Environment number                                              1
 GLOBAL| Basis set file name                                     ../../BASIS_SET
 GLOBAL| Geminal file name                                         BASIS_GEMINAL
 GLOBAL| Potential file name                                ../../GTH_POTENTIALS
 GLOBAL| MM Potential file name                                     MM_POTENTIAL
 GLOBAL| Coordinate file name                                      __STD_INPUT__
 GLOBAL| Method name                                                        CP2K
 GLOBAL| Project name                                                    si_bulk
 GLOBAL| Preferred FFT library                                             FFTW3
 GLOBAL| Preferred diagonalization lib.                                       SL
 GLOBAL| Run type                                                   ENERGY_FORCE
 GLOBAL| All-to-all communication in single precision                          F
 GLOBAL| FFTs using library dependent lengths                                  F
 GLOBAL| Global print level                                               MEDIUM
 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              8070380       8070380       8070380       8070380
 MEMORY| MemFree               1311556       1311556       1311556       1311556
 MEMORY| Buffers               1081548       1081548       1081548       1081548
 MEMORY| Cached                2514572       2514572       2514572       2514572
 MEMORY| Slab                   545032        545032        545032        545032
 MEMORY| SReclaimable           499260        499260        499260        499260
 MEMORY| MemLikelyFree         5406936       5406936       5406936       5406936


 *** Fundamental physical constants (SI units) ***

 *** Literature: B. J. Mohr and B. N. Taylor,
 ***             CODATA recommended values of the fundamental physical
 ***             constants: 2006, Web Version 5.1
 ***             http://physics.nist.gov/constants

 Speed of light in vacuum [m/s]                             2.99792458000000E+08
 Magnetic constant or permeability of vacuum [N/A**2]       1.25663706143592E-06
 Electric constant or permittivity of vacuum [F/m]          8.85418781762039E-12
 Planck constant (h) [J*s]                                  6.62606896000000E-34
 Planck constant (h-bar) [J*s]                              1.05457162825177E-34
 Elementary charge [C]                                      1.60217648700000E-19
 Electron mass [kg]                                         9.10938215000000E-31
 Electron g factor [ ]                                     -2.00231930436220E+00
 Proton mass [kg]                                           1.67262163700000E-27
 Fine-structure constant                                    7.29735253760000E-03
 Rydberg constant [1/m]                                     1.09737315685270E+07
 Avogadro constant [1/mol]                                  6.02214179000000E+23
 Boltzmann constant [J/K]                                   1.38065040000000E-23
 Atomic mass unit [kg]                                      1.66053878200000E-27
 Bohr radius [m]                                            5.29177208590000E-11

 *** Conversion factors ***

 [u] -> [a.u.]                                              1.82288848426455E+03
 [Angstrom] -> [Bohr] = [a.u.]                              1.88972613288564E+00
 [a.u.] = [Bohr] -> [Angstrom]                              5.29177208590000E-01
 [a.u.] -> [s]                                              2.41888432650478E-17
 [a.u.] -> [fs]                                             2.41888432650478E-02
 [a.u.] -> [J]                                              4.35974393937059E-18
 [a.u.] -> [N]                                              8.23872205491840E-08
 [a.u.] -> [K]                                              3.15774647902944E+05
 [a.u.] -> [kJ/mol]                                         2.62549961709828E+03
 [a.u.] -> [kcal/mol]                                       6.27509468713739E+02
 [a.u.] -> [Pa]                                             2.94210107994716E+13
 [a.u.] -> [bar]                                            2.94210107994716E+08
 [a.u.] -> [atm]                                            2.90362800883016E+08
 [a.u.] -> [eV]                                             2.72113838565563E+01
 [a.u.] -> [Hz]                                             6.57968392072181E+15
 [a.u.] -> [1/cm] (wave numbers)                            2.19474631370540E+05
 [a.u./Bohr**2] -> [1/cm]                                   5.14048714338585E+03
 

 CELL_TOP| Volume [angstrom^3]:                                          160.165
 CELL_TOP| Vector a [angstrom     5.431     0.000     0.000    |a| =       5.431
 CELL_TOP| Vector b [angstrom     0.000     5.431     0.000    |b| =       5.431
 CELL_TOP| Vector c [angstrom     0.000     0.000     5.431    |c| =       5.431
 CELL_TOP| Angle (b,c), alpha [degree]:                                   90.000
 CELL_TOP| Angle (a,c), beta  [degree]:                                   90.000
 CELL_TOP| Angle (a,b), gamma [degree]:                                   90.000
 CELL_TOP| Numerically orthorhombic:                                         YES

 GENERATE|  Preliminary Number of Bonds generated:                             0
 GENERATE|  Achieved consistency in connectivity generation.

 CELL| Volume [angstrom^3]:                                              160.165
 CELL| Vector a [angstrom]:       5.431     0.000     0.000    |a| =       5.431
 CELL| Vector b [angstrom]:       0.000     5.431     0.000    |b| =       5.431
 CELL| Vector c [angstrom]:       0.000     0.000     5.431    |c| =       5.431
 CELL| Angle (b,c), alpha [degree]:                                       90.000
 CELL| Angle (a,c), beta  [degree]:                                       90.000
 CELL| Angle (a,b), gamma [degree]:                                       90.000
 CELL| Numerically orthorhombic:                                             YES

 CELL_REF| Volume [angstrom^3]:                                          160.165
 CELL_REF| Vector a [angstrom     5.431     0.000     0.000    |a| =       5.431
 CELL_REF| Vector b [angstrom     0.000     5.431     0.000    |b| =       5.431
 CELL_REF| Vector c [angstrom     0.000     0.000     5.431    |c| =       5.431
 CELL_REF| Angle (b,c), alpha [degree]:                                   90.000
 CELL_REF| Angle (a,c), beta  [degree]:                                   90.000
 CELL_REF| Angle (a,b), gamma [degree]:                                   90.000
 CELL_REF| Numerically orthorhombic:                                         YES

 *******************************************************************************
 *******************************************************************************
 **                                                                           **
 **     #####                         ##              ##                      **
 **    ##   ##            ##          ##              ##                      **
 **   ##     ##                       ##            ######                    **
 **   ##     ##  ##   ##  ##   #####  ##  ##   ####   ##    #####    #####    **
 **   ##     ##  ##   ##  ##  ##      ## ##   ##      ##   ##   ##  ##   ##   **
 **   ##  ## ##  ##   ##  ##  ##      ####     ###    ##   ######   ######    **
 **    ##  ###   ##   ##  ##  ##      ## ##      ##   ##   ##       ##        **
 **     #######   #####   ##   #####  ##  ##  ####    ##    #####   ##        **
 **           ##                                                    ##        **
 **                                                                           **
 **                                                ... make the atoms dance   **
 **                                                                           **
 **            Copyright (C) by CP2K Developers Group (2000 - 2014)           **
 **                                                                           **
 *******************************************************************************

 DFT| Spin restricted Kohn-Sham (RKS) calculation                            RKS
 DFT| Multiplicity                                                             1
 DFT| Number of spin states                                                    1
 DFT| Charge                                                                   0
 DFT| Self-interaction correction (SIC)                             US/MAURI SIC
 DFT| SIC scaling parameter a                                       1.000000E+00
 DFT| SIC scaling parameter b                                       1.000000E+00
 DFT| Cutoffs: density                                              1.000000E-10
 DFT|          gradient                                             1.000000E-10
 DFT|          tau                                                  1.000000E-10
 DFT|          cutoff_smoothing_range                               0.000000E+00
 DFT| XC density smoothing                                                  NONE
 DFT| XC derivatives                                                          PW
 FUNCTIONAL| ROUTINE=NEW
 FUNCTIONAL| PADE:
 FUNCTIONAL| S. Goedecker, M. Teter and J. Hutter, Phys. Rev. B 54, 1703 (1996)

 QS| Method:                                                                 GPW
 QS| Density plane wave grid type                        NON-SPHERICAL FULLSPACE
 QS| Number of grid levels:                                                    2
 QS| Density cutoff [a.u.]:                                                 25.0
 QS| Multi grid cutoff [a.u.]: 1) grid level                                25.0
 QS|                           2) grid level                                 8.3
 QS| Grid level progression factor:                                          3.0
 QS| Relative density cutoff [a.u.]:                                        12.5
 QS| Consistent realspace mapping and integration 
 QS| Interaction thresholds: eps_pgf_orb:                                3.2E-02
 QS|                         eps_filter_matrix:                          0.0E+00
 QS|                         eps_core_charge:                            1.0E-05
 QS|                         eps_rho_gspace:                             1.0E-03
 QS|                         eps_rho_rspace:                             1.0E-03
 QS|                         eps_gvg_rspace:                             3.2E-02
 QS|                         eps_ppl:                                    1.0E-02
 QS|                         eps_ppnl:                                   3.2E-04


 ATOMIC KIND INFORMATION

  1. Atomic kind: Si                                    Number of atoms:       1

     Orbital Basis Set                                             DZVP-GTH-PADE

       Number of orbital shell sets:                                           2
       Number of orbital shells:                                               5
       Number of primitive Cartesian functions:                                5
       Number of Cartesian basis functions:                                   14
       Number of spherical basis functions:                                   13
       Norm type:                                                              2

       Normalised Cartesian orbitals:

                        Set   Shell   Orbital            Exponent    Coefficient

                          1       1    3s                1.203242       0.269412
                                                         0.468841      -0.102290
                                                         0.167986      -0.147195
                                                         0.057562      -0.015996

                          1       2    4s                1.203242       0.000000
                                                         0.468841       0.000000
                                                         0.167986       0.000000
                                                         0.057562       0.083755

                          1       3    4px               1.203242       0.085242
                                                         0.468841      -0.143473
                                                         0.167986      -0.083408
                                                         0.057562      -0.014565
                          1       3    4py               1.203242       0.085242
                                                         0.468841      -0.143473
                                                         0.167986      -0.083408
                                                         0.057562      -0.014565
                          1       3    4pz               1.203242       0.085242
                                                         0.468841      -0.143473
                                                         0.167986      -0.083408
                                                         0.057562      -0.014565

                          1       4    5px               1.203242       0.000000
                                                         0.468841       0.000000
                                                         0.167986       0.000000
                                                         0.057562       0.040189
                          1       4    5py               1.203242       0.000000
                                                         0.468841       0.000000
                                                         0.167986       0.000000
                                                         0.057562       0.040189
                          1       4    5pz               1.203242       0.000000
                                                         0.468841       0.000000
                                                         0.167986       0.000000
                                                         0.057562       0.040189

                          2       1    3dx2              0.450000       0.406941
                          2       1    3dxy              0.450000       0.704842
                          2       1    3dxz              0.450000       0.704842
                          2       1    3dy2              0.450000       0.406941
                          2       1    3dyz              0.450000       0.704842
                          2       1    3dz2              0.450000       0.406941

     Potential information for                                       GTH-PADE-q4

       Description:                       Goedecker-Teter-Hutter pseudopotential
                                           Goedecker et al., PRB 54, 1703 (1996)
                                          Hartwigsen et al., PRB 58, 3641 (1998)
                                                      Krack, TCA 114, 145 (2005)

       Gaussian exponent of the core charge distribution:               2.582645
       Electronic configuration (s p d ...):                               2   2

       Parameters of the local part of the GTH pseudopotential:

                          rloc        C1          C2          C3          C4
                        0.440000   -7.336103

       Parameters of the non-local part of the GTH pseudopotential:

                   l      r(l)      h(i,j,l)

                   0    0.422738    5.906928   -1.261894
                                   -1.261894    3.258196
                   1    0.484278    2.727013


 MOLECULE KIND INFORMATION


 All atoms are their own molecule, skipping detailed information


 TOTAL NUMBERS AND MAXIMUM NUMBERS

  Total number of            - Atomic kinds:                                   1
                             - Atoms:                                          1
                             - Shell sets:                                     2
                             - Shells:                                         5
                             - Primitive Cartesian functions:                  5
                             - Cartesian basis functions:                     14
                             - Spherical basis functions:                     13

  Maximum angular momentum of- Orbital basis functions:                        2
                             - Local part of the GTH pseudopotential:          0
                             - Non-local part of the GTH pseudopotential:      2


 MODULE QUICKSTEP:  ATOMIC COORDINATES IN angstrom

  Atom  Kind  Element       X           Y           Z          Z(eff)       Mass

       1     1 Si  14    0.000000    0.000000    0.000000      4.00      28.0855




 SCF PARAMETERS         Density guess:                                    ATOMIC
                        --------------------------------------------------------
                        max_scf:                                             100
                        max_scf_history:                                       0
                        max_diis:                                              4
                        --------------------------------------------------------
                        eps_scf:                                        1.00E-03
                        eps_scf_history:                                0.00E+00
                        eps_diis:                                       1.00E-01
                        eps_eigval:                                     1.00E-05
                        --------------------------------------------------------
                        level_shift [a.u.]:                                 0.00
                        --------------------------------------------------------
                        Mixing method:                           DIRECT_P_MIXING
                        --------------------------------------------------------
                        No outer SCF

 PW_GRID| Information for grid number                                          1
 PW_GRID| Cutoff [a.u.]                                                     25.0
 PW_GRID| spherical cutoff:                                                   NO
 PW_GRID|   Bounds   1            -12      11                Points:          24
 PW_GRID|   Bounds   2            -12      11                Points:          24
 PW_GRID|   Bounds   3            -12      11                Points:          24
 PW_GRID| Volume element (a.u.^3)  0.7819E-01     Volume (a.u.^3)      1080.8451
 PW_GRID| Grid span                                                    FULLSPACE

 PW_GRID| Information for grid number                                          2
 PW_GRID| Cutoff [a.u.]                                                      8.3
 PW_GRID| spherical cutoff:                                                   NO
 PW_GRID|   Bounds   1             -7       7                Points:          15
 PW_GRID|   Bounds   2             -7       7                Points:          15
 PW_GRID|   Bounds   3             -7       7                Points:          15
 PW_GRID| Volume element (a.u.^3)  0.3203         Volume (a.u.^3)      1080.8451
 PW_GRID| Grid span                                                    FULLSPACE

 POISSON| Solver                                                        PERIODIC
 POISSON| Periodicity                                                        XYZ

 RS_GRID| Information for grid number                                          1
 RS_GRID|   Bounds   1            -12      11                Points:          24
 RS_GRID|   Bounds   2            -12      11                Points:          24
 RS_GRID|   Bounds   3            -12      11                Points:          24

 RS_GRID| Information for grid number                                          2
 RS_GRID|   Bounds   1             -7       7                Points:          15
 RS_GRID|   Bounds   2             -7       7                Points:          15
 RS_GRID|   Bounds   3             -7       7                Points:          15

 DISTRIBUTION OF THE PARTICLES (ROWS)
              Process row      Number of particles         Number of matrix rows
                        0                        1                            -1
                      Sum                        1                            -1

 DISTRIBUTION OF THE PARTICLES (COLUMNS)
              Process col      Number of particles      Number of matrix columns
                        0                        1                            -1
                      Sum                        1                            -1

 DISTRIBUTION OF THE NEIGHBOR LISTS
              Total number of particle pairs:                                  7
              Total number of matrix elements:                              1183
              Average number of particle pairs:                                7
              Maximum number of particle pairs:                                7
              Average number of matrix element:                             1183
              Maximum number of matrix elements:                            1183


 DISTRIBUTION OF THE OVERLAP MATRIX
              Number  of non-zero blocks:                                      1
              Percentage non-zero blocks:                                 100.00
              Average number of blocks per CPU:                                1
              Maximum number of blocks per CPU:                                1
              Average number of matrix elements per CPU:                     179
              Maximum number of matrix elements per CPU:                     179

 Number of electrons:                                                          4
 Number of occupied orbitals:                                                  2
 Number of molecular orbitals:                                                 2

 Number of orbital functions:                                                 13
 Number of independent orbital functions:                                     13

 Extrapolation method: initial_guess

 Atomic guess: The first density matrix is obtained in terms of atomic orbitals
               and electronic configurations assigned to each atomic kind

 Guess for atomic kind: Si

 Electronic structure
    Total number of core electrons                                         10.00
    Total number of valence electrons                                       4.00
    Total number of electrons                                              14.00
    Multiplicity                                                   not specified
    S   [  2.00  2.00] 2.00
    P   [  6.00] 2.00


 *******************************************************************************
                  Iteration          Convergence                     Energy [au]
 *******************************************************************************
                          1        0.191310                      -3.618313869735
                          2        0.731569E-01                  -3.691159009622
                          3        0.405574E-02                  -3.699900512584
                          4        0.328704E-02                  -3.699908407293
                          5        0.320845E-02                  -3.699909118998
                          6        0.316809E-02                  -3.699909477757
                          7        0.331859E-05                  -3.699923449535
                          8        0.110258E-06                  -3.699923449550

 Energy components [Hartree]           Total Energy ::           -3.699923449550
                                        Band Energy ::           -1.012729790251
                                     Kinetic Energy ::            1.397012768229
                                   Potential Energy ::           -5.096936217779
                                      Virial (-V/T) ::            3.648453567279
                                        Core Energy ::           -5.703543362687
                                          XC Energy ::           -0.980691562795
                                     Coulomb Energy ::            2.984311475932
                       Total Pseudopotential Energy ::           -7.145739758818
                       Local Pseudopotential Energy ::           -7.987908627736
                    Nonlocal Pseudopotential Energy ::            0.842168868918
                                        Confinement ::            0.451836279031

 Orbital energies  State     L     Occupation   Energy[a.u.]          Energy[eV]

                       1     0          2.000      -0.378230          -10.292155

                       1     1          2.000      -0.128135           -3.486734

 Re-scaling the density matrix to get the right number of electrons
                  # Electrons              Trace(P)               Scaling factor
                            4                 3.957                        1.011


 SCF WAVEFUNCTION OPTIMIZATION

  Step     Update method      Time    Convergence         Total energy    Change
  ------------------------------------------------------------------------------

  Trace(PS):                                    4.0000000000
  Electronic density on regular grids:         -3.9998453719        0.0001546281
  Core density on regular grids:                3.9961648024       -0.0038351976
  Total charge density on r-space grids:       -0.0036805695
  Total charge density g-space grids:          -0.0036805695


 *************************************************************************
 *** 16:17:33 ERRORL2 in qs_ks_utils:calc_v_sic_rspace processor 0  :: ***
 *** err=-300 condition FAILED at line 764                             ***
 *************************************************************************


 ===== Routine Calling Stack ===== 

            8 qs_ks_build_kohn_sham_matrix
            7 rebuild_ks_matrix
            6 qs_ks_update_qs_env
            5 scf_env_do_scf_inner_loop
            4 scf_env_do_scf
            3 qs_energies_scf
            2 qs_forces
            1 CP2K