
  ___ ___ ___ _ _ _  
 |   |   |_  | | | | 
 | | | | | . | | | | 
 |__ |  _|___|_____|  23.6.2b1
 |___|_|             

User:   aes38@compute-25
Date:   Thu Sep 28 11:55:24 2023
Arch:   x86_64
Pid:    1050453
CWD:    /users/home/aes38/Rydberg/new/ethylene/sicfromreal/seventhext
Python: 3.9.6
gpaw:   /users/home/aes38/new-gmf-gpaw/gpaw (7a0e975e9f)
_gpaw:  /users/home/aes38/new-gmf-gpaw/build/lib.linux-x86_64-3.9/
        _gpaw.cpython-39-x86_64-linux-gnu.so (7a0e975e9f)
ase:    /hpcapps/source/ase_dev/ase (version 3.23.0b1-d4c62f2276)
numpy:  /hpcapps/libsci-gcc/software/Python/3.9.6-GCCcore-11.2.0/lib/python3.9/site-packages/numpy (version 1.22.4)
scipy:  /hpcapps/libsci-gcc/software/Python/3.9.6-GCCcore-11.2.0/lib/python3.9/site-packages/scipy (version 1.8.1)
libxc:  5.1.6
units:  Angstrom and eV
cores: 48
OpenMP: False
OMP_NUM_THREADS: 1

Reading from /users/home/aes38/Rydberg/new/ethylene/Real/firstext/triplet.gpw
Read dipole, energy, free_energy, magmom, magmoms
Reading input parameters:
  basis: aug-cc-pVDZ_PBE.sz
  eigensolver: {blocksize: 1,
                converge_unocc: False,
                excited_state: True,
                functional: {'name': 'ks'},
                grad_tol_inner_loop: 0.0005,
                grad_tol_pz_localization: 0.0005,
                linesearch_algo: {'name': 'max-step', 'max_step': 0.2},
                localization_tol: None,
                localizationseed: None,
                localizationtype: None,
                max_step_inner_loop: 0.2,
                maxiter_inner_loop: 100,
                maxiter_pz_localization: 50,
                maxiter_unocc: 333,
                momevery: 10,
                name: etdm-fdpw,
                need_init_orbs: False,
                need_localization: False,
                printinnerloop: False,
                restart_canonical: True,
                restartevery_iloop_notconverged: 3,
                searchdir_algo: {'name': 'l-bfgs', 'memory': 3},
                use_prec: True}
  h: 0.15
  maxiter: 1000
  mixer: {backend: no-mixing}
  mode: {interpolation: 3,
         name: fd,
         nn: 3}
  nbands: -10
  occupations: {name: mom,
                numbers: [[1. 1. 1. ... 0. 0. 0.]
 [1. 1. 1. ... 0. 0. 0.]],
                update_numbers: True,
                use_fixed_occupations: False,
                use_projections: True}
  spinpol: True
  symmetry: off
  xc: PBE

Initialize ...

species:
  C:
    name: Carbon
    id: 4aa54d4b901d75f77cc0ea3eec22967b
    Z: 6.0
    valence: 4
    core: 2
    charge: 0.0
    file: /users/home/share/gpaw-setups-0.9.20000/C.PBE.gz
    compensation charges: {type: gauss,
                           rc: 0.20,
                           lmax: 2}
    cutoffs: {filter: 1.14,
              core: 1.14}
    valence states:
      #              energy  rcut
      - 2s(2.00)   -13.751   0.635
      - 2p(2.00)    -5.284   0.635
      -  s          13.461   0.635
      -  p          21.927   0.635
      -  d           0.000   0.635
  
    # LCAO basis set for C:
    #   Name: aug-cc-pVDZ_PBE.sz
    #   File: /users/home/aes38/basissets/C.aug-cc-pVDZ_PBE.sz.basis
    #   Number of radial functions: 8
    #   Number of spherical harmonics: 22
    #     l=0, rc=5.5938 Bohr: 2s-sz confined orbital
    #     l=1, rc=6.8750 Bohr: 2p-sz confined orbital
    #     l=0, rc=7.5625 Bohr: NGTO(s,0.1596000)
    #     l=0, rc=13.1094 Bohr: NGTO(s,0.0469000)
    #     l=1, rc=7.7344 Bohr: NGTO(p,0.1517000)
    #     l=1, rc=13.9531 Bohr: NGTO(p,0.0404100)
    #     l=2, rc=4.4688 Bohr: NGTO(d,0.5500000)
    #     l=2, rc=8.0469 Bohr: NGTO(d,0.1510000)
    #   Number of RI-basis functions 0

  H:
    name: Hydrogen
    id: d65de229564ff8ea4db303e23b6d1ecf
    Z: 1.0
    valence: 1
    core: 0
    charge: 0.0
    file: /users/home/share/gpaw-setups-0.9.20000/H.PBE.gz
    compensation charges: {type: gauss,
                           rc: 0.15,
                           lmax: 2}
    cutoffs: {filter: 0.85,
              core: 0.53}
    valence states:
      #              energy  rcut
      - 1s(1.00)    -6.494   0.476
      -  s          20.717   0.476
      -  p           0.000   0.476
  
    # LCAO basis set for H:
    #   Name: aug-cc-pVDZ_PBE.sz
    #   File: /users/home/aes38/basissets/H.aug-cc-pVDZ_PBE.sz.basis
    #   Number of radial functions: 5
    #   Number of spherical harmonics: 9
    #     l=0, rc=6.6406 Bohr: 1s-sz confined orbital
    #     l=0, rc=8.5469 Bohr: NGTO(s,0.1220000)
    #     l=0, rc=16.0469 Bohr: NGTO(s,0.0297400)
    #     l=1, rc=3.7969 Bohr: NGTO(p,0.7270000)
    #     l=1, rc=7.9844 Bohr: NGTO(p,0.1410000)
    #   Number of RI-basis functions 0

Reference energy: -2105.216880  # eV

Spin-polarized calculation.
Magnetic moment: 0.000000

Convergence criteria:
 Maximum [total energy] change in last 3 cyles: 0.0005 eV / valence electron
 Maximum integral of absolute [dens]ity change: 0.0001 electrons / valence electron
 Maximum integral of absolute [eigenst]ate change: 4e-08 eV^2 / valence electron
 Maximum number of scf [iter]ations: 1000
 (Square brackets indicate name in SCF output, whereas a 'c' in
 the SCF output indicates the quantity has converged.)

Symmetries present (total): 1

  ( 1  0  0)
  ( 0  1  0)
  ( 0  0  1)

1 k-point (Gamma)
1 k-point in the irreducible part of the Brillouin zone
       k-points in crystal coordinates                weights
   0:     0.00000000    0.00000000    0.00000000          1/1

Wave functions: Uniform real-space grid
  Kinetic energy operator: 6*3+1=19 point O(h^6) finite-difference Laplacian
  ScaLapack parameters: grid=1x1, blocksize=None
  Wavefunction extrapolation:
    Improved wavefunction reuse through dual PAW basis 

Excited-state calculation with Maximum Overlap Method
  Gaussian smearing of holes and excited electrons: off

Eigensolver
   Direct minimisation using exponential transformation.
       Search direction: L-BFGS algorithm
       Line search: step size equals one
       Preconditioning: True
       Orbital-density self-interaction corrections: ks
       WARNING: do not use it for metals as occupation numbers are
                not found variationally
 

Densities:
  Coarse grid: 132*144*148 grid
  Fine grid: 264*288*296 grid
  Total Charge: 0.000000 

Density mixing:
  Method: difference
  Backend: no-mixing
  Linear mixing parameter: 0
  old densities: 0
  Damping of long wavelength oscillations: 0 

Hamiltonian:
  XC and Coulomb potentials evaluated on a 264*288*296 grid
  Using the PBE Exchange-Correlation functional
  Interpolation: tri-quintic (5. degree polynomial)
  Poisson solver: FastPoissonSolver using
    Stencil: 6*3+1=19 point O(h^6) finite-difference Laplacian
    FFT axes: []
    FST axes: [0, 1, 2]
 

XC parameters: PBE with 2 nearest neighbor stencil

Memory estimate:
  Process memory now: 187.32 MiB
  Calculator: 41.57 MiB
    Density: 15.11 MiB
      Arrays: 15.00 MiB
      Localized functions: 0.10 MiB
      Mixer: 0.00 MiB
    Hamiltonian: 11.16 MiB
      Arrays: 11.15 MiB
      XC: 0.00 MiB
      Poisson: 0.00 MiB
      vbar: 0.01 MiB
    Wavefunctions: 15.30 MiB
      Arrays psit_nG: 13.23 MiB
      Eigensolver: 2.07 MiB
      Projections: 0.00 MiB
      Projectors: 0.01 MiB

Total number of cores used: 48
Domain decomposition: 3 x 4 x 4

Number of atoms: 6
Number of atomic orbitals: 80
Number of bands in calculation: 16
Number of valence electrons: 12
Bands to converge: occupied

... initialized

Input parameters:
  convergence: {eigenstates: 1e-06}
  eigensolver: {blocksize: 1,
                converge_unocc: False,
                excited_state: True,
                functional: {'name': 'PZ-SIC'},
                grad_tol_inner_loop: 0.0005,
                grad_tol_pz_localization: 0.0005,
                linesearch_algo: {'name': 'max-step', 'max_step': 0.2},
                localization_tol: None,
                localizationseed: None,
                localizationtype: PM,
                max_step_inner_loop: 0.25,
                maxiter_inner_loop: 200,
                maxiter_pz_localization: 50,
                maxiter_unocc: 333,
                momevery: 10,
                name: etdm-fdpw,
                need_init_orbs: False,
                need_localization: True,
                printinnerloop: True,
                restart_canonical: True,
                restartevery_iloop_notconverged: 3,
                searchdir_algo: {'name': 'l-bfgs', 'memory': 3},
                use_prec: True}
  mode: {force_complex_dtype: True,
         interpolation: 3,
         name: fd,
         nn: 3}

Initialize ...

species:
  C:
    name: Carbon
    id: 4aa54d4b901d75f77cc0ea3eec22967b
    Z: 6.0
    valence: 4
    core: 2
    charge: 0.0
    file: /users/home/share/gpaw-setups-0.9.20000/C.PBE.gz
    compensation charges: {type: gauss,
                           rc: 0.20,
                           lmax: 2}
    cutoffs: {filter: 1.14,
              core: 1.14}
    valence states:
      #              energy  rcut
      - 2s(2.00)   -13.751   0.635
      - 2p(2.00)    -5.284   0.635
      -  s          13.461   0.635
      -  p          21.927   0.635
      -  d           0.000   0.635
  
    # LCAO basis set for C:
    #   Name: aug-cc-pVDZ_PBE.sz
    #   File: /users/home/aes38/basissets/C.aug-cc-pVDZ_PBE.sz.basis
    #   Number of radial functions: 8
    #   Number of spherical harmonics: 22
    #     l=0, rc=5.5938 Bohr: 2s-sz confined orbital
    #     l=1, rc=6.8750 Bohr: 2p-sz confined orbital
    #     l=0, rc=7.5625 Bohr: NGTO(s,0.1596000)
    #     l=0, rc=13.1094 Bohr: NGTO(s,0.0469000)
    #     l=1, rc=7.7344 Bohr: NGTO(p,0.1517000)
    #     l=1, rc=13.9531 Bohr: NGTO(p,0.0404100)
    #     l=2, rc=4.4688 Bohr: NGTO(d,0.5500000)
    #     l=2, rc=8.0469 Bohr: NGTO(d,0.1510000)
    #   Number of RI-basis functions 0

  H:
    name: Hydrogen
    id: d65de229564ff8ea4db303e23b6d1ecf
    Z: 1.0
    valence: 1
    core: 0
    charge: 0.0
    file: /users/home/share/gpaw-setups-0.9.20000/H.PBE.gz
    compensation charges: {type: gauss,
                           rc: 0.15,
                           lmax: 2}
    cutoffs: {filter: 0.85,
              core: 0.53}
    valence states:
      #              energy  rcut
      - 1s(1.00)    -6.494   0.476
      -  s          20.717   0.476
      -  p           0.000   0.476
  
    # LCAO basis set for H:
    #   Name: aug-cc-pVDZ_PBE.sz
    #   File: /users/home/aes38/basissets/H.aug-cc-pVDZ_PBE.sz.basis
    #   Number of radial functions: 5
    #   Number of spherical harmonics: 9
    #     l=0, rc=6.6406 Bohr: 1s-sz confined orbital
    #     l=0, rc=8.5469 Bohr: NGTO(s,0.1220000)
    #     l=0, rc=16.0469 Bohr: NGTO(s,0.0297400)
    #     l=1, rc=3.7969 Bohr: NGTO(p,0.7270000)
    #     l=1, rc=7.9844 Bohr: NGTO(p,0.1410000)
    #   Number of RI-basis functions 0

Reference energy: -2105.216880  # eV

Spin-polarized calculation.
Magnetic moment: 0.000000

Convergence criteria:
 Maximum [total energy] change in last 3 cyles: 0.0005 eV / valence electron
 Maximum integral of absolute [dens]ity change: 0.0001 electrons / valence electron
 Maximum integral of absolute [eigenst]ate change: 1e-06 eV^2 / valence electron
 Maximum number of scf [iter]ations: 1000
 (Square brackets indicate name in SCF output, whereas a 'c' in
 the SCF output indicates the quantity has converged.)

Symmetries present (total): 1

  ( 1  0  0)
  ( 0  1  0)
  ( 0  0  1)

1 k-point (Gamma)
1 k-point in the irreducible part of the Brillouin zone
       k-points in crystal coordinates                weights
   0:     0.00000000    0.00000000    0.00000000          1/1

Wave functions: Uniform real-space grid
  Kinetic energy operator: 6*3+1=19 point O(h^6) finite-difference Laplacian
  ScaLapack parameters: grid=1x1, blocksize=None
  Wavefunction extrapolation:
    Improved wavefunction reuse through dual PAW basis 

Excited-state calculation with Maximum Overlap Method
  Gaussian smearing of holes and excited electrons: off

Eigensolver
   Direct minimisation using exponential transformation.
       Search direction: L-BFGS algorithm
       Line search: step size equals one
       Preconditioning: True
       Orbital-density self-interaction corrections: PZ-SIC
       WARNING: do not use it for metals as occupation numbers are
                not found variationally
 

Densities:
  Coarse grid: 132*144*148 grid
  Fine grid: 264*288*296 grid
  Total Charge: 0.000000 

Density mixing:
  Method: difference
  Backend: no-mixing
  Linear mixing parameter: 0
  old densities: 0
  Damping of long wavelength oscillations: 0 

Hamiltonian:
  XC and Coulomb potentials evaluated on a 264*288*296 grid
  Using the PBE Exchange-Correlation functional
  Interpolation: tri-quintic (5. degree polynomial)
  Poisson solver: FastPoissonSolver using
    Stencil: 6*3+1=19 point O(h^6) finite-difference Laplacian
    FFT axes: []
    FST axes: [0, 1, 2]
 

XC parameters: PBE with 2 nearest neighbor stencil

Memory estimate:
  Process memory now: 240.31 MiB
  Calculator: 56.86 MiB
    Density: 15.11 MiB
      Arrays: 15.00 MiB
      Localized functions: 0.10 MiB
      Mixer: 0.00 MiB
    Hamiltonian: 11.16 MiB
      Arrays: 11.15 MiB
      XC: 0.00 MiB
      Poisson: 0.00 MiB
      vbar: 0.01 MiB
    Wavefunctions: 30.60 MiB
      Arrays psit_nG: 26.46 MiB
      Eigensolver: 4.13 MiB
      Projections: 0.00 MiB
      Projectors: 0.01 MiB

Total number of cores used: 48
Domain decomposition: 3 x 4 x 4

Number of atoms: 6
Number of atomic orbitals: 80
Number of bands in calculation: 16
Number of valence electrons: 12
Bands to converge: occupied

... initialized

Initializing position-dependent things.

Density initialized from atomic densities
Creating initial wave functions:
  16 bands from LCAO basis set

              .------------------------------------------------.  
             /|                                                |  
            / |                                                |  
           /  |                                                |  
          /   |                                                |  
         /    |                                                |  
        /     |                                                |  
       /      |                                                |  
      /       |                                                |  
     /        |                                                |  
    /         |                                                |  
   /          |                                                |  
  /           |                                                |  
 *            |                                                |  
 |            |                                                |  
 |            |                                                |  
 |            |                                                |  
 |            |                                                |  
 |            |                  H                             |  
 |            |                                                |  
 |            |                HC                              |  
 |            |                 CH                             |  
 |            |                                                |  
 |            |                H                               |  
 |            |                                                |  
 |            |                                                |  
 |            |                                                |  
 |            |                                                |  
 |            .------------------------------------------------.  
 |           /                                                /   
 |          /                                                /    
 |         /                                                /     
 |        /                                                /      
 |       /                                                /       
 |      /                                                /        
 |     /                                                /         
 |    /                                                /          
 |   /                                                /           
 |  /                                                /            
 | /                                                /             
 |/                                                /              
 *------------------------------------------------*               

Positions:
   0 C     10.000000   10.923542   11.892732    ( 0.0000,  0.0000,  0.0000)
   1 C     10.000000   10.923542   10.563732    ( 0.0000,  0.0000,  0.0000)
   2 H     10.000000   11.847084   12.456465    ( 0.0000,  0.0000,  0.0000)
   3 H     10.000000   10.000000   12.456465    ( 0.0000,  0.0000,  0.0000)
   4 H     10.000000   11.847084   10.000000    ( 0.0000,  0.0000,  0.0000)
   5 H     10.000000   10.000000   10.000000    ( 0.0000,  0.0000,  0.0000)

Unit cell:
           periodic     x           y           z      points  spacing
  1. axis:    no    20.000000    0.000000    0.000000   132     0.1515
  2. axis:    no     0.000000   21.847084    0.000000   144     0.1517
  3. axis:    no     0.000000    0.000000   22.456465   148     0.1517

  Lengths:  20.000000  21.847084  22.456465
  Angles:   90.000000  90.000000  90.000000

Effective grid spacing dv^(1/3) = 0.1517

Initial localization: ...
Pipek-Mezey localization started
