<html 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=utf-8">
<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:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0in;
margin-bottom:.0001pt;
font-size:12.0pt;
font-family:"Calibri",sans-serif;}
span.EmailStyle17
{mso-style-type:personal-compose;
font-family:"Calibri",sans-serif;
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
font-size:12.0pt;
font-family:"Calibri",sans-serif;}
@page WordSection1
{size:8.5in 11.0in;
margin:1.0in 1.0in 1.0in 1.0in;}
div.WordSection1
{page:WordSection1;}
--></style>
</head>
<body lang="EN-US" link="#0563C1" vlink="#954F72">
<div class="WordSection1">
<p class="MsoNormal"><span style="font-size:11.0pt">Hello,<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt">I am attempting to calculate the q=0 phonon modes for a system comprising a CO molecule above a Au(111) surface. The input files for both the pw.x (slab1.pwi) and the ph.x (slab1.phi) calculations as well
as their corresponding output files (slab1.pwo and slab1.pho). These calculations were run with QE 6.5 compiled with the Intel compiler stack (ifort/MKL/IMPI), and the jobs were executed on 5 compute nodes of 36 CPU cores each. Both pw.x and ph.x were run
with the “-nk 5” argument to enable k-point parallelization. Note that the pw.x input file was generated with ASE, and some of the settings like ecutrho and ecutwfc are much rounder numbers when represented in eV rather than Rydberg.<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt">The problem I am encountering is that the ph.x calculation does not converge even for the first representation/mode. In fact, the calculation *<b>diverges</b>*, with |ddv_scf|^2 = 2.464E+18 at the 100<sup>th</sup>
iteration. I am wondering whether I have done something wrong, or if there are some keywords that should be tweaked to improve convergence. I have seen advice on this mailing list to lower the value of alpha_mix(1) from its default value of 0.7, but if I understand
the manual correctly this would result in using more new density, which I imagine would make convergence even worse in my case. I haven’t been able to find any advice online for improving convergence of metallic systems.<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt">Regards,<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt">-- <o:p></o:p></span></p>
<div>
<p class="MsoNormal"><span style="font-size:11.0pt">Eric Hermes<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt">Postdoctoral Researcher<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt">Sandia National Laboratories</span><o:p></o:p></p>
</div>
</div>
</body>
</html>