← Back to branch summary

~javier-junquera/siesta/lda+u+so

« back to all changes in this revision

Viewing changes to Tests/Reference/h2o_dipole2.out

  • Committer: Nick Papior
  • Date: 2018-11-09 09:36:46 UTC
  • mfrom: (560.1.476 4.1)
  • Revision ID: nickpapior@gmail.com-20181109093646-b9rx59dx9vyjict2
Merged r1024-1036, GR-pulay, OMM-reuse and test updates

Fixed GR-pulay, OMM-psi reuse.

Updated all tests from 4.1 and fixed input options for nearly
all tests.

Show diffs side-by-side

added added

removed removed

Lines of Context:
Tests/Reference/h2o_dipole2.out
1
 
 
2
 
                           ***********************       
3
 
                           *  WELCOME TO SIESTA  *       
4
 
                           ***********************       
5
 
 
6
 
reinit: Reading from ../h2o_dipole2.fdf
7
 
 
8
 
reinit: -----------------------------------------------------------------------
9
 
reinit: System Name: Water molecule
10
 
reinit: -----------------------------------------------------------------------
11
 
reinit: System Label: h2o_dipole2
12
 
reinit: -----------------------------------------------------------------------
13
 
Siesta Version: siesta-4.1--731
14
 
Architecture  : x86_64-linux-n-62-18-14
15
 
Compiler flags: mpifort -m64 -fPIC -O3 -march=native -ftree-vectorize -fexpensive-optimizatioons -funroll-loops -fprefetch-loop-arrays -fno-second-underscore  -flto -fuse-linker-plugin
16
 
PP flags      : -I/zdata/groups/common/nicpa/2016-oct/generic/build-tools/1.0/include -I/zdatta/groups/common/nicpa/2016-oct/generic/gmp/6.1.0/include -I/zdata/groups/common/nicpa/2016-oct/generic/mpfr/3.1.4/include -I/zdata/groups/common/nicpa/2016-oct/generic/mpc/1.0.3/include -I/zdata/groups/common/nicpa/2016-oct/generic/isl/0.16.1/include -I/zdata/groups/common/nicpa/2016-oct/generic/gcc/6.2.0/include -I/zdata/groups/common/nicpa/2016-oct/generic/numactl/2.0.11/include -I/zdata/groups/common/nicpa/2016-oct/XeonE5-2665/zlib/1.2.8/gnu-6.2.0/include -I/zdata/groups/common/nicpa/2016-oct/XeonE5-2665/libxml2/2.9.4/gnu-6.2.0/include -I/zdata/groups/common/nicpa/2016-oct/XeonE5-2665/hwloc/1.11.4/gnu-6.2.0/include -I/zdata/groups/common/nicpa/2016-oct/XeonE5-2665/openmpi/2.0.1/gnu-6.2.0/include -I/zdata/groups/common/nicpa/2016-oct/XeonE5-2665/hdf5/1.8.17/gnu-6.2.0/include -I/zdata/groups/common/nicpa/2016-oct/XeonE5-2665/pnetcdf/1.7.0/gnu-6.2.0/include -I/zdata/groups/common/nicpa/2016-oct/XeonE5-2665/netcdf/4.4.1/gnu-6.2.0/include -I/zdata/groups/common/nicpa/2016-oct/XeonE5-2665/flook/0.7.0/gnu-6.2.0/include -I/zdata/groups/common/nicpa/2016-oct/XeonE5-2665/parmetis/4.0.3/gnu-6.2.0/include -I/zdata/groups/common/nicpa/2016-oct/XeonE5-2665/scalapack/204/gnu-6.2.0/include -I/zdata/groups/common/nicpa/2016-oct/XeonE5-2665/openblas/0.2.19/gnu-6.2.0/include -I/zdata/groups/common/nicpa/2016-oct/XeonE5-2665/mumps/5.0.2/gnu-6.2.0/include -I/zdata/groups/common/nicpa/2016-oct/XeonE5-2665/fftw/3.3.5/gnu-6.2.0/include -DMPI -DFC_HAVE_FLUSH -DFC_HAVE_ABORT -DCDF -DCDF4 -DSIESTA__METIS -DSIESTA__MUMPS -DTS_NOCHECKS -DNCDF -DNCDF_4 -DNCDF_PARALLEL -DSIESTA__FLOOK
17
 
Libraries     : -lnetcdff -lnetcdf -lpnetcdf -lhdf5_hl -lhdf5 -lz  -lzmumps -lmumps_common -llpord -lparmetis -lmetis -lscalapack -lopenblas  -lmetis -flto -fuse-linker-plugin  -L/zdata/groups/common/nicpa/2016-oct/XeonE5-2665/flook/0.7.0/gnu-6.2.0/lib -Wl,-rpath=/zdata/groups/common/nicpa/2016-oct/XeonE5-2665/flook/0.7.0/gnu-6.2.0/lib -lflookall -ldl
 
1
Siesta Version  : siesta-4.1-1033
 
2
Architecture    : x86_64-linux-n-62-18-18
 
3
Compiler version: GNU Fortran (GCC) 7.3.0
 
4
Compiler flags  : mpifort -m64 -fPIC -O3 -march=native -ftree-vectorize -fexpensive-optimizatioons -funroll-loops -fprefetch-loop-arrays -fno-second-underscore  -flto
 
5
PP flags        : -I/zdata/groups/common/nicpa/2018-feb/generic/build-tools/1.0/include -I/zdatta/groups/common/nicpa/2018-feb/generic/gcc/7.3.0/include -I/zdata/groups/common/nicpa/2018-feb/generic/gcc/7.3.0/include -I/zdata/groups/common/nicpa/2018-feb/XeonE5-2665/zlib/1.2.11/gnu-7.3.0/include -I/zdata/groups/common/nicpa/2018-feb/generic/numactl/2.0.11/include -I/zdata/groups/common/nicpa/2018-feb/XeonE5-2665/libxml2/2.9.7/gnu-7.3.0/include -I/zdata/groups/common/nicpa/2018-feb/XeonE5-2665/hwloc/1.11.9/gnu-7.3.0/include -I/zdata/groups/common/nicpa/2018-feb/XeonE5-2665/openmpi/3.0.0/gnu-7.3.0/include -I/zdata/groups/common/nicpa/2018-feb/XeonE5-2665/szip/2.1.1/gnu-7.3.0/include -I/zdata/groups/common/nicpa/2018-feb/XeonE5-2665/hdf5/1.8.18/gnu-7.3.0/include -I/zdata/groups/common/nicpa/2018-feb/XeonE5-2665/pnetcdf/1.8.1/gnu-7.3.0/include -I/zdata/groups/common/nicpa/2018-feb/XeonE5-2665/netcdf/4.6.0/gnu-7.3.0/include -I/zdata/groups/common/nicpa/2018-feb/XeonE5-2665/flook/0.7.0/gnu-7.3.0/include -I/zdata/groups/common/nicpa/2018-feb/XeonE5-2665/parmetis/4.0.3/gnu-7.3.0/include -I/zdata/groups/common/nicpa/2018-feb/XeonE5-2665/scalapack/204/gnu-7.3.0/include -I/zdata/groups/common/nicpa/2018-feb/XeonE5-2665/openblas/0.2.20/gnu-7.3.0/include -I/zdata/groups/common/nicpa/2018-feb/XeonE5-2665/scotch/6.0.4/gnu-7.3.0/include -I/zdata/groups/common/nicpa/2018-feb/XeonE5-2665/mumps/5.1.2/gnu-7.3.0/include -I/zdata/groups/common/nicpa/2018-feb/XeonE5-2665/fftw/3.3.7/gnu-7.3.0/include -I/zdata/groups/common/nicpa/2018-feb/XeonE5-2665/elpa/2017.05.003/gnu-7.3.0/include/elpa -DSIESTA__ELPA -DMPI -D1 -DFC_HAVE_ABORT -DCDF -DCDF4 -DSIESTA__FLOOK  -DNCDF -DNCDF_4 -DNCDF_PARALLEL -DSIESTA__METIS -DSIESTA__MUMPS -DTS_NOCHECKS -DSIESTA__MRRR
 
6
Libraries       : libncdf.a libfdict.a -lnetcdff -lnetcdf -lpnetcdf -lhdf5_hl -lhdf5 -lz  -lzmuumps -lmumps_common -lesmumps -lscotch -lscotcherr -lpord -lparmetis -lmetis -lelpa -lscalapack -lopenblas  -L/zdata/groups/common/nicpa/2018-feb/XeonE5-2665/flook/0.7.0/gnu-7.3.0/lib -Wl,-rpath=/zdata/groups/common/nicpa/2018-feb/XeonE5-2665/flook/0.7.0/gnu-7.3.0/lib -lflookall -ldl  -flto -fuse-linker-plugin  -lmetis
18
7
PARALLEL version
19
8
NetCDF support
20
9
NetCDF-4 support
22
11
METIS ordering support
23
12
 
24
13
* Running on 8 nodes in parallel
25
 
>> Start of run:   2-JUL-2017  11:40:41
 
14
>> Start of run:   5-NOV-2018  10:49:29
 
15
 
 
16
                           ***********************       
 
17
                           *  WELCOME TO SIESTA  *       
 
18
                           ***********************       
 
19
 
 
20
reinit: Reading from ../h2o_dipole2.fdf
 
21
 
 
22
reinit: -----------------------------------------------------------------------
 
23
reinit: System Name: Water molecule
 
24
reinit: -----------------------------------------------------------------------
 
25
reinit: System Label: h2o_dipole2
 
26
reinit: -----------------------------------------------------------------------
26
27
 
27
28
initatom: Reading input for the pseudopotentials and atomic orbitals ----------
28
29
Species number:   1 Atomic number:    8 Label: O
341
342
redata: Number of Atomic Species                    =        2
342
343
redata: Charge density info will appear in .RHO file
343
344
redata: Write Mulliken Pop.                         = NO
 
345
redata: Matel table size (NRTAB)                    =     1024
344
346
redata: Mesh Cutoff                                 =   200.0000 Ry
345
347
redata: Net charge of the system                    =     0.0000 |e|
346
348
redata: Min. number of SCF Iter                     =        0
347
 
redata: Max. number of SCF Iter                     =       50
 
349
redata: Max. number of SCF Iter                     =     1000
 
350
redata: SCF convergence failure will abort job
348
351
redata: SCF mix quantity                            = Hamiltonian
349
352
redata: Mix DM or H after convergence               =   F
350
353
redata: Recompute H after scf cycle                 =   F
357
360
redata: Require Harris convergence for SCF          =   F
358
361
redata: Harris energy tolerance for SCF             =     0.000100 eV
359
362
redata: Require DM convergence for SCF              =   T
360
 
redata: DM tolerance for SCF                        =     0.0001
 
363
redata: DM tolerance for SCF                        =     0.000100
361
364
redata: Require EDM convergence for SCF             =   F
362
365
redata: EDM tolerance for SCF                       =     0.001000 eV
363
366
redata: Require H convergence for SCF               =   T
368
371
redata: Use continuation files for DM               =   F
369
372
redata: Neglect nonoverlap interactions             =   F
370
373
redata: Method of Calculation                       = Diagonalization
371
 
redata: Divide and Conquer                          =   T
372
374
redata: Electronic Temperature                      =   299.9869 K
373
375
redata: Fix the spin of the system                  =   F
374
376
redata: Dynamics option                             = Single-point calculation
415
417
siesta: k-grid:    0   1   0      0.000
416
418
siesta: k-grid:    0   0   1      0.000
417
419
 
 
420
diag: Algorithm                                     = D&C
 
421
diag: Parallel over k                               =   F
 
422
diag: Use parallel 2D distribution                  =   T
 
423
diag: Parallel block-size                           = 6
 
424
diag: Parallel distribution                         =     2 x     4
 
425
diag: Used triangular part                          = Lower
 
426
diag: Absolute tolerance                            =  0.100E-15
 
427
diag: Orthogonalization factor                      =  0.100E-05
 
428
diag: Memory factor                                 =  1.0000
 
429
 
 
430
 
 
431
ts: **************************************************************
 
432
ts: Save H and S matrices                           =    F
 
433
ts: Save DM and EDM matrices                        =    F
 
434
ts: Fix Hartree potential                           =    F
 
435
ts: Only save the overlap matrix S                  =    F
 
436
ts: **************************************************************
 
437
 
 
438
************************ Begin: TS CHECKS AND WARNINGS ************************
 
439
************************ End: TS CHECKS AND WARNINGS **************************
 
440
 
418
441
 
419
442
                     ====================================
420
443
                        Single-point calculation
428
451
outcell: Cell vector modules (Ang)   :    3.000000   12.000000    5.000000
429
452
outcell: Cell angles (23,13,12) (deg):     90.0000     90.0000     90.0000
430
453
outcell: Cell volume (Ang**3)        :    180.0000
431
 
 Folding of H and S is implicitly performed
 
454
Gamma-point calculation with interaction between periodic images
 
455
Some features might not work optimally:
 
456
e.g. DM initialization from atomic data
432
457
<dSpData1D:S at geom step 0
433
458
  <sparsity:sparsity for geom step 0
434
459
    nrows_g=46 nrows=6 sparsity=.0652 nnzs=138, refcount: 7>
442
467
    nrows_g=46 nrows=6 sparsity=.0652 nnzs=138, refcount: 8>
443
468
  <dData2D:DM n=138 m=1, refcount: 1>
444
469
refcount: 1>
445
 
 
446
 
Note: For starting DM, Qtot, Tr[D*S] =         16.00000000         16.18353818
447
470
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps:       3      45
448
471
New grid distribution:   1
449
472
           1       1:   15    1:   27    1:    6
493
516
Electric field for dipole correction =    0.000000   -0.000000    0.000000 eV/Ang/e
494
517
 
495
518
siesta: Program's energy decomposition (eV):
496
 
siesta: Ebs     =      -162.566207
 
519
siesta: Ebs     =      -179.640045
497
520
siesta: Eions   =      1631.708957
498
 
siesta: Ena     =       349.727602
499
 
siesta: Ekin    =       775.005215
500
 
siesta: Enl     =      -144.531527
 
521
siesta: Ena     =       349.730334
 
522
siesta: Ekin    =       756.447513
 
523
siesta: Enl     =      -139.624038
501
524
siesta: Eso     =         0.000000
502
525
siesta: Eldau   =         0.000000
503
 
siesta: DEna    =       -43.526448
504
 
siesta: DUscf   =         6.468271
 
526
siesta: DEna    =       -32.649064
 
527
siesta: DUscf   =         4.641921
505
528
siesta: DUext   =         0.000000
506
 
siesta: Exc     =      -237.279293
 
529
siesta: Enegf   =         0.000000
 
530
siesta: Exc     =      -234.322959
507
531
siesta: eta*DQ  =         0.000000
508
532
siesta: Emadel  =         0.000000
509
533
siesta: Emeta   =         0.000000
510
534
siesta: Emolmec =         0.000000
511
535
siesta: Ekinion =         0.000000
512
 
siesta: Eharris =      -917.542515
513
 
siesta: Etot    =      -925.845138
514
 
siesta: FreeEng =      -925.845138
 
536
siesta: Eharris =      -923.358033
 
537
siesta: Etot    =      -927.485251
 
538
siesta: FreeEng =      -927.485251
515
539
 
516
540
        iscf     Eharris(eV)        E_KS(eV)     FreeEng(eV)     dDmax    Ef(eV) dHmax(eV)
517
 
   scf:    1     -917.542515     -925.845138     -925.845138  1.426479 -3.165748  8.445312
518
 
timer: Routine,Calls,Time,% = IterSCF        1       0.079   5.74
519
 
Dipole moment in unit cell   =      0.0000     -0.0000      0.0000 D
520
 
Electric field for dipole correction =   -0.000000    0.000000   -0.000000 eV/Ang/e
521
 
   scf:    2     -863.633215     -908.882108     -908.882108  0.299113 -2.847056 18.264538
522
 
Dipole moment in unit cell   =     -0.0000      0.0000     -0.0000 D
523
 
Electric field for dipole correction =    0.000000   -0.000000    0.000000 eV/Ang/e
524
 
   scf:    3     -933.417105     -929.208963     -929.208963  0.296710 -5.728868  1.167038
525
 
Dipole moment in unit cell   =     -0.0000      0.0000     -0.0000 D
526
 
Electric field for dipole correction =    0.000000   -0.000000    0.000000 eV/Ang/e
527
 
   scf:    4     -929.321121     -929.282340     -929.282340  0.011650 -5.530886  0.272993
528
 
Dipole moment in unit cell   =     -0.0000      0.0000     -0.0000 D
529
 
Electric field for dipole correction =    0.000000   -0.000000    0.000000 eV/Ang/e
530
 
   scf:    5     -929.292844     -929.288281     -929.288281  0.012652 -5.479732  0.127470
531
 
Dipole moment in unit cell   =     -0.0000      0.0000     -0.0000 D
532
 
Electric field for dipole correction =    0.000000   -0.000000    0.000000 eV/Ang/e
533
 
   scf:    6     -929.292279     -929.290368     -929.290368  0.009946 -5.456889  0.072816
534
 
Dipole moment in unit cell   =     -0.0000      0.0000     -0.0000 D
535
 
Electric field for dipole correction =    0.000000   -0.000000    0.000000 eV/Ang/e
536
 
   scf:    7     -929.291606     -929.291035     -929.291035  0.008408 -5.430449  0.024878
537
 
Dipole moment in unit cell   =     -0.0000      0.0000     -0.0000 D
538
 
Electric field for dipole correction =    0.000000   -0.000000    0.000000 eV/Ang/e
539
 
   scf:    8     -929.291112     -929.291075     -929.291075  0.002163 -5.416896  0.011514
540
 
Dipole moment in unit cell   =     -0.0000      0.0000     -0.0000 D
541
 
Electric field for dipole correction =    0.000000   -0.000000    0.000000 eV/Ang/e
542
 
   scf:    9     -929.291091     -929.291084     -929.291084  0.000711 -5.407729  0.007066
543
 
Dipole moment in unit cell   =     -0.0000      0.0000     -0.0000 D
544
 
Electric field for dipole correction =    0.000000   -0.000000    0.000000 eV/Ang/e
545
 
   scf:   10     -929.291099     -929.291092     -929.291092  0.000388 -5.407718  0.005467
546
 
Dipole moment in unit cell   =     -0.0000      0.0000     -0.0000 D
547
 
Electric field for dipole correction =    0.000000   -0.000000    0.000000 eV/Ang/e
548
 
   scf:   11     -929.291101     -929.291097     -929.291097  0.000545 -5.408396  0.002918
549
 
Dipole moment in unit cell   =     -0.0000      0.0000     -0.0000 D
550
 
Electric field for dipole correction =    0.000000   -0.000000    0.000000 eV/Ang/e
551
 
   scf:   12     -929.291097     -929.291097     -929.291097  0.000138 -5.409054  0.001783
552
 
Dipole moment in unit cell   =     -0.0000      0.0000     -0.0000 D
553
 
Electric field for dipole correction =    0.000000   -0.000000    0.000000 eV/Ang/e
554
 
   scf:   13     -929.291097     -929.291097     -929.291097  0.000175 -5.410233  0.000518
555
 
Dipole moment in unit cell   =     -0.0000      0.0000     -0.0000 D
556
 
Electric field for dipole correction =    0.000000   -0.000000    0.000000 eV/Ang/e
557
 
   scf:   14     -929.291097     -929.291097     -929.291097  0.000010 -5.410468  0.000324
 
541
   scf:    1     -923.358033     -927.485251     -927.485251  1.435109 -2.746971  6.169673
 
542
timer: Routine,Calls,Time,% = IterSCF        1       0.066   3.29
 
543
Dipole moment in unit cell   =     -0.0000      0.0000     -0.0000 D
 
544
Electric field for dipole correction =    0.000000   -0.000000    0.000000 eV/Ang/e
 
545
   scf:    2     -930.156854     -929.119461     -929.119461  0.041016 -5.958812  1.893694
 
546
Dipole moment in unit cell   =     -0.0000      0.0000     -0.0000 D
 
547
Electric field for dipole correction =    0.000000   -0.000000    0.000000 eV/Ang/e
 
548
   scf:    3     -929.315633     -929.279298     -929.279298  0.015906 -5.556998  0.212174
 
549
Dipole moment in unit cell   =     -0.0000      0.0000     -0.0000 D
 
550
Electric field for dipole correction =    0.000000   -0.000000    0.000000 eV/Ang/e
 
551
   scf:    4     -929.289579     -929.284805     -929.284805  0.009538 -5.536502  0.188299
 
552
Dipole moment in unit cell   =     -0.0000      0.0000     -0.0000 D
 
553
Electric field for dipole correction =    0.000000   -0.000000    0.000000 eV/Ang/e
 
554
   scf:    5     -929.288797     -929.286889     -929.286889  0.007090 -5.506711  0.145216
 
555
Dipole moment in unit cell   =     -0.0000      0.0000     -0.0000 D
 
556
Electric field for dipole correction =    0.000000   -0.000000    0.000000 eV/Ang/e
 
557
   scf:    6     -929.289280     -929.288212     -929.288212  0.013811 -5.423575  0.038152
 
558
Dipole moment in unit cell   =     -0.0000      0.0000     -0.0000 D
 
559
Electric field for dipole correction =    0.000000   -0.000000    0.000000 eV/Ang/e
 
560
   scf:    7     -929.288263     -929.288244     -929.288244  0.001113 -5.405724  0.017906
 
561
Dipole moment in unit cell   =     -0.0000      0.0000     -0.0000 D
 
562
Electric field for dipole correction =    0.000000   -0.000000    0.000000 eV/Ang/e
 
563
   scf:    8     -929.288308     -929.288277     -929.288277  0.000684 -5.405544  0.011004
 
564
Dipole moment in unit cell   =     -0.0000      0.0000     -0.0000 D
 
565
Electric field for dipole correction =    0.000000   -0.000000    0.000000 eV/Ang/e
 
566
   scf:    9     -929.288311     -929.288296     -929.288296  0.000822 -5.406694  0.005690
 
567
Dipole moment in unit cell   =     -0.0000      0.0000     -0.0000 D
 
568
Electric field for dipole correction =    0.000000   -0.000000    0.000000 eV/Ang/e
 
569
   scf:   10     -929.288299     -929.288298     -929.288298  0.000274 -5.408054  0.003093
 
570
Dipole moment in unit cell   =     -0.0000      0.0000     -0.0000 D
 
571
Electric field for dipole correction =    0.000000   -0.000000    0.000000 eV/Ang/e
 
572
   scf:   11     -929.288298     -929.288298     -929.288298  0.000193 -5.409685  0.001173
 
573
Dipole moment in unit cell   =     -0.0000      0.0000     -0.0000 D
 
574
Electric field for dipole correction =    0.000000   -0.000000    0.000000 eV/Ang/e
 
575
   scf:   12     -929.288298     -929.288298     -929.288298  0.000026 -5.410388  0.000820
558
576
 
559
577
SCF Convergence by DM+H criterion
560
 
max |DM_out - DM_in|         :     0.0000095834
561
 
max |H_out - H_in|      (eV) :     0.0003244212
562
 
SCF cycle converged after 14 iterations
 
578
max |DM_out - DM_in|         :     0.0000263046
 
579
max |H_out - H_in|      (eV) :     0.0008196475
 
580
SCF cycle converged after 12 iterations
563
581
 
564
582
Using DM_out to compute the final energy and forces
565
583
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps:       3      45
566
584
Dipole moment in unit cell   =     -0.0000      0.0000     -0.0000 D
567
585
Electric field for dipole correction =    0.000000   -0.000000    0.000000 eV/Ang/e
568
586
 
569
 
siesta: E_KS(eV) =             -929.2911
 
587
siesta: E_KS(eV) =             -929.2883
570
588
 
571
 
siesta: E_KS - E_eggbox =      -929.2911
 
589
siesta: E_KS - E_eggbox =      -929.2883
572
590
 
573
591
siesta: Atomic forces (eV/Ang):
574
592
----------------------------------------
575
 
   Tot   -0.000000    0.000000   -0.000000
576
 
----------------------------------------
577
 
   Max    1.048693
578
 
   Res    0.603481    sqrt( Sum f_i^2 / 3N )
579
 
----------------------------------------
580
 
   Max    1.048693    constrained
 
593
   Tot    0.000000    0.000000    0.000000
 
594
----------------------------------------
 
595
   Max    1.076128
 
596
   Res    0.623852    sqrt( Sum f_i^2 / 3N )
 
597
----------------------------------------
 
598
   Max    1.076128    constrained
581
599
 
582
 
Stress-tensor-Voigt (kbar):      -79.18        8.75       -1.16        0.00        0.00       -0.00
583
 
(Free)E + p*V (eV/cell)     -926.6104
584
 
Target enthalpy (eV/cell)     -929.2911
 
600
Stress-tensor-Voigt (kbar):      -78.42        9.19       -1.16       -0.00       -0.00        0.00
 
601
(Free)E + p*V (eV/cell)     -926.6522
 
602
Target enthalpy (eV/cell)     -929.2883
585
603
 
586
604
siesta: Program's energy decomposition (eV):
587
 
siesta: Ebs     =      -219.474983
 
605
siesta: Ebs     =      -219.476393
588
606
siesta: Eions   =      1631.708957
589
 
siesta: Ena     =       349.727602
590
 
siesta: Ekin    =       711.300185
591
 
siesta: Enl     =      -127.077889
 
607
siesta: Ena     =       349.730334
 
608
siesta: Ekin    =       711.300778
 
609
siesta: Enl     =      -127.078745
592
610
siesta: Eso     =         0.000000
593
611
siesta: Eldau   =         0.000000
594
 
siesta: DEna    =        -6.219236
595
 
siesta: DUscf   =         1.647317
 
612
siesta: DEna    =        -6.219103
 
613
siesta: DUscf   =         1.647223
596
614
siesta: DUext   =         0.000000
597
 
siesta: Exc     =      -226.960120
 
615
siesta: Enegf   =         0.000000
 
616
siesta: Exc     =      -226.959827
598
617
siesta: eta*DQ  =         0.000000
599
618
siesta: Emadel  =         0.000000
600
619
siesta: Emeta   =         0.000000
601
620
siesta: Emolmec =         0.000000
602
621
siesta: Ekinion =         0.000000
603
 
siesta: Eharris =      -929.291097
604
 
siesta: Etot    =      -929.291097
605
 
siesta: FreeEng =      -929.291097
 
622
siesta: Eharris =      -929.288298
 
623
siesta: Etot    =      -929.288298
 
624
siesta: FreeEng =      -929.288298
606
625
 
607
626
siesta: Final energy (eV):
608
 
siesta:  Band Struct. =    -219.474983
609
 
siesta:       Kinetic =     711.300185
610
 
siesta:       Hartree =     709.568373
 
627
siesta:  Band Struct. =    -219.476393
 
628
siesta:       Kinetic =     711.300778
 
629
siesta:       Hartree =     709.566801
611
630
siesta:       Eldau   =       0.000000
612
631
siesta:       Eso     =       0.000000
613
632
siesta:    Ext. field =       0.000000
614
 
siesta:   Exch.-corr. =    -226.960120
615
 
siesta:  Ion-electron =   -2052.537441
616
 
siesta:       Ion-ion =     -70.662094
 
633
siesta:       Enegf   =       0.000000
 
634
siesta:   Exch.-corr. =    -226.959827
 
635
siesta:  Ion-electron =   -2052.536687
 
636
siesta:       Ion-ion =     -70.659362
617
637
siesta:       Ekinion =       0.000000
618
 
siesta:         Total =    -929.291097
619
 
siesta:         Fermi =      -5.410468
 
638
siesta:         Total =    -929.288298
 
639
siesta:         Fermi =      -5.410388
620
640
 
621
641
siesta: Atomic forces (eV/Ang):
622
 
siesta:      1    0.000000    0.844799    0.000000
623
 
siesta:      2   -1.048693   -0.426910   -0.000000
624
 
siesta:      3    1.048693   -0.426910   -0.000000
625
 
siesta:      4   -0.000000   -0.844799    0.000000
626
 
siesta:      5   -1.048693    0.426910   -0.000000
627
 
siesta:      6    1.048693    0.426910   -0.000000
 
642
siesta:      1   -0.000000    0.886408    0.000000
 
643
siesta:      2   -1.076128   -0.447716   -0.000000
 
644
siesta:      3    1.076128   -0.447716   -0.000000
 
645
siesta:      4    0.000000   -0.886408    0.000000
 
646
siesta:      5   -1.076128    0.447716   -0.000000
 
647
siesta:      6    1.076128    0.447716   -0.000000
628
648
siesta: ----------------------------------------
629
 
siesta:    Tot   -0.000000    0.000000   -0.000000
 
649
siesta:    Tot    0.000000    0.000000    0.000000
630
650
 
631
651
siesta: Stress tensor (static) (eV/Ang**3):
632
 
siesta:    -0.049420    0.000000    0.000000
633
 
siesta:     0.000000    0.005463    0.000000
634
 
siesta:    -0.000000    0.000000   -0.000721
 
652
siesta:    -0.048948   -0.000000    0.000000
 
653
siesta:    -0.000000    0.005734   -0.000000
 
654
siesta:     0.000000   -0.000000   -0.000721
635
655
 
636
656
siesta: Cell volume =        180.000000 Ang**3
637
657
 
638
658
siesta: Pressure (static):
639
659
siesta:                Solid            Molecule  Units
640
 
siesta:           0.00016220          0.00024534  Ry/Bohr**3
641
 
siesta:           0.01489263          0.02252596  eV/Ang**3
642
 
siesta:          23.86087381         36.09095324  kBar
643
 
(Free)E+ p_basis*V_orbitals  =        -928.120763
644
 
(Free)Eharris+ p_basis*V_orbitals  =        -928.120763
 
660
siesta:           0.00015950          0.00024530  Ry/Bohr**3
 
661
siesta:           0.01464509          0.02252254  eV/Ang**3
 
662
siesta:          23.46427262         36.08547403  kBar
 
663
(Free)E+ p_basis*V_orbitals  =        -928.117963
 
664
(Free)Eharris+ p_basis*V_orbitals  =        -928.117963
645
665
 
646
666
siesta: Electric dipole (a.u.)  =   -0.000000    0.000000   -0.000000
647
667
siesta: Electric dipole (Debye) =   -0.000000    0.000000   -0.000000
648
668
Dipole moment in unit cell   =     -0.0000      0.0000     -0.0000 D
649
669
Electric field for dipole correction =    0.000000   -0.000000    0.000000 eV/Ang/e
650
670
 
651
 
dhscf: Vacuum level (max, mean) =    1.118146    0.439745 eV
652
 
>> End of run:   2-JUL-2017  11:40:43
 
671
dhscf: Vacuum level (max, mean) =    1.118332    0.439821 eV
 
672
>> End of run:   5-NOV-2018  10:49:31
653
673
Job completed