Density functional theory (DFT) underestimates the experimental band gap—the infamous band gap problem. As the band gap defines the energy scale of defect levels, this complicates computation of charge transition energies for atomic defects. In the extreme case of narrow-gap semiconductors, the DFT band gap collapses to zero, seemingly precluding quantitative predictions of defect levels. We present a band-avoiding occupation-constrained DFT (ba-occ-DFT) approach that prevents spurious occupation of band-edge states and enables reliable total energy calculations of atomic defects. Application to indium arsenide (InAs) shows that ba-occ-DFT circumvents the band gap problem, separates band-edge errors from defect level calculations, and enables rigorous defect level predictions in a narrow-gap semiconductor despite a zero DFT band gap.
Semiconductor materials play a central role in current and future electronics technologies. From microprocessors and advanced computers to optical components, device functionality is dependent on the creation and control of point defects in semiconductors. Incorporating a fundamental understanding of defect kinetics, including formation, migration, and chemistry, is essential for advancing materials science, assessing their device impact, ensuring the reliability of modern electronics, and leveraging new materials for next-generation technologies. This article explores the kinetics of point defects from experimental observations and atomistic modeling, to dynamical multiscale descriptions of defect kinetics. A survey of experiments reveals the important role of kinetics in defect behavior during synthesis, implantation doping, radiation exposure, and long-term defect evolution, while highlighting the impact of evolving material properties on device performance. Atomistic modeling, including molecular dynamics and density functional theory, is surveyed emphasizing its ability to describe dynamical behavior and predict kinetic pathways that govern defect evolution in semiconductors. Dynamical and multiscale modeling methods that integrate experimental and atomistic defect properties into continuum-scale codes are examined for their role in bridging atomic-scale defect behavior to device-level performance. By addressing critical challenges and revealing the inherent difficulties in modeling and experimental validation, this article aims to advance the understanding of defect kinetics and provide insights into the short-term and long-term reliability of materials and devices.
Irradiation of gallium arsenide (GaAs) produces immobile vacancies and mobile interstitials. Yet, after decades of experimental investigation, the immobile Ga vacancy continues to evade detection, raising the question: where is the Ga vacancy? Static first-principles calculations predict a Ga vacancy should be readily observed. We find that short-time dynamical evolution of primary defects is the key to explaining this conundrum. Using a dynamical multiscale atomistically informed device engineering (AIDE) method, we discover that during the initial displacement damage, the Ga vacancy (3-/2-) defect level pins the Fermi level near the midgap, producing oppositely charged vacancies and interstitials. Driven by Coulomb attraction, fast As interstitials preferentially annihilate Ga vacancies. The Ga vacancy population plummets below detectable limits—and the now unpinned Fermi level recovers—before being experimentally observed. This dynamical model solves the mystery of the missing Ga vacancy and reveals the importance of a multiscale approach to explore the dynamical chemical behavior in experimentally inaccessible short-time regimes.
While first-principles density functional theory modeling has become a vital tool to investigate defect properties in semiconductors, the lack of crystalline periodicity in pseudobinary random composition alloys, such as In1-xGaxAs, complicates such analyses. We present a simulation strategy to systematically take into account the variability in the local defect environment in order to predict statistical properties of neutral intrinsic defects in In1-xGaxAs. We use a comprehensive sampling from a modest-sized 64-atom special quasirandom structure (SQS) to define a statistically representative set of defects, and use a 512-atom hypercell, a 2x2x2 supercell of SQS supercells, to achieve cell-size convergence. We articulate an equivalent site principle and describe how it constrains atomic chemical reference energies in computation of defect formation energies in pseudobinary alloys. A simple protocol for estimating reference energies for the Ga and In atoms sharing the III site succeeds in obtaining the equivalence of defects at Ga-sites and In sites in the SQS supercell, (<30 meV differences in average formation energies). For III-site defects, such as the As antisite As-III, the statistical variability in formation energies is modest, approximate to 0.1-0.2 eV. The variability in formation energy at As-site defects, such as the As vacancy v(As), can be much larger, >1 eV. The As antisite is shown to be a low-energy defect and the most likely to be present in as-grown materials, just as in GaAs. All other defects are higher-energy defects unlikely to be important in native material, but potentially important in radiation-damaged material. With a strong variability in defect energies, especially on the As-site, explicit consideration of statistical variability due to compositional randomness will be imperative for meaningful and quantitative comparisons to experiment.
Radiation-induced defects can have a significant impact on the longevity and performance of semiconductor devices. We present an Atomistically Informed Device Engineering (AIDE) method that integrates first-principles defect properties and experimentally measured parameters into a device model to dynamically simulate the defect chemistry in semiconductors. For a silicon-doped gallium arsenide (GaAs) material, we showcase three capabilities: (i) Fermi level $E_F$ movement including its component electron and hole Fermi levels, (ii) dynamical charge equilibration with the arsenic vacancy serving as an example, and a (iii) diffusion-driven reaction between Coulomb attracted gallium interstitial ($Ga_i$) and arsenic vacancy ($v_{As}$). Governed by charge carrier reactions, the electron and hole Fermi levels remained dissimilar until equilibrium was achieved at $E_F\approx1.32$ eV. The equilibrium Fermi level was verified by successfully identifying $v_{As}^{3-}$ as the most populated charge state within the arsenic vacancy defect. Lastly, a Coulomb attraction, created by the shifted Fermi level and the charge equilibration process, between $Ga_i^{1+}$ and $v_{As}^{3-}$ resulted in the formation of a doubly negative gallium antisite ($Ga_{As}^{2-}$). The AIDE method can access experimentally inaccessible short-time and low-concentration regimes, is generalizable to other more complex systems (e.g., indium gallium arsenide), and, after solving open problems in GaAs, will serve as a virtual experiment to bound estimates for difficult-to-measure physical quantities.
X-ray nanotomography is a powerful tool for the characterization of nanoscale materials and structures, but it is difficult to implement due to the competing requirements of X-ray flux and spot size. Due to this constraint, state-of-the-art nanotomography is predominantly performed at large synchrotron facilities. We present a laboratory-scale nanotomography instrument that achieves nanoscale spatial resolution while addressing the limitations of conventional tomography tools. The instrument combines the electron beam of a scanning electron microscope (SEM) with the precise, broadband X-ray detection of a superconducting transition-edge sensor (TES) microcalorimeter. The electron beam generates a highly focused X-ray spot on a metal target held micrometers away from the sample of interest, while the TES spectrometer isolates target photons with a high signal-to-noise ratio. This combination of a focused X-ray spot, energy-resolved X-ray detection, and unique system geometry enables nanoscale, element-specific X-ray imaging in a compact footprint. The proof of concept for this approach to X-ray nanotomography is demonstrated by imaging 160 nm features in three dimensions in six layers of a Cu-SiO2 integrated circuit, and a path toward finer resolution and enhanced imaging capabilities is discussed.
Using ground state density functional theory (DFT) and implementing an occupation-constrained DFT (occ-DFT) for self-consistent excited state calculations, we decipher the electronic structure of the Mn dopant and other 3d defects in GaN across the band gap. Our analysis, validated with broad agreement with defect levels (ground-state calculations) and photoluminescence data (excited-state calculations), mandates reinterpretation and reassignment of 3d defect data in GaN. The Mn-Ga defect is determined to span stable charge states from (1-) in n-type GaN through (2+) in p-type GaN. The Mn(2+) is predicted to be a d(2) ground state spin triplet defect with a singlet excited state, isoelectronic with the defect associated with the 1.19 eV photoluminescence in n-type GaN. The combined analysis of defect levels and excited states invites reassessment of all d(2)-capable dopants in GaN. We demonstrate that the 1.19 eV defect, a candidate defect for optically controlled quantum applications, cannot be the Cr(1+) assumed in literature and instead must be the V(0). The combined ground-state/excited-state DFT analysis is shown to be able to chemically fingerprint defects.
This chapter describes the design, simulation, fabrication, and characterization of a solid-state DC circuit breaker based on a normally-off, light-triggered, gallium nitride photoconductive switch combined with a cascaded, normally-on, silicon carbide junction field-effect transistor circuit leg. This design provides a very fast response time to fault events. Simulations of the various parts of the breaker and their predicted behavior in system designs have guided a first hardware demonstration. Circuit breaker voltage and current timing diagrams illustrate the interplay between different parts of the breaker and the sensitivity of the timing. The experimental performance closely matches the predicted behavior, allowing voltage and current scaling for future designs.
Using first-principles density functional theory (DFT) methods and size-converged supercell models, we analyze the electronic and atomic structure of magnetic 3d transition metal dopants in cubic gallium nitride (c-GaN). All stable defect charge states for Fermi levels across the full experimental gap are computed using a method that correctly resolves the boundary condition problem (without a jellium approximation) and eliminates finite-size errors. The resulting computed defect levels are not impacted by the DFT band-gap problem, they span a width consistent with the experimental gap rather than being limited to the single-particle DFT gap. All defects with electronically degenerate (half-metal) Td ground states are found to have significant distortions, relaxing to D2d structures driven by the Jahn-Teller instability. This leads to insulating ground states for all substitutional 3d dopants, refuting claims in the literature that +U or hybrid functional methods are required to avoid artificial half-metal results. Interpreting the dn atomic occupations within a crystal-field model and exchange splittings, we identify a systematic trend across the 3d transition metal series. Approaches to estimate excited-state energies as observed in photoluminescence from defect centers are assessed, ranging from a Koopmans-type single-particle energy interpretation to relaxed total energy differences in fully self-consistent DFT. The single-particle interpretations are found to be qualitatively predictive and the calculations are consistent with the limited available experimental data across the 3d dopant series. These results provide a baseline understanding to guide future studies and a conceptual framework within which to interpret new results.
The E3 transition in irradiated GaAs observed in deep level transient spectroscopy (DLTS) was recently discovered in Laplace-DLTS to encompass three distinct components. The component designated E3c was found to be metastable, reversibly bleached under minority carrier (hole) injection, with an introduction rate dependent upon Si doping density. It is shown through first-principles modeling that the E3c must be the intimate Sivacancy pair, best described as a Si sitting in a divacancy Sivv. The bleached metastable state is enabled by a doubly site-shifting mechanism: Upon recharging, the defect undergoes a second site shift rather returning to its original E3c-active configuration via reversing the first site shift. Identification of this defect offers insights into the short-time annealing kinetics in irradiated GaAs.
We report the first nonjellium, systematic, density functional theory (DFT) study of intrinsic and extrinsic defects and defect levels in zinc-blende (cubic) gallium nitride. We use the local moment counter charge (LMCC) method, the standard Perdew-Becke-Ernzerhoff (PBE) exchange-correlation potential, and two pseudopotentials, where the Ga 3d orbitals are either in the core (d(0)) or explicitly in the valence set (d(10)). We studied 64, 216, 512, and 1000 atom supercells, and demonstrated convergence to the infinite limit, crucial for delineating deep from shallow states near band edges, and for demonstrating the elimination of finite cell-size errors. Contrary to common claims, we find that exact exchange is not required to obtain defect levels across the experimental band gap. As was true in silicon, silicon carbide, and gallium arsenide, the extremal LMCC defect levels of the aggregate of defects yield an effective LMCC defect band gap that is within 10% of the experimental gap (3.3 eV) for both pseudopotentials. We demonstrate that the gallium vacancy is more complicated than previously reported. There is dramatic metastability-a nearest-neighbor nitrogen atom shifts into the gallium site, forming an antisite, nitrogen vacancy pair, which is more stable than the simple vacancy for positive charge states. Our assessment of the d(0) and d(10) pseudopotentials yields minimal differences in defect structures and defect levels. The better agreement of the d(0) lattice constant with experiment suggests that the more computationally economical d(0) pseudopotentials are sufficient to achieve the fidelity possible within the physical accuracy of DFT, and thereby enable calculations in larger supercells necessary to demonstrate convergence with respect to finite size supercell errors.
Structural disorder causes materials’ surface electronic properties, e.g., work function (ϕ), to vary spatially, yet it is challenging to prove exact causal relationships to underlying ensemble disorder, e.g., roughness or granularity. For polycrystalline Pt, nanoscale resolution photoemission threshold mapping reveals a spatially varying ϕ=5.70±0.03 eV over a distribution of (111) vicinal grain surfaces prepared by sputter deposition and annealing. With regard to field emission and related phenomena, e.g., vacuum arc initiation, a salient feature of the ϕ distribution is that it is skewed with a long tail to values down to 5.4 eV, i.e., far below the mean, which is exponentially impactful to field emission via the Fowler–Nordheim relation. We show that the ϕ spatial variation and distribution can be explained by ensemble variations of granular tilts and surface slopes via a Smoluchowski smoothing model wherein local ϕ variations result from spatially varying densities of electric dipole moments, intrinsic to atomic steps, that locally modify ϕ. Atomic step-terrace structure is confirmed with scanning tunneling microscopy (STM) at several locations on our surfaces, and prior works showed STM evidence for atomic step dipoles at various metal surfaces. From our model, we find an atomic step edge dipole μ=0.12 D/edge atom, which is comparable to values reported in studies that utilized other methods and materials. Our results elucidate a connection between macroscopic ϕ and the nanostructure that may contribute to the spread of reported ϕ for Pt and other surfaces and may be useful toward more complete descriptions of polycrystalline metals in the models of field emission and other related vacuum electronics phenomena, e.g., arc initiation.
The first-principles computation of the surfaces of metals is typically accomplished through slab calculations of finite thickness. The extraction of a convergent surface formation energy from slab calculations is dependent upon defining an appropriate bulk reference energy. I describe a method for an independently computed, slab-consistent bulk reference that leads to convergent surface formation energies from slab calculations that also provides realistic uncertainties for the magnitude of unavoidable nonlinear divergence in the surface formation energy with slab thickness. The accuracy is demonstrated on relaxed, unreconstructed low-index aluminum surfaces with slabs with up to 35 layers.
The stability of low-index platinum surfaces and their electronic properties is investigated with density functional theory, toward the goal of understanding the surface structure and electron emission, and identifying precursors to electrical breakdown, on nonideal platinum surfaces. Propensity for electron emission can be related to a local work function, which, in turn, is intimately dependent on the local surface structure. The $(1\ifmmode\times\else\texttimes\fi{}N)$ missing row reconstruction of the Pt(110) surface is systematically examined. The $(1\ifmmode\times\else\texttimes\fi{}3)$ missing row reconstruction is found to be the lowest in energy, with the $(1\ifmmode\times\else\texttimes\fi{}2)$ and $(1\ifmmode\times\else\texttimes\fi{}4)$ slightly less stable. In the limit of large $(1\ifmmode\times\else\texttimes\fi{}N)$ with wider (111) nanoterraces, the energy accurately approaches the asymptotic limit of the infinite Pt(111) surface. This suggests a local energetic stability of narrow (111) nanoterraces on free Pt surfaces that could be a common structural feature in the complex surface morphologies, leading to work functions consistent with those on thermally grown Pt substrates.
be amorphized by displacing C atoms exclusively and suggest that short-range disorder provides the driving force for amorphization of SiC.
of alternative UWBG semiconductor materials with improved doping and transport properties. As an example, we discuss the case of rutile GeO2 (r-GeO2), a water-insoluble GeO2 polytype, which is theoretically predicted to combine an ultra-wide gap with ambipolar dopability, high carrier mobilities, and a higher thermal conductivity than β-Ga2O3. The subsequent realization of single-crystalline r-GeO2 thin films by molecular beam epitaxy provides the opportunity to realize r-GeO2 for electronic applications. Future efforts toward the predictive discovery and design of new UWBG semiconductors include advances in first-principles theory and high-performance computing software, as well as the demonstration of controlled doping in high-quality thin films with lower dislocation densities and optimized film properties.