Dear pwforum,<br><br>I am calculating scf for a ibrav=4 structure with espresso-4.3.2.<br>But it can not find the symmetry. In the output file, it report "no symmetry found".<br>I am also calculation with ph.x with the same. the phonon mode is not identified. <br clear="all">
<br>I confirmed the atomic positions in the input file.<br>It is consistent with 32 point symmetry (space group 154). <br><br>I tried ibrav=2 and 12, they can find symmetry and phonon modes, even with <br>atomic positions slightly different from point symmety formula..<br>
<br>So I am wondering: what kind of situation, the code can not find symmtery?<br>Is there any suggestion to make it find symmetry?<br><br>I attached the input file and output file of pw.x followiing.<br><br><br>Thanks very much for you attention and helps.<br>
<br>Best withes.<br><br>WANG Riping<br>2011.8.11 <br><br><br><br><br><br><br><br><br> &CONTROL<br> calculation = scf ,<br> restart_mode = 'from_scratch' ,<br> prefix = alpo4-154.0_0 ,<br> pseudo_dir = '~/pseudo',<br>
outdir = outdir ,<br> nstep = 500 ,<br> max_seconds = 360000 ,<br> etot_conv_thr = 5.0D-6 ,<br> forc_conv_thr = 1.0D10 ,<br> tprnfor = .true.<br> /<br> &system <br> ibrav = 4 ,<br> celldm(1) = 9.471757419 , <br>
celldm(3) = 2.21251394 , <br> nat = 18 , <br> ntyp = 3 ,<br> ecutwfc = 50.0 ,<br> ecutrho = 400.0 , <br> /<br> &electrons <br> conv_thr = 1.0d-10 ,<br> /<br>&ions<br> ion_dynamics = none ,<br>
/<br> &CELL<br> cell_dynamics = none ,<br> press = 0 , <br> press_conv_thr = 0.05D0 ,<br>/<br>ATOMIC_SPECIES<br> O 15.9994 O.pbe-van_ak.UPF<br> Al 26.981538 Al.pbe-n-van.UPF<br> P 30.973761 P.pbe-van_ak.UPF<br>
<br>ATOMIC_POSITIONS crystal<br>Al 0.4692 0.0000 0.1667 1 0 0<br>Al 0.0000 0.4692 0.8333 0 1 0<br>Al 0.5308 0.5308 0.5000 1 1 0<br>P 0.4694 0.0000 0.6667 1 0 0<br>
P 0.0000 0.4694 0.3333 0 1 0<br>P 0.5306 0.5306 0.0000 1 1 0<br>O 0.4160 0.2860 0.1011<br>O 0.7140 0.1300 0.7678<br>O 0.8700 0.5840 0.4344<br>O 0.2860 0.4160 0.8989<br>
O 0.1300 0.7140 0.2322<br>O 0.5840 0.8700 0.5655<br>O 0.4150 0.2538 0.6140<br>O 0.7462 0.1612 0.2807<br>O 0.8388 0.5850 0.9473<br>O 0.2538 0.4150 0.3860<br>
O 0.1612 0.7462 0.7193<br>O 0.5850 0.8388 0.0526<br><br>K_POINTS automatic<br>4 4 3 0 0 0<br><br><br><br><br><br><br><br><br><br><br><br><br><br><br><br><br><br><br><br><br><br><br><br><br><br> Program PWSCF v.4.3.2 starts on 11Aug2011 at 11:58: 4 <br>
This program is part of the open-source Quantum ESPRESSO suite<br> for quantum simulation of materials; please cite<br> "P. Giannozzi et al., J. Phys.:Condens. Matter 21 395502 (2009);<br> URL <a href="http://www.quantum-espresso.org">http://www.quantum-espresso.org</a>", <br>
in publications or presentations arising from this work. More details at<br> <a href="http://www.quantum-espresso.org/wiki/index.php/Citing_Quantum-ESPRESSO">http://www.quantum-espresso.org/wiki/index.php/Citing_Quantum-ESPRESSO</a><br>
Parallel version (MPI), running on 3 processors<br> R & G space division: proc/pool = 3<br> Current dimensions of program PWSCF are:<br> Max number of different atomic species (ntypx) = 10<br>
Max number of k-points (npk) = 40000<br> Max angular momentum in pseudopotentials (lmaxx) = 3<br> Waiting for input...<br> Reading input from stdin<br>Warning: card &IONS ignored<br>Warning: card ION_DYNAMICS = NONE , ignored<br>
Warning: card / ignored<br>Warning: card &CELL ignored<br>Warning: card CELL_DYNAMICS = NONE , ignored<br>Warning: card PRESS = 0 , ignored<br>Warning: card PRESS_CONV_THR = 0.05D0 , ignored<br>Warning: card / ignored<br>
Subspace diagonalization in iterative solution of the eigenvalue problem:<br> a serial algorithm will be used<br> Parallelization info<br> --------------------<br> sticks: dense smooth PW G-vecs: dense smooth PW<br>
Min 826 412 116 73322 25972 3822<br> Max 827 413 117 73331 25986 3827<br> Sum 2479 1237 349 219979 77939 11473<br>
bravais-lattice index = 4<br> lattice parameter (alat) = 9.4718 a.u.<br> unit-cell volume = 1628.2022 (a.u.)^3<br> number of atoms/cell = 18<br> number of atomic types = 3<br>
number of electrons = 96.00<br> number of Kohn-Sham states= 48<br> kinetic-energy cutoff = 50.0000 Ry<br> charge density cutoff = 400.0000 Ry<br> convergence threshold = 1.0E-10<br>
mixing beta = 0.7000<br> number of iterations used = 8 plain mixing<br> Exchange-correlation = SLA PW PBE PBE ( 1 4 3 4 0)<br> EXX-fraction = 0.00<br>
celldm(1)= 9.471757 celldm(2)= 0.000000 celldm(3)= 2.212514<br> celldm(4)= 0.000000 celldm(5)= 0.000000 celldm(6)= 0.000000<br> crystal axes: (cart. coord. in units of alat)<br> a(1) = ( 1.000000 0.000000 0.000000 ) <br>
a(2) = ( -0.500000 0.866025 0.000000 ) <br> a(3) = ( 0.000000 0.000000 2.212514 ) <br> reciprocal axes: (cart. coord. in units 2 pi/alat)<br> b(1) = ( 1.000000 0.577350 0.000000 ) <br>
b(2) = ( 0.000000 1.154701 0.000000 ) <br> b(3) = ( 0.000000 0.000000 0.451975 ) <br> PseudoPot. # 1 for O read from file:<br> ~/pseudo/O.pbe-van_ak.UPF<br> MD5 check sum: <br>
Pseudo is Ultrasoft, Zval = 6.0<br> Generated by new atomic code, or converted to UPF format<br> Using radial grid of 737 points, 4 beta functions with: <br> l(1) = 0<br> l(2) = 0<br>
l(3) = 1<br> l(4) = 1<br> Q(r) pseudized with 8 coefficients, rinner = 0.800 0.800 0.800<br> PseudoPot. # 2 for Al read from file:<br> ~/pseudo/Al.pbe-n-van.UPF<br> MD5 check sum: <br>
Pseudo is Ultrasoft + core correction, Zval = 3.0<br> Generated by new atomic code, or converted to UPF format<br> Using radial grid of 893 points, 2 beta functions with: <br> l(1) = 0<br>
l(2) = 1<br> Q(r) pseudized with 8 coefficients, rinner = 0.900 0.900 0.900<br> PseudoPot. # 3 for P read from file:<br> ~/pseudo/P.pbe-van_ak.UPF<br> MD5 check sum: <br> Pseudo is Ultrasoft, Zval = 5.0<br>
Generated by new atomic code, or converted to UPF format<br> Using radial grid of 807 points, 6 beta functions with: <br> l(1) = 0<br> l(2) = 0<br> l(3) = 1<br>
l(4) = 1<br> l(5) = 2<br> l(6) = 2<br> Q(r) pseudized with 8 coefficients, rinner = 1.100 1.100 1.100<br> 1.100 1.100<br>
atomic species valence mass pseudopotential<br> O 6.00 15.99940 O ( 1.00)<br> Al 3.00 26.98154 Al( 1.00)<br> P 5.00 30.97376 P ( 1.00)<br>
No symmetry found<br> Cartesian axes<br> site n. atom positions (alat units)<br> 1 Al tau( 1) = ( 0.4692000 0.0000000 0.3688261 )<br> 2 Al tau( 2) = ( -0.2346000 0.4063391 1.8436879 )<br>
3 Al tau( 3) = ( 0.2654000 0.4596863 1.1062570 )<br> 4 P tau( 4) = ( 0.4694000 0.0000000 1.4750830 )<br> 5 P tau( 5) = ( -0.2347000 0.4065123 0.7374309 )<br>
6 P tau( 6) = ( 0.2653000 0.4595131 0.0000000 )<br> 7 O tau( 7) = ( 0.2730000 0.2476833 0.2236852 )<br> 8 O tau( 8) = ( 0.6490000 0.1125833 1.6987682 )<br>
9 O tau( 9) = ( 0.5780000 0.5057588 0.9611161 )<br> 10 O tau( 10) = ( 0.0780000 0.3602666 1.9888288 )<br> 11 O tau( 11) = ( -0.2270000 0.6183421 0.5137457 )<br>
12 O tau( 12) = ( 0.1490000 0.7534421 1.2511766 )<br> 13 O tau( 13) = ( 0.2881000 0.2197972 1.3584836 )<br> 14 O tau( 14) = ( 0.6656000 0.1396033 0.6210527 )<br>
15 O tau( 15) = ( 0.5463000 0.5066249 2.0959145 )<br> 16 O tau( 16) = ( 0.0463000 0.3594005 0.8540304 )<br> 17 O tau( 17) = ( -0.2119000 0.6462282 1.5914613 )<br>
18 O tau( 18) = ( 0.1656000 0.7264221 0.1163782 )<br> number of k points= 26<br> cart. coord. in units 2pi/alat<br> k( 1) = ( 0.0000000 0.0000000 0.0000000), wk = 0.0416667<br>
k( 2) = ( 0.0000000 0.0000000 0.1506582), wk = 0.0833333<br> k( 3) = ( 0.0000000 0.2886751 0.0000000), wk = 0.0833333<br> k( 4) = ( 0.0000000 0.2886751 0.1506582), wk = 0.0833333<br>
k( 5) = ( 0.0000000 -0.5773503 0.0000000), wk = 0.0416667<br> k( 6) = ( 0.0000000 -0.5773503 0.1506582), wk = 0.0833333<br> k( 7) = ( 0.2500000 0.4330127 0.0000000), wk = 0.0833333<br>
k( 8) = ( 0.2500000 0.4330127 0.1506582), wk = 0.0833333<br> k( 9) = ( -0.2500000 0.1443376 0.0000000), wk = 0.0833333<br> k( 10) = ( 0.2500000 0.1443376 0.0000000), wk = 0.0833333<br>
k( 11) = ( 0.0000000 -0.2886751 0.1506582), wk = 0.0833333<br> k( 12) = ( -0.2500000 0.1443376 0.1506582), wk = 0.0833333<br> k( 13) = ( 0.2500000 0.1443376 0.1506582), wk = 0.0833333<br>
k( 14) = ( -0.2500000 -0.1443376 0.1506582), wk = 0.0833333<br> k( 15) = ( 0.2500000 -0.1443376 0.1506582), wk = 0.0833333<br> k( 16) = ( 0.5000000 -0.2886751 0.0000000), wk = 0.0416667<br>
k( 17) = ( -0.5000000 -0.2886751 0.0000000), wk = 0.0416667<br> k( 18) = ( 0.5000000 -0.2886751 0.1506582), wk = 0.0833333<br> k( 19) = ( -0.5000000 -0.2886751 0.1506582), wk = 0.0833333<br>
k( 20) = ( -0.2500000 0.4330127 0.0000000), wk = 0.0833333<br> k( 21) = ( 0.5000000 0.0000000 0.0000000), wk = 0.0833333<br> k( 22) = ( -0.2500000 -0.4330127 0.1506582), wk = 0.0833333<br>
k( 23) = ( -0.2500000 0.4330127 -0.1506582), wk = 0.0833333<br> k( 24) = ( 0.2500000 -0.4330127 -0.1506582), wk = 0.0833333<br> k( 25) = ( 0.5000000 0.0000000 0.1506582), wk = 0.0833333<br>
k( 26) = ( -0.5000000 0.0000000 0.1506582), wk = 0.0833333<br> Dense grid: 219979 G-vectors FFT dimensions: ( 64, 64, 135)<br> Smooth grid: 77939 G-vectors FFT dimensions: ( 45, 45, 96)<br>
Largest allocated arrays est. size (Mb) dimensions<br> Kohn-Sham Wavefunctions 2.39 Mb ( 3264, 48)<br> NL pseudopotentials 8.07 Mb ( 3264, 162)<br> Each V/rho on FFT grid 2.81 Mb ( 184320)<br>
Each G-vector array 0.56 Mb ( 73326)<br> G-vector shells 0.07 Mb ( 8914)<br> Largest temporary arrays est. size (Mb) dimensions<br> Auxiliary wavefunctions 9.56 Mb ( 3264, 192)<br>
Each subspace H/S matrix 0.56 Mb ( 192, 192)<br> Each <psi_i|beta_j> matrix 0.12 Mb ( 162, 48)<br> Arrays for rho mixing 22.50 Mb ( 184320, 8)<br> Initial potential from superposition of free atoms<br>
starting charge 95.99220, renormalised to 96.00000<br> Starting wfc are 72 atomic wfcs<br> total cpu time spent up to now is 30.1 secs<br> per-process dynamical memory: 60.5 Mb<br> Self-consistent Calculation<br>
iteration # 1 ecut= 50.00 Ry beta=0.70<br> Davidson diagonalization with overlap<br> ethr = 1.00E-02, avg # of iterations = 2.0<br> total cpu time spent up to now is 98.3 secs<br> total energy = -461.81705211 Ry<br>
Harris-Foulkes estimate = -465.06353120 Ry<br> estimated scf accuracy < 4.61917014 Ry<br> iteration # 2 ecut= 50.00 Ry beta=0.70<br> Davidson diagonalization with overlap<br> ethr = 4.81E-03, avg # of iterations = 3.0<br>
total cpu time spent up to now is 183.2 secs<br> total energy = -463.04537545 Ry<br> Harris-Foulkes estimate = -464.89517769 Ry<br> estimated scf accuracy < 3.74928200 Ry<br>
iteration # 3 ecut= 50.00 Ry beta=0.70<br> Davidson diagonalization with overlap<br> ethr = 3.91E-03, avg # of iterations = 2.2<br> total cpu time spent up to now is 248.6 secs<br> total energy = -463.78132866 Ry<br>
Harris-Foulkes estimate = -463.78151842 Ry<br> estimated scf accuracy < 0.01633610 Ry<br> iteration # 4 ecut= 50.00 Ry beta=0.70<br> Davidson diagonalization with overlap<br> ethr = 1.70E-05, avg # of iterations = 3.0<br>
total cpu time spent up to now is 334.9 secs<br> total energy = -463.79565591 Ry<br> Harris-Foulkes estimate = -463.79691483 Ry<br> estimated scf accuracy < 0.00321029 Ry<br>
iteration # 5 ecut= 50.00 Ry beta=0.70<br> Davidson diagonalization with overlap<br> ethr = 3.34E-06, avg # of iterations = 2.0<br> total cpu time spent up to now is 396.0 secs<br> total energy = -463.79574724 Ry<br>
Harris-Foulkes estimate = -463.79592447 Ry<br> estimated scf accuracy < 0.00052331 Ry<br> iteration # 6 ecut= 50.00 Ry beta=0.70<br> Davidson diagonalization with overlap<br> ethr = 5.45E-07, avg # of iterations = 3.0<br>
total cpu time spent up to now is 471.5 secs<br> total energy = -463.79586898 Ry<br> Harris-Foulkes estimate = -463.79591371 Ry<br> estimated scf accuracy < 0.00010453 Ry<br>
iteration # 7 ecut= 50.00 Ry beta=0.70<br> Davidson diagonalization with overlap<br> ethr = 1.09E-07, avg # of iterations = 3.0<br> total cpu time spent up to now is 536.1 secs<br> total energy = -463.79589022 Ry<br>
Harris-Foulkes estimate = -463.79589247 Ry<br> estimated scf accuracy < 0.00000826 Ry<br> iteration # 8 ecut= 50.00 Ry beta=0.70<br> Davidson diagonalization with overlap<br> ethr = 8.61E-09, avg # of iterations = 2.0<br>
total cpu time spent up to now is 591.2 secs<br> total energy = -463.79589151 Ry<br> Harris-Foulkes estimate = -463.79589156 Ry<br> estimated scf accuracy < 0.00000067 Ry<br>
iteration # 9 ecut= 50.00 Ry beta=0.70<br> Davidson diagonalization with overlap<br> ethr = 6.95E-10, avg # of iterations = 2.0<br> total cpu time spent up to now is 655.0 secs<br> total energy = -463.79589164 Ry<br>
Harris-Foulkes estimate = -463.79589171 Ry<br> estimated scf accuracy < 0.00000014 Ry<br> iteration # 10 ecut= 50.00 Ry beta=0.70<br> Davidson diagonalization with overlap<br> ethr = 1.44E-10, avg # of iterations = 2.0<br>
total cpu time spent up to now is 712.5 secs<br> total energy = -463.79589166 Ry<br> Harris-Foulkes estimate = -463.79589166 Ry<br> estimated scf accuracy < 4.4E-09 Ry<br>
iteration # 11 ecut= 50.00 Ry beta=0.70<br> Davidson diagonalization with overlap<br> ethr = 4.56E-12, avg # of iterations = 3.0<br> total cpu time spent up to now is 791.6 secs<br> total energy = -463.79589166 Ry<br>
Harris-Foulkes estimate = -463.79589166 Ry<br> estimated scf accuracy < 1.3E-09 Ry<br> iteration # 12 ecut= 50.00 Ry beta=0.70<br> Davidson diagonalization with overlap<br> ethr = 1.33E-12, avg # of iterations = 2.0<br>
total cpu time spent up to now is 850.1 secs<br> total energy = -463.79589166 Ry<br> Harris-Foulkes estimate = -463.79589166 Ry<br> estimated scf accuracy < 2.2E-10 Ry<br>
iteration # 13 ecut= 50.00 Ry beta=0.70<br> Davidson diagonalization with overlap<br> ethr = 2.28E-13, avg # of iterations = 2.8<br> total cpu time spent up to now is 913.2 secs<br> End of self-consistent calculation<br>
k = 0.0000 0.0000 0.0000 ( 9719 PWs) bands (ev):<br> -18.2892 -18.0881 -18.0867 -15.9173 -15.9156 -15.6028 -15.5923 -15.5905<br> -15.5665 -15.5576 -15.5529 -15.5503 -7.2679 -6.4153 -6.4127 -4.2540<br>
-4.2515 -3.7637 -3.7619 -3.6202 -3.5868 -3.4079 -3.4061 -3.0504<br> -0.9591 -0.9561 -0.7449 -0.7435 -0.7311 -0.6583 -0.4266 -0.1731<br> -0.1724 0.0091 0.0100 0.3102 0.3274 0.3422 0.3426 0.9407<br>
1.0990 1.1019 1.3139 1.6728 1.6740 1.7044 1.7065 1.8290<br> k = 0.0000 0.0000 0.1507 ( 9791 PWs) bands (ev):<br> -18.2606 -18.1839 -18.0176 -16.0194 -15.7681 -15.6458 -15.5963 -15.5839<br>
-15.5741 -15.5569 -15.5532 -15.5501 -7.1578 -6.8497 -6.0144 -4.5881<br> -3.9280 -3.7872 -3.7167 -3.6867 -3.6010 -3.5203 -3.2467 -3.1067<br> -1.0906 -0.8931 -0.7996 -0.7174 -0.6652 -0.6287 -0.3709 -0.2407<br>
-0.0692 -0.0402 0.1734 0.1855 0.2979 0.3232 0.3432 0.9522<br> 0.9906 1.1059 1.3723 1.5241 1.6702 1.7786 1.7914 1.8207<br> k = 0.0000 0.2887 0.0000 ( 9730 PWs) bands (ev):<br>
-18.2028 -18.0635 -18.0394 -15.9736 -15.8951 -15.7409 -15.6275 -15.6122<br> -15.5962 -15.5852 -15.5731 -15.5553 -6.9070 -6.4340 -6.1449 -4.5657<br> -4.3553 -3.7248 -3.6578 -3.5693 -3.5436 -3.4765 -3.2220 -3.1165<br>
-1.1419 -0.9691 -0.9373 -0.8182 -0.7052 -0.5524 -0.3345 -0.3186<br> -0.1866 -0.1364 0.1072 0.1527 0.2525 0.3475 0.4208 1.0729<br> 1.1139 1.1541 1.2769 1.5174 1.5556 1.7135 1.7762 1.7995<br>
k = 0.0000 0.2887 0.1507 ( 9717 PWs) bands (ev):<br> -18.1814 -18.1226 -18.0012 -16.0039 -15.8329 -15.7774 -15.6227 -15.6133<br> -15.5969 -15.5823 -15.5724 -15.5580 -6.8334 -6.6295 -5.9790 -4.7206<br>
-4.1941 -3.8758 -3.6399 -3.5932 -3.4333 -3.4090 -3.2624 -3.1490<br> -1.1648 -1.0560 -0.8086 -0.7578 -0.6559 -0.6201 -0.4106 -0.3292<br> -0.1756 -0.1214 0.0719 0.1708 0.2954 0.3514 0.3973 1.0566<br>
1.1598 1.1796 1.2866 1.3636 1.5813 1.6770 1.8160 1.8562<br> k = 0.0000-0.5774 0.0000 ( 9752 PWs) bands (ev):<br> -18.0901 -18.0392 -18.0053 -16.0145 -15.9451 -15.8680 -15.6520 -15.6364<br>
-15.6106 -15.6021 -15.5819 -15.5705 -6.4558 -6.3777 -5.9407 -4.8536<br> -4.4030 -3.9749 -3.6302 -3.5706 -3.5100 -3.2781 -3.1998 -3.0009<br> -1.2918 -1.0762 -1.0683 -0.8035 -0.7343 -0.5019 -0.3689 -0.3009<br>
-0.2164 -0.1559 -0.0653 0.0759 0.2251 0.3694 0.4393 1.1205<br> 1.1660 1.2487 1.3104 1.3631 1.4223 1.6839 1.7546 1.8830<br> k = 0.0000-0.5774 0.1507 ( 9714 PWs) bands (ev):<br>
-18.0893 -18.0400 -18.0054 -16.0159 -15.9320 -15.8819 -15.6455 -15.6368<br> -15.6205 -15.5973 -15.5805 -15.5707 -6.4371 -6.3980 -5.9395 -4.8554<br> -4.3601 -4.0707 -3.6306 -3.5432 -3.4594 -3.2811 -3.1677 -3.0497<br>
-1.2868 -1.1071 -1.0153 -0.8777 -0.6611 -0.5268 -0.4198 -0.3497<br> -0.1828 -0.1036 -0.0209 0.1102 0.1965 0.2851 0.4885 1.0658<br> 1.1184 1.3026 1.3539 1.3990 1.4486 1.6281 1.7850 1.8486<br>
k = 0.2500 0.4330 0.0000 ( 9737 PWs) bands (ev):<br> -18.0894 -18.0357 -18.0102 -16.0145 -15.9215 -15.8859 -15.6482 -15.6262<br> -15.6169 -15.6034 -15.5916 -15.5730 -6.4747 -6.3693 -5.8972 -4.7356<br>
-4.5108 -4.1846 -3.5807 -3.4958 -3.4199 -3.3437 -3.2277 -2.9957<br> -1.2180 -1.0435 -0.9521 -0.8658 -0.7706 -0.6387 -0.3718 -0.2706<br> -0.2577 -0.0745 0.0016 0.0988 0.1388 0.3787 0.4323 1.0934<br>
1.2037 1.2586 1.3412 1.4066 1.4469 1.6178 1.6624 1.8936<br> k = 0.2500 0.4330 0.1507 ( 9719 PWs) bands (ev):<br> -18.0808 -18.0535 -18.0009 -16.0154 -15.9125 -15.8940 -15.6410 -15.6292<br>
-15.6179 -15.6087 -15.5968 -15.5656 -6.4513 -6.3968 -5.8929 -4.7509<br> -4.4560 -4.2412 -3.6092 -3.4868 -3.3978 -3.2902 -3.2025 -3.0608<br> -1.1946 -1.0517 -1.0015 -0.8280 -0.7645 -0.6548 -0.3837 -0.3060<br>
-0.1825 -0.0975 -0.0051 0.0633 0.2266 0.3345 0.4507 1.0731<br> 1.1491 1.2658 1.3164 1.4230 1.5119 1.6146 1.7064 1.8484<br> k =-0.2500 0.1443 0.0000 ( 9730 PWs) bands (ev):<br>
-18.2029 -18.0641 -18.0391 -15.9727 -15.8938 -15.7377 -15.6268 -15.6091<br> -15.5967 -15.5860 -15.5783 -15.5578 -6.9074 -6.4361 -6.1428 -4.5649<br> -4.3540 -3.7244 -3.6592 -3.5703 -3.5452 -3.4751 -3.2216 -3.1163<br>
-1.1419 -0.9683 -0.9364 -0.8197 -0.7046 -0.5528 -0.3343 -0.3185<br> -0.1881 -0.1364 0.1085 0.1517 0.2534 0.3486 0.4212 1.0719<br> 1.1155 1.1533 1.2771 1.5161 1.5568 1.7140 1.7744 1.7983<br>
k = 0.2500 0.1443 0.0000 ( 9730 PWs) bands (ev):<br> -18.2027 -18.0629 -18.0402 -15.9735 -15.8950 -15.7406 -15.6273 -15.6118<br> -15.5964 -15.5857 -15.5733 -15.5555 -6.9070 -6.4342 -6.1445 -4.5654<br>
-4.3556 -3.7248 -3.6580 -3.5696 -3.5437 -3.4761 -3.2220 -3.1165<br> -1.1419 -0.9690 -0.9374 -0.8184 -0.7052 -0.5525 -0.3345 -0.3185<br> -0.1865 -0.1364 0.1075 0.1521 0.2526 0.3478 0.4209 1.0725<br>
1.1145 1.1538 1.2771 1.5167 1.5559 1.7134 1.7759 1.7993<br> k = 0.0000-0.2887 0.1507 ( 9717 PWs) bands (ev):<br> -18.1813 -18.1228 -18.0011 -16.0018 -15.8464 -15.7627 -15.6256 -15.6094<br>
-15.6046 -15.5813 -15.5681 -15.5599 -6.8328 -6.6303 -5.9784 -4.7241<br> -4.1774 -3.9030 -3.6458 -3.5767 -3.4484 -3.3871 -3.2483 -3.1661<br> -1.1648 -1.0323 -0.8494 -0.7277 -0.6909 -0.5991 -0.4087 -0.3360<br>
-0.1734 -0.1418 0.0787 0.1775 0.3029 0.3452 0.4092 1.0430<br> 1.1107 1.2436 1.2826 1.3812 1.5887 1.7005 1.7490 1.8735<br> k =-0.2500 0.1443 0.1507 ( 9717 PWs) bands (ev):<br>
-18.1814 -18.1229 -18.0011 -16.0007 -15.8458 -15.7591 -15.6251 -15.6069<br> -15.6036 -15.5845 -15.5711 -15.5628 -6.8334 -6.6308 -5.9777 -4.7226<br> -4.1773 -3.9019 -3.6477 -3.5786 -3.4473 -3.3866 -3.2482 -3.1660<br>
-1.1645 -1.0316 -0.8498 -0.7290 -0.6901 -0.5991 -0.4065 -0.3373<br> -0.1744 -0.1432 0.0799 0.1781 0.3025 0.3467 0.4097 1.0439<br> 1.1099 1.2427 1.2827 1.3824 1.5886 1.6997 1.7475 1.8731<br>
k = 0.2500 0.1443 0.1507 ( 9717 PWs) bands (ev):<br> -18.1811 -18.1227 -18.0013 -16.0016 -15.8467 -15.7621 -15.6255 -15.6092<br> -15.6043 -15.5822 -15.5680 -15.5602 -6.8328 -6.6301 -5.9782 -4.7238<br>
-4.1778 -3.9027 -3.6461 -3.5769 -3.4482 -3.3870 -3.2483 -3.1661<br> -1.1649 -1.0322 -0.8496 -0.7278 -0.6910 -0.5992 -0.4083 -0.3362<br> -0.1734 -0.1421 0.0790 0.1774 0.3027 0.3455 0.4095 1.0434<br>
1.1103 1.2433 1.2827 1.3815 1.5885 1.7004 1.7486 1.8734<br> k =-0.2500-0.1443 0.1507 ( 9717 PWs) bands (ev):<br> -18.1812 -18.1226 -18.0014 -16.0040 -15.8323 -15.7774 -15.6227 -15.6132<br>
-15.5965 -15.5826 -15.5731 -15.5580 -6.8332 -6.6297 -5.9785 -4.7206<br> -4.1940 -3.8761 -3.6398 -3.5936 -3.4333 -3.4088 -3.2624 -3.1489<br> -1.1649 -1.0559 -0.8088 -0.7579 -0.6559 -0.6200 -0.4103 -0.3295<br>
-0.1755 -0.1216 0.0720 0.1707 0.2955 0.3513 0.3978 1.0567<br> 1.1599 1.1794 1.2862 1.3638 1.5813 1.6768 1.8155 1.8562<br> k = 0.2500-0.1443 0.1507 ( 9717 PWs) bands (ev):<br>
-18.1815 -18.1228 -18.0012 -16.0033 -15.8302 -15.7757 -15.6228 -15.6105<br> -15.5961 -15.5827 -15.5765 -15.5618 -6.8336 -6.6307 -5.9778 -4.7200<br> -4.1923 -3.8763 -3.6407 -3.5952 -3.4330 -3.4085 -3.2624 -3.1485<br>
-1.1646 -1.0553 -0.8083 -0.7586 -0.6555 -0.6203 -0.4093 -0.3312<br> -0.1761 -0.1215 0.0720 0.1714 0.2964 0.3510 0.3986 1.0568<br> 1.1600 1.1792 1.2864 1.3641 1.5814 1.6762 1.8149 1.8554<br>
k = 0.5000-0.2887 0.0000 ( 9752 PWs) bands (ev):<br> -18.0905 -18.0398 -18.0051 -16.0101 -15.9437 -15.8668 -15.6498 -15.6335<br> -15.6109 -15.6026 -15.5883 -15.5750 -6.4579 -6.3801 -5.9375 -4.8501<br>
-4.4021 -3.9738 -3.6320 -3.5722 -3.5135 -3.2777 -3.1990 -2.9994<br> -1.2914 -1.0771 -1.0687 -0.8009 -0.7340 -0.5005 -0.3703 -0.3014<br> -0.2178 -0.1562 -0.0646 0.0745 0.2273 0.3707 0.4398 1.1219<br>
1.1656 1.2463 1.3118 1.3633 1.4226 1.6832 1.7520 1.8818<br> k =-0.5000-0.2887 0.0000 ( 9752 PWs) bands (ev):<br> -18.0895 -18.0386 -18.0066 -16.0140 -15.9449 -15.8682 -15.6514 -15.6360<br>
-15.6106 -15.6028 -15.5825 -15.5706 -6.4561 -6.3779 -5.9394 -4.8531<br> -4.4036 -3.9751 -3.6302 -3.5711 -3.5105 -3.2781 -3.1995 -3.0005<br> -1.2920 -1.0763 -1.0683 -0.8032 -0.7343 -0.5016 -0.3691 -0.3012<br>
-0.2168 -0.1559 -0.0653 0.0758 0.2254 0.3696 0.4394 1.1209<br> 1.1660 1.2477 1.3109 1.3629 1.4222 1.6836 1.7542 1.8827<br> k = 0.5000-0.2887 0.1507 ( 9714 PWs) bands (ev):<br>
-18.0897 -18.0406 -18.0051 -16.0117 -15.9301 -15.8809 -15.6441 -15.6333<br> -15.6202 -15.5980 -15.5867 -15.5754 -6.4392 -6.4003 -5.9363 -4.8519<br> -4.3590 -4.0701 -3.6325 -3.5461 -3.4606 -3.2820 -3.1662 -3.0480<br>
-1.2867 -1.1082 -1.0154 -0.8760 -0.6600 -0.5251 -0.4208 -0.3508<br> -0.1826 -0.1050 -0.0212 0.1104 0.1988 0.2846 0.4900 1.0639<br> 1.1199 1.3032 1.3547 1.3987 1.4480 1.6273 1.7829 1.8470<br>
k =-0.5000-0.2887 0.1507 ( 9714 PWs) bands (ev):<br> -18.0887 -18.0394 -18.0066 -16.0155 -15.9317 -15.8821 -15.6451 -15.6363<br> -15.6203 -15.5982 -15.5815 -15.5703 -6.4373 -6.3982 -5.9382 -4.8549<br>
-4.3606 -4.0710 -3.6306 -3.5436 -3.4597 -3.2813 -3.1672 -3.0492<br> -1.2870 -1.1071 -1.0154 -0.8774 -0.6610 -0.5267 -0.4198 -0.3498<br> -0.1829 -0.1039 -0.0211 0.1103 0.1970 0.2847 0.4890 1.0651<br>
1.1190 1.3027 1.3541 1.3987 1.4482 1.6278 1.7847 1.8482<br> k =-0.2500 0.4330 0.0000 ( 9737 PWs) bands (ev):<br> -18.0900 -18.0362 -18.0095 -16.0125 -15.9205 -15.8852 -15.6482 -15.6246<br>
-15.6156 -15.6072 -15.5928 -15.5743 -6.4755 -6.3704 -5.8962 -4.7338<br> -4.5101 -4.1856 -3.5821 -3.4945 -3.4210 -3.3432 -3.2276 -2.9957<br> -1.2180 -1.0429 -0.9516 -0.8653 -0.7708 -0.6390 -0.3723 -0.2703<br>
-0.2588 -0.0751 0.0008 0.1004 0.1395 0.3792 0.4330 1.0935<br> 1.2036 1.2592 1.3414 1.4061 1.4461 1.6169 1.6626 1.8931<br> k = 0.5000 0.0000 0.0000 ( 9737 PWs) bands (ev):<br>
-18.0896 -18.0360 -18.0101 -16.0124 -15.9201 -15.8853 -15.6482 -15.6241<br> -15.6154 -15.6076 -15.5933 -15.5743 -6.4755 -6.3706 -5.8957 -4.7338<br> -4.5102 -4.1857 -3.5824 -3.4942 -3.4213 -3.3429 -3.2276 -2.9957<br>
-1.2180 -1.0430 -0.9517 -0.8653 -0.7708 -0.6391 -0.3726 -0.2699<br> -0.2590 -0.0748 0.0005 0.1006 0.1393 0.3789 0.4336 1.0933<br> 1.2036 1.2597 1.3412 1.4057 1.4459 1.6168 1.6626 1.8930<br>
k =-0.2500-0.4330 0.1507 ( 9719 PWs) bands (ev):<br> -18.0808 -18.0535 -18.0009 -16.0154 -15.9126 -15.8940 -15.6411 -15.6293<br> -15.6167 -15.6097 -15.5968 -15.5655 -6.4514 -6.3967 -5.8929 -4.7510<br>
-4.4558 -4.2414 -3.6088 -3.4870 -3.3981 -3.2904 -3.2022 -3.0608<br> -1.1946 -1.0515 -1.0016 -0.8282 -0.7647 -0.6544 -0.3839 -0.3060<br> -0.1823 -0.0974 -0.0050 0.0630 0.2262 0.3347 0.4508 1.0734<br>
1.1488 1.2658 1.3165 1.4224 1.5123 1.6147 1.7067 1.8482<br> k =-0.2500 0.4330-0.1507 ( 9719 PWs) bands (ev):<br> -18.0815 -18.0536 -18.0005 -16.0135 -15.9116 -15.8933 -15.6412 -15.6280<br>
-15.6156 -15.6118 -15.5990 -15.5669 -6.4523 -6.3977 -5.8919 -4.7494<br> -4.4547 -4.2426 -3.6098 -3.4876 -3.3970 -3.2906 -3.2020 -3.0609<br> -1.1947 -1.0511 -1.0011 -0.8280 -0.7646 -0.6546 -0.3844 -0.3065<br>
-0.1828 -0.0977 -0.0054 0.0642 0.2274 0.3343 0.4519 1.0738<br> 1.1483 1.2655 1.3166 1.4228 1.5119 1.6147 1.7061 1.8475<br> k = 0.2500-0.4330-0.1507 ( 9719 PWs) bands (ev):<br>
-18.0816 -18.0535 -18.0006 -16.0136 -15.9109 -15.8939 -15.6409 -15.6278<br> -15.6186 -15.6081 -15.6004 -15.5667 -6.4520 -6.3980 -5.8919 -4.7491<br> -4.4554 -4.2421 -3.6096 -3.4878 -3.3976 -3.2893 -3.2028 -3.0608<br>
-1.1945 -1.0510 -1.0014 -0.8280 -0.7644 -0.6548 -0.3845 -0.3062<br> -0.1826 -0.0978 -0.0058 0.0645 0.2268 0.3354 0.4511 1.0722<br> 1.1498 1.2657 1.3170 1.4234 1.5109 1.6144 1.7057 1.8480<br>
k = 0.5000 0.0000 0.1507 ( 9719 PWs) bands (ev):<br> -18.0810 -18.0537 -18.0009 -16.0134 -15.9109 -15.8937 -15.6413 -15.6276<br> -15.6157 -15.6115 -15.5997 -15.5670 -6.4520 -6.3981 -5.8913 -4.7493<br>
-4.4550 -4.2426 -3.6096 -3.4880 -3.3973 -3.2902 -3.2020 -3.0609<br> -1.1947 -1.0510 -1.0013 -0.8279 -0.7644 -0.6548 -0.3851 -0.3064<br> -0.1826 -0.0973 -0.0054 0.0642 0.2271 0.3346 0.4520 1.0732<br>
1.1487 1.2656 1.3168 1.4228 1.5113 1.6145 1.7059 1.8475<br> k =-0.5000 0.0000 0.1507 ( 9719 PWs) bands (ev):<br> -18.0810 -18.0537 -18.0009 -16.0134 -15.9111 -15.8935 -15.6406 -15.6274<br>
-15.6186 -15.6089 -15.6007 -15.5667 -6.4519 -6.3982 -5.8913 -4.7491<br> -4.4553 -4.2423 -3.6100 -3.4879 -3.3971 -3.2895 -3.2028 -3.0607<br> -1.1945 -1.0511 -1.0014 -0.8279 -0.7646 -0.6547 -0.3844 -0.3068<br>
-0.1824 -0.0979 -0.0055 0.0646 0.2269 0.3348 0.4518 1.0724<br> 1.1494 1.2657 1.3169 1.4237 1.5106 1.6143 1.7058 1.8477<br>! total energy = -463.79589166 Ry<br> Harris-Foulkes estimate = -463.79589166 Ry<br>
estimated scf accuracy < 3.8E-12 Ry<br> The total energy is the sum of the following terms:<br> one-electron contribution = -239.66672104 Ry<br> hartree contribution = 168.32448929 Ry<br>
xc contribution = -126.75238928 Ry<br> ewald contribution = -265.70127063 Ry<br> convergence has been achieved in 13 iterations<br> Forces acting on atoms (Ry/au):<br> atom 1 type 2 force = 0.00036827 0.00002292 -0.00059820<br>
atom 2 type 2 force = 0.00013705 0.00020556 0.00038929<br> atom 3 type 2 force = 0.00009817 -0.00015317 -0.00016344<br> atom 4 type 3 force = 0.00071556 0.00002817 -0.00137366<br>
atom 5 type 3 force = 0.00022662 0.00033906 0.00068870<br> atom 6 type 3 force = 0.00019071 -0.00028185 -0.00040866<br> atom 7 type 1 force = 0.00033604 0.00002755 0.00006190<br>
atom 8 type 1 force = -0.00030375 -0.00013590 -0.00032493<br> atom 9 type 1 force = -0.00014438 -0.00019788 -0.00001937<br> atom 10 type 1 force = -0.00019983 0.00016599 -0.00010529<br>
atom 11 type 1 force = 0.00018213 -0.00023653 0.00007906<br> atom 12 type 1 force = -0.00080929 0.00045370 0.00115239<br> atom 13 type 1 force = 0.00060665 -0.00040630 0.00039096<br>
atom 14 type 1 force = -0.00011932 0.00013773 -0.00021200<br> atom 15 type 1 force = -0.00058975 -0.00029335 0.00029344<br> atom 16 type 1 force = -0.00064759 0.00026839 -0.00041331<br>
atom 17 type 1 force = 0.00030691 0.00002954 0.00004003<br> atom 18 type 1 force = -0.00035419 0.00002637 0.00052308<br> Total force = 0.002481 Total SCF correction = 0.000004<br>
Writing output data file alpo4-154.0_0.save<br> init_run : 26.25s CPU 29.95s WALL ( 1 calls)<br> electrons : 703.87s CPU 883.12s WALL ( 1 calls)<br> forces : 11.24s CPU 15.92s WALL ( 1 calls)<br>
Called by init_run:<br> wfcinit : 21.86s CPU 24.92s WALL ( 1 calls)<br> potinit : 0.80s CPU 0.86s WALL ( 1 calls)<br> Called by electrons:<br> c_bands : 571.02s CPU 720.39s WALL ( 13 calls)<br>
sum_band : 104.28s CPU 127.95s WALL ( 13 calls)<br> v_of_rho : 3.04s CPU 3.64s WALL ( 14 calls)<br> newd : 26.68s CPU 32.12s WALL ( 14 calls)<br> mix_rho : 0.52s CPU 0.63s WALL ( 13 calls)<br>
Called by c_bands:<br> init_us_2 : 4.41s CPU 4.67s WALL ( 728 calls)<br> cegterg : 554.44s CPU 701.58s WALL ( 338 calls)<br> Called by *egterg:<br> h_psi : 420.10s CPU 522.77s WALL ( 1196 calls)<br>
s_psi : 40.40s CPU 42.33s WALL ( 1196 calls)<br> g_psi : 3.46s CPU 3.58s WALL ( 832 calls)<br> cdiaghg : 12.32s CPU 12.98s WALL ( 1170 calls)<br> Called by h_psi:<br>
add_vuspsi : 40.62s CPU 42.71s WALL ( 1196 calls)<br> General routines<br> calbec : 58.93s CPU 75.22s WALL ( 1560 calls)<br> fft : 3.46s CPU 8.48s WALL ( 231 calls)<br>
ffts : 0.09s CPU 0.12s WALL ( 27 calls)<br> fftw : 315.38s CPU 414.00s WALL ( 117346 calls)<br> interpolate : 0.45s CPU 0.74s WALL ( 27 calls)<br> davcio : 0.05s CPU 2.03s WALL ( 1066 calls)<br>
Parallel routines<br> fft_scatter : 74.75s CPU 165.11s WALL ( 117604 calls)<br> PWSCF : 12m21.62s CPU 15m29.43s WALL<br> This run was terminated on: 12:13:34 11Aug2011 <br>=------------------------------------------------------------------------------=<br>
JOB DONE.<br>=------------------------------------------------------------------------------=<br><br><br>-- <br><div>******************************************************************************<br>WANG Riping<br>Ph.D student, </div>
<div>Institute for Study of the Earth's Interior,Okayama University,<br>827 Yamada, Misasa, Tottori-ken 682-0193, Japan<br>Tel: +81-858-43-3739(Office), +81-858-43-1215(Inst)</div>
<div>E-mail: <a href="mailto:wang.riping.81@gmail.com" target="_blank">wang.riping.81@gmail.com</a><br>******************************************************************************<br><br><br></div><br>