Applications of polarizable continuum models to determine accurate solution-phase thermochemical values across a broad range of cation charge - the case of U(III-VI).

Contributing factors to the solution-phase correction to the free energy of the molecular clusters U(H2O)n(3+/4+) and UO2(H2O)m(1+/2+) (n = 8, 9, 30, 41, 77; m = 4, 5, 30, 41, 77) have been examined as a function of cavity type in the integrated-equation-formalism-protocol (IEF) and SMD polarizable continuum models (PCMs). It is observed that the free energy correction, Gcorr, does not smoothly converge to zero as the number of explicitly solvating water molecules approaches the bulk limit, and the convergence behavior varies significantly with cavity and model. The rates of convergence of the gas-phase hydration energy, ΔGhyd, wherein the bare metal ion is inserted into a molecular water cluster and ΔGcorr for the reaction exhibit wide variations as a function of ion charge, cavity, and model. This is the likely source of previously reported discrepancies in predicted free energies of solvation for metal ions when using different PCM cavities and/or models. The cancellation of errors in ΔGhyd and ΔGcorr is optimal for clusters consisting of only a second solvation shell of explicit water molecules (n = m = 30). The UFF cavity within IEF, in particular, exhibits the most consistent cancellation of errors when using a molecular cluster consisting of a second shell of solvating water for all oxidation states of uranium, leading to accurate free energies of solvation ΔGsolv for these species.

Multireference configuration interaction calculations of the first six ionization potentials of the uranium atom.

The first 6 ionization potentials (IPs) of the uranium atom have been calculated using multireference configuration interaction (MRCI+Q) with extrapolations to the complete basis set limit using new all-electron correlation consistent basis sets. The latter was carried out with the third-order Douglas-Kroll-Hess Hamiltonian. Correlation down through the 5s5p5d electrons has been taken into account, as well as contributions to the IPs due to the Lamb shift. Spin-orbit coupling contributions calculated at the 4-component Kramers restricted configuration interaction level, as well as the Gaunt term computed at the Dirac-Hartree-Fock level, were added to the best scalar relativistic results. The final ionization potentials are expected to be accurate to at least 5 kcal/mol (0.2 eV) and thus more reliable than the current experimental values of IP3 through IP6.

Solutes at the liquid:liquid phase boundary--Solubility and solvent conformational response alter interfacial microsolvation.

A detailed understanding of solvent structure and dynamics at liquid:liquid interfaces is a necessary precursor for control and manipulation of these phase boundaries. Experimentally, amphiphilic solutes are often used to alter transport properties across water:organic interfaces; however, a fundamental model for the mechanism of this action has not been determined. This work compares the solvation profiles of ampiphilic solutes that traverse the phase boundary in binary water:n-hexane, and the individual microsolvation processes for interfacial water and hexane molecules therein. Microsolvation is defined as the rare event where one solvent molecule temporarily penetrates the co-solvent phases and is fully solvated therein. The solutes tri-butyl phosphate (TBP), hydrogen di-butyl phosphate, and di-hydrogen mono-butyl phosphate have been examined as they exhibit a systematic increase in aqueous solubility and selectively partition to the interfacial region at the infinite dilution limit. The relationship between adopted configurations of the solute, orientation of the solvent, and the ability of the solute to enhance microsolvation, specifically the ability of n-hexane to penetrate the aqueous phase, is demonstrated within a 20 Å radius of TBP.

Static electric dipole polarizabilities of An(5+/6+) and AnO2 (+/2+) (An = U, Np, and Pu) ions.

The parallel components of static electric dipole polarizabilities have been calculated for the lowest lying spin-orbit states of the penta- and hexavalent oxidation states of the actinides (An) U, Np, and Pu, in both their atomic and molecular diyl ion forms (An(5+/6+) and AnO2 (+/2+)) using the numerical finite-field technique within a four-component relativistic framework. The four-component Dirac-Hartree-Fock method formed the reference for MP2 and CCSD(T) calculations, while multireference Fock space coupled-cluster (FSCC), intermediate Hamiltonian Fock space coupled-cluster (IH-FSCC) and Kramers restricted configuration interaction (KRCI) methods were used to incorporate additional electron correlation. It is observed that electron correlation has significant (∼5 a.u.(3)) impact upon the parallel component of the polarizabilities of the diyls. To the best of our knowledge, these quantities have not been previously reported and they can serve as reference values in the determination of various electronic and response properties (for example intermolecular forces, optical properties, etc.) relevant to the nuclear fuel cycle and material science applications. The highest quality numbers for the parallel components (αzz) of the polarizability for the lowest Ω levels corresponding to the ground electronic states are (in a.u.(3)) 44.15 and 41.17 for UO2 (+) and UO2 (2+), respectively, 45.64 and 41.42 for NpO2 (+) and NpO2 (2+), respectively, and 47.15 for the PuO2 (+) ion.

Static electric dipole polarizabilities of tri- and tetravalent U, Np, and Pu ions.

High-quality static electric dipole polarizabilities have been determined for the ground states of the hard-sphere cations of U, Np, and Pu in the III and IV oxidation states. The polarizabilities have been calculated using the numerical finite field technique in a four-component relativistic framework. Methods including Fock-space coupled cluster (FSCC) and Kramers-restricted configuration interaction (KRCI) have been performed in order to account for electron correlation effects. Comparisons between polarizabilities calculated using Dirac-Hartree-Fock (DHF), FSCC, and KRCI methods have been made using both triple- and quadruple-ζ basis sets for U(4+). In addition to the ground state, this study also reports the polarizability data for the first two excited states of U(3+/4+), Np(3+/4+), and Pu(3+/4+) ions at different levels of theory. The values reported in this work are the most accurate to date calculations for the dipole polarizabilities of the hard-sphere tri- and tetravalent actinide ions and may serve as reference values, aiding in the calculation of various electronic and response properties (for example, intermolecular forces, optical properties, etc.) relevant to the nuclear fuel cycle and material science applications.