Phyllite/bentonite mixture—an alternative effective buffer material for a geological disposal of radioactive waste

The use of phyllite (Phy) instead of quartz in mixtures with bentonite (B) is recommended as a buffer material for engineering barriers in a geological repository of nuclear waste. The recommendation is based on experimentally determined sorption properties of various Phy/B mixtures. The adsorption capacity of Phy/B mixtures (Phy/B: 75/25, 50/50, and 25/75), the removal efficacy of Eu(III) ions (an analog for fissiongenic lanthanides and actinides), and the rate of their binding reaction were studied using the batch adsorption equilibrium and kinetic experiments at different Eu(III) initial concentrations, solution pH, and solution to adsorbent (L/S) ratio. The adsorption capacity of the Phy/B mixtures increased with the increased bentonite content in the mixture depending on the L/S ratio and solution pH. The highest increase in the adsorption capacity of the Phy/B mixtures compared to phyllite was observed for the Phy/B proportions of 25/75 and 50/50. The rate of the Eu(III) adsorption was the best fitted by the pseudo-second-order kinetic model indicating that the adsorption rate was controlled by chemisorption. The Sips model provided the best correlation of the adsorption experimental data, indicative of more than one adsorption site. The results of this study show the advantage of the Phy/B mixtures in immobilizing Eu and certain fission products by combining adsorption properties of the materials. Supplementary Information The online version contains supplementary material available at 10.1007/s11356-023-31102-6.


Introduction
One of the specific goals of radioactive waste disposal is to "inhibit, reduce, and delay the migration of radionuclides at any time from the waste to the accessible biosphere" (IAEA Safety Standards 2011).To achieve this goal, the concept of a multibarrier disposal system in an underground (geological) repository has been developed.The system is made up of engineered and geological (host rock) barriers.The engineered barriers comprise the packaging with vitrified high-level nuclear waste (HLNW) or spent nuclear fuel (SNF), the clay barrier (buffor), and the backfill (Kale and Ravi 2021).The long-term safety of a geological repository for nuclear waste relies on the performance of each of these barriers.However, the role of the clay barrier is crucial in the repository safety assessment (e.g., Itälä 2009; Yang et al. 2014).The principal function of a clay barrier material around waste packages is to protect them from the intrusion of groundwater and, if such an event occurred, to prevent or limit the migration of released radionuclides.
Responsible Editor: Georg Steinhauser Highlights • Eu(III) adsorption on phyllite/bentonite mixtures as a function of pH, solution to rock ratio, and of initial concentration.
• The pseudo-second-order kinetic model as the best fit for all Phy/B mixtures indicating the chemisorption as controlling adsorption rate.
• The Sips isotherm as the best adsorption experimental model suggests more than one adsorption site.
• Very good effect of bentonite on phyllite Eu(III) adsorption capacity for the 50Phy/50B mixture.
• Chemisorption controls adsorption rate of Eu(III) on all Phy/B mixtures.
• Hydrogen bond formation and precipitation of Eu(OH) 3 are predominant mechanisms of Eu binding by phyllite/bentonite mixtures.
• Very good fit of Sips isotherms to experimental data allows estimation of coefficient of determination values.
Extended author information available on the last page of the article One of the possible scenarios considered in the longterm performance assessment of a radioactive waste repository is the release of radionuclides from breached waste containers and their migration into groundwater through the barriers.The distribution and mobility of the released radionuclides are controlled by the adsorption and ion exchangeable properties of the clay barrier and backfill materials.Moreover, adsorption is important in retarding the potential migration of radionuclides to the aquatic systems and biosphere.Understanding the kinetic process and removal efficacy/adsorption capacity is fundamental for the prediction of radionuclide migration in the surroundings of deep geological disposal sites of high-level nuclear wastes.The studies mainly concern the adsorption capacity of radionuclides on pure minerals or their modified forms as well as on composite materials (Bradbury and Baeyens 2005;García et al. 2019;Hanza et al. 2018;Sharma and Tomar 2011;Verma et al. 2019;Virtanen et al. 2018).To design a highly effective geochemical reactive barrier to immobilize radionuclides in a waste repository, different materials must be parallel studied to analyze their adsorption capacity and adsorption kinetics.However, there have been only a few studies on the adsorption kinetics in this respect (e.g., Baumer and Hixon 2018;Wang 2018).
Currently, bentonite or bentonite mixed with quartz sand is considered to be used as a barrier and backfill material in the deep geological repository of HLNW and SNF because of their low permeability, plasticity (ability to swell), very good adsorption/ion exchange properties, and high thermal conductivity (Itälä 2009;Kónya et al. 2005;Pusch 1992Pusch , 2006;;Sellin and Leupin 2014;Yang et al. 2014).However, bentonites undergo dehydration, mechanical degradation, and structural transformation decreasing their adsorption capacity at a temperature above 100 °C expected at the contact with waste containers, and other rocks of high thermal resistance and high adsorption capacity for radionuclides are being sought (Xiao et al. 2013).Previous studies show that argillaceous phyllites, while less efficient adsorbents than bentonites, have high adsorption capacity for low initial concentrations of actinides (< 5•10 −5 M) and much higher adsorption capacity than quartz, especially at low pH (Kar et al. 2011;Kyzioł-Komosińska et al. 2019a).Moreover, they are thermally stable up to 400 °C.Therefore, argillaceous phyllite can be considered as a component of a mixture with bentonite instead of quartz to take advantage of phyllite adsorption capacity and proper mechanical properties (compensation for volume changes in swelling/shrinking bentonite), particularly at elevated temperatures.Phyllite clays have recently been used as barrier material in municipal waste landfills due to their low permeability (Garzón et al 2016).
The main purpose of this study is to determine the removal efficacy of Eu(III) ions by phyllite/bentonite (Phy/B) mixtures.Europium is an analog for fissiongenic lanthanides and actinides: Am and Cm; present in HLNW and SNF (NRC Regulation Part 63 2016;Sun et al 2020).Various Phy/B proportions were used in the adsorption experiments to determine the most efficacious mixture for immobilization Eu(III) ions.The rate of the Eu(III) ion removal from solution, the adsorption capacity of the Phy/B mixtures, and mechanisms of adsorption were investigated.

Adsorbents
Phyllite used in this study was from the Dewon-Pokrzywna phyllite deposit in the northern foothills of the Opava Mountains in the Eastern Polish Sudetes (Sawicka et al. 2018).Bentonite was from the Miocene Kopernica deposit in Slovakia (Górniak et al. 2016).Mineral and chemical compositions together with textural and physical properties of these minerals are given in Kyzioł-Komosińska et al. (2019a).Powdered phyllite (Phy) and bentonite (B), both with particle size < 10 µm, were thoroughly mixed in weight proportions of 75/25, 50/50, and 25/75.The properties of phyllite, bentonite, and their mixtures relevant for the evaluation of their adsorption effectiveness are given in Table 1.

Adsorption experiments
The batch equilibration method was applied to examine the adsorption of Eu(III) from EuCl 3 •6H 2 O solution by the Ph/B mixtures under atmospheric pressure and at Moreover, a high initial Eu(III) concentration may reflect the worst-case scenario of mass radionuclide release from the waste package after some catastrophic event.
Solution initial pH was adjusted to 4.5 and 7.0 by drop-wise addition of appropriate amounts of 0.01 M HCl or 0.01 M NaOH, respectively.The ionic strength was 0.010 M NaCl.The contact time of adsorbents and solutions was 360 min, i.e., more than necessary to achieve equilibrium inferred from the adsorption kinetics determination.The initial (C 0 ) and equilibrium (C eq ) concentrations of Eu(III) in supernatant solutions were measured after centrifugation and filtration either by the ICP-OES or by ICP-MS spectrometry depending on the Eu(III) concentrations.The removal efficacy (RE) of Eu(III) ions was calculated from the equation: Results obtained during this study allowed for the determination of the maximum initial concentrations (C 0max ) for which Eu(III) ions were fully adsorbed, i.e., RE = 99.9% and C eq ≤ 0.001 mg/L.
Adsorption capacity (q) was calculated using the equation: where m mineral mass (g); V solution volume (L).
Duplicate samples were measured, and the standard error in the readings was less than 4%.
Two-parameter isotherm models of Freundlich (Freundlich 1906), Langmuir (Langmuir 1916), and Dubinin-Raduskevich (Dubinin 1960) and the three-parameter Sips isotherm model (Sips 1948) were used to estimate maximum adsorption capacity of the Phy/B mixtures and to determine the affinity between Eu(III) ions and adsorption centers in Phy and B and the bond energy (Table S1 in Supporting Information).The isotherm parameters were determined by non-linear regression using algorithms based on the Levenberg-Marquardt method (Statistic ver.9.0 software).Two non-linear error functions (Foo and Hameed 2010) were examined in addition to the determination coefficient (R 2 ) for fitting the calculated isotherm to the experimental data S2 in Supporting Information).The correlation of experimental results to adsorption isotherms can help to understand the adsorption mechanisms and the degree of the heterogeneity of the adsorbent surface.

Adsorption kinetics
To understand adsorption/desorption processes, both thermodynamic equilibrium and kinetics must be considered (Wang 2018).While thermodynamic data provide information only on the final state of a system, the study of kinetics reveals changes in chemical properties in time and is concerned with the rate of these changes (Azizian 2004).
The Phy/B mixtures were used in the batch kinetic experiments.The effect of the contact time was examined within 2.5-1440-min time-frame for C 0 of 10 and 200 mg/L at pH 4.5 and 7.0 for L/S of 100:1.To determine the impact of the L:S on the adsorption rate, tests were carried out for L/S of 100:1 and 500:1 at C 0 of 10 and 200 mg/L, and pH 7.0.Time elapsed to reach equilibrium in the soluteadsorbent system was further considered in the evaluation of the experimental results.Two kinetic-based reactions, i.e., Lagergren pseudo-first order (PFO) (Begg et al. 2018) and pseudo-second order (PSO) (Blanchard et al. 1984;Ho and McKay 1999), were applied to find the rate-determining step (Table S3 in Supporting Information).The kinetic parameters were estimated by non-linear regression.Tran et al. (2017) recommended the use of nonlinear forms of PFO and PSO to describe the kinetics of adsorption.Due to the transformation of nonlinear to linear forms, the units of the Y and X axes were changed.Two non-linear error functions (Table S2 in Supporting Information) were examined in addition to the coefficient (R 2 ) to evaluate the fit of the equations to the experimental data (Foo and Hameed 2010).

Auxilliary methods: scanning electron microscopy (SEM), attenuated total reflection Fourier transform infrared spectroscopy (ATR-FTIR)
Operating conditions for these methods are given in the Supporting Information.

Adsorbents properties
The Phy/B mixtures differed significantly in physicochemical properties depending on the B content.The addition of 25% B into the Phy/B mixture caused a sixfold increase in CEC up to 21.6 cmol + /kg, i.e., to the value observed for kaolinite (Ma and Eggleton 1999).The CEC increased 12 times to 42.11 cmol + /kg in mixtures with 50% B. A similar effect of the B content was observed for the mixtures' SSA (Table 1).Values of pH and pH PZC of the Phy/B mixtures suspension increased from 6.87 to 7.15 and from 6.14 to 6.01, respectively, with an increase in B content.

Europium precipitates
Scanning electron microscopy (SEM)/EDS study of bentonite, phyllite, and the Phy/B mixtures (Table S4 in Supporting Information) after adsorption experiments revealed precipitates of the Eu-phase (Fig. 1).Elongated crystals of the Eu-phase were 5 μm long and < 4 μm wide.The rounded edges of those crystals suggest partial dissolution at the crystal-solution boundary layer.From semi-quantitative EDS analyses, the Eu contents in the crystals precipitated on bentonite and phyllite are on average 69.31 and 68.74 wt.%, respectively.These values are closer to the Eu(OH) 3 stoichiometry (74.86 wt.% Eu) than to Eu 2 O 3 (86.36wt.% Eu) (Table S5 in Supporting Information).
The Eu precipitation was pH-dependent and occurred only at pH 7.0.The bulk contents of Eu adsorbed by montmorillonite in bentonite and by and chloritized biotite in phyllite were about 4.5 and 2.5 wt.%, respectively, regardless of solution pH.

Infrared spectroscopy
Spectra of phyllite collected before and after Eu adsorption are not significantly different.The expected band shifts due to the formation of hydrogen bonds caused by the Eu(OH) 2+ bound to -Al-OH were not observed in accord with suggestions by Bradbury and Baeyens (2005).The only significant change seen in the IR spectrum of phyllite after the adsorption experiment is the weakening of the 750 cm −1 band assigned to the SiAl-O-OH libration of chlorite hydroxyl sites.FT-IR spectra of bentonite after Eu adsorption show more distinct differences compared to phyllite.The 3620 cm −1 band is prominent, and its intensity is much higher than neighboring bands.Russell et al. (1970) interpreted the 3624 cm −1 montmorillonite band as derived from OH stretching vibrations in Al 2 OH and AlMgOH.Bending vibrations of these groups give bands at 915 cm −1 (Al 2 OH) and much weaker at 843 cm −1 (AlMgOH).After Eu adsorption, the intensity of the 915 cm −1 band increased slightly, possibly as a result of Eu binding to OH groups.Weakening of bands intensity after Eu adsorption occurred below 900 cm −1 , i.e., in the region of OH deformation vibrations in aluminol and silanol groups (Vaculíková and Plevová 2005).The description of the Phy/B mixture spectra is presented in the Supporting Information.

Adsorption kinetics
Adsorption kinetics was examined for 1400 min; however, the experimental results showed that 180 min was sufficient for reaching an equilibrium state in all Eu-adsorbent systems.There was a rapid and efficient one-step adsorption for C 0 of 10 mg/L, whereas for C 0 of 200 mg/L, a two-step reaction was observed with  2).The two-step process is attributed to a fast, diffusion-controlled surface reaction followed by a rate-limiting step that can be explained by a variety of processes including surface precipitation, For the PFO and PSO models, the curve fitting parameters (adsorption rate constants, k 1 , and k 2 ) and the equilibrium adsorbed concentration, q, were calculated (Table 2).High correlation coefficients (R 2 > 0.95) were obtained for the non-linear plots of the pseudofirst Lagergren and pseudo-second-order rate equations.Statistical analysis showed that the pseudo-second-order model resulted in better correlations compared to the pseudo-first-order model for C 0 of 200 mg/L and all pH values and L/S ratios.The pseudo-first-order equation fitted well experimental data for C 0 of 200 mg/L for the first 20-30 min of the adsorption process (Fig. 2).According to Ho and McKay (1998), the pseudo-firstorder model applies only to a limited fraction of the reaction range.For this reason, the pseudo-second-order model was used to determine the rate constant, k 2 , and the Eu(III) equilibrium adsorbed concentration, q, with respect to C 0 , L/S, and pH.Estimated q was similar to or a little higher than experimentally determined adsorption capacity (Table 2) and depended on the B content in the Phy/B mixtures for C 0 of 200 mg/L.At C 0 of 10 mg/L, the q was constant regardless of the B content in the mixtures (Fig. 2a).A small difference in k 2 determined for various Phy/B mixtures suggests that B content in the mixtures was less important than other experimental conditions.The k 2 strongly depended on the initial Eu(III) concentration, pH, and L/S ratio.It was higher for L:S of 100:1 than for 500:1 for both initial Eu concentrations and all of the Phy/B mixtures.However, the effect of L/S was higher for C 0 of 10 than 200 mg/L (Table 2).The increase in k 2 with increasing content of adsorbents presumably was due to an elevated number of total adsorption sites, including strong adsorption sites in montmorillonite (Bradbury andBaeyens 2002, 2005;Bradbury et al. 2005).
The Eu(III) adsorption rate increased with increasing pH (Fig. 2b) at C 0 of 10 mg/L.This observation can be explained by the mineral surface becoming increasingly negatively charged and, hence, attracting the Eu cations due to electrostatic interactions.There is no correlation between k 2 and pH for C 0 of 200 mg/L, apparently, because the electrostatic interactions had a small share in bonding Eu(III) at the high initial concentration.
The PSO rate constant strongly depended on the initial Eu(III) concentration (Fig. 2a, b).Adsorption rate constant increased with decreasing Eu(III) concentration because there was a larger fraction of binding sites in the adsorbent available for the Eu(III) ions.The higher the initial Eu(III) concentration, the longer it takes to reach equilibrium (Plazinski et al. 2009).
The observed effect of the initial Eu(III) concentration on k 2 was approximately 2 orders of magnitude higher at a lower L/S (Table 2).According to Al-Degs et al. (2006), an inverse relationship between initial Eu(III) concentration and k 2 in adsorption by a natural adsorbent is due to the effect of increased competition for the limited number of active adsorption sites by the adsorbate at higher initial concentrations of heavy metals.The observed decrease in k 2 with increasing Eu(III) concentrations can be explained by the occurrence of strong and weak adsorption sites.At high concentrations, Eu(III) ions tend to preferentially occupy strong binding sites and then weaker ones in clay minerals of the Phy/B mixtures.
The pseudo-second-order kinetic model assumes that the adsorption rate is controlled by chemisorption (He et al 2020), whereas the adsorption capacity is controlled by the number of active adsorption sites (Robati 2013).

Removal efficacy of Eu(III) by the Phy/B mixtures
The maximum initial concentration, C 0max , depended on the B content in the Phy/B mixtures (Fig. 3a).The highest impact of the B content on C 0max was observed for the 50Phy/50B mixture.The highest C 0max and the most favorable adsorption conditions were at pH 7 and L:S 100:1.No effect of the initial solution pH on C 0max was observed at L:S 500:1.The full adsorption of the 50Phy/50B mixture was close to the full adsorption of B at L:S of 500:1 for both pH values, and it was 0.68 and 0.79 of B adsorption capacity at L:S of 100:1 and pH 4.5 and 7.0, respectively (Fig. 3b).
Changes in RE of Eu(III) with increasing Eu(III) initial concentrations ranging from 10 to 200 mg/L were investigated for four pH-L/S systems and six concentrations (Fig. 4).Removal efficacy decreased as the initial Eu concentration increased, and it depended not only on Phy/B but also on L/S and the initial solution pH.
For the Eu(III) initial concentrations of 10-50 mg/L at pH 7.0 and L/S 100:1, Eu(III) ions were removed from the solution by the adsorption on the Phy/B mixtures as efficiently as on B. A stronger impact of pH and L/S on RE was observed at lower contents of B in the mixtures, whereas a higher impact of L/S than pH was observed in all Phy/B mixtures.The effect of pH was usually higher at L/S 100:1 than at 500:1 (Fig. 4).
Multi-way analysis of variance was used to study the effect of pH, L/S ratio, proportions of Phy to B in the mixtures, and Eu(III) initial concentration on adsorption efficacy of the Phy/B mixtures (Supporting Information).

Adsorption capacity of the Phy/B mixtures and mechanisms of Eu(III) bonding
The adsorption capacity of the Phy/B mixtures for Eu(III) increased with the increase in B content depending on the L/S and pH (Fig. 5).The adsorption system with the highest Eu(III) initial concentration did not reach the complete saturation for the Phy/B mixtures of 50/50 and 25/75 at L/S 100:1, regardless of pH as indicated by the lack of plateau on the adsorption curves (Fig. 5(a, b, c)).The adsorption capacities of the Phy/B mixtures for Eu(III) expressed in cmol + /kg ranged from 0.51 to 0.65 of their CECs.High adsorbent concentration at L:S of 100:1 and a high number of total adsorption sites resulted in the partial occupation of the adsorption centers at C 0 of 200 mg/L because of the high cation exchange capacity of montmorillonite in B (Table 1).The 75Phy/25B mixtures adsorbed 6.6 mg/g of Eu(III) at pH 4.5 and 12.14 mg/g at pH 7.0, i.e., 0.60 and 1.1 of CEC, respectively.High concentrations of the Phy/B mixture in solution enhance the possibility of collision between the suspended solid particles of the adsorbents leading to their aggregation.The aggregated particles have lower total surface area which results in a decrease in their adsorption capacity (Yu et al. 2015).
With an increase in the L/S to 500:1, i.e., a decrease in the adsorbent concentration, the adsorption capacity of the Phy/B mixtures for Eu(III) ions increased (Fig. 5(c,  d)).
The maximum adsorption capacity of all Phy/B mixtures was higher than their CEC and increased 1.19-1.62-foldat the solution initial pH 4.5 and 1.7-2.4-foldat pH 7.0.Therefore, a mechanism other than ion exchange must have been responsible for the Eu-binding.
The adsorption capacity of Phy and Phy/B mixtures was higher at the solution initial pH of 7.0 than 4.5, probably due to the precipitation of Eu(OH) 3 at pH 7.0, both on the surface of Phy/B mixtures and Phy (Fig. 1).The comparison of the maximum adsorption capacity of different Phy/B mixtures reveals the highest increase in adsorption capacity for 75Phy/25B relative to Phy (1.33-2.26-foldincrease) and for 50Phy/50B relative to 75Phy/25B (1.53-2.0-foldincrease) (Fig. 6).
Speciation of Eu(III), surface properties of the mixtures constituents, and protonation/deprotonation reactions of amphoteric surface functional groups (Si-OH, Al 2 -OH) in clay minerals were all affected by solution pH.
(3) (4) Europium may occur in solution at different valencies: Eu 3+ and EuCl 2+ at pH < 5 and Eu 2+ in hydrolyzed species, i.e., Eu(OH) 2+ and Eu(OH) 2 + at pH > 5.5 (Fig. 7).According to Bradbury and Baeyens (2002), hydrolysis of Eu does not occur until pH ~ 6 because of the large Eu ionic radius (0.95 Å).The presence of Eu(OH) 2+  and Eu(OH) 2 + ions in solution suggests that they can be bound by hydrogen bonds to the Si-OH and/ or Al 2 -OH surface functional groups of clay minerals in the Phy/B mixtures according to the reactions: The hydrolyzed species are more hydrophobic than trivalent cations.From the Fig. 7, it can be inferred that Eu(OH) 3 precipitates at pH > 6.0.However, high adsorbent concentration in the solution with the L:S ratio of 100:1 and high initial concentration of Eu(III) (100-200 mg/L) enables precipitation of Eu(OH) 3 even at pH ≤ 7.0 as it was the case in this study (Fig. 1).
The equilibrium pH of the Ph/B mixtures suspension decreased as the concentration of Eu(III) in the initial solution increased (Fig. 8) and ranged between 6.89 and 5.25 for L/S 100:1 at initial pH 4.5 and 7.0 and L/S 500:1 at pH 4.5.The equilibrium pH ranged from 7.53 to 5.87 at L/S 500:1 and pH 7.0 (Fig. 8).Equilibrium pH higher than the initial one results in good buffer properties of bentonite and removal of Ca 2+ and Mg 2+ ions from montmorillonite (Aitken et al. 1990).
The surface functional groups of minerals are protonated/ deprotonated depending on the pH, according to the following equations: At the initial Eu(III) concentrations between 0.01 and 35 mg/L, the equilibrium solution pH > pH PZC (Table 1), resulting in a negative surface charge of the mixtures-constituting minerals.Europium ions can be electrostatically bound to the negatively charged mineral surfaces according to the reaction: At C 0 > 35 mg/L, the surface charge of the Phy/B mixtures was positive due to the protonation reaction of the surface hydroxyl groups; hence, repulsive forces occur between Eu(III) and mineral surfaces.
As a consequence of the layered structure of clay minerals in the Phy/B mixture, Eu(III) ions are adsorbed by: (5) 1. Forming the outer-sphere complexes via cation exchange within the interlayer space in montmorillonite (Liu et al. 2023) and by the attachment to the permanently negatively charged sites on external basal surfaces of montmorillonite at low pH: 2. The formation of the Eu hydroxyl surface inner-sphere complex composed of the cation and silanol (Si-OH) or aluminol (Al 2 -OH) groups at the crystal edges in both bentonite and phyllite: The inner-sphere complexation at the edges occurs with increased pH (Geckeis et al. 2013).
According to the two-site protolysis non-electrostatic surface complexation and cation exchange (2SPNE SC/ CE) model of Bradbury and Baeyens (2002), the Ln/An uptake by montmorillonite can be explained by (i) cation exchange on planar sites, (ii) surface complexation on strong sites (≡S S OH) with high affinity, but low capacity, and (iii) surface complexation on weak sites (≡S W1,2 OH) with high capacity and low affinity.According to Gao et al. (2021), Am(III) adsorption on Gaomiaozi bentonite may be described by the surface complexation model including cation exchange species ((≡X) 3 Am) and three surface complexes (≡S S OAm 2+ , ≡S S OAm(OH) + , and ≡S S OAm(OH) 2 ).Fernandes et al. (2016) observed that Ln/ An 3+ strong inner-sphere complexes were formed at low loadings through binding to three Al(O,OH) 6 octahedra, most likely by occupying vacant sites in the octahedral layers of montmorillonite, which are exposed on {010} and {110} edge faces.At higher loadings, Ln/An 3+ were bonded only via one Al octahedron, forming a weaker, edge-sharing surface complex.
Iron oxyhydroxides in phyllite are additional active centers for the binding of Eu(III).
This study revealed the dependence of the edge site adsorption capacity of the Phy/B mixtures on the surface area.The 25Phy/75B mixture had the highest specific surface area and the highest adsorption capacity.Moreover, adsorption capacity increased with increasing CEC of the Phy/B mixtures.The contribution of individual mechanisms to the binding of Eu(III) by the Phy/B mixtures depends on the amount of bentonite. (9)

Adsorption isotherms
The parameters of isotherm models determined by the nonlinear regression analysis, values of determination coefficient, and non-linear error functions are summarized in Table 3, and the fits of experimental data with four models are shown in Fig. 9.The isotherm equations described the Eu(III) adsorption on the Phy/B mixtures very well, as suggested by the regression coefficients (R 2 ) ranging from 0.9161 to 0.9877 for the Freundlich model and from 0.9748 to 0.9962 for the Sips model (Table 3).Values of 1/n estimated from the Freundlich isotherm (Eq.(S1) -Table S1) were below 1 (0.1991-0.3281) and did not depend on experimental conditions.They were indicative of chemical adsorption of Eu(III) ions on all mixtures.The good adsorption capacity of the Phy/B mixtures was indicated by n > 2. Values of K F were higher for mixtures with a high content of bentonite and indicated a higher adsorption intensity at pH 7.0 and L/S 500:1 than at pH 4.5 and L/S 100:1.The fitting of adsorption data to the Freundlich isotherm assumes surface heterogeneity and confirms the high porosity of the Phy/B mixture suggested by the heterogeneous adsorption not limited to the monolayer adsorption (Sharma and Tomar 2011).K F ((mg/g)•(L/mg) 1/n ); q L (mg/g); K L (L/mg); q D (mmol/g); β (mol 2 /J 2 ); E (kJ/mol); q S (mg/g); K S (L/mg) m Results of the error analysis showed that the Langmuir model better simulated the adsorption isotherms than the Freundlich model (Table 3) (He et al 2020).The Langmuir constant, K L (Eq. (S2) -Table S1) varied depending on the adsorbent surface chemical composition and experimental conditions and was higher for the Phy/B mixture of 25/75 than 75/25.In addition, the effect of the experimental conditions on the affinity of the Phy/B mixture surface for Eu(III) and the bonding energy was similar to that of the K F values, i.e., it increased with the increase in pH and L:S.The values of the maximum adsorption capacity (q L ) were close to or higher than the theoretical (q max ) values for all Eu(III)-Phy/B mixture systems.This observation suggests a monolayer Eu(III) adsorption model associated with ion exchange, in accord with previously obtained results for minerals used for Eu(III) removal from solutions (Hu et al 2010;Kautenburger and Beck 2010;Songsheng et al 2012).
Mean free energy values, E (Eq.(S4)-Table S1), for the Dubinin-Radushkevich isotherm was above 16 kJ/mol ranging from 16.07 to 19.62 kJ/mol as a result of a strong covalent interaction between Eu(III) and the silanol -Si-OH or aluminol -Al-OH groups on bentonite and phyllite surfaces.
The Sips model is a combination of two other models and represents systems in which one adsorbed molecule can occupy more than one adsorption site (Eq.(S5) -Table S1).The Sips model improved the fit of experimental data in the zone of higher curvature of the adsorption isotherms (Fig. 9).Using a third parameter, it logically improves the quality of the mathematical fit.The maximum adsorption capacity values (q S ) for Eu(III) bonding obtained from the Sips model were higher than those derived from the Langmuir isotherm and than the experimentally determined.The parameter K S changed in the same manner as the constants K L of the Langmuir model.Differences between experimental and calculated values of adsorption capacity using the Langmuir and Freundlich models showed that at the lower initial concentration, the Freundlich model was a better choice for the characterization of adsorption but at higher initial concentrations; adsorption followed the Langmuir model.
The initial slope of the adsorption isotherm curve provides important information.A curve with a steep initial slope indicates that the adsorbent has a high affinity for Eu(III), and experimental conditions are most favorable.This affinity is indicated by the affinity constant derived from the Sips model; the highest affinity constant (K S = 1.563) was obtained for the adsorption of Eu(III) on the 25Phy/75B mixture for the L:S 100:1 and pH 7.0.

Distribution coefficient values for Eu(III) adsorbed on Phy/B mixtures
The laboratory batch sorption equilibrium test allowed the estimation of Eu(III) distribution coefficients (K d ).This is an important parameter in transport modelling and one of the key parameters used in the safety analysis of spent fuel repositories (Lehto et al 2019).Distribution coefficients are relevant to the specific geochemical conditions of each system under study and give an indication of the extent of Eu(III) partitioning to the solid phase in different systems and are useful in assessing the migration of Eu(III) in the barrier.It can be defined as the ratio of the solid-phase concentration (q) to the aqueous phase concentration at equilibrium and can be expressed as: If the adsorption is described by a non-linear isotherm, the partition coefficient is given by the equation: and for the Sips isotherm, which best describes the adsorption of Eu(III) ions by the equation: The distribution coefficient was calculated for the Phy50/ B50 mixture for the following conditions: initial concentration of Eu(III) 0.01-200 mg/L, L/S 100:1 and 500:1, initial concentration of Eu(III) 0.01-200 mg/L, pH 4.5 and 7.0, and is presented in Table 4, and the plots of the logK d values for The distribution coefficient values for Eu(III) are strongly dependent on the initial Eu(II) concentration, the pH of the solution, the L/S ratio, and the mineral composition of the solid phase.
The results show that the distribution coefficient values was: -Higher for bentonite than for phyllite under all conditions and increased with the percentage of bentonite in the mixtures -Higher for L/S of 100:1 than 500:1 for any initial concentration of Eu(III) ions -Higher for pH 7.0 than 4.5 -Decreased with increasing initial Eu concentration The lowest K d value was observed for Eu(III) adsorption on Phy, and it increased with increasing bentonite in the mixture.Furthermore, as shown in Fig. 10, there were greater differences in K d values for L/S 500:1 and 100:1 at pH 7.0 than at pH 4.5, and at pH 4.5 and 7.0 for L/S 100:1 than for 500:1.The results for Eu(III) adsorbed on bentonite at L:S 500:1, pH 4.5, and an initial concentration of 40 mg/L indicated that the K d value was twice as high as the K d value obtained by Pshinko et al. (2004) for Eu adsorption on Cherkasky montmorillonite under similar conditions.

Fig. 1
Fig. 1 BSE images and EDS spectra of Eu precipitates on: a bentonite; b biotite partially replaced by Fe-chlorite in phyllite; close-ups (c, d) show platy crystals of Eu-phase with rounded edges; EDS spectra of the Eu-phase (e, f).Peaks of Al and Si are from the background minerals and C from the carboncoating

Fig. 2
Fig. 2 Isothermal (23 ± 2 °C) time-dependent Eu(III) adsorption onto Phy/B mixtures for different L/S ratio at solution pH 7.0 (a) and at different initial solution pH for L:S 100:1 (b); exp, experimental data

Fig. 3
Fig. 3 Maximum initial Eu(III) concentration in solution, C 0max (a) and full adsorption of Eu(III) ions by the Phy/B mixtures (b) for the total removal efficacy RE > 99.99% and C eq < 0.0001 mg/L

Fig. 4
Fig. 4 Removal efficacy of Eu(III) by the Phy/B mixtures as a function of initial Eu(III) concentrations

Fig. 5 Fig. 6
Fig. 5 The adsorption capacity of the Phy/B mixtures for Eu(III) ions

Fig. 7 Fig. 8
Fig. 7 Distribution of hydrolysed products of Eu(III) at 25 °C in aqueous solution

Fig. 9
Fig. 9 Experimental data and predicted adsorption isotherms of Eu(III) for the Phy/B mixture of 75/25 Eu(III)  on Phy/B of 50/50 against the initial concentration are shown in Fig.10.

1.
Batch experiments showed a rapid and efficient Eu(III) adsorption on all of the Phy/B mixtures investigated.The adsorption rate constant was strongly dependent on the initial Eu(III) concentration, pH, and L/S ratio.A high correlation of pseudo-second order kinetic model with experimental data showed that the adsorption rate was controlled by chemisorption, while the adsorption capacity was controlled by the number of active adsorption sites.2. The adsorption capacity of the Phy/B mixtures for Eu(III) increased with the increase of the bentonite content in the Phy/B mixtures, but was dependent on both L/S and solution pH.The addition of bentonite to the mixtures had a strong effect on their CEC and specific surface area.3. The highest increase in adsorption capacity of Phy/B mixtures compared to phyllite was observed for phyllite/ bentonite ratios of 25/75 and 50/50.A greater effect of L/S than of solution pH on both the adsorption capacity and removal efficacy of Eu(III) was observed.4. Europium can be adsorbed by (i) the outer-sphere complexes via cation exchange in the interlayer space of montmorillonite; (ii) attachment to the permanently negatively charged sites on the outer basal surfaces of montmorillonite; (iii) the inner-sphere complex of the hydroxyl surface composed of the cation and silanol (Si-OH) or aluminol (Al 2 -OH) groups at the crystal edges in both bentonite and phyllite; (iv) hydrogen bonding and attraction to a negatively charged mineral surface; and (v) precipitation as Eu(OH) 3 .5. The Sips model showed the best fit to the experimental data, particularly in the zone of higher curvature on the adsorption isotherms indicating that adsorbed Eu(III) can occupy more than one adsorption site.6.The results of this study show the advantage of the Phy/B mixtures in immobilizing actinides and certain fission products by combining the adsorption properties of montmorillonite and chlorite.

Table 1
Physicochemical properties of phyllite (Phy), bentonite (B), and the Phy/B mixtures While in most disposal scenarios concentrations of released radionuclides are much lower than considered in this study (in the low ppm to ppb range), the wide range of initial Eu(III) concentrations used in the adsorption experiments allowed for the more precise measurements.

Table 3
The isotherm constants for the adsorption of Eu(III) ions onto Phy/B mixtures