<div>Dear users,</div><div><br></div><div>I&#39;m computing the anharmonic scattering time for bulk Silicon. While pw and ph work as expected, d3 crashes. It seems there is some missing parameter. Is &quot;fild0rho&quot; necessary? Below are the scripts. </div>
<div><br></div><div>Any help is greatly appreciated. </div><div><br></div><div>Thanks,</div><div><br></div><div>Giuseppe</div><div><br></div><div><b><a href="http://si.pw.in">si.pw.in</a>:</b></div><div><br></div><div> &amp;system</div>
<div>    ibrav=2, celldm(1)=10.20, nat=2, ntyp=1,</div><div>    ecutwfc =16.0</div><div> /</div><div> &amp;electrons</div><div>    conv_thr =  1.0d-8</div><div>    mixing_beta = 0.7</div><div> /</div><div>ATOMIC_SPECIES</div>
<div> Si  28.0855 Si.vbc.UPF</div><div>ATOMIC_POSITIONS</div><div> Si 0.00 0.00 0.00</div><div> Si 0.25 0.25 0.25</div><div>K_POINTS</div><div>  10</div><div>   0.1250000  0.1250000  0.1250000   1.00</div><div>   0.1250000  0.1250000  0.3750000   3.00</div>
<div>   0.1250000  0.1250000  0.6250000   3.00</div><div>   0.1250000  0.1250000  0.8750000   3.00</div><div>   0.1250000  0.3750000  0.3750000   3.00</div><div>   0.1250000  0.3750000  0.6250000   6.00</div><div>   0.1250000  0.3750000  0.8750000   6.00</div>
<div>   0.1250000  0.6250000  0.6250000   3.00</div><div>   0.3750000  0.3750000  0.3750000   1.00</div><div>   0.3750000  0.3750000  0.6250000   3.00</div><div><br></div><div><b><a href="http://si.ph.in">si.ph.in</a>:</b></div>
<div><br></div><div><div>phonons of Si</div><div> &amp;inputph</div><div>  tr2_ph=1.0d-14,</div><div>  ldisp=.true.,</div><div>  nq1=4, nq2=4, nq3=4</div><div>  amass(1)=28.0855,</div><div>  amass(2)=28.0855,</div><div>  prefix=&#39;PH&#39;,</div>
<div>  outdir=&#39;tmp/&#39;,</div><div>  fildyn=&#39;si.dyn&#39;,</div><div>  fildrho = &#39;si.drho&#39;,</div><div> /</div></div><div><br></div><div><b><a href="http://si.d3.in">si.d3.in</a>:</b></div><div><br></div><div>
<div>anharm at the X-point</div><div>&amp;inputph</div><div>   prefix = &#39;PH&#39;,</div><div>   fildrho = &#39;si.din&#39;,</div><div>   amass(1) = 28.0855,</div><div>   amass(2) = 28.0855,</div><div>   outdir = &#39;tmp/&#39;,</div>
<div>   fildyn = &#39;si.anh&#39;,</div><div> /</div><div>    </div></div><div><br></div><div><br></div><div><br></div><div><br></div>