Siesta Version  : 4.1.5
Architecture    : intel-mpi-mkl
Compiler version: ifort (IFORT) 2021.6.0 20220226
Compiler flags  : mpiifort -O2 -fPIC -fp-model source
PP flags        : -DMPI -DFC_HAVE_ABORT
Libraries       :  -L/opt/intel/oneapi/mkl/2022.1.0/lib/intel64 -lmkl_scalapack_lp64 -lmkl_blaccs_intelmpi_lp64 -Wl,--start-group -lmkl_intel_lp64 -lmkl_sequential -lmkl_core -Wl,--end-group -lpthread -lm -ldl
PARALLEL version

* Running on 8 nodes in parallel
>> Start of run:   7-OCT-2026  23:08:35

                           ***********************       
                           *  WELCOME TO SIESTA  *       
                           ***********************       

reinit: Reading from standard input
reinit: Dumped input in INPUT_TMP.74413
************************** Dump of input data file ****************************
SystemName       GNR_6unit       # Descriptive name of the system
SystemLabel      scat            # Short name for naming files
%include STRUCT.fdf
%include BASIS.fdf
%include TS.fdf
%include KPT.fdf
#%include ELEC.fdf
#%include TBTRANS.fdf
%block kgrid_Monkhorst_Pack
 1    0    0    0.0
 0   30    0    0.0
 0    0    1   0.0
%endblock kgrid_Monkhorst_Pack
SolutionMethod      transiesta
XC.functional       GGA
XC.authors          PBE
MeshCutoff          100.0 Ry
SaveHS              T
MaxSCFIterations      500           # Default: 50
DM.MixingWeight       0.10          # Default: 0.25
DM.NumberPulay        3             # Default: 0
DM.PulayOnFile        F             # SystemLabel.P1, SystemLabel.P2
DM.Tolerance          1.d-4         # Default: 1.d-4
DM.UseSaveDM          T             # because of the bug
Diag.ParallelOverK    T
TS.SCF.DM.Tolerance 1.d-4
TS.SCF.dQ.Tolerance 1.d-0
************************** End of input data file *****************************

reinit: -----------------------------------------------------------------------
reinit: System Name: GNR_6unit
reinit: -----------------------------------------------------------------------
reinit: System Label: scat
reinit: -----------------------------------------------------------------------

initatom: Reading input for the pseudopotentials and atomic orbitals ----------
Species number:   1 Atomic number:    6 Label: C
 
Ground state valence configuration:   2s02  2p02
Reading pseudopotential information in formatted form from C.psf

Valence configuration for pseudopotential generation:
2s( 2.00) rc: 1.54
2p( 2.00) rc: 1.54
3d( 0.00) rc: 1.54
4f( 0.00) rc: 1.54
For C, standard SIESTA heuristics set lmxkb to 2
 (one more than the basis l, including polarization orbitals).
Use PS.lmax or PS.KBprojectors blocks to override.

<basis_specs>
===============================================================================
C                    Z=   6    Mass=  12.010        Charge= 0.17977+309
Lmxo=1 Lmxkb= 2    BasisType=split      Semic=F
L=0  Nsemic=0  Cnfigmx=2
          i=1  nzeta=1  polorb=0  (2s)
            splnorm:   0.15000    
               vcte:    0.0000    
               rinn:    0.0000    
               qcoe:    0.0000    
               qyuk:    0.0000    
               qwid:   0.10000E-01
                rcs:    0.0000    
            lambdas:    1.0000    
L=1  Nsemic=0  Cnfigmx=2
          i=1  nzeta=1  polorb=0  (2p)
            splnorm:   0.15000    
               vcte:    0.0000    
               rinn:    0.0000    
               qcoe:    0.0000    
               qyuk:    0.0000    
               qwid:   0.10000E-01
                rcs:    0.0000    
            lambdas:    1.0000    
-------------------------------------------------------------------------------
L=0  Nkbl=1  erefs: 0.17977+309
L=1  Nkbl=1  erefs: 0.17977+309
L=2  Nkbl=1  erefs: 0.17977+309
===============================================================================
</basis_specs>

atom: Called for C                     (Z =   6)

read_vps: Pseudopotential generation method:
read_vps: ATM3      Troullier-Martins                       
Valence charge for ps generation:    4.00000

read_vps: Pseudopotential includes a core correction:
read_vps: Pseudo-core for xc-correction

xc_check: Exchange-correlation functional:
xc_check: GGA Perdew, Burke & Ernzerhof 1996
V l=0 = -2*Zval/r beyond r=  1.5227
V l=1 = -2*Zval/r beyond r=  1.5227
V l=2 = -2*Zval/r beyond r=  1.5227
All V_l potentials equal beyond r=  1.5227
This should be close to max(r_c) in ps generation
All pots = -2*Zval/r beyond r=  1.5227

VLOCAL1: 99.0% of the norm of Vloc inside     18.722 Ry
VLOCAL1: 99.9% of the norm of Vloc inside     42.668 Ry
atom: Maximum radius for 4*pi*r*r*local-pseudopot. charge    1.83678
atom: Maximum radius for r*vlocal+2*Zval:    1.58088
GHOST: No ghost state for L =  0
GHOST: No ghost state for L =  1
GHOST: No ghost state for L =  2

KBgen: Kleinman-Bylander projectors: 
   l= 0   rc=  1.704034   el= -1.009801   Ekb=  5.335044   kbcos=  0.316824
   l= 1   rc=  1.704034   el= -0.388707   Ekb= -3.785862   kbcos= -0.366464
   l= 2   rc=  1.791422   el=  0.001971   Ekb= -0.970471   kbcos= -0.009212

KBgen: Total number of Kleinman-Bylander projectors:    9
atom: -------------------------------------------------------------------------

atom: SANKEY-TYPE ORBITALS:

SPLIT: Orbitals with angular momentum L= 0

SPLIT: Basis orbitals for state 2s

SPLIT: PAO cut-off radius determined from an
SPLIT: energy shift=  0.020000 Ry

   izeta = 1
                 lambda =    1.000000
                     rc =    4.088342
                 energy =   -0.989105
                kinetic =    0.904439
    potential(screened) =   -1.893544
       potential(ionic) =   -5.484581

SPLIT: Orbitals with angular momentum L= 1

SPLIT: Basis orbitals for state 2p

SPLIT: PAO cut-off radius determined from an
SPLIT: energy shift=  0.020000 Ry

   izeta = 1
                 lambda =    1.000000
                     rc =    4.870301
                 energy =   -0.368906
                kinetic =    2.520369
    potential(screened) =   -2.889275
       potential(ionic) =   -6.370528
atom: Total number of Sankey-type orbitals:  4

atm_pop: Valence configuration (for local Pseudopot. screening):
 2s( 2.00)                                                            
 2p( 2.00)                                                            
Vna: chval, zval:    4.00000   4.00000

Vna:  Cut-off radius for the neutral-atom potential:   4.870301
comcore: Pseudo-core radius Rcore=  1.791422

atom: _________________________________________________________________________

prinput: Basis input ----------------------------------------------------------

PAO.BasisType split     

%block ChemicalSpeciesLabel
    1    6 C                       # Species index, atomic number, species label
%endblock ChemicalSpeciesLabel

%block PAO.Basis                 # Define Basis set
C                     2                    # Species label, number of l-shells
 n=2   0   1                         # n, l, Nzeta 
   4.088   
   1.000   
 n=2   1   1                         # n, l, Nzeta 
   4.870   
   1.000   
%endblock PAO.Basis

prinput: ----------------------------------------------------------------------

coor:   Atomic-coordinates input format  =     Cartesian coordinates
coor:                                          (in Angstroms)

siesta: Atomic coordinates (Bohr) and species
siesta:     18.89727   0.00000   0.00000  1        1
siesta:     18.89727   4.63144   0.00000  1        2
siesta:     18.89727   0.00000   2.67396  1        3
siesta:     18.89727   4.63144   2.67396  1        4
siesta:     18.89727   2.31572   4.01095  1        5
siesta:     18.89727   6.94716   4.01095  1        6
siesta:     18.89727   2.31572   6.68491  1        7
siesta:     18.89727   6.94716   6.68491  1        8
siesta:     18.89727   0.00000   8.02189  1        9
siesta:     18.89727   4.63144   8.02189  1       10
siesta:     18.89727   0.00000  10.69585  1       11
siesta:     18.89727   4.63144  10.69585  1       12
siesta:     18.89727   2.31572  12.03284  1       13
siesta:     18.89727   6.94716  12.03284  1       14
siesta:     18.89727   2.31572  14.70680  1       15
siesta:     18.89727   6.94716  14.70680  1       16
siesta:     18.89727   0.00000  16.04378  1       17
siesta:     18.89727   4.63144  16.04378  1       18
siesta:     18.89727   0.00000  18.71774  1       19
siesta:     18.89727   4.63144  18.71774  1       20
siesta:     18.89727   2.31572  20.05473  1       21
siesta:     18.89727   6.94716  20.05473  1       22
siesta:     18.89727   2.31572  22.72869  1       23
siesta:     18.89727   6.94716  22.72869  1       24
siesta:     18.89727   0.00000  24.06567  1       25
siesta:     18.89727   4.63144  24.06567  1       26
siesta:     18.89727   0.00000  26.73964  1       27
siesta:     18.89727   4.63144  26.73964  1       28
siesta:     18.89727   2.31572  28.07662  1       29
siesta:     18.89727   6.94716  28.07662  1       30
siesta:     18.89727   2.31572  30.75058  1       31
siesta:     18.89727   6.94716  30.75058  1       32
siesta:     18.89727   0.00000  32.08756  1       33
siesta:     18.89727   4.63144  32.08756  1       34
siesta:     18.89727   0.00000  34.76153  1       35
siesta:     18.89727   4.63144  34.76153  1       36
siesta:     18.89727   2.31572  36.09851  1       37
siesta:     18.89727   6.94716  36.09851  1       38
siesta:     18.89727   2.31572  38.77247  1       39
siesta:     18.89727   6.94716  38.77247  1       40
siesta:     18.89727   0.00000  40.10945  1       41
siesta:     18.89727   4.63144  40.10945  1       42
siesta:     18.89727   0.00000  42.78342  1       43
siesta:     18.89727   4.63144  42.78342  1       44
siesta:     18.89727   2.31572  44.12040  1       45
siesta:     18.89727   6.94716  44.12040  1       46
siesta:     18.89727   2.31572  46.79436  1       47
siesta:     18.89727   6.94716  46.79436  1       48

siesta: System type = slab      

initatomlists: Number of atoms, orbitals, and projectors:     48   192   432

siesta: ******************** Simulation parameters ****************************
siesta:
siesta: The following are some of the parameters of the simulation.
siesta: A complete list of the parameters used, including default values,
siesta: can be found in file out.fdf
siesta:
redata: Spin configuration                          = none
redata: Number of spin components                   = 1
redata: Time-Reversal Symmetry                      = T
redata: Spin spiral                                 = F
redata: Long output                                 =   F
redata: Number of Atomic Species                    =        1
redata: Charge density info will appear in .RHO file
redata: Write Mulliken Pop.                         = NO
redata: Matel table size (NRTAB)                    =     1024
redata: Mesh Cutoff                                 =   100.0000 Ry
redata: Net charge of the system                    =     0.0000 |e|
redata: Min. number of SCF Iter                     =        0
redata: Max. number of SCF Iter                     =      500
redata: SCF convergence failure will abort job
redata: SCF mix quantity                            = Hamiltonian
redata: Mix DM or H after convergence               =   F
redata: Recompute H after scf cycle                 =   F
redata: Mix DM in first SCF step                    =   T
redata: Write Pulay info on disk                    =   F
redata: New DM Occupancy tolerance                  = 0.000000000001
redata: No kicks to SCF
redata: DM Mixing Weight for Kicks                  =     0.5000
redata: Require Harris convergence for SCF          =   F
redata: Harris energy tolerance for SCF             =     0.000100 eV
redata: Require DM convergence for SCF              =   T
redata: DM tolerance for SCF                        =     0.000100
redata: Require EDM convergence for SCF             =   F
redata: EDM tolerance for SCF                       =     0.001000 eV
redata: Require H convergence for SCF               =   T
redata: Hamiltonian tolerance for SCF               =     0.001000 eV
redata: Require (free) Energy convergence for SCF   =   F
redata: (free) Energy tolerance for SCF             =     0.000100 eV
redata: Using Saved Data (generic)                  =   F
redata: Use continuation files for DM               =   T
redata: Neglect nonoverlap interactions             =   F
redata: Method of Calculation                       = Transiesta
redata: Fix the spin of the system                  =   F
redata: Dynamics option                             = Single-point calculation
mix.SCF: Pulay mixing                            = Pulay
mix.SCF:    Variant                              = stable
mix.SCF:    History steps                        = 3
mix.SCF:    Linear mixing weight                 =     0.100000
mix.SCF:    Mixing weight                        =     0.100000
mix.SCF:    SVD condition                        = 0.1000E-07
redata: ***********************************************************************

%block SCF.Mixers
  Pulay
%endblock SCF.Mixers

%block SCF.Mixer.Pulay
  # Mixing method
  method pulay
  variant stable

  # Mixing options
  weight 0.1000
  weight.linear 0.1000
  history 3
%endblock SCF.Mixer.Pulay
 
DM_history_depth set to one: no extrapolation allowed by default for geometry relaxation
Size of DM history Fstack: 1
Total number of electrons:   192.000000
Total ionic charge:   192.000000

* ProcessorY, Blocksize:    2  25


* Orbital distribution balance (max,min):    25    17

k-point displ. along   2 input, could be:     0.00    0.50
 Kpoints in:           16 . Kpoints trimmed:           16

siesta: k-grid: Number of k-points =    16
siesta: k-grid: Cutoff (effective) =    10.000 Ang
siesta: k-grid: Supercell and displacements
siesta: k-grid:    1   0   0      0.000
siesta: k-grid:    0  30   0      0.000
siesta: k-grid:    0   0   1      0.000
 
diag: Algorithm                                     = D&C
diag: Parallel over k                               =   T
diag: Use parallel 2D distribution                  =   T
diag: Parallel block-size                           = 25
diag: Parallel distribution                         =     2 x     4
diag: Used triangular part                          = Lower
diag: Absolute tolerance                            =  0.100E-15
diag: Orthogonalization factor                      =  0.100E-05
diag: Memory factor                                 =  1.0000

superc: Internal auxiliary supercell:     1 x     3 x     3  =       9
superc: Number of atoms, orbitals, and projectors:    432   1728   3888
 
k-point displ. along   2 input, could be:     0.00    0.50
 Kpoints in:           16 . Kpoints trimmed:           16

transiesta: k-grid: Number of Green function k-points =    16
transiesta: k-grid: Supercell and displacements
transiesta: k-grid:    1   0   0      0.000
transiesta: k-grid:    0  30   0      0.000
transiesta: k-grid:    0   0   1      0.000
 
ts: **************************************************************
ts: Voltage                                         =    0.0000 Volts
ts: Save H and S matrices                           =    T
ts: Electronic temperature                          =  299.986869 K
ts: Transport along unit-cell vector                =    A3
ts: Transport along Cartesian vector                =    Z
ts: Fixing Hartree potential at cell boundary       =    -C
ts: Fix Hartree potential fraction                  =    1.0000
ts: Hartree potential offset                        =    0.000000 eV
ts: Solution method                                 =    BTD
ts: BTD pivoting method                             =    atom+Left
ts: BTD creation algorithm                          =    speed
ts: SCF DM tolerance                                = 0.100E-03
ts: SCF Hamiltonian tolerance                       =    0.001000 eV
ts: SCF converge charge                             =    T
ts: SCF charge tolerance                            = 0.100E+01
ts: Initialize DM by                                =    diagon
ts: *** TranSiesta will NOT update forces ***
ts: Will not correct charge fluctuations
ts: TS.SCF mixing options same as SCF
ts:           >> Electrodes << 
ts: >> Left
ts:   Electrode cell pivoting: E1, E2, E3           = A1, A2, A3
ts:   GF file                                       = scat.TSGFLeft
ts:   Reuse existing GF-file                        =    T
ts:   Electrode TSHS file                           = elec.TSHS
ts:   # atoms used in electrode                     =    8
ts:   Electrode Bloch unity [E1 x E2 x E3]          = 1 x 1 x 1
ts:   Position in geometry                          = 1 -- 8
ts:   Semi-infinite direction for electrode         = negative wrt. E3
ts:   Chemical shift                                =    0.000000 eV
ts:   Electronic temperature                        =  299.986869 K
ts:   Gamma-only electrode                          =    F
ts:   Bulk H, S in electrode region                 =    T
ts:   Initial DM for electrodes                     = diagon
ts:   Cross-terms are updated
ts:   Manual delta-Ef shift                         =    0.000000 eV
ts:   Electrode self-energy imaginary Eta           =  0.1000E-02 eV
ts:   Electrode self-energy accuracy                =  0.1000E-12 eV
ts:   Electrode inter-layer distance (semi-inf)     =    0.7075 Ang
ts: >> Right
ts:   Electrode cell pivoting: E1, E2, E3           = A1, A2, A3
ts:   GF file                                       = scat.TSGFRight
ts:   Reuse existing GF-file                        =    T
ts:   Electrode TSHS file                           = elec.TSHS
ts:   # atoms used in electrode                     =    8
ts:   Electrode Bloch unity [E1 x E2 x E3]          = 1 x 1 x 1
ts:   Position in geometry                          = 41 -- 48
ts:   Semi-infinite direction for electrode         = positive wrt. E3
ts:   Chemical shift                                =    0.000000 eV
ts:   Electronic temperature                        =  299.986869 K
ts:   Gamma-only electrode                          =    F
ts:   Bulk H, S in electrode region                 =    T
ts:   Initial DM for electrodes                     = diagon
ts:   Cross-terms are updated
ts:   Manual delta-Ef shift                         =    0.000000 eV
ts:   Electrode self-energy imaginary Eta           =  0.1000E-02 eV
ts:   Electrode self-energy accuracy                =  0.1000E-12 eV
ts:   Electrode inter-layer distance (semi-inf)     =    0.7075 Ang
ts:  ----------------- Contour ----------------- 
ts:            >> Residual contour << 
ts: Pole chemical potential                         =    0.0000 eV
ts:    Chemical potential temperature               =  299.9869 K
ts:    Number of poles                              =   16
ts:    Top energy point                             =    2.5988 eV
ts:          >> Equilibrium contour << 
ts: Contour name                                    =    TS.Contour.c-Left
ts:   circle contour E_min                          =  -40.0000 eV
ts:   circle contour E_max                          =   -0.2585 eV
ts:   circle contour points                         =   25
ts:   circle contour method                         =    Gauss-Legendre
ts: Contour name                                    =    TS.Contour.t-Left
ts:   tail contour E_min                            =   -0.2585 eV
ts:   tail contour E_max                            =     Infinity
ts:   tail contour points                           =   10
ts:   tail contour method                           =    Gauss-Fermi (0kT)
ts: **************************************************************
 
************************ Begin: TS CHECKS AND WARNINGS ************************
Without loosing performance you can increase the equilibrium integration precision.
You can add 5 more energy points in the equilibrium contours, for FREE!
This is 5 more energy points per chemical potential.
ts: ***       TranSiesta will NOT update forces       ***
ts: *** ALL FORCES AFTER TRANSIESTA HAS RUN ARE WRONG ***
************************ End: TS CHECKS AND WARNINGS **************************


>>> TranSiesta block information for FDF-file START <<<

%block TS.ChemPots
  Left
%endblock TS.ChemPots

%block TS.ChemPot.Left
  mu  V/2
  contour.eq
   begin
    c-Left
    t-Left
   end
%endblock TS.ChemPot.Left
%block TS.Contour.c-Left
  part circle
     from -40.00000 eV + V/2 to -10. kT + V/2
      points 25
        method g-legendre
%endblock TS.Contour.c-Left
%block TS.Contour.t-Left
  part tail
     from prev to + inf
      points 10
        method g-fermi
%endblock TS.Contour.t-Left

>>> TranSiesta block information for FDF-file END <<<


transiesta: Regions of atoms:
   ## (8): Elec.Left
     [ 1 -- 8 ]
   ## (8): Elec.Right
     [ 41 -- 48 ]
   -- (32): Device
     [ 9 -- 40 ]
 
  <sparsity:elec.TSHS
    nrows_g=32 nrows=32 sparsity=4.6250 nnzs=4736, refcount: 3>
  <sparsity:(TM [--,--, 0] of: elec.TSHS)
    nrows_g=32 nrows=32 sparsity=2.2500 nnzs=2304, refcount: 3>
  <sparsity:(TM [--,--,-1] of: elec.TSHS)
    nrows_g=32 nrows=32 sparsity=1.1250 nnzs=1152, refcount: 3>
 Left principal cell is extending out with 64 elements:
    Atom 1 connects with atom 7
    Orbital 3 connects with orbital 27
    Hamiltonian value: |H(3,315)|@R=-2 =  0.282E-03 eV
    Overlap          :  S(3,315)|@R=-2 = -0.122E-07

Calculating all surface Green functions for: Left
 Fermi level shift in electrode (chemical potential) :        0.00000  eV
 Atoms available    / used atoms   :      8 /      8
 Orbitals available / used orbitals:     32 /     32
 Total self-energy calculations: 816
 Saving surface Green functions in: scat.TSGFLeft
 Estimated file size:     13.250 MB
 Lopez Sancho, Lopez Sancho & Rubio recursive surface self-energy calculation...
 Lopez Sancho, Lopez Sancho & Rubio: Mean/std iterations:     4.7181 /     0.7774
 Lopez Sancho, Lopez Sancho & Rubio: Min/Max iterations :          4 /          9
Done creating 'scat.TSGFLeft'.

  <sparsity:elec.TSHS
    nrows_g=32 nrows=32 sparsity=4.6250 nnzs=4736, refcount: 3>
  <sparsity:(TM [--,--, 0] of: elec.TSHS)
    nrows_g=32 nrows=32 sparsity=2.2500 nnzs=2304, refcount: 3>
  <sparsity:(TM [--,--, 1] of: elec.TSHS)
    nrows_g=32 nrows=32 sparsity=1.1250 nnzs=1152, refcount: 3>
 Right principal cell is extending out with 64 elements:
    Atom 7 connects with atom 1
    Orbital 27 connects with orbital 3
    Hamiltonian value: |H(27,195)|@R=2 =  0.282E-03 eV
    Overlap          :  S(27,195)|@R=2 = -0.122E-07

Calculating all surface Green functions for: Right
 Fermi level shift in electrode (chemical potential) :        0.00000  eV
 Atoms available    / used atoms   :      8 /      8
 Orbitals available / used orbitals:     32 /     32
 Total self-energy calculations: 816
 Saving surface Green functions in: scat.TSGFRight
 Estimated file size:     13.250 MB
 Lopez Sancho, Lopez Sancho & Rubio recursive surface self-energy calculation...
 Lopez Sancho, Lopez Sancho & Rubio: Mean/std iterations:     4.7181 /     0.7774
 Lopez Sancho, Lopez Sancho & Rubio: Min/Max iterations :          4 /          9
Done creating 'scat.TSGFRight'.


                     ====================================
                        Single-point calculation
                     ====================================

superc: Internal auxiliary supercell:     1 x     3 x     3  =       9
superc: Number of atoms, orbitals, and projectors:    432   1728   3888

outcell: Unit cell vectors (Ang):
       20.000000    0.000000    0.000000
        0.000000    4.901704    0.000000
        0.000000    0.000000   25.470000

outcell: Cell vector modules (Ang)   :   20.000000    4.901704   25.470000
outcell: Cell angles (23,13,12) (deg):     90.0000     90.0000     90.0000
outcell: Cell volume (Ang**3)        :   2496.9280
<dSpData1D:S at geom step 0
  <sparsity:sparsity for geom step 0
    nrows_g=192 nrows=25 sparsity=.1004 nnzs=3700, refcount: 7>
  <dData1D:(new from dSpData1D) n=3700, refcount: 1>
refcount: 1>
new_DM -- step:     1
Initializing Density Matrix...
 
Attempting to read DM, EDM from TSDE file... Failed...
Attempting to read DM from file... Failed...
DM filled with atomic data:
<dSpData2D:DM initialized from atoms
  <sparsity:sparsity for geom step 0
    nrows_g=192 nrows=25 sparsity=.1004 nnzs=3700, refcount: 8>
  <dData2D:DM n=3700 m=1, refcount: 1>
refcount: 1>

transiesta: Initialization run using siesta

No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps:      30      52
New grid distribution:   1
           1       1:   64    1:    8    1:   20
           2       1:   64    1:    8   21:   40
           3       1:   64    1:    8   41:   60
           4       1:   64    1:    8   61:   80
           5       1:   64    9:   15    1:   20
           6       1:   64    9:   15   21:   40
           7       1:   64    9:   15   41:   60
           8       1:   64    9:   15   61:   80

InitMesh: MESH =   128 x    30 x   160 =      614400
InitMesh: (bp) =    64 x    15 x    80 =       76800
InitMesh: Mesh cutoff (required, used) =   100.000   103.526 Ry
ExtMesh (bp) on 0 =   100 x    40 x    56 =      224000
New grid distribution:   2
           1      33:   64    1:   15    1:   20
           2       1:   32    1:   15   21:   40
           3      33:   64    1:   15   41:   60
           4       1:   32    1:   15   61:   80
           5       1:   32    1:   15    1:   20
           6      33:   64    1:   15   21:   40
           7       1:   32    1:   15   41:   60
           8      33:   64    1:   15   61:   80
New grid distribution:   3
           1      33:   64    1:   15    1:   20
           2       1:   32    1:   15   21:   40
           3      33:   64    1:   15   41:   60
           4       1:   32    1:   15   61:   80
           5       1:   32    1:   15    1:   20
           6      33:   64    1:   15   21:   40
           7       1:   32    1:   15   41:   60
           8      33:   64    1:   15   61:   80
Setting up quadratic distribution...
ExtMesh (bp) on 0 =    68 x    47 x    56 =      178976
PhiOnMesh: Number of (b)points on node 0 =                 9600
PhiOnMesh: nlist on node 0 =                57055
ts-Vha:  0.23317E-08 eV

stepf: Fermi-Dirac step function
ts-Vha: -0.11812E-01 eV

siesta: Program's energy decomposition (eV):
siesta: Ebs     =     -3012.023557
siesta: Eions   =     12268.674044
siesta: Ena     =      2546.772430
siesta: Ekin    =      4921.651470
siesta: Enl     =      -477.381897
siesta: Eso     =         0.000000
siesta: Edftu   =         0.000000
siesta: DEna    =       168.018395
siesta: DUscf   =         5.889913
siesta: DUext   =         0.000000
siesta: Exc     =     -2621.541644
siesta: eta*DQ  =         0.000000
siesta: Emadel  =         0.000000
siesta: Emeta   =         0.000000
siesta: Emolmec =         0.000000
siesta: Ekinion =         0.000000
siesta: Eharris =     -7660.656105
siesta: Etot    =     -7725.265376
siesta: FreeEng =     -7725.270167

        iscf     Eharris(eV)        E_KS(eV)     FreeEng(eV)     dDmax    Ef(eV) dHmax(eV)
   scf:    1    -7660.656105    -7725.265376    -7725.270167  1.441105 -4.572003  1.845079
timer: Routine,Calls,Time,% = IterSCF        1       0.532  39.58
ts-Vha: -0.11728E-01 eV
   scf:    2    -7725.273454    -7725.269417    -7725.274208  0.000826 -4.722409  1.658468
ts-Vha: -0.11400E-01 eV
   scf:    3    -7725.303183    -7725.286432    -7725.291226  0.007301 -6.068174  0.009014
ts-Vha: -0.12162E-01 eV
   scf:    4    -7725.286533    -7725.286498    -7725.291291  0.000403 -6.067210  0.003403
ts-Vha: -0.12022E-01 eV
   scf:    5    -7725.286500    -7725.286499    -7725.291293  0.000029 -6.067556  0.002367
ts-Vha: -0.11795E-01 eV
   scf:    6    -7725.286501    -7725.286500    -7725.291294  0.000061 -6.068176  0.000528

SCF Convergence by DM+H criterion
max |DM_out - DM_in|         :     0.0000611262
max |H_out - H_in|      (eV) :     0.0005278234
SCF cycle converged after 6 iterations

                     ***************************
                     *  WELCOME TO TRANSIESTA  *
                     ***************************

ts-Vha: -0.11795E-01 eV

transiesta: created H and S sparsity pattern:
  <sparsity:(M of: G ((UC of: T T T sparsity for geom step 0)))
    nrows_g=192 nrows=192 sparsity=.3160 nnzs=11648, refcount: 1>

transiesta: update sparsity pattern same as H and S.

transiesta: Determining an optimal tri-matrix using: atom+Left
transiesta: Established a near-optimal partition for the tri-diagonal matrix.
BTD partitions (7): 
  [ 8, [32] * 5, 24 ]
transiesta: Matrix elements in % of full matrix:  43.40278 %

transiesta: mem of electrodes (static):                     0.03MB
transiesta: mem of global update arrays (static):           0.53MB
transiesta: mem of tri-diagonal matrices:                   0.49MB
transiesta: Total memory usage:                             1.05MB

transiesta: Charge distribution, target =    192.00000
Total charge                  [Q]  :    192.00000
Device                        [D]  :    124.54319
Left                          [E1] :     28.54336
Left                 / device [C1] :      3.45678
Right                         [E2] :     28.54230
Right                / device [C2] :      3.45717
Other                         [O]  :      3.45720
Excess charge                [dQ]  : -1.13687E-13

transiesta: *** The forces are NOT updated ***
transiesta: ***  Will set the energy-density matrix to 0!  ***
ts-q:         D        E1        C1        E2        C2        dQ
ts-q:   124.544    28.543     3.455    28.542     3.456 -.275E-02
ts-Vha: -0.14213E-01 eV

siesta: Program's energy decomposition (eV):
siesta: Ebs     =     -3010.794534
siesta: Eions   =     12268.674044
siesta: Ena     =      2546.772430
siesta: Ekin    =      4927.029594
siesta: Enl     =      -481.602698
siesta: Eso     =         0.000000
siesta: Edftu   =         0.000000
siesta: DEna    =       167.042595
siesta: DUscf   =         5.877187
siesta: DUext   =         0.000000
siesta: Exc     =     -2621.738202
siesta: eta*DQ  =         0.000000
siesta: Emadel  =         0.000000
siesta: Emeta   =         0.000000
siesta: Emolmec =         0.000000
siesta: Ekinion =         0.000000
siesta: Eharris =     -7725.299786
siesta: Etot    =     -7725.293138
siesta: FreeEng =     -7725.293138
ts: E-dN      =         0.000000
ts: Band Stru =     -2076.612141
ts: Kinetic   =      3342.281401
ts: Ion-elec  =    -78717.625893
ts: Eharris   =     -9148.536540
ts: Etot      =     -9148.529892
ts: FreeEng   =     -9148.529892

        iscf     Eharris(eV)        E_KS(eV)     FreeEng(eV)     dDmax    Ef(eV) dHmax(eV)
ts-scf:    1    -9148.536540    -9148.529892    -9148.529892  0.002557 -6.068176  0.016469
timer: Routine,Calls,Time,% = TS             1       0.551  13.72
ts-q:         D        E1        C1        E2        C2        dQ
ts-q:   124.545    28.543     3.455    28.542     3.456 -.105E-02
ts-Vha: -0.12509E-01 eV
ts-scf:    2    -9148.578554    -9148.554240    -9148.554240  0.000201 -6.068176  0.004062
ts-q:         D        E1        C1        E2        C2        dQ
ts-q:   124.545    28.543     3.455    28.542     3.456 -.666E-03
ts-Vha: -0.12291E-01 eV
ts-scf:    3    -9148.565930    -9148.560087    -9148.560087  0.000028 -6.068176  0.001576
ts-q:         D        E1        C1        E2        C2        dQ
ts-q:   124.545    28.543     3.455    28.542     3.456 -.582E-03
ts-Vha: -0.12355E-01 eV
ts-scf:    4    -9148.563923    -9148.562006    -9148.562006  0.000013 -6.068176  0.001193
ts-q:         D        E1        C1        E2        C2        dQ
ts-q:   124.545    28.543     3.455    28.542     3.456 -.521E-03
ts-Vha: -0.12397E-01 eV
ts-scf:    5    -9148.565531    -9148.563769    -9148.563769  0.000013 -6.068176  0.000622

SCF Convergence by DM+H+dQ criterion
max |DM_out - DM_in|         :     0.0000129419
max |H_out - H_in|      (eV) :     0.0006216528
     Q(TS) - Q(dev)          :    -0.0005206015
SCF cycle converged after 5 iterations

Using DM_out to compute the final energy and forces
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps:      30      52
ts-Vha: -0.12397E-01 eV

siesta: E_KS(eV) =            -7725.3270

siesta: E_KS - E_eggbox =     -7725.3270

siesta: Atomic forces (eV/Ang):
----------------------------------------
   Tot   -0.000000    0.000000   -0.117792
----------------------------------------
   Max    0.066150
   Res    0.021535    sqrt( Sum f_i^2 / 3N )
----------------------------------------
   Max    0.066150    constrained

Stress tensor Voigt[x,y,z,yz,xz,xy] (kbar):        0.08      106.20      106.52        0.00       -0.00        0.00
(Free)E + p*V (eV/cell)    -7835.8729
Target enthalpy (eV/cell)    -7725.3270

transiesta: Charge distribution, target =    192.00000
Total charge                  [Q]  :    191.99948
Device                        [D]  :    124.54534
Left                          [E1] :     28.54336
Left                 / device [C1] :      3.45545
Right                         [E2] :     28.54230
Right                / device [C2] :      3.45584
Other                         [O]  :      3.45720
Excess charge                [dQ]  : -5.20601E-04


siesta: Program's energy decomposition (eV):
siesta: Ebs     =     -3010.329434
siesta: Eions   =     12268.674044
siesta: Ena     =      2546.772430
siesta: Ekin    =      4927.089357
siesta: Enl     =      -481.608433
siesta: Eso     =         0.000000
siesta: Edftu   =         0.000000
siesta: DEna    =       166.992867
siesta: DUscf   =         5.873238
siesta: DUext   =         0.000000
siesta: Exc     =     -2621.772428
siesta: eta*DQ  =         0.000000
siesta: Emadel  =         0.000000
siesta: Emeta   =         0.000000
siesta: Emolmec =         0.000000
siesta: Ekinion =         0.000000
siesta: Eharris =     -7725.328776
siesta: Etot    =     -7725.327014
siesta: FreeEng =     -7725.327014
ts: E-dN      =         0.000000
ts: Band Stru =     -2076.540231
ts: Kinetic   =      3342.341164
ts: Ion-elec  =    -78718.194500
ts: Eharris   =     -9148.565531
ts: Etot      =     -9148.563769
ts: FreeEng   =     -9148.563769

siesta: Final energy (eV):
siesta:  Band Struct. =   -3010.329434
siesta:       Kinetic =    4927.089357
siesta:       Hartree =   36941.311130
siesta:       Edftu   =       0.000000
siesta:       Eso     =       0.000000
siesta:    Ext. field =       0.000000
siesta:   Exch.-corr. =   -2621.772428
siesta:  Ion-electron =  -78879.705938
siesta:       Ion-ion =   31907.750866
siesta:       Ekinion =       0.000000
siesta:         Total =   -7725.327014
ts:         E-dN  =       0.000000
ts:     Band Stru =   -2076.540231
ts:       Kinetic =    3342.341164
ts:      Ion-elec =  -78718.194500
ts:         Total =   -9148.563769
siesta:         Fermi =      -6.068176

siesta: Atomic forces (eV/Ang):
siesta:      1   -0.000000   -0.000000   -0.040293
siesta:      2   -0.000000    0.000000   -0.040293
siesta:      3   -0.000000   -0.000000    0.023803
siesta:      4   -0.000000    0.000000    0.023803
siesta:      5   -0.000000    0.000000   -0.017661
siesta:      6   -0.000000   -0.000000   -0.017661
siesta:      7    0.000000    0.000000    0.057571
siesta:      8    0.000000   -0.000000    0.057571
siesta:      9   -0.000000   -0.000000   -0.058715
siesta:     10   -0.000000    0.000000   -0.058715
siesta:     11   -0.000000    0.000000   -0.047746
siesta:     12   -0.000000   -0.000000   -0.047746
siesta:     13   -0.000000   -0.000000   -0.007938
siesta:     14   -0.000000    0.000000   -0.007938
siesta:     15   -0.000000    0.000000   -0.004497
siesta:     16   -0.000000   -0.000000   -0.004497
siesta:     17   -0.000000    0.000000    0.066150
siesta:     18   -0.000000   -0.000000    0.066150
siesta:     19   -0.000000    0.000000    0.030820
siesta:     20   -0.000000   -0.000000    0.030820
siesta:     21   -0.000000   -0.000000   -0.017693
siesta:     22   -0.000000    0.000000   -0.017693
siesta:     23   -0.000000    0.000000   -0.024436
siesta:     24   -0.000000   -0.000000   -0.024436
siesta:     25   -0.000000    0.000000   -0.000912
siesta:     26   -0.000000    0.000000   -0.000912
siesta:     27   -0.000000    0.000000   -0.013781
siesta:     28   -0.000000    0.000000   -0.013781
siesta:     29   -0.000000   -0.000000    0.021328
siesta:     30   -0.000000    0.000000    0.021328
siesta:     31   -0.000000    0.000000    0.022921
siesta:     32   -0.000000   -0.000000    0.022921
siesta:     33   -0.000000    0.000000   -0.053182
siesta:     34   -0.000000    0.000000   -0.053182
siesta:     35   -0.000000   -0.000000   -0.048242
siesta:     36   -0.000000    0.000000   -0.048242
siesta:     37   -0.000000    0.000000   -0.008330
siesta:     38   -0.000000   -0.000000   -0.008330
siesta:     39   -0.000000    0.000000    0.001933
siesta:     40   -0.000000   -0.000000    0.001933
siesta:     41    0.000000    0.000000    0.038329
siesta:     42    0.000000   -0.000000    0.038329
siesta:     43   -0.000000   -0.000000    0.064114
siesta:     44   -0.000000    0.000000    0.064114
siesta:     45   -0.000000    0.000000   -0.057610
siesta:     46   -0.000000   -0.000000   -0.057610
siesta:     47   -0.000000    0.000000    0.015169
siesta:     48   -0.000000   -0.000000    0.015169
siesta: ----------------------------------------
siesta:    Tot   -0.000000    0.000000   -0.117792

siesta: Stress tensor (static) (eV/Ang**3):
siesta:     0.000051    0.000000   -0.000000
siesta:     0.000000    0.066286   -0.000000
siesta:    -0.000000    0.000000    0.066481

siesta: Cell volume =       2496.928018 Ang**3

siesta: Pressure (static):
siesta:                Solid            Molecule  Units
siesta:          -0.00048219         -0.00048387  Ry/Bohr**3
siesta:          -0.04427274         -0.04442713  eV/Ang**3
siesta:         -70.93351333        -71.18088188  kBar
(Free)E+ p_basis*V_orbitals  =       -7718.489017
(Free)Eharris+ p_basis*V_orbitals  =       -7718.490779

siesta: Electric dipole (a.u.)  =   -0.000000    0.000000    0.000000
siesta: Electric dipole (Debye) =   -0.000000    0.000000    0.000000

cite: Please see "scat.bib" for an exhaustive BiBTeX file.
cite: Please clearly indicate Siesta version in published work: 4.1.5
cite: This calculation has made use of the following articles
cite: which are encouraged to be cited in a published work.
        Primary SIESTA paper
          DOI: www.doi.org/10.1088/0953-8984/14/11/302
        Primary TranSiesta paper
          DOI: www.doi.org/10.1103/PhysRevB.65.165401
        TranSiesta N-electrode
          DOI: www.doi.org/10.1016/j.cpc.2016.09.022
 
>> End of run:   7-OCT-2026  23:08:42
Job completed
