Analysis of the First Ion Coordination Sphere: A Toolkit to Analyze the Coordination Sphere of Ions.

Rapid and accurate approaches to characterizing the coordination structure of an ion are important for designing ligands and quantifying structure-property trends. Here, we introduce AFICS (Analysis of the First Ion Coordination Sphere), a tool written in Python 3 for analyzing the structural and geometric features of the first coordination sphere of an ion over the course of molecular dynamics simulations. The principal feature of AFICS is its ability to quantify the distortion a coordination geometry undergoes compared to uniform polyhedra. This work applies the toolkit to analyze molecular dynamics simulations of the well-defined coordination structure of aqueous Cr3+ along with the more ambiguous structure of aqueous Eu3+ chelated to ethylenediaminetetraacetic acid. The tool is targeted for analyzing ions with fluxional or irregular coordination structures (e.g., solution structures of f-block elements) but is generalized such that it may be applied to other systems.

Water Defect Stabilizes the Bi3+ Lone-Pair Electronic State Leading to an Unusual Aqueous Hydration Structure.

The aqueous hydration structure of the Bi3+ ion is probed using a combination of extended X-ray absorption fine structure (EXAFS) spectroscopy and density functional theory (DFT) simulations of ion-water clusters and condensed-phase solutions. Anomalous features in the EXAFS spectra are found to be associated with a highly asymmetric first-solvent water shell. The aqueous chemistry and structure of the Bi3+ ion are dramatically controlled by the water stabilization of a lone-pair electronic state involving the mixed 6s and 6p orbitals. This leads to a distinct multimodal distribution of water molecules in the first shell that are separated by about 0.2 Å. The lone-pair structure is stabilized by a collective response of multiple waters that are localized near the lone-pair anti-bonding site. The findings indicate that the lone-pair stereochemistry of aqueous Bi3+ ions plays a major role in the binding of water and ligands in aqueous solutions.

Ionic Contraction across the Lanthanide Series Decreases the Temperature-Induced Disorder of the Water Coordination Sphere.

In liquid, temperature affects the structures of lanthanide complexes in multiple ways that depend upon complex interactions between ligands, anions, and solvent molecules. The relative simplicity of lanthanide aqua ions (Ln3+) make them well suited to determine how temperature induces structural changes in lanthanide complexes. We performed a combination of ab initio molecular dynamics (AIMD) simulations and extended X-ray absorption fine structure (EXAFS) measurements, both at 25 and 90 °C, to determine how temperature affects the first- and second-coordination spheres of three Ln3+ (Ce3+, Sm3+, and Lu3+) aqua ions. AIMD simulations show first lanthanide coordination spheres that are similar at 25 and 90 °C, more so for the Lu3+ ion that remains as eight-coordinate than for the Ce3+ and Sm3+ ions that change their preferred coordination number from nine (at 25 °C) to eight (at 90 °C). The measured EXAFS spectra are very similar at 25 and 90 °C, for the Ce3+, Sm3+, and Lu3+ ions, suggesting that the dynamical disorder of the Ln3+ ions in liquid water is sufficient such that temperature-induced changes do not clearly manifest changes in the structure of the three ions. Both AIMD simulations and EXAFS measurements show very similar structures of the first coordination sphere of the Lu3+ ion at 25 and 90 °C.

Computational Prediction of All Lanthanide Aqua Ion Acidity Constants.

The protonation state of lanthanide-ligand complexes, or lanthanide-containing porous materials, with many Brønsted acid sites can change due to proton loss/gain reactions with water or other heteroatom-containing compounds. Consequently, variations in the protonation state of lanthanide-containing species affect their molecular structure and desired properties. Lanthanide(III) aqua ions undergo hydrolysis and form hydroxides; they are the best characterized lanthanide-containing species with multiple Brønsted acid sites. We employed constrained ab initio molecular dynamics simulations and electronic structure calculations to determine all acidity constants of the lanthanide(III) aqua ions solely from computation. The first, second, and third acidity constants of lanthanide(III) aqua ions were predicted, on average, within 1.2, 2.5, and 4.7 absolute pKa units from experiment, respectively. A table includes our predicted pKa values alongside most experimentally measured pKa values known to date. The approach presented is particularly suitable to determine the Brønsted acidity of lanthanide-containing systems with multiple acidic sites, including those whose measured acidity constants cannot be linked to specific acid sites.

Coordination Sphere of Lanthanide Aqua Ions Resolved with Ab Initio Molecular Dynamics and X-ray Absorption Spectroscopy.

To resolve the fleeting structures of lanthanide Ln3+ aqua ions in solution, we (i) performed the first ab initio molecular dynamics (AIMD) simulations of the entire series of Ln3+ aqua ions in explicit water solvent using pseudopotentials and basis sets recently optimized for lanthanides and (ii) measured the symmetry of the hydrating waters about Ln3+ ions (Nd3+, Dy3+, Er3+, Lu3+) for the first time with extended X-ray absorption fine structure (EXAFS). EXAFS spectra were measured experimentally and generated from AIMD trajectories to directly compare simulation, which concurrently considers the electronic structure and the atomic dynamics in solution, with experiment. We performed a comprehensive evaluation of EXAFS multiple-scattering analysis (up to 6.5 Å) to measure Ln-O distances and angular correlations (i.e., symmetry) and elucidate the molecular geometry of the first hydration shell. This evaluation, in combination with symmetry-dependent L3- and L1-edge spectral analysis, shows that the AIMD simulations remarkably reproduces the experimental EXAFS data. The error in the predicted Ln-O distances is less than 0.07 Å for the later lanthanides, while we observed excellent agreement with predicted distances within experimental uncertainty for the early lanthanides. Our analysis revealed a dynamic, symmetrically disordered first coordination shell, which does not conform to a single molecular geometry for most lanthanides. This work sheds critical light on the highly elusive coordination geometry of the Ln3+ aqua ions.