quantum-espresso/PW/examples/ESM_example/reference/H2O_pbc.out

312 lines
12 KiB
Plaintext

Program PWSCF v.6.1 (svn rev. 13369) starts on 21Aug2017 at 23: 6: 1
This program is part of the open-source Quantum ESPRESSO suite
for quantum simulation of materials; please cite
"P. Giannozzi et al., J. Phys.:Condens. Matter 21 395502 (2009);
URL http://www.quantum-espresso.org",
in publications or presentations arising from this work. More details at
http://www.quantum-espresso.org/quote
Parallel version (MPI & OpenMP), running on 8 processor cores
Number of MPI processes: 4
Threads/MPI process: 2
R & G space division: proc/nbgrp/npool/nimage = 4
Waiting for input...
Reading input from standard input
Current dimensions of program PWSCF are:
Max number of different atomic species (ntypx) = 10
Max number of k-points (npk) = 40000
Max angular momentum in pseudopotentials (lmaxx) = 3
gamma-point specific algorithms are used
Subspace diagonalization in iterative solution of the eigenvalue problem:
a serial algorithm will be used
Parallelization info
--------------------
sticks: dense smooth PW G-vecs: dense smooth PW
Min 1591 794 198 114643 40525 5074
Max 1594 797 199 114648 40530 5081
Sum 6369 3181 793 458581 162113 20303
bravais-lattice index = 0
lattice parameter (alat) = 20.0000 a.u.
unit-cell volume = 9599.9886 (a.u.)^3
number of atoms/cell = 3
number of atomic types = 2
number of electrons = 8.00
number of Kohn-Sham states= 8
kinetic-energy cutoff = 25.0000 Ry
charge density cutoff = 200.0000 Ry
convergence threshold = 1.0E-06
mixing beta = 0.3000
number of iterations used = 8 plain mixing
Exchange-correlation = SLA PW PBE PBE ( 1 4 3 4 0 0)
Effective Screening Medium Method
=================================
Ordinary Periodic Boundary Conditions
grid points for fit at edges = 4
celldm(1)= 19.999992 celldm(2)= 0.000000 celldm(3)= 0.000000
celldm(4)= 0.000000 celldm(5)= 0.000000 celldm(6)= 0.000000
crystal axes: (cart. coord. in units of alat)
a(1) = ( 1.000000 0.000000 0.000000 )
a(2) = ( 0.000000 1.000000 0.000000 )
a(3) = ( 0.000000 0.000000 1.200000 )
reciprocal axes: (cart. coord. in units 2 pi/alat)
b(1) = ( 1.000000 0.000000 0.000000 )
b(2) = ( 0.000000 1.000000 0.000000 )
b(3) = ( 0.000000 0.000000 0.833333 )
PseudoPot. # 1 for H read from file:
/Users/otani/ESPRESSO/qe-6.1_stress/pseudo/H.pbe-van_ak.UPF
MD5 check sum: 077eb6d537518a38cb46c6de387227b7
Pseudo is Ultrasoft, Zval = 1.0
Generated by new atomic code, or converted to UPF format
Using radial grid of 615 points, 1 beta functions with:
l(1) = 0
Q(r) pseudized with 8 coefficients, rinner = 0.800
PseudoPot. # 2 for O read from file:
/Users/otani/ESPRESSO/qe-6.1_stress/pseudo/O.pbe-van_ak.UPF
MD5 check sum: 78c19812648e8988a0fe644b0f25305e
Pseudo is Ultrasoft, Zval = 6.0
Generated by new atomic code, or converted to UPF format
Using radial grid of 737 points, 4 beta functions with:
l(1) = 0
l(2) = 0
l(3) = 1
l(4) = 1
Q(r) pseudized with 8 coefficients, rinner = 0.800 0.800 0.800
atomic species valence mass pseudopotential
H 1.00 1.00794 H ( 1.00)
O 6.00 55.84700 O ( 1.00)
No symmetry found
Cartesian axes
site n. atom positions (alat units)
1 O tau( 1) = ( 0.0000000 0.5000000 0.0000000 )
2 H tau( 2) = ( 0.0431388 0.4310286 0.0430783 )
3 H tau( 3) = ( 0.0366354 0.5764064 0.0359492 )
number of k points= 1 Methfessel-Paxton smearing, width (Ry)= 0.0300
cart. coord. in units 2pi/alat
k( 1) = ( 0.0000000 0.0000000 0.0000000), wk = 2.0000000
Dense grid: 229291 G-vectors FFT dimensions: ( 96, 96, 120)
Smooth grid: 81057 G-vectors FFT dimensions: ( 64, 64, 80)
Estimated max dynamical RAM per process > 45.72MB
Estimated total allocated dynamical RAM > 182.89MB
Initial potential from superposition of free atoms
Check: negative starting charge= -0.002904
starting charge 7.99998, renormalised to 8.00000
negative rho (up, down): 2.904E-03 0.000E+00
Starting wfc are 6 randomized atomic wfcs + 2 random wfc
total cpu time spent up to now is 1.1 secs
Self-consistent Calculation
iteration # 1 ecut= 25.00 Ry beta=0.30
Davidson diagonalization with overlap
ethr = 1.00E-02, avg # of iterations = 7.0
negative rho (up, down): 2.997E-03 0.000E+00
total cpu time spent up to now is 1.8 secs
total energy = -34.38954834 Ry
Harris-Foulkes estimate = -34.57127197 Ry
estimated scf accuracy < 0.44378693 Ry
iteration # 2 ecut= 25.00 Ry beta=0.30
Davidson diagonalization with overlap
ethr = 5.55E-03, avg # of iterations = 2.0
negative rho (up, down): 4.701E-03 0.000E+00
total cpu time spent up to now is 2.5 secs
total energy = -34.42662395 Ry
Harris-Foulkes estimate = -34.44037537 Ry
estimated scf accuracy < 0.04574427 Ry
iteration # 3 ecut= 25.00 Ry beta=0.30
Davidson diagonalization with overlap
ethr = 5.72E-04, avg # of iterations = 3.0
negative rho (up, down): 9.317E-03 0.000E+00
total cpu time spent up to now is 3.2 secs
total energy = -34.42836850 Ry
Harris-Foulkes estimate = -34.43238564 Ry
estimated scf accuracy < 0.01327854 Ry
iteration # 4 ecut= 25.00 Ry beta=0.30
Davidson diagonalization with overlap
ethr = 1.66E-04, avg # of iterations = 2.0
negative rho (up, down): 1.083E-02 0.000E+00
total cpu time spent up to now is 3.9 secs
total energy = -34.42865574 Ry
Harris-Foulkes estimate = -34.42992413 Ry
estimated scf accuracy < 0.00292845 Ry
iteration # 5 ecut= 25.00 Ry beta=0.30
Davidson diagonalization with overlap
ethr = 3.66E-05, avg # of iterations = 4.0
negative rho (up, down): 9.200E-03 0.000E+00
total cpu time spent up to now is 4.7 secs
total energy = -34.42913698 Ry
Harris-Foulkes estimate = -34.42915694 Ry
estimated scf accuracy < 0.00003482 Ry
iteration # 6 ecut= 25.00 Ry beta=0.30
Davidson diagonalization with overlap
ethr = 4.35E-07, avg # of iterations = 8.0
negative rho (up, down): 8.950E-03 0.000E+00
total cpu time spent up to now is 5.5 secs
total energy = -34.42914393 Ry
Harris-Foulkes estimate = -34.42914460 Ry
estimated scf accuracy < 0.00001397 Ry
iteration # 7 ecut= 25.00 Ry beta=0.30
Davidson diagonalization with overlap
ethr = 1.75E-07, avg # of iterations = 2.0
negative rho (up, down): 8.735E-03 0.000E+00
total cpu time spent up to now is 6.2 secs
total energy = -34.42914495 Ry
Harris-Foulkes estimate = -34.42914565 Ry
estimated scf accuracy < 0.00000102 Ry
iteration # 8 ecut= 25.00 Ry beta=0.30
Davidson diagonalization with overlap
ethr = 1.27E-08, avg # of iterations = 2.0
negative rho (up, down): 8.989E-03 0.000E+00
total cpu time spent up to now is 6.8 secs
End of self-consistent calculation
k = 0.0000 0.0000 0.0000 ( 10152 PWs) bands (ev):
-25.2127 -12.9940 -9.1815 -7.1440 -1.0677 -0.0518 0.3743 0.4918
the Fermi energy is -6.8004 ev
! total energy = -34.42914446 Ry
Harris-Foulkes estimate = -34.42914511 Ry
estimated scf accuracy < 0.00000018 Ry
The total energy is the sum of the following terms:
one-electron contribution = -65.54740095 Ry
hartree contribution = 34.21803356 Ry
xc contribution = -8.40337372 Ry
ewald contribution = 5.30185727 Ry
smearing contrib. (-TS) = 0.00173937 Ry
convergence has been achieved in 8 iterations
Forces acting on atoms (cartesian axes, Ry/au):
atom 1 type 2 force = -0.00183981 -0.00032426 -0.00154706
atom 2 type 1 force = 0.00089110 -0.00025910 0.00075452
atom 3 type 1 force = 0.00094871 0.00058337 0.00079255
Total force = 0.003030 Total SCF correction = 0.001216
SCF correction compared to forces is large: reduce conv_thr to get better values
Writing output data file H2O_pbc.save
init_run : 0.99s CPU 1.01s WALL ( 1 calls)
electrons : 5.68s CPU 5.85s WALL ( 1 calls)
forces : 0.37s CPU 0.40s WALL ( 1 calls)
Called by init_run:
wfcinit : 0.03s CPU 0.03s WALL ( 1 calls)
potinit : 0.52s CPU 0.51s WALL ( 1 calls)
Called by electrons:
c_bands : 0.70s CPU 0.71s WALL ( 8 calls)
sum_band : 1.53s CPU 1.48s WALL ( 8 calls)
v_of_rho : 2.41s CPU 2.40s WALL ( 9 calls)
newd : 0.72s CPU 0.84s WALL ( 9 calls)
mix_rho : 0.31s CPU 0.39s WALL ( 8 calls)
Called by c_bands:
init_us_2 : 0.01s CPU 0.02s WALL ( 17 calls)
regterg : 0.68s CPU 0.70s WALL ( 8 calls)
Called by sum_band:
sum_band:bec : 0.00s CPU 0.00s WALL ( 8 calls)
addusdens : 0.92s CPU 0.85s WALL ( 8 calls)
Called by *egterg:
h_psi : 0.63s CPU 0.64s WALL ( 39 calls)
s_psi : 0.01s CPU 0.00s WALL ( 39 calls)
g_psi : 0.00s CPU 0.00s WALL ( 30 calls)
rdiaghg : 0.01s CPU 0.01s WALL ( 38 calls)
Called by h_psi:
h_psi:pot : 0.63s CPU 0.64s WALL ( 39 calls)
h_psi:calbec : 0.01s CPU 0.01s WALL ( 39 calls)
vloc_psi : 0.61s CPU 0.62s WALL ( 39 calls)
add_vuspsi : 0.01s CPU 0.01s WALL ( 39 calls)
General routines
calbec : 0.01s CPU 0.01s WALL ( 51 calls)
fft : 2.62s CPU 2.44s WALL ( 143 calls)
ffts : 0.08s CPU 0.09s WALL ( 17 calls)
fftw : 0.60s CPU 0.61s WALL ( 250 calls)
interpolate : 0.42s CPU 0.42s WALL ( 17 calls)
davcio : 0.00s CPU 0.00s WALL ( 1 calls)
Parallel routines
fft_scatter : 1.25s CPU 1.33s WALL ( 410 calls)
PWSCF : 7.14s CPU 7.39s WALL
This run was terminated on: 23: 6: 9 21Aug2017
=------------------------------------------------------------------------------=
JOB DONE.
=------------------------------------------------------------------------------=