<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40"><head><META HTTP-EQUIV="Content-Type" CONTENT="text/html; charset=us-ascii"><meta name=Generator content="Microsoft Word 15 (filtered medium)"><style><!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:"\B9D1\C740 \ACE0\B515";
panose-1:2 11 5 3 2 0 0 2 0 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
{font-family:"\@\B9D1\C740 \ACE0\B515";}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0in;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri",sans-serif;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:#0563C1;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:#954F72;
text-decoration:underline;}
span.EmailStyle17
{mso-style-type:personal-compose;
font-family:"Calibri",sans-serif;
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
font-family:"Calibri",sans-serif;}
@page WordSection1
{size:8.5in 11.0in;
margin:85.05pt 1.0in 1.0in 1.0in;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]--></head><body lang=EN-US link="#0563C1" vlink="#954F72"><div class=WordSection1><p class=MsoNormal>Dear Members,<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>I am trying a simple HSE calculation but it stop progressing after some time where output file says “convergence has been achieved in 10 iterations” without showing any further error. CPU usage still remains high but nothing is further printed in the output file even if I wait for hours. Here is my input file:<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>&CONTROL<o:p></o:p></p><p class=MsoNormal> calculation = "scf"<o:p></o:p></p><p class=MsoNormal> max_seconds = 8.64000e+04<o:p></o:p></p><p class=MsoNormal> outdir = "./"<o:p></o:p></p><p class=MsoNormal> prefix = "espresso"<o:p></o:p></p><p class=MsoNormal> pseudo_dir = " \.pseudopot"<o:p></o:p></p><p class=MsoNormal> title = "LiFePO4 (SCF)"<o:p></o:p></p><p class=MsoNormal> wf_collect = .TRUE.<o:p></o:p></p><p class=MsoNormal> wfcdir = "./"<o:p></o:p></p><p class=MsoNormal> tstress = .TRUE.<o:p></o:p></p><p class=MsoNormal> tprnfor = .true.<o:p></o:p></p><p class=MsoNormal> verbosity = 'high'<o:p></o:p></p><p class=MsoNormal>/<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>&SYSTEM<o:p></o:p></p><p class=MsoNormal> a = 4.75196e+00<o:p></o:p></p><p class=MsoNormal> b = 6.08735e+00<o:p></o:p></p><p class=MsoNormal> c = 1.04507e+01<o:p></o:p></p><p class=MsoNormal> degauss = 1.00000e-02<o:p></o:p></p><p class=MsoNormal> ecutrho = 3.50000e+02<o:p></o:p></p><p class=MsoNormal> ecutwfc = 3.50000e+01<o:p></o:p></p><p class=MsoNormal> exxdiv_treatment = "gygi-baldereschi"<o:p></o:p></p><p class=MsoNormal> ibrav = 8<o:p></o:p></p><p class=MsoNormal> input_dft = "hse"<o:p></o:p></p><p class=MsoNormal> nat = 28<o:p></o:p></p><p class=MsoNormal> nbnd = 132<o:p></o:p></p><p class=MsoNormal> nqx1 = 1<o:p></o:p></p><p class=MsoNormal> nqx2 = 1<o:p></o:p></p><p class=MsoNormal> nqx3 = 1<o:p></o:p></p><p class=MsoNormal> exx_fraction = 0.25<o:p></o:p></p><p class=MsoNormal> ntyp = 4<o:p></o:p></p><p class=MsoNormal> occupations = "smearing"<o:p></o:p></p><p class=MsoNormal> smearing = "gaussian"<o:p></o:p></p><p class=MsoNormal> x_gamma_extrapolation = .TRUE.<o:p></o:p></p><p class=MsoNormal>/<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>&ELECTRONS<o:p></o:p></p><p class=MsoNormal> conv_thr = 1.00000e-06<o:p></o:p></p><p class=MsoNormal> electron_maxstep = 200<o:p></o:p></p><p class=MsoNormal> adaptive_thr = .TRUE.<o:p></o:p></p><p class=MsoNormal> mixing_beta = 4.00000e-01<o:p></o:p></p><p class=MsoNormal> startingpot = "atomic"<o:p></o:p></p><p class=MsoNormal> startingwfc = "atomic+random"<o:p></o:p></p><p class=MsoNormal>/<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>K_POINTS {gamma}<o:p></o:p></p><p class=MsoNormal> <o:p></o:p></p><p class=MsoNormal>ATOMIC_SPECIES<o:p></o:p></p><p class=MsoNormal>Li 6.94100 Li.pbe-n-van.UPF<o:p></o:p></p><p class=MsoNormal>Fe 55.84500 Fe.pbe-sp-van.UPF<o:p></o:p></p><p class=MsoNormal>P 30.97376 P.pbe-n-van.UPF<o:p></o:p></p><p class=MsoNormal>O 15.99940 O.pbe-rrkjus.UPF<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>ATOMIC_POSITIONS {angstrom}<o:p></o:p></p><p class=MsoNormal>Li 0.000000 0.000000 0.000000<o:p></o:p></p><p class=MsoNormal>Li 2.375978 3.043675 5.225334<o:p></o:p></p><p class=MsoNormal>Li 2.375978 0.000000 5.225334<o:p></o:p></p><p class=MsoNormal>Li 0.000000 3.043675 0.000000<o:p></o:p></p><p class=MsoNormal>Fe 2.263451 4.565513 8.168085<o:p></o:p></p><p class=MsoNormal>Fe 4.639429 1.521838 7.507916<o:p></o:p></p><p class=MsoNormal>Fe 0.112526 4.565513 2.942751<o:p></o:p></p><p class=MsoNormal>Fe 2.488504 1.521838 2.282583<o:p></o:p></p><p class=MsoNormal>P 2.760801 4.565513 0.989699<o:p></o:p></p><p class=MsoNormal>P 0.384823 1.521838 4.235635<o:p></o:p></p><p class=MsoNormal>P 4.367133 4.565513 6.215033<o:p></o:p></p><p class=MsoNormal>P 1.991155 1.521838 9.460968<o:p></o:p></p><p class=MsoNormal>O 1.223358 4.565513 1.010109<o:p></o:p></p><p class=MsoNormal>O 1.359820 2.756426 8.719891<o:p></o:p></p><p class=MsoNormal>O 1.359820 0.287250 8.719891<o:p></o:p></p><p class=MsoNormal>O 3.528598 1.521838 9.440558<o:p></o:p></p><p class=MsoNormal>O 1.396395 1.521838 0.446421<o:p></o:p></p><p class=MsoNormal>O 3.599335 1.521838 4.215225<o:p></o:p></p><p class=MsoNormal>O 3.392136 5.800101 1.730777<o:p></o:p></p><p class=MsoNormal>O 3.735797 3.330925 6.956111<o:p></o:p></p><p class=MsoNormal>O 1.016158 2.756426 3.494557<o:p></o:p></p><p class=MsoNormal>O 1.016158 0.287250 3.494557<o:p></o:p></p><p class=MsoNormal>O 3.355560 4.565513 10.004246<o:p></o:p></p><p class=MsoNormal>O 3.772373 4.565513 4.778913<o:p></o:p></p><p class=MsoNormal>O 0.979582 1.521838 5.671755<o:p></o:p></p><p class=MsoNormal>O 1.152620 4.565513 6.235443<o:p></o:p></p><p class=MsoNormal>O 3.735797 5.800101 6.956111<o:p></o:p></p><p class=MsoNormal>O 3.392136 3.330925 1.730777<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>And last few lines of output file are<o:p></o:p></p><p class=MsoNormal>End of self-consistent calculation<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal> k = 0.0000 0.0000 0.0000 ( 3575 PWs) bands (ev):<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal> -79.0030 -79.0023 -79.0019 -78.9839 -45.2255 -45.2252 -45.2234 -45.2059<o:p></o:p></p><p class=MsoNormal> -45.1487 -45.1485 -45.1473 -45.1309 -44.9335 -44.9334 -44.9323 -44.9147<o:p></o:p></p><p class=MsoNormal> -15.3106 -15.1682 -15.0711 -14.9777 -12.6234 -12.6211 -12.5757 -12.5494<o:p></o:p></p><p class=MsoNormal> -12.5403 -12.5249 -12.5217 -12.4742 -12.4125 -12.3868 -12.3624 -12.3553<o:p></o:p></p><p class=MsoNormal> -2.5810 -2.1773 -2.1375 -2.0024 -0.9665 -0.6145 -0.5687 -0.4843<o:p></o:p></p><p class=MsoNormal> -0.1764 -0.1015 0.1339 0.2764 0.3243 0.3686 0.4349 0.4662<o:p></o:p></p><p class=MsoNormal> 1.5656 1.9325 2.0392 2.0667 2.2272 2.2489 2.2587 2.2709<o:p></o:p></p><p class=MsoNormal> 2.3867 2.5635 2.7745 2.9206 3.0322 3.0398 3.1022 3.1224<o:p></o:p></p><p class=MsoNormal> 3.1659 3.1803 3.3798 3.4103 3.4174 3.5710 3.5829 3.5833<o:p></o:p></p><p class=MsoNormal> 3.7764 3.8840 4.0086 4.0279 4.1655 4.3439 4.7040 4.8756<o:p></o:p></p><p class=MsoNormal> 7.1039 7.1654 7.2249 7.2788 7.2957 7.3083 7.3196 7.3530<o:p></o:p></p><p class=MsoNormal> 7.4510 7.4565 7.4624 7.4743 8.0415 8.0544 8.0944 8.3315<o:p></o:p></p><p class=MsoNormal> 8.4876 8.5099 8.6253 8.6437 10.2682 12.0817 12.4471 13.1204<o:p></o:p></p><p class=MsoNormal> 13.3931 13.9195 13.9632 13.9830 14.1605 14.3201 14.5312 14.7573<o:p></o:p></p><p class=MsoNormal> 14.8969 14.9143 15.1052 15.6954 15.7497 16.0537 16.3593 16.5595<o:p></o:p></p><p class=MsoNormal> 16.5733 16.6684 16.6782 16.7533 16.9513 17.0567 17.2068 17.6553<o:p></o:p></p><p class=MsoNormal> 17.7728 17.7858 18.0596 18.5210<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal> occupation numbers<o:p></o:p></p><p class=MsoNormal> 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000<o:p></o:p></p><p class=MsoNormal> 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000<o:p></o:p></p><p class=MsoNormal> 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000<o:p></o:p></p><p class=MsoNormal> 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000<o:p></o:p></p><p class=MsoNormal> 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000<o:p></o:p></p><p class=MsoNormal> 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000<o:p></o:p></p><p class=MsoNormal> 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000<o:p></o:p></p><p class=MsoNormal> 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000<o:p></o:p></p><p class=MsoNormal> 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000<o:p></o:p></p><p class=MsoNormal> 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000<o:p></o:p></p><p class=MsoNormal> 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000<o:p></o:p></p><p class=MsoNormal> 0.9994 0.9993 0.9991 0.9987 0.0019 0.0013 0.0003 0.0000<o:p></o:p></p><p class=MsoNormal> 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000<o:p></o:p></p><p class=MsoNormal> 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000<o:p></o:p></p><p class=MsoNormal> 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000<o:p></o:p></p><p class=MsoNormal> 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000<o:p></o:p></p><p class=MsoNormal> 0.0000 0.0000 0.0000 0.0000<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal> the Fermi energy is 7.7636 ev<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>! total energy = -1586.80562557 Ry<o:p></o:p></p><p class=MsoNormal> Harris-Foulkes estimate = -1586.80581500 Ry<o:p></o:p></p><p class=MsoNormal> estimated scf accuracy < 0.00024789 Ry<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal> convergence has been achieved in 10 iterations<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>Shoaib Muhammad, PhD<o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>Department of Energy Science<o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>Sungkyunkwan University<o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>South Korea<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p></div></body></html>