Ab Initio Molecular Dynamics Studies of Ionic Dissolution and Precipitation of Sodium Chloride and Silver Chloride in Water Clusters, NaCl(H2O)n and AgCl(H2O)n,n = 6, 10, and 14

Master equation modeling of water dissociation in small ionic water clusters: Ag+(H2O) n , n = 4-6

We model temperature-dependent blackbody infrared radiative dissociation (BIRD) rate coefficients of Ag+(H2O) n , n = 4-6, a system with loosely bound water molecules. We employ a master equation modeling (MEM) approach with consideration of absorption and emission of blackbody radiation, comparing single and multiple-well descriptions. The unimolecular dissociation rate coefficients are obtained using the Rice-Ramsperger-Kassel-Marcus (RRKM) theory, employing two approaches to model the sum of states in the transition state, the rigid activated complex (RAC) and the phase space limit (PSL) approach. A genetic algorithm is used to find structures of low-lying isomers for the kinetic modeling. We show that the multiple-well MEM approach with PSL RRKM in the All Wells and Transition States Are Relevant (AWATAR) variant provides a reliable description of Ag+(H2O) n BIRD, in agreement with previously published experimental data. Higher-lying isomers contribute significantly to the overall dissociation rate coefficient, underlying the importance of the multiple-well ansatz in which all isomers are treated on the same footing.

Structures and bonding characteristics of KCl(H2O)n clusters with n = 1-10 based on density functional theory

Aqueous inorganic salt solutions play a prominent role in both physiological and chemical experiments, and significant attention has been directed toward understanding the mechanisms underlying salt dissolution. In our effort to elucidate the hydration process of potassium chloride, we employed a comprehensive genetic algorithm to explore the structures of KCl(H2O)n (n = 1-10). A series of stable structures were identified by high-level ab initio optimization and single-point energy calculations with a zero-point energy correction. An analysis of the probability distribution of KCl(H2O)1-10 revealed that clusters with high probability at low temperatures exhibit reduced probabilities at higher temperatures, while others become more prevalent. When n = 1-9, the contact ion pair configurations or partially dissociated structures dominate in the system, and the probability distribution plot shows that the proportion of the solvent-separated ion pair (SSIP) structures of KCl(H2O)n is very small, while the SSIP configuration in KCl(H2O)10 becomes a stable structure with increasing temperature. The results from natural bond orbital analysis reveal a stronger interaction between chloride ions and water molecules. These findings provide valuable insights for a more comprehensive understanding of the intricacies of potassium chloride dissolution in water.

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.

Surface or Internal Hydration - Does It Really Matter?

The precise location of an ion or electron, whether it is internally solvated or residing on the surface of a water cluster, remains an intriguing question. Subtle differences in the hydrogen bonding network may lead to a preference for one or the other. Here we discuss spectroscopic probes of the structure of gas-phase hydrated ions in combination with quantum chemistry, as well as H/D exchange as a means of structure elucidation. With the help of nanocalorimetry, we look for thermochemical signatures of surface vs internal solvation. Examples of strongly size-dependent reactivity are reviewed which illustrate the influence of surface vs internal solvation on unimolecular rearrangements of the cluster, as well as on the rate and product distribution of ion-molecule reactions.

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.

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₂.

Generating Excess Protons in Microsolvated Acid Clusters under Ambient Conditions: An Issue of Configurational Entropy versus Internal Energy

Acid dissociation, and thus liberation of excess protons in small water droplets, impacts on diverse fields such as interstellar, atmospheric or environmental chemistry. At cryogenic temperatures below 1 K, it is now well established that as few as four water molecules suffice to dissociate the generic strong acid HCl, yet temperature-driven recombination sets in simply upon heating that cluster. Here, the fundamental question is posed of how many more water molecules are required to stabilize a hydrated excess proton at room temperature. Ab initio path integral simulations disclose that not five, but six water molecules are needed at 300 K to allow for HCl dissociation independently from nuclear quantum effects. In order to provide the molecular underpinnings of these observations, the classical and quantum free energy profiles were decomposed along the dissociation coordinate in terms of the corresponding internal energy and entropy profiles. What decides in the end about acid dissociation, and thus ion pair formation, in a specific microsolvated water cluster at room temperature is found to be a fierce competition between classical configurational entropy and internal energy, where the former stabilizes the undissociated state whereas the latter favors dissociation. It is expected that these are generic findings with broad implications on acid-base chemistry depending on temperature in small water assemblies.

Emergence of Solvent-Separated Na+-Cl- Ion Pair in Salt Water: Photoelectron Spectroscopy and Theoretical Calculations

Solvation of salts in water is a fundamental physical chemical process, but the underlying mechanism remains unclear. We investigated the contact ion pair (CIP) to solvent-separated ion pair (SSIP) transition in NaCl(H₂O)_n clusters with anion photoelectron spectroscopy and ab initio calculations. It is found that the SSIP type of structures show up at n = 2 for NaCl^-(H₂O)_n anions. For neutral NaCl(H₂O)_n, the CIP structures are dominant at n < 9. At n = 9-12, the CIP structures and SSIP structures of NaCl(H₂O)_n are nearly degenerate in energy, coincident to the H₂O:NaCl molar ratio of NaCl saturated solution and implying that the CIP and SSIP structures can coexist in concentrated solutions. These results are useful for understanding the solvation of salts at the molecular level.

How Many Water Molecules Does it Take to Dissociate HCl?

The potential energy surfaces of the HCl(H2O)n (n is the number of water molecules) clusters are systematically explored using density functional theory and high-level ab initio computations. On the basis of electronic energies, the number of water molecules needed for HCl dissociation is four as reported by some experimental groups. However, this number is five owing to the inclusion of entropic factors. Wiberg bond indices are calculated and analyzed, and the results provide a quadratic correlation and classification of clusters according to the nondissociated, partially dissociated, and fully dissociated character of the H-Cl bond. Our computations show that if temperature is not controlled during the experiment, the values obtained for the dipole moment (or for any measurable property) are susceptible to change, providing a different picture of the number of water molecules needed for HCl dissociation in a nanoscopic droplet.

Physicochemical Properties of Ion Pairs of Biological Macromolecules

Ion pairs (also known as salt bridges) of electrostatically interacting cationic and anionic moieties are important for proteins and nucleic acids to perform their function. Although numerous three-dimensional structures show ion pairs at functionally important sites of biological macromolecules and their complexes, the physicochemical properties of the ion pairs are not well understood. Crystal structures typically show a single state for each ion pair. However, recent studies have revealed the dynamic nature of the ion pairs of the biological macromolecules. Biomolecular ion pairs undergo dynamic transitions between distinct states in which the charged moieties are either in direct contact or separated by water. This dynamic behavior is reasonable in light of the fundamental concepts that were established for small ions over the last century. In this review, we introduce the physicochemical concepts relevant to the ion pairs and provide an overview of the recent advancement in biophysical research on the ion pairs of biological macromolecules.

Microsolvation of LiBO2 in water: anion photoelectron spectroscopy and ab initio calculations

The microsolvation of LiBO2 in water was investigated by conducting anion photoelectron spectroscopy and ab initio studies on the LiBO2(H2O)n(-) (n = 0-5) clusters. By comparing calculations with experiments, the structures of these clusters and their corresponding neutrals were assigned, and their structural evolutions were revealed. During the anionic structural evolution with n increasing to 5, hydroxyborate and metaborate channels were identified and the metaborate channel is more favorable. For the hydroxyborate structures, the anionic Li(+)-BO2(-) ion pair reacts with a water molecule to produce the LiBO(OH)2(-) moiety and three water molecules tend to dissolve this moiety. In the metaborate channel, two types of solvent-separated ion pair (SSIP) geometries were determined as the ring-type and linear-type. The transition from the contact ion pair (CIP) to the ring-type of SSIP starts at n = 3, while that to the linear-type of SSIP occurs at n = 4. In neutral LiBO2(H2O)n clusters, the first water molecule prefers to react with the Li(+)-BO2(-) ion pair to generate the LiBO(OH)2 moiety, analogous to the bulk crystal phase of α-LiBO2 with two O atoms substituted by two OH groups. The Li-O distance in the LiBO(OH)2 moiety increases with the increasing number of water molecules and elongates abruptly at n = 4. Our studies provide new insight into the initial dissolution of LiBO2 salt in water at the molecular level and may be correlated to the bulk state.

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.

Molecular dynamic simulations of ocular tablet dissolution

Small tablets for implantation into the subconjunctival space in the eye are being developed to inhibit scarring after glaucoma filtration surgery (GFS). There is a need to evaluate drug dissolution at the molecular level to determine how the chemical structure of the active may correlate with dissolution in the nonsink conditions of the conjunctival space. We conducted molecular dynamics simulations to study the dissolution process of tablets derived from two drugs that can inhibit fibrosis after GFS, 5-fluorouracil (5-FU) and the matrix metalloprotease inhibitor (MMPi), ilomastat. The dissolution was simulated in the presence of simple point charge (SPC) water molecules, and the liquid turnover of the aqueous humor in the subconjunctival space was simulated by removal of the dissolved drug molecules at regular intervals and replacement by new water molecules. At the end of the simulation, the total molecular solvent accessible surface area of 5-FU tablets increased by 60 times more than that of ilomastat as a result of tablet swelling and release of molecules into solution. The tablet dissolution pattern shown in our molecular dynamic simulations tends to correlate with experimental release profiles. This work indicates that a series of molecular dynamic simulations can be used to predict the influence of the molecular properties of a drug on its dissolution profile and could be useful during preformulation where sufficient amounts of the drug are not always available to perform dissolution studies.

Hydrated metal ions in the gas phase

Studying metal ion solvation, especially hydration, in the gas phase has developed into a field that is dominated by a tight interaction between experiment and theory. Since the studied species carry charge, mass spectrometry is an indispensable tool in all experiments. Whereas gas-phase coordination chemistry and reactions of bare metal ions are reasonably well understood, systems containing a larger number of solvent molecules are still difficult to understand. This review focuses on the rich chemistry of hydrated metal ions in the gas phase, covering coordination chemistry, charge separation in multiply charged systems, as well as intracluster and ion-molecule reactions. Key ideas of metal ion solvation in the gas phase are illustrated with rare-gas solvated metal ions.