INVESTIGATING THE SPECIATION AND EXTRACTION OF COMMERCIALLY
IMPORTANT METALS FROM SPENT NUCLEAR FUEL RAFFINATES: AN
INTEGRATED COMPUTATIONAL AND EXPERIMENTAL APPROACH
By
ALEX CHRISTOPHER SAMUELS
A dissertation submitted in partial fulfillment of the requirements for the degree of
DOCTOR OF PHILOSOPHY
WASHINGTON STATE UNIVERSITY
Program in Materials Science and Engineering
August 2014
© Copyright by ALEX CHRISTOPHER SAMUELS, 2014 All Rights Reserved
© Copyright by ALEX CHRISTOPHER SAMUELS, 2014 All Rights Reserved
To the Faculty of Washington State University:
The members of the Committee appointed to examine the dissertation of ALEX CHRISTOPHER SAMUELS find it satisfactory and recommend that it be accepted.
______
Aurora E. Clark, Ph.D., Co-Chair
______
Nathalie A. Wall, Ph.D., Co-Chair
______
Kirk A. Peterson, Ph.D.
______
Kenneth L. Nash, Ph.D.
iii
Acknowledgments
This work could not have been completed without the support of my family, advisors, committee members, and colleagues. Without their help I would not have found my true passion and enjoyed my time at Washington State University.
I would like to thank my wife for being there for me during these past few years.
Without her I am not sure if I would have remembered to take a break from research and
writing to eat a meal. Most importantly you always believed in me and pushed me to be
my best.
My advisors were key in helping me to find my true passion in chemistry. I struggled my first few years trying to find exactly what kind of project I wanted to work on. Aurora Clark and Nathalie Wall, along with my committee, saw my strengths and helped me to my passion of using a combined theoretical and experimental approach to solving chemical problems. This change in research direction reinvigorated my passion
for science. Aurora and Nathalie always pushed my to do by best and would not accept
anything less. I know the knowledge and skills they have passed on to me will help to me
grow in my career. Most importantly, Aurora and Nathalie made me at home at WSU and
like an important member or their respective groups.
My committee was always there to challenge me and make me better. I will never
forget what Dr. Nash told me during my defense, “If you want to make it as a chemist,
you need to be a badass. You need to have passion for what you do and be confident
about it.” That quote has followed me throughout the last few years. When I found my
passion I was able to gain a new confidence that allowed me to grow as a professional.
iv
Additionally, Dr. Peterson helped me to find my footing in an area of chemistry where I was the weakest: quantum mechanics. Dr. Peterson always welcomed my questions and did not turn me away even when I visited his office several times a day during my first few years.
My colleagues became more than my co-workers and more like family over my time at WSU. From Yasi Bross helping me with my research, to Aaron Johnson teaching me practical aspects of solvent extraction, I was well supported by those I worked with
v
INVESTIGATING THE SPECIATION AND EXTRACTION OF COMMERCIALLY
IMPORTANT METALS FROM SPENT NUCLEAR FUEL RAFFINATES: AN
INTEGRATED COMPUTATIONAL AND EXPERIMENTAL APPROACH
Abstract
by Alex Christopher Samuels, Ph.D. Washington State University August 2014
Co-Chairs: Aurora E. Clark and Nathalie A. Wall
Understanding the solution chemistry of fission products in solution is central to the operation of the nuclear fuel cycle and to the environmental remediation of any byproducts that
migrate away from waste disposal facilities. Actinides can be separated from fission products and can be reprocessed into mixed oxide fuel (MOX). The actinides also pose an hazard if
released into the environment. Additionally other fission products, such as the platinum group
metals, could potentially be extracted from spent fuel raffiniates and spent nuclear fuel (SNF)
provide a new domestic feedstock for these commercially important metals.
Actinides that leach from waste tanks can exist in various oxidation states in solution,
thereby complicating the chemistry. Computational studies can supplement experimental work in
this area and can help to create a holistic understanding of actinides in solution. A starting point
for computational studies is to gain fundamental insight into their static geometries and
electronic structure. In the past, several computational studies have examined different actinides in solution, but a systematic study of the actinide series, including various oxidation states, does
vi
not exist. In this study uranium, neptunium, and plutonium are examined in various oxidation
states.
The separation of rhodium (III) from other platinum group metals (PGM) continues to be relevant to modern separations chemistry, as natural deposits become depleted and SNF is being considered a potential feedstock. Rhodium, as well as other PGMs are produced by 235U fission,
and in a commercial light water reactor around 4 kg of PGMs can be produced per ton of waste.
This potential source could provide this much needed element used primarily in automotive
catalytic converters for years to come.
A combination of computational and experimental techniques have been applied to
understand the aqueous behavior of selected fission products. Using computational calculations
at the quantum mechanical level the hydration numbers of U, Pu, and Np in various oxidation states been determined. A combined approach was utilized to understand the speciation of
Rh(III) in acidic media along with solvent exaction of Rh(III) from nitric acid was developed.
vii
TABLE OF CONTENTS
ACKNOWLEDGEMENTS ...... iv ABSTRACT ...... vi-vii LIST OF TABLES ...... xi LIST OF FIGURES ...... xiv CHAPTER 1: INTRODUCTION ...... 1 CHAPTER 2: COMPUTATIONAL METHODS ...... 7 INTRODUCTION ...... 7 INTRODUCTION TO DENSITY FUNCTIONAL THEORY ...... 7 BASIS SET SELECTION ...... 9 GEOMETRY OPTIMIZATION AND FREQUENCY CALCULATIONS ...... 10 SOLUTION PHASE FREE ENERGY AND CONTINUUM MODELS ...... 14 SIMULATION OF UV-VIS SPECTRA ...... 17 REFERENCES ...... 18
CHAPTER 3: EXPERIMENTAL METHODS...... 20 INTRODUCTION ...... 20 EXTRACTION EQUILIBRIA AND FREE ENERGY ...... 20 SLOPE ANALYSES ...... 21 EXTRACTANT SELECTION ...... 21 LABORATORY MATERIALS AND METHODS ...... 23 CONTROLLING PH OF EXTRACTIONS ...... 24 ION COUPLED PLASMA – OPTICAL EMISSION SPECTROSCOPY (ICP-OES) DETECTION OF RH(III) ...... 24 REFERENCES ...... 25
CHAPTER 4: THERMODYNAMIC AND SPECTROSCOPIC ASSIGNMENT OF THE AQUEOUS SOLVATION ENVIRONMENTS OF TRI- TO HEXAVALENT U, NP, AND PU USING LARGE HYDRATED CLUSTERS CALCULATED BY DENSITY FUNCTIONAL THEORY ...... 26 INTRODUCTION ...... 26
viii
COMPUTATIONAL METHODS ...... 28 RESULTS AND DISCUSSION ...... 31 CONCLUSIONS ...... 55 ACKNOWLEDGEMENTS ...... 57 REFERENCES ...... 57 CHAPTER 5: An Integrated Computational and Experimental Protocol for Understanding Rh(III) Speciation in Hydrochloric and Nitric Acid Solutions: ...... 60 INTRODUCTION ...... 60 EXPERIMENTAL AND COMPUTATIONAL METHODS ...... 62 RESULTS AND DISCUSSION ...... 67 CONCLUSIONS ...... 85 SUPPLEMENTARY INFORMATION ...... 86 ACKNOWLEDGEMENTS ...... 89 REFERENCES ...... 90 CHAPTER 6: Rh(III) Extraction by Phosphinic Acids: A Combined Experimental and Computational Protocol: ...... 92 INTRODUCTION ...... 92 EXPERIMENTAL AND COMPUTATIONAL METHODS ...... 93 RESULTS AND DISCUSSION ...... 97 CONCLUSIONS ...... 106 CONTINUING WORK AND FUTURE DIRECTIONS ...... 107 SUPPLEMENTARY INFORMATION ...... 110 REFERENCES ...... 111 CHAPTER 7: Conclusions ...... 113
APPENDIX A: Applications of Polarizable Continuum Models to Determine Predictive Solution Phase Thermochemical Properties Across A Broad Range of Cation Charge – The Case of U(III-VI) : ...... 115 INTRODUCTION ...... 115 COMPUTATIONAL METHODOLOGY ...... 118 RESULTS AND DISCUSSION ...... 121 CONCLUSIONS ...... 136
ix
SUPPLEMENTARY INFORMATION ...... 137 ACKNOWLEDGEMENTS ...... 140 REFERENCES ...... 141
APPENDIX B: Modulation of Hydride Formation Energies in Transition Metal Doped Mg by Alteration of Spin State: ...... 144 INTRODUCTION ...... 144 COMPUTATIONAL METHODS ...... 147 RESULTS AND DISCUSSION ...... 148 SUPPLEMENTARY INFORMATION ...... 156 ACKNOWLEDGEMENTS ...... 163 REFERENCES ...... 163
x
List of Tables
Table 4.1. UB3LYP calculated ΔGhyd values (in kcal/mol) for tetravalent actinides using reaction r2 calculated in aqueous solutions using different dielectric continuum models. The free energy
for water addition, ΔGadd, using reaction r6 is also presented...... 39
4+ Table 4.2. Average UB3LYP
Table 4.3. UB3LYP calculated ΔGhyd values (in kcal/mol) for trivalent actinides using reactions r2 calculated in aqueous solutions using different cavity models. The free energy for water addition to the first solvation shell ΔGadd using reaction r6 is also presented...... 47
Table 4.4. UB3LYP calculated ΔGhyd values (in kcal/mol) for hexavalent actinides using reaction r4 calculated in aqueous solutions using different cavity models. The free energy for water
addition to the first solvation shell ΔGadd using reaction r8 is also presented...... 50
Table 4.5. Average UB3LYP equatorial, rAn-OH2, and axial metal-oxygen bond lengths, rAn=O, in 2+ + (Å) for AnO2(H2O)4,5(H2O)26,25 and AnO2(H2O)4,5(H2O)26,25 in comparison with experimental values (in parentheses). Average charges (q) are reported, determined by NPA for the metal center, actinyl oxo-atom, O-atoms associated with the inner-sphere (IS) and outer-sphere (OS) water molecules...... 53
Table 4.6. UB3LYP calculated ΔGhyd values (in kcal/mol) for pentavalent actinides using reaction r4 calculated in aqueous solutions using different cavity models. The free energy for
water addition to the first solvation shell ΔGadd using reaction r8 is also presented...... 55
Table 5.1. Concentrations of Rh (in M) used for each portion of the UV-Vis spectrum as a
function of concentration of added HCl or HNO3...... 63
Table 5.2. Calculated solution phase free energies of the successive nitrate addition reactions to 3+ the initial Rh(H2O)6 species (in kcal/mol), along with the TDDFT calculated λmax (in nm) of the 3-x Rh(NO3)x(H2O)y (x = 0 – 3; y = 6 – 2x) products between 180 and 800 nm...... 68
Table 5.3. Calculated excited state transitions and their corresponding oscillator strength...... 73
xi
Table 5.4. Literature UV-Vis λmax for individual Rh species and UV-Vis absorption maxima as a function of the HCl concentration of the Rh chloride salt solutions. The observed λmax values from Figure 5.3 are also presented ...... 75
Table 5.5. Solution phase B3LYP/cc-pVDZ predicted ΔGrxn (kcal/mol) of chloride complexation reactions in aqueous solution...... 81
Table 5.6. Calculated speciation of Rh(III) complexes in various concentrations of HCl with ranges that provide a NRMSD ≤ 5%. For HCl concentrations above 2 M a a fit with a NRMSD ≤ 5% was not found. Fits with minimized NRMSD are presented in Table S5.3 in Supplementary Information...... 82
Table S5.1. Average bond lengths (Å) of Rh(III) nitrate species ...... 86
Table S5.2. Molar absorptivity of λmax (202 nm) for Rh(NO3)3 in various concentrations of HNO3...... 86
Table S5.3. Calculated speciation of Rh(III) complexes in various concentrations of HCl with error expressed as NRMSD of the UV-Vis spectra fit to experiment...... 87
Table S5.4 Predicted speciation of Rh(III) species in various concentrations of HCl from CZE...... 87
Table 6.1. Calculated free energy of extraction using proposed equations in kJ/mol compared to experimentally determined free energy...... 104
Table E6.1. Calculated ΔGD (kcal/mol) for acac distribution between water and hexane. DZ and TZ refer to the aug-cc-pVDZ and aug-cc-pVTZ basis sets, respectively...... 108
Table A.1. Gas-phase ΔGhyd values (in kcal/mol) for actinide hydration (reactions r1 and r3), where RSC(1994)-C and RSC(1994)-U refer to the contracted and uncontracted forms of the RSC(1994) basis set, and so forth...... 123
Table A.2. Gas-phase ΔGadd values (in kcal/mol) for reactions r5 and r7, where RSC(1994)-C and RSC(1994)-U refer to the contracted and uncontracted forms of the RSC(1994) basis set, and so forth...... 125
Table A.3. UB3LYP calculated ΔGsolv values (in kcal/mol) for U(III-VI) ...... 136
Table SA1. Gas-phase UB3LYP calculated ΔErxn values (in kcal/mol) for actinide hydration for the second solvation shell, where RSC(1994)-C and RSC(1994)-U refer to the contracted and uncontracted forms of the RSC(1994) basis set, and so forth...... 138
xii
Table SA2. UB3LYP calculated ΔGsol values (in kcal/mol) for actinide hydration using IEF- PCM...... 140
Table SA3. UB3LYP calculated ΔGsol values (in kcal/mol) for water addition reactions for actinide ions using IEF-PCM...... 140
Table SB1: NPA charges for Mg-TM (TM = Ti, V, Fe) clusters calculated at the CCSD(T) optimized geometries, at the DFT/B3LYP level of theory. LS = low-spin, IS = intermediate-spin, HS = high-spin. Mg is always positive, while the TM is always negatively charged...... 156
Table SB2: Mg-TM (TM = Ti, V, Fe) bond lengths (in Å) calculated at the DFT (B3LYP), MP2, and CCSD(T) levels...... 157