Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Semiconductor Boltzmann equation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

A fast implicit solver for semiconductor models in one space dimension

Several different approaches are proposed for solving fully implicit discretizations of a simplified Boltzmann-Poisson system with a linear relaxation-type collision kernel. This system models the evolution of free electrons in semiconductor devices under a low-density assumption. At each implicit time step, the discretized system is formulated as a fixed-point problem, which can then be solved with a variety of methods. A key algorithmic component in all the approaches considered here is a recently developed sweeping algorithm for Vlasov-Poisson systems. A synthetic acceleration scheme has been implemented to accelerate the convergence of iterative solvers by using the solution to a drift-diffusion equation as a preconditioner. The performance of four iterative solvers and their accelerated variants has been compared on problems modeling semiconductor devices with various electron mean-free-path.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2001

A Hydrodynamic Theory for Spatially Inhomogeneous Semiconductor Lasers: Microscopic Approach

Starting from the microscopic semiconductor Bloch equations (SBEs) including the Boltzmann transport terms in the distribution function equations for electrons and holes, we derived a closed set of diffusion equations for carrier densities and temperatures with self-consistent coupling to Maxwell's equation and to an effective optical polarization equation. The coherent many-body effects are included within the screened Hartree-Fock approximation, while scatterings are treated within the second Born approximation including both the in- and out-scatterings. Microscopic expressions for electron-hole (e-h) and carrier-LO (c-LO) phonon scatterings are directly used to derive the momentum and energy relaxation rates. These rates expressed as functions of temperatures and densities lead to microscopic expressions for self- and mutual-diffusion coefficients in the coupled density-temperature diffusion equations. Approximations for reducing the general two-component description of the electron-hole plasma (EHP) to a single-component one are discussed. In particular, we show that a special single-component reduction is possible when e-h scattering dominates over c-LO phonon scattering. The ambipolar diffusion approximation is also discussed and we show that the ambipolar diffusion coefficients are independent of e-h scattering, even though the diffusion coefficients of individual components depend sensitively on the e-h scattering rates. Our discussions lead to new perspectives into the roles played in the single-component reduction by the electron-hole correlation in momentum space induced by scatterings and the electron-hole correlation in real space via internal static electrical field. Finally, the theory is completed by coupling the diffusion equations to the lattice temperature equation and to the effective optical polarization which in turn couples to the laser field.

Li, Jianzhong↗
NASA NTRS2002

Microscopic Foundation and Simulation of Coupled Carrier-Temperature Diffusions in Semiconductor Lasers

A typical semiconductor-based optoelectronic device, such as a diode laser, consists of three subsystems: an optical field, an electron-hole plasma (EHP), and a host crystal lattice. The physics of such a device involves the interplay of optical, electrical and thermal processes. A proper description of such a device requires that all three processes are treated on equal footing and in a self-consistent fashion. Furthermore, since a semiconductor laser has intrinsic spatial inhomogeneity, such a self-consistency naturally leads to a set of partial differential equations in space and time. There is a significant lacking of research interest and results on the transport aspects of optical devices in the literature with only a few exceptions. Even the most important carrier diffusion coefficient has not been properly derived and studied so far for optically excited plasma, while most of the work adopted results from electronics community where heavily doped semiconductors with mainly one type of carriers are dealt with. The corresponding transport equation for plasma energy or temperature has received even less attention. In this talk we describe our recent results on such a self-consistent derivation of temperature and carrier-density diffusion equations coupled with the lasing process. Starting from the microscopic semiconductor Bloch equations (SBEs) including the Boltzmann transport terms in the distribution function equations for electrons and holes, we derived a closed set of diffusion equations for carrier densities and temperatures with self-consistent coupling to Maxwell's equation and to an effective optical polarization equation. The coherent many-body effects are included within the screened Hartree-Fock approximation, while scatterings are treated within the second Born approximation including both the in- and out-scatterings. Microscopic expressions for electron-hole (e-h) and carrier-LO (c-LO) phonon scatterings are directly used to derive the momentum and energy relaxation rates. These rates expressed as functions of temperatures and densities lead to microscopic expressions for self- and mutual-diffusion coefficients in the coupled density-temperature diffusion equations. Approximations for reducing the general two-component description of the electron-hole plasma (EHP) to a single-component one are discussed. In particular, we show that a special single-component reduction is possible when e-h scattering dominates over c-LO phonon scattering. The ambipolar diffusion approximation is also discussed and we show that the ambipolar diffusion coefficients are independent of e-h scattering, even though the diffusion coefficients of individual components depend sensitively on the e-h scattering rates. Our discussions lead to new perspectives into the roles played in the single-component reduction by the electron-hole correlation in momentum space induced by scatterings and the electron-hole correlation in real space via internal static electrical field. Finally, the theory is completed by coupling the diffusion equations to the lattice temperature equation and to the effective optical polarization which in turn couples to the laser field. The equations derived above are implemented in various limiting cases to a typical diode laser to study the consequences of nonlinear diffusion and the cross diffusion terms on laser behavior, especially the dynamic behavior of a diode laser under modulation. Detailed results will be presented by comparing with the standard rate equation results.

Li, J.↗
NASA NTRS1991

Energy Models for One-Carrier Transport in Semiconductor Devices

Moment models of carrier transport, derived from the Boltzmann equation, made possible the simulation of certain key effects through such realistic assumptions as energy dependent mobility functions. This type of global dependence permits the observation of velocity overshoot in the vicinity of device junctions, not discerned via classical drift-diffusion models, which are primarily local in nature. It was found that a critical role is played in the hydrodynamic model by the heat conduction term. When ignored, the overshoot is inappropriately damped. When the standard choice of the Wiedemann-Franz law is made for the conductivity, spurious overshoot is observed. Agreement with Monte-Carlo simulation in this regime required empirical modification of this law, or nonstandard choices. Simulations of the hydrodynamic model in one and two dimensions, as well as simulations of a newly developed energy model, the RT model, are presented. The RT model, intermediate between the hydrodynamic and drift-diffusion model, was developed to eliminate the parabolic energy band and Maxwellian distribution assumptions, and to reduce the spurious overshoot with physically consistent assumptions. The algorithms employed for both models are the essentially non-oscillatory shock capturing algorithms. Some mathematical results are presented and contrasted with the highly developed state of the drift-diffusion model.

Jerome, Joseph W.↗
DOE OSTI2021

Analysis of a New Implicit Solver for a Semiconductor Model

Here, we present and analyze a new iterative solver for implicit discretizations of a simplified Boltzmann--Poisson system. The algorithm builds on recent work that incorporated a sweeping algorithm for the Vlasov--Poisson equations as part of nested inner-outer iterative solvers for the Boltzmann--Poisson equations. The new method eliminates the need for nesting and requires only one transport sweep per iteration. It arises as a new fixed-point formulation of the discretized system which we prove to be contractive for a given electric potential. We also derive an accelerator to improve the convergence rate for systems in the drift-diffusion regime. We numerically compare the efficiency of the new solver, with and without acceleration, against a recently developed nested iterative solver.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Magnetotransport in semiconductors and two-dimensional materials from first principles

Here we demonstrate a first-principles method to study magnetotransport in materials by solving the Boltzmann transport equation (BTE) in the presence of an external magnetic field. Our approach employs ab initio electron-phonon interactions and takes spin-orbit coupling into account. We apply our method to various semiconductors (Si and GaAs) and two-dimensional (2D) materials (graphene) as representative case studies. The magnetoresistance, Hall mobility, and Hall factor in Si and GaAs are in very good agreement with experiments. In graphene, our method predicts a large magnetoresistance, consistent with experiments. Analysis of the steady-state electron occupations in graphene shows the dominant role of optical phonon scattering and the breaking of the relaxation time approximation. Our paper provides a detailed understanding of the microscopic mechanisms governing magnetotransport coefficients, establishing the BTE in a magnetic field as a broadly applicable first-principles tool to investigate transport in semiconductors and 2D materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

High-throughput screening of 2D materials identifies p-type monolayer WS2 as potential ultra-high mobility semiconductor

Abstract 2D semiconductors offer a promising pathway to replace silicon in next-generation electronics. Among their many advantages, 2D materials possess atomically-sharp surfaces and enable scaling the channel thickness down to the monolayer limit. However, these materials exhibit comparatively lower charge carrier mobility and higher contact resistance than 3D semiconductors, making it challenging to realize high-performance devices at scale. In this work, we search for high-mobility 2D materials by combining a high-throughput screening strategy with state-of-the-art calculations based on the ab initio Boltzmann transport equation. Our analysis singles out a known transition metal dichalcogenide, monolayer WS 2 , as the most promising 2D semiconductor, with the potential to reach ultra-high room-temperature hole mobilities in excess of 1300 cm 2 /Vs should Ohmic contacts and low defect densities be achieved. Our work also highlights the importance of performing full-blown ab initio transport calculations to achieve predictive accuracy, including spin–orbital couplings, quasiparticle corrections, dipole and quadrupole long-range electron–phonon interactions, as well as scattering by point defects and extended defects.

Chemistry↗
NASA NTRS1977

Quantum statistical theory of semiconductor junctions in thermal equilibrium

Free carrier and electric field distributions of one-dimensional semiconductor junctions are evaluated using a quantum mechanical phase-space distribution and its corresponding Boltzmann equation. Attention is given to quantum and exchange corrections in cases of high doping concentrations when carrier densities become degenerate. Quantitative differences between degenerate and classical junction characteristics, e.g., maximum electric field and built-in voltage and carrier concentration within the transition region, are evaluated numerically.

Von Roos, O.↗
DOE OSTI2021

Origin of the background absorption in carbon nanotubes: Phonon-assisted excitonic continuum

Excitonic effects in 1D semiconductors can be qualitatively different from those in higher dimensions. In particular, the Sommerfeld factor, the ratio of the above-band-edge excitonic continuum absorption to free electron-hole pair generation, has been shown to be less than 1 (suppressed) in 1D systems while it is larger than 1 (enhanced) in 2D and 3D systems. Strong continuum suppression indeed exists in semiconducting single-wall carbon nanotubes, a prototypical 1D semiconductor. However, absorption spectra for carbon nanotubes are typically fit with a combination of Lorentzians and a polynomial background baseline with little physical meaning. We performed absorption measurements in aligned single-chirality (6,5) carbon nanotube films. The obtained spectra were fit with our theoretical model obtained by solving the Boltzmann scattering equation involving fifty-nine different types of transitions among three different types of quasiparticles. Specifically, we took into account microscopic interactions between photons, phonons, and excitons, including their dispersions, which unambiguously demonstrated that the background absorption is due to phonon-assisted transitions from the semiconductor vacuum to finite-momentum continuum states of excitons. The excellent agreement we obtained between experiment and theory suggests that our numerical technique can be seamlessly extended to compute strongly out-of-equilibrium many-body dynamics and time-resolved spectra in low-dimensional materials.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Anisotropic-strain-enhanced hole mobility in GaN by lattice matching to ZnGeN 2 and MgSiN 2

The key obstacle toward realizing integrated gallium nitride (GaN) electronics is its low hole mobility. Here, we explore the possibility of improving the hole mobility of GaN via epitaxial matching to II–IV nitride materials that have recently become available, namely, ZnGeN 2 and MgSiN 2 . We perform state-of-the-art calculations of the hole mobility of GaN using the ab initio Boltzmann transport equation. We show that effective uniaxial compressive strain of GaN along the [$1\bar{1}00$] by lattice matching to ZnGeN 2 and MgSiN 2 results in the inversion of the heavy hole band and split-off hole band, thereby lowering the effective hole mass in the compression direction. We find that lattice matching to ZnGeN 2 and MgSiN 2 induces an increase in the room-temperature hole mobility by 50% and 260% as compared to unstrained GaN, respectively. Further, examining the trends as a function of strain, we find that the variation in mobility is highly nonlinear; lattice matching to a hypothetical solid solution of Zn 0.75 Ge 0.75 Mg 0.25 Si 0.25 N 2 would already increase the hole mobility by 160%.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Electron-phonon drag enhancement of transport properties from a fully coupled ab initio Boltzmann formalism

We present a combined treatment of the nonequilibrium dynamics and transport of electrons and phonons by carrying out ab initio calculations of the fully coupled electron and phonon Boltzmann transport equations. We find that the presence of mutual drag between the two carriers causes the thermopower to be enhanced and dominated by the transport of phonons, rather than electrons as in the traditional semiconductor picture. Drag also strongly boosts the intrinsic electron mobility, thermal conductivity and the Lorenz number. Impurity scattering is seen to suppress the drag enhancement of the thermal and electrical conductivities, while having weak effects on the enhancement of the Lorenz number and thermopower. We demonstrate these effects in n-doped 3C-SiC at room temperature, and explain their origins. This work establishes the roles of microscopic scattering mechanisms in the emergence of strong drag effects in the transport of the interacting electron-phonon gas.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Electron mobility of SnO 2 from first principles

The transparent conducting oxide SnO 2 is a wide bandgap semiconductor that is easily n-type doped and widely used in various electronic and optoelectronic applications. Experimental reports of the electron mobility of this material vary widely depending on the growth conditions and doping concentrations. In this work, we calculate the electron mobility of SnO 2 from first principles to examine the temperature and doping concentration dependence and to elucidate the scattering mechanisms that limit transport. We include both electron–phonon scattering and electron-ionized impurity scattering to accurately model scattering in a doped semiconductor. Furthermore, we find a strongly anisotropic mobility that favors transport in the direction parallel to the c-axis. At room temperature and intrinsic carrier concentrations, the low-energy polar-optical phonon modes dominate scattering, while ionized-impurity scattering dominates above 10 18 cm −3 .

36 MATERIALS SCIENCE↗
DOE OSTI2022

First-principles ionized-impurity scattering and charge transport in doped materials

Scattering of carriers with ionized impurities governs charge transport in doped semiconductors. However, electron interactions with ionized impurities cannot be fully described with quantitative first-principles calculations, so their understanding relies primarily on simplified models. Here we show an ab initio approach to compute the interactions between electrons and ionized impurities or other charged defects. It includes the short- and long-range electron-defect (e-d) interactions on equal footing and allows for efficient interpolation of the e-d matrix elements. Here we combine the e-d and electron-phonon interactions in the Boltzmann transport equation to compute the carrier mobilities in doped silicon over a wide range of temperature and doping concentrations, seamlessly spanning the defect- and phonon-limited transport regimes. The individual contributions of the defect- and phonon-scattering mechanisms to the carrier relaxation times and mean-free paths are analyzed. Our method provides a powerful tool to study electronic interactions in doped materials. It broadens the scope of first-principles transport calculations, enabling studies of a wide range of doped semiconductors and oxides with application to electronics, energy, and quantum technologies.

36 MATERIALS SCIENCE↗
DOE OSTI2024

High Field Transport in (Ultra) Wide Bandgap Semiconductors: Diamond Versus Cubic GaN

We provide an analysis of nonlinear transport in diamond and cubic GaN (c-GaN) with emphasis on the different types of optical phonon scattering [i.e., optical deformation potential (ODP) scattering versus polar optic phonon (POP) scattering] limiting the carrier velocity. Both types of carrier mobilities and carrier saturation (peak) velocities in diamond and c-GaN as functions of different doping types and concentrations are obtained by directly solving the Boltzmann equation. Our model indicates that the nonrandomizing nature of POP scattering causes carrier temperature cooling, resulting in higher carrier drift velocity than with ODP scattering. This effect, in addition to the small carrier effective masses and large optical phonon energies, is responsible for the higher peak velocities in c-GaN, compared to carrier drift velocity in diamond.

42 ENGINEERING↗
DOE OSTI2021

Accelerated discovery of a large family of quaternary chalcogenides with very low lattice thermal conductivity

The development of efficient thermal energy management devices such as thermoelectrics and barrier coatings often relies on compounds having low lattice thermal conductivity (κ l ). Here, we present the computational discovery of a large family of 628 thermodynamically stable quaternary chalcogenides, AMM'Q 3 (A = alkali/alkaline earth/post-transition metals; M/M' = transition metals, lanthanides; Q = chalcogens) using high-throughput density functional theory (DFT) calculations. We validate the presence of low κ l in these materials by calculating κ l of several predicted stable compounds using the Peierls–Boltzmann transport equation. Our analysis reveals that the low κ l originates from the presence of either a strong lattice anharmonicity that enhances the phonon-scatterings or rattler cations that lead to multiple scattering channels in their crystal structures. Our thermoelectric calculations indicate that some of the predicted semiconductors may possess high energy conversion efficiency with their figure-of-merits exceeding 1 near 600 K. Our predictions suggest experimental research opportunities in the synthesis and characterization of these stable, low κ l compounds.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Intermediate polaronic charge transport in organic crystals from a many-body first-principles approach

Abstract Charge transport in organic molecular crystals (OMCs) is conventionally categorized into two limiting regimes − band transport, characterized by weak electron-phonon (e-ph) interactions, and charge hopping due to localized polarons formed by strong e-ph interactions. However, between these two limiting cases there is a less well understood intermediate regime where polarons are present but transport does not occur via hopping. Here we show a many-body first-principles approach that can accurately predict the carrier mobility in this intermediate regime and shed light on its microscopic origin. Our approach combines a finite-temperature cumulant method to describe strong e-ph interactions with Green-Kubo transport calculations. We apply this parameter-free framework to naphthalene crystal, demonstrating electron mobility predictions within a factor of 1.5−2 of experiment between 100 and 300 K. Our analysis reveals the formation of a broad polaron satellite peak in the electron spectral function and the failure of the Boltzmann equation in the intermediate regime.

36 MATERIALS SCIENCE↗
DOE OSTI2024

An Asymptotic Preserving Discontinuous Galerkin Method for a Linear Boltzmann Semiconductor Model

A key property of the linear Boltzmann semiconductor model is that as the collision frequency tends to infinity, the phase space density $f$ = $f$ ($x, v, t$) converges to an isotropic function $M (v)$$ρ$$(x, t)$, called the drift-diffusion limit, where $M$ is a Maxwellian and the physical density $ρ$ satisfies a second-order parabolic PDE known as the drift-diffusion equation. Numerical approximations that mirror this property are said to be asymptotic preserving. In this paper we build a discontinuous Galerkin method to the semiconductor model, and we show this scheme is both uniformly stable in $ε$, where 1/$ε$ is the scale of the collision frequency, and asymptotic preserving. Here in particular, we discuss what properties the discrete Maxwellian must satisfy in order for the schemes to converge in $ε$ to an accurate $h$-approximation of the drift-diffusion limit. Discrete versions of the drift-diffusion equation and error estimates in several norms with respect to $ε$ and the spacial resolution are also included.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Ab initio calculation of carrier mobility in semiconductors including ionized-impurity scattering

The past decade has seen the emergence of ab initio computational methods for calculating phonon-limited carrier mobilities in semiconductors with predictive accuracy. More realistic calculations ought to take into account additional scattering mechanisms such as, for example, impurity and grain-boundary scattering. Here, in this paper, we investigate the effect of ionized-impurity scattering on the carrier mobility. We model the analytical impurity potential parameterized from first principles by a collection of randomly distributed Coulomb scattering centers, and we include this relaxation channel into the ab initio Boltzmann transport equation, as implemented in the EPW code. We demonstrate this methodology by considering silicon, silicon carbide, and gallium phosphide, for which detailed experimental data are available. Our calculations agree well with experiments over a broad range of temperatures and impurity concentrations. For each compound investigated here, we compare the relative importance of electron-phonon scattering and ionized-impurity scattering, and we critically assess the reliability of Matthiessen's rule. We also show that an accurate description of dielectric screening and carrier effective masses can improve quantitative agreement with experiments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗