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