quantum-espresso/PW/examples/vdwDF_example/reference/water.scf.out

629 lines
25 KiB
Plaintext

Program PWSCF v.6.4.1 starts on 23Sep2019 at 12:39:59
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);
"P. Giannozzi et al., J. Phys.:Condens. Matter 29 465901 (2017);
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), running on 4 processors
MPI processes distributed on 1 nodes
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
IMPORTANT: XC functional enforced from input :
Exchange-correlation= VDW-DF
( 1 4 4 0 1 0 0)
Any further DFT definition will be discarded
Please, verify this is what you really want
file O.pbe-rrkjus.UPF: wavefunction(s) 2S renormalized
file H.pbe-rrkjus.UPF: wavefunction(s) 1S renormalized
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 768 512 127 40311 21940 2739
Max 770 513 130 40318 21949 2746
Sum 3077 2049 515 161263 87777 10971
bravais-lattice index = 8
lattice parameter (alat) = 15.0000 a.u.
unit-cell volume = 3953.7707 (a.u.)^3
number of atoms/cell = 6
number of atomic types = 2
number of electrons = 16.00
number of Kohn-Sham states= 8
kinetic-energy cutoff = 30.0000 Ry
charge density cutoff = 180.0000 Ry
convergence threshold = 1.0E-08
mixing beta = 0.7000
number of iterations used = 8 plain mixing
Exchange-correlation= VDW-DF
( 1 4 4 0 1 0 0)
celldm(1)= 15.000000 celldm(2)= 0.954545 celldm(3)= 1.227273
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 0.954545 0.000000 )
a(3) = ( 0.000000 0.000000 1.227273 )
reciprocal axes: (cart. coord. in units 2 pi/alat)
b(1) = ( 1.000000 0.000000 0.000000 )
b(2) = ( 0.000000 1.047619 0.000000 )
b(3) = ( 0.000000 0.000000 0.814815 )
PseudoPot. # 1 for O read from file:
/benchmarks/pseudo/O.pbe-rrkjus.UPF
MD5 check sum: fadcf19ee70a498d3030e2e79cf929a4
Pseudo is Ultrasoft, Zval = 6.0
Generated by new atomic code, or converted to UPF format
Using radial grid of 1269 points, 4 beta functions with:
l(1) = 0
l(2) = 0
l(3) = 1
l(4) = 1
Q(r) pseudized with 0 coefficients
PseudoPot. # 2 for H read from file:
/benchmarks/pseudo/H.pbe-rrkjus.UPF
MD5 check sum: 2d52a7f45632fd764c17cc6affed14c8
Pseudo is Ultrasoft, Zval = 1.0
Generated by new atomic code, or converted to UPF format
Using radial grid of 1061 points, 2 beta functions with:
l(1) = 0
l(2) = 0
Q(r) pseudized with 0 coefficients
atomic species valence mass pseudopotential
O 6.00 15.99940 O ( 1.00)
H 1.00 1.00794 H ( 1.00)
2 Sym. Ops. (no inversion) found
s frac. trans.
isym = 1 identity
cryst. s( 1) = ( 1 0 0 )
( 0 1 0 )
( 0 0 1 )
cart. s( 1) = ( 1.0000000 0.0000000 0.0000000 )
( 0.0000000 1.0000000 0.0000000 )
( 0.0000000 0.0000000 1.0000000 )
isym = 2 inv. 180 deg rotation - cart. axis [1,0,0]
cryst. s( 2) = ( -1 0 0 )
( 0 1 0 )
( 0 0 1 )
cart. s( 2) = ( -1.0000000 0.0000000 0.0000000 )
( 0.0000000 1.0000000 0.0000000 )
( 0.0000000 0.0000000 1.0000000 )
point group C_s (m)
there are 2 classes
the character table:
E s
A' 1.00 1.00
A'' 1.00 -1.00
the symmetry operations in each class and the name of the first element:
E 1
identity
s 2
inv. 180 deg rotation - cart. axis [1,0,0]
Cartesian axes
site n. atom positions (alat units)
1 O tau( 1) = ( 0.0000000 0.0016540 -0.0072484 )
2 H tau( 2) = ( 0.0000000 0.0981485 -0.0826521 )
3 H tau( 3) = ( 0.0000000 0.0490883 0.1065556 )
4 O tau( 4) = ( 0.0000000 0.1117595 0.3550478 )
5 H tau( 5) = ( -0.0975766 0.0656956 0.4133167 )
6 H tau( 6) = ( 0.0975766 0.0656956 0.4133167 )
Crystallographic axes
site n. atom positions (cryst. coord.)
1 O tau( 1) = ( 0.0000000 0.0017328 -0.0059061 )
2 H tau( 2) = ( 0.0000000 0.1028222 -0.0673461 )
3 H tau( 3) = ( 0.0000000 0.0514258 0.0868231 )
4 O tau( 4) = ( 0.0000000 0.1170814 0.2892982 )
5 H tau( 5) = ( -0.0975766 0.0688239 0.3367766 )
6 H tau( 6) = ( 0.0975766 0.0688239 0.3367766 )
number of k points= 1
cart. coord. in units 2pi/alat
k( 1) = ( 0.0000000 0.0000000 0.0000000), wk = 2.0000000
cryst. coord.
k( 1) = ( 0.0000000 0.0000000 0.0000000), wk = 2.0000000
Dense grid: 80632 G-vectors FFT dimensions: ( 72, 64, 80)
Smooth grid: 43889 G-vectors FFT dimensions: ( 54, 50, 72)
Dynamical RAM for wfc: 0.17 MB
Dynamical RAM for wfc (w. buffer): 0.17 MB
Dynamical RAM for str. fact: 0.62 MB
Dynamical RAM for local pot: 0.00 MB
Dynamical RAM for nlocal pot: 0.50 MB
Dynamical RAM for qrad: 1.47 MB
Dynamical RAM for rho,v,vnew: 3.03 MB
Dynamical RAM for rhoin: 1.01 MB
Dynamical RAM for rho*nmix: 4.92 MB
Dynamical RAM for G-vectors: 1.24 MB
Dynamical RAM for h,s,v(r/c): 0.02 MB
Dynamical RAM for <psi|beta>: 0.00 MB
Dynamical RAM for psi: 0.67 MB
Dynamical RAM for hpsi: 0.67 MB
Dynamical RAM for spsi: 0.67 MB
Dynamical RAM for wfcinit/wfcrot: 0.51 MB
Dynamical RAM for addusdens: 14.46 MB
Dynamical RAM for addusforce: 18.31 MB
Dynamical RAM for addusstress: 15.38 MB
Estimated static dynamical RAM per process > 17.39 MB
Estimated max dynamical RAM per process > 35.70 MB
Estimated total dynamical RAM > 142.79 MB
Initial potential from superposition of free atoms
starting charge 15.61518, renormalised to 16.00000
negative rho (up, down): 2.810E-05 0.000E+00
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% %
% You are using vdW-DF, which was implemented by the Thonhauser group. %
% Please cite the following two papers that made this development %
% possible and the two reviews that describe the various versions: %
% %
% T. Thonhauser et al., PRL 115, 136402 (2015). %
% T. Thonhauser et al., PRB 76, 125112 (2007). %
% K. Berland et al., Rep. Prog. Phys. 78, 066501 (2015). %
% D.C. Langreth et al., J. Phys.: Condens. Matter 21, 084203 (2009). %
% %
% %
% If you are calculating the stress with vdW-DF, please also cite: %
% %
% R. Sabatini et al., J. Phys.: Condens. Matter 24, 424209 (2012). %
% %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Carrying out vdW-DF run using the following parameters:
Nqs = 20 Npoints = 1024 r_max = 100.000
q_mesh = 0.00001000 0.04494208 0.09755937 0.15916263
0.23128650 0.31572767 0.41458969 0.53033537
0.66584808 0.82450364 1.01025438 1.22772762
1.48234092 1.78043706 2.12944203 2.53805004
3.01644009 3.57652955 4.23227104 5.00000000
-----------------------------------------------
Non-local corr. energy = 0.27494246 Ry
-----------------------------------------------
Starting wfcs are 12 randomized atomic wfcs
total cpu time spent up to now is 26.5 secs
Self-consistent Calculation
iteration # 1 ecut= 30.00 Ry beta= 0.70
Davidson diagonalization with overlap
ethr = 1.00E-02, avg # of iterations = 2.0
negative rho (up, down): 2.476E-05 0.000E+00
-----------------------------------------------
Non-local corr. energy = 0.30229411 Ry
-----------------------------------------------
total cpu time spent up to now is 26.9 secs
total energy = -68.61255685 Ry
Harris-Foulkes estimate = -69.74358395 Ry
estimated scf accuracy < 1.47053641 Ry
iteration # 2 ecut= 30.00 Ry beta= 0.70
Davidson diagonalization with overlap
ethr = 9.19E-03, avg # of iterations = 2.0
negative rho (up, down): 8.259E-05 0.000E+00
-----------------------------------------------
Non-local corr. energy = 0.29919547 Ry
-----------------------------------------------
total cpu time spent up to now is 27.3 secs
total energy = -68.82976620 Ry
Harris-Foulkes estimate = -69.34187945 Ry
estimated scf accuracy < 0.99498446 Ry
iteration # 3 ecut= 30.00 Ry beta= 0.70
Davidson diagonalization with overlap
ethr = 6.22E-03, avg # of iterations = 2.0
negative rho (up, down): 7.154E-04 0.000E+00
-----------------------------------------------
Non-local corr. energy = 0.30435425 Ry
-----------------------------------------------
total cpu time spent up to now is 27.7 secs
total energy = -69.04870183 Ry
Harris-Foulkes estimate = -69.06820230 Ry
estimated scf accuracy < 0.03427255 Ry
iteration # 4 ecut= 30.00 Ry beta= 0.70
Davidson diagonalization with overlap
ethr = 2.14E-04, avg # of iterations = 2.0
negative rho (up, down): 6.256E-04 0.000E+00
-----------------------------------------------
Non-local corr. energy = 0.30390405 Ry
-----------------------------------------------
total cpu time spent up to now is 28.1 secs
total energy = -69.05570611 Ry
Harris-Foulkes estimate = -69.05607760 Ry
estimated scf accuracy < 0.00064847 Ry
iteration # 5 ecut= 30.00 Ry beta= 0.70
Davidson diagonalization with overlap
ethr = 4.05E-06, avg # of iterations = 2.0
negative rho (up, down): 1.768E-03 0.000E+00
-----------------------------------------------
Non-local corr. energy = 0.30409047 Ry
-----------------------------------------------
total cpu time spent up to now is 28.5 secs
total energy = -69.05590281 Ry
Harris-Foulkes estimate = -69.05588163 Ry
estimated scf accuracy < 0.00003152 Ry
iteration # 6 ecut= 30.00 Ry beta= 0.70
Davidson diagonalization with overlap
ethr = 1.97E-07, avg # of iterations = 2.0
negative rho (up, down): 1.877E-03 0.000E+00
-----------------------------------------------
Non-local corr. energy = 0.30407639 Ry
-----------------------------------------------
total cpu time spent up to now is 28.9 secs
total energy = -69.05590617 Ry
Harris-Foulkes estimate = -69.05590851 Ry
estimated scf accuracy < 0.00000286 Ry
iteration # 7 ecut= 30.00 Ry beta= 0.70
Davidson diagonalization with overlap
ethr = 1.79E-08, avg # of iterations = 2.0
negative rho (up, down): 1.890E-03 0.000E+00
-----------------------------------------------
Non-local corr. energy = 0.30407454 Ry
-----------------------------------------------
total cpu time spent up to now is 29.3 secs
total energy = -69.05590651 Ry
Harris-Foulkes estimate = -69.05590663 Ry
estimated scf accuracy < 0.00000007 Ry
iteration # 8 ecut= 30.00 Ry beta= 0.70
Davidson diagonalization with overlap
ethr = 4.43E-10, avg # of iterations = 2.0
negative rho (up, down): 1.893E-03 0.000E+00
-----------------------------------------------
Non-local corr. energy = 0.30407472 Ry
-----------------------------------------------
total cpu time spent up to now is 29.7 secs
total energy = -69.05590648 Ry
Harris-Foulkes estimate = -69.05590653 Ry
estimated scf accuracy < 0.00000003 Ry
iteration # 9 ecut= 30.00 Ry beta= 0.70
Davidson diagonalization with overlap
ethr = 1.69E-10, avg # of iterations = 2.0
negative rho (up, down): 1.891E-03 0.000E+00
-----------------------------------------------
Non-local corr. energy = 0.30407408 Ry
-----------------------------------------------
total cpu time spent up to now is 30.1 secs
End of self-consistent calculation
k = 0.0000 0.0000 0.0000 ( 5486 PWs) bands (ev):
-25.6440 -24.2672 -13.5672 -12.2977 -9.7600 -8.3490 -7.6830 -6.4199
occupation numbers
1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000
highest occupied level (ev): -6.4199
! total energy = -69.05590650 Ry
Harris-Foulkes estimate = -69.05590649 Ry
estimated scf accuracy < 8.8E-10 Ry
The total energy is the sum of the following terms:
one-electron contribution = -122.46807785 Ry
hartree contribution = 64.43497759 Ry
xc contribution = -17.37527098 Ry
ewald contribution = 6.35246473 Ry
convergence has been achieved in 9 iterations
Forces acting on atoms (cartesian axes, Ry/au):
atom 1 type 1 force = 0.00000000 0.00095625 -0.00217434
atom 2 type 2 force = 0.00000000 0.00013616 0.00067786
atom 3 type 2 force = 0.00000000 -0.00035421 0.00013294
atom 4 type 1 force = 0.00000000 0.00087773 0.00122786
atom 5 type 2 force = 0.00001972 -0.00080797 0.00006784
atom 6 type 2 force = -0.00001972 -0.00080797 0.00006784
The non-local contrib. to forces
atom 1 type 1 force = 0.00000000 2.37277510 0.53165394
atom 2 type 2 force = 0.00000000 -0.15344377 0.12907039
atom 3 type 2 force = 0.00000000 -0.06954767 -0.18712626
atom 4 type 1 force = 0.00000000 -1.58418865 1.86394271
atom 5 type 2 force = 0.16527621 0.07169165 -0.09125198
atom 6 type 2 force = -0.16527621 0.07169165 -0.09125198
The ionic contribution to forces
atom 1 type 1 force = -0.00000000 -4.70114709 -3.35460747
atom 2 type 2 force = -0.00000000 2.82811782 -2.66092597
atom 3 type 2 force = 0.00000000 1.08890147 2.53590863
atom 4 type 1 force = 0.00000000 3.35642550 -0.50250843
atom 5 type 2 force = -3.10291310 -1.28614885 1.99106662
atom 6 type 2 force = 3.10291310 -1.28614885 1.99106662
The local contribution to forces
atom 1 type 1 force = 0.00000011 2.32987984 2.82104555
atom 2 type 2 force = 0.00000001 -2.67394764 2.53280069
atom 3 type 2 force = 0.00000001 -1.01911854 -2.34837722
atom 4 type 1 force = 0.00000001 -1.77074534 -1.35995075
atom 5 type 2 force = 2.93765576 1.21423850 -1.89947756
atom 6 type 2 force = -2.93765576 1.21423851 -1.89947757
The core correction contribution to forces
atom 1 type 1 force = 0.00000000 0.00000000 0.00000000
atom 2 type 2 force = 0.00000000 0.00000000 0.00000000
atom 3 type 2 force = 0.00000000 0.00000000 0.00000000
atom 4 type 1 force = 0.00000000 0.00000000 0.00000000
atom 5 type 2 force = 0.00000000 0.00000000 0.00000000
atom 6 type 2 force = 0.00000000 0.00000000 0.00000000
The Hubbard contrib. to forces
atom 1 type 1 force = 0.00000000 0.00000000 0.00000000
atom 2 type 2 force = 0.00000000 0.00000000 0.00000000
atom 3 type 2 force = 0.00000000 0.00000000 0.00000000
atom 4 type 1 force = 0.00000000 0.00000000 0.00000000
atom 5 type 2 force = 0.00000000 0.00000000 0.00000000
atom 6 type 2 force = 0.00000000 0.00000000 0.00000000
The SCF correction term to forces
atom 1 type 1 force = -0.00000013 0.00003751 0.00000265
atom 2 type 2 force = -0.00000000 -0.00000113 0.00000176
atom 3 type 2 force = -0.00000001 -0.00000035 -0.00000321
atom 4 type 1 force = -0.00000005 -0.00002466 0.00001333
atom 5 type 2 force = 0.00000084 -0.00000016 -0.00000022
atom 6 type 2 force = -0.00000085 -0.00000016 -0.00000022
Total force = 0.003140 Total SCF correction = 0.000047
Computing stress (Cartesian axis) and pressure
total stress (Ry/bohr**3) (kbar) P= -1.70
-0.00001157 0.00000000 0.00000000 -1.70 0.00 0.00
0.00000000 -0.00001113 -0.00000082 0.00 -1.64 -0.12
0.00000000 -0.00000082 -0.00001194 0.00 -0.12 -1.76
kinetic stress (kbar) 773.91 0.00 0.00
0.00 795.27 35.09
0.00 35.09 748.19
local stress (kbar) -1260.29 0.00 -0.00
0.00 -1418.75 -162.43
-0.00 -162.43 -2492.38
nonloc. stress (kbar) 278.32 0.00 0.00
0.00 280.60 11.20
0.00 11.20 269.13
hartree stress (kbar) 584.70 -0.00 0.00
-0.00 685.81 102.02
0.00 102.02 1126.88
exc-cor stress (kbar) -207.38 -0.00 0.00
-0.00 -207.95 -1.03
0.00 -1.03 -205.76
corecor stress (kbar) 0.00 0.00 0.00
0.00 0.00 0.00
0.00 0.00 0.00
ewald stress (kbar) -173.91 -0.00 0.00
-0.00 -139.58 15.06
0.00 15.06 549.84
hubbard stress (kbar) 0.00 0.00 0.00
0.00 0.00 0.00
0.00 0.00 0.00
london stress (kbar) 0.00 0.00 0.00
0.00 0.00 0.00
0.00 0.00 0.00
DFT-D3 stress (kbar) 0.00 0.00 0.00
0.00 0.00 0.00
0.00 0.00 0.00
XDM stress (kbar) 0.00 0.00 0.00
0.00 0.00 0.00
0.00 0.00 0.00
dft-nl stress (kbar) 2.95 0.00 -0.00
0.00 2.97 -0.04
-0.00 -0.04 2.35
TS-vdW stress (kbar) 0.00 0.00 0.00
0.00 0.00 0.00
0.00 0.00 0.00
Writing output data file /benchmarks/tempdir/water_vdw.save/
init_run : 25.98s CPU 26.03s WALL ( 1 calls)
electrons : 3.38s CPU 3.72s WALL ( 1 calls)
forces : 0.20s CPU 0.21s WALL ( 1 calls)
stress : 0.75s CPU 0.77s WALL ( 1 calls)
Called by init_run:
wfcinit : 0.01s CPU 0.01s WALL ( 1 calls)
wfcinit:atom : 0.00s CPU 0.00s WALL ( 1 calls)
wfcinit:wfcr : 0.01s CPU 0.01s WALL ( 1 calls)
potinit : 25.73s CPU 25.78s WALL ( 1 calls)
hinit0 : 0.20s CPU 0.20s WALL ( 1 calls)
Called by electrons:
c_bands : 0.20s CPU 0.20s WALL ( 9 calls)
sum_band : 0.20s CPU 0.21s WALL ( 9 calls)
v_of_rho : 28.29s CPU 28.50s WALL ( 10 calls)
v_h : 0.03s CPU 0.04s WALL ( 10 calls)
v_xc : 28.25s CPU 28.46s WALL ( 10 calls)
newd : 0.17s CPU 0.19s WALL ( 10 calls)
mix_rho : 0.05s CPU 0.06s WALL ( 9 calls)
vdW_kernel : 25.24s CPU 25.25s WALL ( 1 calls)
Called by c_bands:
init_us_2 : 0.01s CPU 0.01s WALL ( 19 calls)
regterg : 0.19s CPU 0.19s WALL ( 9 calls)
Called by sum_band:
sum_band:bec : 0.00s CPU 0.00s WALL ( 9 calls)
addusdens : 0.12s CPU 0.12s WALL ( 9 calls)
Called by *egterg:
h_psi : 0.17s CPU 0.18s WALL ( 28 calls)
s_psi : 0.00s CPU 0.00s WALL ( 28 calls)
g_psi : 0.00s CPU 0.00s WALL ( 18 calls)
rdiaghg : 0.01s CPU 0.01s WALL ( 27 calls)
regterg:over : 0.00s CPU 0.00s WALL ( 18 calls)
regterg:upda : 0.00s CPU 0.00s WALL ( 18 calls)
regterg:last : 0.00s CPU 0.00s WALL ( 9 calls)
Called by h_psi:
h_psi:calbec : 0.01s CPU 0.01s WALL ( 28 calls)
vloc_psi : 0.16s CPU 0.17s WALL ( 28 calls)
add_vuspsi : 0.00s CPU 0.00s WALL ( 28 calls)
General routines
calbec : 0.01s CPU 0.01s WALL ( 42 calls)
fft : 1.39s CPU 1.47s WALL ( 659 calls)
ffts : 0.02s CPU 0.02s WALL ( 19 calls)
fftw : 0.16s CPU 0.17s WALL ( 256 calls)
interpolate : 0.04s CPU 0.04s WALL ( 10 calls)
Parallel routines
fft_scatt_xy : 0.17s CPU 0.17s WALL ( 934 calls)
fft_scatt_yz : 0.56s CPU 0.61s WALL ( 934 calls)
PWSCF : 30.79s CPU 31.25s WALL
This run was terminated on: 12:40:30 23Sep2019
=------------------------------------------------------------------------------=
JOB DONE.
=------------------------------------------------------------------------------=