Re: Pb with spin orbit calculation
Posted: Mon Mar 06, 2017 8:53 am
Sorry to keep bothering you, but I am still having difficulties. I am getting reasonable values of Tc, 3.5 K to 4.2 K, however the lambda values at the end are too large:
lambda : 2.5557500
lambda_tr : 2.3351044
Estimated Allen-Dynes Tc
logavg = 0.0001387 l_a2F = 2.7644099
mu = 0.10 Tc = 3.795693352742 K
mu = 0.12 Tc = 3.665318799993 K
mu = 0.14 Tc = 3.533794724152 K
mu = 0.16 Tc = 3.401194883342 K
mu = 0.18 Tc = 3.267604244930 K
mu = 0.20 Tc = 3.133120164925 K
What would cause this? The input file is below.
--
&inputepw
prefix = 'pb'
amass(1) = 207.2
outdir = './'
elph = .true.
kmaps = .false.
epbwrite = .true.
epbread = .false.
epwwrite = .true.
epwread = .false.
nbndsub = 8
nbndskip = 10
wannierize = .true.
num_iter = 300
dis_win_max = 21
dis_win_min = -3
dis_froz_min= -3
dis_froz_max= 13.5
proj(1) = 'Pb:sp3'
iverbosity = 0
elecselfen = .true.
phonselfen = .true.
parallel_k = .true.
parallel_q = .false.
fsthick = 6 ! eV
eptemp = 0.075 ! K
degaussw = 0.05 ! eV
degaussq = 0.5
a2f = .true.
dvscf_dir = './save'
nkf1 = 30
nkf2 = 30
nkf3 = 30
nqf1 = 30
nqf2 = 30
nqf3 = 30
nk1 = 6
nk2 = 6
nk3 = 6
nq1 = 6
nq2 = 6
nq3 = 6
/
16 cartesian
0.000000000000 0.000000000000 0.000000000000 0.0092593
-0.166666666667 0.166666666667 -0.166666666667 0.0740741
-0.333333333333 0.333333333333 -0.333333333333 0.0740741
0.500000000000 -0.500000000000 0.500000000000 0.0370370
0.000000000000 0.333333333333 0.000000000000 0.0555556
-0.166666666667 0.500000000000 -0.166666666667 0.2222222
0.666666666667 -0.333333333333 0.666666666667 0.2222222
0.500000000000 -0.166666666667 0.500000000000 0.2222222
0.333333333333 0.000000000000 0.333333333333 0.1111111
0.000000000000 0.666666666667 0.000000000000 0.0555556
0.833333333333 -0.166666666667 0.833333333333 0.2222222
0.666666666667 0.000000000000 0.666666666667 0.1111111
0.000000000000 -1.000000000000 0.000000000000 0.0277778
0.666666666667 -0.333333333333 1.000000000000 0.2222222
0.500000000000 -0.166666666667 0.833333333333 0.2222222
-0.333333333333 -1.000000000000 0.000000000000 0.1111111
lambda : 2.5557500
lambda_tr : 2.3351044
Estimated Allen-Dynes Tc
logavg = 0.0001387 l_a2F = 2.7644099
mu = 0.10 Tc = 3.795693352742 K
mu = 0.12 Tc = 3.665318799993 K
mu = 0.14 Tc = 3.533794724152 K
mu = 0.16 Tc = 3.401194883342 K
mu = 0.18 Tc = 3.267604244930 K
mu = 0.20 Tc = 3.133120164925 K
What would cause this? The input file is below.
--
&inputepw
prefix = 'pb'
amass(1) = 207.2
outdir = './'
elph = .true.
kmaps = .false.
epbwrite = .true.
epbread = .false.
epwwrite = .true.
epwread = .false.
nbndsub = 8
nbndskip = 10
wannierize = .true.
num_iter = 300
dis_win_max = 21
dis_win_min = -3
dis_froz_min= -3
dis_froz_max= 13.5
proj(1) = 'Pb:sp3'
iverbosity = 0
elecselfen = .true.
phonselfen = .true.
parallel_k = .true.
parallel_q = .false.
fsthick = 6 ! eV
eptemp = 0.075 ! K
degaussw = 0.05 ! eV
degaussq = 0.5
a2f = .true.
dvscf_dir = './save'
nkf1 = 30
nkf2 = 30
nkf3 = 30
nqf1 = 30
nqf2 = 30
nqf3 = 30
nk1 = 6
nk2 = 6
nk3 = 6
nq1 = 6
nq2 = 6
nq3 = 6
/
16 cartesian
0.000000000000 0.000000000000 0.000000000000 0.0092593
-0.166666666667 0.166666666667 -0.166666666667 0.0740741
-0.333333333333 0.333333333333 -0.333333333333 0.0740741
0.500000000000 -0.500000000000 0.500000000000 0.0370370
0.000000000000 0.333333333333 0.000000000000 0.0555556
-0.166666666667 0.500000000000 -0.166666666667 0.2222222
0.666666666667 -0.333333333333 0.666666666667 0.2222222
0.500000000000 -0.166666666667 0.500000000000 0.2222222
0.333333333333 0.000000000000 0.333333333333 0.1111111
0.000000000000 0.666666666667 0.000000000000 0.0555556
0.833333333333 -0.166666666667 0.833333333333 0.2222222
0.666666666667 0.000000000000 0.666666666667 0.1111111
0.000000000000 -1.000000000000 0.000000000000 0.0277778
0.666666666667 -0.333333333333 1.000000000000 0.2222222
0.500000000000 -0.166666666667 0.833333333333 0.2222222
-0.333333333333 -1.000000000000 0.000000000000 0.1111111