=------------------------------------------------------------------------------= CP: variable-cell Car-Parrinello molecular dynamics using norm-conserving and ultrasoft Vanderbilt pseudopotentials Version: 4.0 - Mon Apr 28 15:32:33 CEST 2008 Authors: Alfredo Pasquarello, Kari Laasonen, Andrea Trave, Roberto Car, Paolo Giannozzi, Nicola Marzari, Carlo Cavazzoni, Guido Chiarotti, Sandro Scandolo, Paolo Focher, Gerardo Ballabio, and others =------------------------------------------------------------------------------= This run was started on: 16:30:16 28Apr2008 Serial Build Warning: card &CELL ignored Warning: card CELL_DYNAMICS = 'NONE', ignored Warning: card / ignored Job Title: MD Simulation Atomic Pseudopotentials Parameters ---------------------------------- Reading pseudopotential for specie # 1 from file : /home/giannozz/espresso/pseudo/O.pz-rrkjus.UPF file type is 20: UPF Reading pseudopotential for specie # 2 from file : /home/giannozz/espresso/pseudo/Mg.pz-n-vbc.UPF file type is 20: UPF Main Simulation Parameters (from input) --------------------------------------- Restart Mode = 1 restart Number of MD Steps = 1 Print out every 1 MD Steps Reads from unit = 50 Writes to unit = 50 MD Simulation time step = 1.00 Electronic fictitious mass (emass) = 1000.00 emass cut-off = 4.00 Simulation Cell Parameters (from input) external pressure = 0.00 [GPa] wmass (calculated) = 160952.54 [AU] ibrav = 1 alat = 15.91200000 a1 = 15.91200000 0.00000000 0.00000000 a2 = 0.00000000 15.91200000 0.00000000 a3 = 0.00000000 0.00000000 15.91200000 b1 = 0.06284565 0.00000000 0.00000000 b2 = 0.00000000 0.06284565 0.00000000 b3 = 0.00000000 0.00000000 0.06284565 omega = 4028.78703053 ======================================== | CONJUGATE GRADIENT | ======================================== | iterations = 250 | | conv_thr = 0.00000001000 a.u. | | passop = 0.30000 a.u. | | niter_cg_restart = 20 | ======================================== Energy Cut-offs --------------- Ecutwfc = 20.0 Ry, Ecutrho = 160.0 Ry, Ecuts = 80.0 Ry Gcutwfc = 11.3 , Gcutrho = 32.0 Gcuts = 22.7 NOTA BENE: refg, mmx = 0.050000 3840 Orthog. with Gram-Schmidt Electron dynamics with newton equations Electron dynamics : the temperature is not controlled initial random displacement of el. coordinates with amplitude= 0.020000 Electronic states ----------------- Number of Electron = 256, of States = 128 Occupation numbers : 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 Exchange and correlations functionals ------------------------------------- Using Local Density Approximation with Exchange functional: SLATER Correlation functional: PERDEW AND ZUNGER Exchange-correlation = SLA PZ NOGX NOGC (1100) Ions Simulation Parameters -------------------------- Ions are not allowed to move Ionic position (from input) sorted by specie, and converted to real a.u. coordinates Species 1 atoms = 32 mass = 21892.89 (a.u.), 12.01 (amu) rcmax = 0.80 (a.u.) 0.066300 0.066300 4.044300 4.044300 4.044300 4.044300 4.044300 0.066300 8.022300 0.066300 4.044300 8.022300 0.066300 0.066300 12.000300 4.044300 4.044300 12.000300 4.044300 0.066300 15.978300 0.066300 4.044300 15.978300 0.066300 8.022300 4.044300 4.044300 12.000300 4.044300 4.044300 8.022300 8.022300 0.066300 12.000300 8.022300 0.066300 8.022300 12.000300 4.044300 12.000300 12.000300 4.044300 8.022300 15.978300 0.066300 12.000300 15.978300 8.022300 0.066300 4.044300 12.000300 4.044300 4.044300 12.000300 0.066300 8.022300 8.022300 4.044300 8.022300 8.022300 0.066300 12.000300 12.000300 4.044300 12.000300 12.000300 0.066300 15.978300 8.022300 4.044300 15.978300 8.022300 8.022300 4.044300 12.000300 12.000300 4.044300 12.000300 8.022300 8.022300 8.022300 12.000300 8.022300 8.022300 8.022300 12.000300 12.000300 12.000300 12.000300 12.000300 8.022300 15.978300 8.022300 12.000300 15.978300 Species 2 atoms = 32 mass = 44296.19 (a.u.), 24.30 (amu) rcmax = 0.50 (a.u.) 0.066300 0.066300 0.066300 4.044300 4.044300 0.066300 4.044300 0.066300 4.044300 0.066300 4.044300 4.044300 0.066300 0.066300 8.022300 4.044300 4.044300 8.022300 4.044300 0.066300 12.000300 0.066300 4.044300 12.000300 0.066300 8.022300 0.066300 4.044300 12.000300 0.066300 4.044300 8.022300 4.044300 0.066300 12.000300 4.044300 0.066300 8.022300 8.022300 4.044300 12.000300 8.022300 4.044300 8.022300 12.000300 0.066300 12.000300 12.000300 8.022300 0.066300 0.066300 12.000300 4.044300 0.066300 12.000300 0.066300 4.044300 8.022300 4.044300 4.044300 8.022300 0.066300 8.022300 12.000300 4.044300 8.022300 12.000300 0.066300 12.000300 8.022300 4.044300 12.000300 8.022300 8.022300 0.066300 12.000300 12.000300 0.066300 12.000300 8.022300 4.044300 8.022300 12.000300 4.044300 8.022300 8.022300 8.022300 12.000300 12.000300 8.022300 12.000300 8.022300 12.000300 8.022300 12.000300 12.000300 Ionic position read from input file Cell Dynamics Parameters (from STDIN) ------------------------------------- Starting cell generated from CELLDM Cell parameters will be re-read from restart file Constant VOLUME Molecular dynamics cell parameters are not allowed to move ===================================== | BERRY PHASE ELECTRIC FIELD 1 ===================================== | direction = 3 | intensity = 0.00100 a.u. ===================================== Verbosity: iprsta = 1 Simulation dimensions initialization ------------------------------------ unit vectors of full simulation cell in real space: in reciprocal space (units 2pi/alat): 1 15.9120 0.0000 0.0000 1.0000 0.0000 0.0000 2 0.0000 15.9120 0.0000 0.0000 1.0000 0.0000 3 0.0000 0.0000 15.9120 0.0000 0.0000 1.0000 Stick Mesh ---------- nst = 1617, nstw = 203, nsts = 805 PEs n.st n.stw n.sts n.g n.gw n.gs 1 3233 405 1609 137929 6043 48693 0 3233 405 1609 137929 6043 48693 Real Mesh --------- Global Dimensions Local Dimensions Processor Grid .X. .Y. .Z. .X. .Y. .Z. .X. .Y. .Z. 72 72 72 72 72 72 1 1 1 Array leading dimensions ( nr1x, nr2x, nr3x ) = 72 72 72 Local number of cell to store the grid ( nnrx ) = 373248 Number of x-y planes for each processors: nr3l = 72 Smooth Real Mesh ---------------- Global Dimensions Local Dimensions Processor Grid .X. .Y. .Z. .X. .Y. .Z. .X. .Y. .Z. 45 45 45 45 45 45 1 1 1 Array leading dimensions ( nr1x, nr2x, nr3x ) = 45 45 45 Local number of cell to store the grid ( nnrx ) = 91125 Number of x-y planes for each processors: nr3sl = 45 Small Box Real Mesh ------------------- Global Dimensions Local Dimensions Processor Grid .X. .Y. .Z. .X. .Y. .Z. .X. .Y. .Z. 30 30 30 30 30 30 1 1 1 Array leading dimensions ( nr1x, nr2x, nr3x ) = 30 30 30 Local number of cell to store the grid ( nnrx ) = 27000 unit vectors of box grid cell in real space: in reciprocal space: 6.6300 0.0000 0.0000 1.0000 0.0000 0.0000 0.0000 6.6300 0.0000 0.0000 1.0000 0.0000 0.0000 0.0000 6.6300 0.0000 0.0000 1.0000 Reciprocal Space Mesh --------------------- Large Mesh PE Global(ngmt) Local(ngm) MaxLocal(ngmx) 1 68965 68965 68965 Smooth Mesh PE Global(ngst) Local(ngs) MaxLocal(ngsx) 1 24347 24347 24347 Wave function Mesh PE Global(ngwt) Local(ngw) MaxLocal(ngwx) 1 3022 3022 3022 Small box Mesh ngb = 4970 not distributed to processors System geometry initialization ------------------------------ Scaled positions from standard input O 0.416667E-02 0.416667E-02 0.254167E+00 O 0.254167E+00 0.254167E+00 0.254167E+00 O 0.254167E+00 0.416667E-02 0.504167E+00 O 0.416667E-02 0.254167E+00 0.504167E+00 O 0.416667E-02 0.416667E-02 0.754167E+00 O 0.254167E+00 0.254167E+00 0.754167E+00 O 0.254167E+00 0.416667E-02 0.100417E+01 O 0.416667E-02 0.254167E+00 0.100417E+01 O 0.416667E-02 0.504167E+00 0.254167E+00 O 0.254167E+00 0.754167E+00 0.254167E+00 O 0.254167E+00 0.504167E+00 0.504167E+00 O 0.416667E-02 0.754167E+00 0.504167E+00 O 0.416667E-02 0.504167E+00 0.754167E+00 O 0.254167E+00 0.754167E+00 0.754167E+00 O 0.254167E+00 0.504167E+00 0.100417E+01 O 0.416667E-02 0.754167E+00 0.100417E+01 O 0.504167E+00 0.416667E-02 0.254167E+00 O 0.754167E+00 0.254167E+00 0.254167E+00 O 0.754167E+00 0.416667E-02 0.504167E+00 O 0.504167E+00 0.254167E+00 0.504167E+00 O 0.504167E+00 0.416667E-02 0.754167E+00 O 0.754167E+00 0.254167E+00 0.754167E+00 O 0.754167E+00 0.416667E-02 0.100417E+01 O 0.504167E+00 0.254167E+00 0.100417E+01 O 0.504167E+00 0.504167E+00 0.254167E+00 O 0.754167E+00 0.754167E+00 0.254167E+00 O 0.754167E+00 0.504167E+00 0.504167E+00 O 0.504167E+00 0.754167E+00 0.504167E+00 O 0.504167E+00 0.504167E+00 0.754167E+00 O 0.754167E+00 0.754167E+00 0.754167E+00 O 0.754167E+00 0.504167E+00 0.100417E+01 O 0.504167E+00 0.754167E+00 0.100417E+01 Mg 0.416667E-02 0.416667E-02 0.416667E-02 Mg 0.254167E+00 0.254167E+00 0.416667E-02 Mg 0.254167E+00 0.416667E-02 0.254167E+00 Mg 0.416667E-02 0.254167E+00 0.254167E+00 Mg 0.416667E-02 0.416667E-02 0.504167E+00 Mg 0.254167E+00 0.254167E+00 0.504167E+00 Mg 0.254167E+00 0.416667E-02 0.754167E+00 Mg 0.416667E-02 0.254167E+00 0.754167E+00 Mg 0.416667E-02 0.504167E+00 0.416667E-02 Mg 0.254167E+00 0.754167E+00 0.416667E-02 Mg 0.254167E+00 0.504167E+00 0.254167E+00 Mg 0.416667E-02 0.754167E+00 0.254167E+00 Mg 0.416667E-02 0.504167E+00 0.504167E+00 Mg 0.254167E+00 0.754167E+00 0.504167E+00 Mg 0.254167E+00 0.504167E+00 0.754167E+00 Mg 0.416667E-02 0.754167E+00 0.754167E+00 Mg 0.504167E+00 0.416667E-02 0.416667E-02 Mg 0.754167E+00 0.254167E+00 0.416667E-02 Mg 0.754167E+00 0.416667E-02 0.254167E+00 Mg 0.504167E+00 0.254167E+00 0.254167E+00 Mg 0.504167E+00 0.416667E-02 0.504167E+00 Mg 0.754167E+00 0.254167E+00 0.504167E+00 Mg 0.754167E+00 0.416667E-02 0.754167E+00 Mg 0.504167E+00 0.254167E+00 0.754167E+00 Mg 0.504167E+00 0.504167E+00 0.416667E-02 Mg 0.754167E+00 0.754167E+00 0.416667E-02 Mg 0.754167E+00 0.504167E+00 0.254167E+00 Mg 0.504167E+00 0.754167E+00 0.254167E+00 Mg 0.504167E+00 0.504167E+00 0.504167E+00 Mg 0.754167E+00 0.754167E+00 0.504167E+00 Mg 0.754167E+00 0.504167E+00 0.754167E+00 Mg 0.504167E+00 0.754167E+00 0.754167E+00 ibrav = 1 cell parameters 15.91200 0.00000 0.00000 0.00000 15.91200 0.00000 0.00000 0.00000 15.91200 Pseudopotentials initialization ------------------------------- nlinit nh(is), ngb, is, kkbeta, lmaxq = 8 4970 1 865 3 qqq -0.0987 0.4865 0.0000 0.0000 0.4865 -2.1787 0.0000 0.0000 0.0000 0.0000 0.2330 0.2950 0.0000 0.0000 0.2950 0.3737 Common initialization Specie: 1 1 indv= 1 ang. mom= 0 2 indv= 2 ang. mom= 0 3 indv= 3 ang. mom= 1 4 indv= 3 ang. mom= 1 5 indv= 3 ang. mom= 1 6 indv= 4 ang. mom= 1 7 indv= 4 ang. mom= 1 8 indv= 4 ang. mom= 1 dion 0.4817 -1.2813 0.0000 0.0000 -1.2813 2.3075 0.0000 0.0000 0.0000 0.0000 0.6338 0.8752 0.0000 0.0000 0.8752 1.2039 Specie: 2 1 indv= 1 ang. mom= 0 2 indv= 2 ang. mom= 1 3 indv= 2 ang. mom= 1 4 indv= 2 ang. mom= 1 dion 3.2181 0.0000 0.0000 12.0575 Short Legend and Physical Units in the Output --------------------------------------------- NFI [int] - step index EKINC [HARTREE A.U.] - kinetic energy of the fictitious electronic dynamics TEMPH [K] - Temperature of the fictitious cell dynamics TEMP [K] - Ionic temperature ETOT [HARTREE A.U.] - Scf total energy (Kohn-Sham hamiltonian) ENTHAL [HARTREE A.U.] - Enthalpy ( ETOT + P * V ) ECONS [HARTREE A.U.] - Enthalpy + kinetic energy of ions and cell ECONT [HARTREE A.U.] - Constant of motion for the CP lagrangian reading restart file: /mgo_50.save restart file read in 0.091 sec. nprint_nfi= -2 nprint_nfi= 1 formf: eself= 676.60611 formf: vps(g=0)= -0.0020774 rhops(g=0)= -0.0014893 formf: sum_g vps(g)= -1.4455266 sum_g rhops(g)= -1.0530008 formf: vps(g=0)= 0.0012432 rhops(g=0)= -0.0004964 formf: sum_g vps(g)= 1.6241241 sum_g rhops(g)= -1.4102948 Delta V(G=0): 0.108296Ry, 2.946897eV Initial random displacement of el. coordinates Amplitude = 0.020000 Initialize Berry phase electric field #not found, gtable: 405.000000000000 #not found, gtabin: 405.000000000000 * Physical Quantities at step: 2 PERFORMING CONJUGATE GRADIENT MINIMIZATION OF EL. STATES total energy = -547.28628 Hartree a.u. kinetic energy = 235.43095 Hartree a.u. electrostatic energy = -588.81023 Hartree a.u. esr = 0.00000 Hartree a.u. eself = 676.60611 Hartree a.u. pseudopotential energy = -167.87987 Hartree a.u. n-l pseudopotential energy = 108.84068 Hartree a.u. exchange-correlation energy = -134.85134 Hartree a.u. average potential = 0.00000 Hartree a.u. Eigenvalues (eV), kp = 1 , spin = 1 -14.79 -14.01 -14.01 -14.01 -14.01 -14.01 -14.01 -13.53 -13.53 -13.53 -13.53 -13.53 -13.53 -13.53 -13.53 -13.53 -13.53 -13.53 -13.53 -13.47 -13.47 -13.47 -13.47 -13.14 -13.14 -13.14 -13.11 -13.11 -13.11 -13.11 -13.11 -13.11 -2.27 -2.27 -2.27 -2.27 -1.83 -1.83 -1.83 -1.15 -1.15 -1.15 -1.15 -1.15 -1.15 -1.15 -1.15 -1.15 -1.15 -1.15 -1.15 -0.69 -0.69 -0.69 -0.69 -0.69 -0.69 -0.69 -0.69 -0.69 -0.69 -0.69 -0.69 -0.21 -0.21 -0.21 -0.21 -0.21 -0.21 0.33 0.33 0.33 0.33 0.33 0.33 0.33 0.33 0.33 0.33 0.33 0.33 0.63 0.63 0.63 0.63 0.63 0.63 1.02 1.02 1.02 1.02 1.02 1.02 1.75 1.75 1.75 1.75 1.75 1.75 1.75 1.75 1.75 1.75 1.75 1.75 1.87 1.87 1.87 1.87 1.87 1.87 1.87 1.87 1.87 1.87 1.87 1.87 1.92 1.92 1.92 1.92 1.92 1.92 1.92 1.92 2.46 2.46 2.46 Allocated memory (kb) = 119232 CELL_PARAMETERS 15.91200000 0.00000000 0.00000000 0.00000000 15.91200000 0.00000000 0.00000000 0.00000000 15.91200000 System Density [g/cm^3] : 3.2319 Center of mass square displacement (a.u.): 0.000000 ATOMIC_POSITIONS Mg 0.663000E-01 0.663000E-01 0.663000E-01 Mg 0.404430E+01 0.404430E+01 0.663000E-01 Mg 0.404430E+01 0.663000E-01 0.404430E+01 Mg 0.663000E-01 0.404430E+01 0.404430E+01 Mg 0.663000E-01 0.663000E-01 0.802230E+01 Mg 0.404430E+01 0.404430E+01 0.802230E+01 Mg 0.404430E+01 0.663000E-01 0.120003E+02 Mg 0.663000E-01 0.404430E+01 0.120003E+02 Mg 0.663000E-01 0.802230E+01 0.663000E-01 Mg 0.404430E+01 0.120003E+02 0.663000E-01 Mg 0.404430E+01 0.802230E+01 0.404430E+01 Mg 0.663000E-01 0.120003E+02 0.404430E+01 Mg 0.663000E-01 0.802230E+01 0.802230E+01 Mg 0.404430E+01 0.120003E+02 0.802230E+01 Mg 0.404430E+01 0.802230E+01 0.120003E+02 Mg 0.663000E-01 0.120003E+02 0.120003E+02 Mg 0.802230E+01 0.663000E-01 0.663000E-01 Mg 0.120003E+02 0.404430E+01 0.663000E-01 Mg 0.120003E+02 0.663000E-01 0.404430E+01 Mg 0.802230E+01 0.404430E+01 0.404430E+01 Mg 0.802230E+01 0.663000E-01 0.802230E+01 Mg 0.120003E+02 0.404430E+01 0.802230E+01 Mg 0.120003E+02 0.663000E-01 0.120003E+02 Mg 0.802230E+01 0.404430E+01 0.120003E+02 Mg 0.802230E+01 0.802230E+01 0.663000E-01 Mg 0.120003E+02 0.120003E+02 0.663000E-01 Mg 0.120003E+02 0.802230E+01 0.404430E+01 Mg 0.802230E+01 0.120003E+02 0.404430E+01 Mg 0.802230E+01 0.802230E+01 0.802230E+01 Mg 0.120003E+02 0.120003E+02 0.802230E+01 Mg 0.120003E+02 0.802230E+01 0.120003E+02 Mg 0.802230E+01 0.120003E+02 0.120003E+02 O 0.663000E-01 0.663000E-01 0.404430E+01 O 0.404430E+01 0.404430E+01 0.404430E+01 O 0.404430E+01 0.663000E-01 0.802230E+01 O 0.663000E-01 0.404430E+01 0.802230E+01 O 0.663000E-01 0.663000E-01 0.120003E+02 O 0.404430E+01 0.404430E+01 0.120003E+02 O 0.404430E+01 0.663000E-01 0.159783E+02 O 0.663000E-01 0.404430E+01 0.159783E+02 O 0.663000E-01 0.802230E+01 0.404430E+01 O 0.404430E+01 0.120003E+02 0.404430E+01 O 0.404430E+01 0.802230E+01 0.802230E+01 O 0.663000E-01 0.120003E+02 0.802230E+01 O 0.663000E-01 0.802230E+01 0.120003E+02 O 0.404430E+01 0.120003E+02 0.120003E+02 O 0.404430E+01 0.802230E+01 0.159783E+02 O 0.663000E-01 0.120003E+02 0.159783E+02 O 0.802230E+01 0.663000E-01 0.404430E+01 O 0.120003E+02 0.404430E+01 0.404430E+01 O 0.120003E+02 0.663000E-01 0.802230E+01 O 0.802230E+01 0.404430E+01 0.802230E+01 O 0.802230E+01 0.663000E-01 0.120003E+02 O 0.120003E+02 0.404430E+01 0.120003E+02 O 0.120003E+02 0.663000E-01 0.159783E+02 O 0.802230E+01 0.404430E+01 0.159783E+02 O 0.802230E+01 0.802230E+01 0.404430E+01 O 0.120003E+02 0.120003E+02 0.404430E+01 O 0.120003E+02 0.802230E+01 0.802230E+01 O 0.802230E+01 0.120003E+02 0.802230E+01 O 0.802230E+01 0.802230E+01 0.120003E+02 O 0.120003E+02 0.120003E+02 0.120003E+02 O 0.120003E+02 0.802230E+01 0.159783E+02 O 0.802230E+01 0.120003E+02 0.159783E+02 ATOMIC_VELOCITIES Mg 0.000000E+00 0.000000E+00 0.000000E+00 Mg 0.000000E+00 0.000000E+00 0.000000E+00 Mg 0.000000E+00 0.000000E+00 0.000000E+00 Mg 0.000000E+00 0.000000E+00 0.000000E+00 Mg 0.000000E+00 0.000000E+00 0.000000E+00 Mg 0.000000E+00 0.000000E+00 0.000000E+00 Mg 0.000000E+00 0.000000E+00 0.000000E+00 Mg 0.000000E+00 0.000000E+00 0.000000E+00 Mg 0.000000E+00 0.000000E+00 0.000000E+00 Mg 0.000000E+00 0.000000E+00 0.000000E+00 Mg 0.000000E+00 0.000000E+00 0.000000E+00 Mg 0.000000E+00 0.000000E+00 0.000000E+00 Mg 0.000000E+00 0.000000E+00 0.000000E+00 Mg 0.000000E+00 0.000000E+00 0.000000E+00 Mg 0.000000E+00 0.000000E+00 0.000000E+00 Mg 0.000000E+00 0.000000E+00 0.000000E+00 Mg 0.000000E+00 0.000000E+00 0.000000E+00 Mg 0.000000E+00 0.000000E+00 0.000000E+00 Mg 0.000000E+00 0.000000E+00 0.000000E+00 Mg 0.000000E+00 0.000000E+00 0.000000E+00 Mg 0.000000E+00 0.000000E+00 0.000000E+00 Mg 0.000000E+00 0.000000E+00 0.000000E+00 Mg 0.000000E+00 0.000000E+00 0.000000E+00 Mg 0.000000E+00 0.000000E+00 0.000000E+00 Mg 0.000000E+00 0.000000E+00 0.000000E+00 Mg 0.000000E+00 0.000000E+00 0.000000E+00 Mg 0.000000E+00 0.000000E+00 0.000000E+00 Mg 0.000000E+00 0.000000E+00 0.000000E+00 Mg 0.000000E+00 0.000000E+00 0.000000E+00 Mg 0.000000E+00 0.000000E+00 0.000000E+00 Mg 0.000000E+00 0.000000E+00 0.000000E+00 Mg 0.000000E+00 0.000000E+00 0.000000E+00 O 0.000000E+00 0.000000E+00 0.000000E+00 O 0.000000E+00 0.000000E+00 0.000000E+00 O 0.000000E+00 0.000000E+00 0.000000E+00 O 0.000000E+00 0.000000E+00 0.000000E+00 O 0.000000E+00 0.000000E+00 0.000000E+00 O 0.000000E+00 0.000000E+00 0.000000E+00 O 0.000000E+00 0.000000E+00 0.000000E+00 O 0.000000E+00 0.000000E+00 0.000000E+00 O 0.000000E+00 0.000000E+00 0.000000E+00 O 0.000000E+00 0.000000E+00 0.000000E+00 O 0.000000E+00 0.000000E+00 0.000000E+00 O 0.000000E+00 0.000000E+00 0.000000E+00 O 0.000000E+00 0.000000E+00 0.000000E+00 O 0.000000E+00 0.000000E+00 0.000000E+00 O 0.000000E+00 0.000000E+00 0.000000E+00 O 0.000000E+00 0.000000E+00 0.000000E+00 O 0.000000E+00 0.000000E+00 0.000000E+00 O 0.000000E+00 0.000000E+00 0.000000E+00 O 0.000000E+00 0.000000E+00 0.000000E+00 O 0.000000E+00 0.000000E+00 0.000000E+00 O 0.000000E+00 0.000000E+00 0.000000E+00 O 0.000000E+00 0.000000E+00 0.000000E+00 O 0.000000E+00 0.000000E+00 0.000000E+00 O 0.000000E+00 0.000000E+00 0.000000E+00 O 0.000000E+00 0.000000E+00 0.000000E+00 O 0.000000E+00 0.000000E+00 0.000000E+00 O 0.000000E+00 0.000000E+00 0.000000E+00 O 0.000000E+00 0.000000E+00 0.000000E+00 O 0.000000E+00 0.000000E+00 0.000000E+00 O 0.000000E+00 0.000000E+00 0.000000E+00 O 0.000000E+00 0.000000E+00 0.000000E+00 O 0.000000E+00 0.000000E+00 0.000000E+00 Forces acting on atoms (au): Mg -0.649832E-05 -0.628798E-05 0.195571E-02 Mg -0.630989E-05 -0.645165E-05 0.195565E-02 Mg -0.653127E-05 -0.649702E-05 0.195561E-02 Mg -0.637109E-05 -0.657283E-05 0.195564E-02 Mg -0.651775E-05 -0.627099E-05 0.195564E-02 Mg -0.664243E-05 -0.651217E-05 0.195581E-02 Mg -0.636087E-05 -0.640417E-05 0.195556E-02 Mg -0.637183E-05 -0.641185E-05 0.195556E-02 Mg -0.642386E-05 -0.648346E-05 0.195559E-02 Mg -0.650301E-05 -0.640031E-05 0.195551E-02 Mg -0.659315E-05 -0.657916E-05 0.195562E-02 Mg -0.656403E-05 -0.629709E-05 0.195549E-02 Mg -0.647184E-05 -0.646337E-05 0.195575E-02 Mg -0.649410E-05 -0.641195E-05 0.195564E-02 Mg -0.640628E-05 -0.658838E-05 0.195562E-02 Mg -0.643621E-05 -0.656990E-05 0.195557E-02 Mg -0.635782E-05 -0.635625E-05 0.195564E-02 Mg -0.631268E-05 -0.623755E-05 0.195569E-02 Mg -0.648748E-05 -0.638138E-05 0.195564E-02 Mg -0.633606E-05 -0.650411E-05 0.195551E-02 Mg -0.635631E-05 -0.653026E-05 0.195567E-02 Mg -0.629815E-05 -0.646136E-05 0.195563E-02 Mg -0.654692E-05 -0.629069E-05 0.195566E-02 Mg -0.645116E-05 -0.646623E-05 0.195567E-02 Mg -0.638365E-05 -0.644341E-05 0.195558E-02 Mg -0.647191E-05 -0.639863E-05 0.195553E-02 Mg -0.643854E-05 -0.633238E-05 0.195555E-02 Mg -0.646869E-05 -0.642115E-05 0.195550E-02 Mg -0.633885E-05 -0.638086E-05 0.195564E-02 Mg -0.642009E-05 -0.653014E-05 0.195564E-02 Mg -0.631911E-05 -0.641705E-05 0.195555E-02 Mg -0.641786E-05 -0.643383E-05 0.195547E-02 O -0.771989E-04 -0.780493E-04 -0.204013E-02 O -0.784301E-04 -0.779739E-04 -0.203989E-02 O -0.776212E-04 -0.787391E-04 -0.203973E-02 O -0.777262E-04 -0.782079E-04 -0.203983E-02 O -0.782545E-04 -0.777601E-04 -0.204004E-02 O -0.771697E-04 -0.786720E-04 -0.204023E-02 O -0.776499E-04 -0.781590E-04 -0.203907E-02 O -0.785355E-04 -0.781766E-04 -0.204020E-02 O -0.778548E-04 -0.785407E-04 -0.204076E-02 O -0.787347E-04 -0.784281E-04 -0.203964E-02 O -0.774420E-04 -0.787714E-04 -0.203966E-02 O -0.793503E-04 -0.784549E-04 -0.203928E-02 O -0.773382E-04 -0.778549E-04 -0.204029E-02 O -0.786335E-04 -0.779923E-04 -0.203940E-02 O -0.779783E-04 -0.784798E-04 -0.204032E-02 O -0.780622E-04 -0.782377E-04 -0.204013E-02 O -0.782255E-04 -0.786974E-04 -0.203988E-02 O -0.778113E-04 -0.784045E-04 -0.203978E-02 O -0.784041E-04 -0.782231E-04 -0.204003E-02 O -0.776132E-04 -0.781861E-04 -0.203986E-02 O -0.788691E-04 -0.786263E-04 -0.203972E-02 O -0.774608E-04 -0.786141E-04 -0.203932E-02 O -0.786061E-04 -0.779920E-04 -0.204014E-02 O -0.781905E-04 -0.773731E-04 -0.204037E-02 O -0.773647E-04 -0.783484E-04 -0.204021E-02 O -0.791976E-04 -0.782701E-04 -0.204051E-02 O -0.780675E-04 -0.777471E-04 -0.203984E-02 O -0.778008E-04 -0.779806E-04 -0.203978E-02 O -0.781168E-04 -0.787094E-04 -0.203980E-02 O -0.778605E-04 -0.783572E-04 -0.204067E-02 O -0.779739E-04 -0.778978E-04 -0.203909E-02 O -0.778582E-04 -0.781345E-04 -0.204018E-02 Partial temperatures (for each ionic specie) Species Temp (K) Mean Square Displacement (a.u.) 1 0.00 0.0000 2 0.00 0.0000 nfi ekinc temph tempp etot enthal econs econt vnhh xnhh0 vnhp xnhp0 nfi tempp E -T.S-mu.nbsp +K_p #Iter Step 2 0 -547.28628 -547.28628 -547.28628 27 Elct. dipole 1 15.415051 Ionic dipole 1 1.060800 writing restart file: /mgo_50.save restart file written in 0.175 sec. Averaged Physical Quantities accomulated this run ekinc : 0.00000 0.00000 (AU) ekin : 235.43101 235.43095 (AU) epot : -891.54169 -891.54144 (AU) total energy : -547.27818 -547.28628 (AU) temperature : 0.00000 0.00000 (K ) enthalpy : 0.00000 0.00000 (AU) econs : -547.27818 -547.28628 (AU) pressure : -25.94390 0.00000 (Gpa) volume : 4028.78703 4028.78703 (AU) initialize : 3.30s CPU total_time : 267.86s CPU formf : 0.06s CPU rhoofr : 18.59s CPU ( 54 calls, 0.344 s avg) vofrho : 12.60s CPU ( 53 calls, 0.238 s avg) dforce : 29.29s CPU ( 1728 calls, 0.017 s avg) calphi : 1.38s CPU ( 27 calls, 0.051 s avg) gram : 42.23s CPU ( 54 calls, 0.782 s avg) newd : 26.63s CPU ( 27 calls, 0.986 s avg) calbec : 16.87s CPU ( 186 calls, 0.091 s avg) prefor : 0.38s CPU ( 29 calls, 0.013 s avg) strucf : 0.01s CPU nlfl : 0.01s CPU nlfq : 0.27s CPU set_cc : 1.98s CPU ( 54 calls, 0.037 s avg) rhov : 4.99s CPU ( 54 calls, 0.092 s avg) nlsm1 : 16.87s CPU ( 186 calls, 0.091 s avg) nlsm2 : 0.27s CPU forcecc : 3.58s CPU ( 53 calls, 0.068 s avg) fft : 5.90s CPU ( 267 calls, 0.022 s avg) ffts : 0.48s CPU ( 107 calls, 0.004 s avg) fftw : 17.80s CPU ( 6912 calls, 0.003 s avg) fftb : 20.06s CPU ( 21120 calls, 0.001 s avg) runcg_uspp : 267.76s CPU calcmt : 0.00s CPU pc2 : 7.18s CPU ( 78 calls, 0.092 s avg) pcdaga2 : 17.60s CPU ( 26 calls, 0.677 s avg) set_x_minus1 : 1.82s CPU ( 2 calls, 0.912 s avg) xminus1 : 15.60s CPU ( 52 calls, 0.300 s avg) emass_p_tpa : 0.00s CPU CP : 4m31.33s CPU time, 4m35.07s wall time This run was terminated on: 16:34:51 28Apr2008 =------------------------------------------------------------------------------= JOB DONE. =------------------------------------------------------------------------------=