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:  11-OCT-2026  23:30:48

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

reinit: Reading from standard input
reinit: Dumped input in INPUT_TMP.45287
************************** Dump of input data file ****************************
#(1) General system descriptors
SystemName       Right           # Descriptive name of the system
SystemLabel      Right           # Short name for naming files
%include STRUCT.fdf
%include KPT.fdf
%include BASIS.fdf
#(4) DFT, Grid, SCF
XC.functional         GGA            # LDA or GGA (default = LDA)
XC.authors            PBE            # CA (Ceperley-Aldr) = PZ
MeshCutoff            100.000000    Ry      # Default: 50.0 Ry ~ 0.444 Bohr
                                    #         100.0 Ry ~ 0.314 Bohr
MaxSCFIterations      500           # Default: 50
DM.MixingWeight       0.05          # 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          .true.        # because of the bug
SCFMustConverge       .true.
NeglNonOverlapInt     F             # Default: F
#(5) Eigenvalue problem: order-N or diagonalization
SolutionMethod        Diagon
ElectronicTemperature 300.0 K       # Default: 300.0 K
Diag.ParallelOverK    true
TS.HS.Save .true.
TS.DE.Save .true.
************************** End of input data file *****************************

reinit: -----------------------------------------------------------------------
reinit: System Name: Right
reinit: -----------------------------------------------------------------------
reinit: System Label: Right
reinit: -----------------------------------------------------------------------

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

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

<basis_specs>
===============================================================================
Si                   Z=  14    Mass=  28.090        Charge= 0.17977+309
Lmxo=1 Lmxkb= 3    BasisType=split      Semic=F
L=0  Nsemic=0  Cnfigmx=3
          i=1  nzeta=1  polorb=0  (3s)
            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=3
          i=1  nzeta=1  polorb=1  (3p)
            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
L=3  Nkbl=1  erefs: 0.17977+309
===============================================================================
</basis_specs>

atom: Called for Si                    (Z =  14)

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=  2.0104
V l=1 = -2*Zval/r beyond r=  2.0873
V l=2 = -2*Zval/r beyond r=  2.1670
V l=3 = -2*Zval/r beyond r=  2.1670
All V_l potentials equal beyond r=  2.0613
This should be close to max(r_c) in ps generation
All pots = -2*Zval/r beyond r=  2.1670

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

KBgen: Kleinman-Bylander projectors: 
   l= 0   rc=  2.306811   el= -0.791441   Ekb=  5.289775   kbcos=  0.291595
   l= 1   rc=  2.394965   el= -0.300632   Ekb=  1.944256   kbcos=  0.287657
   l= 2   rc=  2.549439   el=  0.001958   Ekb= -1.788485   kbcos= -0.056492
   l= 3   rc=  2.782581   el=  0.003046   Ekb= -0.472620   kbcos= -0.006203

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

atom: SANKEY-TYPE ORBITALS:

SPLIT: Orbitals with angular momentum L= 0

SPLIT: Basis orbitals for state 3s

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

   izeta = 1
                 lambda =    1.000000
                     rc =    5.007352
                 energy =   -0.772735
                kinetic =    0.581259
    potential(screened) =   -1.353994
       potential(ionic) =   -3.769451

SPLIT: Orbitals with angular momentum L= 1

SPLIT: Basis orbitals for state 3p

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

   izeta = 1
                 lambda =    1.000000
                     rc =    6.270866
                 energy =   -0.282590
                kinetic =    0.876421
    potential(screened) =   -1.159012
       potential(ionic) =   -3.390735

POLgen: Perturbative polarization orbital with L=  2

POLgen: Polarization orbital for state 3p

   izeta = 1
                     rc =    6.270866
                 energy =    0.439433
                kinetic =    1.304939
    potential(screened) =   -0.865506
       potential(ionic) =   -2.929821
atom: Total number of Sankey-type orbitals:  9

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

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

atom: _________________________________________________________________________

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

PAO.BasisType split     

%block ChemicalSpeciesLabel
    1   14 Si                      # Species index, atomic number, species label
%endblock ChemicalSpeciesLabel

%block PAO.Basis                 # Define Basis set
Si                    2                    # Species label, number of l-shells
 n=3   0   1                         # n, l, Nzeta 
   5.007   
   1.000   
 n=3   1   1 P   1                   # n, l, Nzeta, Polarization, NzetaPol
   6.271   
   1.000   
%endblock PAO.Basis

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

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

siesta: Atomic coordinates (Bohr) and species
siesta:      1.76951   5.64968  89.49900  1        1
siesta:      9.17761  13.05779  89.49900  1        2
siesta:      9.18405   5.64484  89.50321  1        3
siesta:      1.77594  13.05295  89.50321  1        4
siesta:      1.79262   1.94838  92.14047  1        5
siesta:      9.20072   9.35649  92.14047  1        6
siesta:      9.18351   1.94517  92.14628  1        7
siesta:      1.77541   9.35328  92.14628  1        8
siesta:     12.90106   1.95838  94.77598  1        9
siesta:      5.49296   9.36649  94.77598  1       10
siesta:      5.49006   1.96153  94.78969  1       11
siesta:     12.89816   9.36963  94.78969  1       12
siesta:     12.91111   5.66602  97.41986  1       13
siesta:      5.50300  13.07412  97.41986  1       14
siesta:      5.50344   5.67050  97.42266  1       15
siesta:     12.91154  13.07861  97.42266  1       16
siesta:      1.79895   5.66601 100.03902  1       17
siesta:      9.20706  13.07412 100.03902  1       18
siesta:      9.20706   5.66602 100.03902  1       19
siesta:      1.79895  13.07412 100.03902  1       20
siesta:      1.79895   1.96196 102.65818  1       21
siesta:      9.20706   9.37007 102.65818  1       22
siesta:      9.20706   1.96196 102.65818  1       23
siesta:      1.79895   9.37007 102.65818  1       24
siesta:      5.50300   1.96196 105.27734  1       25
siesta:     12.91111   9.37007 105.27734  1       26
siesta:     12.91111   1.96196 105.27734  1       27
siesta:      5.50300   9.37007 105.27734  1       28
siesta:     12.91111   5.66601 107.89650  1       29
siesta:      5.50300  13.07412 107.89650  1       30
siesta:      5.50300   5.66602 107.89650  1       31
siesta:     12.91111  13.07412 107.89650  1       32

siesta: System type = bulk      

initatomlists: Number of atoms, orbitals, and projectors:     32   288   512

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                       = Diagonalization
redata: Electronic Temperature                      =   299.9978 K
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.050000
mix.SCF:    Mixing weight                        =     0.050000
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.0500
  weight.linear 0.0500
  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:   128.000000
Total ionic charge:   128.000000

* ProcessorY, Blocksize:    2  24


* Orbital distribution balance (max,min):    48    24

k-point displ. along   1 input, could be:     0.00    0.50
k-point displ. along   2 input, could be:     0.00    0.50
k-point displ. along   3 input, could be:     0.00    0.50
 Kpoints in:          144 . Kpoints trimmed:          132

siesta: k-grid: Number of k-points =   132
siesta: k-grid: Cutoff (effective) =    23.520 Ang
siesta: k-grid: Supercell and displacements
siesta: k-grid:    6   0   0      0.000
siesta: k-grid:    0   6   0      0.000
siesta: k-grid:    0   0   6      0.000
 
diag: Algorithm                                     = D&C
diag: Parallel over k                               =   T
diag: Use parallel 2D distribution                  =   T
diag: Parallel block-size                           = 24
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:     3 x     3 x     3  =      27
superc: Number of atoms, orbitals, and projectors:    864   7776  13824
 
 
ts: **************************************************************
ts: Save H and S matrices                           =    T
ts: Save DM and EDM matrices                        =    T
ts: Only save the overlap matrix S                  =    F
ts: **************************************************************
 
************************ Begin: TS CHECKS AND WARNINGS ************************
************************ End: TS CHECKS AND WARNINGS **************************


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

superc: Internal auxiliary supercell:     3 x     3 x     3  =      27
superc: Number of atoms, orbitals, and projectors:    864   7776  13824

outcell: Unit cell vectors (Ang):
        7.840000    0.000000    0.000000
        0.000000    7.840000    0.000000
        0.000000    0.000000   11.088000

outcell: Cell vector modules (Ang)   :    7.840000    7.840000   11.088000
outcell: Cell angles (23,13,12) (deg):     90.0000     90.0000     90.0000
outcell: Cell volume (Ang**3)        :    681.5306
<dSpData1D:S at geom step 0
  <sparsity:sparsity for geom step 0
    nrows_g=288 nrows=48 sparsity=.7378 nnzs=61192, refcount: 7>
  <dData1D:(new from dSpData1D) n=61192, refcount: 1>
refcount: 1>
new_DM -- step:     1
Initializing Density Matrix...
 
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=288 nrows=48 sparsity=.7378 nnzs=61192, refcount: 8>
  <dData2D:DM n=61192 m=1, refcount: 1>
refcount: 1>
No. of atoms with KB's overlaping orbs in proc 0. Max # of overlaps:     106     249
New grid distribution:   1
           1       1:   24    1:   12    1:    9
           2       1:   24    1:   12   10:   18
           3       1:   24    1:   12   19:   27
           4       1:   24    1:   12   28:   36
           5       1:   24   13:   24    1:    9
           6       1:   24   13:   24   10:   18
           7       1:   24   13:   24   19:   27
           8       1:   24   13:   24   28:   36

InitMesh: MESH =    48 x    48 x    72 =      165888
InitMesh: (bp) =    24 x    24 x    36 =       20736
InitMesh: Mesh cutoff (required, used) =   100.000   103.598 Ry
ExtMesh (bp) on 0 =    68 x    56 x    53 =      201824
New grid distribution:   2
           1       1:   12   13:   24   19:   36
           2      13:   24    1:   12    1:   18
           3      13:   24    1:   12   19:   36
           4      13:   24   13:   24    1:   18
           5       1:   12    1:   12   19:   36
           6       1:   12    1:   12    1:   18
           7      13:   24   13:   24   19:   36
           8       1:   12   13:   24    1:   18
New grid distribution:   3
           1       1:   12   13:   24   19:   36
           2      13:   24    1:   12    1:   18
           3      13:   24    1:   12   19:   36
           4      13:   24   13:   24    1:   18
           5       1:   12    1:   12   19:   36
           6       1:   12    1:   12    1:   18
           7      13:   24   13:   24   19:   36
           8       1:   12   13:   24    1:   18
Setting up quadratic distribution...
ExtMesh (bp) on 0 =    56 x    56 x    62 =      194432
PhiOnMesh: Number of (b)points on node 0 =                 2592
PhiOnMesh: nlist on node 0 =               176250

stepf: Fermi-Dirac step function

siesta: Program's energy decomposition (eV):
siesta: Ebs     =     -1181.418898
siesta: Eions   =      6041.767191
siesta: Ena     =      1837.320519
siesta: Ekin    =      1217.166550
siesta: Enl     =       512.258306
siesta: Eso     =         0.000000
siesta: Edftu   =         0.000000
siesta: DEna    =       126.850565
siesta: DUscf   =         4.982550
siesta: DUext   =         0.000000
siesta: Exc     =     -1766.104768
siesta: eta*DQ  =         0.000000
siesta: Emadel  =         0.000000
siesta: Emeta   =         0.000000
siesta: Emolmec =         0.000000
siesta: Ekinion =         0.000000
siesta: Eharris =     -4121.263759
siesta: Etot    =     -4109.293470
siesta: FreeEng =     -4109.293470

        iscf     Eharris(eV)        E_KS(eV)     FreeEng(eV)     dDmax    Ef(eV) dHmax(eV)
   scf:    1    -4121.263759    -4109.293470    -4109.293470  1.080176 -3.739755  0.296819
timer: Routine,Calls,Time,% = IterSCF        1       5.127  69.08
   scf:    2    -4109.309908    -4109.301743    -4109.301743  0.000384 -3.747277  0.264980
   scf:    3    -4109.385342    -4109.345606    -4109.345606  0.002429 -3.805400  0.173385
   scf:    4    -4109.379853    -4109.364091    -4109.364091  0.002015 -3.868654  0.014450
   scf:    5    -4109.364227    -4109.364163    -4109.364163  0.000150 -3.868791  0.011608
   scf:    6    -4109.364281    -4109.364227    -4109.364227  0.000192 -3.867887  0.011841
   scf:    7    -4109.364268    -4109.364257    -4109.364257  0.000228 -3.866204  0.008842
   scf:    8    -4109.364276    -4109.364266    -4109.364267  0.000032 -3.866244  0.007572
   scf:    9    -4109.364369    -4109.364319    -4109.364319  0.000071 -3.867486  0.004472
   scf:   10    -4109.364342    -4109.364331    -4109.364331  0.000025 -3.867977  0.003190
   scf:   11    -4109.364348    -4109.364339    -4109.364339  0.000046 -3.868611  0.002020
   scf:   12    -4109.364341    -4109.364340    -4109.364340  0.000015 -3.868512  0.001733
   scf:   13    -4109.364341    -4109.364341    -4109.364341  0.000042 -3.868059  0.001153
   scf:   14    -4109.364341    -4109.364341    -4109.364341  0.000018 -3.868194  0.000882

SCF Convergence by DM+H criterion
max |DM_out - DM_in|         :     0.0000175018
max |H_out - H_in|      (eV) :     0.0008821509
SCF cycle converged after 14 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:     106     249

siesta: E_KS(eV) =            -4109.3643

siesta: E_KS - E_eggbox =     -4109.3643

siesta: Atomic forces (eV/Ang):
----------------------------------------
   Tot    0.136493   -0.251937    0.035109
----------------------------------------
   Max    0.366944
   Res    0.105689    sqrt( Sum f_i^2 / 3N )
----------------------------------------
   Max    0.366944    constrained

Stress tensor Voigt[x,y,z,yz,xz,xy] (kbar):      -17.43      -11.30      -14.49       -0.43       -1.02       -0.06
(Free)E + p*V (eV/cell)    -4103.2345
Target enthalpy (eV/cell)    -4109.3643

siesta: Program's energy decomposition (eV):
siesta: Ebs     =     -1178.115728
siesta: Eions   =      6041.767191
siesta: Ena     =      1837.320519
siesta: Ekin    =      1221.764267
siesta: Enl     =       514.319729
siesta: Eso     =         0.000000
siesta: Edftu   =         0.000000
siesta: DEna    =       121.914214
siesta: DUscf   =         4.908786
siesta: DUext   =         0.000000
siesta: Exc     =     -1767.824664
siesta: eta*DQ  =         0.000000
siesta: Emadel  =         0.000000
siesta: Emeta   =         0.000000
siesta: Emolmec =         0.000000
siesta: Ekinion =         0.000000
siesta: Eharris =     -4109.364341
siesta: Etot    =     -4109.364341
siesta: FreeEng =     -4109.364341

siesta: Final energy (eV):
siesta:  Band Struct. =   -1178.115728
siesta:       Kinetic =    1221.764267
siesta:       Hartree =     210.835867
siesta:       Edftu   =       0.000000
siesta:       Eso     =       0.000000
siesta:    Ext. field =       0.000000
siesta:   Exch.-corr. =   -1767.824664
siesta:  Ion-electron =   -1447.289831
siesta:       Ion-ion =   -2326.849980
siesta:       Ekinion =       0.000000
siesta:         Total =   -4109.364341
siesta:         Fermi =      -3.868194

siesta: Atomic forces (eV/Ang):
siesta:      1    0.215021    0.015361    0.366944
siesta:      2    0.212131    0.012308    0.362783
siesta:      3    0.165071    0.023395    0.346911
siesta:      4    0.167482    0.021621    0.347363
siesta:      5   -0.015343   -0.020383    0.025569
siesta:      6   -0.017178   -0.022400    0.025664
siesta:      7    0.057292   -0.012451    0.021984
siesta:      8    0.059279   -0.014643    0.018682
siesta:      9   -0.005018    0.013275   -0.008526
siesta:     10   -0.003206    0.011323   -0.008506
siesta:     11   -0.007776    0.002790   -0.013213
siesta:     12   -0.010524    0.000772   -0.010578
siesta:     13   -0.012491   -0.036228   -0.090715
siesta:     14   -0.009524   -0.038468   -0.086703
siesta:     15   -0.012612   -0.013831   -0.095287
siesta:     16   -0.015064   -0.016958   -0.096452
siesta:     17    0.003887   -0.017484    0.000555
siesta:     18    0.001066   -0.020572   -0.002512
siesta:     19    0.011744   -0.019602    0.000673
siesta:     20    0.014173   -0.022710    0.001763
siesta:     21    0.003259   -0.000964    0.000173
siesta:     22    0.001457   -0.003091    0.000596
siesta:     23    0.000208   -0.000948    0.002763
siesta:     24    0.002091   -0.003068    0.000203
siesta:     25    0.002889   -0.010792   -0.020901
siesta:     26   -0.000158   -0.012997   -0.019132
siesta:     27    0.004158   -0.013049   -0.023362
siesta:     28    0.006949   -0.015382   -0.022925
siesta:     29   -0.182482   -0.015643   -0.263272
siesta:     30   -0.179282   -0.017909   -0.259075
siesta:     31   -0.158949   -0.000793   -0.232752
siesta:     32   -0.162056   -0.002414   -0.233603
siesta: ----------------------------------------
siesta:    Tot    0.136493   -0.251937    0.035109

siesta: Stress tensor (static) (eV/Ang**3):
siesta:    -0.010881   -0.000036   -0.000635
siesta:    -0.000036   -0.007056   -0.000269
siesta:    -0.000635   -0.000269   -0.009046

siesta: Cell volume =        681.530573 Ang**3

siesta: Pressure (static):
siesta:                Solid            Molecule  Units
siesta:           0.00009796          0.00017161  Ry/Bohr**3
siesta:           0.00899418          0.01575633  eV/Ang**3
siesta:          14.41042187         25.24469699  kBar
(Free)E+ p_basis*V_orbitals  =       -4094.022949
(Free)Eharris+ p_basis*V_orbitals  =       -4094.022950

cite: Please see "Right.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
 
>> End of run:  11-OCT-2026  23:31:59
Job completed
