Effects of Electronic Structure on the Hydration of Pbno 3 + and Srno 3 + Ion Pairs †

Hydration of PbNO 3 + and SrNO 3 + with up to 30 water molecules was investigated with infrared photodissociation (IRPD) spectroscopy and with theory. These ions are the same size, yet the IRPD spectra of these ion pairs for n = 2–8 are significantly different. Bands in the bonded O–H region (B3000–3550 cm À1) indicate that the onset of a second hydration shell begins at n = 5 for PbNO 3 + and n = 6 for SrNO 3 +. Spectra for [PbNO 3 ] + (H 2 O) 2–5 and [SrNO 3 ] + (H 2 O) 3–6 indicate that the structures of clusters with Pb(II) are hemidirected with a void in the coordinate sphere. A natural bond orbital analysis of [PbNO 3 ] + (H 2 O) 5 indicates that the anisotropic solvation of the ion is due to a region of asymmetric electron density on Pb(II) that can be explained by charge transfer from the nitrate and water ligands into unoccupied p-orbitals on Pb(II). There are differences in the IRPD spectra of PbNO 3 + and SrNO 3 + with up to 25 water molecules attached. IR intensity in the bonded O–H region is blue-shifted by B50 cm À1 in nanodrops containing SrNO 3 + compared to those containing PbNO 3 + , indicative of a greater perturbation of the water H-bond network by strontium. The free O–H stretches of surface water molecules in nanodrops containing 10, 15, 20, and 25 water molecules are red-shifted by B3–8 cm À1 for PbNO 3 + compared to those for SrNO 3 + , consistent with more charge transfer between water molecules and Pb(II). These results demonstrate that the different electronic structure of these ions significantly influences how they are solvated.


Introduction
2][3] Once ingested, lead deposits into tissues in the central nervous system 4,5 and various organs throughout the body, 6 disrupting their function.Divalent lead in particular inhibits the function of regulatory metalloproteins by acting as a heavy-metal substitute for metal ion cofactors.Pb(II) can bind to the Ca(II)-regulated calmodulin messenger protein 7,8 thereby falsely activating cell signaling processes. 9,10ther calcium-dependent enzymes are affected by Pb(II), including protein kinase C, whose activity can be stimulated by picomolar concentrations of lead. 11Similarly, Pb(II) has been found to inhibit the function of Zn(II) and Fe(II) metalloenzymes involved in heme biosynthesis, 2,12 resulting in anemia.Lead binding depends both on its specific interaction with proteins but also its interactions with water.Knowledge of the coordination chemistry of divalent lead is thus important for understanding its biological activity.The stability of lead in an oxidation state two less than the group valency is a characteristic shared amongst the heavier p-block elements, the so-called ''inert pair effect''. 13This property has been explained by a relativistic contraction of the 6s orbital, which stabilizes the electrons against oxidation and confers chemical inertness. 14,15Although these inert electrons do not participate in chemical bonding, crystal structures of Pb(II)containing complexes indicate that they can be categorized into two broad families depending on the apparent stereochemical activity of the inert pair. 16,17In holodirected complexes, electron density is evenly distributed around the metal ion, as are ligands in the coordination sphere.Just the opposite is the case in hemidirected Pb(II)-complexes, where asymmetric electron density on the ion results in voids in the coordination sphere.Shimoni-Livny and co-workers conducted a survey of the binding preferences of lead by analyzing 329 crystal structures of lead-containing compounds. 16They reported that Pb(II) can adopt coordination numbers (CN) ranging from 2-10, with the exclusive This journal is © the Owner Societies 2015 formation of hemidirected structures for CN = 2-5, the exclusive formation of holodirected structures for CN = 9-10, and the presence of both structures for CN = 6-8.It was initially postulated that the asymmetric electron density on Pb(II) arises from a mixing of the occupied atomic 6s orbital with unoccupied 6p orbitals, creating directional s-p hybrid orbitals. 18Ab initio calculations on solid-state Pb(II) complexes indicate that the 6s electrons remain largely unhybridized, and that asymmetric electron density on the ion is due to charge transfer between the ligands and vacant 6p orbitals on Pb(II). 19These findings better explain the presence of both holo-and hemidirected structures for the intermediate CNs of Pb(II), where there is a dependence on the identity of the ligand.
Although there is much data for solid-state Pb(II) complexes, there is little experimental data on the hydration of lead in aqueous solution. 17An early proton magnetic resonance measurement indicated that the primary hydration number for Pb(II) is 5.7, 20 although this value has been called into question. 21Persson et al. measured an extended X-ray absorption fine structure spectrum of a 0.1 M solution of Pb(ClO 4 ) 2 , and from the broad distribution of Pb-O bond distances concluded that Pb(II) is hemidirected in aqueous solution. 17Ab initio calculations on Pb 2+ (H 2 O) n clusters indicate that hemidirected structures are favored for n o 6, with holodirected structures becoming energetically competitive at n = 6. 22,23In contrast, molecular dynamics simulations of aqueous Pb(II) indicate predominantly holodirected structures. 21,24Even less is known about solvated lead complexes.
7][48] Stace and co-workers showed that many complexes of Pb(II) with organic ligands are stable in the gas phase, 49,50 but that gaseous clusters of Pb 2+ (H 2 O) n are unstable to hydrolysis for n o 11 and consequently are not present in electrospray ionization (ESI) mass spectra.51 However, they report that hydrates of the contact ion pair [PbNO 3 ] + (H 2 O) n are stable at small cluster sizes, and Raman measurements indicate that ion pairing between metal ions and nitrate can occur in solution.52 IRPD studies of ion pair hydrates have so far been limited to two reports on the first hydration shell of metal hydroxides, 53,54 and a report on 55 How the hydration of an ion pair evolves past the completion of an inner shell around the metal ion has hitherto not been studied.
Here, the hydration of [ .We attribute this structural preference to charge transfer between the nitrate and water ligands and unoccupied 6p orbitals on Pb(II).IRPD spectra of larger hydrates of these ion pairs show that differences in their hydration persist up to n = 25, well past the first hydration shell.

IRPD Spectroscopy
IRPD spectra of hydrated PbNO 3 + and SrNO 3 + were acquired using a 7.0 T Fourier-transform ion cyclotron resonance mass spectrometer coupled to a tunable OPO/OPA tabletop laser system.A detailed description of the apparatus and experiment can be found elsewhere. 56Briefly, hydrated ions are generated from nanoelectrospray ionization of 3-5 mM aqueous solutions of Pb(NO 3 ) 2 (J.T. Baker, Phillipsburg, NJ) and Sr(NO 3 ) 2 (Mallinckrodt, Paris, KY) using borosilicate glass capillaries that are pulled to an inner tip diameter of B1 mm.A platinum wire is inserted into the capillary tube so that it is in contact with the solution, and a B +600 V potential relative to the heated capillary entrance of the mass spectrometer is applied to the wire to initiate electrospray.Ions are directed from atmosphere through five stages of differential pumping into the ion cell using electrostatic lenses.The cell is surrounded by a temperature-regulated copper jacket 57 that is equilibrated to 133 K for 8 h prior to experiments.A pulse of dry nitrogen gas is introduced into the vacuum chamber at a pressure of B10 À6 Torr for B5 s to aid with thermalizing and trapping the ions.This is followed by a B7 s delay to reduce the pressure in the cell to o10 À8 Torr.
Precursor ions are selected using stored waveform inverse Fourier transforms.Mass-selected clusters are photodissociated at specific frequencies using an OPO/OPA laser system (LaserVision, Bellevue, WA) pumped by the 1064 nm fundamental of a Nd:YAG laser (Continuum Surelight I-10, Santa Clara, CA) at a 10 Hz repetition rate with pulse energies between B1-4 mJ between 3000-3800 cm À1 .Irradiation times between 1 to 10 s are used to produce substantial, but not complete, fragmentation of the isolated precursor ions which dissociate by the sequential loss of water molecules.First-order photodissociation rate constants are determined from the relative abundances of the precursor and fragment ions, and are corrected for frequency-dependent variations in laser power as well as blackbody infrared radiative dissociation (BIRD), which occurs as a result of the absorption of blackbody photons from the 133 K ion cell and cell jacket.

Computational chemistry
Low-energy structures of PbNO

Results and discussion
IRPD spectra of small PbNO 3 + and SrNO 3 + hydrates The structures of [PbNO 3 ] + (H 2 O) n and [SrNO 3 ] + (H 2 O) n with n between 2 and 8 were investigated by measuring IRPD spectra in the hydrogen stretch region (B3000-3800 cm À1 ) at 133 K (Fig. 1).Vibrational resonances in this region arise from both hydrogen-bonded (H-bonded) and ''free'' O-H stretches of water molecules.5][36][37][38] The overlaid spectra of PbNO 3 + and SrNO 3 + hydrates shown in Fig. 1 can be delineated into two principal spectral regions.Resonances in the higher-frequency free O-H region (B3550-3800 cm À1 ) arise from O-H oscillators that are not involved in H-bonding and are relatively sharp bands.The broad bands of dissociation observed in the lower-frequency region (B3000-3600 cm À1 ) are a result of bonded O-H stretches, including both water-water and water-nitrate hydrogen bonds.
The IRPD spectra of the size selected PbNO 3 + and SrNO 3 + clusters are significantly different, especially in the bonded O-H region, indicating that the solvation of these two ion pairs by water is also significantly different.1 In the IRPD spectrum of ) n spectra also provide structural information. The fee O-H region can be subdivided into two regions arising from the antisymmetric and symmetric stretches of water molecules, which appear between 3650-3750 cm À1 and 3550-3650 cm À1 , respectively (Fig. 1).The onset of H-bonding in the bonded O-H region is accompanied by changes in the free O-H band structure.There are two distinct symmetric stretching bands in the spectrum of [PbNO 3 ] + (H 2 O) 5 , whereas the symmetric stretching bands near 3600 cm À1 are degenerate for the smaller PbNO    To obtain more detailed information about structures, the IRPD spectra are compared to calculated spectra of low-energy isomers.
The and 3756 cm À1 , respectively. 63 The IRPD spectrum of [SrNO 3 ] + (H 2 O) 4 along with calculated infrared spectra of four low-energy isomers are shown in Fig. 3.
The most intense band in the experimental spectrum is the antisymmetric stretch at 3708 cm À1 that is closely reproduced in the calculated spectra of isomers 4a-4c.The positions of water molecules around the strontium ion in these three isomers differ substantially, yet the antisymmetric stretching bands are similar, indicating that the antisymmetric stretching frequencies of these inner shell water molecules are not sensitive to their position around the ion pair.The two well-defined bands in this region near 3585 cm À1 and 3635 cm À1 do not unequivocally correspond to the symmetric stretching bands in the This journal is © the Owner Societies 2015 calculated spectrum of any single isomer.The spectrum of isomer 4c (+2.5 kJ mol À1 ) is arguably the best fit with experiment, reproducing the measured number of resonances, although the lowest-energy symmetric stretch is calculated to be B60 cm À1 lower than that in the IRPD spectrum.These calculations indicate that the frequencies of the symmetric stretching bands are sensitive to the positions of water molecules within the cluster, and discrepancies between experiment and theory could arise from uncertainties in harmonic oscillator frequencies and intensities 64 as well as the presence of energetically competitive isomers.For ions where there are many interconverting isomers or that have vibrations that are significantly anharmonic, calculated harmonic spectra of low-energy isomers can differ from the experimental infrared spectra. 65,66ncertainties as a result of the harmonic frequency approximation can also affect calculated Gibbs free energies, which may affect the relative energetic ordering of isomers presented here.Structure 4d (+25.1 kJ mol À1 ) is the lowest-energy isomer identified that includes a hydrogen bond and is significantly higher in energy than all the isomers without H-bonding (4a-4c).The poor match between its calculated spectrum and the IRPD spectrum, most notably the calculated H-bonded stretch centered on 3182 cm À1 that is absent in the experiment, indicate that it does not constitute a significant fraction of the experimental population.
The experimental and calculated infrared spectra for [PbNO 3 ] + (H 2 O) 4 are shown in Fig. 4. The lowest-energy isomer (4a) reproduces the water antisymmetric and symmetric stretching bands in the IRPD spectrum near 3705 cm À1 and 3612 cm À1 , respectively, but there are additional features in the experimental spectrum of [PbNO 3 ] + (H 2 O) 4 that cannot be explained by isomer 4a alone.An intense band at 3684 cm À1 is present in the antisymmetric stretching region, as well as two broad regions of low intensity in the bonded O-H region near 3460 cm À1 and 3100 cm À1 .Collectively, these resonances may arise from a minor population of structures in which there is hydrogen bonding, such as isomer 4b (+8.1 kJ mol À1 ).This structure has a resonance at 3686 cm À1 corresponding to the antisymmetric free O-H stretch of a water molecule donating one H-bond.The bonded O-H stretch of this water molecule is calculated to occur at 3052 cm À1 , which could account for the small recorded resonance near 3100 cm À1 .The experimentally observed dissociation between 3300-3500 cm À1 is not reproduced in any of the calculated spectra, and we postulate that this band may arise from an O-H oscillator forming a weak H-bond with nitrate.Isomers 4c and 4d have calculated spectra in poor agreement with the IRPD spectrum in the free O-H region, and it is unlikely that they contribute significantly to the experimental population.

Structures of [MNO 3 ] + (H 2 O) 5
The IRPD and calculated spectra of [SrNO 3 ] + (H 2 O) 5 are shown in Fig. 5.There is good agreement between the IRPD spectrum and the calculated spectrum for lowest-energy isomer 5a.In this structure, there is no H-bonding and water molecules spherically solvate the strontium ion.There is no spectroscopic evidence for a substantial population of isomers with H-bonded structures, such as isomers 5b-5e, which are 414 kJ mol À1 higher in energy.
The IRPD spectrum of [PbNO 3 ] + (H 2 O) 5 (Fig. 6) differs significantly from that of [SrNO 3 ] + (H 2 O) 5 (Fig. 5).In striking contrast to [SrNO 3 ] + (H 2 O) 5 , there is substantial dissociation in the   S1, ESI †), and at n = 5, the water molecule binding energies of these two ions are within 20 kJ mol À1 .The maximum infrared rate constant for SrNO 3 + is within a factor of B2 of that for PbNO 3 + at n = 5.This indicates that the differences in the relative binding energies of water to these ions has a relatively small effect on the ability to observe  photodissociation, and that the differences in the spectra of these two ions in the H-bonding region are significant.The best matches between experiment and theory for [PbNO 3 ] + (H 2 O) 5 are isomer 5b (+1.7 kJ mol À1 ) and lowestenergy isomer 5c, indicating the presence of both isomers in the experiment (Fig. 6).The broad antisymmetric stretching band between 3660-3740 cm À1 in the IRPD spectrum is consistent with a combination of the antisymmetric stretches calculated for 5b and 5c.There are two free O-H symmetric stretching bands centered on 3620 cm À1 and 3585 cm À1 in the IRPD spectrum.These resonances are closely matched by the calculated symmetric stretches of isomer 5b at 3611 cm À1 and 3570 cm À1 .The greater relative intensity of the higher-energy symmetric stretch in the experimental spectrum is consistent with the presence of isomer 5c, which is calculated to have a symmetric stretching band centered on 3605 cm À1 .
Resonances in the bonded O-H region of the [PbNO 3 ] + (H 2 O) 5 spectrum can also be attributed to a combination of isomers 5b and 5c.The experimental band that has a maximum near 3475 cm À1 corresponds to the O-H stretch of the second shell water molecule in isomer 5c that weakly H-bonds with nitrate and is predicted to occur at 3504 cm À1 .The broad region of dissociation in the lowest-energy region of the experimental spectrum can be attributed to isomer 5b, which has a waterwater H-bond calculated to have a resonance near 3134 cm À1 .In general, the agreement between experimental and calculated frequencies and intensities is poorer in the bonded O-H region.8][69] and can be attributed to uncertainties associated with using the harmonic approximation, as well as effects of basis set and level of theory used. 64IRPD action spectra can also differ from computed spectra owing to fragmentation efficiencies that depend on the photon energy. 70,71or larger hydrates of SrNO 3 + and PbNO 3 + , a comparison between experiment and theory is more challenging owing to the increasing number of possible isomers and increasing computational difficulty with increasing cluster size.Our results for [SrNO 3 ] + (H 2 O) 6 , however, indicate that the outer shell water molecule accepts a H-bond from an inner shell water molecule and donates a weak H-bond to nitrate, thus forming a bridge between the hydration shell around Sr(II) and nitrate (Fig. S4, ESI †).The effects of basis set and type of theory on the energetic ordering of isomers was evaluated from single point calculations on optimized structures of the hydrated ions with n = 4 and 5 using two types of theory (MP2 and DFT/ B3LYP) in conjunction with two different basis sets (aug-cc-pVDZ and 6-31++G**).There is no significant energetic reordering of the isomers for any of the ions investigated, and in all cases, the same lowest-energy structure was identified (Tables S2  and S3, ESI †).

NBO analysis of [PbNO 3 ] + (H 2 O) 5
A NBO analysis 72 of the lowest-energy isomer of [PbNO 3 ] + (H 2 O) 5 was performed to elucidate the relationship between the apparent hemidirected nature of Pb(II) in this complex and its electronic structure.The filled NBOs f i that represent the natural Lewis structure of a molecule commonly account for 499.9% of the total electron density. 73The remaining electron density is described by non-Lewis NBOs f j of lower occupancy that describe portions of the molecular valence space not associated with the localized Lewis structure.Occupancy of these non-Lewis NBOs thus accounts for irreducible delocalization effects in the description of the total electron density and leads to a lowering of the electronic energy.The resonance delocalization of electron density from filled (donor) f i NBOs to unfilled (acceptor) f j NBOs constitutes charge transfer between the orbitals, and can be calculated from second order perturbation theory. 74Analyses of these stabilizing donor-acceptor interactions have led to chemically intuitive descriptions of hyperconjugative effects, 75 H-bonded systems, 76,77 and ion-molecule complexes. 73esults from an NBO analysis of [PbNO 3 ] + (H 2 O) 5 indicates that the electron density is strongly resonance delocalized, with the occupancy of Lewis NBOs accounting for only 99.16% of the total electron density.Shown in Fig. 7a is a set of three unfilled, non-Lewis NBOs that are primarily responsible for this delocalization, accepting charge from filled Lewis NBOs on the water and nitrate ligands (Fig. 7b).These acceptor orbitals are localized 6p atomic orbitals on Pb(II).The orientation of the 6p 1 orbital in particular is consistent with a region of repulsive electron density that could partially explain the hemidirected character of Pb(II).
Calculated stabilization energies due to charge transfer between donor NBOs on the water and nitrate ligands and the acceptor 6p NBOs on Pb are shown in Table 1.The most stabilizing interaction (91.6 kJ mol À1 ) occurs between a lone pair on an oxygen atom of nitrate (LP 1 ) and the 6p 1 orbital on Pb.These orbitals are relatively close in energy and have good spatial overlap (Fig. 7) making this the most favorable charge transfer interaction.The next most stabilizing interaction occurs between a lone pair on the second coordinating oxygen atom of the nitrate ligand (LP 2 ) and the 6p 1 orbital on Pb(II).Significant charge transfer also occurs between the four inner shell water molecules and the 6p 2 and 6p 3 NBOs on Pb(II).The high occupancies of these 6p orbitals (Table 1), especially the 6p 1 orbital that accepts charge from nitrate, indicate that they play an essential role in accurately describing the electronic structure of [PbNO 3 ] + (H 2 O) 5 .Charge transfer from nitrate to lead slightly shortens the average nitrate-cation distance to 2.56 Å in this complex as compared to 2.63 Å for [SrNO 3 ] + (H 2 O) 5 .Conversely, average cation-oxygen bond lengths for inner shell water molecules are somewhat greater for lead (2.68 Å) than strontium (2.62 Å) for n = 5 as a result of the reduced positive charge on the lead cation (Table 2).An analysis of the net atomic Mulliken charges on Pb and Sr in the lowest-energy structures for the n = 2-10 hydrates (Table 2) shows that the net atomic charge on lead is B +0.5 less than strontium for these cluster sizes, indicating that the charge transfer efficiency does not depend strongly on the extent of hydration within this range of cluster sizes.
Interestingly, the steric influence of lead's 6s electrons (i.e., the ''inert pair'') appears to be minor.valence NBO associated with these electrons is largely unhybridized, having predominately s character -the amount of p character is calculated to be only 1.47%.These results suggest that the hemidirected nature of the [PbNO 3 ] + (H 2 O) 5 complex can mostly be attributed to a charge transfer effect from the nitrate ligand into an unoccupied 6p orbital on lead.This is consistent with ab initio calculations performed on structures of lead solvated by water 22 as well as other ligands. 78Several groups have found that the calculated hybridization of the inert pair is sensitive to the size of the effective core potential used to describe lead, 22,23 potentially explaining an earlier finding where the 6s electrons were calculated to be more significantly hybridized. 16ructural signatures in IRPD spectra of larger PbNO + as a result of the onset of H-bonding at smaller cluster size.At this cluster size, the differences in water molecule binding energies between these clusters is calculated to be B6 kJ mol À1 (Table S1, ESI †) and this difference should only have a minor effect on relative dissociation efficiencies in this spectral region.The intensity of this band becomes similar for both ions with increasing cluster size, indicating that the  There are resonances at this frequency in IRPD spectra of hydrated ions that form clathrate structures with B20 water molecules attached, including the hydrated proton, [43][44][45] alkali metal ions, 82 ammonium, 83,84 and alkylammonium ions 84 corresponding to water molecules that accept one H-bond and donate two H-bonds (''ADD'' water molecules) that are integral to forming cage-like, three-dimensional H-bond networks.Results from Fujii and co-workers indicate that an intense ADD band is a general feature of small water clusters (n B 20-50), where the dominant structural motif is calculated to be ''centered cage'' structures wherein several water molecules are entrained in a three-dimensional water cage. 85The larger size of these ion pairs precludes structures where an ion pair is completely encapsulated by a water cage with only 20 water molecules or so, but it is possible that partial cage-like structures with a high percentage of three-coordinate water molecules can form around the ion pair.Indeed, the lowest-energy structures of [SrNO 3 ] + (H 2 O) 2-10 in Fig. 2   The frequencies of free O-H oscillators depend on the charge of the ion and, to a lesser extent, the size of the ion. 32Sr(II) and Pb(II) have the same ionic radii 62 indicating that the red-shift in the frequency of surface water stretches in hydrates of PbNO 3 + is due to higher charge transfer between water molecules and lead.This is consistent with findings by Duncan and co-workers, who reported that charge transfer can red-shift the O-H stretch frequencies in IRPD spectra of water molecules bound to singly charged cations. 86As our NBO analysis of [PbNO 3 ] + (H 2 O) 5 suggests, the unfilled 6p orbitals of Pb(II) can readily accept charge from donating water and nitrate ligands.Charge transfer from water molecules that are located in outer solvation shells will decrease with increasing cluster size because of poorer overlap with the atomic orbitals on the metal ion.The number of water molecules in the second and higher hydration shells increases relative to the number of water molecules in the first hydration shell and the IRPD spectra increasingly reflect absorption of these outer shell water molecules with increasing cluster size.Thus, the IRPD spectra of these clusters are expected to converge at some cluster size where the contributions to the spectra from water molecules that are significantly perturbed by these ions becomes negligible.

Conclusions
The structures of [

3 +.
photodissociation (IRPD) spectroscopy and with theory.These ions are the same size, yet the IRPD spectra of these ion pairs for n = 2-8 are significantly different.Bands in the bonded O-H region (B3000-3550 cm À1 ) indicate that the onset of a second hydration shell begins at n = 5 for PbNO 3 + and n = 6 for SrNO 3 + .Spectra for [PbNO 3 ] + (H 2 O) 2-5 and [SrNO 3 ] + (H 2 O) 3-6 indicate that the structures of clusters with Pb(II) are hemidirected with a void in the coordinate sphere.A natural bond orbital analysis of [PbNO 3 ] + (H 2 O) 5 indicates that the anisotropic solvation of the ion is due to a region of asymmetric electron density on Pb(II) that can be explained by charge transfer from the nitrate and water ligands into unoccupied p-orbitals on Pb(II).There are differences in the IRPD spectra of PbNO 3 + and SrNO 3 + with up to 25 water molecules attached.IR intensity in the bonded O-H region is blue-shifted by B50 cm À1 in nanodrops containing SrNO 3 + compared to those containing PbNO 3 3 + hydrates.This indicates that water molecules in [PbNO 3 ] + (H 2 O) 5 reside in different environments, as would be expected of inner shell versus outer shell water molecules.Changes in the free O-H band structure of [SrNO 3 ] + (H 2 O) n spectra also accompany the onset of H-bonding in the n = 6 cluster.Most apparent is the broadening of the antisymmetric stretching band centered on 3715 cm À1 in IRPD spectrum of [SrNO 3 ] + (H 2 O) 5 upon the addition of a water molecule, indicating the greater diversity of water molecule binding sites in [SrNO 3 ] + (H 2 O) 6 .Calculated lowest-energy structures for [MNO 3 ] + (H 2 O) n Calculated lowest-energy isomers of [PbNO 3 ] + (H 2 O) n and [SrNO 3 ] + (H 2 O) n for n = 2-10 at the B3LYP/aug-cc-pVDZ level of
IRPD spectra of [PbNO 3 ] + (H 2 O) 2-3 and [SrNO 3 ] + (H 2 O) 3 shown in Fig. 1 have simple band structures, and a detailed comparison between experiment and theory can be found in the ESI.† An IRPD spectrum of [SrNO 3 ] + (H 2 O) 2 was not acquired at 133 K because photodissociation for this stable complex was not observed for up to 60 s of laser irradiation.Calculations indicate that the water molecule binding energy for [SrNO 3 ] + (H 2 O) 2 (B103 kJ mol À1) is significantly higher than that for [PbNO 3 ] + (H 2 O) 2 (B77 kJ mol À1 ), consistent with the absence of observable photodissociation for the former ion under these conditions (TableS1, ESI †).The two bands appearing at 3595 cm À1 and 3685 cm À1 in the IRPD spectrum of [PbNO 3 ] + (H 2 O) 2 correspond to the symmetric and antisymmetric stretches of water molecules, which are substantially redshifted from their neutral values in water vapor of 3657 cm À1

Fig. 4
Fig. 4 IRPD spectrum of [PbNO 3 ] + (H 2 O) 4 and calculated spectra of low-energy isomers at the B3LYP/aug-cc-pVDZ level of theory, along with their relative Gibbs free energies in kJ mol À1 at 133 K.
bonded O-H region of [PbNO 3 ] + (H 2 O) 5 signifying the onset of H-bonding in this cluster.The water molecule binding energies of these ions differ, and this difference affects the relative dissociation efficiencies of these ions which is a complicating factor in comparing spectral intensities in the bonded O-H region.The water molecule binding energies decrease with increasing cluster size (Table

Fig. 5
Fig. 5 IRPD spectrum of [SrNO 3 ] + (H 2 O) 5 and calculated spectra of low-energy isomers at the B3LYP/aug-cc-pVDZ level of theory, along with their relative Gibbs free energies in kJ mol À1 at 133 K.

Fig. 6
Fig. 6 IRPD spectrum of [PbNO 3 ] + (H 2 O) 5 and calculated spectra of low-energy isomers at the B3LYP/aug-cc-pVDZ level of theory, along with their relative Gibbs free energies in kJ mol À1 at 133 K.
Fig. 7c shows that the This journal is © the Owner Societies 2015 Phys.Chem.Chem.Phys., 2015, 17, 15963--15975 | 15971 The extent to which the structural differences in nascent hydrates of PbNO 3 + and SrNO 3 + persist at larger cluster sizes was investigated by measuring IRPD spectra of [PbNO 3 ] + (H 2 O) n and [SrNO 3 ] + (H 2 O) n , where n = 10, 15, 20, 25, 30 (Fig. 8).The IRPD spectra of both ion pairs with up to 25 water molecules attached are distinguishable in the free O-H and bonded O-H regions but are indistinguishable at n = 30 (Fig. 8).Low-energy structures of these ion pairs with 25 water molecules attached generated using molecular mechanics indicate that there are several water molecules in the third solvation shell (Fig. S5, ESI †).This indicates that differences in the electronic structures of PbNO 3 + and SrNO 3 + influence their hydration past the first and even the second solvation shell.The overall intensity of the bonded O-H stretch region is greater for [PbNO 3 ] + (H 2 O) 10 than [SrNO 3 ] + (H 2 O) 10 , which reflects the more extensive H-bond network in hydrates of PbNO 3

Ion n = 2 n 10 [
Fig. 7 (a) Unfilled NBOs centered on Pb(II) that accept charge.(b) Filled NBOs on the nitrate and water ligands that donate charge.(c) NBO of valence 6s electrons on Pb(II) that has almost pure s character.All orbitals are plotted using an isosurface value of 0.05 e Å À3 .
suggests such a hydration motif, where water molecules first nucleate around Sr(II) and subsequently align along the approximately dipolar electric field of SrNO 3 + , forming H-bonds with each other and the nitrate anion.This leads to the beginnings of a cage-like hydration network in the calculated lowest-energy structure of [SrNO 3 ] + (H 2 O) 10 .Two distinct bands appear in the free O-H region of these larger PbNO 3 + and SrNO 3 + clusters (Fig. 8): an intense band near 3700 cm À1 associated with water molecules accepting two H-bonds and donating one H-bond (AAD water molecule), and a small band near 3720 cm À1 arising from water molecules that accept and donate only one H-bond (AD water molecule).The free O-H stretches of water molecules in [PbNO 3 ] + (H 2 O) n are red-shifted by 5.6-7.6 cm À1 compared to those in [SrNO 3 ] + -(H 2 O) n for n = 10, 15 and 20, but they are closer in frequency at n = 25 (3.3 cm À1 ) and are indistinguishable at n = 30.The centroid frequencies of the AAD free O-H bands were determined by fitting each peak with a Gaussian function, and the frequencies for all cluster sizes are shown in Fig. S6 (ESI †).
PbNO 3 ] + (H 2 O) n and [SrNO 3 ] + (H 2 O) n were investigated with IRPD spectroscopy and with theory.Although these ions are the same size, they differ significantly in their hydration.H-bonded stretches in the IRPD spectra indicate the onset of second shell formation at n = 5 for PbNO 3 + (and a minor population of second shell structures at n = 4) and n = 6 for SrNO 3 + .Calculated structures of PbNO 3 + and SrNO 3 + with 10 and fewer water molecules indicate that hemidirected structures are favored for Pb(II) where there is a noticeable void in its coordination sphere, and these results are consistent with the IRPD spectra.The asymmetric solvation for hydrates of PbNO 3 + is due to the asymmetric electron density on Pb(II) that repels water molecules away from the portion of its coordination sphere opposite to the nitrate ligand, constraining the available binding site for water molecules.Results from an NBO analysis of [PbNO 3 ] + (H 2 O) 5 indicate that the asymmetric charge density on Pb(II) is a result of charge transfer from the nitrate and water ligands into vacant 6p orbitals on the ion.In contrast, Sr(II) in SrNO 3 + is isotropically solvated by water.The IRPD spectra of PbNO 3 + and SrNO 3 + with up to 25 water molecules attached differ.The maximum intensity of the H-bonded stretches of SrNO 3 + are blue-shifted compared to those of PbNO 3 + , indicating a greater perturbation of the water H-bond network by strontium than lead.This is consistent with the higher net atomic Mulliken charges of [SrNO 3 ] + (H 2 O) n compared to [PbNO 3 ] + (H 2 O) n .The intense band at 3550 cm À1 in the spectra of SrNO 3 + with 15, 20, and 25 waters attached suggests that partial clathrate-like structures may form around this ion pair.The free O-H stretches of AAD water molecules in hydrates of PbNO 3 + are red-shifted by up to B8 cm À1 compared to those in hydrates of SrNO 3 + for n = 10, 15, 20 but are the same by n = 30.The red-shifting of surface water molecules in [PbNO 3 ] + (H 2 O) n is consistent with the transfer of electron density from their O-H bonds to the lead ion.Collectively, the spectral differences between the larger hydrates of PbNO 3 + and SrNO 3 + demonstrate that the electronic structure of an ion can play an important role in how it is solvated past the first and even second solvent shells.
PbNO 3 ] + (H 2 O) n and [SrNO 3 ] + (H 2 O) n is probed with IRPD spectroscopy, with the latter species serving as a reference for an ion pair incorporating a metal ion with a closed shell of electrons.IRPD spectra of [MNO 3 ] + (H 2 O) n with 2 r n r 8 show that the onset of hydrogen bonding (H-bonding) begins at n = 5 for PbNO 3 + and n = 6 for SrNO 3 + indicating differences in their inner hydration shells.Comparisons between IRPD and calculated spectra of [PbNO 3 ] + (H 2 O) 2-5 and [SrNO 3 ] + (H 2 O) 3-6 show that hemidirected structures are energetically favorable for PbNO 3 + Pb(II) in PbNO 3 + does not form in this range of cluster sizes, leading to structures where the ion pair is not isotropically solvated by water.The directional distribution of electron density on Pb establishes a ''hydrophobic'' surface that repels water molecules away from part of the ion so they preferentially form water-water H-bonds, eschewing short-range ion-dipole interactions with the cation.It is interesting to note that for both [SrNO 3 ] + (H 2 O) n and [PbNO 3 ] + (H 2 O) n , these calculations predict that the nitrate ligand transitions from being bidentate to monodentate between the n = 8 and 10 cluster sizes, presumably to maximize the number of H-bonds formed between water molecules and nitrate.
This journal is © the Owner Societies 2015 Phys.Chem.Chem.Phys., 2015, 17, 15963--15975 | 15967 around The calculated spectrum of the lowest-energy isomer of [PbNO 3 ] + (H 2 O) 2 matches closely with the IRPD spectrum (Fig.S1, ESI †).Similarly, the IRPD spectra of [PbNO 3 ] + (H 2 O) 3 and [SrNO 3 ] + (H 2 O) 3 have resonances associated with the symmetric and antisymmetric stretching motions of water molecules, and match most closely with the calculated spectra of the lowest-energy isomers (Fig.S2 and S3, ESI †).The broadened symmetric stretching band in the IRPD spectrum of [SrNO 3 ] + (H 2 O) 3 may arise from the existence of low-lying isomers where water molecules are distributed more evenly around the strontium ion (Fig.S3, ESI †).The maximum infrared photodissociation rate constant for [SrNO 3 ] + (H 2 O) 3 is about four times lower than that for [PbNO 3 ] + (H 2 O) 3 indicating that water molecules in this complex are more tightly bound and that the absorption of multiple photons may be required to induce observable laser-induced photodissociation on the timescale of the experiment.

Table 1
Resonance delocalization energies between donor NBO f i and acceptor NBO f j calculated by second order perturbation theory.The occupancies of the acceptor orbitals as well as the acceptor (e j ) to donor (e i ) orbital energy difference are also shown