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 cm s] | 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 cmseV] | 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 Vcm] | 2.1254 | 6.3369 | 0.6651 | 4.2133 | 5.4664 | |
[10 VcmeV] | ||||||
[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 cmseV] | 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 Vcm] | 7.6075 | 3.7613 | 5.8913 | 0.25071 | 1.1382 | |
[10 VcmeV] | 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 cms eV] | 4.0975 | 2.8429 | 3.53 | 3.4858 | |
[eV] | |||||
[eV] | |||||
[1] | 0.75896 | 0.76209 | 2.2891 | ||
[10 Vcm] | 3.9982 | 5.9335 | 4.223 | 6.5366 | |
[10 VcmeV] | 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.