It was discussed in Section 5.1 that the enhancement of the bulk low-field electron mobility saturates at around 1.7 [Dhar05b]. In order to maintain the desired mobility enhancement, the g-type and f-type coupling constants had to be adjusted to the values stated in Table 5.1. In addition, it was required to adjust the acoustic deformation potential from its original value of 8.9 eV [Jungemann03a] for analytical band Monte Carlo to 8.5 eV for full-band Monte Carlo simulation. The effect of impact ionization was neglected for the field regime investigated.
Fig. 5.11 presents the velocity-field characteristics for unstrained and strained Si for different field directions as obtained from Monte Carlo simulations. Also displayed are the results from Bufler [Bufler97], Canali [Canali71], Smith [Smith80], Fischer [Fischer99], and Ismail [Ismail93]. The simulation results agree well with measured data from Smith for the [111] field direction and with Canali for the [100] field direction for the unstrained case and with Jungemann [Jungemann03a], Bufler, and Ismail for the strained case.
The velocity-field characteristics for field along the [001] direction exhibit
an untypical form for high strain levels. This phenomenon can be explained as
follows. For field along [001] direction the -valleys are lowered in
energy with increasing strain and have the longitudinal mass in the field
direction. These valleys are located at a scaled distance of 0.85 and 1.15 from
the center of the first Brillouin zone and are separated by an energy barrier
of 129 meV at the X-point (Fig. 5.13). The average velocity in
the left and right valley and also the average of these velocities are shown
in Fig. 5.14. For low-fields, electrons in both valleys
are slightly displaced with respect to the valley minima. This results in the
initial velocity increase for both valleys shown in
Fig. 5.14. However, as the field increases, electrons in
both valleys gain energy, and electrons from one valley can surpass the energy
barrier and drift to the valley in the next Brillouin zone. As sketched in
Fig. 5.13, there are more electrons populating the right side
of the double valley than the left side, giving rise to a slight increase in
average velocity. If only the left valley is considered, there are more
electrons populating the left edge of the single valley resulting in a negative
valley velocity, as shown in Fig. 5.14.
Parameter | Units | E [100] | E [110] | E [11![]() |
E [001] | E [101] | |||
![]() |
[10![]() ![]() |
1.026 | 1.058 | 1.042 |
![]() |
[1] | 1.085 | 1.2475 | 1.273 |
Parameter | Units | E [100] | E [001] | E[110] | E [101] | E [11![]() |
![]() |
[10![]() ![]() ![]() |
![]() |
33.731 | 1.4988 | 11.739 | 11.067 |
![]() |
[eV![]() |
![]() |
![]() |
0.22885 | ![]() |
![]() |
![]() |
[eV] | 0.37994 | ![]() |
0.45615 | ![]() |
![]() |
![]() |
[1] | 1.6239 | 1.0333 | 1.5468 | 6.3401 | 4.7611 |
![]() |
[10![]() ![]() |
2.1254 | 6.3369 | 0.6651 | 4.2133 | 5.4664 |
![]() |
[10![]() ![]() ![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
[1] | 1.3707 | 2.6051 | 1.3869 | 2.4453 | 3.4612 |
![]() |
[eV![]() |
![]() |
![]() |
0.61215 | ![]() |
![]() |
Parameter | Units | E [100] | E [001] | E [110] | E [101] | E [11![]() |
![]() |
[10![]() ![]() ![]() |
![]() |
10.822 | ![]() |
2.2239 | 3.5825 |
![]() |
[eV![]() |
0.472 | 0.5135 | 0.21785 | ![]() |
![]() |
![]() |
[eV] | 0.47701 | ![]() |
1.0876 | ![]() |
1.9381 |
![]() |
[1] | 3.1569 | 2.1639 | ![]() |
![]() |
5.5323 |
![]() |
[10![]() ![]() |
7.6075 | 3.7613 | 5.8913 | 0.25071 | 1.1382 |
![]() |
[10![]() ![]() ![]() |
1.471 | 2.7214 | 1.3928 | 0.92226 | ![]() |
![]() |
[1] | 3.815 | 1.163 | 4.7754 | 1.4471 | 0.7351 |
![]() |
[eV![]() |
2.9118 | 4.8595 | 5.2425 | 0.14618 | 5.2995 |
Parameter | Units | E [101] | E [101] | E [11![]() |
E [11![]() |
![]() |
![]() |
![]() |
![]() |
||
![]() |
[10![]() ![]() ![]() |
4.0975 | 2.8429 | 3.53 | 3.4858 |
![]() |
[eV![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
[eV] | ![]() |
![]() |
![]() |
![]() |
![]() |
[1] | 0.75896 | ![]() |
0.76209 | 2.2891 |
![]() |
[10![]() ![]() |
3.9982 | 5.9335 | 4.223 | 6.5366 |
![]() |
[10![]() ![]() ![]() |
![]() |
1.9583 | ![]() |
1.8587 |
![]() |
[1] | 1.8204 | 1.9209 | 2.0921 | 1.668 |
![]() |
[eV![]() |
![]() |
3.5664 | ![]() |
5.4713 |
Fig. 5.16b shows a comparison of the velocity components and the total velocity as obtained from the interpolation and Monte Carlo simulations for field along the [111] direction for uniaxial tensile stressed Si. In this case the perpendicular velocity is negative.