Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “first principles calculation”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Thermodynamic modeling of the Pd–Zn system with uncertainty quantification and its implication to tailor catalysts

Pd–Zn intermetallic catalysts show encouraging combinations of activity and selectivity on well-defined active site ensembles. Thermodynamic description of the Pd–Zn system, delineating phase boundaries and enumerating site occupancies within intermediate alloy phases, is essential to determining the ensembles of Pd–Zn atoms as a function of composition and temperature. Combining the present extensive first-principles calculations based on density functional theory (DFT) and available experimental data, the Pd–Zn system was remodeled using the CALculation of PHAse Diagrams (CALPHAD) approach. High throughput modeling tools with uncertainty quantification, i.e., ESPEI and PyCalphad, were incorporated in the phase analysis. Here, the site occupancies across the γ phase composition region were given special attention. A four-sublattice model was used for the γ phase owing to its four Wyckoff positions, i.e., the outer tetrahedral (OT) site 8c, the inner tetrahedral (IT) site 8c, the octahedral (OH) site 12e, and the cuboctahedral (CO) site 24g. The site fractions of Pd and Zn calculated from the present thermodynamic model show the occupancy preference of Pd in the OT and OH sublattices in agreement with experimental observations. The force constants obtained from DFT-based phonon calculations further supports the tendency of Pd occupying the OH sublattice compared with the IT and CO sublattice. The catalytic ensembles changing from Pd monomers (Pd 1 ) to trimers (Pd 3 ) on the surface of γ phase are attributed to the increase of Pd occupancy in the OH sublattice.

36 MATERIALS SCIENCE↗
DOE OSTI2022

High thermal conductivity and ultrahigh thermal boundary conductance of homoepitaxial AlN thin films

Wurtzite aluminum nitride (AlN) has attracted increasing attention for high-power and high-temperature operations due to its high piezoelectricity, ultrawide-bandgap, and large thermal conductivity k. The k of epitaxially grown AlN on foreign substrates has been investigated; however, no thermal studies have been conducted on homoepitaxially grown AlN. In this study, the thickness dependent k and thermal boundary conductance G of homoepitaxial AlN thin films were systematically studied using the optical pump–probe method of frequency-domain thermoreflectance. Our results show that k increases with the thickness and k values are among the highest reported for film thicknesses of 200 nm, 500 nm, and 1 μm, with values of 71.95, 152.04, and 195.71 W/(mK), respectively. Our first-principles calculations show good agreement with our measured data. Remarkably, the G between the epilayer and the substrate reported high values of 328, 477, 1180, and 2590 MW/(m 2 K) for sample thicknesses of 200 nm, 500 nm, 1 μm, and 3 μm, respectively. The high k and ultrahigh G of homoepitaxially grown AlN are very promising for efficient heat dissipation, which helps in device design and has advanced applications in micro-electromechanical systems, ultraviolet photonics, and high-power electronics.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Spin-Triplet Excitonic Insulator: The Case of Semihydrogenated Graphene

While various excitonic insulators have been studied in the literature, due to the perceived too-small spin splitting, spin-triplet excitonic insulator is rare. In two-dimensional systems such as a semihydrogenated graphene (known as graphone), however, it is possible, as revealed by first-principles calculations coupled with Bethe-Salpeter equation. Here, the critical temperature, given by an effective Hamiltonian, is 11.5 K. While detecting excitonic insulators is still a daunting challenge, the condensation of triplet excitons will result in spin superfluidity, which can be directly measured by a transport experiment. Nonlocal dielectric screening also leads to an unexpected phenomenon, namely, an indirect-to-direct transition crossover between single-particle band and exciton dispersion in the semihydrogenated graphene, which offers yet another test by experiment.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Control of polymorphism during epitaxial growth of hyperferroelectric candidate LiZnSb on GaSb (111)B

A major challenge for ferroelectric devices is the depolarization field, which competes with and often destroys long-range polar order in the limit of ultrathin films. Recent theoretical predictions suggest a new class of materials, termed hyperferroelectics, that should be robust against the depolarization field and enable ferroelectricity down to the monolayer limit. Here, the authors demonstrate the epitaxial growth of hexagonal LiZnSb, one of the hyperferroelectric candidate materials, by molecular-beam epitaxy on GaSb (111)B substrates. Due to the high volatility of all three atomic species, they find that LiZnSb can be grown in an adsorption-controlled window, using an excess zinc flux. Within this window, the desired polar hexagonal phase is stabilized with respect to a competing cubic polymorph, as revealed by x-ray diffraction and transmission electron microscopy measurements. First-principles calculations show that for moderate amounts of epitaxial strain and moderate concentrations of Li vacancies, the cubic LiZnSb phase is lower in formation energy than the hexagonal phase, but only by a few millielectronvolts per formula unit. Therefore, they suggest that kinetics plays a role in stabilizing the desired hexagonal phase at low temperatures. Their results provide a path toward experimentally demonstrating ferroelectricity and hyperferroelectricity in a new class of ternary intermetallic compounds.

42 ENGINEERING↗
DOE OSTI2024

Heterostructural Alloy Phase Diagram for (Cd 1-x Zn x ) 3 As 2

Alloying the topological semimetal Cd 3 As 2 with Zn 3 As 2 provides a potential route for controlling the electronic properties. We predict the alloy phase diagram from first-principles calculations, considering that both end members have a crystal structure derived from the antifluorite lattice, but with different arrangements of the unoccupied cation sites. To overcome the limitations of the regular solution approximation and to include short-range order effects, we perform Monte Carlo simulations, parameterize the temperature dependence of the mixing enthalpy ΔH m , and perform thermodynamic integration of the free energy. The resulting phase diagram exhibits features that are unique to heterostructural alloy systems and provides computational predictions of solubility limits and composition ranges that are stable against spinodal decomposition.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Iodine Vacancies do not Cause Nonradiative Recombination in Halide Perovskites

The iodine vacancy (V I ) has frequently been discussed as a strong nonradiative recombination center in halide perovskites. This proposition was mainly based on the presence of charge-state transition levels in the band gap, as found in early first-principles calculations. In this work, we perform accurate hybrid-density-functional calculations for V I in CsPbI 3 , CsSnI 3 , and CsGeI 3 and find that V I does not have any transition levels in the band gap in CsPbI 3 , in contrast to the results from calculations based on semilocal functionals. The iodine vacancy V I does introduce levels in the band gap in CsSnI 3 and CsGeI 3 , but our explicitly computed nonradiative capture coefficients demonstrate that V I has a negligible impact on nonradiative recombination. Our study corrects a misunderstanding of the role of V I in the iodide-based perovskites, and shifts the focus toward identifying and mitigating actual recombination centers in order to further improve the optoelectronic performance.

14 SOLAR ENERGY↗
DOE OSTI2020

Prospects for n -type doping of (Al x Ga 1– x ) 2 O 3 alloys

Here we systematically explore the properties of group-IV (C, Si, Ge, and Sn) and transition metal (Hf, Zr, and Ta) dopants substituting on the cation site in (Al x Ga 1– x ) 2 O 3 (AlGO) alloys using first-principles calculations with a hybrid functional. In Ga 2 O 3 , each of these dopants acts as a shallow donor. In Al 2 O 3 , they are deep defects characterized by the formation of either DX centers or positive- U (+/0) levels. Combining our calculations of dopant charge-state transition levels with information of the AlGO alloy band structure, we estimate the critical Al composition at which each dopant transitions from being a shallow to a deep donor. We identify Si to be the most efficient dopant to achieve n -type conductivity in high Al-content AlGO alloys, acting as a shallow donor over the entire predicted stability range for AlGO solid solution alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Giant linear magnetoelectric effect at the morphotropic phase boundary of epitaxial Sr 0.5 Ba 0.5 MnO 3 films

First-principles calculations are conducted to compute linear magnetoelectric coupling coefficients in epitaxial (001) Sr 0.5 Ba 0.5 MnO 3 films. A large enhancement of different linear magnetoelectric elements is found in a strained-induced morphotropic phase boundary region. Such enhancement is demonstrated to originate from the behavior of the dielectric susceptibility, thanks to a simple phenomenological model that is presently shown to be relevant and accurate. Finally, this work can thus provide a promising approach towards designing highly desired single-phase multiferroic with a colossal magnetoelectric conversion.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1991

Global Expression for Representing Diatomic Potential-Energy Curves

A three-parameter expression that gives an accurate fit to diatomic potential curves over the entire range of separation for charge transfers between 0 and 1. It is based on a generalization of the universal binding-energy relation of Smith et al. (1989) with a modification that describes the crossover from a partially ionic state to the neutral state at large separations. The expression is tested by comparison with first-principles calculations of the potential curves ranging from covalently bonded to ionically bonded. The expression is also used to calculate spectroscopic constants form a curve fit to the first-principles curves. A comparison is made with experimental values of the spectroscopic constants.

Ferrante, John↗
DOE OSTI2020

Hidden role of Bi incorporation in nonradiative recombination in methylammonium lead iodide

Bismuth incorporation has been shown to significantly decrease the performance of perovskite solar cells, which has led to the common belief that Bi Pb is a strong nonradiative recombination center in hybrid perovskites. Using rigorous first-principles calculations, we demonstrate that Bi Pb is not an efficient recombination center. However, Bi acts as a donor, and pushes the Fermi level closer to the conduction band; this shift promotes the formation of iodine interstitials, which are the actual nonradiative recombination centers. Here, these insights explain why Bi incorporation is detrimental for the photovoltaic performance. More generally, it draws attention to the detrimental role unintentionally incorporated impurities can play, not by acting as nonradiative recombination centers themselves, but by shifting the Fermi level and thereby promoting the formation of efficiency-killing defects.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Inert Gas Element as Active Infrared‐Absorption Source and Donor in Silicon for Forbidden‐Wavelength Sensing

Intrinsic silicon (Si) is forbidden for infrared (IR) sensing at the communication wavelength like 1.31 or 1.55 µm due to the well‐known bandgap limitation. In this work, an unexpected physical picture of using argon (Ar) is identified, which is usually inert to the surrounding chemical environment and used as a protective agent in semiconductor processing, to overcome the IR‐sensing‐forbidden problem in Si. Here, in this work, it is shown by an analysis of a dynamic secondary ion mass spectrometer that such a Si, when exposed to laser pulse in Ar gas, can contain a very high dose of Ar up to 10 20 cm −3 even after 1300 days. First‐principles calculations, molecular dynamics, and Hall effect measurements reveal that, due to both steric and dynamic repulsions by Ar orbitals to Si dangling bonds, the Ar‐filled‐vacancy produces a much wider defect band inside the gap, which is not only responsible for strong infrared absorption, but also causes a significant increase in n‐type conductivity, both in line with experiments. The study proves that originally inert elements in fact can act as active impurities in semiconductors for advanced applications, which updates the current knowledge of defect physics.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Structural disorder-driven topological phase transition in noncentrosymmetric BiTeI

In this work, we investigate using local structural disorder to induce a topologically nontrivial phase in a solid state system. Using first-principles calculations, we introduce structural disorder in the trivial insulator BiTeI and observe the emergence of a topological insulating phase. By modifying the bonding environments, the crystal-field splitting is enhanced, with spin-orbit interactions producing a band inversion in the bulk electronic structure. Analysis of the Wannier charge centers and the surface electronic structure reveals a strong topological insulator with Dirac surface states. Finally, we propose a prescription for inducing topological states from disorder in crystalline materials. Understanding how local environments produce topological phases is a key step for predicting disordered and amorphous topological materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Observation of paramagnetic spin-degeneracy lifting in EuZn 2 Sb 2

Taken together, time-reversal and spatial inversion symmetries impose a twofold spin degeneracy of the electronic states in crystals. In centrosymmetric materials, this degeneracy can be lifted by introducing magnetism, either via an externally applied field or through internal magnetization. However, a correlated alignment of spins, even in the paramagnetic phase, can lift the spin degeneracy of electronic states. Here, we report an in-depth study of the electronic band structure of the Eu-ternary pnictide EuZn 2 Sb 2 through a combination of high-resolution angle-resolved photoemission spectroscopy measurements and first-principles calculations. An analysis of the photoemission line shapes over a range of incident photon energies and sample temperatures is shown to reveal the presence of band spin-degeneracy lifting in the paramagnetic phase. Our angle-resolved photoemission spectroscopy results are in good agreement with theoretical ferromagnetic-phase calculations, which indicates the importance of ferromagnetic fluctuations in the system. Through our calculations, we predict that spin-polarized bands in EuZn 2 Sb 2 generate a single pair of Weyl nodes. Our observation of band splittingin EuZn 2 Sb 2 provides a key step toward realizing time-reversal symmetry breaking physics in the absence of long-range magnetic order.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Convert Widespread Paraelectric Perovskite to Ferroelectrics

While nature provides a plethora of perovskite materials, only a few exhibit large ferroelectricity and possibly multiferroicity. The majority of perovskite materials have the nonpolar CaTiO 3 (CTO) structure, limiting the scope of their applications. Based on the effective Hamiltonian model as well as first-principles calculations, we propose a general thin-film design method to stabilize the functional BiFeO 3 (BFO)- type structure, which is a common metastable structure in widespread CTO-type perovskite oxides. It is found that the improper antiferroelectricity in CTO-type perovskite and ferroelectricity in BFO-type perovskite have distinct dependences on mechanical and electric boundary conditions, both of which involve oxygen octahedral rotation and tilt. Thus the above difference can be used to stabilize the highly polar BFO-type structure in many CTO-type perovskite materials.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Anisotropy of the magnetic and transport properties of EuZn 2 As 2

Several recent studies have shown that the anisotropy in the magnetic structure of EuCd 2 As 2 plays a significant role in stabilizing the Weyl nodes. Therefore, to investigate the relationship between magnetic anisotropy and Weyl physics, we present a comparative study between EuZn 2 As 2 and EuCd 2 As 2 that are isostructural but with different magnetic anisotropy. We performed structural analysis, electronic transport, and magnetization experiments on millimeter-sized single crystals of EuZn 2 As 2 , and compared the results to those of EuCd 2 As 2 . By combining the first principle calculations and neutron diffraction experiment, we identify the magnetic ground state of EuZn 2 As 2 as A-type antiferromagnetic order with a transition temperature (T_N = 19.6 K) twice that of EuCd 2 As 2 . Like EuCd 2 As 2 , the negative magnetoresistance of EuZn 2 As 2 is observed after suppressing the resistivity peak at T_N with increasing fields. However, the anisotropy in both transport and magnetization are much reduced in EuZn 2 As 2 . The difference could be ascribed to the weaker spin-orbit coupling, more localized d orbitals, and a larger contribution from the Eu s orbitals in the zinc compound, as suggested by the electronic band calculations. The same band structure effect could be also responsible for the observation of a smaller nonlinear anomalous Hall effect in EuZn 2 As 2 compared to EuCd 2 As 2 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Design of High-Performance Lead-Free Quaternary Antiperovskites for Photovoltaics via Ion Type Inversion and Anion Ordering

The emergence of halide double perovskites significantly increases the compositional space for lead-free and air-stable photovoltaic absorbers compared to halide perovskites. Nevertheless, most halide double perovskites exhibit oversized band gaps (>1.9 eV) or dipole-forbidden optical transition, which are unfavorable for efficient single-junction solar cell applications. The current device performance of halide double perovskite is still inferior to that of lead-based halide perovskites, such as CH 3 NH 3 PbI 3 (MAPbI 3 ). Here, by ion type inversion and anion ordering on perovskite lattice sites, two new classes of pnictogen-based quaternary antiperovskites with the formula of X 6 B 2 AA' and X 6 BB'A 2 are designed. Phase stability and tunable band gaps in these quaternary antiperovskites are demonstrated based on first-principles calculations. Further photovoltaic-functionality-directed screening of these materials leads to the discovery of 5 stable compounds (Ca 6 N 2 AsSb, Ca 6 N 2 PSb, Sr 6 N 2 AsSb, Sr 6 N 2 PSb, and Ca 6 NPSb 2 ) with suitable direct band gaps, small carrier effective masses and low exciton binding energies, and dipole-allowed strong optical absorption, which are favorable properties for a photovoltaic absorber material. The calculated theoretical maximum solar cell efficiencies based on these five compounds are all larger than 29%, comparable to or even higher than that of the MAPbI 3 based solar cell. Finally, our work reveals the huge potential of quaternary antiperovskites in the optoelectronic field and provides a new strategy to design lead-free and air-stable perovskite-based photovoltaic absorber materials.

three-dimensional electronic transport↗
DOE OSTI2021

Phonon- and defect-limited electron and hole mobility of diamond and cubic boron nitride: A critical comparison

Diamond and cBN are two of the most promising ultra-wide bandgap semiconductors for applications in high-power high-frequency electronic devices. Despite extensive studies on carrier transport in these materials, there are large discrepancies in their reported carrier mobilities. In this work, we investigate the phonon- and dopant-limited electron and hole mobilities of cBN and diamond with atomistic first-principles calculations in order to understand their fundamental upper bounds to carrier transport. Our results show that although the phonon-limited electron mobilities are comparable between cBN and diamond, the hole mobility is significantly lower in cBN due to its heavier hole effective mass. Moreover, although lattice scattering dominates the mobility at low doping, neutral impurity scattering becomes the dominant scattering mechanism at higher dopant concentrations due to the high dopant ionization energies. Our analysis provides critical insights and reveals the intrinsic upper limits to the carrier mobilities of diamond and cBN as a function of doping and temperature for applications in high-power electronic devices.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Strain effects on Auger–Meitner recombination in silicon

We study the effects of compressive and tensile biaxial strain on direct and phonon-assisted Auger–Meitner recombination (AMR) in silicon using first-principles calculations. We find that the application of strain has a non-trivial effect on the AMR rate. For most AMR processes, the application of strain increases the AMR rate. However, the recombination rate for the AMR process involving two holes and one electron is suppressed by 38% under tensile strain. We further analyze the specific phonon contributions that mediate the phonon-assisted AMR mechanism, demonstrating the increased anisotropy under strain. In conclusion, our results indicate that the application of tensile strain increases the lifetime of minority electron carriers in p-type silicon and can be leveraged to improve the efficiency of silicon devices.

36 MATERIALS SCIENCE↗