Temperature dependence of water-water and ion-water correlations in bulk water and electrolyte solutions probed by femtosecond elastic second harmonic scattering

The temperature dependence of the femtosecond elastic second harmonic scattering (fs-ESHS) response of bulk light and heavy water and their electrolyte solutions is presented. We observe clear temperature dependent changes in the hydrogen (H)-bond network of water that show a decrease in the orientational order of water with increasing temperature. Although D2O has a more structured H-bond network (giving rise to more fs-ESHS intensity), the relative temperature dependence is larger in H2O. The changes are interpreted in terms of the symmetry of H-bonds and are indicators of nuclear quantum effects. Increasing the temperature in electrolyte solutions decreases the influence of the total electrostatic field from ions on the water-water correlations, as expected from Debye-Hückel theory, since the Debye length becomes longer. The effects are, however, 1.9 times (6.3 times) larger than those predicted for H2O (D2O). Since fs-ESHS responses can be computed from known molecular coordinates, our observations provide a unique opportunity to refine quantum mechanical models of water.


INTRODUCTION
Water in its liquid form is the medium of life. In organisms, the structural building blocks of cells (lipid membranes, proteins, DNA, etc.) are embedded in electrolyte solutions. The structure of water and the hydration of ions have continuously been of high interest as they are involved in numerous physiological, medical, biological, and chemical processes. 1 Examples are protein folding/precipitation, enzyme, ion channel, and ion pump activity, action and membrane potential generation, transport across membranes, self-assembly, interfacial charging, and aerosol formation. [1][2][3][4][5][6] Most studies of ion-water interactions have focused on concentrations above 0.1 M, and from these studies it has emerged that ion specific water-ion interactions occur over up to 3 hydration shells. [7][8][9][10][11][12][13][14][15][16][17] Recently, femtosecond elastic second harmonic scattering (fs-ESHS) experiments, in combination with reflection second harmonic generation and surface tension measurements, 18 have shown that there is also an interaction that ranges over tens of hydrations shells and displays significant nuclear quantum effects. [19][20][21][22][23][24][25][26] This effect involves the interaction of the combined electrostatic field of all the ions with the hydrogen (H)-bond network in solution. 18,21,27 The orientational correlations or orientational order between the water molecules increases in response to an ionic electrostatic field, leading to an increase in the fs-ESHS response (as in a) Y. Chen and N. Dupertuis contributed equally to this work. b) sylvie.roke@epfl.ch Fig. 1), a decrease in the surface tension, and an increase in the dielectric energy of the bulk medium. 28 Although the qualitative ingredients of this behavior are understood, there are still many open questions. In the present work, we provide a further experimental investigation into the temperature dependence of the fs-ESHS response of bulk light and heavy water and their electrolyte solutions.
With fs-ESHS, we probe the second harmonic (SH) response of liquid molecules. The measured SH intensity originates from electronic anisotropy within every single nonspherical molecule and is modified by orientational correlations of the non-spherical molecules. Randomly distributed liquid molecules as in an ideal gas are uncorrelated with each other; second harmonic scattering (SHS) from these uncorrelated molecules is an incoherent sum of the SH response of individual molecules, which is also referred to as hyper-Rayleigh scattering (HRS). Molecular correlations can be induced by H-bonding, Coulombic interaction, dipole-dipole interaction, or other interactions. For water, the orientational correlation of a pair of water molecules can be changed by breaking an H-bond through a rotational or bending/libration motion, but not by stretching (see the inset in Fig. 1). The orientational correlation modifies the coherent contribution of the correlated molecules to the SH intensity. Note that femtosecond laser pulses that are shorter than the orientational relaxation time of water are necessary to measure snapshots of the static collective orientational correlations in water. As shown in Ref. 18, information on the two origins of SH responses can be obtained from the SH intensity measured in different polarization combinations of the incident and scattered light. The polarization combinations are indicated with a three letter code, from left to right, representing the polarization of the scattered SH light and the incident fundamental light with P(S) referring to linearly polarized light parallel (perpendicular) to the scattering plane. The PPP and PSS combinations mainly report on the collective orientational order, and the SH intensity consists of both coherent contributions from correlated molecules and incoherent contributions from uncorrelated molecules. The SSS and SPP combinations report on the electronic structure of individual molecules, and the SH intensity is dominated by the incoherent contributions from uncorrelated molecules. According to the selection rule of SHS, molecular correlations of spherical symmetry only contribute to the PPP and PSS SH intensities. 29 Figure 1(a) shows the relative change of the PPP fs-ESHS intensity (at scattering angle θ = 90 • ) for 21 different electrolyte solutions in the ionic strength range of 1 µM-0.1M, as reproduced from Ref. 18. The SH intensity was recorded at room temperature (296.15 K, 23 • C). The 21 electrolytes include a monovalent Cl cation series (NH 4 Cl, KCl, NaCl, LiCl, and HCl), an Na + anion series (NaH 2 PO 4 , NaOH, NaCl, NaBr, NaI, NaNO 3 , NaClO 4 , and NaSCN), divalent cationic-monovalent anionic pairs (MgCl 2 , CaCl 2 , SrCl 2 , and BaCl 2 ), monovalent cationic-divalent anionic pairs (Na 2 S 2 O 3 and Na 2 SO 4 ), and divalent cationic-anionic pairs (BeSO 4 , MgSO 4 , and CaSO 4 ). For all electrolyte solutions, the same trend in the SH intensity is observed: A continuous rise up to a plateau with a saturated increase by 30% ± 3%. The 50% half-saturation concentration (c * ) is at 55 ± 5 µM. This non-ion-specific, saturated increase in the PPP fs-ESHS intensity indicates that electrolytes induce orientational order in the H-bond network of bulk water that extends over tens of hydration shells. Figure 1(b) shows a comparison of the PPP fs-ESHS intensity of NaCl dissolved in H 2 O and NaCl, KCl, DCl, and NaOD dissolved in D 2 O. A similar trend in the SH intensity is observed for the D 2 O solution, but the intensity changes up to the plateau are only 9% ± 3% with c * = 310 ± 71 µM. Additional electrolyte solutions in D 2 O display the same behavior as NaCl solutions [ Fig. 1(b)]. Provided that the differences in the linear dielectric properties of H 2 O and D 2 O are minimal, such as permanent dipole moment and dielectric constant, this significant difference in the SH response is a manifestation of stronger H-bonding in D 2 O than in H 2 O, which arises from nuclear quantum effects. [19][20][21][22]25,26,[30][31][32][33][34][35][36] It also indicates the importance of collective H-bonding. 18 Since these studies clearly show the importance of H-bonding and since temperature can modify the probability that H-bonds are formed and broken, 37 it is of interest to investigate the above behavior for different temperatures.
Here we investigate the temperature dependence of the fs-ESHS response of bulk H 2 O and D 2 O as well as NaCl solutions in light and heavy water. We find that bulk light and heavy water both display a reduction in their intermolecular coherent responses as the temperature is increased from 263.35 K (9.8 • C) to 321.45 K (48.3 • C). However, both liquids have distinct temperature behaviors with a larger intermolecular coherent contribution in D 2 O arising from stronger spherically symmetric orientational correlations in the H-bond network. The relative temperature dependent change is bigger for H 2 O than for D 2 O. For electrolyte solutions, the intensity vs temperature curves, such as recorded in Fig. 1(b), shift to higher ionic strengths due to a decreased screening of the electrostatic field at higher temperatures as explained by the Debye-Hückel theory. However, the half-saturation concentrations change more than the model prediction, with D 2 O deviating by a factor of 6.3 and H 2 O by a factor of 1.9. These measurements display the complexity of liquid water and provide a challenge for future models of water and electrolyte solutions.

Femtosecond elastic second harmonic scattering
Laser pulses (190 fs) centered at 1028 nm with a 200 kHz repetition rate were used as the light source for femtosecond elastic second harmonic scattering measurements. The polarization of the input pulses was controlled by a Glan-Taylor polarizer (GT10-B, Thorlabs) and a zero-order half wave plate (WPH05M-1030). The incident laser pulses, filtered by a long pass filter (FEL0750, Thorlabs), with a pulse energy of 0.3 µJ (incident laser power P = 60 mW) were focused into a cylindrical glass sample cell (4.2 mm inner diameter, LS instruments) with a waist diameter of ∼35 µm and a Rayleigh length of 0.94 mm. The sample cell was placed in a customized temperature controller (Quantum Northwest) that provided a precise control of the temperature of the sample. The temperature can be tuned in the range of 253.15 K (20 • C)-423.15 K (150 • C) with a precision of ±0.1 K. The scattered SH light was collected with a plano-convex lens (f = 5 cm) and then filtered by a bandpass filter (ET525/50, Chroma). A Glan-Taylor polarizer (GT10-A, Thorlabs) was used for the polarization analysis of the SH light. In the end, the SH light was focused into a gated photomultiplier tube (PMT) (H7421-40, Hamamatsu). The detection angle was set to 90 • with an acceptance angle of 11.4 • . For every data point, 3 to 11 measurements were performed and the averaged results are shown. In every measurement, the SH intensity was acquired with 50 × 1 s acquisition time (i.e., using 50 × 2·10 5 pulses in total) and a gate width of 10 ns. The SH intensity of pure H 2 O or D 2 O (Figs. 1 and 3) was measured between every two samples and was used as a reference. The reproducibility of the fs-ESHS measurements is 1%-3%. Samples were stored and measured in sealed glass sample cells. The obtained relative intensities as a function of ionic strength were fit with a Langmuir-type equation derived from Debye-Hückel theory, whose functional form is given by 1+(c / c * ) , with c as the ion concentration and A as a fitting parameter. The 50% onset concentration c * was extracted from the fitting. More details of fs-ESHS can be found in Ref. 38  Fig. 1(a) and to possibly assign it to changes in the orientational correlations in the H-bond network, we need to consider the sources of intensity change.

Temperature dependence of molecular correlations in water
The SH intensity has two contributions: an intermolecular coherent and a molecular incoherent one. The temperature dependence of the incoherent response arises from the temperature dependence of the molecular hyperpolarizability and the dielectric constant that is connected to the local field factors. Thus, with increasing temperature, we may expect that the dielectric constant, the second-order hyperpolarizability, or the square of the second-order hyperpolarizability decreases or that the orientational correlations in the liquid become less probable. The dielectric constants of H 2 O and D 2 O differ by 0.5% at different temperatures, 39 and the hyperpolarizability is a reflection of the electronic structure of individual molecules and has been reported to be almost the same. 40 Therefore, the difference in Fig. 1(a) stems from different degrees of intermolecular correlations between H 2 O and D 2 O. D 2 O has more intermolecular correlations (or correlations that live longer), but changes less with temperature. This is consistent with the notion that D 2 O has stronger H-bonds than H 2 O.
In order to separate the molecular incoherent and intermolecular coherent contributions, we compare the response for PPP polarized light (that contains coherent and incoherent contributions for spherically symmetric systems) to the response for SPP polarized light (which contains only incoherent contributions). The incoherent contribution to the SH PPP and SPP intensity can be computed considering random uncorrelated point dipoles (with p (2) i = β (2) : E(ω) E(ω) the molecular second-order dipole moment for molecule i) 41,42 where c is the velocity of light in vacuum, n is the refractive index of air, k 0 represents the magnitude of the wavevector of the emitted SH light from the induced molecular dipole, E x is the component of the electric field along the horizontal x-axis of the lab frame, orthogonal to the forward direction of incident laser beam (z-axis), θ is the scattering angle between the observation direction and the forward direction, N m is the number density of the molecule, V is the volume in which molecules contribute to the SHS intensity and typically corresponds to the focal volume of the incident laser beam, R is the distance between the scattering point and the observation point, and β (2) ijk corresponds to the values of molecular second-order hyperpolarizability tensor for different polarization combinations (ijk). Similar expressions can be found for the other polarization combinations. 42 For a molecular distribution with spherical symmetry, the SPP polarization combination contains only molecular contributions. The PPP polarized light generally contains both molecular and intermolecular contributions. Since at the used scattering angle (θ = 90 • ) Eq. (1) is identical for both polarization combinations, dividing these at θ = 90 • will thus remove the molecular incoherent contribution. The caveat here is that this works only for nonlinear optical scattering theory using spherical symmetry. This spherical symmetry refers to the spatial orientational distribution of H-bonded water molecules. A water molecule with 2 accepting and 2 donating H-bonds is spherically symmetric. A water molecule with a single accepting H bond is not spherically symmetric just like other H-bond-missing configurations. To go beyond this approximation requires extensive path integral molecular dynamics. 40 Figure 2(b) shows the ratio of the measured PPP and SPP SH intensity at different temperatures. The temperature dependence of this intensity ratio is remarkably different compared to the PPP intensity shown in Fig. 2 This different temperature dependence of the intensity ratio indicates a more complicated picture of the bulk water structure evolving with temperature than a mere change in the H-bond number. Breaking down the SH intensities, the SPP intensity consists primarily of an incoherent contribution from uncorrelated water molecules. We note that it could also contain a coherent contribution from orientationally correlated molecules that are distributed non-spherically as remarked above (as this is an approximation in the nonlinear optical theory 29 ). By contrast, the PPP intensity contains a coherent contribution from spherically correlated molecules (and of non-spherically correlated molecules if the spatial distribution is non-spherical). We thus have two possible ratios (R = I PPP /I SPP ) under the measurement condition, where incoh in the subscripts denotes an incoherent contribution, coh, nonsph and coh, sph denote coherent contributions from spherically correlated molecules and from non-spherically correlated molecules, respectively. The excess amount of PPP SH intensity over SPP SH intensity [ Fig. 2(b)] can then be attributed to two possible factors. From Eq. (2) we derive that the intermolecular orientational correlations should increase with increasing temperature. This, however, is counter-intuitive because there are less H-bonds and intermolecular correlations with increasing temperature [ Fig. 2(a)]. Then, using Eq. (3) an increase in the ratio R with increasing temperature can arise from a reduction in the amount of non-spherical orientational correlations combined with a relative increase in the amount of spherical orientational correlations. This makes sense because the H-bond network of water loses structure when temperature is increased. [43][44][45][46] The larger ratio R of D 2 O thus indicates relatively more spherical intermolecular correlations than H 2 O.

Temperature dependence of ion-water interactions
In a second set of experiments, we measured the temperature dependence of the water-water orientational correlations induced by ions in NaCl solution in H 2 O and D 2 O. Figure 3 shows the normalized fs-ESHS intensity as a function of ionic strength in H 2 O [ Fig. 3 Fig. 3(b)], comparing the effects of Na + and Cl ions at the temperatures of maximum density for H 2 O (277.15 K, 4 • C) and D 2 O (284.15 K, 11 • C) to a fixed high temperature (321.15 K, 87 • C). As can be seen, there is a significant temperature effect in both solvents: Increasing the solution temperature shifts the increase in the SH intensity to higher ionic strength. The halfsaturation points (c * ) in these curves change from c * (4 • C) = 44 ± 5 µM → c * (87 • C) = 69 ± 5 µM (∆c * = 25 ± 7 µM) for H 2 O and c * (4 • C) = 300 ± 51 µM → c * (87 • C) = 805 ± 130 µM (∆c * = 505 ± 140 µM) for D 2 O. Note that the c * values at room temperature are at 55 ± 5 µM and 310 ± 71 µM. The shift in the SH intensity curves indicate a weakened influence of the electrostatic fields of the ions on the water structure induced by increasing the temperature. The temperature dependent change in concentration can be qualitatively explained by the Debye-Hückel model. 18 With increasing tem-

(a)] and D 2 O [
is the Boltzmann constant, T the temperature in the unit K, N A the Avogadro number, and I the ionic strength) increases, leading to a reduction of the influence of ions on the orientational correlations in the water. Additionally the ion induced dipole-dipole correlations have a temperature dependence of the form 27 C (r) ∼ 1/T 2 exp (−κ D r) /r. Here, an increase of temperature results in a reduction of the dipole-dipole correlations. Together, these effects result in 18 c * = q 2 ε 0 εk B T 2(Ze) 2 , where q is the magnitude of the scattering vector and Z is the valence of the considered ion.
As such, a shifted intensity curve is expected. However, the model predicts ∆c * = ∼13 µM and ∆c * = ∼80 µM shifts in the ionic strength for H 2 O and D 2 O, respectively, using the measured values for c * at the maximum density temperatures as input values. This qualitative agreement shows that the Debye-Hückel model can partially capture the temperature dependence of the ions' influence on the water structure. Although the direction of the shift in c * is in the correct direction, for both H 2 O and D 2 O, the degree of change is significantly different. The screening capacity of the liquid is larger than that predicted for both liquids and relatively more for D 2 O than for H 2 O. This is surprising and somewhat counter-intuitive as one might expect that with higher temperatures H-bonds become weaker and so the influence of the electrostatic field on the orientational correlations would be reduced. That we observe the opposite suggests that there are more complicated interactions than the ones considered here. Indeed, the Debye-Hückel theory is not sophisticated enough to achieve a quantitative agreement. The underlying reason for that is that the model neglects any form of H-bonds. Also, nuclear quantum effects are not captured by the Debye-Hückel theory. 18 As recent computational studies show the importance of instantaneous fluctuations in the local H bond structure and since minute millidegree changes in the orientational order of water can already significantly impact the SH intensity, 18 these measurements provide a future opportunity for water models that should include both quantum effects and detailed non-local structures. 51

CONCLUSIONS
In summary, using fs-ESHS we observed clear temperature dependent changes in the molecular structure of the water H-bond network. The ensemble orientational order of water molecules decreases as the temperature increases. H 2 O and D 2 O show differences in their temperature dependence of the molecular structure of the H-bond network. With stronger H-bonding, D 2 O has a more structured H-bond network in which molecular correlations exhibit relatively more spherical symmetry compared to H 2 O. An increase in the temperature leads to a transformation of molecular correlations in pure water from correlations with non-spherical symmetry (e.g., single donor species) toward correlations with spherical symmetry (e.g., double donor species). For electrolyte solutions, we measured temperature dependent SH intensity versus concentration curves that show a decrease in the influence of the combined electrostatic field in the solution on the water-water correlations. The trend qualitatively follows the prediction from the Debye-Hückel theory, but the effects are larger than predicted: ∆c * = 25 µM for H 2 O and ∆c * = 505 µM for D 2 O were measured while the Debye-Hückel theory predicts ∆c * = 13 µM for H 2 O and ∆c * = 80 µM for D 2 O. Since the SH response can be computed from molecular dynamics simulations, the present observations provide a unique opportunity to refine quantum mechanical models of water.