Program PWSCF v.> 4.2 starts on 25Jan2011 at 15:38:33 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/wiki/index.php/Citing_Quantum-ESPRESSO Parallel version (MPI), running on 4 processors R & G space division: proc/pool = 4 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 Waiting for input... XC functional enforced from input : Exchange-correlation = VDW-DF (1449) EXX-fraction = 0.00 !!! 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 Stick Mesh ---------- nst = 1539, nstw = 258, nsts = 1025 n.st n.stw n.sts n.g n.gw n.gs min 768 127 512 40311 2739 21932 max 770 130 513 40318 2746 21968 3077 515 2049 161263 10971 87777 bravais-lattice index = 8 lattice parameter (a_0) = 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 (1449) EXX-fraction = 0.00 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 a_0) 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/a_0) 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 O.pbe-rrkjus.UPF MD5 check sum: 390ba29e75625707450f3bd3f0eb6be9 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 H.pbe-rrkjus.UPF MD5 check sum: 7cc9d459525c9a0585f487a71c3c9563 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 vdW kernel table read from file vdW_kernel_table MD5 check sum: fe4853e4b29e331a1c05f2446fb42441 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) s frac. trans. isym = 1 identity cryst. s( 1) = ( 1 0 0 ) f =( 0.0000000 ) ( 0 1 0 ) ( 0.0000000 ) ( 0 0 1 ) ( 0.0000000 ) cart. s( 1) = ( 1.0000000 0.0000000 0.0000000 ) f =( 0.0000000 ) ( 0.0000000 1.0000000 0.0000000 ) ( 0.0000000 ) ( 0.0000000 0.0000000 1.0000000 ) ( 0.0000000 ) isym = 2 inv. 180 deg rotation - cart. axis [1,0,0] cryst. s( 2) = ( -1 0 0 ) f =( 0.0000000 ) ( 0 1 0 ) ( 0.0000000 ) ( 0 0 1 ) ( 0.0000000 ) cart. s( 2) = ( -1.0000000 0.0000000 0.0000000 ) f =( 0.0000000 ) ( 0.0000000 1.0000000 0.0000000 ) ( 0.0000000 ) ( 0.0000000 0.0000000 1.0000000 ) ( 0.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: E 1 s 2 Cartesian axes site n. atom positions (a_0 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/a_0 k( 1) = ( 0.0000000 0.0000000 0.0000000), wk = 2.0000000 cryst. coord. k( 1) = ( 0.0000000 0.0000000 0.0000000), wk = 2.0000000 G cutoff = 1025.8770 ( 80632 G-vectors) FFT grid: ( 72, 64, 80) G cutoff = 683.9180 ( 43889 G-vectors) smooth grid: ( 54, 50, 72) Largest allocated arrays est. size (Mb) dimensions Kohn-Sham Wavefunctions 0.17 Mb ( 1373, 8) NL pseudopotentials 0.50 Mb ( 1373, 24) Each V/rho on FFT grid 1.41 Mb ( 92160) Each G-vector array 0.15 Mb ( 20158) G-vector shells 0.07 Mb ( 9015) Largest temporary arrays est. size (Mb) dimensions Auxiliary wavefunctions 0.34 Mb ( 1373, 32) Each subspace H/S matrix 0.01 Mb ( 32, 32) Each matrix 0.00 Mb ( 24, 8) Arrays for rho mixing 11.25 Mb ( 92160, 8) Initial potential from superposition of free atoms starting charge 15.61518, renormalised to 16.00000 negative rho (up, down): 0.281E-04 0.000E+00 --------------------------------------------------------------------------------- Carrying out vdW-DF run using the following parameters: Nqs = 20 Nr_points = 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 correlation energy = 0.274904696116047 ---------------------------------------------------------------- Starting wfc are 12 atomic wfcs total cpu time spent up to now is 3.83 secs per-process dynamical memory: 43.6 Mb 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): 0.249E-04 0.000E+00 ---------------------------------------------------------------- Non-local correlation energy = 0.301859606238180 ---------------------------------------------------------------- total cpu time spent up to now is 6.04 secs total energy = -68.62054345 Ry Harris-Foulkes estimate = -69.74338030 Ry estimated scf accuracy < 1.46176767 Ry iteration # 2 ecut= 30.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 9.14E-03, avg # of iterations = 2.0 negative rho (up, down): 0.806E-04 0.000E+00 ---------------------------------------------------------------- Non-local correlation energy = 0.299176558067486 ---------------------------------------------------------------- total cpu time spent up to now is 8.21 secs total energy = -68.83477871 Ry Harris-Foulkes estimate = -69.33634384 Ry estimated scf accuracy < 0.97176223 Ry iteration # 3 ecut= 30.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 6.07E-03, avg # of iterations = 2.0 negative rho (up, down): 0.687E-03 0.000E+00 ---------------------------------------------------------------- Non-local correlation energy = 0.304378686108314 ---------------------------------------------------------------- total cpu time spent up to now is 10.30 secs total energy = -69.04855398 Ry Harris-Foulkes estimate = -69.06844693 Ry estimated scf accuracy < 0.03452251 Ry iteration # 4 ecut= 30.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 2.16E-04, avg # of iterations = 2.0 negative rho (up, down): 0.639E-03 0.000E+00 ---------------------------------------------------------------- Non-local correlation energy = 0.303938430184660 ---------------------------------------------------------------- total cpu time spent up to now is 12.39 secs total energy = -69.05564659 Ry Harris-Foulkes estimate = -69.05601332 Ry estimated scf accuracy < 0.00065530 Ry iteration # 5 ecut= 30.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 4.10E-06, avg # of iterations = 2.0 negative rho (up, down): 0.177E-02 0.000E+00 ---------------------------------------------------------------- Non-local correlation energy = 0.304125165708814 ---------------------------------------------------------------- total cpu time spent up to now is 14.47 secs total energy = -69.05584470 Ry Harris-Foulkes estimate = -69.05582135 Ry estimated scf accuracy < 0.00002804 Ry iteration # 6 ecut= 30.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.75E-07, avg # of iterations = 2.0 negative rho (up, down): 0.187E-02 0.000E+00 ---------------------------------------------------------------- Non-local correlation energy = 0.304109872827670 ---------------------------------------------------------------- total cpu time spent up to now is 16.53 secs total energy = -69.05584769 Ry Harris-Foulkes estimate = -69.05584984 Ry estimated scf accuracy < 0.00000254 Ry iteration # 7 ecut= 30.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.59E-08, avg # of iterations = 2.0 negative rho (up, down): 0.189E-02 0.000E+00 ---------------------------------------------------------------- Non-local correlation energy = 0.304110383940431 ---------------------------------------------------------------- total cpu time spent up to now is 18.59 secs total energy = -69.05584799 Ry Harris-Foulkes estimate = -69.05584814 Ry estimated scf accuracy < 0.00000005 Ry iteration # 8 ecut= 30.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 2.93E-10, avg # of iterations = 2.0 negative rho (up, down): 0.189E-02 0.000E+00 ---------------------------------------------------------------- Non-local correlation energy = 0.304109884755319 ---------------------------------------------------------------- total cpu time spent up to now is 20.65 secs total energy = -69.05584798 Ry Harris-Foulkes estimate = -69.05584800 Ry estimated scf accuracy < 0.00000001 Ry iteration # 9 ecut= 30.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 6.39E-11, avg # of iterations = 2.0 negative rho (up, down): 0.189E-02 0.000E+00 ---------------------------------------------------------------- Non-local correlation energy = 0.304109294756646 ---------------------------------------------------------------- total cpu time spent up to now is 22.61 secs End of self-consistent calculation k = 0.0000 0.0000 0.0000 ( 5486 PWs) bands (ev): -25.6442 -24.2676 -13.5673 -12.2979 -9.7601 -8.3492 -7.6831 -6.4201 occupation numbers 2.0000 2.0000 2.0000 2.0000 2.0000 2.0000 2.0000 2.0000 ! total energy = -69.05584800 Ry Harris-Foulkes estimate = -69.05584798 Ry estimated scf accuracy < 7.8E-10 Ry The total energy is the sum of the following terms: one-electron contribution = -122.46795392 Ry hartree contribution = 64.43483108 Ry xc contribution = -17.37518989 Ry ewald contribution = 6.35246473 Ry convergence has been achieved in 9 iterations Forces acting on atoms (Ry/au): atom 1 type 1 force = 0.00000000 0.00094057 -0.00216126 atom 2 type 2 force = 0.00000000 0.00013998 0.00068311 atom 3 type 2 force = 0.00000000 -0.00035102 0.00014004 atom 4 type 1 force = 0.00000000 0.00088157 0.00119056 atom 5 type 2 force = 0.00002030 -0.00080555 0.00007378 atom 6 type 2 force = -0.00002030 -0.00080555 0.00007378 Total force = 0.003113 Total SCF correction = 0.000038 entering subroutine stress ... VDW GRADIENT stress 0.00000366 0.00000000 0.00000000 0.00000000 0.00000346 0.00000000 0.00000000 -0.00000002 0.00000315 VDW KERNEL stress -0.00002787 0.00000000 0.00000000 0.00000000 -0.00002821 0.00000000 0.00000000 -0.00000003 -0.00002322 VDW ALL stress 0.00002421 0.00000000 0.00000000 0.00000000 0.00002475 0.00000004 0.00000000 0.00000004 0.00002007 total stress (Ry/bohr**3) (kbar) P= -1.07 -0.00000740 0.00000000 0.00000000 -1.09 0.00 0.00 0.00000000 -0.00000655 -0.00000048 0.00 -0.96 -0.07 0.00000000 -0.00000048 -0.00000784 0.00 -0.07 -1.15 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.33 0.00 0.00 0.00 280.61 11.21 0.00 11.21 269.14 hartree stress (kbar) 584.70 0.00 0.00 0.00 685.80 102.02 0.00 102.02 1126.88 exc-cor stress (kbar) -207.39 0.00 0.00 0.00 -207.96 -1.03 0.00 -1.03 -205.77 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-nl stress (kbar) 3.56 0.00 0.00 0.00 3.64 0.01 0.00 0.01 2.95 Writing output data file water_vdw.save init_run : 3.38s CPU 3.47s WALL ( 1 calls) electrons : 18.42s CPU 18.77s WALL ( 1 calls) forces : 0.92s CPU 0.92s WALL ( 1 calls) stress : 3.55s CPU 3.60s WALL ( 1 calls) Called by init_run: wfcinit : 0.06s CPU 0.07s WALL ( 1 calls) potinit : 2.29s CPU 2.36s WALL ( 1 calls) Called by electrons: c_bands : 1.09s CPU 1.09s WALL ( 9 calls) sum_band : 1.57s CPU 1.58s WALL ( 9 calls) v_of_rho : 15.52s CPU 15.68s WALL ( 10 calls) v_h : 0.20s CPU 0.20s WALL ( 10 calls) v_xc : 15.32s CPU 15.47s WALL ( 10 calls) newd : 0.98s CPU 0.99s WALL ( 10 calls) mix_rho : 0.37s CPU 0.37s WALL ( 9 calls) vdW_energy : 1.76s CPU 1.76s WALL ( 10 calls) vdW_ffts : 7.28s CPU 7.32s WALL ( 22 calls) vdW_v : 2.87s CPU 2.88s WALL ( 10 calls) Called by c_bands: init_us_2 : 0.02s CPU 0.02s WALL ( 19 calls) regterg : 1.06s CPU 1.06s WALL ( 9 calls) Called by sum_band: sum_band:bec : 0.00s CPU 0.00s WALL ( 9 calls) addusdens : 0.74s CPU 0.74s WALL ( 9 calls) Called by *egterg: h_psi : 1.02s CPU 1.02s WALL ( 28 calls) s_psi : 0.01s CPU 0.02s WALL ( 28 calls) g_psi : 0.00s CPU 0.01s WALL ( 18 calls) rdiaghg : 0.01s CPU 0.01s WALL ( 27 calls) regterg:over : 0.02s CPU 0.02s WALL ( 18 calls) regterg:upda : 0.02s CPU 0.01s WALL ( 18 calls) regterg:last : 0.00s CPU 0.01s WALL ( 9 calls) Called by h_psi: h_psi:vloc : 0.97s CPU 0.97s WALL ( 28 calls) h_psi:vnl : 0.04s CPU 0.05s WALL ( 28 calls) add_vuspsi : 0.02s CPU 0.02s WALL ( 28 calls) General routines calbec : 0.04s CPU 0.04s WALL ( 42 calls) fft : 9.98s CPU 10.05s WALL ( 605 calls) ffts : 0.14s CPU 0.15s WALL ( 19 calls) fftw : 0.97s CPU 0.98s WALL ( 254 calls) interpolate : 0.50s CPU 0.51s WALL ( 19 calls) davcio : 0.00s CPU 0.01s WALL ( 9 calls) Parallel routines fft_scatter : 3.23s CPU 3.27s WALL ( 878 calls) PWSCF : 26.67s CPU 27.31s WALL This run was terminated on: 15:39: 1 25Jan2011 =------------------------------------------------------------------------------= JOB DONE. =------------------------------------------------------------------------------=