Abstract
Smectites are a group of minerals traditionally analyzed by thermal methods due to their exceptionally large adsorbed-water contents and the presence of OH groups, which makes them unique among all common soil- and rock-forming minerals. The dehydration reaction of smectite is a low-temperature endothermic effect that ends typically below 200°C. Although the removal of bulk interlayer water requires activation energy (Ea) of just above 30 kJ/mol, the removal of the last few H2O molecules attached strongly to interlayer cations requires Ea > 100 kJ/mol. Dehydroxylation is the loss of structural OH groups that proceeds as evolution of H2O molecules out of the smectite structure and occurs in the 300–900°C range. In trioctahedral species, dehydroxylation is combined with recrystallization and proceeds usually at > 700°C. In dioctahedral species, the temperature of dehydroxylation is controlled by the type of octahedral vacancy, having trans-vacant and cis-vacant distinguished by the boundary at ~ 600°C, and by the octahedral cation–OH bond strength, following the order Mg > Al > Fe. The Ea of dehydroxylation correlates linearly with the temperature of maximum dehydroxylation; from > 170 kJ/mol for Cs+-exchanged beidellite and nontronite, through ~ 300 kJ/mol in Mg-rich montmorillonite, to > 500 kJ/mol in trioctahedral saponite. Dehydration and dehydroxylation of smectites can be accompanied by a number of other phenomena, such as dehydrogenation or defluorination. At high temperatures, smectite amorphization and recrystallization occurs. Unless amorphized and/or recrystallized, smectites can undergo rehydration and rehydroxylation, which are opposite reactions to dehydration and dehydroxylation, respectively. This review discusses the details of the above-mentioned thermal reactions of smectites, focusing on thermogravimetric methods, evolved gas analysis, and structural alterations. Factors affecting the accuracy and precision of thermal analysis of smectite are discussed along with examples of best laboratory practices. The paper also provides the most recent description and critical evaluation of smectite reaction kinetics.
Introduction
Thermal analysis is the oldest instrumental method applied to clays and clay minerals. The first application of thermal analysis to classify and identify smectitic species (Le Chatelier, 1887) happened long before the discovery of X-ray radiation and the identification of the actual smectite structure using X-ray diffractometry (XRD). Right after the description of the smectite structure by Hofman et al. (1933), the already-developed differential thermal analysis (DTA) methodology quickly become a major tool in studying smectites, along with their classification and transformations (Grim & Bradley, 1940; Grim & Rowland, 1942a, b). Already in the early years of investigation, the nature of the main DTA peaks in a smectite heating curve had been identified correctly: the low-temperature endothermic DTA peak is caused by dehydration and the high-temperature endothermic DTA peak is caused by dehydroxylation (Fig. 1). These two reactions are followed by structure decomposition and recrystallization.
Thermal Analysis Methodology
Thermal analysis is a ubiquitous and widely used general term that refers to a number of analytical methods monitoring the change in a material’s property with temperature, usually measured over small intervals of time or temperature change. The principles and various methods have been described in numerous reviews and books (Byrn et al., 2017; Gabbott, 2008; Plante et al., 2009), including excellent, in-depth reviews devoted specifically to clay minerals (Emmerich, 2010; Stucki et al., 1990). Therefore, the present paper focuses only on the application of thermal analysis to smectites (also applicable to vermiculites and smectite-interstratified phases) and only to smectites in a form corresponding to their most common natural occurrence in rocks and sediments, i.e. having inorganic interlayer cations. Thermal analysis of smectites with adsorbed organic compounds is a field in and of itself with dedicated reviews (e.g. Xi et al., 2005; Plante et al., 2009; Yariv et al., 2011).
Among many thermal analysis methods (see Emmerich, 2010), a few have been used widely for smectites and related clay minerals, and are the most meaningful for their characterization: thermogravimetric analysis (TGA), differential thermal analysis (DTA), and its modern successor – differential scanning calorimetry (DSC). If a sample is analyzed simultaneously using more than one of the techniques (e.g. TGA and DTA), such a mode of analysis is often referred to as Simultaneous Thermal Analysis (STA). The evolved gas analysis (EGA), which allows determination of the composition of gases and vapors evolved from a sample during heating in a TGA or DTA/DSC analyzer, is a form of STA. Among the different means of analyzing gas evolved from clay samples, infrared spectroscopy (IR) and mass spectrometry (MS; also coupled with gas chromatography – GC-MS) are the most commonly used.
Thermogravimetric analysis presents weight change with temperature or time. In most cases this corresponds to the extent of a clay’s weight loss (in mg or weight %) during heating in an inert gas atmosphere, such as N2, Ar, or He. However, under specific circumstances, mass gain can occur during TGA experiments, i.e. when an oxidizing agent or water molecules are present in the sample’s environment. Common practice when reporting TGA results is to present a first derivative of the original mass change curve over time or temperature. Such a derivative is referred to as the DTG curve. Differential thermal analysis returns the temperature difference between the sample and an inert reference material under the same heating conditions. Similar to DTA, DSC measures the heat flow (J K–1 mol–1) needed to maintain the same temperature of the sample and the inert reference material, which enables the heat capacity to be obtained (Giese, 1990). All the methods mentioned require some kind of purge gas which, in general, is meant to remove the thermal reaction products. A moisture analyzer is a low-temperature alternative to the TGA instrument, capable of measuring weight change (usually due to dehydration) from room temperature typically up to 200 or 250°C, with no purge gas flow.
Any thermal analysis can be performed in two modes depending on the rate of heating: (1) isothermal – temperature is increased quickly to the desired value and kept constant for the remaining part of the experiment; (2) non-isothermal – temperature varies during the experiment. This mode is represented by, for example, heating or cooling with a constant rate to the final temperature. The most common way of performing thermal analysis is to use a predefined (isothermal or non-isothermal) temperature program. Controlled-rate thermal analysis (CRTA) is a specific case of non-isothermal or quasi-isothermal conditions in modern TGA instruments, where the sample’s reaction rate (e.g. weight loss per time) is set as constant, but the heating rate varies (Ding & Frost, 2002; Rouquerol et al., 2013). As the temperature controller keeps the reaction rate constant during CRTA, temperature change with time for a particular experiment is not known beforehand. Whereas DTA and DSC usually need the linear heating/cooling mode (with exceptions, see Boller et al., 1994), TGA can operate in any isothermal or non-isothermal mode, or a combination of both.
Measurement Accuracy and Precision
Thermal-reaction temperatures are commonly reported based on the position of the corresponding DTG or DTA peak. A number of factors affect the position of the peak, which can hamper straightforward comparisons among reaction temperatures reported in various studies (Emmerich, 2010). First, a positive correlation occurs between the DTG or DTA peak temperatures and the heating rate in non-isothermal conditions due to kinetic effects. Second, the peak temperature is dependent on the particular thermal analyzer’s temperature calibration, which should be based on at least two points in the range corresponding to the tested reaction. Curie temperature reference materials are commonly used for TGA instruments, whereas temperature and enthalpy of fusion are used in DTA and DSC instruments (Della Gatta et al., 2006). Theoretically, ~1 K (1°C) should be considered as a standard error in the thermal analysis of smectite in well-calibrated instruments. In practice, however, several factors may cause the actual sample temperature to deviate from the measured value even more (see Emmerich, 2010 for comprehensive review). As a result of a poorly crystalline structure, heterogeneous cation distribution, the particularly small size of its crystallites (both planar and crystallite thickness), and their diverse means of agglomeration, smectite samples do not return sharp, well-defined thermal events. These features, along with operator-dependent choices during data mathematical processing, all result in the real accuracy and precision of determining reaction temperature being in the range of a few °C.
Accuracy of the mass-loss measurement can be tested using a homogeneous material with a well-constrained mass loss in a given temperature interval, such as Ca-oxalate monohydrate (Bish & Duffy, 1990; Plante et al., 2009). Thermal drift (i.e. the profile of mass or heat flow change with temperature in the absence of any reactive material) should be measured for each thermal analysis instrument and a mode of analysis. Thermal drift can be determined on an empty pan or an unreactive sample with a well-calibrated instrument under conditions corresponding to the experimental sample analysis and the test results should be applied to correct an experimental sample curve. A well-calibrated TGA instrument and method should return the mass-loss precision in the range of 0.1%.
Thermal Reactions of Smectites
Dehydration
Dehydration is the loss of molecular water (H2O) adsorbed on or within clay-mineral particles. Upon drying, H2O molecules are removed sequentially following the strength of their interactions with the smectite surface. The entire physisorbed water condensed in meso- and micropores of particles is removed at room temperature under the flow of dry purge gas. Bulk (bound via H-bonds to other H2O molecules) smectite interlayer water is in equilibrium with the external moisture and, therefore, it is removed along with the physisorbed water albeit at a lower rate (Emmerich et al., 2018; Ferrage et al., 2010; Salles et al., 2009). Therefore, distinction between the bulk interlayer and the external, physisorbed water is hardly possible in thermal analysis as no temperature–time cutoffs are applicable. The removal of physisorbed and bulk interlayer water is completed much below 100°C. Due to the difference in solidification and fusion temperatures of physisorbed water, the DSC-based thermoporometry (freezing conditions, below 0°C) is the best temperature-differentiation method that can be applied to address the distinction between total interlayer and micropore capillary water (Grekov et al., 2019; Salles et al., 2010).
In a fully hydrated smectite, some water molecules form a hydration shell around the smectite’s interlayer cations. Once the bulk interlayer water is removed during dehydration, the remaining H2O molecules become organized around the interlayer cations in a fashion and number dependent on the cation’s hydration energy (Salles et al., 2015), expressed as the cation’s hydration enthalpy (Table 1). The larger the absolute value of the cation hydration enthalpy (see Table 1 for an explanation of ΔHo), the lower the probability for the cation to bind directly to the smectite’s siloxane surface, the stronger the interaction between water molecules and the interlayer cation, and, therefore, the greater the energy that is required to remove any remaining strongly bound water (Webb et al., 1986).
The cation’s hydration enthalpy is reflected by the pattern of smectite dehydration. Except for Li+-saturated smectites, the interlayer’s monovalent cations produce a single DTG dehydration peak below 100°C, whereas two DTG peaks or a high-temperature shoulder (up to ~200°C) usually correspond to high-enthalpy cations (polyvalent and Li+; Fig. 2; El-Barawy et al., 1986). The fewer the H2O molecules around interlayer cations, the greater the activation energy (Ea) of dehydration and the higher the temperature needed to remove the remaining H2O. In extreme cases, H2O dissociation leads to the hydrolysis of interlayer cations by OH that can remain in the smectite interlayer, even up to 500°C, which affects the course of dehydroxylation (Kuligiewicz & Derkowski, 2017).
NH4+-exchanged smectites show very low mass loss due to dehydration as the ammonium cation has a low hydration enthalpy (Table 1); however, the smectite undergoes an exothermic deammonization reaction at 300–400°C resulting in weight loss and leaving smectite in H+ form (Mackenzie & Bishui, 1958; Seiffarth & Kaps, 2009).
Unless a sample contains an impurity that provides mass loss at low temperature (< 200°C, e.g. gypsum, zeolites, halloysite, and soil organic matter; see the Impurities section below), the quantity of H2O removed from a clay during thermal dehydration corresponds to the smectite content in a sample (both as a discrete mineral and as a fraction of smectite-interstratified phases), and also depends on the hydration enthalpy of the interlayer cation, and – to a minor extent – the smectite layer charge (Środoń & McCarty, 2008; Środoń et al., 2009; Webb et al., 1986). If used for a semi-quantitative estimation of smectite content, a sample needs to be pre-equilibrated at a selected relative humidity (RH) that produces a particular degree of interlayer H2O filling, the latter also being dependent on the interlayer cation type (Berend et al., 1995; Cases et al., 1997). Despite this ambiguity, the measured fraction of mass loss due to dehydration is good enough for a relative comparison among samples tested. For example, pure Ca2+-saturated smectite equilibrated at ~50% RH contains ~16.5 wt.% of H2O (as determined by drying at 200°C; Środoń & McCarty, 2008); the corresponding mass loss upon dehydration decreases proportionally to illite interstratification in illite-smectite (Środoń et al., 2009) and smectite mineral content in a sample (Środoń, 2009). When equilibrated at ~50% RH, dehydration would correspond to 6% mass loss in K+-saturated smectite, ~10% in Na+- and Li+-saturated smectite, and almost 20% in Mg2+-saturated smectite (Berend et al., 1995; Cases et al., 1997). The above-mentioned semi-quantitative measurement should not be carried out under any dry gas purge (common in TG instruments) as it triggers dehydration during the sample-weight stabilization under dry gas flow; a moisture analyser is thus preferred (Środoń et al., 2009).
Dehydroxylation
The actual thermogravimetric manifestation of the reaction is the loss of H2O after its migration through the interlayer and out of the crystallite. The theoretical mass loss upon dehydroxylation is calculated as the molecular weight of H2O per a smectite’s formula unit molecular weight. Dehydroxylation produces an evolved gas pattern of H2O exactly following the DTG curve and TG mass loss, corresponding quantitatively to Eq. 1 (Fig. 1). In Ca2+- or Na+-exchanged Fe-poor smectites mass loss corresponding to dehydroxylation is 4.7–4.9%, whereas in nontronite it can be as low as 4.3%. Exchanging a smectite with a heavy cation such as Cs+ can decrease the theoretical mass loss significantly. The smectite’s mass loss due exclusively to dehydroxylation may not be clearly distinguishable in TGA analysis, and thus, may not correspond to the theoretical value. In numerous dioctahedral smectites, especially those having an interlayer cation of high hydration enthalpy (Ca2+, Mg2+, etc.), the onset of dehydroxylation overlaps with the dehydration of residual, strongly bound water (Bray & Redfern, 2000; Kuligiewicz & Derkowski, 2017).
The first control on the temperature of dehydroxylation is the octahedral occupancy. In dioctahedral 2:1 layer phyllosilicates the proton from one OH group migrates either along the shared edge of the octahedron or jumps across the octahedral vacancy to join the adjacent OH group (Guggenheim et al., 1987; Ogloza & Malhotra, 1989; Muñoz-Santiburcio et al., 2016). In dioctahedral smectites, dehydroxylation transforms the octahedral structural arrangement into a five-fold coordinated bipyramidal structure and the previously octahedral cations remain in a similar arrangement to the pristine structure (for details see Drits et al., 1995). The 2:1 layer thus remains at least partially intact. In the case of trioctahedral phyllosilicates, due to the lack of an octahedral vacancy where the residual oxygen could reside, OH groups cannot be removed without a reconstruction of the octahedral sheet; their dehydroxylation is accompanied by recrystallization and occurs at higher temperatures than in dioctahedral species, therefore (Fig. 1). The high-temperature DTG peak (usually >700°C) can serve as a distinction between trioctahedral layers or domains and dioctahedral counterparts in natural smectites. Synthetic saponites, however, showed strong dependence of the dehydroxylation temperature on the octahedral sheet cation composition: Mg > Ni > Co > Zn (Vogels et al., 2005). Vogels et al. (2005) proposed that the smaller the octahedral cations, the greater the corrugation of the tetrahedral sheet, which results in greater thermal stability.
For dioctahedral smectites, the next factor controlling the temperature of dehydroxylation is octahedral symmetry. Drits et al. (1995, 1998) showed that the cis-vacant (cv) structure has a significantly higher dehydroxylation temperature (>600°C) than the trans-vacant (tv) layer structure (< 600°C; Fig. 3). Those authors assigned that difference to the longer distance between non-shared OH edges in the cv structure than the shared OH edge of the octahedron in the tv coordination, resulting in the less probable and longer H+ jump between the OH groups in the cv coordination. Besides, in contrast to one-step dehydroxylation of tv layers, dehydroxylation of aluminous cv smectite is a two-step process that involves the Al cation migration to cis-sites and results in the pseudo-tv structure of the dehydroxylate (Drits et al., 1995; Stackhouse et al., 2004), the same as that of the dehydroxylated pristine tv smectite. In tv Mg-Fe-rich smectites (e.g. nontronite), however, some cations migrate to trans-positions during dehydroxylation in order to redistribute the Fe–O bonds, resulting in a cv dehydroxylated structure (Muller et al., 2010a, b Drits & Zviagina, 2009).
Following Drits et al. (1998), Wolters and Emmerich (2007) and Emmerich et al. (2009) adopted 600°C as the dehydroxylation temperature boundary and suggested using it for the determination of cv and tv compounds in smectites.
Heterogeneity of the octahedral sheet’s cation composition contributes to the broad range of smectite dehydroxylation temperature; in smectite that has one cation species dominating the octahedral sheets the temperature range of dehydroxylation is narrower than those observed in smectites with heterogeneous octahedral composition (Drits et al., 2012a). In trioctahedral smectites the octahedral sheet is rarely fully occupied by cations. Locally occurring vacancies, e.g. in stevensite, produce dioctahedral domains which dehydroxylate at temperatures corresponding to the dioctahedral smectite (i.e. montmorillonite).
The interlayer cations’ hydration enthalpy controls the temperature of final dehydration and, thus, the degree of overlap for the dehydration and dehydroxylation reactions in some dioctahedral smectites (Kuligiewicz & Derkowski, 2017; Fig. 4). Also, the interlayer cation size and its potential for migration into the octahedral sheet play a role in the pathway of smectite dehydroxylation. Early DTA studies demonstrated that the interlayer cation influences the temperature of maximum dehydroxylation (El-Barawy et al., 1986; Mackenzie & Bishui, 1958). Later, Bray and Redfern (2000) and Kuligiewicz and Derkowski (2017) found differences in the kinetics of the dehydroxylation reaction in a smectite containing monovalent and divalent interlayer cations. Cs+ is the largest inorganic cation that can pillar the interlayer (Table 1), allowing H2O molecules generated during dehydroxylation to diffuse easily, resulting in a decrease in the temperature of dehydroxylation in dioctahedral smectites (Fig. 4). The more beidellitic the smectite structure is, the greater the effect of Cs+-pillaring (Fig. 4: beidellite > beidellitic montmorillonite and nontronite > montmorillonite) because of the repulsion forces between adjacent siloxane surfaces created by Al substitution in the tetrahedral sheet and also the potential for the Cs+ cation to be located above the Al-substituted tetrahedra instead of the ditrigonal cavity (Cuadros, 2002; Ngouana & Kalinichev, 2014). The cation size, however, seems to result in an opposite trend in trioctahedral smectites; the larger the interlayer cation, the higher the temperature of dehydroxylation (Fig. 4). The ability to block efficiently the ditrigonal cavity by a large cation, which suppresses H+ migration into the adjacent OH group, is likely the reason for this effect.
In montmorillonites, where layer charge is generated in the octahedral sheet, dehydration enables small cations to migrate into the ditrigonal cavity of the tetrahedral sheet and even further into the octahedral vacancy, forming pseudo-trioctahedral domains (Emmerich et al., 1999; Hofmann & Klemen, 1950; Komadel et al., 2005, and references therein; Skoubris et al., 2013). The Hofmann-Klemen effect affects the pattern of dehydroxylation, forming an additional DTG peak at temperatures corresponding to trioctahedral domains. Not only Li+, commonly used in reducing layer charge, but also other small cations can produce such thermal effects. Koster van Groos and Guggenheim (1989) explained differences observed between the dehydroxylation of Mg2+- and Ca2+-exchanged montmorillonite as resulting from the different abilities of Mg2+ and Ca2+ to migrate toward the residual oxygen during the dehydroxylation reaction. The ability to form pseudo-trioctahedral domains must then depend on the kinetics of dehydration and small-cation migration through the ditrigonal cavity. Indeed, as presented in Fig. 5, Li+ migration (Li+ –ΔHo = 519 kJ/mol; Table 1) forms the high-temperature DTG peak weakly dependent on the dehydration temperature prior to dehydroxylation, whereas Mg2+-montmorillonite (Mg2+ –ΔHo = 1921 kJ/mol; Table 1) produces a pseudo-trioctahedral DTG signature at ~900°C, which is proportional to the temperature and time of isothermal dehydration (Fig. 5). Fast dehydroxylation during ramp-heating itself may not be sufficient to quickly and completely dehydrate Mg2+ and fix the cation in the ditrigonal cavity or the octahedral sheet.
Dehydrogenation
Dehydrogenation proceeds during thermal treatment of Fe(II)-bearing phyllosilicates, resulting in the removal of one hydrogen atom while simultaneously oxidizing one Fe(II) cation that is octahedrally coordinated with the OH group (Farmer et al. 1971; MacKenzie & Berezowski, 1984; Rancourt et al., 2001; Sanz et al., 1983; Wicks & Ramik, 1990).
The presence of oxygen enhances dehydrogenation, shifting the reaction to lower temperatures. The MS curve of m/z = 18 of evolved gas in the reaction from Eq. 3 is excessive to the DTG curve, because H2O is not produced within the heated crystallite and thus does not produce the corresponding weight loss (Lempart et al., 2020).
Because of structural constraints, Fe(II) in natural smectites occurs in numerous trioctahedral species (Kohyama et al., 1973; Stucki, 1988); therefore, these phases are the most prone to dehydrogenation. In nature, dehydrogenation may be responsible for the depletion of OH groups in so-called oxy-smectites (Dainyak et al., 1981). Dehydrogenation following Eq. 2 or Eq. 3 is expected upon heating Fe-reduced nontronites (Heller-Kallai, 1997a; Lear & Stucki, 1985).
Smectite Structural Transformation upon Heating
Amorphization and recrystallization are the reactions that occur during thermal destabilization of the dehydroxylation and dehydrogenation products (Grim & Rowland, 1942a,b; Meyers & Speyer, 2003). At temperatures higher than dehydroxylation, dioctahedral smectites undergo amorphization and then recrystallize into a topotaxic set of mullite, spinels, quartz or cristobalite, and pyroxenes (Bradley & Grim, 1951; Brett et al., 1970). The degree and temperature of amorphization and transformation depends on the octahedral-sheet composition. Beidellite dehydroxylate is the most stable high-temperature smectite phase, the irreversible amorphization of which has not been observed below 900°C (Derkowski et al., 2012a; Derkowski & Kuligiewicz, 2017; Kawano & Tomita, 1991b), whereas recrystallization occurs above 1200°C (Bradley & Grim, 1951). Sequentially heated Mg-rich dioctahedral smectites show an onset of amorphization at the end of dehydroxylation, whereas maximum amorphization is reached by ~800°C (Garg & Skibsted, 2014; Derkowski & Kuligiewicz, 2017), which provides their pozzolanic properties, but irreversibly removes their adsorption potential. In smectites rich in redox-active metals, i.e. nontronite, the temperature of transformation and the composition of transformation products depend heavily on the type of ambient gas (MacKenzie & Rogers, 1977).
In trioctahedral smectites, transformation proceeds along with dehydroxylation; after initial amorphization, a pyroxene-like structure and an olivine-like structure are formed (Vogels et al., 2005).
In contrast to dehydroxylation and dehydration, phase transformation is manifested in the DTA or DSC curves and not in TGA and TGA-MS, because it is not associated with a mass change or evolution of volatiles. However, note that at high temperature (>800°C), especially for isothermal experiments, some alkali cations can be removed due to vaporization, which may be responsible for occasionally observed minor and slow weight loss over a broad temperature range.
Rehydration and Rehydroxylation
Rehydration is the adsorption of H2O molecules by a dried material. Rehydroxylation is the reconstruction of the structural OH groups, i.e. reversal of the reaction in Eq. 1, in a dehydroxylated smectite structure that has not undergone phase transition (amorphization or recrystallization). Whereas rehydroxylation in general must be preceded by H2O diffusion through the dehydroxylated structure, advanced rehydroxylation requires specifically the diffusion within the interlayer in order to access all potential reaction sites available in the 2:1 layer. Both rehydration and rehydroxylation can be observed in situ as mass gain once wet gas is delivered directly to the TGA chamber (Derkowski & Kuligiewicz, 2017). When exposed to air, dried smectite having divalent cations in the interlayer starts resorbing H2O from ambient moisture within seconds (Fig. 6; Kanik et al., 2022). For this reason, in practice, smectite material for rehydration and rehydroxylation tests also needs to be pre-dehydrated or pre-dehydroxylated in situ. The rate and extent of rehydration in pristine, dried (i.e. dried below the temperature of dehydroxylation) smectite depends heavily on the pre-drying conditions and, thus, the extent of interlayer cation fixation, and the temperature and RH during rehydration (Fig. 6; Derkowski et al., 2012b, Gailhanou et al., 2017). Structurally, the extent of rehydration is controlled by the bonding between the siloxane surface and the interlayer cation; the greater the layer charge and the lower the cation hydration enthalpy, the stronger the bonding. K+-exchanged high-layer charge smectites would not rehydrate extensively keeping their interlayer closed, whereas Ca2+-exchanged smectite is prone to fast rehydration (Derkowski et al., 2012b; Kawano & Tomita, 1991a, b). Once a small interlayer cation, like Li+ is dehydrated and migrates through the ditrigonal cavity in montmorillonite (Hofmann-Klemen effect; Fig. 5), however, the interlayer is fixed permanently, prohibiting rehydration (Derkowski & Kuligiewicz, 2017; Komadel et al., 2002). A similar effect occurs even if the interlayer cations do not migrate into the 2:1 layer: if layer charge is located in the octahedral sheet, dehydrated interlayer cations are bonded strongly to the 2:1 layer, resulting in rehydration slower than if the charge came from tetrahedral cation substitution (Kawano & Tomita, 1991a,b).
Rehydroxylation occurs mostly in the intact dehydroxylated phase that preserves the layer structure and can undergo the reverse of the Eq. 1 reaction (Derkowski et al., 2012a). Once smectite crystallites reach the high-temperature amorphization or recrystallization, the structural integrity is lost. Therefore, dehydroxylated trioctahedral smectite, Mg-rich montmorillonite, and nontronite do not have great potential for rehydroxylation, whereas beidellite is the most prone to rehydroxylation because its layer structure remains intact at high temperatures (Derkowski et al., 2012a).
Because dehydroxylation of originally tv and cv aluminum smectites produces a similar tv phase, when rehydroxylated it resembles a tv structure (Drits et al., 1995; Muller et al., 2000a, b). Indeed, when analyzed with thermal methods, such a rehydroxylated phase dehydroxylates at the temperature range typical for tv beidellite (<500°C; irrespective of the temperature of dehydroxylation of its pristine structure (Derkowski & Kuligiewicz, 2017; Fig. 7). A tv Fe,Mg-rich smectite that adopts a cv character after dehydroxylation can probably return to tv when rehydroxylated (Drits & Zviagina, 2009; Muller et al., 2000b). The most intense rehydroxylation occurs under H2O vapor at temperatures of 250–350°C, even at extremely low H2O activity. At such conditions, up to tens of % of the dehydroxylated dioctahedral smectite structure can be rehydroxylated within an hour of reaction (Fig. 7; Derkowski & Kuligiewicz, 2017). The mass gain observed during isothermal treatment of the prior dehydroxylated phase has been proven to correspond to true rehydroxylation by XRD and infrared spectroscopy (Kawano & Tomita, 1989, 1991b; Muller et al., 2000b; Derkowski et al., 2012a; Derkowski & Kuligiewicz, 2017). At reaction temperatures providing intense rehydroxylation, the interlayer’s H2O diffusion does not result in permanent rehydration (Derkowski et al., 2012b).
Rehydroxylation is site-specific and is influenced by several factors. If the layer charge is located in the tetrahedral sheet, the opening of the interlayer space during rehydroxylation is relatively easy due to large electrostatic repulsion between adjacent layers. Interlayer cation type also affects the rehydroxylation behavior of a smectite (Kawano & Tomita, 1989, 1991b; Emmerich et al., 2001; Derkowski et al., 2012a; Derkowski & Kuligiewicz, 2017). Interlayer cations with relatively large radii (>0.13 nm in the eight-fold coordination; Shannon, 1976) act as pillars, which leads to a relatively large interlayer space of the dehydroxylate and enhances its rehydroxylation potential. Cations with large hydration potential also increase H2O migration in the interlayer region, which is the prerequisite for further H2O migration into the octahedral sheet and reformation of OH groups (Derkowski et al., 2012a,b). If an anhydrous cation has a small ionic radius (< 0.08 nm, Li+, Mg2+), it can migrate into a pseudohexagonal cavity or into a vacant position of the former octahedral sheet and physically block access to the residual oxygen, which results in a reduced ability to rehydroxylate. Such cation migration results in irreversible collapse of the interlayer space (Emmerich, 2000; Emmerich et al., 2001; Derkowski et al., 2012a). In practice, rehydroxylation in smectites is controlled by the fraction of AlAlOH environment in the octahedral sheet, with no Al substitution in the neighboring tetrahedral sites (“pyrophyllite-like” domains), resulting in the lack of charge-compensating interlayer cations blocking the ditrigonal cavity for H2O migration. Another factor is the valence of the exchangeable cation – the greater the valence, the fewer cations are present in the interlayer space and the lower the chance of blocking the pathway of potential rehydroxylation (Derkowski et al., 2012a).
Experimental Factors Affecting the Results of Thermal Analysis of Smectites
Carrier Gas Composition and Pressure
Smectite thermal analysis is usually performed under dry gas purge in order to remove the products of dehydroxylation or dehydration. The presence of water vapor in the ambient gas changes the reaction pattern, shifting the reaction to higher temperatures following the H2O content in the purge gas (Emmerich et al. 2017; Koster van Groos & Guggenheim, 1987, 1989; Levy & Hurst, 1993). Smectite dehydration and dehydroxylation reactions under high H2O activity and high pressure were the subject of dedicated experiments enabling the prediction of smectite behavior under conditions relevant to geologic systems (Koster van Groos & Guggenheim, 1990, and references therein). Vacuum conditions decrease the smectite reaction temperatures with respect to ≥1 bar conditions (Wicks & Ramik, 1990).
An analysis of redox-sensitive minerals requires use of a virtually oxygen-free purge gas (Bottom, 2008, Bish & Duffy, 1990). The presence of oxygen in the purge gas can be tested in such cases using CuC2O4 (Lempart et al., 2018). The purge gas is considered virtually oxygen-free if no mass gain is detected during the experiment with pure CuC2O4.
Sample Form and Packing
With some exceptions, smectite samples are typically analyzed in a powder form. Before a thermal experiment, powdered samples should be homogenized and preferably sifted through a 100 or 200 μm mesh sieve, which allows for the reduction of the size of smectite aggregates present in the sample. In the case of montmorillonite that has a temperature of amorphization close to that of dehydroxylation, the presence of large aggregates leads to the formation of macroscopic bubbles and hollow spheres, especially for samples with exchangeable cations of high hydration enthalpy, such as Ca2+, Mg2+, or Li+ (Fig. 8; Emmerich, 2010; Emmerich et al., 2017). The spheres are probably produced by slow heat transfer inward and slow diffusion of H2O outward from the aggregates. Interestingly, the formation of such spheres has virtually no effect on the DTG curve or mass loss during a TG experiment (Fig. 8).
The way powder is loaded into a holder can also affect the reaction temperature (Emmerich, 2010); a pressed sample delays the diffusion of the reaction volatile products out of the holder due to decreased permeability. Although the effect is not greater than a few °C, and thus is within the temperature error (Schilling, 1990), use of a loose powder or exactly the same mode of pressing for all samples being compared is recommended. Also, the experiments involving surface reaction with ambient gas (i.e. dehydrogenation under oxidizing conditions) may be strongly affected by the sample permeability to gases.
The maximum reaction temperature correlates positively with the sample mass due to the heat capacity and the diffusion rate of volatile reaction products (Schilling, 1990; Wolters & Emmerich, 2007). However, as found by Derkowski and Kuligiewicz (unpublished data), up to 20% relative of smectite sample weight difference has a negligible effect on the reaction temperature and remains within the analytical error; 20% would thus be the maximum difference in the sample weights recommended for relative comparison. The greater the sample weight, the lower the relative influence of total measurement errors (weighing imprecision, thermal drift) that are usually expressed in absolute quantities, i.e. are independent of the sample size. However, upon heating, a smectite sample can change its bulk volume (due to recrystallization, amorphization) or electrostatic properties (due to dehydration, cation fixation), resulting in the holder overflow and the drop of a sample portion off the holder, which is detrimental to the analysis. Fortunately, a material drop is usually abrupt enough to be distinguishable from a slower thermal reaction of smectites, allowing the results to be rejected from interpretation. Such a sample drop is typical for Li+-montmorillonites at the temperature of dehydroxylation, when a sample is blown off the holder by purge gas (Fig. 9). Using a holder with a lid may prevent that problem (e.g. Emmerich et al., 2017; Wolters & Emmerich, 2007), but this would also induce a delay in removal of volatiles from the sample. The delayed removal of the reaction products by the gas purge results in momentary production of a water vapor atmosphere around the sample.
Factors Affecting Evolved Gas Analysis (EGA)
Isothermal reactions are usually not recorded in EGA. The EGA requires a high heating rate and high carrier gas flow, unless a vacuum system is used (Wicks & Ramik, 1990 cf. Lempart et al., 2020). High heating rates (and possibly large sample masses) are usually needed to keep the evolved gas in a concentration range that is easily determined by the applied spectrometric method.
In the MS signal, H2O is represented by four spectral masses, the dominant m/z = 18, and the minor m/z = 17 (OH+), whereas, the H3O+ (m/z = 19) and O+ (m/z = 16) can be detected in orders of magnitude lower intensity (Meija et al., 2006; Emmerich, 2010). With pure dehydration and dehydroxylation reactions, the m/z 17 follows precisely m/z 18, and both signals overlap the DTG curve (Fig. 1). The NH3 evolution from NH4+-smectites would produce excess m/z 17 whereas dehydrogenation under O2 (or air) would produce an additional H2O signal distinguishable as an excess to the DTG curve (Heller-Kallai, 1997b, Fig. 5 in Lempart et al., 2020; Fig. 1). Dehydrogenation that produces H2 (m/z = 2) is difficult to detect and requires a fast sample-heating rate (≥20°C/min; Lempart et al., 2020) in MS or employing a GC–MS system.
Halogen elements present in trioctahedral smectites substituting OH groups (mostly F, also Cl, B) evolve from a smectite sample as hydrogen halides or radicals forming gaseous diatomic molecules (Heller-Kallai, 1997b; Wicks & Ramik, 1990). Strong reactivity of H-halide molecules with metals makes them react with any metal-based elements of the analyzer setup, however, and so ceramic pans and skimmer TGA-MS interfaces are suggested.
Impurities
The mass loss recorded with TGA in a pure smectite should correspond to the cation-dependent dehydration and the theoretical mass loss corresponding to dehydroxylation (Emmerich, 2010; Kuligiewicz & Derkowski, 2017). In nature, smectite never occurs as a pure mineral, even in bentonite beds. Also, even a thorough chemical treatment aimed at removing impurities and grain-size fractionation (Jackson, 1969), usually leaves some impurity in the studied smectite fraction. Some impurities are not reactive in the temperature range corresponding to major smectite reactions, e.g. the presence of quartz and feldspar only reduces the theoretical mass loss due to dehydroxylation. Other impurities, however, have their own thermal reactions that can strongly affect the thermal analysis pattern of smectite samples (Emmerich, 2010). The presence of discrete or interstratified kaolinite frequently coexisting with aluminous smectites (such as beidellites), increases the sample’s mass loss in the dehydroxylation range, as kaolinite has three times more OH groups per formula unit to form H2O than smectite. The above-mentioned impurities can be detected by XRD and for kaolinite also by infrared spectroscopy.
During smectite cation exchange with divalent cations, atmospheric CO2 dissolved in the exchange solution binds with the cations forming carbonate impurities. In TGA tests they provide a significant mass loss corresponding to 44–52% of the carbonate content depending on the type of carbonate; the same mass loss would come from naturally occurring carbonate minerals present in a clayey sample (Emmerich & Smykatz-Kloss, 2002; Emmerich, 2010), if they were not carefully removed during the chemical pretreatment (Jackson, 1969). Therefore, even minute quantities of carbonate minerals, untraceable using other methods (e.g. XRD), can affect the smectite mass-loss determination. Thermal decarbonatization resulting in CO2 evolved usually at temperatures of 700–900°C can be observed with EGA, however, e.g. by monitoring m/z = 44 in the TGA-MS setup (Fig. 1), in order to apply appropriate mass-loss corrections.
Sedimentary or soil organic matter (SOM) can be bound strongly to smectite surfaces and is difficult to remove even by extensive chemical treatment (Mikutta et al., 2005; Szymański et al., 2022). Thermal decomposition of recalcitrant SOM under neutral gas results in various evolved molecules, usually below 500°C, identifiable by EGA-MS as a series of various signals with m/z > 50 (Emmerich, 2010), along with those belonging to hydrocarbons (methane with m/z = 16), and CO2 (m/z = 44; Fig. 1). The SOM-evolved H2O (m/z = 18) and NH3 (m/z = 17) overlap with the smectite-derived H2O signal. Under neutral gas, only a portion of SOM decomposes, leaving charcoal as a high-temperature product. If run under oxidizing gas (i.e. air), SOM oxidizes completely, usually below 400°C, resulting in evolved H2O and CO2 identifiable by EGA, and the excess mass loss equal the SOM content. The difference in mass loss and evolved gas composition (especially CO2) between the experiments under neutral and oxidizing gas can reveal the SOM presence. However, structural Fe(II) present in the smectite or associated phases can decrease mass loss under oxidizing conditions, canceling the SOM effect, leaving CO2 (m/z = 44) as the most reliable proxy.
Hydrous sulfate minerals (e.g. gypsum) are the first to be removed during smectite sample pretreatment. However, if present in a sample even in minute amounts, their mass loss upon dehydration at 100–200°C (for example, 21% in gypsum) increases smectite dehydration. Desulfurization occurring at >700°C – in the range of trioctahedral smectite dehydroxylation – results in a loss of almost 50% of gypsum mass, which can be detected easily using EGA-MS (m/z: 64, 48, 32; Emmerich, 2010; Dietel et al., 2017). If analyzed under oxidizing conditions, trace sulfides, e.g. pyrite, reacting with O2 also release SO2 albeit at lower temperatures (<500°C; Dietel et al., 2017). When reacting with H2O coming from a smectite’s dehydroxylation at 500–800°C, sulfides produce H2S molecules (Levy & White, 1988).
Kinetic Calculations of Smectite Thermal Reactions
Theoretical Background
As in any other reaction, de- and re-hydration and de- and re-hydroxylation can be tracked using common kinetic approaches. The comprehensive review of methodology applicable in solid-state kinetics studies of thermally stimulated processes can be found elsewhere (e.g. Khawam & Flanagan, 2006a,b; Vyazovkin et al., 2011, 2020). A brief introduction of the key methodological points is presented below, with special emphasis on methods applicable to the thermal reactions of smectite.
Finding the triplet f(α), Ea, and A allows for the full kinetic description of the reaction studied. However, many kinetic studies present only an estimation of Ea as the most relevant reaction-specific parameter. Depending on the mode in which the thermal analysis experiment was performed, Ea can be calculated using isothermal or non-isothermal data, the latter most commonly coming from linear heating experiments (Fig. 10). In any case, using at least three temperature programs is strongly advised (Vyazovkin et al., 2011; Figs. 10, 11b), and at least five points are required to detect any non-Arrhenius behavior of the reaction, which manifests itself by non-linearity of Eqs. 6 and 8 (Vyazovkin et al., 2020).
Two general types of kinetic calculations are used for thermal analysis data. The model-fitting approach is based on fitting one of the reaction models available in the literature to Eq. 5 or its integral form and finding the best-fitting model, simply by comparing its linearity with experimental data. The reaction mechanism may or may not include adjustable parameters, which can later be interpreted in terms of a rate-controlling model (Bray & Redfern, 2000; Giese, 1990; Vyazovkin et al., 2011). In smectites, the isothermal model-fitting approach has been proved useful for determining the kinetics of dehydration (Bray & Redfern, 1999), dehydroxylation (Bray & Redfern, 2000), and rehydroxylation (Derkowski & Kuligiewicz, 2017).
The second type of calculations are those using the isoconversional methods (Criado et al., 2008, Perez-Rodriguez et al., 2010; Vyazovkin, 2015, Vyazovkin, 2021), in which Ea and A are calculated without any assumptions of a particular kinetic model or the reaction order. The basic principle behind the isoconversional methods is that the value of f(α) in Eq. 5 is independent of T. In other words, the value of f(α) is constant for a given α, irrespective of the temperature program. In its simplest form for non-isothermal conditions (the Friedman method), calculating isoconversional Ea involves plotting the values of ln(dα/dt) against 1/T for a given α. The isoconversional approach enables changes in Ea to be detected during the reaction progress with no prior assumptions regarding the reaction mechanism. Such changes can, for example, reflect the changes in the bond strength with the reaction progress (Kuligiewicz & Derkowski, 2017).
Note that some studies apply kinetics analysis to data obtained from a single temperature program only (i.e. only one isothermal heating temperature or only one linear heating rate), which is especially common in the older literature. While this approach was justifiable because of hardware and computational limitations in those times, it is now generally considered unreliable (Vyazovkin et al., 2011).
Kinetics of Smectite Thermal Reactions
Dehydration
Despite the large number of studies of smectite–water interactions, kinetics studies of smectite dehydration employing thermal analysis methods do not present a consistent picture. Besides the dependence of dehydration energy on the hydration enthalpy of the interlayer cations (Webb et al., 1986), it has been demonstrated that a change occurs in the reaction mechanism along the course of dehydration (Koster van Groos & Guggenhem, 1990; Bray & Redfern, 1999; Kuligiewicz & Derkowski, 2017; Poinsignon et al., 1982; Zabat & Van Damme, 2000). As already mentioned, physisorbed water in the micro- and mesopores and bulk interlayer water represent the ‘weakly-bound water’ (WBW). WBW is removed preferentially during dehydration because of its weaker interactions with the smectite surface. Change in the reaction mechanism during dehydration and the associated change in Ea (Table 2) correspond to the completion of weakly bound water removal and the onset of strongly bound water removal. The Ea of dehydration of WBW is in the range 35–79 kJ/mol for smectites having monovalent and divalent interlayer cations commonly found in nature (Ca2+, Mg2+, K+, Na+, Li+; Table 2). As deduced from Table 2, Ea of WBW dehydration in general does not depend on the exchangeable cation type or chemistry of the 2:1 layer. Once weakly bound water is removed, removal of the sequential portions of strongly bound water requires increasingly higher energy, as is manifested by the continuous increase in Ea with reaction progress, up to 190 kJ/mol (Kuligiewicz & Derkowski, 2017). This was interpreted as the result of bond reconfiguration upon removal of water molecules from hydration spheres of interlayer cations, in agreement with theoretical calculations (Pavlov et al., 1998).
Dehydroxylation
In the majority of dioctahedral smectites studied, Ea of dehydroxylation does not vary much with the reaction progress (Fig. 11A,C,D). If the initial (α < 0.2) and final (α > 0.8) portions of the data, corresponding also to lower R2 coefficients because of an overlap with other thermal effects are excluded, the Ea values remain similar over a broad range of α for beidellite SBId-1, cv beidellitic montmorillonite SWy-1, and ferruginous smectite SWa-1 (an approximation for nontronite) (Fig. 11). In Mg-rich montmorillonites that have charge generated exclusively in the octahedral sheet, however, Ea increases significantly with α, especially in the samples saturated with divalent cations (Fig. 11E,F). The increase of Ea in montmorillonite remains in line with the interpretation of Koster van Groos and Guggenheim (1989) who suggested the increase in the residual OH bond strength upon sequential removal of OH groups in heterogeneous structure.
The kinetics of smectite dehydroxylation is controlled by the same structural factors as the temperature of dehydroxylation (Figs. 1 and 3). Activation energy of dehydroxylation as computed from linear heating experiments with isoconversional methods correlates linearly with the temperature of maximum dehydroxylation, based on the DTG peak (Fig. 12, Table 3). The observed relationship implies that the rate-controlling mechanism of dehydroxylation is the same in all smectites regardless of their crystal-chemical composition.
Another factor influencing the kinetics of smectite dehydroxylation is the interlayer cation type. The application of the model-fitting approach to isothermal data allowed demonstration that the interlayer cation type can modify the dominant reaction model, which provides the best description of the process. Bray and Redfern (2000) suggested that dehydroxylation in smectite saturated with divalent cations (Ca2+ or Mg2+) is diffusion-controlled, while dehydroxylation in smectites saturated with monovalent cations (Na+ or K+) follows the first-order reaction kinetics. The influence of the interlayer cation on the dehydroxylation Ea value is less clear (Bray & Redfern, 2000; Kuligiewicz & Derkowski, 2017) but it seems that Mg2+-exchanged smectites have significantly lower Ea values than Ca2+- or Na+-exchanged counterparts (Bray & Redfern, 2000; Kuligiewicz & Derkowski, 2017). A speculative explanation may be that Mg2+, which is capable of retaining dissociated H2O (Mg-OH) up to temperatures of dehydroxylation, acts as a flux and lowers the Ea of dehydroxylation.
All of the above-mentioned results and interpretations come from studies using simple (isothermal or linear heating rate) thermal experiments. An alternative approach based on multiple heating–cooling cycles was employed by Drits et al. (2012a) to determine the kinetics of partial dehydroxylation in dioctahedral smectites and other 2:1 layer minerals. The protocol involves heating the samples with a constant heating rate and subsequent cooling by increasing each following cycle’s maximum temperature in small increments (i.e. 25°C). If a temperature increment between cycles is sufficiently small, the increase of α within one and the same cycle is negligible; therefore, α can be considered as constant for the cycle, meeting the requirement of the isoconversional kinetic approach. Within the cycle, the majority of the ln(dα/dt) over 1/T returns a linear fit, which Drits et al. (2012a) interpreted as the zeroth-order kinetics of partial dehydroxylation. Such a heating–cooling cycle approach apparently allows tracking of Ea in different stages of the reaction or portions of sample reacted at subsequently higher temperatures. However, for a portion of the sample reacting with constant f(α), the Arrhenius relationship represents the population that has already been treated thermally in the preceding cycle, hence, not corresponding to the actual rate-controlling mechanism in full reaction (Muñoz-Santiburcio et al., 2012; Ogloza & Malhotra, 1989). The difficulty in explaining the concave shape of Ea evolution over α (Fig. 13) and the maximum Ea being much lower than that for full reaction (Figs. 11 and 12) implies the limited use of the heating–cooling protocol and the need for specific interpretation of the kinetic data.
Rehydration and Rehydroxylation
The kinetics of smectite rehydration is methodologically challenging to study because the maximum amount of water which can be adsorbed by the smectite’s structure is dependent on temperature; for a fixed H2O molar content in the carrier gas, the lower the reaction temperature, the greater the rehydration (Fig. 6), but for constant RH the degree of rehydration is likely to be similar in spite of different temperatures (Vieillard et al., 2019). This renders non-isothermal experiments inapplicable. The rehydration process can be most conveniently tracked using the TGA experimental setup, allowing in situ dehydration and rehydration (Derkowski & Kuligiewicz, 2017). The mechanism of rehydration at temperatures <100°C cannot be described well by a commonly used solid-state reaction model (such as those presented in Table 1 of Khawam & Flanagan, 2006a). Models developed specifically for adsorption kinetics, such as a pseudo nth order model, usually provide much better description of the process (see Largitte & Pasquier, 2016, for a comprehensive review of models used in the adsorption kinetics). Rehydration kinetics are influenced by the exchangeable cation type; cations with low hydration energy (such as Cs+) adsorb quickly the small amount of water which they are capable of holding. For smectite saturated with cations with high hydration energy (such as Mg2+), the adsorption process continues over longer timescales (Fig. 6).
Isothermal rehydroxylation experiments under H2O vapor at temperatures above those typical for rehydration have shown that kinetics of smectite rehydroxylation can be described by the time-to-the-quarter (TTTQ) model expressed as relative mass gain equal to the temperature-dependent rate constant multiplied by the 4th root of reaction time (Hamilton & Hall, 2012; Derkowski & Kuligiewicz, 2017). Explained as corresponding to restricted, single-file diffusion in fired clay minerals, the TTTQ model has been applied to clay-based traditional ceramics and used to predict the time elapsed since firing, hence for archeological dating (Hare et al., 2016; Wilson et al., 2009). Both the root coefficient of the TTTQ (Barret, 2017) and the overall validity of TTTQ model for kinetic calculations have been criticized strongly (Kuligiewicz & Derkowski, 2021).
Studies of Ea of smectite rehydroxylation are still rare and the Ea values calculated based on the TTTQ model remain within an unexpectedly broad range, between 19 and 173 kJ/mol for experiments performed at 110–400°C (Derkowski & Kuligiewicz, 2017; Stevenson & Gurnick, 2016). However, Kuligiewicz and Derkowski (2021) recently proved that the TTTQ model describes well the progress of various reactions governed by unexpectedly different mechanisms. Ea obtained with the TTTQ model for those reactions was grossly incorrect despite an apparent agreement with simulated data, making the TTTQ model unreliable. Therefore, although the TTTQ model fits well the smectite isothermal rehydroxylation in some cases and has served for occasional successful dating of clay ceramics, its link to kinetics formalism and the mechanism of rehydroxylation are still unclear and a need exists to develop an alternative approach.
Acknowledgements
The authors thank Katja Emmerich, Georgios Chryssikos, Marek Szczerba, Małgorzata Lempart-Drozd, and three anonymous reviewers for their valuable comments. Nadine Kanik kindly provided corrections of the English language. The work was financed by the National Science Centre, Poland; grant OPUS-13 UMO-2017/25/B/ST10/01675.
Declarations
Conflict of Interest
The authors declare that they have no conflict of interest.