Gabriela A Picayo1, Brian D Etz1, Shubham Vyas1, Mark P Jensen1,2. 1. Chemistry Department, Colorado School of Mines, 1012 14th St, Golden, Colorado 80401, United States. 2. Nuclear Science and Engineering Program, Colorado School of Mines, 920 15th St, Golden, Colorado 80401, United States.
Abstract
We have determined the identity of the complexes extracted into the ALSEP process solvent from solutions of nitric acid. The ALSEP process is a new solvent extraction separation designed to separate americium and curium from trivalent lanthanides in irradiated nuclear fuel. ALSEP employs a mixture of two extractants, 2-ethylhexyl phosphonic acid mono-2-ethylhexyl ester (HEH[EHP]) and N,N,N',N'-tetra(2-ethylhexyl)diglycolamide (TEHDGA) in n-dodecane, which makes it difficult to ascertain the nature of the extracted metal complexes. It is often asserted that the weak acid extractant HEH[EHP] does not participate in the extracted complex under ALSEP extraction conditions (2-4 M HNO3). However, the analysis of the Am extraction equilibria, Nd absorption spectra, and Eu fluorescence emission spectra of metal-loaded organic phases argues for the participation of HEH[EHP] in the extracted complex despite the high acidity of the aqueous phases. The extracted complex was determined to contain fully protonated molecules of HEH[EHP] with an overall stoichiometry of M(TEHDGA)2(HEH[EHP])2·3NO3. Computations also demonstrate that replacing one TEHDGA molecule with one (HEH[EHP])2 dimer is likely energetically favorable compared to Eu(TEHDGA)3·3NO3, whether the HEH[EHP] dimer is monodentate or bidentate.
We have determined the identity of the complexes extracted into theALSEP process solvent from solutions of nitric acid. TheALSEP process is a new solvent extraction separation designed to separate americium and curium from trivalent lanthanides in irradiated nuclear fuel. ALSEP employs a mixture of two extractants, 2-ethylhexyl phosphonic acid mono-2-ethylhexyl ester (HEH[EHP]) and N,N,N',N'-tetra(2-ethylhexyl)diglycolamide (TEHDGA) in n-dodecane, which makes it difficult to ascertain the nature of the extracted metal complexes. It is often asserted that theweak acid extractant HEH[EHP] does not participate in the extracted complex under ALSEP extraction conditions (2-4 M HNO3). However, the analysis of the Am extraction equilibria, Nd absorption spectra, and Eu fluorescence emission spectra of metal-loaded organic phases argues for the participation of HEH[EHP] in the extracted complex despite the high acidity of the aqueous phases. The extracted complex was determined to contain fully protonated molecules of HEH[EHP] with an overall stoichiometry of M(TEHDGA)2(HEH[EHP])2·3NO3. Computations also demonstrate that replacing one TEHDGA molecule with one (HEH[EHP])2 dimer is likely energetically favorable compared to Eu(TEHDGA)3·3NO3, whether theHEH[EHP] dimer is monodentate or bidentate.
Implementation of advanced
nuclear fuel cycles is critical to plans
for sustainable use of nuclear energy. Proposed closed nuclear fuel
cycles based on partitioning and transmutation will recycle an class="Chemical">uranium
and n class="Chemical">plutonium from used nuclear fuel and then separate and transmute
minor actinide elements to short-lived nuclides in fast reactors.
The separation and recycle of uranium and plutonium by theplutoniumuranium reduction extraction (PUREX) process is well-studied, but
efficient chemical separation of thetrivalent minor actinides, americium
and curium, from thelanthanide-rich matrix the PUREX process leaves
behind is difficult due to the physical and chemical similarities
of thetrivalent lanthanide (Ln) and actinide (An) ions.[1,2]
In the past 60 years, a number of solvent extraction systems
for
actinide–lanthanide separations have been proposed.[3−9] Recently, a three-step solvent extraction separation, theALSEP
(actinidelanthanide separation) process, has been proposed to simplify
isolation of americium and curium from thelanthanides. In ALSEP,
theAn(III) and Ln(III) cations are extracted together from aqueous
nitric acid solutions into ann-dodecane organic
phase containing two potential ligands, an acidic dialkylorganophosphorus
extractant, and a neutral diglycolamide extractant.[10,11] Although several ALSEP formulations have been considered, the most
effective formulation uses a mixture of the extractants 2-ethylhexylphosphonic acid mono-2-ethylhexyl ester (HEH[EHP]) and N,N,N′,N′-tetra(2-ethylhexyl)diglycolamide (TEHDGA) in n-dodecane (Figure ).[12] Following extraction of thetrivalent
lanthanide and actinide cations from ca. 3 M HNO3, the
organic phase is contacted with a scrub solution to adjust the acidity
and remove minor impurities. Then, theamericium and curium are selectively
stripped from the scrubbed organic phase by the addition of a polyaminocarboxylic
acid such as diethylenetriaminepentaacetic acid (DTPA) or N-(2-hydroxyethyl)ethylenediaminetriacetic acid (HEDTA)
to the aqueous phase. TheALSEP process greatly simplifies the separation
of trivalent actinides by combining the partitioning of trivalent
f-elements and theAn/Ln separation processes into a single separation
cycle, allowing direct use of PUREX raffinate solutions, exhibiting
fast extraction rates, and performing robustly under a broad range
of process conditions.[13,14]
Figure 1
Chemical structures for (A) 2-ethylhexyl
phosphonic acid mono-2-ethylhexyl
ester (HEH[EHP]) and (B) N,N,N’,N’-tetra(2-ethylhexyl)
diglycolamide (TEHDGA), the two extractants used in combination in
the ALSEP organic phase. Structures of additional extractants used
in the computational analyses of this study and in related liquid–liquid
extraction processes are presented in Figure S1.
Chemical structures for (A) 2-ethylhexylphosphonic acid mono-2-ethylhexyl
ester (HEH[EHP]) and (B) N,N,N’,N’-tetra(2-ethylhexyl)
diglycolamide (TEHDGA), the two extractants used in combination in
theALSEP organic phase. Structures of additional extractants used
in the computational analyses of this study and in related liquid–liquid
extraction processes are presented in Figure S1.Our research aims to develop a
kinetic model of the processes that
underlie americium and curium stripping in ALSEP, but probing the
kinetics of this complex multiextractant system requires anunderstanding
of themetal complexes present in the bulk phases at equilibrium.
This work begins that process by dissecting the speciation of the
organic phase metal complexes in theALSEP extraction step.Tan class="Chemical">he complexes formed in solutions containing only a n class="Chemical">diglycolamide
or an acidic dialkylorganophosphorus extractant have been studied
extensively,[12,15−29] but only a few published works have investigated the bulk-phase
chemical speciation of trivalent actinide or lanthanide complexes
in the presence of both diglycolamide and acidic dialkylorganophosphorus
extractants.[12,18,30−32] These studies suggest ternary metal–extractant
complexes are present in the organic phase under certain conditions,
but they either employ different extractants (for example TODGA rather
than TEHDGA or thephosphoric acid HDEHP rather than the phosphonic
acid HEH[EHP]), examine pH ranges that are not relevant to extracting
conditions, or simulate organic phase conditions using solutions that
are not in equilibrium with an aqueous phase. Unsurprisingly, these
studies do not entirely agree. Some hypothesized that a ternary metal–diglycolamide–organophosphorus
complex is present under all conditions[33,34] while others
reported that a mixed-extractant complex exists only under specific
conditions.[12,31,32]
an class="Gene">Understanding tn class="Chemical">he speciation of theALSEP process extraction
step
is important to the accurate interpretation of kinetic data and thus
to pinpointing theunderlying processes that lead to extraction and
ultimately limit theactinide stripping rate. To resolve the ambiguities
of the earlier reports, we explore the composition and stoichiometry
of theequilibrium f-element complexes extracted into TEHDGA/HEH[EHP]/n-dodecane from nitric acid solutions under conditions directly
relevant to theALSEP process. Theequilibrium organic phase species
are characterized through equilibrium analysis, UV–vis and
time-resolved laser-induced fluorescence spectroscopy (TRLFS), and
density functional theory (DFT) calculations to probe themetal–ligand
geometries and orientation and identify the thermodynamically stable
metal-containing species present at equilibrium in theALSEPmetal
extraction step.
Results
Optical Spectroscopy of
Loaded Organic Phases
The coordination
environment of the organic phase complexes formed by extraction into
theALSEP process solvent was probed by optical spectroscopy. Lanthanides
were loaded into organic phases for spectroscopic measurements from
aqueous phases containing 0.01 M Nd or Eu. Neodymium is a particularly
useful surrogate for Am because the ionic radii of Am(III) and Nd(III)
are nearly identical,[35] and thehypersensitive
optical transitions of Nd can be useful for identifying changes in
themetal coordination sphere.[36] Europium(III),
on the other hand, is isoelectronic with Am(III) (4f6 vs
5f6) and displays strong and sensitive fluorescence in
the visible region.[37] Extractions were
performed as closely as possible to ALSEP process relevant conditions
(3–4 M HNO3 and 0.05 M HEH[EHP]/0.75 M HEH[EHP]),
but the need to balance adequate extraction while avoiding third-phase
formation meant using acid concentrations lower than the ideal extraction
conditions for ALSEP in some measurements.Neodymium absorption
spectra for the 4I9/2 → 4G5/2, 2G7/2 transition between 560 and
620 nm, the 4I9/2 → 4S3/2, 4F7/2 transition between 720 and
770 nm, and the 4I9/2 → 4F5/2, 2H9/2 transition between 775 and
825 nm are summarized in Figure . The spectra represent Nd extracted by 0.75 M HEH[EHP]/n-dodecane, 0.1 M TEHDGA/n-dodecane, or
theALSEP organic phase, 0.05 M TEHDGA/0.75 M HEH[EHP]/n-dodecane. The extraction of Nd into theALSEP organic phase from
2 M HNO3 (not shown), 3 M HNO3 (not shown),
and 4 M HNO3 produced identical spectra, but these spectra
of theALSEP organic phase are different from those of the organic
phases containing only 0.1 M TEHDGA/n-dodecane or
0.75 M HEH[EHP]/n-dodecane. The most notable differences
between the Nd-TEHDGA and the Nd-ALSEP spectra are diminished absorbances
at 575 and 591.5 nm in theTEHDGA sample relative to theALSEP sample
and enhanced absorbance at 735 nm for theTEHDGA sample. Furthermore,
a general narrowing of the peaks occurs in theTEHDGA sample, particularly
those centered around 802 and 587 nm, which is accompanied by a slight
blue shift at 582 and 587 nm.
Figure 2
UV–vis spectra of the neodymium 4I9/2→ 4G5/2, 2G7/2 (560–620 nm), 4I9/2 → 4S3/2, 4F7/2, and 2H9/2, 4F5/2 (700–850
nm) transitions
after extraction into 0.75 M HEH[EHP]/n-dodecane
from 0.001 M HNO3/1 M NaNO3 (short dash), 0.1
M TEHDGA/n-dodecane from 3.5 M HNO3 (short
dot), and 0.75 M HEH[EHP]/0.05 M TEHDGA/n-dodecane
(ALSEP) from 4 M HNO3 (solid line).
UV–vis spectra of theneodymium 4I9/2→ 4G5/2, 2G7/2 (560–620 nm), 4I9/2 → 4S3/2, 4F7/2, and 2H9/2, 4F5/2 (700–850
nm) transitions
after extraction into 0.75 M HEH[EHP]/n-dodecane
from 0.001 M HNO3/1 M NaNO3 (short dash), 0.1
M TEHDGA/n-dodecane from 3.5 M HNO3 (short
dot), and 0.75 M HEH[EHP]/0.05 M TEHDGA/n-dodecane
(ALSEP) from 4 M HNO3 (solid line).While the 0.1 M TEHDGA and 0.05 M TEHDGA/0.75 M HEH[EHP] spectra
are similar, each organic phase composition investigated exhibits
a unique spectrum, indicating that different average Nd coordination
environments are encountered in each organic phase. The principal
component analysis of the spectra and target transformation using
the program SIXpack[38] indicate the Nd absorption
spectrum of the two-extractant TEHDGA/HEH[EHP] ALSEP system does not
arise from a linear combination of the Nd spectra observed for the
single-extractant TEHDGA/n-dodecane and HEH[EHP]/n-dodecane systems. Moreover, peaks at 570 and 606 nm are
absent in the spectrum of theALSEP system, indicating that Nd{H(EH[EHP])2}3 is not present. Although the presence of anequilibrium mixture of multiple Nd complexes in theALSEP system cannot
be excluded based on these absorption spectra, the results of the
principal component analysis imply the presence of at least one Nd
complex in theALSEP organic phase that is not found in either of
the two single-extractant systems.The TRLFS emission spectra
of Eu extracted in the three different
extraction systems (Figure ) more clearly demonstrate that themetal coordination environment
of Eu in theALSEP system is distinctly different from the complexes
in the solutions containing only TEHDGA or HEH[EHP]. Each solution
displays distinctly different fluorescence spectra from the other
two systems, and the spectrum of theALSEP system is obviously not
related to the spectra of either of the single extractant systems.
Figure 3
TRLFS
spectra of europium extracted from 0.001 M HNO3/1 M NaNO3 into 0.75 M HEH[EHP]/n-dodecane
(dashed line), from 3.5 M HNO3 into 0.1 M TEHDGA/n-dodecane (dotted line), and from 4 M HNO3 into
0.05 M TEHDGA/0.75 M HEH[EHP]/n-dodecane (solid line),
respectively.
TRLFS
spectra of europium extracted from 0.001 M HNO3/1 M NaNO3 into 0.75 M HEH[EHP]/n-dodecane
(dashed line), from 3.5 M HNO3 into 0.1 M TEHDGA/n-dodecane (dotted line), and from 4 M HNO3 into
0.05 M TEHDGA/0.75 M HEH[EHP]/n-dodecane (solid line),
respectively.The emission spectra of all three
complexes were collected over
a wavelength range spanning the 5D0 → 7F transitions (J = 0–4), and three key differences are apparent in the spectra.
First, a significant intensity for the 5D0 → 7F0 transition at 579.8 nm is only observed for
the Eu-ALSEP system, and only one peak is seen for this nondegenerate
transition. Second, while the 5D0 → 7F1 transitions at 593 nm are similar for both theALSEP and TEHDGA/n-dodecane systems, as expected
for this magnetic dipole transition, the intensities of thehypersensitive 5D0 → 7F2 manifold
(618 nm) are substantially different for each of the solutions, with
the Eu-ALSEP system displaying a more intense 5D0 → 7F2 emission than the Eu-TEHDGA/n-dodecane system. Third, the band shapes and barycenters
of the 5D0 → 7F4 transition at ca. 700 nm are distinct for each organic phase composition.
Together the spectral differences in the TRLFS data indicate that
the Eu inner-sphere coordination is different in each extraction system
and that theALSEP organic phase contains a single Eu coordination
environment that does not match the Eu coordination in the organic
phases containing only TEHDGA or HEH[EHP].
Probing Water Coordination
by TRLFS
Water molecules
have been proposed as inner-sphere ligands in the organic phase complexes
of Eu in a mixed malonamide/HDEHP extraction system.[39] The possible presence of water in the inner coordination
sphere of the extracted complexes in theALSEP process was investigated
by measuring the lifetimes of Eu3+ samples prepared in
light water and identical heavy water samples made with deuterated
reagents. The presence of OH oscillators directly coordinated with
europium decreases the emission decay lifetime (τ = 0.11 ms,
aqueous Eu solution) significantly when compared to heavier OD oscillators
(τ = 3.3 ms).[40] The fluorescence
lifetime of the major fluorescence peaks in each solution could be
fit to a single exponential decay with χ2 < 1
and fluorescence lifetimes on the order of milliseconds (Table ). As expected, Eu3+ in the aqueous stock solutions was found to contain nine
water molecules directly coordinated in the inner coordination sphere.[41] In contrast, all the organic-phase Eu–ligand
complexes we studied exhibited minimal inner-sphere water molecules
as shown in Table .
Table 1
Fluorescence Lifetimes in H2O and D2O for Eu Complexesa
complex
τ
in H2O (ms)
τ in D2O (ms)
number of H2O
aqueous Eu solution
0.11
3.30
9.20
0.05 M TEHDGA/0.75 M HEH[EHP]
2.17
2.26
0.02
0.1 M TEHDGA
2.10
2.35
0.05
0.75 M HEH[EHP]
3.03
3.27
0.03
Samples prepared by extracting
0.01 M Eu from 0.001 M HNO3/1 M NaNO3, 3.5 M
HNO3, and 4 M HNO3 into 0.75 M HEH[EHP]/n-dodecane, 0.1 M TEHDGA/n-dodecane, and
0.05 M TEHDGA/0.75 M HEH[EHP]/n-dodecane, respectively.
Number of inner-sphere water molecules calculated based on eq , with an absolute uncertainty
of ±0.5 H2O.
Samples prepared by extracting
0.01 M Eu from 0.001 M HNO3/1 M NaNO3, 3.5 M
HNO3, and 4 M HNO3 into 0.75 M HEH[EHP]/n-dodecane, 0.1 M TEHDGA/n-dodecane, and
0.05 M TEHDGA/0.75 M HEH[EHP]/n-dodecane, respectively.
Number of inner-sphere water molecules calculated based on eq , with an absolute uncertainty
of ±0.5 H2O.
Liquid–Liquid Extraction Equilibria
Ananalysis
of the partitioning behavior of Am3+ or Nd3+ between aqueous solutions of nitric acid and organic phases consisting
of TEHDGA, HEH[EHP], and n-dodecane was used to define
the stoichiometries of the organic-phase complexes in theALSEP system.
TEHDGA Dependence
Anamericium extraction from either
2 or 4 M HNO3 into organic phases containing 0.75 M HEH[EHP]
and varying amounts of TEHDGA in n-dodecane was examined.
For comparison, additional measurements of the Am extraction from
2 M HNO3 into n-dodecane containing only
variable concentrations of TEHDGA were also evaluated. The extraction
of Am from 4 M HNO3 into TEHDGA/n-dodecane
in the absence of HEH[EHP] could not be studied due to the formation
of a third phase under these conditions.[42]The results of the three TEHDGA dependence experiments are
summarized in Figure and Table . One
set of extractions was conducted in the absence of HEH[EHP], while
two sets were conducted in the presence of 0.75 M HEH[EHP]. The addition
of 0.75 M HEH[EHP] to theTEHDGA solution caused a substantial increase
in Am extraction, and the highest DAm values
were always observed in the 4 M HNO3–TEHDGA/0.75
M HEH[EHP]/n-dodecane system. In the solutions containing
only TEHDGA, the slope analysis of the Am extraction from 2 M HNO3 indicates a third power dependence on theTEHDGA concentration
(slope = 2.95 ± 0.09, Table ), implying the formation of a 1:3 Am:TEHDGA complex
in the organic phase, as previously reported for the extraction of
Eu by TEHDGA[20] and An3+/Ln3+ extraction by TODGA.[19,28] TheALSEP organic phases,
on the other hand, contain 0.75 M HEH[EHP] in addition to TEHDGA.
In theALSEP organic phases, the slopes of the extraction curves decrease
substantially from theTEHDGA-only case, giving second-power dependence
on theTEHDGA concentration at both acidities (2.11 ± 0.07 for
2 M HNO3 and 2.1 ± 0.1 for 4 M HNO3, Table ). This implies an
average 1:2 Am:TEHDGA stoichiometry for the extracted complexes at
both acidities when 0.75 M HEH[EHP] is also present in the organic
phase. Consequently, one TEHDGA molecule is displaced from theequilibrated
extracted complex when 0.75 M HEH[EHP] is added to theTEHDGA/n-dodecane organic phase to create theALSEP extraction
system.
Figure 4
TEHDGA dependence for the extraction of Am3+ from 4
M HNO3 into 0.03–0.075 M TEHDGA/0.75 M HEH[EHP]/n-dodecane (filled square), from 2 M HNO3 into
0.02–0.07 M TEHDGA/0.75 M HEH[EHP]/n-dodecane
(empty square), and from 2 M HNO3 into 0.02–0.2
M TEHDGA/n-dodecane (filled diamond).
Table 2
Stoichiometric Coefficients Determined
by Linear Regression Analysis of Logarithmic Am Extraction Data
aqueous phase
organic phase [TEHDGA]
organic phase [HEH[EHP]]
slope
intercept
2 M HNO3
0.02–0.2 M
0 M
2.95 ± 0.09
3.5 ± 0.1
2 M HNO3
0.02–0.07 M
0.75 M
2.11 ± 0.07
3.6 ± 0.1
4 M HNO3
0.03–0.075 M
0.75 M
2.1 ± 0.1
4.4 ± 0.1
2 M HNO3
0.05 M
0.005–0.75 M
0.43 ± 0.02
0.80 ± 0.02
1–5 M HNO3
0.05 M
0.75 M
3.10 ± 0.10a
4.35 ± 0.05
Slope = 2n –
3 with n = 3.05 ± 0.05 (see eq S15).
TEHDGA dependence for the extraction of Am3+ from 4
M HNO3 into 0.03–0.075 M TEHDGA/0.75 M HEH[EHP]/n-dodecane (filled square), from 2 M HNO3 into
0.02–0.07 M TEHDGA/0.75 M HEH[EHP]/n-dodecane
(empty square), and from 2 M HNO3 into 0.02–0.2
M TEHDGA/n-dodecane (filled diamond).Slope = 2n –
3 with n = 3.05 ± 0.05 (see an class="Chemical">eq S15).
HEH[EHP] Dependence
Attempts to study Am extraction
from 4 M HNO3 by 0.05 M TEHDGA/n-dodecane
with variable HEH[EHP] concentrations resulted in the formation of
a third phase at lower HEH[EHP] concentrations, and these investigations
were not pursued. However, third phases were not observed during contact
with 2 M HNO3, yielding distribution data across theHEH[EHP]
concentration range of 0–0.75 M (Figure ). The strong propensity of HEH[EHP] to dimerize
in n-dodecane means that the dimer (HEH[EHP])2 is actually the active form of the extractant in our experiments.[20] TheHEH[EHP] dependence of the Am extraction derived from linear regression
analysis was 0.43 ± 0.02 (Table ), similar to a previous report.[31]
Figure 5
HEH[EHP] dependence for the extraction of 241Am3+ from 2 M HNO3 into 0.05 M TEHDGA/0–0.75
M HEH[EHP]/n-dodecane. The 0 M HEH[EHP] data is represented
by a dashed line.
HEH[EHP] dependence for the extraction of 241Am3+ from 2 M HNO3 into 0.05 M TEHDGA/0–0.75
M HEH[EHP]/n-dodecane. The 0 M HEH[EHP] data is represented
by a dashed line.The simplest interpretations
of this result suggest an approximate
2:1 Am:(HEH[EHP])2 stoichiometry in the extracted complex.
This condition could be met either by the disruption of anHEH[EHP]
dimer, (HEH[EHP])2,[21] accompanied
by the incorporation of one HEH[EHP] molecule into a mononuclear extracted
complex, or by the formation of a polynuclear organic phase complex
with a 2:1 Am:(HEH[EHP])2 stoichiometry. Neither explanation
is satisfying, however. First, the ease of HEH[EHP] monomer formation
will decrease with increasing HEH[EHP] concentration and would cause
noticeable curvature in the distribution data over the wide range
of HEH[EHP] concentrations studied. Second, extraction experiments
at higher Am concentrations demonstrate that the distribution ratio
is independent of [Am] across the range of Am concentrations studied.
This observation is not consistent with the presence of a 2:1 Am:(HEH[EHP])2 complex in the organic phase, as DAm is independent of themetal concentration only if each extracted
complex contains a single Am cation.A reasonable alternate
explanation for the nonintegral HEH[EHP]
dependence is that HEH[EHP] acts as a phase modifier as well as an
extractant. The addition of HEH[EHP] substantially alters the polarity
of theTEHDGA/n-dodecane organic phases, as suggested
by its ability to suppress third-phase formation in theALSEP process.[10,11] This, in turn, will affect the extractability of the organic-phase
complex with the effect that theequilibrium constant for the extraction
will vary with the organic phase HEH[EHP] concentration, a condition
that is incompatible with using the distribution ratio to probe complexation
equilibria.In light of this complication, a different approach
for probing
theHEH[EHP] stoichiometry of the extracted complexes was devised.
The speciation of the extracted complex was studied by UV–vis
spectrophotometry in organic phases consisting of 0.05 M TEHDGA in n-dodecane and 0–0.075 M HEH[EHP]. The organic phases
were contacted with aqueous phases containing 0.01 M Nd in 2 M HNO3 at 35 °C. The higher temperature was used to discourage
third phase formation at these metal concentrations. After measuring
the amount of Nd extracted, the resulting organic phases were analyzed
by spectrophotometry between 550 and 620 nm to determine changes in
the Nd coordination environment at different HEH[EHP] concentrations.
The series of organic phase spectra obtained from these experiments
are displayed in Figure .
Figure 6
(a) UV–vis spectra of neodymium extraction at 35 °C
from 2 M HNO3 into 0–0.075 M HEH[EHP]/0.05 M TEHDGA/n-dodecane. (b) Speciation of extracted Nd calculated from eq , with m = 0.9 ± 0.1.
(a) UV–vis spectra of neodymium extraction at 35 °C
from 2 M HNO3 into 0–0.075 M HEH[EHP]/0.05 M TEHDGA/n-dodecane. (b) Speciation of extracted Nd calculated from eq , with m = 0.9 ± 0.1.Principal component analysis
conducted using both tan class="Chemical">he programs
SIXpack[38] and GlobalWorks (Online Instrument
Systems, Inc.) indicated that tn class="Chemical">he set of spectra for all 10 HEH[EHP]
concentrations studied was composed of two unique spectral components
(Figure S2). The spectrum of the first
species determined from the model-linked singular value decomposition
in GlobalWorks matched the spectrum of Nd(TEHDGA)3·3NO3 extracted from 2 M HNO3. The spectrum of the second
species matched that of the Nd-loaded ALSEP organic phase depicted
in Figure . Using
the spectra of these two Nd complexes, the ratio of the two complexes
in the organic phase could be calculated for each solution and used
to derive the number of HEH[EHP] dimers, m, coordinated
in theALSEP system for the following generalized equilibrium for
trivalent metal, M3+.where
species in the organic phase are indicated
by overbars, h represents the number of acidic hydrogens
in each metal-complexed HEH[EHP] dimer (h = 0, 1,
or 2), and n represents the number of nitrate ions
in the product complex, with n = 3 – m(2 – h). Defining theequilibrium
constant of this reaction to be Kexchange, theequilibrium constant expression can be written in the following
form:with b = log Kexchange + m(2 – h) log [HNO3]. Between 0.005 and 0.050 M HEH[EHP], an analysis
of the spectroscopic data with eq yields a value of m = 0.9 ±
0.1 (Figure ), implying
that under these conditions one HEH[EHP] dimer displaces one TEHDGA
molecule from M(TEHDGA)3·3NO3 to form the
complex M(TEHDGA)2(H(EH[EHP])2)·nNO3 in theALSEP organic
phase.
Nitric Acid Dependence
Tan class="Chemical">he stoichiometry of n class="Chemical">nitrate
and hydrogen ions in the extracted complex was probed by examining
Am3+ extraction into the standard ALSEP organic phase (0.05
M TEHDGA/0.75 M HEH[EHP]/n-dodecane) as a function
of theequilibrium aqueous concentration of nitric acid (Figure ). As previously
described, a minimum in the Am extraction is observed at approximately
0.5 M HNO3 due to a change in the organic-phase metal speciation
from complexes containing TEHDGA and nitrate at high acidity to homoleptic
complexes of monodeprotonated HEH[EHP] dimers, which show an inverse
dependence of themetal distribution ratio on the aqueous acidity
at low acidities.[11] However, even at aqueous
acidities higher than 1 M HNO3, direct slope analysis of
the concentration dependence of the distribution ratio is not a reliable
indication of the proton or nitrate stoichiometry of themetal complex
extracted in theALSEP process. The large change in ionic strength
across the range of aqueous nitric acid concentrations studied causes
substantial changes in the activity coefficients of the aqueous solutes,
the activity of water, the fraction of aqueous Am present as nitrate
complexes, the degree of nitric acid dissociation, the amount of nitric
acid extracted into the organic phase, and the fraction of each extractant
available to interact with Am. Consequently, a multiequilibrium thermodynamic
model for Am extraction from nitric acid into theALSEP organic phase
was developed to interpret the effect of nitric acid on Am extraction
and thenitrate and hydrogen ion stoichiometry of the extracted complex
using the following general extraction equilibrium:where h = n – 1 and h = 0, 1, or 2. Details of the model
and the relationship between the slope of thenitrate dependence and n, the number of nitrate ions incorporated into the extracted
complex, are described in the Supporting Information. The evaluation of n, the number of nitrate anions
extracted with each Am3+ cation, and Kex, theequilibrium constant for the extraction reaction,
between 1 and 5 M HNO3 gave n = 3.05 ±
0.06 and log Kex = 4.35 ± 0.05 (Figure ).
Figure 7
Effect of aqueous nitric
acid on Am extraction into 0.05 M TEHDGA/0.75
M HEH[EHP]/n-dodecane. (a) Experimental extraction
data (squares) fit to a third-degree polynomial (short dashed line).
The filled squares represent data used in the equilibrium activity
modeling shown in panel b. (b) Determination of the nitrate stoichiometry, n, in the extracted complex from the correlation of the
Am distribution ratio corrected for variations in aqueous activity
coefficients, aqueous nitrate complexation, and extraction of nitric
acid to the activity of nitrate ions in the aqueous phase using the
model described in the Supporting Information.
Effect of aqueous nitric
acid on Am extraction into 0.05 M TEHDGA/0.75
M HEH[EHP]/n-dodecane. (a) Experimental extraction
data (squares) fit to a third-degree polynomial (short dashed line).
The filled squares represent data used in theequilibrium activity
modeling shown in panel b. (b) Determination of thenitrate stoichiometry, n, in the extracted complex from the correlation of the
Am distribution ratio corrected for variations in aqueous activity
coefficients, aqueous nitrate complexation, and extraction of nitric
acid to the activity of nitrate ions in the aqueous phase using the
model described in the Supporting Information.
Computational Data
Further insight into tn class="Chemical">he inner-sphere
coordination of trivalent f-element cations with TEHDGA, HEH[EHP],
nitrate anions, and water molecules was gained by comparing the complexation
energies of europium–ligand complexes potentially present in
theALSEP organic phase. Many different coordination complexes are
possible in the organic phase given the flexibility of theAn3+ and Ln3+ coordination spheres and the presence
of multiple potential ligands. However, the experimental data rule
out many possibilities for the complexes extracted from 2–4
M HNO3, such as complexes containing only singly deprotonated
HEH[EHP] dimers (e.g., M{H(EH[EHP])2}3, complexes
with inner-sphere water molecules, or complexes containing a single
TEHDGA extractant. Consequently, we calculated energies for the following
complexes and explored the conformational space of each: TEHDGA only
(Eu(TEHDGA)3·3NO3), 2:1 ratio of TEHDGA:deprotonated
HEH[EHP] dimer (Eu(TEHDGA)2(H(EH[EHP])2)·2NO3), and 2:1 ratio of TEHDGA:protonated HEH[EHP] dimer (Eu(TEHDGA)2(HEH[EHP])2·3NO3). To enable tractable
calculations on these complexes, theethylhexyl substituents were
replaced with ethyl groups, with N,N,N′,N′-tetraethyldiglycolamide
(TEDGA) replacing TEHDGA and ethyl phosphonic acid monoethyl ester
(HE[EP]) replacing HEH[EHP] for the calculations. The structures of
TEDGA and HE[EP] are summarized in Figure S1.
Tan class="Chemical">he minimum energy geometries calculated for these three
stoichiometries of theTEDGA and HE[EP] complexes are shown in Figure . In the Eu(TEDGA)3·3NO3 complex, theeuropium ion is coordinated
only by tridentate TEDGA ligands, yielding a complex with D3 symmetry, a coordination number of 9, and
a calculated minimum energy of −9.7 kcal/mol. The three nitrate
anions are positioned in clefts between the alkyl-chains of adjacent
TEDGA ligands to balance the charge of europium and produce a neutral
complex, as found by Brigham et al. for the closely related extractant
TODGA.[19] Replacing one of theTEDGA ligands
with a singly deprotonated dialkylphoshonic acid dimer, H(E[EP])2–, yields Eu(TEDGA)2(H(E[EP])2)·2NO3. This complex was observed to have
two possible minimum energy geometries within 0.2 kcal/mol (−11.3
and −11.1 kcal/mol). The two geometries differ primarily in
the location of thenitrate anions: geometry A contains both nitrate
anions in the outer sphere creating an 8-coordinate Eu complex with C2 site symmetry, while geometry B contains one
monodentate nitrate anion in the Eu inner coordination sphere, giving
europium a coordination number of 9 and C1 symmetry. In geometry A, neither nitrate remains in the cleft between
theTEDGA ligands; instead, each nitrate is associated with one of
theTEDGA ligands (Figure ). In geometry B, the two nitrates are found in spaces between
theTEDGA molecules and theHE[EP] dimer. The movement of thenitrate
counteranions can be explained by the increase in inner-sphere free
space when one tridentate TEDGA is replaced by a bidentate H(E[EP])2– anion. The extra room allows the remaining
two TEDGA ligands to shift, removing the original clefts where thenitrate ions were positioned. Consequently, thenitrate ions reposition,
either outer sphere or inner sphere, to stabilize the newly formed
Eu complex. The energetic similarity of these complexes is unusual
because 9-coordinate Eu complexes are generally expected to be energetically
more stable than the 8-coordinate europium complex. However, an analysis
of theeuropium–ligand coordination distances reveals that
inner-sphere coordination of thenitrate anion in geometry B causes
the Eu–O bond lengths of theTEDGA ligands to increase (Table S2). As a result, the energies of the 8-
and 9-coordinate complexes differ by only 0.2 kcal/mol, which suggests
that if Eu(TEHDGA)2(H(EH[EHP])2)·2NO3 were to form in the organic phase, both complexes A and B
could coexist in solution.
Figure 8
Optimized Eu complexes, Eu(TEDGA)3·3NO3, Eu(TEDGA)2(H(E[EP])2)·2NO3 (A) and (B), and
Eu(TEDGA)2(HE[EP])2·3NO3 (C)
and (D), investigated investigated as possible species existing during
ALSEP extraction. Geometries A and B refer to two different Eu complexes
containing deprotonated HE[EP] dimers and geometries C and D refer
to two different Eu complexes containing protonated HE[EP] dimers.
See the text for further details of the complexes. Europium is teal,
phosphorus atoms are orange, oxygen atoms are red, nitrogen atoms
are blue, carbon atoms are gray, and hydrogen atoms are white. Dashed
lines from Eu highlight the chelating oxygen in every complex. Hydrogen
atoms outlined with black circles correspond to the acidic HE[EP]
hydrogen. Black dashed lines represent hydrogen bonding while solid
black lines represent bonds.
Optimized Eu complexes, Eu(TEDGA)3·3NO3, Eu(TEDGA)2(H(E[EP])2)·2NO3 (A) and (B), and
Eu(TEDGA)2(HE[EP])2·3NO3 (C)
and (D), investigated investigated as possible species existing during
ALSEP extraction. Geometries A and B refer to two different Eu complexes
containing deprotonated HE[EP] dimers and geometries C and D refer
to two different Eu complexes containing protonated HE[EP] dimers.
See the text for further details of the complexes. Europium is teal,
phosphorus atoms are orange, oxygen atoms are red, nitrogen atoms
are blue, carbon atoms are gray, and hydrogen atoms are white. Dashed
lines from Eu highlight the chelating oxygen in every complex. Hydrogen
atoms outlined with black circles correspond to the acidic HE[EP]
hydrogen. Black dashed lines represent hydrogen bonding while solid
black lines represent bonds.Five possible minima were located for the Eu(TEDGA)2(HE[EP])2·3NO3 complex, which forms from
a neutral, fully protonated (HE[EP])2 dimer (i.e., H2(E[EP])2) without the loss of a hydrogen ion to
the aqueous phase. Two possible complexes were close in energy (−25.7
and −22.9 kcal/mol), while the other three complexes (−12.8,
−11.1, and −5.9 kcal/mol) were substantially less stable.
Only the two most stable configurations, geometries C and D, are considered
here.In both of these complexes, one P=O--H–O
hydrogen
bond in the (HE[EP])2 dimer is disrupted. The critical
difference between the two lowest energy conformations of Eu(TEDGA)2(HE[EP])2·3NO3 is the interaction
between the (HE[EP])2 dimer and the Eu ion. In geometry
C of the Eu(TEDGA)2(HE[EP])2·3NO3 complex (Figure ), theeuropium sits in a C2 symmetry
site, all thenitrate ions remain in the outer coordination sphere,
a hydrogen bond forms between one HE[EP] and a nitrate anion, and
bidentate coordination between theHE[EP] dimer and Eu is observed.
The energy of this complex is calculated to be −22.9 kcal/mol.
The PO–H bond of theHE[EP] hydrogen bonding to nitrate elongates
from 1.03 to 1.57 Å, while the H--ON hydrogen bond between this
HE[EP] and thenitrate anion is 1.02 Å. This indicates that such
a complex would take on the character of a deprotonated HE[EP] dimer
hydrogen bonded to nitric acid. Alternately, the Eu(TEDGA)2(HE[EP])2·3NO3 complex can distort, leading
to the C1 symmetry complex shown in geometry
D. (Figure ) This
is the most stable Eu(TEDGA)2(HE[EP])2·3NO3 complex studied at −25.7 kcal/mol. In geometry D,
the (HE[EP])2 dimer becomes monodentate, as one Eu–OP
bond and one P=O--H–O–P hydrogen bond are cleaved,
and nitrate anions reposition to stabilize the entire complex. One
nitrate anion moves into the Eu inner coordination sphere with a Eu–ON
bond length of 2.34 Å while a second nitrate anionhydrogen bonds
to thehydrogen on the noncoordinated HE[EP] in the outer coordination
sphere. The noncoordinated HE[EP] PO–H bond length remains
1.04 Å while the H--ON hydrogen bond between this HE[EP] and
thenitrate anion is 1.48 Å. This implies that complex D would
have the character of a fully protonated HE[EP] dimer hydrogen bonded
to a nitrate anion. This asymmetric complex featuring inner- and outer-sphere
nitrate ions is 2.8 kcal/mol more stable than the Eu(TEDGA)2(HE[EP])2·3NO3 complex with all three
nitrates in the outer sphere and 19.0 kcal/mol more stable than Eu(TEDGA)3·3NO3.
Discussion
Together,
theequilibrium partitioning experiments, optical spectroscopy,
and modeling clearly demonstrate the formation of a new extracted
complex when HEH[EHP] is added to solutions of TEHDGAunder metal
extraction conditions in acidic biphasic systems. Previous evidence
of mixed TEHDGA–HEH[EHP] complexes under ALSEP extraction conditions
has been inconsistent. The original work on theALSEP process postulated
that ternary metal–TEHDGA–HDEHP (HDEHP = bis(2-ethylhexyl)phosphoric
acid) complexes were likely involved in An3+ and Ln3+ extraction under acidic conditions but also concluded that
it was unclear if metal–TEHDGA–HEH[EHP] complexes form.[11] Subsequent EXAFS studies of Nd3+ extracted
into mixtures of TEHDGA and HEH[EHP] from 1 M HNO3 suggested
the possibility of a Nd–O–P scattering path consistent
with “relatively weak interactions between Nd–HEH[EHP]–T2EHGDA
when extracted from a mildly acidic aqueous environment”.[12] In contrast to both these studies, a later EXAFS
study of Eu–TODGA–HDEHP complexes concluded that only
1:3:3 Eu:TODGA:nitrate complexes are extracted into 0.05 M TODGA/0.75
M HDEHP/n-dodecane from 3 M HNO3.[32] On the other hand, from FTIR experiments, Rama
Swami et al. proposed simultaneous coordination of TEHDGA and HEH[EHP]
or HDEHP when Eu3+ is extracted into mixtures of these
extractants from 3 M HNO3.[18] The observation of multiextractant organic-phase complexes in theALSEP system is also similar to reports for An3+ and Ln3+ extraction by mixtures of N,N′-dimethyl-N,N′-dioctylhexylethoxymalonamide
and dialkylphosphoric acid extractants.[39,43,44] Our experiments and modeling confirm the formation
and identity of the mixed metal–TEHDGA–HEH[EHP] complex
formed in theALSEP organic phase under acidic extracting conditions.Tan class="Chemical">he presence of neutral and acidic coextractants at variable concentrations
and high aqueous acidities makes tn class="Chemical">he determination of the organic-phase
speciation in theALSEP process difficult. At low acidity, theALSEP
process’ acidic extractant, HEH[EHP], forms the M{H(EH[EHP])2}3 complex,[12] liberating
hydrogen ions to the aqueous phase according to the following equilibrium.Eq is clearly opposed by high acidities, and
the extraction
of An3+ or lanthanides lighter than terbium by HEH[EHP]
is negligible at aqueous acidities of 1 M or greater.[45] By itself, the neutral extractant TEHDGA forms complexes
in the organic phase according the following equilibrium.[19,20]If thenitrate anions are
supplied by nitric
acid, eq will be favored
by high nitric acid concentrations and negligible at the low acidities
favored by eq . At high
nitric acid concentrations, experiments also suggest coextraction
of nitric acid, depending on the conditions,[20] with the following equilibrium.Each of these
equilibria has a characteristic
dependence on the extractant, nitrate, and hydrogen ion activities
(Table ). However,
in the mixed HEH[EHP]/TEHDGAALSEP system, a distinctly different
stoichiometry is observed, 1 An3+/Ln3+:2.1 TEHDGA:0.9
(HEH[EHP])2:3.05 NO3–, with
no aqueous hydrogen ions as reactant or product and no inner-sphere-coordinated
water. Furthermore, the optical spectroscopy of the extracted Nd and
Eu complexes demonstrate the stoichiometry observed in theALSEP system
solvent does not result from mixtures of the products of eq and eq or eq because the distinct optical signatures of M{H(EH[EHP])2}3, the absorption peaks at 570 and 606 nm in the
Nd spectrum (Figure ) and the Eu fluorescence emission band at 611 nm (Figure ), are absent. This disappearance
of the homoleptic Nd–HEH[EHP] complex under ALSEP extracting
conditions is also consistent with experiments by Gullekson et al.,
who titrated TEHDGA·HNO3/n-dodecane
into isolated organic phases initially containing Nd{H(EH[EHP])2}3.[12]
Table 3
Expected Dependence of Distribution
Ratios on the Organic-Phase Concentrations of Extractants and the
Activity of Hydrogen and Nitrate Ions in Aqueous Nitric Acid for Proposed
Extraction Equilibria
slope
slope
stoichiometry
slopea
equilibrium
extracted
complex
[TEHDGA]
[(HEH[EHP])2]
{H+}
{NO3–}
{NO3–}
3
M(TEHDGA)2(Hh(EH[EHP])2)·nNO3
2
1
n – 3
n
2n – 3
4
M{H(EH[EHP])2}3
0
3
–3
–3
5
M(TEHDGA)3·3NO3
3
0
0
3
3
6
M(TEHDGA)3·3NO3·pHNO3
3
0
p
p + 3
2p + 3
Slope expected
in nitric acid solutions
when [H+] ≈ [NO3–]
(eqs S7c and S15).
Slope expected
in nitric acid solutions
when [H+] ≈ [NO3–]
(eqs S7c and S15).Tan class="Chemical">he change in tn class="Chemical">he metal ion site symmetry of the extracted
complex
between organic phases containing only TEHDGA and those containing
mixtures of TEHDGA and HEH[EHP] is also apparent in the changes in
the optical spectra (Figures and 6). The 5D0 → 7F0 transition of Eu3+, which occurs near 580 nm in the condensed phases, is particularly
diagnostic. This transition is only allowed for a subset of noncubic
space groups without an inversion center, C, C, and C.[46] In other symmetries, for example D or any space group with an inversion
center, this transition is forbidden and is generally very weak if
observed at all. The extraction of Eu3+ into 0.75 M HEH[EHP]/n-dodecane or 0.2 M TEHDGA/n-dodecane yielded
solutions with barely discernible 5D0 → 7F0 emissions (Figure and Figure S6), in agreement with previous reports.[21,47] In contrast,
the extraction of Eu3+ from 4 M HNO3 into 0−0.75
M HEH[EHP]/0.05 M TEHDGA/n-dodecane gives a solution
with an obvious 5D0 → 7F0 emission band at 579.7 nm. The absence of a significant 5D0 → 7F0 emission
for the organic-phase Eu complexes containing only HEH[EHP] or TEHDGA
is readily understood based on the Eu site symmetry. Eu{H(EH[EHP])2}3 is believed to possess approximately octahedral
symmetry with an inversion center.[21,29] The Eu(TEHDGA)3·3NO3 complex extracted into 0.05 M TEHDGA
will have D3 site symmetry,[19] similar to theeuropiumtrisoxydiacetato and
trisdipicolinato complexes, which also show exceedingly weak 5D0 → 7F0 emission.[48,49] However, replacing one TEHDGA molecule with one mono- or bidentate
(HEH[EHP])2 will necessarily yield a lower symmetry complex
(C2, C2, or C1) with an allowed 5D0 → 7F0 transition,
as we observe. The lowest energy conformations we found in our theoretical
calculations (geometries C and D) show these lowered symmetries.
The changes in the 5D0 → 7F4 emission band at approximately 700 nm are also consistent
with a decrease in symmetry of the Eu site in theALSEP system compared
to the organic phase containing only TEHDGA. TheEu3+7F4 state is split into 4 levels in crystal fields
of D3 symmetry (i.e., for Eu(TEHDGA)3·3NO3), while it splits into 7 levels for C2 symmetry and 9 levels for C2 or C1 symmetry.[50] The broadening of the 5D0 → 7F4 transition in the Eu-ALSEP emission
spectrum compared to the Eu-TEHDGA spectrum is accompanied by a transition
from 4 discernible bands in the Eu-TEHDGA spectrum to the presence
of 6 discernible shoulders on the main peak for the Eu-ALSEP spectrum
(Figure S6), suggesting that the 7F4 state splits into at least 7 discrete levels. This
change in site symmetry of the extracted complex is consistent with
direct coordination of HEH[EHP] in theALSEP system.With tan class="Chemical">he
n class="Chemical">metal:ligand stoichiometry established by equilibrium
partitioning and spectroscopic experiments, the protonation state
of HEH[EHP] in the extracted complex remains key to understanding
the nature of the extracted complex. The extracted complex must be
charge neutral. The charge on the extracted An3+ and Ln3+ cations can be balanced by coextracted nitrate anions (eq ), deprotonated HEH[EHP]
molecules (e.g., H(EH[EHP])2– in eq ), or a combination of
the two. In addition, HNO3 may be incorporated into the
outer sphere of the extracted complex (eq ). The aqueous acidity disfavors the formation
and complexation of theH(EH[EHP])2– anion,
as the pKa of HEH[EHP] is ca. 4.1,[51] and our nitrate dependence experiments (Figure ) clearly show the
coextraction of 3 nitrates per trivalent cation (n = 3) without the loss of hydrogen ions from HEH[EHP] to the aqueous
phase according to eq . If H(EH[EHP])2– or 2 equiv of EH[EHP]− was to form by the release of H+ from theHEH[EHP] dimers to the aqueous phase, the experimental nitrate activity
dependence in Figure would have a slope of 1 (n = 2) or −1 (n = 1), respectively (see eq , eq S15, and the Supporting Information). Instead, the experimentally
verified lack of H+ release to the aqueous phase suggests
either that HEH[EHP] coordinates to themetal center as a neutral,
protonated species[45] or that HEH[EHP] deprotonates
during the formation of the extracted complex with the released hydrogen
ion remaining in the organic phase and forming a molecule of nitric
acid by reacting with one nitrate anion.
Computational modeling
of tan class="Chemical">he possible extracted complexes provides
further insight into HEH[EHP] coordination and the likely protonation
state of these complexes. The complexation energy of Eu(TEDGA)3·3NO3 was calculated to be −9.7 kcal/mol
compared to −11.3/–11.1 and −22.9/–25.7
kcal/mol for Eu(TEDGA)2(H(E[EP])2)·2NO3 geometries A/B and Eu(TEDGA)2(HE[EP])2·3NO3 geometries C/D, respectively (Figure and Table ). Although the formation of M(TEDGA)3·3NO3 is exothermic, the substantially greater
stability of the ternary metal–diglycolamide–phosphonic
acid complexes suggest that fully protonated HEH[EHP] can readily
displace a TEHDGA molecule, since HEH[EHP] is available at 15 times
the concentration of TEHDGA in theALSEP organic phase. Indeed, the
biphasic spectrophotometric titrations (Figure ) indicate that an appreciable concentration
of a 1:2:1 Nd:TEHDGA:(HEH[EHP])2 complex is already present
with 0.015 M concentrations of HEH[EHP] and 0.05 M TEHDGA. Consequently,
M(TEHDGA)3·3NO3 is unlikely to be present
during ALSEP extractions.
Table 4
Calculated Gibbs
Free Energies of
Complexation for the Optimal Geometries of Eu(TEDGA)·3NO3, Eu(TEDGA)2(H(E[EP])2)·2NO3, and Eu(TEDGA)2(HE[EP])2·3NO3
complexa
complexation energy (kcal/mol)
Eu(TEDGA)·3NO3
–9.7
Eu(TEDGA)2(H(E[EP])2)·2NO3 (A)
–11.3
Eu(TEDGA)2(H(E[EP])2)·2NO3 (B)
–11.1
Eu(TEDGA)2(HE[EP])2·3NO3 (C)
–22.9
Eu(TEDGA)2(HE[EP])2·3NO3 (D)
–25.7
Extractants truncated from 2-ethylhexyl
to ethyl chains for computational feasibility.
Extractants truncated from an class="Chemical">2-ethylhexyl
to an class="Chemical">ethyl chains for computational feasibility.
Substitutions of an class="Chemical">HEH[EHP] for n class="Chemical">TEHDGA
in the extracted complex of
theALSEP system also appear to affect thenitrate coordination environment
profoundly. In the M(TEHDGA)3·3NO3 complex,
thenitrate anions required to balance the positive charge of themetal cation sit in the outer coordination sphere in clefts between
the coordinated TEHDGA ligands,[19,52] likely interacting
with coextracted water molecules present in the outer coordination
sphere.[53] Two of these 3-fold interligand
nitrate binding clefts are lost when one TEHDGA molecule is replaced
with HEH[EHP] to form theALSEP system’s extracted complex.
Consider the hypothetical ALSEP extraction reaction:It would display a nitrate
activity dependence slope of 1 with n = 2 (eq S15 and the Supporting Information) and produce complexes with geometries A or B.
For these complexes, the computations indicate that replacing one
TEHDGA with H(EH[EHP])2– will cause both
nitrate anions to migrate out of the one remaining cleft between TEHDGA
molecules. In geometry A, thenitrates occupy spaces between theTEHDGA
and HEH[EHP] molecules. In the case of geometry B, one of thenitrate
anions repositions further to form an inner-sphere complex with the
Eu ion (Figure ).
The changes in thenitrate environment are even more profound for
the most stable complexes studied: geometries C and D of Eu(TEDGA)2(HE[EP])2·3NO3. Unlike the complexes
represented by geometries A and B, the stoichiometry of these complexes
C and D matches the experimentally determined stoichiometry and the
overall extraction equilibrium.In both geometries C and D,
the computations suggest that one nitrate anion will remain in the
cleft between the two TEHDGA molecules and a second nitrate will interact
with the acidic hydrogen of one of the neutral HEH[EHP] molecules.
In geometry C the third nitrate anion resides in the outer coordination
sphere between one TEHDGA and one HEH[EHP] molecule. However, the
bond lengths suggest that the third nitrate reacts with the acidic
hydrogen from the other HEH[EHP] molecule in the complex to form a
molecule of HNO3 that hydrogen bonds to the coordinated
oxygen of theHEH[EHP] (Figure ). In complex C, one inter-HEH[EHP]hydrogen bond remains
intact, and the (HEH[EHP])2 dimer is coordinated to themetal in a bidentate fashion. Theequilibrium for this process can
be represented asIn geometry D, similar to geometry C, one
nitrate remains in the inter-TEHDGA cleft, a second nitratehydrogen
bonds to one fully protonated HEH[EHP] molecule, and one of thehydrogen
bonds between HEH[EHP] monomers remains intact. However, in complex
D, theHEH[EHP] dimer becomes monodentate. TheHEH[EHP] monomer that
hydrogen bonds to nitrate is rotated away from the Eu and does not
coordinate to themetal center. Instead, the third nitrate anion moves
into the inner coordination sphere to give an 8-coordinate complex.
Highlighting this change in nitrate coordination, theequilibrium
for the formation of complex D can be written asFunctionally, the low-energy complexes with
geometries C and D form by sharing an acidic hydrogen between a HEH[EHP]
molecule and a nitrate ion in the organic phase. Stoichiometrically,
complexes C and D are indistinguishable, they match our experimentally
determined speciation of the organic-phase complex, and they present
similarly favorable complexation energies compared to the other complexes
studied.Complexation-driven an class="Chemical">HEH[EHP] proton transfer to create
n class="Chemical">hydrogen-bonded
HNO3 in the organic phase is an intriguing reaction mechanism
for trivalent f-element extraction in theALSEP process. It explains
the observed synergistic enhancement to the extraction caused by HEH[EHP]
addition, the observed stoichiometry under common extracting conditions,
and the eventual loss of HEH[EHP] dependence at very high aqueous
acidities,[31] and it is further supported
by FTIR studies of Eu extracted from 3 M HNO3 into TEHDGA/HEH[EHP]/n-dodecane.[18] Moreover, the predicted
stabilization of HNO3 in the organic phase is not surprising.
The ability of HEH[EHP] to extract nitric acid on its own is known
(Figure S3), and the coextraction of nitric
acid with An(NO3)3 or Ln(NO3)3 has been reported for bifunctional extractants such as octyl(phenyl)-N,N-diisopropylcarbamoylmethanephosphine
oxide (CMPO),[54,55] malonamides,[56] and diglycolamides TODGA and TEHDGA.[20,57] In theALSEP process, the nondissociated acidic hydrogens of HEH[EHP]
appear to exert a stabilizing influence on the extracted complexes,
enabling the formation of f-element-TEHDGA-HEH[EHP] complexes under
conditions where HEH[EHP]’s conventional metal–proton
exchange equilibrium is strongly hindered.
Conclusions
Our
studies conclusively demonstrate tan class="Chemical">he formation of a previously
suspected mixed trivalent f-element–TEHDGA/HEH[EHP] complex
in the organic phase under the acidic extracting conditions of theALSEP process (2–4 M HNO3). The mutual presence
of a strong neutral extractant and a weakly acidic extractant in theALSEP organic phase provides unique opportunities for the synergistic
extraction of An3+ and Ln3+ from aqueous molar
nitric acid solutions. As a more basic extractant, the dialkylphosphonic
acid HEH[EHP] in theALSEP system displays different behavior than
thedialkylphosphoric acidHDEHP in similar mixed-extractant systems.[11,43,58] HEH[EHP] does not appear to undergo
deprotonation when it cooperates with TEHDGA to extract trivalent
lanthanide or actinide nitratesunder theALSEP process extracting
conditions. While the strongly acidic aqueous phases inhibit the deprotonation
of HEH[EHP], the relatively high concentration of HEH[EHP], strong
PO–M bonds, a smaller coordination footprint (mono or bidentate
(HEH[EHP])2 vs tridentate TEHDGA), and an ability to hydrogen
bond with nitrate anions work together to boost the stability of the
mixed M(TEHDGA)2(HEHEHP)2·3NO3 complex in theALSEP organic phase.
Methods
Materials
Chemicals were purchased from Sigma-Aldrich
except where otherwise indicated. The extractant HEH[EHP] was obtained
from BOC Sciences (95%) and purified to ≥98% by the third-phase
method,[59] as confirmed by 31P NMR and acid–base titration in anethanol–water mixture
(80:20). Stock extractant solutions were made by combining weighed
amounts of TEHDGA (Eichrom Technologies, > 99%) and the purified
HEH[EHP]
in n-dodecane (anhydrous, 99%) and then diluting
to a known volume.Aqueous phases were made from nitric acid
(Baker, ULTREX II) and standardized by titration with sodium hydroxide.
Solutions of lanthanidenitrates were prepared at the desired acidities
from stock solutions of neodymium(III)nitratehexahydrate (99.9%
trace metals basis) or a europium nitrate solution prepared by dissolving
a weighed amount of europium oxide (Treibacher Industrie AG, 99.99%
REE) in nitric acid. Theneodymium stock solution was standardized
by titration with a standard EDTA solution (Fisher Scientific), xylenol
orange indicator, and a saturated solution of hexamethylenetetramine
(>99%) to buffer the titration at pH 5–5.5. The pH was adjusted
with 4 M HNO3 and monitored throughout the titration with
a ThermoOrion Ross semimicro pH electrode. Aqueous solutions were
prepared using Millipore 18 MΩ cm deionized water, while deuterium
oxide (99.9 atom % D) and 65 wt % nitric acid-d in
D2O (99 atom %) were used for theD2O luminescence
experiments.
Absorption Spectroscopy
Aqueous
phases of neodymium
were made to 0.01 M and diluted with the relevant titrated nitric
acid solutions. Acid concentrations were chosen to optimize extraction
of Nd while avoiding the formation of a third phase; therefore, the
acid concentrations at ALSEP process relevant conditions were not
always possible to use. Organic-phase solutions were prepared by weight
and pre-equilibrated through contact with twice their volume of the
appropriate metal-free aqueous phase and then vortexing for 30 s,
followed by centrifugation for phase separation. This process was
repeated twice to complete pre-equilibration. Metal extractions were
performed by vortexing equal volumes of pre-equilibrated organic phases
with aqueous phases containing 0.001–0.01 M Nd(NO3)3 at the desired nitric acid concentration for 5 min
followed by 5 min of centrifugation before the phases were separated.Spectra of both the aqueous and organic phases were obtained on
a Varian Cary 300 Spectrophotometer in 1.000 cm quartz cuvettes from
480 to 850 at 0.2 nm resolution. Theequilibrium concentrations of
neodymium in the organic phases were calculated from the total amount
of Nd in the system and the difference in the initial and final spectra
of the aqueous phase.The influence of HEH[EHP] on the stoichiometry
of the organic-phase
Nd complexes also was studied spectroscopically. Organic phases containing
0.05 M TEHDGA and varying amounts of HEH[EHP] were made by combining
aliquots of 0.2 M HEH[EHP]/n-dodecane, 0.2 M TEHDGA/n-dodecane, and n-dodecane, all of which
had been pre-equilibrated with 2 M HNO3. These solutions
were pipetted in appropriate ratios to vary theHEH[EHP] concentration
only. In addition, 0.05 M HEH[EHP]/n-dodecane and
0.05 M TEHDGA/n-dodecane organic phases were also
tested. Culture tubes containing equal volumes organic phase and 0.01
M neodymium in 2 M HNO3 were initially placed into a 35
°C water bath for 30 min, and the samples were then vortexed
for 30 s and placed back in thewater bath for several minutes. These
intervals were repeated until a total vortexing time of 5 min was
reached. Phases were separated, and spectra were collected between
490 and 610 nm in a 1.000 cm jacketed cuvette held at 35 °C by
a recirculating water bath.
TRLFS measurements were
performed on an class="Chemical">europium-loaded organic solutions
prn class="Chemical">epared by extraction from light and heavy water. Deuterated samples
were made entirely with deuterated chemicals except for theeuropiumnitrate stock solution, which had been prepared in light water and
contributed no more than 2 atom % 1H to the aqueous solutions.
Organic solutions were pre-equilibrated with appropriate D2O/DNO3 or H2O/HNO3 aqueous phases
and then contacted with a fresh aqueous phase containing europium
similar to the procedure for the absorption spectroscopy experiments.
Fluorescence spectra were obtained using an Edinburgh Instruments
LP980 spectrometer with a Continuum Surelite nanosecond Nd:YAG laser
as the excitation source. The excitation wavelength was fixed at 394
nm, while the emission spectrum was recorded in the range of 300 to
800 nm. Fluorescence lifetimes were obtained by monitoring the lifetimes
of theeuropium 5D0 → 7F transitions (J = 0–4)
at the wavelengths of interest. Kinetic traces at each wavelength
were fit as a single exponential decay, and thegoodness of fit was
evaluated using the reduced chi-squared (χ2). The
coordination environment of europium was elucidated based on the fluorescence
emission spectra and fluorescence lifetime of the complexes, and the
change in lifetime between nondeuterated and deuterated samples. In
addition, the number of inner-sphere-coordinated water molecules (NH2O) was calculated using Horrock’s equation:where kH and kD are the fluorescence
decay rate constants in ms–1 for samples prepared
in H2O and D2O, respectively.[60−63] Theequation calculates the number
of water molecules to within 0.5 molecules.
Extraction Experiments
The procedures for studying
the extraction of nitric acid are described in the Supporting Information. All radiotracer experiments were performed
using radiochemically pure 241Am (Eckert and Ziegler) dissolved
in 2 M HNO3. Extractions were performed by pipetting equal
volumes of pre-equilibrated organic and aqueous phases containing
the desired concentrations of acid and extractant into glass culture
tubes, spiking with 2 μL of the 241Am solution, and
vortexing for 5 min, followed by 5 min of centrifugation. The phases
were separated, and aliquots of each phase were taken for liquid scintillation
counting on a Packard 2500 TR liquid scintillation analyzer using
an Ecoscint liquid scintillation cocktail. Distribution ratios (D) were calculated from the
following equation:assuming equal counting efficiency of the
aqueous and organic phases. All experimental uncertainties are reported
at two standard deviations, and error bars are not shown in figures
if the uncertainty is smaller than the data points.
Theoretical
Calculations
All calculations were performed
using Gaussian 09 software.[64] Tan class="Chemical">he alkyl
chains of tn class="Chemical">he ligands under investigation were truncated from 2-ethylhexyl
chains to ethyl chains for ease of calculation, as previously reported.[65] Therefore, the computational analysis replaces
TEHDGA with TEDGA and HEH[EHP] with HE[EP] (see Figure S1 for structures of the truncated ligands). Theeuropium
was treated as a trivalent atom containing a +3 charge and a multiplicity
of seven. The geometries of the hypothesized europium–ligand
complexes were optimized using Becke’s three parameter exchange
functional with Lee–Yang–Parr’s correlation functional
(B3LYP) and the 6-31G(d,p) basis sets for all second and third row
elements.[66,67] A relativistic effective core pseudopotential
and corresponding basis set were employed (Stuttgart RSC 1997 ECP)
for treatment of theeuropium cation.[68] These calculations were performed without explicit or implicit solvation
since theeuropium extraction complexes under investigation will be
present in thedodecane organic phase, which has a dielectric constant
near that of the gas phase. Calculated stationary points were identified
as a minimum by verifying the lack of any imaginary frequencies. Complexation
energies were calculated from the difference of the Gibbs Free Energies
of the products and reactants for each complex. Even with truncation
of the alkyl chains, the number of degrees of freedom of each complex
is large. Therefore, the exploration of the conformational spaces
was necessarily limited to dozens of conformations for each complex.