EPW Forum

Forum for EPW code
https://forum.epw-code.org/

Si transport problem with carrier concentration

https://forum.epw-code.org/viewtopic.php?t=1402

Page 1 of 1

Si transport problem with carrier concentration

Posted: Tue Mar 23, 2021 1:50 am
by weihua-xiao
Dear all
Recently, I plan to calculate mobility of Si with different carrier concentration, but the result seem wrong and inconsistent with this paper https://journals.aps.org/prb/abstract/1 ... .97.121201
I use EPW v 5.2.0
following is my input file, how should i modify it.

scf.in

Code: Select all

&control
  calculation = 'scf'
  prefix = 'si'
  outdir = './'
  pseudo_dir = '/home/kaike/test/QE/2-si/pp'
  etot_conv_thr = 1.0d-10
  forc_conv_thr = 1.0d-8
/
&system
  ibrav = 0
  celldm(1) = 10.3527
  nat = 2, ntyp = 1
  ecutwfc = 100
  occupations = 'smearing', smearing = 'gauss', degauss = 0.0001
  nbnd = 8
/
&electrons
  conv_thr = 1d-12
  mixing_beta=0.3
  diagonalization = 'david'
/

ATOMIC_SPECIES
  Si  28.0855 Si.SG15.PBE.UPF

CELL_PARAMETERS (alat)
  -0.499620504   0.000000000   0.499620504
  -0.000000000   0.499620504   0.499620504
  -0.499620504   0.499620504   0.000000000

ATOMIC_POSITIONS {crystal}
  Si   0.0000000000d0   0.0000000000d0   0.0000000000d0
  Si   0.2500000000d0   0.2500000000d0   0.2500000000d0

K_POINTS {automatic}
  15 15 15 0 0 0
ph.in

Code: Select all

--
&inputph
  prefix = 'si'
  epsil = .true.
  fildyn = 'si.dyn'
  ldisp = .true.
  fildvscf = 'dvscf'
  nq1 = 6, nq2 = 6, nq3 = 6
  tr2_ph = 1.0d-12
/
nscf.in

Code: Select all

&control
  calculation = 'nscf'
  prefix = 'si'
  outdir = './'
  pseudo_dir = '/home/kaike/test/QE/2-si/pp'
  etot_conv_thr = 1.0d-10
  forc_conv_thr = 1.0d-8
/
&system
  ibrav = 0
  celldm(1) = 10.3527
  nat = 2, ntyp = 1
  ecutwfc = 100
  occupations = 'smearing', smearing = 'gauss', degauss = 0.0001
  nbnd = 8
/
&electrons
  conv_thr = 1d-12
  mixing_beta=0.3
  diagonalization = 'david'
/

ATOMIC_SPECIES
  Si  28.0855 Si.SG15.PBE.UPF

CELL_PARAMETERS (alat)
  -0.499620504   0.000000000   0.499620504
  -0.000000000   0.499620504   0.499620504
  -0.499620504   0.499620504   0.000000000

ATOMIC_POSITIONS {crystal}
  Si   0.0000000000d0   0.0000000000d0   0.0000000000d0
  Si   0.2500000000d0   0.2500000000d0   0.2500000000d0

K_POINTS crystal
216
  0.00000000  0.00000000  0.00000000  4.629630e-03
  0.00000000  0.00000000  0.16666667  4.629630e-03
  0.00000000  0.00000000  0.33333333  4.629630e-03
  0.00000000  0.00000000  0.50000000  4.629630e-03
  0.00000000  0.00000000  0.66666667  4.629630e-03
  0.00000000  0.00000000  0.83333333  4.629630e-03
  0.00000000  0.16666667  0.00000000  4.629630e-03
  0.00000000  0.16666667  0.16666667  4.629630e-03
  0.00000000  0.16666667  0.33333333  4.629630e-03
  0.00000000  0.16666667  0.50000000  4.629630e-03
  0.00000000  0.16666667  0.66666667  4.629630e-03
  0.00000000  0.16666667  0.83333333  4.629630e-03
  0.00000000  0.33333333  0.00000000  4.629630e-03
  0.00000000  0.33333333  0.16666667  4.629630e-03
  0.00000000  0.33333333  0.33333333  4.629630e-03
  0.00000000  0.33333333  0.50000000  4.629630e-03
  0.00000000  0.33333333  0.66666667  4.629630e-03
  0.00000000  0.33333333  0.83333333  4.629630e-03
  0.00000000  0.50000000  0.00000000  4.629630e-03
  0.00000000  0.50000000  0.16666667  4.629630e-03
  0.00000000  0.50000000  0.33333333  4.629630e-03
  0.00000000  0.50000000  0.50000000  4.629630e-03
  0.00000000  0.50000000  0.66666667  4.629630e-03
  0.00000000  0.50000000  0.83333333  4.629630e-03
  0.00000000  0.66666667  0.00000000  4.629630e-03
  0.00000000  0.66666667  0.16666667  4.629630e-03
  0.00000000  0.66666667  0.33333333  4.629630e-03
  0.00000000  0.66666667  0.50000000  4.629630e-03
  0.00000000  0.66666667  0.66666667  4.629630e-03
  0.00000000  0.66666667  0.83333333  4.629630e-03
  0.00000000  0.83333333  0.00000000  4.629630e-03
  0.00000000  0.83333333  0.16666667  4.629630e-03
  0.00000000  0.83333333  0.33333333  4.629630e-03
  0.00000000  0.83333333  0.50000000  4.629630e-03
  0.00000000  0.83333333  0.66666667  4.629630e-03
  0.00000000  0.83333333  0.83333333  4.629630e-03
  0.16666667  0.00000000  0.00000000  4.629630e-03
  0.16666667  0.00000000  0.16666667  4.629630e-03
  0.16666667  0.00000000  0.33333333  4.629630e-03
  0.16666667  0.00000000  0.50000000  4.629630e-03
  0.16666667  0.00000000  0.66666667  4.629630e-03
  0.16666667  0.00000000  0.83333333  4.629630e-03
  0.16666667  0.16666667  0.00000000  4.629630e-03
  0.16666667  0.16666667  0.16666667  4.629630e-03
  0.16666667  0.16666667  0.33333333  4.629630e-03
  0.16666667  0.16666667  0.50000000  4.629630e-03
  0.16666667  0.16666667  0.66666667  4.629630e-03
  0.16666667  0.16666667  0.83333333  4.629630e-03
  0.16666667  0.33333333  0.00000000  4.629630e-03
  0.16666667  0.33333333  0.16666667  4.629630e-03
  0.16666667  0.33333333  0.33333333  4.629630e-03
  0.16666667  0.33333333  0.50000000  4.629630e-03
  0.16666667  0.33333333  0.66666667  4.629630e-03
  0.16666667  0.33333333  0.83333333  4.629630e-03
  0.16666667  0.50000000  0.00000000  4.629630e-03
  0.16666667  0.50000000  0.16666667  4.629630e-03
  0.16666667  0.50000000  0.33333333  4.629630e-03
  0.16666667  0.50000000  0.50000000  4.629630e-03
  0.16666667  0.50000000  0.66666667  4.629630e-03
  0.16666667  0.50000000  0.83333333  4.629630e-03
  0.16666667  0.66666667  0.00000000  4.629630e-03
  0.16666667  0.66666667  0.16666667  4.629630e-03
  0.16666667  0.66666667  0.33333333  4.629630e-03
  0.16666667  0.66666667  0.50000000  4.629630e-03
  0.16666667  0.66666667  0.66666667  4.629630e-03
  0.16666667  0.66666667  0.83333333  4.629630e-03
  0.16666667  0.83333333  0.00000000  4.629630e-03
  0.16666667  0.83333333  0.16666667  4.629630e-03
  0.16666667  0.83333333  0.33333333  4.629630e-03
  0.16666667  0.83333333  0.50000000  4.629630e-03
  0.16666667  0.83333333  0.66666667  4.629630e-03
  0.16666667  0.83333333  0.83333333  4.629630e-03
  0.33333333  0.00000000  0.00000000  4.629630e-03
  0.33333333  0.00000000  0.16666667  4.629630e-03
  0.33333333  0.00000000  0.33333333  4.629630e-03
  0.33333333  0.00000000  0.50000000  4.629630e-03
  0.33333333  0.00000000  0.66666667  4.629630e-03
  0.33333333  0.00000000  0.83333333  4.629630e-03
  0.33333333  0.16666667  0.00000000  4.629630e-03
  0.33333333  0.16666667  0.16666667  4.629630e-03
  0.33333333  0.16666667  0.33333333  4.629630e-03
  0.33333333  0.16666667  0.50000000  4.629630e-03
  0.33333333  0.16666667  0.66666667  4.629630e-03
  0.33333333  0.16666667  0.83333333  4.629630e-03
  0.33333333  0.33333333  0.00000000  4.629630e-03
  0.33333333  0.33333333  0.16666667  4.629630e-03
  0.33333333  0.33333333  0.33333333  4.629630e-03
  0.33333333  0.33333333  0.50000000  4.629630e-03
  0.33333333  0.33333333  0.66666667  4.629630e-03
  0.33333333  0.33333333  0.83333333  4.629630e-03
  0.33333333  0.50000000  0.00000000  4.629630e-03
  0.33333333  0.50000000  0.16666667  4.629630e-03
  0.33333333  0.50000000  0.33333333  4.629630e-03
  0.33333333  0.50000000  0.50000000  4.629630e-03
  0.33333333  0.50000000  0.66666667  4.629630e-03
  0.33333333  0.50000000  0.83333333  4.629630e-03
  0.33333333  0.66666667  0.00000000  4.629630e-03
  0.33333333  0.66666667  0.16666667  4.629630e-03
  0.33333333  0.66666667  0.33333333  4.629630e-03
  0.33333333  0.66666667  0.50000000  4.629630e-03
  0.33333333  0.66666667  0.66666667  4.629630e-03
  0.33333333  0.66666667  0.83333333  4.629630e-03
  0.33333333  0.83333333  0.00000000  4.629630e-03
  0.33333333  0.83333333  0.16666667  4.629630e-03
  0.33333333  0.83333333  0.33333333  4.629630e-03
  0.33333333  0.83333333  0.50000000  4.629630e-03
  0.33333333  0.83333333  0.66666667  4.629630e-03
  0.33333333  0.83333333  0.83333333  4.629630e-03
  0.50000000  0.00000000  0.00000000  4.629630e-03
  0.50000000  0.00000000  0.16666667  4.629630e-03
  0.50000000  0.00000000  0.33333333  4.629630e-03
  0.50000000  0.00000000  0.50000000  4.629630e-03
  0.50000000  0.00000000  0.66666667  4.629630e-03
  0.50000000  0.00000000  0.83333333  4.629630e-03
  0.50000000  0.16666667  0.00000000  4.629630e-03
  0.50000000  0.16666667  0.16666667  4.629630e-03
  0.50000000  0.16666667  0.33333333  4.629630e-03
  0.50000000  0.16666667  0.50000000  4.629630e-03
  0.50000000  0.16666667  0.66666667  4.629630e-03
  0.50000000  0.16666667  0.83333333  4.629630e-03
  0.50000000  0.33333333  0.00000000  4.629630e-03
  0.50000000  0.33333333  0.16666667  4.629630e-03
  0.50000000  0.33333333  0.33333333  4.629630e-03
  0.50000000  0.33333333  0.50000000  4.629630e-03
  0.50000000  0.33333333  0.66666667  4.629630e-03
  0.50000000  0.33333333  0.83333333  4.629630e-03
  0.50000000  0.50000000  0.00000000  4.629630e-03
  0.50000000  0.50000000  0.16666667  4.629630e-03
  0.50000000  0.50000000  0.33333333  4.629630e-03
  0.50000000  0.50000000  0.50000000  4.629630e-03
  0.50000000  0.50000000  0.66666667  4.629630e-03
  0.50000000  0.50000000  0.83333333  4.629630e-03
  0.50000000  0.66666667  0.00000000  4.629630e-03
  0.50000000  0.66666667  0.16666667  4.629630e-03
  0.50000000  0.66666667  0.33333333  4.629630e-03
  0.50000000  0.66666667  0.50000000  4.629630e-03
  0.50000000  0.66666667  0.66666667  4.629630e-03
  0.50000000  0.66666667  0.83333333  4.629630e-03
  0.50000000  0.83333333  0.00000000  4.629630e-03
  0.50000000  0.83333333  0.16666667  4.629630e-03
  0.50000000  0.83333333  0.33333333  4.629630e-03
  0.50000000  0.83333333  0.50000000  4.629630e-03
  0.50000000  0.83333333  0.66666667  4.629630e-03
  0.50000000  0.83333333  0.83333333  4.629630e-03
  0.66666667  0.00000000  0.00000000  4.629630e-03
  0.66666667  0.00000000  0.16666667  4.629630e-03
  0.66666667  0.00000000  0.33333333  4.629630e-03
  0.66666667  0.00000000  0.50000000  4.629630e-03
  0.66666667  0.00000000  0.66666667  4.629630e-03
  0.66666667  0.00000000  0.83333333  4.629630e-03
  0.66666667  0.16666667  0.00000000  4.629630e-03
  0.66666667  0.16666667  0.16666667  4.629630e-03
  0.66666667  0.16666667  0.33333333  4.629630e-03
  0.66666667  0.16666667  0.50000000  4.629630e-03
  0.66666667  0.16666667  0.66666667  4.629630e-03
  0.66666667  0.16666667  0.83333333  4.629630e-03
  0.66666667  0.33333333  0.00000000  4.629630e-03
  0.66666667  0.33333333  0.16666667  4.629630e-03
  0.66666667  0.33333333  0.33333333  4.629630e-03
  0.66666667  0.33333333  0.50000000  4.629630e-03
  0.66666667  0.33333333  0.66666667  4.629630e-03
  0.66666667  0.33333333  0.83333333  4.629630e-03
  0.66666667  0.50000000  0.00000000  4.629630e-03
  0.66666667  0.50000000  0.16666667  4.629630e-03
  0.66666667  0.50000000  0.33333333  4.629630e-03
  0.66666667  0.50000000  0.50000000  4.629630e-03
  0.66666667  0.50000000  0.66666667  4.629630e-03
  0.66666667  0.50000000  0.83333333  4.629630e-03
  0.66666667  0.66666667  0.00000000  4.629630e-03
  0.66666667  0.66666667  0.16666667  4.629630e-03
  0.66666667  0.66666667  0.33333333  4.629630e-03
  0.66666667  0.66666667  0.50000000  4.629630e-03
  0.66666667  0.66666667  0.66666667  4.629630e-03
  0.66666667  0.66666667  0.83333333  4.629630e-03
  0.66666667  0.83333333  0.00000000  4.629630e-03
  0.66666667  0.83333333  0.16666667  4.629630e-03
  0.66666667  0.83333333  0.33333333  4.629630e-03
  0.66666667  0.83333333  0.50000000  4.629630e-03
  0.66666667  0.83333333  0.66666667  4.629630e-03
  0.66666667  0.83333333  0.83333333  4.629630e-03
  0.83333333  0.00000000  0.00000000  4.629630e-03
  0.83333333  0.00000000  0.16666667  4.629630e-03
  0.83333333  0.00000000  0.33333333  4.629630e-03
  0.83333333  0.00000000  0.50000000  4.629630e-03
  0.83333333  0.00000000  0.66666667  4.629630e-03
  0.83333333  0.00000000  0.83333333  4.629630e-03
  0.83333333  0.16666667  0.00000000  4.629630e-03
  0.83333333  0.16666667  0.16666667  4.629630e-03
  0.83333333  0.16666667  0.33333333  4.629630e-03
  0.83333333  0.16666667  0.50000000  4.629630e-03
  0.83333333  0.16666667  0.66666667  4.629630e-03
  0.83333333  0.16666667  0.83333333  4.629630e-03
  0.83333333  0.33333333  0.00000000  4.629630e-03
  0.83333333  0.33333333  0.16666667  4.629630e-03
  0.83333333  0.33333333  0.33333333  4.629630e-03
  0.83333333  0.33333333  0.50000000  4.629630e-03
  0.83333333  0.33333333  0.66666667  4.629630e-03
  0.83333333  0.33333333  0.83333333  4.629630e-03
  0.83333333  0.50000000  0.00000000  4.629630e-03
  0.83333333  0.50000000  0.16666667  4.629630e-03
  0.83333333  0.50000000  0.33333333  4.629630e-03
  0.83333333  0.50000000  0.50000000  4.629630e-03
  0.83333333  0.50000000  0.66666667  4.629630e-03
  0.83333333  0.50000000  0.83333333  4.629630e-03
  0.83333333  0.66666667  0.00000000  4.629630e-03
  0.83333333  0.66666667  0.16666667  4.629630e-03
  0.83333333  0.66666667  0.33333333  4.629630e-03
  0.83333333  0.66666667  0.50000000  4.629630e-03
  0.83333333  0.66666667  0.66666667  4.629630e-03
  0.83333333  0.66666667  0.83333333  4.629630e-03
  0.83333333  0.83333333  0.00000000  4.629630e-03
  0.83333333  0.83333333  0.16666667  4.629630e-03
  0.83333333  0.83333333  0.33333333  4.629630e-03
  0.83333333  0.83333333  0.50000000  4.629630e-03
  0.83333333  0.83333333  0.66666667  4.629630e-03
  0.83333333  0.83333333  0.83333333  4.629630e-03
epw.in

Code: Select all

--
&inputepw
  prefix      = 'si'
  amass(1)    = 28.0855
  outdir      = './'
  iverbosity  = 0

  elph        = .true.
  epbwrite    = .false.
  epbread     = .true.

  epwwrite    = .true.
  epwread     = .false.

  nbndsub     =  8
  nbndskip    =  0

  wannierize  = .true.
  num_iter    = 300
  iprint      = 2
  dis_win_max = 12
  dis_froz_max= 8
  dis_froz_min= -8
  proj(1)     = 'f=0,0,0:l=-3'
  proj(2)     = 'f=0.25,0.25,0.25:l=-3'
!  band_plot   = .true.
!  filkf       = './Kpath.dat'
!  filqf       = './Kpath.dat'

  elecselfen  = .false.
  phonselfen  = .false.
  a2f         = .false.
  fsthick     = 0.4 ! eV
  degaussw    = 0.001 ! eV
  lpolar      = .true.
  iterative_bte = .true.
  mob_maxiter = 300
  efermi_read = .true.
  fermi_energy = 6.2492
  epmatkqread = .false.
  scattering  = .true.
  scattering_serta = .true.
  nstemp = 3
  tempsmin = 200
  tempsmax = 400
  int_mob = .true.
  carrier = .true.
  ncarrier = -1E17
  asr_typ = 'simple'
!  mp_mesh_k = .true.

  dvscf_dir   = '/home/kaike/test/QE/2-si/phonons/725823.ph/save'
  filukk      = './si.ukk'
  nkf1        = 66
  nkf2        = 66
  nkf3        = 66

  nqf1        = 66
  nqf2        = 66
  nqf3        = 66

  nk1         = 6
  nk2         = 6
  nk3         = 6

  nq1         = 6
  nq2         = 6
  nq3         = 6
/
      16 cartesian
   0.000000000000000E+00   0.000000000000000E+00   0.000000000000000E+00
  -0.166793261417737E+00   0.166793261417737E+00  -0.166793261417737E+00
  -0.333586522835473E+00   0.333586522835473E+00  -0.333586522835473E+00
   0.500379784253210E+00  -0.500379784253210E+00   0.500379784253210E+00
   0.925888595891470E-17   0.333586522835473E+00   0.925888595891470E-17
  -0.166793261417737E+00   0.500379784253210E+00  -0.166793261417737E+00
   0.667173045670947E+00  -0.333586522835473E+00   0.667173045670947E+00
   0.500379784253210E+00  -0.166793261417737E+00   0.500379784253210E+00
   0.333586522835473E+00   0.370355438356588E-16   0.333586522835473E+00
   0.185177719178294E-16   0.667173045670947E+00   0.185177719178294E-16
   0.833966307088683E+00  -0.166793261417737E+00   0.833966307088683E+00
   0.667173045670947E+00  -0.370355438356588E-16   0.667173045670947E+00
   0.000000000000000E+00  -0.100075956850642E+01   0.000000000000000E+00
   0.667173045670947E+00  -0.333586522835473E+00   0.100075956850642E+01
   0.500379784253210E+00  -0.166793261417737E+00   0.833966307088683E+00
  -0.333586522835473E+00  -0.100075956850642E+01   0.370355438356588E-16
following is my ouput file

ncarrier = -1E15

Code: Select all

     =============================================================================================
     BTE in the self-energy relaxation time approximation (SERTA)
     =============================================================================================

     =============================================================================================
       Temp     Fermi   Elec density  Population SR                  Elec mobility
        [K]      [eV]     [cm^-3]      [e per cell]                    [cm^2/Vs]
     =============================================================================================

      200.000   6.2975   0.10000E+16   0.61550E-03    0.327202E+04    0.188211E+03    0.246924E+03
                                                                              0.188211E+03    0.361755E+04   -0.215613E+03
                                                                               0.246924E+03   -0.215613E+03    0.330406E+04
      300.000   6.1987   0.99999E+15   0.78252E-04    0.132193E+04    0.521989E+02    0.735535E+02
                                                                              0.521989E+02    0.140901E+04   -0.536017E+02
                                                                               0.735535E+02   -0.536017E+02    0.134374E+04
      400.000   6.0954   0.99989E+15   0.15429E-04    0.721963E+03    0.257305E+02    0.338835E+02
                                                                              0.257305E+02    0.758366E+03   -0.251827E+02
                                                                              0.338835E+02   -0.251827E+02    0.732404E+03

     =============================================================================================
     Start solving iterative Boltzmann Transport Equation
     =============================================================================================

     Iteration number:         1

     =============================================================================================
       Temp     Fermi   Elec density  Population SR                  Elec mobility
        [K]      [eV]     [cm^-3]      [e per cell]                    [cm^2/Vs]
     =============================================================================================

      200.000   6.2975   0.10000E+16   0.83926E-03    0.344161E+04    0.147672E+03    0.273901E+03
                                                                              0.164882E+03    0.370528E+04   -0.210492E+03
                                                                               0.280233E+03   -0.186536E+03    0.333477E+04
      300.000   6.1987   0.99999E+15   0.14562E-03    0.139229E+04    0.471541E+02    0.908287E+02
                                                                              0.508595E+02    0.145293E+04   -0.560100E+02
                                                                              0.925457E+02   -0.506374E+02    0.137158E+04
      400.000   6.0954   0.99989E+15   0.43290E-04    0.756508E+03    0.260658E+02    0.436382E+02
                                                                              0.275699E+02    0.781424E+03   -0.277408E+02
                                                                              0.443082E+02   -0.256764E+02    0.749712E+03
ncarrier = -1E17

Code: Select all

     =============================================================================================
     BTE in the self-energy relaxation time approximation (SERTA)
     =============================================================================================

     =============================================================================================
       Temp     Fermi   Elec density  Population SR                  Elec mobility
        [K]      [eV]     [cm^-3]      [e per cell]                    [cm^2/Vs]
     =============================================================================================

      200.000   6.3769   0.10000E+18   0.60870E-01    0.326529E+04    0.186969E+03    0.245763E+03
                                                                              0.186969E+03    0.360866E+04   -0.213953E+03
                                                                              0.245763E+03   -0.213953E+03    0.329758E+04
      300.000   6.3178   0.99999E+17   0.77441E-02    0.132057E+04    0.520750E+02    0.733779E+02
                                                                              0.520750E+02    0.140740E+04   -0.534463E+02
                                                                              0.733779E+02   -0.534463E+02    0.134237E+04
      400.000   6.2542   0.99990E+17   0.15273E-02    0.721499E+03    0.257031E+02    0.338368E+02
                                                                              0.257031E+02    0.757839E+03   -0.251501E+02
                                                                               0.338368E+02   -0.251501E+02    0.731930E+03

     =============================================================================================
     Start solving iterative Boltzmann Transport Equation
     =============================================================================================

     Iteration number:         1

     =============================================================================================
       Temp     Fermi   Elec density  Population SR                  Elec mobility
        [K]      [eV]     [cm^-3]      [e per cell]                    [cm^2/Vs]
     =============================================================================================

      200.000   6.3769   0.10000E+18   0.83201E-01    0.343536E+04    0.146756E+03    0.272980E+03
                                                                              0.163842E+03    0.369685E+04   -0.208995E+03
                                                                              0.279290E+03   -0.185195E+03    0.332900E+04
      300.000   6.3178   0.99999E+17   0.14463E-01    0.139090E+04    0.470758E+02    0.906660E+02
                                                                              0.507719E+02    0.145131E+04   -0.558685E+02
                                                                              0.923774E+02   -0.505138E+02    0.137026E+04
      400.000   6.2542   0.99990E+17   0.43074E-02    0.756021E+03    0.260503E+02    0.435902E+02
                                                                              0.275522E+02    0.780885E+03   -0.277119E+02
                                                                              0.442587E+02   -0.256515E+02    0.749242E+03
the mobility almost no change

Re: Si transport problem with carrier concentration

Posted: Tue Mar 30, 2021 8:17 am
by maxwhite
hi,
i notice there are some difference in epw1.in

for Si,it should be not polar materials.

So lpolar should not be counted

In epw.in "lpolar =.false."

Not sure it will help or not :?:

Re: Si transport problem with carrier concentration

Posted: Thu Apr 08, 2021 1:25 am
by weihua-xiao
my mistake
Si is not a polar material .
I had modified epwinput, But the result is almost the same as before.

Re: Si transport problem with carrier concentration

Posted: Thu Apr 08, 2021 2:18 am
by maxwhite
weihua-xiao wrote: ↑Thu Apr 08, 2021 1:25 am my mistake
Si is not a polar material .
I had modified epwinput, But the result is almost the same as before.
Hi，i find IBET in your caculation only has 1 interation,
as far as i konw ,you should change input:
"! scattering_serta = .true
scattering = .true.
iterative_bte = .true."

in that way i suppose ibte calculation will begin to iterate
there will be interations in the output;

Besides i want to know if you set "nbndskip = 0" in your epw.in
does the output call out error?
Removal of nbndskip tag has been made in EPW v5.3,

should i input "bands_skipped = 'exclude_bands =0:0' or just “！” it?
Need @hlee help here!!

Re: Si transport problem with carrier concentration

Posted: Thu Apr 08, 2021 12:11 pm
by weihua-xiao
thank you for your reply
scattering_serta = .true.
it will computes scattering rates in the self-energy relaxation time approximation.
iterative_bte = .true.
it will compute the iterative Boltzmann Transport Equation. i didn't show the complete output, EPW compute the iterative Boltzmann Transport Equation.
These two keywords are true, output will contain the calculation results of these two methods.
the output doesn't call out error
All times are UTC
Page 1 of 1

Powered by phpBB® Forum Software © phpBB Limited