Abstract

This work reports characterization of the possible effects that might distress the hydration properties of Mg-exchanged low-charge montmorillonite (SWy-2) when it undergoes external environmental solicitation. This perturbation was created by an alteration of relative humidity rates (i.e., RH%) over two hydration-dehydration cycles with different sequence orientations. Structural characterization is mainly based on the X-ray diffraction (XRD) profile-modeling approach achieved by comparing the “in situ” obtained experimental 00l reflections with other ones calculated from theoretical models. This method allows assessing the evolution of the interlayer water retention mechanism and the progress of diverse hydration state’s contributions versus external strain. Obtained results prove that the hydration behavior of the studied materials is strongly dependent on the RH sequence orientation which varied over cycles. The interlayer organization of Mg-exchanged montmorillonite (i.e., SWy-2-Mg) is characterized by a heterogeneous hydration behavior, which is systematically observed at different stages of both cycles. By comparing the interlayer water process evolution of Mg-exchanged montmorillonite with the observed SWy-2-Ni sample hydration behaviors, a same hysteresis thickness characterized by obvious fluctuations of interlayer water molecule abundances is observed. Nevertheless, in the case of Hg and Ba-saturated montmorillonite, the retention water process versus the applied cycles was steadier comparing with Mg ions.

1. Introduction

Smectites are swelling clay minerals that belong to the family of the phyllosilicates 2 : 1, which naturally occur in both terrestrial and marine environments [1] where they often represent the most effective components. In fact, by dint of its significant intrinsic physical and chemical properties, mainly the high specific surface areas (up to 760 m2·g−1), the cation-exchange capacities as well as the high sorption efficiency of cations, and its strong mechanical stability [2], these materials were used in many extensive exploitation. Indeed, it plays an important role in many geological processes, the phenomena of the petroleum migration, greenhouse gas sequestration, oilfield [3], and in engineering. Moreover, for last decades, the smectites were widely used as crucial components for elaboration of natural barriers to isolate the hazardous wastes and for the removal of heavy metal cations from various effluents of industrial and wastewater treatment [412] and also proposed as geotechnical barriers in many nuclear waste disposal concepts in order to retard the potential transport of radionuclides towards the biosphere [1315].

Notwithstanding the diverse beneficial effects of smectites, the hydration behavior of these mineral is very sensitive to the change of the environmental surrounding conditions (i.e., temperature, pressure, and RH), which may influence the stability of the clay microstructure therefore the stability of the geotechnical barrier. Hence, the well understanding of swelling smectites properties (the hydration and the dehydration mechanism) is of paramount importance for many natural processes and for such applications.

The water-metal-smectites system was the subject of diverse studies, and it was extensively studied at different scales using several characterization methods [1626]. At nanometer scale, the swelling process corresponds to the presence of hydrated chemical species on the interlayer space. This process depends on several factors including composition such as the layer charge amount, the charge location and distribution [2730], the interlayer cation type, the interlayer cation valence, and the hydration energy [31] and also on the environment conditions such as the relative humidity, temperature, and H2O pressure [3237]. The pioneering studies focusing the crystalline swelling was performed using X-ray diffraction technique. By following the evolution of d001 basal-spacing value, different hydration states were defined with the insertion of 0, 1, 2, or 3 planes of H2O molecules in the interlamellar spaces leading to growth of the spacing values as function of RH% [28,3840].

Furthermore, the development of the XRD profile-modeling procedures and the methodology proposed by [41] allows the characterization of the structural modifications that occur during smectite swelling process. In addition, by using the XRD modeling approach, the quantification of hydration heterogeneity becomes easier especially when different layer hydration state types coexist in the same smectite structure. Several studies used this modeling approach tools to study the hydration properties of different smectite sample types and demonstrate systematic hydration heterogeneity whatever the interlayer cation, the RH%, the amount, the layer charge location, and the surrounding temperature [40,4246].

According to this method, the present work focuses the damage that may affect the structural properties of the host materials and the hydrous behavior of an Mg-saturated smectite (i.e., montmorillonite: SWy-2) when it was submitted to a continuous variation of an environmental surroundings condition (i.e, variable RH%). This environmental solicitation is performed by varying “in situ” the RH% in reverse sequence orientation upon two hydration/dehydration cycles. Upon the applied cycles and under controlled atmosphere, a complex progress of the interlamellar space configuration of SWy-2-Mg sample is followed and quantified. After that the obtained results are compared with the earlier studies related to the same studied specimen which is saturated with other bivalent cations (Hg, Ba, and Ni) and submitted to the same climatic changes over the hydration-dehydration cycles for the purpose of discriminating the effect of the ionic potential on the crystalline swelling.

2. Materials and Methods

2.1. Host Materials

A dioctahedral smectite SWy-2 originated from bentonites of Wyoming (USA) is selected for the present study. Clay fraction is supplied by the Source Clay Minerals Repository Collection of the Clay Minerals Society and characterized by the following half-cell structural formula [47]:

The cation-exchange capacity (CEC) of this smectite is 101 meq/100 g where the charge deficit is majority resulting from cationic substitutions in the octahedral sheet and extremely limited tetrahedral ones.

2.2. Sample Treatments

To saturate all exchangeable sites by homoionic cations (Na+) and to guarantee better colloidal dispersion, a pretreatment of the host material is carried out in order to prepare a Na-rich montmorillonite suspension. This aims is based on a classical protocol of an exchange process [48] which consists of dispersing ∼20 g of solid in ∼200 ml of NaCl solution (1 M) and stirring mechanically for 24 h. A SIGMA laboratory centrifuge is used for the separation of the solid fraction at 4000 rpm speed. These steps were repeated five times to ensure saturation of all exchangeable sites by Na+ cations. Excess chloride was removed by washing with distilled water five times, and the separation of the solid-liquid was performed by centrifugation at 8000 rpm speed. The same ionic exchange procedure was followed to prepare Mg-rich montmorillonite suspension using MgCl2 solution (1M), and the final obtained clay suspension was labeled SWy-2-Mg.

Two oriented slides were prepared for the obtained samples (SWy-2-Mg) to be analyzed by XRD technique. For that the specimen suspension was deposed on a glass slide and then dried at room temperature for 24 hours to obtain an air-dried preparation [49].

2.3. “In Situ” XRD Analysis

The “in situ” XRD patterns produced by SWy-2-Mg sample are obtained upon two overturn hydration-dehydration cycles. These cycles were created by varying gradually the RH% in reverse orientation with a 10% step. For that an Anton Paar TTK 450 chamber coupled with a D8 Advance Brüker installation (Cu-Kα radiation) equipped with solid-state detector and operating at 40 KV and 30 mA was used. The reflection setting diffractometer installation is equipped with an Ansyco rh-plus 2250 humidity control device used to vary manually the RH%. The hydration-dehydration cycle orientations were detailed in Figure 1 and can be resumed as follows.

Both applied cycles were performed in three successive processes starting from 40% RH, which coincides with the relative humidity value of the room condition (297 K and ∼40% RH).

The first cycle is divided into three domains summarized as follows.

The first one (I) consists of a hydration process realized by increasing the RH% from 40 to almost saturated condition (80%) followed by a dehydration procedure (II) assured by the decrease of the RH% towards extremely dry condition (10%). Finally, a rehydration process (III) was performed to a second return to 40% RH (Figure 1).

The sequential RH orientation was accomplished in an inverse way for the second cycle. Indeed, starting from 40% RH, a dehydration process was performed decreasing the RH% to 10%, followed by a hydration procedure reaching 80% RH and to finish with a second dehydration process by lessening the RH% to 40% (Figure 1).

Over these applied cycles, experimental XRD patterns were registered, in situ, every 10% at the fixed relative humidity condition values where the usual scanning parameters were 0.04°2θ as step size and 6 s as counting time per step over the angular range 2–40°2θ. In total, sixteen experimental patterns were recorded per cycle, and for all obtained 001 reflection, quantitative and quanlitative XRD analyses are performed.

2.3.1. Qualitative XRD Pattern Investigations

A primeval interpretation about the hydration states of the studied samples can be obtained from a qualitative XRD analysis of the experimental profile. This information was deduced through the d001 basal-spacing values and a description of the 001 reflection profile geometry (pics symmetry and/or asymmetry). In addition, the correlation between the calculated parameters including the full width at half maximum intensity (FWHM) for the 001 reflections and ξ parameters [42, 50] can supply information about the hydration character of the studied samples (homogenous or interstratified). Nevertheless, the qualitative examination cannot provide information about structural transformation related to the position and organization of H2O molecule and the exchangeable cations along the axis. In addition, it is impossible to distinguish the nature and the relative contributions of different hydration phases at different RH values varied over the cycles. Thus, a quantitative evaluation of diverse changes which occurs within the smectite structure was required to accomplish the aims of the studies.

2.3.2. XRD Profile-Modeling Approach

The quantitative analysis is based on the XRD profile-modeling approach in order to propose theoretical structural models estimating, respectively, the gradual evolution of the interlamellar space content versus the hydration sequences, the nature of the different layer types coexisting within crystallites, and their proportion and their structural composition at different stages of both the applied cycles.

This method consists of adjusting the experimental XRD patterns (001 reflection) to theoretical ones where the calculated intensity is based on the algorithms developed by [41]. The theoretical matrix formalism was detailed by [41]. The used Z atomic coordinates of the interlayer space correspond to those proposed by [41]. The abundances of the different types of layers (Wi), the mode of stacking of the different kinds of layers, and the mean number of layers per coherent scattering domain (CSD) [51] can be determined also through XRD profile-modeling approach. The layer-type stacking is described by a set of junction probabilities (Pij) where the relationships between these probabilities and the abundances Wi of the different types of layers were given by [41]. A detailed description of the fitting strategy was detailed in the work of [37, 40]. Indeed, XRD pattern-modeling was performed assuming the possible presence of different layer types. These different layer types correspond to the different hydration states commonly reported in smectites as a function of relative humidity. In the fitting process, we have introduced dehydrated layers (0W layers, layer thickness at 9.6–10.0 Å), monohydrated layers with one plane of H2O molecules in the interlayer (1W layers at 11.5–13.0 Å), bihydrated layers with two planes of H2O molecules in the interlayer (2W layers at 13.9–15.8 Å), and trihydrated layers (3W layers at 18.0–18.5 Å).

3. Results

3.1. Quantifying Interlamellar Space Content during the First Cycle
3.1.1. Qualitative XRD Analysis

All experimental XRD patterns produced by the SWy-2-Mg sample during the first cycle are presented in Figure 2 with the calculated profile, obtained using the corresponding contributions of the various mixed-layer structures (MLSs).

A qualitative investigation was performed providing preliminary information about the hydration property evolution. In fact, a homogenous hydration character is observed at the highest RH range extending between 70% (hydration process) ≤ RH ≤ 60% (dehydration process). This description is confirmed by the low value of the calculated FWHM, the ξ parameter (Table 1) of the 001 reflection positions over the RH fields. On the contrary, a heterogeneous hydration behavior was deduced at the lowest RH range justified by the irrationality, for all measurable reflection positions, characterized by high ξ parameter values (Table 1).

The gradual evolution of the d001 basal-spacing values as function of RH% (Figure 3(a)) shows a clear hysteresis between 40% ≤ RH ≤ 10%. For this RH range, the structure is dominated by an interstratified hydration behavior between 1W and 2W layer types. However, the variation of the d001 spacing between 80 and 40% RH tails the same roads which can be interpreted by a hydrous stability on the interlayer spaces.

3.1.2. Theoretical Models and the Hydration Properties

The structural parameters used to reproduce experimental patterns of SWy-2-Mg as a function of RH% were regrouped in Table 2. Results derived from the quantitative XRD investigation demonstrate that the studied structure (i.e., SWy-2-Mg) changes, all over the cycle, with an interstratified hydration behavior. Indeed, all proposed models are described by a main structure composed by diverse hydration state contributions (0W, 1W, 2W, and 3W) at different RH values (Figure 2).

These heterogeneous hydration behaviors can be explained by sequential transitions between different hydration states induced by the continuous variation of the RH%. The evolution of different relative layer-type proportions as a function of the RH% (Figure 4) shows a continuous diffusion of H2O molecules in the interlamellar spaces during the hydration process leading to progressive and continuous 1W → 2W and 2W → 3W transitions.

In fact, at the beginning of the cycle (40% RH), the studied sample (SWy-2-Mg) was characterized by a main structure composed by 21% of 1W layer types with major contribution of the bihydrated ones (79%). By reaching 80% RH, a different configuration composed by 7%, 65%, and 28% attributed, respectively, to 1W, 2W, and 3W hydration states is manifested (Figure 4). On the contrary, along the dehydration process, the decrease of the RH values from 80% to 30% leads to a gradual emptying of the microscopic porosity corresponding to the reduction of the interlamellar water molecule abundances. Indeed, successive 3W → 2W and 2W → 1W transitions are observed over this RH range, where the structure is composed at 30% RH by 19%, 75%, and 6%, respectively, for 1W, 2W, and 3W layer types. An obvious transformation on the structural composition and a notable change on the hydration behavior were distinguished towards the lowest RH domain starting from 20% RH (rehydration process) to the end of the cycle at 40% RH (rehydration process). In fact, fast emptying of the interlamellar spaces from water planes was noted by decreasing the % RH rates during the dehydration procedure where the 2W layer-type contribution decreases rapidly with a complete disappearance of the 3W layer from the structure (Figure 4). Over this RH range, fast increasing of the 1W layer-type contribution was noted, which dominates the structure until the end of the cycle (Figure 4).

3.2. Quantifying Interlamellar Space Content during the Second Cycle
3.2.1. Qualitative XRD Investigation

Experimental XRD patterns obtained under controlled RH condition along the second cycle are represented with the respective contributions of the various MLSs used to calculate profiles in Figure 5.

The correlation between FWHM and the ξ parameter (Table 3) proposes that over a large RH range extending between 30% (dehydration process) ≤ RH ≤ 70% RH (hydration procedure), the structure is characterized by an interstratified hydration character. However, the study suggests homogenous hydration behaviors with decreasing the RH values from 80% to 40% over the second dehydration procedure (Table 3).

Evolution of the layer thickness (d001) as function of the RH%, along the second cycle, was characterized by the appearance of a clear hysteresis at the RH range extending between 40% and 10% (Figure 3). Indeed, a fast shift of the d001 spacing values from 15.39 Å at 40% RH to 12.90 Å at 10% RH was observed. However, a slow d001 value progress is noted by inversing the RH% orientation. Over these RH domains, the d001 is attributed to an interstratified 2W-1W hydration state.

3.2.2. Theoretical Models and the Hydration Properties

Optimum structural parameters used to fit experimental XRD patterns in the case of the second cycle are summarized in Table 4. The quantitative XRD investigation demonstrates that all calculated theoretical models which allow reproducing the 001 reflections are characterized by heterogeneous hydration states composed of two different MLSs with diverse relative proportions of layer types (Figures 5(a) and 5(b)).

Evolution of different layer type contributions as a function of RH is detailed in Figure 6. With decreasing the RH% during the first dehydration procedure, a slow transition from the 2W to 1W state was observed between 40% RH and 20% RH where the smectite crystallite presents major contributions of the 2W layer types.

A notable increase of the 1W layer proportion, towards the lowest RH condition (i.e.,10% RH), is detected. At this RH rate, the structure is reproduced by 70% of the monohydrated states (1W) and 30% of the bihydrated ones (Figure 6).

Progressive intercalation of water planes is noted during the hydration procedure. This accelerated intercalation is achieved by a continuous transition from the 1W to the 2W layer types between 20 and 40% RH. Reaching 50% RH, the insertion of the H2O molecules becomes more easier; indeed, the 3W state appears at this stage and persists in the structure for a wide RH field extending between 50 and 60% RH (dehydration procedure) (Figure 6). Over these RH domains, the structure grows gradually containing different contributions of the three layer types including 1W, 2W, and 3W with a clear dominance of the bihydrated phases (2W). During the second dehydration procedure performed by decreasing the RH% from the highest RH value (i.e., 80% RH) to the starting point (i.e., 40% RH), a gradual emptying of the interlamellar spaces from water is noted. Theoretical models propose structures characterized by the major proportion of the 2W layer types which kept the highest proportion until the end of the cycle (Figure 6).

4. Discussion

4.1. RH, Interlamellar Water Amounts, and Ionic Radius

The continuous variation of the RH values along both cycles automatically brings a sequential transition between different hydration states in the interlayer spaces of SWy-2-Mg complex, which logically leads to a continuous change on the interlayer water amount and distribution.

Figure 7 represents a comparison between the interlamellar water molecule content evolution along the first and the second cycle, which establishes the dependence of the hydration behavior progression on the sequence variation and orientation of the RH values clearly. In fact, a dissimilar progression of the water amount was perceived in three principal RH fields: Over the first RH section spreading between 50% RH and 80% RH, labeled (I), the interlayer water amounts retained in the structure were more important in the case of first hydration-dehydration cycle than in the second one. An analogous hydration behavior of the studied complex is observed in the short RH domain over the second domain (II) between 50 and 40% RH of both cycles. However, at the lower RH fields (III), quantitative results prove that the calculated interlamellar water amounts were more important in the case of the second cycle than the first one.

Such obtained result is in concordance with the previous work [5254] where the same studied montmoriollonite was saturated with other bivalent cations characterized by different ionic radius (Ni2+, Hg2+, and Ba2+) and have undergone the same “in situ” adsorption/desorption sequences. Results derived from modeling of X-ray diffraction (XRD) patterns related to the cited study have proven that the hydration mechanism at crystal scale of diverse samples was found to evolve gradually in different ways as a function of the applied cycles, thus proving the dependence of the progress of the swelling property of the montmorillonite on the sequence orientation, which means problem related to how the RH rates varied over the cycle.

On the contrary, the comparison between the development of the interlayer H2O molecule retained in the interlamellar spaces for the montmorillonite saturated by diverse bivalent cations (Ba2+, Ni2+, Hg2+, and Mg2+) upon the first cycle and the second one (Figure 8) shows an analogous hydration performance between the SWy-2-Ni and SWy-2-Mg. In fact, a clear hysteresis is observed in both cases indicating irreversible interlayer water progress during the hydration-dehydration cycles. However, location of Ba2+ and Hg2+ in exchangeable sites of the same matrix leads to more stability in interlamellar spaces of both structures; thus, the evolution of the interlayer water contents was more respected during the applied cycles. The assessment of these results revealed that the environmental solicitation performed by the continuous change of the RH% affects deeply the hydration performance of the host materials (SWy-2) especially in presence of bivalent cations with the lowest ionic radius in their exchangeable sites. In fact, unlike the cations with little size (Mg2+ and Ni2+), the huge size of Ba2+ and Hg2+ favors a decrease of the hydration heterogeneities degree and thus establish more stability of the crystal structure. Thus, the ionic radius of the bivalent compensator ions represents an intrinsic parameter which has important impact on the evolution of the hydration behaviors of the studied smectite.

The hydrous disruption that appears on the interlamellar spaces can be interpreted by the appearance of structural perturbation and new organization at the crystal scale versus the applied hydrous strain. Indeed, the hydration heterogeneity was explained in many studies by the disorder distribution of the surface charge sites which leads to heterogeneous structural composition responsible for such behaviors [55, 56].

Such structural perturbation is justified in the present work through quantitative XRD investigation. Indeed, all theoretical structural models suggest the coexistence of more than one mixed-layer structure (MLS) and propose continuous variations of the interlamellar spaces configurations whatever the interlayer bivalent cations nature (i.e., Mg2+, Hg2+, Ba2+,and Ni2+) and whatever the relative humidity rates that varied upon the applied cycles [5254].

However, the unsteadiness of the material behavior face to the environmental surroundings changes represent inconvenience when using such mineral as natural barriers in industrial wastes and radioactive treatment specially in long-live storage application.

Indeed, [45] demonstrates that the cation-exchange process of Na-rich montmorillonite in contact with solution containing bivalent cations (Cu2+ and Co2+) was affected when applied 15 hydration-dehydration cycles, created by continuous variation of RH% rates. In this case, results derivatived from quantitative XRD analysis prove that the hydrous strain was accompanied by an obvious structural change characterized by a decrease in the amount of exchangeable sites which affect the selective exchange process.

In this regard, taking into account the change of surrounding environmental condition is crucial to avoid its consequence on the microstructure stability of the geotechnical barrier especially in long term.

4.2. Hydration Hysteresis

The comparison between the average hydration hysteresis and the evolution of the calculated standard deviation (SD) for all specimens along both cases of cycles were reported in Figure 9. The SD is generally used to quantify the amount of variation or dispersion of a set of data values. The comparison shows that every sample was characterized by different average hydration hysteresis as well as their associated SD as function of the applied cycle. This result confirms the impact of the RH sequence orientation on the hydration behavior of the studied montmorillonite whatever the nature of the bivalent cation present in their exchangeable sites.

Moreover, the evolution of the SD parameter versus RH% over both cycles (Figure 10) indicate that for all studied complexes, the calculated SD values were more important for the major parts of the RH field in the case of the second cycle than in the first one. The high SD values can be interpreted by the appearance of important structural fluctuations and perturbation on the hydration properties over the second cycle comparing with the first one. However, when focusing the progress of the SD values’ evolution during the first cycle, the obtained results show that the calculated SD decreases gradually for different samples when increasing RH values, which indicates more structural stability trends for all studied complexes at the highest RH fields.

On the contrary, according to the literature [20, 28, 33], the classification of the divalent cations used in this work (based on their ionic potential) should respect the following order: Mg2+ > Ni2+ > Hg2+ > Ba2+. However, the exploitation of Figure 9, by following the evolution of different average hydration hysteresis as function of the %RH, indicates that this classification was low respected in major parts of the explored RH fields and in both cycle types. Indeed, at the exception of narrow RH fields extending between 10% and 30% during the second cycle (Figure 9), the classification varied arbitrarily between different cations and changed randomly from an RH field to others ones.

As a consequence of this last result, the ionic potential cannot explain or justify the hydration behavior evolution of the montmorillonite saturated with bivalent cations when submitted to continuous changes of the RH%. In this regard, the water retention mechanism and the crystalline swelling of this smectite became more complex phenomenon under this type of environmental surrounding condition change. In fact, several factors contribute simultaneously to govern the interlamellar hydration processes, and their combined impacts complicate the well understanding of the interlayer swelling process.

5. Conclusion

This work focuses on the detailed hydration behavior response of an Mg-rich montmorillonite when it undergoes an environmental solicitation created by continuous variation of the RH% along two different hydration-dehydration cycles. Quantitative analysis is mainly based on XRD-modeling approach. The evolution of the hydration properties of the studied sample is quantified and followed every 10% RH upon both cycles. The obtained results are compared with the hydration behaviors of Hg-, Ba-, and Ni-rich montmorillonite studied in the same environmental condition changes.

The main results obtained through quantitative XRD investigation shows the following:(i)The hydration behavior of the studied sample (i.e., SWy-2-Mg) strongly depended on the sequence orientation of the RH that varied over cycles.(ii)The proposed theoretical models describing the evolution of the structural properties suggests the coexistence of more than one MLS indicating the hydration heterogeneity character for the SWy-2-Mg whatever the RH% sequence orientation.(iii)The montmorillonite’s interlamellar water content growth was dependent on the nature of the bivalent exchangeable cations. In fact, the presence of Mg2+ as well as Ni2+ ions in the structure leads an irreversible interlayer water content process confirmed by the appearance of a hydration hysteresis. However, the location of cations with largest ionic radius like Ba2+ and Hg2+, in exchangeable sites, was accompanied by more orderliness of systems and decrease in the water-content fluctuation.(iv)The effect of the ionic potential parameter on the interlayer water retention mechanism, under controlled atmosphere (i.e, variable RH), is not justified and cannot give a logical explanation on the progress of hydration behaviors of different complex.

Data Availability

No data were used to support this study.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The results presented are a part of the Ph.D. thesis of Marwa Ammar realized at the PMLNMH (UR05/13-01) Faculty of Science of Bizerte, Tunisia, cosupervised by Dr. Walid Oueslati and Prof. Abdesslem Ben Haj Amara. Marwa Ammar acknowledges Dr. W. Oueslati for the fruitful discussions about smectite hydration, her main contribution in the XRD modeling approach, and the proof reading of the manuscript.