Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “first principles density functional theory”

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 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Cluster expansion by transfer learning for phase stability predictions

Recent progress towards universal machine-learned interatomic potentials holds considerable promise for materials discovery. Yet the accuracy of these potentials for predicting phase stability may still be limited. In contrast, cluster expansions provide accurate phase stability predictions but are computationally demanding to parameterize from first principles, especially for structures of low dimension or with a large number of components, such as interfaces or multimetal catalysts. We overcome this trade-off via transfer learning. Using Bayesian inference, we incorporate prior statistical knowledge from machine-learned and physics-based potentials, enabling us to sample the most informative configurations and to efficiently fit first-principles cluster expansions. Furthermore, this algorithm is tested on Pt:Ni, showing robust convergence of the mixing energies as a function of sample size with reduced statistical fluctuations.

36 MATERIALS SCIENCE↗
DOE OSTI2024

High-Temperature Observation of Intralayer, Interlayer, and Rydberg Excitons in Bulk Van Der Waals Alloy Single Crystals

Transition metal dichalcogenides exhibit remarkable optical properties due to the diverse number of strongly bound excitons, which can be fine-tuned by alloying. Despite a flurry of research activity in characterizing these excitons, a comprehensive and profound understanding of their behavior with temperature is lacking. Here, we report the rich spectrum of excitonic features within bulk van der Waals alloy Mo0.5 W0.5 S2 and Mo0.5 W0.5 Se2 single crystals through temperature-dependent reflectance spectroscopy and first-principles calculations. We observed Rydberg excitons and interlayer excitons in both the single crystals. Notably, we provide the first experimental evidence of highly energetic A' and B' excitons in Mo0.5 W0.5 S2 at room temperature. The strong carrier-phonon scattering significantly broadens the A', B', and interlayer excitons at room temperature in bulk Mo0.5 W0.5 S2 single crystal compared to its selenide. Our findings, supported by density functional theory and Bethe-Salpeter equation calculations, signify the crucial role of carrier-phonon interactions. These results open pathways for next-generation optoelectronic devices and quantum technologies operating at high temperature.

excitons↗
DOE OSTI2025

First-Principles Assessment of ZnTe and CdSe as Prospective Tunnel Barriers at the InAs/Al Interface

Majorana zero modes are predicted to emerge in semiconductor/ superconductor interfaces, such as InAs/Al. Majorana modes could be utilized for fault tolerant topological qubits. However, their realization is hindered by materials challenges. The coupling between the superconductor and the semiconductor may be too strong for Majorana modes to emerge, due to effective doping of the semiconductor by the metallic contact. This could be mediated by adding a tunnel barrier of controlled thickness. We use density functional theory (DFT) with Hubbard U corrections, whose values are machine-learned via Bayesian optimization (BO), to assess ZnTe and CdSe as prospective tunnel barriers for the InAs/Al interface. The results of DFT +U(BO) for ZnTe are validated by comparison to angle resolved photoemission spectroscopy (ARPES). We then study bilayer interfaces of the three semiconductors with each other and with Al, as well as trilayer interfaces with a varying number of ZnTe or CdSe layers inserted between InAs and Al. We find that 16 atomic layers of either material completely insulate the InAs from metal induced gap states (MIGS). However, ZnTe and CdSe differ significantly in their band alignment, such that ZnTe forms an effective barrier for electrons, whereas CdSe forms a barrier for holes. Because of Fermi level pinning in the conduction band at the interface, only electron transport is relevant for InAs-based Majorana devices. Therefore, ZnTe is the better choice. Based on the results of our simulations, we suggest conducting experiments with ZnTe barriers in the thickness range of 6–18 atomic layers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Defect modeling in semiconductors: the role of first principles simulations and machine learning

Abstract Point defects in semiconductors dictate their electronic and optical properties. Vacancies, interstitials, substitutional defects, and defect complexes can form in the semiconductor lattice and significantly impact its performance in applications such as solar absorption, light emission, electronics, and catalysis. Understanding the nature and energetics of point defects is essential for the design and optimization of next-generation semiconductor technologies. Here, we provide a comprehensive overview of the current state of research on point defects in semiconductors, focusing on the application of density functional theory (DFT) and machine learning (ML) in accelerating the prediction and understanding of defect properties. DFT has been instrumental in accurately calculating defect formation energies, charge transition levels, and other defect-related properties such as carrier recombination rates and lifetimes, and ion migration barriers. ML techniques, particularly neural networks, have emerged as powerful tools for enabling rapid prediction of defect properties at DFT-accuracy in order to overcome the expense of using large supercells and advanced functionals. We begin this article with a discussion of different types of point defects and complexes, their impact on semiconductor properties, and the experimental and DFT approaches typically used for their characterization. Through multiple case studies, we explore how DFT has been successfully applied to understand defect behavior across a variety of semiconductors, and how ML approaches integrated with DFT can efficiently predict defect properties and facilitate the discovery of new materials with tailored defect behavior. Overall, the advent of ‘DFT+ML’ promises to drive advancements in semiconductor technology, catalysis, and renewable energy applications, paving the way for the development of high-performance semiconductors which are defect-tolerant or have desirable dopability.

Rahman, Md Habibur (ORCID:000000027705984X)↗
DOE OSTI2024

Gold–Thiolate Nanocluster Dynamics and Intercluster Reactions Enabled by a Machine Learned Interatomic Potential

Monolayer protected metal clusters comprise a rich class of molecular systems and are promising candidate materials for a variety of applications. While a growing number of protected nanoclusters have been synthesized and characterized in crystalline forms, their dynamical behavior in solution, including prenucleation cluster formation, is not well understood due to limitations both in characterization and first-principles modeling techniques. Recent advancements in machine-learned interatomic potentials are rapidly enabling the study of complex interactions such as dynamical behavior and reactivity on the nanoscale. Here, we develop an Au-S-C-H atomic cluster expansion (ACE) interatomic potential for efficient and accurate molecular dynamics simulations of thiolate-protected gold nanoclusters (Au n (SCH 3 ) m ). Trained on more than 30,000 density functional theory calculations of gold nanoclusters, the interatomic potential exhibits ab initio level accuracy in energies and forces and replicates nanocluster dynamics including thermal vibration and chiral inversion. Long dynamics simulations (up to 0.1 μs time scale) reveal a mechanism explaining the thermal instability of neutral Au 25 (SR) 18 clusters. Specifically, we observe multiple stages of isomerization of the Au 25 (SR) 18 cluster, including a chiral isomer. Additionally, we simulate coalescence of two Au 25 (SR) 18 clusters and observe series of clusters where the formation mechanisms are critically mediated by ligand exchange in the form of [Au-S] n rings.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Electronic structure of Humble defects in Ge and Ge 0.8⁢ Si 0.2

The group-IV diamond-structure elements are known to host a variety of planar defects, including {001} planar defects in C and {001}, {111}, and {113} planar defects in Si and Ge. Among the {001} planar defects, the Humble defect, known for some time to occur in Ge, has recently also been observed in Si-Ge alloys, but the details of its electronic structure remain poorly understood. Here, in this work, we perform first-principles density-functional calculations to study Humble defects in both Ge and Ge 0.8⁢ Si 0.2 . We also measure the Si L 2,3 -edge electron energy-loss spectra both at the defect and in a bulk-like region far from the defect and compare with theoretical calculations on corresponding Si sites in our first-principles calculations. We find that inclusion of core-hole effects in the theory is essential for reproducing the observed L 2,3 edge spectra, and that once they are included, the results provide a set of fingerprints for different types of local atomic bonding environments in Ge 0.8⁢ Si 0.2 . Our first-principles calculations reveal that the Humble defects have a tendency to enlarge the electronic band gap, which may have potential uses in band engineering. The use of hybrid functionals for an improved description of the band gap in these systems is also discussed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

First-Principles Study of Interaction between Molecules and Lewis Acid Zeolites Manipulated by Injection of Energized Charge Carriers

Metal dopants in Lewis acid zeolites can be applied to manipulate the interaction between the active sites and molecular reactants. However, only a few elements can be introduced into the zeolite framework to generate strong Lewis acid sites. In this study, we propose an alternative approach to tune the interaction between zeolite and molecules by introducing energized charge carriers that are photoexcited in a complex of plasmonic metal and zeolites. Here, through density functional theory (DFT) calculations, we show that although the energized electrons do not impact the adsorption of ammonia, which is used as a probe molecule, on Zr- and Si-BEA, ammonia’s heat of adsorption can be reduced up to 60% on Sn-BEA and 20% on Ti-BEA, respectively, because of the localized population of the electrons at the Sn and Ti center. This most pronounced effect in Sn-BEA is attributed to the Sn-induced midgap states and the strong electronegativity of Sn atoms. In contrast, the energized holes can strengthen the interaction between zeolites and ammonia by at least 10%. Particularly, if an energized hole is present in the system, the heat of adsorption of ammonia within Si-BEA can increase 3-fold, shifting from weak to moderate interaction. We find that, similar to closed sites, the effect of energized charge carriers can also be extended to open sites of Lewis acid zeolites. Our study provides fundamental insights for unlocking the potentials of integrating plasmonic metals with zeolites and brings an additional degree of control to tune the interfacial interaction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Benchmarking boron carbide equation of state using computation and experiment

Boron carbide ( B 4 C ) is of both fundamental scientific and practical interest due to its structural complexity and how it changes upon compression, as well as its many industrial uses and potential for use in inertial confinement fusion (ICF) and high-energy density physics experiments. Here, we report the results of a comprehensive computational study of the equation of state (EOS) of B 4 C in the liquid, warm dense matter, and plasma phases. Our calculations are cross-validated by comparisons with Hugoniot measurements up to 61 megabar from planar shock experiments performed at the National Ignition Facility (NIF). Our computational methods include path integral Monte Carlo, activity expansion, as well as all-electron Green's function Korringa-Kohn-Rostoker and molecular dynamics that are both based on density functional theory. We calculate the pressure-internal energy EOS of B 4 C over a broad range of temperatures ( ~ 6 × 10 3 – 5 × 10 8 K) and densities (0.025–50 g / cm 3 ). We assess that the largest discrepancies between theoretical predictions are ≲ 5 % near the compression maximum at 1– 2 × 10 6 K. This is the warm-dense state in which the K shell significantly ionizes and has posed grand challenges to theory and experiment. By comparing with different EOS models, we find a Purgatorio model (LEOS 2122) that agrees with our calculations. The maximum discrepancies in pressure between our first-principles predictions and LEOS 2122 are ~ 18 % and occur at temperatures between 6 × 10 3 – 2 × 10 5 K, which we believe originate from differences in the ion thermal term and the cold curve that are modeled in LEOS 2122 in comparison with our first-principles calculations. To account for potential differences in the ion thermal term, we have developed three new equation-of-state models that are consistent with theoretical calculations and experiment. We apply these new models to 1D hydrodynamic simulations of a polar direct-drive NIF implosion, demonstrating that these new models are now available for future ICF design studies.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Identification of novel organic polar materials: A machine learning study with importance sampling

Recent advances in the synthesis of polar molecular materials have produced practical alternatives to ferroelectric ceramics, opening up exciting new avenues for their incorporation into modern electronic devices. However, in order to realize the full potential of polar polymer and molecular crystals for modern technological applications, it is paramount to assemble and evaluate all the available data for such compounds, identifying descriptors that could be associated with an emergence of ferroelectricity. In this paper, we utilized data-driven approaches to judiciously shortlist candidate materials from a wide chemical space that could possess ferroelectric functionalities. A machine learning study with importance sampling was employed to address the challenge of having a limited amount of available data on already-known organic ferroelectrics. Sets of molecular- and crystal-level descriptors were combined with a Random Forest Regression algorithm in order to predict the spontaneous polarization of the shortlisted compounds. First-principles simulations were performed to further validate the predictions obtained from the machine learning model.

36 MATERIALS SCIENCE↗
DOE OSTI2025

First-Principles insights into group-V impurities and their impact on germanium detector performance

The outstanding properties of high-purity germanium (HPGe) detectors, such as excellent energy resolution, high energy sensitivity, and a low background-to-signal ratio, make them essential and ideal candidates for detecting particle signatures in nuclear processes such as neutrino-less double beta decay (0νββ). However, the presence of defects and impurities in HPGe crystals can lead to charge trapping, which affects carrier mobility and results in significant energy resolution degradation. In this work, we employ density functional theory with a hybrid functional to study the energetics of possible point defects in Ge. Our findings indicate that group-V impurities form more readily in Ge compared to vacancy and interstitial of Ge. Unlike N dopants, which yield deep trap states, P, As, and Sb create shallow traps close to the conduction band edge of Ge. Furthermore, we predict that group-V defects can condense into defect complexes with Ge vacancies. These vacancy-impurity complexes form deep traps in Ge, similar to Ge vacancies, suggesting that both vacancies and vacancy-impurity complexes contribute to charge trapping in these detectors, thereby diminishing their performance.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Interspecies radiative transition in warm and superdense plasma mixtures

Superdense plasmas widely exist in planetary interiors and astrophysical objects such as brown-dwarf cores and white dwarfs. How atoms behave under such extreme-density conditions is not yet well understood, even in single-species plasmas. Here, we have applied thermal density functional theory to investigate the radiation spectra of superdense iron–zinc plasma mixtures at mass densities of ρ = 250 to 2000 g cm -3 and temperatures of kT = 50 to 100 eV, accessible by double-shell–target implosions. Our ab initio calculations reveal two extreme atomic-physics phenomena—firstly, an interspecies radiative transition; and, secondly, the breaking down of the dipole-selection rule for radiative transitions in isolated atoms. Our first-principles calculations predict that for superdense plasma mixtures, both interatomic radiative transitions and dipole-forbidden transitions can become comparable to the normal intra-atomic Kα-emission signal. These physics phenomena were not previously considered for extreme high-density plasma mixtures at super-high energy densities.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Interplay of broken symmetry and delocalized excitations in the insulating state of 1 T –Ta⁢S 2

Coexistence of localized and extended excitations is central to the macroscopic properties of correlated materials. For 5⁢d transition-metal compounds, electron correlations alone generally do not lead to a metal-insulator (Mott) transition, with insulating behavior usually resulting from their coupling with magnetic ordering and/or structural distortions. 1⁢T–Ta⁢S 2 is a prototypical example of such correlated insulating behavior, with a high-symmetry metallic phase transforming into a distorted, charge-density wave (CDW) insulating state at low temperatures. The nature of the insulating phase as well as the existence and relevance of the localized electron physics remains debated. We resolve this standing controversy in 1⁢T–Ta⁢S 2 combining resonant inelastic x-ray spectroscopy and first-principles calculations. We observe five electronic excitations arising from the interband transitions of the Ta 5d orbitals and the S 3⁢p ligand state, with none of the excitations on the order of the Mott gap. These excitations cannot be explained within the framework of standard multiplet calculations that assume a localized wave function, but instead, are captured by a band-theory framework accounting for the low symmetry of the crystal field in the CDW state. Our findings suggest that the electronic properties of 1T–Ta⁢S 2 are dominated in the visible by both plasmonic quasiparticles and interband transitions, with no resonance associated with a putative Mott transition observed in the 0–3 eV energy range. Finally, our discovery provides insights into the electron localization and the Mott vs band-insulator debate in 1T–Ta⁢S 2 and other transition-metal materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

DFT-Based Calculation of Dissolution Activation Energy and Kinetics of Ni–Cr Alloys

A density functional theory investigation of the process of dissolution of Ni and Cr adatoms from model Ni–Cr(111) alloy surfaces is presented, both in vacuum and with explicit water molecules. The goal is to understand how the electronic structures solved using DFT can provide insights as to changes in valence, energy and coordination of the adatoms in the process of dissolving from the alloy surface. It is found that nearby Cr solute atoms increase the dissolution activation energy of Ni. Cr adatoms have a similar dissolution activation energy as Ni adatoms, except for one particular surface configuration that has a much smaller dissolution activation energy, which might promote selective dissolution of Cr in Ni–Cr alloys. To interpret the first-principles modeling results for the potential energy trajectory along the dissolution reaction coordinate, we provide a thermodynamic breakdown of the various terms that contribute to the dissolution activation and total reaction energy. The role of addition of Cr is explored on the dissolution kinetics of a Ni–Cr alloys of varying composition by constructing a model for the corrosion current density with DFT-based activation energies. Lastly, dissolution resistance index is proposed as a quantifiable descriptor of the dissolution resistance of corrosion resistant alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Precipitation of dopants on acceptor-doped LaMnO3±δ revealed by defect chemistry from first principles

Perovskite oxides degrade at elevated temperatures while precipitating dopant-rich particles on the surface. A knowledge-based improvement of surface stability requires a fundamental and quantitative understanding of the dopant precipitation mechanism on these materials. We propose that dopant precipitation is a consequence of the variation of dopant solubility between calcination and operating conditions in solid oxide fuel cells (SOFCs) and electrolyzer cells (SOECs). To study dopant precipitation, we use 20% (D = Ca, Sr, Ba)-doped LaMnO3+δ (LDM20) as a model system. We employ a defect model taking input from density functional theory calculations. The defect model considers the equilibration of LDM20 with a reservoir consisting of dopant oxide (DO), peroxide (DO2), and O2 in the gas phase. The equilibrated non-stoichiometry of the A-site and B-site as a function of temperature, T, and oxygen partial pressure, p(O2), reveals three regimes for LDM20: A-site deficient (oxidizing conditions), A-site rich (atmospheric conditions), and near-stoichiometric (reducing conditions). Assuming an initial A/B non-stoichiometry, we compute the dopant precipitation boundaries in a p–T phase diagram. Our model predicts precipitation both under reducing (DO) and under highly oxidizing conditions (DO2). We found precipitation under anodic, SOEC conditions to be promoted by large dopant size, while under cathodic, SOFC conditions precipitation is promoted by initial A-site excess. The main driving forces for precipitation are oxygen uptake by the condensed phase under oxidizing conditions and oxygen release assisted by B-site vacancies under reducing conditions. Possible strategies for mitigating dopant precipitation under in electrolytic and fuel cell conditions are discussed.

Hess, Franziska (ORCID:000000024460643X)↗
DOE OSTI2024

The Interplay of Solvation and Polarization Effects on Ion Pairing in Nanoconfined Electrolytes

The nature of ion–ion interactions in electrolytes confined to nanoscale pores has important implications for energy storage and separation technologies. However, the physical effects dictating the structure of nanoconfined electrolytes remain debated. Here we employ machine-learning-based molecular dynamics simulations to investigate ion–ion interactions with density functional theory level accuracy in a prototypical confined electrolyte, aqueous NaCl within graphene slit pores. We find that the free energy of ion pairing in highly confined electrolytes deviates substantially from that in bulk solutions, observing a decrease in contact ion pairing but an increase in solvent-separated ion pairing. These changes arise from an interplay of ion solvation effects and graphene’s electronic structure. Notably, the behavior observed from our first-principles-level simulations is not reproduced even qualitatively with the classical force fields conventionally used to model these systems. The insight provided in this work opens new avenues for predicting and controlling the structure of nanoconfined electrolytes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Finite-element-based simulations of electrodes for CO 2 cascade reduction reactions

The multielectron reduction of CO 2 to liquid fuels could be a path to scalable energy storage, but reaching this goal requires major advances in catalysis and systems engineering. Cascade catalysis, which couples sequential reactions without isolating intermediates, has emerged as a promising route to enhance selectivity and efficiency in CO 2 reduction (CO 2 R). In this review, we examine how finite-element-based simulations of continuum model [finite element method (FEM)] approaches are being used to analyze and guide CO 2 R cascade systems. We first outline the fundamentals of cascade catalysis and recent advances in catalytic materials (metallic, molecular, and hybrid architectures). We then focus on FEM developments at the electrode and device scales, emphasizing how these models capture transport phenomena, local microenvironments, and geometry-dependent effects. To clarify design principles, we present case studies of cascade electrodes organized in systems without and with integrated semiconductors. We further emphasize the integration of FEM with multiscale frameworks (density functional theory, molecular dynamics, kinetic Monte Carlo) and its role in bridging atomic-level insights with device-level performance. Finally, we identify current limitations and future prospects, including improved boundary conditions, coupling with operando experiments, and machine learning-accelerated model development. Together, these insights provide design principles for next-generation CO 2 R cascade systems for efficient solar fuel production.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

First Principles Determination of the Potential-of-Zero-Charge in an Alumina-Coated Aluminum/Water Interface Model for Corrosion Applications

The surfaces of most metals immersed in aqueous electrolytes have a several-nanometer-thick oxide/hydroxide surface layer. This gives rise to the existence of both metal|oxide and oxide|liquid electrotlyte interfaces, and makes it challenging to correlate atomic length-scale structures with electrochemical properties such the potential-of-zero-charge (PZC). The PZC has been shown to be correlated the onset potential for pitting corrosion. Here, we conduct large-scale Density Functional Theory and ab initio molecular dynamics to calculate the PZC of a Al(111)| γ -Al 2 O 3 (110)| water double-interface model within the context of aluminum corrosion. By partitioning the multiple interfaces involved into binary components with additive contributions to the overall work function and voltage, we predict the PZC to be -1.53 V vs SHE for this model. Furthermore, we calculate the orbital energy levels of defects like oxygen vacancies in the oxide, which are critical parameters in theories associated with pitting corrosion. We predict that the Fermi level at the PZC lies above the impurity defect levels of the oxygen vacancies, which are therefore uncharged at the PZC. From the PZC estimate, we predict the voltage needed to create oxygen vacancies with net positive charges within a flatband approximation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Octahedral rotations in Ruddlesden-Popper layered oxides under pressure from first principles

The combination of reduced dimensionality and tunable structural distortions in layered perovskite oxides makes these materials ideal platforms for designing novel properties and functionalities. One example is hybrid improper ferroelectricity in n = 2 Ruddlesden-Popper oxides, where the combination of a layered crystal structure and rotations of the metal-oxide octahedra break symmetry and induce a polarization. Precisely controlling the octahedral rotation distortions, for example by the application of hydrostatic pressure, provides a pathway to tune and optimize the properties of these materials. Here we combine group theoretic methods, density functional theory calculations, and Landau theory analysis to investigate how octahedral rotations respond to pressure in the hybrid improper ferroelectrics Sr 3 Zr 2 O 7 , Ca 3 Ti 2 O 7 , and Sr 3 Sn 2 O 7 . We find that factors that are known to control the pressure response of ABO 3 perovskites—the formal charge of the A- and B-site cations, tolerance factor, and B-site chemistry—also impact the pressure response of these layered perovskites. We also show that coupling between the octahedral rotation and strain order parameters plays a key role in determining the overall pressure response. Despite some similarities, we find that these layered perovskites display a distinct pressure response compared to their ABO 3 perovskite analogs. By identifying trends and underlying mechanisms that control octahedral rotations in Ruddlesden-Popper oxides under pressure, this work lays the foundation for tailoring the structure and properties of these materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗