Using synchrotron X-ray and neutron methods to investigate structural aspects of metal ion solvation and solution structure: An approach using empirical potential structure refinement

On the Elusive Links between Solution Microstructure, Dynamics, and Solvation Thermodynamics: Demystifying the Path through a Bridge over Troubled Conjectures and Misinterpretations

We build a fundamentally based bridge between the solute-induced microstructural perturbation of the species environment and the dynamic as well as thermodynamic responses of the fluid system, regardless of the state conditions, composition, nature of the solvent, and either the magnitude or the type of solute-solvent intermolecular-interaction asymmetries. For that purpose, we advance a fluctuation-based solvation formalism of fluid mixtures to provide meaningful descriptors of solvation phenomena, the microstructural signatures of their solute-solvent intermolecular interaction asymmetry, and the thermodynamic manifestations linked to the solution nonideality. The rigorous foundations afford us to address some crucial issues frequently invoked in the literature including the microstructural perturbation domain, its proper identification and molecular-based meaning toward the interpretation of the solvation process, and the potential impact of the local differential behavior between anions and cations on the actual salt-induced perturbation of the solvent microstructure. Indeed, we link the precisely characterized species solvation behavior to fundamental thermodynamic residual-property relations, and the dynamics associated with either the viscous flow or diffusive behavior of the solvent, to finally illustrate their outcome with experimental data of aqueous electrolyte solutions from the available literature. Ultimately, this effort provides a highly desirable unambiguous identification of the cause-effect connections between the microstructurally perturbed domains and the experimentally measured macroscopic solvation properties, including their effect on the dynamics of the solvent environment. More importantly, it lends a well-established solvation framework to bridge rigorously the microstructural details of the mixture, its dynamics, and its solvation thermodynamics to enhance our understanding of well-defined ranked Hofmeister series, i.e., by avoiding ad hoc conjectures and unsupported microscopic interpretations of solvation phenomena.

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 Cr³⁺ along with the more ambiguous structure of aqueous Eu³⁺ 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.

Insights into the solution structure of the hydrated uranyl ion from neutron scattering and EXAFS experiments

The solution structure of 1.0 M Uranyl Chloride has been determined by the EPSR modelling of a combination of neutron scattering and EXAFS data. The experimental data show an equilibrium in solution between [UO₂(H₂O)₅]²⁺ and [UO₂Cl(H₂O)₄]⁺ with a stability constant of 0.23 ± 0.03 mol^-1 dm^-3. A much smaller fraction of the neutral [UO₂Cl₂(H₂O)₃] ion is also observed. The data also show, for the first time in solution, that the uranyl ion is a very poor hydrogen bond acceptor, but the coordinated waters show enhanced hydrogen bond ability compared to the bulk water.

EXAFS investigations of temperature-dependent structure in cobalt-59 molecular NMR thermometers

Cobalt-59 nuclei are known for extremely thermally sensitive chemical shifts (δ), which in the long term could yield novel magnetic resonance thermometers for bioimaging applications. In this manuscript, we apply extended X-ray absorption fine structure (EXAFS) spectroscopy for the first time to probe the exact variations in physical structure that produce the exceptional thermal sensitivity of the 59Co NMR chemical shift. We apply this spectroscopic technique to five Co(iii) complexes: [Co(NH3)6]Cl3 (1), [Co(en)3]Cl3 (2) (en = ethylenediamine), [Co(tn)3]Cl3 (3) (tn = trimethylenediamine), [Co(tame)2]Cl3 (4) (tame = 1,1,1-tris(aminomethyl)ethane), and [Co(diNOsar)]Cl3 (5) (diNOsar = dinitrosarcophagine). The solution-phase EXAFS data reveal increasing Co-N bond distances for these aqueous complexes over a ∼50 °C temperature window, expanding by Δr(Co-N) = 0.0256(6) Å, 0.0020(5) Å, 0.0084(5) Å, 0.0006(5) Å, and 0.0075(6) Å for 1-5, respectively. Computational analyses of the structural changes reveal that increased connectivity between the donor atoms encourages complex structural variations. These results imply that rich temperature-dependent structural variations define 59Co NMR thermometry in macrocyclic complexes.

The structural elucidation of aqueous H3BO3 solutions by DFT and neutron scattering studies

The micro-structure of aqueous boric acid (H3BO3) solutions is of broad interest in earth sciences, geochemistry, material science, as well as chemical engineering. In the present study, the structure of aqueous H3BO3 solutions was studied via neutron scattering with 2H and 11B isotope labelling combined with empirical potential structure refinement (EPSR) modelling. In aqueous H3BO3 solutions, B(OH)3 is the dominant borate species. Density function theory (DFT) calculations show that the boron hydroxyl has a lower electrostatic potential (ESP), which makes B(OH)3 a relatively weakly hydrated, compared with the bulk water. In the 0.95 mol L-1 H3BO3 solution at 298 K (saturated), ∼18 water molecules enter the hydration sphere of B(OH)3 with the hydration distance (B-O(W)) of 3.75 Å, while only 4.23 of them hydrate with H3BO3 as the hydrogen bond (H-bond) acceptor or H-bond donor. Both neutron scattering and DFT calculations for 2B(OH)3·6H2O clusters at the ωB97XD/6-311++g(3df,3pd) basis level show that B(OH)3 forms molecular clusters in bidentate contact molecular pairs (BCMP), mono-dentate molecular pairs (MCMP), solvent-shared molecular pairs (SMP), and parallel solvent-shared molecular pairs (PSMP) in aqueous solutions. Their relative contents are both concentration- and temperature-sensitive. BCMP with the B-B distance of ∼4.1 Å is the dominant molecular pair in the aqueous solutions. Relatively less content and van der Waals interactions stabilized PSMP, with a B-B distance of ∼3.6 Å between the two parallel layers, which is a crucial species for the crystallization of H3BO3 from aqueous solution.

Ion hydration and association in aqueous potassium tetrahydroxyborate solutions

Potassium tetrahydroxyborate solution is a significant material in the borate solution family, but there is limited knowledge about hydration structures and interactions of K⁺, [B(OH)₄^-], and water. In this study, the X-ray diffraction measurements of potassium tetrahydroxyborate solutions have been made. The experimental structure factors are subjected to Empirical Potential Structure Refinement (EPSR) modeling to reveal the details of ion hydration and association in the aqueous solutions. This study shows that the O(W)-O(W) distance of water in the studied solutions ranges from 2.82 to 2.76 Å with a coordination number that ranges from 4.7 ± 1.4 to 3.1 ± 1.3 when the value of the water-salt molar ratio (WSR) is decreased from 30 to 6. The addition of ions slightly affects the tetrahedral structure of water even when the concentration of ions is high. The first hydration distance of K⁺ remained at ∼2.67 Å, whereas the value of the coordination number (CN) decreased from 5.4 ± 1.3 to 3.9 ± 1.5 when the concentration of the borate solution was increased. The hydration ability of K⁺ was weak and almost did not have a fixed local hydration structure. The pair distribution function (PDF) of g_B-O(W)(r) shows that [B(OH)₄^-] has a broad hydration distance from 2.9 to 5.4 Å because of the complex interactive relationship between K⁺, [B(OH)₄^-] and water. There is a competitive hydration between K⁺ and [B(OH)₄^-]. Both the X-ray diffraction and DFT-based calculations confirm that the main species is monodentate contact ion pairs when WSR = 30, bidentate contact ion pairs when WSR = 14, and triple contact ion pairs when WSR = 6. These results will provide a new understanding about potassium tetrahydroxyborate solution.

Ion Hydration and Ion Pairing as Probed by THz Spectroscopy

Ion hydration is of pivotal importance for many fundamental processes. Various spectroscopic methods are used to study the retardation of the hydration bond dynamics in the vicinity of anions and cations. Here we introduce THz-FTIR spectroscopy as a powerful method to answer the open questions. We show through dissection of THz spectra that we can pinpoint characteristic absorption features that can be attributed to the rattling modes of strongly hydrating ions within their hydration cages as well as vibrationally induced charge fluctuations in the case of weakly hydrating ions. Further analysis yields information on anion-cation cooperativity, the size of the dynamic hydration shell, as well as the lifetimes of these collective ion-hydration water modes and their connecting thermal bath states. Our study provides evidence for a non-additive behavior, thus questioning the simplified Hofmeister model. THz spectroscopy enables ion pairing to be observed and quantified at a high salt concentration.

Using EXAFS data to improve atomistic structural models of glasses

Quantitative characterization of the atomic structure of multi-component glasses is a long-standing scientific challenge. This is because in most cases no single experimental technique is capable of completely resolving all aspects of a disordered system's structure. In this situation, the most practical solution for the materials scientist is to apply multiple experimental probes offering differing degrees of insight into a material's properties. This powerful and widely adopted approach does, however, transfer the characterization challenge to the task of developing a coherent data analysis framework that can appropriately combine the diverse experimental insight into a single, data-consistent, structural model. Here, taking a terbium metaphosphate glass as an example system, it is illustrated how this can be achieved for X-ray diffraction and extended X-ray absorption fine-structure (EXAFS) spectroscopy data, using an empirical potential structure refinement approach. This methodology is based on performing a Monte Carlo simulation of the structure of a disordered material that is guided to a solution consistent with the provided experimental data, by a series of pairwise perturbation potentials operating on a classical reference potential foundation. For multi-component glasses the incorporation of EXAFS data into the resulting bulk structural models is shown to make a critical contribution that is required to properly account for the increase in local structural order that can develop in the melt-quench process of glass formation.

Using neutrons, X-rays and nuclear magnetism to determine the role of transition metal oxide inclusions on both glass structure and stability in automotive glass enamels

Low levels of transition metal oxides in alkali borosilicate glass systems can drastically influence crystallisation and phase separation properties. We investigated the non-monotonic effect of manganese doping on suppressing crystallisation, and the influence on optical properties by iron oxide doping, in terms of local atomic structure. Structural models based on empirical potential structure refinement were generated from neutron and X-ray scattering data, and compared against multinuclear solid-state NMR. This revealed that a 2.5% manganese doping had a disruptive effect on the entire glass network, supressing crystallisation of an undesired bismuth silicate phase, and that iron species preferentially locate near borate tetrahedra. Preventing phase separation and controlling crystallisation behaviour of glass are critical to the ultimate properties of automotive glass enamels.

A carbon nanopore model to quantify structure and kinetics of ion electrosorption with in situ small-angle X-ray scattering

In situ small-angle X-ray scattering was carried out on a custom-built supercapacitor cell and is presented together with a novel data analysis strategy to study the structure and kinetics of ion electrosorption in a nanoporous carbon electrode.

Rotation of X-ray polarization in the glitches of a silicon crystal monochromator

EXAFS studies on dilute samples are usually carried out by collecting the fluorescence yield using a large-area multi-element detector. This method is susceptible to the 'glitches' produced by all single-crystal monochromators. Glitches are sharp dips or spikes in the diffracted intensity at specific crystal orientations. If incorrectly compensated, they degrade the spectroscopic data. Normalization of the fluorescence signal by the incident flux alone is sometimes insufficient to compensate for the glitches. Measurements performed at the state-of-the-art wiggler beamline I20-scanning at Diamond Light Source have shown that the glitches alter the spatial distribution of the sample's quasi-elastic X-ray scattering. Because glitches result from additional Bragg reflections, multiple-beam dynamical diffraction theory is necessary to understand their effects. Here, the glitches of the Si(111) four-bounce monochromator of I20-scanning just above the Ni K edge are associated with their Bragg reflections. A fitting procedure that treats coherent and Compton scattering is developed and applied to a sample of an extremely dilute (100 micromolal) aqueous solution of Ni(NO3)2. The depolarization of the wiggler X-ray beam out of the electron orbit is modeled. The fits achieve good agreement with the sample's quasi-elastic scattering with just a few parameters. The X-ray polarization is rotated up to ±4.3° within the glitches, as predicted by dynamical diffraction. These results will help users normalize EXAFS data at glitches.