<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: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: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:70.85pt 70.85pt 56.7pt 70.85pt;}
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 All, <o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">I am a new member of this community, trying to make my way to performing calculations with QE (V.6.7).
<br>
I successfully managed to optimize a Pd metallic supercell, construct and optimize a slab and run a SCF calculation on it (Input below). As I want to check if I reached a good minimum before continuing with adsorbing molecules on the surface, I am trying to
run phonon calculations to calculate the vibrational frequencies. Nevertheless, I invariably end in a phonon calculation that reaches no convergence (Input below).
<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Could you please guide me on whether I am making mistakes with the calculations parameters? I would highly appreciate it.
<o:p></o:p></p>
<p class="MsoNormal"> <br>
Thank you very much! <br>
<br>
Natalia Levin Rojas <o:p></o:p></p>
<p class="MsoNormal">Postdoc<br>
Max Planck Institute for Chemical Energy Conversion <o:p></o:p></p>
<p class="MsoNormal">Muelheim an der Ruhr, Germany<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><b>SCF Input<o:p></o:p></b></p>
<p class="MsoNormal">&CONTROL<o:p></o:p></p>
<p class="MsoNormal"> calculation = "scf"<o:p></o:p></p>
<p class="MsoNormal"> <span lang="ES-AR">max_seconds = 1.3e+07<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="ES-AR"> pseudo_dir = "/home/levin-rojas/pseudopot"<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="DE">/<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="DE"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="DE">&SYSTEM<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="DE"> a = 5.5513e+00<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="DE"> c = 3.69057e+01<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="DE"> </span><span lang="ES-AR">degauss = 1.00000e-02<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="ES-AR"> ecutrho = 4.344e+02<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="ES-AR"> ecutwfc = 3.62136e+01<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="ES-AR"> </span>ibrav = 4<o:p></o:p></p>
<p class="MsoNormal"> nat = 16<o:p></o:p></p>
<p class="MsoNormal"> nspin = 2<o:p></o:p></p>
<p class="MsoNormal"> ntyp = 1<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"> starting_magnetization(1) = 2.00000e-01<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-10 ! Reduced to make the phonon calc more precise<o:p></o:p></p>
<p class="MsoNormal"> electron_maxstep = 200<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 {automatic}<o:p></o:p></p>
<p class="MsoNormal">3 3 1 1 1 1<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">Pd 106.42000 Pd.pbe-n-kjpaw_psl.1.0.0.UPF.txt<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><span lang="DE">ATOMIC_POSITIONS {angstrom}<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="DE">Pd 1.3878326176 0.8012653980 21.8246640127<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="DE">Pd 4.1634673824 0.8012653980 21.8246640127<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="DE">Pd -0.0000000000 3.2050450000 21.8246701760<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="DE">Pd 2.7756500000 3.2050362040 21.8246640127<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="DE">Pd 0.0000000000 1.6025200243 19.5362847798<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="DE">Pd 2.7756500000 1.6025220000 19.5362955020<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="DE">Pd -1.3878267110 4.0063069878 19.5362847798<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="DE">Pd 1.3878267110 4.0063069878 19.5362847798<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="DE">Pd 0.0000000000 0.0000000000 17.2663090000
<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="DE">Pd 2.7756500000 0.0000000000 17.2663090000
<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="DE">Pd 1.3878250000 2.4037830000 17.2663090000
<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="DE">Pd -1.3878250000 2.4037830000 17.2663090000
<o:p></o:p></span></p>
<p class="MsoNormal">Pd 1.3878250000 0.8012610000 15.0000000000
<o:p></o:p></p>
<p class="MsoNormal">Pd 4.1634750000 0.8012610000 15.0000000000
<o:p></o:p></p>
<p class="MsoNormal">Pd -0.0000000000 3.2050450000 15.0000000000
<o:p></o:p></p>
<p class="MsoNormal">Pd 2.7756500000 3.2050450000 15.0000000000
<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><b>Phonon Input<o:p></o:p></b></p>
<p class="MsoNormal">&inputph<o:p></o:p></p>
<p class="MsoNormal"> outdir = '/home/levin-rojas/Pd/012_Pd_slab_SCF_33248'<o:p></o:p></p>
<p class="MsoNormal"> <span lang="ES-AR">prefix = 'pwscf'<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="ES-AR"> tr2_ph = 1.0e-14<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="ES-AR"> </span>amass(1) = 106.42000<o:p></o:p></p>
<p class="MsoNormal"> trans = .true.<o:p></o:p></p>
<p class="MsoNormal"> ldisp = .true.<o:p></o:p></p>
<p class="MsoNormal"> nq1 = 12<o:p></o:p></p>
<p class="MsoNormal"> nq2 = 12<o:p></o:p></p>
<p class="MsoNormal"> nq3 = 12<o:p></o:p></p>
<p class="MsoNormal"> fildyn = 'pd.dyn'<o:p></o:p></p>
<p class="MsoNormal">/<o:p></o:p></p>
</div>
</body>
</html>