Hello, everyone,<br><br>I am trying splitting q points for electron-phonon coupling calculations on a computer cluster. Though jobs were finished after successful running, I was not given by file elph$qpoint . I did these calculations on example07 as follows:<br clear="all">
1) a electronic scf calculation on a very dense k mesh with la2f=.true.,<br>2) a scf calculation on a medium k mesh without la2f,<br>3) split q point using following INPUT.<br><br>Electron-phonon coefficients for Al<br> &inputph<br>
  tr2_ph=1.0d-12,<br>  prefix='al',<br>  amass(1)=26.98,<br>  outdir='$OUTDIR',<br>  fildyn='al.dyn',<br>  ldisp=.true.,<br>  nq1=4, nq2=4, nq3=4<br>  start_irr=0,<br>  last_irr=0,<br>/<br><br>4) copy file $prefix.a2Fsave, $prefix.wfc*, $prefix.dyn0, and _ph0$prefix.phsave(including both data-file.xml and data-file.xml.$qpoint) to corresponding sub-job. Then a ph.x run for electron-phonon coupling using following INPUT.<br>
<br>phonons of al<br> &inputph<br>  tr2_ph=1.0d-12,<br>  nq1=4, nq2=4, nq3=4<br>  amass(1)=26.98,<br>  prefix='$prefix'<br>  start_q=$i,<br>  last_q=$i,<br>  recover=.true.,<br>  outdir='$INDIR/$i',<br>
  fildyn='$prefix.dyn',<br>  fildvscf='$prefix.dv',<br>  elph=.true.,<br>  trans=.true.,<br>  ldisp=.true.,<br>/<br><br>5) a run of ph.x by setting trans=.false.<br><br>Though the job runs successfully, it does not give me the file elph$qpoint. Having browsed the source code, the PH should produce files elph$qpoint when elph=.true.<br>
Any suggestion would be appreciated.<br><br>Thanks a lot.<br><br>Mithra Chan<br><br>-- <br>Dr. Mithra Chan,<br>Department of Physics, National University of Singapore, Singapore<br><a href="mailto:chan.mithra@gmail.com">chan.mithra@gmail.com</a><br>
<br>