The impact ionization rate in silicon is numerically derived from wave functions and energy band structure based on an empirical pseudopotential method. The calculated impact ionization rate is well fitted to an analytical formula with a power exponent of 4.6, indicating soft threshold of impact ionization rate, which originates from the complexity of the Si band structure. The calculated impact ionization rate shows strong anisotropy at low electron energy (ε<3 eV), while it becomes isotropic at higher energy. Numerical calculation also reveals that the average energy of secondary generated carriers depends linearly on the primary electron energy at the moment of their generation. A full band Monte Carlo simulation using the newly derived impact ionization rate demonstrates that calculated quantum yield and ionization coefficient agree well with reported experimental data.
The physics of electron transport in bulk silicon is investigated by using a newly developed Monte Carlo simulator which improves the state-of-the-art treatment of hot carrier transport. (1) The full band structure of the semiconductor was computed by using an empirical-pseudopotential method. (2) A phonon dispersion curve was obtained from an adiabatic bond-charge model. (3) Electron-phonon scattering was computed by using a rigid pseudo-ion model. The calculated scattering rate is consistent with the full band structure and the phonon dispersion curve of silicon, thus leaving no adjustable parameters such as deformation potential coefficients. (4) The impact-ionization rate was calculated by using Fermi’s golden rule directly from the full band structure. We took into account the dielectric function depending on both wave vector and transition energy in the numerical calculation of the rate. The impact-ionization rate obtained in the present study strongly depends on both wave vector and band index of the conduction electron, which is ignored by the traditional Keldysh formula. (5) In the simulator, the final state of a scattering electron is determined in such a way as to conserve both energy and momentum in scattering processes. The simulated results, under the steady-state conditions as well as under the nonequilibrium conditions, are presented and compared with experimental results. Special attention is focused on anisotropic transport during velocity overshoot. Quantitative agreement between calculated and experimental results confirms the validity of the newly developed Monte Carlo simulator and the physical models that were used.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.