<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:"Yu Gothic";
panose-1:2 11 4 0 0 0 0 0 0 0;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
{font-family:"\@Yu Gothic";
panose-1:2 11 4 0 0 0 0 0 0 0;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0in;
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:8.5in 11.0in;
margin:1.0in 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" style="word-wrap:break-word">
<div class="WordSection1">
<p class="MsoNormal"><b>Hello,<o:p></o:p></b></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><b>For gfortran(10.3.0)-compiled QE 6.8, I’m getting the following error for jobs that include smearing:<o:p></o:p></b></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><b>Program received signal SIGSEGV: Segmentation fault - invalid memory reference.<o:p></o:p></b></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Backtrace for this error:<o:p></o:p></p>
<p class="MsoNormal">#0 0x146941220d21 in ???<o:p></o:p></p>
<p class="MsoNormal">#1 0x14694121fef5 in ???<o:p></o:p></p>
<p class="MsoNormal">#2 0x1469414dd3bf in ???<o:p></o:p></p>
<p class="MsoNormal">#3 0x7ffd48ddd9cc in ???<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><b>Job terminates right after this line:<o:p></o:p></b></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"> iteration # 1 ecut= 35.00 Ry beta= 0.70<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><b>This does not occur with the latest Intel compiler. It also does not occur with either compiler on QE 6.7.<o:p></o:p></b></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><b>For completeness, here are a few jobs that fail:<o:p></o:p></b></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><b>1.<o:p></o:p></b></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">&CONTROL<o:p></o:p></p>
<p class="MsoNormal"> !<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"> tprnfor = .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">&SYSTEM<o:p></o:p></p>
<p class="MsoNormal"> !<o:p></o:p></p>
<p class="MsoNormal"> ibrav = 8<o:p></o:p></p>
<p class="MsoNormal"> celldm(1) = 10.6881<o:p></o:p></p>
<p class="MsoNormal"> celldm(2) = 0.866025<o:p></o:p></p>
<p class="MsoNormal"> celldm(3) = 3.95422<o:p></o:p></p>
<p class="MsoNormal"> nat = 10<o:p></o:p></p>
<p class="MsoNormal"> ntyp = 3<o:p></o:p></p>
<p class="MsoNormal"> ecutwfc = 35<o:p></o:p></p>
<p class="MsoNormal"> ecutrho = 300<o:p></o:p></p>
<p class="MsoNormal"> occupations = 'smearing'<o:p></o:p></p>
<p class="MsoNormal"> degauss = 0.03<o:p></o:p></p>
<p class="MsoNormal"> smearing = 'mv'<o:p></o:p></p>
<p class="MsoNormal"> nbnd = 80<o:p></o:p></p>
<p class="MsoNormal"> tot_charge = 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">&ELECTRONS<o:p></o:p></p>
<p class="MsoNormal"> !<o:p></o:p></p>
<p class="MsoNormal"> conv_thr = 5.D-6<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">ATOMIC_SPECIES<o:p></o:p></p>
<p class="MsoNormal">C 1 C.pbe-rrkjus.UPF<o:p></o:p></p>
<p class="MsoNormal">O 1 O.pbe-rrkjus.UPF<o:p></o:p></p>
<p class="MsoNormal">Pt 1 Pt.pbe-n-rrkjus_psl.0.1.UPF<o:p></o:p></p>
<p class="MsoNormal">ATOMIC_POSITIONS (bohr)<o:p></o:p></p>
<p class="MsoNormal">C 5.335084148 4.646723426 12.901029877<o:p></o:p></p>
<p class="MsoNormal">O 5.335009643 4.619623254 15.079854269<o:p></o:p></p>
<p class="MsoNormal">Pt 8.061327071 0.098057998 8.992142901<o:p></o:p></p>
<p class="MsoNormal">Pt 2.608989366 0.098058283 8.992140585<o:p></o:p></p>
<p class="MsoNormal">Pt 0.000036609 4.720846294 8.968756935<o:p></o:p></p>
<p class="MsoNormal">Pt 5.335159557 4.721612729 9.380196435<o:p></o:p></p>
<p class="MsoNormal">Pt 0.000041121 7.802951963 4.604626508<o:p></o:p></p>
<p class="MsoNormal">Pt 5.335161233 7.697749113 4.753489408<o:p></o:p></p>
<p class="MsoNormal">Pt 2.697860636 3.152173889 4.688412329<o:p></o:p></p>
<p class="MsoNormal">Pt 7.972463687 3.152174491 4.688415209<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><b>2.<o:p></o:p></b></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">&CONTROL<o:p></o:p></p>
<p class="MsoNormal"> !<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"> tprnfor = .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">&SYSTEM<o:p></o:p></p>
<p class="MsoNormal"> !<o:p></o:p></p>
<p class="MsoNormal"> ecutrho = 300<o:p></o:p></p>
<p class="MsoNormal"> ecutwfc = 30<o:p></o:p></p>
<p class="MsoNormal"> ibrav = 1<o:p></o:p></p>
<p class="MsoNormal"> celldm( 1 ) = 20<o:p></o:p></p>
<p class="MsoNormal"> nat = 2<o:p></o:p></p>
<p class="MsoNormal"> ntyp = 2<o:p></o:p></p>
<p class="MsoNormal"> nbnd = 6<o:p></o:p></p>
<p class="MsoNormal"> tot_charge = 0<o:p></o:p></p>
<p class="MsoNormal"> nspin = 2<o:p></o:p></p>
<p class="MsoNormal"> starting_magnetization(1) = 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"> degauss = 0.01<o:p></o:p></p>
<p class="MsoNormal"> nosym = .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"> !<o:p></o:p></p>
<p class="MsoNormal"> conv_thr = 5.D-3<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">ATOMIC_SPECIES <o:p></o:p></p>
<p class="MsoNormal">C 12 C.pbe-rrkjus.UPF <o:p></o:p></p>
<p class="MsoNormal">H 1 H.pbe-rrkjus.UPF <o:p></o:p></p>
<p class="MsoNormal">ATOMIC_POSITIONS (angstrom)<o:p></o:p></p>
<p class="MsoNormal">C 0.00000000 0.00000000 0.16007400 <o:p></o:p></p>
<p class="MsoNormal">H 0.00000000 0.00000000 -0.96044600 <o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><b>3.<o:p></o:p></b></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"> tprnfor = .true.,<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"> nspin = 1,<o:p></o:p></p>
<p class="MsoNormal"> ecutwfc = 35,<o:p></o:p></p>
<p class="MsoNormal"> nat = 2,<o:p></o:p></p>
<p class="MsoNormal"> occupations = 'smearing',<o:p></o:p></p>
<p class="MsoNormal"> ibrav = 0,<o:p></o:p></p>
<p class="MsoNormal"> degauss = 0.01,<o:p></o:p></p>
<p class="MsoNormal"> smearing = 'marzari-vanderbilt',<o:p></o:p></p>
<p class="MsoNormal"> tot_charge = 0,<o:p></o:p></p>
<p class="MsoNormal"> ntyp = 1,<o:p></o:p></p>
<p class="MsoNormal"> ecutrho = 350,<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"> mixing_beta = 0.3<o:p></o:p></p>
<p class="MsoNormal"> conv_thr = 1e-06,<o:p></o:p></p>
<p class="MsoNormal">/<o:p></o:p></p>
<p class="MsoNormal">&IONS<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">Ag 107.868000 Ag.pbe-n-rrkjus_psl.1.0.0.UPF<o:p></o:p></p>
<p class="MsoNormal">ATOMIC_POSITIONS {bohr}<o:p></o:p></p>
<p class="MsoNormal">Ag 2.771763 2.771763 18.897260 <o:p></o:p></p>
<p class="MsoNormal"> Ag 0.000000 0.000000 22.740340 <o:p></o:p></p>
<p class="MsoNormal">K_POINTS {automatic}<o:p></o:p></p>
<p class="MsoNormal"> 1 1 1 0 0 0<o:p></o:p></p>
<p class="MsoNormal">CELL_PARAMETERS {bohr}<o:p></o:p></p>
<p class="MsoNormal">5.543526 0.000000 0.000000<o:p></o:p></p>
<p class="MsoNormal">0.000000 5.543526 0.000000<o:p></o:p></p>
<p class="MsoNormal">0.000000 0.000000 65.08548<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><b><o:p> </o:p></b></p>
<p class="MsoNormal"><b>Thanks in advance for any info,<o:p></o:p></b></p>
<p class="MsoNormal"><b><o:p> </o:p></b></p>
<p class="MsoNormal"><b>Dr. Edan Bainglass<o:p></o:p></b></p>
<p class="MsoNormal"><b>Postdoctoral Research Assistant<o:p></o:p></b></p>
<p class="MsoNormal"><b>Environ Developer<o:p></o:p></b></p>
<p class="MsoNormal"><b>MATERIALab Group<o:p></o:p></b></p>
<p class="MsoNormal"><b>Physics Department<o:p></o:p></b></p>
<p class="MsoNormal"><b>University of North Texas<o:p></o:p></b></p>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
</body>
</html>