nextnano.com
 nextnano³
 Download | SearchCopyright | Publications
 * password protected

 

nextnano³ software

 

  1D doped semiconductors

 

 

 
Up
 

nextnano3 - Tutorial

next generation 3D nano device simulator

1D Tutorial

Electron concentration in doped semiconductors (Si, Ge, GaAs) - Compensated semiconductors

Author: Stefan Birner

If you want to obtain the input files that are used within this tutorial, please check if you can find them in the installation directory.
If you cannot find them, please submit a Support Ticket.
1) -> 1Dn_doped_bulk_Si.in
2) -> 1Dn_doped_bulk_GaAs.in
3) -> 1Dn_doped_bulk_Ge_a.in
3) -> 1Dn_doped_bulk_Ge_b.in
4) -> 1Dn_doped_bulk_Ge_a_compensated.in
5) -> 1Dn_doped_bulk_Ge_a_compensated_two_donors.in
 


1) Effective density of states Nc(T) of the conduction band in Si and GaAs

The effective density of states Nc(T) of a conduction band is defined as
 Nc(T) = 2 (2 pi me kBT / h² )3/2

            2: spin degeneracy factor

For an effective mass of me = 1 at a temperature of 300 K, the value of Nc = 25.09*1018 cm-3.

The same equation holds for the effective density of states Nv(T) of a valence band where the valence band mass mh has to be used.

Note that some of the conduction bands are degenerate (valley degeneracy), e.g.
- the X       conduction band minimum in a zinc blende semiconductor is three-fold degenerate, so N'c(T) = 3 Nc(T). (Note: Each valley counts only 1/2 because they are at the boundary of the Brillouin zone.)
- the Delta conduction band minimum in silicon                                  is    six-fold degenerate, so N'c(T) = 6 Nc(T). (Note: Each valley counts because they are not at the boundary of the Brillouin zone.)
- the L       conduction band minimum in germanium                           is   four-fold degenerate, so N'c(T) = 4 Nc(T). (Note: Each valley counts only 1/2 because they are at the boundary of the Brillouin zone.)

 

1a) Si

In silicon, the electron mass at the X valley is an ellipsoid with a longitudinal and two transverse masses:
  conduction-band-masses = 0.156d0  0.156d0  0.156d0 ! GAMMA: m,m,m
                           1.420d0  0.130d0  0.130d0 ! L: ml,mt,mt
                           0.916d0  0.190d0  0.190d0 ! X: ml,mt,mt
The density of states mass is calculated as follows:
 
meX = me*DOS = (ml·mt·mt)1/3 = (0.916 · 0.1902)1/3 m0 = 0.321 m0
 

The value of kBT at different temperatures reads:

  • T =  10 K:       kBT =  0.86173 meV
  • T = 100 K:      kBT =  8.6173   meV
  • T = 300 K:      kBT = 25.852    meV

 

Thus for Si we obtain the following values for the effective density of states Nc(T) at the sixfold X (i.e. Delta) conduction band minimum:

  • T = 10K:       Nc = 6 * 2 (2 pi me kBT / h² )3/2 = 2 (0.321  * 0.8617*10-3 * 2.0886*1014)3/2 =   0.166626 * 1018 cm-3
  • T = 100K:      Nc = 6 * 2 (2 pi me kBT / h² )3/2 = 2 (0.321 * 8.617*10-3   * 2.0886*1014)3/2 =   5.26919   * 1018 cm-3
  • T = 300K:      Nc = 6 * 2 (2 pi me kBT / h² )3/2 = 2 (0.321 * 25.85*10-3   * 2.0886*1014)3/2 = 27.3795     * 1018 cm-3

The output for Nc is contained in the following file: Nc_cb1D.dat
The units are 1 * 1018 cm-3.
The first column refers to Nc of the conduction band at the Gamma point, the second at the L point, and the third at the X (or Delta) point.

For the valence bands, the effective density of states of the heavy (hh), light (lh) and split-off (so) hole bands are contained in the file:
Nv_vb1D.dat: Nv(hh), Nv(lh), Nv(so)

 

To vary the temperature, the following specifier has to be varied:

 $global-parameters
 !lattice-temperature =  10.0d0  ! 
10 Kelvin
 !lattice-temperature = 100.0d0  !
100 Kelvin
  lattice-temperature = 300.0d0  !
300 Kelvin
 $end_global-parameters

 

 

1b) GaAs

For GaAs where the electron mass at the Gamma point is me = 0.067 m0 we obtain:

  • T = 10K:       Nc = 2 (2 pi me kBT / h² )3/2 = 2 (0.067  * 0.8617*10-3 * 2.0886*1014)3/2 = 0.00264852 * 1018 cm-3
  • T = 100K:      Nc =2 (2 pi me kBT / h² )3/2 = 2 (0.067  * 8.617*10-3    * 2.0886*1014)3/2 = 0.0837535   * 1018 cm-3
  • T = 300K:      Nc = 2 (2 pi me kBT / h² )3/2 = 2 (0.067 * 25.85*10-3   * 2.0886*1014)3/2 = 0.435196     * 1018 cm-3

 

 

2) Electron concentration in doped semiconductors and position of Fermi level

The electron concentration is equivalent to the concentration of the ionized donors.

The Fermi level is assumed to be constant and equal to 0 eV in an equilibrium nextnano³ simulation (EF = 0).

 

2a) n-doped Si

The Si layer is n-type doped with a donor concentration of ND = 1*1017 cm-3.
The energy level of the donor is set to 100 meV below the conduction band edge.

 $doping-function
  doping-function-number = 1
  impurity-number        = 1        !
properties of this impurity type have to be specified below
  doping-concentration   = 0.1d0    ! 0.1
* 1018 cm-3 = 1 * 1017 cm-3
  only-region            = 0d0 20d0 !
xmin = 0 nm, xmax = 20 nm
 $end_doping-function

 $impurity-parameters
  impurity-number             =       !
impurity numbers labelled in doping-function
  impurity-type               = n-type   ! n-type
, p-type
  number-of-energy-levels     = 1        !
number of energy levels of this impurity (only 1 is currently allowed)
  energy-levels-relative      = 0.100d0  !
100 meV below conduction band
 !energy-levels-relative      = 0.049d0  ! n-As-in-Si
, energy relative to 'nearest' band edge (n-type -> conduction band
                                                                                                                                                           p-type -> valence band)
  degeneracy-of-energy-levels = 2        !
degeneracy of energy levels, 2 for n-type, 4 for p-type
 $end_impurity-parameters

 

At T = 0 K, the donors are neutral: ND+ = ND0

For T > 0 K, the electrons that are bound to the donors can be thermally activated into the conduction band.
It holds:

                      ND+ / ND0 = 1/g exp[ (ED-EF) / kBT] =  1/g exp[ (Ec-ED'-EF) / kBT]

where ED' is the donor level relative to the conduction band edge and g is the degeneracy factor (g = 2 for donors).
Note that the degeneracy of neutral donors is 2, whereas the degeneracy of the ionized donors is 1.

 

If the position of the Fermi level relative to the conduction band edge is known, one can calculate the electron density analytically using:

        ND = ND0 + ND+

        n = p + ND+ ~= ND+     (charge neutrality)

==> n ~= ND / ( 1 + 2 exp[ (EF-ED) / (kBT) ] ) =
             = ND / ( 1 + 2 exp[ (EF-(Ec-ED')) / (kBT) ] )                                 (Eq. 1)

 

 

The calculated electron density can be found in this file: densities/density1Del.dat
This is the total electron density, i.e. the sum over all conduction band minima at the Gamma, L and X (or Delta) points.

The file densities/density1DGamma_L_X.dat
contains the electron density for each valley (Gamma, L and X).

 

  • At T =  10 K, the electron density is 0.0000372311 * 1018 cm-3 = 3.72311 * 1013 cm-3, i.e. it is practially zero.
  • At T = 100 K, the electron density is    0.00154008 * 1018 cm-3 = 1.54008 * 1015 cm-3.
  • At T = 300 K, the electron density is      0.0784660 * 1018 cm-3 = 7.84660 * 1016 cm-3.

 

This can be approximated analytically for small temperatures, leading to a formula which is independent of the position of the Fermi level:

n ~= ( 0.5 * Nc * ND)1/2 exp[ (ED-Ec) / (2kBT) ] = ...

  • T =  100 K: n ~= ( 0.5 * 5.26919*1018 cm-3 * 1*1017 cm-3)1/2 exp[ (-100 meV) / (2 * 8.6173 meV) ] =
                          =  1.550 * 1015 cm-3.

 

2b) n-doped GaAs

For GaAs the electron densities at the Gamma conduction band minimum are as follows:

  • At T =  10 K, the electron density is 0.0000323902 * 1018 cm-3 = 3.23902 * 1013 cm-3, i.e. it is practially zero.
  • At T = 100 K, the electron density is  0.000217719 * 1018 cm-3 = 2.17719 * 1014 cm-3.
  • At T = 300 K, the electron density is     0.0190349 * 1018 cm-3 =  1.90349 * 1016 cm-3.

 

This can be approximated analytically for small temperatures, leading to a formula which is independent of the position of the Fermi level:

n ~= ( 0.5 * Nc * ND)1/2 exp[ (ED-Ec) / (2kBT) ] = ...

  • T =  100 K: n ~= ( 0.5 * 0.0837535*1018 cm-3 * 1*1017 cm-3)1/2 exp[ (-100 meV) / (2 * 8.6173 meV) ] =
                          =  1.955 * 1014 cm-3.

 

Position of the Fermi level

For small temperatures, the following relationship holds for the Fermi level:

      EF(T) = 1/2 ( ED + Ec )   + 1/2 kBT ln[ 1/2 ND/Nc ] =
              = 1/2 ( 2 Ec - ED' ) + 1/2 kBT ln[ 1/2 ND/Nc ]

where ED' is the donor level relative to the conduction band edge.

Thus for small temperatures, the distance of the Fermi level to the conduction band edge reads:

 Ec - EF(T) = 1/2 ED' - 1/2 kBT ln[ 1/2 ND/Nc ]

 

3a) The conduction band edge in Si is thus given by:

  • T =    10 K, EF= 0 eV:  Ec (10 K)   = 100 / 2 meV - 0.5 * 0.86173 meV  ln[ 1/2 1*1017 cm-3 /  0.166626*1018 cm-3 ] =
                                                      =        50 meV + 0.52 meV = 0.05052 eV
    nextnano³ result: Ec (X,  10 K) = 0.00724398 eV
    (The result differs from the analytical formuala because the electron density is practically zero, and thus of similar magnitude than the hole density!)
  • T =  100 K, EF= 0 eV:  Ec (100 K) = 100 / 2 meV - 0.5 * 8.6173 meV    ln[ 1/2 1*1017 cm-3 / 5.26919*1018 cm-3 ] =
                                                      =       50 meV + 20.07 meV = 0.07007 eV
    nextnano³ result: Ec (X, 100 K) = 0.0701253 eV

 

3b) The conduction band edge in GaAs is thus given by:

  • T =    10 K, EF= 0 eV:  Ec (10 K)   = 100 / 2 meV - 0.5 * 0.86173 meV  ln[ 1/2 1*1017 cm-3 /  0.00264852*1018 cm-3 ] =
                                                      =        50 meV - 1.27 meV = 0.04873 eV
    nextnano³ result: Ec (Gamma,  10 K) = 0.00379127 eV
    (The result differs from the analytical formuala because the electron density is practically zero, and thus of similar magnitude than the hole density!)
  • T =  100 K, EF= 0 eV:  Ec (100 K) = 100 / 2 meV - 0.5 * 8.6173 meV    ln[ 1/2 1*1017 cm-3 / 0.0837535*1018 cm-3 ] =
                                                      =       50 meV + 2.22 meV = 0.05222 eV
    nextnano³ result: Ec (Gamma, 100 K) = 0.0512862 eV

 

2a) n-doped Si

Using Eq. (1), one can calculate the electron density analytically:

  • At T =  10 K, the electron density is

    n ~= 1*1017 cm-3 / ( 1 + 2 exp[ 49.48 meV / 0.86173 meV ] ) ~ 0

    nextnano³ result: 3.72311 * 1013 cm-3, i.e. it is practially zero.
  • At T =  100 K, the electron density is

    n ~= 1*1017 cm-3 / ( 1 + 2 exp[ 29.93 meV / 8.6173 meV ] ) = 1.53 * 1015 cm-3

    nextnano³ result: 1.54008 * 1015 cm-3

 

2b) n-doped GaAs

  • At T =  10 K, the electron density is

    n ~= 1*1017 cm-3 / ( 1 + 2 exp[ 51.27 meV / 0.86173 meV ] ) ~ 0

    nextnano³ result: 3.23902 * 1013 cm-3, i.e. it is practially zero.
  • At T =  100 K, the electron density is

    n ~= 1*1017 cm-3 / ( 1 + 2 exp[ 47.78 meV / 8.6173 meV ] ) = 1.95 * 1014 cm-3

    nextnano³ result: 2.17719 * 1014 cm-3

     

Numerical details

The distance of the conduction band from the Fermi level is calculated in the following way inside nextnano³:

For details, see also this tutorial: I-V characteristics of an n-doped Si structure

The Semiconductor equation: np = ni2 = Nc Nv exp( - Egap/kBT)

Maxwell-Boltzmann Fermi-Dirac
n (T) = Nc(T) exp( (EF - Ec) / (kBT) ) n (T) = Nc(T) F1/2( (EF - Ec) / (kBT) )
p (T) = Nv(T) exp( (Ev - EF) / (kBT) ) p (T) = Nv(T) F1/2( (Ev - EF) / (kBT) )
                     F1/2 = Fermi-Dirac integral of order 1/2 multiplied by 2/SQRT(pi) (i.e. F1/2 includes the Gamma prefactor)

When using the Maxwell-Boltzmann statistics as an approximation we obtain for Ec:
Ec = kBT ln (Nc / n)
=> Note: Inside the code we make use of the Fermi-Dirac integrals (Fermi-Dirac statistics).

The distance of the valence band from the Fermi level can be calculated in the following way:
Ev = - kBT ln (Nv / p)         (Maxwell-Boltzmann statistics approximation)
=> Note: Inside the code we make use of the Fermi-Dirac integrals (Fermi-Dirac statistics).

 

 

3) Electron and hole concentration in doped Ge

This part of the tutorial is based on Fig. 5.5 and Fig. 5.6 in

      Holger T. Grahn, Introduction to Semiconductor Physics, World Scientific (2001)

 

3a) Electron and hole concentration in doped Ge as a function of temperature T

   ==> 1Dn_doped_bulk_Ge_a.in

   flow-scheme = 13: vary temperature as T

The Ge sample is n-type doped with a concentration of ND = 1 * 1017 cm-3. The ionization energy of the donor ED is set to 10 meV.
(Note that the melting point of germanium is at 1211.40 K.)

The following figure shows the electron density as a function of temperature for a Ge sample.

In the saturation regime, the charge carrier density is equal to the doping density (between 150 K and 600 K), i.e. n =~ 1 * 1017 cm-3.

For large temperatures, the intrinsic carrier density dominates, i.e. here the electron density is equal to the hole density.

The figure also shows the band gap of the Gamma, L and Delta conduction band minima as a function of temperature. The valence band edges are also shown (heavy, light and split-off hole).
One can clearly see that the band gaps decrease with increasing temperature.
This temperature dependence can be described using the Varshni formula.

For very small temperatures, the position of the Fermi level (chemical potential) is close to the L conduction band edge, i.e. at EF = EgapL - ED/2 = EcL - ED/2.
For temperatures between 600 and 900 K, the Fermi level is around the middle of the band gap.



In this example, the effective density of states at 300 K was calculated as follows:

    Nc (300 K) = 61.068 * 1018 cm-3
    Nv (300 K) =   5.411 * 1018 cm-3


This figure shows the same results as above, but this time the Fermi level is set to a constant value EF = 0 eV.
Consequently, the conduction and valence band edges are plotted with respect to the constant Fermi level.
This figure shows more clearly how the band edges align with respect to the Fermi level.

The data that has been plotted in the above shown figures can be found in this file: sweep_temperature1D.dat

 

3b) Electron and hole concentration in doped Ge as a function of inverse temperature 1/T

   ==> 1Dn_doped_bulk_Ge_b.in

   flow-scheme = 14: vary temperature as 1000/T

The same data as above is plotted but this time as a function of inverse temperature in units of 1000/T.

The slope of the electron density approaches at high temperatures (i.e. left side of the figure) the value of Egap/2 (if it holds Nc = Nv which is not the case in our example).
At 300 K,  Egap/2 = 0.332 eV.
At 600 K,  Egap/2 = 0.267 eV.
At 900 K,  Egap/2 = 0.200 eV.
At low temperatures the slope is ED/2 = 5 meV = 0.0005 eV.


 

This figure shows the same results as above, but this time the Fermi level is set to a constant value EF = 0 eV.
Consequently, the conduction and valence band edges are plotted with respect to the constant Fermi level.
This figure shows more clearly how the band edges align with respect to the Fermi level.

 

 

4) Compensated semiconductors - one donor level and one acceptor level

   ==> 1Dn_doped_bulk_Ge_a_compensated.in

   flow-scheme = 13: vary temperature as T

  • An interesting situation occurs in real semiconductors that have both, n-type and p-type doping. It is rather impossible to fabricate devices that are merely n- or p-type because unintentional doping through impurities is always present.
  • We assume the same Ge sample as above, i.e. the Ge sample is n-type doped with a concentration of ND = 1 * 1017 cm-3.
    The ionization energy of the donor ED is set to 10 meV.
    In contrast to the above shown figures, also acceptors are present with varying acceptor concentration ranging from
    - NA = 1 * 1014 cm-3    (NA << ND)  ==>  ND - NA =~ ND
    - NA = 1 * 1016 cm-3    (NA  <  ND)  ==>  ND - NA = 9*1016 cm-3
    - NA = 5 * 1016 cm-3    (NA  <  ND)  ==>  ND - NA = 5*1016 cm-3
    - NA = 7 * 1016 cm-3    (NA  <  ND)  ==>  ND - NA = 3*1016 cm-3
    - NA = 9 * 1016 cm-3    (NA  <  ND)  ==>  ND - NA = 1*1016 cm-3
    - NA = 1 * 1017 cm-3    (NA  =  ND)  ==>  ND - NA = 0  ==>  intrinisc density is recovered
    The ionization energy of the acceptor EA is set to 10 meV.
  • The following figure shows the electron density vs. temperature for different acceptor concentrations NA.
    Note that the donor concentration ND is kept constant.



    If NA = ND, then the intrinsic behavior of the electron density is recovered. In any case, the intrinsic density dominates at high temperatures.
    In the saturation regime (between 100 K and 400 K), the nominally dopant density (ND - NA) determines the electron density.

    The variation of the acceptor concentration is specified here:

      $doping-function
       doping-function-number = 2
       impurity-number        = 2         !
    p-type impurity
       doping-concentration   = 0.0001d0  ! ND >> NA ! 0.0001 * 1018 / cm3 = 1 * 1014 / cm3
     ! doping-concentration   = 0.01d0    ! ND >  NA ! 0.01   * 1018 / cm3 = 1 * 1016 / cm3
     ! doping-concentration   = 0.05d0    ! ND >  NA ! 0.05   * 1018 / cm3 = 5 * 1016 / cm3
     ! doping-concentration   = 0.07d0    ! ND >  NA ! 0.07   * 1018 / cm3 = 7 * 1016 / cm3
     ! doping-concentration   = 0.09d0    ! ND >  NA ! 0.09   * 1018 / cm3 = 9 * 1016 / cm3
     ! doping-concentration   = 0.1d0     ! ND =  NA ! 0.1    * 1018 / cm3 = 1 * 1017 / cm3
     
  • The following figure shows the Fermi level as a function of temperature for various acceptor densities NA.
    Note that the indirect band gap at the L point decreases with temperature.
    As expected, for the "intrinisic case" (ND - NA = 0  ==>  intrinisc density is recovered), the Fermi level at 0 Kelvin lies in the middle of the band gap at EF = 0.74 eV / 2 = 0.37 eV.

 

 

5) Compensated semiconductors - two donor levels and one acceptor level

   ==> 1Dn_doped_bulk_Ge_a_compensated_two_donors.in

   flow-scheme = 13: vary temperature as T

We assume a Ge sample which has two donor levels and one acceptor level.

The acceptor concentration is NA = 4 * 1016 cm-3.
The ionization energy of the acceptor EA is set to 10 meV.

The donor concentrations ND,1 and ND,2 are varied.
The ionization energy of the upper donor level ED,1 is set to 10 meV.
The ionization energy of the lower donor level ED,2 is set to 100 meV.
 

We distinguish between three cases:

  • 5a)
        ND,1 = ND,2 = 5 * 1016 cm-3  >  NA
     => ND,1 + ND,2  >  NA
     
  • 5b)
        ND,1 = 3 * 1016 cm-3  <  NA
        ND,2 = 7 * 1016 cm-3  >  NA
        ND,1 + ND,2  >  NA
     
  • 5c)   (one donor level and one acceptor level)
        ND,1 = 10 * 1016 cm-3  >  NA
        ND,2 = 0
     
  • The following figure shows the result for these three cases. 5c) is similar to the figures that have been discussed already.
    At the dopant density of ND,1 - NA = 6 * 1016 cm-3 the saturation regime forms.

    For case 5b) only the ionization of the second donor level ED,2 is relevant at low temperatures.

    For case 5a) two plateaus form:
    - for the ionization of the upper donor level ED,1: ND,1               - NA = 1 * 1016 cm-3
    - for the ionization of the lower donor level ED,2: ND,1 + ND,2 - NA = 6 * 1016 cm-3




    For comparison, for case 5a) the energetic position of the second donor level ED,2 has been varied between 70 meV, 100 meV and 130 meV.
    The energetic position of the upper donor level is at 10 meV.
     
  • The following figure shows the Fermi level for the cases 5a), 5b) and 5c).
    The valence band edge is fixed at 0 eV. The indirect band gap at the L point decreases with increasing temperature.

    For case 5b), where only the ionization of the second donor level ED,2 is relevant at low temperatures, one can clearly see the position of the Fermi level at 100 meV below the conduction band edge: EF (T = 0 K) = EC - ED,2

    For case 5a) the Fermi level moves from the upper donor level ED,1 at low temperatures to the lower donor level ED,2 at around 200 K and start to ionize the second donor level, eventually forming the second "plateau" in the electron density (see figure above for the electron density).

    At high temperatures, the Fermi levels moves approximately into the middle of the band gap.