Accurate Predictions of Volatile Plutonium Thermodynamic Properties

Abstract

The ability to predict the nature and amounts of plutonium emissions in industrial accidents, such as in solvent fires at PUREX nuclear reprocessing facilities, is a key concern of nuclear safety agencies. In accident conditions and in the presence of oxygen and water vapor, plutonium is expected to form the three major volatile species $\mathrm{Pu}{\mathrm{O}}_2$, $\mathrm{Pu}{\mathrm{O}}_3$, and $\mathrm{Pu}{\mathrm{O}}_2(\mathrm{OH}{)}_2$, for which the thermodynamic data necessary for predictions (enthalpies of formation and heat capacities) presently show either large uncertainties or are lacking. In this work we aim to alleviate such shortcomings by obtaining the aforementioned data via relativistic correlated electronic structure calculations employing a two-step multi-reference approach (MS-CASPT2 with SO-RASSI), which is able to describe the multireference character of the ground-state wave functions of $\mathrm{Pu}{\mathrm{O}}_3$ and $\mathrm{Pu}{\mathrm{O}}_2(\mathrm{OH}{)}_2$. We benchmark this approach by comparing it to relativistic coupled cluster calculations for the ground, ionized, and excited states of $\mathrm{Pu}{\mathrm{O}}_2$. Our results allow us to predict enthalpies of formation ${\Delta }_f{H}^{\ominus }(298.15\mathrm{K})$ of $\mathrm{Pu}{\mathrm{O}}_2$, $\mathrm{Pu}{\mathrm{O}}_3$ and $\mathrm{Pu}{\mathrm{O}}_2(\mathrm{OH}{)}_2$ to be $-449.5\pm 8.8$, $-553.2\pm 27.5$, and $-1012.6\pm 38.1\mathrm{kJ}\mathrm{mo}{\mathrm{l}}^{-1}$, respectively, which confirm the predominance of plutonium dioxide, but also reveal the existence of plutonium trioxide in the gaseous phase under oxidative conditions, though the partial pressures of $\mathrm{Pu}{\mathrm{O}}_3$ and $\mathrm{Pu}{\mathrm{O}}_2(\mathrm{OH}{)}_2$ are nonetheless always rather low under a wet atmosphere. Our calculations also permit us to reassess prior results for $\mathrm{Pu}{\mathrm{O}}_2$, establishing that the ground state of the $\mathrm{Pu}{\mathrm{O}}_2$ molecule is mainly of ${}^5{\mathrm{\Sigma }}_{\mathrm{g}}^{+}$ character, as well as to confirm the experimental value for the adiabatic ionization energy of $\mathrm{Pu}{\mathrm{O}}_2$.

arXiv comment17 pages, 5 figures, 9 tables and supplementary information

References in corpus (4)

• Comment on: "Estimating the Hartree-Fock limit from finite basis set calculations" [Jensen F (2005) Theor Chem Acc 113:267]
• Equation-of-Motion Coupled-Cluster Theory based on the 4-component Dirac-Coulomb(-Gaunt) Hamiltonian. Energies for single electron detachment, attachment and electronically excited states
• On the Multi-Reference Nature of Plutonium Oxides: PuO${}_2^{2+}$, PuO${}_2$, PuO${}_3$ and PuO${}_2$(OH)${}_2$
• Accurate Predictions of Volatile Plutonium Thermodynamic Properties

Cited by in corpus (7)

• The DIRAC code for relativistic molecular calculations
• Implementation of relativistic coupled cluster theory for massively parallel GPU-accelerated computing architectures
• Accurate Predictions of Volatile Plutonium Thermodynamic Properties
• Assessing MP2 frozen natural orbitals in relativistic correlated electronic structure calculations
• Formulation and Implementation of Frequency-Dependent Linear Response Properties with Relativistic Coupled Cluster Theory for GPU-accelerated Computer Architectures
• Electronic Spectra of Ytterbium Fluoride from Relativistic Electronic Structure Calculations
• Reassessing the potential of TlCl for laser cooling experiments via four-component correlated electronic structure calculations