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

Fundamental Thermodynamic, Kinetic, and Mechanical Properties of Lithium and Its Alloys

Lithium alloying reactions are beneficial in promoting uniform plating and stripping of lithium metal in all-solid-state batteries. First-principles calculations are performed to predict thermodynamic, kinetic, and mechanical properties of lithium and several important Li–M alloys (M = Mg, Ag, Zn, Al, Ga, In, Sn, Sb, and Bi). While the Li–Mg binary system forms a solid solution, most other lithium–metal alloys prefer stoichiometric intermetallic compounds with common local motifs that enable fast Li diffusion. Lithium and Li-rich alloys exhibit an unusually flat energy landscape along paths that connect BCC to close-packed structures like FCC and HCP, with important implications for mechanical properties. Very low migration barriers for Li diffusion that rival those of superion conductors are predicted, both in pure Li and in Li–M intermetallics. However, vacancy concentration, which is crucial for substitutional diffusion, is predicted to be low in metallic Li and most Li–M intermetallics. Compounds such as B32 LiAl and LiGa as well as D0 3 Li 3 Sb and Li 3 Bi exhibit structural vacancies at higher ends of their voltage windows, which together with low migration barriers leads to exceptionally high Li mobilities. In the Li–Mg solid solution, the addition of Mg is found to decrease the vacancy tracer diffusion coefficient by an order of magnitude.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Ligand Hole Driven Metal–Insulator Transition Exemplified in a Layered Transition Metal Oxide

The interplay of cooperative Jahn–Teller (JT) distortions and charge-disproportionation (CD) with a strong electronic correlation in transition metal oxides leads to structural symmetry breaking. Both JT and CD often manifest in the form of significant modifications in electronic and structural properties such as band splitting, metal–insulator transitions (MIT), and enhanced electron lattice interactions. Notably, the charge-disproportionation is a key electronic feature that drives the MIT. Here, we demonstrate and quantify it using first-principles calculations combining density-functional theory, dynamical mean-field theory, and spin–lattice Monte Carlo simulations. Taking Ca 2 FeMnO 6 as a prototypical example of a correlated oxide, our ab initio study shows that MIT in Ca 2 FeMnO 6 arises from the partial localization of oxygen ligand holes at alternate Fe sites that control both charge and magnetic ordering. Interestingly, the band gap was found to be fundamentally controlled by the strength of the charge-transfer energy and not by the Mott–Hubbard interactions. The novel physics and insights presented in this work reveal promising routes for tuning the electronic functionality in transition-metal oxides.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Leveraging Polymorphism in YbCuBi to Map Transport and Elastic Properties

AMX Zintl compounds with the hexagonal ZrBeSi structure have gained significant attention for their remarkable vacancy tolerance and low thermal conductivity. Their 2D honeycomb sublattice, composed of M–X covalent bonds, is believed to contribute to high anharmonicity and unusual thermal transport properties. In this study, we explore the temperature-dependent polymorphism of YbCuBi as a model system to investigate the relationship between the structure and elastic and thermal transport properties in AMX Zintls. YbCuBi undergoes a structural transition from the “flat” Cu–Bi layers in the ZrBeSi structure to corrugated layers in the LiGaGe structure below 410 K, resulting in a distortion of its centrosymmetric structure. To probe the effects of this crystallographic transition, we employ inelastic neutron scattering and temperature-dependent resonant ultrasound spectroscopy. These experimental findings, coupled with first-principles calculations and thermal conductivity measurements, allow us to elucidate a direct relationship between corrugation of the honeycomb lattice and the observed changes in elastic and thermal transport properties. These insights can be extended to other Zintl phases with similar structure types, providing a platform for the rational design of functional materials with tailored thermal properties.

Chemical structure↗
DOE OSTI2026

Anisotropic Ferromagnetism in CrAu 3 Sb 6

The crystal structure and properties of CrAu 3 Sb 6 are presented, determined by measurements on single crystal and polycrystalline samples and first-principles calculations. The trigonal structure (space group P$\overline{3}$1m) comprises a CdI 2 -like sublattice of AuSb 2 with Cr occupying octahedral holes in a fully ordered triangular array. It can be viewed as a variation of the interesting and well-known families of partially intercalated transition metal dichalcogenides, but with stronger interactions along the stacking direction evidenced by short Cr–Au distances. The compound is metallic and ferromagnetic with a Curie temperature of 164 K. A strong anomalous contribution to the Hall effect is seen in the ferromagnetic state, and quantum oscillations are observed in magnetization at 2 K. Magnetization measurements show that the ordered moments of 1.5 μ B per Cr are oriented along the c-axis with relatively strong magnetocrystalline anisotropy. Electronic structure calculations confirm this uniaxial anisotropy and the important role of spin–orbit coupling in CrAu 3 Sb 6 and reveal strongly favored ferromagnetic ground state consistent with the measured Curie temperature. Through combined experiment and theory, this work provides a detailed picture of the basic properties and behaviors of this uniquely structured, Cr-based, anisotropic ferromagnet.

Crystal structure↗
DOE OSTI2024

Proton Conducting Neuromorphic Materials and Devices

Neuromorphic computing and artificial intelligence hardware generally aims to emulate features found in biological neural circuit components and to enable the development of energy-efficient machines. In the biological brain, ionic currents and temporal concentration gradients control information flow and storage. It is therefore of interest to examine materials and devices for neuromorphic computing wherein ionic and electronic currents can propagate. Protons being mobile under an external electric field offers a compelling avenue for facilitating biological functionalities in artificial synapses and neurons. In this review, we first highlight the interesting biological analog of protons as neurotransmitters in various animals. We then discuss the experimental approaches and mechanisms of proton doping in various classes of inorganic and organic proton-conducting materials for the advancement of neuromorphic architectures. Since hydrogen is among the lightest of elements, characterization in a solid matrix requires advanced techniques. We review powerful synchrotron-based spectroscopic techniques for characterizing hydrogen doping in various materials as well as complementary scattering techniques to detect hydrogen. First-principles calculations are then discussed as they help provide an understanding of proton migration and electronic structure modification. Outstanding scientific challenges to further our understanding of proton doping and its use in emerging neuromorphic electronics are pointed out.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Iridate Li8IrO 6 : An Antiferromagnetic Insulator

In this work, a polycrystalline iridate Li 8 IrO 6 material was prepared via heating Li 2 O and IrO 2 starting materials in a sealed quartz tube at 650 °C for 48 h. The structure was determined from Rietveld refinement of room-temperature powder neutron diffraction data. Li 8 IrO 6 adopts the nonpolar space group $R\bar{3}$ with Li atoms occupying the tetrahedral and octahedral sites, which is supported by the electron diffraction and solid-state 7Li NMR. This results in a crystal structure consisting of LiO 4 tetrahedral layers alternating with mixed IrO 6 and LiO 6 octahedral layers along the crystallographic c -axis. The +4 oxidation state of Ir 4+ was confirmed by near-edge X-ray absorption spectroscopy. An in situ synchrotron X-ray diffraction study of Li 8 IrO 6 indicates that the sample is stable up to 1000 °C and exhibits no structural transitions. Magnetic measurements suggest long-range antiferromagnetic ordering with a Néel temperature (TN) of 4 K, which is corroborated by heat capacity measurements. The localized effective moment $μ_{\text{eff}}$ (Ir) = 1.73 $μ_{\text{B}}$ and insulating character indicate that Li 8 IrO 6 is a correlated insulator. First-principles calculations support the nonpolar crystal structure and reveal the insulating behavior both in paramagnetic and antiferromagnetic states.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Theoretical Evaluation of Metal–Ligand Bonding in Neptunium Compounds in Relation to 237 Np Mössbauer Spectroscopy

The 237 Np Mössbauer isomer shift and quadrupole splitting (QS) are powerful probes for the metal–ligand bonding of neptunium, a 5f-element of vital importance in the nuclear fuel cycle. A large set of Np compounds with different oxidation states (III) to (VII) is studied to investigate, by first-principles calculations, isomer shifts and the QS trends in relation to the Np oxidation state. Natural Bond Orbital analysis reveals that in addition to donation bonding to the 5f shell, participation of the 6d and 7s neptunium shells in covalent (donation) bonding substantially impacts the isomer shifts. The isomer shift cannot be interpreted solely by the 5f shell electron count. The isomer shift for Np(II) compounds is estimated to be in the range of 31–34 mm/s, less positive than for Np(III) compounds. For the QS, density functional calculations fail to reproduce the quadrupole splitting for some Np(VI) ionic solids. A multiconfigurational wave function approach reproduces the observed QS trends. The calculations give a semiquantitative interpretation of the trends for Np oxidation states (V) to (VII). Furthermore, the contrasting QS for standard and “reverse” neptunyl(VI), at the opposite extremes of the observed QS scale, arises predominantly from the different crystal environments.

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

Two Distinct Cu(II)-V(IV) Superexchange Interactions with Similar Bond Angles in a Triangular ?CuV2? Fragment

The strength and sign of superexchange interactions are often predicted on the basis of the bond angles between magnetic ions, but complications may arise in situations with a nontrivial arrangement of the magnetic orbitals. We report on a novel molecular tetramer compound [Cu(H2O)dmbpy]2[V2O2F8] (dmbpy = 4,4 '-dimethyl-2,2 ' bipyridyl) that is composed of triangular "CuV2" fragments and displays a spin gap behavior. By combining first-principles calculations and electronic models, we reveal that superexchange Cu-V interactions carry drastically diHerent coupling strengths along two Cu-F-V pathways with comparable bond angles in the triangular "CuV2" fragment. Counterintuitively, their strong disparity is found to originate from the restricted symmetry of the half-filled Cu dx2-y2 orbital stabilized by the crystal field, leading to one dominating antiferromagnetic Cu-V coupling in each fragment. We revisit the magnetic properties of the reported spin-gapped chain compound [enH2]Cu(H2O)2[V2O2F8] (enH2 = ethylene diammonium) containing similar triangular "CuV2" fragments, and the magnetic behavior of the molecular tetramer and the chain compounds is rationalized as that of weakly coupled spin dimers and spin trimers, respectively. This work demonstrates that fundamentally diHerent magnetic couplings can be observed between magnetic ions with similar bond angles in a single spin motif, thus providing a strategy to introduce various exchange interactions combined with low dimensionality in heterometallic Cu(II)-V(IV) compounds.

Complexes↗
DOE OSTI2022

Stabilization of Dinuclear Rhodium and Iridium Clusters on Layered Titanate and Niobate Supports

Atomically dispersed organometallic clusters can provide well-defined nuclearity of active sites for both fundamental studies as well as new regimes of activity and selectivity in chemical transformations. More recently, dinuclear clusters adsorbed onto solid surfaces have shown novel catalytic properties resulting from the synergistic effect of two metal centers to anchor different reactant species. Difficulty in synthesizing, stabilizing, and characterizing isolated atoms and clusters without agglomeration challenges allocating catalytic performance to atomic structure. Here, we explore the stability of dinuclear rhodium and iridium clusters adsorbed onto layered titanate and niobate supports using molecular precursors. Both systems maintain their nuclearity when Statistical analysis of HAADF-STEM images revealed that rhodium and iridium dimers had mean cluster-to-cluster distances very similar to what is expected from a random distribution of atoms over a large area, indicating that they are dispersed without aggregation. The stability of dinuclear rhodium clusters supported on titanate nanosheets was also investigated by X-ray absorption fine structure (EXAFS), DRIFTS, and first-principles calculations. Both X-ray absorption spectroscopy and HAADF-STEM simulations, guided by density functional theory (DFT)-optimized structure models, suggested that rhodium dimers adsorb onto the nanosheets in an end-on binding mode that is stable up to 100 degrees C under reducing conditions. Finally, this study highlights that crystalline nanosheets derived from layered metal oxides can be used as model supports to selectively stabilize dinuclear clusters, which could have implications for heterogeneous catalysis.

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

Electronic and Thermal Properties of the Cation Substitution-Derived Quaternary Chalcogenide CuInSnSe 4

Quaternary chalcogenides continue to be of interest for a variety of technological applications, with physical properties stemming from their structural complexity and stoichiometric variation. In certain structure types, partial vacancies on specific lattice positions present an opportunity to investigate electrical and thermal properties in light of these lattice defects. In this work, we investigated the structural, thermal, and electronic properties of CuInSnSe 4 , a material that belongs to a relatively unexplored class of quaternary chalcogenides with a defect adamantine crystal structure. First-principles calculations together with experimental measurements revealed a chalcopyrite-like structure with inherent vacancies and characteristic s–p and p–d orbital hybridizations in the electronic structure of the material. Cation disorder and lattice anharmonicity result in very low thermal conductivity with values significantly lower than those for related compositions. In conclusion, this work reveals the fundamental physical properties of a previously uninvestigated quaternary chalcogenide and may aid investigations of similar as well as other quaternary chalcogenide compositions.

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

Incommensurate Spiral Spin Order in CaMn 2 Bi 2 Observed via High-Pressure Neutron Diffraction

Here, high-pressure neutron diffraction is employed to investigate the magnetic behavior of CaMn 2 Bi 2 in extreme conditions. In contrast to antiferromagnetic ordering on Mn atoms reported at ambient pressure, our results reveal that at high pressure, incommensurate spiral spin order emerges due to the interplay between magnetism on the Mn atoms and strong spin–orbit coupling on the Bi atoms: sinusoidal spin order is observed at pressures as high as 7.4 GPa. First-principles calculations with a noncollinear spin orientation demonstrate band crossing behavior near the Fermi level as a result of strong hybridization between the d orbitals of Mn and the p orbitals of Bi atoms. Competing antiferromagnetic order is observed at different temperatures in the partially frustrated lattice. Theoretical models have been developed to investigate spin dynamics. This research provides a unique toolbox for conducting experimental and theoretical magnetic and spin dynamics studies of magnetic quantum materials via high-pressure neutron diffraction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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↗