<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 >XXX.scf.in <<EOF</span><br><span style="color: rgb(0, 0, 255);">&control</span><br><span style="color: rgb(0, 0, 255);">     calculation = 'scf' ,</span><br><span style="color: rgb(0, 0, 255);">     prefix = 'XXX' ,</span><br><span style="color: rgb(0, 0, 255);">     outdir= '$TempDIR/' ,</span><br><span style="color: rgb(0, 0, 255);">     pseudo_dir = '$PseudoDIR/' ,</span><br><span style="color: rgb(0, 0, 255);">     etot_conv_thr = 1.0d-9 ,</span><br><span style="color: rgb(0, 0, 255);">     restart_mode = 'from_scratch' ,</span><br><span style="color: rgb(0, 0, 255);">     disk_io = 'low' ,</span><br><span style="color: rgb(0, 0, 255);">/</span><br><span style="color: rgb(0, 0, 255);">&system</span><br><span style="color: rgb(0, 0, 255);">     ibrav = 2 ,</span><br><span style="color: rgb(0, 0, 255);">     celldm(1) = *** ,</span><br><span style="color: rgb(0, 0, 255);">     nat = 2 ,</span><br><span style="color: rgb(0, 0, 255);">     ntyp = 2 ,</span><br><span style="color: rgb(0, 0, 255);">     ecutwfc = 42.5 ,</span><br><span style="color: rgb(0, 0, 255);">     ecutrho = 400 ,</span><br><span style="color: rgb(0, 0, 255);">     occupations = 'smearing' ,</span><br><span style="color: rgb(0, 0, 255);">     smearing = 'mp' ,</span><br><span style="color: rgb(0, 0, 255);">     degauss = 0.015 ,</span><br><span style="color: rgb(0, 0, 255);">/</span><br><span style="color: rgb(0, 0, 255);">&electrons</span><br><span style="color: rgb(0, 0, 255);">     conv_thr = 1.0d-9 ,</span><br><span style="color: rgb(0, 0, 255);">     startingpot = 'atomic' ,</span><br><span style="color: rgb(0, 0, 255);">     startingwfc = 'atomic+random' ,</span><br><span style="color: rgb(0, 0, 255);">     diagonalization = 'david' ,</span><br><span style="color: rgb(0, 0, 255);">     mixing_mode = 'plain' ,</span><br><span style="color: rgb(0, 0, 255);">     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);"> *******</span><br><span style="color: rgb(0, 0, 255);">ATOMIC_POSITIONS crystal</span><br><span style="color: rgb(0, 0, 255);">  *******</span><br><span style="color: rgb(0, 0, 255);">K_POINTS</span><br><span style="color: rgb(0, 0, 255);"> 10</span><br><span style="color: rgb(0, 0, 255);">  0.0000000   0.0000000   0.0000000  0.0160000</span><br><span style="color: rgb(0, 0, 255);">***********</span><br><span style="color: rgb(0, 0, 255);">  0.6000000  -0.2000000   1.0000000  0.3840000</span><br><span style="color: rgb(0, 0, 255);">EOF</span><br><span style="color: rgb(0, 0, 255);">pw.x <XXX.scf.in >XXX.scf.out</span><br><span style="color: rgb(0, 0, 255);">cat >XXX.band.in <<EOF</span><br><span style="color: rgb(0, 0, 255);">&control</span><br><span style="color: rgb(0, 0, 255);">     calculation = 'bands' ,</span><br><span style="color: rgb(0, 0, 255);">     prefix = 'XXX' ,</span><br><span style="color: rgb(0, 0, 255);">     outdir= '$TempDIR/' ,</span><br><span style="color: rgb(0, 0, 255);">     pseudo_dir = '$PseudoDIR/' ,</span><br><span style="color: rgb(0, 0, 255);">     etot_conv_thr = 1.0d-9 ,</span><br><span style="color: rgb(0, 0, 255);">     disk_io = 'low' ,</span><br><span style="color: rgb(0, 0, 255);">/</span><br><span style="color: rgb(0, 0, 255);">&system</span><br><span style="color: rgb(0, 0, 255);">     ibrav = 2 ,</span><br><span style="color: rgb(0, 0, 255);">     celldm(1) = *** ,</span><br><span style="color: rgb(0, 0, 255);">     nat = 2 ,</span><br><span style="color: rgb(0, 0, 255);">     ntyp = 2 ,</span><br><span style="color: rgb(0, 0, 255);">     ecutwfc = 42.5 ,</span><br><span style="color: rgb(0, 0, 255);">     ecutrho = 400 ,</span><br><span style="color: rgb(0, 0, 255);">     nbnd = 12 ,</span><br><span style="color: rgb(0, 0, 255);">/</span><br><span style="color: rgb(0, 0, 255);">&electrons</span><br><span style="color: rgb(0, 0, 255);">     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);"> ******</span><br><span style="color: rgb(0, 0, 255);">ATOMIC_POSITIONS crystal</span><br><span style="color: rgb(0, 0, 255);">  ******</span><br><span style="color: rgb(0, 0, 255);">K_POINTS</span><br><span style="color: rgb(0, 0, 255);"> 28</span><br><span style="color: rgb(0, 0, 255);">   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);">   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 <XXX.band.in >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);">        Each <psi_i|beta_j> matrix      0.01 Mb     (     52,  12)</span><br><span style="color: rgb(255, 0, 0);"> %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%</span><br><span style="color: rgb(255, 0, 0);">     from davcio : error #        99</span><br><span style="color: rgb(255, 0, 0);">     error while reading from file</span><br><span style="color: rgb(255, 0, 0);"> %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%</span><br><span style="color: rgb(255, 0, 0);">     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>