Search NASA⌕ Search

SEARCH · Search NASA

Results for “molecular decomposition”

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 OSTI2025

A Combined Crossed Molecular Beam and Theoretical Investigation of the Elementary Reaction of Tricarbon (C 3 (X 1 Σ g + )) with Diacetylene (C 4 H 2 (X 1 Σ g + )): Gas Phase Formation of the Heptatriynylidyne Radical ( l -C 7 H(X 2 Π))

An elucidation of the underlying formation pathways to acyclic hydrocarbons such as polyynes (C n H 2 ), cumulenes (C n H 2 ), and linear resonantly stabilized linear radicals (l-C n H) is indispensable to understand the hydrocarbon chemistry in extreme low and high temperature environments. In this study, we exploited the crossed molecular beam technique to investigate the reaction of tricarbon C 3 (X 1 Σ g + ) with diacetylene (butadiyne; HCCCCH; X 1 Σ g + ) at a collision energy of 47 ± 1 kJ mol ⁻1 . The experimental data were merged with ab initio calculations of the singlet C 7 H 2 potential energy surface (PES) revealing that the reaction is initiated via the formation of an initial van der Waals reactant complex in the entrance channel. Subsequent rearrangements lead to various carbene-type and cyclic intermediates via ring-opening, ring-closure, and hydrogen migration processes eventually forming acyclic C 7 H 2 isomers prior to their barrierless unimolecular decomposition to the most stable linear isomer, heptatriynylidyne (C 7 H, X 2 Π) in an overall endoergic reaction (+57 kJ mol ⁻1 ). The reaction exhibits strong similarities to the tricarbon – acetylene (C 3 – C 2 H 2 ). Furthermore, the significant energy threshold suggests that the tricarbon reaction with (poly)acetylenes forming resonantly stabilized linear radicals are open in high-temperature environments such as combustion flames and circumstellar envelopes of carbon stars and planetary nebulae as their descendants; however, these reactions are closed in low-temperature environments as in cold molecular clouds and hydrocarbon-rich atmospheres of planets and their moons such as in Titan.

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

Investigation of the isoscalar response of 24 Mg to 6 Li scattering

24 Mg is a strongly deformed nucleus in the ground state. Deformation effects can be observed in the structure of the isoscalar giant monopole and quadrupole resonances. 24 Mg is also a nucleus that is well known to present different types of cluster-oscillation modes. Both giant resonances and cluster states are strongly populated by isoscalar transitions. To extract the E0, E1, and E2 transition strengths via 6 Li scattering. The 6 Li probe is a powerful tool for investigating the isoscalar nuclear response with a very favorable ratio of resonance-to-continuum background. Double-differential cross sections of 6 Li inelastic scattering, at the beam energy of 100 MeV/u, were measured in the excitation-energy range 10 – 40 MeV and scattering angles 0 - 3°. A multipole-decomposition analysis was performed for extracting the isoscalar E0 , E1, and E2 strength distributions. Results: The extracted multipole strengths were compared with predictions from consistent quasiparticle random phase approximation calculations. The theoretical predictions are in fair agreement with the experimental data. The E0 strength was also compared with results from antisymmetrized molecular dynamics calculations found in the literature. A few peaks in the experimental data might be associated with clustering in 24 Mg. Ground-state deformation effects were observed in the isoscalar giant monopole resonance (ISGMR) and isoscalar giant quadrupole resonance (ISGQR) distributions. The ISGMR strength is split in two peaks around 19 and 28 MeV. The ISGQR exhibits a pronounced peak at 20 MeV with a broadening at the low-energy region, similar to predictions from microscopic calculations. Signatures of excitation of cluster states were observed in the E0 response. Further studies including particle-decay measurements will be required to confirm the nature of the observed peaks.

20 ≤ A ≤ 38↗
DOE OSTI2025

Infrared spectral signatures of interfacial water at TiO 2 –electrolyte interfaces from deep potential molecular dynamics

Vibrational spectroscopy is a powerful tool for probing water at oxide–electrolyte interfaces, but its molecular interpretation can be challenging. Here, we employ deep potential long-range molecular dynamics simulations with layer-resolved spectral analysis to investigate the microscopic origins of the infrared (IR) response of water at the interface with anatase TiO 2 (101), a prototypical oxide surface. The calculated interfacial spectra exhibit characteristic modifications compared to bulk water IR spectra, including enhanced intensities, a red shifted and broadened stretching band, and a higher-frequency shoulder, in qualitative agreement with experiments. Spectral decomposition shows that these signatures originate mainly from the first interfacial water layer, dominated by surface-bound H 2 O at Ti 5C sites, with secondary contributions from the second layer. A moderate salt concentration (0.4 M NaCl) leaves both the interfacial structure and the spectra essentially unchanged, while tuning the pH strongly modulates the spectral intensity. We establish a scaling relation linking the spectral intensity to the surface water dissociation fraction and the dipole moment, both governed by interfacial electric fields. These findings provide a microscopic framework for interpreting IR spectra of oxide–electrolyte interfaces.

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

Porting fragmentation methods to GPUs using an OpenMP API: Offloading the resolution-of-the-identity second-order Møller–Plesset perturbation method

Here, using an OpenMP Application Programming Interface, the resolution-of-the-identity second-order Møller–Plesset perturbation (RI-MP2) method has been off-loaded onto graphical processing units (GPUs), both as a standalone method in the GAMESS electronic structure program and as an electron correlation energy component in the effective fragment molecular orbital (EFMO) framework. First, a new scheme has been proposed to maximize data digestion on GPUs that subsequently linearizes data transfer from central processing units (CPUs) to GPUs. Second, the GAMESS Fortran code has been interfaced with GPU numerical libraries (e.g., NVIDIA cuBLAS and cuSOLVER) for efficient matrix operations (e.g., matrix multiplication, matrix decomposition, and matrix inversion). The standalone GPU RI-MP2 code shows an increasing speedup of up to 7.5× using one NVIDIA V100 GPU with one IBM 42-core P9 CPU for calculations on fullerenes of increasing size from 40 to 260 carbon atoms using the 6-31G(d)/cc-pVDZ-RI basis sets. A single Summit node with six V100s can compute the RI-MP2 correlation energy of a cluster of 175 water molecules using the correlation consistent basis sets cc-pVDZ/cc-pVDZ-RI containing 4375 atomic orbitals and 14 700 auxiliary basis functions in ~0.85 h. In the EFMO framework, the GPU RI-MP2 component shows near linear scaling for a large number of V100s when computing the energy of an 1800-atom mesoporous silica nanoparticle in a bath of 4000 water molecules. The parallel efficiencies of the GPU RI-MP2 component with 2304 and 4608 V100s are 98.0% and 96.1%, respectively.

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

Heat transport in liquid water from first-principles and deep neural network simulations

In this work, we compute the thermal conductivity of water within linear response theory from equilibrium molecular dynamics simulations, by adopting two different approaches. In one, the potential energy surface (PES) is derived on the fly from the electronic ground state of density functional theory (DFT) and the corresponding analytical expression is used for the energy flux. In the other, the PES is represented by a deep neural network (DNN) trained on DFT data, whereby the PES has an explicit local decomposition and the energy flux takes a particularly simple expression. By virtue of a gauge invariance principle, established by Marcolongo, Umari, and Baroni, the two approaches should be equivalent if the PES were reproduced accurately by the DNN model. We test this hypothesis by calculating the thermal conductivity, at the GGA (PBE) level of theory, using the direct formulation and its DNN proxy, finding that both approaches yield the same conductivity, in excess of the experimental value by approximately 60%. Besides being numerically much more efficient than its direct DFT counterpart, the DNN scheme has the advantage of being easily applicable to more sophisticated DFT approximations, such as meta-GGA and hybrid functionals, for which it would be hard to derive analytically the expression of the energy flux. We find in this way that a DNN model, trained on meta-GGA (SCAN) data, reduces the deviation from experiment of the predicted thermal conductivity by about 50%, leaving the question open as to whether the residual error is due to deficiencies of the functional, to a neglect of nuclear quantum effects in the atomic dynamics, or, likely, to a combination of the two.

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

Dynamics of poly[n]catenane melts

Inspired by advances in the chemical synthesis of interlocking polymer architectures, extensive molecular dynamics simulations have been conducted to study the dynamical properties of poly[n]catenanes-polymers composed entirely of interlocking rings-in the melt state. Both the degree of polymerization (number of links) and the number of beads per ring are systematically varied, and the results are compared to linear and ring polymers. Herein, a simple Rouse-like model is presented, and its analytical solution suggests a decomposition of the dynamics into "ring-like" and "linear-like" regimes at short and long times, respectively. In agreement with this picture, multiple sub-diffusive regimes are observed in the monomer mean-squared-displacements even though interchain entanglement is not prevalent in the system. However, the Rouse-type model does not account for the topological effects of the mechanical bonds, which significantly alter the dynamics at intermediate length scales both within the rings and at the chain segment scales. The stress relaxation in the system is extremely rapid and may be conveniently separated into ring-like and linear-like contributions, again in agreement with the Rouse picture. However, the viscosity has a non-monotonic dependence on the ring size for long chains, which disagrees strongly with theoretical predictions. This unexpected observation cannot be explained in terms of chain disentanglement and is inconsistent with other measures of polymer relaxation. Possible mechanisms for this behavior are proposed and implications for materials design are discussed.

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

Internally Catalyzed Hydrogen Atom Transfer (I-CHAT)—A New Class of Reactions in Combustion Chemistry

The current paradigm of low-T combustion and autoignition of hydrocarbons is based on the sequential two-step oxygenation of fuel radicals. The key chain-branching occurs when the second oxygenation adduct (OOQOOH) is isomerized releasing an OH radical and a key ketohydroperoxide (KHP) intermediate. The subsequent homolytic dissociation of relatively weak O–O bonds in KHP generates two more radicals in the oxidation chain leading to ignition. Based on the recently introduced intramolecular “catalytic hydrogen atom transfer” mechanism (J. Phys. Chem. 2024, 128, 2169), abbreviated here as I-CHAT, we have identified a novel unimolecular decomposition channel for KHPs to form their classical isomers—enol hydroperoxides (EHP). The uncertainty in the contribution of enols is typically due to the high computed barriers for conventional (“direct”) keto–enol tautomerization. Remarkably, the I-CHAT dramatically reduces such barriers. The novel mechanism can be regarded as an intramolecular version of the intermolecular relay transfer of H-atoms mediated by an external molecule following the general classification of such processes (Catal. Rev.-Sci. Eng. 2014, 56, 403). Here, we present a detailed mechanistic and kinetic analysis of the I-CHAT-facilitated pathways applied to n-hexane, n-heptane, and n-pentane models as prototype molecules for gasoline, diesel, and hybrid rocket fuels. We particularly examined the formation kinetics and subsequent dissociation of the γ-enol-hydroperoxide isomer of the most abundant pentane-derived isomer γ-C5-KHP observed experimentally. To gain molecular-level insight into the I-CHAT catalysis, we have also explored the role of the internal catalyst moieties using truncated models. All applied models demonstrated a significant reduction in the isomerization barriers, primarily due to the decreased ring strain in transition states. In addition, the longer-range and sequential H-migration processes were also identified and illustrated via a combined double keto–enol conversion of heptane-2,6-diketo-4-hydroperoxide as a potential chain-branching model. To assess the possible impact of the I-CHAT channels on global fuel combustion characteristics, we performed a detailed kinetic analysis of the isomerization and decomposition of γ-C5-KHP comparing I-CHAT with key alternative reactions—direct dissociation and Korcek channels. Calculated rate parameters were implemented into a modified version of the n-pentane kinetic model developed earlier using RMG automated model generation tools (ACS Omega, 2023, 8, 4908). Simulations of ignition delay times revealed the significant effect of the new pathways, suggesting an important role of the I-CHAT pathways in the low-T combustion of large alkanes.

Biochemistry & Molecular Biology↗
DOE OSTI2026

Drought shifts dissolved organic matter sources from above- to belowground and stress-induced processes in Amazon white-sand forests

White-sand forests contribute significantly to dissolved organic matter (DOM) production in the central Amazon, forming blackwater rivers that dominate organic matter export from the Amazon basin to the ocean. Despite their importance in controlling DOM export, white-sand forests are understudied, and it remains unclear whether systematic changes in the formation of blackwater DOM occur and how seasonal variations and extremes like El Niño-associated droughts impact them. We collected soil porewater from two central Amazon white-sand forests for two years, spanning a wet La Niña year followed by an El Niño drought year. The molecular composition of DOM was analyzed using high-resolution mass spectrometry, and correlation network analysis was employed to identify ecologically meaningful DOM subsets. Using additional chemical characterization, database annotations, correlation with 14C-age of DOM and climatic variables, and ecological null modeling, we propose five distinct DOM sources: plant litter and throughfall, soil organic matter (SOM) decomposition, root exudation, and two drought response subsets of likely microbial and plant origin. During drought conditions, aboveground plant-derived compounds decreased, while SOM products, root exudates, and drought response compounds increased. These drought responses were qualitatively similar in both years but notably amplified in the drier El Niño year. Drought amplified deterministic control over DOM composition, indicating that DOM reflected directed biological responses and that future droughts are likely to generate similar shifts. Overall, drought substantially altered belowground carbon cycling by shifting DOM sources and inducing stress responses, effects expected to recur and potentially intensify under future climate scenarios.

Lange, Dan F.↗
DOE OSTI2024

High-temperature structure, elasticity, and thermal expansion of ε-ZrH 1.8

Zirconium hydride is a promising candidate material for nuclear microreactor applications as a solid-state moderator component, owing to its favorable neutronics properties and good thermal stability over other metal hydrides. Here, in the present work, the crystal structure, thermal expansion, and elastic properties of the hydrogen-rich ε phase hydride were measured at elevated temperatures in the range 300–900 K. Samples were prepared by direct hydriding Zircaloy-4 metal – a nuclear-grade zirconium alloy. Room-temperature lattice parameters agree well with those reported from literature for unalloyed zirconium hydride and fall within an observed quadratic H-content dependence. The coefficients of thermal expansion, determined from lattice expansion and dilatometry, agree well within our work but were about 30 % lower than those reported by others for unalloyed hydrides. Density functional theory-based molecular dynamics simulations were used to compare with thermal expansion and elasticity measurements. Results showed lattice parameter temperature dependence and slope of thermal expansion align with those from measurements. Based on diffraction scans at select temperatures, ε phase remained stable in air up to at least 770 K. Likewise, dilatometry showed smooth thermal expansion up to the thermal decomposition temperature around 950 K. The precise decomposition temperature was not determined via diffraction due to sparse scanning. The complete elastic property measurements were gathered for ε-phase Ziracloy-4 hydride for the first time. Young's modulus was lower compared to the metal and δ hydride phases. High-temperature elasticity measurements were limited to <350 K due to acoustic dissipation effects.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Gas-Phase Study of the Elementary Reaction of the D1-Ethynyl Radical (C 2 D; X 2 Σ + ) with Propylene (C 3 H 6 ; X 1 A') under Single-Collision Conditions

The bimolecular gas-phase reactions of the D1-ethynyl radical (C 2 D; X 2 Σ + ) with propylene (C 3 H 6 ; X 1 A’) and partially substituted D3-3,3,3-propylene (C 2 H 3 CD 3 ; X 1 A’) were studied under single collision conditions utilizing the crossed molecular beams technique. Combining our laboratory data with electronic structure and statistical calculations, the D1-ethynyl radical is found to add without barrier to the C1 and C2 carbons of the propylene reactant, resulting in doublet C 5 H 6 D intermediate(s) with lifetime(s) longer than their rotational period(s). These intermediates undergo isomerization and unimolecular decomposition via atomic hydrogen loss through tight exit transition states forming predominantly cis/trans-3-penten-1-yne ((HCC)CH=CH(CH 3 )) and to a minor amount 3-methyl-3-buten-1-yne ((HCC)C(CH 3 )=CH 2 ) via overall exoergic reactions. Although the title reaction does not lead to the cyclopentadiene molecule (c-C 5 H 6 , X 1 A 1 ), high temperature environments can convert the identified acyclic C 5 H 6 isomers through hydrogen atom assisted isomerization to cyclopentadiene (c-C 5 H 6 , X 1 A 1 ). Since both the ethynyl radical and propylene reactants have been observed in cold interstellar environments such as TMC-1 and the reaction is exoergic and all barriers lie below the energy of the separated reactants, the these C 5 H 6 product isomers are predicted to form in those low temperature regions.

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

Root exudate lipids: Uncovering chemodiversity and carbon stability potential

Root-derived carbon has been shown to contribute more to soil carbon stocks than aboveground litter. Yet the molecular chemodiversity of root exudates remains poorly understood due to limited characterization and annotation. In this study, we characterized the molecular chemodiversity and production of metabolites and lipids in root exudates from field grown mature tall wheatgrass (Thinopyrum ponticum). We discovered a diversity of lipids, including substantial levels of triacylglycerols (∼19 μg/g fresh root per min), fatty acyls, sphingolipids, sterol lipids, and glycerophospholipids, some of which have not been previously documented in root exudates. By integrating tandem mass spectral library searching and deep learning-based chemical class assignment, our metabo-lipidomics approach significantly expanded the known molecular diversity of root exudates. Rates of lipid derived carbon production were approximately double that of polar metabolites (lipids: 81.52 ± 13.81 vs polar metabolites: 38.41 ± 5.93 μg C g −1 fresh root mass min −1 ) with an order of magnitude higher carbon to nitrogen ratios (lipids: 459 ± 90 vs polar metabolites: 14.40 ± 0.58). Exudate lipids displayed highly negative nominal oxidation state of carbon (−1.182 to −1.909), indicating that these compounds may be less favorable for microbial decomposition. Together our results suggest the potential of root exudate lipids to contribute to stable carbon pools in soil, supporting long-term carbon storage. This work advances understanding of plant-derived lipid inputs to soil and underscores the need for future studies on the functional roles of lipids in shaping root-microbe-soil interactions, microbial activity, soil structure, and nutrient availability – contributing to soil health.

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

Code modernization strategies for short-range non-bonded molecular dynamics simulations

Modern HPC systems are increasingly relying on greater core counts and wider vector registers. Thus, applications need to be adapted to fully utilize these hardware capabilities. One class of applications that can benefit from this increase in parallelism are molecular dynamics simulations. In this paper, we describe our efforts at modernizing the ESPResSo++ simulation package for molecular dynamics by restructuring its particle data layout for efficient memory accesses and applying vectorization techniques to benefit the calculation of short-range non-bonded forces, which results in an overall three times speedup and serves as a baseline for further optimizations. We also implement fine-grained parallelism for multi-core CPUs through HPX, a C++ runtime system which uses lightweight threads and an asynchronous many-task approach to maximize concurrency. Our goal is to evaluate the performance of an HPX-based approach compared to the bulk-synchronous MPI-based implementation. This requires the introduction of an additional layer to the domain decomposition scheme that defines the task granularity. On spatially inhomogeneous systems, which impose a corresponding load-imbalance in traditional MPI-based approaches, we demonstrate that by choosing an optimal task size, the efficient work-stealing mechanisms of HPX can overcome the overhead of communication resulting in an overall 1.4 times speedup compared to the baseline MPI version.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Molecular structures of dissolved and colloidal As V –Fe III complexes and their roles in the mobilization of As V under strongly acidic conditions

The effect of high concentration of iron (Fe III ) on the speciation and mobility of arsenic (As) under strongly acidic conditions remains unclear. This work studied the redistribution and speciation of As V and Fe III at Fe/As molar ratio of 1–14 and pH 1.5–2.0 in the dissolved, colloidal, and solid phases. Results showed that the elevated Fe III induced the decomposition of the precipitated poorly crystalline ferric arsenate by forming dissolved (< 3 kDa) and colloidal (3 kDa–0.1 µm) As–Fe complexes. The fraction of particulate As (> 0.1 µm) decreased from 70–90% to less than 20% when the Fe/As molar ratio increased from 1 to 14. The particle size of the bulk samples decreased significantly with the increase of Fe III concentration. The FTIR results suggested that AsV in dissolved/colloidal As–Fe complexes dominantly occurred as HAsO 4 2– species. The EXAFS results indicated that each HAsO 4 2– coordinated with approximately two Fe atoms in dissolved/colloidal As–Fe complexes at Fe/As ≥ 2. Furthermore, the findings suggest that high aqueous Fe III concentration can promote the mobility of As by forming dissolved/colloidal Fe–As complexes in acidic waters, potentially accelerating As transport from source to downstream in acid mine drainage systems.

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

Confocal Raman Microscopy as a Probe of Material Deconstruction in Processed Low-Density Polyethylene Particles

Confocal Raman microscopy was applied to detect structural change within individual particles of low-density polyethylene (LDPE) following chemical and electrochemical processing steps that aimed to facilitate material decomposition. A high numerical aperture (NA) oil-immersion objective enabled depth-profiling through the near surface region (20 μm–40 μm) of irregularly shaped particles with an axial spatial resolution < 2 μm estimated from measurements of instrument detection efficiency profiles. Changes in vibrational bands sensitive to polyethylene crystallinity were evident following treatments and linked to the release of low molecular weight compounds present as additives and products of processing. Effects of processing were probed by monitoring the rise of Raman scattering intensity in vibrational modes associated with polyethylene chains in a zig-zag (trans) conformation near 1128 cm –1 , 1294 cm –1 , and 1418 cm –1 , signaling chain clustering and development of organized, crystalline-like assemblies. Pristine LDPE particles displayed a uniform structure across the near surface region, while particles treated initially with chemical extractant and then further processed displayed increasingly enhanced crystallinity up to the maximum depth probed (40 μm). As a step toward measurements on ensembles of particles, least squares modeling was adapted to derive pure component spectra reflecting crystallinity change within spectral datasets. The work demonstrates high spatial resolution Raman depth-profiling for the characterization of processed polymers using a high NA immersion objective to overcome the limitations of air-objectives often used for confocal Raman microscopy.

Wahiduzzaman, Md. [Department of Chemistry and Bio↗
DOE OSTI2022

SEI formation mechanisms and Li + dissolution in lithium metal anodes: Impact of the electrolyte composition and the electrolyte-to-anode ratio

The lithium metal battery is one of today's most promising high-energy-density storage devices. Its full-scale implementation depends on solving operational and safety issues intrinsic to the Li metal high reactivity leading to uncontrolled electrolyte decomposition and uneven Li deposition. In this work, we study the spontaneous formation of the solid electrolyte interphase (SEI) upon contact of Li metal with the electrolyte and describe the heterogeneous SEI morphological features. Multiple electrolyte formulations based on lithium bis(fluorosulfonyl)imide (LiFSI), dimethoxyethane (DME), dimethyl carbonate (DMC), 1,2,2-tetrafluoroethyl-2,2,3,3-tetrafluoropropyl ether (TTE) and bis(2,2,2-trifluoroethyl) ether (BTFE) are used. Findings include the description of the SEI evolution from dispersed LiO, LiS, LiN, and LiF clusters to a continuous and compact inorganic phase in which the LiO and LiF content depend on the presence of fluorine diluents. Further, the role of the DME ether solvent helping the growth of a “wet-SEI” is compared to that of the highly unstable carbonate DMC, which decomposes into complex radical oligomers that might contribute to further electrolyte decomposition. The impact of the electrolyte-to-anode ratio on LiFSI decomposition is highlighted. Finally, we suggest the existence of a critical LiFSI concentration and electrolyte-to-anode ratio that could potentially balance the rate of electrolyte depletion and lithium consumption.

25 ENERGY STORAGE↗
DOE OSTI2026

The Interplay of Pauli Repulsion, Electrostatics, and Field Inhomogeneity for Blueshifting and Redshifting Vibrational Probe Molecules

Many molecules' vibrational frequencies are sensitive to intermolecular electric fields, enabling them to probe the field in complex molecular environments. However, it is often unclear whether the probe is responding to the local electric field or other types of intermolecular interactions, inhibiting interpretation of the frequency and effectiveness as probes. This is especially true for molecules whose vibrational frequencies blueshift instead of the more typical redshift in hydrogen bonding configurations. In this work, we computationally investigate the causes of redshifting versus blueshifting over a range of vibrational reporters. First, we apply adiabatic energy decomposition analysis to a paradigmatic set of probes, finding that redshifting only occurs when electrostatic interactions are strong enough to overcome the dominant and large blueshifting contribution of Pauli repulsion. Furthermore, we demonstrate that field inhomogeneity can further shift the frequency of many probes substantially to either reinforce or counteract the shift expected from a homogeneous field. We find that redshifting is reinforced by electric field inhomogeneity, otherwise field inhomogeneity further weakens the electrostatic contribution relative to Pauli repulsion, leading to blueshifting. Further calculations indicate that the probe's response to field inhomogeneity can be understood by considering the mass of the atoms involved in the stretching mode and sign of the electric field. In explaining the interplay of different intermolecular interactions and field inhomogeneity for many probes, our results should enable the use and interpretation of spectroscopic probes and their connection to electric fields in more complex systems.

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

Exchange Coupling Determines Metal-Dependent Efficiency for Iron- and Cobalt-Catalyzed Photochemical CO 2 Reduction

Catalysts promoting multielectron charge delocalization offer selectivity for the CO 2 reduction reaction (CO 2 RR) over the competing hydrogen evolution reaction. Here, we show metal–ligand exchange coupling as an example of charge delocalization that can determine the efficiency for photocatalytic CO 2 RR. A comparative evaluation of iron and cobalt complexes supported by the redox-active ligand tpyPY2Me establishes that the two-electron reduction of [Co(tpyPY2Me)] 2+ ([Co] 2+ ) occurs at potentials 770 mV more negative than the [Fe(tpyPY2Me)] 2+ ([Fe] 2+ ) analogue by maximizing the exchange coupling in the latter compound. The positive shift in the reduction potential promoted by metal–ligand exchange coupling drives [Fe] 2+ to be among the most active and selective molecular catalysts for photochemical CO 2 RR reported to date, maintaining up to 99% CO product selectivity with total turnover numbers (TONs) and initial turnover frequencies exceeding 30,000 and 900 min –1 , respectively. In contrast, [Co] 2+ shows much lower CO 2 RR activity, reaching only ca. 600 TON at 83% CO product selectivity under similar conditions accompanied by rapid catalyst decomposition. Further, the spin density plots of the two-electron reduced [Co] 0 complex implicate a paramagnetic open-shell doublet ground state compared to the diamagnetic open-shell singlet ground state of reduced [Fe] 0 , rationalizing the observed negative shift in two-electron reduction potentials from the [M] 2+ species and lowered CO 2 RR efficiency for the cobalt complex relative to its iron congener. This work emphasizes the contributions of multielectron metal–ligand exchange coupling in promoting effective CO 2 RR and provides a starting point for the broader incorporation of this strategy in catalyst design.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗