<div style="line-height:1.7;color:#000000;font-size:14px;font-family:arial"><div style="line-height: 1.7; color: rgb(0, 0, 0); font-size: 14px; font-family: arial;">Dear QE Developer and Users:<br>I met a problem while calculating band structure. My in put file is like this:<br><i></i><span style="color: rgb(0, 255, 255);"><i><span style="color: rgb(0, 255, 255);"></span></i></span><span style="color: rgb(0, 0, 255);"><span style="color: rgb(0, 255, 255);"><i><span style="color: rgb(0, 0, 255);">cat &gt;XXX.scf.in &lt;&lt;EOF</span><br><span style="color: rgb(0, 0, 255);">&amp;control</span><br><span style="color: rgb(0, 0, 255);">&nbsp;&nbsp;&nbsp;&nbsp; calculation = 'scf' ,</span><br><span style="color: rgb(0, 0, 255);">&nbsp;&nbsp;&nbsp;&nbsp; prefix = 'XXX' ,</span><br><span style="color: rgb(0, 0, 255);">&nbsp;&nbsp;&nbsp;&nbsp; outdir= '$TempDIR/' ,</span><br><span style="color: rgb(0, 0, 255);">&nbsp;&nbsp;&nbsp;&nbsp; pseudo_dir = '$PseudoDIR/' ,</span><br><span style="color: rgb(0, 0, 255);">&nbsp;&nbsp;&nbsp;&nbsp; etot_conv_thr = 1.0d-9 ,</span><br><span style="color: rgb(0, 0, 255);">&nbsp;&nbsp;&nbsp;&nbsp; restart_mode = 'from_scratch' ,</span><br><span style="color: rgb(0, 0, 255);">&nbsp;&nbsp;&nbsp;&nbsp; disk_io = 'low' ,</span><br><span style="color: rgb(0, 0, 255);">/</span><br><span style="color: rgb(0, 0, 255);">&amp;system</span><br><span style="color: rgb(0, 0, 255);">&nbsp;&nbsp;&nbsp;&nbsp; ibrav = 2 ,</span><br><span style="color: rgb(0, 0, 255);">&nbsp;&nbsp;&nbsp;&nbsp; celldm(1) = *** ,</span><br><span style="color: rgb(0, 0, 255);">&nbsp;&nbsp;&nbsp;&nbsp; nat = 2 ,</span><br><span style="color: rgb(0, 0, 255);">&nbsp;&nbsp;&nbsp;&nbsp; ntyp = 2 ,</span><br><span style="color: rgb(0, 0, 255);">&nbsp;&nbsp;&nbsp;&nbsp; ecutwfc = 42.5 ,</span><br><span style="color: rgb(0, 0, 255);">&nbsp;&nbsp;&nbsp;&nbsp; ecutrho = 400 ,</span><br><span style="color: rgb(0, 0, 255);">&nbsp;&nbsp;&nbsp;&nbsp; occupations = 'smearing' ,</span><br><span style="color: rgb(0, 0, 255);">&nbsp;&nbsp;&nbsp;&nbsp; smearing = 'mp' ,</span><br><span style="color: rgb(0, 0, 255);">&nbsp;&nbsp;&nbsp;&nbsp; degauss = 0.015 ,</span><br><span style="color: rgb(0, 0, 255);">/</span><br><span style="color: rgb(0, 0, 255);">&amp;electrons</span><br><span style="color: rgb(0, 0, 255);">&nbsp;&nbsp;&nbsp;&nbsp; conv_thr = 1.0d-9 ,</span><br><span style="color: rgb(0, 0, 255);">&nbsp;&nbsp;&nbsp;&nbsp; startingpot = 'atomic' ,</span><br><span style="color: rgb(0, 0, 255);">&nbsp;&nbsp;&nbsp;&nbsp; startingwfc = 'atomic+random' ,</span><br><span style="color: rgb(0, 0, 255);">&nbsp;&nbsp;&nbsp;&nbsp; diagonalization = 'david' ,</span><br><span style="color: rgb(0, 0, 255);">&nbsp;&nbsp;&nbsp;&nbsp; mixing_mode = 'plain' ,</span><br><span style="color: rgb(0, 0, 255);">&nbsp;&nbsp;&nbsp;&nbsp; mixing_beta = 0.7 ,</span><br><span style="color: rgb(0, 0, 255);">/</span><br><span style="color: rgb(0, 0, 255);">ATOMIC_SPECIES</span><br><span style="color: rgb(0, 0, 255);">&nbsp;*******</span><br><span style="color: rgb(0, 0, 255);">ATOMIC_POSITIONS crystal</span><br><span style="color: rgb(0, 0, 255);">&nbsp; *******</span><br><span style="color: rgb(0, 0, 255);">K_POINTS</span><br><span style="color: rgb(0, 0, 255);">&nbsp;10</span><br><span style="color: rgb(0, 0, 255);">&nbsp; 0.0000000&nbsp;&nbsp; 0.0000000&nbsp;&nbsp; 0.0000000&nbsp; 0.0160000</span><br><span style="color: rgb(0, 0, 255);">***********</span><br><span style="color: rgb(0, 0, 255);">&nbsp; 0.6000000&nbsp; -0.2000000&nbsp;&nbsp; 1.0000000&nbsp; 0.3840000</span><br><span style="color: rgb(0, 0, 255);">EOF</span><br><span style="color: rgb(0, 0, 255);">pw.x &lt;XXX.scf.in &gt;XXX.scf.out</span><br><span style="color: rgb(0, 0, 255);">cat &gt;XXX.band.in &lt;&lt;EOF</span><br><span style="color: rgb(0, 0, 255);">&amp;control</span><br><span style="color: rgb(0, 0, 255);">&nbsp;&nbsp;&nbsp;&nbsp; calculation = 'bands' ,</span><br><span style="color: rgb(0, 0, 255);">&nbsp;&nbsp;&nbsp;&nbsp; prefix = 'XXX' ,</span><br><span style="color: rgb(0, 0, 255);">&nbsp;&nbsp;&nbsp;&nbsp; outdir= '$TempDIR/' ,</span><br><span style="color: rgb(0, 0, 255);">&nbsp;&nbsp;&nbsp;&nbsp; pseudo_dir = '$PseudoDIR/' ,</span><br><span style="color: rgb(0, 0, 255);">&nbsp;&nbsp;&nbsp;&nbsp; etot_conv_thr = 1.0d-9 ,</span><br><span style="color: rgb(0, 0, 255);">&nbsp;&nbsp;&nbsp;&nbsp; disk_io = 'low' ,</span><br><span style="color: rgb(0, 0, 255);">/</span><br><span style="color: rgb(0, 0, 255);">&amp;system</span><br><span style="color: rgb(0, 0, 255);">&nbsp;&nbsp;&nbsp;&nbsp; ibrav = 2 ,</span><br><span style="color: rgb(0, 0, 255);">&nbsp;&nbsp;&nbsp;&nbsp; celldm(1) = *** ,</span><br><span style="color: rgb(0, 0, 255);">&nbsp;&nbsp;&nbsp;&nbsp; nat = 2 ,</span><br><span style="color: rgb(0, 0, 255);">&nbsp;&nbsp;&nbsp;&nbsp; ntyp = 2 ,</span><br><span style="color: rgb(0, 0, 255);">&nbsp;&nbsp;&nbsp;&nbsp; ecutwfc = 42.5 ,</span><br><span style="color: rgb(0, 0, 255);">&nbsp;&nbsp;&nbsp;&nbsp; ecutrho = 400 ,</span><br><span style="color: rgb(0, 0, 255);">&nbsp;&nbsp;&nbsp;&nbsp; nbnd = 12 ,</span><br><span style="color: rgb(0, 0, 255);">/</span><br><span style="color: rgb(0, 0, 255);">&amp;electrons</span><br><span style="color: rgb(0, 0, 255);">&nbsp;&nbsp;&nbsp;&nbsp; diagonalization = 'david' ,</span><br><span style="color: rgb(0, 0, 255);">/</span><br><span style="color: rgb(0, 0, 255);">ATOMIC_SPECIES</span><br><span style="color: rgb(0, 0, 255);">&nbsp;******</span><br><span style="color: rgb(0, 0, 255);">ATOMIC_POSITIONS crystal</span><br><span style="color: rgb(0, 0, 255);">&nbsp; ******</span><br><span style="color: rgb(0, 0, 255);">K_POINTS</span><br><span style="color: rgb(0, 0, 255);">&nbsp;28</span><br><span style="color: rgb(0, 0, 255);">&nbsp;&nbsp; 0.0 0.0 0.0 1.0</span><br><span style="color: rgb(0, 0, 255);">*******</span><br><span style="color: rgb(0, 0, 255);">&nbsp;&nbsp; 0.5 0.5 0.5 1.0</span><br><span style="color: rgb(0, 0, 255);">EOF</span><br><span style="color: rgb(0, 0, 255);">pw.x &lt;XXX.band.in &gt;XXX.band.out</span></i></span></span><br>Then, I got an error from <i>pw.x</i> and process stopped. The error was:<br><i></i><span style="color: rgb(255, 0, 0);"><i><span style="color: rgb(255, 0, 0);">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Each &lt;psi_i|beta_j&gt; matrix&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.01 Mb&nbsp;&nbsp;&nbsp;&nbsp; (&nbsp;&nbsp;&nbsp;&nbsp; 52,&nbsp; 12)</span><br><span style="color: rgb(255, 0, 0);">&nbsp;%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%</span><br><span style="color: rgb(255, 0, 0);">&nbsp;&nbsp;&nbsp;&nbsp; from davcio : error #&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 99</span><br><span style="color: rgb(255, 0, 0);">&nbsp;&nbsp;&nbsp;&nbsp; error while reading from file</span><br><span style="color: rgb(255, 0, 0);">&nbsp;%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%</span><br><span style="color: rgb(255, 0, 0);">&nbsp;&nbsp;&nbsp;&nbsp; stopping ...</span></i></span><br>However, after I removed <i>disk_io = 'low'</i> in SCF calculation, non-SCF calculation can run without any problem.<br>I checked <i>XXX.save/charge-density.dat</i>, the files' size seemed same when disk_io equal to 'low' and 'default'. Is this phenomenon means <i>disk_io = 'low'</i> cannot be used during <i>calculation='bands'</i> and<i> 'nscf'</i>? I am afraid <i>disk_io = 'default'</i> would lower speed when model was huge~<br><div>--<br>GAO Zhe<br>CMC Lab, MSE, SNU, Seoul, S.Korea<br>        
</div></div></div><br><br><span title="neteasefooter"><span id="netease_mail_footer"></span></span>