hello siesta user, i am running siesta in parallel AIX machine. I use the Diag.ParallelOverK .true. and then run the siesta... i use two processor to run siesta... but i saw that the siesta run will take same time, when i dont use the above option... can any body suggest me how to reduce the computational time..
--------------------------------- Heres a new way to find what you're looking for - Yahoo! Answers
Sender: LSF System <[EMAIL PROTECTED]> Subject: Job 161390: </afs/enea.it/software/bin/poe.lsf shell.sh -procs 2 -sharedmemory yes> Done Job </afs/enea.it/software/bin/poe.lsf shell.sh -procs 2 -sharedmemory yes> was submitted from host <sp5-1.frascati.enea.it> by user <rakshit>. Job was executed on host(s) <2*sp4-3.frascati.enea.it>, in queue <large_72h>, as user <rakshit>. </afs/enea.it/fra/user/rakshit> was used as the home directory. </afs/enea.it/fra/user/rakshit/aix/siesta-2.0/ptn/ptn_LDA_RS> was used as the working directory. Started at Thu Mar 22 12:42:38 2007 Results reported at Thu Mar 22 13:26:10 2007 Your job looked like: ------------------------------------------------------------ # LSBATCH: User input /afs/enea.it/software/bin/poe.lsf shell.sh -procs 2 -sharedmemory yes ------------------------------------------------------------ Successfully completed. Resource usage summary: CPU time : 0.36 sec. Max Memory : 4 MB Max Swap : 6 MB Max Processes : 4 Max Threads : 4 The output (if any) follows: Siesta Version: siesta-2.0-release Architecture : powerpc-ibm-aix5.2.0.0--Xlf Compiler flags: mpxlf95 -O3 -qarch=auto -qtune=auto -qcache=auto -q64 -qfixed -qnolm -qzerosize -qstrict -qmaxmem=-1 PARALLEL version * Running on 2 nodes in parallel >> Start of run: 22-MAR-2007 12:42:40 *********************** * WELCOME TO SIESTA * *********************** reinit: Reading from standard input ************************** Dump of input data file **************************** # $Id: ptn.fdf,v 1.1 1999/04/20 14:43:44 emilio Exp $ # ----------------------------------------------------------------------------- # FDF fo # # E. Artacho, April 1999 # ----------------------------------------------------------------------------- SystemName ptnRS SystemLabel ptnRS NumberOfAtoms 2 NumberOfSpecies 2 %block ChemicalSpeciesLabel 1 78 Pt 2 7 N %endblock ChemicalSpeciesLabel PAO.BasisType split #PAO.BasisSize DZP PAO.EnergyShift 0.1 eV PAO.SplitNorm 0.2000 %block PAO.Basis # Define Basis set Pt 2 # Species label, number of l-shells n=6 0 2 P 1 # n, l, Nzeta, Polarization, NzetaPol 6.982 5.645 1.000 1.000 n=5 2 2 # n, l, Nzeta 5.044 2.803 1.000 1.000 N 3 -0.00139 n=2 0 2 E 65.50216 4.29661 5.64483 3.02914 1.00 1.00 n=2 1 2 E 30.54417 5.81284 7.25855 2.85547 1.00 1.00 n=3 2 1 E 59.15335 0.14049 3.65788 1.00 %endblock PAO.Basis LatticeConstant 4.8041 Ang #%block LatticeParameters # 3.92 3.92 3.92 60.0 60.0 60.0 #%endblock LatticeParameters %block LatticeVectors 0.5 0.5 0.0 0.5 0.0 0.5 0.0 0.5 0.5 %endblock LatticeVectors MeshCutoff 250.00 Ry # SCF options MaxSCFIterations 200 # Maximum number of SCF iter DM.MixingWeight 0.3 # New DM amount for next SCF cycle DM.Tolerance 1.d-4 # Tolerance in maximum difference DM.NumberPulay 8 # Number of pulay mixing steps DM.UseSaveDM .false. # tells if already existing density matrix is to be used or not WriteCoorXmol WriteMullikenPop 1 WriteForces .true. ElectronicTemperature 30 meV xc.functional LDA xc.authors CA # WriteCoorStep .true. #AtomCoorFormatOut Ang SolutionMethod Diagon # OrderN or Diagon AtomicCoordinatesFormat Fractional %block AtomicCoordinatesAndAtomicSpecies 0.0000 0.0000 0.0000 1 0.5000 0.5000 0.5000 2 %endblock AtomicCoordinatesAndAtomicSpecies MD.TypeOfRun CG # Type of dynamics: MD.NumCGsteps 180 # Number of CG steps for MD.MaxCGDispl 0.4 Ang # Maximum atomic displacement MD.MaxForceTol 0.01 eV/Ang # Tolerance in the maximum MD.MaxStressTol 0.1 GPa MD.VariableCell .true. Diag.PatallelOverK .true. %block kgrid_Monkhorst_Pack 12 0 0 0.0 0 12 0 0.0 0 0 12 0.0 %endblock kgrid_Monkhorst_Pack ************************** End of input data file ***************************** reinit: ----------------------------------------------------------------------- reinit: System Name: ptnRS reinit: ----------------------------------------------------------------------- reinit: System Label: ptnRS reinit: ----------------------------------------------------------------------- initatom: Reading input for the pseudopotentials and atomic orbitals ---------- Species number: 1 Label: Pt Atomic number: 78 Species number: 2 Label: N Atomic number: 7 Ground state valence configuration: 5d10 Reading pseudopotential information in formatted form from Pt.psf Ground state valence configuration: 2s02 2p03 Reading pseudopotential information in formatted form from N.psf For Pt, 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... For N, 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> =============================================================================== Pt Z= 78 Mass= 195.08 Charge= 0.0000 Lmxo=2 Lmxkb=3 BasisType=split Semic=F L=0 Nsemic=0 Cnfigmx=6 n=1 nzeta=2 polorb=1 vcte: 0.0000 rinn: 0.0000 rcs: 6.9820 5.6450 lambdas: 1.0000 1.0000 L=1 Nsemic=0 Cnfigmx=6 L=2 Nsemic=0 Cnfigmx=5 n=1 nzeta=2 polorb=0 vcte: 0.0000 rinn: 0.0000 rcs: 5.0440 2.8030 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 Pt (Z = 78) read_vps: Pseudopotential generation method: read_vps: ATM 3.2.2 Troullier-Martins read_vps: Pseudopotential generated from a relativistic atomic calculation read_vps: There are spin-orbit pseudopotentials available read_vps: Spin-orbit interaction is not included in this calculation read_vps: Valence configuration (pseudopotential and basis set generation): 6s( 1.00) rc: 2.44 6p( 0.00) rc: 2.73 5d( 9.00) rc: 2.84 5f( 0.00) rc: 2.44 Total valence charge: 10.00000 xc_check: Exchange-correlation functional: xc_check: Ceperley-Alder V l=0 = -2*Zval/r beyond r= 3.8798 V l=1 = -2*Zval/r beyond r= 3.8798 V l=2 = -2*Zval/r beyond r= 3.8798 V l=3 = -2*Zval/r beyond r= 3.8798 All V_l potentials equal beyond r= 2.8033 This should be close to max(r_c) in ps generation All pots = -2*Zval/r beyond r= 3.8798 Using large-core scheme for Vlocal atom: Estimated core radius 3.87980 atom: Including non-local core corrections could be a good idea atom: Maximum radius for 4*pi*r*r*local-pseudopot. charge 4.23458 atom: Maximum radius for r*vlocal+2*Zval: 3.92860 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.910370 el= -0.436769 Ekb= 4.269138 kbcos= 0.200135 l= 1 rc= 2.910370 el= -0.072130 Ekb= 1.476375 kbcos= 0.245690 l= 2 rc= 2.874217 el= -0.471638 Ekb= -5.183329 kbcos= -0.722504 l= 3 rc= 2.910370 el= 0.003524 Ekb= -1.958140 kbcos= -0.022829 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 6s izeta = 1 lambda = 1.000000 rc = 6.981659 energy = -0.429931 kinetic = 0.358719 potential(screened) = -0.788650 potential(ionic) = -7.462152 izeta = 2 rmatch = 5.716093 splitnorm = 0.182532 energy = -0.343112 kinetic = 0.625502 potential(screened) = -0.968614 potential(ionic) = -8.134943 SPLIT: Orbitals with angular momentum L= 2 SPLIT: Basis orbitals for state 5d izeta = 1 lambda = 1.000000 rc = 5.044431 energy = -0.464168 kinetic = 4.254891 potential(screened) = -4.719060 potential(ionic) = -13.831024 izeta = 2 rmatch = 2.838512 splitnorm = 0.190633 energy = -0.091487 kinetic = 6.301344 potential(screened) = -6.392830 potential(ionic) = -16.358932 POLgen: Perturbative polarization orbital with L= 1 POLgen: Polarization orbital for state 6s izeta = 1 rc = 6.981659 energy = -0.106926 kinetic = 0.637884 potential(screened) = -0.744810 potential(ionic) = -7.029134 atom: Total number of Sankey-type orbitals: 15 atm_pop: Valence configuration(local Pseudopot. screening): 6s( 0.00) 6p( 0.00) 5d(10.00) Vna: chval, zval: 10.00000 10.00000 Vna: Cut-off radius for the neutral-atom potential: 5.044431 atom: _________________________________________________________________________ <basis_specs> =============================================================================== N Z= 7 Mass= 14.010 Charge=-0.13900E-02 Lmxo=2 Lmxkb=3 BasisType=split Semic=F L=0 Nsemic=0 Cnfigmx=2 n=1 nzeta=2 polorb=0 vcte: 65.502 rinn: 4.2966 rcs: 5.6448 3.0291 lambdas: 1.0000 1.0000 L=1 Nsemic=0 Cnfigmx=2 n=1 nzeta=2 polorb=0 vcte: 30.544 rinn: 5.8128 rcs: 7.2585 2.8555 lambdas: 1.0000 1.0000 L=2 Nsemic=0 Cnfigmx=3 n=1 nzeta=1 polorb=0 vcte: 59.153 rinn: 0.14049 rcs: 3.6579 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 N (Z = 7) read_vps: Pseudopotential generation method: read_vps: ATM 3.2.2 Troullier-Martins read_vps: Pseudopotential generated from a relativistic atomic calculation read_vps: There are spin-orbit pseudopotentials available read_vps: Spin-orbit interaction is not included in this calculation read_vps: Valence configuration (pseudopotential and basis set generation): 2s( 2.00) rc: 0.99 2p( 3.00) rc: 0.99 3d( 0.00) rc: 0.99 4f( 0.00) rc: 0.99 Total valence charge: 5.00000 xc_check: Exchange-correlation functional: xc_check: Ceperley-Alder V l=0 = -2*Zval/r beyond r= 1.1662 V l=1 = -2*Zval/r beyond r= 1.1662 V l=2 = -2*Zval/r beyond r= 1.1662 V l=3 = -2*Zval/r beyond r= 1.1662 All V_l potentials equal beyond r= 0.9790 This should be close to max(r_c) in ps generation All pots = -2*Zval/r beyond r= 1.1662 VLOCAL1: 99.0% of the norm of Vloc inside 45.294 Ry VLOCAL1: 99.9% of the norm of Vloc inside 103.226 Ry atom: Maximum radius for 4*pi*r*r*local-pseudopot. charge 1.19577 atom: Maximum radius for r*vlocal+2*Zval: 1.02916 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.372085 el= -1.353310 Ekb= 12.245155 kbcos= 0.275311 l= 1 rc= 1.389349 el= -0.532335 Ekb= -7.229919 kbcos= -0.278452 l= 2 rc= 1.675947 el= 0.002346 Ekb= -2.063346 kbcos= -0.002640 l= 3 rc= 1.947230 el= 0.003449 Ekb= -1.004017 kbcos= -0.000198 KBgen: Total number of Kleinman-Bylander projectors: 16 atom: ------------------------------------------------------------------------- atom: SANKEY-TYPE ORBITALS: atom: Selected multiple-zeta basis: split atom: basis set generated (by rescaling the valence charge) atom: for an anion of charge -0.0014 SPLIT: Orbitals with angular momentum L= 0 SPLIT: Basis orbitals for state 2s izeta = 1 lambda = 1.000000 rc = 5.635168 energy = -1.352726 kinetic = 1.209558 potential(screened) = -2.562284 potential(ionic) = -8.066652 izeta = 2 rmatch = 3.016124 splitnorm = 0.141580 energy = -1.172872 kinetic = 1.899409 potential(screened) = -3.072281 potential(ionic) = -9.001375 SPLIT: Orbitals with angular momentum L= 1 SPLIT: Basis orbitals for state 2p izeta = 1 lambda = 1.000000 rc = 7.235800 energy = -0.531938 kinetic = 3.571754 potential(screened) = -4.103692 potential(ionic) = -9.418660 izeta = 2 rmatch = 2.869008 splitnorm = 0.205560 energy = -0.267318 kinetic = 5.900435 potential(screened) = -6.167752 potential(ionic) = -12.282486 SPLIT: Orbitals with angular momentum L= 2 SPLIT: Basis orbitals for state 3d izeta = 1 lambda = 1.000000 rc = 3.683980 energy = 3.372639 kinetic = 6.107250 potential(screened) = -2.734611 potential(ionic) = -8.093890 atom: Total number of Sankey-type orbitals: 13 atm_pop: Valence configuration(local Pseudopot. screening): 2s( 2.00) 2p( 3.00) 3d( 0.00) Vna: chval, zval: 5.00000 5.00000 Vna: Cut-off radius for the neutral-atom potential: 7.235800 atom: _________________________________________________________________________ prinput: Basis input ---------------------------------------------------------- PAO.BasisType split %block ChemicalSpeciesLabel 1 78 Pt # Species index, atomic number, species label 2 7 N # Species index, atomic number, species label %endblock ChemicalSpeciesLabel %block PAO.Basis # Define Basis set Pt 2 # Species label, number of l-shells n=6 0 2 P 1 # n, l, Nzeta, Polarization, NzetaPol 6.982 5.716 1.000 1.000 n=5 2 2 # n, l, Nzeta 5.044 2.839 1.000 1.000 N 3 -0.001 # Label, l-shells, ionic net charge n=2 0 2 # n, l, Nzeta 5.635 3.016 1.000 1.000 n=2 1 2 # n, l, Nzeta 7.236 2.869 1.000 1.000 n=3 2 1 # n, l, Nzeta 3.684 1.000 %endblock PAO.Basis prinput: ---------------------------------------------------------------------- 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: coor: Atomic-coordinates input format = Fractional 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 = 250.0000 Ry redata: Net charge of the system = 0.0000 |e| redata: Max. number of SCF Iter = 200 redata: Performing Pulay mixing using = 8 iterations redata: Mix DM in first SCF step ? = F redata: Write Pulay info on disk? = F redata: New DM Mixing Weight = 0.3000 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.000100 redata: Require Energy convergence for SCF = F redata: DM 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 = T redata: Electronic Temperature = 0.0022 Ry redata: Fix the spin of the system = F redata: Dynamics option = CG coord. optimization redata: Variable cell = T redata: Use continuation files for CG = F redata: Max atomic displ per move = 0.7559 Bohr redata: Maximum number of CG moves = 180 redata: Force tolerance = 0.0004 Ry/Bohr redata: Stress tolerance = 0.1000 GPa redata: *********************************************************************** siesta: Atomic coordinates (Bohr) and species siesta: 0.00000 0.00000 0.00000 1 1 siesta: 4.53922 4.53922 4.53922 2 2 initatomlists: Number of atoms, orbitals, and projectors: 2 28 32 siesta: System type = bulk * ProcessorY, Blocksize: 1 14 siesta: k-grid: Number of k-points = 1008 siesta: k-grid: Cutoff = 20.382 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 superc: Internal auxiliary supercell: 7 x 7 x 7 = 343 superc: Number of atoms, orbitals, and projectors: 686 9604 10976 * Maximum dynamic memory allocated = 1 MB siesta: ============================== Begin CG move = 0 ============================== superc: Internal auxiliary supercell: 7 x 7 x 7 = 343 superc: Number of atoms, orbitals, and projectors: 686 9604 10976 outcell: Unit cell vectors (Ang): 2.402050 2.402050 0.000000 2.402050 0.000000 2.402050 0.000000 2.402050 2.402050 outcell: Cell vector modules (Ang) : 3.397012 3.397012 3.397012 outcell: Cell angles (23,13,12) (deg): 60.0000 60.0000 60.0000 outcell: Cell volume (Ang**3) : 27.7189 InitMesh: MESH = 30 x 30 x 30 = 27000 InitMesh: Mesh cutoff (required, used) = 250.000 323.327 Ry * Maximum dynamic memory allocated = 30 MB stepf: Fermi-Dirac step function siesta: Program's energy decomposition (eV): siesta: Eions = 1401.175675 siesta: Ena = 242.829444 siesta: Ekin = 757.621504 siesta: Enl = -374.462935 siesta: DEna = 0.000000 siesta: DUscf = 0.000000 siesta: DUext = 0.000000 siesta: Exc = -189.856809 siesta: eta*DQ = 0.000000 siesta: Emadel = 0.000000 siesta: Ekinion = 0.000000 siesta: Eharris = -986.005833 siesta: Etot = -965.044470 siesta: FreeEng = -965.044470 siesta: iscf Eharris(eV) E_KS(eV) FreeEng(eV) dDmax Ef(eV) siesta: 1 -986.0058 -965.0445 -965.0445 1.3929 0.2663 timer: Routine,Calls,Time,% = IterSCF 1 69.940 88.32 elaps: Routine,Calls,Wall,% = IterSCF 1 33.420 85.50 siesta: 2 -1065.2342 -937.5528 -937.5665 9.3320 1.2211 siesta: 3 -976.4276 -965.6816 -965.6825 1.2871 -1.9678 siesta: 4 -974.1184 -967.3168 -967.3224 0.7175 -2.5677 siesta: 5 -974.0781 -968.5805 -968.5860 0.7283 -2.8238 siesta: 6 -973.9512 -971.2537 -971.2590 0.5139 -3.1146 siesta: 7 -973.9598 -972.2561 -972.2623 0.3870 -3.1078 siesta: 8 -973.9388 -973.2159 -973.2220 0.2271 -3.1011 siesta: 9 -973.9380 -973.2897 -973.2956 0.2031 -3.1002 siesta: 10 -973.9339 -973.8802 -973.8861 0.0101 -3.1018 siesta: 11 -973.9339 -973.8914 -973.8970 0.0050 -3.1015 siesta: 12 -973.9339 -973.9413 -973.9469 0.0012 -3.1021 siesta: 13 -973.9339 -973.9409 -973.9465 0.0011 -3.1022 siesta: 14 -973.9339 -973.9286 -973.9341 0.0007 -3.1014 siesta: 15 -973.9339 -973.9331 -973.9387 0.0001 -3.1014 siesta: E_KS(eV) = -973.9339 siesta: E_KS - E_eggbox = -973.9339 siesta: Atomic forces (eV/Ang): 1 0.000002 -0.000004 -0.000006 2 -0.000186 -0.000185 -0.000188 ---------------------------------------- Tot -0.000184 -0.000189 -0.000194 ---------------------------------------- Max 0.000188 Res 0.000132 sqrt( Sum f_i^2 / 3N ) ---------------------------------------- Max 0.000188 constrained Stress-tensor-Voigt (kbar): 352.03 352.03 352.03 -0.00 -0.00 0.00 Target enthalpy (eV/cell) -973.9394 siesta: Stress tensor (static) (eV/Ang**3): 0.219718 -0.000000 0.000000 -0.000000 0.219718 -0.000000 0.000000 -0.000000 0.219717 siesta: Pressure (static): -352.03010112 kBar siesta: Stress tensor (total) (eV/Ang**3): 0.219718 -0.000000 0.000000 -0.000000 0.219718 -0.000000 0.000000 -0.000000 0.219717 siesta: Pressure (total): -352.03010112 kBar mulliken: Atomic and Orbital Populations: Species: Pt Atom Qatom Qorb 6s 6s 5dxy 5dyz 5dz2 5dxz 5dx2-y2 5dxy 5dyz 5dz2 5dxz 5dx2-y2 6Ppy 6Ppz 6Ppx 1 9.263 -0.125 0.615 1.906 1.906 1.286 1.906 1.286 0.039 0.039 0.029 0.039 0.029 0.102 0.102 0.102 Species: N Atom Qatom Qorb 2s 2s 2py 2pz 2px 2py 2pz 2px 3dxy 3dyz 3dz2 3dxz 3dx2-y2 2 5.737 2.124 -0.257 1.410 1.410 1.410 -0.123 -0.123 -0.123 0.002 0.002 0.002 0.002 0.002 mulliken: Qtot = 15.000 cgvc: No target stress found, assuming hydrostatic MD.TargetPressure. cgvc: Target stress (kBar) cgvc: -0.000 0.000 0.000 cgvc: 0.000 -0.000 0.000 cgvc: 0.000 0.000 -0.000 * Maximum dynamic memory allocated = 30 MB siesta: ============================== Begin CG move = 1 ============================== siesta: k-grid: Number of k-points = 1008 siesta: k-grid: Cutoff = 19.441 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 superc: Internal auxiliary supercell: 7 x 7 x 7 = 343 superc: Number of atoms, orbitals, and projectors: 686 9604 10976 outcell: Unit cell vectors (Ang): 2.291104 2.291104 -0.000000 2.291104 0.000000 2.291104 -0.000000 2.291104 2.291104 outcell: Cell vector modules (Ang) : 3.240110 3.240110 3.240110 outcell: Cell angles (23,13,12) (deg): 60.0000 60.0000 60.0000 outcell: Cell volume (Ang**3) : 24.0527 InitMesh: MESH = 30 x 30 x 30 = 27000 InitMesh: Mesh cutoff (required, used) = 250.000 355.399 Ry siesta: iscf Eharris(eV) E_KS(eV) FreeEng(eV) dDmax Ef(eV) siesta: 1 -974.6874 -976.2344 -976.2399 0.3607 -2.4287 siesta: 2 -975.5534 -974.3392 -974.3429 0.7283 -2.0519 siesta: 3 -974.6168 -975.1079 -975.1108 0.1510 -2.2994 siesta: 4 -974.6064 -974.3879 -974.3921 0.0201 -2.2446 siesta: 5 -974.6052 -974.5583 -974.5629 0.0051 -2.2828 siesta: 6 -974.6052 -974.6020 -974.6066 0.0033 -2.2893 siesta: 7 -974.6051 -974.5921 -974.5967 0.0009 -2.2816 siesta: 8 -974.6051 -974.6044 -974.6090 0.0001 -2.2818 siesta: E_KS(eV) = -974.6058 siesta: Atomic forces (eV/Ang): 1 -0.000007 -0.000019 -0.000024 2 -0.000233 -0.000231 -0.000189 ---------------------------------------- Tot -0.000240 -0.000250 -0.000214 ---------------------------------------- Max 0.000233 Res 0.000155 sqrt( Sum f_i^2 / 3N ) ---------------------------------------- Max 0.000233 constrained Stress-tensor-Voigt (kbar): 209.01 209.01 209.01 0.00 0.00 -0.00 Target enthalpy (eV/cell) -974.6104 siesta: Stress tensor (static) (eV/Ang**3): 0.130451 0.000000 -0.000000 0.000000 0.130451 0.000000 -0.000000 0.000000 0.130452 siesta: Pressure (static): -209.00794191 kBar siesta: Stress tensor (total) (eV/Ang**3): 0.130451 0.000000 -0.000000 0.000000 0.130451 0.000000 -0.000000 0.000000 0.130452 siesta: Pressure (total): -209.00794191 kBar mulliken: Atomic and Orbital Populations: Species: Pt Atom Qatom Qorb 6s 6s 5dxy 5dyz 5dz2 5dxz 5dx2-y2 5dxy 5dyz 5dz2 5dxz 5dx2-y2 6Ppy 6Ppz 6Ppx 1 9.150 -0.407 0.836 1.931 1.931 1.270 1.931 1.270 0.020 0.020 0.020 0.020 0.020 0.096 0.096 0.096 Species: N Atom Qatom Qorb 2s 2s 2py 2pz 2px 2py 2pz 2px 3dxy 3dyz 3dz2 3dxz 3dx2-y2 2 5.850 2.200 -0.356 1.485 1.485 1.485 -0.155 -0.155 -0.155 0.003 0.003 0.004 0.003 0.004 mulliken: Qtot = 15.000 * Maximum dynamic memory allocated = 31 MB siesta: ============================== Begin CG move = 2 ============================== siesta: k-grid: Number of k-points = 1008 siesta: k-grid: Cutoff = 17.934 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 superc: Internal auxiliary supercell: 8 x 8 x 8 = 512 superc: Number of atoms, orbitals, and projectors: 1024 14336 16384 outcell: Unit cell vectors (Ang): 2.113590 2.113590 -0.000000 2.113590 0.000000 2.113591 -0.000000 2.113590 2.113591 outcell: Cell vector modules (Ang) : 2.989068 2.989069 2.989069 outcell: Cell angles (23,13,12) (deg): 60.0000 60.0000 60.0000 outcell: Cell volume (Ang**3) : 18.8839 InitMesh: MESH = 24 x 24 x 24 = 13824 InitMesh: Mesh cutoff (required, used) = 250.000 267.266 Ry siesta: iscf Eharris(eV) E_KS(eV) FreeEng(eV) dDmax Ef(eV) siesta: 1 -974.2032 -976.1242 -976.1288 1.2301 -0.2736 siesta: 2 -977.6756 -972.8164 -972.8178 1.3444 -0.1606 siesta: 3 -974.0662 -973.9468 -973.9510 0.2227 -0.4806 siesta: 4 -974.1918 -973.6159 -973.6188 0.0841 -0.5165 siesta: 5 -974.1949 -973.9562 -973.9593 0.0263 -0.4999 siesta: 6 -974.1998 -974.1210 -974.1241 0.0066 -0.4953 siesta: 7 -974.1998 -974.1979 -974.2010 0.0003 -0.4978 siesta: 8 -974.1998 -974.2026 -974.2057 0.0002 -0.4979 siesta: 9 -974.1998 -974.2010 -974.2041 0.0001 -0.4979 siesta: E_KS(eV) = -974.1998 siesta: Atomic forces (eV/Ang): 1 -0.002340 -0.002335 -0.002319 2 0.002011 0.001997 0.001986 ---------------------------------------- Tot -0.000329 -0.000338 -0.000333 ---------------------------------------- Max 0.002340 Res 0.002171 sqrt( Sum f_i^2 / 3N ) ---------------------------------------- Max 0.002340 constrained Stress-tensor-Voigt (kbar): -655.86 -655.86 -655.86 0.00 0.00 0.00 Target enthalpy (eV/cell) -974.2029 siesta: Stress tensor (static) (eV/Ang**3): -0.409352 0.000000 0.000001 0.000000 -0.409351 0.000000 0.000001 0.000000 -0.409351 siesta: Pressure (static): 655.86026584 kBar siesta: Stress tensor (total) (eV/Ang**3): -0.409352 0.000000 0.000001 0.000000 -0.409351 0.000000 0.000001 0.000000 -0.409351 siesta: Pressure (total): 655.86026584 kBar mulliken: Atomic and Orbital Populations: Species: Pt Atom Qatom Qorb 6s 6s 5dxy 5dyz 5dz2 5dxz 5dx2-y2 5dxy 5dyz 5dz2 5dxz 5dx2-y2 6Ppy 6Ppz 6Ppx 1 9.051 -0.759 1.180 2.006 2.006 1.258 2.006 1.258 -0.045 -0.045 0.006 -0.045 0.006 0.073 0.073 0.073 Species: N Atom Qatom Qorb 2s 2s 2py 2pz 2px 2py 2pz 2px 3dxy 3dyz 3dz2 3dxz 3dx2-y2 2 5.949 2.353 -0.548 1.546 1.546 1.546 -0.176 -0.176 -0.176 0.005 0.005 0.009 0.005 0.009 mulliken: Qtot = 15.000 * Maximum dynamic memory allocated = 35 MB siesta: ============================== Begin CG move = 3 ============================== siesta: k-grid: Number of k-points = 1008 siesta: k-grid: Cutoff = 19.006 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 superc: Internal auxiliary supercell: 7 x 7 x 7 = 343 superc: Number of atoms, orbitals, and projectors: 686 9604 10976 outcell: Unit cell vectors (Ang): 2.239853 2.239853 -0.000000 2.239853 0.000000 2.239854 -0.000000 2.239853 2.239854 outcell: Cell vector modules (Ang) : 3.167631 3.167631 3.167631 outcell: Cell angles (23,13,12) (deg): 60.0000 60.0000 60.0000 outcell: Cell volume (Ang**3) : 22.4744 InitMesh: MESH = 30 x 30 x 30 = 27000 InitMesh: Mesh cutoff (required, used) = 250.000 371.849 Ry siesta: iscf Eharris(eV) E_KS(eV) FreeEng(eV) dDmax Ef(eV) siesta: 1 -973.9187 -972.1939 -972.1969 0.8458 -1.8908 siesta: 2 -976.1450 -974.2839 -974.2873 0.6746 -1.6376 siesta: 3 -974.6628 -974.0948 -974.0972 0.1678 -1.7804 siesta: 4 -974.7406 -974.4597 -974.4637 0.0254 -1.8147 siesta: 5 -974.7406 -974.4695 -974.4729 0.0240 -1.8162 siesta: 6 -974.7365 -974.6674 -974.6708 0.0083 -1.8440 siesta: 7 -974.7366 -974.7304 -974.7338 0.0008 -1.8439 siesta: 8 -974.7366 -974.7372 -974.7406 0.0001 -1.8439 siesta: 9 -974.7366 -974.7367 -974.7400 0.0000 -1.8439 siesta: E_KS(eV) = -974.7365 siesta: Atomic forces (eV/Ang): 1 -0.000275 -0.000279 -0.000268 2 0.000157 0.000154 0.000159 ---------------------------------------- Tot -0.000118 -0.000126 -0.000109 ---------------------------------------- Max 0.000279 Res 0.000223 sqrt( Sum f_i^2 / 3N ) ---------------------------------------- Max 0.000279 constrained Stress-tensor-Voigt (kbar): 51.94 51.94 51.94 -0.00 -0.00 0.00 Target enthalpy (eV/cell) -974.7399 siesta: Stress tensor (static) (eV/Ang**3): 0.032416 -0.000001 0.000000 -0.000001 0.032415 -0.000000 0.000000 -0.000000 0.032416 siesta: Pressure (static): -51.93603020 kBar siesta: Stress tensor (total) (eV/Ang**3): 0.032416 -0.000001 0.000000 -0.000001 0.032415 -0.000000 0.000000 -0.000000 0.032416 siesta: Pressure (total): -51.93603020 kBar mulliken: Atomic and Orbital Populations: Species: Pt Atom Qatom Qorb 6s 6s 5dxy 5dyz 5dz2 5dxz 5dx2-y2 5dxy 5dyz 5dz2 5dxz 5dx2-y2 6Ppy 6Ppz 6Ppx 1 9.101 -0.539 0.949 1.946 1.946 1.264 1.946 1.264 0.006 0.006 0.017 0.006 0.017 0.090 0.090 0.090 Species: N Atom Qatom Qorb 2s 2s 2py 2pz 2px 2py 2pz 2px 3dxy 3dyz 3dz2 3dxz 3dx2-y2 2 5.899 2.251 -0.414 1.515 1.515 1.515 -0.168 -0.168 -0.168 0.003 0.003 0.005 0.003 0.005 mulliken: Qtot = 15.000 cgvc: Finished line minimization 1. Mean atomic displacement = 0.7814 * Maximum dynamic memory allocated = 35 MB siesta: ============================== Begin CG move = 4 ============================== siesta: k-grid: Number of k-points = 1008 siesta: k-grid: Cutoff = 17.630 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 superc: Internal auxiliary supercell: 8 x 8 x 8 = 512 superc: Number of atoms, orbitals, and projectors: 1024 14336 16384 outcell: Unit cell vectors (Ang): 2.077738 2.077754 0.000019 2.077698 0.000059 2.077704 0.000038 2.077735 2.077725 outcell: Cell vector modules (Ang) : 2.938376 2.938313 2.938354 outcell: Cell angles (23,13,12) (deg): 59.9985 59.9989 59.9990 outcell: Cell volume (Ang**3) : 17.9384 InitMesh: MESH = 24 x 24 x 24 = 13824 InitMesh: Mesh cutoff (required, used) = 250.000 276.567 Ry siesta: iscf Eharris(eV) E_KS(eV) FreeEng(eV) dDmax Ef(eV) siesta: 1 -972.7317 -976.0091 -976.0125 1.0572 -0.0069 siesta: 2 -976.7356 -972.4322 -972.4335 1.2481 0.4298 siesta: 3 -973.4453 -973.4470 -973.4515 0.2002 -0.0291 siesta: 4 -973.7299 -973.1142 -973.1163 0.0570 -0.0393 siesta: 5 -973.7075 -973.4521 -973.4543 0.0203 -0.0453 siesta: 6 -973.7101 -973.6054 -973.6075 0.0093 -0.0472 siesta: 7 -973.7137 -973.7184 -973.7205 0.0004 -0.0466 siesta: 8 -973.7137 -973.7160 -973.7181 0.0002 -0.0465 siesta: 9 -973.7137 -973.7143 -973.7165 0.0001 -0.0465 siesta: E_KS(eV) = -973.7137 siesta: Atomic forces (eV/Ang): 1 0.284198 0.279964 0.278990 2 -0.285124 -0.285033 -0.284175 ---------------------------------------- Tot -0.000926 -0.005069 -0.005185 ---------------------------------------- Max 0.285124 Res 0.282925 sqrt( Sum f_i^2 / 3N ) ---------------------------------------- Max 0.285124 constrained Stress-tensor-Voigt (kbar): -992.02 -991.99 -992.16 -2.45 -1.53 -2.41 Target enthalpy (eV/cell) -973.7159 siesta: Stress tensor (static) (eV/Ang**3): -0.619166 -0.001529 -0.001502 -0.001529 -0.619146 -0.000953 -0.001502 -0.000953 -0.619251 siesta: Pressure (static): 992.05889402 kBar siesta: Stress tensor (total) (eV/Ang**3): -0.619166 -0.001529 -0.001502 -0.001529 -0.619146 -0.000953 -0.001502 -0.000953 -0.619251 siesta: Pressure (total): 992.05889402 kBar mulliken: Atomic and Orbital Populations: Species: Pt Atom Qatom Qorb 6s 6s 5dxy 5dyz 5dz2 5dxz 5dx2-y2 5dxy 5dyz 5dz2 5dxz 5dx2-y2 6Ppy 6Ppz 6Ppx 1 9.093 -0.769 1.218 2.029 2.029 1.263 2.029 1.263 -0.066 -0.066 0.002 -0.066 0.002 0.075 0.075 0.075 Species: N Atom Qatom Qorb 2s 2s 2py 2pz 2px 2py 2pz 2px 3dxy 3dyz 3dz2 3dxz 3dx2-y2 2 5.907 2.313 -0.541 1.533 1.533 1.533 -0.168 -0.168 -0.168 0.006 0.006 0.011 0.006 0.011 mulliken: Qtot = 15.000 * Maximum dynamic memory allocated = 36 MB siesta: ============================== Begin CG move = 5 ============================== siesta: k-grid: Number of k-points = 1008 siesta: k-grid: Cutoff = 18.921 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 superc: Internal auxiliary supercell: 7 x 7 x 7 = 343 superc: Number of atoms, orbitals, and projectors: 686 9604 10976 outcell: Unit cell vectors (Ang): 2.229924 2.229925 0.000001 2.229921 0.000004 2.229922 0.000002 2.229923 2.229923 outcell: Cell vector modules (Ang) : 3.153589 3.153585 3.153588 outcell: Cell angles (23,13,12) (deg): 59.9999 59.9999 59.9999 outcell: Cell volume (Ang**3) : 22.1768 InitMesh: MESH = 30 x 30 x 30 = 27000 InitMesh: Mesh cutoff (required, used) = 250.000 375.168 Ry siesta: iscf Eharris(eV) E_KS(eV) FreeEng(eV) dDmax Ef(eV) siesta: 1 -971.5134 -969.4902 -969.4923 0.9659 -1.9551 siesta: 2 -975.5636 -974.4558 -974.4594 0.4704 -1.5119 siesta: 3 -974.5019 -974.0990 -974.1021 0.1149 -1.6596 siesta: 4 -974.7537 -974.5773 -974.5810 0.0383 -1.6977 siesta: 5 -974.7515 -974.6526 -974.6559 0.0115 -1.7352 siesta: 6 -974.7440 -974.7009 -974.7042 0.0056 -1.7527 siesta: 7 -974.7439 -974.7429 -974.7462 0.0002 -1.7537 siesta: 8 -974.7438 -974.7436 -974.7468 0.0001 -1.7537 siesta: 9 -974.7438 -974.7438 -974.7471 0.0000 -1.7538 siesta: E_KS(eV) = -974.7439 siesta: Atomic forces (eV/Ang): 1 0.005218 0.005183 0.006368 2 -0.003412 -0.003394 -0.003347 ---------------------------------------- Tot 0.001805 0.001788 0.003021 ---------------------------------------- Max 0.006368 Res 0.004637 sqrt( Sum f_i^2 / 3N ) ---------------------------------------- Max 0.006368 constrained Stress-tensor-Voigt (kbar): 27.37 27.36 27.38 -0.42 -0.24 -0.24 Target enthalpy (eV/cell) -974.7471 siesta: Stress tensor (static) (eV/Ang**3): 0.017081 -0.000263 -0.000148 -0.000263 0.017079 -0.000149 -0.000148 -0.000148 0.017090 siesta: Pressure (static): -27.37068337 kBar siesta: Stress tensor (total) (eV/Ang**3): 0.017081 -0.000263 -0.000148 -0.000263 0.017079 -0.000149 -0.000148 -0.000148 0.017090 siesta: Pressure (total): -27.37068337 kBar mulliken: Atomic and Orbital Populations: Species: Pt Atom Qatom Qorb 6s 6s 5dxy 5dyz 5dz2 5dxz 5dx2-y2 5dxy 5dyz 5dz2 5dxz 5dx2-y2 6Ppy 6Ppz 6Ppx 1 9.094 -0.556 0.965 1.950 1.950 1.263 1.950 1.263 0.003 0.003 0.016 0.003 0.016 0.089 0.089 0.089 Species: N Atom Qatom Qorb 2s 2s 2py 2pz 2px 2py 2pz 2px 3dxy 3dyz 3dz2 3dxz 3dx2-y2 2 5.906 2.257 -0.423 1.519 1.519 1.519 -0.169 -0.169 -0.169 0.003 0.003 0.006 0.003 0.006 mulliken: Qtot = 15.000 cgvc: Finished line minimization 2. Mean atomic displacement = 0.0479 * Maximum dynamic memory allocated = 36 MB siesta: ============================== Begin CG move = 6 ============================== siesta: k-grid: Number of k-points = 1008 siesta: k-grid: Cutoff = 18.836 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 superc: Internal auxiliary supercell: 7 x 7 x 7 = 343 superc: Number of atoms, orbitals, and projectors: 686 9604 10976 outcell: Unit cell vectors (Ang): 2.220225 2.220226 0.000183 2.220152 0.000255 2.220147 0.000254 2.220155 2.220149 outcell: Cell vector modules (Ang) : 3.139873 3.139766 3.139769 outcell: Cell angles (23,13,12) (deg): 59.9925 59.9934 59.9935 outcell: Cell volume (Ang**3) : 21.8839 InitMesh: MESH = 30 x 30 x 30 = 27000 InitMesh: Mesh cutoff (required, used) = 250.000 378.430 Ry siesta: iscf Eharris(eV) E_KS(eV) FreeEng(eV) dDmax Ef(eV) siesta: 1 -974.7473 -974.9508 -974.9541 0.0653 -1.6760 siesta: 2 -974.7548 -974.7434 -974.7465 0.0750 -1.6413 siesta: 3 -974.7463 -974.8061 -974.8095 0.0221 -1.6621 siesta: 4 -974.7461 -974.7364 -974.7396 0.0011 -1.6576 siesta: 5 -974.7461 -974.7459 -974.7490 0.0001 -1.6608 siesta: E_KS(eV) = -974.7458 siesta: Atomic forces (eV/Ang): 1 -0.022381 -0.022223 -0.022800 2 0.018547 0.018505 0.019236 ---------------------------------------- Tot -0.003834 -0.003718 -0.003564 ---------------------------------------- Max 0.022800 Res 0.020700 sqrt( Sum f_i^2 / 3N ) ---------------------------------------- Max 0.022800 constrained Stress-tensor-Voigt (kbar): -6.43 -6.46 -6.52 -0.20 -0.27 -0.29 Target enthalpy (eV/cell) -974.7490 siesta: Stress tensor (static) (eV/Ang**3): -0.004013 -0.000126 -0.000178 -0.000126 -0.004034 -0.000171 -0.000178 -0.000171 -0.004072 siesta: Pressure (static): 6.47269288 kBar siesta: Stress tensor (total) (eV/Ang**3): -0.004013 -0.000126 -0.000178 -0.000126 -0.004034 -0.000171 -0.000178 -0.000171 -0.004072 siesta: Pressure (total): 6.47269288 kBar mulliken: Atomic and Orbital Populations: Species: Pt Atom Qatom Qorb 6s 6s 5dxy 5dyz 5dz2 5dxz 5dx2-y2 5dxy 5dyz 5dz2 5dxz 5dx2-y2 6Ppy 6Ppz 6Ppx 1 9.086 -0.581 0.988 1.953 1.953 1.263 1.953 1.263 0.000 0.000 0.015 0.000 0.015 0.088 0.088 0.088 Species: N Atom Qatom Qorb 2s 2s 2py 2pz 2px 2py 2pz 2px 3dxy 3dyz 3dz2 3dxz 3dx2-y2 2 5.914 2.269 -0.436 1.523 1.523 1.523 -0.171 -0.171 -0.171 0.004 0.004 0.006 0.004 0.006 mulliken: Qtot = 15.000 * Maximum dynamic memory allocated = 36 MB siesta: ============================== Begin CG move = 7 ============================== siesta: k-grid: Number of k-points = 1008 siesta: k-grid: Cutoff = 18.856 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 superc: Internal auxiliary supercell: 7 x 7 x 7 = 343 superc: Number of atoms, orbitals, and projectors: 686 9604 10976 outcell: Unit cell vectors (Ang): 2.222404 2.222406 0.000142 2.222348 0.000199 2.222344 0.000197 2.222350 2.222345 outcell: Cell vector modules (Ang) : 3.142955 3.142872 3.142875 outcell: Cell angles (23,13,12) (deg): 59.9942 59.9949 59.9949 outcell: Cell volume (Ang**3) : 21.9495 InitMesh: MESH = 30 x 30 x 30 = 27000 InitMesh: Mesh cutoff (required, used) = 250.000 377.693 Ry siesta: iscf Eharris(eV) E_KS(eV) FreeEng(eV) dDmax Ef(eV) siesta: 1 -974.7463 -974.6994 -974.7026 0.0156 -1.6785 siesta: 2 -974.7468 -974.7461 -974.7493 0.0171 -1.6851 siesta: 3 -974.7462 -974.7325 -974.7357 0.0052 -1.6817 siesta: 4 -974.7462 -974.7478 -974.7510 0.0002 -1.6827 siesta: 5 -974.7462 -974.7462 -974.7494 0.0000 -1.6820 siesta: E_KS(eV) = -974.7462 siesta: Atomic forces (eV/Ang): 1 -0.025394 -0.025407 -0.024897 2 0.012831 0.012810 0.013304 ---------------------------------------- Tot -0.012563 -0.012596 -0.011594 ---------------------------------------- Max 0.025407 Res 0.020066 sqrt( Sum f_i^2 / 3N ) ---------------------------------------- Max 0.025407 constrained Stress-tensor-Voigt (kbar): 3.07 3.07 3.03 -0.08 -0.08 -0.08 Target enthalpy (eV/cell) -974.7494 siesta: Stress tensor (static) (eV/Ang**3): 0.001918 -0.000050 -0.000050 -0.000050 0.001917 -0.000049 -0.000050 -0.000049 0.001892 siesta: Pressure (static): -3.05834420 kBar siesta: Stress tensor (total) (eV/Ang**3): 0.001918 -0.000050 -0.000050 -0.000050 0.001917 -0.000049 -0.000050 -0.000049 0.001892 siesta: Pressure (total): -3.05834420 kBar mulliken: Atomic and Orbital Populations: Species: Pt Atom Qatom Qorb 6s 6s 5dxy 5dyz 5dz2 5dxz 5dx2-y2 5dxy 5dyz 5dz2 5dxz 5dx2-y2 6Ppy 6Ppz 6Ppx 1 9.088 -0.574 0.982 1.953 1.953 1.263 1.953 1.263 0.001 0.001 0.015 0.001 0.015 0.088 0.088 0.088 Species: N Atom Qatom Qorb 2s 2s 2py 2pz 2px 2py 2pz 2px 3dxy 3dyz 3dz2 3dxz 3dx2-y2 2 5.912 2.266 -0.433 1.523 1.523 1.523 -0.170 -0.170 -0.170 0.003 0.003 0.006 0.003 0.006 mulliken: Qtot = 15.000 cgvc: Finished line minimization 3. Mean atomic displacement = 0.0371 * Maximum dynamic memory allocated = 36 MB siesta: ============================== Begin CG move = 8 ============================== siesta: k-grid: Number of k-points = 1008 siesta: k-grid: Cutoff = 18.820 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 superc: Internal auxiliary supercell: 7 x 7 x 7 = 343 superc: Number of atoms, orbitals, and projectors: 686 9604 10976 outcell: Unit cell vectors (Ang): 2.218380 2.218382 0.000269 2.218305 0.000344 2.218321 0.000343 2.218308 2.218322 outcell: Cell vector modules (Ang) : 3.137265 3.137168 3.137171 outcell: Cell angles (23,13,12) (deg): 59.9895 59.9910 59.9910 outcell: Cell volume (Ang**3) : 21.8282 InitMesh: MESH = 30 x 30 x 30 = 27000 InitMesh: Mesh cutoff (required, used) = 250.000 379.043 Ry siesta: iscf Eharris(eV) E_KS(eV) FreeEng(eV) dDmax Ef(eV) siesta: 1 -974.7424 -974.8279 -974.8311 0.0503 -1.6373 siesta: 2 -974.7435 -974.7408 -974.7439 0.0344 -1.6230 siesta: 3 -974.7418 -974.7429 -974.7462 0.0065 -1.6299 siesta: 4 -974.7417 -974.7395 -974.7426 0.0031 -1.6295 siesta: 5 -974.7416 -974.7414 -974.7445 0.0001 -1.6307 siesta: 6 -974.7416 -974.7415 -974.7447 0.0001 -1.6307 siesta: E_KS(eV) = -974.7416 siesta: Atomic forces (eV/Ang): 1 0.134452 0.134417 0.134762 2 -0.111413 -0.111457 -0.110681 ---------------------------------------- Tot 0.023039 0.022960 0.024081 ---------------------------------------- Max 0.134762 Res 0.123418 sqrt( Sum f_i^2 / 3N ) ---------------------------------------- Max 0.134762 constrained Stress-tensor-Voigt (kbar): -24.19 -24.19 -24.24 -1.57 -1.60 -1.60 Target enthalpy (eV/cell) -974.7448 siesta: Stress tensor (static) (eV/Ang**3): -0.015095 -0.000978 -0.001001 -0.000978 -0.015097 -0.001001 -0.001001 -0.001001 -0.015131 siesta: Pressure (static): 24.20553657 kBar siesta: Stress tensor (total) (eV/Ang**3): -0.015095 -0.000978 -0.001001 -0.000978 -0.015097 -0.001001 -0.001001 -0.001001 -0.015131 siesta: Pressure (total): 24.20553657 kBar mulliken: Atomic and Orbital Populations: Species: Pt Atom Qatom Qorb 6s 6s 5dxy 5dyz 5dz2 5dxz 5dx2-y2 5dxy 5dyz 5dz2 5dxz 5dx2-y2 6Ppy 6Ppz 6Ppx 1 9.085 -0.589 0.995 1.953 1.953 1.264 1.953 1.264 -0.000 -0.000 0.015 -0.000 0.015 0.088 0.088 0.088 Species: N Atom Qatom Qorb 2s 2s 2py 2pz 2px 2py 2pz 2px 3dxy 3dyz 3dz2 3dxz 3dx2-y2 2 5.915 2.272 -0.439 1.524 1.524 1.524 -0.171 -0.171 -0.171 0.004 0.004 0.006 0.004 0.006 mulliken: Qtot = 15.000 * Maximum dynamic memory allocated = 36 MB siesta: ============================== Begin CG move = 9 ============================== siesta: k-grid: Number of k-points = 1008 siesta: k-grid: Cutoff = 18.851 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 superc: Internal auxiliary supercell: 7 x 7 x 7 = 343 superc: Number of atoms, orbitals, and projectors: 686 9604 10976 outcell: Unit cell vectors (Ang): 2.221847 2.221849 0.000160 2.221788 0.000219 2.221787 0.000218 2.221791 2.221788 outcell: Cell vector modules (Ang) : 3.142168 3.142082 3.142085 outcell: Cell angles (23,13,12) (deg): 59.9935 59.9944 59.9944 outcell: Cell volume (Ang**3) : 21.9327 InitMesh: MESH = 30 x 30 x 30 = 27000 InitMesh: Mesh cutoff (required, used) = 250.000 377.880 Ry siesta: iscf Eharris(eV) E_KS(eV) FreeEng(eV) dDmax Ef(eV) siesta: 1 -974.7469 -974.6617 -974.6649 0.0432 -1.6699 siesta: 2 -974.7479 -974.7457 -974.7489 0.0294 -1.6821 siesta: 3 -974.7464 -974.7443 -974.7474 0.0054 -1.6777 siesta: 4 -974.7464 -974.7476 -974.7508 0.0027 -1.6780 siesta: 5 -974.7464 -974.7463 -974.7495 0.0001 -1.6768 siesta: 6 -974.7464 -974.7463 -974.7495 0.0001 -1.6768 siesta: E_KS(eV) = -974.7464 siesta: Atomic forces (eV/Ang): 1 -0.001162 -0.001184 -0.000959 2 -0.003546 -0.003531 -0.002912 ---------------------------------------- Tot -0.004709 -0.004714 -0.003871 ---------------------------------------- Max 0.003546 Res 0.002490 sqrt( Sum f_i^2 / 3N ) ---------------------------------------- Max 0.003546 constrained Stress-tensor-Voigt (kbar): 1.07 1.07 1.08 -0.04 -0.06 -0.06 Target enthalpy (eV/cell) -974.7496 siesta: Stress tensor (static) (eV/Ang**3): 0.000669 -0.000025 -0.000038 -0.000025 0.000668 -0.000038 -0.000038 -0.000038 0.000674 siesta: Pressure (static): -1.07403989 kBar siesta: Stress tensor (total) (eV/Ang**3): 0.000669 -0.000025 -0.000038 -0.000025 0.000668 -0.000038 -0.000038 -0.000038 0.000674 siesta: Pressure (total): -1.07403989 kBar mulliken: Atomic and Orbital Populations: Species: Pt Atom Qatom Qorb 6s 6s 5dxy 5dyz 5dz2 5dxz 5dx2-y2 5dxy 5dyz 5dz2 5dxz 5dx2-y2 6Ppy 6Ppz 6Ppx 1 9.088 -0.576 0.983 1.953 1.953 1.263 1.953 1.263 0.001 0.001 0.015 0.001 0.015 0.088 0.088 0.088 Species: N Atom Qatom Qorb 2s 2s 2py 2pz 2px 2py 2pz 2px 3dxy 3dyz 3dz2 3dxz 3dx2-y2 2 5.912 2.267 -0.434 1.523 1.523 1.523 -0.170 -0.170 -0.170 0.003 0.003 0.006 0.003 0.006 mulliken: Qtot = 15.000 cgvc: Finished line minimization 4. Mean atomic displacement = 0.0051 * Maximum dynamic memory allocated = 36 MB siesta: ============================== Begin CG move = 10 ============================== siesta: k-grid: Number of k-points = 1008 siesta: k-grid: Cutoff = 18.845 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 superc: Internal auxiliary supercell: 7 x 7 x 7 = 343 superc: Number of atoms, orbitals, and projectors: 686 9604 10976 outcell: Unit cell vectors (Ang): 2.221230 2.221232 0.000225 2.221181 0.000274 2.221175 0.000273 2.221184 2.221176 outcell: Cell vector modules (Ang) : 3.141295 3.141220 3.141223 outcell: Cell angles (23,13,12) (deg): 59.9919 59.9925 59.9925 outcell: Cell volume (Ang**3) : 21.9137 InitMesh: MESH = 30 x 30 x 30 = 27000 InitMesh: Mesh cutoff (required, used) = 250.000 378.078 Ry siesta: iscf Eharris(eV) E_KS(eV) FreeEng(eV) dDmax Ef(eV) siesta: 1 -974.7463 -974.7597 -974.7629 0.0045 -1.6717 siesta: 2 -974.7463 -974.7463 -974.7495 0.0050 -1.6697 siesta: 3 -974.7463 -974.7492 -974.7524 0.0013 -1.6707 siesta: 4 -974.7463 -974.7459 -974.7491 0.0001 -1.6705 siesta: 5 -974.7463 -974.7463 -974.7495 0.0000 -1.6707 siesta: E_KS(eV) = -974.7463 siesta: Atomic forces (eV/Ang): 1 0.004082 0.004118 0.004900 2 0.001609 0.001576 0.001819 ---------------------------------------- Tot 0.005690 0.005694 0.006719 ---------------------------------------- Max 0.004900 Res 0.003317 sqrt( Sum f_i^2 / 3N ) ---------------------------------------- Max 0.004900 constrained Stress-tensor-Voigt (kbar): -3.04 -3.04 -3.03 -0.22 -0.24 -0.23 Target enthalpy (eV/cell) -974.7495 siesta: Stress tensor (static) (eV/Ang**3): -0.001896 -0.000135 -0.000145 -0.000135 -0.001894 -0.000147 -0.000145 -0.000147 -0.001893 siesta: Pressure (static): 3.03543681 kBar siesta: Stress tensor (total) (eV/Ang**3): -0.001896 -0.000135 -0.000145 -0.000135 -0.001894 -0.000147 -0.000145 -0.000147 -0.001893 siesta: Pressure (total): 3.03543681 kBar mulliken: Atomic and Orbital Populations: Species: Pt Atom Qatom Qorb 6s 6s 5dxy 5dyz 5dz2 5dxz 5dx2-y2 5dxy 5dyz 5dz2 5dxz 5dx2-y2 6Ppy 6Ppz 6Ppx 1 9.087 -0.578 0.985 1.953 1.953 1.263 1.953 1.263 0.001 0.001 0.015 0.001 0.015 0.088 0.088 0.088 Species: N Atom Qatom Qorb 2s 2s 2py 2pz 2px 2py 2pz 2px 3dxy 3dyz 3dz2 3dxz 3dx2-y2 2 5.913 2.268 -0.435 1.523 1.523 1.523 -0.171 -0.171 -0.171 0.004 0.004 0.006 0.004 0.006 mulliken: Qtot = 15.000 * Maximum dynamic memory allocated = 36 MB siesta: ============================== Begin CG move = 11 ============================== siesta: k-grid: Number of k-points = 1008 siesta: k-grid: Cutoff = 18.848 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 superc: Internal auxiliary supercell: 7 x 7 x 7 = 343 superc: Number of atoms, orbitals, and projectors: 686 9604 10976 outcell: Unit cell vectors (Ang): 2.221610 2.221612 0.000185 2.221555 0.000240 2.221552 0.000239 2.221558 2.221553 outcell: Cell vector modules (Ang) : 3.141833 3.141751 3.141754 outcell: Cell angles (23,13,12) (deg): 59.9929 59.9936 59.9937 outcell: Cell volume (Ang**3) : 21.9254 InitMesh: MESH = 30 x 30 x 30 = 27000 InitMesh: Mesh cutoff (required, used) = 250.000 377.956 Ry siesta: iscf Eharris(eV) E_KS(eV) FreeEng(eV) dDmax Ef(eV) siesta: 1 -974.7463 -974.7380 -974.7412 0.0028 -1.6738 siesta: 2 -974.7464 -974.7463 -974.7495 0.0031 -1.6751 siesta: 3 -974.7463 -974.7446 -974.7478 0.0008 -1.6744 siesta: 4 -974.7463 -974.7466 -974.7497 0.0001 -1.6746 siesta: E_KS(eV) = -974.7463 siesta: Atomic forces (eV/Ang): 1 0.005784 0.005767 0.005588 2 -0.002033 -0.002060 -0.001531 ---------------------------------------- Tot 0.003750 0.003707 0.004058 ---------------------------------------- Max 0.005784 Res 0.004256 sqrt( Sum f_i^2 / 3N ) ---------------------------------------- Max 0.005784 constrained Stress-tensor-Voigt (kbar): -1.28 -1.28 -1.27 0.04 0.05 0.05 Target enthalpy (eV/cell) -974.7495 siesta: Stress tensor (static) (eV/Ang**3): -0.000800 0.000026 0.000029 0.000026 -0.000801 0.000029 0.000029 0.000029 -0.000794 siesta: Pressure (static): 1.27946245 kBar siesta: Stress tensor (total) (eV/Ang**3): -0.000800 0.000026 0.000029 0.000026 -0.000801 0.000029 0.000029 0.000029 -0.000794 siesta: Pressure (total): 1.27946245 kBar mulliken: Atomic and Orbital Populations: Species: Pt Atom Qatom Qorb 6s 6s 5dxy 5dyz 5dz2 5dxz 5dx2-y2 5dxy 5dyz 5dz2 5dxz 5dx2-y2 6Ppy 6Ppz 6Ppx 1 9.088 -0.577 0.984 1.953 1.953 1.263 1.953 1.263 0.001 0.001 0.015 0.001 0.015 0.088 0.088 0.088 Species: N Atom Qatom Qorb 2s 2s 2py 2pz 2px 2py 2pz 2px 3dxy 3dyz 3dz2 3dxz 3dx2-y2 2 5.912 2.267 -0.434 1.523 1.523 1.523 -0.171 -0.171 -0.171 0.004 0.004 0.006 0.004 0.006 mulliken: Qtot = 15.000 cgvc: Finished line minimization 5. Mean atomic displacement = 0.0020 * Maximum dynamic memory allocated = 36 MB siesta: ============================== Begin CG move = 12 ============================== siesta: k-grid: Number of k-points = 1008 siesta: k-grid: Cutoff = 18.847 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 superc: Internal auxiliary supercell: 7 x 7 x 7 = 343 superc: Number of atoms, orbitals, and projectors: 686 9604 10976 outcell: Unit cell vectors (Ang): 2.221450 2.221452 0.000205 2.221399 0.000256 2.221394 0.000255 2.221402 2.221395 outcell: Cell vector modules (Ang) : 3.141606 3.141529 3.141531 outcell: Cell angles (23,13,12) (deg): 59.9925 59.9931 59.9931 outcell: Cell volume (Ang**3) : 21.9204 InitMesh: MESH = 30 x 30 x 30 = 27000 InitMesh: Mesh cutoff (required, used) = 250.000 378.007 Ry siesta: iscf Eharris(eV) E_KS(eV) FreeEng(eV) dDmax Ef(eV) siesta: 1 -974.7463 -974.7498 -974.7530 0.0023 -1.6731 siesta: 2 -974.7463 -974.7463 -974.7495 0.0013 -1.6726 siesta: 3 -974.7463 -974.7463 -974.7495 0.0002 -1.6728 siesta: 4 -974.7463 -974.7462 -974.7494 0.0001 -1.6728 siesta: 5 -974.7463 -974.7463 -974.7495 0.0000 -1.6729 siesta: E_KS(eV) = -974.7463 siesta: Atomic forces (eV/Ang): 1 -0.000375 -0.000331 -0.000108 2 0.004105 0.004067 0.004336 ---------------------------------------- Tot 0.003730 0.003736 0.004228 ---------------------------------------- Max 0.004336 Res 0.002957 sqrt( Sum f_i^2 / 3N ) ---------------------------------------- Max 0.004336 constrained Stress-tensor-Voigt (kbar): -1.99 -2.00 -1.89 -0.02 -0.03 -0.03 Target enthalpy (eV/cell) -974.7495 siesta: Stress tensor (static) (eV/Ang**3): -0.001245 -0.000014 -0.000021 -0.000014 -0.001247 -0.000019 -0.000021 -0.000019 -0.001181 siesta: Pressure (static): 1.96179475 kBar siesta: Stress tensor (total) (eV/Ang**3): -0.001245 -0.000014 -0.000021 -0.000014 -0.001247 -0.000019 -0.000021 -0.000019 -0.001181 siesta: Pressure (total): 1.96179475 kBar mulliken: Atomic and Orbital Populations: Species: Pt Atom Qatom Qorb 6s 6s 5dxy 5dyz 5dz2 5dxz 5dx2-y2 5dxy 5dyz 5dz2 5dxz 5dx2-y2 6Ppy 6Ppz 6Ppx 1 9.087 -0.577 0.984 1.953 1.953 1.263 1.953 1.263 0.001 0.001 0.015 0.001 0.015 0.088 0.088 0.088 Species: N Atom Qatom Qorb 2s 2s 2py 2pz 2px 2py 2pz 2px 3dxy 3dyz 3dz2 3dxz 3dx2-y2 2 5.913 2.267 -0.434 1.523 1.523 1.523 -0.171 -0.171 -0.171 0.004 0.004 0.006 0.004 0.006 mulliken: Qtot = 15.000 * Maximum dynamic memory allocated = 36 MB siesta: ============================== Begin CG move = 13 ============================== siesta: k-grid: Number of k-points = 1008 siesta: k-grid: Cutoff = 18.848 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 superc: Internal auxiliary supercell: 7 x 7 x 7 = 343 superc: Number of atoms, orbitals, and projectors: 686 9604 10976 outcell: Unit cell vectors (Ang): 2.221582 2.221584 0.000188 2.221528 0.000243 2.221525 0.000242 2.221531 2.221526 outcell: Cell vector modules (Ang) : 3.141793 3.141713 3.141716 outcell: Cell angles (23,13,12) (deg): 59.9928 59.9935 59.9936 outcell: Cell volume (Ang**3) : 21.9245 InitMesh: MESH = 30 x 30 x 30 = 27000 InitMesh: Mesh cutoff (required, used) = 250.000 377.964 Ry siesta: iscf Eharris(eV) E_KS(eV) FreeEng(eV) dDmax Ef(eV) siesta: 1 -974.7463 -974.7434 -974.7466 0.0019 -1.6740 siesta: 2 -974.7463 -974.7463 -974.7495 0.0011 -1.6744 siesta: 3 -974.7463 -974.7463 -974.7495 0.0002 -1.6742 siesta: 4 -974.7463 -974.7464 -974.7496 0.0001 -1.6742 siesta: 5 -974.7463 -974.7463 -974.7495 0.0000 -1.6742 siesta: E_KS(eV) = -974.7463 siesta: Atomic forces (eV/Ang): 1 0.004512 0.004672 0.004045 2 -0.000938 -0.000948 -0.000436 ---------------------------------------- Tot 0.003574 0.003723 0.003609 ---------------------------------------- Max 0.004672 Res 0.003176 sqrt( Sum f_i^2 / 3N ) ---------------------------------------- Max 0.004672 constrained Stress-tensor-Voigt (kbar): -1.50 -1.46 -1.48 0.09 0.01 -0.01 Target enthalpy (eV/cell) -974.7495 siesta: Stress tensor (static) (eV/Ang**3): -0.000933 0.000054 -0.000008 0.000054 -0.000914 0.000005 -0.000008 0.000005 -0.000925 siesta: Pressure (static): 1.48065367 kBar siesta: Stress tensor (total) (eV/Ang**3): -0.000933 0.000054 -0.000008 0.000054 -0.000914 0.000005 -0.000008 0.000005 -0.000925 siesta: Pressure (total): 1.48065367 kBar mulliken: Atomic and Orbital Populations: Species: Pt Atom Qatom Qorb 6s 6s 5dxy 5dyz 5dz2 5dxz 5dx2-y2 5dxy 5dyz 5dz2 5dxz 5dx2-y2 6Ppy 6Ppz 6Ppx 1 9.088 -0.577 0.984 1.953 1.953 1.263 1.953 1.263 0.001 0.001 0.015 0.001 0.015 0.088 0.088 0.088 Species: N Atom Qatom Qorb 2s 2s 2py 2pz 2px 2py 2pz 2px 3dxy 3dyz 3dz2 3dxz 3dx2-y2 2 5.912 2.267 -0.434 1.523 1.523 1.523 -0.171 -0.171 -0.171 0.004 0.004 0.006 0.004 0.006 mulliken: Qtot = 15.000 cgvc: Finished line minimization 6. Mean atomic displacement = 0.0003 * Maximum dynamic memory allocated = 36 MB siesta: ============================== Begin CG move = 14 ============================== siesta: k-grid: Number of k-points = 1008 siesta: k-grid: Cutoff = 18.849 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 superc: Internal auxiliary supercell: 7 x 7 x 7 = 343 superc: Number of atoms, orbitals, and projectors: 686 9604 10976 outcell: Unit cell vectors (Ang): 2.221626 2.221627 0.000188 2.221574 0.000240 2.221571 0.000239 2.221576 2.221571 outcell: Cell vector modules (Ang) : 3.141855 3.141778 3.141780 outcell: Cell angles (23,13,12) (deg): 59.9929 59.9936 59.9936 outcell: Cell volume (Ang**3) : 21.9259 InitMesh: MESH = 30 x 30 x 30 = 27000 InitMesh: Mesh cutoff (required, used) = 250.000 377.950 Ry siesta: iscf Eharris(eV) E_KS(eV) FreeEng(eV) dDmax Ef(eV) siesta: 1 -974.7463 -974.7454 -974.7486 0.0003 -1.6746 siesta: 2 -974.7463 -974.7463 -974.7495 0.0004 -1.6747 siesta: 3 -974.7463 -974.7462 -974.7494 0.0001 -1.6746 siesta: E_KS(eV) = -974.7463 siesta: Atomic forces (eV/Ang): 1 0.003750 0.003729 0.004111 2 -0.000278 -0.000304 0.000155 ---------------------------------------- Tot 0.003473 0.003425 0.004267 ---------------------------------------- Max 0.004111 Res 0.002741 sqrt( Sum f_i^2 / 3N ) ---------------------------------------- Max 0.004111 constrained Stress-tensor-Voigt (kbar): -1.38 -1.38 -1.31 -0.03 -0.09 -0.09 Target enthalpy (eV/cell) -974.7495 siesta: Stress tensor (static) (eV/Ang**3): -0.000863 -0.000016 -0.000056 -0.000016 -0.000862 -0.000057 -0.000056 -0.000057 -0.000816 siesta: Pressure (static): 1.35756227 kBar siesta: Stress tensor (total) (eV/Ang**3): -0.000863 -0.000016 -0.000056 -0.000016 -0.000862 -0.000057 -0.000056 -0.000057 -0.000816 siesta: Pressure (total): 1.35756227 kBar mulliken: Atomic and Orbital Populations: Species: Pt Atom Qatom Qorb 6s 6s 5dxy 5dyz 5dz2 5dxz 5dx2-y2 5dxy 5dyz 5dz2 5dxz 5dx2-y2 6Ppy 6Ppz 6Ppx 1 9.088 -0.577 0.984 1.953 1.953 1.263 1.953 1.263 0.001 0.001 0.015 0.001 0.015 0.088 0.088 0.088 Species: N Atom Qatom Qorb 2s 2s 2py 2pz 2px 2py 2pz 2px 3dxy 3dyz 3dz2 3dxz 3dx2-y2 2 5.912 2.267 -0.434 1.523 1.523 1.523 -0.171 -0.171 -0.171 0.004 0.004 0.006 0.004 0.006 mulliken: Qtot = 15.000 * Maximum dynamic memory allocated = 36 MB siesta: ============================== Begin CG move = 15 ============================== siesta: k-grid: Number of k-points = 1008 siesta: k-grid: Cutoff = 18.849 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 superc: Internal auxiliary supercell: 7 x 7 x 7 = 343 superc: Number of atoms, orbitals, and projectors: 686 9604 10976 outcell: Unit cell vectors (Ang): 2.221696 2.221696 0.000186 2.221648 0.000235 2.221645 0.000235 2.221648 2.221644 outcell: Cell vector modules (Ang) : 3.141953 3.141882 3.141882 outcell: Cell angles (23,13,12) (deg): 59.9931 59.9937 59.9937 outcell: Cell volume (Ang**3) : 21.9281 InitMesh: MESH = 30 x 30 x 30 = 27000 InitMesh: Mesh cutoff (required, used) = 250.000 377.926 Ry siesta: iscf Eharris(eV) E_KS(eV) FreeEng(eV) dDmax Ef(eV) siesta: 1 -974.7463 -974.7448 -974.7480 0.0005 -1.6752 siesta: 2 -974.7463 -974.7463 -974.7495 0.0006 -1.6754 siesta: 3 -974.7463 -974.7462 -974.7494 0.0001 -1.6753 siesta: 4 -974.7463 -974.7464 -974.7495 0.0000 -1.6753 siesta: E_KS(eV) = -974.7463 siesta: Atomic forces (eV/Ang): 1 0.004187 0.004175 -0.003995 2 0.000847 0.000853 0.001196 ---------------------------------------- Tot 0.005034 0.005029 -0.002798 ---------------------------------------- Max 0.004187 Res 0.002994 sqrt( Sum f_i^2 / 3N ) ---------------------------------------- Max 0.004187 constrained Stress-tensor-Voigt (kbar): -1.21 -1.21 0.56 -0.07 -0.04 -0.04 Target enthalpy (eV/cell) -974.7495 siesta: Stress tensor (static) (eV/Ang**3): -0.000756 -0.000044 -0.000026 -0.000044 -0.000756 -0.000026 -0.000026 -0.000026 0.000352 siesta: Pressure (static): 0.61939932 kBar siesta: Stress tensor (total) (eV/Ang**3): -0.000756 -0.000044 -0.000026 -0.000044 -0.000756 -0.000026 -0.000026 -0.000026 0.000352 siesta: Pressure (total): 0.61939932 kBar mulliken: Atomic and Orbital Populations: Species: Pt Atom Qatom Qorb 6s 6s 5dxy 5dyz 5dz2 5dxz 5dx2-y2 5dxy 5dyz 5dz2 5dxz 5dx2-y2 6Ppy 6Ppz 6Ppx 1 9.088 -0.576 0.984 1.953 1.953 1.263 1.953 1.263 0.001 0.001 0.015 0.001 0.015 0.088 0.088 0.088 Species: N Atom Qatom Qorb 2s 2s 2py 2pz 2px 2py 2pz 2px 3dxy 3dyz 3dz2 3dxz 3dx2-y2 2 5.912 2.267 -0.434 1.523 1.523 1.523 -0.171 -0.171 -0.171 0.004 0.004 0.006 0.004 0.006 mulliken: Qtot = 15.000 * Maximum dynamic memory allocated = 36 MB siesta: ============================== Begin CG move = 16 ============================== siesta: k-grid: Number of k-points = 1008 siesta: k-grid: Cutoff = 18.850 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 superc: Internal auxiliary supercell: 7 x 7 x 7 = 343 superc: Number of atoms, orbitals, and projectors: 686 9604 10976 outcell: Unit cell vectors (Ang): 2.221809 2.221806 0.000185 2.221767 0.000227 2.221762 0.000228 2.221763 2.221761 outcell: Cell vector modules (Ang) : 3.142110 3.142050 3.142045 outcell: Cell angles (23,13,12) (deg): 59.9933 59.9939 59.9938 outcell: Cell volume (Ang**3) : 21.9316 InitMesh: MESH = 30 x 30 x 30 = 27000 InitMesh: Mesh cutoff (required, used) = 250.000 377.887 Ry siesta: iscf Eharris(eV) E_KS(eV) FreeEng(eV) dDmax Ef(eV) siesta: 1 -974.7463 -974.7439 -974.7471 0.0008 -1.6763 siesta: 2 -974.7464 -974.7464 -974.7496 0.0009 -1.6766 siesta: 3 -974.7464 -974.7461 -974.7493 0.0002 -1.6765 siesta: 4 -974.7464 -974.7464 -974.7496 0.0000 -1.6765 siesta: E_KS(eV) = -974.7464 siesta: Atomic forces (eV/Ang): 1 -0.004719 -0.004723 -0.003994 2 0.002587 0.002589 0.002849 ---------------------------------------- Tot -0.002131 -0.002134 -0.001145 ---------------------------------------- Max 0.004723 Res 0.003698 sqrt( Sum f_i^2 / 3N ) ---------------------------------------- Max 0.004723 constrained Stress-tensor-Voigt (kbar): 1.12 1.12 1.16 0.02 0.03 0.03 Target enthalpy (eV/cell) -974.7496 siesta: Stress tensor (static) (eV/Ang**3): 0.000697 0.000011 0.000021 0.000011 0.000698 0.000021 0.000021 0.000021 0.000723 siesta: Pressure (static): -1.13056893 kBar siesta: Stress tensor (total) (eV/Ang**3): 0.000697 0.000011 0.000021 0.000011 0.000698 0.000021 0.000021 0.000021 0.000723 siesta: Pressure (total): -1.13056893 kBar mulliken: Atomic and Orbital Populations: Species: Pt Atom Qatom Qorb 6s 6s 5dxy 5dyz 5dz2 5dxz 5dx2-y2 5dxy 5dyz 5dz2 5dxz 5dx2-y2 6Ppy 6Ppz 6Ppx 1 9.088 -0.576 0.983 1.953 1.953 1.263 1.953 1.263 0.001 0.001 0.015 0.001 0.015 0.088 0.088 0.088 Species: N Atom Qatom Qorb 2s 2s 2py 2pz 2px 2py 2pz 2px 3dxy 3dyz 3dz2 3dxz 3dx2-y2 2 5.912 2.267 -0.434 1.523 1.523 1.523 -0.170 -0.170 -0.170 0.004 0.004 0.006 0.004 0.006 mulliken: Qtot = 15.000 * Maximum dynamic memory allocated = 36 MB siesta: ============================== Begin CG move = 17 ============================== siesta: k-grid: Number of k-points = 1008 siesta: k-grid: Cutoff = 18.850 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 superc: Internal auxiliary supercell: 7 x 7 x 7 = 343 superc: Number of atoms, orbitals, and projectors: 686 9604 10976 outcell: Unit cell vectors (Ang): 2.221731 2.221730 0.000186 2.221685 0.000232 2.221681 0.000233 2.221683 2.221680 outcell: Cell vector modules (Ang) : 3.142001 3.141934 3.141932 outcell: Cell angles (23,13,12) (deg): 59.9931 59.9937 59.9937 outcell: Cell volume (Ang**3) : 21.9291 InitMesh: MESH = 30 x 30 x 30 = 27000 InitMesh: Mesh cutoff (required, used) = 250.000 377.914 Ry siesta: iscf Eharris(eV) E_KS(eV) FreeEng(eV) dDmax Ef(eV) siesta: 1 -974.7464 -974.7481 -974.7512 0.0006 -1.6758 siesta: 2 -974.7464 -974.7464 -974.7495 0.0006 -1.6755 siesta: 3 -974.7464 -974.7465 -974.7497 0.0001 -1.6757 siesta: 4 -974.7464 -974.7463 -974.7495 0.0000 -1.6756 siesta: E_KS(eV) = -974.7464 siesta: Atomic forces (eV/Ang): 1 -0.005633 -0.005642 -0.004869 2 0.001392 0.001391 0.001659 ---------------------------------------- Tot -0.004241 -0.004251 -0.003209 ---------------------------------------- Max 0.005642 Res 0.003956 sqrt( Sum f_i^2 / 3N ) ---------------------------------------- Max 0.005642 constrained Stress-tensor-Voigt (kbar): 0.79 0.79 0.83 -0.04 -0.01 -0.01 Target enthalpy (eV/cell) -974.7495 siesta: Stress tensor (static) (eV/Ang**3): 0.000493 -0.000023 -0.000007 -0.000023 0.000494 -0.000007 -0.000007 -0.000007 0.000517 siesta: Pressure (static): -0.80316114 kBar siesta: Stress tensor (total) (eV/Ang**3): 0.000493 -0.000023 -0.000007 -0.000023 0.000494 -0.000007 -0.000007 -0.000007 0.000517 siesta: Pressure (total): -0.80316114 kBar mulliken: Atomic and Orbital Populations: Species: Pt Atom Qatom Qorb 6s 6s 5dxy 5dyz 5dz2 5dxz 5dx2-y2 5dxy 5dyz 5dz2 5dxz 5dx2-y2 6Ppy 6Ppz 6Ppx 1 9.088 -0.576 0.984 1.953 1.953 1.263 1.953 1.263 0.001 0.001 0.015 0.001 0.015 0.088 0.088 0.088 Species: N Atom Qatom Qorb 2s 2s 2py 2pz 2px 2py 2pz 2px 3dxy 3dyz 3dz2 3dxz 3dx2-y2 2 5.912 2.267 -0.434 1.523 1.523 1.523 -0.171 -0.171 -0.171 0.004 0.004 0.006 0.004 0.006 mulliken: Qtot = 15.000 * Maximum dynamic memory allocated = 36 MB outcoor: Relaxed atomic coordinates (fractional): -0.00015043 -0.00009060 -0.00009192 1 Pt 1 0.49979152 0.49982518 0.49982538 2 N 2 outcell: Unit cell vectors (Ang): 2.221731 2.221730 0.000186 2.221685 0.000232 2.221681 0.000233 2.221683 2.221680 outcell: Cell vector modules (Ang) : 3.142001 3.141934 3.141932 outcell: Cell angles (23,13,12) (deg): 59.9931 59.9937 59.9937 outcell: Cell volume (Ang**3) : 21.9291 coxmol: Writing XMOL coordinates into file ptnRS.xyz siesta: Program's energy decomposition (eV): siesta: Eions = 1401.175675 siesta: Ena = 241.332615 siesta: Ekin = 691.485265 siesta: Enl = -317.830910 siesta: DEna = -4.953655 siesta: DUscf = 4.165897 siesta: DUext = 0.000000 siesta: Exc = -187.769894 siesta: eta*DQ = 0.000000 siesta: Emadel = 0.000000 siesta: Ekinion = 0.000000 siesta: Eharris = -974.746366 siesta: Etot = -974.746358 siesta: FreeEng = -974.749547 siesta: Final energy (eV): siesta: Kinetic = 691.485265 siesta: Hartree = 82.317724 siesta: Ext. field = 0.000000 siesta: Exch.-corr. = -187.769894 siesta: Ion-electron = -719.815005 siesta: Ion-ion = -840.964447 siesta: Ekinion = 0.000000 siesta: Total = -974.746358 siesta: Stress tensor (static) (eV/Ang**3): siesta: 0.000493 -0.000023 -0.000007 siesta: -0.000023 0.000494 -0.000007 siesta: -0.000007 -0.000007 0.000517 siesta: Constrained stress tensor (static) (eV/Ang**3): siesta: 0.000493 -0.000023 -0.000007 siesta: -0.000023 0.000494 -0.000007 siesta: -0.000007 -0.000007 0.000517 siesta: Cell volume = 21.929137 Ang**3 siesta: Pressure (static): siesta: Solid Molecule Units siesta: -0.00000546 -0.00000925 Ry/Bohr**3 siesta: -0.00050129 -0.00084885 eV/Ang**3 siesta: -0.80316114 -1.36002229 kBar * Maximum dynamic memory allocated : Node 0 = 36 MB * Maximum dynamic memory allocated : Node 1 = 35 MB * Maximum memory occured during ranger timer: CPU execution times: timer: Routine Calls Time/call Tot.time % timer: siesta 1 5494.790 5494.790 100.00 timer: Setup 1 2.630 2.630 0.05 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 18 305.072 5491.290 99.94 timer: hsparse 22 2.139 47.060 0.86 timer: overfsm 36 0.158 5.691 0.10 timer: IterSCF 133 40.670 5409.089 98.44 timer: kinefsm 36 0.155 5.591 0.10 timer: nlefsm 36 10.420 375.131 6.83 timer: DHSCF 133 5.329 708.719 12.90 timer: DHSCF1 18 0.809 14.559 0.26 timer: DHSCF2 18 14.373 258.721 4.71 timer: REORD 834 0.000 0.090 0.00 timer: POISON 151 0.009 1.329 0.02 timer: DHSCF3 133 2.787 370.689 6.75 timer: rhoofd 133 1.244 165.440 3.01 timer: cellXC 133 0.016 2.131 0.04 timer: vmat 133 1.450 192.878 3.51 timer: diagon 115 36.170 4159.502 75.70 timer: cdiag 231840 0.012 2752.704 50.10 timer: cdiag1 231840 0.001 154.803 2.82 timer: cdiag2 231840 0.001 216.046 3.93 timer: cdiag3 231840 0.009 2132.258 38.81 timer: cdiag4 231840 0.000 62.868 1.14 timer: DHSCF4 18 3.587 64.570 1.18 timer: dfscf 18 2.737 49.260 0.90 timer: optical 1 0.000 0.000 0.00 elaps: ELAPSED times: elaps: Routine Calls Time/call Tot.time % elaps: siesta 1 2608.570 2608.570 100.00 elaps: Setup 1 2.220 2.220 0.09 elaps: bands 1 0.000 0.000 0.00 elaps: writewave 1 0.010 0.010 0.00 elaps: KSV_init 1 0.000 0.000 0.00 elaps: IterMD 18 144.772 2605.890 99.90 elaps: hsparse 22 1.091 24.000 0.92 elaps: overfsm 36 0.087 3.130 0.12 elaps: IterSCF 133 19.283 2564.690 98.32 elaps: kinefsm 36 0.083 2.990 0.11 elaps: nlefsm 36 4.182 150.540 5.77 elaps: DHSCF 133 2.671 355.260 13.62 elaps: DHSCF1 18 0.406 7.300 0.28 elaps: DHSCF2 18 7.183 129.290 4.96 elaps: REORD 834 0.000 0.060 0.00 elaps: POISON 151 0.004 0.610 0.02 elaps: DHSCF3 133 1.391 184.980 7.09 elaps: rhoofd 133 0.652 86.670 3.32 elaps: cellXC 133 0.009 1.150 0.04 elaps: vmat 133 0.724 96.260 3.69 elaps: diagon 115 16.521 1899.880 72.83 elaps: cdiag 231840 0.005 1206.230 46.24 elaps: cdiag1 231840 0.000 69.680 2.67 elaps: cdiag2 231840 0.000 80.290 3.08 elaps: cdiag3 231840 0.004 934.660 35.83 elaps: cdiag4 231840 0.000 26.350 1.01 elaps: DHSCF4 18 1.869 33.650 1.29 elaps: dfscf 18 1.441 25.930 0.99 elaps: optical 1 0.000 0.000 0.00 >> End of run: 22-MAR-2007 13:26:09 PS: Read file <err.file> for stderr output of this job.