2.1. X-ray Standing Wave Studies
Additional XRSW measurements were carried out to examine the adsorption behavior of
Ln3+ ions at solution surface, in particular possible enrichment of
Ln3+ ions at the air/liquid interface. In XRSW experiments the intensity of characteristic fluorescence exited by the incident X-ray beam is recorded as a function of the incident angle θ. The key idea of XRSW method is as follows: the angular dependence of fluorescence yield is highly sensitive to the position of atoms in the direction normal to the sample surface. Thus, XRSW measurements offer an opportunity to locate the atoms directly from the analysis of the corresponding fluorescence curve [
27].
Experimental XRSW data, obtained for studied
Ln3+ aqueous solutions, exhibited essentially similar behavior, as a typical example
Figure 1 shows angular dependence of Ce L3-fluorescence from Се(NO
3)
3 · 6H
2O salt solution. This section may be divided by subheadings. It should provide a concise and precise description of the experimental results, their interpretation, as well as the experimental conclusions that can be drawn.
Generally, two factors should be kept in mind, when analyzing XRSW data collected under TER. The first factor is the dramatic changes in electric field intensity above the reflecting surface as the incidence angle is scanned within the TER region. As a direct result of these changes large modulations arise in angular dependence of fluorescence yield from atoms located in the near-surface region; calculated angular dependence, presented in
Figure 1, is a perfect illustration of this characteristic feature. The calculations have been performed using the recursion formalism developed by Parratt [
28]. As can be seen in
Figure 1 angular dependence of fluorescence yield for near-surface distribution of atoms drastically increases from zero at θ = 0° and reaches the maximum value in the vicinity of the critical angle θ
C.
In marked contrast, angular dependence of fluorescence yield from the atoms, which are present in the bulk liquid subphase, exhibits absolutely different behavior. In this case the most important factor is the changes in the penetration depth of the electric field into the subphase. As the incident angle increases above the critical angle θ
С the penetration depth abruptly rises from several nanometers (for θ<θ
С) to several hundreds of microns (for θ > θ
С), that in turn results in the sharp increase in fluorescence signal. From
Figure 1 it is clear that Ce-fluorescence curve, recorded in our experiments, corresponds to this type of angular dependence. These observations evidence, that no enhancement of
Ln3+ concentration in the near-surface region occurs in the dilute aqueous solutions of
Ln(NO
3)
3 · xH
2O salts and the distribution of
Ln3+ ions can be considered as homogeneous.
2.2. X-ray Absorption Spectroscopy Studies
As a reference point in studying
Ln(NO
3)
3 · xH
2O aqueous solutions (c=20mM) we performed X-ray absorption spectroscopy measurements also for
Ln(NO
3)
3 · xH
2O microcrystalline powder samples (
Table S1). For appropriate fitting of the first intense peak in Fourier-transformed spectrum one should take into account two coordination shells around metal ion: oxygen atoms (in crystalline
Ln(NO
3)
3 · xH
2O: СN
Ln=4O
H2O+6O
NO3 except CNCe=11 - 5O
H2O+6O
NO3, CNYb(Lu)=9 - 3O
H2O+6O
NO3;
Ln bond is monodentate with water oxygen and bidentate with oxygen of nitrate groups), nitrogen atoms (in crystalline
Ln(NO
3)
3 ·
LH
2O coordination number is 3: CN
Ln=3N
NO3). Results of such fitting presented in
Table S1 and in
Figure S1.
These fits performed with fixed coordination numbers. Such fitting allows to estimate the possibilities of EXAFS technique to determine structural parameters. Observation results presented in
Table S1 demonstrate that only first coordination shell radius could be determined this appropriate precision.
Further, such analysis applied for experimental EXAFS data for solutions. But in case of solutions coordination numbers for oxygen and nitrogen coordination shells were treated like free parameters, but Debay-Waller parameter for nitrogen coordination shell was fixed to decrease the number of independent parameters during the fitting. Results also presented in
Table S1 and in
Figure S2. In case of solutions the only parameter which could be determined from EXAFS fitting with appropriate precision is an interatomic distance
Ln-O for the first coordination shell.
Figure 2 contains the dependence of the
Ln–O bond distances over the
Ln series.
Also plotted error bars associated with each data point are based on EXAFS fitting uncertainty. It’s clearly seen from presented
Ln-O distances that
Ln-O for
Ln = Tb, Dy, Ho, Er, Tm ions demonstrate significant contraction of oxygen shell in solutions. Such observation is also in agreement with measured XANES spectra, which are plotted in
Figure S3. Vertical lines show the position of the second maximum in XANES spectra, shifting these lines towards higher energies indicates shrinking of the first coordination shell.
Concluding the results obtained by examination of XANES spectra and fitting of the EXAFS ones, the contraction of the first coordination shell for metals Ln = Tb, Dy, Ho, Er, Tm is observed. Stability of the nitrate complexes couldn’t be determined by simple EXAFS fitting procedure.
It should be noted that nitrogen coordination number in aquatic solutions has been a subject of disagreement. Questions such as whether an eventual coordination change occurs within each series or whether inner or outer sphere complexation occurs for specific ligands have not been completely resolved. Thus, according to [
19] the «light» lanthanides in aqueous solutions are nine-coordinated, whereas the “heavy” ones are eight-coordinated. According to [
29] La-Nd has larger coordination numbers than Tb-Lu, and in the region between them (
Ln=Sm-Gd) there are transitional structures or a mixture of structures. The authors of [
30] believe that the coordination numbers change in the Nd-Dy region. It is interesting to note that the authors of the review [
31] do not exclude fractional coordination (between 8 and 9)
Ln=Gd-Ho in
Ln(NO
3) solutions. According to [
32] the exact coordination numbers of rare-earth ions in the nitrate complexes of solutions are unknown. The latter statement seems to be correct, since CN
Ln in aquatic solutions depends on the concentration: in highly concentrated solutions, coordination with fewer solvent molecules can be realized than in diluted ones, and even more so in significantly diluted ones, due to a shortage of water molecules.
Analysis of
Figure 2 and
Table S1 data allows distinguishing two
Ln(NO
3)
3 · xH
2O regions:
Ln=Ce-Eu (region 1) and
Ln=Eu-Lu (region 2). Moreover, it is not excluded that in region 1 of
Ln(NO
3)
3 · xH
2O solutions both bidentate and monodentate bonds of
Ln with nitrate groups are realized, and in region 2 only monodentate bonds occur.
According to Raman spectroscopy data [
33], in a wide range of Се(NO
3)
3 concentrations in aqueous solution, nitrate ions are bound to Ce both monodentately and bidentately, which is consistent with our data (
Table S1,
Figure 2). Complexation of Nd
3+ and Eu
3+ ions with nitrate group was studied by spectrophotometry and microcalorimetry [
20]. The authors are of the opinion that both inner-sphere and outer-sphere nitrate complexes of
Ln3+ ions exist in solutions. Based on thermodynamic and spectroscopic data, it is assumed that the weak complex of Nd
3+ with nitrate in solution forms an inner sphere, and the nature of complex formation increases with an increase in temperature. In addition, it is possible that nitrate binds Eu
3+ and possibly also Nd
3+ bidentately in aqueous solutions.
Bonal S. et al. [
19] used microcalorimetry to determine the stability constant, Gibbs energy, enthalpies and entropies to analyze the very weak complexation of
Ln=La-Lu nitrate anion in dilute aqueous solutions of
Ln(NO
3)
3 · xH
2O at room temperature (298.15 K) with analysis of changes in the thermodynamic properties of (
LnNO
3)
2+ across the lanthanide series. With a decrease in the ionic radius of
Ln3+, it is more difficult for the nitrate anion with a slightly larger size of the NO
3 group than a water molecule to “penetrate” into the inner sphere of the
Ln3+ cation. Therefore, for Tm-Lu the inner sphere consists exclusively of water molecules, which was confirmed by studying aqueous solutions of rare earth metal nitrates in a glassy state by Raman spectroscopy [
32]. Moreover, as the cation becomes smaller, the preference for monodentate nitrate binding increases, due to avoided repulsions in the first coordination sphere. These literature data do not contradict our analysis of EXAFS results (
Table S1,
Figure 2).
Dobler et al. [
16] showed by quantum chemical investigation that an increase in the number of water molecules in the first coordination sphere of
Ln3+ promotes monodentate coordination of nitrate. The change from bidentate to monodentate coordination is also observed before the salt dissociation, i.e. in an aqueous solution, nitrates on the way to dissociation pass from bidentate to monodentate coordination. The thermodynamic explanation for this process is interesting. In solution, the enthalpy (ΔH) and entropy (ΔS) energy components are antagonists. Bidentate binding may be promoted by entropy, which reduces the number of “frozen” water molecules coordinated with
Ln3+ ions. On the other hand, monodentate bonding may be preferable from an enthalpy point of view, since up to six hydrogen bonds can be formed with the oxygens of the three monodentate nitrates (O
NO3) instead of three with the bidentate nitrate groups. If we compare the course of interatomic distances (
Table S1,
Figure 2) according to EXAFS data with the curves ΔS of complexes (
LnNO
3)
2+ and ΔH from
Ln [
16], then we can detect similarities for
Ln=Ce-Eu with ΔS curve. This may be indirect confirmation of the presence of a certain amount of nitrate groups bidentately bound to
Ln=Ce-Eu.
The authors of [
17] used MD with explicit polarization and UV-visible spectroscopy to study solutions of Nd
3+ and Dy
3+ nitrates from a “highly diluted” solution to experimental saturation. It has been established that the bidentate mode is somewhat more stable for Nd
3+ than for Dy
3+; at the end of the lanthanide series, the ratio between the bidentate and monodentate conformations decreases, while only the monodentate mode is present in Lu.
Moreover, in solution, hydrated nitrate complexes can exhibit an equilibrium between several polyhydrate forms involving different types of nitrate binding modes [
34].
Analysis of the limited literature data on the study of dilute solutions of Ln(NO3)3 · xH2O salts confirms the different structural behavior of Ln3+ ions in them, and different results are observed depending on the calculation methods and experimental conditions. It should be noted that experimental concentrations are typically far from standard and cannot be correctly extrapolated to infinite dilution to obtain a thermodynamic equilibrium constant that is valid only for a particular medium and concentration range (ionic strength).
2.3. X-ray Diffraction
According to the structural analysis of
Ln(NO
3)
3 · xH
2O salts [
35,
36,
37,
38,
39,
40,
41,
42,
43,
44,
45], with an increase in the
Ln atomic number, the total content of water molecules decreases: x=6 for Ce - Sm, x=6 and 5 for Eu - Tb, x=5 for Dy-Yb (except for Tm with x=6 and 5), x=4 and 3 for Lu. In the inner sphere (in square brackets) of coordination compounds of the form [
Ln(O
2NO
3)
3(OH
2)
n] ·
(x-n)H
2O, the number of water molecules also decreases along the
Ln series (n=5 for Ce, n=4 for Pr-Yb, n=3 for Lu) connected monodentantly to
Ln. In the outer sphere (outside square brackets), the number of water molecules
(x-n) for all
Ln except Ce, and the content of (NO
3)
1- groups bidentantly coordinated with
Ln ions remains constant (3 NO
3).
It should be noted that in the
Ln(NO
3)
3 · 6H
2O structures, starting from Gd, the length of one
Ln-O
NO3 bond is greatly increased compared to others (bond asymmetry: Δ=202Å). The bond asymmetry increases in the structure with Tb (Δ=220Å) and is maximum (Δ=523Å) in the structure with Tm [
36]. This gave reason to the authors of [
36] to assume one monodentate and two bidentate bonds of NO
3 with Tm
3+ ions (there is no asymmetry of bonds in the Tm(NO
3)
3 · 5H
2O structure). This asymmetric bonding seems to be associated with a steric effect of the coordinating water molecules. The increasing asymmetry in the binding mode of one nitrate group is related to the decrease in the ionic radius of the
Ln.
An analysis of the interatomic distances given in the literature showed that for the same
Ln, for which x=5 and x=6 in the
Ln(NO
3)
3 · xH
2O composition, with an increase in molecules of crystallization water, the average
Ln-O
H2O distance decreases (Δ
Ln-O(H2O) ~0.032Å), and
Ln-O
NO3, on the contrary, increases, reaching a maximum value (Δ
Ln-O(NO3) =0.058Å) for Tm [
36]. In this case, the average
Ln-O interatomic distance in the first coordination sphere increases upon transition from x=5 to x=6. However, for Lu, in the transition from x=3 to x=4, the value of Δ
Ln-O(NO3) =0.004Å, which contributes to an increase in
Ln-O (R,Å) in the Lu(NO
3)
3 · 3H
2O structure compared to Lu(NO
3)
3 · 4H
2O Lu(NO
3)
3 · 3H
2O (
Table S3,
Figure 3).
The X-ray diffraction study of commercial
Ln(NO
3)
3 · xH
2O samples [
46] indicates single-phase
Ln(NO
3)
3 · 6H
2O with Ce (structure 1; x=6), Pr-Tb (structure 2; x=6), Tm (structure 3; x=5), two-phase with Dy-Er (structure 2 + structure 3), Yb (structure 3 + unknown structure 4) and uncertainty with Lu (unknown structure 5 or non-single-phase sample). Although the crystal structures of 1-3 differ from each other (the closest are structures 2 and 3), and the structures of 4 and 5 are unknown, their main structural fragments are the same.
Table S3 and
Figure 3a-c show the average interatomic distances
Ln-O (R, Å) without separation into
Ln-O
H2O and
Ln-O
NO3, calculated based on the
Ln coordination in the inner sphere, taking into account the quantitative analysis of non-single-phase
Ln(NO
3)
3 · xH
2O with Ln = Dy, Ho, Er (
Figure 3a, c, d; red dots). The structural parameters are consistent with the data given in [
16] for
Ln(NO
3)
3 · xH
2O with
Ln=Ce-Sm (x=6; CN=11 for Ce, CN=10 for Nd-Sm), but differ for Eu (x=6 according to our data, x=5 according to [
16]) and interatomic Lu-O distances (
Figure 3a, blue double dots), which we either calculated or took from structural data [
43].
X-ray diffraction patterns of
Ln(NO
3)
3 · xH
2O solutions (
Figure 4) are represented by two pronounced halos with interplanar distances in the intervals
d = ~7.16- ~6.34 Å (1st peak) and
d = 3.423-3.119 Å (2nd peak) and diffuse peak with a maximum at
d ~ 2.15Å (3rd peak).
A comparison of diffuse reflections of
Ln(NO
3)
3 · xH
2O solutions with
Ln=Ce, Sm, Gd, Yb (
Figure 4) with similar reflections of a dilute aqueous solution of lanthanum nitrate [
47] shows their similarity. The main diffuse peaks for water (
Figure 4, 1st and 2nd peaks) occur at 14.17 and 27.92°(2θ) (
d = 6.24 and 3.193 Å) and a very weak and diffuse peaks at ~39.85°(2θ) (
d ~ 2.26 Å) (
Figure 4, 3rd peak) are present in the diffraction patterns of all solutions.
If we take the intensities of the 1st (
I(1)
rel) and 2nd (
I(2)
rel) diffuse reflection for water as reference, then the intensity
I(1)
rel>
I(2)
rel for all concentrations of
Ln(NO
3)
3 · xH
2O solutions (
Figure 4). Moreover,
I(1)
rel>
I(2)
rel with c=134 Mm for
Ln=Ce, Sm (for Sm to a greater extent), for Yb with c=50 and 134 mM the values of
I(1)
rel and
I(2)
rel are almost identical, and for Gd the 1st and 2nd diffuse reflections with c = 50 mM are most pronounced (
Figure 4). This pattern of changes in the reflection intensity of
Ln(NO
3)
3 · xH
2O (
Ln = Ce, Sm, Gd, Yb) solutions corresponds to the characteristic structural features of commercial samples of these compositions.
When comparing diffraction patterns of aqueous solutions [
31,
47,
48] and
Ln(NO
3)
3 · xH
2O salts [
35,
37,
38,
43], the position of the 1st diffuse peak generally coincides with the region of the most intense Bragg reflections in the range of 2θ ~12÷~15°, in particular, 1
0 with
d~6.7 Å, caused by
Ln ions. These values are in the regions of increased electron density observed in the radial distribution functions of aqueous solutions.
It is necessary to pay attention to the background level in the diffraction patterns of solutions and water (
Figure 4), in particular, in the low-angle region adjacent to 2θ=2°, which logically increases with an increase in the
Ln(NO
3)
3 · xH
2O concentration from 0 (water) to c = 134 mM. This may be an indirect indication of the presence of clusters in solutions, the size of which increases with an increase in the solution concentration, since it is known that a system consisting of weakly bound atoms, ions or molecules forms a so-called cluster ion. It occupies an intermediate position between van der Waals molecules and molecular formations with a covalent chemical bond, and the components of the cluster retain their specific individuality.
2.4. FT-IR Spectroscopy Data
Figure 5a shows the IR spectra of solutions of commercial
Ln(NO
3)
3 · xH
2O salts (c=20 mM and 50 mM) with
Ln = Ce, Sm, Gd, Yb, and the bandwidth correspondence is presented in
Table S4.
The band shift ~1044 cm
-1 is associated with the bond covalency degree of nitrate ions with
Ln ones, which increases along the
Ln series. Based upon this shift in the 1050 cm
-1 region, the following series may be set up for increasing covalency of metal-nitrate bond: Ce<La<Gd=Y [
49] This should lead to a general tendency for interatomic distances and coordination numbers to decrease due to a decrease in the size of
Ln, which is confirmed by the data in
Figure 2. A very weak band at 1146 cm
-1 (
Figure 5b) may correspond to bending vibrations of the hydronium ion δ (H
3O
+), which can be present in a solution with pH<7 [
50]. The intensity of this band increases with an increase in
Ln concentration in
Ln(NO
3)
3 · xH
2O solutions with
Ln=Gd and Yb (
Figure 5b).
The splitting of the ~1400 cm
-1 band into two (~1345 and ~1400 cm
-1) (
Figure 5c) indicates an asymmetric vibration of the uncoordinated hydrated nitrate ion. The difference (Δν~125 cm
-1) between the bands ~1345 and ~1470 cm
-1 indicates partial coordination of the nitrate ions to
Ln [
51]. With an increase in the concentration of the
Ln(NO
3)
3 · xH
2O solutions (from 20 mM to 100 mM), the ~1470 cm
-1 band becomes more clearly defined (
Figure 5c). This means an increase in the number of bound nitrate ions and their gradual dominance over water molecules in competition to enter the inner coordination sphere. The band at 1470 cm
-1 is quite wide and, in accordance with [
12], apparently results from the superposition of bands of monodentate and bidentate coordinated nitrate ions [
52,
53].
Strong band at ~1635 см
-1 is the result of superposition of water scissoring bending and N=O stretching: the vibrational peak for N=O stretching of bidentate coordinated nitrate ions usually appears at 1630-1788 cm
-1 (
Figure 5d). The band at 1640 cm
–1 is related to the intermolecular bending oscillation ν
2 of the water molecule, which is due to the change in the angle between OH bands [
54]. From the above results it can be concluded that all bands in the 1625-1524 cm
-1 region characterize isolated N=O bonds (consistent with bidentate structure of nitrates) and the bands below 1520 cm
-1 are coupled vibrations, consistent with monodentate (but not only) nitrate structure [
55].
The weak band in the region of ~1750 cm-1 (Figures 5a, 5e) corresponds to the combination band of the symmetrical stretch and in-plane bending of the nitrate-ion. With an increase in the Ln atomic number, the position of this band shifts (Figures 5a, 5e) to the long-wavelength region (1745 cm-1 for Ce, 1749.6 cm-1 for Sm, 1752.3 cm-1 for Yb; it was not possible to detect it for Gd), due to a weakening of the O–N bond in the nitrate ligand owing to an increase in electron transfer to the Ln3+ metal ion with an increase in the charge density.
Bands above 1600 cm
-1 are usually attributed to N=O bond and suggest bidentate coordination. In contrast, when the highest frequency band is below ca. 1500 cm
-1, this is an indication of monodentate coordination [
55].
The position and intensity of the transmission bands corresponding to vibrations of nitrate ions in these spectra are very close for all samples, and the intensity of the bands increases with an increase in the concentration of solutions. It should be noted that from the IR spectra of Ln(NO3)3 · xH2O solutions with low-intensity, overlapping or even partially absent transmission bands, it is quite difficult to identify nitrate ions coordinated monodentately by Ln3+ ions. However, judging by the spectra, part (from ~10% to ~50%) of nitrate ions is coordinated by Ln3+ ions mono- or bidentately (mainly monodentately), and the remaining nitrate groups are uncoordinated.
A wide weak band at ~2114 cm
-1 (
Figure 5a) corresponds to the composite vibration of water molecules: bending vibration together with stretching, intermolecular, due to the rotation of the water molecule. The broad strong band ~3314 cm
-1 corresponds to the stretching of free, hydrogen-bonded and coordinated water molecules. The band is greatly broadened, which indicates the implementation of strong hydrogen bonding in
Ln(NO
3)
3 · xH
2O solutions (
Figure 5a).
The wide intensive band at 2700–3700 cm
–1 corresponds to three oscillations of the water molecule: asymmetric stretching oscillation (3490 cm
–1) symmetric stretching oscillation (3280 cm
–1, and the overtone of the bending oscillation (3250 cm
–1). The intensity of the maximum of the bending band, on the contrary, decreases with a reduction of the temperature in the same range [
54].
A wide intense band below ~800 cm
-1, responsible for stretching vibrations of water molecules, overlaps the bands of bending vibrations of nitrate ions lying in this spectrum region (
Figure 5a).
The IR spectra of
Ln(NO
3)
3 · xH
2O salts contain transmission bands both observed for solutions (marked with * in
Table S4 and in the text) and those belonging only to the salts. Thus, bidentately coordinated nitrate ions correspond to intense bands at ~1460 cm
-1* and 1280 cm
-1, medium-intensity bands at ~1660 cm
-1* (in solutions, this band is present for Gd and Yb) (
Figure 5d), and 1044 cm
-1*. The difference between the band values at 1280 and 1044 cm
-1 (Δ=236 cm-1) indicates bidentate chelating of nitrate ions [
55].
The ν(H2O) band splits into ~3560 cm-1, 3480 cm-1, ~3200 cm-1, with the exception of the salt with Ce, for which a broad band is observed at ~3500 cm-1 due to hydrolysis in air, characteristic of salts with Ln=La-Sm (mostly for La). The intense band at 1635 cm-1* refers to bending vibrations of water molecules.
Weak transmission bands at ~810 and 750 cm-1 are attributed to stretching and bending vibrations of (NO3); very weak bands at 675 and 610 cm-1 belong to Ln-О (νLn-O) vibrations.
Thus, the results of IR spectroscopy of commercial
Ln(NO
3)
3 · xH
2O salts with
Ln=Ce, Sm, Gd, Yb confirm the structural analysis data known from the literature on the implementation of only bidentately coordinated nitrate groups in them. As for
Ln(NO
3)
3 · xH
2O solutions, according to IR spectroscopy, in dilute solutions of 20 mM there are
Ln3+ ions coordinated by nitrate groups, but in very small quantities, bidentantly and monodentantly (mainly) linked with
Ln3+ ions. The number of monodentate nitrate groups increases from Ce to Yb: the ratio of the intensity of the band (shoulder) at 1750 cm
-1 to the intensity of the water bending vibration band at 1635 cm
-1 increases from Ce to Yb (
Figure 5e).
It should be noted that the peculiarity of the diffraction patterns of
Ln(NO
3)
3 · xH
2O solutions (
Ln=Ce, Sm, Gd, Yb) with c=50 and 134 mM (
Figure 4) is also preserved in the IR spectra of solutions with c=20 and 50 mM (
Figures 5b and 5d) and is consistent with the EXAFS data (
Table S1,
Figure 3).