problems calculating the superconducting gap and Tc using The Eliashberg formalism
Posted: Sun Apr 11, 2021 2:14 pm
dear EPW user
I am still having problems calculating the superconducting gap and Tc using The Eliashberg formalism.but getting an error:
Reading supplied temperature list.
Reading xml data from directory:
./NbS2.save/
Message from routine qes_read:magnetizationType:
do_magnetization: wrong number of occurrences
Message from routine qexsd_readschema :
error reading output_obj of xsd data file
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Error in routine read_xml_file (4):
fatal error reading xml file
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
stopping ...
Any help is appreciated. the input file is :
&inputepw
prefix = 'NbS2',
amass(1) = 92.9064,
amass(2) = 32.065,
outdir = './'
ep_coupling = .true.
elph = .true.
kmaps = .false.
epbwrite = .true.
epbread = .false.
epwwrite = .true.
epwread = .false.
etf_mem = 1
nbndsub = 16,
wannierize = .true.
num_iter = 300
proj(1) = 'S:pz'
wdata(1) = 'bands_plot = .true.'
wdata(2) = 'begin kpoint_path'
wdata(3) = 'G 0.00 0.00 0.00 M 0.50 0.00 0.00'
wdata(4) = 'X 0.50 0.00 0.00 K 0.33 0.33 0.00'
wdata(5) = 'end kpoint_path'
wdata(6) = 'bands_plot_format = gnuplot'
iverbosity = 2
eps_acustic = 2.0 ! Lowest boundary for the phonon frequency
ephwrite = .true. ! Writes .ephmat files used when Eliasberg = .true.
fsthick = 0.4 ! eV
degaussw = 0.10 ! eV
nsmear = 1
delta_smear = 0.04 ! eV
degaussq = 0.5 ! meV
nqstep = 500
eliashberg = .true.
laniso = .true.
limag = .true.
lpade = .true.
conv_thr_iaxis = 1.0d-4
wscut = 1.0 ! eV Upper limit over frequency integration/summation in the Elisashberg eq
nstemp = 1 ! Nr. of temps
temps = 2.00 ! K provide list of temperetures OR (nstemp and temps = tempsmin tempsmax for even space mode)
nsiter = 500
muc = 0.16
dvscf_dir = '../Phonon/save'
nk1 = 6
nk2 = 6
nk3 = 6
nq1 = 6
nq2 = 6
nq3 = 6
mp_mesh_k = .true.
nkf1 = 20
nkf2 = 20
nkf3 = 20
nqf1 = 20
nqf2 = 20
nqf3 = 20
/
i have check the inputs or outputs of scf or nscf calculations. it's all done and there's no error. out is
total cpu time spent up to now is 3555.3 secs
End of self-consistent calculation
Number of k-points >= 100: set verbosity='high' to print the bands.
the Fermi energy is 0.3018 ev
! total energy = -320.14512623 Ry
estimated scf accuracy < 4.5E-16 Ry
smearing contrib. (-TS) = -0.47134806 Ry
internal energy E=F+TS = -319.67377817 Ry
The total energy is F=E-TS. E is the sum of the following terms:
one-electron contribution = -631.20182508 Ry
hartree contribution = 319.30152387 Ry
xc contribution = -57.71365110 Ry
ewald contribution = 49.94017414 Ry
convergence has been achieved in 35 iterations
negative rho (up, down): 1.472E-05 0.000E+00
Forces acting on atoms (cartesian axes, Ry/au):
atom 1 type 1 force = -0.04711203 0.08157864 0.85420280
atom 2 type 1 force = 0.10521036 0.00913841 -0.13833855
atom 3 type 2 force = 0.00006968 0.00000702 -0.00227483
atom 4 type 2 force = -0.00001038 0.00002062 0.00244982
atom 5 type 2 force = 0.04710880 -0.08159195 -0.85691639
atom 6 type 2 force = -0.10526642 -0.00915275 0.14087716
Total force = 1.242188 Total SCF correction = 0.000000
Computing stress (Cartesian axis) and pressure
negative rho (up, down): 1.472E-05 0.000E+00
total stress (Ry/bohr**3) (kbar) P= 60.84
-0.00010851 0.00000080 -0.00023507 -15.96 0.12 -34.58
0.00000080 -0.00015862 0.00011190 0.12 -23.33 16.46
-0.00023507 0.00011190 0.00150790 -34.58 16.46 221.82
Writing output data file ./NbS2.save/
init_run : 31.31s CPU 32.69s WALL ( 1 calls)
electrons : 3376.99s CPU 3522.51s WALL ( 1 calls)
forces : 8.51s CPU 8.85s WALL ( 1 calls)
stress : 27.28s CPU 27.81s WALL ( 1 calls)
Called by init_run:
wfcinit : 29.09s CPU 30.42s WALL ( 1 calls)
potinit : 2.03s CPU 2.05s WALL ( 1 calls)
hinit0 : 0.08s CPU 0.09s WALL ( 1 calls)
Called by electrons:
c_bands : 2967.32s CPU 3093.11s WALL ( 35 calls)
sum_band : 406.29s CPU 425.90s WALL ( 35 calls)
v_of_rho : 3.88s CPU 3.95s WALL ( 36 calls)
newd : 1.41s CPU 1.44s WALL ( 36 calls)
mix_rho : 0.08s CPU 0.09s WALL ( 35 calls)
vdW_kernel : 1.95s CPU 1.96s WALL ( 1 calls)
Called by c_bands:
init_us_2 : 41.51s CPU 49.48s WALL ( 120450 calls)
cegterg : 2772.69s CPU 2871.42s WALL ( 57750 calls)
Called by *egterg:
cdiaghg : 400.01s CPU 404.66s WALL ( 241585 calls)
h_psi : 1993.25s CPU 2077.74s WALL ( 243235 calls)
s_psi : 68.35s CPU 69.55s WALL ( 243235 calls)
g_psi : 5.59s CPU 5.66s WALL ( 183835 calls)
Called by h_psi:
h_psi:calbec : 116.78s CPU 123.97s WALL ( 243235 calls)
vloc_psi : 1786.74s CPU 1860.90s WALL ( 243235 calls)
add_vuspsi : 73.23s CPU 76.02s WALL ( 243235 calls)
General routines
calbec : 153.61s CPU 163.68s WALL ( 309235 calls)
fft : 1.20s CPU 1.25s WALL ( 2435 calls)
ffts : 0.02s CPU 0.02s WALL ( 71 calls)
fftw : 1989.43s CPU 2073.34s WALL (14243148 calls)
interpolate : 0.03s CPU 0.03s WALL ( 36 calls)
Parallel routines
fft_scatt_xy : 177.32s CPU 184.51s WALL (14245654 calls)
fft_scatt_yz : 733.79s CPU 766.10s WALL (14245654 calls)
PWSCF : 57m34.03s CPU 1h 0m WALL
This run was terminated on: 21:32: 5 11Apr2021
=------------------------------------------------------------------------------=
JOB DONE.
=------------------------------------------------------------------------------=
Is it because I didn't set the magnetic parameters ?
&control
prefix = 'NbS2'
calculation = 'scf'
restart_mode = 'from_scratch'
tstress = .true.
tprnfor = .true.
outdir = './'
wfcdir = './wfcdir'
pseudo_dir = '../'
nstep = 1000 ,
! etot_conv_thr = 3.0E-9 ,
forc_conv_thr = 1.0D-5 ,
/
&system
ibrav = 4
A = 3.482607807
C = 30
nat = 6
ntyp = 2
ecutwfc=50,
ecutrho=300,
input_dft='vdw-DF',
occupations='smearing',
smearing='gaussian',
degauss=0.1,
nspin = 1,
nbnd = 40,
lda_plus_u = .false. ,
/
&electrons
electron_maxstep = 100
conv_thr = 1.0d-14
mixing_mode = 'plain'
/
&IONS
ion_dynamics='bfgs'
/
&CELL
cell_dynamics = 'bfgs' ,
press = 0.00 ,
cell_dofree=2Dxy,
/
ATOMIC_SPECIES
Nb 92.91 Nb.pbe-nsp-van.UPF
S 32.065 S.pbe-n-rrkjus_psl.1.0.0.UPF
ATOMIC_POSITIONS crystal
Nb 0 0 9.93675
Nb 0 0 3.31225
S 1.68138 0.970743 11.5033
S 0 1.94149 1.74571
S 0 1.94149 4.87878
S 1.68138 0.970743 8.37021
K_POINTS crystal
324
0.00000000 0.00000000 0.00000000 0.00308642
0.00000000 0.05555556 0.00000000 0.00308642
0.00000000 0.11111111 0.00000000 0.00308642
I am still having problems calculating the superconducting gap and Tc using The Eliashberg formalism.but getting an error:
Reading supplied temperature list.
Reading xml data from directory:
./NbS2.save/
Message from routine qes_read:magnetizationType:
do_magnetization: wrong number of occurrences
Message from routine qexsd_readschema :
error reading output_obj of xsd data file
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Error in routine read_xml_file (4):
fatal error reading xml file
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
stopping ...
Any help is appreciated. the input file is :
&inputepw
prefix = 'NbS2',
amass(1) = 92.9064,
amass(2) = 32.065,
outdir = './'
ep_coupling = .true.
elph = .true.
kmaps = .false.
epbwrite = .true.
epbread = .false.
epwwrite = .true.
epwread = .false.
etf_mem = 1
nbndsub = 16,
wannierize = .true.
num_iter = 300
proj(1) = 'S:pz'
wdata(1) = 'bands_plot = .true.'
wdata(2) = 'begin kpoint_path'
wdata(3) = 'G 0.00 0.00 0.00 M 0.50 0.00 0.00'
wdata(4) = 'X 0.50 0.00 0.00 K 0.33 0.33 0.00'
wdata(5) = 'end kpoint_path'
wdata(6) = 'bands_plot_format = gnuplot'
iverbosity = 2
eps_acustic = 2.0 ! Lowest boundary for the phonon frequency
ephwrite = .true. ! Writes .ephmat files used when Eliasberg = .true.
fsthick = 0.4 ! eV
degaussw = 0.10 ! eV
nsmear = 1
delta_smear = 0.04 ! eV
degaussq = 0.5 ! meV
nqstep = 500
eliashberg = .true.
laniso = .true.
limag = .true.
lpade = .true.
conv_thr_iaxis = 1.0d-4
wscut = 1.0 ! eV Upper limit over frequency integration/summation in the Elisashberg eq
nstemp = 1 ! Nr. of temps
temps = 2.00 ! K provide list of temperetures OR (nstemp and temps = tempsmin tempsmax for even space mode)
nsiter = 500
muc = 0.16
dvscf_dir = '../Phonon/save'
nk1 = 6
nk2 = 6
nk3 = 6
nq1 = 6
nq2 = 6
nq3 = 6
mp_mesh_k = .true.
nkf1 = 20
nkf2 = 20
nkf3 = 20
nqf1 = 20
nqf2 = 20
nqf3 = 20
/
i have check the inputs or outputs of scf or nscf calculations. it's all done and there's no error. out is
total cpu time spent up to now is 3555.3 secs
End of self-consistent calculation
Number of k-points >= 100: set verbosity='high' to print the bands.
the Fermi energy is 0.3018 ev
! total energy = -320.14512623 Ry
estimated scf accuracy < 4.5E-16 Ry
smearing contrib. (-TS) = -0.47134806 Ry
internal energy E=F+TS = -319.67377817 Ry
The total energy is F=E-TS. E is the sum of the following terms:
one-electron contribution = -631.20182508 Ry
hartree contribution = 319.30152387 Ry
xc contribution = -57.71365110 Ry
ewald contribution = 49.94017414 Ry
convergence has been achieved in 35 iterations
negative rho (up, down): 1.472E-05 0.000E+00
Forces acting on atoms (cartesian axes, Ry/au):
atom 1 type 1 force = -0.04711203 0.08157864 0.85420280
atom 2 type 1 force = 0.10521036 0.00913841 -0.13833855
atom 3 type 2 force = 0.00006968 0.00000702 -0.00227483
atom 4 type 2 force = -0.00001038 0.00002062 0.00244982
atom 5 type 2 force = 0.04710880 -0.08159195 -0.85691639
atom 6 type 2 force = -0.10526642 -0.00915275 0.14087716
Total force = 1.242188 Total SCF correction = 0.000000
Computing stress (Cartesian axis) and pressure
negative rho (up, down): 1.472E-05 0.000E+00
total stress (Ry/bohr**3) (kbar) P= 60.84
-0.00010851 0.00000080 -0.00023507 -15.96 0.12 -34.58
0.00000080 -0.00015862 0.00011190 0.12 -23.33 16.46
-0.00023507 0.00011190 0.00150790 -34.58 16.46 221.82
Writing output data file ./NbS2.save/
init_run : 31.31s CPU 32.69s WALL ( 1 calls)
electrons : 3376.99s CPU 3522.51s WALL ( 1 calls)
forces : 8.51s CPU 8.85s WALL ( 1 calls)
stress : 27.28s CPU 27.81s WALL ( 1 calls)
Called by init_run:
wfcinit : 29.09s CPU 30.42s WALL ( 1 calls)
potinit : 2.03s CPU 2.05s WALL ( 1 calls)
hinit0 : 0.08s CPU 0.09s WALL ( 1 calls)
Called by electrons:
c_bands : 2967.32s CPU 3093.11s WALL ( 35 calls)
sum_band : 406.29s CPU 425.90s WALL ( 35 calls)
v_of_rho : 3.88s CPU 3.95s WALL ( 36 calls)
newd : 1.41s CPU 1.44s WALL ( 36 calls)
mix_rho : 0.08s CPU 0.09s WALL ( 35 calls)
vdW_kernel : 1.95s CPU 1.96s WALL ( 1 calls)
Called by c_bands:
init_us_2 : 41.51s CPU 49.48s WALL ( 120450 calls)
cegterg : 2772.69s CPU 2871.42s WALL ( 57750 calls)
Called by *egterg:
cdiaghg : 400.01s CPU 404.66s WALL ( 241585 calls)
h_psi : 1993.25s CPU 2077.74s WALL ( 243235 calls)
s_psi : 68.35s CPU 69.55s WALL ( 243235 calls)
g_psi : 5.59s CPU 5.66s WALL ( 183835 calls)
Called by h_psi:
h_psi:calbec : 116.78s CPU 123.97s WALL ( 243235 calls)
vloc_psi : 1786.74s CPU 1860.90s WALL ( 243235 calls)
add_vuspsi : 73.23s CPU 76.02s WALL ( 243235 calls)
General routines
calbec : 153.61s CPU 163.68s WALL ( 309235 calls)
fft : 1.20s CPU 1.25s WALL ( 2435 calls)
ffts : 0.02s CPU 0.02s WALL ( 71 calls)
fftw : 1989.43s CPU 2073.34s WALL (14243148 calls)
interpolate : 0.03s CPU 0.03s WALL ( 36 calls)
Parallel routines
fft_scatt_xy : 177.32s CPU 184.51s WALL (14245654 calls)
fft_scatt_yz : 733.79s CPU 766.10s WALL (14245654 calls)
PWSCF : 57m34.03s CPU 1h 0m WALL
This run was terminated on: 21:32: 5 11Apr2021
=------------------------------------------------------------------------------=
JOB DONE.
=------------------------------------------------------------------------------=
Is it because I didn't set the magnetic parameters ?
&control
prefix = 'NbS2'
calculation = 'scf'
restart_mode = 'from_scratch'
tstress = .true.
tprnfor = .true.
outdir = './'
wfcdir = './wfcdir'
pseudo_dir = '../'
nstep = 1000 ,
! etot_conv_thr = 3.0E-9 ,
forc_conv_thr = 1.0D-5 ,
/
&system
ibrav = 4
A = 3.482607807
C = 30
nat = 6
ntyp = 2
ecutwfc=50,
ecutrho=300,
input_dft='vdw-DF',
occupations='smearing',
smearing='gaussian',
degauss=0.1,
nspin = 1,
nbnd = 40,
lda_plus_u = .false. ,
/
&electrons
electron_maxstep = 100
conv_thr = 1.0d-14
mixing_mode = 'plain'
/
&IONS
ion_dynamics='bfgs'
/
&CELL
cell_dynamics = 'bfgs' ,
press = 0.00 ,
cell_dofree=2Dxy,
/
ATOMIC_SPECIES
Nb 92.91 Nb.pbe-nsp-van.UPF
S 32.065 S.pbe-n-rrkjus_psl.1.0.0.UPF
ATOMIC_POSITIONS crystal
Nb 0 0 9.93675
Nb 0 0 3.31225
S 1.68138 0.970743 11.5033
S 0 1.94149 1.74571
S 0 1.94149 4.87878
S 1.68138 0.970743 8.37021
K_POINTS crystal
324
0.00000000 0.00000000 0.00000000 0.00308642
0.00000000 0.05555556 0.00000000 0.00308642
0.00000000 0.11111111 0.00000000 0.00308642