On the Hydration of the Rare Earth Ions in Aqueous Solution

The totally symmetric stretching mode ν1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\nu_{1}$$\end{document} Ln–(OH2) of the first hydration shells of all the rare earth (RE) ions across the series from lanthanum to lutetium has been measured on dilute aqueous perchlorate solutions at room temperature. An S-shaped relationship has been found between the ν1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\nu_{1}$$\end{document} Ln–(OH2) peak positions and the Ln–(OH2) bond distances of the lanthanide(III) aqua ions. While the light rare earth ions form nona-hydrates, the heavy ones form octa-hydrates and the rare earth ions in the middle of the series show non integer hydration numbers between 9 and 8. A relationship between wavenumber positions ν1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\nu_{1}$$\end{document} Ln–(OH2) and the Ln–(OH2) bond distances of the RE hydrates has been given. Recent quantum mechanical calculations support the given interpretation.


Introduction
The trivalent lanthanides (Ln 3+ ) are part of the longest series of chemically similar metal ions marked by progressive filling of the 4f orbitals going from La 3+ to Lu 3+ . A subsequent decrease of the ionic radii along the series is observed due to the increase of the nuclear electric field. The decrease of the ionic radii of the rare earth ions with increasing atomic number, known as lanthanide contraction [1,2], was first coined by Goldschmidt [2]. A smooth plot of thermodynamic properties of lanthanide salts in aqueous solution would, therefore, be expected. This is especially true because the f-electrons do not take part in the bonding (negligible ligand field effect).
Various physical and chemical properties of the series of rare earth ions in aqueous solution have been studied over the years (apparent molal volumes, relative viscosities, 1 3 apparent molal heat capacities, heats of dilution and electrical conductance [3][4][5][6]). However, many thermodynamic and transport properties of rare earth ions in aqueous solution actually follow an irregular or S-shaped behavior when plotted against the atomic number or the ionic radius [3][4][5][6][7][8][9]. This behavior was seen as the result of the coordination number change of the first sphere in the middle of the lanthanide series following neodymium to terbium (see [7] and references therein). From X-ray and neutron diffraction measurements [10][11][12][13][14][15][16][17], it was deduced that the light rare earth ions exist as nona-hydrates, the heavy ones octhydrates and the ions in the middle of the series are mixtures of both (see ref [7].). Recent extended X-ray absorption fine structure (EXAFS) techniques confirmed these coordination structures [18,19]. Theoretical simulations [20][21][22][23] on the rare earth ion hydrates have also been carried out and supported the picture given by experimental structure techniques.
The tetrad effect [24] and the gadolinium break (half-filled f electron shell) [25] leads to small deviations from the expected smooth curvature. These effects are caused by the quantum mechanical inter electronic repulsion energy of the q electrons in the 4 f q electronic configuration [26]. The tetrad effect is much smaller than the effect caused by the change in coordination number across the whole lanthanide(III) series and it should be mentioned that the tetrad effect is not always observable [27].
In recent years, we characterized the hydrates of light rare earth ions and the members of the heavy ones in aqueous solution using Raman spectroscopy [28][29][30]. Raman spectroscopic data in conjunction with Ln-O bond distances are helpful to further shed light on the structure of the rare earth ions. If the totally symmetric stretching modes of these ions, 1 Ln-(OH 2 ) were to follow a rise with atomic number or decreasing ionic radius of these ions (lanthanide contraction) then the symmetric stretching modes 1 Ln-OH 2 would monotonically increase.
As an extension of our previous work on rare earth ions in solution, the Raman spectra were measured in the middle of the series of the lanthanide(III) hydrates in order to complete the data for the whole series. The Raman spectra of perchlorate solutions of all rare earth ions, Ln(ClO 4 ) 3 (aq), except promethium, are presented and discussed in this paper. Different concentrations of these rare earth ion solutions were measured, including fairly dilute solutions because the 1 Ln-(OH 2 ) modes showed slight changes with concentration [28][29][30]. Therefore it is crucial to consider only dilute solutions when comparing the 1 Ln-(OH 2 ) modes. Viewing our earlier results [28][29][30], it is apparent that the totally symmetric stretching modes, 1 Ln-(OH 2 ) for these rare earth ions, plot on two different straight lines with respect to the atomic number z, or the Ln-O bond distances. The light rare earth ions forming nona-hydrates in aqueous solution plot on an almost straight line and the one for the heavy rare earth octa-hydrates on another. These spectroscopic results, including the ones for the trivalent rare earth (RE) ions in the middle of the series, will be discussed in detail, taking into account recent results concerning thermodynamic and transport properties, structural results as well as quantum mechanical calculations on these lanthanide hydrates.
Raman spectroscopic measurements have been described on many occasions in detail and therefore only a brief discussion may be given. The scattering geometries, excitation conditions and Raman equipment were the same as reported in earlier publications and for details see [32,33]. Briefly, Raman spectra were measured in the macro chamber of the T 64,000 Raman spectrometer from Jobin Yvon in a 90° scattering geometry at 22 °C. A quartz cuvette was used (Hellma Analytics, Müllheim, Germany) with a 10 mm path length and a volume at 1000 µL. The spectra were excited with the 487.987 nm or the 514.532 nm line of an Ar + laser at a power level of 1000 mW at the sample. The Eu(ClO 4 ) 3 , Tb(ClO 4 ) 3 and Dy(ClO 4 ) 3 solutions were measured with the 514.532 nm Ar + line while Gd(ClO 4 ) 3 solutions with the one at 487.987 nm. After passing the spectrometer in subtractive mode, with gratings of 1800 grooves/mm, the scattered light was detected with a cooled CCD detector. The scattering geometries I VV = (X[ZZ]Y) and I VH = (X[ZX]Y) are defined as follows: the propagation (wave vector direction) of the exciting laser beam is in X direction and the propagation of the observed scattered light is in Y direction, the 90° geometry. The polarization (electrical field vector) of the laser beam is fixed in the Z direction (vertical) and the observed scattered light is polarized in he Z direction (vertical) for the I VV scattering geometry and polarized in X direction (horizontal) for the I VH scattering geometry. Thus, for the two scattering geometries it follows: and The symbols ′ and ′ denote the isotropic and the anisotropic invariant of the Raman polarizability tensor [31]. The isotropic spectrum, I iso was constructed according to Eq. 3: The polarization degree of the Raman bands, ρ (ρ = I VH /I VV ) was determined using an analyzer and adjusted, if necessary, before each measuring cycle using CCl 4 [32,33].
In order to characterize the spectral features in the low wavenumber region, the Raman spectra in I-format were reduced and the R-spectra obtained in order to obtain spectra defined as R(̃) which are independent of the excitation wavenumber, ̃L ; the measured Stokes intensity should be corrected for the scattering factor (̃L −̃) 3 in the case of applying counting methods [32,33]. The spectra were further corrected for the Bose-Einstein temperature factor, B = [1 − exp(−h̃c∕kT)] and the frequency factor, ̃ , to give the socalled reduced spectrum, R(̃) and again details were given in earlier publications [32,33]. (1) The isotropic spectrum in R-format follows from the corrected R VV and R VH spectra in the same manner as those in I-format (see Eq. 3). In the low wavenumber region, the I(̃) and R(̃) spectra are quite different and allowed the observation of an almost flat baseline [33,34].

Results
To complete the characterization of the RE(III) ions for the whole lanthanide series (except promethium), it is necessary to consider dilute aqueous perchlorate solutions of europium, gadolinium, terbium and dysprosium. The peak positions of the entire series of the modes 1 Ln-(OH) 2 for the rare earth ions, Ln 3+ (aq) including the ones in the middle of the series and those measured in previous studies [28][29][30] are given in Table 1. In addition to the 1 Ln-OH 2 values, the Ln-O bond distances were also included in Table 1 for all the RE ions. In Figs. 1 and 2, the representative Raman spectra of dilute Eu(ClO 4 ) 3 and Gd(ClO 4 ) 3 solutions are presented. The ClO − 4 (aq) modes have been given in recent publications (for details see refs [28,34].). Table 1 The symmetric stretching band, 1 Ln-O n of the hydrated [Ln(OH 2 ) n ] 3+ species (n = 9, 8 and 9/8) for dilute aqueous perchlorate solutions (C T < 0.35 mol·L -1 ) at 22 °C Given are the band positions, the force constants, the fwhh and the scattering intensities, S h of these bands. Furthermore, the Ln 3+ -O bond distances obtained by different methods (ND [17] and EXAFS [37] are presented as well as the rate constants for sulfate complex formation (k f ) [38] and the directly measured water exchange constants (k ex ) [39,40]   is depolarized and appears as a broad band at 1112 cm -1 . The inset gives the isotropic scattering profile in the terahertz region in more detail. The isotropic band at 374 cm -1 is assigned to the Gd-OH 2 symmetric stretching mode It is important to note that the 1 Ln-OH 2 mode is also called the symmetric breathing mode 1 since the oxygen atoms of the first hydration shell vibrate along the Ln-O bonds, stretching them in-phase. The molecule ion retains its symmetry and, therefore, the resulting dipole moment is zero. This normal mode depends on only one geometrical parameter: the radius r, the Ln-(OH 2 ) bond distance. The breathing mode 1 Ln-OH 2 of the rare earth ion hydrates was found to be almost completely polarized (ρ ~ 0) in solution, hence isotropic in character.
The peak positions of the totally symmetric stretching mode, ν 1 Ln-(OH) 2 , for the rare earth ions in the middle of the series show a slight concentration effect which has been observed also for the light and heavy rare earth ions even though perchlorate is known as a weak complex forming anion [28][29][30]35]. The concentration dependence of the 1 Ln-(OH 2 ) mode of the RE ions in the middle of the series may be illustrated by the spectra of Gd(ClO 4 ) 3 (aq) and Dy(ClO 4 ) 3 (aq). This phenomenon was described for a variety of lanthanide solutions in our previous contributions [28][29][30]35]. These earlier Raman spectroscopic contributions characterizing rare earth ions in aqueous perchlorate solutions showed that the solute concentration must be sufficiently small (~ 0.3 mol·L -1 , which translates to a solute to water ration at ca. 1:182) in order to avoid ion pair formation. The concentration effect expresses itself in concentrated solutions as a slight downshift of the peak position on the 1 Ln-(OH 2 ) mode and the subsequent broadening of its band width (full width at half height, fwhh). In concentrated solutions of Gd(ClO 4 ) 3 (aq) larger than 2 mol·L -1 , for instance, evidence for ion pair formation is found, although perchlorate is a weak complex forming anion. In these concentrated solutions there is simply not enough water to form complete hydration spheres of all ions (see [41] and discussions therein). The Raman spectrum of a 2.418 mol·L -1 Gd(ClO 4 ) 3 solution given in Fig. S1 shows the isotropic 1 Gd-OH 2 mode which is slightly shifted to a lower peak position at 371 cm -1 compared with the one in a dilute solution where it appears at 374 cm -1 (Table 1). For the Dy(ClO 4 ) 3 (aq) solution at 2.554 mol·L -1 (Fig. S2), a similar situation was observed, namely a shift of the ν 1 Dy-(OH 2 ) mode to slightly lower wavenumbers with a peak position at 381 cm -1 compared with the one for a 0.364 mol·L -1 Dy(ClO 4 ) 3 (aq) with a 1 Dy-(OH 2 ) mode at 384.5 cm -1 (Table 1). These peak shifts are also accompanied with a broadening of the band by ca. 15 cm -1 . Furthermore, for an undisturbed observation of the 1 Ln-(OH 2 ) mode, the nature of the counter ion is of importance. It may either act as weak complex forming anions such as perchlorate or trifluoromethansulfonate (triflate) or the ones with the tendency to form complexes such as chloride, nitrate and sulfate. In aqueous La 2 (SO 4 ) 3 solutions, which are quite dilute because of the sparingly soluble solute, the sulfato complex formation could, nevertheless, be detected in solutions as dilute as 10 -3 mol·L -1 . The breathing mode of the La 3+ of the sulfato-complex, [La(OH 2 ) n SO 4 ] + , is down shifted by 31 cm -1 compared with the breathing mode of the non-complexed, fully hydrated [La(OH 2 ) 9 ] 3+ species with a 1 value at 343 cm -1 [42].
When following the breathing modes 1 Ln-OH 2 for the RE ions of the whole series as a function of the atomic number z, an S-shaped curve results (Fig. S3). However, when giving the 1 Ln-OH 2 breathing mode dependence as a function of their corresponding Ln-O bond distances an inverse shaped function results. This is simply due to the fact that with increasing atomic number of the RE ions the Ln-O bond distances of the RE ion hydrates decrease. The 1 Ln-OH 2 breathing modes for the RE ions of the whole series as a function 1 3 of the Ln-O bond distance, r is given in Fig. 3 and the RE ions are given on top of the abscissa of the graph. The rare earth ions Ln 3+ (aq) La 3+ , Ce 3+ , Pr 3+ and Nd 3+ form nonahydrates and show the Ln-OH 2 breathing modes at 343 cm -1 , 344 cm -1 , 347 cm -1 and 351 cm -1 , respectively [30]. However, the 1 breathing mode for Sm 3+ (aq) appears at higher wavenumbers, namely at 362 cm -1 and it is thought to have a coordination number smaller than 9 namely 8.54 [15]. In addition, the Ln-OH 2 breathing modes of the heavy rare earth ions Ho 3+ , Er 3+ , Tm 3+ , Yb 3+ and Lu 3+ were also characterized by Raman spectroscopy and these ions are known to form octa-hydrates [15][16][17] and their Ln-(OH 2 ) breathing modes appeared at 387 cm -1 , 389 cm -1 , 391 cm -1 , 394 cm -1 , and 396 cm -1 , respectively [31]. It is of interest that Y 3+ (aq), having properties similar to the heavy rare earth ions dysprosium and holmium, is coordinated by eight water molecules in the first coordination sphere forming a square antiprismatic coordination polyhedron. The 1 Y-(OH 2 ) breathing mode for [Y(OH 2 ) 8 ] 3+ appears at 384 cm −1 [36] and with the ionic radius of Y 3+ (aq) at 2.36 Å [36] plots close to the one for 1 Dy-OH 2 (see Fig. 3).
The individual band parameters for all the 1 Ln-(OH 2 ) modes of the Ln 3+ (aq) ions are given in Table 1, specifically the scattering intensity of these bands, fwhh and the force constants of the modes. Additionally, the Ln-OH 2 bond distances (r) and kinetic parameters regarding the H 2 O exchange of the water molecules in the first sphere are presented. Our spectroscopic results will now be discussed in detail together with these parameters on Ln 3+ (aq) ions. High level quantum mechanical calculations in recent works and our own  [36]. The Ln-OH 2 bond distances including the one for Y-OH 2 were taken from ref [37] calculations will be taken into account to shed light upon the structure of the rare earth ion hydrates in solution.

Discussion of the Raman Spectroscopic Results
The peak position of the 1 Ln-(OH 2 ) modes of the rare earth ions given as a function of the Ln-O bond distances, r, 2) (Table 1) may be expressed by Eq. 4: and the following parameters were obtained: a = 380.5 cm -1 , b = 32.35 cm -1 , c = -8.97 Å -1 and r m = 2.435 Å. As previously mentioned, an inverse S-shaped dependence appears as a function of their Ln-OH 2 bond distances given in Fig. 3. The 1 Ln-(OH 2 ) breathing modes show a drastic change between samarium and terbium. However, then the 1 Ln-(OH 2 ) frequencies depend upon the atomic number z of the RE(III) ions following an S-shaped curve (Fig. S3). The light rare earth ions forming nona-hydrates in aqueous solution plot on an almost straight line beneath the one for the heavy rare earth octa-hydrates. The difference in wavenumbers between LRE ions and HRE ions is ca. 18 cm -1 and the ions in the middle of the series are displayed between these two lines. The first derivative of the Eq. 4 describing the relationship of 1 Ln-(OH 2 ) mode as a function of Ln-O bond distance, r, was taken and shows how much an error in bond distance with ± 0.02 Å affects the frequencies of the modes of the rare earth ion hydrates, 1 Ln-(OH 2 ). The derivative of Eq. 4 leads to Eq. 5: and the relation Δ 1 (r) ≈ d 1 (r) dr ⋅ Δr allows one to estimate the uncertainty Δ 1 (r) and to determine the frequency of the Ln-O mode when the uncertainty Δr in determining the Ln-O bond distances is given. For the light and heavy rare earth ions, an uncertainty of ± ~ 3 cm -1 follows and for the intermediate ions ± 6 cm -1 . The uncertainties in determining the peak positions of the 1 Ln-(OH 2 ) of the rare earth ions are ± 1 cm -1 and therefore are much smaller than the above calculated uncertainties. In other words, the symmetric breathing modes reflect the change in the first coordination numbers across the lanthanide series much more sensitively than methods determining the Ln-(OH 2 ) bond distances. The reason for that is the Raman frequencies depend not only on the bond distance but also on the reduced mass of the ligand.
The scattering intensity (relative integrated molar scattering intensity, S h ) of the 1 Ln-(OH 2 ) modes is very small and with their broad band profiles makes them difficult to detect by Raman spectroscopy ( Table 1). The relative integrated molar scattering intensity, S h , is defined by Eq. 6 as: 2 The published Ln-OH 2 bond distances vary several tenths of an Å which is outside the expected accuracy at ± 0.02 Ẳ [17]. Measurements via the ND first order technique [15][16][17] may serve as benchmark. The nature of the counter ions may be a problem especially in concentrated solutions because the frequently used anions such as Cl -, NO − 3 and SO 2− 4 form complexes with RE ions. Recent Raman spectroscopic measurements on a variety of RE ions in ClO − 4 solutions confirmed this point [28][29][30]. Furthermore, method specific problems may also contribute to the uncertainty of the Ln-OH 2 bond distances.
where A Ln 3+ and A ClO − 4 are integrated band intensities from R-spectra for the 1 Ln-OH 2 bands and the 1 ClO − 4 band respectively. C Ln 3+ and C ClO − 4 are the molar concentrations of the Ln 3+ ions and the perchlorate ions and C instr is a constant which corrects for the instrumental factors. The S h values for the lanthanide ion hydrates are small compared to softer more polarizable ions such as In 3+ [43]. It is known that the lanthanide ions are hard cations and therefore less polarizable ions [44].
The bandwidths of the isotropic 1 Ln-OH 2 bands are quite broad and are inhomogeneously broadened which is known as environmental broadening. In solution, the local environment varies significantly and slightly different spectrum will result and so the overall spectral shape is significantly broadened. Molecular dynamics (MD) results [45,46] demonstrated the flexibility of the lanthanide hydrate structures such as minor ultrafast geometric modifications within the hydrate or from pseudorotations, requiring substantially longer times in the scale at ca. 1 ps for establishment of an altered geometry [45]. The vibrational periods of the 1 Ln-OH 2 vibrations range from 0.097 ps for La 3+ to 0.082 ps for Lu 3+ [30,31,35] so that the modes vibrate several hundred thousand times before the hydrated Ln 3+ species experiences a water exchange. (The establishment of slightly altered geometry according to Rode et al. [45,46] is of the order of a ps and allows only a few vibrational periods to occur.) Water exchange reactions for the [Ln(OH 2 ) n ] 3+ species in aqueous solution (see Table 1) were reported from ultrasound absorption measurements [38] and more recently from NMR studies [39,40] (Table 1). From water exchange rates k 298 i of all the Ln 3+ (aq) ions (interchange rate of an inner-sphere water molecule by a SO 2− 4 ion [38]), it follows that the exchange rates k 298 i start at smaller values with La 3+ at 2.1 × 10 8 s -1 , reaching a maximum at Sm 3+ (7.4 × 10 8 s -1 ), remaining at high values in going to Tb 3+ (ions in the middle of the series) and then falling rapidly reaching Lu 3+ with a value at 6 × 10 7 s -1 . (More recent data on the rates of water exchange, k 298 ex measured directly by NMR [39,40] correlate closely with the k 298 i values for the RE ions reassuring the validity of the data in [38]). It is worthwhile to mention that recently a plethora of MD results (see for instance [45,46] and references therein) have been published giving structural data for the hydrated rare earth ions and also dynamic parameters. However, the water exchange rates put forward by these theoretical simulations are very short; actually too short to justify the aforementioned experimental values. It is clear that the water exchange kinetic in aqueous solution between first sphere water and second sphere water is a diffusion controlled process and therefore of limited rate (see Eigen [47]).
The force constant calculations for the 1 Ln-OH 2 breathing modes of the Ln(III) hydrate species (Table 1), applying a simple model, taking only one ligand into account, have been carried out according to Eq. 7: with c, the velocity of light, ̃i the wavenumber of the mode i, N the Avogadro constant and A L the molecular weight of the ligand, in our case, water. Such a simplification is justifiable because the totally symmetric character of the 1 Ln-OH 2 mode dictates that the central cation, Ln 3+ remains stationary and only the water molecules are involved in the breathing motion along the bonds without disturbing the symmetry and therefore the mode is totally polarized. The force constants, k Ln-O , calculated for the measured 1 Ln-OH 2 breathing modes, and the corresponding Ln 3+ -O bond distances [37] are given in Table 1.
The force constants increase in going from La 3+ to Lu 3+ form an S-shaped curve plotted against the atomic number of the rare earth ion which is shown in Fig. S3. The ions of the light rare earths (LRE) plot on as almost linear curve as do the heavy rare earth (HRE) ions which are displayed also on a linear curve offset to the former and with the ions in the middle of the series (Sm-Tb) appearing between the two lines (see Fig. 4). The LRE ions forming nona-hydrates, the heavy ones octa-hydrates and the RE ions in the middle show an overall hydration number between 9 and 8, with non integer coordination numbers. It should be mentioned that this S-shape behavior has been observed for a variety of thermodynamic properties of these ions and has been published for quite some time [7][8][9]. The assumed break of the properties between the light and heavy rare earths does not mean that the hydration number drops suddenly in the middle of the series but a dynamic hydration change from 9 to 8 may be envisioned. Such a dynamic structure has been presented recently by a quantum chemical (QC) study [22] taking into account two aspects, energy and wave function.  [37]. and the black filled squares refs [15][16][17]. from the ND first order technique The QC analysis by Zhang et al. [22] contributed greatly to the understanding of the structure and kinetics of the Ln 3+ (aq) ions across the whole lanthanide series. From their wave function analysis including the atoms in molecules approach by Bader (see ref. in [22]) it was shown that the nature of the Ln 3+ -oxygen bonds are electrostatic with very small covalent contributions. A specific behavior was found regarding the capping Ln−O bonds for the RE ions across the whole lanthanide series. The bond strength for the tri-capped water molecules is strong at the beginning of the series for the nonahydrates form lanthanum to samarium. But going from samarium to lutetium, the capping Ln−O bonds in the nona-aqua lanthanide hydrates become weaker and, unexpectedly, shorten. Finally, as the capping Ln−O bonds are getting easier to disrupt, heavier lanthanides will prefer a CN between 9 and 8. 3 This model suggests that the capping Ln−O bonds of moderate strength, which occur for intermediate lanthanides, are advantageous for the formation of an intermediate bi-capped trigonal prism (BTP) during water exchange and explains the kinetic of the water exchange rate reaching a maximum for the intermediate region of the rare earth ions but is lowest at the beginning and end of the series. This behavior of the capping Ln−O bonds explains not only the kinetics of the Ln 3+ hydrates across the series but also the structures and energetics of the Ln 3+ hydrates [22]. As a result, the hydrates of the lanthanide ions across the whole series form stable nona-hydrates at the beginning of the series (lanthanum to promethium), the heavy ones form octa-hydrates while the intermediate rare earth ions have a non integer coordination number between 9 and 8 according to the ratio of nona-and octaaqua hydrates. Such a peculiar behavior of capping Ln−O bonds put forward by Zhang and Dolg [22] determining the hydration behavior of Ln 3+ ions across the whole series is also portrayed by the course of the totally symmetric stretching mode, 1 Ln-OH 2 and their respective force constants across the series (Table 1 and Fig. S3) and reinforces these quantum mechanical model findings.
In the following, we will compare our experimental results on the 1 Ln-OH 2 breathing modes across the series (Figs. 4,5) with results from DFT calculations. However, the large number of electrons, partly filled 4f shells, and relativistic effects of Ln 3+ ions make quantum mechanical calculations difficult [21][22][23]49]. It should be pointed out that Zhang et al. [22] used different methods such as B3LYP, (SCS)-MP2, and CCSD(T) for the water cluster calculations in the inner shell and COSMO model for simulating the outer water shells. Reliable first principle values on hydration Gibbs energies were reported which correspond well with the experimental data and the corresponding coordination numbers 9, 8 and intermediate ones for light, heavy and intermediate rare earth ions. However, the geometry optimizations were restricted to gas phase clusters where the MP2 method was used but frequencies of the normal vibrations of the clusters were, unfortunately, not reported. The Ln-OH 2 bond distances of naked clusters are systematically overestimated in comparison with the experimental ones due to the neglect of the surrounding water molecules beyond the first hydration shell. As a result, vibrational frequencies are, therefore, underestimated. Due to many saddle points on the potential energy surface and convergence problems [49] it proved difficult to perform Ln-OH 2 cluster geometry optimizations within continuum models while taking higher hydration shells into account.
We  [50] with self-consistent field (SCF) procedure and starting without geometrical constraints. The presence of the tricapped trigonal structure (TTS) and square Given are the calculated bond distances for the nona-and octa-hydrates using the corresponding effective ionic radii for coordination 9 and 8, respectively [51], and adding the radius for the water molecules [52]. Upper (broken) line and red small data points correspond to the data of the nona-hydrates and below the small blue data points and the black line correspond to the octa-hydrates. The experimental Ln-OH 2 bond distances of the [Ln(OH 2 ) n ] 3+ (n = 9 or 8) hydrates are also given: filled triangle EXAFS [37], filled square XRD [10][11][12], inverted triangle ND [13] + ND [14], x ND [16]. Furthermore, the DFT Ln-OH 2 bond distances ( Table 2) are given for the nona-hydrate of La 3+ and its (hypothetical) octa-hydrate as well as for the octa-hydrate of Lu 3+ and its (hypothetical) nona-hydrate (Color figure online) antiprism structure (SAP) at the beginning and the end of the Ln 3+ series, respectively, was established and this is generally accepted. The B3LYP functional and Stuttgart/Dresden basis set SDD with effective core potential were used [51]. In Table 2 (Table 2). It should be pointed out that all normal vibrations of the naked clusters were without imaginary frequencies. With respect to the Ln-O bond distances and the breathing frequencies, two effects must be taken into account which we estimated and introduced in Table 2 and Fig. 4. Beginning with cerium and the elements that follow with increasing atomic number z, the 4f shell is successively filled with electrons shielding the central positive charge and diminishing the lanthanide ion radius by about 0.02 Å between neighboring ions for both cluster series, namely the nona-coordinated hydrates and the octa-coordinated ones (Fig. 5). This effect is known as the lanthanide contraction. Taking Shannon's effective radii for the coordination numbers 9 and 8 [52] and adding the mean water molecule radius at 1.34 Å [53] results in the Ln-OH 2 bond distances for the whole series (Fig. 5). The experimental bond distances appear at or between these two lines with the light rare earth ion hydrates at the line for the nona-hydrates, the heavy rare earths on the one for octa-hydrates and the intermediate ones between the two lines taking the considerable uncertainty of the bond distances into account. Secondly, a change from coordination number 9 to CN 8 results in a reduced Ln 3+ -O bond distance by ca. 0.04 Å for La 3+ and by 0.05 Å for Lu 3+ (Table 2). Decreasing the Ln-O bond distance, however, is accompanied by a frequency increase of the breathing vibration, 1 Ln-(OH 2 ). Changing the CN from 9 to 8 results in a frequency increase of 1 by ca. 14 cm -1 in going from [La(OH 2 ) 9 ] 3+ to the (hypothetical) [La(OH 2 ) 8 ] 3+ and 20 cm -1 frequency decrease in going from [Lu(OH 2 ) 8 ] 3+ to the (hypothetical) [Lu(OH 2 ) 9 ] 3+ . Along the two straight lines in Fig. 5 for the nona-hydrates and the octa-hydrates, this corresponds to a ca. 2.4 cm -1 frequency difference Δ for the LRE and a ca. 2.8 cm -1 Δ for the HRE between neighboring ions. It should be pointed out that this frequency difference between nona-and octa-hydrates is much larger than the uncertainty in determining 1 with an accuracy of ± 1 cm -1 [30,31,35] and, therefore, reflects the coordination change well.

Conclusion
The isotropic Raman bands of the 1 Ln-OH 2 modes of the rare earth ions in aqueous solution, [Ln(OH 2 ) n ] 3+ (n = 9,8 and 9/8) as a function of the Ln-OH 2 bond distance across the whole series show an S-shaped behavior. The light rare earth ions Ln 3+ (aq) (Ln = La 3+ to Nd 3+ ) plot on a straight line as do the heavy rare earth ions (Dy 3+ , Ho 3+ , Er 3+ , Tm 3+ , Yb 3+ and Lu 3+ ) which plot on a straight line above the one for the LRE ions. The middle rare earth ions (Sm 3+ , Eu 3+ , Gd 3+ and Tb 3+ ) plot inbetween. For the LRE ions which are known to be nona-hydrates, the oxygen atoms form a tricapped trigonal prismatic structure and the isotropic bands of these hydrates increase in an almost linear fashion with increasing atomic number, whilst the octa-hydrates of the HRE ions plot on a much higher straight line and the wavenumber difference between the nona-and octa-hydrates is ca. 18 cm -1 . Again, the 1 Ln-O frequencies of the Raman bands of the MRE ions plot in between these two straight lines. The S-shaped behavior found for the 1 Ln-OH 2 modes plotted against the atomic number for the hydrates of the rare earth ions across the whole series is in agreement with data on apparent molal volumes, relative viscosities, and molal heat capacities. No abrupt gadolinium break could be observed but a change of the coordination number from 9 to 8 in the middle of the series was detected. Recent QC results [22] explain the Raman data across the lanthanide series in a qualitative way but contradict the view of the specific asymmetry of the capping waters in the TTP structure of the nona-hydrate structure across the lanthanide series with water deficient configurations for the smallest ions (see [18] p. 3061-3062).