Siesta Version: siesta-3.2-pl-5
Architecture : i686-apple-darwin15.0.0--Gfortran
Compiler flags: gfortran -g -O2
SERIAL version
* Running in serial mode
>> Start of run: 30-JAN-2017 5:53:15
***********************
* WELCOME TO SIESTA *
***********************
reinit: Reading from standard input
************************** Dump of input data file ****************************
# Created by GDIS version 0.90.0
#
SystemLabel Zn1Al_fcc_FCC
NumberOfAtoms 4
NumberOfSpecies 2
%block ChemicalSpeciesLabel
1 13 Al
2 30 Zn
%endblock ChemicalSpeciesLabel
LatticeConstant 3.99 Ang
%block LatticeParameters
0.99 1.02 0.99 90.000000 90.000000 90.000000
%endblock LatticeParameters
AtomicCoordinatesFormat Fractional
%block AtomicCoordinatesAndAtomicSpecies
0 0 0 1
0.5 0.5 0 2
0.5 0 0.5 2
0 0.5 0.5 2
%endblock AtomicCoordinatesAndAtomicSpecies
%block kgrid_Monkhorst_Pack # k points
12 0 0 0.000
0 12 0 0.000
0 0 12 0.000
%endblock kgrid_Monkhorst_Pack
MeshCutoff 250.0 Ry
%block PAO.Basis # Define Basis set
Al 2
n=3 0 2
5.81289 3.04130
1.00000 1.00000
n=3 1 2 P 1
6.44483 5.01313
1.00000 1.00000
Zn 2 # Species label, number of l-shells
n=4 0 2 P 1 # n, l, Nzeta, Polarization, NzetaPol
6.847 5.072
1.000 1.000
n=3 2 2 # n, l, Nzeta
3.665 1.547
1.000 1.000
%endblock PAO.Basis
MaxSCFIterations 300
DM.Tolerance 5.d-4
DM.NumberPulay 10
DM.MixingWeight 0.01
Diag.DivideAndConquer .false.
SolutionMethod diagon
XC.authors PBE # Echange-correlation
XC.functional GGA
###################################################
# Simulation Type
MD.TypeOfRun Broyden
MD.NumCGsteps 0
MD.MaxCGDispl 0.1 Ang
MD.MaxForceTol 0.04 eV/Ang
MD.MaxStressTol 0.001000 GPa
MD.VariableCell .true.
MD.RelaxCellOnly .true.
SlabDipoleCorrection .true.
###########################
WriteForces .true.
************************** End of input data file *****************************
reinit: -----------------------------------------------------------------------
reinit: System Name:
reinit: -----------------------------------------------------------------------
reinit: System Label: Zn1Al_fcc_FCC
reinit: -----------------------------------------------------------------------
initatom: Reading input for the pseudopotentials and atomic orbitals ----------
Species number: 1 Label: Al Atomic number: 13
Species number: 2 Label: Zn Atomic number: 30
Ground state valence configuration: 3s02 3p01
Reading pseudopotential information in formatted form from Al.psf
Pseudopotential generated from a relativistic atomic calculation
There are spin-orbit pseudopotentials available
Spin-orbit interaction is not included in this calculation
Valence configuration for pseudopotential generation:
3s( 2.00) rc: 2.28
3p( 1.00) rc: 2.28
3d( 0.00) rc: 2.28
4f( 0.00) rc: 2.28
Ground state valence configuration: 4s02 3d10
Reading pseudopotential information in formatted form from Zn.psf
Valence configuration for pseudopotential generation:
4s( 2.00) rc: 2.14
4p( 0.00) rc: 2.14
3d(10.00) rc: 2.14
4f( 0.00) rc: 2.14
For Al, 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.
For Zn, 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.
Warning: Empty PAO shell. l = 1
Will have a KB projector anyway...
===============================================================================
Al Z= 13 Mass= 26.980 Charge= 0.17977+309
Lmxo=1 Lmxkb= 3 BasisType=split Semic=F
L=0 Nsemic=0 Cnfigmx=3
n=1 nzeta=2 polorb=0
splnorm: 0.15000
vcte: 0.0000
rinn: 0.0000
rcs: 5.8129 3.0413
lambdas: 1.0000 1.0000
L=1 Nsemic=0 Cnfigmx=3
n=1 nzeta=2 polorb=1
splnorm: 0.15000
vcte: 0.0000
rinn: 0.0000
rcs: 6.4448 5.0131
lambdas: 1.0000 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
===============================================================================
atom: Called for Al (Z = 13)
read_vps: Pseudopotential generation method:
read_vps: ATM 3.2.2 Troullier-Martins
Total valence charge: 3.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.4534
V l=1 = -2*Zval/r beyond r= 2.4229
V l=2 = -2*Zval/r beyond r= 2.4229
V l=3 = -2*Zval/r beyond r= 2.4229
All V_l potentials equal beyond r= 2.2478
This should be close to max(r_c) in ps generation
All pots = -2*Zval/r beyond r= 2.4534
VLOCAL1: 99.0% of the norm of Vloc inside 8.591 Ry
VLOCAL1: 99.9% of the norm of Vloc inside 19.579 Ry
atom: Maximum radius for 4*pi*r*r*local-pseudopot. charge 2.67776
atom: Maximum radius for r*vlocal+2*Zval: 2.30474
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.677757 el= -0.569464 Ekb= 3.981824 kbcos= 0.292046
l= 1 rc= 2.745549 el= -0.199425 Ekb= 1.608990 kbcos= 0.255417
l= 2 rc= 2.886326 el= 0.001961 Ekb= -1.154048 kbcos= -0.055097
l= 3 rc= 3.111140 el= 0.003052 Ekb= -0.339449 kbcos= -0.007888
KBgen: Total number of Kleinman-Bylander projectors: 16
atom: -------------------------------------------------------------------------
atom: SANKEY-TYPE ORBITALS:
atom: Selected multiple-zeta basis: split
SPLIT: Orbitals with angular momentum L= 0
SPLIT: Basis orbitals for state 3s
izeta = 1
lambda = 1.000000
rc = 5.812541
energy = -0.554173
kinetic = 0.409002
potential(screened) = -0.963175
potential(ionic) = -2.405428
izeta = 2
rmatch = 3.034321
splitnorm = 0.891294
energy = 3.817364
kinetic = 4.676562
potential(screened) = -0.859199
potential(ionic) = -2.231100
SPLIT: Orbitals with angular momentum L= 1
SPLIT: Basis orbitals for state 3p
izeta = 1
lambda = 1.000000
rc = 6.423871
energy = -0.162864
kinetic = 0.675764
potential(screened) = -0.838628
potential(ionic) = -2.200402
izeta = 2
rmatch = 5.002874
splitnorm = 0.266355
energy = -0.011336
kinetic = 1.103145
potential(screened) = -1.114481
potential(ionic) = -2.589707
POLgen: Perturbative polarization orbital with L= 2
POLgen: Polarization orbital for state 3p
izeta = 1
rc = 6.423871
energy = 0.396031
kinetic = 1.061177
potential(screened) = -0.665146
potential(ionic) = -1.955564
atom: Total number of Sankey-type orbitals: 13
atm_pop: Valence configuration (for local Pseudopot. screening):
3s( 2.00)
3p( 1.00)
Vna: chval, zval: 3.00000 3.00000
Vna: Cut-off radius for the neutral-atom potential: 6.423871
comcore: Pseudo-core radius Rcore= 3.395648
atom: _________________________________________________________________________
===============================================================================
Zn Z= 30 Mass= 65.390 Charge= 0.17977+309
Lmxo=2 Lmxkb= 3 BasisType=split Semic=F
L=0 Nsemic=0 Cnfigmx=4
n=1 nzeta=2 polorb=1
splnorm: 0.15000
vcte: 0.0000
rinn: 0.0000
rcs: 6.8470 5.0720
lambdas: 1.0000 1.0000
L=1 Nsemic=0 Cnfigmx=4
L=2 Nsemic=0 Cnfigmx=3
n=1 nzeta=2 polorb=0
splnorm: 0.15000
vcte: 0.0000
rinn: 0.0000
rcs: 3.6650 1.5470
lambdas: 1.0000 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
===============================================================================
atom: Called for Zn (Z = 30)
read_vps: Pseudopotential generation method:
read_vps: ATM3 Troullier-Martins
Total valence charge: 12.00000
xc_check: Exchange-correlation functional:
xc_check: GGA Perdew, Burke & Ernzerhof 1996
V l=0 = -2*Zval/r beyond r= 3.0764
V l=1 = -2*Zval/r beyond r= 3.0764
V l=2 = -2*Zval/r beyond r= 3.0764
V l=3 = -2*Zval/r beyond r= 3.0764
All V_l potentials equal beyond r= 2.1144
This should be close to max(r_c) in ps generation
All pots = -2*Zval/r beyond r= 3.0764
Using large-core scheme for Vlocal
atom: Estimated core radius 3.07645
atom: Including non-local core corrections could be a good idea
atom: Maximum radius for 4*pi*r*r*local-pseudopot. charge 3.35777
atom: Maximum radius for r*vlocal+2*Zval: 3.11515
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.195182 el= -0.428170 Ekb= 3.641528 kbcos= 0.209834
l= 1 rc= 2.195182 el= -0.085221 Ekb= 1.185033 kbcos= 0.193783
l= 2 rc= 2.167912 el= -0.776629 Ekb=-16.402020 kbcos= -0.773612
l= 3 rc= 2.195182 el= 0.003130 Ekb= -1.611899 kbcos= -0.007520
KBgen: Total number of Kleinman-Bylander projectors: 16
atom: -------------------------------------------------------------------------
atom: SANKEY-TYPE ORBITALS:
atom: Selected multiple-zeta basis: split
SPLIT: Orbitals with angular momentum L= 0
SPLIT: Basis orbitals for state 4s
izeta = 1
lambda = 1.000000
rc = 6.846865
energy = -0.421513
kinetic = 0.364327
potential(screened) = -0.785840
potential(ionic) = -10.198130
izeta = 2
rmatch = 5.072261
splitnorm = 0.250691
energy = -0.300715
kinetic = 0.698757
potential(screened) = -0.999473
potential(ionic) = -11.519293
SPLIT: Orbitals with angular momentum L= 2
SPLIT: Basis orbitals for state 3d
izeta = 1
lambda = 1.000000
rc = 3.664824
energy = -0.769709
kinetic = 11.344692
potential(screened) = -12.114401
potential(ionic) = -29.325261
izeta = 2
rmatch = 1.566353
splitnorm = 0.239282
energy = 0.782278
kinetic = 17.889091
potential(screened) = -17.106813
potential(ionic) = -35.977884
POLgen: Perturbative polarization orbital with L= 1
POLgen: Polarization orbital for state 4s
izeta = 1
rc = 6.846865
energy = -0.028511
kinetic = 0.688555
potential(screened) = -0.717066
potential(ionic) = -9.212158
atom: Total number of Sankey-type orbitals: 15
atm_pop: Valence configuration (for local Pseudopot. screening):
4s( 2.00)
4p( 0.00)
3d(10.00)
Vna: chval, zval: 12.00000 12.00000
Vna: Cut-off radius for the neutral-atom potential: 6.846865
atom: _________________________________________________________________________
prinput: Basis input ----------------------------------------------------------
PAO.BasisType split
%block ChemicalSpeciesLabel
1 13 Al # Species index, atomic number, species label
2 30 Zn # Species index, atomic number, species label
%endblock ChemicalSpeciesLabel
%block PAO.Basis # Define Basis set
Al 2 # Species label, number of l-shells
n=3 0 2 # n, l, Nzeta
5.813 3.034
1.000 1.000
n=3 1 2 P 1 # n, l, Nzeta, Polarization, NzetaPol
6.424 5.003
1.000 1.000
Zn 2 # Species label, number of l-shells
n=4 0 2 P 1 # n, l, Nzeta, Polarization, NzetaPol
6.847 5.072
1.000 1.000
n=3 2 2 # n, l, Nzeta
3.665 1.566
1.000 1.000
%endblock PAO.Basis
prinput: ----------------------------------------------------------------------
coor: Atomic-coordinates input format = Fractional
siesta: Atomic coordinates (Bohr) and species
siesta: 0.00000 0.00000 0.00000 1 1
siesta: 3.73231 3.84541 0.00000 2 2
siesta: 3.73231 0.00000 3.73231 2 3
siesta: 0.00000 3.84541 3.73231 2 4
siesta: System type = bulk
initatomlists: Number of atoms, orbitals, and projectors: 4 58 64
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: Non-Collinear-spin run = F
redata: SpinPolarized (Up/Down) run = F
redata: Number of spin components = 1
redata: Long output = F
redata: Number of Atomic Species = 2
redata: Charge density info will appear in .RHO file
redata: Write Mulliken Pop. = NO
redata: Mesh Cutoff = 250.0000 Ry
redata: Net charge of the system = 0.0000 |e|
redata: Max. number of SCF Iter = 300
redata: Performing Pulay mixing using = 10 iterations
redata: Mix DM in first SCF step ? = F
redata: Write Pulay info on disk? = F
redata: Discard 1st Pulay DM after kick = F
redata: New DM Mixing Weight = 0.0100
redata: New DM Occupancy tolerance = 0.000000000001
redata: No kicks to SCF
redata: DM Mixing Weight for Kicks = 0.5000
redata: DM Tolerance for SCF = 0.000500
redata: Require Energy convergence for SCF = F
redata: DM Energy tolerance for SCF = 0.000100 eV
redata: Require Harris convergence for SCF = F
redata: DM Harris energy tolerance for SCF = 0.000100 eV
redata: Using Saved Data (generic) = F
redata: Use continuation files for DM = F
redata: Neglect nonoverlap interactions = F
redata: Method of Calculation = Diagonalization
redata: Divide and Conquer = F
redata: Electronic Temperature = 0.0019 Ry
redata: Fix the spin of the system = F
redata: Dynamics option = Broyden coord. optimization
redata: Variable cell = T
redata: Maximum number of CG moves = 0
redata: Force tolerance = 0.0016 Ry/Bohr
redata: Stress tolerance = 0.0010 GPa
redata: ***********************************************************************
Total number of electrons: 39.000000
Total ionic charge: 39.000000
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: 1008 . Kpoints trimmed: 948
siesta: k-grid: Number of k-points = 948
siesta: k-grid: Cutoff (effective) = 23.701 Ang
siesta: k-grid: Supercell and displacements
siesta: k-grid: 12 0 0 0.000
siesta: k-grid: 0 12 0 0.000
siesta: k-grid: 0 0 12 0.000
Naive supercell factors: 6 6 6
superc: Internal auxiliary supercell: 6 x 6 x 6 = 216
superc: Number of atoms, orbitals, and projectors: 864 12528 13824
* Maximum dynamic memory allocated = 2 MB
siesta: ==============================
Begin CG move = 0
==============================
superc: Internal auxiliary supercell: 6 x 6 x 6 = 216
superc: Number of atoms, orbitals, and projectors: 864 12528 13824
outcell: Unit cell vectors (Ang):
3.950100 0.000000 0.000000
0.000000 4.069800 0.000000
0.000000 0.000000 3.950100
outcell: Cell vector modules (Ang) : 3.950100 4.069800 3.950100
outcell: Cell angles (23,13,12) (deg): 90.0000 90.0000 90.0000
outcell: Cell volume (Ang**3) : 63.5023
New_DM. Step: 1
Initializing Density Matrix...
InitMesh: MESH = 40 x 40 x 40 = 64000
InitMesh: Mesh cutoff (required, used) = 250.000 266.978 Ry
efield: SlabDipoleCorrection = .true. in input file
efield: A dipole layer will be introduced in the vacuum
efield: region to compensate the system dipole
Dipole moment in unit cell = 0.0000 0.0000 0.0000 D
Electric field for dipole correction = -0.000000 -0.000000 -0.000000 Ry/Bohr/e
* Maximum dynamic memory allocated = 62 MB
stepf: Fermi-Dirac step function
siesta: Program's energy decomposition (eV):
siesta: Ebs = -323.576504
siesta: Eions = 4532.180669
siesta: Ena = 77.880204
siesta: Ekin = 4680.770977
siesta: Enl = -4037.603429
siesta: DEna = -0.000019
siesta: DUscf = 0.000000
siesta: DUext = -0.000000
siesta: Exc = -748.736176
siesta: eta*DQ = 0.000000
siesta: Emadel = 0.000000
siesta: Emeta = 0.000000
siesta: Emolmec = 0.000000
siesta: Ekinion = 0.000000
siesta: Eharris = -4553.513465
siesta: Etot = -4559.869113
siesta: FreeEng = -4559.869113
siesta: iscf Eharris(eV) E_KS(eV) FreeEng(eV) dDmax Ef(eV)
siesta: 1 -4553.5135 -4559.8691 -4559.8691 1.3785 -2.7447
timer: Routine,Calls,Time,% = IterSCF 1 10.112 94.60
elaps: Routine,Calls,Wall,% = IterSCF 1 10.144 94.59
Dipole moment in unit cell = 0.0000 0.0000 0.0000 D
Electric field for dipole correction = -0.000000 -0.000000 -0.000000 Ry/Bohr/e
siesta: 2 -4553.4385 -4550.9574 -4550.9601 0.3115 -3.1547
Dipole moment in unit cell = 0.0000 0.0000 0.0000 D
Electric field for dipole correction = -0.000000 -0.000000 -0.000000 Ry/Bohr/e
siesta: 3 -4552.4516 -4554.1402 -4554.1430 0.6117 -3.0548
Dipole moment in unit cell = 0.0000 0.0000 0.0000 D
Electric field for dipole correction = -0.000000 -0.000000 -0.000000 Ry/Bohr/e
siesta: 4 -4552.5461 -4553.9643 -4553.9668 0.6079 -3.0634
Dipole moment in unit cell = 0.0000 0.0000 0.0000 D
Electric field for dipole correction = -0.000000 -0.000000 -0.000000 Ry/Bohr/e
siesta: 5 -4552.2428 -4552.1861 -4552.1888 0.1308 -3.0825
Dipole moment in unit cell = 0.0000 0.0000 0.0000 D
Electric field for dipole correction = -0.000000 -0.000000 -0.000000 Ry/Bohr/e
siesta: 6 -4552.2427 -4552.1858 -4552.1885 0.1271 -3.0818
Dipole moment in unit cell = 0.0000 0.0000 0.0000 D
Electric field for dipole correction = -0.000000 -0.000000 -0.000000 Ry/Bohr/e
siesta: 7 -4552.2280 -4552.1414 -4552.1442 0.0085 -3.0588
Dipole moment in unit cell = 0.0000 0.0000 0.0000 D
Electric field for dipole correction = -0.000000 -0.000000 -0.000000 Ry/Bohr/e
siesta: 8 -4552.2280 -4552.1363 -4552.1389 0.0112 -3.0597
Dipole moment in unit cell = 0.0000 0.0000 0.0000 D
Electric field for dipole correction = -0.000000 -0.000000 -0.000000 Ry/Bohr/e
siesta: 9 -4552.2279 -4552.1898 -4552.1924 0.0034 -3.0590
Dipole moment in unit cell = 0.0000 0.0000 0.0000 D
Electric field for dipole correction = -0.000000 -0.000000 -0.000000 Ry/Bohr/e
siesta: 10 -4552.2279 -4552.1939 -4552.1965 0.0036 -3.0591
Dipole moment in unit cell = 0.0000 0.0000 0.0000 D
Electric field for dipole correction = -0.000000 -0.000000 -0.000000 Ry/Bohr/e
siesta: 11 -4552.2278 -4552.2246 -4552.2272 0.0002 -3.0586
Dipole moment in unit cell = 0.0000 0.0000 0.0000 D
Electric field for dipole correction = -0.000000 -0.000000 -0.000000 Ry/Bohr/e
siesta: E_KS(eV) = -4552.2248
siesta: E_KS - E_eggbox = -4552.2248
siesta: Atomic forces (eV/Ang):
1 -0.000001 -0.000001 -0.000001
2 -0.000002 -0.000002 -0.000005
3 -0.000000 -0.000002 0.000001
4 -0.000004 -0.000002 -0.000001
----------------------------------------
Tot -0.000006 -0.000007 -0.000006
----------------------------------------
Max 0.000005
Res 0.000002 sqrt( Sum f_i^2 / 3N )
----------------------------------------
Max 0.000005 constrained
Stress-tensor-Voigt (kbar): -10.62 17.28 -10.62 -0.00 0.00 0.00
(Free)E + p*V (eV/cell) -4552.1750
Target enthalpy (eV/cell) -4552.2274
siesta: Stress tensor (static) (eV/Ang**3):
-0.006631 -0.000000 0.000000
-0.000000 0.010785 0.000000
0.000000 0.000000 -0.006631
siesta: Pressure (static): 1.32248946 kBar
siesta: Stress tensor (total) (eV/Ang**3):
-0.006631 -0.000000 0.000000
-0.000000 0.010785 0.000000
0.000000 0.000000 -0.006631
siesta: Pressure (total): 1.32248946 kBar
* Maximum dynamic memory allocated = 63 MB
outcell: Unit cell vectors (Ang):
3.950100 0.000000 0.000000
0.000000 4.069800 0.000000
0.000000 0.000000 3.950100
outcell: Cell vector modules (Ang) : 3.950100 4.069800 3.950100
outcell: Cell angles (23,13,12) (deg): 90.0000 90.0000 90.0000
outcell: Cell volume (Ang**3) : 63.5023
siesta: Program's energy decomposition (eV):
siesta: Ebs = -388.916289
siesta: Eions = 4532.180669
siesta: Ena = 77.880204
siesta: Ekin = 4611.476154
siesta: Enl = -3974.257248
siesta: DEna = 0.661444
siesta: DUscf = 0.974099
siesta: DUext = -0.000000
siesta: Exc = -736.778752
siesta: eta*DQ = 0.000000
siesta: Emadel = 0.000000
siesta: Emeta = 0.000000
siesta: Emolmec = 0.000000
siesta: Ekinion = 0.000000
siesta: Eharris = -4552.227793
siesta: Etot = -4552.224769
siesta: FreeEng = -4552.227395
siesta: Final energy (eV):
siesta: Band Struct. = -388.916289
siesta: Kinetic = 4611.476154
siesta: Hartree = 1481.936175
siesta: Ext. field = -0.000000
siesta: Exch.-corr. = -736.778752
siesta: Ion-electron = -7069.314934
siesta: Ion-ion = -2839.543412
siesta: Ekinion = 0.000000
siesta: Total = -4552.224769
siesta: Stress tensor (static) (eV/Ang**3):
siesta: -0.006631 -0.000000 0.000000
siesta: -0.000000 0.010785 0.000000
siesta: 0.000000 0.000000 -0.006631
siesta: Cell volume = 63.502270 Ang**3
siesta: Pressure (static):
siesta: Solid Molecule Units
siesta: 0.00000899 0.00000899 Ry/Bohr**3
siesta: 0.00082542 0.00082539 eV/Ang**3
siesta: 1.32248946 1.32243991 kBar
(Free)E+ p_basis*V_orbitals = -4549.626858
(Free)Eharris+ p_basis*V_orbitals = -4549.630035
* Maximum dynamic memory allocated : Node 0 = 63 MB
* Maximum memory occured during ranger
timer: CPU execution times:
timer: Routine Calls Time/call Tot.time %
timer: siesta 1 77.384 77.384 100.00
timer: Setup 1 0.142 0.142 0.18
timer: bands 1 0.000 0.000 0.00
timer: writewave 1 0.000 0.000 0.00
timer: KSV_init 1 0.000 0.000 0.00
timer: IterMD 1 77.122 77.122 99.66
timer: hsparse 1 0.023 0.023 0.03
timer: overlap 1 0.393 0.393 0.51
timer: IterSCF 11 6.778 74.558 96.35
timer: kinefsm 2 0.198 0.396 0.51
timer: nlefsm 2 0.744 1.489 1.92
timer: DHSCF 12 1.163 13.954 18.03
timer: DHSCF1 1 0.058 0.058 0.08
timer: DHSCF2 1 2.307 2.307 2.98
timer: REORD 98 0.000 0.006 0.01
timer: POISON 13 0.004 0.046 0.06
timer: DHSCF3 12 0.896 10.752 13.89
timer: rhoofd 12 0.515 6.176 7.98
timer: cellXC 12 0.056 0.675 0.87
timer: vmat 12 0.320 3.841 4.96
timer: MolMec 2 0.000 0.000 0.00
timer: diagon 11 5.524 60.761 78.52
timer: c-eigval 10428 0.002 20.758 26.83
timer: c-buildHS 10428 0.001 10.845 14.01
timer: cdiag 20856 0.001 18.239 23.57
timer: cdiag1 20856 0.000 1.039 1.34
timer: cdiag2 20856 0.000 2.729 3.53
timer: cdiag3 20856 0.001 13.924 17.99
timer: cdiag4 20856 0.000 0.342 0.44
timer: c-eigvec 10428 0.002 19.287 24.92
timer: c-buildD 10428 0.002 20.485 26.47
timer: DHSCF4 1 0.835 0.835 1.08
timer: dfscf 1 0.705 0.705 0.91
timer: overfsm 1 0.009 0.009 0.01
timer: optical 1 0.000 0.000 0.00
elaps: ELAPSED times:
elaps: Routine Calls Time/call Tot.time %
elaps: siesta 1 77.795 77.795 100.00
elaps: Setup 1 0.143 0.143 0.18
elaps: bands 1 0.000 0.000 0.00
elaps: writewave 1 0.000 0.000 0.00
elaps: KSV_init 1 0.000 0.000 0.00
elaps: IterMD 1 77.532 77.532 99.66
elaps: hsparse 1 0.022 0.022 0.03
elaps: overlap 1 0.394 0.394 0.51
elaps: IterSCF 11 6.815 74.961 96.36
elaps: kinefsm 2 0.198 0.396 0.51
elaps: nlefsm 2 0.749 1.499 1.93
elaps: DHSCF 12 1.167 14.003 18.00
elaps: DHSCF1 1 0.059 0.059 0.08
elaps: DHSCF2 1 2.317 2.317 2.98
elaps: REORD 98 0.000 0.008 0.01
elaps: POISON 13 0.004 0.049 0.06
elaps: DHSCF3 12 0.899 10.788 13.87
elaps: rhoofd 12 0.516 6.197 7.97
elaps: cellXC 12 0.056 0.674 0.87
elaps: vmat 12 0.321 3.851 4.95
elaps: MolMec 2 0.000 0.000 0.00
elaps: diagon 11 5.550 61.055 78.48
elaps: c-eigval 10428 0.002 20.851 26.80
elaps: c-buildHS 10428 0.001 10.900 14.01
elaps: cdiag 20856 0.001 18.295 23.52
elaps: cdiag1 20856 0.000 1.073 1.38
elaps: cdiag2 20856 0.000 2.763 3.55
elaps: cdiag3 20856 0.001 13.937 17.92
elaps: cdiag4 20856 0.000 0.328 0.42
elaps: c-eigvec 10428 0.002 19.351 24.87
elaps: c-buildD 10428 0.002 20.627 26.51
elaps: DHSCF4 1 0.837 0.837 1.08
elaps: dfscf 1 0.706 0.706 0.91
elaps: overfsm 1 0.009 0.009 0.01
elaps: optical 1 0.000 0.000 0.00
>> End of run: 30-JAN-2017 5:54:33