Ab initio study on the structures of Th(IV) hydrate and its hydrolysis products in aqueous solution

Combined Thermodynamic, Theoretical, and Biological Study for Investigating N-(2-Acetamido)iminodiacetic Acid as a Potential Thorium Decorporation Agent

The most effective approach to mitigate the toxic effects of internal exposure of radiometals to humans is metal-ligand (ML) chelation therapy. Thorium (Th)-induced carcinogenesis as well as other health hazards to humans as a result of chronic internal exposure necessitates the development of efficient Th-decorporating agents. In this regard, chemical and biological studies were carried out to evaluate N-(2-Acetamido)iminodiacetic acid (ADA), a comparatively cost-effective, readily available, and biologically safe complexing agent for Th decorporation. In the present work, detailed thermodynamic studies for complexation of ADA with Th(IV) have been carried out to understand Th-ADA interaction, using potentiometry, calorimetry, electrospray ionization mass spectrometry, and theoretical studies, followed by its biological assessment for Th decorporation. Thermodynamic studies revealed the formation of strong Th-ADA complexes, which are enthalpically as well as entropically favored. Interestingly, density functional theory calculations, to obtain a thermodynamically favored mode of coordination, showed the uncommon trend of lower denticity of ADA in ML than in ML2, which has been explained on the basis of stabilization of ML by hydrogen bonding. The same was also reflected in the unusual trend of enthalpy for Th-ADA complexes. Biological experiments using human erythrocytes, whole human blood, and lung cells showed good cytocompatibility and ability of ADA to significantly prevent Th-induced hemolysis. Th removal of ADA from erythrocytes, human blood, and normal lung cells was found to be comparable with that of diethylenetriamine pentaacetate (DTPA), an FDA approved decorporating agent. The present study contributed significant data about Th complexation chemistry of ADA and its Th decorporation efficacy from human erythrocytes, blood, and lung cells.

Electronic Structure and Chemical Bonding of [AmO2(H2O) n ]2+/1

Systematic americyl-hydration cations were investigated theoretically to understand the electronic structures and bonding in [(AmO2)(H2O) n ]2+/1+ (n = 1-6). We obtained the binding energy using density functional theory methods with scalar relativistic and spin-orbit coupling effects. The geometric structures of these species have been investigated in aqueous solution via an implicit solvation model. Computational results reveal that the complexes of five equatorial water molecules coordinated to americyl ions are the most stable due to the enhanced ionic interactions between the AmO2 2+/1+ cation and multiple oxygen atoms as electron donors. As expected, Am-Owater bonds in such series are electrostatic in nature and contain a generally decreasing covalent character when hydration number increases.

Complexation of thorium with pyridine monocarboxylate-N-oxides: Thermodynamic and computational studies

The feed wastes and waste water treatment plants are the major sources for the entry of N-oxides into the soils then to aquatic life. The complexation of actinides with potentially stable anthropogenic ligands facilitate the transportation and migration of the actinides from the source confinement. The present study describes the determination of thermodynamic parameters for the complexation of Th(IV) with the three isomeric pyridine monocarboxylates (PCNO) namely picolinic acid-N-oxide (PANO), nicotinic acid-N-oxide (NANO) and isonicotinic acid-N-oxide (IANO). The potentiometric and isothermal calorimetric titrations were carried out to determine the stability and enthalpy of the formations for all the Th(IV)-PCNO complexes. Th-PANO complexes are more stable than Th-NANO and Th-IANO complexes which can be attributed to chelate formation in the former complexes. Formation of all the Th-PCNO complexes are endothermic and are entropy driven. The geometries for all the predicted complexes are optimized the energies, bond distances and charges on individual atoms are obtained using TURBOMOLE software. The theoretical calculation corroborated the experimental determinations.

Computational study of Th(4+) and Np(4+) hydration and hydrolysis of Th(4+) from first principles

The aqueous solvation of Th and Np in the IV oxidation state was examined using cluster models generated by Monte Carlo simulations and density functional theory embedded within the COSMO continuum model to approximate the effect of bulk water. Our results suggest that the coordination number (CN) for both Th(IV) and NP(IV) should be 9, in accordance to some experimental and theoretical results from the literature. The structural values for average oxygen-metal distances are within 0.01 Å compared to experimental data, and also within the experimental error. The calculated ΔG _Sol⁰ are in very good agreement with experimental reported values, with deviations at CN = 9 lower than 1% for both Th(IV) and Np(IV). The hydrolysis constants are also in very good agreement with experimental values. Finally, this [corrected] methodology has the advantage of using a GGA functional (BP86) that not only makes the calculations more affordable computationally than hybrid functional or ab initio molecular dynamics simulations (Car-Parrinello) calculations, but also opens the perspective to use resolution of identity (RI) calculations for more extended systems.

Structures and formation mechanisms of aquo/hydroxo oligomeric beryllium in aqueous solution: a density functional theory study

The structures and formation mechanisms of a wide variety of aquo/hydroxo oligomeric beryllium clusters were investigated using density functional theory. The structural parameters of beryllium clusters were found to vary regularly with the stepwise substitution of bound water molecules in the inner coordination sphere by hydroxyl groups. According to the Gibbs free energies deduced from SMD solvation model computations, unhydrolyzed oligomeric beryllium species are the most favorable products of polymerization, independent of the degrees of hydrolysis of the reactants. Simulation of the formation processes of oligomeric beryllium showed that polymerization, in essence, involves the nucleophilic attack of a terminal hydroxyl group in one BeO4 tetrahedron on the beryllium center in another BeO4 tetrahedron, leading to the bridging of two BeO4 tetrahedrons by a hydroxyl group.

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.

Solution structures, stabilities, kinetics, and dynamics of DO3A and DO3A-sulphonamide complexes

The Gd(3+)-DO3A-arylsulphonamide (DO3A-SA) complex is a promising pH-sensitive MRI agent. The stability constants of the DO3A-SA and DO3A complexes formed with Mg(2+), Ca(2+), Mn(2+), Zn(2+), and Cu(2+) ions are similar, whereas the logKLnL values of Ln(DO3A-SA) complexes are 2 orders of magnitude higher than those of the Ln(DO3A) complexes. The protonation constant (log KMHL) of the sulphonamide nitrogen in the Mg(2+), Ca(2+), Mn(2+), Zn(2+), and Cu(2+) complexes is very similar to that of the free ligand, whereas the logKLnHL values of the Ln(DO3A-SA) complexes are lower by about 4 logK units, indicating a strong interaction between the Ln(3+) ions and the sulphonamide N atom. The Ln(HDO3A-SA) complexes are formed via triprotonated *Ln(H3DO3A-SA) intermediates which rearrange to the final complex in an OH(-)-assisted deprotonation process. The transmetalation reaction of Gd(HDO3A-SA) with Cu(2+) is very slow (t1/2 = 5.6 × 10(3) h at pH = 7.4), and it mainly occurs through proton-assisted dissociation of the complex. The (1)H and (13)C NMR spectra of the La-, Eu-, Y-, and Lu(DO3A-SA) complexes have been assigned using 2D correlation spectroscopy (COSY, EXSY, HSQC). Two sets of signals are observed for Eu-, Y-, and Lu(DO3A-SA), showing two coordination isomers in solution, that is, square antiprismatic (SAP) and twisted square antiprismatic (TSAP) geometries with ratios of 86-14, 93-7, and 94-6%, respectively. Line shape analysis of the (13)C NMR spectra of La-, Y- , and Lu(DO3A-SA) gives higher rates and lower activation entropy values compared to Ln(DOTA) for the arm rotation, which indicates that the Ln(DO3A-SA) complexes are less rigid due to the larger flexibility of the ethylene group in the sulphonamide pendant arm. The fast isomerization and the lower activation parameters of Ln(DO3A-SA) have been confirmed by theoretical calculations in vacuo and by using the polarizable continuum model. The solid state X-ray structure of Cu(H2DO3A-SA) shows distorted octahedral coordination. The coordination sites of Cu(2+) are occupied by two ring N- and two carboxylate O-atoms in equatorial position. The other two ring N-atoms complete the coordination sphere in axial positions. The solid state structure also indicates that a carboxylate O atom and the sulphonamide nitrogen are protonated and noncoordinated.

Understanding the structure and electronic properties of Th4+-water complexes

We present a systematic quantum chemistry study of the [Th(H2O)n]4+ (n = 1 to 10) complexes to gain insight into their electronic structure and properties: the effect of the ligand distribution on the valence shells of the thorium(IV) ion is studied by means of the electron localization function (ELF) topological analysis. Particular care is given to the study of the mono-aqua complex both at its equilibrium geometry, using various tools such as energy decomposition analyses (EDA), and along its dissociation pathway. Indeed, as several electronic states cross the Th4 +-H2O0 ground state along the minimum energy path, we demonstrate that the diabatic representation implemented in MOLPRO is able to generate reference potential energy surfaces that will lead to the evaluation of diabatic dissociation curves. The calculated diabatic interaction energy curve will allow for a consistent parameterization of new generation force fields dedicated to heavy metals based on quantum chemistry.

Structure and dynamics of the U4+ ion in aqueous solution: an ab initio quantum mechanical charge field molecular dynamics study

The structure and dynamics of the stable four-times positively charged uranium(IV) cation in aqueous solution have been investigated by ab initio quantum mechanical charge field (QMCF) molecular dynamics (MD) simulation at the Hartree-Fock double-zeta quantum mechanical level. The QMCF-MD approach enables investigations with the accuracy of a quantum mechanics/molecular mechanics approach without the need for the construction of solute-solvent potentials. Angular distribution functions; radial distribution functions; coordination numbers of the first, second, and third shell (9, 19, and 44, respectively); coordination number distribution functions; tilt- and Theta-angle distribution functions; as well as local density corrected triangle distribution functions have been employed for the evaluation of the hydrated ion's structure. Special attention was paid to the determination of the geometry of the first hydration layer, and the results were compared to experimental large-angle X-ray scattering and extended X-ray absorption fine structure data. The solvent dynamics around the ion were also investigated using mean ligand residence times and related data and, resulting from the unavailability of any experimental data, were compared to ions with similar properties.

Sorption of Th(IV) onto iron corrosion products: EXAFS study

Long-term performance assessment of nuclear waste repositories is affected by the ability of the outer barrier systems to retain radionuclides after possible corrosive leakage of waste containers. The mobility of the radionuclides released from the spent fuel depends strongly on the processes that take place in the backfill material. The interaction of steel corrosion products and radionuclides is part of such a scenario. In this work, the sorption of Th(IV) onto 2-line-ferrihydrite (FeOOH x H2O) and magnetite (Fe3O4), used as models for steel corrosion products, has been studied using EXAFS spectroscopy. Sorption samples were prepared in 0.1 M NaClO4 solutions at acidic pH (initial pH values in the range 3.0-4.2) either from undersaturation and supersaturation conditions with respect to amorphous ThO2. Two oxygen subshells, one at 2.37 A and another at 2.54 A, were observed in the first hydration sphere of Th in the case of the ferrihydrite samples. Th-Fe distances for the different ferrihydrite samples are approximately 3.60 A. These results indicate a corner sharing surface complex of Th(IV) ion onto the ferrihydrite surface where the Th atom shares one O atom with each of two coordinated octahedra. The longer Th-O distance accounts for coordinated water molecules. No significant changes in the structural environment of Th in terms of coordination numbers and distances were detected as a function of Th(IV) concentration. Magnetite samples sorbing Th(IV) also showed also a strong distortion of the O shell, but in contrast to ferrihydrite, two types of nearest Fe atoms were detected at 3.50 A and 3.70 A. These results indicate that Th(IV) ion sorbs onto the magnetite surface as bidentate-corner sharing arrangements to [FeO6] octahedra and [FeO4] tetrahedra.

EXAFS Investigation of U(VI), U(IV), and Th(IV) Sulfato Complexes in Aqueous Solution

The local structure of U(VI), U(IV), and Th(IV) sulfato complexes in aqueous solution was investigated by U-L(3) and Th-L(3) EXAFS spectroscopy for total sulfate concentrations 0.05 < or = [SO(4)(2-)] < or = 3 M and 1.0 < or = pH < or = 2.6. The sulfate coordination was derived from U-S and Th-S distances and coordination numbers. The spectroscopic results were combined with thermodynamic speciation and density functional theory (DFT) calculations. In equimolar [SO(4)(2-)]/[UO(2)(2+)] solution, a U-S distance of 3.57 +/- 0.02 Angstrom suggests monodentate coordination, in line with UO(2)SO(4)(aq) as the dominant species. With increasing [SO(4)(2-)]/[UO(2)(2+)] ratio, an additional U-S distance of 3.11 +/- 0.02 Angstrom appears, suggesting bidentate coordination in line with the predominance of the UO(2)(SO(4))(2)(2-) species. The sulfate coordination of Th(IV) and U(IV) was investigated at [SO(4)(2-)]/[M(IV)] ratios > or = 8. The Th(IV) sulfato complex comprises both, monodentate and bidentate coordination, with Th-S distances of 3.81 +/- 0.02 and 3.14 +/- 0.02 Angstrom, respectively. A similar coordination is obtained for U(IV) sulfato complexes at pH 1 with monodentate and bidentate U-S distances of 3.67 +/- 0.02 and 3.08 +/- 0.02 Angstrom, respectively. By increasing the pH value to 2, a U(IV) sulfate precipitates. This precipitate shows only a U-S distance of 3.67 +/- 0.02 Angstrom in line with a monodentate linkage between U(IV) and sulfate. Previous controversially discussed observations of either monodentate or bidentate sulfate coordination in aqueous solutions can now be explained by differences of the [SO(4)(2-)]/[M] ratio. At low [SO(4)(2-)]/[M] ratios, the monodentate coordination prevails, and bidentate coordination becomes important only at higher ratios.

Cyclic voltammetry modeling, geometries, and electronic properties for metallofullerene complexes withμ3-η2:η2:η2-C60bonding mode

Reduction potential (E(red)) values have been calculated and compared with available cyclic voltammetry (CV) data for 10 metallofullerene complexes with the mu(3)-eta(2):eta(2):eta(2)-C(60) (M(3)-C(6)[C(60)]) bonding mode. Consideration of bulk solvent effects is essential for the calculation of the E(red) values. Scaling factors for the electrostatic terms of the solvation energies have been introduced to fully describe the experimental cyclic voltammograms with a small mean deviation of 0.07 V. Multiple electron reductions induce movement of the metal cluster moieties on the C(60) surface, which is accompanied with the changes in some M-C[C(60)] bonds from pi-type to sigma-type mode. However, the changes in M(3)-C(60) distances, as well as the geometric changes of M(3) and C(60), are small for the reductions, which is in harmony with the high chemical and electrochemical stability of the metallofullerenes. Our population analyses reveal that the added electrons are not localized at the C(60) moieties, and electron population in the metal clusters is significant, more than 20% (av. 37%), for all the reductions. Furthermore, we demonstrated that the two close one-electron redox waves in CV diagrams are strongly correlated with significant electron delocalization, about 40-80%, to the metal-cluster moieties in these metallofullerene complexes.

Theoretical Study on Molecular Property of Protactinium(V) and Uranium(VI) Oxocations: Why Does Protactinium(V) Form Monooxo Cations in Aqueous Solution?

The stability of Pa(V) and U(VI) oxocations in aqueous solution were theoretically investigated by means of density functional theory calculations. As a result, the present calculations clearly supported an experimental result from an energetic point of view that monooxo protactinyl cation, PaO3+, is a preferable species for Pa(V) in aqueous solution, although dioxo protactinyl cation, PaO2+, is not a feasible form. By an analysis of molecular orbitals, we revealed that 6d orbitals of Pa(V) destabilize the pi orbitals of PaO2+, because 6d-2p antibonding orbital conflicts with another 5f-2p bonding orbital. For stable dioxo uranyl cation, UO2(2+), we found that 6d orbitals of U(VI), in contrast, form a bonding orbital with the 2p orbitals, and this bonding orbital coexists at an angle with the 5f-2p bonding orbital due to an electron correlation.

Speciation of the Curium(III) Ion in Aqueous Solution: A Combined Study by Quantum Chemistry and Molecular Dynamics Simulation

The structures of aquo complexes of the curium(III) ion have been systematically studied using quantum chemical and molecular dynamics (MD) methods. The first hydration shell of the Cm3+ ion has been calculated using density functional theory (DFT), with and without inclusion of the conductor-like polarizable continuum medium (CPCM) model of solvation. The calculated results indicate that the primary hydration number of Cm3+ is nine, with a Cm-O bond distance of 2.47-2.48 A. The calculated bond distances and the hydration number are in excellent agreement with available experimental data. The inclusion of a complete second hydration shell of Cm3+ has been investigated using both DFT and MD methods. The presence of the second hydration shell has significant effects on the primary coordination sphere, suggesting that the explicit inclusion of second-shell effects is important for understanding the nature of the first shell. The calculated results indicate that 21 water molecules can be coordinated in the second hydration shell of the Cm3+ ion. MD simulations within the hydrated-ion model suggest that the second-shell water molecules exchange with the bulk solvent with a lifetime of 161 ps.

Pyridine and phosphonate containing ligands for stable lanthanide complexation. An experimental and theoretical study to assess the solution structure

We report an experimental and theoretical study of the stability and solution structure of lanthanide complexes with two novel ligands containing pyridine units and phosphonate pendant arms on either ethane-1,2-diamine (L2) or cyclohexane-1,2-diamine (L3) backbones. Potentiometric studies have been carried out to determine the protonation constants of the ligands and the stability constants of the complexes with Gd(III) and the endogenous metal ions Zn(II) and Cu(II). While the stability constant of the GdL2 complex is too high to be determined by direct pH-potentiometric titrations, the cyclohexyl derivative GdL3 has a lower and assessable stability (log K(GdL3)=17.62). Due to the presence of the phosphonate groups, various protonated species can be detected up to pH approximately 8 for both ligands and all metal ions studied. The molecular clusters [Ln(L)(H2O)](3-).19H2O (Ln=La, Nd, Ho or Lu; L=L2 or L3) were characterized by theoretical calculations at the HF level. Our calculations provide two minimum energy geometries where the ligand adopts different conformations: twist-wrap (tw), in which the ligand wraps around the metal ion by twisting the pyridyl units relative to each other, and twist-fold (tf), where the slight twisting of the pyridyl units is accompanied by an overall folding of the two pyridine units towards one of the phosphonate groups. The relative free energies of the tw and tf conformations of [Ln(L)(H2O)]3- (L=L2, L3) complexes calculated in aqueous solution (C-PCM) by using the B3LYP model indicate that the tw form is the most stable one along the whole lanthanide series for the complexes of L3, while for those of L2 only the Gd(III) complex is more stable in the tf conformation by ca. 0.5 kcal mol-1. 1H NMR studies of the Eu(III) complex of L3 show the initial formation of the tf complex in aqueous solution, which slowly converts to the thermodynamically stable tw form. The structures calculated for the Nd(III) complexes are in reasonably good agreement with the experimental solution structures, as demonstrated by Nd(III)-induced relaxation rate enhancement effects in the 1H NMR spectra.

The highest water exchange rate ever measured for a Gd(III) chelate

The complex [Gd(L)(H2O)]3- (H(6)L =N,N'-bis(6-carboxy-2-pyridylmethyl)ethylenediamine-N,N'-methylenephosphonic acid) displays the highest water exchange rate ever measured for a Gd(III) chelate (k(298)(ex)= 8.8 x 10(8) s(-1)), which is attributed to the flexibility of the metal coordination environment.