Structures and Energetics of Erbium Chloride Complexes in Aqueous Solution

Bridging molecular-scale interfacial science with continuum-scale models

Solid-water interfaces are crucial for clean water, conventional and renewable energy, and effective nuclear waste management. However, reflecting the complexity of reactive interfaces in continuum-scale models is a challenge, leading to oversimplified representations that often fail to predict real-world behavior. This is because these models use fixed parameters derived by averaging across a wide physicochemical range observed at the molecular scale. Recent studies have revealed the stochastic nature of molecular-level surface sites that define a variety of reaction mechanisms, rates, and products even across a single surface. To bridge the molecular knowledge and predictive continuum-scale models, we propose to represent surface properties with probability distributions rather than with discrete constant values derived by averaging across a heterogeneous surface. This conceptual shift in continuum-scale modeling requires exponentially rising computational power. By incorporating our molecular-scale understanding of solid-water interfaces into continuum-scale models we can pave the way for next generation critical technologies and novel environmental solutions.

Direct Observation of Contact Ion-Pair Formation in La3+ Methanol Solution

An approach combining molecular dynamics (MD) simulations and X-ray absorption spectroscopy (XAS) has been used to carry out a comparative study about the solvation properties of dilute La(NO₃)₃ solutions in water and methanol, with the aim of elucidating the still elusive coordination of the La³⁺ ion in the latter medium. The comparison between these two systems enlightened a different behavior of the nitrate counterions in the two environments: while in water the La(NO₃)₃ salt is fully dissociated and the La³⁺ ion is coordinated by water molecules only, the nitrate anions are able to enter the metal first solvation shell to form inner-sphere complexes in methanol solution. The speciation of the formed complexes showed that the 10-fold coordination is preferential in methanol solution, where the nitrate anions coordinate the La³⁺ cations in a monodentate fashion and the methanol molecules complete the solvation shell to form an overall bicapped square antiprism geometry. This is at variance with the aqueous solution where a more balanced situation is observed between the 9- and 10-fold coordination. An experimental confirmation of the MD results was obtained by La K-edge XAS measurements carried out on 0.1 M La(NO₃)₃ solutions in the two solvents, showing the distinct presence of the nitrate counterions in the La³⁺ ion first solvation sphere of the methanol solution. The analysis of the extended X-ray absorption fine structure (EXAFS) part of the absorption spectrum collected on the methanol solution was carried out starting from the MD results and confirmed the structural arrangement observed by the simulations.

The formation of contact ion pairs between the La³⁺ ions and the nitrate anions has been demonstrated in diluted methanol solution using a combined approach using Molecular Dynamics simulations and X-ray absorption spectroscpy.

Hydration of ferric chloride and nitrate in aqueous solutions: water-mediated ion pairing revealed by Raman spectroscopy

Iron is the most abundant transition metal in the earth's crust and is important for the proper functioning of many technological and natural processes. Despite the importance, a complete microscopic understanding of the hydration of ferric ions and water mediated ion pairing has not been realized. Hydrated Fe(iii) is difficult to study due to the process of complexation to the anion and hydrolysis of the hydrating water molecules leading to a heterogeneous solution with diverse speciation. Here, ferric chloride and nitrate aqueous solutions are studied using polarized Raman spectroscopy as a function of concentration and referenced to their respective sodium salt or mineral acid. Perturbed water spectra (PWS) were generated using multivariate curve resolution-alternating least squares (MCR-ALS) to show the residual spectral response uniquely attributable to the hydration of ferric speciation. The hydrogen bonding network associated with the hydrating water molecules in ferric chloride solutions are found to be more similar to hydrochloric acid solutions, whereas in ferric nitrate solutions, the network behaves more similar to sodium nitrate, despite increased acidity. Thus, in the FeNO3 and FeCl3 solutions, ion pairing and coordination, respectively, are significantly influencing the hydration spectra signature. These results further reveal concentration dependent changes to the hydrogen bonding network, hydrating water symmetry, and changes to the relative abundance of solvent shared ion pairs that are governed primarily by the ferric salt identity.

Ion Association in Lanthanide Chloride Solutions

A better understanding of the solution chemistry of the lanthanide (Ln) salts in water would have wide ranging implications in materials processing, waste management, element tracing, medicine and many more fields. This is particularly true for minerals processing, given governmental concerns about lanthanide security of supply and the drive to identify environmentally sustainable processing routes. Despite much effort, even in simple systems, the mechanisms and thermodynamics of LnIII association with small anions remain unclear. In the present study, molecular dynamics (MD), using a newly developed force field, provide new insights into LnCl3 (aq) solutions. The force field accurately reproduces the structure and dynamics of Nd3+ , Gd3+ and Er3+ in water when compared to calculations using density functional theory (DFT). Adaptive-bias MD simulations show that the mechanisms for ion pairing change from dissociative to associative exchange depending upon cation size. Thermodynamics of association reveal that whereas ion pairing is favourable, the equilibrium distribution of species at low concentration is dominated by weakly bound solvent-shared and solvent-separated ion pairs, rather than contact ion pairs, reconciling a number of contrasting observations of LnIII -Cl association in the literature. In addition, we show that the thermodynamic stabilities of a range of inner sphere and outer sphere LnCl x ( 3 - x ) + coordination complexes are comparable and that the kinetics of anion binding to cations may control solution speciation distributions beyond ion pairs. The techniques adopted in this work provide a framework with which to investigate more complex solution chemistries of cations in water.

Molecular Hydroxo-Bridged Dimers of Uranium(VI), Neptunium(VI), and Plutonium(VI): [Me4N]2[(AnO2)2(OH)2(NO3)4]

The synthesis of a series of molecular actinyl(VI), namely, uranium(VI), neptunium(VI), and plutonium(VI), hydroxo-bridged dimers is reported. These complexes were isolated from an aqueous nitrate solution by titration with tetramethylammonium hydroxide. The solid-state structures were determined using single-crystal X-ray diffraction, revealing molecular complexes with the formula [Me₄N]₂[(AnO₂)₂(μ₂-OH)₂(NO₃)₄], where An = U^VI, Np^VI, and Pu^VI. Spectroscopic data-UV-vis-near-IR absorption, IR, and Raman-were collected on the solutions and solid-state complexes where available and compared to those of the aqueous solutions from which the crystals formed. These data provide structural evidence for higher-order polynuclear complexes of actinyl(VI) complexes upon a pH increase in the aqueous solution, confirming earlier thermodynamic models.

Anion dependent ion pairing in concentrated ytterbium halide solutions

We have studied ion pairing of ytterbium halide solutions. THz spectra (30-400 cm^-1) of aqueous YbCl₃ and YbBr₃ solutions reveal fundamental differences in the hydration structures of YbCl₃ and YbBr₃ at high salt concentrations: While for YbBr₃ no indications for a changing local hydration environment of the ions were experimentally observed within the measured concentration range, the spectra of YbCl₃ pointed towards formation of weak contact ion pairs. The proposed anion specificity for ion pairing was confirmed by supplementary Raman measurements.

Linking Solution Structures and Energetics: Thorium Nitrate Complexes

Seeking predictive insights into how metal-ion speciation impacts solution chemistry as well as the composition and structure of solid-precipitates, thorium correlations, with both solvent and other solute ions, were quantitatively probed in a series of acidic, nitrate/perchlorate solutions held at constant ionic strength. Difference pair-distribution functions (dPDF), obtained from high-energy X-ray scattering (HEXS) data, provide unprecedented structural information on the number of Th ligating ions in solution and how they change with increasing nitrate concentration. A fit of the end member solution, Th (4 m perchloric acid and no nitrate), reveals a homoleptic Th aqua ion with 10 waters in its first coordination shell. Analyses of the acidic solutions containing nitrate reveal exclusively bidentate NO₃^- complexation with Th, consistent with published solid-state M^IV nitrate structures, where M^IV = Ce, Th, U, Np, Pu. Metrical fits of Th coordination as a function of nitrate concentration are used to calculate Th-NO₃ stability constants, information important to a molecular-scale description of reaction energetics. The coordination environments of Th in solution were compared with single-crystal structures obtained from their precipitates, Th(NO₃)₄(H₂O)₄ and Th(NO₃)₄(H₂O)₃·(H₂O)₂. Relative stabilities of the solid-state compounds, assessed based on the results of molecular quantum chemical calculations, reveal the importance of including an accurate description of complexed waters when predicting relative energetics of dissolved ions in aqueous solution.

Effects of multivalent hexacyanoferrates and their ion pairs on water molecule dynamics measured with terahertz spectroscopy

Broadband Fourier transform terahertz spectroscopy reveals that dynamical perturbations to the low-frequency dynamics of water molecules by multivalent hexacyanoferrate salts extend beyond the primary solvation shell.

Understanding the role of aqueous solution speciation and its application to the directed syntheses of complex oxidic Zr chlorides and sulfates

The lack of an in-depth understanding of solution-phase speciation and its relationship to solid-state phase formation is a grand challenge in synthesis science. It has severely limited the ability of inorganic chemists to predict or rationalize the formation of compounds from solutions. The need to investigate mechanisms that underlie self-assembly has motivated this study of aqueous Zr-sulfate chemistry as a model system, with the goal of understanding the structures of oligomeric clusters present in solution. We used high-energy X-ray scattering (HEXS) data to quantify Zr correlations in a series of solutions as a function of sulfate concentration. The pair distribution function (PDF) from the sulfate-free sample reveals that the average oligomeric Zr moiety is larger than the tetrameric building unit, [Zr4(OH)8(H2O)16](8+), generally understood to dominate its solution speciation. At sulfate concentrations greater than 1 m (molal), bidentate sulfate is observed, a coordination not seen in Zr(SO4)2·4H2O (2), which forms upon evaporation. Also seen in solution are correlations consistent with sulfate-bridged Zr dimers and the higher-order oligomers seen in 2. At intermediate sulfate concentrations there are correlations consistent with large Zr hydroxo-/oxo-bridged clusters. Crystals of [Zr18(OH)26O20(H2O)23.2(SO4)12.7]Cl0.6·nH2O (3) precipitate from these solutions. The Raman spectrum of 3 has a peak at 1017 cm(-1) that can be used as a signature for its presence in solution. Raman studies on deuterated solutions point to the important role of sulfate in the crystallization process. These solution results emphasize the presence of well-defined prenucleation correlations on length scales of <1 nm, often considered to be within the structurally amorphous regime.

Short-range and long-range order of phyllomanganate nanoparticles determined using high-energy X-ray scattering

High-energy X-ray scattering (HEXS) is used to explore the pH-dependent structure of randomly stacked manganese oxide nanosheets of nominal formula δ-MnO2. Data are simulated in real space by pair distribution function (PDF) analysis and in reciprocal space by both the Bragg-rod method and the Debye equation in order to maximize the information gained from the total scattering measurements. The essential new features of this triple-analysis approach are (1) the use of a two-dimensional supercell in PDF modeling to describe local distortions around Mn layer vacancies, (2) the implementation in Bragg-rod calculations of a lognormal crystal size distribution in the layer plane and an empirical function for the effect of strain, and (3) the incorporation into the model used with the Debye equation of an explicit elastic deformation of the two-dimensional nanocrystals. The PDF analysis reveals steady migration at acidic pH of the Mn atoms from layer to interlayer sites, either above or below the Mn layer vacancies, and important displacement of the remaining in-layer Mn atoms toward vacancies. The increased density of the vacancy–interlayer Mn pairs at low pH causes their mutual repulsion and results in short-range ordering. The layer microstructure, responsible for the long-range lateral disorder, is modeled with spherically and cylindrically bent crystallites having volume-averaged radii of 20–40 Å. Thebunit-cell parameter from the hexagonal layer has different values in PDF, Bragg-rod and Debye equation modeling, because of the use of different weighting contributions from long-range and short-range distances in each method. The PDFbparameter is in effect a measure of the average inlayer Mn...Mn distance and consistently deviates from the average structure value determined by the Bragg-rod method by 0.02 Å at low pH, as a result of the local relaxation induced by vacancies. The layer curvature increases the Bragg-rod value by 0.01–0.02 Å with the cylindrical model and as much as 0.04–0.05 Å with the spherical model. Therefore, in principle, the diffraction alone can unambiguously determine with good accuracy only a volume-averaged apparent layer dimension of the manganese oxide nanosheets. Thebparameter is model dependent and has no single straightforward interpretation, so comparison ofbbetween different samples only makes sense if done in the context of a single specified model.

A Review of High-Energy X-Ray Diffraction from Glasses and Liquids

This paper summarizes the scientific trends associated with the rapid development of the technique of high-energy X-ray diffraction over the past decade pertaining to the field of liquids, glasses, and amorphous materials. The measurement of high-quality X-ray structure factors out to large momentum transfers leads to high-resolution pair distribution functions which can be directly compared to theory or combined with data from other experimental techniques. The advantages of combining highly penetrating radiation with low angle scattering are outlined together with the data analysis procedure and formalism. Also included are advances in high-energy synchrotron beamline instrumentation, sample environment equipment, and an overview of the role of simulation and modeling for interpreting data from disordered materials. Several examples of recent trends in glass and liquid research are described. Finally, directions for future research are considered within the context of past and current developments in the field.

Hydration of lanthanide chloride salts: a quantum chemical and classical molecular dynamics simulation study

We present the results of a quantum chemical and classical molecular dynamics simulation study of some solutions containing chloride salts of La(3+), Gd(3+), and Er(3+) at various concentrations (from 0.05 to 5 M), with the purpose of understanding their structure and dynamics and analyzing how the coordination varies along the lanthanide series. In the La-Cl case, nine water molecules surround the central La(3+) cation in the first solvation shell, and chloride is present only in the second shell for all solutions but the most concentrated one (5 M). In the Gd(3+) case, the coordination number is ∼8.6 for the two lowest concentrations (0.05 and 0.1 M), and then it decreases rapidly. In the Er(3+) case, the coordination number is 7.4 for the two lowest concentrations (0.05 and 0.1 M), and then it decreases. The counterion Cl(-) is not present in the first solvation shell in the La(3+) case for most of the solutions, but it becomes progressively closer to the central cation in the Gd(3+) and Er(3+) cases, even at low concentrations.

Interactions of Al(III), La(III), Gd(III), and Lu(III) with the fused silica/water interface studied by second harmonic generation

The interactions of the trivalent metal cations Al(III), La(III), Gd(III), and Lu(III) with the silica/water interface were studied using the nonlinear optical technique of second harmonic generation (SHG). Specifically, the Eisenthal chi(3) technique was used to quantify the thermodynamics of trivalent ion adsorption to the bare fused silica surface. SHG adsorption isotherms were measured and fit with the triple layer surface complexation model to obtain adsorption free energies, binding constants, and interfacial charge densities. The adsorption free energy for Al(III) was found to be -37.2(5) kJ/mol, while the adsorption free energies for the three trivalent lanthanide cations ranged from -29.9(9) to -32.2(7) kJ/mol. Despite identical ionic charges, the metals under investigation displayed different affinities for the fused silica/water interface, and this finding is analyzed and interpreted in the context of size-dependent metal cation properties and metal ion speciation. The thermodynamic results from this work are valuable benchmarks for computer simulations of trivalent metal transport in the environment.

Activation of a PARACEST agent for MRI through selective outersphere interactions with phosphate diesters

Ln(S-THP)(3+) complexes are paramagnetic chemical exchange saturation transfer (PARACEST) agents for magnetic resonance imaging (MRI; S-THP = (1S,4S,7S,10S)-1,4,7,10-tetrakis(2-hydroxypropyl)-1,4,7,10-tetraazacyclododecane, Ln(III) = Ce(III), Eu(III), Yb(III)). CEST spectra at 11.7 T show that the PARACEST effect of these complexes is enhanced at neutral pH in buffered solutions containing 100 mM NaCl upon the addition of 1-2 equiv of diethylphosphate (DEP). CEST images of phantoms at 4.7 T confirm that DEP enhances the properties of Yb(S-THP)(3+) as a PARACEST MRI agent in buffered solutions at neutral pH and 100 mM NaCl. Studies using (1)H NMR, direct excitation Eu(III) luminescence spectroscopy, and UV-visible spectroscopy show that DEP is an outersphere ligand. Dissociation constants for [Ln(S-THP)(OH(2))](DEP) are 1.9 mM and 2.8 mM for Ln(III) = Yb(III) at pH 7.0 and Eu(III) at pH 7.4. Related ligands including phosphorothioic acid, O,O-diethylester, ethyl methylphosphonate, O-(4-nitrophenylphosphoryl)choline, and cyclic 3,5-adenosine monophosphate do not activate PARACEST. BNPP (bis(4-nitrophenyl phosphate) activates PARACEST of Ln(S-THP)(3+) (Ln(III) = Eu(III), Yb(III)), albeit less effectively than does DEP. These data show that binding through second coordination sphere interactions is selective for phosphate diesters with two terminal oxygens and two identical ester groups. A crystal structure of [Eu(S-THP)(OH(2))]((O(2)NPhO)(2)PO(2))(2)(CF(3)SO(3)) x 2 H(2)O x iPrOH has two outersphere BNPP anions that form hydrogen bonds to the alcohol groups of the macrocycle and the bound water ligand. This structure supports (1)H NMR spectroscopy studies showing that outersphere interactions of the phosphate diester with the alcohol protons modulate the rate of alcohol proton exchange to influence the PARACEST properties of the complex. Further, DEP interacts only with the nonionized form of the complex, Ln(S-THP)(OH(2))(3+) contributing to the pH dependence of the PARACEST effect.