The Role of the Solvent in the Condensed-Phase Dynamics and Identity of Chemical Bonds: The Case of the Sodium Dimer Cation in THF

How Solvation Alters the Thermodynamics of Asymmetric Bond-Breaking: Quantum Simulation of NaK+ in Liquid Tetrahydrofuran

Gas-phase potential energy surfaces (PESs) are often used to provide an intuitive understanding of molecular chemical reactivity. Most chemical reactions, however, take place in solution, and it is unclear whether gas-phase PESs accurately represent chemical processes in solvent environments. In this work we use quantum simulations to investigate the dissociation energetics of NaK+ in liquid tetrahydrofuran (THF) to understand the degree to which solvent interactions alter the gas-phase picture. Using umbrella sampling and thermodynamic integration techniques, we construct condensed-phase free energy surfaces of NaK+ on THF in both the ground and electronic excited states. We find that solvation by THF completely alters the nature of the NaK+ bond by reordering the thermodynamic dissociation products. Reaching the thermodynamic dissociation limit in THF also requires a long-range charge transfer process that has no counterpart in the gas phase. Gas-phase PESs, even with perturbations, cannot adequately describe the reactivity of simple asymmetric molecules in solution.

Using Machine Learning to Understand the Causes of Quantum Decoherence in Solution-Phase Bond-Breaking Reactions

Decoherence is a fundamental phenomenon that occurs when an entangled quantum state interacts with its environment, leading to collapse of the wave function. The inevitability of decoherence provides one of the most intrinsic limits of quantum computing. However, there has been little study of the precise chemical motions from the environment that cause decoherence. Here, we use quantum molecular dynamics simulations to explore the photodissociation of Na2+ in liquid Ar, in which solvent fluctuations induce decoherence and thus determine the products of chemical bond breaking. We use machine learning to characterize the solute-solvent environment as a high-dimensional feature space that allows us to predict when and onto which photofragment the bonding electron will localize. We find that reaching a requisite photofragment separation and experiencing out-of-phase solvent collisions underlie decoherence during chemical bond breaking. Our work highlights the utility of machine learning for interpreting complex solution-phase chemical processes as well as identifies the molecular underpinnings of decoherence.

Electron Presolvation in Tetrahydrofuran-Incorporated Supramolecular Sodium Entities

Alkali metal atoms can repopulate their valence electrons toward solvation due to impact from solvents or microsurroundings and provide the remaining alkali metal cations for coordinating with a variety of specific solvents, forming various electron-expanded complexes or solvated ionic pairs with special interactions. Such special solute-solvent interactions not only affect their electronic structures but also enable the formation of entirely new species. Taking Na(THF)_n (n = 1-6, THF = tetrahydrofuran) and Na₂@THF complexes as typical representatives, density functional theory calculations are carried out to explore the solvation of a sodium atom and its dimer in THF and characterize their complexes as solvent-incorporated supramolecular entities and particularly valence electron presolvation due to their interaction with solvent THF. Electron presolvation is caused by the Pauli repulsion between THF containing a coordinating O atom with a lone pair of electrons and the alkali metal Na or Na₂ containing valence electrons, and THF coordination to them forces their valence electrons to redistribute, which can be easily realized in such solvents. Compared with strongly bound valance electrons of alkali metal atoms, THF coordination enables Na or Na₂ electrons to exhibit much more active states (i.e., the presolvated states) featuring small vertical detachment energies of electrons and distorted diffuse distributions in the frames of the generally structured metal cation complexes, acting as the electron-expanded chemical entities. Furthermore, the degree of electron diffusion and the polarity of the Na-Na bond are proportional to the coordination number (n) and the coordination number difference (Δn) between two Na centers in Na₂@THF. The unique properties of such entities are also discussed. This work offers a theoretical support to the supramolecular entities formed by alkali-metal atoms or their dimers with ligands containing O or N and uncovers the unique electron presolvation phenomena and also enriches our understanding of the novel metal atom complexes.

Solvent Control of Chemical Identity Can Change Photodissociation into Photoisomerization

In solution-phase chemistry, the solvent is often considered to be merely a medium that allows reacting solutes to encounter each other. In this work, however, we show that moderate locally specific solute-solvent interactions can affect not only the nature of the solute but also the types of reactive chemistry. We use quantum simulation methods to explore how solvent participation in solute chemical identity alters reactions involving the breaking of chemical bonds. In particular, we explore the photoexcitation dynamics of Na₂⁺ dissolved in liquid tetrahydrofuran. In the gas phase, excitation of Na₂⁺ directly leads to dissociation, but in solution, photoexcitation leads to an isomerization reaction involving rearrangement of the first-shell solvent molecules; this isomerization must go to completion before the solute can dissociate. Despite the complexity, the solution-phase reaction dynamics can be captured by a two-dimensional energy surface where one dimension involves only the isomerization of the first-shell solvent molecules.

Bond-Breaking Reactions Encounter Distinct Solvent Environments Causing Breakdown of Linear Response

Solvent effects are important for understanding solution-phase chemical reactions. Surprisingly, very few studies have explored how solvent dynamics change during the course of a reaction with solutes that encounter a wide range of configurations. Here, we use quantum simulation methods to explore the solvent dynamics during a solution-phase bond-breaking reaction: the photodissociation of Na₂⁺ in liquid Ar. We find that the solute experiences a small number of distinct solvent environments that change in a discrete fashion as the bond lengthens. In characterizing the solvent environments, we show also that linear response fails by all measures, even when nonstationarity of solvent dynamics is considered. This observation of distinct solvent response environments with a solvent that can undergo only translational motions highlights the complexity of solute-solvent interactions, but that there are only a few environments gives hope to the idea that solvation dynamics can be understood for solution-phase reactions that explore a wide configuration space.

How Water-Ion Interactions Control the Formation of Hydrated Electron:Sodium Cation Contact Pairs

Although solvated electrons are a perennial subject of interest, relatively little attention has been paid to the way they behave in aqueous electrolytes. Experimentally, it is known that the hydrated electron's (e_aq^-) absorption spectrum shifts to the blue in the presence of salts, and the magnitude of the shift depends on the ion concentration and the identities of both the cation and anion. Does the blue-shift result from some type of dielectric effect from the bulk electrolyte, or are there specific interactions between the hydrated electron and ions in solution? Previous work has suggested that e_aq^- forms contact pairs with aqueous ions such as Na⁺, leading to the question of what controls the stability of such contact pairs and their possible connection to the observed spectroscopy. In this work, we use mixed quantum/classical simulations to examine the nature of Na⁺:e^- contact pairs in water, using a novel method for quantum umbrella sampling to construct e_aq^--ion potentials of mean force (PMF). We find that the nature of the contact pair PMF depends sensitively on the choice of the classical interactions used to describe the Na⁺-water interactions. When the ion-water interactions are slightly stronger, the corresponding cation:e^- contact pairs form at longer distances and become free energetically less stable. We show that this is because there is a delicate balance between solvation of the cation, solvation of e_aq^- and the direct electronic interaction between the cation and the electron, so that small changes in this balance lead to large changes in the formation and stability of e^--ion contact pairs. In particular, strengthening the ion-water interactions helps to maintain a favorable local solvation environment around Na⁺, which in turn forces water molecules in the first solvation shell of the cation to be unfavorably oriented toward the electron in a contact pair; stronger solvation of the cation also reduces the electronic overlap of e_aq^- with Na⁺. We also find that the calculated spectra of different models of Na⁺:e^- contact pairs do not shift monotonically with cation-electron distance, and that the calculated spectral shifts are about an order of magnitude larger than experiment, suggesting that isolated contact pairs are not the sole explanation for the blue-shift of the hydrated electron's spectrum in the presence of electrolytes.

Nonequilibrium Solvent Effects during Photodissociation in Liquids: Dynamical Energy Surfaces, Caging, and Chemical Identity

In the gas phase, potential energy surfaces can be used to provide insight into the details of photochemical reaction dynamics. In solution, however, it is unclear what potential energy surfaces, if any, can be used to describe even simple chemical reactions such as the photodissociation of a diatomic solute. In this paper, we use mixed quantum/classical (MQC) molecular dynamics (MD) to study the photodissociation of Na₂⁺ in both liquid Ar and liquid tetrahydrofuran (THF). We examine both the gas-phase potential surfaces and potentials of mean force (PMF), which assume that the solvent remains at equilibrium with the solute throughout the photodissociation process and show that neither resemble a nonequilibrium dynamical energy surface that is generated by taking the time integral of work. For the photodissociation of Na₂⁺ in liquid Ar, the dynamical energy surface shows clear signatures of solvent caging, and the degree of caging is directly related to the mass of the solvent atoms. For Na₂⁺ in liquid THF, local specific interactions between the solute and solvent lead to changes in chemical identity that create a kinetic trap that effectively prevents the molecule from dissociating. The results show that nonequilibrium effects play an important role even in simple solution-phase reactions, requiring the use of dynamical energy surface to understand such chemical events.