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.

SourceResearch connectionFrom yearThrough year

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

At least 595 records · Page 33

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Search for Stable and Low-Energy Ce–Co–Cu Ternary Compounds Using Machine Learning

Cerium-based intermetallics have garnered significant research attention as potential new permanent magnets. In this study, we explore the compositional and structural landscape of Ce−Co−Cu ternary compounds using a machine learning (ML)- guided framework integrated with first-principles calculations. We employ a crystal graph convolutional neural network (CGCNN), which enables efficient screening for promising candidates, significantly accelerating the material discovery process. With this approach, we predict five stable compounds, Ce 3 Co 3 Cu, CeCoCu 2 , Ce 12 Co 7 Cu, Ce 11 Co 9 Cu, and Ce 10 Co 11 Cu 4 , with formation energies below the convex hull, along with hundreds of low-energy (possibly metastable) Ce−Co−Cu ternary compounds. Firstprinciples calculations reveal that several structures are both energetically and dynamically stable. Notably, two Co-rich low-energy compounds, Ce 4 Co 33 Cu and Ce 4 Co 31 Cu 3 , are predicted to have high magnetizations.

Chemical structure↗

DOE OSTI2025

La 3 CuTe 5 : A Narrow-Gap Semiconductor with Indirect Gap and Dual-Regime Thermally Activated Transport

Metal-chalcogenide systems remain a long-standing research topic because of their structural diversity and potential to host emergent phenomena. Here, we report a new compound, La 3 CuTe 5 , synthesized from the halide-flux method. Single-crystal X-ray diffraction studies indicate the structure to be unique among reported ones. Here, the compound crystallizes in a novel structure type adopting the orthorhombic space group Pnma and a unit cell of a = 24.3947(14) Å, b = 4.4232(2) Å, and c = 10.2142(5) Å. The tetrahedral [CuTe 4 ] building blocks form chains along [010] by corner sharing and link [LaTe 7 ] and [LaTe 8 ] polyhedra via edge sharing, resulting in a three-dimensional bulk structure. Thermal analysis results indicate that the material remains stable with a temperature up to 950 °C and decomposable at 1400 °C. First-principles calculations reveal an indirect electronic band gap and flat valence bands dominated by Te p and Cu d states. Optical absorption measurements yield a band gap of ∼0.65 eV, consistent with semiconducting behavior observed in transport measurements. Fittings to the temperature-dependent resistivity reveal two thermally activated regimes associated with Arrhenius-type conduction and three-dimensional variable range hopping, respectively.

Chalcogenides↗

DOE OSTI2020

Perspectives on van der Waals Density Functionals: The Case of TiS 2

The van der Waals interaction is of foundational importance for a wide variety of physical systems. In particular, van der Waals forces lie at the heart of potential device technologies that may be realized from the functional organization of layered two-dimensional (2D) nanomaterials. Furthermore, for intermediate to large-scale applications modeling, van der Waals density functionals have become the de facto choice for first-principles calculations. In particular, the vdW-DF family of functionals have provided a systematic approach to this theoretically challenging problem. While much progress has been made, there remains room for improvement in the microscopic description of vdW forces from these density functionals. In this work, we compute benchmark results for the binding energy and the electronic density response to binding in TiS2 via accurate diffusion quantum Monte Carlo calculations. We compare these benchmark data to results obtained from local, semilocal, and van der Waals functionals. In particular, we gauge the quality of the original vdW-DF/vdW-DF2 functionals, as well as updated variants such as vdW-DF-C09, vdW-DF-optB88, vdW-DF-optB86b, and vdW-DF2-B86R. We find a close relationship between the accuracy of predicted interlayer separation distances and binding energies for TiS 2 , with the vdW-DF-optB88 functional performing very well in terms of both quantities. In general, the more recently developed functionals are systematic improvements over older ones. However, when considering the response of the electron density to binding, we find that local-density approximation (LDA) and PBEsol generally outperform the vdW-DF functionals in describing the interlayer charge accumulation with vdW-DF-C09 variants performing the best overall.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

High-Pressure Equation of State of 1,3,5-triamino-2,4,6-trinitrobenzene: Insights into the Monoclinic Phase Transition, Hydrogen Bonding, and Anharmonicity

The high-pressure equation of state (EOS) of energetic materials (EMs) is important for continuum and mesoscale models of detonation performance and initiation safety. Additionally, obtaining a high-fidelity EOS of the insensitive EM 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) has proven to be difficult because of challenges in experimental characterization at high pressures (HPs). In this work, powder X-ray diffraction patterns were fitted using the recently discovered monoclinic I2/a phase above 4 GPa, which shows that TATB is less compressible than when indexed with the triclinic $P\bar{1}$ phase. First-principles calculations were performed with Perdew–Burke–Ernzerhof (PBE) and PBE0 functionals including thermal effects using the $P\bar{1}$ phase. PBE0 improves the description of hydrogen bonding and thus predicts accurate planar a and b lattice parameters under ambient conditions. However, discrepancies in the predicted lattice parameters above 4–10 GPa compared with experimental measurements indexed with $P\bar{1}$ are further evidence of a structural modification at high pressure. Layer sliding defects are formed during molecular dynamics simulations, which induces an anharmonic effect on the thermal expansion of the c lattice parameter. In short, the results provide several insights into determining high-fidelity EOS parameters for TATB and other molecular crystals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Electronic Structure and Anion Photoelectron Spectroscopy of Uranium–Gold Clusters UAu n – , n = 3–7

A collaborative effort between experiment and theory towards elucidating the electronic and molecular structures of uranium-gold clusters is presented. Anion photoelectron spectra of UAu n – (n=3-7) were taken at third (355 nm) and fourth (266 nm) harmonics of a Nd:YAG laser, as well as excimer (ArF 193nm) photon energies, where the experimental adiabatic electron affinities (AEA) and vertical detachment energies (VDE) values were measured. Complementary first principles calculations were subsequently carried out to corroborate experimentally determined electron detachment energies and to determine the geometry and electronic structure for each cluster. Except for the ring-like neutral isomer of UAu 6 where one unpaired electron is spread over the Au atoms, all other neutral and anionic UAu n clusters (n = 3 - 7) were calculated to possess open-shell electrons with the unpaired electrons localized on the central U atom. The smaller clusters closely resemble the analogous UF n species, but significant deviations are seen starting with UAu 5 where a competition between U-Au and Au-Au bonding begins to become apparent. The UAu 6 system appears to mark a transition where Au-Au interactions begin to dominate, where both a ring-like and two heavily distorted octahedral structures around the central U atom are to be nearly isoenergetic. With UAu 7 , only ring-like structures are calculated. Altogether, the calculated electron detachment energies are in good agreement with the experimental values

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Ionic Liquids under Radiation and the Dimer Radical Dicyanamide Anion

Ionic liquids (ILs) used in solar, electrochemical, nuclear, or even space exploration applications will necessarily be exposed to conditions in which electron detachment may occur and transient charge-depleted radical species may form. In a recent article (J. Am. Chem. Soc. 2025, 147, 23395–23398), we studied the different possible kinetic fates of an excess electron in the popular 1-butyl-1-methylpyrrolidinium dicyanamide IL but did not address the nature of electron-deficient transient radical species these leaving electrons produce, the so-called holes. Using pulse radiolysis measurements and first-principles calculations in the bulk and gas phases to interpret these, we show that the most likely hole species in this IL is a dimer radical anion that, somewhat surprisingly, has similar spectral features to the excess electron, including a broad near-infrared absorption band.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Methane Activation on Metal-Doped (111) and (100) Ceria Surfaces with Charge-Compensating Oxygen Vacancies

Reducing the temperature for methane activation is an important objective that would benefit many technological applications. In this work we explore the possibility to achieve this goal using single atom catalysts (SACs) on ceria surfaces. In particular, we focus on Ag SACs, which have recently been suggested to be promising catalysts for both the oxidation of molecular hydrogen and CH 4 activation. The rate limiting step of methane activation on ceria is the breaking of the first C-H bond. Using first principles calculations, we investigate this reaction on CeO 2 (111) and (100) surfaces, which are both frequently exposed on ceria nanoparticles. For undoped ceria, our results show that the dissociation of CH 4 to CH 3 + H has a lower barrier on the (001) surface than on the (111) one. On both surfaces, the presence of Ag enhances the binding of CH 4 to the surface and reduces the activation energy for its dissociation. Interestingly, the formation of methanol via the Mars-van Krevelen mechanism has a slightly lower energy barrier than the dissociation of CH 4 to CH 3 + H on Ag-doped CeO 2 (111), suggesting that methanol is a likely product of methane activation on this surface. Here, a novel aspect of our work relative to previous studies on CH 4 activation over transition metal doped ceria is that we mainly focus on stable surface structures where each Ag dopant substituting a Ce atom forms a complex with a surface oxygen vacancy (VO). Our results show that these Ag-VO complexes play a critical role and accounting for their presence is essential for a proper description of the surface reactivity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

RbB 3 Si 3 : An Alkali Metal Borosilicide that is Metastable and Superconducting at 1 atm

The stability, electronic structure and potential superconductivity in AB 3 Si 3 (A = Na, K, Rb, and Cs) compounds that assume the clathrate-based sodalite structure whose frameworks consists of covalent B-Si bonds is investigated via first-principles calculations. This structure type has recently been predicted in a number of high temperature superconducting hydrides, but these are only stable under megabar pressures. Herein, we predict a novel superconducting phase, RbB 3 Si 3 , that could be synthesized under pressures that are a factor of ten smaller, ~10 GPa, and quenched to atmospheric conditions. Electron-phonon coupling calculations predict that RbB 3 Si 3 possesses a superconducting critical temperature, T c , of 14 K at 1 atm. Here, the dynamic stability of RbB 3 Si 3 and CsB 3 Si 3 at ambient pressure can be explained by considering the chemical pressure exerted on the B-Si framework that is caused by the size effect of the alkali metal atom.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Intercalation of Lithium into Graphite: Insights from First-Principles Simulations

Understanding ion intercalation at electrode–electrolyte interfaces is key to the development of energy storage and water desalination. Here, we investigate Li + kinetics at a prototypical interface between graphite anodes and an organic electrolyte, and we elucidate key factors that determine ion transport, using first-principles methodology coupling ab initio molecular dynamics simulations with a solvation model. We show that surface chemical composition significantly influences the kinetics of ion intercalation from the liquid into graphite. We find that this is partly related to the ion desolvation process, which varies notably for different graphite surface chemical terminations. In addition, interfacial polarization is found to play an important role in determining energy barriers for ion transfer. We also discuss the impact of electrode potentials, which is often neglected in conventional first-principles calculations despite being a key factor in device configurations. Our study provides insights into the coupling of electronic and ionic effects of interfacial chemistry on ion transport at complex electrode–electrolyte interfaces.

25 ENERGY STORAGE↗

DOE OSTI2020

Electronic Structure of Double-Layer Epitaxial Graphene on SiC(0001) Modified by Gd Intercalation

In this work, we systematically study the effects of Gd adsorption and intercalation on the electronic band structure of double-layer epitaxial graphene on Si-terminated SiC(0001) by first-principles calculations. We show that Gd adsorption and intercalation exhibit strong effects on the coupling between the graphene layers and between the buffer layer and substrate. Different adsorption/intercalation geometries can result in very different electron band structures. The number of Dirac cones and the positions of the Dirac cones relative to the Fermi level can be effectively manipulated through controlling the Gd adsorption/intercalation geometries. Our calculations provide useful insights to guide the experimental design of graphene-based materials with desirable functionalities for applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Mechanism of Metal Intercalation under Graphene through Small Vacancy Defects

Metal intercalation under graphene has attracted extensive experimental and theoretical research because of its capability to manipulate the electronic structure and properties of graphene. However, the pathways and mechanisms of intercalation are still not well understood. Here, we systematically investigate the intercalation process of metal atoms through graphene vacancies using first-principles calculations. We show that the energy barrier for metal atom penetration through the vacancies in graphene is small as long as the size of the vacancy is larger than a mono-vacancy. However, metal atoms are strongly bonded to the vacancy so that the detachment energy of a metal atom from the vacancy is extremely high. This inhibits the diffusion of the metal atom into the gallery beneath the surface to complete the intercalation process. On the other hand, our calculation results show that the detachment energy of a metal atom from a metal dimer at small vacancy defects is significantly reduced, making intercalation much easier. Therefore, the key step limiting the intercalation process is the detachment of the metal atoms from vacancy defects. This finding from our study provides useful insight into the defect-assisted intercalation mechanism.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Correlation between Spin and Orbital Dynamics during Laser-Induced Femtosecond Demagnetization

Spin and orbital angular momenta are two intrinsic properties of an electron and are responsible for the physics of a solid. How the spin and orbital evolve with respect to each other on several hundred femtoseconds is largely unknown, but it is at the center of laser-induced ultrafast demagnetization. In this paper, we introduce a concept of the spin–orbital correlation diagram, where spin angular momentum is plotted against orbital angular momentum, much like the position-velocity phase diagram in classical mechanics. We use four sets of highly accurate time-resolved X-ray magnetic circular dichroism data to construct four correlation diagrams for iron and cobalt. To our surprise, a pattern emerges. The trace on the correlation diagram for iron is an arc, and at the end of demagnetization, it has a pronounced cusp. The correlation diagram for cobalt is different and appears more linear but with kinks. We carry out first-principles calculations with two different methods: time-dependent density functional theory (TDDFT) and time-dependent Liouville density functional theory. These two methods agree that the experimental findings for both Fe and Co are not due to experimental errors. It is the spin–orbit coupling that correlates the spin dynamics to the orbital dynamics. Microscopically, Fe and Co have different orbital occupations, which leads to distinctive correlation diagrams. We believe that this correlation diagram presents a useful tool to better understand spin and orbital dynamics on an ultrafast time scale. A brief discussion on the magnetic anisotropy energy is also provided.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Computational Investigation of the Catalytic Hydrodeoxygenation of Propanoic Acid over a Cu(111) Surface

Cu-based alloy catalysts have recently been investigated experimentally for the hydrodeoxygenation (HDO) of biomass-derived organic acids. Here, the HDO of propanoic acid (PAc) has been studied over Cu(111) by mean-field microkinetic modeling based on parameters obtained from first-principles calculations. Models were developed for the gas- and liquid-phase HDO in condensed water and 1,4-dioxane. In agreement with experimental observations, the gas-phase PAc conversion rate is low at 573 K and increases in liquid water by 1 order of magnitude. In all reaction environments, the decarboxylation mechanism is dominant at low hydrogen partial pressures less than 0.1 bar, and the C–COO bond dissociation is the rate-controlling elementary step. This observation contrasts with the rate-controlling step identified over most group VIII metal surfaces, which is the C–OH bond dissociation in the decarbonylation mechanism. At high hydrogen (H 2 ) partial pressures greater than 10 bar, the HDO of PAc produces propionaldehyde that can readsorb and further react through decarbonylation to produce C 2 alkane products, which is conceptually different from the low H 2 partial pressure scenario. At high H 2 partial pressures, the initial hydrogenation at the carbonyl carbon of PAc becomes the rate-controlling elementary step.

Catalysts↗

DOE OSTI2021

Quantifying Temperature Dependence of Electronic Band Gaps and Optical Properties in SnO 2 and SnO via First-Principles Simulations

Tin metal oxides SnO x ( x = 1, 2) have gained interest as gas-sensing materials. For their applications as high-temperature sensors, however, a better understanding of their temperature dependence sensing responses is needed. Here, we comparatively quantify the temperature-dependent electronic band gaps and optical properties of SnO 2 and SnO using first-principles calculations. Without considering the temperature effect, SnO 2 and SnO are predicted to have direct and indirect band gaps of 2.18 and 1.66 eV, respectively, at the PBE + U -GGA level that we employed. The temperature effect on the electronic and optical properties is captured by taking account of the electron–phonon interaction. Band gap renormalization with temperature is calculated via the Allen–Heine–Cardona theory. For both oxides, we find a monotonic decrease in the electronic band gap such that renormalization at zero point (0 K) is ~–0.17 and ~–0.52 eV at 1000 K. These results are also analyzed by employing an analytical equation that helps characterize the band gap shift with temperature. In addition, the optical properties at finite temperatures are simulated using the frozen-phonon method that combines electron–phonon coupling with the momentum matrix. As temperature increases, the optical property spectra are smoothed because of the smearing effect, which diminishes optical constants at shorter wavelengths. Our results are of interest for high-temperature functional materials in applications of optical detection.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Straintronic Effect on Phonon-Mediated Superconductivity of Nb 2 CT 2 (T = O, S, Se, or Te) MXenes

Here, the electronic structures, phonon dispersions, and electron–phonon coupling of Nb 2 CT 2 (T = O, S, Se, or Te) MXenes were investigated via first-principles calculations. Different models of Nb 2 CT 2 were constructed, and the results show that the low-energy models of Nb 2 CT 2 are intrinsic phonon-mediated superconductors. Of the four Nb 2 CT 2 MXenes, Nb 2 CO 2 MXene exhibits the largest superconducting critical temperature ($T_c$) of 14.43 K. The existence of soft modes induced by Kohn anomalies and the contribution of Nb atoms to the Fermi level lead to strong electron–phonon coupling (λ = 0.92) in Nb 2 CO 2 MXene. The $T_c$ of Nb 2 CO 2 is further enhanced by biaxial tensile strain and reaches up to 18.28 K under 4% tensile strain. The predicted $T_c$ of Nb 2 CS 2 is 4.5 K, which is comparable with experimental data. These findings will further stimulate the search for superconducting MXenes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Effect of Multiple Oxygen Vacancies on the Optical and Thermodynamic Properties of La 0.75 Sr 0.25 Co 0.25 Fe 0.75 O 3-δ Perovskite

We report the oxygen vacancy concentration of perovskite changes with the work condition in applications such as oxygen carriers and oxygen sensors. To explore the effect of multiple oxygen vacancies on the optical and thermodynamic properties, we perform first-principles calculations to investigate the formation of multiple oxygen vacancies in the paramagnetic La 0.75 Sr 0.25 Co 0.25 Fe 0.75 O 3-δ . Our results show that the higher VO concentration in La 0.75 Sr 0.25 Co 0.25 Fe 0.75 O 3-δ occurs at the conditions of higher temperatures and lower oxygen partial pressures. The electrons liberated upon oxygen vacancy formation are partially localized on the nearby Fe/Co sites and partially delocalized through the lattice, with the result that the unoccupied Fe/Co states near the Fermi level decrease and the Fermi level shows an upward shift. From the results of optical absorption, the peaks from 250 to 600 nm are weakened and blue-shifted with increase of oxygen vacancy concentration, due to the decrease of empty Fe states and an upward shift of the Fermi level. The obvious change of optical absorption with oxygen vacancy concentration indicates that La 0.75 Sr 0.25 Co 0.25 Fe 0.75 O 3-δ could be a good candidate for optical oxygen sensors in the visible spectrum.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Unraveling 1-Hexene Hydrogenation over Dilute Pd-in-Au Alloys

Dilute Pd-in-Au alloys are valuable materials for selectively hydrogenating alkynes and isomerizing alkenes. By diluting Pd in a Au host, the selectivity toward semihydrogenated alkene isomers can be significantly enhanced and the unfavorable overhydrogenation to alkanes is suppressed. However, a detailed mechanistic study on the origin of the enhanced alkene selectivity over dilute alloy catalysts is still missing. Here, we combine experiment and theory to unravel the reaction mechanism, identifying rate-limiting and selectivity-controlling steps in 1-hexene hydrogenation over dilute Pd-in-Au catalysts. Using isotope-exchange hydrogenation experiments, we show that 1-hexene and hydrogen over a bimetallic Pd4Au96 in silica catalyst preferentially form 1-hexene isomers, (trans and cis) 2- and 3-hexene and only a small amounts of hexane. The reaction is consistent with a Horiuti–Polanyi mechanism, similar to a monometallic Pd nanoparticle catalyst. Computation of the free-energy profiles for 1-hexene hydrogenation and isomerization over a single Pd atom in a Au surface using first principles calculations indicated that the isomerization of 1-hexene to 2-hexene is energetically favorable due to the relatively large barrier for H2 dissociation preventing hydrogenation to n-hexane. Microkinetic modeling established that H2 dissociation on the single-atom Pd sites and H spillover from these sites onto the Au host are rate-limiting and key in steering the selectivity of dilute Pd-in-Au alloys toward the hexene isomers. Furthermore, the mechanistic insights from this study contribute to the rational design of optimized dilute alloy catalysts for selective alkene isomerization.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Computational Screening of Cationic Zeolites for n -Butane/Methane Separations Using Quantitatively Accurate First-Principles-Derived Force Fields

We developed force fields for linear alkanes in Na and Ca-exchanged zeolites based on periodic first-principles calculations with coupled cluster corrections. These force fields were validated by comparing the simulated adsorption properties with extensive experimental data. As an example of using these force fields, we screened silica and aluminosilicate zeolites for separation of n-butane/methane mixtures in several different separation processes. A key step in this screening was the use of an interpolation method to create accurate zeolite structures that take account of the lattice constant changes for different Si/Al ratios for a given zeolite topology. We focused on zeolite structures that have known synthesis routes, so the materials selected from our calculations can be tested experimentally. Several promising materials that show good separation performance were selected for more detailed simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗