  __  _  _
 | _ |_)|_||  |
 |__||  | ||/\| - 26.7.0

User:   root@login01
Date:   Wed Oct 28 00:00:36 2026
Arch:   x86_64
Pid:    836198
CWD:    /public/.ys-full-function-20261027/wave24-gpaw/water-fd
units:  Angstrom and eV
cores:  1
OpenMP: True

Versions:
  Python: 3.11.16
  ASE:    3.29.0
  Numpy:  2.3.5
  Scipy:  1.16.3
  Libxc:  7.0.0

Paths:
  gpaw:  /public/toolkit/miniconda3/envs/gpaw-26.7.0-8173m/lib/python3.11/site-packages/gpaw
  _gpaw: /public/toolkit/miniconda3/envs/gpaw-26.7.0-8173m/lib/python3.11/site-packages/
        _gpaw.cpython-311-x86_64-linux-gnu.so
  ase:   /public/toolkit/miniconda3/envs/gpaw-26.7.0-8173m/lib/python3.11/site-packages/ase
  numpy: /public/toolkit/miniconda3/envs/gpaw-26.7.0-8173m/lib/python3.11/site-packages/numpy
  scipy: /public/toolkit/miniconda3/envs/gpaw-26.7.0-8173m/lib/python3.11/site-packages/scipy

Git-hashes:

OMP_NUM_THREADS: 1


Environment variables:
  GPAW_CPUPY:          False
  GPAW_DEBUG:          False
  GPAW_MPI_BACKEND:    serial
  GPAW_MPI_OPTIONS:    None
  GPAW_NO_C_EXTENSION: False
  GPAW_NO_GPU_MPI:     False
  GPAW_SETUP_PATH:     /public/software/gpaw/datasets/26.8.0
  GPAW_TRACE:          False
  GPAW_USE_GPUS:       False

Input parameters:
  mode=FD(),
  convergence={'energy': 1e-06, 'density': 1e-05, 'eigenstates': 1e-08},
  h=0.25,
  kpts=MonkhorstPack(size=array([1, 1, 1])),
  maxiter=200,
  nbands=12,
  occupations=Occupations(name='fermi-dirac', width=0.01, fixmagmom=False),
  xc=XC(name='PBE')

Diagonalizing LCAO Hamiltonian
Converting LCAO to fd mode
      .--------------.  
     /|              |  
    / |              |  
   /  |              |  
  /   |              |  
 *    |              |  
 |    |    OH        |  
 |    |   H          |  
 |    .--------------.  
 |   /              /   
 |  /              /    
 | /              /     
 |/              /      
 *--------------*       

Atoms  # Å, Bohr magnetons
|   | symbol |        x |        y |        z |    initial magnetic moments |
|---|--------|----------|----------|----------|-----------------------------|
| 0 |      O | 3.000000 | 3.763239 | 3.596309 | ( 0.0000,  0.0000,  0.0000) |
| 1 |      H | 3.000000 | 4.526478 | 3.000000 | ( 0.0000,  0.0000,  0.0000) |
| 2 |      H | 3.000000 | 3.000000 | 3.000000 | ( 0.0000,  0.0000,  0.0000) |


Unit cell  # Å
| periodic |        x |        y |        z |  lengths |    angles |
|----------|----------|----------|----------|----------|-----------|
|       no | 6.000000 | 0.000000 | 0.000000 | 6.000000 | 90.000000 |
|       no | 0.000000 | 7.526478 | 0.000000 | 7.526478 | 90.000000 |
|       no | 0.000000 | 0.000000 | 6.596309 | 6.596309 | 90.000000 |


Number of symmetries: 4
Symmetry operations  # r_c * U_cc + t_c -> r_c
|   | kind |            matrix            |
|---|------|------------------------------|
| 0 |  E   |  1, 0, 0,  0, 1, 0,  0, 0, 1 |
| 1 |  σ   |  1, 0, 0,  0,-1, 0,  0, 0, 1 |
| 2 |  σ   | -1, 0, 0,  0, 1, 0,  0, 0, 1 |
| 3 |  C2  | -1, 0, 0,  0,-1, 0,  0, 0, 1 |

BZ-sampling:
  Number of BZ points: 1
  Number of IBZ points: 1

  Monkhorst-Pack size: [1, 1, 1]
  Monkhorst-Pack shift: [0, 0, 0]

K-points  # in reciprocal-cell coordinates
|   |                                coordinates |     weight |
|---|--------------------------------------------|------------|
| 0 | (  0.00000000,   0.00000000,   0.00000000) | 1.00000000 |

Uniform wave function grid shape: (np.int64(23), np.int64(31), np.int64(27))
Spin-components: 1 (collinear spins)
Bands:           12
Projectors:      11
Spin-degeneracy: 2
Datatype:        float64

Memory:
  Storage: CPU
  Wave functions: 1_848_096 bytes (1_848_096 per core)

Parallelization:
  kpt:    1
  domain: 1
  band:   1

Density:
  Valence electrons: 8.0
  Components: 1
  Grid points: [24 32 28]
  Charge: 0.0  # |e|

Potential:
  Grid points: [24 32 28]

species:
  O:
    name:    Oxygen
    id:      32ecd46bf208036f09c70b0ec9a88b78
    Z:       8.0
    valence: 6
    core:    2
    charge:  0.0  # [|e|]
    file:    /public/software/gpaw/datasets/26.8.0/O.PBE.gz
    compensation charges:
      type: gauss
      rc: 0.21  # [Å]
      lmax: 2
    cutoffs:
      filter: 1.17  # [Å]
      core: 0.83  # [Å]
    projectors
    | kind |    occ | energy [eV] | rcut [Å] |
    |------|--------|-------------|----------|
    |   2s |  2.000 |     -23.961 |    0.688 |
    |   2p |  4.000 |      -9.029 |    0.598 |
    |    s |        |       3.251 |    0.688 |
    |    p |        |      18.182 |    0.598 |
    |    d |        |       0.000 |    0.619 |
  
    Using partial waves for O as LCAO basis

  H:
    name:    Hydrogen
    id:      d65de229564ff8ea4db303e23b6d1ecf
    Z:       1.0
    valence: 1
    core:    0
    charge:  0.0  # [|e|]
    file:    /public/software/gpaw/datasets/26.8.0/H.PBE.gz
    compensation charges:
      type: gauss
      rc: 0.15  # [Å]
      lmax: 2
    cutoffs:
      filter: 0.85  # [Å]
      core: 0.53  # [Å]
    projectors
    | kind |    occ | energy [eV] | rcut [Å] |
    |------|--------|-------------|----------|
    |   1s |  1.000 |      -6.494 |    0.476 |
    |    s |        |      20.717 |    0.476 |
    |    p |        |       0.000 |    0.476 |
  
    Using partial waves for H as LCAO basis

Reference energy: -2065.832626  # eV

eigensolver:
  PPCG
  Converge bands: occupied
  niter=2
  
density mixing:
  MSR1-Mixer
  Pratt-step length: 0.05
  Maximum mixing coefficient: 0.75
  Old densities: 10
  SVD Regularization: 0.005
  Good broyden scale: 1.0
  Trust region scale: 1.0
  Backstep threshold: 4.0
  History deletion threshold: 6.0
  Fourier metric: weight=20.0, sigma=0.015
  
occupation numbers:
  Fermi-Dirac:
    width: 0.0100  # eV
  

FastPoissonSolver using
    Stencil: 6*3+1=19 point O(h^6) finite-difference Laplacian
    FFT axes: []
    FST axes: [np.int64(0), np.int64(2), np.int64(1)]

xc functional:
  name: PBE with 2 nearest neighbor stencil
interpolation: tri-quintic # 5. degree polynomial

Initialization done in: 0.799 s

convergence criteria:
- Maximum [total energy] change in last 3 cyles: 1e-06 eV / valence electron
- Maximum integral of absolute [dens]ity change: 1e-05 electrons / valence electron
- Maximum integral of absolute [eigenst]ate change: 1e-08 eV^2 / valence electron
maximum number of iterations: 200

SCF iterations
|iter     |     time |      energy | eigst |  dens |
|---------|----------|-------------|-------|-------|
|iter:   1| 00:00:37 |  -16.898226 | +1.08 |       |
|iter:   2| 00:00:37 |  -16.438257 | -0.06 | -0.65 |
|iter:   3| 00:00:38 |  -14.872580 | +0.11 | -0.70 |
|iter:   4| 00:00:38 |  -14.858903 | -0.90 | -1.20 |
|iter:   5| 00:00:38 |  -14.854956 | -1.75 | -1.51 |
|iter:   6| 00:00:38 |  -14.855032 | -2.62 | -2.48 |
|iter:   7| 00:00:39 |  -14.855059 | -3.56 | -2.93 |
|iter:   8| 00:00:39 |  -14.855060 | -4.40 | -3.58 |
|iter:   9| 00:00:39 |  -14.855061c| -5.19 | -4.11 |
|iter:  10| 00:00:40 |  -14.855061c| -5.99 | -4.58 |
|iter:  11| 00:00:40 |  -14.855061c| -6.77 | -5.26c|
|iter:  12| 00:00:40 |  -14.855061c| -7.59 | -5.31c|
|iter:  13| 00:00:40 |  -14.855061c| -8.18c| -6.21c|

Converged in 13 steps
SCF loop duration: 3.472 s (0.267 s/step)
Energy contributions relative to reference atoms: (reference = -2065.832626)

Energy contributions:
  Kinetic:          17.826630
  Coulomb:         -20.673742
  Zero:              0.919128
  External:          0.000000
  Xc:              -12.927077
  Entropy:          -0.000000
  Spinorbit:         0.000000
  ----------------------------
  Free energy:     -14.855061
  Extrapolated:    -14.855061

dipole moment: (0.000000, -0.000000, -0.386923)  # |e|*Ang


Fermi level: -5.491
Eigenvalues  # (at kpt = [0.000, 0.000, 0.000])
| band index |      eig [eV] | occ [0-2] |
|------------|---------------|-----------|
|          0 |       -25.110 |     2.000 |
|          1 |       -13.035 |     2.000 |
|          2 |        -9.303 |     2.000 |
|          3 |        -7.237 |     2.000 |
|          4 |        -0.128 |     0.000 |
|          5 |         1.869 |     0.000 |
|          6 |         3.709 |     0.000 |
|          7 |         4.114 |     0.000 |
|          8 |         4.502 |     0.000 |
|          9 |         5.369 |     0.000 |
|         10 |         5.568 |     0.000 |
|         11 |         7.502 |     0.000 |


Gap: 7.109 eV
Transition (v -> c):
  (s=0, k=0, n=3, [0.00, 0.00, 0.00]) -> (s=0, k=0, n=4, [0.00, 0.00, 0.00])
No difference between direct/indirect transitions
vacuum-level: 0.000  # V

Forces  # eV/Ang
|       x |        y |        z |
|---------|----------|----------|
| 0.00000 |  0.00000 | -0.35639 |  # O    0
| 0.00000 |  0.17271 | -0.10953 |  # H    1
| 0.00000 | -0.17271 | -0.10953 |  # H    2


Force computation time: 0.012  # s
# Writing to state.gpw (mode='all')


timing:  # [seconds]
  Total:  4.352  # ━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━
  Init:   0.798  # ━━━━━━━━━━╸
  SCF:    3.472  # ━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━╸
  Forces: 0.012  # 

Max RSS: 167.223  # MiB
