Electronic structure of hydrated Ce4+ ions in solution: An x-ray-absorption study

Structures and Free Energies of Cerium Ions in Acidic Electrolytes

The Ce³⁺/Ce⁴⁺ redox potential changes with the electrolyte, which could be due to unequal anion complexation free energies between Ce³⁺ and Ce⁴⁺ or a change in the solvent electrostatic screening. Ce complexation with anions and solvent screening also affect the solubility of Ce and charge transfer kinetics for electrochemical reactions involving waste remediation and energy storage. We report the structures and free energies of cerium complexes in seven acidic electrolytes based on Extended X-ray Absorption Fine Structure, UV-vis, and Density Functional Theory calculations. Ce³⁺ coordinates with nine water molecules as [Ce(H₂O)₉]³⁺ in all studied electrolytes. However, Ce⁴⁺ complexes with anions in all electrolytes except HClO₄. Thus, our results suggest that Ce⁴⁺-anion complexation leads to the large shifts in standard redox potential. Long range screening effects are smaller than the anion complexation energies but could be responsible for changes in the Ce solubility with acid.

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.

Aqueous chemistry of Ce(iv): estimations using actinide analogues

The prediction of cerium (Ce) aqueous speciation is relevant in many research fields. Indeed, Ce compounds are used for many industrial applications, which may require the control of Ce aqueous chemistry for their synthesis. The aquatic geochemistry of Ce is also of interest. Due to its growing industrial use and its release into the environment, Ce is now considered as an emerging contaminant. Cerium is also used as a proxy of (paleo)redox conditions due to the Ce(iv)/Ce(iii) redox transition. Finally, Ce(iv) is often presented as a relevant analogue of tetravalent actinides (An(iv)). In the present study, quantum chemical calculations were conducted to highlight the similarities between the structures of Ce(iv) and tetravalent actinide (An(iv); An = Th, Pa, U, Np, Pu) aqua-ions, especially Pu(iv). The current knowledge of An(iv) hydrolysis, solubility and colloid formation in water was briefly reviewed but important discrepancies were observed in the available data for Ce(iv). Therefore, new estimations of the hydrolysis constants of Ce(iv) and the solubility of Ce(iv)-(hydr)oxides are proposed, by analogy with Pu(iv). By plotting pH-Eh (Pourbaix) diagrams, we showed that the pH values corresponding to the onset of Ce(iv) species formation (i.e. Ce(iv)-(hydr)oxide or dissolved Ce(iv)) agreed with various experimental results. Although further experimental studies are required to obtain a more accurate thermodynamic database, the present work might yet help to predict more accurately the Ce chemical behavior in aqueous solution.

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.

A QMCF-MD investigation of the structure and dynamics of Ce4+ in aqueous solution

A quantum-mechanical charge-field molecular dynamics simulation has been performed for a tetravalent Ce ion in aqueous solution. In this framework, the complete first and second hydration spheres are treated by ab initio quantum mechanics supplemented by an electrostatic embedding technique, making the construction of non-Coulombic solute-solvent potentials unnecessary. During the 10 ps of simulation time, the structural aspects of the solution were analyzed by various methods. Experimental results such as the mean Ce-O bond distance and the predicted first-shell coordination number were compared to the results obtained from the simulation resolving some ambiguities in the literature. The dynamics of the system were characterized by mean ligand residence times and frequency/force constant calculations. Furthermore, Ce-O and Ce-H angular radial distribution plots were employed, yielding deeper insight into the structural and dynamical aspects of the system.

Hydration of trivalent lanthanum revisited - An ab initio QMCF-MD approach

The previously investigated La³⁺-hydrate has been re-evaluated by means of the quantum mechanical charge field (QMCF) molecular dynamics (MD) approach. Improved description of the hydration characteristics has been realised by including the full second hydration shell into the quantum mechanically treated region and by introducing the influence of the surrounding bulk via an electrostatic embedding technique. Analytical tools such as the ligand angular radial distribution analysis have been employed to gain deeper insight into the structural features of the hydrate. La³⁺ simultaneously forms nona- and decahydrates with capped trigonal and quadratic prismatic structure, besides small amounts of an octahydrate.