Binding Energies of Hexahydrated Alkaline Earth Metal Ions, M2+(H2O)6, M = Mg, Ca, Sr, Ba: Evidence of Isomeric Structures for Magnesium

One strategy that has commonly been used to gain a more detailed understanding of ion chemistry in bulk solution is to investigate the structure and energetics of solvated ions in the gas phase. While solvation of singly charged metal ions has been studied extensively during the last two decades,1 significantly less information about doubly and triply charged metal ions is known. Solvated divalent metal ions can now be readily produced using electrospray ionization (ESI).2–5 This has led to a significant interest in obtaining thermochemical information about these ions with both experiment3–5 and theory.6–8 Kebarle and co-workers measured the Gibbs free energies of hydration for water molecules located in the second solvation shell of a variety of divalent metal ions using equilibrium experiments.3 More recently, Posey and co-workers combined ESI with laser photofragmentation mass spectrometry to study divalent transition metal–ligand complexes with methanol in the second solvation shell.4 Binding energies of inner solvent shell water molecules around Ni2+ and Ca2+ ions have been determined from blackbody infrared radiative dissociation (BIRD) experiments.5 Here, BIRD kinetics of the hexahydrated alkaline earth metal ions, Mg2+, Ca2+, Sr2+, and Ba2+, are presented. At low temperatures, binding energies obtained from these kinetic data are directly correlated with the radii of the metal ions (Mg2+ > Ca2+ > Sr2+ > Ba2+). In contrast, the binding energies at higher temperatures follow the trend Ca2+ > Mg2+ > Sr2+ > Ba2+. This is the first direct evidence for two distinct gas-phase structures for a hydrated divalent metal ion. Experiments were performed using an external electrospray ionization source Fourier transform mass spectrometer that has been described previously.9 Hydrated alkaline earth M2+(H2O)n ions were generated from ~10−4 M aqueous solutions of the metal chloride salts using nanoelectrospray. Ions are loaded into the cell for 5 s during which time N2 gas (10−6 Torr) is introduced. Dissociation kinetics of the mass-selected ion are investigated at pressures <10−8 Torr. At these pressures, ion activation occurs by absorption of blackbody photons emitted by the heated vacuum chamber walls.9–11 Unimolecular dissociation rate constants in the zero-pressure limit are measured as a function of temperature (22–140 °C) from which Arrhenius activation parameters are obtained. Under these conditions, the ion population does not have a Boltzmann distribution of internal energies.10 Measured Arrhenius preexponentials and activation energies are smaller than those measured in the rapid energy exchange (REX) limit. To determine the threshold dissociation energies, master equation modeling of the kinetic data is necessary.12 The application of the master equation model to obtain threshold dissociation energies has been described in detail previously.11 The reverse activation barriers for these dissociation reactions should be negligible. Thus, the threshold dissociation energy should be equal to the binding energy. For all metal ion hydrates, dissociation occurs exclusively by loss of one water molecule. The kinetics are first order at all temperatures. Figure 1 shows the Arrhenius plots for dissociation of M2+(H2O)6 and M2+(H2O)7. For the heptahydrated clusters, the dissociation rate constants follow the trend Ba2+ > Sr2+ > Ca2+ > Mg2+. This trend is consistent with the expected reactivity based on the ionic radii of the metals (water molecules attached to larger metal ions are expected to be less strongly bound). For the hexahydrated ions, this same ordering of rate constants is observed at low temperatures. In striking contrast, the order is Ba2+ > Sr2+ > Mg2+ > Ca2+ at higher temperatures; the slope of the Arrhenius plot for Mg2+(H2O)6 changes around 80°C! BIRD rate constants for the pentahydrated species have the same trend as that observed for the hexahydrated species at high temperature (Ba2+ > Sr2+ > Mg2+ > Ca2+).13 Thus, these dissociation kinetics do not follow the trend expected based on ionic radii. The change in the slope of the Arrhenius plot for Mg2+-(H2O)6 shows the presence of two distinct gas-phase structuresfor this ion. The first-order behavior (>80% precursor ion depletion) at all temperatures indicates the presence of one structure or several structures that interconvert. Figure 1 Arrhenius plots obtained from blackbody infrared radiative dissociation data of the hexa- and heptahydrated alkaline earth metal ions; Mg2+ (○), Ca2+ (•), Sr2+ (*), and Ba2+ (▵). The measured zero-pressure limit Arrhenius parameters and the binding energies (Eo) obtained from master equation modeling of the kinetic data are given in Table 1. The low-temperature Eo’s consistently increase with decreasing cation size. At high temperatures, the binding energy of the sixth water molecule to Mg2+ (21.1 kcal/mol) is less than that at low temperatures (23.5 kcal/mol). It is also slightly less than that for the larger calcium cation (21.6 kcal/mol). For Mg2+(H2O)6, anomalously large radiative rate constants were used in the modeling in order to fit the kinetic data at low temperature. The reason for this is not clear, but this suggests that the error in Eo at low temperature may be larger than reported. The unexpected ordering of binding energies at high temperature is consistent with ΔH values recently reported by Kebarle and co-workers from equilibrium experiments at 170–225°C.3a Recent density functional theory (DFT) calculations by Pavlov et al.8 have been used to evaluate structures for hydrated divalent magnesium and calcium ions, M2+(H2O)m-(H2O)n, where m and n correspond to the number of water molecules in the first and second shell, respectively. The successive water binding energies (ΔEbinding) reported for these structures are also included in Table 1. Table 1 Measured Zero-Pressure Limit Arrhenius Parameters and EO Values Obtained from Master Equation Modeling for Loss of Water from M2+(H2O)6 Ionsa To further investigate this unexpected reactivity of Mg2+(H2O)6, we also performed DFT calculations on all four hexahydrated divalent metals.14 Three different types of structures were studied in which the six water molecules are either all in the first shell or distributed between the first and second solvation shells around the metals. By the use of the nomenclature of Pavlov et al., the structures investigated were M2+(H2O)6, M2+(H2O)5(H2O),15 and M2+(H2O)4(H2O)2. Calculations indicate that the order of stability is M2+(H2O)6 > M2+(H2O)5(H2O) > M2+(H2O)4(H2O)2 for all of the metals. The energy differences between each of the higher energy structures and the M2+(H2O)6 conformation are reported in Table 2. Also included in the table are the values of Pavlov et al.8 for magnesium and calcium. At this level of theory, all three hexahydrated Mg2+ structures have energies within 5 kcal/mol. Several similar energy structures for hexahydrated barium ions were also found. However, no evidence of multiple gas-phase hexahydrated barium structures was observed over the temperature range of these experiments Table 2 B3LYP Energy Differences between M2+(H2O)5(H2O) and M2+(H2O)4(H2O)2 and the Most Stable M2+(H2O)6 Structurea Although the experimental evidence for two distinct structures for Mg2+(H2O)6 is clear, the exact nature of these isomers is not. One possibility is that the isomers have different number of water molecules in the first and second shell. The results of DFT calculations show that all three hexahydrated Mg2+ structures investigated have 0 K energies within 5 kcal/mol (Table 2). The entropies of Mg2+(H2O)5(H2O) and Mg 2+(H2O)4(H2O)2 under standard conditions are higher than that of Mg2+(H2O)6. The opposite is true for the Ca2+, Sr2+, and Ba2+ hexahydrates. The anomalous entropy for Mg2+(H2O)6 is presumably due to the crowding of the six water molecules around the much smaller metal ion. This would suggest that the higher temperature structure corresponds to a two-shell structure. Another possible explanation is that one of the isomers is a salt-bridge structure (Scheme 1). Such a salt-bridge structure was found to be stable at the semiempirical PM3 level of theory, but not with the density functional levels of theory applied. The proton-transfer process necessary for the formation of a salt-bridge is expected to occur more readily for hydrated Mg than for the larger metal hydrates.3b Scheme 1 In conclusion, evidence for two distinct isomeric structures of Mg2+(H2O)6 is presented. This is the first experimental evidence for two isomeric structures of hydrated metal dications. Binding energies of water to Mg2+(H2O)7 and Mg2+(H2O)6 (low temperature) follow the expected trend in cation size. These values for Mg2+(H2O)6 (high temperature) do not. Calculations suggest that the most stable structure for M2+(H2O)6 is one in which all of the water molecules are in the inner shell. However, further experiments and higher level calculations are necessary to definitively elucidate these structures.