Non-adiabatic molecular dynamics of excited Na2+ solvated in Ar17 clusters

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.

Structure, energetics, and spectroscopy of the K2+(X2Σ+g) interacting with the noble gas atoms Ar, Kr and Xe

The structure, energetic, and spectroscopy properties of the ionic system K₂⁺(X²Σ⁺_g) interacting with the noble gas atoms Argon, Krypton and Xenon are studied. The computations are done by an accurate ab initio approach based on the pseudo-potential technique, Gaussian basis sets, parameterized l-dependent polarization potentials and an analytic potential form for the K⁺Ar, K⁺Kr and K⁺Xe interactions. These interactions are added via the CCSD(T) potential taken from literature and fitted applying the analytical expression of Tang and Toennies. The application of the pseudo-potential approach reduces the number of active electrons of each complex to only one electron. The potential energy surfaces are analyzed on a large range of the Jacobi coordinates, R and θ. By the general interpolation outline based on the RKHS (Reproducing Kernel Hilbert Space) procedure, we have reproduced for each complex from our ab initio results the two-dimensional contour plots of an analytical potential. To evaluate the stability of each complex, we have determined from the potential energy surfaces the equilibrium distance (R_e), the well depth (D_e), the quantum energy (D₀), the zero-point-energy (ZPE) and the ZPE%. The results showed that the linear configurations, where the noble gas atom connected to the K₂⁺(X²Σ⁺_g) system, are the more stable.

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.

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.

Spectroscopy and Dynamics of K Atoms on Argon Clusters

We present a combined experimental and simulation study of the 4s → 4p photoexcitation of the K atom trapped at the surface of ArN clusters made of a few hundred Ar atoms. Our experimental method based on photoelectron spectroscopy allows us to firmly establish that one single K atom is trapped at the surface of the cluster. The absorption spectrum is characterized by the splitting of the atomic absorption line into two broad bands, a Π band associated with p orbitals parallel to the cluster surface and a Σ band associated with the perpendicular orientation. The spectrum is consistent with observations reported for K atoms trapped on lighter inert gas clusters, but the splitting between the Π and Σ bands is significantly larger. We show that a large amount of K atoms are transiently stuck and eventually lost by the Ar cluster, in contrast with previous observations reported for alkaline earth metal systems. The excitation in the Σ band leads systematically to the ejection of the K atom from the Ar cluster. On the contrary, excitation in the Π band leads to the formation of a bound state. In this case, the analysis of the experimental photoelectron spectrum by means of nonadiabatic molecular dynamics simulation shows that the relaxation drives the system toward a basin where the coordination of the K atom is 2.2 Ar atoms on the average, in a poorly structured surface.

Nonadiabatic molecular dynamics of photoexcited Li2(+) Ne(n) clusters

We investigate the relaxation of photoexcited Li(2)(+) chromophores solvated in Ne(n) clusters (n = 2-22) by means of molecular dynamics with surface hopping. The simplicity of the electronic structure of these ideal systems is exploited to design an accurate and computationally efficient model. These systems present two series of conical intersections between the states correlated with the Li+Li(2s) and Li+Li(2p) dissociation limits of the Li(2)(+) molecule. Frank-Condon transition from the ground state to one of the three lowest excited states, hereafter indexed by ascending energy from 1 to 3, quickly drives the system toward the first series of conical intersections, which have a tremendous influence on the issue of the dynamics. The states 1 and 2, which originate in the Frank-Condon area from the degenerated nondissociative 1(2)Π(u) states of the bare Li(2)(+) molecule, relax mainly to Li+Li(2s) with a complete atomization of the clusters in the whole range of size n investigated here. The third state, which originates in the Frank-Condon area from the dissociative 1(2)Σ(u)(+) state of the bare Li(2)(+) molecule, exhibits a richer relaxation dynamics. Contrary to intuition, excitation into state 3 leads to less molecular dissociation, though the amount of energy deposited in the cluster by the excitation process is larger than for excitation into state 1 and 2. This extra amount of energy allows the system to reach the second series of conical intersections so that approximately 20% of the clusters are stabilized in the 2(2)Σ(g)(+) state potential well for cluster sizes n larger than 6.