Americium is a highly radioactive actinide element (e.g., t1/2,Am−241 = 432.6 a) that is formed in nuclear reactors and, thus, is present in used nuclear fuel as a minor actinide. Americium has also been introduced to the environment through nuclear weapons testing and aging legacy nuclear waste infrastructure. The emplacement of used nuclear fuel in deep underground repositories is being proposed as the best strategy for its long-term disposal [1–6]; commercially-available bentonite clay, which is mostly composed of the aluminum phyllosilicate mineral montmorillonite, will be used as a backfill material. Aluminum (hydr)oxide minerals like corundum (α-Al2O3), γ-alumina (γ-Al2O3), and gibbsite (γ-Al(OH)3) are ubiquitous in the subsurface environment and may influence the fate and transport of actinides (including americium). Furthermore, the aluminol sites present in aluminum (hydr)oxide minerals are also present in bentonite [7, 8]. Therefore, adsorption of americium on aluminum (hydr)oxide minerals is an important phenomenon to study.
Surface complexation modeling is a predictive tool that is often used to elucidate the adsorption of metal cations on minerals [9–12]. Surface complexes are analogous to aqueous complexes; however, unlike the aqueous complexation reaction, surface electrostatic effects, which are dependent on surface potential, need to be considered for the formation of surface complexes [10, 13–16]. Different types of electrostatic model frameworks can be used to develop surface complexation models for a given metal-mineral system and vary from one another in how they treat charge distribution at the mineral surface. Diffuse double layer (DDL) and charge distribution multi-site complexation (CD-MUSIC) are widely used electrostatic modeling frameworks for surface complexation modeling [9–11, 17]. In the CD-MUSIC framework, the charge on the mineral surface is distributed in three planes, which is a more realistic depiction of the charge distribution as compared to the DDL framework, where a point charge distribution is assumed. Although the CD-MUSIC framework is more realistic than DDL framework, a larger number of parameters are required for defining the surface complexation modeling under the CD-MUSIC framework. These parameters are the inner and the outer layer capacitances and the charge distribution coefficients for the three planes (i.e., ΔZ0, ΔZ1, ΔZ2). The capacitances denote the rate of change of the surface potential as a function of distance from the surface. In the DDL framework, the surface potential remains constant between the naught plane and the d-plane and then decreases exponentially beyond the d-plane. Apart from the DDL and the CD-MUSIC framework, constant capacitance model (CCM) is also used for elucidating the surface electrostatics. In CCM, the surface potential decreases linearly from its maximum value on the naught plane to zero on the d-plane, and no surface effect is present beyond the d-plane. The details of different electrostatic models and their fundamentals have been explained in literature [18].
Europium is a lanthanide metal that resembles americium in its ionic size (1.066 Å and 1.09 Å for Eu and Am, respectively, in 8-fold coordination), oxidation state, coordination number, and the properties of its coordination complexes [19, 20]. Therefore, europium is widely used as a chemical analog for americium. Surface complexation models have been developed for various metal ion sorption onto iron oxide minerals for a large and varied dataset sourced from multiple different studies which denote large variation in the input conditions [10, 21–25]. Many surface complexation models for Eu(III)-γ-alumina system have been developed since 2000 [26–28]. Rabung et al. [28] report the results of parameter optimization for Eu(III) adsorption to γ-alumina as a function of pH using a CCM electrostatic framework and two site types (strong and weak). The log K values for the same surface complex at different pH values varied, but were all close to 2.5; no global optimization was completed in order to obtain a unique log K for this system. Almost a decade later, two more surface complexation modeling studies [26, 27] were published for the Eu(III)-γ-alumina system. Kumar et al. [29] assume the presence of only a single surface site on γ-alumina and both monodentate and bidentate surface complexation of Eu(III) on the sorbent surface. The average log K values for the monodentate and the bidentate surface complexes are 2.22 ± 0.35 and − 4.99 ± 0.04, respectively. Although sorption edge data were collected as a function of Eu(III) concentration (0.1–100 µM), no global optimization was reported. This could be because of the limitations created by the use of FITEQL as the optimization tool, where only a single set of input data can be introduced for a given optimization problem. Morel et al. [27] generate a surface complexation model based off a single set of adsorption edge data. They assume the presence of only one site and surface complex for optimization. The log K of the surface complex is -1.2, which was almost six times weaker than the log K of the same surface complex optimized by Kumar et al. [26]. To the best of our knowledge, no surface complexation models have been developed for the Eu(III)-corundum and Eu(III)-gibbsite systems.
The objectives of this work are to compile all the available data for Eu(III) adsorption onto corundum, γ-alumina, and gibbsite, use DDL and CD-MUSIC electrostatic frameworks to develop surface complexation models describing Eu(III) adsorption to these minerals, and determine whether the resulting surface complexation models could also describe the corresponding Am(III)-corundum, Am(III)-γ-alumina and Am(III)-gibbsite systems.