<div dir="ltr">Dear QE users,<div><br></div><div>As in the attached portion of pw.x output, I'm getting NaN forces. My cluster is Cray XE (<a href="https://bluewaters.ncsa.illinois.edu/">https://bluewaters.ncsa.illinois.edu/</a>, compiled after ./install/configure ARCH=crayxt --enable-openmp --enable-parallel --with-scalapack --with-elpa)</div><div><br></div><div>This only happens if I activate lda_plus_u= .true. I've tried without elpa and opnemp, but it didn't change.</div><div><br></div><div>Has anyone encountered this before ?</div><div><br></div><div>Any suggestion is appreciated.</div><div><br></div><div>Thanks.</div><div><br></div><div><br></div><div>-------------</div><div><br></div><div><div>&control</div><div>calculation ='scf',</div><div>tstress = .true.,</div><div>tprnfor = .true.,</div><div>pseudo_dir ='./',</div><div>outdir='./tmp',</div><div>disk_io = 'none'</div><div>/</div><div><br></div><div>&system</div><div>ibrav=0,</div><div>celldm(1)=1.889725,</div><div>nat = 40,</div><div>ntyp = 4,</div><div>ecutwfc = 50.0, ecutrho = 200.0,</div><div>nspin= 2,</div><div>starting_magnetization(1)= 0.25,</div><div>occupations= 'smearing', smearing= 'm-p', degauss= 0.005,</div><div>lda_plus_u= .true.</div><div>Hubbard_U(1)= 3.1</div><div>/</div><div><br></div><div>&electrons</div><div> conv_thr = 1.0D-8,</div><div> mixing_mode = 'plain',</div><div> mixing_beta = 0.3,</div><div>/</div><div><br></div><div>&ions</div><div>/</div><div><br></div><div>&cell</div><div>/</div><div><br></div><div>ATOMIC_SPECIES</div><div>Fe 1.00 fe_lda_v1.2.uspp.F.UPF</div><div>Mg 1.00 mg_lda_v1.uspp.F.UPF</div><div>Ge 1.00 ge_lda_v1.uspp.F.UPF</div><div>O 1.00 o_lda_v1.2.uspp.F.UPF</div><div><br></div><div>K_POINTS {automatic}</div><div>2 2 2 1 1 1</div><div><br></div><div>CELL_PARAMETERS</div><div> 5.1049880824 4.9672323236 0.0000000000</div><div> -4.7912675141 5.2733257368 -0.0000000000</div><div> -0.0000000000 0.0000000000 7.2032409922</div><div> </div><div>ATOMIC_POSITIONS angstrom</div><div> Fe 0.0640761766 5.4327565722 1.8008102480 </div><div> Mg -2.2988644562 2.9398516435 1.8008102480 </div><div> Mg 0.2241918836 4.8130480598 5.4024307441 </div><div> Mg 2.6306507132 7.3060445959 5.4024307441 </div><div> Mg 0.2406224949 10.5429355305 1.8008102480 </div><div> Mg -2.1553695100 8.0528043965 1.8008102480 </div><div> Mg 0.3775479659 9.9332538514 5.4024307441 </div><div> Mg -2.3207876196 7.4557439726 5.4024307441 </div><div> Ge 0.2387789083 7.6817096680 -0.0003661912 </div><div> Ge 2.6317776574 5.0436724103 0.0005773902 </div><div> Ge 0.2387789083 7.6817096680 3.6019866873 </div><div> Ge 2.6317776574 5.0436724103 3.6010431059 </div><div> Ge 0.0847798457 2.5566994682 0.0002034988 </div><div> Ge -2.3230677758 5.1976250811 0.0020760533 </div><div> Ge 0.0847798457 2.5566994682 3.6014169973 </div><div> Ge -2.3230677758 5.1976250811 3.5995444428 </div><div> O -0.8125582446 4.1476235583 0.4523679571 </div><div> O -0.8125582446 4.1476235583 3.1492525390 </div><div> O 2.0480850306 4.7754591319 1.8008102480 </div><div> O 0.8365953645 7.4123247931 1.8008102480 </div><div> O -1.3429791197 6.7381611727 0.4456839473 </div><div> O -1.3429791197 6.7381611727 3.1559365488 </div><div> O 1.1226742607 6.0903006692 -0.4575712326 </div><div> O 1.1226742607 6.0903006692 4.0591917287 </div><div> O -2.9299032201 4.9437783857 1.8008102480 </div><div> O 3.6144962241 6.5785289829 0.4534922390 </div><div> O 1.8071826416 8.6303200619 -0.4579511820 </div><div> O 3.6144962241 6.5785289829 3.1481282571 </div><div> O 1.8071826416 8.6303200619 4.0595716781 </div><div> O -0.5067970200 2.8487235832 5.4024307441 </div><div> O 3.2347126384 5.2970177159 5.4024307441 </div><div> O -3.8257139978 6.2420006036 -0.4587575992 </div><div> O 1.6485222082 3.5119019151 -0.4562875575 </div><div> O -0.6456714568 9.2738705978 0.4483821013 </div><div> O -3.8257139978 6.2420006036 4.0603780953 </div><div> O 1.6485222082 3.5119019151 4.0579080536 </div><div> O -0.6456714568 9.2738705978 3.1532383948 </div><div> O 0.6676050140 2.2631227748 1.8008102480 </div><div> O -0.3514415341 7.9676535020 5.4024307441 </div><div> O -1.7112214092 5.4465820124 5.4024307441 </div></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><div> Program PWSCF v.5.1 starts on 8Sep2014 at 13:28:24 </div><div><br></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><br></div><div> Parallel version (MPI & OpenMP), running on 32 processor cores</div><div> Number of MPI processes: 32</div><div> Threads/MPI process: 1</div><div> K-points division: npool = 2</div><div> R & G space division: proc/nbgrp/npool/nimage = 16</div><div> Waiting for input...</div><div> Reading input from standard input</div><div>Warning: card &IONS ignored</div><div>Warning: card / ignored</div><div>Warning: card &CELL ignored</div><div>Warning: card / ignored</div><div> Message from routine read_cards :</div><div> DEPRECATED: no units specified in CELL_PARAMETERS card</div><div><br></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><br></div><div> Subspace diagonalization in iterative solution of the eigenvalue problem:</div><div> ELPA distributed-memory algorithm (size of sub-group: 2* 2 procs)</div><div><br></div><div> </div><div> Parallelization info</div><div> --------------------</div><div> sticks: dense smooth PW G-vecs: dense smooth PW</div><div> Min 180 180 47 7358 7358 1001</div><div> Max 181 181 48 7361 7361 1006</div><div> Sum 2885 2885 767 117745 117745 16061</div><div> </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.4125 (alat units) 0.7795 (a.u.) for type 3</div><div> new r_m : 0.4125 (alat units) 0.7795 (a.u.) for type 4</div><div><br></div><div><br></div><div> bravais-lattice index = 0</div><div> lattice parameter (alat) = 1.8897 a.u.</div><div> unit-cell volume = 2465.4694 (a.u.)^3</div><div> number of atoms/cell = 40</div><div> number of atomic types = 4</div><div> number of electrons = 342.00</div><div> number of Kohn-Sham states= 205</div><div> kinetic-energy cutoff = 50.0000 Ry</div><div> charge density cutoff = 200.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)</div><div><br></div><div> celldm(1)= 1.889725 celldm(2)= 0.000000 celldm(3)= 0.000000</div><div> celldm(4)= 0.000000 celldm(5)= 0.000000 celldm(6)= 0.000000</div><div><br></div><div> crystal axes: (cart. coord. in units of alat)</div><div> a(1) = ( 5.104988 4.967232 0.000000 ) </div><div> a(2) = ( -4.791268 5.273326 0.000000 ) </div><div> a(3) = ( 0.000000 0.000000 7.203241 ) </div><div><br></div><div> reciprocal axes: (cart. coord. in units 2 pi/alat)</div><div> b(1) = ( 0.103970 0.094466 0.000000 ) </div><div> b(2) = ( -0.097935 0.100651 0.000000 ) </div><div> b(3) = ( 0.000000 0.000000 0.138826 ) </div><div><br></div><div><br></div><div> PseudoPot. # 1 for Fe read from file:</div><div> ./fe_lda_v1.2.uspp.F.UPF</div><div> MD5 check sum: 32cbb033663340320c88e12ca63afbaf</div><div> Pseudo is Ultrasoft + core correction, Zval = 16.0</div><div> Generated by new atomic code, or converted to UPF format</div><div> Using radial grid of 861 points, 6 beta functions with: </div><div> l(1) = 0</div><div> l(2) = 0</div><div> l(3) = 1</div><div> l(4) = 1</div><div> l(5) = 2</div><div> l(6) = 2</div><div> Q(r) pseudized with 8 coefficients, rinner = 1.200 1.200 1.200</div><div> 1.200 1.200</div><div><br></div><div> PseudoPot. # 2 for Mg read from file:</div><div> ./mg_lda_v1.uspp.F.UPF</div><div> MD5 check sum: 3e2f2ee1761bb9e932e4e49f523ff5f9</div><div> Pseudo is Ultrasoft, Zval = 10.0</div><div> Generated by new atomic code, or converted to UPF format</div><div> Using radial grid of 821 points, 7 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> l(7) = 2</div><div> Q(r) pseudized with 8 coefficients, rinner = 1.200 1.200 1.200</div><div> 1.200 1.200</div><div><br></div><div> PseudoPot. # 3 for Ge read from file:</div><div> ./ge_lda_v1.uspp.F.UPF</div><div> MD5 check sum: 15ad864a5cf8594a2bbd64283e74b2b4</div><div> Pseudo is Ultrasoft + core correction, Zval = 14.0</div><div> Generated by new atomic code, or converted to UPF format</div><div> Using radial grid of 873 points, 6 beta functions with: </div><div> l(1) = 0</div><div> l(2) = 0</div><div> l(3) = 1</div><div> l(4) = 1</div><div> l(5) = 2</div><div> l(6) = 2</div><div> Q(r) pseudized with 8 coefficients, rinner = 1.400 1.400 1.400</div><div> 1.400 1.400</div><div><br></div><div> PseudoPot. # 4 for O read from file:</div><div> ./o_lda_v1.2.uspp.F.UPF</div><div> MD5 check sum: 194b818514d2a7dbb603bb78519a4f40</div><div> Pseudo is Ultrasoft + core correction, Zval = 6.0</div><div> Generated by new atomic code, or converted to UPF format</div><div> Using radial grid of 737 points, 5 beta functions with: </div><div> l(1) = 0</div><div> l(2) = 0</div><div> l(3) = 1</div><div> l(4) = 1</div><div> l(5) = 2</div><div> Q(r) pseudized with 8 coefficients, rinner = 0.900 0.900 0.900</div><div> 0.900 0.900</div><div><br></div><div> atomic species valence mass pseudopotential</div><div> Fe 16.00 1.00000 Fe( 1.00)</div><div> Mg 10.00 1.00000 Mg( 1.00)</div><div> Ge 14.00 1.00000 Ge( 1.00)</div><div> O 6.00 1.00000 O ( 1.00)</div><div><br></div><div> Starting magnetic structure </div><div> atomic species magnetization</div><div> Fe 0.250</div><div> Mg 0.000</div><div> Ge 0.000</div><div> O 0.000</div><div><br></div><div><br></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 3.1000 0.0000 0.0000 0.0000</div><div><br></div><div><br></div><div><br></div><div> 2 Sym. Ops. (no inversion) found ( 1 have fractional translation)</div><div><br></div><div><br></div><div><br></div><div> Cartesian axes</div><div><br></div><div> site n. atom positions (alat units)</div><div> 1 Fe tau( 1) = ( 0.0640762 5.4327598 1.8008113 )</div><div> 2 Mg tau( 2) = ( -2.2988658 2.9398534 1.8008113 )</div><div> 3 Mg tau( 3) = ( 0.2241920 4.8130509 5.4024340 )</div><div> 4 Mg tau( 4) = ( 2.6306523 7.3060490 5.4024340 )</div><div> 5 Mg tau( 5) = ( 0.2406226 10.5429419 1.8008113 )</div><div> 6 Mg tau( 6) = ( -2.1553708 8.0528092 1.8008113 )</div><div> 7 Mg tau( 7) = ( 0.3775482 9.9332598 5.4024340 )</div><div> 8 Mg tau( 8) = ( -2.3207890 7.4557484 5.4024340 )</div><div> 9 Ge tau( 9) = ( 0.2387791 7.6817143 -0.0003662 )</div><div> 10 Ge tau( 10) = ( 2.6317792 5.0436754 0.0005774 )</div><div> 11 Ge tau( 11) = ( 0.2387791 7.6817143 3.6019888 )</div><div> 12 Ge tau( 12) = ( 2.6317792 5.0436754 3.6010453 )</div><div> 13 Ge tau( 13) = ( 0.0847799 2.5567010 0.0002035 )</div><div> 14 Ge tau( 14) = ( -2.3230692 5.1976282 0.0020761 )</div><div> 15 Ge tau( 15) = ( 0.0847799 2.5567010 3.6014192 )</div><div> 16 Ge tau( 16) = ( -2.3230692 5.1976282 3.5995466 )</div><div> 17 O tau( 17) = ( -0.8125587 4.1476260 0.4523682 )</div><div> 18 O tau( 18) = ( -0.8125587 4.1476260 3.1492544 )</div><div> 19 O tau( 19) = ( 2.0480863 4.7754620 1.8008113 )</div><div> 20 O tau( 20) = ( 0.8365959 7.4123292 1.8008113 )</div><div> 21 O tau( 21) = ( -1.3429799 6.7381652 0.4456842 )</div><div> 22 O tau( 22) = ( -1.3429799 6.7381652 3.1559384 )</div><div> 23 O tau( 23) = ( 1.1226749 6.0903043 -0.4575715 )</div><div> 24 O tau( 24) = ( 1.1226749 6.0903043 4.0591942 )</div><div> 25 O tau( 25) = ( -2.9299050 4.9437813 1.8008113 )</div><div> 26 O tau( 26) = ( 3.6144984 6.5785329 0.4534925 )</div><div> 27 O tau( 27) = ( 1.8071837 8.6303252 -0.4579515 )</div><div> 28 O tau( 28) = ( 3.6144984 6.5785329 3.1481301 )</div><div> 29 O tau( 29) = ( 1.8071837 8.6303252 4.0595741 )</div><div> 30 O tau( 30) = ( -0.5067973 2.8487253 5.4024340 )</div><div> 31 O tau( 31) = ( 3.2347146 5.2970209 5.4024340 )</div><div> 32 O tau( 32) = ( -3.8257163 6.2420043 -0.4587579 )</div><div> 33 O tau( 33) = ( 1.6485232 3.5119040 -0.4562878 )</div><div> 34 O tau( 34) = ( -0.6456718 9.2738762 0.4483824 )</div><div> 35 O tau( 35) = ( -3.8257163 6.2420043 4.0603805 )</div><div> 36 O tau( 36) = ( 1.6485232 3.5119040 4.0579105 )</div><div> 37 O tau( 37) = ( -0.6456718 9.2738762 3.1532403 )</div><div> 38 O tau( 38) = ( 0.6676054 2.2631241 1.8008113 )</div><div> 39 O tau( 39) = ( -0.3514417 7.9676583 5.4024340 )</div><div> 40 O tau( 40) = ( -1.7112224 5.4465853 5.4024340 )</div><div><br></div><div> number of k points= 4 Methfessel-Paxton smearing, width (Ry)= 0.0050</div><div> cart. coord. in units 2pi/alat</div><div> k( 1) = ( 0.0015088 0.0487792 0.0347066), wk = 0.5000000</div><div> k( 2) = ( 0.0504763 -0.0015463 0.0347066), wk = 0.5000000</div><div> k( 3) = ( 0.0015088 0.0487792 0.0347066), wk = 0.5000000</div><div> k( 4) = ( 0.0504763 -0.0015463 0.0347066), wk = 0.5000000</div><div><br></div><div> Dense grid: 117745 G-vectors FFT dimensions: ( 64, 64, 64)</div><div><br></div><div> Largest allocated arrays est. size (Mb) dimensions</div><div> Kohn-Sham Wavefunctions 2.87 Mb ( 919, 205)</div><div> Atomic Hubbard wavefuncts 0.07 Mb ( 919, 5)</div><div> NL pseudopotentials 8.51 Mb ( 919, 607)</div><div> Each V/rho on FFT grid 0.50 Mb ( 16384, 2)</div><div> Each G-vector array 0.06 Mb ( 7361)</div><div> G-vector shells 0.03 Mb ( 3673)</div><div> Largest temporary arrays est. size (Mb) dimensions</div><div> Auxiliary wavefunctions 11.50 Mb ( 919, 820)</div><div> Each subspace H/S matrix 2.57 Mb ( 410, 410)</div><div> Each <psi_i|beta_j> matrix 1.90 Mb ( 607, 205)</div><div> Arrays for rho mixing 2.00 Mb ( 16384, 8)</div><div><br></div><div> Check: negative/imaginary core charge= -0.000017 0.000000</div><div><br></div><div> Initial potential from superposition of free atoms</div><div><br></div><div> starting charge 338.89287, renormalised to 342.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) = 3.10000000</div><div>alpha( 1) = 0.00000000</div><div>atom 1 Tr[ns(na)] (up, down, total) = 5.00000 1.00000 6.00000</div><div> spin 1</div><div> eigenvalues: </div><div> 1.000 1.000 1.000 1.000 1.000</div><div> eigenvectors:</div><div> 1.000 0.000 0.000 0.000 0.000</div><div> 0.000 1.000 0.000 0.000 0.000</div><div> 0.000 0.000 1.000 0.000 0.000</div><div> 0.000 0.000 0.000 1.000 0.000</div><div> 0.000 0.000 0.000 0.000 1.000</div><div> occupations:</div><div> 1.000 0.000 0.000 0.000 0.000</div><div> 0.000 1.000 0.000 0.000 0.000</div><div> 0.000 0.000 1.000 0.000 0.000</div><div> 0.000 0.000 0.000 1.000 0.000</div><div> 0.000 0.000 0.000 0.000 1.000</div><div> spin 2</div><div> eigenvalues: </div><div> 0.200 0.200 0.200 0.200 0.200</div><div> eigenvectors:</div><div> 1.000 0.000 0.000 0.000 0.000</div><div> 0.000 1.000 0.000 0.000 0.000</div><div> 0.000 0.000 1.000 0.000 0.000</div><div> 0.000 0.000 0.000 1.000 0.000</div><div> 0.000 0.000 0.000 0.000 1.000</div><div> occupations:</div><div> 0.200 0.000 0.000 0.000 0.000</div><div> 0.000 0.200 0.000 0.000 0.000</div><div> 0.000 0.000 0.200 0.000 0.000</div><div> 0.000 0.000 0.000 0.200 0.000</div><div> 0.000 0.000 0.000 0.000 0.200</div><div>atomic mag. moment = 4.000000</div><div>N of occupied +U levels = 6.000000</div><div> --- exit write_ns ---</div><div> Atomic wfc used for LDA+U Projector are NOT orthogonalized</div><div> Starting wfc are 216 randomized atomic wfcs</div><div><br></div><div> total cpu time spent up to now is 4.2 secs</div><div><br></div><div> per-process dynamical memory: 33.4 Mb</div><div><br></div><div> Self-consistent Calculation</div><div><br></div><div> iteration # 1 ecut= 50.00 Ry beta=0.30</div><div> Davidson diagonalization with overlap</div><div> ethr = 1.00E-02, avg # of iterations = 4.5</div><div> --- enter write_ns ---</div><div> LDA+U parameters:</div><div>U( 1) = 3.10000000</div><div>alpha( 1) = 0.00000000</div><div>atom 1 Tr[ns(na)] (up, down, total) = 4.99609 2.35813 7.35422</div><div> spin 1</div><div> eigenvalues: </div><div> 0.999 0.999 0.999 0.999 0.999</div><div> eigenvectors:</div><div> -0.473 0.000 0.881 -0.015 0.000</div><div> 0.000 0.747 0.000 0.000 -0.664</div><div> 0.000 0.664 0.000 0.000 0.747</div><div> 0.862 0.000 0.459 -0.214 0.000</div><div> 0.181 0.000 0.114 0.977 0.000</div><div> occupations:</div><div> 0.999 0.000 0.000 0.000 0.000</div><div> 0.000 0.999 0.000 0.000 0.000</div><div> 0.000 0.000 0.999 0.000 0.000</div><div> 0.000 0.000 0.000 0.999 0.000</div><div> 0.000 0.000 0.000 0.000 0.999</div><div> spin 2</div><div> eigenvalues: </div><div> 0.316 0.334 0.365 0.369 0.973</div><div> eigenvectors:</div><div> -0.989 0.000 0.000 -0.143 -0.031</div><div> 0.000 0.968 -0.250 0.000 0.000</div><div> 0.000 -0.250 -0.968 0.000 0.000</div><div> -0.109 0.000 0.000 0.579 0.808</div><div> 0.097 0.000 0.000 -0.803 0.589</div><div> occupations:</div><div> 0.318 0.000 0.000 -0.021 -0.006</div><div> 0.000 0.336 0.007 0.000 0.000</div><div> 0.000 0.007 0.363 0.000 0.000</div><div> -0.021 0.000 0.000 0.763 0.287</div><div> -0.006 0.000 0.000 0.287 0.578</div><div>atomic mag. moment = 2.637957</div><div>N of occupied +U levels = 7.354220</div><div> --- exit write_ns ---</div><div> </div><div> Magnetic moment per site:</div><div> atom: 1 charge: 7.3095 magn: 1.6861 constr: 0.0000</div><div> atom: 2 charge: 5.8857 magn: 0.0000 constr: 0.0000</div><div> atom: 3 charge: 5.8864 magn: 0.0000 constr: 0.0000</div><div> atom: 4 charge: 5.8820 magn: 0.0000 constr: 0.0000</div><div> atom: 5 charge: 5.8824 magn: 0.0000 constr: 0.0000</div><div> atom: 6 charge: 5.8814 magn: 0.0000 constr: 0.0000</div><div> atom: 7 charge: 5.8860 magn: 0.0000 constr: 0.0000</div><div> atom: 8 charge: 5.8852 magn: 0.0000 constr: 0.0000</div><div> atom: 9 charge: 6.8132 magn: 0.0002 constr: 0.0000</div><div> atom: 10 charge: 6.8140 magn: 0.0002 constr: 0.0000</div><div> atom: 11 charge: 6.8132 magn: 0.0002 constr: 0.0000</div><div> atom: 12 charge: 6.8140 magn: 0.0002 constr: 0.0000</div><div> atom: 13 charge: 6.8108 magn: 0.0002 constr: 0.0000</div><div> atom: 14 charge: 6.8109 magn: 0.0002 constr: 0.0000</div><div> atom: 15 charge: 6.8108 magn: 0.0002 constr: 0.0000</div><div> atom: 16 charge: 6.8109 magn: 0.0002 constr: 0.0000</div><div> atom: 17 charge: 2.0927 magn: 0.0171 constr: 0.0000</div><div> atom: 18 charge: 2.0927 magn: 0.0171 constr: 0.0000</div><div> atom: 19 charge: 2.0974 magn: 0.0144 constr: 0.0000</div><div> atom: 20 charge: 2.1066 magn: 0.0143 constr: 0.0000</div><div> atom: 21 charge: 2.0916 magn: 0.0136 constr: 0.0000</div><div> atom: 22 charge: 2.0916 magn: 0.0136 constr: 0.0000</div><div> atom: 23 charge: 2.0975 magn: 0.0092 constr: 0.0000</div><div> atom: 24 charge: 2.0975 magn: 0.0092 constr: 0.0000</div><div> atom: 25 charge: 2.1065 magn: 0.0011 constr: 0.0000</div><div> atom: 26 charge: 2.1003 magn: 0.0014 constr: 0.0000</div><div> atom: 27 charge: 2.0993 magn: 0.0023 constr: 0.0000</div><div> atom: 28 charge: 2.1003 magn: 0.0014 constr: 0.0000</div><div> atom: 29 charge: 2.0993 magn: 0.0023 constr: 0.0000</div><div> atom: 30 charge: 2.1039 magn: 0.0021 constr: 0.0000</div><div> atom: 31 charge: 2.1029 magn: 0.0021 constr: 0.0000</div><div> atom: 32 charge: 2.1006 magn: 0.0018 constr: 0.0000</div><div> atom: 33 charge: 2.1023 magn: 0.0011 constr: 0.0000</div><div> atom: 34 charge: 2.0971 magn: 0.0026 constr: 0.0000</div><div> atom: 35 charge: 2.1006 magn: 0.0018 constr: 0.0000</div><div> atom: 36 charge: 2.1023 magn: 0.0011 constr: 0.0000</div><div> atom: 37 charge: 2.0971 magn: 0.0026 constr: 0.0000</div><div> atom: 38 charge: 2.1058 magn: 0.0002 constr: 0.0000</div><div> atom: 39 charge: 2.1053 magn: 0.0012 constr: 0.0000</div><div> atom: 40 charge: 2.1065 magn: 0.0014 constr: 0.0000</div><div><br></div><div> total cpu time spent up to now is 13.4 secs</div><div><br></div><div> total energy = -3560.92084665 Ry</div><div> Harris-Foulkes estimate = -3564.91160194 Ry</div><div> estimated scf accuracy < 9.18723739 Ry</div><div><br></div><div> total magnetization = 4.03 Bohr mag/cell</div><div> absolute magnetization = 4.03 Bohr mag/cell</div><div><br></div><div> iteration # 2 ecut= 50.00 Ry beta=0.30</div><div> Davidson diagonalization with overlap</div><div> ethr = 2.69E-03, avg # of iterations = 3.5</div><div> </div><div> Magnetic moment per site:</div><div> atom: 1 charge: 7.1778 magn: 1.6765 constr: 0.0000</div><div> atom: 2 charge: 5.8732 magn: 0.0001 constr: 0.0000</div><div> atom: 3 charge: 5.8737 magn: 0.0000 constr: 0.0000</div><div> atom: 4 charge: 5.8695 magn: 0.0000 constr: 0.0000</div><div> atom: 5 charge: 5.8700 magn: 0.0000 constr: 0.0000</div><div> atom: 6 charge: 5.8690 magn: 0.0000 constr: 0.0000</div><div> atom: 7 charge: 5.8735 magn: 0.0000 constr: 0.0000</div><div> atom: 8 charge: 5.8728 magn: 0.0000 constr: 0.0000</div><div> atom: 9 charge: 6.8057 magn: 0.0001 constr: 0.0000</div><div> atom: 10 charge: 6.8068 magn: 0.0002 constr: 0.0000</div><div> atom: 11 charge: 6.8057 magn: 0.0001 constr: 0.0000</div><div> atom: 12 charge: 6.8068 magn: 0.0002 constr: 0.0000</div><div> atom: 13 charge: 6.8038 magn: 0.0001 constr: 0.0000</div><div> atom: 14 charge: 6.8037 magn: 0.0001 constr: 0.0000</div><div> atom: 15 charge: 6.8038 magn: 0.0001 constr: 0.0000</div><div> atom: 16 charge: 6.8037 magn: 0.0001 constr: 0.0000</div><div> atom: 17 charge: 2.1006 magn: 0.0167 constr: 0.0000</div><div> atom: 18 charge: 2.1006 magn: 0.0167 constr: 0.0000</div><div> atom: 19 charge: 2.1022 magn: 0.0140 constr: 0.0000</div><div> atom: 20 charge: 2.1130 magn: 0.0139 constr: 0.0000</div><div> atom: 21 charge: 2.1016 magn: 0.0135 constr: 0.0000</div><div> atom: 22 charge: 2.1016 magn: 0.0135 constr: 0.0000</div><div> atom: 23 charge: 2.1036 magn: 0.0090 constr: 0.0000</div><div> atom: 24 charge: 2.1036 magn: 0.0090 constr: 0.0000</div><div> atom: 25 charge: 2.1040 magn: 0.0009 constr: 0.0000</div><div> atom: 26 charge: 2.1016 magn: 0.0012 constr: 0.0000</div><div> atom: 27 charge: 2.1036 magn: 0.0019 constr: 0.0000</div><div> atom: 28 charge: 2.1016 magn: 0.0012 constr: 0.0000</div><div> atom: 29 charge: 2.1036 magn: 0.0019 constr: 0.0000</div><div> atom: 30 charge: 2.1064 magn: 0.0015 constr: 0.0000</div><div> atom: 31 charge: 2.1055 magn: 0.0016 constr: 0.0000</div><div> atom: 32 charge: 2.1035 magn: 0.0012 constr: 0.0000</div><div> atom: 33 charge: 2.1022 magn: 0.0008 constr: 0.0000</div><div> atom: 34 charge: 2.1014 magn: 0.0022 constr: 0.0000</div><div> atom: 35 charge: 2.1035 magn: 0.0012 constr: 0.0000</div><div> atom: 36 charge: 2.1022 magn: 0.0008 constr: 0.0000</div><div> atom: 37 charge: 2.1014 magn: 0.0022 constr: 0.0000</div><div> atom: 38 charge: 2.1010 magn: 0.0001 constr: 0.0000</div><div> atom: 39 charge: 2.1062 magn: 0.0010 constr: 0.0000</div><div> atom: 40 charge: 2.1060 magn: 0.0010 constr: 0.0000</div><div><br></div><div> total cpu time spent up to now is 21.7 secs</div><div><br></div><div> total energy = -3562.05386923 Ry</div><div> Harris-Foulkes estimate = -3562.38770592 Ry</div><div> estimated scf accuracy < 0.92880777 Ry</div><div><br></div><div> total magnetization = 4.02 Bohr mag/cell</div><div> absolute magnetization = 4.07 Bohr mag/cell</div><div><br></div><div> iteration # 3 ecut= 50.00 Ry beta=0.30</div><div> Davidson diagonalization with overlap</div><div> ethr = 2.72E-04, avg # of iterations = 3.5</div><div> </div><div> Magnetic moment per site:</div><div> atom: 1 charge: 7.1892 magn: 1.6799 constr: 0.0000</div><div> atom: 2 charge: 5.8681 magn: 0.0000 constr: 0.0000</div><div> atom: 3 charge: 5.8687 magn: -0.0001 constr: 0.0000</div><div> atom: 4 charge: 5.8645 magn: 0.0000 constr: 0.0000</div><div> atom: 5 charge: 5.8650 magn: -0.0001 constr: 0.0000</div><div> atom: 6 charge: 5.8641 magn: -0.0001 constr: 0.0000</div><div> atom: 7 charge: 5.8686 magn: 0.0000 constr: 0.0000</div><div> atom: 8 charge: 5.8678 magn: 0.0000 constr: 0.0000</div><div> atom: 9 charge: 6.8067 magn: 0.0001 constr: 0.0000</div><div> atom: 10 charge: 6.8078 magn: 0.0002 constr: 0.0000</div><div> atom: 11 charge: 6.8067 magn: 0.0001 constr: 0.0000</div><div> atom: 12 charge: 6.8078 magn: 0.0002 constr: 0.0000</div><div> atom: 13 charge: 6.8047 magn: 0.0001 constr: 0.0000</div><div> atom: 14 charge: 6.8046 magn: 0.0001 constr: 0.0000</div><div> atom: 15 charge: 6.8047 magn: 0.0001 constr: 0.0000</div><div> atom: 16 charge: 6.8046 magn: 0.0001 constr: 0.0000</div><div> atom: 17 charge: 2.0873 magn: 0.0190 constr: 0.0000</div><div> atom: 18 charge: 2.0873 magn: 0.0190 constr: 0.0000</div><div> atom: 19 charge: 2.0887 magn: 0.0158 constr: 0.0000</div><div> atom: 20 charge: 2.0994 magn: 0.0157 constr: 0.0000</div><div> atom: 21 charge: 2.0880 magn: 0.0148 constr: 0.0000</div><div> atom: 22 charge: 2.0880 magn: 0.0148 constr: 0.0000</div><div> atom: 23 charge: 2.0910 magn: 0.0098 constr: 0.0000</div><div> atom: 24 charge: 2.0910 magn: 0.0098 constr: 0.0000</div><div> atom: 25 charge: 2.0925 magn: 0.0008 constr: 0.0000</div><div> atom: 26 charge: 2.0898 magn: 0.0012 constr: 0.0000</div><div> atom: 27 charge: 2.0911 magn: 0.0018 constr: 0.0000</div><div> atom: 28 charge: 2.0898 magn: 0.0012 constr: 0.0000</div><div> atom: 29 charge: 2.0911 magn: 0.0018 constr: 0.0000</div><div> atom: 30 charge: 2.0929 magn: 0.0015 constr: 0.0000</div><div> atom: 31 charge: 2.0929 magn: 0.0013 constr: 0.0000</div><div> atom: 32 charge: 2.0917 magn: 0.0010 constr: 0.0000</div><div> atom: 33 charge: 2.0912 magn: 0.0006 constr: 0.0000</div><div> atom: 34 charge: 2.0891 magn: 0.0020 constr: 0.0000</div><div> atom: 35 charge: 2.0917 magn: 0.0010 constr: 0.0000</div><div> atom: 36 charge: 2.0912 magn: 0.0006 constr: 0.0000</div><div> atom: 37 charge: 2.0891 magn: 0.0020 constr: 0.0000</div><div> atom: 38 charge: 2.0902 magn: 0.0000 constr: 0.0000</div><div> atom: 39 charge: 2.0941 magn: 0.0009 constr: 0.0000</div><div> atom: 40 charge: 2.0937 magn: 0.0010 constr: 0.0000</div></div><div><br></div><div><br clear="all"><div>...........</div><div><br></div><div><div><br></div><div> -72.3758 -62.5598 -62.5537 -62.5511 -62.5480 -62.5425 -62.5369 -62.5357</div><div> -39.4541 -38.8156 -38.6976 -30.0229 -30.0158 -30.0130 -30.0103 -30.0084</div><div> -30.0040 -30.0035 -30.0009 -29.9992 -29.9984 -29.9982 -29.9924 -29.9876</div><div> -29.9864 -29.9812 -29.9751 -29.9728 -29.9689 -29.9646 -29.9590 -29.9563</div><div> -14.3294 -14.3240 -14.3033 -14.3024 -14.2913 -14.2772 -14.2751 -14.2435</div><div> -14.2228 -14.2223 -14.2170 -14.1807 -14.1794 -14.1656 -14.1642 -14.1571</div><div> -14.0171 -14.0153 -14.0145 -14.0131 -14.0086 -14.0081 -14.0069 -14.0058</div><div> -14.0043 -14.0027 -14.0024 -14.0008 -13.9995 -13.9982 -13.9976 -13.9961</div><div> -13.9930 -13.9927 -13.9912 -13.9910 -13.9868 -13.9859 -13.9824 -13.9792</div><div> -9.7843 -9.0442 -9.0383 -9.0323 -8.3879 -8.3807 -8.3738 -7.9140</div><div> -7.4971 -7.4713 -7.4614 -7.4010 -7.3331 -7.3190 -7.2015 -7.1094</div><div> -6.9305 -6.9052 -6.8962 -6.4064 -6.3960 -6.3038 -6.2927 -6.1830</div><div> -0.3517 0.8496 1.1505 1.1568 1.9540 1.9947 2.0011 2.2299</div><div> 3.5932 3.7287 3.9068 3.9226 4.2995 4.3190 4.4402 4.4506</div><div> 4.7519 4.8447 4.8627 4.9689 5.1655 5.2280 5.2371 5.6173</div><div> 5.6668 5.7574 5.7815 5.7941 5.8425 5.9005 5.9230 6.6227</div><div> 6.6506 6.6777 6.7071 6.7640 6.8696 6.8963 6.9263 7.0988</div><div> 7.4069 7.4359 7.4917 7.6701 7.7331 7.7787 7.8363 7.9296</div><div> 7.9556 7.9730 8.0333 8.1732 8.2016 8.2275 8.3257 8.3813</div><div> 8.4031 8.4117 8.4690 8.5685 8.6023 8.6228 8.6804 8.7394</div><div> 8.7624 8.7812 8.7938 8.8262 8.9081 8.9152 8.9798 9.0363</div><div> 11.4012 13.3197 13.7582 13.8042 14.1482 14.5752 14.6947 14.7607</div><div> 14.9438 15.5053 15.7670 15.8208 16.5742 17.0272 18.0198 18.1416</div><div> 18.1597 18.4134 18.4539 18.4813 18.5170 18.7278 18.8137 18.8752</div><div> 19.0233 19.2133 19.4970 19.5531 19.5742 19.8056 20.2091 20.2829</div><div> 20.3267 20.4359 20.7334 20.8263 20.8985</div><div><br></div><div> the Fermi energy is 11.9044 ev</div><div><br></div><div>! total energy = -3562.13880898 Ry</div><div> Harris-Foulkes estimate = -3562.13880897 Ry</div><div> estimated scf accuracy < 1.5E-09 Ry</div><div><br></div><div> The total energy is the sum of the following terms:</div><div><br></div><div> one-electron contribution = -1662.94477982 Ry</div><div> hartree contribution = 1067.83968469 Ry</div><div> xc contribution = -752.14181346 Ry</div><div> ewald contribution = -2214.94330371 Ry</div><div> Hubbard energy = 0.05140332 Ry</div><div> smearing contrib. (-TS) = 0.00000000 Ry</div><div><br></div><div> total magnetization = 4.00 Bohr mag/cell</div><div> absolute magnetization = 4.03 Bohr mag/cell</div><div><br></div><div> convergence has been achieved in 19 iterations</div><div><br></div><div> Forces acting on atoms (Ry/au):</div><div><br></div><div><b><font color="#ff0000"> atom 1 type 1 force = NaN NaN NaN</font></b></div><div><b><font color="#ff0000"> atom 2 type 2 force = NaN NaN NaN</font></b></div><div><b><font color="#ff0000"> atom 3 type 2 force = NaN NaN NaN</font></b></div><div><b><font color="#ff0000"> atom 4 type 2 force = NaN NaN NaN</font></b></div><div><b><font color="#ff0000"> atom 5 type 2 force = NaN NaN NaN</font></b></div><div><b><font color="#ff0000"> atom 6 type 2 force = NaN NaN NaN</font></b></div><div><b><font color="#ff0000"> atom 7 type 2 force = NaN NaN NaN</font></b></div><div><b><font color="#ff0000"> atom 8 type 2 force = NaN NaN NaN</font></b></div><div><b><font color="#ff0000"> atom 9 type 3 force = NaN NaN NaN</font></b></div><div><b><font color="#ff0000"> atom 10 type 3 force = NaN NaN NaN</font></b></div><div><b><font color="#ff0000"> atom 11 type 3 force = NaN NaN NaN</font></b></div><div><b><font color="#ff0000"> atom 12 type 3 force = NaN NaN NaN</font></b></div><div><b><font color="#ff0000"> atom 13 type 3 force = NaN NaN NaN</font></b></div><div><b><font color="#ff0000"> atom 14 type 3 force = NaN NaN NaN</font></b></div><div><b><font color="#ff0000"> atom 15 type 3 force = NaN NaN NaN</font></b></div><div><b><font color="#ff0000"> atom 16 type 3 force = NaN NaN NaN</font></b></div><div><b><font color="#ff0000"> atom 17 type 4 force = NaN NaN NaN</font></b></div><div><b><font color="#ff0000"> atom 18 type 4 force = NaN NaN NaN</font></b></div><div><b><font color="#ff0000"> atom 19 type 4 force = NaN NaN NaN</font></b></div><div><b><font color="#ff0000"> atom 20 type 4 force = NaN NaN NaN</font></b></div><div><b><font color="#ff0000"> atom 21 type 4 force = NaN NaN NaN</font></b></div><div><b><font color="#ff0000"> atom 22 type 4 force = NaN NaN NaN</font></b></div><div><b><font color="#ff0000"> atom 23 type 4 force = NaN NaN NaN</font></b></div><div><b><font color="#ff0000"> atom 24 type 4 force = NaN NaN NaN</font></b></div><div><b><font color="#ff0000"> atom 25 type 4 force = NaN NaN NaN</font></b></div><div><b><font color="#ff0000"> atom 26 type 4 force = NaN NaN NaN</font></b></div><div><b><font color="#ff0000"> atom 27 type 4 force = NaN NaN NaN</font></b></div><div><b><font color="#ff0000"> atom 28 type 4 force = NaN NaN NaN</font></b></div><div><b><font color="#ff0000"> atom 29 type 4 force = NaN NaN NaN</font></b></div><div><b><font color="#ff0000"> atom 30 type 4 force = NaN NaN NaN</font></b></div><div><b><font color="#ff0000"> atom 31 type 4 force = NaN NaN NaN</font></b></div><div><b><font color="#ff0000"> atom 32 type 4 force = NaN NaN NaN</font></b></div><div><b><font color="#ff0000"> atom 33 type 4 force = NaN NaN NaN</font></b></div><div><b><font color="#ff0000"> atom 34 type 4 force = NaN NaN NaN</font></b></div><div><b><font color="#ff0000"> atom 35 type 4 force = NaN NaN NaN</font></b></div><div><b><font color="#ff0000"> atom 36 type 4 force = NaN NaN NaN</font></b></div><div><b><font color="#ff0000"> atom 37 type 4 force = NaN NaN NaN</font></b></div><div><b><font color="#ff0000"> atom 38 type 4 force = NaN NaN NaN</font></b></div><div><b><font color="#ff0000"> atom 39 type 4 force = NaN NaN NaN</font></b></div><div><b><font color="#ff0000"> atom 40 type 4 force = NaN NaN NaN</font></b></div><div><b><font color="#ff0000"><br></font></b></div><div><b><font color="#ff0000"> Total force = NaN Total SCF correction = 0.000083</font></b></div><div><br></div><div><br></div><div> entering subroutine stress ...</div><div><br></div><div> total stress (Ry/bohr**3) (kbar) P= -21.21</div><div> -0.00015382 0.00000344 0.00000000 -22.63 0.51 0.00</div><div> 0.00000344 -0.00015120 0.00000000 0.51 -22.24 0.00</div><div> 0.00000000 0.00000000 -0.00012761 0.00 0.00 -18.77</div><div><br></div><div> </div><div> init_run : 3.43s CPU 3.56s WALL ( 1 calls)</div><div> electrons : 149.61s CPU 150.70s WALL ( 1 calls)</div><div> forces : 4.01s CPU 4.01s WALL ( 1 calls)</div><div> stress : 18.35s CPU 18.36s WALL ( 1 calls)</div><div><br></div><div> Called by init_run:</div><div> wfcinit : 2.00s CPU 2.06s WALL ( 1 calls)</div><div> potinit : 0.28s CPU 0.28s WALL ( 1 calls)</div><div><br></div><div> Called by electrons:</div><div> c_bands : 108.27s CPU 109.25s WALL ( 19 calls)</div><div> sum_band : 28.75s CPU 28.82s WALL ( 19 calls)</div><div> v_of_rho : 0.28s CPU 0.28s WALL ( 20 calls)</div><div> newd : 12.50s CPU 12.51s WALL ( 20 calls)</div><div> mix_rho : 0.31s CPU 0.31s WALL ( 19 calls)</div><div><br></div><div> Called by c_bands:</div><div> init_us_2 : 0.46s CPU 0.58s WALL ( 104 calls)</div><div> cegterg : 99.68s CPU 100.60s WALL ( 38 calls)</div><div><br></div><div> Called by *egterg:</div><div> h_psi : 48.38s CPU 48.61s WALL ( 158 calls)</div><div> s_psi : 9.59s CPU 9.59s WALL ( 180 calls)</div><div> g_psi : 0.17s CPU 0.17s WALL ( 118 calls)</div><div> cdiaghg : 28.96s CPU 29.32s WALL ( 156 calls)</div><div><br></div><div> Called by h_psi:</div><div> add_vuspsi : 8.32s CPU 8.55s WALL ( 158 calls)</div><div> vhpsi : 2.63s CPU 2.60s WALL ( 158 calls)</div><div><br></div><div> General routines</div><div> calbec : 20.62s CPU 19.83s WALL ( 3336 calls)</div><div> fft : 0.50s CPU 0.42s WALL ( 226 calls)</div><div> fftw : 32.89s CPU 31.71s WALL ( 44122 calls)</div><div> </div><div> Parallel routines</div><div> fft_scatter : 19.21s CPU 12.04s WALL ( 44348 calls)</div><div> Hubbard U routines</div><div> new_ns : 0.27s CPU 0.30s WALL ( 19 calls)</div><div> vhpsi : 2.63s CPU 2.60s WALL ( 158 calls)</div><div> force_hub : 1.19s CPU 1.20s WALL ( 1 calls)</div><div> stres_hub : 6.66s CPU 6.68s WALL ( 1 calls)</div><div> </div><div> PWSCF : 2m56.00s CPU 2m57.40s WALL</div><div><br></div><div> </div><div> This run was terminated on: 13:31:21 8Sep2014 </div><div><br></div><div>=------------------------------------------------------------------------------=</div><div> JOB DONE.</div><div>=------------------------------------------------------------------------------=</div><div>Application 5692467 resources: utime ~5600s, stime ~86s, Rss ~108008, inblocks ~646892, outblocks ~142424</div></div><div><br></div><div><br></div><div><br></div><div><br></div>-- <br><div dir="ltr"><div><div>Mehmet Topsakal</div><div>Chemical Engineering and Materials Science,</div><div>University of Minnesota, Postdoctoral Associate,</div><div><a href="https://www.researchgate.net/profile/Mehmet_Topsakal/" style="color:rgb(17,85,204)" target="_blank">www.researchgate.net/profile/Mehmet_Topsakal</a></div></div></div>
</div></div>