We
propose the statistical thermodynamic model for the prediction
of the liquid–liquid extraction efficiency in the case of rare-earth
metal cations using the common bis(2-ethyl-hexyl)phosphoric acid (HDEHP)
extractant. In this soft matter-based approach, the solutes are modeled
as colloids. The leading terms in free-energy representation account
for: the complexation, the formation of a highly curved extractant
film, lateral interactions between the different extractant head groups
in the film, configurational entropy of ions and water molecules,
the dimerization, and the acidity of the HDEHP extractant. We provided
a full framework for the multicomponent study of extraction systems.
By taking into account these different contributions, we are able
to establish the relation between the extraction and general complexation
at any pH in the system. This further allowed us to rationalize the
well-defined optimum in the extraction engineering design. Calculations
show that there are multiple extraction regimes even in the case of
lanthanide/acid system only. Each of these regimes is controlled by
the formation of different species in the solvent phase, ranging from
multiple metal cation-filled aggregates (at the low acid concentrations
in the aqueous phase), to the pure acid-filled aggregates (at the
high acid concentrations in the aqueous phase). These results are
contrary to a long-standing opinion that liquid–liquid extraction
can be modeled with only a few species. Therefore, a traditional multiple
equilibria approach is abandoned in favor of polydisperse spherical
aggregate formations, which are in dynamic equilibrium.