In standard simulators for semiconductor devices the electrical behavior of the
devices is calculated using the drift-diffusion model. The development of more
complex and smaller devices demands for more exact simulation methods, such as
the Monte Carlo method examined within this work.
An important application of the Monte Carlo method is the study of the electrical
behavior of Silicon and Germanium. Mechanical strain can raise significantly
the carrier mobility in a semiconductor.
This effect has been utilized over the last few years to enhance performance of
CMOS technology.
This work starts with an introduction of the stress- and strain tensors. Then the symmetry properties of the band structure of the relaxed and the strained diamond lattice are presented. The irreducible domains of the Brillouin zone for band structure calculation are derived for important strain configurations. The band structure is calculated using the pseudo potential method and discretized on a mesh.
The main part of this work is about the simulation of carrier transport using the full-band Monte Carlo method. Important numerical algorithms and scattering models are presented. The scope is on new algorithms, which reduce simulation times. In this regard the generation of locally refined meshes for the Brillouin zone and the efficient implementation of rejection algorithms are explored.
Results from band structure and carrier transport calculations are discussed in the context of recent theoretical findings. In strained Silicon the degeneration of the X-valleys is lifted. As a consequence the valleys lower in energy are higher populated. Band structure calculations show, that in addition to the shift of the conduction band valleys relative to each other, the effective electron masses can be changed. This is caused by shear strain, which leads to an deformation of the valley minima. An analysis of the valence bands of a strained crystal shows that the degeneration of the heavy hole and light hole bands at the -point is lifted. Shear strain increases the hole mobility along certain directions, an effect caused by a deformation of the valence bands.
Mass manufactured CMOS transistors feature a uniaxially tensile strained channel in [110] direction. Simulations at 1.5GPa for this strain configuration show a mobility gain by a factor of 1.68 to 2410cm /Vs for electrons in bulk Silicon. For compressively strained Germanium hole mobility is raised by a factor of 2.55 up to 4790cm /Vs at 1.5GPa stress in [110] direction.
The final part of this work deals with the simulation of blocked impurity band photo detectors. These devices operate in the long wave length infrared range at temperatures below K. A photon is detected by lifting a carrier from a heavy doped impurity band by optical excitation to the conduction band. The simulator is extended by an inelastic scattering model for acoustic phonons, which is appropriate for simulations at low device temperature. A model for individual simulation of every carrier from an avalanche caused by impact ionization is also implemented. The distribution in energy and arrival time of the carrier avalanche is calculated.