<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=ks_c_5601-1987">
<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:"Malgun Gothic";
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:"\@Malgun Gothic";
panose-1:2 11 5 3 2 0 0 2 0 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
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;}
p.MsoListParagraph, li.MsoListParagraph, div.MsoListParagraph
{mso-style-priority:34;
margin-top:0cm;
margin-right:0cm;
margin-bottom:0cm;
margin-left:36.0pt;
margin-bottom:.0001pt;
font-size:11.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-family:"Calibri",sans-serif;}
@page WordSection1
{size:612.0pt 792.0pt;
margin:72.0pt 72.0pt 72.0pt 72.0pt;}
div.WordSection1
{page:WordSection1;}
/* List Definitions */
@list l0
{mso-list-id:1394310577;
mso-list-type:hybrid;
mso-list-template-ids:-852177868 67698705 67698713 67698715 67698703 67698713 67698715 67698703 67698713 67698715;}
@list l0:level1
{mso-level-text:"%1\)";
mso-level-tab-stop:none;
mso-level-number-position:left;
text-indent:-18.0pt;}
@list l0:level2
{mso-level-number-format:alpha-lower;
mso-level-tab-stop:none;
mso-level-number-position:left;
text-indent:-18.0pt;}
@list l0:level3
{mso-level-number-format:roman-lower;
mso-level-tab-stop:none;
mso-level-number-position:right;
text-indent:-9.0pt;}
@list l0:level4
{mso-level-tab-stop:none;
mso-level-number-position:left;
text-indent:-18.0pt;}
@list l0:level5
{mso-level-number-format:alpha-lower;
mso-level-tab-stop:none;
mso-level-number-position:left;
text-indent:-18.0pt;}
@list l0:level6
{mso-level-number-format:roman-lower;
mso-level-tab-stop:none;
mso-level-number-position:right;
text-indent:-9.0pt;}
@list l0:level7
{mso-level-tab-stop:none;
mso-level-number-position:left;
text-indent:-18.0pt;}
@list l0:level8
{mso-level-number-format:alpha-lower;
mso-level-tab-stop:none;
mso-level-number-position:left;
text-indent:-18.0pt;}
@list l0:level9
{mso-level-number-format:roman-lower;
mso-level-tab-stop:none;
mso-level-number-position:right;
text-indent:-9.0pt;}
ol
{margin-bottom:0cm;}
ul
{margin-bottom:0cm;}
--></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 trying to calculate the band structure of a (tiny) simple cubic unit cell containing 12 atoms of 5 species. I imported to structure from a CIF file via VESTA.<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoListParagraph" style="text-indent:-18.0pt;mso-list:l0 level1 lfo1"><![if !supportLists]><span style="mso-list:Ignore">1)<span style="font:7.0pt "Times New Roman"">
</span></span><![endif]>First I am running an &scf and it seems to converge.<o:p></o:p></p>
<p class="MsoListParagraph" style="text-indent:-18.0pt;mso-list:l0 level1 lfo1"><![if !supportLists]><span style="mso-list:Ignore">2)<span style="font:7.0pt "Times New Roman"">
</span></span><![endif]>Then I am running &bands and it finishes fine without errors<o:p></o:p></p>
<p class="MsoListParagraph" style="text-indent:-18.0pt;mso-list:l0 level1 lfo1"><![if !supportLists]><span style="mso-list:Ignore">3)<span style="font:7.0pt "Times New Roman"">
</span></span><![endif]>Then I am running bands.x and generate an output file which seems fine, no NaN there<o:p></o:p></p>
<p class="MsoListParagraph" style="text-indent:-18.0pt;mso-list:l0 level1 lfo1"><![if !supportLists]><span style="mso-list:Ignore">4)<span style="font:7.0pt "Times New Roman"">
</span></span><![endif]>When I use plotband.x (interactive mode) the created ps file contains NaN and I don¡¯t understand where they come from<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">The output of the bands looks something like this:<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">&plot nbnd= 24, nks= 4 /<o:p></o:p></p>
<p class="MsoNormal"> 0.000000 0.000000 0.000000<o:p></o:p></p>
<p class="MsoNormal">-18.411 -12.728 -12.287 -11.961 -10.967 -7.359 -7.181 -5.332 -4.940 -3.029<o:p></o:p></p>
<p class="MsoNormal"> -3.013 -2.097 -2.045 -1.873 -0.112 -0.047 0.005 0.136 0.448 0.470<o:p></o:p></p>
<p class="MsoNormal"> 5.704 7.022 7.067 7.347<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">But then the ps file contains lots of:<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"> NaN 0.000 NaN 284.530 riga<o:p></o:p></p>
<p class="MsoNormal"> NaN 11.303 dot<o:p></o:p></p>
<p class="MsoNormal"> NaN 51.728 dot<o:p></o:p></p>
<p class="MsoNormal"> NaN 54.865 dot<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">My input files are as follows (maybe someone can spot the error). Version is 5.0.2 on ubuntu<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Any help is greatly appreciated!<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Yours,<o:p></o:p></p>
<p class="MsoNormal">Chris <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"> restart_mode='from_scratch',<o:p></o:p></p>
<p class="MsoNormal"> pseudo_dir = '/usr/share/espresso/pseudo/',<o:p></o:p></p>
<p class="MsoNormal"> outdir='/usr/share/doc/quantum-espresso/examples/MAPbBr3/tmp/'<o:p></o:p></p>
<p class="MsoNormal"> prefix='mapbbr3'<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"> ibrav=1, a=5.9200000763, nat=12, ntyp=5, ecutwfc = 24.0, nbnd = 24<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">/<o:p></o:p></p>
<p class="MsoNormal">ATOMIC_SPECIES<o:p></o:p></p>
<p class="MsoNormal"> C 12.01 C.pbe-hgh.UPF <o:p></o:p></p>
<p class="MsoNormal"> N 14.01 N.pbe-hgh.UPF<o:p></o:p></p>
<p class="MsoNormal"> H 1 H.pbe-hgh.UPF<o:p></o:p></p>
<p class="MsoNormal"> Pb 207.19 Pb.pbe-hgh.UPF<o:p></o:p></p>
<p class="MsoNormal"> Br 79.9 Br.pbe-hgh.UPF<o:p></o:p></p>
<p class="MsoNormal">ATOMIC_POSITIONS<o:p></o:p></p>
<p class="MsoNormal"> C 0.896270000 0.999980000 0.988510000<o:p></o:p></p>
<p class="MsoNormal"> N 0.142130000 0.999830000 0.032690000<o:p></o:p></p>
<p class="MsoNormal"> H 0.866960000 0.999760000 0.805200000<o:p></o:p></p>
<p class="MsoNormal"> H 0.822470000 0.152470000 0.064660000<o:p></o:p></p>
<p class="MsoNormal"> H 0.822200000 0.847840000 0.065070000<o:p></o:p></p>
<p class="MsoNormal"> H 0.221280000 0.142420000 0.965660000<o:p></o:p></p>
<p class="MsoNormal"> H 0.221080000 0.856980000 0.965950000<o:p></o:p></p>
<p class="MsoNormal"> H 0.180460000 0.999970000 0.205370000<o:p></o:p></p>
<p class="MsoNormal"> Pb 0.475340000 0.500000000 0.477920000<o:p></o:p></p>
<p class="MsoNormal"> Br 0.428780000 0.499880000 0.972690000<o:p></o:p></p>
<p class="MsoNormal"> Br 0.433480000 0.000020000 0.512850000<o:p></o:p></p>
<p class="MsoNormal"> Br 0.967550000 0.500170000 0.439210000<o:p></o:p></p>
<p class="MsoNormal">K_POINTS (automatic)<o:p></o:p></p>
<p class="MsoNormal">10 10 10 0 0 0<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='bands'<o:p></o:p></p>
<p class="MsoNormal"> pseudo_dir = '/usr/share/espresso/pseudo/',<o:p></o:p></p>
<p class="MsoNormal"> outdir='/usr/share/doc/quantum-espresso/examples/MAPbBr3/tmp/'<o:p></o:p></p>
<p class="MsoNormal"> prefix='mapbbr3'<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"> ibrav=1, a=5.9200000763, nat=12, ntyp=5, ecutwfc = 24.0, nbnd = 24<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">/<o:p></o:p></p>
<p class="MsoNormal">ATOMIC_SPECIES<o:p></o:p></p>
<p class="MsoNormal"> C 12.01 C.pbe-hgh.UPF <o:p></o:p></p>
<p class="MsoNormal"> N 14.01 N.pbe-hgh.UPF<o:p></o:p></p>
<p class="MsoNormal"> H 1 H.pbe-hgh.UPF<o:p></o:p></p>
<p class="MsoNormal"> Pb 207.19 Pb.pbe-hgh.UPF<o:p></o:p></p>
<p class="MsoNormal"> Br 79.9 Br.pbe-hgh.UPF<o:p></o:p></p>
<p class="MsoNormal">ATOMIC_POSITIONS<o:p></o:p></p>
<p class="MsoNormal"> C 0.896270000 0.999980000 0.988510000<o:p></o:p></p>
<p class="MsoNormal"> N 0.142130000 0.999830000 0.032690000<o:p></o:p></p>
<p class="MsoNormal"> H 0.866960000 0.999760000 0.805200000<o:p></o:p></p>
<p class="MsoNormal"> H 0.822470000 0.152470000 0.064660000<o:p></o:p></p>
<p class="MsoNormal"> H 0.822200000 0.847840000 0.065070000<o:p></o:p></p>
<p class="MsoNormal"> H 0.221280000 0.142420000 0.965660000<o:p></o:p></p>
<p class="MsoNormal"> H 0.221080000 0.856980000 0.965950000<o:p></o:p></p>
<p class="MsoNormal"> H 0.180460000 0.999970000 0.205370000<o:p></o:p></p>
<p class="MsoNormal"> Pb 0.475340000 0.500000000 0.477920000<o:p></o:p></p>
<p class="MsoNormal"> Br 0.428780000 0.499880000 0.972690000<o:p></o:p></p>
<p class="MsoNormal"> Br 0.433480000 0.000020000 0.512850000<o:p></o:p></p>
<p class="MsoNormal"> Br 0.967550000 0.500170000 0.439210000<o:p></o:p></p>
<p class="MsoNormal">K_POINTS |(crystal)<o:p></o:p></p>
<p class="MsoNormal">4<o:p></o:p></p>
<p class="MsoNormal">0.0 0.0 0.0 20 !G<o:p></o:p></p>
<p class="MsoNormal">0.5 0.5 0.5 20 !R<o:p></o:p></p>
<p class="MsoNormal">0.0 0.5 0.0 20 !X<o:p></o:p></p>
<p class="MsoNormal">0.5 0.5 0.0 20 !M<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">&bands<o:p></o:p></p>
<p class="MsoNormal">outdir='/usr/share/doc/quantum-espresso/examples/MAPbBr3/tmp/'<o:p></o:p></p>
<p class="MsoNormal">prefix='mapbbr3'<o:p></o:p></p>
<p class="MsoNormal">filband='bands.dat'<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"><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>
</div>
</body>
</html>