Exponential Sensitivity and Speciation of Al(III), Sc(III), Y(III), La(III), and Gd(III) at Fused Silica/Water Interfaces

Modeling separation of lanthanides via heterogeneous ligand binding

Individual lanthanide elements have physical/electronic/magnetic properties that make each useful for specific applications. Several of the lanthanides cations (Ln3+) naturally occur together in the same ores. They are notoriously difficult to separate from each other due to their chemical similarity. Predicting the Ln3+ differential binding energies (ΔΔE) or free energies (ΔΔG) at different binding sites, which are key figures of merit for separation applications, will help design of materials with lanthanide selectivity. We apply ab initio molecular dynamics (AIMD) simulations and density functional theory (DFT) to calculate ΔΔG for Ln3+ coordinated to ligands in water and embedded in metal-organic frameworks (MOFs), and ΔΔE for Ln3+ bonded to functionalized silica surfaces, thus circumventing the need for the computational costly absolute binding (free) energies ΔG and ΔE. Perturbative AIMD simulations of water-inundated simulation cells are applied to examine the selectivity of ligands towards adjacent Ln3+ in the periodic table. Static DFT calculations with a full Ln3+ first coordination shell, while less rigorous, show that all ligands examined with net negative charges are more selective towards the heavier lanthanides than a charge-neutral coordination shell made up of water molecules. Amine groups are predicted to be poor ligands for lanthanide-binding. We also address cooperative ion binding, i.e., using different ligands in concert to enhance lanthanide selectivity.

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.

Interplay of physically different properties leading to challenges in separating lanthanide cations - an ab initio molecular dynamics and experimental study

Lanthanide elements have well-documented similarities in their chemical behavior, which make the valuable trivalent lanthanide cations (Ln3+) particularly difficult to separate from each other in water. In this work, we apply ab initio molecular dynamics simulations to compare the free energies (ΔGads) associated with the adsorption of lanthanide cations to silica surfaces at a pH condition where SiO- groups are present. The predicted ΔGads for lutetium (Lu3+) and europium (Eu3+) are similar within statistical uncertainties; this is in qualitative agreement with our batch adsorption measurements on silica. This finding is remarkable because the two cations exhibit hydration free energies (ΔGhyd) that differ by >2 eV, different hydration numbers, and different hydrolysis behavior far from silica surfaces. We observe that the similarity in Lu3+ and Eu3+ ΔGads is the result of a delicate cancellation between the difference in Eu3+ and Lu3+ hydration (ΔGhyd), and their difference in binding energies to silica. We propose that disrupting this cancellation at the two end points, either for adsorbed or completely desorbed lanthanides (e.g., via nanoconfinment or mixed solvents), will lead to effective Ln3+ separation.

Restructuring of Hydration Shell Water due to Solvent-Shared Ion Pairing (SSIP): A Case Study of Aqueous MgCl2 and LaCl3 Solutions

Hydration of ions plays a crucial role in interionic interactions and associated processes in aqueous media, but selective probing of the hydration shell water is nontrivial. Here, we introduce Raman difference with simultaneous curve fitting (RD-SCF) analysis to extract the OH-stretch spectrum of hydration shell water, not only for the fully hydrated ions (Mg²⁺, La³⁺, and Cl^-) but also for the ion pairs. RD-SCF analyses of diluted MgCl₂ (0.18 M) and LaCl₃ (0.12 M) solutions relative to aqueous NaCl of equivalent Cl^- concentrations provide the OH-stretch spectra of water in the hydration shells of fully hydrated Mg²⁺ and La³⁺ cations relative to that of Na⁺. Integrated intensities of the hydration shell spectra of Mg²⁺ and La³⁺ ions increase linearly with the salt concentration (up to 2.0 M MgCl₂ and 1.3 M LaCl₃), which suggests no contact ion pair (CIP) formation in the MgCl₂ and LaCl₃ solutions. Nevertheless, the band shapes of the cation hydration shell spectra show a growing signature of Cl^--associated water with the rising salt concentration, which is a manifestation of the formation of a solvent-shared ion pair (SSIP). The OH-stretch spectrum of the shared/intervening water in the SSIP, retrieved by second-round RD-SCF analysis (2RD-SCF), shows that the average H-bonding of the shared water is weaker than that of the hydration water of the fully hydrated cation (Mg²⁺ or La³⁺) but stronger than that of the anion (Cl^-). The shared water displays an overall second-order dependence on the concentration of the interacting ions, unveiling 1:1 stoichiometry of the SSIP formed between Mg²⁺ and Cl^- as well as La³⁺ and Cl^-.

Interfacial electrostatics of poly(vinylamine hydrochloride), poly(diallyldimethylammonium chloride), poly-l-lysine, and poly-l-arginine interacting with lipid bilayers

Charge densities of cationic polymers adsorbed to lipid bilayers are estimated from SHG spectroscopy and QCM-D measurements.

Geochemical Insight from Nonlinear Optical Studies of Mineral-Water Interfaces

The physics and chemistry of mineral-water interfaces are complex, even in idealized systems. Our need to understand this complexity is driven by both pure and applied sciences, that is, by the need for basic understanding of earth systems and for the knowledge to mitigate our influences upon them. The second-order nonlinear optical techniques of second-harmonic generation and sum-frequency generation spectroscopy have proven adept at probing these types of interfaces. This review focuses on the contributions to geochemistry made by nonlinear optical methods. The types of questions probed have included a basic description of the structure adopted by water molecules at the mineral interface, how flow and porosity affect this structure, adsorption of trace metal and organic species, and dissolution mechanisms. We also discuss directions and challenges that lie ahead and the outlook for the continued use of nonlinear optical methods for studies of mineral-water boundaries.

Precipitates of Al(III), Sc(III), and La(III) at the muscovite-water interface

The interaction of Al(III), Sc(III), and La(III) with muscovite-water interfaces was studied at pH 4 and 10 mM NaCl using second harmonic generation (SHG) and X-ray photoelectron spectroscopy (XPS). SHG data for Sc(III) and La(III) suggest complete and/or partial irreversible adsorption that is attributed by XPS to the growth of Sc(III) and La(III) hydroxides/oxides on the muscovite surface. Al(III) adsorption appears to coincide with the growth of gibbsite (Al(OH)3) deposits on the muscovite surface, as indicated by the magnitude of the interfacial potential computed from the SHG data. This interpretation of the data is consistent with previous studies reporting the epitaxial growth of gibbsite on the muscovite surface under similar conditions. The implication of our findings is that the surface charge density of mica may change (and in the case of Al(III), even flip sign from negative (mica) to positive (gibbsite)) when Al(III), Sc(III), or La(III) is present in aqueous phases in contact with heterogeneous geochemical media rich in mica-class minerals, even at subsaturation conditions.

Y(III) interactions with guanine oligonucleotides covalently attached to aqueous/solid interfaces

The binding of Y(III) ions to surface-immobilized single-stranded 20-mers of guanine was studied using the Eisenthal χ((3)) technique and AFM. The free energy of binding for Y(III) to the G(20) sequence was found to be -39.5(8) kJ/mol. Furthermore, yttrium binds much more strongly to surface-immobilized oligonucleotides than the divalent metals previously reported. At maximum surface coverage, Y(III) ion densities range between one to three ions bound per strand. Comparatively, Mg(II) binds to the G(20)-functionalized interface in much higher ion densities. This result may be explained, in part, by the larger hydration sphere radius of Y(III) compared to that of Mg(II). The ion loading and binding free energy results, in conjunction with other surface and bulk aqueous phase studies, suggest that a fully hydrated +2 or +3 yttrium ion binds to the oligonucleotides through an outer-sphere mechanism. Tapping mode AFM results indicate that oligonucleotide height does not appreciably decrease following Y(III) binding. These results, together with the low ion densities for Y(III) ions, indicate that Y(III) strand loading may not significantly decrease the intrastrand Coulombic repulsions in order to cause a significant decrease in oligomer height.

Uranyl adsorption at the muscovite (mica)/water interface studied by second harmonic generation

Uranyl adsorption at the muscovite (mica)/water interface was studied by second harmonic generation (SHG). Using the nonresonant χ(3) technique and the Gouy-Chapman model, the initial surface charge density of the mica surface was determined to be -0.022(1) C/m(2) at pH 6 and in the presence of dissolved carbonate. Under these same conditions, uranyl adsorption isotherms collected using nonresonant χ(3) experiments and resonantly enhanced SHG experiments that probe the ligand-to-metal charge transfer bands of the uranyl cation yielded a uranyl binding constant of 3(1) × 10(7) M(-1), corresponding to a Gibbs free energy of adsorption of -52.6(8) kJ/mol, and a maximum surface charge density at monolayer uranyl coverage of 0.028(3) C/m(2). These results suggest favorable adsorption of uranyl ions to the mica interface through strong ion-dipole or hydrogen interactions, with a 1:1 uranyl ion to surface site ratio that is indicative of monovalent cations ((UO(2))(3)(OH)(5)(+), (UO(2))(4)(OH)(7)(+), UO(2)OH(+), UO(2)Cl(+), UO(2)(CH(3)COO(-))(+)) binding at the interface, in addition to neutral uranyl species (UO(2)(OH)(2) and UO(2)CO(3)). This work provides benchmark measurements to be used in the improvement of contaminant transport modeling.