The
extraction of lanthanide series by Cyanex301, i.e.,
bis(2,4,4-trimethylpentyl)dithiophosphinic acid (HC301), has been
modeled by density functional theory calculation, taking into account
the formation of both inner- and outer-sphere complexes. The inner-sphere
complex Ln(C301)
3
and the outer-sphere complex Ln(H
2
O)
9
(C301)
3
are optimized, followed by
the analysis of interaction energy, bond length, Laplacian bond orders,
and Mulliken populations. The covalency degree increases in Ln–S
and Ln–O bonds in the inner- and outer-sphere complexes, respectively,
as the lanthanide series is traversed. Mulliken population analysis
indicates the important role of the 5d-orbital participation in bonding
in the formation of inner- and outer-sphere complexes. Two thermodynamic
cycles regarding the formation of inner- and outer-sphere complexes
are established to calculate the extraction Gibbs free energies (Δ
G
extr
), and relaxed potential energy surface
scan is utilized to model the kinetic complexation of C301 anion with
hydrated metal ions. Light lanthanides can form both inner- and outer-sphere
complexes, whereas heavy lanthanides only form outer-sphere complexes
in biphasic extraction. After adopting the data of forming inner-sphere
complex for light Ln(III) and that of forming outer-sphere complexes
for heavy Ln(III), the trend of the calculated −Δ
G
extr
agrees very well with that of the experimental
distribution ratios on crossing the Ln(III) series. Results from this
work help to theoretically understand the extraction behavior of Cyanex301
with respect to different Ln(III).