SEARCH · Engineering Papers

Results for “MOLECULAR DISSOCIATION”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 883 records · Page 49

All sources NASA NTRS DOE OSTI

NASA NTRS2003

Effects of carbonyl bond, metal cluster dissociation, and evaporation rates on predictions of nanotube production in high-pressure carbon monoxide

The high-pressure carbon monoxide (HiPco) process for producing single-wall carbon nanotubes (SWNTs) uses iron pentacarbonyl as the source of iron for catalyzing the Boudouard reaction. Attempts using nickel tetracarbonyl led to no production of SWNTs. This paper discusses simulations at a constant condition of 1300 K and 30 atm in which the chemical rate equations are solved for different reaction schemes. A lumped cluster model is developed to limit the number of species in the models, yet it includes fairly large clusters. Reaction rate coefficients in these schemes are based on bond energies of iron and nickel species and on estimates of chemical rates for formation of SWNTs. SWNT growth is measured by the conformation of CO2. It is shown that the production of CO2 is significantly greater for FeCO because of its lower bond energy as compared with that of NiCO. It is also shown that the dissociation and evaporation rates of atoms from small metal clusters have a significant effect on CO2 production. A high rate of evaporation leads to a smaller number of metal clusters available to catalyze the Boudouard reaction. This suggests that if CO reacts with metal clusters and removes atoms from them by forming MeCO, this has the effect of enhancing the evaporation rate and reducing SWNT production. The study also investigates some other reactions in the model that have a less dramatic influence.

Models, Chemical↗

DOE OSTI2023

Evaluation of an Image-Derived Input Function for Kinetic Modeling of Nicotinic Acetylcholine Receptor-Binding PET Ligands in Mice

Positron emission tomography (PET) radioligands that bind with high-affinity to α4β2-type nicotinic receptors (α4β2Rs) allow for in vivo investigations of the mechanisms underlying nicotine addiction and smoking cessation. Here, we investigate the use of an image-derived arterial input function and the cerebellum for kinetic analysis of radioligand binding in mice. Two radioligands were explored: 2-[ 18 F]FA85380 (2-FA), displaying similar pKa and binding affinity to the smoking cessation drug varenicline (Chantix), and [ 18 F]Nifene, displaying similar pKa and binding affinity to nicotine. Time–activity curves of the left ventricle of the heart displayed similar distribution across wild type mice, mice lacking the β2-subunit for ligand binding, and acute nicotine-treated mice, whereas reference tissue binding displayed high variation between groups. Binding potential estimated from a two-tissue compartment model fit of the data with the image-derived input function were higher than estimates from reference tissue-based estimations. Rate constants of radioligand dissociation were very slow for 2-FA and very fast for Nifene. We conclude that using an image-derived input function for kinetic modeling of nicotinic PET ligands provides suitable results compared to reference tissue-based methods and that the chemical properties of 2-FA and Nifene are suitable to study receptor response to nicotine addiction and smoking cessation therapies.

2-FA85380↗

DOE OSTI2024

Investigating the impact of donor atoms of single-source precursors on f-element nanomaterial formation and mechanisms of initial bond cleavage

Nanotechnology is advancing exponentially, with novel discoveries and applications emerging every day across various fields, including environmental remediation, biomedicine, electronics, and textiles. Nanoparticles (NPs) are defined as particles with sizes ranging from 1 to 100 nm, and their properties are influenced by factors such as composition, phase, size, and shape. Single-source precursors (SSPs) offer a promising method for controlling the synthesis of phase-pure nanoparticles with the desired stoichiometry and have been successfully applied to prepare transition metal and lanthanide NPs. SSPs are molecular compounds that contain all the necessary components to form the final nanomaterial and decompose under mild conditions, such as heating, to create the solid-state product. The objective of this work was to garner a mechanistic understanding of cerium oxide NP formation from cerium SSPs. To achieve these goals, first, ligands with oxygen and sulfur donor atoms were synthesized; among the array of ligands synthesized, some were novel. These ligands were then chelated to cerium ions to synthesize cerium precursors. Initially, known cerium alkoxide precursors were synthesized. Subsequently, these basic cerium complexes were further manipulated through acid-base protonolysis reactions to install a diverse set of ligands. Within this body of work, 10 novel cerium precursors were successfully synthesized. The ligands and the SSPs were analyzed using nuclear magnetic resonance (NMR) and thermogravimetric analysis coupled with differential scanning calorimetry (TGA-DSC). These SSPs were then introduced to a collision-induced dissociation mass spectrometer (CID-MS), and their dissociation pathways were studied.

37 - INORGANIC, ORGANIC, PHYSICAL AND ANALYTICAL C↗

DOE OSTI2022

Ultrafast probing of charge transfer dynamics with an XUV free-electron laser

We study charge transfer in dissociating molecules of 2-iodopropane (C 3 H 7 I) in a UV pump – XUV probe regime at the SACLA free-electron laser. Signatures of charge transfer from the recently dissociated iodine atom to the neutral propyl radical are observed in recorded 3D ion momentum distributions. This signature has evaded observation in previous related experiments. Detailed analysis of these signatures allows the identification of a region of internuclear separation where charge transfer occurs.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2026

Fluid-kinetic modeling of a high power density radio frequency inductively coupled positive hydrogen ion source

High power density radio-frequency (RF) inductively coupled positive ion sources are attractive candidates for next-generation neutral beam injection (NBI) systems, where higher injected power and longer pulse lengths are desired without sacrificing source reliability. Operating at absorbed power densities of order $\gt 1~\mathrm{W\,cm}^{-3}$ places these sources in a regime with stronger gas heating, higher dissociation, and non-Maxwellian electron energy distributions. The Large Uniform Plasma for Ionizing Neutrals (LUPIN) is an RF inductively coupled plasma source designed to explore this high power density regime and to provide guidance for a positive ion source upgrade for the DIII-D NBI system. LUPIN is designed to operate at up to 20 kW of RF power at 2 MHz, coupling energy through a cylindrical quartz vessel to achieve target ion current densities of $2100\,\mathrm{A\,m}^{-2}$ . This paper presents fluid-kinetic modeling of LUPIN using the hybrid plasma equipment model where electrons are treated kinetically, and the simulations reveal that electron energy distribution function transitions from nearly Maxwellian in the core to bi-Maxwellian towards the edge. Parametric simulations investigate the effects of RF power, gas pressure, and frequency on plasma density, ion flux, and uniformity. Parametric sweeps reveal that increasing power shifts the primary ionization channel from molecular to atomic with diminishing flux gains due to skin-depth contraction and gas rarefaction. Higher frequency localizes heating and increases $\mathrm{H}_2^+$ and $\mathrm{H}_3^+$ delivery to the grid, while elevated pressure boosts ionization yet hinders ion transport due to increase in collisionality.

inductively coupled plasma↗

DOE OSTI2023

Radiolysis of Thin Water Ice in Electron Microscopy

Little is known about the radiolysis of water ice, especially as compared to the radiolysis of liquid water. In this study, the radiolytic decomposition of thin water ice films is probed for an initial 80 keV electron beam using electron energy loss spectroscopy, EELS. Pre-peaks in the oxygen K-edge spectra give the relative intensities of most of the oxygen containing species produced by radiolytic water decomposition. Contrary to expectations from the extrapolation of liquid water radiolysis data to high dose rates where significant H 2 O 2 production is expected, the main molecular product observed in the oxygen K-edge EELS spectra of water ice is O 2 . Significant mass loss of water and decomposition of most of its radiolytic products are observed for higher exposures, and the O atom seems to have a major role in the subsequent chemistry. An inverse relationship between the formation of O 2 and of ·OH is observed. A new high-dose reaction scheme is proposed. Furthermore, the significance of our results with respect to the radiolysis of water ice and for electron microscopy studies is discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

One-colour (~220 nm) resonance-enhanced (S 1 –S 0 ) multi-photon dissociation of acetylene: probe of the C 2 A 1 Π u – X 1 Σ + g band by frequency-modulation spectroscopy

In a recent paper, we demonstrated that one-colour (~220 nm), resonance-enhanced (S 1 –S 0 ), photodissociation of acetylene generates strong $C_2$ Swan band ($d^3Π_g–a^3Π_u$) and $C_2$ Deslandres-d'Azambuja band ($C^1Π_g–A^1Π_u$) fluorescence, and long-lived (>3 µs) fluorescence from an electronically-excited $C_2$H * species. It was not known whether the $C_2A^1Π_u$ and $X^1Σ^+_g$ states are also directly populated in this process. In this paper, multiple vibration-rotation transitions between the $C_2A$-state v = 2 and the X-state v = 0 level are examined by time-resolved frequency-modulation (FM) spectroscopy. The photolysis laser wavelength is tuned into resonance at the one-photon level with S 1 –S 0 transitions that populate individual rotational levels of the S 1 $_{trans}$-conformer $3^2,3^3,$ and $3^4$ vibrational states. By comparing the phase of the FM signals from the $C_2A–X$ transitions with that from the Rb D 1 -line absorption transition, we determine that, for all of the probed A–X transitions, the X-state level is more populated than the A-state level. We propose that the acetylene S 1 level is excited by the second photon to an acetylene dissociation precursor state, which undergoes sequential C-H bond-breaking to produce the $C_2X$ state. The dissociation precursor is assigned as the $1^1B_g(trans-bent)/1^1B_1(cis-bent)$ valence state, which correlates to a doubly-excited configuration, $(1π_u)^2(1π_g)^2$, at linear geometry. Based on the rotational distributions of the $C_2X$-state fragments, we believe that at least one of the transition states involved in the photolysis via S 1 $3^4$ has a larger CC-H bend-angle for the departing H-atom than that involved in the S 1 $3^2$ and $3^3$ photolysis.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2024

Diabatic Valence-Hole Concept

A global diabatization scheme, based on the “valence-hole” concept, has been previously applied to model webs of avoided crossings that exist in four electronic-state symmetry manifolds of C 2 ( 1 Π g , 3 Π g , 1 Σ u + , and 3 Σ u + ). Here, this model is extended to the electronically excited states of four more molecules: CN ( 2 Σ + ), N 2 ( 3 Π u ), SiC ( 3 Π), and Si 2 ( 3 Π g ). Many strangenesses in the spectroscopic observations (e.g., energy level structure, predissociation linewidths, and radiative lifetimes) for all four electronic state systems discussed here are accounted for by this unified model. The key concept of the model is valence-hole electron configurations: 3σ 2 4σ 1 1π 4 5σ 2 in CN ( 2 Σ + ), 2σ g 2 2σ u 1 1π u 4 3σ g 2 1π g 1 in N 2 ( 3 Π u ), 5σ 2 6σ 1 7σ 2 2π 3 in SiC ( 3 Π), and 4σ g 2 4σ u 1 5σ g 2 2π u 3 in Si 2 ( 3 Π g ), all of which have a triply occupied “valence-core” (i.e., 2σ g 2 2σ u 1 or the equivalent). These valence-hole configurations have a nominal bond order of three or higher and correlate with high-energy separated-atom limits with an np ← ns (n = 2, 3) promotion in one of the atomic constituents. On its way to dissociation, the strongly bound diabatic valence-hole state crosses multiple weakly bound or repulsive states, which belong to electron configurations with a completely filled valence-core. These curve crossings between diabatic potentials result in a network of many avoided crossings among multiple electronic states, analogous to the well-studied electronic structure landscape of ionic-covalent crossings in strongly ionic molecules. In conclusion, considering the unique role of valence-hole states in shaping the global electronic structure, the valence-hole concept should be added to our intuitive framework of chemical bonding.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Photodissociation and Theory to Investigate Uranium Oxide Cluster Cations

Uranium oxide cluster cations of the form U n O m + are produced by laser vaporization of a depleted uranium rod in a pulsed supersonic expansion. Ions are mass-analyzed and mass-selected with a time-of-flight spectrometer and studied with UV laser multiphoton dissociation. Cations of the stoichiometry UO 2 (UO 3 ) n + were observed as photofragments from all photodissociated cluster cations. (UO 3 ) n + clusters were also observed to result from dissociation of larger (UO 3 ) n + clusters, with UO 3 neutral as a common leaving group. Electronic structure calculations were used to investigate the stability of the prominent uranium oxide cluster cations using density functional theory (DFT) with the hybrid B3LYP exchange-correlation functional and at the CCSD(T) level and cc-pVnZ-PP basis sets (n = D,T), including diffuse orbitals as computational expense and availability permitted. Clustering energies, relative energies and dissociation energies of the cations are reported. The lowest energy neutral (UO 3 ) n clusters up to n = 3 are rings, n = 4 and 5 are chains with very low energy rings, and n = 6 is 3D. Finally, the lowest energy structures for UO 2 (UO 3 ) n + are composed of uranyl-like UO 2 + units bound by bridging oxygens to other UO 2 2+ units for n = 2 and 3, and for n = 4 a more complex 3D structure is predicted.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Revealing the nature of electron correlation in transition metal complexes with symmetry breaking and chemical intuition

Here, in this work, we provide a nuanced view of electron correlation in the context of transition metal complexes, reconciling computational characterization via spin and spatial symmetry breaking in single-reference methods with qualitative concepts from ligand-field and molecular orbital theories. These insights provide the tools to reliably diagnose the multi-reference character, and our analysis reveals that while strong (i.e., static) correlation can be found in linear molecules (e.g., diatomics) and weakly bound and antiferromagnetically coupled (monometal-noninnocent ligand or multi-metal) complexes, it is rarely found in the ground-states of mono-transition-metal complexes. This leads to a picture of static correlation that is no more complex for transition metals than it is, e.g., for organic biradicaloids. In contrast, the ability of organometallic species to form more complex interactions, involving both ligand-to-metal σ-donation and metal-to-ligand π-backdonation, places a larger burden on a theory's treatment of dynamic correlation. We hypothesize that chemical bonds in which inter-electron pair correlation is non-negligible cannot be adequately described by theories using MP2 correlation energies and indeed find large errors vs experiment for carbonyl-dissociation energies from double-hybrid density functionals. A theory's description of dynamic correlation (and to a less important extent, delocalization error), which affects relative spin-state energetics and thus spin symmetry breaking, is found to govern the efficacy of its use to diagnose static correlation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Characterization of plasma catalytic decomposition of methane: role of atomic O and reaction mechanism

In this work, we investigated atmospheric pressure plasma jet (APPJ)-assisted methane oxidation over a Ni-SiO 2 /Al 2 O 3 catalyst. We evaluated possible reaction mechanisms by analyzing the correlation of gas phase, surface and plasma-produced species. Plasma feed gas compositions, plasma powers, and catalyst temperatures were varied to expand the experimental parameters. Real-time Fourier-transform infrared spectroscopy was applied to quantify gas phase species from the reactions. The reactive incident fluxes generated by plasma were measured by molecular beam mass spectroscopy using an identical APPJ operating at the same conditions. A strong correlation of the quantified fluxes of plasma-produced atomic oxygen with that of CH 4 consumption, and CO and CO 2 formation implies that O atoms play an essential role in CH 4 oxidation for the investigated conditions. With the integration of APPJ, the apparent activation energy was lowered and a synergistic effect of 30% was observed. We also performed in-situ diffuse reflectance infrared Fourier-transform spectroscopy to analyze the catalyst surface. The surface analysis showed that surface CO abundance mirrored the surface coverage of CH n at 25 °C. This suggests that CH n adsorbed on the catalyst surface as an intermediate species that was subsequently transformed into surface CO. We observed very little surface CH n absorbance at 500 °C, while a ten-fold increase of surface CO and stronger CO 2 absorption were seen. This indicates that for a nickel catalyst at 500 °C, the dissociation of CH 4 to CHn may be the rate-determining step in the plasma-assisted CH 4 oxidation for our conditions. We also found the CO vibrational frequency changes from 2143 cm –1 for gas phase CO to 2196 cm –1 for CO on a 25 °C catalyst surface, whereas the frequency of CO on a 500 °C catalyst was 2188 cm –1 . Here, the change in CO vibrational frequency may be related to the oxidation of the catalyst.

methane oxidation↗

DOE OSTI2024

Atomistic simulations for investigation of substrate and salt effects on lipid in-source fragmentation in secondary ion mass spectrometry: A follow-up study

In-source fragmentation (ISF) poses a significant challenge in secondary ion mass spectrometry (SIMS). These fragment ions increase the spectral complexity and can lead to incorrect annotation of fragments as intact species. The presence of salt that is ubiquitous in biological samples can influence the fragmentation and ionization of analytes in a significant manner, but their influences on SIMS have not been well characterized. To elucidate the effect of substrates and salt on ISF in SIMS, we have employed experimental SIMS in combination with atomistic simulations of a sphingolipid on a gold surface with various NaCl concentrations as a model system. Our results revealed that a combination of bond dissociation energy and binding energy between N-palmitoyl-sphingomyelin and a gold surface is a good predictor of fragment ion intensities in the absence of salt. However, ion-fragment interactions play a significant role in determining fragment yields in the presence of salt. Additionally, the charge distribution on fragment species may be a major contributor to the varying effects of salt on fragmentation. This study demonstrates that atomistic modeling can help predict ionization potential when salts are present, providing insights for more accurate interpretations of complex biological spectra.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2022

Observation of site-selective chemical bond changes via ultrafast chemical shifts

The concomitant motion of electrons and nuclei on the femtosecond time scale marks the fate of chemical and biological processes. Here we demonstrate the ability to initiate and track the ultrafast electron rearrangement and chemical bond breaking site-specifically in real time for the carbon monoxide diatomic molecule. We employ a local resonant x-ray pump at the oxygen atom and probe the chemical shifts of the carbon core-electron binding energy. We observe charge redistribution accompanying core-excitation followed by Auger decay, eventually leading to dissociation and hole trapping at one site of the molecule. The presented technique is general in nature with sensitivity to chemical environment changes including transient electronic excited state dynamics. This work provides a route to investigate energy and charge transport processes in more complex systems by tracking selective chemical bond changes on their natural timescale.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS2003

Towards the Structure Determination of a Modulated Protein Crystal: The Semicrystalline State of Profilin:Actin

One of the remaining challenges to structural biology is the solution of modulated structures. While small molecule crystallographers have championed this type of structure, to date, no modulated macromolecular structures have been determined. Modulation of the molecular structures within the crystal can produce satellite reflections or a superlattice of reflections in reciprocal space. We have developed the data collection methods and strategies that are needed to collect and analyze these data. If the macromolecule's crystal lattice is composed of physiologically relevant packing contacts, structural changes induced under physiological conditions can cause distortion relevant to the function and biophysical processes of the molecule making up the crystal. By careful measurement of the distortion, and the corresponding three-dimensional structure of the distorted molecule, we will visualize the motion and mechanism of the biological macromolecule(s). We have measured the modulated diffraction pattern produced by the semicrystalline state of profilin:actin crystals using highly parallel and highly monochromatic synchrotron radiation coupled with fine phi slicing (0.001-0.010 degrees) for structure determination. These crystals present these crystals present a unique opportunity to address an important question in structural biology. The modulation is believed to be due to the formation of actin helical filaments from the actin beta ribbon upon the pH-induced dissociation of profilin. To date, the filamentous state of actin has resisted crystallization and no detailed structures are available. The semicrystalline state profilin:actin crystals provides a unique opportunity to understand the many conformational states of actin. This knowledge is essential for understanding the dynamics underlying shape changes and motility of eukaryotic cells. Many essential processes, such as cytokinesis, phagocytosis, and cellular migration depend upon the capacity of the actin microfilament system to be restructured in a controlled manner via polymerization, depolymerization, severing, cross-linking, and anchorage. The structure the semicrystalline state of profilin:actin will challenge and validate current models of muscle contraction and cell motility. The methodology and theory under development will be easily extendable to other systems.

Borgstahl, G.↗

DOE OSTI2023

Dissociation dynamics of ion-pair states accessed by low-energy electron collisions

Abstract Ion-pair (IP) states are the superexcited states of the neutral molecule that dissociate as a pair of positive and negative ions. These states are generally present near the ionization continuum of the molecule. IP states can be accessed by using photoexcitation or by using the electron collision technique. Different photoexcitation techniques are widely used over the years to study the threshold (threshold IP production spectroscopy) and the dynamics (IP imaging spectroscopy) of the IP states. However, the electron collision technique is ignored over the years and only a few studies are available. In this review, we will discuss different experimental techniques to probe IP states by using electron collision and also the dynamics of the IP states that are accessed by the electron collision.

74 ATOMIC AND MOLECULAR PHYSICS↗

NASA NTRS1995

A far-infrared spectral line survey of 23 infrared-bright Galaxy nuclei

We present results from a KAO survey of fine-structure lines observed in 23 infrared-luminous galaxies. One or more of the following lines was observed and/or detected in each galaxy: (S III) 19, 33 microns, (Ne V) 24 microns, (O IV) 26 microns, (Fe II) 26 microns, (Si II) 35 microns, (O III) 52, 88 microns, (0 I) 63, 146 microns, (N III) 57 micro ns, (N II) 122, 205 microns, (C II) 158 microns. The galaxies span a wide range of morphologies (irregular to grand design), have varying metallicities, and include mergers, AGN's, and starburst systems. The observations were made beginning in 1988 using the facility Cryogenic Grating Spectrometer onboard the KAO at a typical resolution of approximately 60-140 km/s and with a 30-44 deg beam. We interpret the (C II) and (O I) fluxes, along with previous measurements of the IR continuum fluxes, in the context of photo dissociation region (PDR) models (Tielens & Hollenbach 1985; Wolfire et al. 1990). With these models, we obtain estimates of the typical interstellar UV fields incident on the line emitting regions (102-104 times the local interstellar radiation field) and the total masses (10(exp 7)-10(exp 8) Solar Mass), densities (10(exp 3)-10(exp 4)/cu cm), and temperatures (100-250 K) of the warm atomic gas. The (O III) (52/88) and (S III) (33/19) line flux ratios constrain the range of electron densities and pressures found within the ionized regions. The (O III) and (S III) lines also provide estimates of the effective temperature of the ionizing stars and elemental abundances within the ionized regions of these galactic nuclei. Our measurements imply typical gas pressures of nT approximately 5 x 10(exp 6)/cu cm K and typical upper mass cutoffs of 25-35 Solar Mass. The low-metallicity systems show high (C II)/CO and (O I)/CO flux ratios, 3-5 times the Milky Way value, indicating that they contain a larger fraction of photodissociated gas relative to the molecular material.

Lord, Steven D.↗

DOE OSTI2023

Feasibility of a Novel Density Functional Method Outside the Kohn-Sham Framework for Modeling Global Potential Energy Surfaces of Molecular Chemical Reactions (Final Technical Report)

The project aimed towards the construction of computational methods capable of modeling the global potential energy surface (PES) of small molecules, molecular ions, and radicals — including the parts of the PES which correspond to chemical reactions, and the reaction paths connecting reactants, intermediates, and reaction products. This goal may seem humble at first glance, but for chemical systems with more than about six atoms in total, at the time the project was proposed, there were no established theoretical methods capable of simulating such systems reliably, not even for small molecules in the gas phase which are electronically benign. This restriction severely hampers our ability to model and control chemical processes under harsh conditions. With the goal of constructing a method capable of modeling such chemical systems, we proposed to pursue a novel approach towards a Multi-Configuration (MC) DFT outside the traditional frameworks of Kohn-Sham theory and other methods of coupling wave function theory with DFT. Rather than being a complete active space (CAS) method, the proposed DOCI-DFT would employ a special restricted form of the active space wave function, called Doubly-Occupied CI (DOCI)—this wave function form is sufficient to describe not only heterolytic, but also homolytic bond dissociation processes at the zeroth order (i.e., as active space wave function). There are no other standard mean field methods which can do so. In the original proposal, we also outlined strong formal and practical arguments speaking for this method. The project proposal was accepted by the Department of Energy and provided two years of seed funding for one graduate student, as well as two years of two weeks PI summer-salary for the PI. However, ultimately the project could not be effectively pursued due to severe interferences outside the context scientific problems and was cancelled. Apart from two weeks of PI salary (and associated fringe & overhead costs issues by the performing institution), no costs were charged against the award, with the entire rest of funding returned to the Department of Energy. This report reiterates the primary background information regarding the project, its original goals, and outlines the preliminary work performed during the two weeks of DOE funding charged.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Potentials of mean force fail to describe chemical bond-breaking in solution

Many liquid phase studies assume that the potential energy surfaces of reacting molecules are the same as in the gas phase, neglecting complex solvent dynamics that can completely alter the nature of chemical reactivity. Even studies that include solvent effects typically only consider them in an average, equilibrium way as part of a potential of mean force (PMF). In this work, we use mixed quantum/classical simulations to compare how equilibrium and non-equilibrium solvent motions affect the photodissociation of a simple diatomic molecule, NaK + , in liquid tetrahydrofuran. A PMF analysis shows that as the excited-state molecule dissociates with the solvent at equilibrium, the bonding electron remains associated with K + at short bond distances but eventually localizes on Na + at the end of dissociation. When we examine non-equilibrium dynamical photodissociation trajectories, however, we find that they fall into three distinct categories: about a quarter of them have the bonding electron mainly associated with Na + , another quarter stay mainly associated with K + , and about half have the bonding electron shared roughly equally between the two ions. The results show that equilibrium PMFs cannot accurately describe the dynamics of bond-breaking chemical reactions in solution because there is insufficient time for the solvent to reach equilibrium on the time scale over which bond dissociation occurs. Furthermore, our analysis shows that the solvent coupling between the electronic energy surfaces is similar at and away from equilibrium, suggesting that other factors, such as solute velocity-driven solvent memory effects, play a more important role in explaining the failure of the equilibrium PMF to predict the non-equilibrium dynamics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗