Acidity of Uranyl(VI) Hydrate Studied with First-Principles Molecular Dynamics Simulations

Norm-Conserving 4f-in-Core Pseudopotentials and Basis Sets Optimized for Trivalent Lanthanides (Ln = Ce-Lu)

We present here a set of scalar-relativistic norm-conserving 4f-in-core pseudopotentials, together with complementary valence-shell Gaussian basis sets, for the lanthanide (Ln) series (Ce-Lu). The Goedecker, Teter, and Hutter (GTH) formalism is adopted with the generalized gradient approximation (GGA) for the exchange-correlation Perdew-Burke-Ernzerhof (PBE) functional. The 4f-in-core pseudopotentials are built through attributing 4f-subconfiguration 4fⁿ (n = 1-14) for Ln (Ln = Ce-Lu) into the atomic core region, making it possible to circumvent the difficulty of the description of the open 4fⁿ valence shell. A wide variety of computational benchmarks and tests have been carried out on lanthanide systems including Ln³⁺-containing molecular complexes, aqueous solutions, and bulk solids to validate the accuracy, reliability, and efficiency of the optimized 4f-in-core GTH pseudopotentials and basis sets. The 4f-in-core GTH pseudopotentials successfully replicate the main features of lanthanide structural chemistry and reaction energetics, particularly for nonredox reactions. The chemical bonding features and solvation shells, hydrolysis energetics, acidity constants, and solid-state properties of selected lanthanide systems are also discussed in detail by utilizing these new 4f-in-core GTH pseudopotentials. This work bridges the idea of keeping highly localized 4f electrons in the atomic core and efficient pseudopotential formalism of GTH, thus providing a highly efficient approach for studying lanthanide chemistry in multi-scale modeling of constituent-wise and structurally complicated systems, including electronic structures of the condensed phase and first-principles molecular dynamics simulations.

Ab Initio Molecular Dynamics Investigation of Beryllium Complexes

Structures of aqueous [Be(H₂O)₄]²⁺, its outer-sphere and inner-sphere complexes with F^-, Cl^-, and SO₄^2-, and dinuclear complexes with a [Be₂(κ-OH)(κ-SO₄)]⁺ core have been studied through Car-Parrinello molecular dynamics (CPMD) simulations with the BLYP functional. According to constrained CPMD/BLYP simulations and pointwise thermodynamic integration, the free energy of deprotonation of [Be(H₂O)₄]²⁺ and its binding free energy with F^- are 9.6 and -6.2 kcal/mol, respectively, in good accord with available experimental data. The computed activation barriers for replacing a water ligand in [Be(H₂O)₄]²⁺ with F^- and SO₄^2-, 10.9 and 13.6 kcal/mol, respectively, are also in good qualitative agreement with available experimental data. These ligand-substitution reactions are indicated to follow associative interchange mechanisms with backside (S_N2-like) attack of the anion relative to the aquo ligand it is displacing. Outperforming static density functional theory computations of the salient kinetic and thermodynamic quantities involving simple polarizable continuum solvent models, CPMD simulations are validated as a promising tool for studying the structures and speciation of beryllium complexes in aqueous solution.

Uranyl Arsenate Complexes in Aqueous Solution: Insights from First-Principles Molecular Dynamics Simulations

In this study, the structures and acidity constants (p K_a's) of uranyl arsenate complexes in solutions have been revealed by using the first principle molecular dynamics technique. The results show that uranyl and arsenate form stable complexes with the U/As ratios of 1:1 and 1:2, and the bidentate complexation between U and As is highly favored. Speciation-pH distributions are derived based on free energy and p K_a calculations, which indicate that for the 1:1 species, UO₂(H₂AsO₄)(H₂O)₃⁺ is the major species at pH < 7, while UO₂(HAsO₄)(H₂O)₃⁰ and UO₂(AsO₄)(H₂O)₃^- dominate in acid-to-alkaline and extreme alkaline pH ranges. For the 1:2 species, UO₂(H₂AsO₄)₂(H₂O)⁰ is dominant under acid-to-neutral pH conditions, while UO₂(HAsO₄)(H₂AsO₄)(H₂O)^-, UO₂(HAsO₄)(HAsO₄)(H₂O)^2-, and UO₂(AsO₄)(HAsO₄)(H₂O)^3- become the major forms in the pH range of 7.2-10.7, 10.7-12.1, and >12.1, respectively.

Uranyl Carbonate Complexes in Aqueous Solution and Their Ligand NMR Chemical Shifts and 17O Quadrupolar Relaxation Studied by ab Initio Molecular Dynamics

Dynamic structural effects, NMR ligand chemical shifts, and ¹⁷O NMR quadrupolar relaxation rates are investigated in the series of complexes UO₂²⁺, UO₂(CO₃)₃^4-, and (UO₂)₃(CO₃)₆^6-. Car-Parrinello molecular dynamics (CPMD) is used to simulate the dynamics of the complexes in water. NMR properties are computed on clusters extracted from the CPMD trajectories. In the UO₂²⁺ complex, coordination at the uranium center by water molecules causes a decrease of around 300 ppm for the uranyl ¹⁷O chemical shift. The final value of this chemical shift is within 40 ppm of the experimental range. The UO₂(CO₃)₃^4- and (UO₂)₃(CO₃)₆^6- complexes show a solvent dependence of the terminal carbonate ¹⁷O and ¹³C chemical shifts that is less pronounced than that for the uranyl oxygen atom. Corrections to the chemical shift from hybrid functionals and spin-orbit coupling improve the accuracy of chemical shifts if the sensitivity of the uranyl chemical shift to the uranyl bond length (estimated at 140 ppm per 0.1 Å from trajectory data) is taken into consideration. The experimentally reported trend in the two unique ¹³C chemical shifts is correctly reproduced for (UO₂)₃(CO₃)₆^6-. NMR relaxation rate data support large ¹⁷O peak widths, but remain below those noted in the experimental literature. Comparison of relaxation data for solvent-including versus solvent-free models suggest that carbonate ligand motion overshadows explicit solvent effects.

Predicting pKa in Implicit Solvents: Current Status and Future Directions

Computational prediction of condensed phase acidity is a topic of much interest in the field today. We introduce the methods available for predicting gas phase acidity and pKas in aqueous and non-aqueous solvents including high-level electronic structure methods, empirical linear free energy relationships (LFERs), implicit solvent methods, explicit solvent statistical free energy methods, and hybrid implicit–explicit approaches. The focus of this paper is on implicit solvent methods, and we review recent developments including new electronic structure methods, cluster-continuum schemes for calculating ionic solvation free energies, as well as address issues relating to the choice of proton solvation free energy to use with implicit solvation models, and whether thermodynamic cycles are necessary for the computation of pKas. A comparison of the scope and accuracy of implicit solvent methods with ab initio molecular dynamics free energy methods is also presented. The present status of the theory and future directions are outlined.

Coordination and hydrolysis of plutonium ions in aqueous solution using Car-Parrinello molecular dynamics free energy simulations

Car-Parrinello molecular dynamics (CPMD) simulations have been used to examine the hydration structures, coordination energetics, and the first hydrolysis constants of Pu(3+), Pu(4+), PuO2(+), and PuO2(2+) ions in aqueous solution at 300 K. The coordination numbers and structural properties of the first shell of these ions are in good agreement with available experimental estimates. The hexavalent PuO2(2+) species is coordinated to five aquo ligands while the pentavalent PuO2(+) complex is coordinated to four aquo ligands. The Pu(3+) and Pu(4+) ions are both coordinated to eight water molecules. The first hydrolysis constants obtained for Pu(3+) and PuO2(2+) are 6.65 and 5.70, respectively, all within 0.3 pH unit of the experimental values (6.90 and 5.50, respectively). The hydrolysis constant of Pu(4+), 0.17, disagrees with the value of -0.60 in the most recent update of the Nuclear Energy Agency Thermochemical Database (NEA-TDB) but supports recent experimental findings. The hydrolysis constant of PuO2(+), 9.51, supports the experimental results of Bennett et al. [Radiochim. Acta 1992, 56, 15]. A correlation between the pKa of the first hydrolysis reaction and the effective charge of the plutonium center was found.

Theoretical study of the microhydration of mononuclear and dinuclear uranium(VI) species derived from solvolysis of uranyl nitrate in water

The structures and energetics of mononuclear and dinuclear uranium species formed upon speciation of uranyl(VI) nitrate, UO(2)(NO(3))(2), in water are investigated by quantum chemistry using density functional theory and the wavefunction-based methods (MP2, CCSD, CCSD(T)). We provide a discussion of the basic coordination patterns of the various mono- and dinuclear uranyl compounds [(UO(2))(m)(X,Y)(2m-1)(H2O)(n)](+) (m = 1, 2; n = 0-4) found in a recent mass spectrometric study (Tsierkezos et al., Inorg Chem 2009, 48, 6287). The energetics of the complexation of the uranyl dication to the counterions OH(-) and NO(3) (-) as well as the degradation of the dinuclear species were studied by reference to a test set of 16 representative molecules with the MP2 method and the B3LYP, M06, M06-HF, and M06-2X DFT functionals. All DFT functionals provide structures and energetics close to MP2 results, with M06 family being slightly superior to the standard B3LYP functional.

Structure and stability of aquotetrafluorouranyl(VI) in the solid state — Density functional study of [UO2F4(H2O)][NMe4]2·2H2O,

Periodic density functional computations have been performed for solid [UO2F4(H2O)][NMe4]2·2H2O at the BLYP level. A model with disordered fluoro and aquo ligands in the uranyl anion is significantly lower in energy than one with a symmetrical assignment of these sites, which was favored in the original X-ray crystallography study. According to optimized energies and Car–Parrinello molecular dynamics (CPMD) simulations, the [UO2F4(H2O)]2− ion in the solid is stable with respect to loss of the coordinated water molecule. In contrast, CPMD simulations had found this ligand to be unbound in aqueous solution. The role of the counterions in stabilizing the higher coordination number in the crystal is discussed.

Binding of Pertechnetate to Uranyl(VI) in Aqueous Solution. A Density Functional Theory Molecular Dynamics Study

According to constrained Car-Parrinello molecular dynamics simulations and thermodynamic integration, the free binding energy between uranyl hydrate and pertechnetate in aqueous solution is significantly lower than that between uranyl and nitrate, namely, by 1.7 kcal mol(-1). This is the first study of the differential binding of these two ligands to uranyl, which can have implications for the separability of uranium and technetium during the reprocessing of nuclear waste.