Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Density functional calculations”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Metal Identity and Coordination Environment Modulate Single-Atom Catalyst Stability During Electrocatalysis

A major hurdle to the implementation of single-atom catalysts (SACs) in real-world systems is a poor understanding of their stability under operating conditions, which is particularly relevant due to the high surface free energy of SACs. Here, we evaluated the aggregation behavior of a suite of SACs varied by metal identity (Fe, Co, Ni, and Cu) during electrocatalytic nitrate reduction using in situ X-ray absorption spectroscopy. The metal center had significant influence on reconstruction, where under identical applied reductive potentials, SACs underwent varying levels of reconstruction, ranging from no discernible change to complete reduction into metallic nanoparticles. Such in situ experiments revealed Cu SACs to be the most susceptible to aggregation, prompting a deeper investigation into how coordination environment (O-, B-, and N-graphene) affected Cu SAC aggregation. We further conducted density functional theory calculations to elucidate the relationship between Cu SAC structure and stability. This work deconvolutes the relationship between SAC architecture and stability, which is essential to evaluate and explain for the realization of SACs for electrocatalysis.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Targeted Chemical Looping Materials Discovery by an Inverse Design

Chemical looping with oxygen uncoupling (CLOU) materials is actively sought for combustion of carbonaceous materials to achieve complete conversion and capture of carbon dioxide. These materials may play a vital role in reducing atmospheric carbon via negative carbon output. However, there is no one‐size‐fits‐all approach as different operating conditions and feedstocks may require different CLOU materials. As a result, the exploration and discovery of high‐performance CLOU materials can be a slow process. To address this challenge, a high‐throughput inverse machine learning workflow that identifies optimum materials from perovskite oxides for a given set of targets is developed—temperature and Gibbs free energy of oxygen formation. The model is trained on high‐throughput density functional theory calculations of CLOU materials and inverts the materials design process using a genetic algorithm to produce realistic substituted SrFeO 3‐δ compositions as output. Using the inverse model, it is able to identify several interesting new families of CLOU materials: Sr 1‐ x A x Fe 1‐ y B y O 3‐δ (e.g., A = Ca or K; B = Mg, Bi, Mn, Ni, Co, Cu, or Zn). These materials have shown promising properties, and some of them even outperform the benchmark material in terms of oxygen release kinetics under relevant CLOU operating conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Heavy boron doping effects on biaxially tensile strained germanium (>1.5%) investigated via structural characterization, effective lifetime assessment and atomistic modeling

Highly tensile strained germanium (ε-Ge) represents an essential material system for emerging electronic and photonics applications. Moreover, adjusting the doping levels to moderate or high concentrations can effectively tailor the properties of ε-Ge for specific applications. This article combines experimental characterization with a theoretical framework to examine the effects of heavy elemental boron (B) doping on pseudomorphic sub-50 nm ε-Ge. High resolution X-ray diffractometry is used to validate tensile strain levels of 1.53% and 1.68% in Ge epilayers, surpassing the indirect-to-direct band gap crossover point at ∼1.5% biaxial tensile strain. Cross-sectional transmission electron microscopy revealed visual evidence of stacking faults and surface roughening in 1.68% ε-Ge, although a coherent and abrupt Ge/III–V heterointerface is observed, devoid of interfacial misfit dislocations. Effective lifetime measurements demonstrated approximately twofold enhancement in 1.53% B-doped ε-Ge (N B ∼7 × 10 19 cm −3 ) compared to its unstrained B-doped counterpart, while no such improvement was observed in 1.68% B-doped ε-Ge. This lack of enhancement is attributed to the presence of stacking faults and surface roughness within the ε-Ge epilayer. Through density functional theory calculations, we independently demonstrate that substitutional B atoms induce local deformation of Ge–Ge bonds in both unstrained Ge and ε-Ge epilayers, resulting in an additive tensile strain. This phenomenon could potentially lead to dynamic reduction and overcoming of the critical layer thickness for the system, facilitating the nucleation and subsequent glide of 90° leading Shockley partial dislocations, thereby generating stacking faults. In essence, these findings establish an upper limit on the B-doping concentration that can be achieved in highly ε-Ge epilayers, and collectively, offer valuable insights into the significance of heavy doping in Ge-based heterostructures. As such, this study delineates a fundamental constraint for integrating heavily doped ε-Ge in high-performance optoelectronic systems, necessitating precise strain-doping co-optimization to avoid performance degradation.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2024

Ultrafast Optically Induced Perturbation of Oxygen Octahedral Rotations in Multiferroic BiFeO 3 Thin Films

The functional properties of complex oxides, including magnetism and ferroelectricity, are closely linked to subtle structural distortions. Ultrafast optical excitations provide the means to manipulate structural features and ultimately to affect the functional properties of complex oxides with picosecond-scale precision. Here, we report that the lattice expansion of multiferroic BiFeO 3 following above-bandgap optical excitation leads to distortion of the oxygen octahedral rotation (OOR) pattern. The continuous coupling between OOR and strain was probed using time-resolved X-ray free-electron laser diffraction with femtosecond time resolution. Density functional theory calculations predict a relationship between the OOR and the elastic strain consistent with the experiment, demonstrating a route to employing this approach in a wider range of systems. Ultrafast control of the functional properties of BiFeO 3 thin films is enabled by this approach because the OOR phenomena are related to ferroelectricity, and via the Fe-O-Fe bond angles, the superexchange interaction between Fe atoms.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Machine learned synthesizability predictions aided by density functional theory

Abstract A grand challenge of materials science is predicting synthesis pathways for novel compounds. Data-driven approaches have made significant progress in predicting a compound’s synthesizability; however, some recent attempts ignore phase stability information. Here, we combine thermodynamic stability calculated using density functional theory with composition-based features to train a machine learning model that predicts a material’s synthesizability. Our model predicts the synthesizability of ternary 1:1:1 compositions in the half-Heusler structure, achieving a cross-validated precision of 0.82 and recall of 0.82. Our model shows improvement in predicting non-half-Heuslers compared to a previous study’s model, and identifies 121 synthesizable candidates out of 4141 unreported ternary compositions. More notably, 39 stable compositions are predicted unsynthesizable while 62 unstable compositions are predicted synthesizable; these findings otherwise cannot be made using density functional theory stability alone. This study presents a new approach for accurately predicting synthesizability, and identifies new half-Heuslers for experimental synthesis.

Lee, Andrew (ORCID:0000000153014295)↗
DOE OSTI2021

Synthesis, structure, linear and nonlinear optical properties of noncentrosymmetric quaternary diamond-like semiconductors, Cu 2 ZnGeSe 4 (CZGSe) and the novel Cu 4 ZnGe 2 Se 7

In this work, two quaternary diamond-like semiconductors (DLSs) in the Cu-Zn-Ge-Se system, the known Cu 2 ZnGeSe 4 (also referred to as CZGSe) and the new Cu 4 ZnGe 2 Se 7 , are compared in terms of their crystal structures, electronic structures and physicochemical properties. Both compounds were prepared by high-temperature, solid-state synthesis at 800 degrees C. Single crystal X-ray diffraction was used to determine the structure of Cu 4 ZnGe 2 Se 7 . The structures of both Cu 2 ZnGeSe 4 and Cu 4 ZnGe 2 Se 7 can be considered as derivatives of cubic diamond. Cu 4 ZnGe 2 Se 7 , with a reduced symmetry due to the cation ordering pattern, adopts the Cu 4 NiSi 2 S 7 structure type, with space group C2. The corner-sharing tetrahedra in Cu 4 ZnGe 2 Se 7 are slightly distorted because the charge for some of the S 2 - anions is not compensated by the first-nearest-neighbor cations. Rietveld refinements using synchrotron X-ray powder diffraction data were used to assess the phase purity of the samples and confirm the bulk structural behavior. Diffuse reflectance UV/Vis/NIR spectroscopy shows that Cu 2 ZnGeSe 4 and Cu 4 ZnGe 2 Se 7 have direct optical bandgaps of 1.38 and 0.91 eV, respectively. Electronic structure calculations implementing density functional theory confirm the direct bandgap for Cu 4 ZnGe 2 Se 7 , with a calculated value of 0.62 eV. Both compounds are air stable, thermally stable up to relatively high temperatures, undergo phase transitions and have wide windows of optical clarity. The Kurtz-Perry powder technique was used to determine the second harmonic generation (SHG) responses using a commercial AgGaSe 2 standard. While Cu 4 ZnGe 2 Se 7 displays a weak SHG response, Cu 2 ZnGeSe 4 exhibits a response that is greater than several benchmark materials, with a large SHG coefficient, χ (2) , of 43 ± 6 pm/V at λ = 2900 nm.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Metal–Organic Framework Integrating Ionic Framework and Bimetallic Coupling Effect for Highly Efficient Oxygen Evolution Reaction

Metal–organic frameworks (MOFs) are recognized as promising electrocatalysts for the oxygen evolution reaction (OER) because of their permanent porosity and rich architectural diversity; however, ionic MOFs enabling fast ions exchange during OER are rarely explored. Here, an ionic MOF (Ni-btz) constructed with an azolate ligand is selected, and continuous 3D bimetallic MOF (NiFe-btz) films deriving from high-degree intergrowth of microsized MOFs particles are fabricated. The as-prepared NiFe-btz/NF-OH electrode exhibits excellent OER performance with a low overpotential of 239 mV at 10 mA cm –2 under alkaline condition. The OER charge transfer process and bimetallic coupling effect in ionic NiFe-btz are probed by density functional theory calculations and confirmed via X-ray photoelectron spectroscopy and in situ Raman measurements. The partial density of states of NiFe-btz indicates that the main contribution for electron density around the Fermi level is from Cl ions clarifying the profitable impact of ionic MOF framework. This work systematically demonstrates the relationship of electronic structure and OER activity in ionic, bimetallic MOFs and expands the scope of 3D MOF films for efficient OER.

36 MATERIALS SCIENCE↗
DOE OSTI2021

The ultrahigh pressure stability of silver: An experimental and theoretical study

We measured the atomic volume of Ag in a toroidal diamond anvil cell to a maximum pressure of 416 GPa and calculated the atomic volume and elastic constants of Ag up to 750 and 460 GPa, respectively. Our density functional theory calculations at 0 K utilize an all-electron fully relativistic method and agree well with our volume measurements, particularly at pressures above ~75 GPa. In this work, we corrected our experimental results for non-hydrostaticity using a line shift analysis, and the resulting Vinet equation of state (EOS) parameters are reported. We find that the uniaxial stress sustained by Ag increases linearly up to 4.5 GPa at a pressure of 416 GPa. Our experimental results indicate that the fcc structure of Ag remains stable to at least 416 GPa at room temperature. Our theoretical results show that C44 increases as pressure increases and reaches a maximum at ~100 GPa above which it begins to decrease, a sign that the fcc structure of Ag is becoming unstable, and at V/V 0 = 0.30, the bcc structure is lower in energy than fcc.

36 MATERIALS SCIENCE↗
NASA NTRS1994

An Ab Initio Study of CuCO

Modified coupled-pair functional (MCPF) calculations and coupled cluster singles and doubles calculations, which include a perturbational estimate of the connected triples [CCSD(T)], yield a bent structure for CuCO, thus, supporting the prediction of a nonlinear structure based on density functional (DF) calculations. Our best estimate for the binding energy is 4.9 +/- 1.4 kcal/mol; this is in better agreement with experiment (6.0 +/- 1.2 kcal/mol) than the DF approach which yields a value (19.6 kcal/mol) significantly larger than experiment.

Bauschlicher, Charles W., Jr.↗
DOE OSTI2024

Digital twins and deep learning segmentation of defects in monolayer MX 2 phases

Developing methods to understand and control defect formation in nanomaterials offers a promising route for materials discovery. Monolayer MX 2 phases represent a particularly compelling case for defect engineering of nanomaterials due to the large variability in their physical properties as different defects are introduced into their structure. However, effective identification and quantification of defects remain a challenge even as high-throughput scanning transmission electron microscopy methods improve. This study highlights the benefits of employing first principles calculations to produce digital twins for training deep learning segmentation models for defect identification in monolayer MX 2 phases. Around 600 defect structures were obtained using density functional theory calculations, with each monolayer MX 2 structure being subjected to multislice simulations for the purpose of generating the digital twins. Several deep learning segmentation architectures were trained on this dataset, and their performances evaluated under a variety of conditions such as recognizing defects in the presence of unidentified impurities, beam damage, grain boundaries, and with reduced image quality from low electron doses. Further, this digital twin approach allows benchmarking different deep learning architectures on a theory dataset, which enables the study of defect classification under a broad array of finely controlled conditions. It thus opens the door to resolving the underpinning physical reasons for model shortcomings and potentially chart paths forward for automated discovery of materials defect phases in experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Cubic Crystal Structure Formation and Optical Properties within the Ag–B II –M IV –X (B II = Sr, Pb; M IV = Si, Ge, Sn; X = S, Se) Family of Semiconductors

Quaternary chalcogenide semiconductors are promising materials for energy conversion and nonlinear optical applications, with properties tunable primarily by varying the elemental composition and crystal structure. Here, we first analyze the connections among several cubic crystal structure types, as well as the orthorhombic Ag 2 PbGeS 4 -type structure, reported for select members within the Ag–B II –M IV –X (B II = Sr, Pb; M IV = Si, Ge, Sn; X = S, Se) compositional space. Focusing on the Ag–Pb–Si–S and Ag–Sr–Sn–S systems, we show that one structure type, with the formulas Ag 2 Pb 3 Si 2 S 8 and Ag 2 Sr 3 Sn 2 S 8 , is favored. We have prepared powder and single-crystal samples of Ag 2 Pb 3 Si 2 S 8 and Ag 2 Sr 3 Sn 2 S 8 , showing that each takes on the noncentrosymmetric cubic space group I$\bar{43}$d and is isostructural to the previously reported compound Ag 2 Sr 3 Ge 2 Se 8 . Through hybrid density functional theory calculations, these cubic compounds are demonstrated to be (quasi-)direct band gap semiconductors with high densities of states at the band maxima. Furthermore, the band-gap energies are measured by reflectance spectroscopy as 1.95(3) and 2.66(4) eV for Ag 2 Pb 3 Si 2 S 8 and Ag 2 Sr 3 Sn 2 S 8 , respectively. We further measure the optical properties and show the electronic band structures of three other isostructural A I –B II –M IV –X-type materials, i.e., Ag 2 Sr 3 Si 2 S 8 , Ag 2 Sr 3 Ge 2 S 8 , and Ag 2 Sr 3 Ge 2 Se 8 , showing that the band gaps can be predictably tuned by element substitution. Detailed visual analyses of the different structures and of their relationships with other members of the Ag–B II –M IV –X compositional family provide a basis for a broader understanding of the structure formation and optoelectronic properties within the quaternary chalcogenide semiconductor family.

14 SOLAR ENERGY↗
DOE OSTI2023

Double-Layer Kagome Metals Pt 3 Tl 2 and Pt 3 In 2

The connectivity and inherent frustration of the kagome lattice can produce interesting electronic structures and behaviors in compounds containing this structural motif. Here we report the properties of Pt 3 X 2 (X = In and Tl) that adopt a double-layer kagome net structure related to that of the topologically nontrivial high-temperature ferromagnet Fe 3 Sn 2 and the density wave hosting compound V 3 Sb 2 . We examined the structural and physical properties of single crystal Pt 3 Tl 2 and polycrystalline Pt 3 In 2 using X-ray and neutron diffraction, magnetic susceptibility, heat capacity, and electrical transport measurements, along with density functional theory calculations of the electronic structure. Our calculations show that Fermi levels lie in pseudogaps in the densities of states with several bands contributing to transport, and this is consistent with our Hall effect, magnetic susceptibility, and heat capacity measurements. Although electronic dispersions, characteristic of simple kagome nets with nearest-neighbor hopping, are not clearly seen, likely due to the extended nature of the Pt 5d states, we do observe moderately large and non-saturating magnetoresistance values and quantum oscillations in the magnetoresistance and magnetization associated with the kagome nets of Pt.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Benchmarking Density Functional Theory Methods for Efficient Calculations of a Strongly Correlated Li 1– x Ni 1– y O 2−δ System

Transition metal oxides (TMOs), such as LiNiO 2 , are promising candidates for energy storage and electronic devices due to their unique electronic properties, exceptional physical and chemical characteristics, and ability to adopt multiple oxidation states. However, accurately predicting their properties using mean-field density functional theory (DFT) is challenging due to the presence of strongly correlated d-electrons and the complex interplay between their structural, electronic, and magnetic responses. These challenges are further exacerbated by the need to model defects, surfaces, and interfaces, which require computationally efficient, large-scale simulations. To address these issues, we carry out a benchmark study on the Li 1–x NiO 2 system, evaluating the performance of several popular functionals. Our findings demonstrate that combining SCAN functional relaxation with single-step HSE calculations provides a practical and scalable computational strategy. This approach balances accuracy and efficiency, enabling high-throughput simulations of strongly correlated TMOs and improved predictive modeling capability of TMOs for practical applications.

25 ENERGY STORAGE↗
DOE OSTI2025

Defluorination Mechanisms and Real-Time Dynamics of Per- and Polyfluoroalkyl Substances on Electrified Surfaces

Per- and polyfluoroalkyl substances (PFAS) are persistent environmental contaminants found in groundwater sources and a wide variety of consumer products. In recent years, electrochemical approaches for the degradation of these harmful contaminants have garnered a significant amount of attention due to their efficiency and chemical-free modular nature. However, these electrochemical processes occur in open, highly non-equilibrium systems, and a detailed understanding of PFAS degradation mechanisms in these promising technologies is still in its infancy. To shed mechanistic insight into these complex processes, we present the first constant-electrode potential (CEP) quantum calculations of PFAS degradation on electrified surfaces. These advanced CEP calculations provide new mechanistic details about the intricate electronic processes that occur during PFAS degradation in the presence of an electrochemical bias, which cannot be gleaned from conventional density functional theory calculations. We complement our CEP calculations with large-scale ab initio molecular dynamics simulations in the presence of an electrochemical bias to provide time scales for PFAS degradation on electrified surfaces. Taken together, our CEP-based quantum calculations provide critical reaction mechanisms for PFAS degradation in open electrochemical systems, which can be used to prescreen candidate material surfaces and optimal electrochemical conditions for remediating PFAS and other environmental contaminants.

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

Defect levels and self-compensation in iodine-doped CdTe single crystals

This study systematically documents defect levels in n-type iodine-doped cadmium telluride (CdTe:I) crystals as measured by thermoelectric effect spectroscopy, Hall-effect measurements, and photoluminescence and calculated by density functional theory. The primary donor, I Te , was identified with an activation energy of ∼0.05 eV (measured)/∼0.13 eV (calculated). Deep acceptor states, V Cd and I Te –V Cd complex, exhibited activation energies of ∼0.1 eV (measured)/∼0.19 eV (calculated) and ∼0.12 eV (measured)/∼0.24 eV (calculated), respectively. Self-compensation and Fermi level pinning were observed at ∼0.7 eV (measured)/∼1.08 eV (calculated), contributing to high resistivity in as-grown samples. Post-growth Cd annealing effectively removes compensation centers, releases iodine donors, and significantly reduces resistivity, resulting in an ideal n-type solar cell material.

Cd annealing↗
DOE OSTI2020

Simulated field-modulated x-ray absorption in titania

Here, we present a method to compute the x-ray absorption near-edge structure (XANES) spectra of solid-state transition metal oxides using real-time time-dependent density functional theory, including spin–orbit coupling effects. This was performed on bulk-mimicking anatase titania (TiO 2 ) clusters, which allows for the use of hybrid functionals and atom-centered all electron basis sets. Furthermore, this method was employed to calculate the shifts in the XANES spectra of the Ti L-edge in the presence of applied electric fields to understand how external fields can modify the electronic structure, and how this can be probed using x-ray absorption spectroscopy. Specifically, the onset of t 2g peaks in the Ti L-edge was observed to red shift and the e g peaks were observed to blue shift with increasing fields, attributed to changes in the hybridization of the conduction band (3d) orbitals.

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

Antimony stable isotope fractionation during adsorption onto birnessite: A molecular perspective from X-ray absorption spectroscopy and density functional theory

Sorption of antimony (Sb) onto birnessite significantly influences the fate of Sb in oceanic and terrestrial environments and fractionates Sb isotopes. Nevertheless, little is known about Sb isotopic fractionation during its adsorption on birnessite. Here, in this study, we show the value of Δ 123 Sb adsorbed-aqueous increases from −0.398 to −0.332 ‰ in 1 h and then decreases and stabilizes at −0.384 ‰ in 72 h. The enrichment of the light Sb isotope is predominantly due to the distortion of the octahedral symmetry. X-ray absorption spectroscopy results indicate Sb first forms a double-corner-sharing complex on birnessite and then transforms to a double-edge-sharing complex during adsorption. The optimized bond distances for double-corner-sharing (3.37 Å) and double-edge-sharing (2.90 Å) complexes calculated using density functional theory (DFT) fits well with the structure (3.41 and 3.00 Å) revealed by X-ray absorption spectroscopy, respectively. The fractionation derived from reduced partition function ratios calculated using DFT aligns well with the experimental results. Therefore, the variation in Sb isotopic fractionation during adsorption is attributed to the evolving structure of Sb complexes on birnessite. Our results demonstrate the isotopic fractionation of Sb during adsorption on birnessite and provide a molecular-scale understanding of Sb behavior, contributing to the correct reconstruction of the Sb isotope composition of ancient seawater using ferromanganese crusts and nodules, and efforts to trace Sb migration in epigenetic mining environments.

Adsorption↗
DOE OSTI2021

P–V–T Equation of State of Iridium Up to 80 GPa and 3100 K

In the present study, the high-pressure high-temperature equation of the state of iridium has been determined through a combination of in situ synchrotron X-ray diffraction experiments using laser-heating diamond-anvil cells (up to 48 GPa and 3100 K) and density-functional theory calculations (up to 80 GPa and 3000 K). The melting temperature of iridium at 40 GPa was also determined experimentally as being 4260 (200) K. The results obtained with the two different methods are fully consistent and agree with previous thermal expansion studies performed at ambient pressure. The resulting thermal equation of state can be described using a third-order Birch–Murnaghan formalism with a Berman thermal-expansion model. The present equation of the state of iridium can be used as a reliable primary pressure standard for static experiments up to 80 GPa and 3100 K. A comparison with gold, copper, platinum, niobium, rhenium, tantalum, and osmium is also presented. On top of that, the radial-distribution function of liquid iridium has been determined from experiments and calculations.

36 MATERIALS SCIENCE↗