Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “DFT 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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

SPARC-X: Quantum simulations at extreme scale - reactive dynamics from first principles

We have developed the massively parallel electronic structure code SPARC-X: a computational framework for performing Kohn-Sham Density Functional Theory (DFT) calculations that can scale linearly with the number of atoms in the system, while being able to leverage petascale and emerging exascale parallel computers to study chemical phenomena at unprecedented length and time scales. SPARC-X exploits a recent breakthrough in electronic structure methodologies: systematically improvable, strictly local, orthonormal, discontinuous real-space bases that efficiently and systematically capture the local chemistry of the system. With further adaptation using new machine-learning techniques and the use of the massively parallel Spectral Quadrature (SQ) electronic structure method, the algorithmic complexity and prefactor associated with DFT calculations involving semilocal as well as hybrid functionals are dramatically reduced. Using petascale computational resources, SPARC-X enables quantum mechanical simulations at length and time scales previously accessible only by empirical approaches, e.g., 1,000,000 atoms for a few picoseconds using semilocal functionals or 1,000 atoms for a few picoseconds using hybrid functionals. Using exascale resources, the sizes and times targeted are two orders of magnitude larger. Such a capability has applications in a wide variety of chemical sciences, including reactive interfaces where large length- and/or long time-scales are needed and traditional force fields fail. This is particularly important in dynamic catalysis, where bond breaking and formation must be understood in detail. We developed, tested, and employed the SPARC-X framework to understand the photocatalytic properties of TiO 2 nanoparticles, revealing finite size effects that cannot be captured with standard model systems or functionals. This integrated development and application strategy ensures that SPARC-X remains a robust, efficient, and scalable software package for quantum simulations on current petascale and emerging exascale computing resources.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Ensemble Effects on Hydroxide Bond Dissociation Free Energies in Polyoxovanadate Clusters

Understanding structure-property relationships is foundational to numerous modern chemistries, such as proton-coupled electron transfer (PCET). However, an experimentally measured property is the result of the behavior from an ensemble of molecules. Neglecting ensemble effects, especially under complex chemical environments, may obfuscate these relationships and lead to discrepancies between theory and experiment. In this work, we demonstrate the impact of configurational entropy and local chemical environments on hydroxide bond dissociation free energies [BDFE- (O−H)] for a set of polyoxovanadate nanoclusters, at ambient conditions. The O−H bond strengths are investigated via density functional theory (DFT) coupled with statistical thermodynamic analysis and bilinear modeling, and compared with previous experimental results on the same systems, namely electrochemical solutions of: [V 6 O 13−x (OH) x (TRIOL R ) 2 ] −2 (x = 2, 4, 6; R = NO 2 , Me) and [V 6 O 11−x (OMe) 2 (OH) x (TRIOL NO 2 ) 2 ] −2 (x = 2, 4). Interestingly, we find that ensemble effects, even at room temperature, can account for a significant portion of the BDFE(O−H) trend with the degree of reduction via H atom binding, which cannot be fully captured by single-structure, static DFT calculations. Moreover, we find that the ensemble effects may be replicated statistically, requiring only enumeration of energetically accessible H-binding sites. With the ensemble effects resolved, we present a simple bilinear model to reconcile remaining biases between experiment and ensemble-informed theory, which corelate with clusterspecific electronic environment differences. The bilinear model achieves outstanding accuracy vs experiments with a root-mean squared error of 0.4 kcal/mol. Finally, based on the physicochemical characteristics of hydrogen interaction with polyoxometalates, we present a simple methodology that captures the BDFE(O−H) trend while dramatically reducing required DFT calculations by 98% and achieving accuracy within 1 kcal/mol. Overall, this work elucidates the roles and structural origins of configurational entropy and chemical effects on polyoxometalate hydroxide bond energies, with potential applicability to various atomically precise metal oxide systems. Importantly, it introduces models for rapid and highly accurate property calculations in connection with experiments.

Cluster chemistry↗
DOE OSTI2025

Carboxylic acid induced restructuring of the Fe 3 O 4 (001) surface

The redox properties of Fe 3 O 4 surfaces are central to many catalytic processes and enable dynamic, reduction-induced morphological restructuring during reactions. Here, we investigate the structural changes of the Fe 3 O 4 (001) surface during the decomposition of formic and acetic acids using scanning tunneling microscopy, X-ray photoelectron spectroscopy, and density functional theory (DFT) calculations. Both acids readily deprotonate, forming ordered carboxylate overlayers on the surface. Product formation pathways involve the removal of lattice oxygen, resulting in extensive surface restructuring. For formic acid, only a modest level of surface oxygen removal (∼3%) is observed, resulting in elongated pits along the octahedral Fe rows and exhibiting an aspect ratio of ∼3. In contrast, acetic acid induces more extensive reduction, with the removal of ∼20% of surface oxygen, yielding significantly larger pits while maintaining a similar aspect ratio. Repeated exposure to acetic acid further enlarges the pits, indicating preferential etching at step edges. DFT calculations reveal a mechanistic sequence in which lattice oxygen removal destabilizes adjacent Fe atoms, promoting their migration into the bulk and subsequent pit propagation and step edge formation. Together, these findings provide atomistic insights into the coupling between carboxylic acid conversion and oxide surface restructuring, underscoring the strong interplay between redox chemistry and morphological changes on catalytically active Fe 3 O 4 surfaces.

Ortiz-Garcia, José J. [Pacific Northwest National ↗
DOE OSTI2025

Machine Learning-Assisted Recovery of Delicate Kinetic Information from Transient Reactor Experiments

Identifying active sites and their roles in chemical reaction steps remains a vital challenge in heterogeneous catalysis. Transient experiments offer a unique way to probe active sites and distinguish subtle kinetic features. Although physics-based analysis methods may be well-developed, they can be highly susceptible to experimental noise, and smoothing methods may erase or even distort important features; a smooth curve is not always the best curve. We demonstrate a new workflow for the direct interpretation of intrinsic kinetic information from exit flux curves measured in transient reactor experiments. This workflow contains three artificial neural networks (ANNs), including a noise reducer, a concentration predictor, and a rate predictor to analyze experimental data, followed by the virtual TAP (VTAP) physics-based reactor model and density functional theory (DFT) calculations of adsorption energies on specific sites. We use this workflow to analyze the data from experiments titrating Pt/Al 2 O 3 and Pt/SiO 2 catalysts with carbon monoxide (CO) in the temporal analysis of products (TAP) reactor. Our workflow separates the time-evolving chemical reaction and mass transfer information contained in the TAP pulse response. The existence of strong- and weak-binding sites on the Pt/Al 2 O 3 catalyst is observed in the catalyst titration experiment in the transient reactor. The structures of the strong- and weak-binding sites are then identified by using DFT calculations. We find that the Pt/SiO 2 catalyst has only strong-binding sites, which aligns with the inactive support effect of SiO 2 . We demonstrate how machine learning methods provide unique insights with high-resolution data analysis that cannot be achieved by using state-of-the-art physics-based methods.

Adsorption↗
DOE OSTI2025

2025 ASMS Investigation of the Collision-Induced Dissociation Mechanism of Protonated TODGA with IRIS

Title (20 words): Investigation of the Collision-Induced Dissociation Mechanism of Protonated TODGA with IRIS Introduction (120 words): One of the challenges facing wide-spread adoption of nuclear power is the development of efficient separation processes for used nuclear fuel. The molecules in separation processes are subjected to an extreme environment due to the high radiation fields from the used fuel and highly acidic media used for fuel dissolution, which results in significant molecular degradation, leading to reduced process efficiency. These degradation products must be identified and studied so mitigation strategies can be developed to maintain process efficiency. However, complex systems can have many degradation products, complicating identification. Untargeted analysis tools could be used to understand radiation chemistry in complex systems. However, this would necessitate improved understanding of the gas-phase fragmentation mechanisms of fuel cycle molecules like tetraoctyldiglycolamide (TODGA). Methods (120 words): The gas-phase fragmentation of protonated TODGA was investigated using collision-induced dissociation (CID), resonance ejection, and infrared ion spectroscopy (IRIS). CID and resonance ejection experiments were conducted using a Bruker Daltonics (Bremen, Gemany) SolariX XR fourier transform ion cyclotron resonance (FT-ICR) mass spectrometer. IRIS spectra of protonated TODGA and its two CID fragmentation products were measured using a modified Bruker amaZon Speed ETD 3D quadrupole ion trap mass spectrometer coupled to the Free Electron Lasers for Infrared eXperiments (FELIX) free electron laser. Measured spectra were compared with density functional theory (DFT) calculations using the Gaussian 16, Revision C.02 software package with the ?B97X-D functional and def2-TZVPP basis sets. Candidate structures were generated using the CREST 3.0 conformational sampling software tool. Preliminary Data (300 words): Collision-induced dissociation of protonated TODGA ([C36H73N2O3]+, m/z=581.562) results two fragment ions, one at m/z=340.285 assigned as [C20H38NO3]+ and the other at m/z=312.290, assigned as [C19H38NO2]+. Based on the assigned formula and the structure of protonated TODGA, the fragment at m/z=340.285 is likely formed from elimination of neutral dioctylamine. Comparison of the IRIS spectrum of m/z=340.285 with DFT predictions suggests it contains a ring structure, and is assigned as N-octyl-N-(6-oxo-1,4-dioxan-2-ylidene)octan-1-aminium. Based on this structure and the structure of protonated TODGA, we hypothesize this fragment formed from elimination of neutral dioctylamine followed by a ring closure mechanism. Comparison of the IRIS spectrum of the fragment at m/z=312.290 with DFT predictions also indicated the presence of a ring structure, assigned as N-(1,3-dioxolan-4-ylidene)-N-octyloctan-1-aminium. This product could be formed from elimination of carbon monoxide from the ring of m/z=340.285 as a sequential fragmentation or formed directly from protonated TODGA via elimination of neutral N,N-dioctylformamide followed by a ring closure. Resonance ejection experiments where m/z=340 was continuously ejected from the IRC cell showed no decrease in intensity of m/z=312.290 across several collision energies, suggesting that the later, direct formation mechanism, dominates. The location of the ionizing proton in protonated TODGA is important for modeling the fragmentation mechanisms. DFT calculations suggested that the position of bands involving the coupled vibrations of the amide C—N and C=O bonds in TODGA are the most sensitive to proton location. Evaluation of the IRIS spectrum of protonated TODGA suggests that the ionizing proton is located between the two amid oxygens. This protonation location was calculated to lie approximately 30 kJ/mol lower in energy than the next lowest energy location, with the proton located solely on one of the amide oxygens. Novel aspect (20 words): Infrared ion spectroscopy combined with resonance ejection experiments and density functional theory to probe the collision-induced dissociation mechanism of tetraoctyldiglycolamide.

37 - INORGANIC, ORGANIC, PHYSICAL AND ANALYTICAL C↗
DOE OSTI2023

Polarity‐Driven Atomic Displacements at the 2D Mg 2 TiO 4 ‐MgO (001) Oxide Interface for Hosting Potential Interlayer Excitons

Abstract Interlayer excitons in solid‐state systems have emerged as candidates for realizing novel platforms ranging from excitonic transistors and optical qubits to exciton condensates. Interlayer excitons have been discovered in 2D transition metal dichalcogenides, with large exciton binding energies and the ability to form various van der Waals heterostructures. Here, an oxide system consisting of a single unit cell of Mg 2 TiO 4 on MgO (001) is proposed as a platform for hosting interlayer excitons. Using a combination of density functional theory (DFT) calculations, molecular beam epitaxy growth, and in situ crystal truncation rod measurements, it is shown that the Mg 2 TiO 4 ‐MgO interface can be precisely controlled to yield an internal electric field suitable for hosting interlayer excitons. The atoms in the polar Mg 2 TiO 4 layers are observed to be displaced to reduce polarity at the interface with the non‐polar MgO (001) surface. Such polarity‐driven atomic displacements strongly affect electrostatics of the film and the interface, resulting in localization of filled and empty band‐edge states in different layers of the Mg 2 TiO 4 film. The DFT calculations suggest that the electronic structure is favorable for localization of photoexcited electrons in the bottom layer and holes in the top layer, which may bind to form interlayer exciton states.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Probing the elastic modulus and hardness of superionic boron cluster solid electrolytes

Polyhedral borane salts are highly tunable compounds with the potential to be used as solid electrolytes due to their high Li and Na superionic conductivity and relatively wide electrochemical stability window. In considering their application to all solid-state batteries, their mechanical properties play a critical role in their operation due to stresses the solid electrolyte is subjected to during battery operation. Density functional theory (DFT) calculations have provided an initial assessment of bulk and elastic moduli for some selected boron cluster solid electrolytes, but direct measurements are still scarce. Here, in this paper, we report the elastic moduli and hardnesses of LiCB 11 H 12 , LiCB 9 H 10 , NaCB 11 H 12 , and NaCB 9 H 10 for the first time using nanoindentation continuous stiffness method (CSM) measurements under inert atmosphere. Experimental modulus values ranging 8.8–12.3 GPa and hardness values ranging 0.15–0.38 GPa are lower than that of other inorganic solid-state electrolytes such as oxide- and sulfide-based solid electrolytes. DFT calculations on the expected moduli of these compounds are also presented and discussed. Analysis of the indentation plasticity index reveals that these compounds have higher plasticity index compared to other common oxide and sulfide solid electrolytes. According to the calculated Pugh's ratio, all compounds in this study except LiCB 11 H 12 are considered ductile.

25 ENERGY STORAGE↗
DOE OSTI2023

Bimetallic-Derived Catalytic Structures for CO 2 -Assisted Ethane Activation

In recent years, the simultaneous upgrading of CO 2 and ethane has emerged as a promising approach for generating valuable gaseous (CO, H 2 , and ethylene) and liquid (aromatics and C 3 oxygenates) chemicals from the greenhouse gas CO 2 and large-reserved shale gas. The key challenges for controlling product selectivity lie in the selective C-H and C-C bond cleavage of ethane with the assistance of CO 2 . Bimetallic-derived catalysts likely undergo alloying or oxygen-induced segregation under reaction conditions, thus providing diverse types of interfacial sites, e.g., metal/support (M/M'O x ) interface and metal oxide/metal (M'O x /M) inverse interface, that are beneficial for selective CO 2 -assisted ethane upgrading. The alloying extent can be initially predicted by cohesive energy and atomic radius (or Wigner-Seitz radius), while the preference for segregation to form the on-top suboxide can be approximated using the work function, electronegativity, and binding strength of adsorbed oxygen. Furthermore, bimetallic-derived catalysts are typically supported on high surface area oxides. Modifying the reducibility and acidity/basicity of the oxide supports and introducing surface defects facilitate CO 2 activation and oxygen supplies for ethane activation. Using in situ synchrotron characterization and density functional theory (DFT) calculations, we found that the electronic properties of oxygen species influence the selective cleavage of C-H/C-C bonds in ethane, with electron-deficient oxygen over the metal (or alloy) surface promoting nonselective bond scission to produce syngas and electron-enriched oxygen over the metal oxide/metal interface enhancing selective C-H scission to yield ethylene. Here, we further demonstrate that the preferred structures of the catalyst surfaces, either alloy surfaces or metal oxide/metal inverse interfaces, can be controlled through the appropriate choice of metal combinations and their atomic ratios. Through a comprehensive comparison of experimental results and DFT calculations, the selectivity of C-C/C-H bond scission is correlated with the thermodynamically favorable bimetallic-derived structures (i.e., alloy surfaces or metal oxide/metal inverse interfaces) under reaction conditions over a wide range of bimetallic catalysts. These findings not only offer structural and mechanistic insights into bimetallic-derived catalysts but also provide design principles for selective catalysts for CO 2 -assisted activation of ethane and other light alkanes. This Account concludes by discussing challenges and opportunities in designing advanced bimetallic-derived catalysts, incorporating new reaction chemistries for other products, employing precise synthesis strategies for well-defined structures with optimized site densities, and leveraging time/spatial/energy-resolved in situ spectroscopy/scattering/microscopy techniques for comprehensive structural analysis. The research methodologies established here are helpful for the investigation of dynamic alloy and interfacial structures and should inspire more efforts toward the simultaneous upgrading of CO 2 and shale gas.

03 NATURAL GAS↗
DOE OSTI2024

Structural Studies of Alloyed and Nanoparticle Transition Metal Dichalcogenides by Selenium-77 Solid-State Nuclear Magnetic Resonance Spectroscopy

Layered transition metal dichalcogenides (TMDCs) such as MoS 2 , MoSe 2 , and WSe 2 are under intense investigation because they are atomically thin semiconductors with photophysical properties that can be tuned by changing their composition or morphology. Mechanochemical processing has been proposed as a method to obtain alloyed TMDCs in the series Mo 1–x W x S y Se 2–y (x = 0–1; y = 0–2). However, elucidating the chemical transformations occurring at the atomic scale following mechanochemical processing can be challenging because the products are often amorphous or microcrystalline. To address this challenge, we probe TMDC mixing and alloying by using a combination of powder X-ray diffraction, Raman spectroscopy, diffuse reflectance spectroscopy, 77 Se solid-state nuclear magnetic resonance (SSNMR) spectroscopy, and planewave density functional theory (DFT) calculations. The nature of the milling material and reaction atmosphere are shown to be essential factors in limiting the formation of undesired oxide byproducts. We demonstrate acquisition of 77 Se SSNMR spectra using different combinations of Carr-Purcell Meiboom-Gill acquisition (CPMG) pulse sequences, magic angle spinning (MAS), and MAS dynamic nuclear polarization. Further, the combination of SSNMR with the other characterization methods clearly demonstrates that high energy impact ball milling induces molecular level alloying of Mo, W and chalcogen atoms in the family Mo 1–x W x S y Se 2–y . Gauge including projector augmented wave DFT calculations yield accurate 77 Se chemical shift (CS) tensor components. 77 Se SSNMR spectroscopy was also applied to study the structure of WSe 2 nanocrystals intercalated with ethylenediamine. The intercalated WSe 2 nanocrystals exhibit a more positive isotropic 77 Se CS as compared to bulk WSe 2 , however, the 77 Se CS anisotropy is the same, confirming the WSe 2 layers have a similar structure as in their bulk counterparts.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Band gap and electronic structure of defects in the ternary nitride BP 3 N 6 : experiment and theory

Recent advances in methods to access nitride systems by a high-pressure high-temperature approach have made possible the one-step synthesis of mixed ternary non-metal nitrides. As a prerequisite to use in a practical device, it is important to understand important bulk electronic properties, such as the band gap, as well as characterizing the presence and effect of defects that are present. In this work, the novel ternary nitride BP 3 N 6 is studied using techniques sensitive to the partial electronic density of states, specifically X-ray absorption spectroscopy and X-ray emission spectroscopy. Complementary full-potential all-electron density functional theory (DFT) calculations allow important bulk electronic parameters, such as the band gap, to be elucidated. The band gap of BP 3 N 6 has been determined to be 3.9 ± 0.2 eV and 4.1 ± 0.4 eV at the B K- and N K-edges, respectively. This is close to a theoretical value of 4.3 eV predicted by the PBEsol exchange–correlation functional and considerably less than a value of 5.8 eV predicted by the modified Becke–Johnson exchange–correlation functional. X-Ray excited optical luminescence (XEOL) measurements are performed to interrogate the presence of point defects in this system. Together with DFT calculations, these measurements reveal the presence of nitrogen vacancies which lead to multiple mid-gap trap states.

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

Controlling Selective C–O and C–H Bond Scission of Methanol by Supporting Pt on TiN and Mo 2 N Model Surfaces and Powder Catalysts

Transition metal nitrides (TMNs) have been explored as effective supports for Pt due to their Pt-like electronic properties. However, there is a lack of fundamental understanding regarding the behavior of Pt on different TMNs (Pt/TMN). Herein two TMNs, Mo 2 N and TiN, were modified with Pt and compared using methanol decomposition as a probe reaction via both ultrahigh vacuum (UHV) studies on thin films and ambient-pressure batch reactor studies of powder catalysts. Temperature-programmed desorption (TPD) and high-resolution electron energy loss spectroscopy (HREELS) measurements were conducted under UHV conditions with Mo 2 N and TiN thin films. Mo 2 N was shown to favor C–H bond scission to form CO with a 56.2% selectivity, while TiN favored C–O bond scission to form CH 4 with a 74.5% selectivity. The addition of 0.9 monolayers (MLs) of Pt increased C–H bond scission selectivity to 89.7% and 49.2% for Mo 2 N and TiN respectively. Density functional theory (DFT) calculations on model surfaces revealed that the binding energy of O (BE *O ) was significantly reduced on Pt/TMNs, from −4.02 eV on Mo 2 N to −1.31 eV on Pt/Mo 2 N and −4.74 eV on TiN to −1.37 eV on Pt/TiN. As a result, C–O bond scission pathways were suppressed, leading to the preferential C–H bond scission that was observed experimentally. The C–O and C–H bond scission trends observed on thin films were then extended to powder catalysts, which demonstrated similar trends toward methanol decomposition. In conclusion, results from the current study establish that by combining UHV studies and DFT calculations over model surfaces, one can effectively predict the catalytic behavior of realistic TMN powder catalysts.

08 HYDROGEN↗
DOE OSTI2024

Structural anisotropy in Sb thin films

Sb thin films have attracted wide interest due to their tunable band structure, topological phases, high electron mobility, and thermoelectric properties. We successfully grow epitaxial Sb thin films on a closely lattice-matched GaSb(001) surface by molecular beam epitaxy. We find a novel anisotropic directional dependence on their structural, morphological, and electronic properties. The origin of the anisotropic features is elucidated using first-principles density functional theory (DFT) calculations. The growth regime of crystalline and amorphous Sb thin films was determined by mapping the surface reconstruction phase diagram of the GaSb(001) surface under Sb 2 flux, with confirmation of structural characterizations. Crystalline Sb thin films show a rhombohedral crystal structure along the rhombohedral (211) surface orientation parallel to the cubic (001) surface orientation of the GaSb substrate. At this coherent interface, Sb atoms are aligned with the GaSb lattice along the [1̄10] crystallographic direction but are not aligned well along the [110] crystallographic direction, which results in anisotropic features in reflection of high-energy electron diffraction patterns, misfit dislocation formation, surface morphology, and transport properties. Our DFT calculations show that the preferential orientation of the rhombohedral Sb (211) plane may originate from the GaSb surface, where Sb atoms align with the Ga and Sb atoms on the reconstructed surface. The formation energy calculations confirm the stability of the experimentally observed structures. Our results provide optimal film growth conditions for further studies of novel properties of Bi 1-x Sb x thin films with similar lattice parameters and an identical crystal structure, as well as functional heterostructures of them with III–V semiconductor layers along the (001) surface orientation, supported by a theoretical understanding of the anisotropic film orientation.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Band Energy Dependence of Defect Formation in the Topological Semimetal Cd3As2

Cadmium Arsenide (Cd3As2) is a prototypical Dirac semimetal that manifests topological properties in a 3D bulk material. In defect-free Cd3As2, the Fermi level lies at a minimum in the density of states at the Dirac point, but experimentally it forms with excess electron carriers and an elevated EF, thereby masking the topological features. To computationally study the self-doping of Cd3As2, we combine density functional theory (DFT) calculations for defect formation energies with quasi-particle self-consistent GW (QSGW) electronic structure calculations. We demonstrate an innate dependence of the point defect formation energies on carrier concentrations and use the QSGW calculated density of states to extrapolate formation energies to arbitrary electron concentrations. This approach allows the quantitative modeling of thermodynamic defect equilibria in topological semimetals and is used to predict how Cd3As2 growth conditions affect the position of EF relative to the Dirac point.

CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS↗
DOE OSTI2021

Relationship between A-site cation and magnetic structure in 3d–5d–4f double perovskite iridates Ln 2 NiIrO 6 ( Ln = La, Pr, Nd)

We report a comprehensive investigation of Ln 2 NiIrO 6 (Ln=La, Pr, Nd) using thermodynamic and transport properties, neutron powder diffraction, resonant inelastic x-ray scattering, and density-functional theory (DFT) calculations to investigate the role of A-site cations on the magnetic interactions in this family of hybrid 3d–5d–4f compositions. Magnetic structure determination using neutron diffraction reveals antiferromagnetism for La 2 NiIrO 6 , a collinear ferrimagnetic Ni and Ir state that is driven to long-range antiferromagnetism upon the onset of Nd ordering in Nd 2 NiIrO 6 , and a noncollinear ferrimagnetic Ni and Ir sublattice interpenetrated by a ferromagnetic Pr lattice for Pr 2 NiIrO 6 . For Pr 2 NiIrO 6 , heat-capacity results reveal the presence of two independent magnetic sublattices, and transport resistivity indicates insulating behavior and a conduction pathway that is thermally mediated. Here, a first principles DFT calculation elucidates the existence of the two independent magnetic sublattices within Pr 2 NiIrO 6 and offers insight into the behavior in La 2 NiIrO 6 and Nd 2 NiIrO 6 . Resonant inelastic x-ray scattering is consistent with spin-orbit coupling splitting the t 2 g manifold of octahedral Ir 4+ into a J eff = 1/2 and J eff = 3/2 state for all members of the series considered.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Mechanisms and Energetics of CO 2 Chemisorption in Choline-Based Eutectic Solvents

Understanding the influences of proton sharing and hydrogen transfer in amine and hydroxyl functionalized eutectic solvents is important for tuning CO 2 binding as they relate to the chemisorption capacity and energetics of solvent regeneration in carbon capture. In this study, we examine how the CO 2 binding mechanism and energetics vary in eutectic solvents composed of choline-based hydrogen-bond acceptors (HBAs) with amine and hydroxyl moieties and ethylene glycol (EG), propylene glycol (PG), 1,2-butanediol (BD), and monoethanolamine (MEA)-type hydrogen-bond donors (HBDs), through absorption–desorption experiments and spectral analysis, supported with density functional theory (DFT) calculations. While the HBD-CO 2 complex is the major and the thermodynamically more favored product of absorption in all of the examined eutectic solvents, this complexation is only possible proceeding a proton transfer from HBD to the HBA in diol-based systems (EG, PG, BD). The interplay between HBA-CO 2 and HBD-CO 2 speciation is governed by the HBD concentration, basicity, and steric effects, with increased HBD-CO 2 interactions in the presence of smaller HBDs (size: BD > PG > EG). On the other hand, in MEA-based systems, CO 2 readily binds to the amine followed by a proton transfer from the amine to the HBA as confirmed by DFT calculations. In this way, the MEA-CO 2 chemisorption capacities exceeded that of the conventional amine (up to 0.93 mol of CO 2 per mol of MEA versus 0.5). Furthermore, the increase in CO 2 capacities and enhanced binding in eutectics including MEA or the more basic HBAs also resulted in regeneration temperatures of 70 to 60 °C with calculated reaction energies of −82 to −57 kJ/mol, compared to 50 °C for others with weaker binding, following CO 2 absorption at 25 °C.

absorption↗
DOE OSTI2024

Defect-induced formation and frustration-driven multiple magnetic transitions in Gd 2 Co 0.90 Si 2.90

A new ternary compound, Gd 2 Co 0.90 Si 2.90 , has been synthesized in a chemically single phase by deliberately introducing lattice vacancies in the Co and Si sites. The system is characterized by DC magnetization, heat capacity, resistivity and density functional theory (DFT) calculations. Here, a detailed experimental study reveals that the system exhibits an antiferromagnetic transition below T N = 11.3 K followed by spin freezing behaviour below T B ~ 4.8 K. An additional high-temperature magnetic transition could also be detected at T H ~ 150 K, which is short range in character and is associated with defect-induced polarization of the conduction electrons. The spin-glass-like state formation in the system is additionally bolstered by the observation of magnetic relaxation and associated aging phenomena as well as the magnetic memory effect. On the basis of non-equilibrium dynamical behaviour, we argue that the glassy state in this compound favours the hierarchical model over the droplet model. DFT calculations and experimental outcomes establish that the ground state is magnetically frustrated due to the small energy difference between distinct spin configurations having closely spaced ground-state energies.

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

Mechanistic origin of solvent-dependent thermal stability in sodiated Sn anodes for sodium-ion batteries

Understanding the thermal stability of high-energy density alloy anodes is critical for the safe deployment of sodium-ion batteries (SIBs). Here, accelerating rate calorimetry (ARC), post-mortem characterizations, and density functional theory (DFT) calculations are combined to understand the thermal reactivity of fully sodiated Sn, Sn-hard carbon (HC) blends, and HC anodes in carbonate- and ether-based electrolytes. ARC measurements show that propylene carbonate (PC) causes earlier self-heating rate (SHR) onset and higher reactivity than tetraethylene glycol dimethyl ether (TEGDME), indicating inferior thermal stability. Sodiated Sn exhibits better thermal stability than sodiated HC, while Sn-HC blends show intermediate behavior that improves with increasing Sn content. Post-ARC analyses reveal desodiation of Na15Sn4 to metallic Sn with particle coalescence, whereas Sn-HC blends and HC retain finer morphologies. PC promotes Sn oxidation to SnO, while TEGDME suppresses oxide formation; NaPF6-containing electrolytes additionally form NaF. DFT calculations show that PC adsorption lowers Na extraction energy and enhances interfacial electronic interactions, facilitating Na release and reductive decomposition. These results establish a direct correlation between solvent-dependent reaction pathways and thermal stability in SIB alloy anodes.

Accelerating rate calorimetry↗
DOE OSTI2025

A Novel Solid Form of Erlotinib: Synthesis by Heterogeneous Complexation and Characterization by NMR Crystallography

We describe the synthesis of a novel complex of the anticancer "active pharmaceutical ingredient erlotinib (ERL) via heterogeneous nucleation on polished zinc tiles. The resulting product, ERL 2 ·ZnCl 2 , is characterized by single-crystal X-ray diffraction, multinuclear solid-state NMR (ssNMR) spectroscopy, and density functional theory (DFT) calculations. Also characterized are the hydrochloride salt (ERL·HCl) and monohydrate free base (ERL·H 2 O) forms of erlotinib. 13 C ssNMR spectroscopy is useful for site-by-site assignment and rapid fingerprinting, while also providing preliminary structural interpretations, such as the number of molecules in the asymmetric unit. 35 Cl ssNMR can readily differentiate between the chloride ions in ERL·HCl and the covalently bonded chlorine in ERL 2 ·ZnCl 2 . 15 N ssNMR proves to be critical here because of the large isotropic chemical shift differences between ERL·H 2 O, ERL·HCl, and ERL 2 ·ZnCl 2 . The 15 N chemical shift tensors are linked directly to differences in structure and bonding with the aid of DFT calculations. Together, these results demonstrate the utility of multinuclear NMR crystallography for the characterization of solid forms of APIs, especially when other analytical techniques face significant challenges.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗