Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- mahmood@cluster:A$ date
- Thu Oct 6 20:33:48 IRST 2016
- mahmood@cluster:A$ ../sie/siesta < A.fdf
- Siesta Version: siesta-3.2-pl-5
- Architecture : x86_64-unknown-linux-gnu--unknown
- Compiler flags: f95 -g -O2
- SERIAL version
- * Running in serial mode
- >> Start of run: 6-OCT-2016 20:33:50
- ***********************
- * WELCOME TO SIESTA *
- ***********************
- reinit: Reading from standard input
- ************************** Dump of input data file ****************************
- SystemLabel A
- NumberOfAtoms 54
- NumberOfSpecies 2
- %block ChemicalSpeciesLabel
- 1 6 C
- 2 1 H
- %endblock ChemicalSpeciesLabel
- PAO.EnergyShift 100 meV
- PAO.BasisType split
- %block PAO.BasisSizes
- C DZP
- H DZP
- %endblock PAO.BasisSizes
- LatticeConstant 1.000 Ang
- #kgrid_cutoff 20. Ang
- %block kgrid_Monkhorst_Pack
- 1 0 0 0.0
- 0 1 0 0.0
- 0 0 3 0.0
- %endblock kgrid_Monkhorst_Pack
- # SPIN options
- xc.functional GGA # Exchange-correlation functional
- xc.authors PBE # Exchange-correlation version
- SpinPolarized F # Logical parameters are: yes or no
- FixSpin F
- TotalSpin 0.0
- NonCollinearSpin F # 'T', 'F'
- MeshCutoff 200. Ry # Equivalent plane wave cutoff for the grid
- # SCF options
- DM.Tolerance 4.d-4 # Tolerance in maximum difference between input and output DM
- MaxSCFIterations 1000 # Maximum number of SCF iter
- DM.UseSaveDM T # to use continuation files
- DM.MixingWeight 0.05 # New DM amount for next SCF cycle
- DM.NumberPulay 10
- DM.MixSCF1 F
- DM.PulayOnFile F # Store in memory ('F') or in files ('T')
- # NeglNonOverlapInt T # 'F'=do not neglect
- SolutionMethod Diagon # OrderN or Diagon
- ElectronicTemperature 300 K # Temp. for Fermi smearing
- # MD options
- MD.TypeOfRun CG # Type of dynamics:
- MD.VariableCell F
- MD.NumCGsteps 700 # Number of CG steps for coordinate optimization
- ZM.UnitsLength Ang #the units of length used during Z-matrix input
- ZM.UnitsAngle deg #the units of angles used during Z-matrix input
- ZM.ForceTolLength 0.02 eV/Ang #controls the convergence with respect to forces on Z-matrix lengths
- ZM.ForceTolAngle 0.00356549 Ry/rad #controls the convergence with respect to forces on Z-matrix angles
- ZM.MaxDisplLength 0.1 Ang # controls the maximum change in a Zmatrix length during an optimisation step
- ZM.MaxDisplAngle 0.003 rad # controls the maximum change in a Z-matrix angle during an optimisation step
- ZM.CalcAllForces T # Default value
- AtomicCoordinatesFormat NotScaledCartesianAng
- %block LatticeVectors
- 10.0 0.0 0.0
- 0.0 20.0 0.0
- 0.0 0.0 20.0
- %endblock LatticeVectors
- AtomicCoordinatesFormat NotScaledCartesianAng
- %block AtomicCoordinatesAndAtomicSpecies
- 5.000000000 7.539607750 4.316560000 1
- 5.000000000 10.000609750 4.316560000 1
- 5.000000000 12.461612750 4.316560000 1
- 5.000000000 5.355612750 4.476490000 2
- 5.000000000 14.618759750 4.491990000 2
- 5.000000000 6.309106750 5.026990000 1
- 5.000000000 8.770108750 5.026990000 1
- 5.000000000 11.231112750 5.026990000 1
- 5.000000000 13.692112750 5.026990000 1
- 5.000000000 6.309106750 6.447850000 1
- 5.000000000 8.770108750 6.447850000 1
- 5.000000000 11.231112750 6.447850000 1
- 5.000000000 13.692112750 6.447850000 1
- 5.000000000 14.618759750 6.982850000 2
- 5.000000000 5.355612750 6.998350000 2
- 5.000000000 7.539607750 7.158280000 1
- 5.000000000 10.000609750 7.158280000 1
- 5.000000000 12.461612750 7.158280000 1
- 5.000000000 7.539607750 8.579140000 1
- 5.000000000 10.000609750 8.579140000 1
- 5.000000000 12.461612750 8.579140000 1
- 5.000000000 5.355612750 8.739070000 2
- 5.000000000 14.618759750 8.754570000 2
- 5.000000000 6.309106750 9.289570000 1
- 5.000000000 8.770108750 9.289570000 1
- 5.000000000 11.231112750 9.289570000 1
- 5.000000000 13.692112750 9.289570000 1
- 5.000000000 6.309106750 10.710430000 1
- 5.000000000 8.770108750 10.710430000 1
- 5.000000000 11.231112750 10.710430000 1
- 5.000000000 13.692112750 10.710430000 1
- 5.000000000 14.618759750 11.245430000 2
- 5.000000000 5.355612750 11.260930000 2
- 5.000000000 7.539607750 11.420860000 1
- 5.000000000 10.000609750 11.420860000 1
- 5.000000000 12.461612750 11.420860000 1
- 5.000000000 7.539607750 12.841720000 1
- 5.000000000 10.000609750 12.841720000 1
- 5.000000000 12.461612750 12.841720000 1
- 5.000000000 5.355612750 13.001650000 2
- 5.000000000 14.618759750 13.017150000 2
- 5.000000000 6.309106750 13.552150000 1
- 5.000000000 8.770108750 13.552150000 1
- 5.000000000 11.231112750 13.552150000 1
- 5.000000000 13.692112750 13.552150000 1
- 5.000000000 6.309106750 14.973010000 1
- 5.000000000 8.770108750 14.973010000 1
- 5.000000000 11.231112750 14.973010000 1
- 5.000000000 13.692112750 14.973010000 1
- 5.000000000 14.618759750 15.508010000 2
- 5.000000000 5.355612750 15.523510000 2
- 5.000000000 7.539607750 15.683440000 1
- 5.000000000 10.000609750 15.683440000 1
- 5.000000000 12.461612750 15.683440000 1
- %endblock AtomicCoordinatesAndAtomicSpecies
- WriteWaveFunctions T # coefficients of the wavefunctions in the basis set orbitals expansion
- #WaveFuncKPointsScale ReciprocalLatticeVectors
- # Output options
- WriteCoorInitial T
- WriteCoorStep T
- WriteCoorXmol T
- WriteForces T
- WriteEigenvalues F # If .false., it writes them in the file Systemlabel.EIG
- WriteMullikenPop 1 # Write Mulliken Population Analysis
- #WriteKpoints T
- #WriteKbands T
- #WriteBands T
- WriteMDCoorXmol T
- WriteMDhistory T
- WriteCoorCerius t
- # Options for saving or reading information
- MD.UseSaveZM T # Use stored positions and velocities
- MD.UseSaveCG T # Use stored positions and velocities
- SaveRho T # Write valence pseudocharge at the mesh
- SaveDeltaRho T # Write RHOscf-RHOatm at the mesh
- SaveElectrostaticPotential T # Write the total elect. pot. at the mesh
- # (local pseudopotential + Hartree)
- SaveTotalPotential T # write the valence total effective local potential
- # (local pseudopotential + Hartree + Vxc)
- WriteSiestaDim T # Write minimum dim to siesta.h and stop
- WriteDenchar T # Write information for DENCHAR
- %block ProjectedDensityOfStates
- -10.00 0.00 0.10 500 eV
- %endblock ProjectedDensityOfStates
- ************************** End of input data file *****************************
- reinit: -----------------------------------------------------------------------
- reinit: System Name:
- reinit: -----------------------------------------------------------------------
- reinit: System Label: A
- reinit: -----------------------------------------------------------------------
- initatom: Reading input for the pseudopotentials and atomic orbitals ----------
- Species number: 1 Label: C Atomic number: 6
- Species number: 2 Label: H Atomic number: 1
- 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
- Ground state valence configuration: 1s01
- Reading pseudopotential information in formatted form from H.psf
- Valence configuration for pseudopotential generation:
- 1s( 1.00) rc: 1.31
- 2p( 0.00) rc: 1.31
- 3d( 0.00) rc: 0.37
- 4f( 0.00) rc: 1.31
- resizes: Read basis size for species C = dzp
- resizes: Read basis size for species H = dzp
- For C, 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 H, 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= 3 BasisType=split Semic=F
- L=0 Nsemic=0 Cnfigmx=2
- n=1 nzeta=2 polorb=0
- splnorm: 0.15000
- vcte: 0.0000
- rinn: 0.0000
- rcs: 0.0000 0.0000
- lambdas: 1.0000 1.0000
- L=1 Nsemic=0 Cnfigmx=2
- n=1 nzeta=2 polorb=1
- splnorm: 0.15000
- vcte: 0.0000
- rinn: 0.0000
- rcs: 0.0000 0.0000
- 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
- ===============================================================================
- </basis_specs>
- atom: Called for C (Z = 6)
- read_vps: Pseudopotential generation method:
- read_vps: ATM3 Troullier-Martins
- Total valence charge: 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
- V l=3 = -2*Zval/r beyond r= 1.5038
- 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
- GHOST: No ghost state for L = 3
- 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
- l= 3 rc= 1.859892 el= 0.003065 Ekb= -0.399525 kbcos= -0.001223
- 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 2s
- SPLIT: PAO cut-off radius determined from an
- SPLIT: energy shift= 0.007350 Ry
- izeta = 1
- lambda = 1.000000
- rc = 4.632753
- energy = -1.002270
- kinetic = 0.855760
- potential(screened) = -1.858031
- potential(ionic) = -5.417575
- izeta = 2
- rmatch = 3.431921
- splitnorm = 0.150000
- energy = -0.867111
- kinetic = 1.300804
- potential(screened) = -2.167915
- potential(ionic) = -5.962684
- 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.007350 Ry
- izeta = 1
- lambda = 1.000000
- rc = 5.658549
- energy = -0.381140
- kinetic = 2.399701
- potential(screened) = -2.780841
- potential(ionic) = -6.206085
- izeta = 2
- rmatch = 3.653288
- splitnorm = 0.150000
- energy = -0.249149
- kinetic = 3.542762
- potential(screened) = -3.791911
- potential(ionic) = -7.581078
- POLgen: Perturbative polarization orbital with L= 2
- POLgen: Polarization orbital for state 2p
- izeta = 1
- rc = 5.658549
- energy = 1.109913
- kinetic = 2.304789
- potential(screened) = -1.194876
- potential(ionic) = -4.013592
- atom: Total number of Sankey-type orbitals: 13
- 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: 5.658549
- comcore: Pseudo-core radius Rcore= 1.791422
- atom: _________________________________________________________________________
- <basis_specs>
- ===============================================================================
- H Z= 1 Mass= 1.0100 Charge= 0.17977+309
- Lmxo=0 Lmxkb= 2 BasisType=split Semic=F
- L=0 Nsemic=0 Cnfigmx=1
- n=1 nzeta=2 polorb=1
- splnorm: 0.15000
- vcte: 0.0000
- rinn: 0.0000
- rcs: 0.0000 0.0000
- 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
- ===============================================================================
- </basis_specs>
- atom: Called for H (Z = 1)
- read_vps: Pseudopotential generation method:
- read_vps: ATM3 Troullier-Martins
- Total valence charge: 1.00000
- xc_check: Exchange-correlation functional:
- xc_check: GGA Perdew, Burke & Ernzerhof 1996
- V l=0 = -2*Zval/r beyond r= 1.2977
- V l=1 = -2*Zval/r beyond r= 1.2815
- V l=2 = -2*Zval/r beyond r= 1.2656
- All V_l potentials equal beyond r= 1.2977
- This should be close to max(r_c) in ps generation
- All pots = -2*Zval/r beyond r= 1.2977
- VLOCAL1: 99.0% of the norm of Vloc inside 25.776 Ry
- VLOCAL1: 99.9% of the norm of Vloc inside 58.745 Ry
- atom: Maximum radius for 4*pi*r*r*local-pseudopot. charge 1.52711
- atom: Maximum radius for r*vlocal+2*Zval: 1.28155
- 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.434438 el= -0.477200 Ekb= -1.957043 kbcos= -0.361614
- l= 1 rc= 1.508111 el= 0.001076 Ekb= -0.424224 kbcos= -0.025441
- l= 2 rc= 1.527112 el= 0.002010 Ekb= -0.378578 kbcos= -0.002010
- KBgen: Total number of Kleinman-Bylander projectors: 9
- atom: -------------------------------------------------------------------------
- atom: SANKEY-TYPE ORBITALS:
- atom: Selected multiple-zeta basis: split
- SPLIT: Orbitals with angular momentum L= 0
- SPLIT: Basis orbitals for state 1s
- SPLIT: PAO cut-off radius determined from an
- SPLIT: energy shift= 0.007350 Ry
- izeta = 1
- lambda = 1.000000
- rc = 5.471469
- energy = -0.469871
- kinetic = 0.914968
- potential(screened) = -1.384840
- potential(ionic) = -1.907761
- izeta = 2
- rmatch = 3.952598
- splitnorm = 0.150000
- energy = -0.377108
- kinetic = 1.448266
- potential(screened) = -1.825374
- potential(ionic) = -2.372743
- POLgen: Perturbative polarization orbital with L= 1
- POLgen: Polarization orbital for state 1s
- izeta = 1
- rc = 5.471469
- energy = 0.628222
- kinetic = 1.275876
- potential(screened) = -0.647654
- potential(ionic) = -1.125590
- atom: Total number of Sankey-type orbitals: 5
- atm_pop: Valence configuration (for local Pseudopot. screening):
- 1s( 1.00)
- Vna: chval, zval: 1.00000 1.00000
- Vna: Cut-off radius for the neutral-atom potential: 5.471469
- atom: _________________________________________________________________________
- prinput: Basis input ----------------------------------------------------------
- PAO.BasisType split
- %block ChemicalSpeciesLabel
- 1 6 C # Species index, atomic number, species label
- 2 1 H # 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 2 # n, l, Nzeta
- 4.633 3.432
- 1.000 1.000
- n=2 1 2 P 1 # n, l, Nzeta, Polarization, NzetaPol
- 5.659 3.653
- 1.000 1.000
- H 1 # Species label, number of l-shells
- n=1 0 2 P 1 # n, l, Nzeta, Polarization, NzetaPol
- 5.471 3.953
- 1.000 1.000
- %endblock PAO.Basis
- prinput: ----------------------------------------------------------------------
- coor: Atomic-coordinates input format = Cartesian coordinates
- coor: (in Angstroms)
- siesta: Atomic coordinates (Bohr) and species
- siesta: 9.44863 14.24780 8.15712 1 1
- siesta: 9.44863 18.89842 8.15712 1 2
- siesta: 9.44863 23.54904 8.15712 1 3
- siesta: 9.44863 10.12065 8.45934 2 4
- siesta: 9.44863 27.62546 8.48863 2 5
- siesta: 9.44863 11.92249 9.49964 1 6
- siesta: 9.44863 16.57311 9.49964 1 7
- siesta: 9.44863 21.22374 9.49964 1 8
- siesta: 9.44863 25.87435 9.49964 1 9
- siesta: 9.44863 11.92249 12.18468 1 10
- siesta: 9.44863 16.57311 12.18468 1 11
- siesta: 9.44863 21.22374 12.18468 1 12
- siesta: 9.44863 25.87435 12.18468 1 13
- siesta: 9.44863 27.62546 13.19568 2 14
- siesta: 9.44863 10.12065 13.22497 2 15
- siesta: 9.44863 14.24780 13.52719 1 16
- siesta: 9.44863 18.89842 13.52719 1 17
- siesta: 9.44863 23.54904 13.52719 1 18
- siesta: 9.44863 14.24780 16.21223 1 19
- siesta: 9.44863 18.89842 16.21223 1 20
- siesta: 9.44863 23.54904 16.21223 1 21
- siesta: 9.44863 10.12065 16.51446 2 22
- siesta: 9.44863 27.62546 16.54375 2 23
- siesta: 9.44863 11.92249 17.55475 1 24
- siesta: 9.44863 16.57311 17.55475 1 25
- siesta: 9.44863 21.22374 17.55475 1 26
- siesta: 9.44863 25.87435 17.55475 1 27
- siesta: 9.44863 11.92249 20.23979 1 28
- siesta: 9.44863 16.57311 20.23979 1 29
- siesta: 9.44863 21.22374 20.23979 1 30
- siesta: 9.44863 25.87435 20.23979 1 31
- siesta: 9.44863 27.62546 21.25079 2 32
- siesta: 9.44863 10.12065 21.28008 2 33
- siesta: 9.44863 14.24780 21.58231 1 34
- siesta: 9.44863 18.89842 21.58231 1 35
- siesta: 9.44863 23.54904 21.58231 1 36
- siesta: 9.44863 14.24780 24.26734 1 37
- siesta: 9.44863 18.89842 24.26734 1 38
- siesta: 9.44863 23.54904 24.26734 1 39
- siesta: 9.44863 10.12065 24.56957 2 40
- siesta: 9.44863 27.62546 24.59886 2 41
- siesta: 9.44863 11.92249 25.60986 1 42
- siesta: 9.44863 16.57311 25.60986 1 43
- siesta: 9.44863 21.22374 25.60986 1 44
- siesta: 9.44863 25.87435 25.60986 1 45
- siesta: 9.44863 11.92249 28.29490 1 46
- siesta: 9.44863 16.57311 28.29490 1 47
- siesta: 9.44863 21.22374 28.29490 1 48
- siesta: 9.44863 25.87435 28.29490 1 49
- siesta: 9.44863 27.62546 29.30590 2 50
- siesta: 9.44863 10.12065 29.33519 2 51
- siesta: 9.44863 14.24780 29.63742 1 52
- siesta: 9.44863 18.89842 29.63742 1 53
- siesta: 9.44863 23.54904 29.63742 1 54
- siesta: System type = molecule
- initatomlists: Number of atoms, orbitals, and projectors: 54 606 780
- 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. = Atomic and Orbital charges
- redata: Mesh Cutoff = 200.0000 Ry
- redata: Net charge of the system = 0.0000 |e|
- redata: Max. number of SCF Iter = 1000
- 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.0500
- 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.000400
- 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 = T
- redata: Neglect nonoverlap interactions = F
- redata: Method of Calculation = Diagonalization
- redata: Divide and Conquer = T
- redata: Electronic Temperature = 0.0019 Ry
- redata: Fix the spin of the system = F
- redata: Dynamics option = CG coord. optimization
- redata: Variable cell = F
- redata: Use continuation files for CG = T
- redata: Max atomic displ per move = 0.2000 Bohr
- redata: Maximum number of CG moves = 700
- redata: Force tolerance = 0.0016 Ry/Bohr
- redata: ***********************************************************************
- Total number of electrons: 180.000000
- Total ionic charge: 180.000000
- Kpoints in: 2 . Kpoints trimmed: 2
- siesta: k-grid: Number of k-points = 2
- siesta: k-grid: Cutoff (effective) = 5.000 Ang
- siesta: k-grid: Supercell and displacements
- siesta: k-grid: 1 0 0 0.000
- siesta: k-grid: 0 1 0 0.000
- siesta: k-grid: 0 0 3 0.000
- Naive supercell factors: 2 1 1
- superc: Internal auxiliary supercell: 2 x 1 x 1 = 2
- superc: Number of atoms, orbitals, and projectors: 108 1212 1560
- * Maximum dynamic memory allocated = 2 MB
- siesta: ==============================
- Begin CG move = 0
- ==============================
- outcoor: Atomic coordinates (Ang):
- 5.00000000 7.53960775 4.31656000 1 1 C
- 5.00000000 10.00060975 4.31656000 1 2 C
- 5.00000000 12.46161275 4.31656000 1 3 C
- 5.00000000 5.35561275 4.47649000 2 4 H
- 5.00000000 14.61875975 4.49199000 2 5 H
- 5.00000000 6.30910675 5.02699000 1 6 C
- 5.00000000 8.77010875 5.02699000 1 7 C
- 5.00000000 11.23111275 5.02699000 1 8 C
- 5.00000000 13.69211275 5.02699000 1 9 C
- 5.00000000 6.30910675 6.44785000 1 10 C
- 5.00000000 8.77010875 6.44785000 1 11 C
- 5.00000000 11.23111275 6.44785000 1 12 C
- 5.00000000 13.69211275 6.44785000 1 13 C
- 5.00000000 14.61875975 6.98285000 2 14 H
- 5.00000000 5.35561275 6.99835000 2 15 H
- 5.00000000 7.53960775 7.15828000 1 16 C
- 5.00000000 10.00060975 7.15828000 1 17 C
- 5.00000000 12.46161275 7.15828000 1 18 C
- 5.00000000 7.53960775 8.57914000 1 19 C
- 5.00000000 10.00060975 8.57914000 1 20 C
- 5.00000000 12.46161275 8.57914000 1 21 C
- 5.00000000 5.35561275 8.73907000 2 22 H
- 5.00000000 14.61875975 8.75457000 2 23 H
- 5.00000000 6.30910675 9.28957000 1 24 C
- 5.00000000 8.77010875 9.28957000 1 25 C
- 5.00000000 11.23111275 9.28957000 1 26 C
- 5.00000000 13.69211275 9.28957000 1 27 C
- 5.00000000 6.30910675 10.71043000 1 28 C
- 5.00000000 8.77010875 10.71043000 1 29 C
- 5.00000000 11.23111275 10.71043000 1 30 C
- 5.00000000 13.69211275 10.71043000 1 31 C
- 5.00000000 14.61875975 11.24543000 2 32 H
- 5.00000000 5.35561275 11.26093000 2 33 H
- 5.00000000 7.53960775 11.42086000 1 34 C
- 5.00000000 10.00060975 11.42086000 1 35 C
- 5.00000000 12.46161275 11.42086000 1 36 C
- 5.00000000 7.53960775 12.84172000 1 37 C
- 5.00000000 10.00060975 12.84172000 1 38 C
- 5.00000000 12.46161275 12.84172000 1 39 C
- 5.00000000 5.35561275 13.00165000 2 40 H
- 5.00000000 14.61875975 13.01715000 2 41 H
- 5.00000000 6.30910675 13.55215000 1 42 C
- 5.00000000 8.77010875 13.55215000 1 43 C
- 5.00000000 11.23111275 13.55215000 1 44 C
- 5.00000000 13.69211275 13.55215000 1 45 C
- 5.00000000 6.30910675 14.97301000 1 46 C
- 5.00000000 8.77010875 14.97301000 1 47 C
- 5.00000000 11.23111275 14.97301000 1 48 C
- 5.00000000 13.69211275 14.97301000 1 49 C
- 5.00000000 14.61875975 15.50801000 2 50 H
- 5.00000000 5.35561275 15.52351000 2 51 H
- 5.00000000 7.53960775 15.68344000 1 52 C
- 5.00000000 10.00060975 15.68344000 1 53 C
- 5.00000000 12.46161275 15.68344000 1 54 C
- superc: Internal auxiliary supercell: 2 x 1 x 1 = 2
- superc: Number of atoms, orbitals, and projectors: 108 1212 1560
- outcell: Unit cell vectors (Ang):
- 10.000000 0.000000 0.000000
- 0.000000 20.000000 0.000000
- 0.000000 0.000000 20.000000
- outcell: Cell vector modules (Ang) : 10.000000 20.000000 20.000000
- outcell: Cell angles (23,13,12) (deg): 90.0000 90.0000 90.0000
- outcell: Cell volume (Ang**3) : 4000.0000
- New_DM. Step: 1
- Initializing Density Matrix...
- InitMesh: MESH = 90 x 180 x 180 = 2916000
- InitMesh: Mesh cutoff (required, used) = 200.000 223.865 Ry
- * Maximum dynamic memory allocated = 363 MB
- stepf: Fermi-Dirac step function
- siesta: Program's energy decomposition (eV):
- siesta: Ebs = -2764.552497
- siesta: Eions = 10960.021970
- siesta: Ena = 2402.139500
- siesta: Ekin = 3869.988559
- siesta: Enl = 66.413028
- siesta: DEna = -0.000008
- siesta: DUscf = 0.000000
- siesta: DUext = 0.000000
- siesta: Exc = -2376.340999
- siesta: eta*DQ = 0.000000
- siesta: Emadel = 0.000000
- siesta: Emeta = 0.000000
- siesta: Emolmec = 0.000000
- siesta: Ekinion = 0.000000
- siesta: Eharris = -6997.285168
- siesta: Etot = -6997.821891
- siesta: FreeEng = -6997.821891
- siesta: iscf Eharris(eV) E_KS(eV) FreeEng(eV) dDmax Ef(eV)
- siesta: 1 -6997.2852 -6997.8219 -6997.8219 1.9478 -7.0550
- timer: Routine,Calls,Time,% = IterSCF 1 41.162 97.32
- elaps: Routine,Calls,Wall,% = IterSCF 1 88.245 98.73
- siesta: 2 -6995.5710 -6979.1304 -6979.2076 1.1332 -3.0040
- siesta: 3 -6988.4334 -6984.7326 -6984.7573 0.8727 -3.5989
- siesta: 4 -6988.3877 -6984.1706 -6984.2373 0.9974 -3.6240
- siesta: 5 -6987.4963 -6984.3120 -6984.3197 0.8479 -3.7722
- siesta: 6 -6989.6064 -6983.4914 -6983.5005 0.9377 -3.7433
- siesta: 7 -6989.7368 -6983.5914 -6983.5950 0.9438 -3.8114
- siesta: 8 -6989.5617 -6983.3162 -6983.3220 1.0130 -4.2295
- siesta: 9 -6989.1057 -6983.3055 -6983.3092 0.6526 -4.1951
- siesta: 10 -6987.7380 -6983.6310 -6983.6501 1.1446 -4.4989
- siesta: 11 -6985.9660 -6983.6660 -6983.6660 0.6916 -4.6073
- siesta: 12 -6986.0314 -6983.7028 -6983.7306 0.6918 -4.5761
- siesta: 13 -6985.7867 -6984.2148 -6984.2417 0.6244 -4.5442
- siesta: 14 -6985.8631 -6984.4144 -6984.4421 0.6860 -4.5003
- siesta: 15 -6986.5332 -6985.7475 -6985.7715 0.4947 -4.3396
- siesta: 16 -6987.1387 -6986.6318 -6986.6515 1.2129 -4.4570
- siesta: 17 -6985.6345 -6984.7151 -6984.7153 0.5210 -4.5470
- siesta: 18 -6985.6483 -6984.6991 -6984.7361 0.5439 -4.5495
- siesta: 19 -6985.6643 -6984.6744 -6984.7113 0.5714 -4.5571
- siesta: 20 -6985.5710 -6984.7380 -6984.7744 0.4909 -4.5627
- siesta: 21 -6985.4536 -6984.9508 -6984.9961 0.3452 -4.5589
- siesta: 22 -6985.4983 -6984.8858 -6984.9605 0.3111 -4.5938
- siesta: 23 -6985.4994 -6984.8765 -6984.9374 0.3309 -4.5808
- siesta: 24 -6985.4630 -6984.9673 -6985.0274 0.3605 -4.5691
- siesta: 25 -6985.4709 -6984.8944 -6984.9682 0.4198 -4.5718
- siesta: 26 -6985.5994 -6984.9023 -6984.9772 0.4844 -4.5809
- siesta: 27 -6985.7736 -6984.7111 -6984.7536 0.5681 -4.5390
- siesta: 28 -6985.7722 -6984.7118 -6984.7431 0.5679 -4.5415
- siesta: 29 -6985.7693 -6984.5999 -6984.6312 0.5676 -4.5211
- siesta: 30 -6985.7596 -6984.5913 -6984.6233 0.5880 -4.5384
- ^C
- mahmood@cluster:A$ date
- Thu Oct 6 21:07:18 IRST 2016
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement