<div dir="ltr"><div>Dear all,</div><div><br></div><div>
I used lambda.x to calculate Tc for MgB2, as it explained in the PHonon/.../example03.I have a question about the q-points and in the <a href="http://lambda.in">lambda.in</a>

 input file.

<br></div><div>First time I used the sequence and weights of q-points calculated by kpoints.x

and I got Tc=9 K which is not correct.</div><div>Second time I used the sequence and weights of q-points calculated

by ph.x and I got 44 for Tc which is almost consistent with literature.</div><div>Can anyone help me find out which one is basically right?</div><div>I wonder why it is needed to reproduce q points by means of kpoints.x while these q points are presented in ph.out and also elph_dir/elph.inp_lambda.* files. Is it related to the algorithm of labmda.x, or the algorithm of kpoints.x?<br></div><div><br></div><div>
Thanks and Regards

<br></div><div><br></div><div>==================================<br></div><div>Input file:</div><div>==================================</div><div>&CONTROL<br>    calculation='scf'<br>    restart_mode='from_scratch',<br>    outdir='tmp',<br>    prefix='mgb2',<br>    pseudo_dir = '/home/pseudo/',<br>    tstress =.TRUE.,<br>    tprnfor= .TRUE. ,<br>    verbosity='high' ,<br>/<br>&SYSTEM<br>    ibrav =4,<br>    celldm(1)= 5.824,<br>    celldm(3)=1.141,<br>    nat = 3,<br>    ntyp = 2,<br>    nspin = 1,<br>    ecutwfc = 120,<br>    occupations = 'smearing',<br>    degauss= 0.01,<br>    smearing= 'mp',<br>    la2f=.true.<br>/<br>&ELECTRONS<br>    conv_thr=1.D-11,<br>    mixing_beta=0.3,<br>/<br>ATOMIC_SPECIES<br>    Mg    24.305   Mg.pbe-hgh.UPF<br>    B     10.811   B.pbe-hgh.UPF<br><br>ATOMIC_POSITIONS (crystal)<br>   Mg       0.000000000   0.000000000   0.000000000<br>   B        0.333333333   0.666666666   0.500000000<br>   B        0.666666666   0.333333333   0.500000000<br><br>K_POINTS (Automatic)<br>   12 12 10  0 0 0<br>==========================================</div><div>electron-phonon coefficients for MgB2<br>&inputph<br>  tr2_ph=1.0d-14,<br>  prefix='mgb2',<br>  fildvscf='mgb2dv',<br>  amass(1)=24.305,<br>  amass(2)=10.811,<br>  outdir='tmp',<br>  fildyn='mgb2.dyn',<br>  electron_phonon='interpolated',<br>  el_ph_sigma=0.005,<br>  el_ph_nsigma=10,<br>  trans=.true.,<br>  ldisp=.true.<br>  nq1=6, nq2=6, nq3=5,<br> /<br>===========================================<br></div><div><br></div><div>
<div>M Moaddeli, PhD</div><div>Pasargad HEI, Shiraz, Iran</div>

<br></div></div>