Si transport problem with carrier concentration

General discussion around the EPW software

Moderator: stiwari

Post Reply
weihua-xiao
Posts: 11
Joined: Sun Jun 28, 2020 6:54 am
Affiliation: Condensed Matter Physics

Si transport problem with carrier concentration

Post 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

maxwhite
Posts: 8
Joined: Sat Feb 06, 2021 7:18 am
Affiliation: CIOMP

Re: Si transport problem with carrier concentration

Post 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 :?:

weihua-xiao
Posts: 11
Joined: Sun Jun 28, 2020 6:54 am
Affiliation: Condensed Matter Physics

Re: Si transport problem with carrier concentration

Post by weihua-xiao »

my mistake
Si is not a polar material .
I had modified epwinput, But the result is almost the same as before.

maxwhite
Posts: 8
Joined: Sat Feb 06, 2021 7:18 am
Affiliation: CIOMP

Re: Si transport problem with carrier concentration

Post 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!!

weihua-xiao
Posts: 11
Joined: Sun Jun 28, 2020 6:54 am
Affiliation: Condensed Matter Physics

Re: Si transport problem with carrier concentration

Post 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

Post Reply