mirror of https://gitlab.com/QEF/q-e.git
319 lines
14 KiB
Plaintext
319 lines
14 KiB
Plaintext
|
|
||
|
Program PWSCF v.5.3.0 (svn rev. 11957) starts on 8Jan2016 at 12:28:52
|
||
|
|
||
|
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
|
||
|
|
||
|
Serial version
|
||
|
Reading input from xdm.in
|
||
|
|
||
|
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
|
||
|
file C.pbe-n-kjpaw_psl.0.1.UPF: wavefunction(s) 2P renormalized
|
||
|
|
||
|
gamma-point specific algorithms are used
|
||
|
|
||
|
G-vector sticks info
|
||
|
--------------------
|
||
|
sticks: dense smooth PW G-vecs: dense smooth PW
|
||
|
Sum 301 109 31 10915 2349 287
|
||
|
Tot 151 55 16
|
||
|
|
||
|
|
||
|
|
||
|
bravais-lattice index = 4
|
||
|
lattice parameter (alat) = 4.6600 a.u.
|
||
|
unit-cell volume = 227.8567 (a.u.)^3
|
||
|
number of atoms/cell = 4
|
||
|
number of atomic types = 1
|
||
|
number of electrons = 16.00
|
||
|
number of Kohn-Sham states= 12
|
||
|
kinetic-energy cutoff = 18.0000 Ry
|
||
|
charge density cutoff = 200.0000 Ry
|
||
|
convergence threshold = 1.0E-06
|
||
|
mixing beta = 0.5000
|
||
|
number of iterations used = 20 plain mixing
|
||
|
Exchange-correlation = SLA PW PBX PBC ( 1 4 3 4 0 0)
|
||
|
|
||
|
celldm(1)= 4.660000 celldm(2)= 0.000000 celldm(3)= 2.600000
|
||
|
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.500000 0.866025 0.000000 )
|
||
|
a(3) = ( 0.000000 0.000000 2.600000 )
|
||
|
|
||
|
reciprocal axes: (cart. coord. in units 2 pi/alat)
|
||
|
b(1) = ( 1.000000 0.577350 -0.000000 )
|
||
|
b(2) = ( 0.000000 1.154701 0.000000 )
|
||
|
b(3) = ( 0.000000 -0.000000 0.384615 )
|
||
|
|
||
|
|
||
|
PseudoPot. # 1 for C read from file:
|
||
|
/home/fs395/scratch/QE/espresso/test-suite/..//pseudo/C.pbe-n-kjpaw_psl.0.1.UPF
|
||
|
MD5 check sum: c849119b86527ef673bbe9c4f1f1b6ab
|
||
|
Pseudo is Projector augmented-wave + core cor, Zval = 4.0
|
||
|
Generated using "atomic" code by A. Dal Corso v.5.0.2 svn rev. 9415
|
||
|
Shape of augmentation charge: BESSEL
|
||
|
Using radial grid of 1073 points, 4 beta functions with:
|
||
|
l(1) = 0
|
||
|
l(2) = 0
|
||
|
l(3) = 1
|
||
|
l(4) = 1
|
||
|
Q(r) pseudized with 0 coefficients
|
||
|
|
||
|
|
||
|
atomic species valence mass pseudopotential
|
||
|
C 4.00 12.00000 C( 1.00)
|
||
|
|
||
|
8 Sym. Ops., with inversion, found ( 4 have fractional translation)
|
||
|
|
||
|
|
||
|
|
||
|
Cartesian axes
|
||
|
|
||
|
site n. atom positions (alat units)
|
||
|
1 C tau( 1) = ( -0.5000000 0.8660254 1.9500000 )
|
||
|
2 C tau( 2) = ( 0.5000050 0.2886722 1.9500000 )
|
||
|
3 C tau( 3) = ( -0.5000000 0.8660254 0.6500000 )
|
||
|
4 C tau( 4) = ( -0.0000050 0.5773532 0.6500000 )
|
||
|
|
||
|
number of k points= 1 Marzari-Vanderbilt smearing, width (Ry)= 0.0200
|
||
|
cart. coord. in units 2pi/alat
|
||
|
k( 1) = ( 0.0000000 0.0000000 0.0000000), wk = 2.0000000
|
||
|
|
||
|
Dense grid: 5458 G-vectors FFT dimensions: ( 24, 24, 60)
|
||
|
|
||
|
Smooth grid: 1175 G-vectors FFT dimensions: ( 15, 15, 36)
|
||
|
|
||
|
Largest allocated arrays est. size (Mb) dimensions
|
||
|
Kohn-Sham Wavefunctions 0.03 Mb ( 144, 12)
|
||
|
NL pseudopotentials 0.07 Mb ( 144, 32)
|
||
|
Each V/rho on FFT grid 0.53 Mb ( 34560)
|
||
|
Each G-vector array 0.04 Mb ( 5458)
|
||
|
G-vector shells 0.00 Mb ( 616)
|
||
|
Largest temporary arrays est. size (Mb) dimensions
|
||
|
Auxiliary wavefunctions 0.05 Mb ( 144, 48)
|
||
|
Each subspace H/S matrix 0.02 Mb ( 48, 48)
|
||
|
Each <psi_i|beta_j> matrix 0.00 Mb ( 32, 12)
|
||
|
Arrays for rho mixing 10.55 Mb ( 34560, 20)
|
||
|
|
||
|
Check: negative/imaginary core charge= -0.000011 0.000000
|
||
|
|
||
|
Initial potential from superposition of free atoms
|
||
|
|
||
|
starting charge 15.99983, renormalised to 16.00000
|
||
|
Starting wfc are 16 randomized atomic wfcs
|
||
|
|
||
|
total cpu time spent up to now is 0.7 secs
|
||
|
|
||
|
per-process dynamical memory: 6.7 Mb
|
||
|
|
||
|
Self-consistent Calculation
|
||
|
|
||
|
iteration # 1 ecut= 18.00 Ry beta=0.50
|
||
|
Davidson diagonalization with overlap
|
||
|
ethr = 1.00E-02, avg # of iterations = 3.0
|
||
|
|
||
|
total cpu time spent up to now is 1.0 secs
|
||
|
|
||
|
total energy = -70.75696246 Ry
|
||
|
Harris-Foulkes estimate = -70.93981621 Ry
|
||
|
estimated scf accuracy < 0.53772235 Ry
|
||
|
|
||
|
iteration # 2 ecut= 18.00 Ry beta=0.50
|
||
|
Davidson diagonalization with overlap
|
||
|
ethr = 3.36E-03, avg # of iterations = 1.0
|
||
|
|
||
|
total cpu time spent up to now is 1.3 secs
|
||
|
|
||
|
total energy = -70.74797408 Ry
|
||
|
Harris-Foulkes estimate = -70.77687693 Ry
|
||
|
estimated scf accuracy < 0.08906399 Ry
|
||
|
|
||
|
iteration # 3 ecut= 18.00 Ry beta=0.50
|
||
|
Davidson diagonalization with overlap
|
||
|
ethr = 5.57E-04, avg # of iterations = 2.0
|
||
|
|
||
|
total cpu time spent up to now is 1.6 secs
|
||
|
|
||
|
total energy = -70.75427887 Ry
|
||
|
Harris-Foulkes estimate = -70.75473766 Ry
|
||
|
estimated scf accuracy < 0.00241122 Ry
|
||
|
|
||
|
iteration # 4 ecut= 18.00 Ry beta=0.50
|
||
|
Davidson diagonalization with overlap
|
||
|
ethr = 1.51E-05, avg # of iterations = 2.0
|
||
|
|
||
|
total cpu time spent up to now is 1.8 secs
|
||
|
|
||
|
total energy = -70.75443837 Ry
|
||
|
Harris-Foulkes estimate = -70.75445713 Ry
|
||
|
estimated scf accuracy < 0.00003810 Ry
|
||
|
|
||
|
iteration # 5 ecut= 18.00 Ry beta=0.50
|
||
|
Davidson diagonalization with overlap
|
||
|
ethr = 2.38E-07, avg # of iterations = 9.0
|
||
|
|
||
|
total cpu time spent up to now is 2.1 secs
|
||
|
|
||
|
total energy = -70.74792325 Ry
|
||
|
Harris-Foulkes estimate = -70.74792324 Ry
|
||
|
estimated scf accuracy < 0.00000142 Ry
|
||
|
|
||
|
iteration # 6 ecut= 18.00 Ry beta=0.50
|
||
|
Davidson diagonalization with overlap
|
||
|
ethr = 8.85E-09, avg # of iterations = 2.0
|
||
|
|
||
|
total cpu time spent up to now is 2.4 secs
|
||
|
|
||
|
End of self-consistent calculation
|
||
|
|
||
|
k = 0.0000 0.0000 0.0000 ( 144 PWs) bands (ev):
|
||
|
|
||
|
-12.2285 -11.8060 -0.1530 1.9646 6.9306 6.9308 7.9258 7.9259
|
||
|
12.2469 18.3083 18.3089 18.9456
|
||
|
|
||
|
the Fermi energy is 8.1765 ev
|
||
|
|
||
|
* XDM dispersion
|
||
|
a1 = 0.327500
|
||
|
a2 (ang) = 2.767300
|
||
|
a2 (bohr) = 5.229439
|
||
|
|
||
|
+ Volumes and moments
|
||
|
# All results in atomic units (Hartree,bohr)
|
||
|
# i V Vfree M1 M2 M3
|
||
|
1 3.502570E+01 3.899308E+01 3.018403E+00 3.883027E+01 7.140367E+02
|
||
|
2 3.502051E+01 3.899308E+01 3.030460E+00 3.886995E+01 7.135457E+02
|
||
|
3 3.502564E+01 3.899308E+01 3.018407E+00 3.883029E+01 7.140367E+02
|
||
|
4 3.502047E+01 3.899308E+01 3.030458E+00 3.886997E+01 7.135459E+02
|
||
|
|
||
|
+ Dispersion coefficients
|
||
|
# All results in atomic units (Hartree,bohr).
|
||
|
# i j C6 C8 C10 Rc Rvdw
|
||
|
1 1 1.610111E+01 6.213985E+02 2.642716E+04 6.366254E+00 7.314387E+00
|
||
|
2 1 1.613200E+01 6.216692E+02 2.640908E+04 6.362124E+00 7.313035E+00
|
||
|
2 2 1.616303E+01 6.219413E+02 2.639104E+04 6.357992E+00 7.311682E+00
|
||
|
3 1 1.610111E+01 6.213981E+02 2.642714E+04 6.366253E+00 7.314387E+00
|
||
|
3 2 1.613200E+01 6.216689E+02 2.640906E+04 6.362123E+00 7.313034E+00
|
||
|
3 3 1.610111E+01 6.213978E+02 2.642712E+04 6.366252E+00 7.314387E+00
|
||
|
4 1 1.613199E+01 6.216690E+02 2.640907E+04 6.362126E+00 7.313035E+00
|
||
|
4 2 1.616301E+01 6.219411E+02 2.639103E+04 6.357994E+00 7.311682E+00
|
||
|
4 3 1.613199E+01 6.216686E+02 2.640905E+04 6.362124E+00 7.313035E+00
|
||
|
4 4 1.616300E+01 6.219408E+02 2.639103E+04 6.357995E+00 7.311682E+00
|
||
|
|
||
|
+ van der Waals energies, forces and stresses (Ry,bohr)
|
||
|
Evdw(total,Ry) = -3.716185694926E-02
|
||
|
Evdw(C6,Ry) = -1.857211334124E-02
|
||
|
Evdw(C8,Ry) = -1.079408362252E-02
|
||
|
Evdw(C10,Ry) = -7.795659985495E-03
|
||
|
Fvdw (001,Ry/bohr) = 1.076050517253E-09 -6.212580410888E-10 3.113826143617E-19
|
||
|
Fvdw (002,Ry/bohr) = -8.314769489572E-10 4.800534384154E-10 2.639110339744E-19
|
||
|
Fvdw (003,Ry/bohr) = -1.008476746592E-09 5.822443125410E-10 -2.149648891964E-20
|
||
|
Fvdw (004,Ry/bohr) = 7.639031773011E-10 -4.410397067912E-10 1.241554582658E-19
|
||
|
sigma_vdw (Ry/bohr**3) = -1.671678619648E-04 3.653079666431E-11 6.819010926240E-22
|
||
|
3.653079666431E-11 -1.671678197827E-04 2.177868983706E-22
|
||
|
6.819010926240E-22 2.177868983706E-22 -1.954498798977E-04
|
||
|
sigma_vdw (GPa) = -2.459123736451E+00 5.373864816645E-07 1.003110970658E-17
|
||
|
5.373864816645E-07 -2.459123115930E+00 3.203755344935E-18
|
||
|
1.003110970658E-17 3.203755344935E-18 -2.875166514028E+00
|
||
|
|
||
|
|
||
|
! total energy = -70.78508555 Ry
|
||
|
Harris-Foulkes estimate = -70.78508533 Ry
|
||
|
estimated scf accuracy < 0.00000024 Ry
|
||
|
|
||
|
total all-electron energy = -301.841492 Ry
|
||
|
|
||
|
The total energy is the sum of the following terms:
|
||
|
|
||
|
one-electron contribution = -4.01450216 Ry
|
||
|
hartree contribution = 12.15493336 Ry
|
||
|
xc contribution = -16.66359122 Ry
|
||
|
ewald contribution = -35.87244982 Ry
|
||
|
Dispersion XDM Correction = -0.03716186 Ry
|
||
|
one-center paw contrib. = -26.35883785 Ry
|
||
|
smearing contrib. (-TS) = 0.00652398 Ry
|
||
|
|
||
|
convergence has been achieved in 6 iterations
|
||
|
|
||
|
Forces acting on atoms (Ry/au):
|
||
|
|
||
|
atom 1 type 1 force = 0.00002656 -0.00001534 0.00000000
|
||
|
atom 2 type 1 force = -0.00002956 0.00001707 0.00000000
|
||
|
atom 3 type 1 force = -0.00002656 0.00001534 0.00000000
|
||
|
atom 4 type 1 force = 0.00002956 -0.00001707 0.00000000
|
||
|
|
||
|
Total force = 0.000065 Total SCF correction = 0.000080
|
||
|
SCF correction compared to forces is large: reduce conv_thr to get better values
|
||
|
|
||
|
|
||
|
entering subroutine stress ...
|
||
|
|
||
|
total stress (Ry/bohr**3) (kbar) P=-2730.14
|
||
|
-0.01813130 -0.00000035 0.00000000 -2667.21 -0.05 0.00
|
||
|
-0.00000035 -0.01813170 0.00000000 -0.05 -2667.26 0.00
|
||
|
0.00000000 0.00000000 -0.01941433 0.00 0.00 -2855.95
|
||
|
|
||
|
|
||
|
Writing output data file pwscf.save
|
||
|
|
||
|
init_run : 0.51s CPU 0.53s WALL ( 1 calls)
|
||
|
electrons : 4.81s CPU 4.82s WALL ( 1 calls)
|
||
|
forces : 0.09s CPU 0.09s WALL ( 1 calls)
|
||
|
stress : 0.18s CPU 0.18s WALL ( 1 calls)
|
||
|
|
||
|
Called by init_run:
|
||
|
wfcinit : 0.00s CPU 0.00s WALL ( 1 calls)
|
||
|
potinit : 0.27s CPU 0.28s WALL ( 1 calls)
|
||
|
|
||
|
Called by electrons:
|
||
|
c_bands : 0.04s CPU 0.04s WALL ( 6 calls)
|
||
|
sum_band : 0.08s CPU 0.09s WALL ( 6 calls)
|
||
|
v_of_rho : 0.21s CPU 0.21s WALL ( 7 calls)
|
||
|
newd : 0.03s CPU 0.04s WALL ( 7 calls)
|
||
|
PAW_pot : 1.56s CPU 1.58s WALL ( 7 calls)
|
||
|
mix_rho : 0.02s CPU 0.02s WALL ( 6 calls)
|
||
|
|
||
|
Called by c_bands:
|
||
|
init_us_2 : 0.00s CPU 0.00s WALL ( 13 calls)
|
||
|
regterg : 0.03s CPU 0.03s WALL ( 6 calls)
|
||
|
|
||
|
Called by sum_band:
|
||
|
sum_band:bec : 0.00s CPU 0.00s WALL ( 6 calls)
|
||
|
addusdens : 0.03s CPU 0.03s WALL ( 6 calls)
|
||
|
|
||
|
Called by *egterg:
|
||
|
h_psi : 0.02s CPU 0.02s WALL ( 26 calls)
|
||
|
s_psi : 0.00s CPU 0.00s WALL ( 26 calls)
|
||
|
g_psi : 0.00s CPU 0.00s WALL ( 19 calls)
|
||
|
rdiaghg : 0.01s CPU 0.01s WALL ( 25 calls)
|
||
|
|
||
|
Called by h_psi:
|
||
|
add_vuspsi : 0.00s CPU 0.00s WALL ( 26 calls)
|
||
|
|
||
|
General routines
|
||
|
calbec : 0.00s CPU 0.00s WALL ( 37 calls)
|
||
|
fft : 0.15s CPU 0.15s WALL ( 171 calls)
|
||
|
ffts : 0.00s CPU 0.00s WALL ( 19 calls)
|
||
|
fftw : 0.03s CPU 0.03s WALL ( 350 calls)
|
||
|
interpolate : 0.02s CPU 0.02s WALL ( 19 calls)
|
||
|
davcio : 0.00s CPU 0.00s WALL ( 1 calls)
|
||
|
|
||
|
|
||
|
PWSCF : 5.77s CPU 5.83s WALL
|
||
|
|
||
|
|
||
|
This run was terminated on: 12:28:58 8Jan2016
|
||
|
|
||
|
=------------------------------------------------------------------------------=
|
||
|
JOB DONE.
|
||
|
=------------------------------------------------------------------------------=
|