Insight into Adsorption Kinetics of Cs + , Rb + , Co 2+ , and Sr 2+ on a Zeolites-Based Composite: Comprehensive Diffusional Explanation and Modelling

: Kaolinite-rich soils were used to prepare zeolite-based composites via alkaline activation. The porous material was characterized by conducting XRD and microporosity measurements, as well as ESEM microscopy. The Weber and Morris (W-M) model was used for studying adsorption kinetics of radioactive cations on synthesized alkali-activated material. These investigations evidenced the effects of pore structure and the importance of the intrinsic characteristics of hydrated cations (ionic potential; hydrated radius; B-viscosity parameter; molar Gibbs energy of hydration of cation) on W-M kinetic rate constants. The application of diffusion-based models permitted us to assess the key diffusion parameters controlling successive diffusion regimes, and to reveal strong contributions of surface diffusion to adsorption kinetics during the course of the second and third kinetics stages of the W-M model. The magnitude of the surface diffusion coefficient was related to the capacity of hydrated cationic species to lose water molecules when penetrating brick pores. The HSDM model were tested for predicting radionuclide adsorption in a fixed-bed column. A breakthrough curve simulation indicated the predominance of the surface diffusion regime, which was in agreement with mathematical analysis of (batch) adsorption kinetics data. Ionic diffusion was linked to the characteristics of capillary porosity and connectivity of capillary pores in the composite, suggesting the generation of hydrated nuclides and their immobilization in the form of outer-sphere complexes.


Introduction
In order to decrease greenhouse gas emissions, it is becoming more and more imperative (i) to preserve/control natural resources by attempting to reuse/recycle them; (ii) to decrease industrial wastes/by-products by transforming them into innovative materials for new applications; (iii) and above all, to limit, at best, energy consumption derived from fossil combustibles.It has, therefore, been necessary to rapidly develop more sustainable alternatives to certain industrial products, like cement materials, which are currently recognized to be responsible for about 8% of anthropogenic CO 2 emissions worldwide [1].In this view, aluminosilicate geopolymer materials have attracted research on their use not only as construction materials, but also as ion exchangers in effective water remediation processes [2,3], since their synthesis is easy and low-cost, necessitating natural and available (waste) materials.They are obtained through alkali activation: (i) either directly from natural minerals (such as kaolin [3] and metakaolin [4][5][6]), or (ii) from industrial products/wastes (such as fly ash and blast furnace slag [7][8][9], red mud [10], and brick [11,12]).
Structurally, aluminosilicate geopolymers are amorphous to semi-crystalline materials.Their three-dimensional networks constitute cross-linked tetrahedral silica and alumina units, forming rings comparable to those in zeolites.These rings contain "negatively charged" aluminium ([AlO4] 5− ) charge-balanced by cations such as Na + or K + , facilitating the ion exchange with "positively charged" counter-ions (such as hazardous/toxic heavy metals) in the aqueous phase through electrostatic attraction, and thereby enabling metal(s) immobilization.Furthermore, aluminosilicate geopolymer materials provide relatively good mechanical properties and thermal stability [33,34], but they possess a weak surface area and a low number of active surface sites in comparison with crystalline aluminosilicates such as zeolites.Some of the zeolite types have great importance in heavy metal removal [35], dye removal [36], or even in agricultural area like soil amending [37].In pure zeolites like NaA and NaP, their frameworks contain the highest atomic Al/Si ratio and thereby possess the highest negative-charge density.However, because of their small pore sizes, the ionexchange kinetics of radioactive cations with large diameters (like Cs + and Sr 2+ , with 4.5 Å and 8.2 Å, respectively) take place more slowly on zeolites [38].In this context, as natural zeolites possess micropores with weak interconnection, their surfaces could be easily obstructed by adsorbed (ionic, gaseous or neutral) species, and hence become inactive with time when they are used as catalysts, ion-exchangers, or adsorbents.Conversely, the use of hierarchical structuring of zeolite or Al-rich zeolites with larger pores sizes should allow for the improvement of "large" cation/Na + exchanges and adsorptive capacity [38][39][40][41].These are obtained with zeolites synthesized via conversion of geopolymers [42][43][44].Indeed, the use of geopolymeric zeolites has been found to be more adapted for practical applications because they possess pores with larger sizes and that are interconnected.These pore characteristics confer geopolyparticularly better adsorptive properties towards larger adsorbates to geopolymeric zeolites, and more generally, they enhance the adsorption of toxic materials such as heavy metal ions, dyes, herbicides, and drugs [42][43][44].
Currently, more emphasis is being given to the synthesis of zeolitic materials from low-cost industrial wastes with the aim to reduce preparation costs, as well as to valorise and recycle industrial by-products [45][46][47][48][49].Among these, brick was recently proposed by our research group as an interesting (green) material, mostly owing to the following reasons: (i) the facile and low-cost production of zeolitic/-geopolymic aluminosilicate brick composites; and (ii) the introduction of meso-porosity in generated zeolites NaA and NaP during chemical activation, enabling better ion-exchange kinetics and higher adsorption capacity toward hydrated (radioactive) cations with large sizes [12,50].In addition, for industrial applications, it is imperative that adsorptive mixtures like "zeolites + geopolymers" be supported and shaped in sufficiently large particles in order to avoid high pressure drops and weak permeability during the filtration procedure.In order to overcome this issue, the use of brick presents another notable advantage: the possibility of grinding, granulating, and sieving this material into suitable particle sizes.The technical capacity for producing approximately cylindrical granules with adequate sizes would facilitate the implementation of large-scale filtration systems in which the inlet liquid would flow more easily through adsorbent-grain interstices (thereby ensuring low pressure drops and better permeability) [11,12].In the lab, brick aggregates were activated with sodium hydroxide with the aim of generating a low-cost, mesopore-rich composite according to a one-step process.Preliminary investigations revealed that the obtained material was a good adsorbent for the removal and immobilization of the Cs + and Rb + radionuclides from water, and was well-adapted for column experiments and future large-scale (industrial) applications owing to the facile granulation of the material [11,12].In the present work, caesium, strontium, and cobalt were chosen because their isotopes, 137 Cs, 90 Sr and 60 Co, are among the most frequent radionuclides encountered in liquid radioactive wastes that can lead to dangerous releases of radioactivity to the environment following nuclear incidents or accidents.These radioactive metals are known to be relatively stable at the human scale, and decay at half-life times of a few years: 60 Co, 5.27 years [51]; 137 Cs, 30.17 years [52]; and 90 Sr, 28.7 years [52].Rubidium was also studied in this article because it was recognized first to have great potential to trace various geological processes [53], and second to be useful as a conservative tracer for effluent pollution/impact diagnostics of municipal wastewater by considering the ratio of Rb concentration to that of strontium in highly urbanized, anthropogenically impacted freshwater systems [54].The Rb element possesses two naturally occurring isotopes: the stable isotope 85 Rb (72.17%) and the radioactive isotope 87 Rb (27.83%).The 87 Rb isotope is very stable, with a half-life of 4.97 × 10 10 years, and it decays into 87 Sr [55].If the isotopes 137 Cs, 87 Rb, 60 Co, and/or 90 Sr were to contaminate soils and water streams, they would be capable, with time, of passing through food chains, and thus would pose relevant risks to human beings.Furthermore, as these radioactive elements are high-energy beta-ray or gamma-ray emitters, it is urgent to remove them from radioactive wastewaters using an effective material and an efficient method [56].
The first objective of this research was to characterize the brick-based composite by conducting X-ray diffraction (XRD) and micro-porosity measurements, as well as environmental scanning electron microscopy (ESEM).
In the lab, the aqueous adsorptive characteristics of the brick composite towards some cations (Cs + , Rb + , Co 2+ , and Sr 2+ ) were examined, with special emphasis on the implication of surface properties on ionic motion within composite frameworks as well as on adsorption advancement.For this purpose, it was necessary to estimate preliminary equilibrium adsorption data and adsorption kinetics parameters: (i) for providing insight into the underlying mechanism; (ii) for predicting the reaction time to accomplish adequately specific radioactive wastewater treatment; and (iii) for optimizing the design of (future) industrial batch/column adsorption systems.These objectives could be accomplished by imperatively identifying and modelling heterogeneous reactions involved in the the brick composite/water interface.However, by applying pseudo-first-and -second-order models, the adsorption kinetics process was over-simplified because intra-particle diffusion and external mass transfer were not taken into consideration in the adsorption kinetics.The knowledge and evidencing of intra-particle diffusion and mass-transfer phenomena at the water-solid interface were essential for us to better understand the heterogeneous kinetics process being studied.Thereby, the second objective of this research was to apply to our system the kinetics-based model proposed by Weber and Morris [57] to reveal reaction pathways and rate-controlling steps.In this step, the effects of pore distribution and the importance of intrinsic properties of hydrated cations on diffusion kinetics were also studied.
The Weber-Morris model is, nevertheless, considered more as a descriptive model than a predictive/theoretical model for quantifying diffusion/mass transfer phenomena.The third objective was, therefore, to use diffusion-based models, enabling the estimation of diffusion or mass transfer parameters, which are essential for designing adsorption processes under static (batch) and dynamic (fixed-bed column) conditions.In these modelling approaches, adsorption kinetics were investigated by assuming the occurrence of the following sequential steps: (i) film diffusion/mass transfer; (ii) intra-particle diffusion/mass transfer; and (iii) adsorption on the pore surface of the adsorbent (surface reaction).The corresponding mathematical equations were applied to our heterogeneous system for the purpose of modelling/interpreting the adsorption kinetics process and evidencing/quantifying rate-limiting steps.
Should obtained (batch) information about diffusion kinetics and thermodynamic (interfacial) equilibria be relevant and valuable for predicting the removal performance of the brick composite towards Cs + , Rb + , Co 2+ , and Sr 2+ cations in fixed-bed columns?To verify this, our fourth objective was to use a mass-transfer model to describe the dynamics of the cation/Na + -exchange process in packed-bed columns.Breakthrough curves were experimentally determined and modelled using the homogeneous surface diffusion model (HSDM).The FAST 2.1 program [58] was used for simulating breakthrough curves and solving HSDM calculations.The involved mass-transfer mechanism was investigated in detail by modelling the adsorptive and diffusional behaviour of the elements (Cs, Rb, Sr, and Co) through the porous frameworks of the brick composite.In this part, the dimensionless analysis of the involved diffusional regimes and the implication of the connectivity of pores on ion diffusivity were also examined.
The originality of this work was that it rationally interpreted the effects of hierarchically mesoporous "zeolitic" structures appearing at the surface of an activated brick on the adsorption kinetics of Cs + , Rb + , Co 2+ , and Sr 2+ cations, which were studied individually, and to distinguish/quantify different diffusion regimes taking place at the solid-water interface under static and dynamic operational conditions.

Preparation of Zeolite-Based Composite
Kaolinite-rich soils from Bangui (Central African Republic) were used for making brick.Briefly, the extracted soils were mixed with water, and the obtained mud was shaped manually.The resulting air-dried (48 h) bricks were placed in efficient stackings with air flows in order to constitute a basic oven (that was built simply on ground).They were then heat-treated with dry wood for a period of about three days at temperatures ranging from 500 • C to 800 • C, and finally cooled progressively up to ambient temperature for two/three days.In order to increase the surface area of the brick material, it was broken into grains manually using a hammer.Brick particles were sieved with mechanical sieves afterwards, and the fraction of particles with sizes varying from 0.7 to 1.0 mm was kept for our experiments.This fraction was washed with Milli-Q water and then decanted; after settling, the water was eliminated and the brick grains were dried at 105 • C.This brick contained mainly quartz (60-65 w%) and metakaolinite (20-26 w%), and, to a lesser extent, iron oxide/hydroxide, illite, and titanium dioxide [59].Afterwards, brick pellets were treated with sodium hydroxide according to the following optimized synthesis conditions [60]: 10 g of Bangui brick reacted in 40 mL of a diluted NaOH solution (0.6 mol•L −1 ) at room temperature for 24 h under slow shaking at a speed of 120 rpm.This procedure was followed by a fixed-temperature increase in the mixture at 90 • C for a constant reaction time of six days.Finally, the recovered grains were rinsed several times with MilliQ water and dried at 90 • C for 24 h.

Chemicals
All chemicals employed in the experiments were of analytical grade.Sodium hydroxide, Co(NO 3 ) 2 .6H 2 O, Sr(NO 3 ) 2 , Cs(NO 3 ), and Rb(NO 3 ) were supplied by DISLAB (Paris, France).It is worth noting that only salts of non-radioactive isotopes were used as model substances for the kinetics experiments.

ICP-AES Analyses
During batch and column studies, the recovered solutions were analysed for most of their element contents using an inductively coupled plasma optical emission spectrometer (ICP-OES, model 5110 VDV, Agilent Technologies, Paris, France).As for caesium, because of the low emission intensity detected by ICP-OES, Cs contents in recovered solutions were determined by means of an ICP-MS spectrometer (Agilent Technologies model 7850).

Kinetics (Batch) Experiments
All kinetics experiments, excluding the pH effect tests, were carried out at the natural pH of the suspension initially measured, i.e., 7.74 for Cs + ; 7.42 for Rb + ; 7.37 for Sr 2+ ; and 7.29 for Co 2+ .These experiments were conducted at 298 K.For this, with a a thermostat, we used a double-walled cell into which water maintained at a constant temperature flowed (model Thermofisher Polystat 71; Paris, France).A temperature sensor (Pt 1000 Metrohm, Paris, France) and a pH Methrom electrode (Metrohm 6.0262.100,Paris, France) were employed to accurately measure the temperature and pH of the suspension.These sensors were connected to an automatic titrator (model Metrohm Titrando 888, Paris, France), and the measuring device was controlled by the software Tiamo (version 2.3).
Kinetics tests were performed at the minimum shaking speed above which the kinetics were independent of the degree of agitation, and therefore not influenced by film diffusion.The synthetic solutions which were employed for kinetics experiments were prepared by diluting 1000 mg•L −1 stock solutions of cations with distilled water.The course of the kinetics process was appraised by following the decrease in adsorbate with time in the liquid phase.For this, an element solution of 100 mL (2 mmol•L −1 ) was put into contact with 0.5 g of alkali-activated brick in a 200 mL glass reactor until an equilibrium state was attained (~240 mn).Then, 0.5 mL of suspension solution was collected at different time intervals (ranged from 0.5 to 240 mn) and diluted with distilled water with the aim to determine the element concentration using either ICP-OES or ICP-MS.
It is worth noting that the effects of element concentration (in the aqueous phase) on cation adsorption (on the brick) ought to be attributable preliminarily to the number of active sites of alkali-activated brick interacting with cationic species (i.e., the mass of brick grains in contact with solution).However, the extent of cationic removal ought also to be related to the affinity, reactivity, and diffusional/complexing characters of the adsorbate towards the brick surface.
Note that the adsorptive behaviour of raw brick was studied in the lab, and data were afterwards compared to those obtained with brick activated either by acidic or basic treatment.Our investigations clearly showed the weak adsorption capacity of raw brick towards cations such as Fe 2+ and Pb 2+ because of the absence or negligible presence of negatively charged sites associated with extra-framework ions [59].

Fixed-Bed Column Experiments
The continuous flow adsorption experiment was conducted in a fixed-bed glass column with an inner diameter of 12.5 mm, a height of 25 cm, and a porosity of 160-250 µm, with a sintered pyrex disk at its bottom to prevent any loss of material.A bed depth of 8.6 cm (10 g) was investigated at a constant flow rate of 5 mL/min.Before being used in the experiments, at least 500 mL of Milli-Q water was passed through the column.The schematic diagram of the fixed-bed column reactor used was previously illustrated [61].
In a single-component system (Co 2+ , Sr 2+ , Cs + , and Rb + ), the concentration of cations in the influent was fixed at about 1 mmol•L −1 .The cation solution was pumped through the column at the desired flow rate using a peristaltic pump (LaboModerne France Type KD1170, Paris, France) in up-flow mode.During the column experiment, the effluent pH was measured continuously using a pH meter that was connected to a computer, and effluent samples exiting the bottom of the column were collected at different time intervals and analysed for element contents using ICP-OES or ICP-MS.The flow to the column continued until the effluent metal concentration at time t (C t ) reached the influent metal concentration (C o ): C t /C o ≈ 0.99.The performance of the packed bed is described in the present work using the concept of the breakthrough curve.

Electron Microscopy
Micrographs of representative specimens of the brick composite were recorded using an environmental scanning electron microscope equipped with an Energy-Dispersive X-Ray Spectrometer EDS X flash 3001 (Thermo Fisher Scientific, Courtaboeuf, France).EDS measurements were carried out at 20 kV at low vacuum (1.00 Torr), and the maximum pulse throughput was 20 kcps.Different surface areas ranging from 0.12 to 0.50 mm 2 were targeted on alkali-brick grains and examined by ESEM/EDS.For this, a narrow beam scanned selected areas of brick pellets for chemical analysis.Atomic quantifications and mathematical treatments were undertaken using QUANTA-400 software (elaborated by Bruker) in order to determine the averaged elemental composition of the surface brick and to detect chemical/elemental variability.

X-ray Diffraction (XRD)
XRD analysis was conducted at room temperature in a Bruker D8 Advance diffractometer (Evry, France) using Ni-filtered CuKα radiation (40 kV, 40 mA).Brick composite samples were scanned with a step size of 0.02 • and a counting time of 0.5 s per step.

Nitrogen Adsorption-Desorption Isotherms
The textural characterization of the activated-brick samples was determined from the N 2 adsorption desorption isotherm at 77 K using the micrometrics model Tristar II 3820.Before analysis, the brick material was heated at 373 K and degassed for 2 h under primary vacuum conditions (pressure: 10 Pa).
The specific superficial area was assessed using the Brunauer-Emmet-Teller (BET) method.The pore volume was estimated using the α-plot method.The total pore volume was evaluated from the desorption isotherm branch at P/Po = 0.98, assuming complete pore saturation.Pore size distribution was calculated using non-local density functional theory (NLDFT) and the Berrett-Joyner-Halenda (BJH) method.

Theoretical Approach of Adsorption Kinetics
Under static (batch) experimental conditions, the Weber-Morris model allowed us to gain evidence and follow up with time film mass transfer and pore/surface diffusion phenomena over the course of the kinetics process.Moreover, the combined use of the Weber-Morris model and diffusion-based models (as described in Section 3.3) allowed us to quantify the different diffusional regimes taking place on the brick surface during the adsorption process.
Pore surface diffusion modelling of the cations' breakthrough curves was performed using the dynamic packed-bed model HSDM (homogeneous surface diffusion model).Using this model, two subsequent steps were performed: (i) the mass transport from the solution bulk to the brick's external surface (which was represented here by a mass transport coefficient, k f ); and (ii) the intra-particle diffusion (which was attributable to both pore volume diffusion, D p , and surface diffusion, D s ).On the basis of one-dimensional solute transport in the fixed bed, a series of fundamental differential convective-diffusive equations for a packed-bed setup was established to quantify interfacial phenomena, and hence, to explain the adsorber dynamics.In the presently discussed work, HSDM simulations were undertaken by means of the FAST 2.1 software developed by Schimmelpfennig and Sperlich [58].
On the other hand, a dimensionless analysis of the HDSM model and its application for the simulation of breakthrough curves of modified brick was undertaken by considering the dimensionless parameters: Biot, Stanton, and solute distribution.
Adsorption kinetics models: Experimental data were described using the following adsorption kinetics equations (Weber-Morris adsorption kinetics model [62]): where Q t (in µmol•g −1 ) is the adsorption capacity at any time t; k ipd (in µmol•g −1 •mn −1 ) corresponds to the intra-particle diffusion rate constant; and C (in µmol•g −1 ) is a coefficient indicative of the film thickness (film-surface diffusion model [63]).
where r o represents the initial adsorption rate; k f is the film mass transfer coefficient; C o is the initial concentration of adsorbate; R b is the brick-grain radius; and ρ b is the brick density.Knowing the r o value, the k f coefficient can then be determined from the following equation: According to Yao and Chen [63], at the early stage of the adsorption process, kinetic data can be fitted with an equation equivalent to that used for pseudo-second-order kinetics: where "r o " and "a" represent two fitting parameters.Pore diffusion coefficient: The pore diffusion coefficient (D p ) was estimated from the ionic diffusivity of the radioactive cation in water (D) using the following equation [64]: where τ is the tortuosity factor, which is defined as: Note that in all kinetics experiments, we used nitrate salts.The ionic diffusivity of each cation could then be calculated from the following equation [65]: where D 1 • and D 2 • represent the infinite dilution diffusion coefficients for radioactive cations and anions (NO 3 − ), respectively.Dimensionless (fixed-bed column) parameters: The Biot number (B i ) corresponds to the ratio between the film mass-transfer rate and the intra-particle surface diffusion, and is defined as: The Stanton number (S t ) is the dimensionless liquid-phase mass transfer coefficient, which is defined as: where m is the mass of the adsorbent and F represents the flow rate.The S t parameter is considered as a good indicator of the influence of film diffusion on the shape of the breakthrough curve.As for the solute distribution parameter (D g ), it is defined as the ratio of cationic species adsorbed onto brick particles to those in the liquid phase: where ρ B and ε B are the bed density and bed porosity, respectively; C o is the influent concentration; and Q o is the equilibrium solid-phase concentration corresponding to the influent concentration, C o (Q o is determined using the Freundlich equation: where K F and n represent Freundlich isotherm coefficients).

XRD Analysis
The XRD patterns of raw brick and alkali-activated brick presented in Figure 1a,b revealed that quartz was the major crystalline phase, as were, to a lesser extent, illite and rutile.In the XRD patterns of alkali-activated brick (see Figure 1b), additional peaks with lower intensities appeared in the diffractogram.

XRD Analysis
The XRD patterns of raw brick and alkali-activated brick presented in Figure 1a,b revealed that quartz was the major crystalline phase, as were, to a lesser extent, illite and rutile.In the XRD patterns of alkali-activated brick (see Figure 1b), additional peaks with lower intensities appeared in the diffractogram.; LTA zeolite, .
It is worth noting that the micro-porosity technique used herein did permit us to assess the pore size contributions provided by LTA and Na-P zeolites to the overall surface area.Indeed, the dimension for nitrogen gas to penetrate the pores at 77 K is limited to 5 Å [67], whereas for pore sizes in the LTA zeolite, it is ≤4.1 Å [68], and dimensions of intersecting channels in the Na-P zeolite (i.e., 0.31 nm × 0.44 nm and 0.26 nm × 0.49 nm) are too weak [69][70][71].
Overall, XRD characterizations performed on alkali-activated brick revealed the appearance of new reflections, which were ascribed to two crystalline zeolitic phases: zeolite Na-A and zeolite Na-P.This study confirmed the occurrence of crystal nucleation of the amorphous geopolymeric (aluminosilicate) gel phase (knowing that this gel was intermediately generated through brick-metakaolinite dehydroxylation under fixed thermal/alkaline conditions).

ESEM/EDS Analysis
ESEM micrographs of the treated brick pellets are shown in Figure 2a-c.As can be seen in these ESEM micrographs, the morphological aspects of this material were evidenced at two magnifications.The ESEM micrograph of raw brick is displayed in Figure 2d.Elemental (ESEM/EDS) analysis of alkali-activated brick permitted us to detect the following elements: (i) silicon (Si), aluminum (Al), and oxygen (O) as major elements; and (ii) iron (Fe), titanium (Ti), magnesium (Mg), potassium (K), and calcium (Ca) as minor elements.In order to determine their atomic surface contents, several spherical and ellipsoidal surface zones were targeted and analyzed using ESM/ED, as shown in Figure 2e.The averaged atomic percentages of elements constituting the alkali-activated brick  134), respectively [66].
It is worth noting that the micro-porosity technique used herein did permit us to assess the pore size contributions provided by LTA and Na-P zeolites to the overall surface area.Indeed, the dimension for nitrogen gas to penetrate the pores at 77 K is limited to 5 Å [67], whereas for pore sizes in the LTA zeolite, it is ≤4.1 Å [68], and dimensions of intersecting channels in the Na-P zeolite (i.e., 0.31 nm × 0.44 nm and 0.26 nm × 0.49 nm) are too weak [69][70][71].
Overall, XRD characterizations performed on alkali-activated brick revealed the appearance of new reflections, which were ascribed to two crystalline zeolitic phases: zeolite Na-A and zeolite Na-P.This study confirmed the occurrence of crystal nucleation of the amorphous geopolymeric (aluminosilicate) gel phase (knowing that this gel was intermediately generated through brick-metakaolinite dehydroxylation under fixed thermal/alkaline conditions).

ESEM/EDS Analysis
ESEM micrographs of the treated brick pellets are shown in Figure 2a-c.As can be seen in these ESEM micrographs, the morphological aspects of this material were evidenced at two magnifications.The ESEM micrograph of raw brick is displayed in Figure 2d.Elemental (ESEM/EDS) analysis of alkali-activated brick permitted us to detect the following elements: (i) silicon (Si), aluminum (Al), and oxygen (O) as major elements; and (ii) iron (Fe), titanium (Ti), magnesium (Mg), potassium (K), and calcium (Ca) as minor elements.In order to determine their atomic surface contents, several spherical and ellipsoidal surface zones were targeted and analyzed using ESM/ED, as shown in Figure 2e.The averaged atomic percentages of elements constituting the alkali-activated brick surface were estimated: 19.00 ± 1.24 at.% for silicon; 10.44 ± 0.76 at.% for aluminium; 60.57 ± 1.22 at.% for oxygen; 0.92 ± 0.33 at.% for iron; 0.42 ± 0.36 at.% for titanium; 0.60 ± 0.03 at.% for magnesium; 0.97 ± 0.49 at.% for potassium; and 0.09 ± 0.02 at.% for calcium.Some of these elements were found to be better localized across the brick surfaces, evidencing the presence of SiO 2 (as quartz); alumino-silicates; and, to a lesser extent, iron oxides/hydroxides and titanium oxide [11,12].Aluminosilicates present in alkali-activated brick were detected as aggregates with either spherical or cubic shapes.Their sizes varied from 8 µm to 20 µm for spherical particles and from 15 µm to 29 µm for cubic particles, as shown in Figure 2f.
Detailed ESEM/EDS analysis was performed on several cubic and spherical specimens targeted hazardously on alkali-activated brick samples, and the atomic percentages of Al, Si, and Na were determined and averaged.As metakaolinite-derived zeolites ought to belong to the GIS family, the chemical formula (based on the crystalline unit cell) should correspond to the crystalline phases: , where m is the valence of alkali ions and n represents the number of water molecules per unit cell.Hence, it was possible to establish structural and chemical formulas and calculate standard deviations after averaging ESEM/EDS analysis data.We found: (Na) 13.86±1.25 (?) z (Al) 16.43±0.99 (Si) 15.89±0.98 (O)  32 •nH 2 O (for cubic particles) (Na) 12.57±1.05 (?)z(Al) 14.92±1.14 (Si) 15.58±1.45 (O)  In these formulas, the Si/Al atomic ratios are close enough to 1 (0.97 for cubic specimens and 1.04 for spherical specimens).Note that the atomic deficit for sodium in the formulas is represented by (?) z .This deficit was assumed to be compensated for by other cations, like K + , Ca 2+ , and/or Mg 2+ , which were detected at trace levels on zeolitic frameworks.
The spherically shaped specimens were morphologically similar to those reported in the literature for zeolite Na-P [72][73][74][75], while the cubic-shaped specimens were comparable with those observed previously for the A-type zeolite [76,77].Note the presence of zeolitic (Na-P) aggregates with "sponge" shapes (see Figure 2c).The occurrence of Na-A and Na-P "zeolitic" crystals in the brick composite was confirmed using XRD, as described in Section 3.1.1(see Figure 1).
Overall, ESEM/EDS characterizations performed on alkali-activated brick revealed the presence of numerous spherical and cubic (sponge) zeolitic crystalline specimens.Their specific macro-, meso-, and microporous characteristics and negatively charged surfaces make it possible to use this brick material as an adsorbent for the elimination of ionic pollutants.

Porosity and Specific Surface Area
To characterize the specific surface area and porosity of the brick-based composite, N 2 adsorption analysis was carried out.Figure 3 displays the N 2 adsorption-desorption isotherm and the corresponding pore size distribution curve for the studied material.
The obtained N 2 adsorption-desorption isotherm corresponded to a type-IV isothermal profile, according to the International Union of Pure and Applied Chemistry (IUPAC) classification [78].The adsorption of P/P 0 increased abruptly between 0 and 0.01 (Figure 3a), which confirmed the existence of micropores in the modified brick.The appearance of the hysteresis loop at the medium pressure area (0.2 < P/P 0 < 0.8; see Figure 3a) proved the abundance of mesopores, and the steep increase in the N 2 capacity in the high-pressure area (P/P 0 > 0.8; see Figure 3a) indicated the presence of macropores.The shape of this hysteresis loop resembled an H3-type hysteresis loop [78], indicating that the brick composite was mainly composed of mesopores and macropores [79,80].The pore size distribution shown in Figure 3(b1) agreed with the analysis described above, and revealed that the final composite had hierarchical pores with micropores, mesopores, and macropores.Compared to the pore distribution of raw brick (see Figure 3(b2)), the profile changes in the pore size distribution of the treated brick were intimately related to the significant increase in mesopores, which appeared mostly at low diameters, <50 Å, and reached a maximum pore volume for pore diameters around 28 ± 1 Å.New mesopores were generated through zeolitization of brick-metakaolinite by alkaline treatment.Disordered hollow-like mesoporous frameworks were created by aggregated zeolitic crystals.As for the macroporous structure of the brick (i.e., pores with diameters of about ≥500 Å), it could be directly observed from the ESEM images of a raw brick sample (see Figure 2d), and above all, on ESEM images of treated brick samples taken at higher magnifications (Figure 2b,c), particularly revealing numerous macropores on the surfaces of zeolitic aggregates (Na-P) with "sponge" morphologies (see Figure 2c).
Although the values of average pore diameter, BET surface area, and pore volume obtained for the activated brick were somewhat low, our recent research has revealed that the generated "hierarchical" mesoporous structure of activated brick was well suited for adsorbing inorganic pollutants into the interior space of the alkaline brick under large-scale dynamic conditions [12,81,82].Thus, our previous investigations, which were performed by our research group at Bangui University (Central African Republic), had permitted us to show that the brick-based composite could be an interesting adsorbent in The analysis of N 2 adsorption-desorption isotherms was undertaken in the following range of pore sizes: 17-3100 Å.This analysis permitted us to determine the average pore diameter (94.5-112.4Å), BET surface area (19.3 m 2 •g −1 ), and pore volume (0.036 cm 3 •g −1 ).Overall, the porosity analysis revealed the generation of hierarchical mesopores with "bimodal" pore size distribution at 29 Å and ~70-150 Å.
Compared to the pore distribution of raw brick (see Figure 3(b2)), the profile changes in the pore size distribution of the treated brick were intimately related to the significant increase in mesopores, which appeared mostly at low diameters, <50 Å, and reached a maximum pore volume for pore diameters around 28 ± 1 Å.New mesopores were generated through zeolitization of brick-metakaolinite by alkaline treatment.Disordered hollow-like mesoporous frameworks were created by aggregated zeolitic crystals.As for the macroporous structure of the brick (i.e., pores with diameters of about ≥500 Å), it could be directly observed from the ESEM images of a raw brick sample (see Figure 2d), and above all, on ESEM images of treated brick samples taken at higher magnifications (Figure 2b,c), particularly revealing numerous macropores on the surfaces of zeolitic aggregates (Na-P) with "sponge" morphologies (see Figure 2c).
Although the values of average pore diameter, BET surface area, and pore volume obtained for the activated brick were somewhat low, our recent research has revealed that the generated "hierarchical" mesoporous structure of activated brick was well suited for adsorbing inorganic pollutants into the interior space of the alkaline brick under large-scale dynamic conditions [12,81,82].Thus, our previous investigations, which were performed by our research group at Bangui University (Central African Republic), had permitted us to show that the brick-based composite could be an interesting adsorbent in the conception of a pilot associated with an iron(II)-contaminated drilling [83,84].Hence, by considering the macro-/mesoporous structure of brick grains and optimizing their sizes to improve flow rates, one could anticipate that such porous characteristics ought to be beneficial to enhance the adsorption efficiency of cationic adsorbates.

Weber-Morris (W-M) Model
In order to evidence and differentiate diffusion regions in adsorption kinetics, we decided to test Weber-Morris equation to qualitatively reveal the effects of the heterogeneous and porous nature of alkali-activated brick on the diffusivity of Cs + , Rb + , Co 2+ , and Sr 2+ ions inside adsorbent frameworks.
It has previously been stated that the Weber-Morris (W-M) model is useful for identifying reaction pathways and predicting rate-controlling steps in porous matrices [85].The model proposed by Weber and Morris [62] was applied in this work to gain deeper knowledge regarding the mechanism of nuclide adsorption onto the porous brick composite in aqueous solutions (see Section 2.9).For each element, the plot of Q t vs. t 1/2 was found to be divided into three stages by two break points (Figure 4).As suggested in previous works regarding adsorption kinetics in relation to the adsorbent pore structures and surface properties of adsorbents [86][87][88][89], it was assumed here that the multi-linearity observed for cation adsorption was due to the effects of the adsorbent porosity on cations.The observed multi-linear plot would then result from the broad pore size distribution of the brick composite, including, essentially, macropores and mesopores and, to a lesser extent, meso/micropores of "hierarchical" zeolites (Figure 3b).The higher value of k ipd (1) in the first stage of the Weber-Morris model compared to k ipd (2) and k ipd (3) revealed the importance of external liquid film diffusion (Table 1).
Table 1.Intra-particle diffusion rate constants (k ipd ) and film thickness (C) determined for the different stages observed in Weber-Morris curves.The first linear plot represented the fastest adsorption stage.In connection with the porous structure properties of the brick composite (as shown in Figure 3b), the liquid film diffusion occurred not only on the external brick surface, but also inside macropores and even wider mesopores.This indicates that, in the first W-M stage, the ions Cs + , Rb + , Co 2+ , and Sr 2+ were easily transported into the broader pores of the brick without significant hindrance.As for the second stage of the Weber-Morris model, the lower slopes of the second linear W-M plots (k ipd (2) < k ipd (1); see Table 1) showed that cationic species diffused more slowly into smaller brick pores.The higher correlation coefficient (R 2 ) obtained in the second W-M stage (see Figure 4) suggests, a priori, that intra-particle diffusion was the rate-limiting step of the adsorption process.But, because none of the second linear W-M plots drawn passed through the origin in Figure 4, it can be concluded that neither the liquid film diffusion nor the intra-particle diffusion was the only ratecontrolling step.Finally, in the third stage of the Weber-Morris model, the adsorption process slowed down significantly ((k ipd (3) << k ipd (2)< k ipd (1); see Table 1) due to two features: (i) diffusion through smaller meso/micropores of "hierarchical" zeolites generated during brick activation and (ii) establishment of adsorption equilibrium.The film diffusion rate constants, k ipd (1), were 54.322, 46.860, 15.127, and 7.209 µmol•g −1 •min −1/2 for Cs + , Rb + , Sr 2+ , and Co 2+ , respectively, while the intra-particle diffusion rate constants, k ipd (2), were lower: 20.715, 21.597, 4.192, and 4.600 µmol•g −1 •min −1/2 .Overall, these values indicate that Cs + and Rb + ions diffused similarly into the composite frameworks, while Sr 2+ and Co 2+ ions diffused more slowly through brick pores before being adsorbed by the material.

Radio-Nuclides
In order to evidence and differentiate diffusion regions in adsorption kinetics, we decided to test Weber-Morris equation to qualitatively reveal the effects of the heterogeneous and porous nature of alkali-activated brick on the diffusivity of Cs + , Rb + , Co 2+ , and Sr 2+ ions inside adsorbent frameworks.
It has previously been stated that the Weber-Morris (W-M) model is useful for identifying reaction pathways and predicting rate-controlling steps in porous matrices [85].The model proposed by Weber and Morris [62] was applied in this work to gain deeper knowledge regarding the mechanism of nuclide adsorption onto the porous brick composite in aqueous solutions (see Section 2.9).For each element, the plot of Qt vs. t 1/2 was found to be divided into three stages by two break points (Figure 4).As suggested in previous works regarding adsorption kinetics in relation to the adsorbent pore structures and surface properties of adsorbents [86][87][88][89], it was assumed here that the multi-linearity observed for cation adsorption was due to the effects of the adsorbent porosity on cations.The observed multi-linear plot would then result from the broad pore size distribution of the brick composite, including, essentially, macropores and mesopores and, to a lesser extent, meso/micropores of "hierarchical" zeolites (Figure 3b).The higher value of kipd (1) in the first stage of the Weber-Morris model compared to kipd (2) and kipd(3) revealed the importance of external liquid film diffusion (Table 1).In the following section, we attempt to identify the various determining factors which both controlled the adsorption of cationic elements onto alkali brick pellets and regulated the interfacial diffusion kinetics.For this purpose, we compiled data from the literature on the thermodynamic and physical properties of the ions Cs + , Rb + , Co 2+ , and Sr 2+ (as well as other metallic cations) in water and correlated them with the experimental results.The W-M kinetics curves obtained for the metallic cations were drawn and are shown in Figure S1 (see Supplementary Materials), and the calculated intra-particle diffusion rate constants, k ipd (1) and k ipd (2), are listed in Table S1 (see Supplementary Materials).Note that the values of ionic potential, Jones-Dole viscosity coefficient (B), and hydrated radius used in the present work are also given in Table S1.
Ionic potential is known to be an indicator of the electrostatic properties of an ion when interacting with other ions, either of similar charge (yielding electrostatic repulsion forces) or of opposite charge (yielding electrostatic attraction forces).Ionic potential is defined as the ratio of the electrical charge of an ion (Z) to its radius (R).It was previously considered as an important criterion for the choice of efficient adsorbents for toxic elements [90].In order to calculate the ionic potential, we fixed Z = +1 for monovalent ions and +2 for divalent ions and took the values of hydrated radii (R H ) proposed in Marcus' works [91].According to this author, a cation with an ionic radius, R i , is surrounded with water in a spherical shell, ∆R H2O , in which several H 2 0 molecules are immobilized by the electrostatic field of the cationic entity.The total distance from the centre of the cation to the edge of the shell is then given by R Total = R H = R i + ∆R H2O .Figure 5 displays the evolution of the intra-particle diffusion rate constants, k ipd (1) and k ipd (2), as a function of the ratio Z/R H .
The k ipd (1) rate constant exhibited a relatively good correlation with the ratio Z/R H , at least for divalent cations (R 2 = 0.9365), indicating a noticeable contribution of electrostatic forces to external liquid film diffusion.Conversely, the weak correlation of the k ipd (2) rate constant with the ratio Z/R H (R 2 = 0.2176) seemed to show a less important electrostatic contribution to intra-particle diffusion in comparison to that observed in the first stage of the Weber-Morris model (Figure 5).
Figure 5 displays the change in the hydrated radius (R H ) in relation to the intraparticle diffusion rate constants k ipd (1) and k ipd (2) for the first and second stages of the Weber-Morris model, respectively.As can be seen in this figure, the k ipd (1) and k ipd (2) rate constants exhibited relatively good correlations with R H (R 2 = 0.9649 and 0.9735, respectively).The observed hydrated radius (R H ) dependence on diffusion rate constants suggested that the size of the hydrated ion ought to be taken into consideration as a noticeable factor in the diffusional approach of the studied adsorption process.
In the diffusion adsorption mechanism, the mobility of the cationic species first toward brick surfaces and second through zeolitic pores/channels should depend on the viscosity of the interfacial medium.The viscosity of the hydrated cation solution (η) was referred to that of the pure water (η o ), and could be determined from the Jones-Dole expression [92]: where C (in mol•L −1 ) is the solute concentration and A (in L 1/2 •mol −1/2 ) and B (in L•mol −1 ) are coefficients depending on the solute, the solvent, and the temperature [92].Here, we consider the Jones-Dole viscosity coefficient (B) as a potential indicator of the degree of water structures around the cationic entity, and thereby of the bonding strength of hydration shells covering the element.The B-coefficients of aqueous inorganic ions at 25 • C in L.mol −1 were previously reported by Jenkins and Marcus [92].Thus, by plotting the intra-particle diffusion rate constants k ipd (1) and k ipd (2) against the B "viscosity" parameter (Figure 5), we noticed that the increase in the medium viscosity contributed to a significant inhibition of cationic adsorption kinetics.This suggests that the B-coefficient also had to be taken into account as a determining parameter regulating cationic diffusion into brick surfaces.
In water, the complexation of a cation with H 2 O molecules results in the formation of a hydration shell.This reaction is exothermic, releasing a certain hydration free energy noted here as "∆G H ". To evaluate the contribution of the molar Gibbs energy of hydration of the cation M + or M 2+ (as defined and determined by Marcus [91]) to the diffusion behaviour of cationic species inside composite frameworks, we plotted the intra-particle diffusion rate constants, k ipd (1) and k ipd (2), against "∆G H " (Figure 5).The values of molar Gibbs energy of hydration used for plotting are listed in Table S3.As can be seen in Figure 5, the k ipd (1) and k ipd (2) constants exhibited good correlations with the molar Gibbs energy of hydration (R 2 = 0.9923 and 0.9852, respectively).These findings suggest that, for Cs + and Rb + , which have much softer hydration shells than those of Ni 2+ , Cu 2+ , Co 2+ , Sr 2+ , and Zn 2+ , the easier rearrangement and loss of the water molecules of hydration around these cations contributed to an enhancement of diffusion kinetics.
Globally, the facile kinetic exchange between Na + ions on the brick surface and Cs + or Rb + ions in the solution could be explained in part by the low differences in the ion hydration energies (−385 kJ•mol −1 for Na + , −285 kJ•mol −1 for Rb + , and −245 kJ•mol −1 for Cs + [91]).These differences were substantially smaller in magnitude than those calculated for the ion exchanges between Na + ions and the cationic metals Sr 2+ and Co 2+ (−385 kJ•mol −1 for Na + , −1468 kJ•mol −1 for Sr 2+ , and −2003 kJ•mol −1 for Co 2+ [91]).But the larger barrier to adsorption kinetics for the more strongly hydrated cations (Co 2+ and Sr 2+ ) could not be understood simply by the very high differences in ion hydration energies.Indeed, partial losses of water molecules in cation hydration shells and in oxygen atoms on the brick surface should be responsible for energy barriers to adsorption and, hence, for the magnitude of diffusion kinetic constants observed.On this view, it could be reasonably stated that the H 2 O losses in the coordinating structures of the cations involved in Cs + and Rb + adsorption and Na + desorption were thermodynamically easier and thereby kinetically more favourable than those in the coordinating structures of Sr 2+ and Co 2+ ions.
To summarize, this study clearly shows that adsorption kinetics took place initially with internal film diffusion involving macropores and wider mesopores, followed by intraparticle diffusion involving smaller mesopores and zeolitic micropores.The Weber-Morris model [62] was used here as a first approach for describing the interfacial diffusion aspects of the adsorption kinetics on the brick composite.However, although the application of the Weber-Morris model evidenced the implication of distinctive diffusion regimes in connection with the porous properties of the adsorbent, it did not allow for the measurement of ionic diffusivities of cations, either in the bulk solution (towards the solid surface) or in the solution confined inside the porous brick matrix.
In the following section, the restricted diffusion of cations in the adsorption on the brick composite is analysed and quantified by relating the successive diffusional regimes observed in W-M adsorption kinetics to the specific pore size distribution of the adsorbent.
Appl.Sci.2024, 14, x FOR PEER REVIEW 15 of 34 water in a spherical shell, ΔRH2O, in which several H20 molecules are immobilized by the electrostatic field of the cationic entity.The total distance from the centre of the cation to the edge of the shell is then given by RTotal = RH = Ri + ΔRH2O.The kipd(1) rate constant exhibited a relatively good correlation with the ratio Z/RH, at least for divalent cations (R 2 = 0.9365), indicating a noticeable contribution of electrostatic forces to external liquid film diffusion.Conversely, the weak correlation of the kipd(2) rate Figure 5. Dependence of the intra-particle diffusion rate constants, k ipd (1) and k ipd (2), upon ionic potential (Z/R H ); hydrated radius (R H ); Jones-Dole viscosity coefficient (B); and molar Gibbs energy of the cation's hydration, M + or M 2+ (∆G H ). (Note that only divalent cations were taken into consideration in the curve: k ipd (1) and k ipd (2) versus ionic potential).

Diffusion-Based Models
In this section, different diffusion models were used for determining the key diffusion parameters governing the diffusion kinetics of the adsorption of elements onto the brick composite.

Evaluation of Film Mass Transfer Coefficient
W-M studies permitted us to show that the external film diffusion took place only at the early stage of the adsorption kinetics.
Considering the film-surface diffusion model proposed by Yao and Chen [63], the initial adsorption rate r o could be determined (see Section 2.9).Plotting t/Q t against time, good straight lines were obtained with correlation coefficients ranging from 0.9988 to 0.9996 for the different ion exchange reactions studied (Figure 6).From the slope and intercept values, the "r o " and "a" constants were determined, and results are listed in Table 2.The k f value was then determined and is given in Table 2.
Appl.Sci.2024, 14, x FOR PEER REVIEW 17 of 34 observed in W-M adsorption kinetics to the specific pore size distribution of the adsorbent.

Diffusion-Based Models
In this section, different diffusion models were used for determining the key diffusion parameters governing the diffusion kinetics of the adsorption of elements onto the brick composite.

Evaluation of Film Mass Transfer Coefficient
W-M studies permitted us to show that the external film diffusion took place only at the early stage of the adsorption kinetics.
Considering the film-surface diffusion model proposed by Yao and Chen [63], the initial adsorption rate ro could be determined (see Section 2.9).Plotting t/Qt against time, good straight lines were obtained with correlation coefficients ranging from 0.9988 to 0.9996 for the different ion exchange reactions studied (Figure 6).From the slope and intercept values, the "ro" and "a" constants were determined, and results are listed in Table 2.The kf value was then determined and is given in Table 2.  To summarize, Equation (3) permitted us to estimate the film mass transfer coefficient (kf) through the initial adsorption rate (ro).This latter was determined at the early adsorption stage (i.e., for a short adsorption time of less than 10 min for Cs + , Rb + , and Sr 2+ and 14 min for Co 2+ , as shown in Figure 6).The kf values will be used in Section 3.4 for modelling breakthrough curves obtained under dynamic (fixed-bed column) operational conditions.To summarize, Equation (3) permitted us to estimate the film mass transfer coefficient (k f ) through the initial adsorption rate (r o ).This latter was determined at the early adsorption stage (i.e., for a short adsorption time of less than 10 min for Cs + , Rb + , and Sr 2+ and 14 min for Co 2+ , as shown in Figure 6).The k f values will be used in Section 3.4 for modelling breakthrough curves obtained under dynamic (fixed-bed column) operational conditions.

Evaluation of Pore and Surface Diffusion Coefficients
In the second steps of the Weber-Morris plots, the (single) diffusion of the cations Cs + , Rb + , Sr 2+ , and Co 2+ through the brick slowed down because the adsorption process became more restricted by brick pores with lower diameters.Under this condition, Fick's law might be applied to diffusion kinetics: where C is the concentration of the adsorbate; D eff is the effective diffusion coefficient of the cation; and r represents the radial coordinate.Assuming that the adsorbent particles are spherical and that the pores are identical and homogeneously distributed on the solid surface, Equation ( 13) could be solved to Equation ( 14) [93] and applied to our system.
where X(t) = Q t /Q e is the fractional attainment at time t; Q e and Q t represent the amounts of cation adsorbed on brick composite at equilibrium and time t (mg•g −1 ), respectively; D eff is the effective diffusion coefficient of the cation in the composite exchanger phase (m 2 •s −1 ); R b is the radius of brick grains, which are assumed to be spherical (m); and "z" is an integer.The following equation (valid for X(t) varying from 0 to 1) was used here as an approximation of Equation ( 14): 1.2 (15) This equation can be transformed into: In this equation, it was necessary to introduce a constant (noted "Cst") in order to take into account the reaction time t o when the adsorption kinetics process began to be controlled by surface diffusion according to the Weber-Morris model.This constant could be determined from the following formula: It should be noted that the "Cst" value differed for each cation studied.Hence, Equation ( 16) ought to cover only the kinetic data points which were used for determining the Weber-Morris rate constant k ipd (2).Its use permitted us to reach the effective particle diffusivity as follows.By plotting −ln [1 − X(t) 2 ] against time, we obtained relatively good straight lines with correlation coefficients varying from 0.9741 to 0.9872 for the different adsorption reactions studied.See Figure 7.
The slope values of these linear relations were employed to assess the effective diffusion coefficients, and the D eff values are listed in Table 3.
Likewise, from the coordinates of kinetic data points representative of the third W-M stage, we plotted −ln [1 − X(t) 2 ] against time and estimated the effective particle diffusivity for each element by employing the same calculation procedure as that described for the second W-M stage (see Figure 7 and Table 3).In this case, Equation ( 16) ought to cover only the kinetic data points which were used for determining the Weber-Morris rate constant k ipd (3).The slope values of these linear relations were employed to assess the effective diffusion coefficients, and the Deff values are listed in Table 3.  6); (b) Deff, calculated from models established by Boyd et al. and discussed in ref. [93] and Vermeulen's approximation also documented in ref. [  Likewise, from the coordinates of kinetic data points representative of the third W-M stage, we plotted −ln [1 − X(t) 2 ] against time and estimated the effective particle diffusivity for each element by employing the same calculation procedure as that described for the second W-M stage (see Figure 7 and Table 3).In this case, Equation ( 16) ought to cover only the kinetic data points which were used for determining the Weber-Morris rate constant kipd (3).

93] at reaction time values of the different experimental points of the second stage and third stage of the Weber-Morris (W-M) kinetics model;
To interpret the results of the diffusion (Deff) measurements of the mass-transfer kinetics involved in the second Weber-Morris step, it was assumed in this work that both pore diffusion and surface diffusion might be responsible for the mass transfer resistance.Hence, in order to better model diffusion in the brick composite based on the step diffusion through mesopores and micropores, we decided to apply a heterogeneous model in which both pore diffusion and surface diffusion were considered [94][95][96].By taking into account the pore diffusion and surface diffusion intrinsically involved in the second Weber-Morris step and by using Fick' theories, the total diffusion flux (J) could be expressed as [97][98][99][100]:   To interpret the results of the diffusion (D eff ) measurements of the mass-transfer kinetics involved in the second Weber-Morris step, it was assumed in this work that both pore diffusion and surface diffusion might be responsible for the mass transfer resistance.Hence, in order to better model diffusion in the brick composite based on the step diffusion through mesopores and micropores, we decided to apply a heterogeneous model in which both pore diffusion and surface diffusion were considered [94][95][96].By taking into account the pore diffusion and surface diffusion intrinsically involved in the second Weber-Morris step and by using Fick' theories, the total diffusion flux (J) could be expressed as [97][98][99][100]:

Cations
where D p and D s are the pore diffusion coefficient and solid-phase diffusion coefficient, respectively; C and C s are the concentrations of adsorbate in bulk external liquid and on the adsorbent surface per unit volume, respectively; and r represents the radial coordinate.We considered the transport of cations and mass balance equation within the pores of a spherical adsorbent by using the following expression [100]: By describing diffusion kinetics from Fick's law according to Equation ( 13), one could establish a mathematical expression evidencing both solid-phase diffusion contribution and pore diffusion contribution according to [100]: where ε bc is the internal porosity of the brick composite; ρ bc is the density of the porous support; C is the liquid-phase concentration of the cation per unit volume of the liquid phase; and Q represents the solid-phase concentration per unit mass of brick.As suggested by Zhu and his co-workers [100], under favourable (and ideal) adsorption conditions, it could be assumed that: (i) the concentration of adsorbates on the surface of the adsorbent (C s ) was close to the product "Q•ρ bc "; and (ii) as the examined diffusion region tended to attain thermodynamic equilibrium, the ∂Q/∂C derivative could be assimilated to the derivative at the equilibrium state, (∂Q/∂C) eq .By taking into consideration these hypotheses and according to some authors [98,100], Equation ( 19) could be rearranged, and the effective diffusion coefficient could be deduced and expressed as a function of both the pore diffusion coefficient (D p ) and the surface diffusion coefficient (D s ) as: On the other hand, from batch experiments carried out on the adsorption of ions (Cs + , Rb + , Sr 2+ , and Co 2+ ) on to the brick composite, the Langmuir [101] and Freundlich [102] isotherm models were fitted to the experimental data.Both models described the equilibrium data quite well, with relatively good correlation coefficients (R 2 > 0.98).However, the Langmuir model provided lower standard errors with higher R 2 coefficients, evidencing preferential monolayer adsorption of cations on brick surfaces.The derivation of the Langmuir equation permitted us to deduce it (∂Q/∂C) and replace it with the following expression: From this equation, and by considering the following calculation as a first approximation, (∂Q t /∂C t ) was determined at the reaction time values of the different experimental points representative of the second stage of the Weber-Morris kinetics model.These calculations permitted us to show that, for each element, the (∂Q t /∂C t ) values remained close enough to (∂Q/∂C) eq , which enabled the application of Equation (20) to the delimited intra-particle diffusion zone (second step) revealed by Weber-Morris kinetics curve.
On the other hand, as shown in Section 2.9, the ionic diffusivity of each cation could be attained by using the values of the infinite dilution diffusion coefficients D 1 • for cations (Cs + , Rb + , Co 2+ , or Sr 2+ ) and D 2 • for anions (NO 3 − ) (D 1 • and D 2 • values are given in Table S2) in Equation (8).The values of D are listed in Table S2.From Equations ( 6) and (7), it was then possible to calculate the D p coefficient for each cation studied.The data are listed in Table 3.
Afterwards, the determination of D p and D eff allowed us to estimate the surface diffusion coefficients (D s ) for cation adsorption on brick composites from aqueous solutions at 298 K using Equations (20) and (21) (Table 3).As seen in this table, the values of the surface diffusion coefficients ranged globally from 3.89 × 10 −12 m 2 •s −1 to 4.22 × 10 −11 m 2 •s −1 , indicating that they were in the same order of magnitude as those of the effective diffusion coefficients (4.79 × 10 −12 < D eff < 1.18 × 10 −11 m 2 •s −1 ; see Table 3); however, they were lower than the pore diffusion coefficients (8.77 × 10 −11 < D p < 1.40 × 10 −10 m 2 •s −1 ; see Table 3).These findings seemed to reveal a predominance of surface diffusion on the (intra-particle diffusion) adsorption kinetics.
The mathematical method described above for the determination of D p , D efff , and D s relative to the second Weber-Morris stage, was also used for the experimental points of the third W-M stage.By plotting −ln [1 − X(t) 2 ] against time, we obtained relatively good straight lines with correlation coefficients varying from 0.9137 to 0.9909 for the different adsorption reactions studied (see Figure 7).
In this case, from the kinetic data points used for determining the Weber-Morris rate constant, k ipd (3), it was noticed that, except for caesium, the D eff value was found to be closer to the D s one for each element studied (see Table 3).Moreover, for long adsorption times, one could approximate Equation ( 14) by retaining only the first term in the series, which led to the following equation [93]: As explained above for Equation (16), it was also necessary to take into account a constant Cst' in Equation (22).This Cst' constant could be assessed from the formula , where the reaction time t o ' must be taken at the beginning of the third W-M stage.As for plot calculation, Equation ( 22) ought to cover only the kinetic data points which were used for determining the Weber-Morris rate constant k ipd (3).
Hence, considering the experimental points of the third W-M stage, the term -Ln(1 − X(t).(π 2 /6) was plotted against the time t (Figure 8), and the value of the slope permitted us to attain the D eff coefficient.All the D eff values evaluated from Equation ( 22) were found to be comparable enough with those calculated from Equation ( 16) (see Table 3).Particularly, for each element, the D eff value was closer to the D s one, suggesting that the pore diffusion contribution was negligible at the end of the kinetics process.On the other hand, in order to evidence the importance of ion dehydration for the diffusivity of cationic species inside composite frameworks, the surface diffusion coefficient (Ds) was plotted against the molar Gibbs energy of hydration of the cation M + or M 2+ (ΔGH).See Figure 9 (data used for plotting Figure 9 are given in Table S3).On the other hand, in order to evidence the importance of ion dehydration for the diffusivity of cationic species inside composite frameworks, the surface diffusion coefficient (D s ) was plotted against the molar Gibbs energy of hydration of the cation M + or M 2+ (∆G H ). See Figure 9 (data used for plotting Figure 9 are given in Table S3).
As seen in this figure, Sr 2+ and Co 2+ ions (and, more generally, heavy metal cations) diffused more slowly through the brick pores, not only because of hindered or restricted effects caused by their larger radius sizes (as shown in Figure 5), but also due to their elevated free energy of dehydration (as also shown in Figure 5).
In the next paragraph, we use batch adsorption data to achieve the modelling of breakthrough profiles for lab-scale fixed-bed adsorption columns by considering a simplified mathematical approach involving interfacial diffusion phenomena.On the other hand, in order to evidence the importance of ion dehydration for the diffusivity of cationic species inside composite frameworks, the surface diffusion coefficient (Ds) was plotted against the molar Gibbs energy of hydration of the cation M + or M 2+ (ΔGH).See Figure 9 (data used for plotting Figure 9 are given in Table S3).As seen in this figure, Sr 2+ and Co 2+ ions (and, more generally, heavy metal cations) diffused more slowly through the brick pores, not only because of hindered or restricted  In order to simulate dynamic fixed-bed adsorption breakthrough curves, we made the following basic assumptions: (i) the operating adsorption process was assumed to be isothermal; (ii) the liquid phase flowed according to an axial dispersion; (iii) the compressibility of the liquid phase was negligible; (iv) the liquid phase inside the pores was considered to be stationary and its flow was assumed not to be altered by movement of the mobile phase; (v) the brick bed was homogeneous and the brick particles were spherical in shape and uniform in size; and (vi) the liquid phase velocity, fluid flow rate, and column void fraction remained constant during the entire column experiment.In the modelling procedure, we further hypothesized that the adsorbate was first adsorbed on the external surface of the brick particle and then diffused inside the particle.The obtained thermodynamic equilibrium and kinetics data were used here from independent batch experiments with the aim of providing predictions of fixed-bed breakthrough behaviour (Table 4).As the shape of the breakthrough curve depended on the adsorption equilibrium, intra-particle mass transfer, and hydrodynamic conditions, it was preliminarily necessary to determine the ideal operational conditions (such as adsorbent particle size, flow rate and input concentration) for optimizing the system and producing steep (or ideal) breakthrough curves for optimal column performance [84].The optimization of breakthrough curves of cations was realized under well-defined experimental conditions, which are documented in Table 4.A mathematical model of fixed-bed adsorption (HSDM) was employed herein for the simulation of experimental breakthrough curves.The homogeneous surface diffusion model (HSDM; not detailed here) has been reported extensively in the past for describing the sorption dynamics of porous materials [103][104][105][106][107]. Simulating calculations were undertaken in the present work by using the fixed-bed adsorption simulation tool (FAST 2.1) [58].Note that, from the FAST 2.1 software, mathematical formulations applying HSDM were based upon driving forces that did not separate pore diffusion from surface diffusion by fitting experimental breakthrough curve data.Figure 10 depicts experimental breakthrough curves corresponding to the dynamic adsorption of cations onto the brick composite.In the same figure, simulated breakthrough curves obtained from the HSDM model using FAST 2.1 software are represented as well.
Globally, the comparison of numerical results from the HSDM model with experimental data revealed relatively good model predictions.The numerical value of the diffusion coefficient obtained for each column study of element is provided in Table 4.
To summarize, the present work permitted us to simulate, in a single system, the adsorption dynamics of the ions Cs + , Rb + , Co 2+ , and Sr 2+ on alkali-activated brick.Breakthrough curves of cationic adsorption experiments were successfully described from the HSDM model.The simulated values of diffusion parameters (k f and D s ) obtained from linear driving force approximation together with the Langmuir isotherm constant (K L ) and maximum adsorption capacity (Q AC ) agreed well with those obtained experimentally under static (batch) operational conditions.The goodness-of-fit parameter then validated the successful prediction of experimental breakthrough curves and evidenced the predominance of surface diffusion for the studied cations.To support this, in the next paragraph, we calculate the different column dimensionless parameters which defined the shape of the breakthrough curve via adsorption dynamics, and use them to verify the importance of surface diffusion phenomena along the packed-brick pellets column.
fixed-bed adsorption simulation tool (FAST 2.1) [58].Note that, from the FAST 2.1 soft-ware, mathematical formulations applying HSDM were based upon driving forces that did not separate pore diffusion from surface diffusion by fitting experimental breakthrough curve data.Figure 10 depicts experimental breakthrough curves corresponding to the dynamic adsorption of cations onto the brick composite.In the same figure, simulated breakthrough curves obtained from the HSDM model using FAST 2.1 software are represented as well.Globally, the comparison of numerical results from the HSDM model with experimental data revealed relatively good model predictions.The numerical value of the diffusion coefficient obtained for each column study of element is provided in Table 4.

Simulation of breakthrough curve (Cs
To summarize, the present work permitted us to simulate, in a single system, the adsorption dynamics of the ions Cs + , Rb + , Co 2+ , and Sr 2+ on alkali-activated brick.Breakthrough curves of cationic adsorption experiments were successfully described from the HSDM model.The simulated values of diffusion parameters (kf and Ds) obtained from linear driving force approximation together with the Langmuir isotherm constant (KL) and maximum adsorption capacity (QAC) agreed well with those obtained experimentally under static (batch) operational conditions.The goodness-of-fit parameter then validated

Dimensionless Analysis of Diffusional Regimes
In the precedent paragraph, we employed the homogeneous surface diffusion model (HSDM) for simulating Na + /M + (or 2Na + /M ++ ) exchanges in the fixed-bed column and predicting the adsorption dynamics of nuclides.In the following section, we interpret dimensionless numbers relative to column operations (Biot, B i ; Stanton, S t ; solute distribution constant, D g ) with the aim of gaining information regarding the shapes of breakthrough curves and breakthrough times, as well as the types of diffusional regimes involved [105,[108][109][110].In this work, the dimensionless numbers B i , S t , and D g (as defined in Section 2.9) were determined from experimental adsorption results for different elements in alkali-activated brick, and their values are listed in Table 4. Globally, characteristic dimensionless numbers ranged from 1091 to 5237 for B i , from 162 to 256 for D g , and from 296 to 778 for S t .Ratio Ed = S t /B i was also considered in the dimensionless analysis, knowing that its implication on the adsorption process was directly related to the individual values of the B i and S t numbers.The very large values of B i revealed that film mass transfer occurred very faster than internal diffusion, and hence, internal surface diffusion dominated the adsorption process.The high Biot number values also permitted us to better understand the influence of D s on the dynamics of breakthrough curves.Indeed, since the obtained B i values were high and D s and B i were inversely proportional, these should lead to low D s values, as evidenced experimentally (D s values ranging from 2.00 × 10 −12 to 6.66 × 10 −12 m 2 •s −1 ).Dimensionless analysis further confirmed the occurrence of slow diffusional regimes, which were characterized by smoothed wavefronts.Moreover, accord-ing to Pérez-Foguet and his co-workers [108], longer computational times are expected to be observed for large values of D g in an elevated adsorbent medium.In the studied adsorption systems, as high time-scales were given by D g , weak internal diffusions took place and significantly controlled the adsorption process.This finding was supported by the weak values of Ed, which also confirmed the involvement of slow adsorption diffusion.

Implication of the Connectivity of Pores on Ion Diffusivity
To evaluate the immobilization capacity of a binder toward nuclides, the physical barrier effect of the pore structure on ionic diffusivity ought to be taken into consideration [115].For cementitious materials like ordinary Portland cement, fly ash-based geopolymers, and slag-blended geopolymers, it has previously been shown that the capillary porosity and connectivity of the capillary pores in porous media play a relevant role in ionic diffusion within these porous materials [115,120,121].The capillary pores with pore diameters rang-ing from 10-50 to 10,000 nm [115,121] were indeed assumed to dominate the ionic transport processes in cementitious materials.
It has previously been suggested that the broad range of D eff observed in previous studies was mostly attributed to the variability of the porosity, capillary pore volume, and critical pore diameter, which control the ion diffusivity in cementitious materials [115].In order to describe the connectivity of such porous solids, it was proposed to use the critical pore diameter as a key parameter representative of the highest proportion of interconnected pores, which would be mainly responsible for ion transport characteristics [115,121,122].The critical pore diameter (or continuous pore diameter) was defined as the maximum of the Log differential pore volume, dV/dlog(diameter), versus the pore diameter curve.With this in mind, some investigations permitted us to correlate the critical pore diameter of the binder with the diffusivity behaviour of water-soluble ions such as caesium and strontium [115].
In order to appraise the physical barrier effect of treated brick on the diffusivity of caesium and strontium, we first attempted to attain the critical pore diameter of the brick composite.For this, the calculated dV/d(log d) term was plotted against the pore diameter, and the maximal dV/d(log d) value was determined, as shown in Figure 3c.We found ~50 nm (±10 nm).This value was comparable with that obtained for fly ash-based geopolymers (50-151 nm) and lower than those for fly ash/fly ash-based geopolymers (227-249 nm) and ordinary Portland cement (284 nm), but was much lower than that for fly ash/blast furnace slag-blended geopolymers (1567 nm) [115].Taking into account the empirical linear equation expressing the critical pore diameter of the binder versus the diffusivity of water-soluble ions, which was established by Jang and his co-workers [115], we then evaluated the diffusion coefficient, D soluble , of water-soluble Cs + and Sr 2+ ions.We found D soluble = ~4 × 10 −12 and ~2 × 10 −12 m 2 •s −1 , respectively.These values were somewhat lower than those calculated herein using different diffusion-based models (see Table 3).It could, however, be noted that, based either on the results from surface diffusion models (see Section 3.3) or those derived from the dependence of the connectivity of pores in cement-based materials upon the ion diffusion coefficient (as shown in this paragraph), the diffusion of strontium ions was found, in all cases, to be more difficult than that of caesium ions, which is in agreement, in most cases, with the findings of previous studies (see Table 5).Despite the slight discrepancy between the D soluble and D eff values determined for Cs + and Sr 2+ , this study nevertheless permitted us to demonstrate the implication of hydrated entities of nuclides in the diffusion stages.
From these investigations, one could, therefore, ascertain that each of these elements still remained barely soluble or chemically unfixed to the brick matrix, meaning that the nuclide could be assumed to be adsorbed via hydration changes in the spherical shell of its cationic entity (as pointed out in Section 3.3) and transformed into an unfixed hydrate/complex.To validate this hypothesis, we recently performed a 133 Cs analysis of Cs-doped composite samples and clearly evidenced the existence of (outer-sphere) hydrated-caesium complexes inside the zeolitic frameworks of the modified brick.

Conclusions
In the present work, the diffusional adsorption kinetics of the cations Cs + , Rb + , Co 2+ , and Sr 2+ were investigated in single system.By applying the Weber-Morris model, mathematical treatment of the kinetics (batch) data evidenced a multi-step process with distinctive diffusion regimes that were linked to the pore size distribution of the material.W-M kinetic rate constants were found to be dependent on the following intrinsic characteristics of the hydrated cation: ionic potential; hydrated radius; B-viscosity parameter; and molar Gibbs energy of hydration of cations.By applying adequate diffusion models, a detailed mathematical analysis of the observed W-M kinetics stages allowed us to obtain the key diffusion parameters governing the successive diffusion regimes observed.The surface diffusion coefficient for each cation (D s ) was estimated from the mathematical expression relating it to both the pore diffusion coefficient (D p ) and effective diffusion coefficients (D eff ).
These calculations revealed the strong contribution of surface diffusion to the adsorption kinetics during the course of the second and third kinetics stages of the W-M model.It was noticed that the magnitude of the surface diffusion coefficient was strongly dependent upon the capacity of hydrated cationic species to lose water molecules when penetrating brick pores.
The HSDM model was used for describing and simulating (single) the adsorption of cations in a fixed-bed column under optimized operational conditions.A breakthrough curve simulation clearly indicated the predominance of surface diffusion in the intraparticle diffusion mechanism under fixed operational conditions, in agreement with a mathematical analysis of the (batch) adsorption kinetics data.The simulated (apparent) surface diffusion coefficients were found to be close to those calculated from diffusion-based models, which permitted us to validate the modelling procedure which was employed.
The capillary porosity and connectivity of the capillary pores in the modified brick were found to be relevant characteristics in the ionic diffusion of cations (Cs + and Sr 2+ ) within the porous structure, evidencing the immobilization of hydrated nuclides in the form of outer-sphere complexes.
Since artisanal gold mining activities are growing significantly in developing countries, mercury is being released at elevated levels to the environment.This led our research group to study mercury adsorption on alkali-activated brick.This work will permit future projects to implement dynamic water treatment systems in poor regions, which would be capable of eliminating this toxic metal from immediate contaminated waters and of obtaining potable water for the local population efficiently and at a low cost.Furthermore, experiments regarding the adsorption capability of alkali-activated brick toward soluble uranium species present in waters contaminated by uranium extraction mines in the Central African Republic are also underway in the laboratory.

Future Prospects
The brick composite represents a novel alkali-activated binder for immobilizing radio elements.Its low Si/Al ratio contributes to an abundance of negative surface charges, facilitating ionic exchanges between Na + ions at the solid surface and radioactive ions (Cs + , Rb + , Co 2+ , and Sr 2+ ) in the aqueous phase.The knowledge obtained herein on the diffusion behaviour of these ions inside brick frameworks will help us to better understand nuclear waste immobilization and nuclear radiation shielding after calcination or vitrification.Studies are underway to assess the state of the immobilized radio elements and to follow their degree of encapsulation inside composite frameworks under the effects of high temperatures.For this purpose, by using time-of-flight secondary ion mass spectrometry (ToF-SIMS) combined with sputter depth profiling, we propose the analysis of encapsulated atoms by establishing depth distribution of radioactive atoms inside brick specimens via ionic Bi 3+ bombardment at different durations.Layer-by-layer sputtering will permit the reconstruction of radioactive element distribution by sample depth, as well as the interpretation of leaching behaviour in aggressive media.Using a similar approach, we also propose the micro-analysis of brick composite surfaces with depth using X-ray photoelectron spectroscopy.At the same time, we propose an examination of the chemical and structural transformations of brick composite matrices under the effects of elevated temperatures using X-ray diffraction, solid nuclear magnetic resonance and vibrational spectroscopies, and scanning electron microscopy.
Studies regarding zeolite-geopolymeric transformations of brick under variable alkaline conditions are also underway in the laboratory.This would permit us to determine the amount of extra-framework cations (Na + ) or "active" negatively charged sites, as well as to examine the evolution of porosity with changes in the alkalinity and synthesis time.The objective of this work would be mainly to produce better green, economical, and hierarchically porous materials from brick, with the aim of facilitating the accessibility of more active surface sites for radionuclides through more interconnected meso-and microporous channels, and thereby to increase adsorption capacity.
and εB are the bed density and bed porosity, respectively; Co is the influent concentration; and Qo is the equilibrium solid-phase concentration corresponding to the influent concentration, Co (Qo is determined using the Freundlich equation: Qo = KF•Co n , where KF and n represent Freundlich isotherm coefficients).

Figure 2 .
Figure 2. Morphological aspects of treated brick pellets observed using ESEM at different magnifications (a-c).ESEM micrograph of raw brick (d).Spherical and ellipsoidal surface zones targeted on alkali-activated brick grains for ESEM/EDS analysis (e).ESEM micrograph of alkali-activated brick surface on which the lengths of cube edges and sphere diameters of zeolitic specimens were

Figure 2 .
Figure 2. Morphological aspects of treated brick pellets observed using ESEM at different magnifications (a-c).ESEM micrograph of raw brick (d).Spherical and ellipsoidal surface zones targeted on alkaliactivated brick grains for ESEM/EDS analysis (e).ESEM micrograph of alkali-activated brick surface on which the lengths of cube edges and sphere diameters of zeolitic specimens were measured (f).Note

Figure 3 .
Figure 3. N2 adsorption (in blue) and desorption (in red) isotherm (a) and corresponding pore size distribution curves for brick composite (in blue; b1) and raw brick (in red; b2).Log differential pore volume dV/d(log d) versus pore diameter (c).

Figure 3 .
Figure 3. N 2 adsorption (in blue) and desorption (in red) isotherm (a) and corresponding pore size distribution curves for brick composite (in blue; b1) and raw brick (in red; b2).Log differential pore volume dV/d(log d) versus pore diameter (c).

1 )Figure 4 .
Figure 4. Plotting of the adsorption capacity (Qt, in µmol•g −1 ) against t 1/2 .Application of Weber−Morris kinetics model to adsorption kinetics of cations (Cs + , Rb + , Co 2+ , and Sr 2+ ) on brick com-Figure4.Plotting of the adsorption capacity (Q t , in µmol•g −1 ) against t 1/2 .Application of Weber-Morris kinetics model to adsorption kinetics of cations (Cs + , Rb + , Co 2+ , and Sr 2+ ) on brick composite in aqueous medium.Note that the first, second and third steps observed from W-M model are represented in black squares, red squares and blue triangles, respectively.

Table 3 .
Surface diffusion coefficients values (Deff, Ds, and Dp; expressed in m 2 •s −1 ) obtained by applying the homogeneous surface diffusion model (HSDM) for the adsorption of cations (Cs + , Rb + , Co 2+ , and Sr 2+ ) on brick composite. (a) Dp, calculated from Equation ( (c) Def,f calculated from models established by Boyd et al. or Crank (as described in ref. [93] for long adsorption times; (d) Ds, evaluated from Equation (20).

Figure 7 .
Figure 7. Plotting of −ln [1 − X(t) 2 ] against time by covering kinetic data points which were used for separately determining Weber-Morris rate constants k ipd (2) and k ipd (3).

Table 3 .
Surface diffusion coefficients values (D eff , D s , and D p ; expressed in m 2 •s −1 ) obtained by applying the homogeneous surface diffusion model (HSDM) for the adsorption of cations (Cs + , Rb + , Co 2+ , and Sr 2+ ) on brick composite. (a) D p , calculated from Equation (6); (b) D eff , calculated from models established by Boyd et al. and discussed in ref. [93] and Vermeulen's approximation also documented in ref. [93] at reaction time values of the different experimental points of the second stage and third stage of the Weber-Morris (W-M) kinetics model; (c) D ef , f calculated from models established by Boyd et al. or Crank (as described in ref. [93] for long adsorption times; (d) D s , evaluated from Equation (20).

Figure 9 .
Figure 9. Dependence of the surface diffusion coefficient (Ds) relative to the adsorption of nuclides (Cs + , Rb + , Co 2+ , and Sr 2+ ), monovalent cations (Li + and K + ), and some metallic cations onto brick composites, with molar Gibbs energy of hydration of the cation M + or M 2+ (ΔGH).

Figure 9 .
Figure 9. Dependence of the surface diffusion coefficient (D s ) relative to the adsorption of nuclides (Cs + , Rb + , Co 2+ , and Sr 2+ ), monovalent cations (Li + and K + ), and some metallic cations onto brick composites, with molar Gibbs energy of hydration of the cation M + or M 2+ (∆G H ).

Figure 10 .
Figure 10.Experimental and simulated breakthrough curves relative to the dynamic adsorption of cations (Cs + , Rb + , Co 2+ and Sr 2+ ) on to brick composite.

Figure 10 .
Figure 10.Experimental and simulated breakthrough curves relative to the dynamic adsorption of cations (Cs + , Rb + , Co 2+ and Sr 2+ ) on to brick composite.

Table 2 .
Film mass transfer coefficient evaluated using the film-surface diffusion model proposed by Yao and Chen [63].

Table 2 .
Film mass transfer coefficient evaluated using the film-surface diffusion model proposed by Yao and Chen [63].

Table 4 .
Operating conditions, adsorbent characteristics, and equilibrium/kinetics parameters used for HSDM modelling of the breakthrough curves relative to the (single) adsorption of Cs + , Rb + , Co 2+ , and Sr 2+ on brick composite.

Table 5 .
Comparison of the D eff values (in m 2 •s −1 ) of Cs + , Rb + , Co 2+ , and Sr 2+ ions individually adsorbed on the brick-based composite with those previously reported, concerning their diffusion behaviour through different cements, clays, and clay rocks/soils. -