<div>i am calculating Tc for YC( ytrium carbide)</div><div>but it gives error NaN</div><div>I have all positive frequency in the phonon spectrum</div><div>why the error comes ? </div><div><br></div><div><br></div>i have given my input file<div>
<br></div><div><br></div><div>1. pw.x input</div><div><div>&CONTROL</div><div> calculation = 'scf' ,</div><div> restart_mode = 'from_scratch' ,</div><div> outdir = '/home/bds/yc2test/' ,</div>
<div> wfcdir = '/home/bds/yc2test/' ,</div><div> pseudo_dir = '/home/bds/bhabya/espresso-4.3.2/pseudo/' ,</div><div> prefix = 'yc1' ,</div>
<div> tstress = .true. ,</div><div> tprnfor = .true. ,</div><div> /</div><div> &SYSTEM</div><div> ibrav = 2,</div><div> celldm(1) = 9.6,</div>
<div> nat = 2,</div><div> ntyp = 2,</div><div> ecutwfc = 60 ,</div><div> ecutrho = 600 ,</div><div> occupations = 'smearing' ,</div>
<div> degauss = 0.02 ,</div><div> smearing = 'methfessel-paxton' , </div><div> la2F = .true.,</div><div>/</div><div> &ELECTRONS</div><div> /</div><div>
ATOMIC_SPECIES</div><div> y 88.90580 Y.pw91-nsp-van.UPF </div><div> c 12.01700 C.pw91-van_ak.UPF </div><div>ATOMIC_POSITIONS crystal </div><div> y 0.000000000 0.000000000 0.000000000 </div><div>
c 0.500000000 0.500000000 0.500000000 </div><div>K_POINTS automatic </div><div> 16 16 16 0 0 0 </div></div><div><br></div><div><br></div><div>2 pw.x input&CONTROL</div><div> calculation = 'scf' ,</div>
<div> restart_mode = 'from_scratch' ,</div><div> outdir = '/home/bds/yc2test/' ,</div><div> wfcdir = '/home/bds/yc2test/' ,</div><div> pseudo_dir = '/home/bds/bhabya/espresso-4.3.2/pseudo/' ,</div>
<div> prefix = 'yc1' ,</div><div> tstress = .true. ,</div><div> tprnfor = .true. ,</div><div> /</div><div> &SYSTEM</div><div> ibrav = 2,</div>
<div> celldm(1) = 9.6,</div><div> nat = 2,</div><div> ntyp = 2,</div><div> ecutwfc = 60 ,</div><div> ecutrho = 600 ,</div>
<div> occupations = 'smearing' ,</div><div> degauss = 0.02 ,</div><div> smearing = 'methfessel-paxton' ,</div><div> /</div><div> &ELECTRONS</div><div>
/</div><div>ATOMIC_SPECIES</div><div> y 88.90580 Y.pw91-nsp-van.UPF </div><div> c 12.01700 C.pw91-van_ak.UPF </div><div>ATOMIC_POSITIONS crystal </div><div> y 0.000000000 0.000000000 0.000000000 </div>
<div> c 0.500000000 0.500000000 0.500000000 </div><div>K_POINTS automatic </div><div> 8 8 8 0 0 0 </div><div><br></div><div><br></div><div>ph.x input</div><div><br></div><div><div>&INPUTPH</div><div>
outdir = '/home/bds/yc2test/' ,</div><div> prefix = 'yc1' ,</div><div> fildyn = 'yc1.dyn' ,</div><div> fildvscf = 'yc1dv' ,</div>
<div> ldisp = .true.,</div><div> nq1 = 4 ,</div><div> nq2 = 4 ,</div><div> nq3 = 4 ,</div><div> epsil = .false.,</div>
<div> elph = .true.,</div><div> fpol = .false. ,</div><div> recover = .false. ,</div><div> tr2_ph=1.0d-12,</div><div> trans=.true.,</div>
<div> amass(1) = 88.9058,</div><div> amass(2) = 12.017,</div><div> /</div></div><div><br></div><div><br></div><div>output error shows </div><div>alha2F.dat</div><div><br></div><div><div>
# E(THz) 0.005 0.010 0.015 0.020 0.025 0.030 0.035 0.040 0.045 0.050</div><div> 0.0000 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 0.0754 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 0.1508 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 0.2261 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 0.3015 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 0.3769 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 0.4523 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 0.5276 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 0.6030 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 0.6784 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 0.7538 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 0.8291 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 0.9045 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 0.9799 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 1.0553 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 1.1307 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 1.2060 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 1.2814 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 1.3568 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 1.4322 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 1.5075 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 1.5829 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 1.6583 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 1.7337 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 1.8090 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 1.8844 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 1.9598 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 2.0352 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 2.1106 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 2.1859 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 2.2613 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 2.3367 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 2.4121 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 2.4874 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 2.5628 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 2.6382 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 2.7136 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 2.7889 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 2.8643 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 2.9397 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 3.0151 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 3.0905 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 3.1658 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 3.2412 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 3.3166 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 3.3920 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 3.4673 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 3.5427 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 3.6181 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 3.6935 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 3.7688 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 3.8442 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 3.9196 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 3.9950 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 4.0704 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 4.1457 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 4.2211 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 4.2965 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 4.3719 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 4.4472 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 4.5226 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 4.5980 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 4.6734 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 4.7487 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 4.8241 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 4.8995 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 4.9749 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 5.0503 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 5.1256 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 5.2010 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 5.2764 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 5.3518 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 5.4271 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 5.5025 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 5.5779 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 5.6533 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 5.7286 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 5.8040 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 5.8794 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 5.9548 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 6.0302 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 6.1055 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 6.1809 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 6.2563 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 6.3317 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 6.4070 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 6.4824 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 6.5578 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 6.6332 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 6.7085 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 6.7839 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 6.8593 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 6.9347 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 7.0101 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 7.0854 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 7.1608 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 7.2362 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 7.3116 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 7.3869 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 7.4623 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 7.5377 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 7.6131 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 7.6884 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 7.7638 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 7.8392 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 7.9146 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 7.9899 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 8.0653 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 8.1407 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 8.2161 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 8.2915 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 8.3668 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 8.4422 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 8.5176 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 8.5930 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 8.6683 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 8.7437 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 8.8191 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 8.8945 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 8.9698 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 9.0452 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 9.1206 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 9.1960 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 9.2714 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 9.3467 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 9.4221 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 9.4975 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 9.5729 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 9.6482 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 9.7236 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 9.7990 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 9.8744 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 9.9497 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 10.0251 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 10.1005 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 10.1759 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 10.2513 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 10.3266 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 10.4020 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 10.4774 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 10.5528 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 10.6281 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 10.7035 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 10.7789 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 10.8543 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 10.9296 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 11.0050 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 11.0804 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 11.1558 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 11.2312 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 11.3065 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 11.3819 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 11.4573 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 11.5327 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 11.6080 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 11.6834 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 11.7588 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 11.8342 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 11.9095 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 11.9849 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 12.0603 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 12.1357 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 12.2111 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 12.2864 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 12.3618 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 12.4372 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 12.5126 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 12.5879 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 12.6633 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 12.7387 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 12.8141 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 12.8894 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 12.9648 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 13.0402 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 13.1156 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 13.1910 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 13.2663 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 13.3417 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 13.4171 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 13.4925 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 13.5678 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 13.6432 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 13.7186 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 13.7940 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 13.8693 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 13.9447 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 14.0201 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 14.0955 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 14.1709 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 14.2462 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 14.3216 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 14.3970 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 14.4724 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 14.5477 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 14.6231 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 14.6985 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 14.7739 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 14.8492 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div><div> 14.9246 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div>
<div> 15.0000 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN</div></div><div><br></div><div><br></div><div><br></div><div><br></div><div>B sahoo</div><div>barc mumbai</div>