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

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.

State averaged CASSCF in AMOEBA polarizable water model for simulating nonadiabatic molecular dynamics with nonequilibrium solvation effects

This paper presents a state-averaged complete active space self-consistent field (SA-CASSCF) in the atomic multipole optimized energetics for biomolecular application (AMOEBA) polarizable water model, which enables rigorous simulation of non-adiabatic molecular dynamics with nonequilibrium solvation effects. The molecular orbital and configuration interaction coefficients of the solute wavefunction, and the induced dipoles on solvent atoms, are solved by minimizing the state averaged energy variationally. In particular, by formulating AMOEBA water models and the polarizable continuum model (PCM) in a unified way, the algorithms developed for computing SA-CASSCF/PCM energies, analytical gradients, and non-adiabatic couplings in our previous work can be generalized to SA-CASSCF/AMOEBA by properly substituting a specific list of variables. Implementation of this method will be discussed with the emphasis on how the calculations of different terms are partitioned between the quantum chemistry and molecular mechanics codes. We will present and discuss results that demonstrate the accuracy and performance of the implementation. Next, we will discuss results that compare three solvent models that work with SA-CASSCF, i.e., PCM, fixed-charge force fields, and the newly implemented AMOEBA. Finally, the new SA-CASSCF/AMOEBA method has been interfaced with the ab initio multiple spawning method to carry out non-adiabatic molecular dynamics simulations. This method is demonstrated by simulating the photodynamics of the model retinal protonated Schiff base molecule in water.

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.