Nonadiabatic Trajectory Studies of NaI(H2O)n Photodissociation Dynamics

Solvation of magnesium chloride dimer in water: The case of anionic and neutral clusters

The structures of magnesium chloride dimer-water clusters, (MgCl2)2(H2O)n-/0, were investigated with size-selected anion photoelectron spectroscopy and theoretical calculations to understand the dissolution of magnesium chloride in water. The most stable structures were confirmed by comparing vertical detachment energies (VDEs) with the experimental measurements. A dramatic drop of VDE at n = 3 has been observed in the experiment, which is in accordance with the structural change of (MgCl2)2(H2O)n-. Compared to the neutral clusters, the excess electron induces two significant phenomena in (MgCl2)2(H2O)n-. First, the planar D2h geometry can be converted into a C3v structure at n = 0, making the Mg-Cl bonds easier to be broken by water molecules. More importantly, a negative charge-transfer-to-solvent process occurs after adding three water molecules (i.e., at n = 3), which leads to an obvious deviation in the evolution of the clusters. Such electron transfer behavior was noticed at n = 1 in monomer MgCl2(H2O)n-, indicating that the dimerization between two MgCl2 molecules can make the cluster more capable of binding electron. In neutral (MgCl2)2(H2O)n, this dimerization provides more sites for the added water molecules, which can stabilize the entire cluster and maintain its initial structure. Specifically, filling the coordination number to be 6 for Mg atoms can be seen as a link between structural preferences in the dissolution of the monomers, dimers, and extended bulk-state of MgCl2. This work represents an important step forward into fully understanding the solvation of MgCl2 crystals and other multivalent salt oligomers.

Comparison of the Microsolvation of CaX2 (X = F, Cl, Br, I) in Water: Size-Selected Anion Photoelectron Spectroscopy and Theoretical Calculations

To understand the microsolvation of alkaline-earth dihalides in water and provide information about the dependence of solvation processes on different halides, we investigated CaBr₂(H₂O)_n^-, CaI₂(H₂O)_n^-, and CaF₂(H₂O)_n^- (n = 0-6) clusters using size-selected anion photoelectron spectroscopy and conducted theoretical calculations on these clusters and their neutrals. The results are compared with those of CaCl₂(H₂O)_n^-/0 clusters reported previously. It is found that the vertical detachment energies (VDEs) of CaCl₂(H₂O)_n^-, CaBr₂(H₂O)_n^-, and CaI₂(H₂O)_n^- show a similar trend with increasing cluster size, while the VDEs of CaF₂(H₂O)_n^- show a different trend. The VDEs of CaF₂(H₂O)_n^- are much lower than those of CaCl₂(H₂O)_n^-, CaBr₂(H₂O)_n^-, and CaI₂(H₂O)_n^-. A detailed probing of the structures shows that a significant increase of the Ca-X distance (separation of Ca²⁺-X^- ion pair) in CaCl₂(H₂O)_n^-/0, CaBr₂(H₂O)_n^-/0, and CaI₂(H₂O)_n^-/0 clusters occurred at about n = 5. However, for CaF₂(H₂O)_n^-/0, no abrupt change of the Ca-F distance with the increasing cluster size has been observed. In CaCl₂(H₂O)₆^-/0, CaBr₂(H₂O)₆^-/0, and CaI₂(H₂O)₆^-/0, the Ca atom coordinates directly with 5 H₂O molecules. However, in CaF₂(H₂O)_n^-/0, the Ca atom coordinates directly with only 2 or 3 H₂O molecules. The similarity or differences in the structures and coordination numbers are consistent with the fact that CaCl₂, CaBr₂, and CaI₂ have similar solubility, while CaF₂ has much lower solubility.

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.

Hydration processes of barium chloride: Size-selected anion photoelectron spectroscopy and theoretical calculations of BaCl2-water clusters

In order to understand the hydration processes of BaCl₂, we investigated BaCl₂(H₂O)_n ^- (n = 0-5) clusters using size-selected anion photoelectron spectroscopy and theoretical calculations. The structures of neutral BaCl₂(H₂O)_n clusters up to n = 8 were also investigated by theoretical calculations. It is found that in BaCl₂(H₂O)_n ^-/0, the Ba-Cl distances increase very slowly with the cluster size. The hydration process is not able to induce the breaking of a Ba-Cl bond in the cluster size range (n = 0-8) studied in this work. In small BaCl₂(H₂O)_n clusters with n ≤ 5, the Ba atom has a coordination number of n + 2; however, in BaCl₂(H₂O)_6-8 clusters, the Ba atom coordinates with two Cl atoms and (n - 1) water molecules, and it has a coordination number of n + 1. Unlike the previously studied MgCl₂(H₂O)_n ^- and CaCl₂(H₂O)_n ^-, negative charge-transfer-to-solvent behavior has not been observed for BaCl₂(H₂O)_n ^-, and the excess electron of BaCl₂(H₂O)_n ^- is mainly localized on the Ba atom rather on the water molecules. No observation of Ba²⁺-Cl^- separation in current work is consistent with the lower solubility of BaCl₂ compared to MgCl₂ and CaCl₂. Considering the BaCl₂/H₂O mole ratio in the saturated solution, one would expect that about 20-30 H₂O molecules are needed to break the first Ba-Cl bond in BaCl₂.

Control of photodissociation of the NaI molecule via pulse chirping

A dynamic pump control scheme is proposed to manipulate the predissociation process of NaI molecules in different reaction channels. A linearly chirped pulse is used to excite the NaI molecule, and a time-delayed infrared pulse is employed to modify the molecular potentials in the coupling zone. The predissociation branching ratio of the product from two channels can be controlled by tuning the chirp rate with a proper range of delay times. Furthermore, an additional ultrafast photoionization step is adopted to monitor the wave packet evolution and probe the possible modifications of the electronic potential under the influence of a chirped pump field to reveal the physical mechanism behind the control. Aulter-Townes splitting is observed at a proper chirp rate, and the dressed-state population can be controlled via pulse chirping.

Initial hydration behavior of sodium iodide dimer: photoelectron spectroscopy and ab initio calculations

We investigated (NaI)₂⁻(H₂O)_n (n = 0–6) clusters to examine the initial solvation process of (NaI)₂ in water, using negative ion photoelectron spectroscopy and theoretical calculations.

On the dissolution of lithium sulfate in water: anion photoelectron spectroscopy and density functional theory calculations

The initial dissolution steps of lithium sulfate (Li2SO4) in water were investigated by performing anion photoelectron spectroscopy and density functional theory calculations on the Li2SO4(H2O)n(-) (n = 0-5) clusters. The plausible structures of these clusters and the corresponding neutral clusters were obtained using LC-ωPBE/6-311++G(d,p) calculations by comparing the experimental and theoretical vertical electron detachment energies. Two types of structures for bare Li2SO4(-/0) were found: a turtle-shaped structure and a propeller-shaped structure. For Li2SO4(H2O)n(-) cluster anions with n = 1-3, two kinds of isomers derived from the turtle-shaped and propeller-shaped structures of bare Li2SO4(-) were identified. For n = 4-5, these two kinds of isomers present similar structural and energetic features and thus are not distinguishable. For the anionic clusters the water molecules prefer to firstly interact with one Li atom until fully coordinating it. While for the neutral clusters, the water molecules interact with the two Li atoms alternately, therefore, showing a pairwise solvation behavior. The Li-S distance increases smoothly upon addition of water molecules one by one. Addition of five water molecules to Li2SO4 cannot induce the dissociation of one Li(+) ion because the water molecules are shared by two Li(+) ions.

QM/MM investigations of organic chemistry oriented questions

About 35 years after its first suggestion, QM/MM became the standard theoretical approach to investigate enzymatic structures and processes. The success is due to the ability of QM/MM to provide an accurate atomistic picture of enzymes and related processes. This picture can even be turned into a movie if nuclei-dynamics is taken into account to describe enzymatic processes. In the field of organic chemistry, QM/MM methods are used to a much lesser extent although almost all relevant processes happen in condensed matter or are influenced by complicated interactions between substrate and catalyst. There is less importance for theoretical organic chemistry since the influence of nonpolar solvents is rather weak and the effect of polar solvents can often be accurately described by continuum approaches. Catalytic processes (homogeneous and heterogeneous) can often be reduced to truncated model systems, which are so small that pure quantum-mechanical approaches can be employed. However, since QM/MM becomes more and more efficient due to the success in software and hardware developments, it is more and more used in theoretical organic chemistry to study effects which result from the molecular nature of the environment. It is shown by many examples discussed in this review that the influence can be tremendous, even for nonpolar reactions. The importance of environmental effects in theoretical spectroscopy was already known. Due to its benefits, QM/MM can be expected to experience ongoing growth for the next decade.In the present chapter we give an overview of QM/MM developments and their importance in theoretical organic chemistry, and review applications which give impressions of the possibilities and the importance of the relevant effects. Since there is already a bunch of excellent reviews dealing with QM/MM, we will discuss fundamental ingredients and developments of QM/MM very briefly with a focus on very recent progress. For the applications we follow a similar strategy.

Importance of polarization in quantum mechanics/molecular mechanics descriptions of electronic excited states: NaI(H2O)n photodissociation dynamics as a case study

Sodium iodide has long been a paradigm for ionic and covalent curve crossing and ultrafast nonadiabatic dynamics, and our interest lies in the influence of solvation on this process. The NaI(H2O)n photodissociation dynamics are simulated with the molecular dynamics with quantum transitions method. A quantum mechanics/molecular mechanics (QM/MM) description is adopted for the NaI(H2O)n electronic states, in which a semiempirical valence bond approach is used to describe the NaI electronic structure, and a polarizable optimized potential for cluster simulations model is used to describe solute-solvent and solvent-solvent interactions. In contrast to previous work with a nonpolarizable MM model [Koch et al., J. Phys. Chem. A, 2006, 110, 1438], this approach predicts that the NaI ionic ground- to covalent first-excited-state Franck-Condon energy gaps reach a plateau by cluster size 16, in relatively good agreement with experiment and electronic structure calculations; this allows us to safely extend our previous simulations to larger cluster sizes, i.e., n > 4. The simulations suggest that the disappearance of the two-photon ionization probe signals observed in femtosecond pump-probe experiments of NaI(H2O)n, n >/= 4, is due to the shift of the NaI curve-crossing region toward larger NaI internuclear separations because of solvent stabilization of the NaI ionic state. Further, the latter causes the adiabatic ground and excited states to acquire pure ionic and covalent character, respectively, by cluster 8, resulting in NaI ionic ground-state recombination or dissociation. To make a connection with electron transfer in solution, free energy curves have been generated as a function of a solvent coordinate similar to that of solution theory. Inspection of the free energy curves together with the results of excited-state simulations reveal that the electron-transfer process in clusters is not governed by the collective motion of the solvent molecules, as in solution, but that it rather proceeds along the NaI internuclear separation coordinate, as in the gas phase. In fact, solvation in small clusters mainly influences the nonadiabatic dynamics by modulating the NaI internuclear separation at which the ionic and covalent curve-crossing region occurs. Furthermore, the simulations show that electron transfer does not occur in the inverted regime, as predicted by the free energy curves, because of the extreme nonequilibrium nature of the NaI(H2O)n photodissociation process, and the rate of electron transfer increases with cluster size and increasing solvation. Overall, this work demonstrates the importance of including polarization in realistic excited-state simulations of NaI(H2O)n relaxation.