quantum-espresso/PW/examples/example04/reference/BP.out

439 lines
19 KiB
Plaintext

Program PWSCF v.6.0 (svn rev. 13188) starts on 6Dec2016 at 15: 0: 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), running on 4 processors
R & G space division: proc/nbgrp/npool/nimage = 4
Waiting for input...
Reading input from standard input
Message from routine read_cards :
DEPRECATED: no units specified in ATOMIC_POSITIONS card
Message from routine read_cards :
ATOMIC_POSITIONS: units set to alat
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
Atomic positions and unit cell read from directory:
/home/pietro/espresso-svn/tempdir/pwscf.save/
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 129 129 40 2231 2231 386
Max 130 130 41 2232 2232 389
Sum 517 517 161 8925 8925 1551
bravais-lattice index = 1
lattice parameter (alat) = 7.3699 a.u.
unit-cell volume = 400.2993 (a.u.)^3
number of atoms/cell = 5
number of atomic types = 3
number of electrons = 44.00
number of Kohn-Sham states= 22
kinetic-energy cutoff = 30.0000 Ry
charge density cutoff = 120.0000 Ry
Exchange-correlation = SLA PZ NOGX NOGC ( 1 1 0 0 0 0)
celldm(1)= 7.369900 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.000000 )
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 1.000000 )
PseudoPot. # 1 for Pb read from file:
/home/pietro/espresso-svn/pseudo/Pb.pz-d-van.UPF
MD5 check sum: 4e1e5920686a026ae26139ac417581ff
Pseudo is Ultrasoft, Zval = 14.0
Generated by new atomic code, or converted to UPF format
Using radial grid of 899 points, 6 beta functions with:
l(1) = 0
l(2) = 0
l(3) = 1
l(4) = 1
l(5) = 2
l(6) = 2
Q(r) pseudized with 8 coefficients, rinner = 1.000 1.000 1.000
1.000 1.000
PseudoPot. # 2 for Ti read from file:
/home/pietro/espresso-svn/pseudo/Ti.pz-sp-van_ak.UPF
MD5 check sum: 545d0e6e05332b8871a8093f427cb0ca
Pseudo is Ultrasoft, Zval = 12.0
Generated by new atomic code, or converted to UPF format
Using radial grid of 851 points, 6 beta functions with:
l(1) = 0
l(2) = 0
l(3) = 1
l(4) = 1
l(5) = 2
l(6) = 2
Q(r) pseudized with 8 coefficients, rinner = 1.000 1.000 1.000
1.000 1.000
PseudoPot. # 3 for O read from file:
/home/pietro/espresso-svn/pseudo/O.pz-van_ak.UPF
MD5 check sum: d814fcb982dd9af4fc6452aae6bb9318
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
Pb 14.00 207.20000 Pb( 1.00)
Ti 12.00 47.86700 Ti( 1.00)
O 6.00 15.99940 O ( 1.00)
No symmetry found
Cartesian axes
site n. atom positions (alat units)
1 Pb tau( 1) = ( 0.0000000 0.0000000 0.0100000 )
2 Ti tau( 2) = ( 0.5000000 0.5000000 0.5000000 )
3 O tau( 3) = ( 0.0000000 0.5000000 0.5000000 )
4 O tau( 4) = ( 0.5000000 0.5000000 0.0000000 )
5 O tau( 5) = ( 0.5000000 0.0000000 0.5000000 )
number of k points= 21
cart. coord. in units 2pi/alat
k( 1) = ( 0.1250000 0.1250000 -0.5000000), wk = 0.0714286
k( 2) = ( 0.1250000 0.1250000 -0.3333333), wk = 0.0714286
k( 3) = ( 0.1250000 0.1250000 -0.1666667), wk = 0.0714286
k( 4) = ( 0.1250000 0.1250000 0.0000000), wk = 0.0714286
k( 5) = ( 0.1250000 0.1250000 0.1666667), wk = 0.0714286
k( 6) = ( 0.1250000 0.1250000 0.3333333), wk = 0.0714286
k( 7) = ( 0.1250000 0.1250000 0.5000000), wk = 0.0714286
k( 8) = ( 0.1250000 0.3750000 -0.5000000), wk = 0.1428571
k( 9) = ( 0.1250000 0.3750000 -0.3333333), wk = 0.1428571
k( 10) = ( 0.1250000 0.3750000 -0.1666667), wk = 0.1428571
k( 11) = ( 0.1250000 0.3750000 0.0000000), wk = 0.1428571
k( 12) = ( 0.1250000 0.3750000 0.1666667), wk = 0.1428571
k( 13) = ( 0.1250000 0.3750000 0.3333333), wk = 0.1428571
k( 14) = ( 0.1250000 0.3750000 0.5000000), wk = 0.1428571
k( 15) = ( 0.3750000 0.3750000 -0.5000000), wk = 0.0714286
k( 16) = ( 0.3750000 0.3750000 -0.3333333), wk = 0.0714286
k( 17) = ( 0.3750000 0.3750000 -0.1666667), wk = 0.0714286
k( 18) = ( 0.3750000 0.3750000 0.0000000), wk = 0.0714286
k( 19) = ( 0.3750000 0.3750000 0.1666667), wk = 0.0714286
k( 20) = ( 0.3750000 0.3750000 0.3333333), wk = 0.0714286
k( 21) = ( 0.3750000 0.3750000 0.5000000), wk = 0.0714286
Dense grid: 8925 G-vectors FFT dimensions: ( 25, 25, 25)
Estimated max dynamical RAM per process > 2.00Mb
Estimated total allocated dynamical RAM > 7.99Mb
The potential is recalculated from file :
/home/pietro/espresso-svn/tempdir/pwscf.save/charge-density.dat
Starting wfc are 31 randomized atomic wfcs
Band Structure Calculation
Davidson diagonalization with overlap
ethr = 2.27E-08, avg # of iterations = 10.0
total cpu time spent up to now is 3.5 secs
End of band structure calculation
k = 0.1250 0.1250-0.5000 ( 1116 PWs) bands (ev):
-45.0286 -21.6879 -21.6220 -21.6218 -6.1326 -5.4260 -5.3629 -4.3696
-4.2933 -4.2820 -4.2463 -4.1042 3.5400 6.3080 7.1160 7.7181
8.0885 8.3302 9.1255 9.6741 9.8694 10.7410
k = 0.1250 0.1250-0.3333 ( 1117 PWs) bands (ev):
-45.0290 -21.6692 -21.6228 -21.6226 -6.2744 -5.4698 -5.4253 -4.3523
-4.3007 -4.2873 -4.2070 -4.1049 3.7323 6.7755 7.1083 7.5501
8.2562 8.5571 9.1384 9.6963 9.8243 10.3939
k = 0.1250 0.1250-0.1667 ( 1105 PWs) bands (ev):
-45.0297 -21.6313 -21.6240 -21.6234 -6.5522 -5.5496 -5.5385 -4.3119
-4.3043 -4.2896 -4.1420 -4.1251 4.2259 7.0456 7.0703 7.8187
8.7404 8.7739 9.5773 9.6672 9.7344 9.8378
k = 0.1250 0.1250 0.0000 ( 1109 PWs) bands (ev):
-45.0301 -21.6248 -21.6242 -21.6124 -6.6871 -5.5919 -5.5864 -4.2995
-4.2894 -4.2884 -4.1431 -4.1150 4.5357 7.0441 7.1578 7.5699
9.0082 9.2275 9.4827 9.6521 9.8947 9.9577
k = 0.1250 0.1250 0.1667 ( 1105 PWs) bands (ev):
-45.0297 -21.6313 -21.6240 -21.6234 -6.5522 -5.5496 -5.5385 -4.3119
-4.3043 -4.2896 -4.1420 -4.1251 4.2259 7.0456 7.0703 7.8187
8.7404 8.7739 9.5773 9.6672 9.7344 9.8378
k = 0.1250 0.1250 0.3333 ( 1117 PWs) bands (ev):
-45.0290 -21.6692 -21.6228 -21.6226 -6.2744 -5.4698 -5.4253 -4.3523
-4.3007 -4.2873 -4.2070 -4.1049 3.7323 6.7755 7.1083 7.5501
8.2562 8.5571 9.1384 9.6963 9.8243 10.3939
k = 0.1250 0.1250 0.5000 ( 1116 PWs) bands (ev):
-45.0286 -21.6879 -21.6220 -21.6218 -6.1326 -5.4260 -5.3629 -4.3696
-4.2933 -4.2820 -4.2463 -4.1042 3.5400 6.3080 7.1160 7.7181
8.0885 8.3302 9.1255 9.6741 9.8694 10.7410
k = 0.1250 0.3750-0.5000 ( 1106 PWs) bands (ev):
-45.0274 -21.6848 -21.6740 -21.6198 -5.9903 -5.2766 -5.2247 -4.3318
-4.2775 -4.2413 -4.1747 -3.9800 3.2291 6.3121 6.6225 6.8824
7.1958 8.4624 9.0001 9.3075 9.4565 10.5801
k = 0.1250 0.3750-0.3333 ( 1110 PWs) bands (ev):
-45.0278 -21.6752 -21.6665 -21.6205 -6.0204 -5.3663 -5.2877 -4.3367
-4.3011 -4.2484 -4.1507 -4.0423 3.3576 6.5828 6.7462 6.9480
7.3953 8.3713 9.1406 9.3200 9.8689 10.3673
k = 0.1250 0.3750-0.1667 ( 1121 PWs) bands (ev):
-45.0287 -21.6774 -21.6296 -21.6223 -6.1701 -5.4490 -5.3866 -4.3561
-4.3015 -4.2768 -4.1934 -4.1069 3.6008 6.5934 7.0114 7.4789
8.0946 8.5355 8.9597 9.7687 9.8160 10.4508
k = 0.1250 0.3750 0.0000 ( 1115 PWs) bands (ev):
-45.0290 -21.6782 -21.6229 -21.6104 -6.2800 -5.4492 -5.4234 -4.3579
-4.3244 -4.2862 -4.2529 -4.0982 3.7068 6.5906 7.4780 7.5766
8.1586 8.4586 9.4691 9.5858 9.8375 10.6712
k = 0.1250 0.3750 0.1667 ( 1121 PWs) bands (ev):
-45.0287 -21.6774 -21.6296 -21.6223 -6.1701 -5.4490 -5.3866 -4.3561
-4.3015 -4.2768 -4.1934 -4.1069 3.6008 6.5934 7.0114 7.4789
8.0946 8.5355 8.9597 9.7687 9.8160 10.4508
k = 0.1250 0.3750 0.3333 ( 1110 PWs) bands (ev):
-45.0278 -21.6752 -21.6665 -21.6205 -6.0204 -5.3663 -5.2877 -4.3367
-4.3011 -4.2484 -4.1507 -4.0423 3.3576 6.5828 6.7462 6.9480
7.3953 8.3713 9.1406 9.3200 9.8689 10.3673
k = 0.1250 0.3750 0.5000 ( 1106 PWs) bands (ev):
-45.0274 -21.6848 -21.6740 -21.6198 -5.9903 -5.2766 -5.2247 -4.3318
-4.2775 -4.2413 -4.1747 -3.9800 3.2291 6.3121 6.6225 6.8824
7.1958 8.4624 9.0001 9.3075 9.4565 10.5801
k = 0.3750 0.3750-0.5000 ( 1092 PWs) bands (ev):
-45.0261 -21.6817 -21.6722 -21.6713 -5.5919 -5.4830 -5.4125 -4.2495
-4.2345 -3.9140 -3.9034 -3.8663 4.0549 5.5694 5.7056 6.2366
6.6837 6.7911 7.1829 10.4380 10.5053 10.5822
k = 0.3750 0.3750-0.3333 ( 1103 PWs) bands (ev):
-45.0266 -21.6728 -21.6727 -21.6641 -5.6600 -5.4502 -5.4291 -4.2652
-4.2556 -4.0158 -3.9508 -3.9327 3.8111 5.8290 5.8825 6.5119
6.8020 6.9700 7.9902 10.3259 10.3320 10.4642
k = 0.3750 0.3750-0.1667 ( 1106 PWs) bands (ev):
-45.0274 -21.6746 -21.6741 -21.6272 -5.9395 -5.3601 -5.3003 -4.3245
-4.2805 -4.1979 -4.1259 -4.0051 3.3786 6.5372 6.5421 6.6751
7.0406 8.1607 9.1944 9.5202 9.6783 10.4549
k = 0.3750 0.3750 0.0000 ( 1103 PWs) bands (ev):
-45.0278 -21.6755 -21.6746 -21.6083 -6.0856 -5.2911 -5.2184 -4.3373
-4.3251 -4.3192 -4.1977 -4.0321 3.1966 6.5595 6.6562 7.8839
7.9881 8.0339 8.3931 9.1645 10.0339 10.4563
k = 0.3750 0.3750 0.1667 ( 1106 PWs) bands (ev):
-45.0274 -21.6746 -21.6741 -21.6272 -5.9395 -5.3601 -5.3003 -4.3245
-4.2805 -4.1979 -4.1259 -4.0051 3.3786 6.5372 6.5421 6.6751
7.0406 8.1607 9.1944 9.5202 9.6783 10.4549
k = 0.3750 0.3750 0.3333 ( 1103 PWs) bands (ev):
-45.0266 -21.6728 -21.6727 -21.6641 -5.6600 -5.4502 -5.4291 -4.2652
-4.2556 -4.0158 -3.9508 -3.9327 3.8111 5.8290 5.8825 6.5119
6.8020 6.9700 7.9902 10.3259 10.3320 10.4642
k = 0.3750 0.3750 0.5000 ( 1092 PWs) bands (ev):
-45.0261 -21.6817 -21.6722 -21.6713 -5.5919 -5.4830 -5.4125 -4.2495
-4.2345 -3.9140 -3.9034 -3.8663 4.0549 5.5694 5.7056 6.2366
6.6837 6.7911 7.1829 10.4380 10.5053 10.5822
highest occupied level (ev): 10.7410
==================================================
POLARIZATION CALCULATION
!!! NOT THOROUGHLY TESTED !!!
--------------------------------------------------
K-POINTS STRINGS USED IN CALCULATIONS
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
G-vector along string (2 pi/a): 0.00000 0.00000 1.00000
Modulus of the vector (1/bohr): 0.85255
Number of k-points per string: 7
Number of different strings : 3
IONIC POLARIZATION
~~~~~~~~~~~~~~~~~~
Note: (mod 1) means that the phases (angles ranging from
-pi to pi) have been mapped to the interval [-1/2,+1/2) by
dividing by 2*pi; (mod 2) refers to the interval [-1,+1)
============================================================================
Ion Species Charge Position Phase
----------------------------------------------------------------------------
1 Pb 14.000 0.0000 0.0000 0.0100 0.14000 (mod 2)
2 Ti 12.000 0.5000 0.5000 0.5000 0.00000 (mod 2)
3 O 6.000 0.0000 0.5000 0.5000 -1.00000 (mod 2)
4 O 6.000 0.5000 0.5000 0.0000 0.00000 (mod 2)
5 O 6.000 0.5000 0.0000 0.5000 -1.00000 (mod 2)
----------------------------------------------------------------------------
IONIC PHASE: 0.14000 (mod 2)
============================================================================
ELECTRONIC POLARIZATION
~~~~~~~~~~~~~~~~~~~~~~~
Note: (mod 1) means that the phases (angles ranging from
-pi to pi) have been mapped to the interval [-1/2,+1/2) by
dividing by 2*pi; (mod 2) refers to the interval [-1,+1)
============================================================================
Spin String Weight First k-point in string Phase
----------------------------------------------------------------------------
up 1 0.250000 0.1250 0.1250 -0.5000 -0.05427 (mod 1)
up 2 0.500000 0.1250 0.3750 -0.5000 -0.04874 (mod 1)
up 3 0.250000 0.3750 0.3750 -0.5000 -0.05046 (mod 1)
----------------------------------------------------------------------------
down 1 0.250000 0.1250 0.1250 -0.5000 -0.05427 (mod 1)
down 2 0.500000 0.1250 0.3750 -0.5000 -0.04874 (mod 1)
down 3 0.250000 0.3750 0.3750 -0.5000 -0.05046 (mod 1)
----------------------------------------------------------------------------
Average phase (up): -0.05055 (mod 1)
Average phase (down): -0.05055 (mod 1)
ELECTRONIC PHASE: -0.10111 (mod 2)
============================================================================
SUMMARY OF PHASES
~~~~~~~~~~~~~~~~~
Ionic Phase: 0.14000 (mod 2)
Electronic Phase: -0.10111 (mod 2)
TOTAL PHASE: 0.03889 (mod 2)
VALUES OF POLARIZATION
~~~~~~~~~~~~~~~~~~~~~~
The calculation of phases done along the direction of vector 3
of the reciprocal lattice gives the following contribution to
the polarization vector (in different units, and being Omega
the volume of the unit cell):
P = 0.2866512 (mod 14.7398000) (e/Omega).bohr
P = 0.0007161 (mod 0.0368220) e/bohr^2
P = 0.0409402 (mod 2.1051744) C/m^2
The polarization direction is: ( 0.00000 , 0.00000 , 1.00000 )
==================================================
Writing output data file pwscf.save
init_run : 0.26s CPU 0.31s WALL ( 1 calls)
electrons : 3.10s CPU 3.20s WALL ( 1 calls)
Called by init_run:
wfcinit : 0.00s CPU 0.00s WALL ( 1 calls)
potinit : 0.00s CPU 0.00s WALL ( 1 calls)
Called by electrons:
c_bands : 2.89s CPU 2.97s WALL ( 1 calls)
v_of_rho : 0.00s CPU 0.00s WALL ( 1 calls)
newd : 0.02s CPU 0.02s WALL ( 1 calls)
Called by c_bands:
init_us_2 : 0.02s CPU 0.02s WALL ( 57 calls)
cegterg : 2.57s CPU 2.64s WALL ( 21 calls)
Called by sum_band:
Called by *egterg:
h_psi : 1.12s CPU 1.09s WALL ( 251 calls)
s_psi : 0.20s CPU 0.24s WALL ( 251 calls)
g_psi : 0.00s CPU 0.01s WALL ( 209 calls)
cdiaghg : 0.56s CPU 0.56s WALL ( 230 calls)
Called by h_psi:
h_psi:pot : 1.11s CPU 1.08s WALL ( 251 calls)
h_psi:calbec : 0.23s CPU 0.27s WALL ( 251 calls)
vloc_psi : 0.58s CPU 0.57s WALL ( 251 calls)
add_vuspsi : 0.30s CPU 0.24s WALL ( 251 calls)
General routines
calbec : 0.28s CPU 0.33s WALL ( 287 calls)
fft : 0.00s CPU 0.00s WALL ( 4 calls)
fftw : 0.52s CPU 0.50s WALL ( 8644 calls)
davcio : 0.00s CPU 0.00s WALL ( 57 calls)
Parallel routines
fft_scatter : 0.15s CPU 0.15s WALL ( 8648 calls)
PWSCF : 3.65s CPU 3.83s WALL
This run was terminated on: 15: 0: 5 6Dec2016
=------------------------------------------------------------------------------=
JOB DONE.
=------------------------------------------------------------------------------=