<HTML><HEAD></HEAD>
<BODY dir=ltr>
<DIV dir=ltr>
<DIV style="FONT-FAMILY: 'Century'; COLOR: #000000; FONT-SIZE: 14pt">
<DIV style="FONT-FAMILY: ; COLOR: ">
<DIV><FONT size=4>Dear pwscf developer,</FONT></DIV>
<DIV><FONT size=4></FONT></DIV>
<DIV><FONT size=4>I have found an error of QE6.0. I am trying to optimize the
structure of scorodite using LDA+U. When the k-mesh is not gamma, everything is
ok. When the k-mesh is set to gamma, I got many error in output
file.</FONT></DIV>
<DIV><FONT size=4></FONT> </DIV>
<DIV><FONT size=4>My input is as follow:</FONT></DIV>
<DIV><FONT size=4></FONT></DIV>
<DIV>&CONTROL</DIV>
<DIV>title = 'geometry optimization of scorodite spin-polarized'</DIV>
<DIV>calculation = 'vc-relax' ,</DIV>
<DIV>outdir = './tmp/' ,</DIV>
<DIV>pseudo_dir = '/export/soft/shaofeng/pseudo/ncpp' ,</DIV>
<DIV>prefix = 'vc' ,</DIV>
<DIV>etot_conv_thr = 1.0D-4 ,</DIV>
<DIV>forc_conv_thr = 1.0d-3 ,</DIV>
<DIV>tstress = .true. ,</DIV>
<DIV>tprnfor = .true. ,</DIV>
<DIV>nstep = 150 ,</DIV>
<DIV>/</DIV>
<DIV>&SYSTEM</DIV>
<DIV>ibrav = 0,</DIV>
<DIV>celldm(1) = 1.889726,</DIV>
<DIV>nat = 96,</DIV>
<DIV>ntyp = 5,</DIV>
<DIV>ecutwfc = 80 ,</DIV>
<DIV>vdw_corr = 'DFT-D' ,</DIV>
<DIV>occupations = 'smearing',</DIV>
<DIV>smearing = 'mv',</DIV>
<DIV>degauss = 0.02,</DIV>
<DIV>nspin = 2 , </DIV>
<DIV>starting_magnetization(1) = 0.5 , </DIV>
<DIV>starting_magnetization(2) = -0.5 , </DIV>
<DIV>lda_plus_u=.true. Hubbard_U(1)=4, Hubbard_U(2)=4, </DIV>
<DIV>/</DIV>
<DIV>&ELECTRONS</DIV>
<DIV>electron_maxstep = 150,</DIV>
<DIV>conv_thr = 1.0d-8 ,</DIV>
<DIV>! diagonalization = 'cg' ,</DIV>
<DIV>mixing_beta = 0.3 ,</DIV>
<DIV>/</DIV>
<DIV>&IONS</DIV>
<DIV>ion_dynamics = 'bfgs' ,</DIV>
<DIV>/</DIV>
<DIV>&CELL</DIV>
<DIV>cell_dynamics = 'bfgs' ,</DIV>
<DIV>cell_dofree = xyz ,</DIV>
<DIV>/</DIV>
<DIV>ATOMIC_SPECIES</DIV>
<DIV>Fe 55.845 Fe_00LDA_OP.ncpp </DIV>
<DIV>Fe1 55.845 Fe_00LDA_OP.ncpp </DIV>
<DIV>O 15.99940 O_00LDA_OP.ncpp </DIV>
<DIV>H 1.00794 H_00LDA_OP.ncpp </DIV>
<DIV>As 74.92160 As.pz-hgh.UPF </DIV>
<DIV>CELL_PARAMETERS (alat= 1.88972613)</DIV>
<DIV>8.936999999999999 0.000000000000000 0.000000000000000</DIV>
<DIV>0.000000000000000 10.277999999999999 0.000000000000000</DIV>
<DIV>0.000000000000000 0.000000000000000 9.996000000000000</DIV>
<DIV></DIV>
<DIV>ATOMIC_POSITIONS (crystal)</DIV>
<DIV>H 0.3770000000000002 0.4760000000000001 0.3970000000000001</DIV>
<DIV>H 0.3010000000000001 0.5340000000000000 0.5010000000000001</DIV>
<DIV>H 0.1570000000000000 0.8540000000000001 0.2970000000000002</DIV>
<DIV>H 0.3100000000000002 0.8100000000000001 0.4180000000000002</DIV>
<DIV>H -0.3770000000000002 0.0239999999999999 0.8970000000000002</DIV>
<DIV>H -0.3010000000000001 -0.0340000000000000 1.0010000000000001</DIV>
<DIV>H -0.1570000000000000 -0.3540000000000001 0.7970000000000002</DIV>
<DIV>H -0.3100000000000002 -0.3100000000000001 0.9180000000000003</DIV>
<DIV>H 0.8770000000000002 -0.4760000000000001 0.1030000000000000</DIV>
<DIV>H 0.8010000000000002 -0.5340000000000000 -0.0010000000000001</DIV>
<DIV>H 0.6570000000000000 -0.8540000000000001 0.2029999999999999</DIV>
<DIV>H 0.8100000000000002 -0.8100000000000001 0.0819999999999999</DIV>
<DIV>H 0.1229999999999999 0.9760000000000000 -0.3970000000000001</DIV>
<DIV>H 0.1990000000000000 1.0340000000000000 -0.5010000000000001</DIV>
<DIV>H 0.3430000000000000 1.3540000000000001 -0.2970000000000002</DIV>
<DIV>H 0.1899999999999999 1.3100000000000001 -0.4180000000000002</DIV>
<DIV>H -0.3770000000000002 -0.4760000000000001 -0.3970000000000001</DIV>
<DIV>H -0.3010000000000001 -0.5340000000000000 -0.5010000000000001</DIV>
<DIV>H -0.1570000000000000 -0.8540000000000001 -0.2970000000000002</DIV>
<DIV>H -0.3100000000000002 -0.8100000000000001 -0.4180000000000002</DIV>
<DIV>H 0.3770000000000002 0.9760000000000000 0.1030000000000000</DIV>
<DIV>H 0.3010000000000001 1.0340000000000000 -0.0010000000000001</DIV>
<DIV>H 0.1570000000000000 1.3540000000000001 0.2029999999999999</DIV>
<DIV>H 0.3100000000000002 1.3100000000000001 0.0819999999999999</DIV>
<DIV>H 0.1229999999999999 0.4760000000000001 0.8970000000000002</DIV>
<DIV>H 0.1990000000000000 0.5340000000000000 1.0010000000000001</DIV>
<DIV>H 0.3430000000000000 0.8540000000000001 0.7970000000000002</DIV>
<DIV>H 0.1899999999999999 0.8100000000000001 0.9180000000000003</DIV>
<DIV>H 0.8770000000000002 0.0239999999999999 0.3970000000000001</DIV>
<DIV>H 0.8010000000000002 -0.0340000000000000 0.5010000000000001</DIV>
<DIV>H 0.6570000000000000 -0.3540000000000001 0.2970000000000002</DIV>
<DIV>H 0.8100000000000002 -0.3100000000000001 0.4180000000000002</DIV>
<DIV>O 0.4955000000000001 0.3012000000000000 0.3044000000000001</DIV>
<DIV>O 0.5831000000000002 0.0508000000000000 0.2830000000000001</DIV>
<DIV>O 0.2859000000000002 0.1102000000000000 0.3329000000000002</DIV>
<DIV>O 0.4951000000000000 0.1426000000000000 0.5293000000000002</DIV>
<DIV>O 0.3282000000000002 0.5543000000000000 0.4305000000000002</DIV>
<DIV>O 0.2277000000000000 0.8006000000000001 0.3836000000000001</DIV>
<DIV>O -0.4955000000000001 0.1988000000000000 0.8044000000000002</DIV>
<DIV>O -0.5831000000000002 0.4492000000000000 0.7830000000000001</DIV>
<DIV>O -0.2859000000000002 0.3898000000000000 0.8329000000000002</DIV>
<DIV>O -0.4951000000000000 0.3574000000000000 1.0293000000000003</DIV>
<DIV>O -0.3282000000000002 -0.0543000000000000 0.9305000000000002</DIV>
<DIV>O -0.2277000000000000 -0.3006000000000000 0.8836000000000002</DIV>
<DIV>O 0.9955000000000002 -0.3012000000000000 0.1956000000000000</DIV>
<DIV>O 1.0831000000000002 -0.0508000000000000 0.2169999999999999</DIV>
<DIV>O 0.7859000000000002 -0.1102000000000000 0.1670999999999999</DIV>
<DIV>O 0.9951000000000000 -0.1426000000000000 -0.0293000000000001</DIV>
<DIV>O 0.8282000000000002 -0.5543000000000000 0.0694999999999999</DIV>
<DIV>O 0.7277000000000000 -0.8006000000000001 0.1163999999999999</DIV>
<DIV>O 0.0045000000000000 0.8011999999999999 -0.3044000000000001</DIV>
<DIV>O -0.0831000000000001 0.5508000000000000 -0.2830000000000001</DIV>
<DIV>O 0.2140999999999999 0.6102000000000001 -0.3329000000000002</DIV>
<DIV>O 0.0049000000000001 0.6426000000000000 -0.5293000000000002</DIV>
<DIV>O 0.1717999999999999 1.0543000000000000 -0.4305000000000002</DIV>
<DIV>O 0.2723000000000000 1.3006000000000000 -0.3836000000000001</DIV>
<DIV>O -0.4955000000000001 -0.3012000000000000 -0.3044000000000001</DIV>
<DIV>O -0.5831000000000002 -0.0508000000000000 -0.2830000000000001</DIV>
<DIV>O -0.2859000000000002 -0.1102000000000000 -0.3329000000000002</DIV>
<DIV>O -0.4951000000000000 -0.1426000000000000 -0.5293000000000002</DIV>
<DIV>O -0.3282000000000002 -0.5543000000000000 -0.4305000000000002</DIV>
<DIV>O -0.2277000000000000 -0.8006000000000001 -0.3836000000000001</DIV>
<DIV>O 0.4955000000000001 0.8011999999999999 0.1956000000000000</DIV>
<DIV>O 0.5831000000000002 0.5508000000000000 0.2169999999999999</DIV>
<DIV>O 0.2859000000000002 0.6102000000000001 0.1670999999999999</DIV>
<DIV>O 0.4951000000000000 0.6426000000000000 -0.0293000000000001</DIV>
<DIV>O 0.3282000000000002 1.0543000000000000 0.0694999999999999</DIV>
<DIV>O 0.2277000000000000 1.3006000000000000 0.1163999999999999</DIV>
<DIV>O 0.0045000000000000 0.3012000000000000 0.8044000000000002</DIV>
<DIV>O -0.0831000000000001 0.0508000000000000 0.7830000000000001</DIV>
<DIV>O 0.2140999999999999 0.1102000000000000 0.8329000000000002</DIV>
<DIV>O 0.0049000000000001 0.1426000000000000 1.0293000000000003</DIV>
<DIV>O 0.1717999999999999 0.5543000000000000 0.9305000000000002</DIV>
<DIV>O 0.2723000000000000 0.8006000000000001 0.8836000000000002</DIV>
<DIV>O 0.9955000000000002 0.1988000000000000 0.3044000000000001</DIV>
<DIV>O 1.0831000000000002 0.4492000000000000 0.2830000000000001</DIV>
<DIV>O 0.7859000000000002 0.3898000000000000 0.3329000000000002</DIV>
<DIV>O 0.9951000000000000 0.3574000000000000 0.5293000000000002</DIV>
<DIV>O 0.8282000000000002 -0.0543000000000000 0.4305000000000002</DIV>
<DIV>O 0.7277000000000000 -0.3006000000000000 0.3836000000000001</DIV>
<DIV>Fe 0.1468000000000000 0.1269000000000000 0.1823000000000001</DIV>
<DIV>Fe1 -0.1468000000000000 0.3731000000000000 0.6823000000000001</DIV>
<DIV>Fe 0.6468000000000000 -0.1269000000000000 0.3177000000000000</DIV>
<DIV>Fe1 0.3532000000000000 0.6269000000000000 -0.1823000000000001</DIV>
<DIV>Fe -0.1468000000000000 -0.1269000000000000 -0.1823000000000001</DIV>
<DIV>Fe1 0.1468000000000000 0.6269000000000000 0.3177000000000000</DIV>
<DIV>Fe 0.3532000000000000 0.1269000000000000 0.6823000000000001</DIV>
<DIV>Fe1 0.6468000000000000 0.3731000000000000 0.1823000000000001</DIV>
<DIV>As 0.4646000000000001 0.1516000000000000 0.3636000000000002</DIV>
<DIV>As -0.4646000000000001 0.3484000000000000 0.8636000000000002</DIV>
<DIV>As 0.9646000000000002 -0.1516000000000000 0.1363999999999999</DIV>
<DIV>As 0.0353999999999999 0.6516000000000001 -0.3636000000000002</DIV>
<DIV>As -0.4646000000000001 -0.1516000000000000 -0.3636000000000002</DIV>
<DIV>As 0.4646000000000001 0.6516000000000001 0.1363999999999999</DIV>
<DIV>As 0.0353999999999999 0.1516000000000000 0.8636000000000002</DIV>
<DIV>As 0.9646000000000002 0.3484000000000000 0.3636000000000002</DIV>
<DIV>K_POINTS gamma</DIV>
<DIV> </DIV>
<DIV><FONT size=4></FONT> </DIV>
<DIV><FONT size=4>The output file is </FONT></DIV>
<DIV><FONT size=4></FONT> </DIV>
<DIV><FONT size=4></FONT></DIV>
<DIV><FONT size=4></FONT></DIV>
<DIV>Program PWSCF v.6.0 (svn rev. 13079) starts on 17Mar2017 at 8:21: 8 </DIV>
<DIV></DIV>
<DIV>This program is part of the open-source Quantum ESPRESSO suite</DIV>
<DIV>for quantum simulation of materials; please cite</DIV>
<DIV>"P. Giannozzi et al., J. Phys.:Condens. Matter 21 395502 (2009);</DIV>
<DIV>URL <A
href='http://www.quantum-espresso.org"/'>http://www.quantum-espresso.org"/</A>,
</DIV>
<DIV>in publications or presentations arising from this work. More details
at</DIV>
<DIV><A
href="http://www.quantum-espresso.org/quote">http://www.quantum-espresso.org/quote</A></DIV>
<DIV></DIV>
<DIV>Parallel version (MPI), running on 10 processors</DIV>
<DIV>R & G space division: proc/nbgrp/npool/nimage = 10</DIV>
<DIV>Waiting for input...</DIV>
<DIV>Reading input from standard input</DIV>
<DIV></DIV>
<DIV>Current dimensions of program PWSCF are:</DIV>
<DIV>Max number of different atomic species (ntypx) = 10</DIV>
<DIV>Max number of k-points (npk) = 40000</DIV>
<DIV>Max angular momentum in pseudopotentials (lmaxx) = 3</DIV>
<DIV></DIV>
<DIV>-------------------------------------</DIV>
<DIV>Parameters for Dispersion Correction:</DIV>
<DIV>-------------------------------------</DIV>
<DIV>atom VdW radius C_6 </DIV>
<DIV></DIV>
<DIV>Fe 2.952 374.666</DIV>
<DIV>Fe1 2.952 374.666</DIV>
<DIV>O 2.536 24.284</DIV>
<DIV>H 1.892 4.857</DIV>
<DIV>As 3.326 567.896</DIV>
<DIV></DIV>
<DIV>gamma-point specific algorithms are used</DIV>
<DIV></DIV>
<DIV>Subspace diagonalization in iterative solution of the eigenvalue
problem:</DIV>
<DIV>one sub-group per band group will be used</DIV>
<DIV>scalapack distributed-memory algorithm (size of sub-group: 2* 2
procs)</DIV>
<DIV></DIV>
<DIV></DIV>
<DIV>Parallelization info</DIV>
<DIV>--------------------</DIV>
<DIV>sticks: dense smooth PW G-vecs: dense smooth PW</DIV>
<DIV>Min 833 833 208 59883 59883 7480</DIV>
<DIV>Max 836 836 210 59892 59892 7486</DIV>
<DIV>Sum 8355 8355 2091 598889 598889 74839</DIV>
<DIV></DIV>
<DIV></DIV>
<DIV>Title: </DIV>
<DIV>geometry optimization of scorodite spin-polarized </DIV>
<DIV></DIV>
<DIV></DIV>
<DIV>bravais-lattice index = 0</DIV>
<DIV>lattice parameter (alat) = 1.8897 a.u.</DIV>
<DIV>unit-cell volume = 6196.1673 (a.u.)^3</DIV>
<DIV>number of atoms/cell = 96</DIV>
<DIV>number of atomic types = 5</DIV>
<DIV>number of electrons = 488.00</DIV>
<DIV>number of Kohn-Sham states= 293</DIV>
<DIV>kinetic-energy cutoff = 80.0000 Ry</DIV>
<DIV>charge density cutoff = 320.0000 Ry</DIV>
<DIV>convergence threshold = 1.0E-08</DIV>
<DIV>mixing beta = 0.3000</DIV>
<DIV>number of iterations used = 8 plain mixing</DIV>
<DIV>Exchange-correlation = SLA-PZ-NOGX-NOGC ( 1 1 0 0 0 0)</DIV>
<DIV>nstep = 150</DIV>
<DIV></DIV>
<DIV></DIV>
<DIV>celldm(1)= 1.889726 celldm(2)= 0.000000 celldm(3)= 0.000000</DIV>
<DIV>celldm(4)= 0.000000 celldm(5)= 0.000000 celldm(6)= 0.000000</DIV>
<DIV></DIV>
<DIV>crystal axes: (cart. coord. in units of alat)</DIV>
<DIV>a(1) = ( 8.937000 0.000000 0.000000 ) </DIV>
<DIV>a(2) = ( 0.000000 10.278000 0.000000 ) </DIV>
<DIV>a(3) = ( 0.000000 0.000000 9.996000 ) </DIV>
<DIV></DIV>
<DIV>reciprocal axes: (cart. coord. in units 2 pi/alat)</DIV>
<DIV>b(1) = ( 0.111894 0.000000 0.000000 ) </DIV>
<DIV>b(2) = ( 0.000000 0.097295 0.000000 ) </DIV>
<DIV>b(3) = ( 0.000000 0.000000 0.100040 ) </DIV>
<DIV></DIV>
<DIV></DIV>
<DIV>PseudoPot. # 1 for Fe read from file:</DIV>
<DIV>/export/soft/shaofeng/pseudo/ncpp/Fe_00LDA_OP.ncpp</DIV>
<DIV>MD5 check sum: d105913d3df1c460515ff904f80bf6fe</DIV>
<DIV>Pseudo is Norm-conserving, Zval = 16.0</DIV>
<DIV>Generated by old ld1 code (numerical format)</DIV>
<DIV>Using radial grid of 1160 points, 2 beta functions with: </DIV>
<DIV>l(1) = 1</DIV>
<DIV>l(2) = 2</DIV>
<DIV></DIV>
<DIV>PseudoPot. # 2 for Fe read from file:</DIV>
<DIV>/export/soft/shaofeng/pseudo/ncpp/Fe_00LDA_OP.ncpp</DIV>
<DIV>MD5 check sum: d105913d3df1c460515ff904f80bf6fe</DIV>
<DIV>Pseudo is Norm-conserving, Zval = 16.0</DIV>
<DIV>Generated by old ld1 code (numerical format)</DIV>
<DIV>Using radial grid of 1160 points, 2 beta functions with: </DIV>
<DIV>l(1) = 1</DIV>
<DIV>l(2) = 2</DIV>
<DIV></DIV>
<DIV>PseudoPot. # 3 for O read from file:</DIV>
<DIV>/export/soft/shaofeng/pseudo/ncpp/O_00LDA_OP.ncpp</DIV>
<DIV>MD5 check sum: 22ef4299c8c00f6acd2c8f40dd97b777</DIV>
<DIV>Pseudo is Norm-conserving, Zval = 6.0</DIV>
<DIV>Generated by old ld1 code (numerical format)</DIV>
<DIV>Using radial grid of 1130 points, 1 beta functions with: </DIV>
<DIV>l(1) = 1</DIV>
<DIV></DIV>
<DIV>PseudoPot. # 4 for H read from file:</DIV>
<DIV>/export/soft/shaofeng/pseudo/ncpp/H_00LDA_OP.ncpp</DIV>
<DIV>MD5 check sum: 544dd84d65d09735d9c3e39f71958d78</DIV>
<DIV>Pseudo is Norm-conserving, Zval = 1.0</DIV>
<DIV>Generated by old ld1 code (numerical format)</DIV>
<DIV>Using radial grid of 1076 points, 0 beta functions with: </DIV>
<DIV></DIV>
<DIV>PseudoPot. # 5 for As read from file:</DIV>
<DIV>/export/soft/shaofeng/pseudo/ncpp/As.pz-hgh.UPF</DIV>
<DIV>MD5 check sum: 19ccbb18ac7fdc8ac5652dbe56319a83</DIV>
<DIV>Pseudo is Norm-conserving, Zval = 5.0</DIV>
<DIV>Generated in analytical, separable form</DIV>
<DIV>Using radial grid of 1209 points, 6 beta functions with: </DIV>
<DIV>l(1) = 0</DIV>
<DIV>l(2) = 0</DIV>
<DIV>l(3) = 0</DIV>
<DIV>l(4) = 1</DIV>
<DIV>l(5) = 1</DIV>
<DIV>l(6) = 2</DIV>
<DIV></DIV>
<DIV>atomic species valence mass pseudopotential</DIV>
<DIV>Fe 16.00 55.84500 Fe( 1.00)</DIV>
<DIV>Fe1 16.00 55.84500 Fe( 1.00)</DIV>
<DIV>O 6.00 15.99940 O ( 1.00)</DIV>
<DIV>H 1.00 1.00794 H ( 1.00)</DIV>
<DIV>As 5.00 74.92160 As( 1.00)</DIV>
<DIV></DIV>
<DIV>Starting magnetic structure </DIV>
<DIV>atomic species magnetization</DIV>
<DIV>Fe 0.500</DIV>
<DIV>Fe1 -0.500</DIV>
<DIV>O 0.000</DIV>
<DIV>H 0.000</DIV>
<DIV>As 0.000</DIV>
<DIV></DIV>
<DIV></DIV>
<DIV>Simplified LDA+U calculation (l_max = 2) with parameters (eV):</DIV>
<DIV>atomic species L U alpha J0 beta</DIV>
<DIV>Fe 2 4.0000 0.0000 0.0000 0.0000</DIV>
<DIV>Fe1 2 4.0000 0.0000 0.0000 0.0000</DIV>
<DIV></DIV>
<DIV></DIV>
<DIV></DIV>
<DIV>4 Sym. Ops., with inversion, found ( 2 have fractional translation)</DIV>
<DIV></DIV>
<DIV></DIV>
<DIV></DIV>
<DIV>Cartesian axes</DIV>
<DIV></DIV>
<DIV>site n. atom positions (alat units)</DIV>
<DIV>1 H tau( 1) = ( 3.3692490 4.8923280 3.9684120 )</DIV>
<DIV>2 H tau( 2) = ( 2.6900370 5.4884520 5.0079960 )</DIV>
<DIV>3 H tau( 3) = ( 1.4031090 8.7774120 2.9688120 )</DIV>
<DIV>4 H tau( 4) = ( 2.7704700 8.3251800 4.1783280 )</DIV>
<DIV>5 H tau( 5) = ( -3.3692490 0.2466720 8.9664120 )</DIV>
<DIV>6 H tau( 6) = ( -2.6900370 -0.3494520 10.0059960 )</DIV>
<DIV>7 H tau( 7) = ( -1.4031090 -3.6384120 7.9668120 )</DIV>
<DIV>8 H tau( 8) = ( -2.7704700 -3.1861800 9.1763280 )</DIV>
<DIV>9 H tau( 9) = ( 7.8377490 -4.8923280 1.0295880 )</DIV>
<DIV>10 H tau( 10) = ( 7.1585370 -5.4884520 -0.0099960 )</DIV>
<DIV>11 H tau( 11) = ( 5.8716090 -8.7774120 2.0291880 )</DIV>
<DIV>12 H tau( 12) = ( 7.2389700 -8.3251800 0.8196720 )</DIV>
<DIV>13 H tau( 13) = ( 1.0992510 10.0313280 -3.9684120 )</DIV>
<DIV>14 H tau( 14) = ( 1.7784630 10.6274520 -5.0079960 )</DIV>
<DIV>15 H tau( 15) = ( 3.0653910 13.9164120 -2.9688120 )</DIV>
<DIV>16 H tau( 16) = ( 1.6980300 13.4641800 -4.1783280 )</DIV>
<DIV>17 H tau( 17) = ( -3.3692490 -4.8923280 -3.9684120 )</DIV>
<DIV>18 H tau( 18) = ( -2.6900370 -5.4884520 -5.0079960 )</DIV>
<DIV>19 H tau( 19) = ( -1.4031090 -8.7774120 -2.9688120 )</DIV>
<DIV>20 H tau( 20) = ( -2.7704700 -8.3251800 -4.1783280 )</DIV>
<DIV>21 H tau( 21) = ( 3.3692490 10.0313280 1.0295880 )</DIV>
<DIV>22 H tau( 22) = ( 2.6900370 10.6274520 -0.0099960 )</DIV>
<DIV>23 H tau( 23) = ( 1.4031090 13.9164120 2.0291880 )</DIV>
<DIV>24 H tau( 24) = ( 2.7704700 13.4641800 0.8196720 )</DIV>
<DIV>25 H tau( 25) = ( 1.0992510 4.8923280 8.9664120 )</DIV>
<DIV>26 H tau( 26) = ( 1.7784630 5.4884520 10.0059960 )</DIV>
<DIV>27 H tau( 27) = ( 3.0653910 8.7774120 7.9668120 )</DIV>
<DIV>28 H tau( 28) = ( 1.6980300 8.3251800 9.1763280 )</DIV>
<DIV>29 H tau( 29) = ( 7.8377490 0.2466720 3.9684120 )</DIV>
<DIV>30 H tau( 30) = ( 7.1585370 -0.3494520 5.0079960 )</DIV>
<DIV>31 H tau( 31) = ( 5.8716090 -3.6384120 2.9688120 )</DIV>
<DIV>32 H tau( 32) = ( 7.2389700 -3.1861800 4.1783280 )</DIV>
<DIV>33 O tau( 33) = ( 4.4282835 3.0957336 3.0427824 )</DIV>
<DIV>34 O tau( 34) = ( 5.2111647 0.5221224 2.8288680 )</DIV>
<DIV>35 O tau( 35) = ( 2.5550883 1.1326356 3.3276684 )</DIV>
<DIV>36 O tau( 36) = ( 4.4247087 1.4656428 5.2908828 )</DIV>
<DIV>37 O tau( 37) = ( 2.9331234 5.6970954 4.3032780 )</DIV>
<DIV>38 O tau( 38) = ( 2.0349549 8.2285668 3.8344656 )</DIV>
<DIV>39 O tau( 39) = ( -4.4282835 2.0432664 8.0407824 )</DIV>
<DIV>40 O tau( 40) = ( -5.2111647 4.6168776 7.8268680 )</DIV>
<DIV>41 O tau( 41) = ( -2.5550883 4.0063644 8.3256684 )</DIV>
<DIV>42 O tau( 42) = ( -4.4247087 3.6733572 10.2888828 )</DIV>
<DIV>43 O tau( 43) = ( -2.9331234 -0.5580954 9.3012780 )</DIV>
<DIV>44 O tau( 44) = ( -2.0349549 -3.0895668 8.8324656 )</DIV>
<DIV>45 O tau( 45) = ( 8.8967835 -3.0957336 1.9552176 )</DIV>
<DIV>46 O tau( 46) = ( 9.6796647 -0.5221224 2.1691320 )</DIV>
<DIV>47 O tau( 47) = ( 7.0235883 -1.1326356 1.6703316 )</DIV>
<DIV>48 O tau( 48) = ( 8.8932087 -1.4656428 -0.2928828 )</DIV>
<DIV>49 O tau( 49) = ( 7.4016234 -5.6970954 0.6947220 )</DIV>
<DIV>50 O tau( 50) = ( 6.5034549 -8.2285668 1.1635344 )</DIV>
<DIV>51 O tau( 51) = ( 0.0402165 8.2347336 -3.0427824 )</DIV>
<DIV>52 O tau( 52) = ( -0.7426647 5.6611224 -2.8288680 )</DIV>
<DIV>53 O tau( 53) = ( 1.9134117 6.2716356 -3.3276684 )</DIV>
<DIV>54 O tau( 54) = ( 0.0437913 6.6046428 -5.2908828 )</DIV>
<DIV>55 O tau( 55) = ( 1.5353766 10.8360954 -4.3032780 )</DIV>
<DIV>56 O tau( 56) = ( 2.4335451 13.3675668 -3.8344656 )</DIV>
<DIV>57 O tau( 57) = ( -4.4282835 -3.0957336 -3.0427824 )</DIV>
<DIV>58 O tau( 58) = ( -5.2111647 -0.5221224 -2.8288680 )</DIV>
<DIV>59 O tau( 59) = ( -2.5550883 -1.1326356 -3.3276684 )</DIV>
<DIV>60 O tau( 60) = ( -4.4247087 -1.4656428 -5.2908828 )</DIV>
<DIV>61 O tau( 61) = ( -2.9331234 -5.6970954 -4.3032780 )</DIV>
<DIV>62 O tau( 62) = ( -2.0349549 -8.2285668 -3.8344656 )</DIV>
<DIV>63 O tau( 63) = ( 4.4282835 8.2347336 1.9552176 )</DIV>
<DIV>64 O tau( 64) = ( 5.2111647 5.6611224 2.1691320 )</DIV>
<DIV>65 O tau( 65) = ( 2.5550883 6.2716356 1.6703316 )</DIV>
<DIV>66 O tau( 66) = ( 4.4247087 6.6046428 -0.2928828 )</DIV>
<DIV>67 O tau( 67) = ( 2.9331234 10.8360954 0.6947220 )</DIV>
<DIV>68 O tau( 68) = ( 2.0349549 13.3675668 1.1635344 )</DIV>
<DIV>69 O tau( 69) = ( 0.0402165 3.0957336 8.0407824 )</DIV>
<DIV>70 O tau( 70) = ( -0.7426647 0.5221224 7.8268680 )</DIV>
<DIV>71 O tau( 71) = ( 1.9134117 1.1326356 8.3256684 )</DIV>
<DIV>72 O tau( 72) = ( 0.0437913 1.4656428 10.2888828 )</DIV>
<DIV>73 O tau( 73) = ( 1.5353766 5.6970954 9.3012780 )</DIV>
<DIV>74 O tau( 74) = ( 2.4335451 8.2285668 8.8324656 )</DIV>
<DIV>75 O tau( 75) = ( 8.8967835 2.0432664 3.0427824 )</DIV>
<DIV>76 O tau( 76) = ( 9.6796647 4.6168776 2.8288680 )</DIV>
<DIV>77 O tau( 77) = ( 7.0235883 4.0063644 3.3276684 )</DIV>
<DIV>78 O tau( 78) = ( 8.8932087 3.6733572 5.2908828 )</DIV>
<DIV>79 O tau( 79) = ( 7.4016234 -0.5580954 4.3032780 )</DIV>
<DIV>80 O tau( 80) = ( 6.5034549 -3.0895668 3.8344656 )</DIV>
<DIV>81 Fe tau( 81) = ( 1.3119516 1.3042782 1.8222708 )</DIV>
<DIV>82 Fe1 tau( 82) = ( -1.3119516 3.8347218 6.8202708 )</DIV>
<DIV>83 Fe tau( 83) = ( 5.7804516 -1.3042782 3.1757292 )</DIV>
<DIV>84 Fe1 tau( 84) = ( 3.1565484 6.4432782 -1.8222708 )</DIV>
<DIV>85 Fe tau( 85) = ( -1.3119516 -1.3042782 -1.8222708 )</DIV>
<DIV>86 Fe1 tau( 86) = ( 1.3119516 6.4432782 3.1757292 )</DIV>
<DIV>87 Fe tau( 87) = ( 3.1565484 1.3042782 6.8202708 )</DIV>
<DIV>88 Fe1 tau( 88) = ( 5.7804516 3.8347218 1.8222708 )</DIV>
<DIV>89 As tau( 89) = ( 4.1521302 1.5581448 3.6345456 )</DIV>
<DIV>90 As tau( 90) = ( -4.1521302 3.5808552 8.6325456 )</DIV>
<DIV>91 As tau( 91) = ( 8.6206302 -1.5581448 1.3634544 )</DIV>
<DIV>92 As tau( 92) = ( 0.3163698 6.6971448 -3.6345456 )</DIV>
<DIV>93 As tau( 93) = ( -4.1521302 -1.5581448 -3.6345456 )</DIV>
<DIV>94 As tau( 94) = ( 4.1521302 6.6971448 1.3634544 )</DIV>
<DIV>95 As tau( 95) = ( 0.3163698 1.5581448 8.6325456 )</DIV>
<DIV>96 As tau( 96) = ( 8.6206302 3.5808552 3.6345456 )</DIV>
<DIV></DIV>
<DIV>number of k points= 1 Marzari-Vanderbilt smearing, width (Ry)= 0.0200</DIV>
<DIV>cart. coord. in units 2pi/alat</DIV>
<DIV>k( 1) = ( 0.0000000 0.0000000 0.0000000), wk = 1.0000000</DIV>
<DIV></DIV>
<DIV>Dense grid: 299445 G-vectors FFT dimensions: ( 100, 120, 108)</DIV>
<DIV></DIV>
<DIV>Estimated max dynamical RAM per process > 180.71Mb</DIV>
<DIV></DIV>
<DIV>Estimated total allocated dynamical RAM > 1807.14Mb</DIV>
<DIV>Generating pointlists ...</DIV>
<DIV>new r_m : 0.4125 (alat units) 0.7795 (a.u.) for type 1</DIV>
<DIV>new r_m : 0.4125 (alat units) 0.7795 (a.u.) for type 2</DIV>
<DIV>new r_m : 0.3193 (alat units) 0.6034 (a.u.) for type 3</DIV>
<DIV>new r_m : 0.3193 (alat units) 0.6034 (a.u.) for type 4</DIV>
<DIV>new r_m : 0.4125 (alat units) 0.7795 (a.u.) for type 5</DIV>
<DIV></DIV>
<DIV>Initial potential from superposition of free atoms</DIV>
<DIV></DIV>
<DIV>starting charge 463.99961, renormalised to 488.00000</DIV>
<DIV>Number of +U iterations with fixed ns = 0</DIV>
<DIV>Starting occupations:</DIV>
<DIV>--- enter write_ns ---</DIV>
<DIV>LDA+U parameters:</DIV>
<DIV>U( 1) = 4.00000000</DIV>
<DIV>alpha( 1) = 0.00000000</DIV>
<DIV>U( 2) = 4.00000000</DIV>
<DIV>alpha( 2) = 0.00000000</DIV>
<DIV>atom 81 Tr[ns(na)] (up, down, total) = 5.00000 1.00000 6.00000</DIV>
<DIV></DIV>
<DIV>....................</DIV>
<DIV>iteration # 17 ecut= 80.00 Ry beta=0.30</DIV>
<DIV>Davidson diagonalization with overlap</DIV>
<DIV>ethr = 6.52E-11, avg # of iterations = 1.5</DIV>
<DIV></DIV>
<DIV>total cpu time spent up to now is 357.8 secs</DIV>
<DIV></DIV>
<DIV>total energy = -3574.32866175 Ry</DIV>
<DIV>Harris-Foulkes estimate = -3574.32866175 Ry</DIV>
<DIV>estimated scf accuracy < 0.00000010 Ry</DIV>
<DIV></DIV>
<DIV>total magnetization = 0.00 Bohr mag/cell</DIV>
<DIV>absolute magnetization = 38.76 Bohr mag/cell</DIV>
<DIV></DIV>
<DIV>iteration # 18 ecut= 80.00 Ry beta=0.30</DIV>
<DIV>Davidson diagonalization with overlap</DIV>
<DIV>ethr = 2.10E-11, avg # of iterations = 1.5</DIV>
<DIV></DIV>
<DIV>total cpu time spent up to now is 374.7 secs</DIV>
<DIV></DIV>
<DIV>total energy = -3574.32866176 Ry</DIV>
<DIV>Harris-Foulkes estimate = -3574.32866177 Ry</DIV>
<DIV>estimated scf accuracy < 0.00000004 Ry</DIV>
<DIV></DIV>
<DIV>total magnetization = -0.00 Bohr mag/cell</DIV>
<DIV>absolute magnetization = 38.76 Bohr mag/cell</DIV>
<DIV></DIV>
<DIV>iteration # 19 ecut= 80.00 Ry beta=0.30</DIV>
<DIV>Davidson diagonalization with overlap</DIV>
<DIV>ethr = 7.92E-12, avg # of iterations = 1.5</DIV>
<DIV></DIV>
<DIV>total cpu time spent up to now is 392.6 secs</DIV>
<DIV></DIV>
<DIV>total energy = -3574.32866236 Ry</DIV>
<DIV>Harris-Foulkes estimate = -3574.32866177 Ry</DIV>
<DIV>estimated scf accuracy < 0.00000002 Ry</DIV>
<DIV></DIV>
<DIV>total magnetization = 0.00 Bohr mag/cell</DIV>
<DIV>absolute magnetization = 38.76 Bohr mag/cell</DIV>
<DIV></DIV>
<DIV>iteration # 20 ecut= 80.00 Ry beta=0.30</DIV>
<DIV>Davidson diagonalization with overlap</DIV>
<DIV>ethr = 3.31E-12, avg # of iterations = 2.0</DIV>
<DIV></DIV>
<DIV>Magnetic moment per site:</DIV>
<DIV>atom: 1 charge: 0.2215 magn: -0.0004 constr: 0.0000</DIV>
<DIV>atom: 2 charge: 0.3887 magn: -0.0004 constr: 0.0000</DIV>
<DIV>atom: 3 charge: 0.1459 magn: 0.0006 constr: 0.0000</DIV>
<DIV>...........</DIV>
<DIV>atom: 93 charge: 0.0564 magn: 0.0004 constr: 0.0000</DIV>
<DIV>atom: 94 charge: 0.0564 magn: -0.0004 constr: 0.0000</DIV>
<DIV>atom: 95 charge: 0.0564 magn: 0.0004 constr: 0.0000</DIV>
<DIV>atom: 96 charge: 0.0564 magn: -0.0004 constr: 0.0000</DIV>
<DIV></DIV>
<DIV>total cpu time spent up to now is 413.1 secs</DIV>
<DIV></DIV>
<DIV>End of self-consistent calculation</DIV>
<DIV>--- enter write_ns ---</DIV>
<DIV>LDA+U parameters:</DIV>
<DIV>U( 1) = 4.00000000</DIV>
<DIV>alpha( 1) = 0.00000000</DIV>
<DIV>U( 2) = 4.00000000</DIV>
<DIV>alpha( 2) = 0.00000000</DIV>
<DIV>atom 81 Tr[ns(na)] (up, down, total) = 4.98183 0.96635 5.94818</DIV>
<DIV>spin 1</DIV>
<DIV>eigenvalues: </DIV>
<DIV>0.985 0.998 0.999 0.999 1.001</DIV>
<DIV>eigenvectors:</DIV>
<DIV>0.205 0.012 0.467 0.309 0.008</DIV>
<DIV>0.575 0.171 0.105 0.006 0.143</DIV>
<DIV>0.006 0.054 0.424 0.410 0.105</DIV>
<DIV>0.184 0.571 0.002 0.214 0.029</DIV>
<DIV>0.030 0.192 0.002 0.060 0.715</DIV>
<DIV>occupations:</DIV>
<DIV>0.996 0.005 0.000 -0.003 0.001</DIV>
<DIV>0.005 0.991 -0.001 0.005 -0.003</DIV>
<DIV>0.000 -0.001 0.999 0.000 0.001</DIV>
<DIV>-0.003 0.005 0.000 0.996 0.000</DIV>
<DIV>0.001 -0.003 0.001 0.000 1.000</DIV>
<DIV>spin 2</DIV>
<DIV>eigenvalues: </DIV>
<DIV>0.092 0.109 0.131 0.300 0.333</DIV>
<DIV>eigenvectors:</DIV>
<DIV>0.006 0.745 0.046 0.012 0.191</DIV>
<DIV>0.124 0.075 0.001 0.287 0.513</DIV>
<DIV>0.086 0.070 0.608 0.230 0.005</DIV>
<DIV>0.571 0.110 0.003 0.092 0.225</DIV>
<DIV>0.214 0.000 0.340 0.380 0.066</DIV>
<DIV>occupations:</DIV>
<DIV>0.155 -0.060 -0.021 0.051 -0.036</DIV>
<DIV>-0.060 0.277 -0.039 -0.050 -0.025</DIV>
<DIV>-0.021 -0.039 0.166 -0.038 0.048</DIV>
<DIV>0.051 -0.050 -0.038 0.167 -0.070</DIV>
<DIV>-0.036 -0.025 0.048 -0.070 0.200</DIV>
<DIV>atomic mag. moment = 4.015476</DIV>
<DIV>atom 82 Tr[ns(na)] (up, down, total) = 0.96635 4.98183 5.94818</DIV>
<DIV>spin 1</DIV>
<DIV>eigenvalues: </DIV>
<DIV>0.092 0.109 0.131 0.300 0.333</DIV>
<DIV>eigenvectors:</DIV>
<DIV>0.006 0.745 0.046 0.012 0.191</DIV>
<DIV>0.124 0.075 0.001 0.287 0.513</DIV>
<DIV>0.086 0.070 0.608 0.230 0.005</DIV>
<DIV>0.571 0.110 0.003 0.092 0.225</DIV>
<DIV>0.214 0.000 0.340 0.380 0.066</DIV>
<DIV>occupations:</DIV>
<DIV>0.155 0.060 0.021 0.051 -0.036</DIV>
<DIV>0.060 0.277 -0.039 0.050 0.025</DIV>
<DIV>0.021 -0.039 0.166 0.038 -0.048</DIV>
<DIV>0.051 0.050 0.038 0.167 -0.070</DIV>
<DIV>-0.036 0.025 -0.048 -0.070 0.200</DIV>
<DIV>spin 2</DIV>
<DIV>eigenvalues: </DIV>
<DIV>0.985 0.998 0.999 0.999 1.001</DIV>
<DIV>.................</DIV>
<DIV>-6.1506 -6.0100 -5.9746 -5.7537 -5.7489 -5.6929 -5.5830 -5.4915</DIV>
<DIV>-5.4638 -5.4261 -5.3641 -5.3519 -5.3515 -5.3329 -5.2918 -5.2795</DIV>
<DIV>-5.2782 -5.2718 -5.1972 -5.1310 -4.9251 -4.9128 -4.8267 -4.8149</DIV>
<DIV>-3.9884 -3.9601 -3.9176 -3.8239 -3.7443 -3.7323 -3.6190 -3.5331</DIV>
<DIV>-3.4503 -3.3435 -3.2390 -3.2152 -3.1793 -3.1334 -3.0571 -3.0071</DIV>
<DIV>-2.9992 -2.9942 -2.9576 -2.9506 -2.9240 -2.8794 -2.8003 -2.6721</DIV>
<DIV>-2.6093 -2.3990 -2.3176 -2.1578 -2.1345 -2.1299 -2.1180 -2.0880</DIV>
<DIV>-2.0019 -1.9958 -1.9670 -1.9442 -1.8988 -1.8461 -1.7678 -1.7658</DIV>
<DIV>-1.6268 -1.6204 -1.4648 -1.4644 -1.3709 -1.3146 -1.2969 -1.2893</DIV>
<DIV>-1.0651 -0.9563 -0.8390 -0.6356 -0.3698 -0.2651 -0.1795 -0.1209</DIV>
<DIV>-0.0979 -0.0824 -0.0816 0.0090 0.0106 0.2117 0.2135 0.2195</DIV>
<DIV>0.2258 0.3881 0.4109 0.4381 0.4452 0.5491 0.5774 0.5855</DIV>
<DIV>0.7251 0.7344 0.7866 0.7983 0.7984 0.8264 0.8884 0.9250</DIV>
<DIV>1.0197 1.0488 1.0646 1.1063 1.1142 1.1153 1.1215 1.1977</DIV>
<DIV>1.2414 1.2450 1.2483 1.2517 1.2550 1.3231 1.3756 1.4053</DIV>
<DIV>1.4274 1.4358 1.4941 1.5206 1.5743 1.6436 1.7523 1.7681</DIV>
<DIV>1.7780 1.8048 1.9017 1.9171 2.0213 2.0253 2.0305 2.0399</DIV>
<DIV>2.0630 2.1295 2.2009 2.2122 2.2299 2.2892 2.2938 2.3264</DIV>
<DIV>2.3487 2.3516 2.3573 2.3775 2.4053 2.4189 2.4447 2.5663</DIV>
<DIV>2.6084 2.6110 2.7082 2.7380 2.7565 2.7582 2.8471 2.8472</DIV>
<DIV>2.8887 2.9200 3.1226 3.1265 5.5252 5.5621 5.5840 5.5862</DIV>
<DIV>5.5908 5.6155 5.6449 5.6461 5.6554 5.6579 5.6796 5.6835</DIV>
<DIV>6.2028 6.2194 6.2727 6.3047 6.3640 6.3754 6.4399 6.4874</DIV>
<DIV>6.8382 8.3784 8.4746 8.7730 8.8824 8.9521 8.9907 9.5771</DIV>
<DIV>10.1566 10.2854 10.3246 10.4288 10.5018 10.5770 10.6200 10.6638</DIV>
<DIV>10.7170 10.7483 10.8761 10.8840 10.9615 11.0368 11.0472 11.1268</DIV>
<DIV>11.2449 11.3329 11.4499 11.5404 11.5712</DIV>
<DIV></DIV>
<DIV>------ SPIN DOWN ----------</DIV>
<DIV></DIV>
<DIV></DIV>
<DIV>k = 0.0000 0.0000 0.0000 ( 37420 PWs) bands (ev):</DIV>
<DIV></DIV>
<DIV>-84.1635 -84.1635 -84.1634 -84.1634 -79.6195 -79.6195 -79.6195
-79.6195</DIV>
<DIV>-52.4085 -52.4085 -52.4085 -52.4085 -52.3860 -52.3860 -52.3860
-52.3860</DIV>
<DIV>-52.3295 -52.3295 -52.3295 -52.3295 -47.6583 -47.6583 -47.6582
-47.6582</DIV>
<DIV>-47.6479 -47.6479 -47.6479 -47.6479 -47.6347 -47.6347 -47.6347
-47.6347</DIV>
<DIV>-17.8187 -17.7820 -17.7663 -17.7536 -17.7150 -17.6785 -17.6639
-17.6483</DIV>
<DIV>-16.6495 -16.5036 -16.4882 -16.4872 -16.4539 -16.4148 -16.3778
-16.3688</DIV>
<DIV>-15.6824 -15.6689 -15.5867 -15.5740 -15.5386 -15.5119 -15.4054
-15.3767</DIV>
<DIV>-14.5944 -14.4562 -14.4204 -14.3954 -14.3936 -14.2819 -14.2445
-14.2387</DIV>
<DIV>-13.9761 -13.9034 -13.8879 -13.8614 -13.8128 -13.7619 -13.7554
-13.7384</DIV>
<DIV>-13.7233 -13.6685 -13.6570 -13.6287 -13.5615 -13.5376 -13.5265
-13.5034</DIV>
<DIV>-6.1507 -6.0100 -5.9746 -5.7537 -5.7489 -5.6929 -5.5830 -5.4915</DIV>
<DIV>-5.4638 -5.4262 -5.3641 -5.3519 -5.3515 -5.3329 -5.2918 -5.2795</DIV>
<DIV>-5.2782 -5.2718 -5.1971 -5.1310 -4.9251 -4.9128 -4.8267 -4.8149</DIV>
<DIV>-3.9885 -3.9601 -3.9176 -3.8239 -3.7443 -3.7323 -3.6190 -3.5331</DIV>
<DIV>-3.4503 -3.3435 -3.2390 -3.2152 -3.1793 -3.1334 -3.0571 -3.0071</DIV>
<DIV>-2.9992 -2.9942 -2.9576 -2.9506 -2.9240 -2.8794 -2.8003 -2.6721</DIV>
<DIV>-2.6093 -2.3990 -2.3176 -2.1578 -2.1345 -2.1299 -2.1180 -2.0880</DIV>
<DIV>-2.0019 -1.9958 -1.9670 -1.9442 -1.8988 -1.8461 -1.7678 -1.7658</DIV>
<DIV>-1.6267 -1.6204 -1.4648 -1.4644 -1.3709 -1.3146 -1.2969 -1.2893</DIV>
<DIV>-1.0651 -0.9563 -0.8390 -0.6356 -0.3698 -0.2651 -0.1795 -0.1209</DIV>
<DIV>-0.0979 -0.0824 -0.0816 0.0090 0.0106 0.2117 0.2135 0.2195</DIV>
<DIV>0.2258 0.3881 0.4109 0.4381 0.4452 0.5491 0.5774 0.5855</DIV>
<DIV>0.7251 0.7344 0.7866 0.7983 0.7984 0.8264 0.8884 0.9250</DIV>
<DIV>1.0197 1.0488 1.0646 1.1063 1.1142 1.1153 1.1215 1.1977</DIV>
<DIV>1.2414 1.2450 1.2483 1.2517 1.2550 1.3231 1.3756 1.4053</DIV>
<DIV>1.4274 1.4358 1.4942 1.5205 1.5743 1.6436 1.7523 1.7681</DIV>
<DIV>1.7780 1.8048 1.9017 1.9171 2.0213 2.0253 2.0305 2.0399</DIV>
<DIV>2.0630 2.1295 2.2009 2.2122 2.2299 2.2892 2.2938 2.3264</DIV>
<DIV>2.3487 2.3516 2.3573 2.3775 2.4053 2.4189 2.4447 2.5663</DIV>
<DIV>2.6084 2.6110 2.7082 2.7380 2.7565 2.7582 2.8471 2.8472</DIV>
<DIV>2.8887 2.9199 3.1226 3.1264 5.5252 5.5621 5.5840 5.5863</DIV>
<DIV>5.5908 5.6155 5.6449 5.6462 5.6555 5.6579 5.6796 5.6836</DIV>
<DIV>6.2028 6.2194 6.2727 6.3047 6.3640 6.3754 6.4399 6.4874</DIV>
<DIV>6.8382 8.3784 8.4746 8.7730 8.8824 8.9521 8.9907 9.5771</DIV>
<DIV>10.1566 10.2854 10.3246 10.4288 10.5018 10.5770 10.6200 10.6639</DIV>
<DIV>10.7170 10.7483 10.8761 10.8840 10.9615 11.0368 11.0472 11.1268</DIV>
<DIV>11.2449 11.3329 11.4499 11.5404 11.5710</DIV>
<DIV></DIV>
<DIV>the Fermi energy is 3.2355 ev</DIV>
<DIV></DIV>
<DIV>! total energy = -3574.32866236 Ry</DIV>
<DIV>Harris-Foulkes estimate = -3574.32866236 Ry</DIV>
<DIV>estimated scf accuracy < 3.0E-09 Ry</DIV>
<DIV></DIV>
<DIV>The total energy is the sum of the following terms:</DIV>
<DIV></DIV>
<DIV>one-electron contribution = -2562.56041226 Ry</DIV>
<DIV>hartree contribution = 1464.78176945 Ry</DIV>
<DIV>xc contribution = -580.13791183 Ry</DIV>
<DIV>ewald contribution = -1896.58030792 Ry</DIV>
<DIV>Dispersion Correction = -0.75249209 Ry</DIV>
<DIV>Hubbard energy = 0.87642701 Ry</DIV>
<DIV>smearing contrib. (-TS) = 0.04426529 Ry</DIV>
<DIV></DIV>
<DIV>total magnetization = -0.00 Bohr mag/cell</DIV>
<DIV>absolute magnetization = 38.76 Bohr mag/cell</DIV>
<DIV></DIV>
<DIV>convergence has been achieved in 20 iterations</DIV>
<DIV></DIV>
<DIV></DIV>
<DIV></DIV>
<DIV></DIV>
<DIV></DIV>
<DIV></DIV>
<DIV></DIV>
<DIV></DIV>
<DIV></DIV>
<DIV></DIV>
<DIV>MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .</DIV>
<DIV></DIV>
<DIV>MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .</DIV>
<DIV></DIV>
<DIV>MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .</DIV>
<DIV></DIV>
<DIV>MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .</DIV>
<DIV></DIV>
<DIV>MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .</DIV>
<DIV></DIV>
<DIV>MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .</DIV>
<DIV></DIV>
<DIV>MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .</DIV>
<DIV></DIV>
<DIV>MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .</DIV>
<DIV></DIV>
<DIV>MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .</DIV>
<DIV></DIV>
<DIV>MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .</DIV>
<DIV></DIV>
<DIV>MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .</DIV>
<DIV></DIV>
<DIV>MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .</DIV>
<DIV></DIV>
<DIV>MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .</DIV>
<DIV></DIV>
<DIV>MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .</DIV>
<DIV></DIV>
<DIV>MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .</DIV>
<DIV></DIV>
<DIV>MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .</DIV>
<DIV></DIV>
<DIV>MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .</DIV>
<DIV></DIV>
<DIV>MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .</DIV>
<DIV></DIV>
<DIV>MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .</DIV>
<DIV></DIV>
<DIV>MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .</DIV>
<DIV></DIV>
<DIV>MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .</DIV>
<DIV></DIV>
<DIV>MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .</DIV>
<DIV></DIV>
<DIV>MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .</DIV>
<DIV></DIV>
<DIV>MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .</DIV>
<DIV></DIV>
<DIV>MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .</DIV>
<DIV></DIV>
<DIV>MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .</DIV>
<DIV></DIV>
<DIV>MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .</DIV>
<DIV>................</DIV>
<DIV> </DIV>
<DIV> </DIV>
<DIV><FONT size=4></FONT></DIV>
<DIV><FONT size=4></FONT></DIV>
<DIV><FONT size=4></FONT></DIV>
<DIV><FONT size=4>By the way, this error only appear for Norm-conserving
pseudopotentials.</FONT></DIV>
<DIV><FONT size=4>Can anyone help me to solve this problem?</FONT></DIV>
<DIV><FONT size=4></FONT></DIV>
<DIV><FONT size=4>Regards,</FONT></DIV>
<DIV><FONT size=4></FONT></DIV>
<DIV><FONT size=4>Shaofeng </FONT></DIV>
<DIV
style="FONT-FAMILY: ; COLOR: ">--------------------------------------<BR>Shaofeng
Wang, Ph.D of Geochemistry<BR>Environmental Molecular Science Group<BR>Institute
of Applied Ecology, Chinese Academy of Sciences<BR>Shenyang, 110016,
China<BR>wangshaofeng@iae.ac.cn<BR>www.iae.cas.cn<BR><BR></DIV></DIV>
<DIV
style="FONT-FAMILY: 'Century'; COLOR: #000000; FONT-SIZE: 14pt">--------------------------------------<BR>Shaofeng
Wang, Ph.D of Geochemistry<BR>Environmental Molecular Science Group<BR>Institute
of Applied Ecology, Chinese Academy of Sciences<BR>Shenyang, 110016,
China<BR>wangshaofeng@iae.ac.cn<BR>www.iae.cas.cn<BR><BR></DIV></DIV></DIV></BODY></HTML>