Search NASA⌕ Search

DOE OSTI · 3009396

Confinement-controlled selective CO 2 insertion into a dicopper dihydride core: A multiscale mechanistic study

Abstract

CO 2 is an abundant C 1 feedstock for fuel and chemical synthesis. We have previously demonstrated experimentally a stepwise insertion of CO 2 into a [Cu 2 H 2 ] core via a solid–gas in crystallo reaction, forming formate species that are unstable and inaccessible under solution-phase conditions. This work elucidates how structural confinement within the crystal lattice enables such selective reactivity. In particular, co-crystallized tetrahydrofuran molecules induce site asymmetry around the [Cu 2 H 2 ] unit, modulating both the local electronic environment and CO 2 diffusion pathways. Using a multiscale computational approach that combines classical molecular mechanics, hybrid quantum mechanics/molecular mechanics molecular dynamics, and enhanced-sampling free energy calculations, we demonstrate how site asymmetry affects CO 2 binding affinities and reaction pathways. These results provide detailed mechanistic insight into CO 2 insertion and hydride transfer, highlighting key differences between crystal- and solution-phase pathways and offering a general framework for understanding how lattice confinement shapes chemical reactivity.

Explore related subjects

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

Fuller III, Jack T. [Pacific Northwest National Laboratory (PNNL), Richland, WA (United States)], Patrick, Evan A. [Pacific Northwest National Laboratory (PNNL), Richland, WA (United States)] (ORCID:0000000183229254), Schenter, Gregory K. [Pacific Northwest National Laboratory (PNNL), Richland, WA (United States)] (ORCID:0000000154445484), Ginovska, Bojana [Pacific Northwest National Laboratory (PNNL), Richland, WA (United States)] (ORCID:0000000301658202), Tran, Ba L. [Pacific Northwest National Laboratory (PNNL), Richland, WA (United States)] (ORCID:0000000176613696), Raugei, Simone [Pacific Northwest National Laboratory (PNNL), Richland, WA (United States)] (ORCID:0000000191188480). 2025-11-21. Confinement-controlled selective CO 2 insertion into a dicopper dihydride core: A multiscale mechanistic study. https://doi.org/10.1063/5.0287973

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related reports

Chemical bonding, phase stability and magnetic property in Sm 2 Fe 17 X 3 (X=H, C, N): A first-principles perspective

As a promising alternative to Nd–Fe–B magnets, the critical rare earth free Sm 2 Fe 17 X 3 (X = C, N) exhibits potential for high-performance magnets. However, their poor phase stability remains a major obstacle to developing bulk magnets. We investigated the phase stability and intrinsic magnetic properties of Sm 2 Fe 17 X 3 (X = H, C, N) using first-principles calculations and chemical bond analysis. The formation energies are negative, while the decomposition energies are −1.53, 0.348, and −0.74 eV per formula unit for X = H, C, and N, respectively, which is responsible for the weak thermal stability. Our chemical bond analysis reveals that the bonding asymmetry between Sm–X and Fe–X interactions creates local structural distortions and degrades the phase stability of Sm 2 Fe 17 X 3 . The project Crystal Orbital Hamilton Population (-pCOHP) analysis indicates that the Sm–X bonding remains positive up to the Fermi level, indicating stable bonding interactions. Here, in contrast, the Fe–X bonding becomes negative near the Fermi level, signifying anti-bonding contributions that reduce structural stability. Interstitial atoms X expand the lattice and enhance Fe magnetic moments, but Fe–X bonding suppresses neighboring Fe moments. Electron transfer from Sm to X modifies the valence state of Sm and the crystal field at the site, contributing to enhanced magnetocrystalline anisotropy in Sm 2 Fe 17 X 3 . Among the interstitial elements, carbon and nitrogen—with their larger atomic radius and higher electronegativity—induce greater lattice expansion and form stronger bonds with neighboring Sm and Fe atoms compared to hydrogen. Consequently, Sm 2 Fe 17 X 3 (X = C and N) exhibits better phase stability and significant improvement in magnetic properties.

Chemical bonding↗

First-principles insights into Si substitution effects in Sm 2 (Fe,Si) 17 C x magnet

The partial substitution of Fe by Si enhances the phase stability of Sm 2 Fe 17 C x magnets with x > 1.0. Here, we elucidate the Si-substitution scheme and its impact on phase stability and magnetic properties in Sm 2 (Fe,Si) 17 C 3 from first-principles calculations and chemical bonding analysis. The calculated substitution energies for Si at various Fe sites are negative, indicating improved phase stability. Si preferentially substitutes Fe atoms at the 9d site in Sm 2 (Fe,Si) 17 C 3 while it tends to enter the Fe 18h site in Sm 2 (Fe,Si) 17 . This difference in site preference is attributed to the distinct chemical environments surrounding the Fe (Si) sites in the two compounds. Si substitution favors the formation of Sm–Si bonds while minimizing the Si–C and Si–Si interactions. Crystal orbital Hamilton populations and crystal orbital bond index calculations indicate that the partial replacement of Fe with Si strengthens the chemical bonding of Sm–Fe 3 (18f) and Sm–Fe 4 (18h) and improves overall phase stability in Sm 2 (Fe,Si) 17 C 3 . Beyond the dilution effect, Si substitution also reduces the magnetic moments of neighboring Fe atoms, a phenomenon linked to the strong Fe–Si bonding. These findings highlight the dual role of Si in modifying both the structural and magnetic characteristics of Sm 2 Fe 17 -based magnetic compounds.

Chemical bonding↗

Machine learning-accelerated path integral molecular dynamics simulations of reactive organic electrolytes

Hydrogen bonded electrolytes that exhibit accelerated proton transport via sequential reactive hops have drawn interest for their promise in clean energy applications. Molecular dynamics simulations of these electrolytes offer the opportunity to uncover microscopic mechanistic details that could be used to design and tune the properties of candidate electrolyte technologies. However, accurately modeling the proton transfer reactions and transport properties that give rise to high charge conductivites in these electrolytes proves computationally challenging because of the need to perform lengthy condensed phase simulations, treating both the electronic and nuclear degrees of freedom quantum mechanically. In this paper, we demonstrate that such a modeling task can be efficiently achieved with the use of density functional theory (DFT)-trained machine learning potentials (MLP) to accelerate path integral molecular dynamics (PIMD) simulations. We highlight the practical utility of this approach by using it to benchmark how closely PIMD simulations employing different DFT exchange–correlation functionals reproduce the composition-dependent densities, diffusion coefficients, and electrical conductivities of mixtures consisting of imidazole and levulinic acid. Even with the speedup afforded by our MLPs, PIMD simulations remain quite expensive. Furthermore, in order to render PIMD more computationally tractable, we introduce and benchmark the accuracy of a ring polymer contraction approach that leverages a computationally efficient short-range MLP to accelerate our PIMD simulations by an additional factor of four.

Chemical bonding↗