Abstract

The use of montmorillonite in the context of engineered barriers makes it possible to minimize the spread of heavy metals from industrial and even radioactive waste. An evaluation of the performance of the mechanisms controlling the clay-environment interaction and predicting the dynamics/configuration of the interlayer space (IS) is required. This work focuses on a quantitative identification of the structural changes and porosity alteration in the case of heavy metal-exchanged montmorillonite samples (Co2+ and Cd2+ cations) undergoing mechanical stresses (uniaxial oedometric test (loading/unloading)). Relationships between mechanical stress strength, intrinsic structural response, ion exchanger efficiency, and adsorption performance fate are investigated. This goal is achieved through the correlation of in situ quantitative X-ray diffraction (XRD) analysis (under an extremely controlled atmosphere reached by varying relative humidity rate %rh) and porosity investigation (assured by combining outcomes from BET (Brunauer–Emmett–Teller) and BJH- (Barrett, Joyner, and Halenda-) PSD (pore size distribution) analysis). Obtained results show an upsurge in the structural heterogeneities accompanying the theoretical increase in the mixed layer structure (MLS) number and developing an unconventional hydration behaviour after stress relaxation regardless of exchangeable cation nature. Experimental XRD patterns are reproduced using MLS, which suggests the coexistence of more than one “crystallite” specie and more than one exchangeable cation indicating a complex cation exchange capacity (CEC) saturation. For extremely low %rh value, a new homogeneous dehydrated state trend is observed in the case of the Co2+ cation. Porosity analysis shows mesopore volume growth for the stressed sample and confirms crystallite exfoliation layer trends, results of the layer cohesion damage, and subsequent constraint strength fluctuations.

1. Introduction

Waste disposal has become one of the most serious modern environmental problems associated with the development of countries. The architecture, design, and operation of the waste disposal sites ensure adequate waste management with respect for human health and the ecosystem. The geological disposal, adopted as a safe solution for the long-term management of radioactive waste, is delivered by a system that (1) isolates waste from the biosphere for extremely prolonged periods of time, (2) ensures that residual radioactive substances reaching the biosphere will be at comparatively insignificant concentrations, and (3) provides reasonable assurance that any risk of inadvertent human intrusion would be extremely low [15]. The most common problems of the landfill site are environmental degradation and groundwater contamination [6].

The use of compacted clays, as part of the multibarrier concept, ensures the previously cited requirements [7, 8]. Indeed, the most used materials in landfill disposal are compacted clay liner (CCL) and geosynthetic clay liner (GCL). The main reason for using such materials is their intrinsic properties, such as the large specific surface areas (SSA), the low permeability, the cation exchange capacity (CEC), the high adsorption ability, and the low hydraulic conductivity [611]. Compacted clay liners are less expensive, and it has good attenuation capacity. Clay membranes can be used as toxic gases (originating from the landfill site) eliminator and as a trap for radionuclide migration/diffusion (resulting from a pipeline leak or damage due to various types of stress affecting installation) [1215]. The clay mineral membranes are exposed to various chemical, biological, and physical stresses, and they are influenced by the subsequent leachate [1621]. To evaluate the consistency of these membranes, it is important to examine their chemical compatibility with the different fluids or the leachate to which they are subjected [22]. Based on the fine particle size, the homogeneous micropores and the high surface charges, the low hydraulic conductivity, the high adsorption capacity, and the low cost compared with geosynthetic clay (GC), natural bentonites are most prized in this context [2329]. Indeed, [30] studied the shear strength of compacted clays as affected by the mineral content and wet-dry cycles and demonstrated that the cohesion and frictional strength properties were determined and linked to the proportion of clay minerals and the number of wetting and drying cycles. On the other hand, [31] demonstrated, in the case of compacted clays mixed with a wide range of bentonite for engineered barriers, that compressibility was affected by the bentonite content on geotechnical characteristics.

The use of bentonite, especially as a buffer in the excavated escape galleries between the waste containers and the tunnel walls, achieves important strategic outcomes. Bentonite membrane is usually affected by the surrounding temperature/humidity gradient fluctuation [32, 33]. However, compacted clays present problems with cracking and/or desiccation, especially those containing an appreciable amount of bentonite [3436]. Numerous studies have shown that compacted clays undergo large physicochemical properties modifications when exposed to a cycle of swelling and/or desiccation-wetting, which constitutes a surrounding variable environmental constraint [2839]. Montmorillonite fraction, which is a dioctahedral species of the smectite mineral group, can constitute up to 95% of bentonite. These functional materials provide restriction of groundwater access to radioactive waste, promote conditions for mass transfer between waste and groundwater by diffusion, overwhelm the migration of radionuclides in colloidal form into groundwater, ensure effective sorption after possible depressurization of radioactive waste containers, and so on [4045]. In addition to their presence in the formations housing the storage facilities and their use as backfill material or as sealing material to isolate parts of disposal facility, montmorillonite is also present as a buffer material in the empty space between the package and the host rock [4648]. The montmorillonite layer structure is constituted by the stacking of two types of sheet: an octahedral sheet O (MO4 (OH)2 with M is a metal cation (Al, Mg)) sandwiched between two tetrahedral sheet T (SiO4) often labelled T-O-T with an average layer thickness around 10 Ǻ. The layer charge is related to existing isomorphic substitutions in the tetrahedral sheet (Al3+, Si4+) and/or octahedral ones (Mg2+/Li+, Al3+/Mg2+). This charge deficit is compensated by exchangeable cations (Exc.Cat), coming from the soil solution and incorporated in IS. The presence of solvent (water molecule) facilitates exchangeable cations insertion and enables, respectively, layer thickness expansion, mineral dispersion (colloidal properties) in the soil solution, and micro/macroscopic clay swelling. The progressive layer thickness expansion of the basal spacing value d (001) is done by discrete hydration state going from the dehydrated state (0 W, d001 ≈ 10 Å) to the strongly hydrated ones (4 W, d001 ≈ 21 Å) passing, respectively, by 1 W (d001 ≈ 12.4 Å), 2 W (d001 ≈ 15.4 Å), and 3 W (d001 ≈ 18.2 Å) hydration state [4954].

Montmorillonite fraction interacts physically and chemically with soil and the disposal packages over time, which results in structural changes and alteration of the intrinsic properties [4860]. Depending on these operating constraints (aggressive variable environment), a thorough predictive study of the behaviour and montmorillonite structural response is required. This response can vary from significant structure transformation and properties to yield affected by the nature of the stress.

The applied constraints (aggressive environment) can be considered individually. In this case, they are divided into three groups: (i) atmospheric (temperature, relative humidity, and pressure) and hydrogeological (soil desiccation/wetting cycles) constraints; (ii) chemical constraints (soil solution composition, water chemistry, and soil solution pH); (iii) mechanical stresses (phenomena of sliding/shearing of soils, tectonics and fracturing, and natural mechanical field) [5053]. In reality, these constraints are often coupled with combinations up to the third coupling order inducing an enormous disruptive stress potential, which affect necessarily the clay fraction efficiency. The performances of the starting materials, according to the initial objectives, are no longer guaranteed; the main problem in this case becomes the effectiveness of these stressed/disturbed materials.

In this work, structural changes, hydration behaviour evolution, CEC alteration, and porosity fate of dioctahedral smectite (Wyoming montmorillonite) subsequent to the applied mechanical stress (uniaxial load/unload) are explored. An experimental investigation protocol, developed at the laboratory scale, based on three main levels is realized. Indeed, mechanical stress based on the loading/unloading test is done to create material stress. An in situ XRD study, under an extremely controlled atmosphere (variable %rh), is directed. The strength of this structural analysis lies in the use of the theoretical 00l diffraction profiles modelling approach to define the relationship between mechanical constraint and microstructure changes. Finally, the adsorption measurement and the PSD analysis are done by targeting the stress damage to the metal-exchanged montmorillonite porosity features.

2. Materials and Methods

2.1. Baseline Sample

Starting sample is a reference montmorillonite specimen (Swy-2) supplied by the Source Clays Repository of the Clay Minerals Society. The intrinsic characteristic properties of this reference sample can be consulted in the literature data [61, 62]. The structural formula per half unit cell is given by [62]

2.2. Experimental
2.2.1. Na-Rich Montmorillonite

Raw material requires NaCl pretreatment before application to ensure maximum layer dispersion. This specific treatment is carried out by respecting a classical protocol detailed by [60, 63] to prepare the Na-rich montmorillonite suspension. Indeed, the <2 μm fraction of Wyoming montmorillonite was extracted by allowing larger particles to settle in water and then decanting. The Na-rich montmorillonite was prepared by equilibrating the colloidal fraction with NaCl solution. This step consists of dispersing 40 g of solid in 500 mL of NaCl solution (1 M) and stirring mechanically for 24 h. A SIGMA laboratory centrifuge at 8000 rpm speed is used for the solid fraction separation. To ensure saturation of all exchangeable sites by Na+ cations, this equilibration was repeated five times. The clay was then washed with distilled water and separated in a high-speed centrifuge until the solution was free of CI. Chloride was presumed to be absent when the supernatant liquid gave a negative test with AgN03.

2.2.2. Metal-Exchanged Montmorillonite

The same ionic exchange procedure (used in the case of Na+ cation) was followed to prepare, respectively, Co-montmorillonite (SWy-Co) and Cd-montmorillonite (SWy-Cd) suspension. This goal is realized using metal chloride solution (1 M) and Na-rich montmorillonite respecting the following chemical reaction:where M2+ denotes Co2+ and Cd2+ cations and SWy-M2+ denotes SWy-Co and SWy-Cd.

2.2.3. Experiments

Before beginning the study of the compaction/reswelling cycle effect on the stressed samples, a study of the unstressed samples under the same controlled atmosphere is essential. The detailed experimental protocol is summarized in Figure 1.

2.3. Mechanical Tests: Materials Fatigue

The oedometer testing device used is a front-loading oedometer (WF 24251) (Figure 2). The cell diameter is 50 mm. The comparator race extends 12 mm with an accuracy of 2 microns. A predetermined piston pressure can be applied to the metal-exchanged montmorillonite sample, causing fluid expulsion through the filter until the system reaches equilibrium. The Linear Voltage Displacement Transducer (LVDT) measures the displacement of the piston and records the change in the system volume during the test. The measured equilibrium displacement permits the determination of the equilibrium void ratio e, which is defined as the ratio of the fluid volume to that of solid particles. The sample was maintained in contact with an external reservoir of deionized water until the experiment ended. All pressures quoted are gauges; that is, they are pressures relative to atmosphere and hence correspond to the differential pressures causing filtration. Each test consisted of a step-by-step loading and unloading (compaction and swelling) of the clay suspensions. The pressure sequences were 5, 10, 15, 20, 25, and 30 bar for the clay suspension, the same pressures being used during compaction and swelling. At the end of each test, the pressure was totally released, and the sample allowed swelling until equilibrium was reached [64, 65]. Then, the final void ratio e was measured by drying the cake at 140∼ overnight. From the known weight of clay fraction M (the amount of suspension introduced in the cell was carefully measured) and the measured piston displacements, the equilibrium void ratio at each pressure can be calculated from the equation

In this equation, the values are taken relative to the final position of the piston and may be positive or negative. The dry clay density was taken as 2.75 g/cm3 as previously calculated. In both pieces of apparatus, the sample cell is of cross-sectional area A = 5.07 cm2.

2.4. Structural Analysis

Experimental XRD patterns are obtained from a Brucker D8 Advance X-ray diffractometer with the following setting and scanning parameters: (i) 40 kV and 20 mA; (ii) Cu_Kα monochromatic radiation λ = 0.15406 nm; (iii) 0.04°2θ as step size; (iv) angular range 3, 5–60°2θ; (v) 6 s as counting time per step. Diffractometer installation is equipped with an advanced Ansyco humidity (rh) generator connected to a CHC + Cryo & Humidity Anton Paar Chamber. To decode the structural impact of the applied mechanical constraint on the possible deformations affecting the IS configuration, a study of unstressed and stressed samples under extreme values of %rh is carried out. This appeal targets the understanding of the interstitial water retention and release mechanisms. The environmental rh variation extends from 3% to 97%. The desired rh rates are maintained by equilibrating samples for 0.5 hours with their environment.

2.5. Qualitative XRD Analysis

The qualitative XRD analysis is achieved by QualX 2.0 program [66]. The commercial PDF-2 database associated with the new freely available database POW_COD is used. This analysis targets to identify all present sample phases and provide information, respectively, about the layer thickness (hydrous indication), the 00l reflections positions, and the diffracted profile geometry (based on the symmetry and/or asymmetry peak observations). Additionally, the calculation of “Full Width at Half Maximum” (FWHM) and the rationality deviation parameters (ξ) related to the 00l reflection (calculated as the standard deviation of the 𝑙 × 𝑑00l values for all 𝑋i measurable reflections over the explored 2θ° angular range) can provide information on the degree of hydration heterogeneity [6670]. However, the qualitative interpretation of the XRD profile cannot provide a detailed insight into the coexistence of the different phases with their associated relative proportions as well as their different hydration states, which may coexist in the structure. On the other hand, it does not make it possible to follow the evolution and/or the fluctuations of certain structural parameters along the axis. For these reasons, it is necessary to introduce the quantitative XRD analysis based on the XRD modelling approach.

2.6. Modelling (00l) Reflection Profiles

The theoretical diffracted XRD intensity by a disordered lamellar structure is based on the Drits and Tchoubar [71] algorithm. Based on the expression of the diffracted theoretical intensity, detailed information going from the crystallite size to the elementary tetrahedral or octahedral sheet and its composition can be reached [71]. By optimizing several structural parameters (i.e., Z coordinates of the IS, exchangeable cation position and abundance, and water molecule distribution) during the modelling process, the agreement quality is improved [7274]. Rwp is used as a confidence factor to control this quality [74]. The existing mathematical arsenal behind the modelling approach as well as the details of the fit strategy adopted is exposed in several earlier works [7577]. It remains to be noted that (per half unit cell) the IS exchangeable cation position respects the provided literature data on similar samples [78, 79].

2.7. Layers Interstratification: Probabilistic Description

The layers stacking within crystallite is defined by the probabilities of succession of the different layer type, which is described according to the Markovian statistics model. The probability of the presence of one layer is only dependent on the types of layers that precede it. The relationship between the presence of a layer and the preceding types of layers corresponds to the degree of order of the interlayer or Reichweite (R). In soils, the most common Reichweite is R0 and R1; there are also interstratified R2 or R3 whose probability of occurrence of a layer depends on the two or three previous layers; the latter are not described in the soils but in the diagenetic series. Equations have been developed [70, 78, 79] to respond to clay interstratification problems. We thus find the probabilistic model where we represent (in a simple case) the stacking of two types of layers A and B of proportions WA and WB with the probabilities that layer A follows layer B (PBA). The fundamental equations governing this interstratification are as follows:

Using a combination of the equation, three essential limits for the layers stacking mode appear:(i)The random stacking R0 with PAA = WA (i.e., no stacking sequence is forbidden; the probability of appearance of a layer in a sequence does not depend on its abundance).(ii)Stacking with maximum order R1-MPDO prohibiting the succession of two layers in the minority.(iii)In the case where it is forbidden for two layers of different nature to follow one another (PAB = PBA = 0), the probability that two layers of the same nature do so is therefore equal to 1 (PAA = PBB = 1), and it is no longer interstratification question but of a physical mixture or total segregation.

It should be noted that the montmorillonite hydration heterogeneities can be treated as layer interstratifications at different hydration states within a clay particle depending on the relative humidity. This can constitute the intermediate hydration state (i.e., 0 W/1 W, 1 W/2 W, and 0 W/1 W/2 W). Hydration heterogeneity is linked to the presence of layer structure heterogeneity (proportion and charge localization of charges) [70, 78, 79].

2.8. Adsorption Measurement and Porosity Investigation

The relationship between mechanical stress and montmorillonite water content affects many geotechnical properties of unsaturated soils, including permeability, volume change, deformation, and shear strength. The main clay properties are mainly controlled by their internal and external surface [80, 81]. The investigation of montmorillonite porosity fate, having undergone structural transformations following the application of mechanical stress, seems essential to confirm quantitative XRD results. Also, the water retention properties are related to the pore structure and the level of external relative humidity experienced. With compaction/reswelling cycles and in situ study by variation of %rh, more layers of water molecules can either be absorbed on the surface of the pores or evacuated by drainage/diffusion by a new transport mechanism.

Adsorption-desorption nitrogen by porous media mainly targets the determination of the surface area (SA) and the pore size distribution (PSD). Montmorillonite powders are obtained by drying each sample. The porosity characteristics are entirely determined by nitrogen sorption [82]. The adsorbed gas amount gives a complete description of the porosity state and even the overall structure [8286]. Quantachrome NOVA 2000e series volumetric gas adsorption instrument, which is a USA automated gas adsorption system using nitrogen as the adsorptive, is used for the BET-specific surface area and PSD measurement. Experience setting conditions are as follows: (i) adsorption isotherms are obtained at 77 K and at P/P0∼0.95 (relative pressure); (ii) the SA and PSD measurements require the removal of adsorbed nitrogen and oxygen. This is carried out under reduced pressure (vacuum) at 100°C for 10 h. The desorption isotherms section, assuming a cylindrical pores model, is the basis of the application of the BJH method to define the PSD [8291].

3. Results and Discussion

3.1. The Oedometric Test: Compaction and Reswelling

Several parameters control the clay material response to mechanical stress. Montmorillonite fractions are characterized by high porosity and a total volume which essentially depends on the water content, the water saturation, and the water retention mechanisms. Other intrinsic/multiscale parameters ranging from the geological genesis to the individual layer nanometric properties influence the possible deformations [9295]. This work investigates the response of metal-exchanged montmorillonite specimens (i.e., exchangeable cations are Co2+ or Cd2+) to uniaxial mechanical stress ensured by an oedometric test (loading/unloading). The applied stress is confirmed by investigating the void ratio e modification giving the applied pressure (Figure 3). The reswelling curve is obtained following a release of the applied pressure (Table 1). At first sight, the compaction/reswelling curves show the deformation irreversibility process and the nonlinear process, whatever the exchanged cation nature, which agrees with [72, 91] works. This behaviour finds an explanation in [64], which attributes the observed nonlinearity of the consolidation characteristics to the correlation between the consolidation stress () and the void ratio (e).

The observed nonlinearity during the loading/unloading process concerns only a single stress cycle which resulted in remarkable fluctuations in the void ratio values. The microobjective here is to create material fatigue to examine its hydric behaviour afterwards by varying %rh under a controlled atmosphere. For the two samples, a gap separating the void ratio values in the absence of applied pressure (pressure = 0 bar) is partly explained by the coexistence of two possible scenarios, respectively, irreversible energy losses and interparticle friction, which is an intrinsic phenomenon in lamellar structures and which occurs during the progressive alignment of the particles/crystallites by increasing constraint strength [87, 88].

To understand the obtained void ratio values fluctuations during the reswelling process, the effective stress parameter, the swelling index (Cs), and the compression index (Cc) are introduced. Indeed, [64] uses the constant rate of strain CRS test as a crucial parameter to explore the unloading behaviour of clay minerals. For that, the pore water pressure excess during the unloading sequence is modelled by an assumption based on a cubic polynomial equation. Indeed, the Cc or Cs variations are very sensitive to the consolidation time (at constant load before the unloading) consideration. Along the compaction/reswelling sequence and the global pressure range (0 to 30 bar), in the case of the SWy-Co sample, . This value decreases in the case of the SWy-Cd sample to .

On the other hand, a maximum gap is observed for an applied pressure of 15 bar and for which jumps to 7.09 for Co2+. For Cd2+ cation, the maximum Cs and Cc fluctuations are observed at 10 bar when jumps to 5.52. For a maximum applied pressure intensity (zone on Figure 3) and whatever the nature of the exchangeable cation, a similar behaviour of Cs and Cc is observed during the loading sequence from 15 bar. During the unload sequence (reswelling) and from 30 bar applied pressure values, a divergence between the evolution of Cc and Cs is observed.

This noted disparity is explained by the process irreversibility at the nanoscopic scale, where the layer organization inside the crystallite depends essentially on the composition of the interlamellar space (IS). The intrinsic IS configurations are affected by the nature of the exchangeable cation, which afterwards will affect the layer stacking and the crystallite geometry. In fact, the diffuse ion layers around the particles lead to interparticle repulsion ascribed to the osmotic activity of the ions. The existing repulsion force decreases when increasing, respectively, the distance between clay particles, electrolyte concentration, and exchangeable cation valence. It should be noted that the obtained data for the “free” swelling state is realized when the sample is constrained only by the friction of the piston.

3.2. XRD Analysis of Unstressed Samples

The best agreement between the experimental and theoretical XRD profiles for the stressed samples is summarized in Figure 4. Qualitatively, a classical homogeneous 1 W layer phase is observed in the case of SWy-Na [67]. This hydration homogeneity is confirmed by the low values of FWHM and ξ parameter (Table 2) [74]. For the SWy-Co and SWy-Cd samples, an interstratified 1 W/2 W hydration character is observed throughout the entire exploited angular range. The initial d001 (Å) value increased (Table 2) towards 13.95 and 13.69, indicating probably the achievement of the cation exchange process and a new IS configuration. For both cases, the XRD profile geometry is characterized by an irrational 00l position (ξ parameter) with a large 001 refection (FWHM (2θ°)) = 1.02 and 1.23). The 001 reflection is accompanied by a shoulder towards 7.48 (2θ°) attributed to 1 W hydration phases (d001 ≈ 12 Å). Generally, the observed asymmetric 001 reflection reflects the coexistence of several hydration states within the stack and/or an incomplete (partial) cation exchange. Two possible interpretations arise in this case. Either the cation exchange capacity (CEC) is fully saturated by the exchangeable cations, or a minor fraction of the starting Na+ cation persists in the IS with partial cation exchange. Compared with the work of [5860] and considering the ionic potential affinity of the two cations with optimized, verified, and reproducible experimental conditions, obtained results in the case of Cd2+ and Co2+ cation are in concordance. More details can be given later by quantitative XRD analysis. Quantitative XRD investigations are the solution to decide with the predictions of the qualitative analysis [70, 76]. Indeed, in the case of interstratified structures, several hypotheses are necessary, whose resolution requires a more in-depth theoretical approach. In the case of the SWy-Na sample, an agreement between the experimental and theoretical profile is obtained with an Rwp = 3.89% (Table 3). The position of the sodium and its hydration sphere at the medium of the IS respects the literature data [57]. On the other hand, the experimental profile is reproduced theoretically through a randomly distributed two phases, qualitatively inapparent, between a major contribution of the 1 W state (≈75%) and a minor contribution of (≈25%) relative to the 0 W layer.

The 00l reflections modelling approach in the case of the SWy-Co sample supposes the coexistence of three mixed layers structure MLS saturated with the Co2+ cation. Theoretical models exhibit variable hydration states (e.g., 1 W, 1 W/2 W, and 1 W/2 W). The used MLS (from the best agreement) are decomposed into several layer populations (Table 3) assembled according to a specific succession probability law. The structural parameters released from theoretical models are summarized in Table 3. The MLS are obtained by weighted layer-type populations, which are expected to have, respectively, identical chemical composition, identical layer thickness, and identical Z coordinates of the atoms [67, 9698].

The best fit of the experimental XRD pattern (Rwp = 4.13%) (Figure 4) in the case of the SWy-Cd sample is obtained through theoretical models mainly composed of an interstratified structure using three MLSs including various relative proportions of the hydrated layer (Table 3). The foremost remark in this case is the presence of a reasonable 1 W phase weight (9.04%) attributable to the saturation of the material CEC by starting Na+ cation. This agrees with the starting assumption concerning the observed experimental 001 reflection shoulder at 7.48 (2θ°). This result calls into question the achievement of the cation exchange process despite compliance with the reproducibility of the experimental cation exchange protocol.

3.3. XRD Analysis of Mechanically Stressed Samples
3.3.1. Qualitative Investigation

The effect of the applied mechanical stress on the hydration properties of Cd2+ and Co2+ exchanged montmorillonite is addressed by varying the %rh rate under a controlled atmosphere. The goal is to see the structural evolution of the studied specimens under extreme %rh values (3%rh and 97%rh) (after environment equilibrium), which then allows us to simulate real IS changes affected using clay as a geological membrane. Experimental XRD patterns obtained in the case of SWy-Cd at 3%rh show a d001 = 12.23 Å indicating probably a 1 W hydration state. An interstratified hydration trend is confirmed by the elevated FWHM and ξ parameter value (Table 2). Indeed, a left-handed profile asymmetry towards the low angles (2θ = 5.84°; d001 = 15.04 Å), which is ascribed to a minor 2 W hydration state phase, was observed despite the conditioned dry environment. The d001 basal spacing value shows conservation of the 1 W hydration state regardless of the experimental equilibrium achieved at 3%rh (one hour of equilibrium with its environment before starting the XRD recording).

For the SWy-Co at 3%rh, the sample seems more sensitive to the %rh level drop, and the d001 = 9.27 Å indicates a basically dehydrated 0 W state. The observed result is a real manifestation not only of the %rh fluctuations effect but also of the mechanical stress already applied through the compaction and reswelling cycle. The conclusion that can be drawn by qualitative comparison between these two samples (having undergone the same stresses) lies in the modification of the IS organization. Indeed, the basal spacing shift is directly related to the new IS organization, which is very sensitive to the electrostatic and chemical bonding forces governing the IS equilibrium.

The IS water molecules insertion/release mechanism is interpreted by a hypothesis based on three freedom degree fluctuations, including, respectively, the positions, the abundances of the IS species, and the associated layer stacking probabilities. All cited structural parameters, treated individually or coupled, complicate and/or facilitate the water molecule release process, which can be established when decreasing the %rh rate. This hypothesis is since interpreted in the case of the SWy-Cd sample by the free water transformation (supporting cation exchange process) from the surrounding environment into structural water (intrinsic layer composition), incorporated in the internal and external pore surfaces, which complicates its loss by dehydration [60, 68]. It is noted that, in the case of the SWy-Co sample, an uncomplicated release of an IS water molecule is observed.

At 97%rh, SWy-Cd and SWy-Co are characterized by an asymmetric 00l reflection profile (Figure 5). Both samples present a high hydration state with d001 equal, respectively, to 17.38 Å and 17.32 Å, indicating probably a transition to the 3 W hydration state. A shoulder towards the small angles (d001 = 20.54 Å), which is attributed to a beginning 3 W⟶4 W hydrous transition state, is observed. An interstratified global character is confirmed by the high values of FWHM and the rationality parameter.

3.3.2. XRD Profile Modelling Approach

The best agreement between theoretical and experimental XRD patterns obtained at a variable %rh rate is confirmed by an acceptable Rwp that did not exceed 5.71%. Optimum structural parameters used to imitate experimental XRD profile are reported in Table 3.

The theoretical IS water molecule distribution and the compensator cation position follow a discrete distribution consistent with [57, 67, 74] works. For each accepted theoretical model, several layer types with variable stacking modes are used to improve the agreement. The partial segregation (R1) is the main MLS distribution adopted for all studied samples. At 3%rh, the experimental profile of the SWy-Cd sample is reproduced by mixing three MLS mainly composed of a major 1 W layer fraction (saturated by Cd2+ cation). For the same sample, a minor proportion of the Na+ saturated layer is used. Generally, the absence of a monohomogeneous 1 W phase at 3%rh and an R1 stacking mode type is the main conclusion. In the case of the SWy-Co sample, which presents a 0 W hydration state at 3%rh, the diffracted intensity is reproduced by two MLS (0 W) mainly saturated by Co2+ and Na+ cations (Table 3). At 97%rh and regardless of the exchangeable cation nature, the 0 W and 1 W phase contribution in the MLS disappears. The experimental profiles are entirely reproduced by a three-MLS combination characterized by variable weighted layer population. For SWy-Co, a minority bihydrated phase related to Na+ cation is introduced, and the major layer fraction saturated by Co2+ cation (extremely hydrated 3 W or 4 W water layers) is used to obtain a good fit.

For SWy-Cd, the structure is completely matched by combining three MLS. Indeed, theoretical models contain a major contribution of 3 W/4 W phase, and the sample retains its interstratified character, which is in accordance with the qualitative description mentioning an interstratified phase. Among the strengths of the modelling approach, the accurate structural heterogeneities identification is related to the coexistence of variable layer populations within the clay particle. Subsequently, it allows us to show the structural changes affecting the IS after applying stresses (Table 3).

3.4. Adsorption and Porosity Features

The clay properties are mainly monitored by their internal and external surface. The total surface includes the external surface, between the clay particles, and the inner surface, relating to the IS. Internal and external layer surfaces are affected during the compaction and reswelling process. The surface modification will probably influence the chemical and physical bonds with the IS content and the adsorption properties. The layer stacking distribution is affected by the applied mechanical constraint (results obtained from XRD analysis), and a possible consequent porosity disturbance may be conducted. In fact, after stress release, the configuration theoretically leans towards a new equilibrium, which differs from the starting configuration [99, 100].

The montmorillonite porosity is still poorly characterized because of the difficulty of visualizing hydrated samples in their original condition. The potential relation between porosity, CEC, crystallite shape/size, and applied mechanical stress is assured by the correlation of BET adsorption measurement and BJH pore size distribution analysis. Single-point and multipoint BET methods (based on the surface area (SA)) are used to assess the possible chemical transformation affecting the external clay layer surface. This purpose is achieved from nitrogen isotherms [90]. To measure the average pore diameter for each sample, a BJH method is applied. The obtained results from adsorption measurement and porosity for all studied samples are summarized in Figures 4 and 7. Results show an increase in the calculated external surface for samples having undergone mechanical stress, whatever the nature of the exchangeable cation and the constraint type (compaction or reswelling).

By comparing SA value evolution, after applying mechanical stress, a logical values boost is observed for the stressed samples, whatever the applied %rh rate. This is probably interpreted by a crystallite exfoliation tendency which results in a decrease in the average number of layers per crystallite following the layer cohesion damage. Also, it can be interpreted by the diffusion coefficient alteration in relation to the packaging density of the montmorillonite. For saturated CEC by Co2+ cations, high SA values are obtained, suggesting the effect of the exchangeable cation nature on the new IS configuration (Figure 6). The average nanopore diameter (Figure 7) respects the same trend with values boost after stress regardless of the cation exchanged.

For all examined samples, an isotherm of nitrogen adsorption- (ads-) desorption (des) was presented (Figure 8). A type II adsorption isotherm [84, 86, 90] is closer to the shape of the obtained curves. This allocation is due to the mesoporous texture described by coexisting of large pores with nanopores, which agrees with the pore sizes distribution for the porous materials. Large pore sizes obtained in the case of stressed SWy-Co sample are justified by the coupling of three essential parameters, respectively, the loading/unloading cycle, the %rh fluctuations (which constitutes a hydrous constraint), and the nature of the exchangeable cation. The Co2+ ions exchange promotes the exfoliation process and thereafter increases the porosity degree.

The SA concept does not permit a full textural description extending from 2D to 3D in the case of lamellar structures. The failures and limitations of SA results interpretation involve improvements established on the PSD analysis to access IS information. The PSD analysis is directed based on several approximations summarized in Table 4. The V-r (cumulative pore volume versus pore radius) plots or the mesopore PSD are given in Figure 9. A considerable disparity between SWy-Co and SWy-Cd sample pore size evolution was obtained. For unstressed and stressed samples, the saturated SWy-Co specimen varies in a very close way to the SWy-Cd sample, provided that the diameter of the pores does not exceed 3 nm. For higher values, a divergence appears with a “gap” in favor of Co2+ cation, which reaches 0.0068 mL/g (Figure 9). The obtained V-r variation proves, initially, the achievement of the cation exchange process, which is confirmed by specific V-r curve trends yet showing the same appearance. In addition, an increase in the V-r “gap” for stressed samples seems like the direct effect of the application of mechanical stresses. The V-r curve radius derivative for all samples is given in Figure 10.

The maximum mesopore volume (0.00620 mL/g) reached at 3.99 nm is obtained for the stressed SWy-Co sample. In this case, the determined V-r values are far greater than those relating to SWy-Cd. The gap between the two stressed samples exaggerates until it reaches gap values of 0.0030. This variation is consistent with the XRD modelling profile, which predicts layer exfoliation and presence of a highly hydrated rate (4 W) for both SWy-Co and SWy-Cd. Commonly, the applied mechanical stress affected the volume of the mesopores, which are in any case greater than those determined for the starting sample. For all samples, the mesopore radii varied between 1.598 and 16.828 nm (Figure 10). The porosity investigation using the BET-BJH method confirms obtained results from XRD modelling approach indicating, respectively, high hydration rate layer exfoliation trends and crystallite size fluctuations induced by stress.

The maximum mesopore volume (0.00620 mL/g) is obtained for the stressed SWy-Co sample. The determined V-r values, in the case of SWy-Co, are far greater than those relating to SWy-Cd. For the stressed SWy-Co sample, the maximum mesopore volume is 0.00620 mL/g and obtained at 3.99 nm. The gap between the two stressed samples exaggerates until it reaches gap values of 0.0030. This variation is consistent with the XRD modelling profile, which predicts layer exfoliation and presence of a highly hydrated rate (4 W) for both SWy-Co and SWy-Cd. Commonly, the applied mechanical stress affected the volume of the mesopores, which are in any case greater than those determined for the starting sample. For all samples, the mesopore radii varied between 1.598 and 16.828 nm (Figure 10). The porosity investigation using the BET-BJH method confirms obtained results from XRD modelling approach indicating, respectively, high hydration rate layer exfoliation trends and crystallite size fluctuations after stress.

4. Conclusion

This work investigates structural impact, IS response, quantitative CEC instabilities, and porosity alteration of metal-exchanged montmorillonite (Co2+ and Cd2+) according to the mechanical constraints created at the laboratory scale. Structures heterogeneities, hydration behaviour, and relationships between macroscopic stresses and microscopic hydration properties have been disclosed through the XRD profile modelling approach and the adsorption measurement. Obtained results demonstrate the following: (i) The void ratio e deviation increases along the compaction/reswelling process regardless of the exchangeable cation nature. (ii) At a low extreme %rh rate (3%), a clear and easy 0 W dehydration process is observed for the Co2+ cation, which is not obvious for the Cd2+ cation, which keeps a 1 W hydration state despite the mechanical disturbance. (iii) A highly 3 W/4 W hydrated state appears for the two studied samples brought to 97%rh, which is explained by a reduction in the average number of layers per crystallite evidence of an exfoliation tendency and the ease of filling of the IS. (iv) At 97%rh, mechanical damage induces, respectively, a clear interstratified character, segregation layer stacking trends, sensitive layer-type junction probabilities, and an increase in the used MLS. (v) The adsorption measurement results confirm the exfoliation layer trends for stressed samples. (vi) The divergence of V-r fluctuations proves the achievement of the cation exchange process, which is confirmed by the specific conduct of the V-r curves despite having the same appearance. (vii) The mechanical stress effect is reflected in the V-r curve by a high “gap” for the stressed samples. (viii) The mesopore V-r shows raised values in the case of Co2+ cation, indicating a new IS arrangement (compared with the starting sample) allowing better water profit for the SWy-Co sample [101103].

Data Availability

The generated and analysed data during the current study are included within the article and can be obtained from the corresponding authors upon reasonable request.

Conflicts of Interest

The authors declare no conflicts of interest.

Authors’ Contributions

Walid Oueslati (WO) was responsible for conceptualization, methodology, investigation, validation, formal analysis, visualization, data curation, resources, supervision, writing, editing, reviewing, and project administration. Chadha Mejri (CM) was responsible for investigation, visualization, and resources. Abdesslem Ben Haj Amara (ABHA) was responsible for methodology and resources.

Acknowledgments

The authors acknowledge the assistance provided by the editor.