Vibrational Raman Spectroscopy of the Hydration Shell of Ions

: Ionic perturbation of water has important implications in various chemical, biological and environmental processes. Previous studies revealed the structural and dynamical perturbation of water in the presence of ions, mainly with concentrated electrolyte solutions having signiﬁcant interionic interactions. These investigations highlighted the need of selective extraction of the hydration shell water from a dilute electrolyte solution that is largely free from interionic interactions. Double-difference infrared (DDIR) and Raman multivariate curve resolution (Raman-MCR), as well as MD simulation, provided valuable insight in this direction, suggesting that the perturbed water mainly resides in the immediate vicinity of the ion, called the hydration shell. Recently, we have introduced Raman difference spectroscopy with simultaneous curve ﬁtting (Raman-DS-SCF) analysis that can quantitatively extract the vibrational response of the perturbed water pertaining to the hydration shell of fully hydrated ions/solute. The DS-SCF analysis revealed novel hydrogen-bond (H-bond) structural features of hydration water, such as the existence of extremely weakly interacting water–OH ( ν max ~ 3600 cm − 1 ) in the hydration shell of high-charge-density metal ions (Mg 2+ , Dy 3+ ). In addition, Raman-DS-SCF retrieves the vibrational response of the shared water in the water–shared-ion pair (WSIP), which is different from the hydration shell water of either the interacting cation and anion. Herein, we discuss the perturbation of water H-bonding in the immediate vicinity of cation, anion, zwitterion and hydrophobes and also the inter-ionic interactions, with a focus on the recent results from our laboratory using Raman-DS-SCF spectroscopy.


Introduction
Water plays pivotal roles in atmospheric, chemical and biological processes [1][2][3][4]. In the real world, water rarely exists in its pure form; ions and amphiphiles are integral parts of natural water. Water as a liquid is highly structured due to its extensive hydrogen (H)bond network, which is easily affected by external stimuli, including temperature, pressure, charge (ions) and other solutes with varying H-bonding abilities [5][6][7][8]. The majority of molecular interactions that dictate cellular processes in living organisms occur in the ion and amphiphile containing aqueous media known as the intra-and extra-cellular fluids. The water molecules in the immediate vicinity of ions and biomolecules not only act as a solvent but also as a 'reactant'. The H-bonding of water is affected in the hydration shell of the ions and amphiphilic/hydrophobic biomolecules, which in turn influences the properties of these ions/molecules that play pivotal roles in biophysical/chemical processes, including the stability of secondary, tertiary and quaternary structures of proteins; selectivity of ion channel; cell-signaling; and the formation and flexibility of lipid membranes [9]. Being strongly associated, water exhibits the rapid delocalization of energy via its intra-and intermolecular vibrational modes; the structural perturbation of water in the ion hydration shell also affects the energy delocalization and reorganization dynamics in the hydration shell [6,[10][11][12][13][14][15]. fails to appreciate the subtle spectral features of water in the hydration shell of an ion. This limitation can be overcome to a great extent by using spectral fitting-guided minimum-area difference spectroscopy, known as difference spectroscopy, with simultaneous curve fitting (DS-SCF) analysis [41]. The application of Raman-DS-SCF revealed new features of water in the hydration shell of ions [42], which were obscured in the Raman-MCR [43], as well as double difference infrared, spectroscopy (DDIR) [44].
This review article describes the perturbation of the H-bonding and the vibrational coupling of water in the hydration shell of different ions, using DS-SCF and MCR methods on the Raman spectra of the corresponding dilute electrolyte solutions. Herein, hydration water refers to water that are in the close vicinity of ions, collectively called hydration shell water, regardless of their positions in the first, second or third hydration spheres that may exist depending on the charge density of the ion. The subsequent discussions are structured as follows. Section 2 provides a comparative discussion of the Raman-DS-SCF method with other hydration shell spectroscopy techniques. Sections 3 and 4 describe the Raman spectrum of water in the OH stretch region and discuss the effect of isotopic dilution and electrolyte on the spectral line shape. Sections 5 and 6 discuss the perturbation of H-bonding and vibrational coupling in the hydration shell of the anion and cation, respectively, with their varying natures such as charge density, valency, atomic/molecular nature, etc. The restructuring of hydration shell water for the formation of WSIP is described in Section 7. Sections 8 and 9 discuss the hydration characteristics of a hydrophobic molecular ion (CH 3 ) 4 N + ;) and zwitterion (CH 3 ) 3 N + -O − ), respectively. Finally, Section 10 provides a future perspective along with concluding remarks.

Raman-DS-SCF Analysis and Ion Hydration Shell
For the DS-SCF analysis, the Raman spectrum of a dilute electrolyte solution (~sub molar concentration) is considered as a linear combination of three spectral components: (1) vibrational response of the solute/ion (if any), (2) spectrum of the water that is perturbed by the ion/solute and (3) spectrum of the unperturbed water that is equivalent to the spectrum of the bulk water. Components 1 and 2 are associated with the solute and its hydration shell water, which is combinedly treated as the ion/solute-perturbed spectrum or the hydration shell spectrum. While subtracting the bulk-like water spectrum (reference spectrum; R(ν)) from the solution spectrum (S(ν)), the emerging perturbed-water spectrum, along with the inherent vibrational bands of the solute (if any), is fitted with a combination of Gaussian bands following Equation (1).
A n e − (ν−νn ) 2 2Γn 2 where f is the "appropriate fraction" of the reference spectrum (R(ν)) subtracted from the dilute solution spectrum (S(ν)). A n , ν n and Γ n are the amplitude, center frequency and fwhm of the n th band in components 1 and 2 (i.e., solute/ion's own vibrational bands and that of the perturbed water). A 0 is a constant that takes care of the difference in constant background (if any) between the solution and reference spectra. Figure 1 depicts the S(ν) for 0.25 M aqueous MgCl 2 , R(ν) corresponding to the 0.5 M aqueous NaCl and the Raman-DS-SCF-extracted hydration shell spectrum of Mg 2+ ion. Instead of the neat water, a 0.5 M NaCl solution was chosen as the reference (R(ν)), because the latter contains the same concentration of Cl − (counter-ion) as in the 0.25 M aqueous MgCl 2 , such that the effect of the counter ion (Cl − ) is nullified in the DS-SCF-extracted spectrum, representing the hydration shell water of Mg 2+ . To be precise, the Mg 2+ hydration shell spectrum shown in Figure 1 is the response of the Mg 2+ -associated water relative to that of Na + ; however, due to the subtle influence of Na + (on water) [45], the perturbed water spectrum predominantly represents the hydration shell response of Mg 2+ . Liquids 2022, 2, FOR PEER REVIEW 4 Figure 1. Input and output spectra corresponding to the extraction of Mg 2+ hydration shell spectrum using Raman-DS-SCF analysis. Input: Raman spectra of 0.25 M MgCl2 (red) as the solution spectrum and 0.5 M NaCl (black) as the reference spectrum. Output: Mg 2+ hydration shell spectra (blue) and spectrum of water that are unperturbed by Mg 2+ ( fR ν ), which is basically the fraction of the reference spectrum removed from the solution spectrum. The coefficients of the Raman-DS-SCF fitting analysis are shown in Table 1 for reference. As can be seen in the hydration shell spectrum of the Mg 2+ ion (blue curve, Figure 1), a weak band in the high-frequency region (~ 3600 cm −1 ) is distinctly visible in the Raman-DS-SCF extracted spectrum. Such subtle spectral features are not so obvious in the corresponding Raman-MCR-extracted Mg 2+ hydration shell spectrum [43]. Moreover, unlike in the case of Raman-MCR, the Raman-DS-SCF extracted spectrum is not area-normalized, and hence retains quantitative information of the perturbed hydration shell water to a great extent [42]. The DDIR approach also applies difference spectroscopy and curve fitting, but the "fitting" analysis is carried out after the completion of the "subtraction" process [46][47][48], while in Raman-DS-SCF, the "subtraction" and "curve fitting" are carried out simultaneously with the gradual (small-step) increments of the 'f' factor while readjusting the Gaussian fitting coefficients. For the electrolytes with the dominant effect of the counter ion, unlike the DDIR, Raman-DS-SCF and Raman-MCR use an appropriate reference spectrum to automatically subtract the counter ion influence so that the Input and output spectra corresponding to the extraction of Mg 2+ hydration shell spectrum using Raman-DS-SCF analysis. Input: Raman spectra of 0.25 M MgCl 2 (red) as the solution spectrum and 0.5 M NaCl (black) as the reference spectrum. Output: Mg 2+ hydration shell spectra (blue) and spectrum of water that are unperturbed by Mg 2+ (− fR(ν)), which is basically the fraction of the reference spectrum removed from the solution spectrum. The coefficients of the Raman-DS-SCF fitting analysis are shown in Table 1 for reference.
As can be seen in the hydration shell spectrum of the Mg 2+ ion (blue curve, Figure 1), a weak band in the high-frequency region (~3600 cm −1 ) is distinctly visible in the Raman-DS-SCF extracted spectrum. Such subtle spectral features are not so obvious in the corresponding Raman-MCR-extracted Mg 2+ hydration shell spectrum [43]. Moreover, unlike in the case of Raman-MCR, the Raman-DS-SCF extracted spectrum is not area-normalized, and hence retains quantitative information of the perturbed hydration shell water to a great extent [42]. The DDIR approach also applies difference spectroscopy and curve fitting, but the "fitting" analysis is carried out after the completion of the "subtraction" process [46][47][48], while in Raman-DS-SCF, the "subtraction" and "curve fitting" are carried out simultaneously with the gradual (small-step) increments of the 'f' factor while readjusting the Gaussian fitting coefficients. For the electrolytes with the dominant effect of the counter ion, unlike the DDIR, Raman-DS-SCF and Raman-MCR use an appropriate reference spectrum to automatically subtract the counter ion influence so that the subtracted spectrum is predominantly from the ion in question. Nevertheless, it is important to note that both Raman-MCR and Raman-DS-SCF are biased towards the perturbed water in the hydration shell, and the hydration water component having bulk-like vibrational response is not retained in the extracted spectrum.

OH Stretch Spectra of Neat Water
The Raman OH stretch (3000-3800 cm −1 ) spectrum of water is affected in the presence of dissolved ions/solutes. These spectral changes are better understood by comparing the solution spectrum with the neat-water spectrum. Therefore, it important to briefly discuss the OH stretch spectrum of neat water. The Raman spectrum of neat water ( Figure 2) shows a broad band in the OH stretch region, with a maximum at ∼3420 cm −1 and shoulders at ∼3250 and ∼3620 cm −1 . The frequency of OH stretch vibration is strongly correlated with the H-bond strength: strengthening of the H-bond makes the parent O-H bond weaker, leading to a shift in the O-H stretch frequency to a lower wavenumber, termed the 'red-shift'. Similarly, the weakening of H-bonding causes a shift in the O-H stretch frequency to a higher wavenumber, termed the 'blue-shift'. However, the intramolecular vibrational coupling between the HOH-bend overtone and O-H stretch fundamental, called Fermi resonance (FR), as well as the intra-and intermolecular vibrational coupling of the OH-stretch, change the line shape of the OH stretch spectrum of water, barring its straightforward interpretation solely on the basis of the strength of H-bonding [37,[49][50][51][52]. To correctly interpret the spectral change due to H-bonding, generally isotopically diluted water (HOD; a mixture of H 2 O and D 2 O) is used instead of isotopically pure H 2 O or D 2 O. An appropriate mixture of H 2 O and D 2 O (e.g., H 2 O:D 2 O = 1:9, v/v) is free from the intraand intermolecular vibrational coupling due to the energetic mismatch of intramolecular vibrational modes (OD and OH stretches) and the bend overtones. As can be seen in Figure 2a, the peak-normalized OH stretch spectrum of HOD has substantially lower intensity around 3200 cm −1 compared to that of H 2 O. This indicates that the 3200 cm −1 band of water is predominantly due to intra-and inter-molecular vibrational coupling and Fermi resonance. Polarized Raman measurement provides the symmetry information of the vibrational bands. The lower-frequency region (<3300 cm −1 ) is dominated in the isotropic Raman spectrum, whereas the higher frequency region (>3300 cm −1 ) is dominated in the anisotropic spectrum ( Figure 2b). The isotropic Raman spectrum of neat water is dominated by the symmetric OH stretch that are coupled in-phase with the symmetric OH stretch modes of neighboring water molecules and that appear towards the lower-frequency region of the OH stretch band [49]. The band shape of the anisotropic Raman spectrum is similar to the unpolarized Raman spectrum of HOD and assignable to the antisymmetric OH stretch mode of water. Thus, the isotopic dilution and polarized Raman measurements show that the spectral response at ∼3250 cm −1 region is dominated by the symmetric OH stretch that are intra-and intermolecularly coupled along with FR. The band at ∼3420 cm −1 represents the average H-bond strength of bulk water (free from coupling effects). The extreme high-frequency shoulder around 3620 cm −1 mainly corresponds to the very weakly H-bonded water or the water-OH with broken H-bonds (free OH) that exists transiently during the switching of the H-bonded partner among water molecules [53,54]. Thus, at any instant, bulk water can be thought of a mixture of water species with continuously varying H-bond strengths.

Effect of Electrolytes on the OH Stretch Raman Cross-Section and Vibrational Coupling of Water
Coupling among the vibrational modes in liquid water has a close relation with its structure and properties. For instance, strong intramolecular vibrational coupling leads to ultrafast energy delocalization among different intramolecular vibrational modes of a water molecule; while strong intermolecular coupling does the same among several water molecules (collective vibration) [51,55]. These ultrafast elementary processes play pivotal roles in water-mediated energy dissipation processes. Figure 3a shows the effect of isotopic dilution on the OD stretch spectrum of D2O. With an increasing percentage of H2O in D2O, the intensity around 2350 cm −1 (corresponding to the 3200 cm −1 for the OH stretch) decreases and that around 2500 cm −1 increases, making the overall spectrum narrower and blue-shifted compared to D2O. These spectral changes suggest the reduction of Fermi resonance and vibrational decoupling in D2O due to the addition of H2O (isotopic dilution). Interestingly, the addition of electrolytes, such as NaBr, has a similar effect on the OD stretch of D2O (Figure 3b). Except NaF, all other Na-halides lead to such a spectral change in the OD (or OH) stretch regions. Thus, the spectral change in Figure 3a,b indicates the reduced coupling (VC and FR) of water in the presence of NaX (X = Cl, Br, I) [37,45]. The reduction of coupling could be due to the change in OH stretch frequency (variation in Hbond strength) and the preferred orientation of water in the hydration shell of the ions. Quantitative estimation, on the basis of integrated Raman intensity, revealed that the Raman cross section increases in the presence of the halide ions, except F − .

Effect of Electrolytes on the OH Stretch Raman Cross-Section and Vibrational Coupling of Water
Coupling among the vibrational modes in liquid water has a close relation with its structure and properties. For instance, strong intramolecular vibrational coupling leads to ultrafast energy delocalization among different intramolecular vibrational modes of a water molecule; while strong intermolecular coupling does the same among several water molecules (collective vibration) [51,55]. These ultrafast elementary processes play pivotal roles in water-mediated energy dissipation processes. Figure 3a shows the effect of isotopic dilution on the OD stretch spectrum of D2O. With an increasing percentage of H2O in D2O, the intensity around 2350 cm −1 (corresponding to the 3200 cm −1 for the OH stretch) decreases and that around 2500 cm −1 increases, making the overall spectrum narrower and blue-shifted compared to D2O. These spectral changes suggest the reduction of Fermi resonance and vibrational decoupling in D2O due to the addition of H2O (isotopic dilution). Interestingly, the addition of electrolytes, such as NaBr, has a similar effect on the OD stretch of D2O ( Figure 3b). Except NaF, all other Na-halides lead to such a spectral change in the OD (or OH) stretch regions. Thus, the spectral change in Figure 3a,b indicates the reduced coupling (VC and FR) of water in the presence of NaX (X = Cl, Br, I) [37,45]. The reduction of coupling could be due to the change in OH stretch frequency (variation in Hbond strength) and the preferred orientation of water in the hydration shell of the ions. Quantitative estimation, on the basis of integrated Raman intensity, revealed that the Raman cross section increases in the presence of the halide ions, except F − . Raman spectra of H2O and HOD (H2O:D2O = 1:9 (v/v)) in the OH stretch region; (b) )and anisotropic (Ianiso) Raman spectra of H2O. All spectra are peak-normalized at 3420 tropic and anisotropic Raman spectra are calculated as Iiso = I ∥ − 4/3I⊥ and Ianiso = I∥ and I⊥ are the parallel and perpendicular spectra recorded by placing a laminated r (in the detection path) in the parallel and perpendicular orientation with respect to ion of the excitation beam.

Electrolytes on the OH Stretch Raman Cross-Section and Vibrational Couater
ing among the vibrational modes in liquid water has a close relation with its d properties. For instance, strong intramolecular vibrational coupling leads to ergy delocalization among different intramolecular vibrational modes of a wae; while strong intermolecular coupling does the same among several water collective vibration) [51,55]. These ultrafast elementary processes play pivotal ter-mediated energy dissipation processes. Figure 3a shows the effect of isoon on the OD stretch spectrum of D2O.
With an increasing percentage of H2O intensity around 2350 cm −1 (corresponding to the 3200 cm −1 for the OH stretch) nd that around 2500 cm −1 increases, making the overall spectrum narrower and compared to D2O. These spectral changes suggest the reduction of Fermi resvibrational decoupling in D2O due to the addition of H2O (isotopic dilution). y, the addition of electrolytes, such as NaBr, has a similar effect on the OD 2O ( Figure 3b). Except NaF, all other Na-halides lead to such a spectral change or OH) stretch regions. Thus, the spectral change in Figure 3a,b indicates the upling (VC and FR) of water in the presence of NaX (X = Cl, Br, I) [37,45]. The f coupling could be due to the change in OH stretch frequency (variation in Hth) and the preferred orientation of water in the hydration shell of the ions. e estimation, on the basis of integrated Raman intensity, revealed that the Raection increases in the presence of the halide ions, except F − . man spectra (OD stretch) of (a) D2O-H2O mixtures with different mole ratios and (b) ferent mole fractions of NaBr. The spectra are normalized for the concentration of OD and I ⊥ are the parallel and perpendicular spectra recorded by placing laminated film polarizer (in the detection path) in the parallel and perpendicular orientation with respect to the polarization of the excitation beam.

Effect of Electrolytes on the OH Stretch Raman Cross-Section and Vibrational Coupling of Water
Coupling among the vibrational modes in liquid water has a close relation with its structure and properties. For instance, strong intramolecular vibrational coupling leads to ultrafast energy delocalization among different intramolecular vibrational modes of a water molecule; while strong intermolecular coupling does the same among several water molecules (collective vibration) [51,55]. These ultrafast elementary processes play pivotal roles in water-mediated energy dissipation processes. Figure 3a shows the effect of isotopic dilution on the OD stretch spectrum of D 2 O. With an increasing percentage of H 2 O in D 2 O, the intensity around 2350 cm −1 (corresponding to the 3200 cm −1 for the OH stretch) decreases and that around 2500 cm −1 increases, making the overall spectrum narrower and blue-shifted compared to D 2 O. These spectral changes suggest the reduction of Fermi resonance and vibrational decoupling in D 2 O due to the addition of H 2 O (isotopic dilution). Interestingly, the addition of electrolytes, such as NaBr, has a similar effect on the OD stretch of D 2 O (Figure 3b). Except NaF, all other Na-halides lead to such a spectral change in the OD (or OH) stretch regions. Thus, the spectral change in Figure 3a,b indicates the reduced coupling (VC and FR) of water in the presence of NaX (X = Cl, Br, I) [37,45]. The reduction of coupling could be due to the change in OH stretch frequency (variation in H-bond strength) and the preferred orientation of water in the hydration shell of the ions. Quantitative estimation, on the basis of integrated Raman intensity, revealed that the Raman cross section increases in the presence of the halide ions, except F − .

Effect of Electrolytes on the OH Stretch Raman Cross-Section and Vibrational Coupling of Water
Coupling among the vibrational modes in liquid water has a close relation with its structure and properties. For instance, strong intramolecular vibrational coupling leads to ultrafast energy delocalization among different intramolecular vibrational modes of a water molecule; while strong intermolecular coupling does the same among several water molecules (collective vibration) [51,55]. These ultrafast elementary processes play pivotal roles in water-mediated energy dissipation processes. Figure 3a shows the effect of isotopic dilution on the OD stretch spectrum of D2O. With an increasing percentage of H2O in D2O, the intensity around 2350 cm −1 (corresponding to the 3200 cm −1 for the OH stretch) decreases and that around 2500 cm −1 increases, making the overall spectrum narrower and blue-shifted compared to D2O. These spectral changes suggest the reduction of Fermi resonance and vibrational decoupling in D2O due to the addition of H2O (isotopic dilution). Interestingly, the addition of electrolytes, such as NaBr, has a similar effect on the OD stretch of D2O ( Figure 3b). Except NaF, all other Na-halides lead to such a spectral change in the OD (or OH) stretch regions. Thus, the spectral change in Figure 3a,b indicates the reduced coupling (VC and FR) of water in the presence of NaX (X = Cl, Br, I) [37,45]. The reduction of coupling could be due to the change in OH stretch frequency (variation in Hbond strength) and the preferred orientation of water in the hydration shell of the ions. Quantitative estimation, on the basis of integrated Raman intensity, revealed that the Raman cross section increases in the presence of the halide ions, except F − .

Water in the Hydration Shell of a Low-Charge-Density Monovalent Anion
Water interacts with anions through its hydrogen atom by forming an H-bond with the latter. The E-field around the ion reorients the water in the hydration shell such that the water-Hs are pointed 'towards' the hydrated anion and 'away' from the hydrated cation, following the law of electrostatics. Figure 4a shows the Raman spectrum of neat water and the MCR-extracted hydration shell water spectra of different monovalent anions. It is important to note that, for the hydration shell of low-charge-density anions, the band shapes of the MCR-extracted spectra are very similar to the Raman-DS-SCF-extracted spectra (not shown here). It is obvious that the spectral responses of hydration water are markedly different from the Raman spectrum of the bulk water. To examine the H-bond strength of water in the vicinity of ions, it is reasonable to compare the OH stretch spectra of hydration water with that of HOD, in which the intermolecular coupling of OH stretch modes is largely suppressed. The maxima of the OH stretch Raman band of the hydration water of the monovalent halide anions (Br − , I − ) lie in the higher-frequency region (blueshift) compare to that of HOD (Figure 4a), indicating weaker average H-bond strength of water in the hydration shell of the halide anions [36,39,45]. A similar spectral change was also observed for Cl − (spectrum not shown here) [39]. Moreover, the extent of the blue-shift of the hydration water spectra (OH stretch) follows the order of the increasing size of the anions, i.e., Cl − < Br − < I − . A triatomic monovalent anion with lower charge density than I − (e.g., I 3 − ) shows a further decrease in the average H-bond strength of the hydration water than that of the I − [56][57][58]. The monovalent NO 3 − anion also weakens the H-bond strength of its hydration shell water (Figure 4a) [36,44,59]. Further analysis of the integrated area of the 3250 cm −1 band with that of the H 2 O and HOD (decoupled OH) suggested a significant reduction (~80%) of the combined coupling effects (intermolecular coupling and Fermi resonance) in the vicinity of halide ions relative to the bulk water [36].

Water in the Hydration Shell of a Low-Charge-Density Monovalent Anion
Water interacts with anions through its hydrogen atom by forming an H-bond w the latter. The E-field around the ion reorients the water in the hydration shell such th the water-Hs are pointed 'towards' the hydrated anion and 'away' from the hydrat cation, following the law of electrostatics. Figure 4a shows the Raman spectrum of n water and the MCR-extracted hydration shell water spectra of different monovalent a ons. It is important to note that, for the hydration shell of low-charge-density anions, band shapes of the MCR-extracted spectra are very similar to the Raman-DS-SCFtracted spectra (not shown here). It is obvious that the spectral responses of hydrati water are markedly different from the Raman spectrum of the bulk water. To examine H-bond strength of water in the vicinity of ions, it is reasonable to compare the OH stre spectra of hydration water with that of HOD, in which the intermolecular coupling of O stretch modes is largely suppressed. The maxima of the OH stretch Raman band of hydration water of the monovalent halide anions (Br − , I − ) lie in the higher-frequency gion (blue-shift) compare to that of HOD (Figure 4a), indicating weaker average H-bo strength of water in the hydration shell of the halide anions [36,39,45]. A similar spect change was also observed for Cl − (spectrum not shown here) [39]. Moreover, the extent the blue-shift of the hydration water spectra (OH stretch) follows the order of the incre ing size of the anions, i.e., Cl − < Br − < I − . A triatomic monovalent anion with lower char density than I − (e.g., I3 − ) shows a further decrease in the average H-bond strength of hydration water than that of the I − [56][57][58]. The monovalent NO3 − anion also weakens H-bond strength of its hydration shell water (Figure 4a) [36,44,59]. Further analysis of integrated area of the 3250 cm −1 band with that of the H2O and HOD (decoupled O suggested a significant reduction (~80%) of the combined coupling effects (intermolecu coupling and Fermi resonance) in the vicinity of halide ions relative to the bulk water [3

Hydration Shell of Monovalent Oxyhalide Anion
Although the halide ions weaken the H-bond strength of hydration water in the order Cl − < Br − < I − , the oxyhalide ions affect the surrounding water differently. Figure 5a-c shows the hydration shell water spectra of ClO 3 − , BrO 3 − and IO 3 − anions in H 2 O. The band maximum is gradually blue-shifted (shifted to the high-frequency side) as the central halogen atom changes from I to Br to Cl. The relative intensity of the 3200 cm −1 shoulder band is decreased due to reduced vibrational coupling and the Fermi resonance of water in the hydration shell of the oxyhalide anions, as has been observed with the halide ions [36,45].
The hydration shell spectra in HOD (Figure 5d-f) show that the OH stretch band maximum is blue-shifted by~100 cm −1 and~70 cm −1 for ClO 3 − and BrO 3 − , respectively, while in the case of IO 3 − , the relative intensity in the low-frequency region (below 3400 cm −1 ) is enhanced. Hence, ClO 3 − and BrO 3 − weaken the H-bond strength of their hydration water following the order ClO 3 − > BrO 3 − , whereas IO 3 − is strengthened the same [60]. Using MD simulation and XAFS spectroscopy, Baer et al. showed that the central iodine atom of IO 3 − possesses a positive charge, while each of the three oxygen atoms bears a negative charge [61], which indicates two distinct hydration environments around the iodate ion corresponding to the cationic iodine and anionic oxygen atoms, respectively. MD simulation also revealed that the vibrational power spectra and distributions of different types of H-bonded water molecules hydrating the iodate-oxygens are very similar to that of bulk water [62]. Hence, the cationic iodine (of iodate) predominantly affects the structure of surrounding water. Thus, unlike the weakly interacting I − anion, in spite of having an overall negative charge, IO 3 − effectively behaves as a hydrated cation in aqueous solution which makes the H-bonds of surrounding water stronger. The asymmetric hydration of IO 3 − in aqueous KIO 3 solution has also been suggested by Raman-MCR study [56]. The MD simulation of BrO 3 − ion in water suggested that anionic oxygens (of BrO 3 − ) affect the hydrating water to a greater extent than the cationic bromine (of BrO 3 − ) and that the average strength of BrO 3 − -water H-bonding is weaker than that of water-water in bulk [63].

Hydration Shell of High-Charge-Density Anions
For the high-charge-density anions, the Na-salt of F − and CO3 2− anions has been lyzed with respect to the bulk water using Raman-DS-SCF analysis. The hydration w spectrum of F − shows enhanced absorption (compared to neat water) in the low-frequ

Hydration Shell of High-Charge-Density Anions
For the high-charge-density anions, the Na-salt of F − and CO 3 2− anions has been analyzed with respect to the bulk water using Raman-DS-SCF analysis. The hydration water spectrum of F − shows enhanced absorption (compared to neat water) in the lowfrequency region of the OH stretch up to 2600 cm −1 (Figure 6) [64]. This indicates the strong H-bonding of water with the F − ion. MD simulation also suggested that the reorientation dynamics of water in the first hydration shell of F − is much slower than that of I − , invoking stronger F − . . . H 2 O H-bonding [65]. Divalent CO 3 2− anion perturbs the H-bond strength of hydration water quite similar to that of F − , although the average number of water molecules perturbed per CO 3 2− anion is greater than that of F − (Figure 6) [64]. In the hydration shell Raman spectrum of CO 3 2− and F − ions in H 2 O, a new band appears in the lower-frequency region (~3050 cm −1 ) suggesting the presence of a very strong H-bond in hydration water [39,64]. Thus, unlike the low-charge-density monovalent ions (Cl − , Br − , I − , ClO 3 − , BrO 3 − ), higher-charge-density anions, such as F − and CO 3 2− , increase the average H-bond strength of water in their hydration shell compared to that of bulk water. The dip feature around 3200 cm −1 in the hydration shell spectra in H 2 O is due to the decoupling of the vibrational modes of water in the hydration shell of the ions; the dip feature is largely reduced in HOD. Hence, the structure of water in the hydration shell of CO 3 2− and F − is more inhomogeneous compared to bulk water. It is believed that two different types of water molecules exist in the hydration shell of the ions: water in the first hydration shell forms a very strong H-bond (~3050 cm −1 ) with the carbonate oxygen (or F − ), and water in the second hydration shell has an H-bond strength weaker than the first hydration shell but stronger than bulk water. In agreement with this, the calculated (MD simulation) stretch frequency of the water-OD, directly H-bonded to the CO 3 2− oxygen, is red-shifted (i.e., shifted towards low frequency region) from that of bulk-water-OD, whereas the stretch frequency of the second hydration shell water lies in between that of the bulk and first hydration shell water [66]. These observations certainly indicate that the CO 3 2− anion has a long-range effect on neighboring water molecules that goes beyond the first hydration shell.
In the case of SO 4 2− , a distinct band around 3250 cm −1 is clearly visible is the hydration shell spectrum (in HOD; Figure 6e), though the band is not as intense as that with CO 3 2− . In agreement with the OH stretch, the analysis of the HOH bending band revealed stronger H-bonding in the hydration shell of SO 4 2− [39]. The difference spectroscopy-based IR methods in the OD stretch region also suggested that water forms H-bonds with sulfate ions stronger than that between two water molecules in bulk water [44,59]. Kameda et al. showed that the intermolecular H-bonded O-O distance between sulfate and water (in aqueous sulfuric acid solution) is shorter than that of pure water, suggesting a strongly H-bonded intermolecular interaction [67].
As expected, the trivalent PO 4 3− anion induces the strong H-bonding of water in the hydration shell [36]. The Raman spectrum of aqueous PO 4 3− solution in the lowwavenumber region (hindered translation of the H-bonded water molecules or the Hbond stretch) also confirmed the formation of strong P-O . . . H bonds in solution [68]. One interesting feature of the hydration shell water spectrum is that, unlike the case of other anions, the shoulder near 3600 cm −1 , which represents weakly interacting or undercoordinated water-OH, is not reduced for PO 4 3− [36]. Furthermore, Raman measurements of aqueous K 3 PO 4 solutions showed a polarized shoulder band at~3620 cm −1 , which increased slightly in intensity with increasing concentration (range 0.7-5.3 M) of the salt [68]. Chandra and co-worker calculated the OD stretch frequencies in the hydration shells of PO 4 3− [69]. They showed that the stretch frequency of OD in the first hydration shell that directly forms an H-bond with the anion is considerably red-shifted compared to the bulk water; whereas the other OD (of the same water forming an H-bond with the anion) stretch is blue-shifted and has significant overlap with the band of non-H-bonded OD in bulk water. On the basis of this observation, the authors speculated the existence of weakly H-bonded (compared to water-water in bulk) water-OD in the first hydration shell of the PO 4 3− anion.
hydration shell has an H-bond strength weaker than the first hydration shell but stronger than bulk water. In agreement with this, the calculated (MD simulation) stretch frequency of the water-OD, directly H-bonded to the CO3 2− oxygen, is red-shifted (i.e., shifted towards low frequency region) from that of bulk-water-OD, whereas the stretch frequency of the second hydration shell water lies in between that of the bulk and first hydration shell water [66]. These observations certainly indicate that the CO3 2− anion has a longrange effect on neighboring water molecules that goes beyond the first hydration shell. In the case of SO4 2− , a distinct band around 3250 cm −1 is clearly visible is the hydration shell spectrum (in HOD; Figure 6e), though the band is not as intense as that with CO3 2− . In agreement with the OH stretch, the analysis of the HOH bending band revealed stronger H-bonding in the hydration shell of SO4 2− [39]. The difference spectroscopy-based IR methods in the OD stretch region also suggested that water forms H-bonds with sulfate ions stronger than that between two water molecules in bulk water [44,59]. Kameda et al. showed that the intermolecular H-bonded O-O distance between sulfate and water (in aqueous sulfuric acid solution) is shorter than that of pure water, suggesting a strongly H-bonded intermolecular interaction [67].
As expected, the trivalent PO4 3− anion induces the strong H-bonding of water in the hydration shell [36]. The Raman spectrum of aqueous PO4 3− solution in the low-wavenumber region (hindered translation of the H-bonded water molecules or the H-bond stretch) also confirmed the formation of strong P-O…H bonds in solution [68]. One interesting feature of the hydration shell water spectrum is that, unlike the case of other anions, the shoulder near 3600 cm −1 , which represents weakly interacting or under-coordinated water-OH, is not reduced for PO4 3− [36]. Furthermore, Raman measurements of aqueous K3PO4 solutions showed a polarized shoulder band at ~3620 cm −1 , which increased slightly in intensity with increasing concentration (range 0.7-5.3 M) of the salt [68]. Chandra and Copyright 2021 American Chemical Society).

HOH Bend and Bend+Libration Combination Band of Water in the Anion Hydration Shell
The restricted rotational motions of water in the H-bonded network are termed as librational motions that play an important role in the energy transfer mechanism in liquid water. It has been suggested that, unlike fast intermolecular stretch-to-stretch energy transfer in water, intermolecular bend-to-bend transfer is much slower [70]. Rather, librational motions randomize the orientation of the HOH bending mode transition dipole moment much faster, and thus, the vibrational energy of the HOH bending mode is primarily released to the librational modes.
In the Raman spectrum of neat water, the HOH bend vibration of water appears at 1640 cm −1 [71]. As the H-bond strength of water increases, the HOH bending motion becomes more restricted. Therefore, the HOH bending band of water is expected to be blueshifted (red-shifted) in the hydration shell of ions, where the water H-bond is strengthened (weakened). Upon isotopic dilution, the HOH bend maximum is blue-shifted by ∼10 cm −1 and the spectral width is decreased due to vibrational decoupling (compare the black dotted and green dashed curves in Figure 7a,b). As can be seen in Figure 7a,b, the hydration shell spectra for Cl − , I − and SO 4 2− ions are blue-shifted and have lower spectral width than that of neat H 2 O, indicating that the HOH bending mode of water in the hydration shell of these anions are decoupled from the neighboring water molecules. On the contrary, the hydration shell spectrum for CO 3 2− showed comparable spectral width to that of bulk H 2 O (Figure 7b). It is likely that the larger spectral width (compared to that of other anions) for CO 3 2− is because of the greater distribution of H-bond strength in the hydration shell.

Water in the Hydration Shell of a Metal Cation
The local electric field around a cation orients hydration water dipole vectors to po radially away from the cation. The monovalent cations with low charge density (e.g., N Cs + ) exhibit a negligible effect on the H-bonded structure of the hydration water compar to the monovalent anions. Although monovalent Na + does not have significant influen on the neighboring water [76], the cations with higher charge density are expected to p turb the water in their hydration shell. In fact, the DDIR study of divalent transition me ions (Mn 2+ , Fe 2+ , Co 2+ ) and trivalent lanthanide ions (La 3+ , Nd 3+ , Dy 3+ , Yb 3+ ) with compara ionic potential revealed quite a similar effect on the OD stretch band [77]. Rudolph a co-workers investigated a series of lanthanide ions (Ln 3+ : La 3+ , Ce 3+ , Pr 3+ , Nd 3+ , Sm 3+ , Ho Er 3+ , Tm 3+ , Yb 3+ and Lu 3+ ) using Raman measurements that enabled determining the fo constant of the Ln-O breathing modes of the aqua complexes [78,79]. It was observed th the strength of the force constant increases with decreasing Ln-O bond distances (bo distances: Recently, we have selectively extracted the OH stretch spectra of water in the hyd tion shell of a series of high-charge-density cations (Li + , Mg 2+ , Ca 2+ , Sr 2+ , Ba 2+ , La 3+ , Gd Dy 3+ and Lu 3+ , as well as H + ) using Raman-DS-SCF spectroscopy [42]. It was observed th these cations strongly interact with their hydration water and profoundly affect the lo water structure. Interestingly, the strong cation-water interaction not only increases t average H-bond strength of hydration water but also produces a small fraction of tremely weakly interacting water-OH in the hydration shell, which increases the stru tural heterogeneity therein (see Figure 8). Figure 8(a) and (b) compare the hydration sh spectra of monovalent Cs + (low charge density) to trivalent Lu 3+ (high charge density), Water librations represent a broad band over 300-1000 cm −1 in the Raman spectrum of neat water [8,72]. In fact, this broad librational response is composed of three overlapped bands (centered at~450,~580 and~720 cm −1 ) arising from water librations about the three axes of rotation [72][73][74]. The librational band of water is red-shifted with increasing temperature signifying increased librational freedom because of the weakened H-bonding of water at elevated temperatures [75]. Hence, it is expected that the ions that induce the weakening (strengthening) of the H-bond strength of neighboring water will increase (decrease) the librational freedom of the hydration water. The bend+libration combination band of hydration water is red-shifted for monovalent Cl − and I − compared to the neat water (~2130 cm −1 ), while it is blue-shifted for divalent SO 4 2− and CO 3 2− (Figure 7c,d). The calculated librational band positions revealed that water molecules in the hydration shell of the monovalent Cl − and I − weakly interact with vicinal water and have higher librational freedom than bulk water [39]. Conversely, water in the immediate vicinity of divalent (strongly interacting) anions such as SO 4 2− and CO 3 2− have more restricted librational mobility compared to bulk water [39].

Water in the Hydration Shell of a Metal Cation
The local electric field around a cation orients hydration water dipole vectors to point radially away from the cation. The monovalent cations with low charge density (e.g., Na + , Cs + ) exhibit a negligible effect on the H-bonded structure of the hydration water compared to the monovalent anions. Although monovalent Na + does not have significant influence on the neighboring water [76], the cations with higher charge density are expected to perturb the water in their hydration shell. In fact, the DDIR study of divalent transition metal ions (Mn 2+ , Fe 2+ , Co 2+ ) and trivalent lanthanide ions (La 3+ , Nd 3+ , Dy 3+ , Yb 3+ ) with comparable ionic potential revealed quite a similar effect on the OD stretch band [77].
Rudolph and co-workers investigated a series of lanthanide ions (Ln 3+ : La 3+ , Ce 3+ , Pr 3+ , Nd 3+ , Sm 3+ , Ho 3+ , Er 3+ , Tm 3+ , Yb 3+ and Lu 3+ ) using Raman measurements that enabled determining the force constant of the Ln-O breathing modes of the aqua complexes [78,79]. It was observed that the strength of the force constant increases with decreasing Ln-O bond distances (bond distances: Recently, we have selectively extracted the OH stretch spectra of water in the hydration shell of a series of high-charge-density cations (Li + , Mg 2+ , Ca 2+ , Sr 2+ , Ba 2+ , La 3+ , Gd 3+ , Dy 3+ and Lu 3+ , as well as H + ) using Raman-DS-SCF spectroscopy [42]. It was observed that these cations strongly interact with their hydration water and profoundly affect the local water structure. Interestingly, the strong cation-water interaction not only increases the average H-bond strength of hydration water but also produces a small fraction of extremely weakly interacting water-OH in the hydration shell, which increases the structural heterogeneity therein (see Figure 8). Figure 8a,b compare the hydration shell spectra of monovalent Cs + (low charge density) to trivalent Lu 3+ (high charge density), as extracted by the Raman-DS-SCF method. The hydration shell spectrum (OH stretch band-shape) of Cs + is qualitatively similar to that of bulk water. However, the hydration shell spectrum of Lu 3+ shows distinct changes compared to that of bulk water: the spectral intensity is enhanced in the low-frequency region (below 3400 cm −1 ; strong H-bonding) and a small band appears at the high-frequency end near 3600 cm −1 (weakly interacting water-OH). The extremely weakly interacting water-OH is believed to have originated from the second hydration shell water-OH that donates an H-bond to the electron-deficient water-oxygen in the first hydration shell of the high-charge-density cations and distributed unevenly within the range of ionic influence (see Figure 9). In fact, this weakly interacting water-OH is different from the broken H-bonded OH transients (usually appearing at ∼3620 cm −1 ) that is uniformly distributed in bulk water [54]. Hence, the high-charge-density metal ions induce the structural perturbation of water beyond their first hydration shell, which was also suggested by MD simulation studies [80,81]. The number of affected water molecules by the metal cations is correlated with the ionic potential of the ions: higher the ionic potential, the greater the number of affected water molecules by the cations. Interestingly, among the lanthanides, the number of affected water molecules marginally decreases with increasing ionic potential. On moving from La 3+ to Lu 3+ , the number of affected water molecules decreases following a similar trend as in the coordination numbers due to lanthanide contraction [82][83][84]. Nevertheless, the effects of cations on water cannot be described solely by the ionic potential. The ionic potential of Mg 2+ is higher than that of La 3+ ; however, the latter affects a greater number of water molecules. In addition to the ionic potential, the absolute charge (i.e., oxidation state) and the size of an ion play roles in the ion's effect on water. Between two ions with the same ionic potential, the larger one affects a greater number of water molecules than the smaller one. Mg 2+ and Dy 3+ have the same ionic potential, but the latter ion affects a greater number of water molecules [42].
The perturbing influence of a proton (H + ) on the structure of water is very similar to that of high-charge-density metal ions [42,43]. Indeed, H + and a high-charge-density metal ion are similar from the perspective of interaction with hydration water: both of them accept electrons from the oxygen of the water in the first hydration shell.
larger one affects a greater number of water molecules than the smaller one. Mg 2+ and Dy 3+ have the same ionic potential, but the latter ion affects a greater number of water molecules [42].
The perturbing influence of a proton (H + ) on the structure of water is very similar to that of high-charge-density metal ions [42,43]. Indeed, H + and a high-charge-density metal ion are similar from the perspective of interaction with hydration water: both of them accept electrons from the oxygen of the water in the first hydration shell.  Liquids 2022, 2, FOR PEER REVIEW 14 Figure 9. Schematic of the hydration shell of a high-charge-density cation: dashed white lines represent a very weak H-bond; the corresponding parent OH appears at ∼3600 cm −1 in the hydration shell spectrum.

Water Shared Ion-Pair Formation and Restructuring of Hydration Shell Water
In a very dilute aqueous solution of an electrolyte, the cation and anion are hydrated separately and mostly do not influence the hydration shell structure of each other. In such a fully hydrated ion, the structural and dynamical properties of the hydrating water are perturbed by the physicochemical properties of the individual hydrated ion. However, with the increasing concentration of the electrolyte, the oppositely charged ions approach closer with their native hydration shells, which can mutually influence their properties [13,85,86]. In such cases, the native hydration shells of the oppositely charged ions overlap, and the intervening (shared) water becomes perturbed by both the interacting ions [87]. These ion-pairs are generally termed "solvent-shared ion-pairs (SSIP)". The structural perturbation of water in hydrated ion-pairs is important in understanding the ionic mobilities and thermodynamics of solutions.
Raman-DS-SCF analysis successfully enabled the identification of the formation of SSIPs and suitably extracted the vibrational response of the shared water. As can be seen in Figure 10a, the band shape of the OH stretch hydration shell water spectra (extracted by Raman-DS-SCF) of Mg 2+ in its aqueous chloride solution undergoes gradual changes with the increasing concentration of MgCl2 ([MgCl2] ≥ 0.37 M). The apparent band maxima of the Mg 2+ hydration shell spectrum is shifted from ~3200 cm −1 to ~3400 cm −1 (as indicated by the gray dashed arrow in Figure 10a). This concentration-dependent spectral change indicates an interionic interaction between Mg 2+ and Cl − , forming SSIP. Further, a secondround DS-SCF analysis (Figure 10b) revealed that the shared water is affected by both Figure 9. Schematic of the hydration shell of a high-charge-density cation: dashed white lines represent a very weak H-bond; the corresponding parent OH appears at ∼3600 cm −1 in the hydration shell spectrum.

Water Shared Ion-Pair Formation and Restructuring of Hydration Shell Water
In a very dilute aqueous solution of an electrolyte, the cation and anion are hydrated separately and mostly do not influence the hydration shell structure of each other. In such a fully hydrated ion, the structural and dynamical properties of the hydrating water are perturbed by the physicochemical properties of the individual hydrated ion. However, with the increasing concentration of the electrolyte, the oppositely charged ions approach closer with their native hydration shells, which can mutually influence their properties [13,85,86]. In such cases, the native hydration shells of the oppositely charged ions overlap, and the intervening (shared) water becomes perturbed by both the interacting ions [87]. These ionpairs are generally termed "solvent-shared ion-pairs (SSIP)". The structural perturbation of water in hydrated ion-pairs is important in understanding the ionic mobilities and thermodynamics of solutions.
Raman-DS-SCF analysis successfully enabled the identification of the formation of SSIPs and suitably extracted the vibrational response of the shared water. As can be seen in Figure 10a, the band shape of the OH stretch hydration shell water spectra (extracted by Raman-DS-SCF) of Mg 2+ in its aqueous chloride solution undergoes gradual changes with the increasing concentration of MgCl 2 ([MgCl 2 ] ≥ 0.37 M). The apparent band maxima of the Mg 2+ hydration shell spectrum is shifted from~3200 cm −1 to~3400 cm −1 (as indicated by the gray dashed arrow in Figure 10a). This concentration-dependent spectral change indicates an interionic interaction between Mg 2+ and Cl − , forming SSIP. Further, a secondround DS-SCF analysis (Figure 10b) revealed that the shared water is affected by both Mg 2+ and Cl − ; the average H-bonding of the shared water is weaker than the hydration shell water of Mg 2+ but stronger than the hydration shell water of Cl − . A similar observation was also observed for the case of an aqueous LaCl 3 solution [88]. The perturbing influences of both the cation and anion on the shared water have also been suggested in the aqueous solutions of nitrate and the chloride salt of Fe(III) [89].
Liquids 2022, 2, FOR PEER REVIEW 15 between Na + and Cl − (above 3M) [88]. On the contrary, in the cases of Mg 2+ -and La 3+ -salts, the predominant formation of SSIP was identified (data not shown here), in agreement with THz and dielectric relaxation spectroscopy studies [29,90]. X-ray absorption fine structure spectroscopy and neutron scattering investigations of divalent Ca 2+ (chloride salt) solvation provided evidence for the predominant formation of Ca 2+ -OH2-Cl − SSIPs compared to Ca 2+ -Cl − CIPs, even at high (~6 M) concentrations [91,92]. Presumably, the strong interaction (i.e., tight binding) of first hydration shell water with high-charge-density cations (e.g., Mg 2+ and La 3+ ) promotes the preferential formation of SSIP, while weaker ion-water interaction favors the formation of CIP in NaCl solution. Ion-pairing in aqueous solution has also been recently detected using Raman-MCR and IR-MCR analyses [93,94]. The analyses of the band shape (OH stretch) and intensity of the Raman-MCR extracted ion-perturbed water spectra of aqueous NaOH and LiOH solutions revealed that SSIPs are predominantly formed in these solutions, while CIPs only become significant at higher concentrations [93]. In fact, the strong H-bonding between OH − and the first hydration shell water (HOH…OH − ) favors the predominant formation of SSIP. Moreover, the ion-pairing in aqueous NaOH and LiOH is enthalpically disfavored and entropically favored (i.e., positive ΔH and ΔS), implying that the thermodynamics of ion-pairing is dominated by ion-water interaction rather the electrostatic attraction between the oppositely charged ions.

Water in the Hydration Shell of Hydrophobic Molecular Cation ((CH3)4N + )
Ions containing hydrophobic molecular groups exhibit distinct local effects on the neighboring water molecules. The hydration shell spectrum (OH stretch) of tetramethylammonium cation ((CH3)4N + ; TMA + ) closely resembles that of tert-butyl alcohol (TBA; an uncharged hydrophobe), showing "hydrophobic hydration" around the tetramethyl group ( Figure 11). Water lacks favorable interactions with non-polar hydrophobes and is unable to form H-bonds. Accordingly, water molecules around a hydrophobe reorganize themselves to optimize the water-water interaction, resulting in increased average Hbond strength (i.e., the enhanced intensity ~3200 cm −1 ; compare the green and black dashed curves of Figure 11a and 11b), along with some non-H-bonded "dangling OH" defects (small band ~3670 cm −1 in case of TBA hydration shell) [95,96]. Notably, the dangling OH defect of hydrophobic hydration in TMA + is largely suppressed due to the presence of a cationic charge. The ion-dipole interaction between the positively charged TMA + The Raman-DS-SCF analysis also enabled the identification of the type of ion-pair formed in the solution from the band shape and the intensity of the hydration shell water spectra. In the formation of SSIP, a fraction of the hydration shell water molecules of an ion is shared by its counterion. Hence, the average number of water molecules affected by the ion is not expected to vary significantly. As a result, the integrated intensity of the hydration shell water spectra that is a quantitative measure of the ion-affected water varies linearly with concentration. On the contrary, a fraction of water molecules is excluded from the hydration shell of an ion during the formation of contact ion pairs (CIP) where the ion directly interacts with its counterion. Thus, the average number of affected water molecules is expected to decrease during CIP formation. The integrated intensity of the Cl − -affected water spectrum (OH stretch) linearly increases with concentrations up to 3 M NaCl while showing negative deviation above 3 M NaCl, indicating the formation of CIP between Na + and Cl − (above 3M) [88]. On the contrary, in the cases of Mg 2+ -and La 3+ -salts, the predominant formation of SSIP was identified (data not shown here), in agreement with THz and dielectric relaxation spectroscopy studies [29,90]. X-ray absorption fine structure spectroscopy and neutron scattering investigations of divalent Ca 2+ (chloride salt) solvation provided evidence for the predominant formation of Ca 2+ -OH 2 -Cl − SSIPs compared to Ca 2+ -Cl − CIPs, even at high (~6 M) concentrations [91,92]. Presumably, the strong interaction (i.e., tight binding) of first hydration shell water with high-charge-density cations (e.g., Mg 2+ and La 3+ ) promotes the preferential formation of SSIP, while weaker ion-water interaction favors the formation of CIP in NaCl solution.
Ion-pairing in aqueous solution has also been recently detected using Raman-MCR and IR-MCR analyses [93,94]. The analyses of the band shape (OH stretch) and intensity of the Raman-MCR extracted ion-perturbed water spectra of aqueous NaOH and LiOH solutions revealed that SSIPs are predominantly formed in these solutions, while CIPs only become significant at higher concentrations [93]. In fact, the strong H-bonding between OH − and the first hydration shell water (HOH . . . OH − ) favors the predominant formation of SSIP. Moreover, the ion-pairing in aqueous NaOH and LiOH is enthalpically disfavored and entropically favored (i.e., positive ∆H and ∆S), implying that the thermodynamics of ion-pairing is dominated by ion-water interaction rather the electrostatic attraction between the oppositely charged ions.

Water in the Hydration Shell of Hydrophobic Molecular Cation ((CH 3 ) 4 N + )
Ions containing hydrophobic molecular groups exhibit distinct local effects on the neighboring water molecules. The hydration shell spectrum (OH stretch) of tetramethylammonium cation ((CH 3 ) 4 N + ; TMA + ) closely resembles that of tert-butyl alcohol (TBA; an uncharged hydrophobe), showing "hydrophobic hydration" around the tetramethyl group ( Figure 11). Water lacks favorable interactions with non-polar hydrophobes and is unable to form H-bonds. Accordingly, water molecules around a hydrophobe reorganize themselves to optimize the water-water interaction, resulting in increased average H-bond strength (i.e., the enhanced intensity~3200 cm −1 ; compare the green and black dashed curves of Figure 11a,b), along with some non-H-bonded "dangling OH" defects (small band~3670 cm −1 in case of TBA hydration shell) [95,96]. Notably, the dangling OH defect of hydrophobic hydration in TMA + is largely suppressed due to the presence of a cationic charge. The ion-dipole interaction between the positively charged TMA + and hydration water preferentially orients the latter, pointing their hydrogens away from the TMA + surface, suppressing the propensity of dangling OH (3670 cm −1 band is absent in TMA + hydration shell).
Liquids 2022, 2, FOR PEER REVIEW 16 and hydration water preferentially orients the latter, pointing their hydrogens away from the TMA + surface, suppressing the propensity of dangling OH (3670 cm −1 band is absent in TMA + hydration shell).

Water in the Hydration Shell of Zwitterion, Trimethylamine N-Oxide (CH3)3N + -O − )
In the case of zwitterions, such as trimethylamine N-oxide (TMAO), which contains a hydrophobic trimethyl and an zwitterionic N-oxide (N + -O − ) group, shows different local hydration characteristics: hydrophobic hydration around the trimethyl group and hydrophilic (ionic) hydration of the N-oxide group (Figure 12) [97]. In addition to hydrophobic hydration, the TMAO-oxygen forms a stronger H-bond with the hydration water compared to that in bulk water (reflected by the markedly enhanced intensity ~3230 cm −1 in the TMAO hydration shell spectrum compared to both the TBA and the TMA + hydration shell; compare the green curves of Figures 11 and 12). In fact, a balance between these differential hydration characteristics (hydrophobic vs. ionic) of TMAO is believed to have implications in its interfacial properties and its interactions with a biological membrane/water interface [98,99]. Moreover, the hydration of such an amphoteric molecule is very sensitive to the local physicochemical environment. In aqueous solution, for example, TMAO experiences van der Waals attraction with TBA and hence, interacts through its trimethyl moiety. This is governed by their mutual hydrophobic interaction with water, which makes the trimethyl group of TMAO less exposed to water [100]. On the other hand, TMAO electrostatically interacts with cationic TMA + through its negatively charged + Figure 11. Hydration shell spectra (OH stretch; retrieved by Raman-MCR) of (a) TBA and (b) TMA + in aqueous solution (0.8 M). The spectrum of bulk water (black dashed curve) is shown in each panel for reference. (Reprinted with permission from ref. [97]. Copyright 2016 American Chemical Society).

Water in the Hydration Shell of Zwitterion, Trimethylamine N-Oxide (CH 3 ) 3 N + -O − )
In the case of zwitterions, such as trimethylamine N-oxide (TMAO), which contains a hydrophobic trimethyl and an zwitterionic N-oxide (N + -O − ) group, shows different local hydration characteristics: hydrophobic hydration around the trimethyl group and hydrophilic (ionic) hydration of the N-oxide group (Figure 12) [97]. In addition to hydrophobic hydration, the TMAO-oxygen forms a stronger H-bond with the hydration water compared to that in bulk water (reflected by the markedly enhanced intensity~3230 cm −1 in the TMAO hydration shell spectrum compared to both the TBA and the TMA + hydration shell; compare the green curves of Figures 11 and 12). In fact, a balance between these differential hydration characteristics (hydrophobic vs. ionic) of TMAO is believed to have implications in its interfacial properties and its interactions with a biological membrane/water interface [98,99]. Moreover, the hydration of such an amphoteric molecule is very sensitive to the local physicochemical environment. In aqueous solution, for example, TMAO experiences van der Waals attraction with TBA and hence, interacts through its trimethyl moiety. This is governed by their mutual hydrophobic interaction with water, which makes the trimethyl group of TMAO less exposed to water [100]. On the other hand, TMAO electrostatically interacts with cationic TMA + through its negatively charged oxygen, which orients the trimethyl group of TMAO away from TMA + and renders them exposed to water [100]. pared to that in bulk water (reflected by the markedly enhanced intensity ~3230 cm in the TMAO hydration shell spectrum compared to both the TBA and the TMA + hydration shell; compare the green curves of Figures 11 and 12). In fact, a balance between these differential hydration characteristics (hydrophobic vs. ionic) of TMAO is believed to have implications in its interfacial properties and its interactions with a biological membrane/water interface [98,99]. Moreover, the hydration of such an amphoteric molecule is very sensitive to the local physicochemical environment. In aqueous solution, for example, TMAO experiences van der Waals attraction with TBA and hence, interacts through its trimethyl moiety. This is governed by their mutual hydrophobic interaction with water, which makes the trimethyl group of TMAO less exposed to water [100]. On the other hand, TMAO electrostatically interacts with cationic TMA + through its negatively charged oxygen, which orients the trimethyl group of TMAO away from TMA + and renders them exposed to water [100].

Conclusions and Future Perspective
Anions and cations exert specific perturbation on the H-bonded structure of water. The strength of H-bonding and water libration are altered in the immediate vicinity of the ions, causing the decoupling of the intra-and intermolecular vibrational modes of water. The ion effect on water increases with increasing charge density for both anions as well as cations. Monovalent anions have relatively stronger influences on the hydration water compared to the monovalent cations with a similar charge density. Nevertheless, the multivalent cations strongly interact with water that not only increases the average H-bond strength of hydration water but also produces very weakly interacting water-OH beyond the first hydration shell. In the case of fully hydrated ions, the structure of hydration water becomes perturbed primarily by the individual hydrated ion, while at higher concentrations, interionic interaction between anions and cations starts to form ion-pairs (SSIP and CIP). In SSIP, the intervening shared water is perturbed by both ions, forming an ion-pair. In fact, the predominant formation of SSIP or CIP is driven by the specific nature of the individual ion-water interaction. Hydrophobic molecular ions exhibit distinct local hydration characteristics.
Ion-induced perturbation modifies coupling between vibrational modes of water. Anion-water interaction decouples the vibrational modes of the hydration water from its H-bonded partner. The vibrational coupling of water plays crucial roles in the biochemical reactions of enzymes [101,102]. Quantification of the vibrational coupling of perturbed water by metal ions with biochemical relevance is expected to provide a better understanding of such systems.
The majority of the IR and Raman spectroscopic investigations have used the OH stretch band of water as a marker of the ion effect on the water structure. Presumably, this is because of the high oscillator strength and the higher sensitivity of the OH stretch towards a structural perturbation of the water than that of its other vibrational modes. However, recent FTIR spectroscopic measurements indicated the potential of HOH bending and bend+libration combination band of water as powerful probes for solute-induced perturbations of a water H-bond network [103][104][105]. Moreover, probing low-frequency intermolecular H-bond stretch (~200 cm −1 ) and librational modes can decipher the structural perturbation of water by ions (data not yet published). The low-frequency vibrational modes of water can be applied as complementary probes to the widely studied OH stretch mode in elucidating the structure of water.
One emerging trend is the combined investigation of ionic hydration in bulk as well as at aqueous interface. It is observed that the ion-induced perturbations of H-bonding structure in the bulk water and at the air/water interface are strongly correlated [56,60,64]. The local hydration characteristics of ions in bulk water dictate their surface affinity [56]. In fact, kosmotropic electrolytes perturb the structure of water at the air/water-electrolyte interface through their strongly solvated anion hydration shell [64]. DS-SCF analysis of the vibrational response from bulk water and the aqueous interface can provide a unified picture of a specific ion effect in bulk water and at the aqueous interfaces.