Structure of the Homoleptic Thorium(IV) Aqua Ion [Th(H2O)10]Br4

Prediction of Redox Potentials for Ac, Th, and Pa in Aqueous Solution

Density functional theory in conjunction with small core pseudopotentials and the associated basis sets was used to calculate potentials for multiple redox couples, covering a range of oxidation states for Ac (0 to III), Th (0 to IV), and Pa (0 to V) in aqueous solution. Solvation effects were incorporated using a supermolecule-continuum approach, with 30 water molecules representing two solvation shells, and the COSMO and SMD implicit solvation models. The calculated geometries for Ac(III), Th(IV), and Pa(V) were in reasonable agreement with the available experimental data. Using the COSMO model with the B3LYP functional, the calculated redox potentials were within ±0.2 V from experiment for most redox couples. Several pathways were explored for the Pa(V/IV) redox couple for different forms of Pa(V) and Pa(IV). Most Pa(V/IV) redox couples have very similar potentials, ranging from 0 to -0.4 V up to a pH of 1.4. At pH = 1.4, the potentials shift to values that are more negative than -0.7 V, reflecting the growing unfavorable nature of the redox process at higher pH levels. The calculated values for An(III/II) potentials were consistent with prior estimates and the available experimental data. The predicted redox potentials for An(II/I) were highly negative, as expected. For An(I/0) potentials, Th and Pa exhibited positive values, contrasting with the negative values calculated for Ac. The An+m/An(0) potentials agreed better with the experimental data when using the COSMO solvation model as compared to the SMD model.

Structures of Th4+ aqueous solutions: insights from AIMD and metadynamics simulations

Solution chemistry of actinide ions is critical to understanding the solvation behaviors and hydrolysis process. Using tetravalent thorium ion Th4+ as a representative example, we investigate the local structures and dynamic behaviors of hydrated Th4+ ions by ab initio molecular dynamics (AIMD) simulations using the recently developed norm-conserving pseudopotentials and basis sets optimized for actinides (J.-B. Lu et al., J. Chem. Theory Comput. 2021, 17, 3360-3371). AIMD simulations reveal two distinct solvation shells, with the first shell comprising 9 water molecules at approximately rTh-O = 2.50 Å and exhibiting a tricapped trigonal prism geometry. These conclusions are confirmed through metadynamics simulations and further structural analysis. AIMD simulations also show the slight effect of temperature and counterions on the structure of the solution. The structured solvation shells of the highly charged Th4+ ion with the specific geometry, distinct from the structure of liquid water, lead to corresponding structural changes in the hydrogen bond network in water. Additionally, beyond the solvent-shared ion pair (SIP) state observed in the unbiased AIMD simulations, the metadynamics simulations reconstruct a two-dimensional free energy surface that clearly indicates the potential stability of the contact ion pair (CIP) state in the system with Cl- as a counterion. The findings in this work provide insights into the solution chemistry of actinides and serve as a reference for studying other actinide solution systems.

Structure and size of complete hydration shells of metal ions and inorganic anions in aqueous solution

The structures of nine hydrated metal ions in aqueous solution have been redetermined by large angle X-ray scattering to obtain experimental data of better quality than those reported 40-50 years ago. Accurate M-OI and M-(OI-H)⋯OII distances and M-OI(H)⋯OII bond angles are reported for the hydrated magnesium(II), aluminium(III), manganese(II), iron(II), iron(III), cobalt(II), nickel(II), copper(II) and zinc(II) ions; the subscripts I and II denote oxygen atoms in the first and second hydration sphere, respectively. Reported structures of hydrated metal ions in aqueous solution are summarized and evaluated with emphasis on a possible relationship between M-OI-OII bond angles and bonding character. Metal ions with high charge density have M-OI-OII bond angles close to 120°, indicative of a mainly electrostatic interaction with the oxygen atom in the water molecule in the first hydration shell. Metal ions forming bonds with a significant covalent contribution, as e.g. mercury(II) and tin(II), have M-OI-OII bond angles close to 109.5°. This implies that they bind to one of the free electron pairs in the water molecule. Comparison of M-O bond distances of hydrated metal ions in the solid state with one hydration shell, and in aqueous solution with in most cases at least two hydration shells, shows no significant differences. On the other hand, the X-O bond distance in hydrated oxoanions increases by ca. 0.02 Å in aqueous solution in comparison with the corresponding X-O distance in the solid state. A linear correlation is observed between volume, calculated from the van der Waals radius of the hydrated ion, and the ionic diffusion coefficient in aqueous solution. This correlation strongly indicates that monovalent metal ions, except lithium and silver(I), and singly-charged monovalent oxoanions have a single hydration shell. Divalent metal ions, bismuth(III) and the lanthanoid(III) and actinoid(III) ions have two hydration shells. Trivalent transition and tetravalent metal ions have two full hydration shells and portion of a third one. Doubly charged oxoanions have one well-defined hydration shell and an ill-defined second one.

Spatiotemporal Studies of Soluble Inorganic Nanostructures with X-rays and Neutrons

This Review addresses the use of X-ray and neutron scattering as well as X-ray absorption to describe how inorganic nanostructured materials assemble, evolve, and function in solution. We first provide an overview of techniques and instrumentation (both large user facilities and benchtop). We review recent studies of soluble inorganic nanostructure assembly, covering the disciplines of materials synthesis, processes in nature, nuclear materials, and the widely applicable fundamental processes of hydrophobic interactions and ion pairing. Reviewed studies cover size regimes and length scales ranging from sub-Ångström (coordination chemistry and ion pairing) to several nanometers (molecular clusters, i.e. polyoxometalates, polyoxocations, and metal-organic polyhedra), to the mesoscale (supramolecular assembly processes). Reviewed studies predominantly exploit 1) SAXS/WAXS/SANS (small- and wide-angle X-ray or neutron scattering), 2) PDF (pair-distribution function analysis of X-ray total scattering), and 3) XANES and EXAFS (X-ray absorption near-edge structure and extended X-ray absorption fine structure, respectively). While the scattering techniques provide structural information, X-ray absorption yields the oxidation state in addition to the local coordination. Our goal for this Review is to provide information and inspiration for the inorganic/materials science communities that may benefit from elucidating the role of solution speciation in natural and synthetic processes.

Stabilization of Pu(IV) in PuBr4(OPCy3)2 and Comparisons with Structurally Similar ThX4(OPR3)2 (R = Cy, Ph) Molecules

The pursuit of a trivalent plutonium halide phosphine oxide compound, e.g., "PuBr3(OPR)3," instead led to the isolation of the tetravalent trans-PuIVBr4(OPCy3)2, PuBr/Cy, compound by spontaneous oxidation of PuIII. The donating nature of phosphine oxides has allowed the isolation and characterization of PuBr/Cy by crystallographic, multinuclear NMR, solid state, and solution phase UV-vis-NIR spectroscopic techniques. The presence of a putative plutonyl(VI) complex formulated as "trans-PuVIO2Br2(OPCy3)2" was also observed spectroscopically and tentatively by single-crystal X-ray diffraction as a cocrystal of PuBr/Cy. A series of trans-ThX4(OPCy3)2 (X = Cl, ThCl/Cy; Br, ThBr/Cy; I, ThI/Cy) complexes were synthesized for comparison to PuBr/Cy. The triphenylphosphine oxide, OPPh3, complexes, trans-AnI4(OPPh3)2 (An = Th, ThI/Ph; U, UI/Ph), were also synthesized for comparison, completing the series trans-UX4(OPPh3)2 (X = Cl, Br, I), UX/Ph. To enable the synthesis of ThI/Cy and ThI/Ph, a new nonaqueous thorium iodide starting material, ThI4(Et2O)2, was synthesized. The syntheses of organic solvent soluble ThI4L2 (L = Et2O, OPCy3, and OPPh3) are the first examples of crystallographically characterized neutral thorium tetraiodide materials beyond binary ThI4. To show the viability of ThI4(Et2O)2 as a starting material for organothorium chemistry, (C5Me4H)3ThI was synthesized and crystallographically characterized.

Impact of Noncovalent Interactions on the Structural Chemistry of Thorium(IV)-Aquo-Chloro Complexes

Five novel tetravalent thorium (Th) compounds that consist of Th(H₂O)_xCl_y structural units were isolated from acidic aqueous solutions using a series of nitrogen-containing heterocyclic hydrogen (H) bond donors. Taken together with three previously reported phases, the compounds provide a series of monomeric Th^IV complexes wherein the effects of noncovalent interactions (and H-bond donor identity) on Th structural chemistry can be examined. Seven distinct structural units of the general formulas [Th(H₂O)_xCl_8-x]^x-4 (x = 2, 4) and [Th(H₂O)_xCl_9-x]^x-5 (x = 5-7) are described. The complexes range from chloride-deficient [Th(H₂O)₇Cl₂]²⁺ to chloride-rich [Th(H₂O)₂Cl₆]^2- species, and theory was used to understand the relative energies that separate complexes within this series via the stepwise chloride addition to an aquated Th cation. Electronic structure theory predicted the reaction energies of chloride addition and release of water through a series of transformations, generally highlighting an energetic driving force for chloride complexation. To probe the role of the counterion in the stabilization of these complexes, electrostatic potential (ESP) surfaces were calculated. The ESP surfaces indicated a dependence of the chloride distribution about the Th metal center on the pK_a of the countercation, highlighting the directing effects of noncovalent interactions (e.g., Hbonding) on Th speciation.

Periodic trends and complexation chemistry of tetravalent actinide ions with a potential actinide decorporation agent 5-LIO(Me-3,2-HOPO): A relativistic density functional theory exploration

A relativistic density functional theory (DFT) study is reported which aims to understand the complexation chemistry of An⁴⁺ ions (An = Th, U, Np, and Pu) with a potential decorporation agent, 5-LIO(Me-3,2-HOPO). The calculations show that the periodic change of the metal binding free energy has an excellent correlation with the ionic radii and such change of ionic radii also leads to the structural modulation of actinide-ligand complexes. The calculated structural and binding parameters agree well with the available experimental data. Atomic charges derived from quantum theory of atoms in molecules (QTAIM) and natural bond order (NBO) analysis shows the major role of ligand-to-metal charge transfer in the stability of the complexes. Energy decomposition analysis, QTAIM, and electron localization function (ELF) predict that the actinide-ligand bond is dominantly ionic, but the contribution of orbital interaction is considerable and increases from Th⁴⁺ to Pu⁴⁺ . A decomposition of orbital contributions applying the extended transition state-natural orbital chemical valence method points out the significant π-donation from the oxygen donor centers to the electron-poor actinide ion. Molecular orbital analysis suggests an increasing trend of orbital mixing in the context of 5f orbital participation across the tetravalent An series (Th-Pu). However, the corresponding overlap integral is found to be smaller than in the case of 6d orbital participation. An analysis of the results from the aforementioned electronic structure methods indicates that such orbital participation possibly arises due to the energy matching of ligand and metal orbitals and carries the signature of near-degeneracy driven covalency.

Properties of the tetravalent actinide series in aqueous phase from a microscopic simulation self-consistent engine

In the context of nuclear fuel recycling and environmental issues, the understanding of the properties of radio-elements with various approaches remains a challenge regarding their dangerousness. Moreover, experimentally, some issues are also of importance; first, it is imperative to work at sufficiently high concentrations to reach the sensitivities of the analytical tools, however this condition often leads to precipitation for some of them; second, stabilizing specific oxidation states of some actinides remains a challenge, thus making it difficult to extract general trends across the actinide series. Complementary to experiments, modeling can be used to unbiasedly probe the actinide's properties in an aquatic environment and offers a predictive tool. We report the first molecular dynamics simulations based on homogeneously built force fields for the whole series of the tetravalent actinides in aqueous phase from ThIV to BkIV and including PuIV. The force fields used to model the interactions among the constituents include polarization and charge donation microscopic effects. They are built from a self-consistent iterative ab initio based engine that can be included in future developments as an element of a potential machine learning procedure devoted to generating accurate force fields. The comparison of our simulated hydrated actinide properties to available experimental data shows the model robustness and the relevance of our parameter assignment engine. Moreover, our simulated structural, dynamical and evolution of the hydration free energy data show that, apart from AmIV and CmIV, the actinide properties change progressively along the series.

Density Functional Theory Investigations on the Mechanism of Formation of Pa(V) Ion in Hydrous Solutions

Due to the enormous threat of protactinium to the environment and human health, its disposal and chemistry have long been important topics in nuclear science. [PaO(H₂O)₆]3+ is proposed as the predominant species in hydrous and acidic solutions, but little is known about its formation mechanism. In this study, density functional theory (DFT) calculations demonstrate a water coordination-proton transfer-water dissociation mechanism for the formation of PaO3+ in hydrous solutions. First, Pa(V) ion preferentially forms hydrated complexes with a coordination number of 10. Through hydrogen bonding, water molecules in the second coordination sphere easily capture two protons on the same coordinated H₂O ligand to form [PaO(H₂O)₉]3+. Water dissociation then occurs to generate the final [PaO(H₂O)₆]3+, which is the thermodynamic product of Pa(V) in hydrous solutions.

The coordination chemistry of lanthanide and actinide metal ions with hydroxypyridinone-based decorporation agents: orbital and density based analyses

In the context of the mitigation of the biological effects of internal radionuclide contamination and for efficient decorporation, the design and development of efficient chelators for lanthanide and actinide metal ions has become a central issue. The pioneering work of Raymond and coworkers (Chem. Rev., 2003, 103, 4207-4282) led to the development of siderophore-related hydroxypyridinonate ligands for possible treatment of internalized radionuclides. However, the structure-function relationship of Ln/An bound to these ligands, particularly the bonding and coordination aspects are not clearly understood at the atomic level. Here, we have investigated the structure, binding and energetics of trivalent and tetravalent Ln/An (Sm3+, Eu3+, Am3+, Cm3+, Th4+, Pu4+) ions with spermine-based octadentate hydroxypyridinonate chelators, namely 3,4,3-LI(1,2-HOPO) and its 3,3,3 variant, using relativistic density functional theory (DFT). Furthermore, we have performed orbital and density based analyses to elucidate the nature of bonding in these complexes. In accordance with the experimental stability constant, we found the maximum binding free energy for An4+ (Pu4+, Th4+) as compared to trivalent metal ions. CDA and ECDA analyses along with orbital-based population analyses confirmed the higher ligand to metal charge transfer for An4+ than for trivalent metal ions. Furthermore, the aromaticity index analysis suggested the presence of crucial chelatoaromatic stabilization for all these metal ions with the maximum for An4+. QTAIM descriptors indicated that the binding of An/Ln with the hard oxygen donor of the ligands is of the donor-acceptor type but a higher degree of covalency exists for actinides as compared to lanthanides. Furthermore, QTAIM and molecular orbital analysis confirmed that such covalency is of the energy-driven type and strictly originates from the orbital mixing event of An-5f orbitals with the ligand orbitals.

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.

Structural and thermodynamic investigation of AnIVLI(O)HOPO

For the first time, capillary electrophoresis coupled with inductively coupled plasma mass spectrometry has been used to determine the stability constants of Pu^IV with the multidentate hydroxypyridinonate chelating agents LI(O)HOPO in 0.1 M NaNO₃ solution, p_cH = 1.395 at 25 °C through competition with the NTA ligand.

Structure, Stability, and Electronic Properties of Dimethyl Sulfoxide and Dimethyl Formammide Clusters Containing Th(4.)

By using accurate density functional theory calculations, we have studied the complexes of Th(4+) with dimethyl-sulfoxide (DMSO) and dimethyl-formammide (DMF) molecules. These solvents are prototypes for oxygen-donor organic environments in which the oxygen atom is connected to S and C atom, respectively. Extended structural, energetic, and electronic structure analysis has been performed to provide a complete picture of the physical properties at the basis of the interaction of Th(4+) with the two solvents. By using a cluster grow approach, we have found that, very likely, the first solvation shell contains nine molecules in the case of DMSO, while it contains eight molecules for DMF. The theoretical results shown here are in agreement with experimental data taken from the literature.

First structural characterization of Pa(iv) in aqueous solution and quantum chemical investigations of the tetravalent actinides up to Bk(IV): the evidence of a curium break

More than a century after its discovery the structure of the Pa(4+) ion in acidic aqueous solution has been investigated for the first time experimentally and by quantum chemistry. The combined results of EXAFS data and quantum chemically optimized structures suggest that the Pa(4+) aqua ion has an average of nine water molecules in its first hydration sphere at a mean Pa-O distance of 2.43 Å. The data available for the early tetravalent actinide (An) elements from Th(4+) to Bk(4+) show that the An-O bonds have a pronounced electrostatic character, with bond distances following the same monotonic decreasing trend as the An(4+) ionic radii, with a decrease of the hydration number from nine to eight for the heaviest ions Cm(4+) and Bk(4+). Being the first open-shell tetravalent actinide, Pa(4+) features a coordination chemistry very similar to its successors. The electronic configuration of all open-shell systems corresponds to occupation of the valence 5f orbitals, without contribution from the 6d orbitals. Our results thus demonstrate that Pa(iv) resembles its early actinide neighbors.

Exploring actinide materials through synchrotron radiation techniques

Synchrotron radiation (SR) based techniques have been utilized with increasing frequency in the past decade to explore the brilliant and challenging sciences of actinide-based materials. This trend is partially driven by the basic needs for multi-scale actinide speciation and bonding information and also the realistic needs for nuclear energy research. In this review, recent research progresses on actinide related materials by means of various SR techniques were selectively highlighted and summarized, with the emphasis on X-ray absorption spectroscopy, X-ray diffraction and scattering spectroscopy, which are powerful tools to characterize actinide materials. In addition, advanced SR techniques for exploring future advanced nuclear fuel cycles dealing with actinides are illustrated as well.

Parameterization of highly charged metal ions using the 12-6-4 LJ-type nonbonded model in explicit water

Highly charged metal ions act as catalytic centers and structural elements in a broad range of chemical complexes. The nonbonded model for metal ions is extensively used in molecular simulations due to its simple form, computational speed, and transferability. We have proposed and parametrized a 12-6-4 LJ (Lennard-Jones)-type nonbonded model for divalent metal ions in previous work, which showed a marked improvement over the 12-6 LJ nonbonded model. In the present study, by treating the experimental hydration free energies and ion–oxygen distances of the first solvation shell as targets for our parametrization, we evaluated 12-6 LJ parameters for 18 M(III) and 6 M(IV) metal ions for three widely used water models (TIP3P, SPC/E, and TIP4P_EW). As expected, the interaction energy underestimation of the 12-6 LJ nonbonded model increases dramatically for the highly charged metal ions. We then parametrized the 12-6-4 LJ-type nonbonded model for these metal ions with the three water models. The final parameters reproduced the target values with good accuracy, which is consistent with our previous experience using this potential. Finally, tests were performed on a protein system, and the obtained results validate the transferability of these nonbonded model parameters.

On the structure of thorium and americium adenosine triphosphate complexes

PURPOSE

The actinides are chemical poisons and radiological hazards. One challenge to better appraise their toxicity and develop countermeasures in case of exposure of living organisms is to better assess pathways of contamination. Because of the high chemical affinity of those actinide elements for phosphate groups and the ubiquity of such chemical functions in biochemistry, nucleotides and in particular adenosine triphosphate nucleotide (ATP) may be considered critical target building blocks for actinides.

MATERIALS AND METHODS

Combinations of spectroscopic techniques (Fourier transformed Infra Red [FTIR], Electrospray Ionization Mass Spectrometry [ESI-MS], and Extended X-ray Absorption Fine Structure [EXAFS]) with quantum chemical calculations have been implemented in order to assess the actinides coordination arrangement with ATP.

RESULTS

We describe and compare herein the interaction of ATP with thorium and americium; thorium(IV) as a representative of actinide(IV) like plutonium(IV) and americium(III) as a representative of all heavier actinides. In the case of thorium, an insoluble complex is readily formed. In the case of americium, a behavior identical to that described previously for lutetium has been observed with insoluble and soluble complexes.

CONCLUSIONS

The comparative study of ATP complexation with Th(IV) and Am(III) shows their ability to form insoluble complexes for which a structural model has been proposed by analogy with previously described Lu(III) complexes.

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.

A new x-ray interface and surface scattering environmental cell design for in situ studies of radioactive and atmosphere-sensitive samples

We present a novel design of a purpose-built, portable sample cell for in situ x-ray scattering experiments of radioactive or atmosphere sensitive samples. The cell has a modular design that includes two independent layers of containment that are used simultaneously to isolate the sensitive samples. Both layers of containment can be flushed with an inert gas, thus serving a double purpose as containment of radiological material (either as a solid sample or as a liquid phase) and in separating reactive samples from the ambient atmosphere. A remote controlled solution flow system is integrated into the containment system that allows sorption experiments to be performed on the diffractometer. The cell's design is discussed in detail and we demonstrate the cell's performance by presenting first results of crystal truncation rod measurements. The results were obtained from muscovite mica single crystals reacted with 1 mM solutions of Th(IV) with 0.1 M NaCl background electrolyte. Data were obtained in specular as well as off-specular geometry.