**Thermodynamics of Hydration in Minerals: How to Predict These Entities**

Philippe Vieillard

338 Thermodynamics – Fundamentals and Its Application in Science

Incl. Phenom. Macrocycl. Chem. 50: 3–11.

Chem. High T. Scr. 10: 350–357.

Cyclodextrins and Bile Salts. Supramol. Chem. 16: 371–379.

Recognition of Bile Salts. Bull. Korean Chem. Soc. 29: 2119–2124.

Molecular Recognition of Bile Salts. J. Solution Chem. 38: 417–428.

Bis(*β*-Cyclodextrin) for Bile Salts. Helv. Chim. Acta 93: 999–1011.

Bridged Bis(*β*-Cyclodextrin)s. J. Org. Chem. 69: 6590–6602.

Linkers. Supramol. Chem. 21: 409–415.

Cyclodextrin Dimers. J. Phys. Chem. B 109: 4129–4134.

1355–1362.

[25] Liu Y, Zhao Y.-L, Yang E.-C, Zhang H.-Y (2004) Molecular Recognition Thermodynamics of Steroids by Novel Oligo(aminoethylamino)-*β*-Cyclodextrins Bearing Anthryl: Enhanced Molecular Binding Ability by Co-inclusion Complexation. J.

[26] Li N, Chen Y, Zhang Y.-M, Li Z.-Q, Liu Y (2010) Binding Thermodynamics of Quinolinyl-Modified *β*-Cyclodextrins with Bile Salts. SCIENTIA SINICA Chimica 40:

[27] Liu Y, Li L, Zhang H.-Y, Yang Y.-W, Ding F (2004) Correlation between Thermodynamic Behavior and Structure in the Complexation of Modified β-

[28] Liu Y, Li L, Chen Y, Yu L, Fan Z, Ding F (2005) Molecular Recognition Thermodynamics of Bile Salts by *β*-Cyclodextrin Dimers: Factors Governing the Cooperative Binding of

[29] Liu Y, Kang S, Chen Y, Cao R, Shi J (2007) Thermodynamics of Molecular Recognition of Bile Salts by 3,6'-(Oligoethylenediamino-Bridged) *β*-Cyclodextrin Dimers. Comb.

[30] Zhao Y, Gu J, Chi S.-M, Yang Y.-C, Zhu H.-Y, Wang Y.-F, Liu J.-H, Huang R (2008) Aromatic Diamino-bridged Bis(*β*-Cyclodextrin) as Fluorescent Sensor for the Molecular

[31] Zhao Y, Gu J, Chi S.-M, Yang Y.-C, Zhu H.-Y, Huang R, Wang Y.-F (2009) Bis(*β*-Cyclodextrin)s Linked with an Aromatic Diamine as Fluorescent Sensor for the

[32] Zhao Y, Yang Y.-C, Chi S.-M, Shi H, Zhao Y, Zhu H.-Y, Li Q.-L, Wang Y.-F (2010) Fluorescence Sensing and Selective Binding of a Novel 4,4'-Sulfonyldianiline-Bridged

[33] Liu Y, Yang Y.-W, Yang E.-C, Guan X.-D (2004) Molecular Recognition Thermodynamics and Structural Elucidation of Interactions between Steroids and

[34] Liu Y, Wu H.-X, Chen Y, Chen G.-S (2009) Molecular Binding Behaviours of Bile Salts by Bridged and Metallobridged Bis(β-Cyclodextrin)s with Naphthalenecarboxyl

[35] Liu Y, Song Y, Wang H, Zhang H.-Y, Wada T, Inoue Y (2003) Selective Binding of Steroids by 2,2'-Biquinoline-4,4'-dicarboxamide-Bridged Bis(*β*-cyclodextrin):

[37] Liu Y, Yu H.-M, Chen Y, Zhao Y.-L (2006) Synthesis of Bridged and Metallobridged Bis(β-cyclodextrin)s Containing Fluorescent Oxamidobisbenzoyl Linkers and Their

[38] Ikeda H, Nakamura M, Ise N, Oguma N, Nakamura A, Ikeda T, Toda F, Ueno A (1996) Fluorescent Cyclodextrins for Molecule Sensing: Fluorescent Properties, NMR Characterization, and Inclusion Phenomena of *N*-Dansylleucine-Modified

Fluorescence Enhancement by Guest Inclusion. J. Org. Chem. 68: 3687–3690. [36] Liu Y, Song Y, Chen Y, Li X.-Q, Ding F, Zhong R.-Q (2004) Biquinolino-Modified β-Cyclodextrin Dimers and Their Metal Complexes as Efficient Fluorescent Sensors for

the Molecular Recognition of Steroids. Chem. Eur. J. 10: 3685–3696.

Selective Binding towards Bile Salts. Chem. Eur. J. 12: 3858–3868.

Cyclodextrins. J. Am. Chem. Soc. 118: 10980–10988.

Additional information is available at the end of the chapter

http://dx.doi.org/10.5772/51567

## **1. Introduction**

Many papers have suggested that several specific "types" of H2O exist. These have been labelled variously as " zeolitic", " loosely held", "structural", "crystal", "tightly bound" and "external" [1-3]. These labels suggest both the differing apparent energies of H2O as well as their differing apparent structural roles. Three useful distinctions can be made about H2O in compounds:


The first type (continuously varying H2O) is characteristic of the zeolites and clay minerals and considers that most of the zeolites and clay minerals lose or gain H2O in response to small changes in temperature and pressure over an extended temperature, relative humidity or the nature of cations in exchanged sites.

The second type of H2O is similar in nature to that found in hydrates of salts and considers that dehydration occurs over narrow temperature intervals in some compounds. Some zeolites, like analcime and laumontite, do not exchange water at room temperature. This type can be called "hydrate" H2O and has a specific position in the crystal structure. Hydroxides are excluded from this chapter because the H2O molecule is not identified but is only virtually present as a hydroxyl OH.

The third type of H2O is externally sorbed to the crystal and may be referred to as external. This type is present in quantities much smaller than the H2O present within the structure of any zeolite grain size. In the clay minerals, whose specific surfaces are greater than those of zeolites, the water located in inter-particle spaces, in inter-aggregate spaces and at the

surface of uncharged clays minerals indicates that a special status of H2O may be different in the regimes of adsorption or of capillary condensation, and this is not considered in this chapter.

### **2. Definition of thermodynamic entities of hydration**

In evaluating the thermodynamic quantities, it is important to define the standard state chosen for a given constituent. In the present study, the standard state for water is chosen as pure water at a temperature of 25°C with a vapour pressure of Po, and for the minerals, the dry compound is at a temperature of 25°C with a vapour pressure of P=0.

The reaction of hydration may be expressed by the following reaction:

$$X\_{\text{(c)}} + n^\* H\_2O\_{\text{(l)}} \to X - n H\_2O\_{\text{(c)}} \tag{1}$$

Thermodynamics of Hydration in Minerals: How to Predict These Entities 341

<sup>0</sup> <sup>Ξ</sup> ,298 2 of

,298 2 ,298 2 ,298 Ξ = Ξ\* + Ξ (5)

*<sup>f</sup>* ( ) () *HO n HO Hyd* ( ) *<sup>f</sup> liq hyd <sup>n</sup>* 0 00

the hydration water varies in all inorganic compounds, such as minerals, salts, zeolites and clays. There are two main questions: 1) Is the thermodynamic entity of the hydration water strongly related to the physical-chemical properties of the inorganic compound or does it remain constant within a structural family of inorganic compounds? and 2) What are the adequate relationships that can predict the thermodynamic entity of the hydration water in

Different methods of measurements of thermodynamic entities of hydration are presented


dehydrated zeolite [4-6] or at 25°C in a HF - HNO3 solution for clay minerals [7]. - Transposed-temperature drop calorimetry. The enthalpy of hydration data are obtained by measurements of the heat evolved from a sample dropped at room temperature into a calorimeter maintained at 700–800°C [8], [9-12]. One of the difficulties associated with this measurement includes the enthalpy effects because of the irreversible structural



of the method is that it is a direct measurement of the enthalpy of hydration.

hydration for zeolites [22], [23] and for montmorillonites [24-26].


The aim of this chapter is to determine how the thermodynamic entity *<sup>f</sup>* ( ) *H OHyd*

**3. Measurements of thermodynamics of hydration in minerals** 

any mineral?

in this section:

changes at 700–800°C.

clay minerals [7, 13].

If we consider any standard thermodynamic entity of hydration *hyd*,298 ° <sup>Ξ</sup> (enthalpy *Hhyd*,298 ° Δ free energy *Ghyd*,298 ° <sup>Δ</sup> , entropy *hyd <sup>S</sup>* ,298 ° Δ , heat capacity Cphyd. and molar volume VHyd.) for a given mineral having n bounded molecules, the entity of hydration is related to the total formation entity *<sup>f</sup>* ,298 ° Ξ (enthalpy *<sup>o</sup> <sup>f</sup>* ,298.15 ΔΗ free energy *<sup>o</sup> Gf* ,298.15 <sup>Δ</sup> , entropy *S*<sup>298</sup> ° , heat capacity *Cp*298.15 and molar volume V298) of the mineral by the following:

$$
\Xi^{o}\_{hyd,298} = \Xi^{o}\_{f,298} \left( X - nH\_2O \right) - \Xi^{0}\_{f,298} \left( X \right) - n^\* \Xi^{o}\_{f,298} \left( H\_2O\_{liq} \right) \tag{2}
$$

where ( ) *<sup>o</sup>* Ξ − *<sup>f</sup>* ,298 *X nH O*2 and ( ) *<sup>f</sup> <sup>X</sup>* <sup>0</sup> <sup>Ξ</sup> ,298 are the thermodynamic entity of formation of the mineral in the hydrated and dehydrated states, respectively, and *<sup>f</sup>* ( ) *H Oliq* <sup>0</sup> <sup>Ξ</sup> ,298 2 stands for the thermodynamic entity of formation of bulk water. The unit of the entity *<sup>o</sup>* <sup>Ξ</sup>*hyd*,298 is based per n moles of bound water or per mole of an anhydrous compound X.

The standard thermodynamic parameters of the hydration water can be calculated if the thermodynamic entities of formation (enthalpy *<sup>o</sup> <sup>f</sup>* ,298.15 ΔΗ free energy <sup>o</sup> <sup>Δ</sup><sup>G</sup> 15.298,f , entropy *S*<sup>298</sup> ° , heat capacity *Cp*298.15 . and molar volume V298.) of the anhydrous and hydrous end members are available in thermochemical tables or calculated from calorimetric measurements. However, the thermodynamic properties of the hydration water can be estimated by considering the hypothetical intra-crystalline reaction:

$$X\_{\text{(c)}} + n \, ^\ast H\_2O\_{\text{(hydr)}} \to \, X - nH\_2O\_{\text{(c)}} \tag{3}$$

where H2O(hyd) represents the hydration water. The thermodynamics of the intra-crystalline reaction can be written as follows:

$$
\Xi\_{r,298}^{o} = \Xi\_{f,298}^{o} \left( X - \mathfrak{n} H\_2O \right) - \Xi\_{f,298}^{0} \left( X \right) - \mathfrak{n}^\* \Xi\_{f,298}^{o} \left( H\_2O\_{hyd.} \right) \tag{4}
$$

By assuming that the standard thermodynamic entities of the intra-crystalline reaction are equal to zero, *i.e.*, *<sup>o</sup> <sup>r</sup>*,298 Ξ = 0 , the standard values of the molar thermodynamic entity of the hydration water, *<sup>f</sup>* ( ) *H OHyd* <sup>0</sup> <sup>Ξ</sup> ,298 2 , is finally obtained from eqs (2) and (4):

$$
\Xi\_{f,298}^{0}\left(H\_2O\_{hyd}\right) = \left(\mathfrak{n}^\* \, \Xi\_{f,298}^{0}\left(H\_2O\_{liq}\right) + \Xi\_{hyd,298}^{0}\right) \Big/ \mathfrak{n} \tag{5}
$$

The aim of this chapter is to determine how the thermodynamic entity *<sup>f</sup>* ( ) *H OHyd* <sup>0</sup> <sup>Ξ</sup> ,298 2 of the hydration water varies in all inorganic compounds, such as minerals, salts, zeolites and clays. There are two main questions: 1) Is the thermodynamic entity of the hydration water strongly related to the physical-chemical properties of the inorganic compound or does it remain constant within a structural family of inorganic compounds? and 2) What are the adequate relationships that can predict the thermodynamic entity of the hydration water in any mineral?

### **3. Measurements of thermodynamics of hydration in minerals**

340 Thermodynamics – Fundamentals and Its Application in Science

**2. Definition of thermodynamic entities of hydration** 

dry compound is at a temperature of 25°C with a vapour pressure of P=0.

The reaction of hydration may be expressed by the following reaction:

If we consider any standard thermodynamic entity of hydration *hyd*,298

capacity *Cp*298.15 and molar volume V298) of the mineral by the following:

per n moles of bound water or per mole of an anhydrous compound X.

estimated by considering the hypothetical intra-crystalline reaction:

° <sup>Δ</sup> , entropy *hyd <sup>S</sup>* ,298

thermodynamic entities of formation (enthalpy *<sup>o</sup>*

reaction can be written as follows:

hydration water, *<sup>f</sup>* ( ) *H OHyd*

equal to zero, *i.e.*, *<sup>o</sup>*

° Ξ (enthalpy *<sup>o</sup>*

chapter.

free energy *Ghyd*,298

*S*<sup>298</sup>

formation entity *<sup>f</sup>* ,298

surface of uncharged clays minerals indicates that a special status of H2O may be different in the regimes of adsorption or of capillary condensation, and this is not considered in this

In evaluating the thermodynamic quantities, it is important to define the standard state chosen for a given constituent. In the present study, the standard state for water is chosen as pure water at a temperature of 25°C with a vapour pressure of Po, and for the minerals, the

given mineral having n bounded molecules, the entity of hydration is related to the total

*hyd f X nH O X n H O <sup>f</sup> <sup>f</sup> liq* 0

where ( ) *<sup>o</sup>* Ξ − *<sup>f</sup>* ,298 *X nH O*2 and ( ) *<sup>f</sup> <sup>X</sup>* <sup>0</sup> <sup>Ξ</sup> ,298 are the thermodynamic entity of formation of the

the thermodynamic entity of formation of bulk water. The unit of the entity *<sup>o</sup>* <sup>Ξ</sup>*hyd*,298 is based

The standard thermodynamic parameters of the hydration water can be calculated if the

° , heat capacity *Cp*298.15 . and molar volume V298.) of the anhydrous and hydrous end members are available in thermochemical tables or calculated from calorimetric measurements. However, the thermodynamic properties of the hydration water can be

where H2O(hyd) represents the hydration water. The thermodynamics of the intra-crystalline

*r f X nH O X n H O <sup>f</sup> <sup>f</sup> hyd* 0

By assuming that the standard thermodynamic entities of the intra-crystalline reaction are

<sup>0</sup> <sup>Ξ</sup> ,298 2 , is finally obtained from eqs (2) and (4):

( ) () *o o <sup>o</sup>*

( ) () *o o <sup>o</sup>*

mineral in the hydrated and dehydrated states, respectively, and *<sup>f</sup>* ( ) *H Oliq*

*X n H O X nH O* ( ) 2 () *cl c* 2 () + →− \* (1)

° Δ , heat capacity Cphyd. and molar volume VHyd.) for a

*<sup>f</sup>* ,298.15 ΔΗ free energy *<sup>o</sup> Gf* ,298.15 <sup>Δ</sup> , entropy *S*<sup>298</sup>

*X n H O X nH O* () 2 ( ) *<sup>c</sup> hydr* 2 () *<sup>c</sup>* + →− \* (3)

,298 ,298 2 ,298 ,298 2 . Ξ =Ξ − −Ξ − Ξ\*( ) (4)

*<sup>r</sup>*,298 Ξ = 0 , the standard values of the molar thermodynamic entity of the

,298 ,298 2 ,298 ,298 2 Ξ =Ξ − −Ξ − Ξ\*() (2)

° <sup>Ξ</sup> (enthalpy *Hhyd*,298

<sup>0</sup> <sup>Ξ</sup> ,298 2 stands for

*<sup>f</sup>* ,298.15 ΔΗ free energy <sup>o</sup> <sup>Δ</sup><sup>G</sup> 15.298,f , entropy

° Δ

° , heat

Different methods of measurements of thermodynamic entities of hydration are presented in this section:



Thermodynamics of Hydration in Minerals: How to Predict These Entities 343

, Cp and V) of hydration water obtained in

Families Ref. Nb. ΔH°f ΔG°f S° Cp V comp (kJ/mole) (kJ/mole) (J/K/mole) (J/K/mole) cm3/mole

sulphites 1 41 -301.4±7.7 -244.4±7.9 41.5±6.0 39±5.2 14.5±1.2

Chlorides 1 34 -308.3±8.6 -251.8±8.5 43.7±5.9 42±2.7 14.1±2.3

Hydrates 3 342 -298.6±20.5 -242.4 40.9 42.8 14.75 Silicates 4 ? 39.3±4.2 13.7±0.5

Silicates 9 7 43.5 14 Silicates 10 7 -306.99±12.8 44.1±2.0 14.6±0.09

References: 1 - [35]; 2 - [38]; 3 - [39]; 4 - [40]; 5 - [41]; 6 - [42]; 7 - [43]; 8 - [44]; 9 - [45]; 10 - [46]; 11 - [47]; 12 - [48]; 13 - [49].

All of these models of predictions are built on a simple additive scheme and require knowledge of the thermodynamic entity of the anhydrous compound, the number of water molecules and the thermodynamic entity of formation of "ice-like" water or of hydration water obtained by statistical regressions. For each thermodynamic entity (enthalpy, free energy, entropy, heat capacity and molar volume), the values are significantly different within the nature of the compounds or minerals. For each prediction model, the nature of family, number of compounds and statistical errors are provided to check the extent and the validity of the model. There appear to be few models of prediction of free energies of hydration, which is the most useful parameter in geochemical modelling. The best way to reach this target is to combine the enthalpy and entropy of hydration. Despite the great diversity of the enthalpy of formation of the hydration water and its great error within different families of compounds, some recent and accurate models of prediction of the enthalpy of hydration are important to introduce for compounds such as salts, zeolites and

To improve the accuracy of prediction, the next section develops several models of prediction based on different kinds of minerals displaying recent measurements of the thermodynamic entities of the hydration water closely related either to the environmental

conditions and or to the nature and chemical formulas of the anhydrous compounds.

Zeolites 2 15 59.1

Silicates 5 3 54.123

RE chlorides 7 8 -297 46.4

Borates 12 23 -290.42 -237.28 54.86

Silicates 11 6 -292.37±4.6 -239.91±5.7

U bearing minerals 13 14 -295.58 -241.1

**Table 2.** Average thermodynamic entities (ΔH°f, ΔG°f, S°

different prediction models.

clays minerals.

Hydrates 6 ? -297.06

Binary comp. of Na 8 42 -297±8.3

Sulphates &

Chlorates &


## **4. Inventory of different values from prediction methods of thermodynamic entities of hydration water**

The thermodynamic data of liquid water and ice are displayed in Table 1. The formation properties of the hydration water may vary depending on the nature of the physicalchemical properties of the anhydrous compound [35]. In any case, when the minerals exist as hydrated phases, the presence of water molecules stabilises them. Thermodynamically, this corresponds to ΔG0hyd < 0 kJ mol-1 or ΔG0f,(H2OHyd) < -237.14 kJ mol-1, the value for bulk water at 250C and 1 bar [36]. Additionally, the maximum entropy is for bulk water and implies that S0(H2OHyd) < 69.95 J K-1 mol-1 [36]. The minimum entropy may be assessed by considering the entropy of Ih ice, as refined by Tardy et al. [37] (44.77 J K-1 mol-1). To give a minimum value for ΔG0f(H2OHyd) is more complicated. Ice cannot be of use, for example, because it is not stable at room temperature, and on the contrary, ΔG0f(Ice Ih) > ΔG0f (bulk water), with –236.59 kJ mol-1. For the enthalpy term ΔH0f(H2OHyd), because ΔG0hyd < 0 kJ.mol-1 and S0(H2OHyd) < S0(bulk H2O) for hydrated minerals, then ΔH0hyd must be negative or ΔH0f(H2OHyd) < ΔH0f(bulk H2O) with -285.83 kJ mol-1 [36]. For the heat capacity of hydration water, very few measurements have been actually performed so far.


References: 1 – [36]; 2 – [37]; 3 – [35].

**Table 1.** Average thermodynamic entities (ΔH° f, ΔG° f, S° , Cp and V) of liquid water and ice.

Among the numerous methods of prediction of thermodynamics entities of compounds (salts, and minerals), only twelve are retained and exhibit different values of the thermodynamic entities of ice-like water obtained from the statistical results from the anhydrous and hydrated compounds (Table 2).





The thermodynamic data of liquid water and ice are displayed in Table 1. The formation properties of the hydration water may vary depending on the nature of the physicalchemical properties of the anhydrous compound [35]. In any case, when the minerals exist as hydrated phases, the presence of water molecules stabilises them. Thermodynamically, this corresponds to ΔG0hyd < 0 kJ mol-1 or ΔG0f,(H2OHyd) < -237.14 kJ mol-1, the value for bulk water at 250C and 1 bar [36]. Additionally, the maximum entropy is for bulk water and implies that S0(H2OHyd) < 69.95 J K-1 mol-1 [36]. The minimum entropy may be assessed by considering the entropy of Ih ice, as refined by Tardy et al. [37] (44.77 J K-1 mol-1). To give a minimum value for ΔG0f(H2OHyd) is more complicated. Ice cannot be of use, for example, because it is not stable at room temperature, and on the contrary, ΔG0f(Ice Ih) > ΔG0f (bulk water), with –236.59 kJ mol-1. For the enthalpy term ΔH0f(H2OHyd), because ΔG0hyd < 0 kJ.mol-1 and S0(H2OHyd) < S0(bulk H2O) for hydrated minerals, then ΔH0hyd must be negative or ΔH0f(H2OHyd) < ΔH0f(bulk H2O) with -285.83 kJ mol-1 [36]. For the heat capacity of hydration water, very few measurements have been actually

f S° Cp V

, Cp and V) of liquid water and ice.

(kJ mol-1) (kJ mol-1) (J K-1 mol-1) (J K-1 mol-1) cm3 mol-1

Liquid water 1 -285.84±0.04 -237.18±0.05 69.923±0.03 75.576 18.064 Ice 2,3 -292.75±0.1 -236.59±0.1 44.771±0.05 42.008 19.228

Among the numerous methods of prediction of thermodynamics entities of compounds (salts, and minerals), only twelve are retained and exhibit different values of the thermodynamic entities of ice-like water obtained from the statistical results from the

the capacity to measure partial molar enthalpies of hydration directly.

developed by [28, 29] and applied on montmorillonite in [30] and [31]

raman spectroscopy has been tested on zeolites [32] and clays [33] [34].

**4. Inventory of different values from prediction methods of** 

**thermodynamic entities of hydration water** 

Families Ref. ΔH°f ΔG°

**Table 1.** Average thermodynamic entities (ΔH°f, ΔG°f, S°

anhydrous and hydrated compounds (Table 2).

performed so far.

References: 1 – [36]; 2 – [37]; 3 – [35].

References: 1 - [35]; 2 - [38]; 3 - [39]; 4 - [40]; 5 - [41]; 6 - [42]; 7 - [43]; 8 - [44]; 9 - [45]; 10 - [46]; 11 - [47]; 12 - [48]; 13 - [49].

**Table 2.** Average thermodynamic entities (ΔH°f, ΔG°f, S° , Cp and V) of hydration water obtained in different prediction models.

All of these models of predictions are built on a simple additive scheme and require knowledge of the thermodynamic entity of the anhydrous compound, the number of water molecules and the thermodynamic entity of formation of "ice-like" water or of hydration water obtained by statistical regressions. For each thermodynamic entity (enthalpy, free energy, entropy, heat capacity and molar volume), the values are significantly different within the nature of the compounds or minerals. For each prediction model, the nature of family, number of compounds and statistical errors are provided to check the extent and the validity of the model. There appear to be few models of prediction of free energies of hydration, which is the most useful parameter in geochemical modelling. The best way to reach this target is to combine the enthalpy and entropy of hydration. Despite the great diversity of the enthalpy of formation of the hydration water and its great error within different families of compounds, some recent and accurate models of prediction of the enthalpy of hydration are important to introduce for compounds such as salts, zeolites and clays minerals.

To improve the accuracy of prediction, the next section develops several models of prediction based on different kinds of minerals displaying recent measurements of the thermodynamic entities of the hydration water closely related either to the environmental conditions and or to the nature and chemical formulas of the anhydrous compounds.

#### **5. Recommended models of prediction of thermodynamic entities**

#### **5.1. Enthalpy of formation**

**Salts:** A model of the prediction of hydration enthalpy has been developed and tested on 349 different hydrates salts [50-52]. By considering the following reaction defining the hydration of a salt *A Bx x* 1 2 :

$$A\_{x\_1}B\_{x\_2} + n^\* \, H\_2O\_{\text{(c)}} \to \, A\_{x\_1}B\_{x\_2} - n \, H\_2O \tag{6}$$

Thermodynamics of Hydration in Minerals: How to Predict These Entities 345

are given in table 3 with the

(kJ/mol H2O)

Δ ). For a given anhydrous salt

° Δ H2O= -309.25 kJ mol-1) to values close to

° Δ = -297.22 kJ mol-1) and are different from the average

+

*<sup>H</sup> B A* 2 1 ( ) − +

Δ ), the enthalpy of hydration per mole of water appears to not be

<sup>+</sup> β Nb. of data % error1 Std dev.

<sup>+</sup> , β, number of data, % error and statistical error per moles of crystalline water.

The values of *<sup>z</sup> <sup>A</sup>*

hydrated salts.

*z A* α

α

(same value of *z z*

*<sup>H</sup> B A* 2 1 ( ) − +

decreases from large values (NaOH. H2O; *H <sup>f</sup>* ,298

those of ice (Na2Se. 16 H2O; *H <sup>f</sup>* ,298

α

value for 61 Na salts ( *H <sup>f</sup>* ,298

Cations *<sup>z</sup> <sup>A</sup>*

hydrated salts.

**Table 3.** Values of *<sup>z</sup> <sup>A</sup>*

α

<sup>+</sup> and β for the different salts of cation *<sup>z</sup> A* <sup>1</sup>

molecules and the nature of the crystalline salts ( *z z*

° Δ = -297.74 kJ mol-1).

respective standard deviations and % errors of the predicted enthalpy of formation of the

An example is given by considering the data of all of the sodium salts displayed in Figure 1A and shows the great variation of the enthalpy of hydration having the same number as the crystalline water. The different straight lines are obtained from Eq. (10) with values of

<sup>+</sup> and β for Na+ and are plotted for the different values of the number of crystalline water

constant but decreases when the number of molecules of the crystalline water increases. Figure 1B displays the enthalpy of crystalline water versus the nature of the anhydrous salt showing the great disparity and shows that the enthalpy of the crystalline water in all of the sodic salts is more negative than those of the hexagonal ice. When the number of crystalline water molecules of any salt increases, the enthalpy of formation of the crystalline water

Li+ -0.277 -0.341 19 0.34 ± 3.58 Na+ -0.294 1.759 61 0.29 ± 3.44 K+ -0.115 -9.168 15 0.66 ± 10.70 Rb+ -0.360 -6.769 6 0.29 ± 3.38 Cs+ -0.577 18.954 4 0.48 ± 6.35 Be+2 -0.485 -3.068 7 0.25 ± 3.48 Mg+2 -0.312 -10.708 21 0.43 ± 5.39 Ca+2 -0.270 -3.892 25 0.44 ± 4.64 Sr+2 -0.349 -1.894 20 0.32 ± 5.93 Ba+2 -0.327 -4.637 19 0.80 ± 5.31 Zn2+ -0.211 -16.461 12 0.44 ± 6.57 Cd2+ -0.322 2.266 8 0.33 ± 4.29 Cu+2 -0.433 -1.791 11 0.65 ± 4.22 Ni2+ -0.176 -18.056 10 0.46 ± 2.81 Co2+ -0.260 -2.296 11 0.36 ± 3.60 Fe2+ -0.136 -17.806 6 0.18 ± 1.58 Mn2+ -0.238 -5.144 15 0.25 ± 6.14 Cr+2 -0.478 31.459 3 0.28 ± 1.77 UO2+2 -0.247 -7.945 18 0.35 ± 4.80 1: the % error weighs the difference (in %) between the measured and predicted enthalpy of the formation of the

the enthalpy of hydration of the reaction (6) becomes

$$
\Delta H^{o}\_{hyd,298} = \Delta H^{o}\_{f,298} \left( A\_{x\_1} B\_{x\_2} - n H\_2 O \right) - \Delta H^{o}\_{f,298} \left( A\_{x\_1} B\_{x\_2} \right) - n^\* \Delta H^{o}\_{f,298} (H\_2 O\_{\{c\}}) \tag{7}
$$

where n is the number of molecules of crystalline water contained in the hydrate under consideration and *H <sup>f</sup>* ,298 ° Δ (H2O)(c) represents the enthalpy of formation of ice. A parameter *z z <sup>H</sup> B A* 2 1 ( ) − + <sup>Δ</sup> characterising the anhydrous salt *A Bx x* 1 2 , was defined as follows [50-52]:

$$
\Delta\_H(B^{z\_2^-})A^{z\_1^+} = \frac{1}{z\_1^{\*\prime}z\_2} \left[ \Delta H^{o}\_{f,298} \left( A\_{x\_1} B\_{x\_2} \right) - x\_1^{\*\prime} \Delta H^{o}\_{f,298} \left( A^{z\_1^+} \right) - x\_2^{\*\prime} \Delta H^{o}\_{f,298} \left( B^{z\_2^-} \right) \right] \tag{8}
$$

This parameter is analogous to the enthalpy of dissolution of an anhydrous salt per equivalent (characterised by the product of charges of the cation and anion). The relationships of the enthalpy of hydration of a salt to the number of the molecules of water and to the nature of the salt is:

$$
\Delta H\_{\text{hyd},298.15} = \alpha\_{A^{\ast \ast}} \, ^\ast \left( n \right)^a \cdot \left[ \Delta\_H (\mathcal{B}^{z\_2^-}) A^{z\_1^+} + b \right] + \beta \tag{9}
$$

 where *<sup>z</sup> <sup>A</sup>* α <sup>+</sup> and β are constants that depend on the nature of the cation *<sup>z</sup> A* <sup>1</sup> + present in the hydrate, n is the number of molecules of the water of crystallisation in the hydrate, and a and b are constants for all hydrates and respectively equal to 0.791 and -80.0 kJ mol-1 [50-52]. Eq. (9) shows that the enthalpy of hydration is closely related to the nature of the cation in the anhydrous salt, to the number of water molecules in the chemical formula and to the nature of the anhydrous salt.

The enthalpy of the formation of ice used in Eq. (7) [50-52] comes from Robie & Waldbaum [53] and is equal to -279.8 kJ mol-1. As this value is very different from those given in table 1, the values of coefficients ( *<sup>z</sup> <sup>A</sup>* α <sup>+</sup> , β, a and b) of Eq. (9) need to be determined with *H <sup>f</sup>* ,298 ° Δ H2O(ice)=-292.75 kJ mol-1 by a minimisation technique of the square difference between the measured and calculated heat of hydration. The values of a and b are assumed constants for all hydrated salts and are a= 0.612 and b= 17.07; Eq (9) becomes:

$$
\Delta H\_{hyd, 298.15} = \alpha\_{A^{z+}} \ast \left( n \right)^{0.612} \ast \left[ \Delta \left( B^{z\_2^{-}} \right) A^{z\_1^{+}} + 17.07 \right] + \beta \tag{10}
$$

The values of *<sup>z</sup> <sup>A</sup>* α <sup>+</sup> and β for the different salts of cation *<sup>z</sup> A* <sup>1</sup> + are given in table 3 with the respective standard deviations and % errors of the predicted enthalpy of formation of the hydrated salts.

344 Thermodynamics – Fundamentals and Its Application in Science

:

the enthalpy of hydration of the reaction (6) becomes

<sup>Δ</sup> characterising the anhydrous salt *A Bx x* 1 2

*z z*

1 2

,298.15

α

,298.15

all hydrated salts and are a= 0.612 and b= 17.07; Eq (9) becomes:

α

+

α

+

**5.1. Enthalpy of formation** 

hydration of a salt *A Bx x* 1 2

consideration and *H <sup>f</sup>* ,298

and to the nature of the salt is:

nature of the anhydrous salt.

the values of coefficients ( *<sup>z</sup> <sup>A</sup>*

*z z <sup>H</sup> B A* 2 1 ( ) − +

 where *<sup>z</sup> <sup>A</sup>* α

**5. Recommended models of prediction of thermodynamic entities** 

**Salts:** A model of the prediction of hydration enthalpy has been developed and tested on 349 different hydrates salts [50-52]. By considering the following reaction defining the

( ) () *o o <sup>o</sup> <sup>o</sup> H H A B nH O H A B n H H O hyd f xx f xx <sup>f</sup> <sup>c</sup>* 1 2 1 2 ,298 ,298 <sup>2</sup> ,298 ,298 2 ( ) Δ =Δ − −Δ − Δ\* () (7)

° Δ (H2O)(c) represents the enthalpy of formation of ice. A parameter

1 2 ,298 1 ,298 2 ,298

where n is the number of molecules of crystalline water contained in the hydrate under

( ) ( ) *z z o oo <sup>z</sup> <sup>z</sup> H f <sup>x</sup> <sup>x</sup> f f B A H AB x H A x H B*

<sup>1</sup> ( ) \* \* () \* − + + − Δ = Δ −Δ −Δ (8)

This parameter is analogous to the enthalpy of dissolution of an anhydrous salt per equivalent (characterised by the product of charges of the cation and anion). The relationships of the enthalpy of hydration of a salt to the number of the molecules of water

*<sup>z</sup>* ( )*<sup>a</sup> z z H nB hyd <sup>A</sup> <sup>H</sup> <sup>A</sup> <sup>b</sup>* 2 1

<sup>+</sup> and β are constants that depend on the nature of the cation *<sup>z</sup> A* <sup>1</sup>

\* \* ()

hydrate, n is the number of molecules of the water of crystallisation in the hydrate, and a and b are constants for all hydrates and respectively equal to 0.791 and -80.0 kJ mol-1 [50-52]. Eq. (9) shows that the enthalpy of hydration is closely related to the nature of the cation in the anhydrous salt, to the number of water molecules in the chemical formula and to the

The enthalpy of the formation of ice used in Eq. (7) [50-52] comes from Robie & Waldbaum [53] and is equal to -279.8 kJ mol-1. As this value is very different from those given in table 1,

H2O(ice)=-292.75 kJ mol-1 by a minimisation technique of the square difference between the measured and calculated heat of hydration. The values of a and b are assumed constants for

*<sup>z</sup>* ( ) ( ) *z z*

\* \* 17.07

*hyd <sup>A</sup> H nB* 2 1 *<sup>A</sup>* 0.612

2 1 1 2

*A B n H O A B nH O x x c xx* 1 2 1 2 2 () <sup>2</sup> + →− \* (6)

, was defined as follows [50-52]:

β

β

+

present in the

° Δ

− +

Δ = Δ ++ (9)

<sup>+</sup> , β, a and b) of Eq. (9) need to be determined with *H <sup>f</sup>* ,298

− +

Δ = Δ ++ (10)

An example is given by considering the data of all of the sodium salts displayed in Figure 1A and shows the great variation of the enthalpy of hydration having the same number as the crystalline water. The different straight lines are obtained from Eq. (10) with values of *z A* α <sup>+</sup> and β for Na+ and are plotted for the different values of the number of crystalline water molecules and the nature of the crystalline salts ( *z z <sup>H</sup> B A* 2 1 ( ) − + Δ ). For a given anhydrous salt (same value of *z z <sup>H</sup> B A* 2 1 ( ) − + Δ ), the enthalpy of hydration per mole of water appears to not be constant but decreases when the number of molecules of the crystalline water increases. Figure 1B displays the enthalpy of crystalline water versus the nature of the anhydrous salt showing the great disparity and shows that the enthalpy of the crystalline water in all of the sodic salts is more negative than those of the hexagonal ice. When the number of crystalline water molecules of any salt increases, the enthalpy of formation of the crystalline water decreases from large values (NaOH. H2O; *H <sup>f</sup>* ,298 ° Δ H2O= -309.25 kJ mol-1) to values close to those of ice (Na2Se. 16 H2O; *H <sup>f</sup>* ,298 ° Δ = -297.22 kJ mol-1) and are different from the average value for 61 Na salts ( *H <sup>f</sup>* ,298 ° Δ = -297.74 kJ mol-1).


1: the % error weighs the difference (in %) between the measured and predicted enthalpy of the formation of the hydrated salts.

**Table 3.** Values of *<sup>z</sup> <sup>A</sup>* α<sup>+</sup> , β, number of data, % error and statistical error per moles of crystalline water.

Thermodynamics of Hydration in Minerals: How to Predict These Entities 347

substitution, the framework density of the anhydrous zeolite, the average cation electronegativity in the exchange site and the intracrystalline water porosity, respectively. All of these parameters can be calculated from the chemical formula and the unit cell

The framework density (FD) [55, 56] represents the number of tetrahedral atoms per 1000 Å3

The electronegativity of site Az+, defined by ΔHO= (site A), represents the weighed average

*i 1*

*<sup>c</sup> i n A i i i 1 x n \*x* =

=

The parameter ΔHO= Mz+ (aq) characterises the electronegativity for a given cation Mz+ and is defined as the difference between the enthalpy of formation of the corresponding oxides

*H MO* ° Δ and the enthalpy of formation of the corresponding aqueous cation

( ) ( ) *i i i*

*x* ( )

*z z H aq f ix f i (c) aq <sup>i</sup> O M H MO - H M*

= °° + + <sup>1</sup> Δ =Δ Δ

where z is the charge of the cation Mz+ and x is the number of oxygen atoms combined with one atom of M in the oxide (x = z/2), such that the difference in Eq. (17) refers to one oxygen atom. A set of internally consistent values was presented by Vieillard & Mathieu [54] and is

The total intracrystalline pore volume, WP, was introduced by Barrer [57] for determining the volume of liquid water that can be recovered thorough the outgassing of the fully hydrated zeolite. Assuming that the unit cell volume of zeolitic water is the same as that of liquid water ( *Vuc H O*<sup>2</sup> . . = 29.89 Å3), during the hydration-dehydration processes, the water porosity (WP) can be calculated as the volume of liquid water in 1Å3 of crystal and is expressed as follows:

( )

*V*

29.89 \*

*WP*

*u c hydr*

*H O*

2 .. .

= (18)

*n*

=

of ΔHO= Mz+(aq) of nc different cations in the exchanged site A:

*H*

*Δ O (site.A)*

<sup>=</sup> <sup>=</sup> =

The number of oxygens balancing site Mz+ (in extra-framework sites) is then:

( ) *anh Al Fe Si FD n n n Unit cell Vol 1* . = ++ / \* 000 (14)

( ) *<sup>c</sup>*

*n \*x . Δ O M*

*A*

*x*

*i n <sup>z</sup> i i H aq*

*i*

<sup>=</sup> (16)

(15)

(17)

= +

volume for any zeolite.

and are obtained as:

( )*<sup>i</sup> f ix (c)*

( ) *<sup>i</sup> <sup>z</sup> f i aq H M*

given in Table 4.

° <sup>+</sup> Δ :

( )

**Figure 1.** A – Hydration enthalpy of sodium salts versus the parameter *z z <sup>H</sup> B A* 2 1 ( ) − + Δ ; B – Enthalpy of the formation of the crystalline water versus the parameter *z z <sup>H</sup> B A* 2 1 ( ) − + Δ .

**Zeolites:** The anhydrous compound X (reaction 1) is a zeolite with the following formula:

( ) *<sup>c</sup> Al Fe Si i n i n nn N i A Al Fe Si O* 1 = = , in which Al, Fe and Si represent the framework with different

channels occupied by various exchangeable cations present, which are represented by Az+. A compilation of the average hydration enthalpies per mole of water of 145 diversely originating zeolites measured using different technical methods [54] provided 76 data points from the transposed-temperature drop calorimetry (TTDC), 57 data points from the immersion calorimetry (IC), 6 data points from the phase equilibriums (PE), 5 data points from the gas-adsorption calorimetry (GAC) and 3 data points from the hydrofluoric acid calorimetry (HF).

From the selected average enthalpies of hydration of the different zeolites with their chemical formulas (145 data points), a regression equation for the enthalpy of hydration per mole of water is proposed [54]:

$$
\Delta H\_{Hydro-W} = -\left(\bigwedge\_{Si}^{Al}\bigvee\nolimits^\*\mathcal{C}^{[\![\![\![\![\![\![\![\![\![\!0.\!]\!]\!]\!]\!]\!]\!]\!]}\right)^\* \tag{11}
$$

where coefficients a and b are

( ) ( ) *anh.* ( ) *<sup>H</sup> a A* 4.664 2.787 \* *<sup>l</sup> Al Si \** 0.056 0 003 *FD . \* <sup>Δ</sup> O site A* <sup>=</sup> = − ++ − (12)

$$a = -1.52 + 0.0047 \ast \left(\Delta\_H \mathcal{O}^{\rightarrow} \text{site} \, A\right) \ast \left(\text{Al/Si}\right) \tag{13}$$

Coefficients of Eqs. (11) (12) and (13), such as the ratio Al/Si, Al/(Al+Si), FDanh., *ΔH ( O (site A) aq)* <sup>=</sup> and WP, represent the extra-framework charge, the ratio of tetrahedral substitution, the framework density of the anhydrous zeolite, the average cation electronegativity in the exchange site and the intracrystalline water porosity, respectively. All of these parameters can be calculated from the chemical formula and the unit cell volume for any zeolite.

346 Thermodynamics – Fundamentals and Its Application in Science

**Figure 1.** A – Hydration enthalpy of sodium salts versus the parameter *z z*

**Zeolites:** The anhydrous compound X (reaction 1) is a zeolite with the following formula:

, in which Al, Fe and Si represent the framework with different

channels occupied by various exchangeable cations present, which are represented by Az+. A compilation of the average hydration enthalpies per mole of water of 145 diversely originating zeolites measured using different technical methods [54] provided 76 data points from the transposed-temperature drop calorimetry (TTDC), 57 data points from the immersion calorimetry (IC), 6 data points from the phase equilibriums (PE), 5 data points from the gas-adsorption calorimetry (GAC) and 3 data points from the hydrofluoric acid

From the selected average enthalpies of hydration of the different zeolites with their chemical formulas (145 data points), a regression equation for the enthalpy of hydration per

> *Hyd W* ( ) *a b WP <sup>H</sup> Al Si <sup>e</sup>* \* \* <sup>−</sup>

*b O* 1.52 0.0047 \* *<sup>H</sup> site <sup>A</sup>* \* ( ) *Al Si* <sup>=</sup>

Coefficients of Eqs. (11) (12) and (13), such as the ratio Al/Si, Al/(Al+Si), FDanh.,

<sup>=</sup> and WP, represent the extra-framework charge, the ratio of tetrahedral

( ) ( ) *anh.* ( ) *<sup>H</sup> a A* 4.664 2.787 \* *<sup>l</sup> Al Si \** 0.056 0 003 *FD . \* <sup>Δ</sup> O site A* <sup>=</sup> = − ++ − (12)

the formation of the crystalline water versus the parameter *z z*

( ) *<sup>c</sup> Al Fe Si*

*A Al Fe Si O*

mole of water is proposed [54]:

where coefficients a and b are

*i n nn N*

*i n*

=

*i*

=

1

calorimetry (HF).

*ΔH ( O (site A) aq)*

  *<sup>H</sup> B A* 2 1 ( ) − +

*<sup>H</sup> B A* 2 1 ( ) − + Δ .

<sup>+</sup> Δ =− (11)

=− + Δ (13)

Δ ; B – Enthalpy of

The framework density (FD) [55, 56] represents the number of tetrahedral atoms per 1000 Å3 and are obtained as:

$$FD\_{amb.} = \left[ \left( n\_{Al} + n\_{Fe} + n\_{Si} \right) / \text{Limit cell Vol} \right] \* 1000 \tag{14}$$

The electronegativity of site Az+, defined by ΔHO= (site A), represents the weighed average of ΔHO= Mz+(aq) of nc different cations in the exchanged site A:

$$
\Delta\_H \mathbf{O}^{\mp}(\text{site}.A) = \frac{\sum\_{i=1}^{i=n\_c} n\_i \text{"\infty} \left(\Delta\_H \mathbf{O}^{\mp} \mathbf{M}\_{\text{aq}}^{z\_i+}\right)}{\mathbf{x}\_A} \tag{15}
$$

The number of oxygens balancing site Mz+ (in extra-framework sites) is then:

$$\mathbf{x}\_A = \sum\_{i=1}^{i=n\_c} \mathbf{n}\_i \mathbf{^\ast x}\_i \tag{16}$$

The parameter ΔHO= Mz+ (aq) characterises the electronegativity for a given cation Mz+ and is defined as the difference between the enthalpy of formation of the corresponding oxides ( )*<sup>i</sup> f ix (c) H MO* ° Δ and the enthalpy of formation of the corresponding aqueous cation ( ) *<sup>i</sup> <sup>z</sup> f i aq H M* ( ) ° <sup>+</sup> Δ :

$$
\Delta\_H O^\equiv M\_{\rm aq}^{z\_{i\cdot}+} = \bigvee\_{\mathbf{X}\_i} \left[ \Delta H\_f^\circ \left( M\_i \mathcal{O}\_{x\_i} \right)\_{\text{(c)}} \cdot \Delta H\_f^\circ \left( M\_i^{z\_{i\cdot}+} \right)\_{\text{(aq)}} \right] \tag{17}
$$

where z is the charge of the cation Mz+ and x is the number of oxygen atoms combined with one atom of M in the oxide (x = z/2), such that the difference in Eq. (17) refers to one oxygen atom. A set of internally consistent values was presented by Vieillard & Mathieu [54] and is given in Table 4.

The total intracrystalline pore volume, WP, was introduced by Barrer [57] for determining the volume of liquid water that can be recovered thorough the outgassing of the fully hydrated zeolite. Assuming that the unit cell volume of zeolitic water is the same as that of liquid water ( *Vuc H O*<sup>2</sup> . . = 29.89 Å3), during the hydration-dehydration processes, the water porosity (WP) can be calculated as the volume of liquid water in 1Å3 of crystal and is expressed as follows:

$$\text{WPP} = \frac{\text{29.89\*} \,\, n\_{H\_2O}}{\left(V\_{u.c.hydr.}\right)}\tag{18}$$

where Vu.c. hyd. represents the unit cell volume of a hydrated zeolite containing nH2O water molecules. Assuming that the number of water molecules nH2O varies from 0 to nH2O max, the unit cell volume of a zeolite can be directly related to the number of water molecules. Although such variations in the unit -cell volume do occur, let us assume a linear variation in the unit-cell volume as a function of the number of water molecules:

$$V\_{u.c.hyd} = V\_{u.c.amh.} + k^\* \ n\_{H\_2O} \tag{19}$$

Thermodynamics of Hydration in Minerals: How to Predict These Entities 349

**Figure 2.** Predicted enthalpy of the hydration of zeolitic water calculated from Eq. (11) versus the

clinoptilolites (-Ca, -Na and –K) with the following chemical formula given by [22]:


A detailed computation of the enthalpy of hydration is given as an example for three natural

The values of the ratio Al/Si yield a constant value of 0.235. The electronegativity of the exchangeable site Az+ for clinoptilolite-Na is, for example, the average electronegativity of

The unit-cell volume of the anhydrous clinoptilolite–K (V u.c.) is 2019.51 Å3 [59]. This means that there are 35.95 tetrahedral atoms (6.95+0.1+29) in 2019.51 Å3; which indicates that FDanh.= 17.8 tetrahedral atoms in 1000Å3. The unit-cell volume of the hydrated clinoptilolite-K (at saturation or fully hydrated) with 18.5 moles of zeolitic waters is 2089.50 Å3 [54], which is slightly greater than the anhydrous form and allows the setting of the relationship of

( ) ( )( ) ( ) 2 2 2.607 0.375 \* 0.12 0.29 \* 3.39 *H aq H aq H aq H aq*

<sup>=</sup> + ++ <sup>=</sup> (21)

*\* Δ O Na Δ OK \* Δ O Ca Δ O Mg*

= + = + = + = +

experimental enthalpy of hydration using different technical measurements.

four cations (Na+, K+, Ca2+ and Mg2+ (table 4)) as follows:

water porosity versus the number of hydration water moles:

*H*

*Δ O (site.A)*

i.e., *ΔHO (site.A)* <sup>=</sup> =51.42 kJ mol-1.

where Vu.c. anh. represents the unit-cell volume of the anhydrous zeolite. Parameter k weights the variation in the unit-cell volume between an anhydrous zeolite and a fully hydrated zeolite per one water molecule. This parameter can be calculated from the available unit-cell volumes of anhydrous and fully hydrated zeolites. Thus, knowledge of the number of water molecules and the unit-cell volume of anhydrous and hydrated zeolites is required for the calculation of the water porosity:


( ) *H O u c anh H O n WP V kn* 2 <sup>2</sup> .. . 29.89 \* \* <sup>=</sup> <sup>+</sup> (20)

**Table 4.** Parameter ΔHO=Mz+ (aq) of selected cations [54].

By eliminating 9 erroneous data points [54], the regression coefficient for 136 data is R2= 0.880, with a standard error of ± 3.46 kJ mol-1 for all data, regardless of the nature of the experimental data. Figure 2 displays the predicted enthalpy of hydration of zeolitic water calculated from Eqs. (11) (12) and (13) versus the experimental enthalpy of hydration coming from the different technical measurements. Regression coefficients and standard errors are different within the two main groups: R2= 0.88 and std. err. = ± 3.41 kJ mol-1 for TTDC data (67 data points), R2= 0.658 and std err. = ± 2.66 kJ mol-1 for the "IC" data (57 data points).

When Al/Si and Al/(Al+Si) = 0, the hypothetical integral enthalpy of water in zeolites obtained by the extrapolation of Eq. (11) is 0, which corresponds to the enthalpy of hydration of water in cordierite (2.2 ± 1.6 kJ mol-1 from Carey [58]). Unlike zeolites, the water in cordierite is not coordinated by cations, so the molecular environment of H2O in cordierite is similar to that of H2O in water.

**Figure 2.** Predicted enthalpy of the hydration of zeolitic water calculated from Eq. (11) versus the experimental enthalpy of hydration using different technical measurements.

A detailed computation of the enthalpy of hydration is given as an example for three natural clinoptilolites (-Ca, -Na and –K) with the following chemical formula given by [22]:


The values of the ratio Al/Si yield a constant value of 0.235. The electronegativity of the exchangeable site Az+ for clinoptilolite-Na is, for example, the average electronegativity of four cations (Na+, K+, Ca2+ and Mg2+ (table 4)) as follows:

$$\Delta\_{\text{II}}O^{\mp}(siteA) = \frac{2.607^{\star} \left(\Delta\_{\text{II}}O^{\mp} \text{M}d\_{\text{aq}}^{\star}\right) + 0.375^{\star} \left(\Delta\_{\text{II}}O^{\mp} \text{K}\_{\text{aq}}^{\star}\right) + 0.12^{\star} \left(\Delta\_{\text{II}}O^{\mp} \text{C}d\_{\text{aq}}^{2+}\right) + 0.29^{\star} \left(\Delta\_{\text{II}}O^{\mp} \text{M}g\_{\text{aq}}^{2+}\right)}{3.39} \tag{21}$$

i.e., *ΔHO (site.A)* <sup>=</sup> =51.42 kJ mol-1.

348 Thermodynamics – Fundamentals and Its Application in Science

calculation of the water porosity:

**Table 4.** Parameter ΔHO=Mz+ (aq) of selected cations [54].

cordierite is similar to that of H2O in water.

where Vu.c. hyd. represents the unit cell volume of a hydrated zeolite containing nH2O water molecules. Assuming that the number of water molecules nH2O varies from 0 to nH2O max, the unit cell volume of a zeolite can be directly related to the number of water molecules. Although such variations in the unit -cell volume do occur, let us assume a linear variation

where Vu.c. anh. represents the unit-cell volume of the anhydrous zeolite. Parameter k weights the variation in the unit-cell volume between an anhydrous zeolite and a fully hydrated zeolite per one water molecule. This parameter can be calculated from the available unit-cell volumes of anhydrous and fully hydrated zeolites. Thus, knowledge of the number of water molecules and the unit-cell volume of anhydrous and hydrated zeolites is required for the

( )

*V kn*

Ions ΔHO=Mz+(aq) Ions ΔHO=Mz+(aq) (kJ mol-1) (kJ mol-1)

Li+(aq) -40.90 Ca+2(aq) -92.1 Na+(aq) 65.80 Mg+2(aq) -134.6 K+(aq) 141.00 Fe+2(aq) -181.58 Rb+(aq) 163.24 Mn+2(aq) -164.4 Cs+(aq) 170.0 Zn+2(aq) -197.07 Ba+2(aq) -15.6 La+3(aq) -124.8 Sr+2(aq) -40.4 H+(aq) -285.83

By eliminating 9 erroneous data points [54], the regression coefficient for 136 data is R2= 0.880, with a standard error of ± 3.46 kJ mol-1 for all data, regardless of the nature of the experimental data. Figure 2 displays the predicted enthalpy of hydration of zeolitic water calculated from Eqs. (11) (12) and (13) versus the experimental enthalpy of hydration coming from the different technical measurements. Regression coefficients and standard errors are different within the two main groups: R2= 0.88 and std. err. = ± 3.41 kJ mol-1 for TTDC data (67

When Al/Si and Al/(Al+Si) = 0, the hypothetical integral enthalpy of water in zeolites obtained by the extrapolation of Eq. (11) is 0, which corresponds to the enthalpy of hydration of water in cordierite (2.2 ± 1.6 kJ mol-1 from Carey [58]). Unlike zeolites, the water in cordierite is not coordinated by cations, so the molecular environment of H2O in

data points), R2= 0.658 and std err. = ± 2.66 kJ mol-1 for the "IC" data (57 data points).

29.89 \* \* <sup>=</sup> <sup>+</sup>

*u c anh H O*

<sup>2</sup> .. .

*n*

*H O*

2

*V V kn u c hyd u c anh H O*<sup>2</sup> .. .. . = + \* (19)

(20)

in the unit-cell volume as a function of the number of water molecules:

*WP*

The unit-cell volume of the anhydrous clinoptilolite–K (V u.c.) is 2019.51 Å3 [59]. This means that there are 35.95 tetrahedral atoms (6.95+0.1+29) in 2019.51 Å3; which indicates that FDanh.= 17.8 tetrahedral atoms in 1000Å3. The unit-cell volume of the hydrated clinoptilolite-K (at saturation or fully hydrated) with 18.5 moles of zeolitic waters is 2089.50 Å3 [54], which is slightly greater than the anhydrous form and allows the setting of the relationship of water porosity versus the number of hydration water moles:

$$\text{WPP} = \frac{29.89 \text{ \* } n\_{H\_2O}}{\left(2019.51 + 3.78 \text{ \* } n\_{H\_2O}\right)}\tag{22}$$

Thermodynamics of Hydration in Minerals: How to Predict These Entities 351

+

 (23)

= 0 (corresponding to an anhydrous zeolite) to a maximum

<sup>−</sup> <sup>=</sup> Δ= Δ (24)

*u c anh H O u c anh H O*

2 2

= (26)

.. .

*V kn bV n*

.. . .. .

( )

<sup>+</sup>

*u c anh H O*

\*

2

2

(25)

The bulk enthalpy of hydration or the integral enthalpy of hydration can be used as the

( ) ( ) ( ) *Hyd Z Hyd W H O H O a b WP H Hn Al <sup>n</sup> Si <sup>e</sup>* 2 2 \* \* \*\* − −

As the integral hydration enthalpy *HHyd Z*<sup>−</sup> <sup>Δ</sup> is the integral function of the partial molar

value of *H O <sup>n</sup>* <sup>2</sup> max (corresponding to the fully hydrated zeolite), the partial molar enthalpy of hydration relative to the liquid water, Δ*<sup>H</sup>* , which is the derivative function of *HHyd Z*<sup>−</sup> <sup>Δ</sup> ,

> *H H dn n*

2

*n*

*H O*

2

<sup>+</sup> <sup>+</sup> Δ =− <sup>+</sup>

*Si V kn*

The fractional water content can be obtained if the maximum number of zeolitic waters is

*n n*

From values of table 5, the partial molar enthalpy of clinoptilolites K, Na and Ca can be calculated versus the fractional water content and plotted in Figure 3 with the experimental partial molar enthalpies measured from immersion calorimetry [62] and from thermogravimetry [22]. The modelled values of the partial molar enthalpy for the three clinoptilolites appear to be closer to the data from the immersion calorimetry than from thermogravimetry. The second important point is the fact that the calculated and

The partial molar enthalpy of hydration per mole of H2O increases smoothly from low water content to high θ and indicates that the H2O in clinoptilolites occupies a continuum of energetic states. This is not the case for chabazite-Ca [23], which exhibits three energetically

This shows that, within different zeolites with the same exchangeable cation, the partial molar enthalpy as a function of the hydration degree may display different functions. This is the reason why the fundamental relationship verifying the enthalpy of hydration by means

θ

experimental partial molar enthalpy of hydration has a similar behaviour.

*H O H O*

2 <sup>2</sup> ,max.

*H O n n*

0 <sup>1</sup> \* <sup>=</sup>

\* 29.89 \* \* \* \* \* 1

2

product of the enthalpy of hydration of the water times the number of zeolitic water:

Δ =Δ = −

*Hyd Z <sup>n</sup> H O*

( ) ( )

29.89 \* \*

*a b*

*e*

known at the saturation state and is equal to:

enthalpy hydration from *H O <sup>n</sup>* <sup>2</sup>

*H Al*

distinct types of H2O.

becomes:

From these previous examples, the parameters requested for the prediction of the enthalpy of hydration are displayed in table 5 for the three clinoptilolites (-Na, -K and –Ca).


1 - [60]; 2 - [59]

**Table 5.** Unit-cell volume of hydrated and anhydrous clinoptilolite, factor k (eq. 20), mean ΔHO=Az+ (site A), framework density and ratio Al/Si for 3 clinoptilolites (-Na, K and –Ca).

From the values given in table 5, the enthalpy of hydration of the zeolitic water, *HHyd W*<sup>−</sup> Δ , can be computed with Eqs. (11), (12) and (13) and are close to experimental values from [22] (Table 6).

The model of the computation of the enthalpy of hydration of the hydration water represents a very useful tool and contributes to the knowledge of enthalpies of formation of hydrated zeolites from anhydrous ones [61].

A useful method for avoiding complications at the outset of a thermodynamic analysis is to emphasise the H2O rather than the zeolite structure. Using such an approach, the equilibriums can be considered as the equilibrium between the H2O in the fluid-vapour phase and the H2O in the zeolite. According to the definition of equilibrium, the chemical potential of H2O in the vapour phase must equal the chemical potential of the H2O component in the zeolite. Any series of measurements of the amount of H2O in zeolite at a known fugacity (or partial pressure of H2O under ideal gas conditions) provides the basis for the thermodynamic description of the system. The most elegant approach to develop a thermodynamic formulation is to know the partial molar enthalpy of hydration, Δ*H* , and the partial molar entropy Δ*S* . A constant partial molar enthalpy of hydration indicates an ideal mixing of the H2O.


**Table 6.** Predicted and experimental enthalpy of hydration, *HHyd W*<sup>−</sup> Δ of clinoptilolites –Na, -K and –Ca.

The bulk enthalpy of hydration or the integral enthalpy of hydration can be used as the product of the enthalpy of hydration of the water times the number of zeolitic water:

350 Thermodynamics – Fundamentals and Its Application in Science

1 - [60]; 2 - [59]

(Table 6).

ideal mixing of the H2O.

Clinoptilolite.-Na -22.20

*WP*

( )

29.89 \* 2019.51 3.78 \* <sup>=</sup> <sup>+</sup>

From these previous examples, the parameters requested for the prediction of the enthalpy

Zeolite minerals V u.c. anh. Vu.c. hydr. k ΔHO= FD. WP Al/Si

Clinoptilolite -Na 2000.00 1 2140.00 1 6.35 51.39 18.02 0.31 0.23 Clinoptilolite -K 2019.51 2 2089.50 2 3.78 109.16 17.80 0.26 0.24 Clinoptilolite -Ca 2019.51 2 2089.50 2 3.20 -64.60 17.78 0.31 0.23

**Table 5.** Unit-cell volume of hydrated and anhydrous clinoptilolite, factor k (eq. 20), mean ΔHO=Az+ (site

From the values given in table 5, the enthalpy of hydration of the zeolitic water, *HHyd W*<sup>−</sup> Δ , can be computed with Eqs. (11), (12) and (13) and are close to experimental values from [22]

The model of the computation of the enthalpy of hydration of the hydration water represents a very useful tool and contributes to the knowledge of enthalpies of formation of

A useful method for avoiding complications at the outset of a thermodynamic analysis is to emphasise the H2O rather than the zeolite structure. Using such an approach, the equilibriums can be considered as the equilibrium between the H2O in the fluid-vapour phase and the H2O in the zeolite. According to the definition of equilibrium, the chemical potential of H2O in the vapour phase must equal the chemical potential of the H2O component in the zeolite. Any series of measurements of the amount of H2O in zeolite at a known fugacity (or partial pressure of H2O under ideal gas conditions) provides the basis for the thermodynamic description of the system. The most elegant approach to develop a thermodynamic formulation is to know the partial molar enthalpy of hydration, Δ*H* , and the partial molar entropy Δ*S* . A constant partial molar enthalpy of hydration indicates an

predicted (kJ mol-1 H2O) measured (kJmol-1 H2O) References

Clinoptilolite -K -20.24 -23.97 [22] Clinoptilolite-Ca -29.03 -32.87 [22]

**Table 6.** Predicted and experimental enthalpy of hydration, *HHyd W*<sup>−</sup> Δ of clinoptilolites –Na, -K and –Ca.


of hydration are displayed in table 5 for the three clinoptilolites (-Na, -K and –Ca).

(Å3) (Å3) Site A Anh.

kJmol-1

A), framework density and ratio Al/Si for 3 clinoptilolites (-Na, K and –Ca).

Zeolite minerals *HHyd W*<sup>−</sup> Δ *HHyd W*<sup>−</sup> Δ

hydrated zeolites from anhydrous ones [61].

*H O*

*n*

*H O*

(22)

2

*n* 2

$$\Delta \tilde{H}\_{Hyd-Z} = \left(\Delta H\_{Hyd-W}\right) \* \left. n\_{H\_2O} \right| = \left\lfloor -\left(\bigvee\_{S\bar{I}} \right) \* \mathcal{C}^{\left(a+b\right.} \mathcal{WPP}\right) \right\rfloor \* \left. n\_{H\_2O} \right\rfloor \tag{23}$$

As the integral hydration enthalpy *HHyd Z*<sup>−</sup> <sup>Δ</sup> is the integral function of the partial molar enthalpy hydration from *H O <sup>n</sup>* <sup>2</sup> = 0 (corresponding to an anhydrous zeolite) to a maximum value of *H O <sup>n</sup>* <sup>2</sup> max (corresponding to the fully hydrated zeolite), the partial molar enthalpy of hydration relative to the liquid water, Δ*<sup>H</sup>* , which is the derivative function of *HHyd Z*<sup>−</sup> <sup>Δ</sup> , becomes:

$$
\Delta \tilde{H}\_{Hyd-Z} = \frac{1}{n\_{H\_2O}} \int\_{n=0}^{n=n\_{H\_2O}} \Delta \overline{H} \, ^\ast \, dn \tag{24}
$$

$$\Delta \overline{H} = \left| - \left( {}^{M}\acute{\bigvee}\_{Si} \right) \* \mathbf{e} \left( \overline{\left( {}^{V\_{\
u,c.mh.}} + {}^{k\*}\boldsymbol{n}\_{H\_2O} \right)} \right) \right| \begin{aligned} \left( {}^{a}a + b \, ^{a}\left( \overline{\left( {}^{V\_{\
u,c.mh.}} + {}^{k\*}\boldsymbol{n}\_{H\_2O} \right)} \right) \right) \Bigg| \begin{aligned} \left( {}^{a}a + b \, ^{a}\boldsymbol{n}\_{H\_2O} \right) \end{aligned} \tag{25}$$

The fractional water content can be obtained if the maximum number of zeolitic waters is known at the saturation state and is equal to:

$$\theta = \frac{\mathfrak{n}\_{H\_2O}}{\mathfrak{n}\_{H\_2O, \text{max.}}} \tag{26}$$

From values of table 5, the partial molar enthalpy of clinoptilolites K, Na and Ca can be calculated versus the fractional water content and plotted in Figure 3 with the experimental partial molar enthalpies measured from immersion calorimetry [62] and from thermogravimetry [22]. The modelled values of the partial molar enthalpy for the three clinoptilolites appear to be closer to the data from the immersion calorimetry than from thermogravimetry. The second important point is the fact that the calculated and experimental partial molar enthalpy of hydration has a similar behaviour.

The partial molar enthalpy of hydration per mole of H2O increases smoothly from low water content to high θ and indicates that the H2O in clinoptilolites occupies a continuum of energetic states. This is not the case for chabazite-Ca [23], which exhibits three energetically distinct types of H2O.

This shows that, within different zeolites with the same exchangeable cation, the partial molar enthalpy as a function of the hydration degree may display different functions. This is the reason why the fundamental relationship verifying the enthalpy of hydration by means

of the chemical composition and the accurate knowledge of the unit-cell parameters of anhydrous and hydrated zeolites needs to be improved with new data of hydration enthalpy of the partially hydrated zeolites on the one hand and a better fit of the effective water porosity as a function of the hydration degree on the other hand. Then, the thermodynamic description of the hydration-dehydration process can be modelled as a function of pressure and temperature with the contribution of predicted enthalpies of hydration and entropies of hydration.

Thermodynamics of Hydration in Minerals: How to Predict These Entities 353

*m l m i Smectite n H O Smectite n H O* 2 () 2 () + →− \* \* (27)

*Gx xW xs as as G* = − \*1 \* ( ) (28)

° <sup>Δ</sup> , is determined from *GHyd*,298

° H2O(i) equal to 55.02 J K-1 mol-1 [67]. From calculated

° Δ by

*G x G x xW Hyd hs Hyd hs hs G* ( ) ,298 \* \* 1 \* ° Δ=Δ + − (29)

*H x H x xW Hyd hs Hyd hs hs H* ( ) ,298 \* \* 1 \* ° Δ =Δ + − (30)

,298 ,298 2 ( ) ,298 2 \* Δ =Δ − Δ (31)

° Δ is the standard free energy of hydration. The previous parameter and WG

are determined only from adsorption isotherms by assuming a maximal number of moles of interlayer water fixed to 4.5 moles. As there are no experimental measurements of heat of

° <sup>Δ</sup> , the standard enthalpy of formation of the interlayer water, *H HO <sup>f</sup>* ,298 2 *<sup>i</sup>* ( ) ° <sup>Δ</sup> , is

Vidal et al. [68-69] have assumed that a smectite could be considered strictly regular solid solutions between four end-members with 0.7 H2O, 2 H2O, 4 H2O, and 7 H2O. Those compositions correspond to four different hydrated states (with 0, 1, 2, and 3 water sheets, respectively). For any strictly regular solid solution, the integral enthalpy of hydration is

*H n H HO H HO hyd m f* ( ) *i f* ( ) *liq*

00 0

hydration-dehydration on SWy-1 [15-17, 64] [18, 19, 30] [20] [34] [24, 26, 65] [25] and on

All recent predictive models of the hydration of smectites are based on an approach that uses the solid solution initiated by Ransom & Helgeson [66]. The hydration of a smectite is considered through the following reaction between hydrated and dehydrated end-members:

where nm represents the maximal number of moles of water that can be included in the smectite on the basis of a half-cell (i.e., O10(OH)2). One considers the interlayer water H2O(i) (where subscript i stands for interlayer) as a variable weighing the hydration ratio. The amount of interlayer water in a smectite is proportional to the mole fraction of hydrated end-member: x\*H2O= nm\*xhs=(1-nm)\*xas; to the mole fraction of the hydrous end-member, xhs or to the anhydrous end-members xas. Ransom & Helgeson [66] considered the system hydrated smectite – anhydrous smectite + bulk water as a strictly regular binary solid solution by considering the free energy only. Thus, the excess free energy of mixing can be

where WG is an excess mixing constant. The integral Gibbs free energy of hydration is:

adsorption, the standard enthalpy of hydration, *HHyd*,298

vermiculites [21].

expressed as follows:

in which *GHyd*,298

*HHyd*,298

with

given in table 7.

assuming a constant value of *S*<sup>298</sup>

expressed as the following:

**Figure 3.** Partial molar enthalpy of hydration of K-, Na- and Ca- clinoptilolites derived from a model in [54], immersion calorimetry [62] and thermogravimetry results [22]

### **Clay minerals:**

Unlike zeolites, smectites are clearly not inert phases, as the particles size increases with the relative humidity [15]. Adsorbed water is distributed throughout the interlayer space, the outer surfaces of particles and the open pores space in the sample. [15-17] provided measured values of the surface covering waters for both the hydration and dehydration reactions of a set of eight homo-ionic SWy-1 montmorillonites saturated by alkali and alkaline-earth cations from their BET specific surface area (Table 1) and their basal spacing variations with relative humidity. To quantify the effective amount of water involved in the hydration reaction, the amount of poral H2O must be quantified and subtracted from the total amount of H2O taken up by the clay sample. The number of external surface water molecules (expressed in mmol/g dry clay) can be expressed as a function of relative humidity [63]. There are numerous papers about the adsorption-desorption isotherms performed on various clays minerals, but few are devoted to the acquisition of enthalpy of hydration-dehydration on SWy-1 [15-17, 64] [18, 19, 30] [20] [34] [24, 26, 65] [25] and on vermiculites [21].

All recent predictive models of the hydration of smectites are based on an approach that uses the solid solution initiated by Ransom & Helgeson [66]. The hydration of a smectite is considered through the following reaction between hydrated and dehydrated end-members:

$$\text{Sametite} + \text{n}\_{\text{m}} \, \, ^\ast H\_2\text{O}\_{\text{(l)}} \to \text{Sametite} - \text{n}\_{\text{m}} \, ^\ast H\_2\text{O}\_{\text{(i)}} \tag{27}$$

where nm represents the maximal number of moles of water that can be included in the smectite on the basis of a half-cell (i.e., O10(OH)2). One considers the interlayer water H2O(i) (where subscript i stands for interlayer) as a variable weighing the hydration ratio. The amount of interlayer water in a smectite is proportional to the mole fraction of hydrated end-member: x\*H2O= nm\*xhs=(1-nm)\*xas; to the mole fraction of the hydrous end-member, xhs or to the anhydrous end-members xas. Ransom & Helgeson [66] considered the system hydrated smectite – anhydrous smectite + bulk water as a strictly regular binary solid solution by considering the free energy only. Thus, the excess free energy of mixing can be expressed as follows:

$$G\_{xs} = \boldsymbol{\pi}\_{as} \, \* \left(1 - \boldsymbol{\pi}\_{as}\right) \, \* \, W\_G \tag{28}$$

where WG is an excess mixing constant. The integral Gibbs free energy of hydration is:

$$
\Delta \tilde{G}\_{Hyd} = \boldsymbol{\chi}\_{hs} \, ^\ast \Delta G\_{Hyd, 298}^\circ + \boldsymbol{\chi}\_{hs} \, ^\ast \left(1 - \boldsymbol{\chi}\_{hs}\right) ^\ast \boldsymbol{W}\_G \tag{29}
$$

in which *GHyd*,298 ° Δ is the standard free energy of hydration. The previous parameter and WG are determined only from adsorption isotherms by assuming a maximal number of moles of interlayer water fixed to 4.5 moles. As there are no experimental measurements of heat of adsorption, the standard enthalpy of hydration, *HHyd*,298 ° <sup>Δ</sup> , is determined from *GHyd*,298 ° Δ by assuming a constant value of *S*<sup>298</sup> ° H2O(i) equal to 55.02 J K-1 mol-1 [67]. From calculated *HHyd*,298 ° <sup>Δ</sup> , the standard enthalpy of formation of the interlayer water, *H HO <sup>f</sup>* ,298 2 *<sup>i</sup>* ( ) ° <sup>Δ</sup> , is given in table 7.

Vidal et al. [68-69] have assumed that a smectite could be considered strictly regular solid solutions between four end-members with 0.7 H2O, 2 H2O, 4 H2O, and 7 H2O. Those compositions correspond to four different hydrated states (with 0, 1, 2, and 3 water sheets, respectively). For any strictly regular solid solution, the integral enthalpy of hydration is expressed as the following:

$$
\Delta \tilde{H}\_{Hyd} = \mathbf{x}\_{hs} \, ^\ast \Delta H\_{Hyd, 298}^\circ + \mathbf{x}\_{hs} \, ^\ast \left(1 - \mathbf{x}\_{hs}\right) ^\ast \, ^\ast W\_H \tag{30}
$$

with

352 Thermodynamics – Fundamentals and Its Application in Science

hydration and entropies of hydration.

of the chemical composition and the accurate knowledge of the unit-cell parameters of anhydrous and hydrated zeolites needs to be improved with new data of hydration enthalpy of the partially hydrated zeolites on the one hand and a better fit of the effective water porosity as a function of the hydration degree on the other hand. Then, the thermodynamic description of the hydration-dehydration process can be modelled as a function of pressure and temperature with the contribution of predicted enthalpies of

**Figure 3.** Partial molar enthalpy of hydration of K-, Na- and Ca- clinoptilolites derived from a model in

Unlike zeolites, smectites are clearly not inert phases, as the particles size increases with the relative humidity [15]. Adsorbed water is distributed throughout the interlayer space, the outer surfaces of particles and the open pores space in the sample. [15-17] provided measured values of the surface covering waters for both the hydration and dehydration reactions of a set of eight homo-ionic SWy-1 montmorillonites saturated by alkali and alkaline-earth cations from their BET specific surface area (Table 1) and their basal spacing variations with relative humidity. To quantify the effective amount of water involved in the hydration reaction, the amount of poral H2O must be quantified and subtracted from the total amount of H2O taken up by the clay sample. The number of external surface water molecules (expressed in mmol/g dry clay) can be expressed as a function of relative humidity [63]. There are numerous papers about the adsorption-desorption isotherms performed on various clays minerals, but few are devoted to the acquisition of enthalpy of

[54], immersion calorimetry [62] and thermogravimetry results [22]

**Clay minerals:** 

$$
\Delta H\_{hyd,298}^{0} = \eta\_m \, ^\ast \left[ \Delta H\_{f,298}^{0} \left( H\_2O\_{\text{(i)}} \right) - \Delta H\_{f,298}^{0} \left( H\_2O\_{\text{liq}} \right) \right] \tag{31}
$$

The hydration enthalpy is retrieved from data of [24] for montmorillonite –Na and –K and derived from differential heats of adsorption obtained from the measurements of the heat of immersion [18, 20] for Na-, Ca- and Mg- montmorillonite. Thus, an integral enthalpy of hydration *H HO <sup>f</sup>* ,298 2 *<sup>i</sup>* ( ) ° <sup>Δ</sup> is provided for each of the four solid solutions with nm equal to 0.7, 2, 4 and 7 and is given in table 7. Vieillard et al. [63] consider the hydration of a smectite to be an asymmetrical regular solid solution between anhydrous and hydrated smectite.

Thermodynamics of Hydration in Minerals: How to Predict These Entities 355

where the coefficient -24.39 is computed from *H HO <sup>f</sup>* ,298 2 *<sup>i</sup>* ( ) ° <sup>Δ</sup> given in Table 7. The curves 2b, 2c and 2d (green dotted lines) displayed in Figure 4A come together from the model of Vidal et al. [68] and are graphically represented by four truncated lines with their respective

( ) hs *H x Hyd*,298 *hs* 0 layer 0 x 0.7 : < < Δ =−9.716 \* (34)

( ) *Hyd* ( ) *H xhs hs hs x x* hs ,298 1 layer 0<x <2.0 : Δ =− + − 20.14 \* \* 1 \* 10.0 (35)

**Figure 4.** A – Integral enthalpy of hydration versus the number of interlayer cation; B – Standard

The curve 3 (black line) in Figure 4A comes from the model of Vieillard et al. [63] and is

It appears that integral enthalpy of hydration modelled by Vieillard et al. [63] (black line) based on experimental measurements of heat of adsorption from [15-17] encompasses nearly all experimental data. As the standard entropy of hydration is assumed constant in [66] and [68] models, the integral hydration enthalpies provided by these two previous models are nearly linear. Those obtained by [68] exhibit three sections of curves corresponding to the

In Figure 4B, a relationship between the standard state enthalpy of interlayer water with the nature of the interlayer cation located in the SWy-1 montmorillonite with a constant layer

enthalpy of formation of interlayer water versus the nature of the interlayer cation.

represented by the following equation:

three states of hydration (1, 2 and 3 water layers).

( ) *Hyd* ( ) *H xhs hs hs x x* hs ,298 2 layers 0 x 4.0 : < < Δ =− + − 32.0 \* \* 1 \* 10.0 (36)

( ) ( ) ,298 3 layers 0 .0 : 45.5 \* \* 1 \* 10.0 *hs Hyd hs hs hs x7 H* < < Δ =− + − *xx x* (37)

*Hyd* ( ) ( ) *H xx x x x hs hs hs hs hs* Δ =− + − 23.0 \* \* 1 \* 8.0 \* 33.0 \* 1 − − (38)

equations:


1 –[66]; 2 – [68]; 3 - [63].

**Table 7.** Standard thermodynamic enthalpy *H HO f i* ,298 2 ( ) ° <sup>Δ</sup> (in kJ mol-1) of interlayer water and Margules interaction parameters for the system H2O – smectite.

The integral hydration enthalpy is obtained by the following relationships:

$$
\Delta \tilde{H}\_{Hyd} = \mathbf{x}\_{hs} \, ^\ast \Delta \boldsymbol{H}\_{Hyd, 298}^\circ + \mathbf{x}\_{hs} \, ^\ast \left(1 - \mathbf{x}\_{hs}\right) ^\ast \left[\boldsymbol{W}\_{H\_1} \, ^\ast \mathbf{x}\_{hs} + \boldsymbol{W}\_{H\_2} \, ^\ast \left(1 - \mathbf{x}\_{hs}\right)\right] \tag{32}
$$

in which *HHyd*,298 ° <sup>Δ</sup> and Margules parameters *WH*<sup>1</sup> and *WH*<sup>2</sup> are determined by a minimisation procedure applied to the difference between computed and calorimetric integral enthalpies on eight homoionic montmorillonite SWy-1 [15-17] and given in table 7. The maximum number of moles of interlayer water was set to 5.5 moles. Figure 4A displays the comparison of integral enthalpy of hydration of a smectite- Na computed from the three models developed previously with different experimental values.

The curve 1 in Figure 4A (yellow dotted line) comes from model of Ransom & Helgeson [66], is linear and is obtained from the following equation of integral enthalpy of hydration:

$$
\Delta \tilde{H}\_{Hyd, 298} = -24.39 \, ^\circ \text{x}\_{hs} \tag{33}
$$

where the coefficient -24.39 is computed from *H HO <sup>f</sup>* ,298 2 *<sup>i</sup>* ( ) ° <sup>Δ</sup> given in Table 7. The curves 2b, 2c and 2d (green dotted lines) displayed in Figure 4A come together from the model of Vidal et al. [68] and are graphically represented by four truncated lines with their respective equations:

354 Thermodynamics – Fundamentals and Its Application in Science

*WH*<sup>1</sup>

*WH*<sup>2</sup>

1 –[66]; 2 – [68]; 3 - [63].

in which *HHyd*,298

The hydration enthalpy is retrieved from data of [24] for montmorillonite –Na and –K and derived from differential heats of adsorption obtained from the measurements of the heat of immersion [18, 20] for Na-, Ca- and Mg- montmorillonite. Thus, an integral enthalpy of hydration *H HO <sup>f</sup>* ,298 2 *<sup>i</sup>* ( ) ° <sup>Δ</sup> is provided for each of the four solid solutions with nm equal to 0.7, 2, 4 and 7 and is given in table 7. Vieillard et al. [63] consider the hydration of a smectite to

be an asymmetrical regular solid solution between anhydrous and hydrated smectite.

Units (kJ mol-1) H2O Cs Rb K Na Ba Sr Ca Mg *H HO <sup>f</sup>* ,298 2 *<sup>i</sup>* ( ) ° <sup>Δ</sup> 1 4.5 -289.93 -290.23 -290.45 -291.25 -294.14 -294.57 -294.86 -295.71

*H HO <sup>f</sup>* ,298 2 *<sup>i</sup>* ( ) ° <sup>Δ</sup> 2 0-0.7 -295.50 -299.71 -301.43 -302.86 WH 2 0. 0. 0. 0. *H HO <sup>f</sup>* ,298 2 *<sup>i</sup>* ( ) ° <sup>Δ</sup> <sup>2</sup> 0-2 -293.58 -295.90 -298.60 -300.05 WH 2 -10. -10. -10. -10. *H HO <sup>f</sup>* ,298 2 *<sup>i</sup>* ( ) ° <sup>Δ</sup> 2 0-4 -291.96 -293.83 -296.25 -297.62 WH 2 -10. -10. -10. -10. *H HO <sup>f</sup>* ,298 2 *<sup>i</sup>* ( ) ° <sup>Δ</sup> 2 0-7 -291.05 -292.33 -293.86 -294.69 WH 2 -10. -10. -10. -10. *H HO <sup>f</sup>* ,298 2 *<sup>i</sup>* ( ) ° <sup>Δ</sup> <sup>3</sup> 5.5 -287.47 -288.01 -288.16 -290.01 -292.10 -292.65 -294.72 -296.56

3 16.0 19.0 15.0 8.0 -60.0 -55.0 -60.0 -49.0

3 -90.0 -95.0 -68.0 -33.0 -75.0 -90.0 -100.0 -129.0

*Hyd hs Hyd hs hs H hs H* ( ) ( ) *H x H x x WxW xhs* 1 2 ,298 \* \* 1\*\* \* <sup>1</sup> ° Δ =Δ + − + − (32)

and *WH*<sup>2</sup>

*H x Hyd*,298 *hs* Δ =−24.39 \* (33)

are determined by a

**Table 7.** Standard thermodynamic enthalpy *H HO f i* ,298 2 ( ) ° <sup>Δ</sup> (in kJ mol-1) of interlayer water and

minimisation procedure applied to the difference between computed and calorimetric integral enthalpies on eight homoionic montmorillonite SWy-1 [15-17] and given in table 7. The maximum number of moles of interlayer water was set to 5.5 moles. Figure 4A displays the comparison of integral enthalpy of hydration of a smectite- Na computed from the three

The curve 1 in Figure 4A (yellow dotted line) comes from model of Ransom & Helgeson [66], is linear and is obtained from the following equation of integral enthalpy of hydration:

The integral hydration enthalpy is obtained by the following relationships:

° <sup>Δ</sup> and Margules parameters *WH*<sup>1</sup>

models developed previously with different experimental values.

Margules interaction parameters for the system H2O – smectite.

$$0 \text{ layer } \left( 0 < \mathbf{x}\_{\text{hs}} < 0.7 \right) \text{: } \Delta \tilde{H}\_{\text{hyd}, 298} = -9.716 \text{ }^{\circ}\mathbf{x}\_{\text{hs}} \tag{34}$$

$$1\text{ layer }\left(0\ltimes\_{\text{hs}}\langle 0\right) \colon \Delta \bar{H}\_{\text{Hyl},298} = -20.14 \, ^\circ \text{x}\_{\text{hs}} + \text{x}\_{\text{hs}} \, ^\circ \left(1 - \text{x}\_{\text{hs}}\right) \text{\*} 10.0\tag{35}$$

$$\text{2 layers } \left( 0 < \mathbf{x}\_{\text{hs}} < 4.0 \right) \text{: } \Delta \tilde{H}\_{\text{Hydro},298} = -\mathbf{3} \text{\textdegree } \mathbf{0} \text{: } \mathbf{x}\_{\text{hs}} + \mathbf{x}\_{\text{hs}} \text{\* } \left( 1 - \mathbf{x}\_{\text{hs}} \right) \text{\* } 10.0 \tag{36}$$

**Figure 4.** A – Integral enthalpy of hydration versus the number of interlayer cation; B – Standard enthalpy of formation of interlayer water versus the nature of the interlayer cation.

$$\text{GB layers } \left( 0 < \mathbf{x}\_{\text{hs}} < 7.0 \right) \text{: } \Delta \tilde{H}\_{Hyd, 298} = -45.5 \text{: } \mathbf{x}\_{\text{hs}} + \mathbf{x}\_{\text{hs}} \text{\* } \left( 1 - \mathbf{x}\_{\text{hs}} \right) \text{\* } 10.0 \tag{37}$$

The curve 3 (black line) in Figure 4A comes from the model of Vieillard et al. [63] and is represented by the following equation:

$$
\Delta \tilde{H}\_{Hyd} = -23.0 \, ^\circ \text{x}\_{hs} + \text{x}\_{hs} \, ^\circ \left(1 - \text{x}\_{hs}\right) \, ^\circ \left[8.0 \, ^\circ \text{x}\_{hs} - 33.0 \, ^\circ \left(1 - \text{x}\_{hs}\right)\right] \tag{38}
$$

It appears that integral enthalpy of hydration modelled by Vieillard et al. [63] (black line) based on experimental measurements of heat of adsorption from [15-17] encompasses nearly all experimental data. As the standard entropy of hydration is assumed constant in [66] and [68] models, the integral hydration enthalpies provided by these two previous models are nearly linear. Those obtained by [68] exhibit three sections of curves corresponding to the three states of hydration (1, 2 and 3 water layers).

In Figure 4B, a relationship between the standard state enthalpy of interlayer water with the nature of the interlayer cation located in the SWy-1 montmorillonite with a constant layer charge of 0.38, characterised by its electronegativity (ΔHO=Mz+ (aq), Table 4) is proposed [63] and given as follows:

$$
\Delta H\_{f,298}^{\circ}(H\_2O)\_i - \Delta H\_{f,298}^{\circ}(H\_2O)\_{\text{(l)}} = 0.02784^{\circ} \,\Delta\_HO^{\circ}M^{z+} - 6.275 \tag{39}
$$

Thermodynamics of Hydration in Minerals: How to Predict These Entities 357

obtained previously in the minimisation

1 *z*

*Hyd aq S M* ( ) ° + 2

(42)

represent the standard entropy of hydration and Margules

*hyd m* ( ) *<sup>i</sup>* ( ) *liq S n S HO S HO* 0 00 with ,298 298 2 ( ) 298 2 \* Δ= − (41)

parameters of the excess entropy of mixing, respectively. The integral entropy is expressed in J K-1 mol-1 of smectite based on O10(OH)2. The determination of these parameters is

procedure of integral hydration enthalpy and from the experimental adsorption-desorption isotherm interlayer water versus relative humidity. This minimisation procedure uses a least

 J/K/mole J/K/mole J/K/mole J/K/mole Cs+ 64.6 163.5 -348.0 -78. Rb+ 63.8 142.0 -355.0 -84. K+ 63.7 104.0 -267.5 -93. Na+ 59.8 49.1 -94.0 -130. Ba2+ 55.7 -116.0 -211.5 -224. Sr2+ 54.3 -117.8 -244.5 -261. Ca2+ 51.5 -93.5 -270.5 -271. Mg2+ 45.4 -42.8 -380.0 -350.

(Table 8).

<sup>1</sup>*WS*<sup>2</sup>

, and *WH*<sup>2</sup>

and *WS*<sup>2</sup>

**Table 8.** Standard thermodynamic enthalpy *Hyd <sup>S</sup>*° <sup>Δ</sup> of interlayer water, Margules interactions parameters and hydration entropy of the ion under its hydrated aqueous state for the eight cations.

equilibrium water – Swy1-Na (black line, N° 3), whose equation is given:

Figure 5A displays the comparison of the integral entropy of hydration of a smectite- Na computed from the three models [67] [68] and [63] developed previously with different

As indicated in the enthalpy section, integral hydration entropies modelled by [66] (yellow dotted line, N°1) and [68] (green dotted line, N°2a, 2b and 2c corresponding to 1st, 2nd and 3rd layer) are linear and merge together. In the model proposed by Vieillard et al. [63], three calculated integral entropies of SWy1-Na –water system were plotted in Figure 5A and correspond to the adsorption (red line), the desorption (blue line) and the theoretical

> ( ) ( ) *Hyd hs hs hs hs hs hs hs hs hs*

55.8 \* \* ln 1 \* ln 1

*S x Rx x x x xx x x*

\* 1 \* 49.05 \* 94.0 \* 1 Δ =− − +− − +− − −

The comparison of hydration-dehydration curves calculated by [63] with those provided by Fu et al. [34] show an opposite interpretation. The data from [34] show that integral entropies during adsorption are less negative than those during desorption. Data from our model show the opposite. This difference is explained by the fact that from experimental

( )( ) ( )

where *Hyd <sup>S</sup>*° <sup>Δ</sup> , *WS*<sup>1</sup>

1 –[63]; 2 - [73].

experimental values.

retrieved from values of *HHyd*

square method and provides *Hyd <sup>S</sup>*° <sup>Δ</sup> , *WS*<sup>1</sup>

*<sup>i</sup> S HO* 298 2 ( )

, and *WS*<sup>2</sup>

° <sup>Δ</sup> , *WH*<sup>1</sup>

° <sup>1</sup>*WS*<sup>1</sup>

Lowering the electronegativity of cation in the interlayer sites stabilises the enthalpy of formation of the interlayer water. It should be kept in mind that Eq. (39) has been settled for the same support, *i.e.,* a constant layer charge. For natural montmorillonite-Na, most enthalpies of hydration of the interlayer water are consistent within a narrow range of interlayer charge, IC (0.31<IC<0.38) [63]. Due to the lack of calorimetric measurements of enthalpy of hydration for high charge (IC>0.5) and low charge (IC<0.3) montmorillonites, the variation of the standard hydration enthalpy as a function of the layer charge remains highly questionable. An important remaining matter of discussion is the extent to which the model is able to predict hydration properties of dioctahedral and trioctahedral smectites with different layer charges and compositions.

#### **5.2. Entropy**

**Salts:** Average values of the entropy of hydration water were obtained [39] and [35] within chlorides, sulphates and salts (table 2). These average values of salts appear very similar to the value of entropy for crystalline water, S°=40.17 J/K-1 Mol-1 [70]. To improve the accuracy of prediction, statistical relationships between entropy and the molar volume of hydration should be investigated.

**Zeolites**: Initially, an average value of entropy of zeolitic water was proposed by Helgeson et al. [70] and was equal to S°H2O (zeol) = 58.99 J/K-1 Mol-1. From very few available calorimetric measurements of anhydrous and hydrated zeolites [4, 6, 71], it has been shown that the entropy of hydration water, at saturation, remains constant at approximately 52.0 J/K-1 Mol-1 [72]. Some authors provide slight deviations of entropy of hydration depending of the nature of cations in clinoptilolite [22].

**Clay minerals**: Entropies of interlayer water in clay minerals have never been measured directly, but have been evaluated from indirect measurements such as a contribution of experimental enthalpies of immersion and adsorption-desorption isotherms [63] [34] or from the Clausius-Clapeyron rules with adsorption isotherms performed at different temperatures [31], [30]. A constant value 55.0 J/K-1 Mol-1 for the entropy of interlayer water has been chosen [67] and [68] regardless of the nature of the interlayer water and the numbers of interlayer waters layers. By considering the hydration of a smectite always as a symmetrical regular solid solution, the integral entropy is expressed as in [63]:

$$\begin{split} \Delta \tilde{\mathbf{S}}\_{Hyd} &= \mathbf{x}\_{hs} \, ^\ast \Delta \mathbf{S}\_{Hyd} - \mathbb{R} \Big[ \mathbf{x}\_{hs} \, ^\ast \ln \Big( \mathbf{x}\_{hs} \Big) + \left( \mathbf{1} - \mathbf{x}\_{hs} \right) ^\ast \ln \Big( \mathbf{1} - \mathbf{x}\_{hs} \Big) \Big] \\ &+ \mathbf{x}\_{hs} \, ^\ast \Big( \mathbf{1} - \mathbf{x}\_{hs} \Big) ^\ast \Big[ \mathcal{W}\_{\mathbf{S}\_1} \, ^\ast \mathbf{x}\_{hs} + \mathcal{W}\_{\mathbf{S}\_2} \, ^\ast \Big( \mathbf{1} - \mathbf{x}\_{hs} \Big) \Big] \end{split} \tag{40}$$

Thermodynamics of Hydration in Minerals: How to Predict These Entities 357

$$\text{with}\quad \Delta S\_{hyd,298}^{0} = n\_m \, ^\ast \left[ S\_{298}^{0} \left( H\_2 O\_{\text{(i)}} \right) - S\_{298}^{0} \left( H\_2 O\_{\text{liq}} \right) \right] \tag{41}$$

where *Hyd <sup>S</sup>*° <sup>Δ</sup> , *WS*<sup>1</sup> , and *WS*<sup>2</sup> represent the standard entropy of hydration and Margules parameters of the excess entropy of mixing, respectively. The integral entropy is expressed in J K-1 mol-1 of smectite based on O10(OH)2. The determination of these parameters is retrieved from values of *HHyd* ° <sup>Δ</sup> , *WH*<sup>1</sup> , and *WH*<sup>2</sup> obtained previously in the minimisation procedure of integral hydration enthalpy and from the experimental adsorption-desorption isotherm interlayer water versus relative humidity. This minimisation procedure uses a least square method and provides *Hyd <sup>S</sup>*° <sup>Δ</sup> , *WS*<sup>1</sup> and *WS*<sup>2</sup> (Table 8).


1 –[63]; 2 - [73].

356 Thermodynamics – Fundamentals and Its Application in Science

with different layer charges and compositions.

of the nature of cations in clinoptilolite [22].

and given as follows:

**5.2. Entropy** 

in [63]:

should be investigated.

charge of 0.38, characterised by its electronegativity (ΔHO=Mz+ (aq), Table 4) is proposed [63]

Lowering the electronegativity of cation in the interlayer sites stabilises the enthalpy of formation of the interlayer water. It should be kept in mind that Eq. (39) has been settled for the same support, *i.e.,* a constant layer charge. For natural montmorillonite-Na, most enthalpies of hydration of the interlayer water are consistent within a narrow range of interlayer charge, IC (0.31<IC<0.38) [63]. Due to the lack of calorimetric measurements of enthalpy of hydration for high charge (IC>0.5) and low charge (IC<0.3) montmorillonites, the variation of the standard hydration enthalpy as a function of the layer charge remains highly questionable. An important remaining matter of discussion is the extent to which the model is able to predict hydration properties of dioctahedral and trioctahedral smectites

**Salts:** Average values of the entropy of hydration water were obtained [39] and [35] within chlorides, sulphates and salts (table 2). These average values of salts appear very similar to the value of entropy for crystalline water, S°=40.17 J/K-1 Mol-1 [70]. To improve the accuracy of prediction, statistical relationships between entropy and the molar volume of hydration

**Zeolites**: Initially, an average value of entropy of zeolitic water was proposed by Helgeson et al. [70] and was equal to S°H2O (zeol) = 58.99 J/K-1 Mol-1. From very few available calorimetric measurements of anhydrous and hydrated zeolites [4, 6, 71], it has been shown that the entropy of hydration water, at saturation, remains constant at approximately 52.0 J/K-1 Mol-1 [72]. Some authors provide slight deviations of entropy of hydration depending

**Clay minerals**: Entropies of interlayer water in clay minerals have never been measured directly, but have been evaluated from indirect measurements such as a contribution of experimental enthalpies of immersion and adsorption-desorption isotherms [63] [34] or from the Clausius-Clapeyron rules with adsorption isotherms performed at different temperatures [31], [30]. A constant value 55.0 J/K-1 Mol-1 for the entropy of interlayer water has been chosen [67] and [68] regardless of the nature of the interlayer water and the numbers of interlayer waters layers. By considering the hydration of a smectite always as a symmetrical regular solid solution, the integral entropy is expressed as

> ( ) ( ) *Hyd hs Hyd hs hs hs hs hs hs S hs S hs*

\* \* ln 1 \* ln 1

*S x S Rx x x x*

° Δ= Δ− +− − +− + −

*x x Wx W x* 1 2

\*1 \* \* \*1

( )( ) ( )

(40)

*<sup>z</sup> H HO H HO f if l* ,298 2 ( ) ( ) 0.02784 \* ,298 2 ( ) *<sup>H</sup>O M* 6.275 ° ° <sup>=</sup> <sup>+</sup> Δ − Δ = Δ− (39)

**Table 8.** Standard thermodynamic enthalpy *Hyd <sup>S</sup>*° <sup>Δ</sup> of interlayer water, Margules interactions parameters and hydration entropy of the ion under its hydrated aqueous state for the eight cations.

Figure 5A displays the comparison of the integral entropy of hydration of a smectite- Na computed from the three models [67] [68] and [63] developed previously with different experimental values.

As indicated in the enthalpy section, integral hydration entropies modelled by [66] (yellow dotted line, N°1) and [68] (green dotted line, N°2a, 2b and 2c corresponding to 1st, 2nd and 3rd layer) are linear and merge together. In the model proposed by Vieillard et al. [63], three calculated integral entropies of SWy1-Na –water system were plotted in Figure 5A and correspond to the adsorption (red line), the desorption (blue line) and the theoretical equilibrium water – Swy1-Na (black line, N° 3), whose equation is given:

$$\begin{aligned} \Delta \tilde{S}\_{Hyd} &= -55.8 \,\mathrm{\*}\,\mathrm{x}\_{hs} - R \Big[ \mathrm{x}\_{hs} \,\mathrm{\*}\ln\left(\mathrm{x}\_{hs}\right) + \left(1 - \mathrm{x}\_{hs}\right) \mathrm{\*}\ln\left(1 - \mathrm{x}\_{hs}\right) \Big] \\ &+ \mathrm{x}\_{hs} \,\mathrm{\*}\left(1 - \mathrm{x}\_{hs}\right) \mathrm{\*}\left[49.05 \,\mathrm{\*}\,\mathrm{x}\_{hs} - 94.0 \,\mathrm{\*}\,\mathrm{(1 - \mathrm{x}\_{hs})}\right] \end{aligned} \tag{42}$$

The comparison of hydration-dehydration curves calculated by [63] with those provided by Fu et al. [34] show an opposite interpretation. The data from [34] show that integral entropies during adsorption are less negative than those during desorption. Data from our model show the opposite. This difference is explained by the fact that from experimental

works of [34], adsorption follows desorption, while with our data, desorption follows adsorption. Thus, a maximum entropy difference between hydration and dehydration functions can be depicted and is equal to 18 J/K/mole. These observations show the importance of movement of exchangeable cations from ditrigonal cavities and the rotation of tetrahedrals in the tetrahedral sheets when dry collapsed layers are progressively exposed to water vapour.

Thermodynamics of Hydration in Minerals: How to Predict These Entities 359

<sup>2</sup> ( ) 56.61 2.64 \* 10 \* <sup>−</sup> = − [41]

<sup>2</sup> ( ) 87.62 7.58 \* 10 \* 0.53 \* 107 \* − − =+ + [74]

<sup>2</sup> ( ) 39.65 32.5 \* 10 \* 4.694 \* 10 \* − − =+ − [45]

2 . ( ) 23.999 66.295 \* 10 \* 0.3948 \* 10 \* − − =+ − [72]

Heat capacities of the hydration in zeolites have been performed by Neuhoff & Wang [75] on three zeolites (analcime, natrolite and wairakite) and exhibit marked variations in the heat capacity of hydration with temperature. Four zeolites (mordenite, wairakite and two different values from different sources for analcime and natrolite) for which the heat capacity of zeolitic water can be obtained by a difference between the heat capacities measured on both the anhydrous and the hydrated phases [72] are represented in Figure 6A. Also displayed are the heat capacity values of water in all its states (ice, liquid and vapour) for comparison [35]. The heat capacity values of zeolitic water obtained by minimisation techniques (Eq. D of table 9, turquish line) matches with the heat capacity function

**Figure 6.** Heat capacity of zeolitic water versus temperature: A – Heat capacities of zeolitic water in zeolites, in different prediction models and of water in its various states (l for liquid, c for ice and g for vapour); B – Experimental heat capacities of zeolitic waters in natural zeolites and predicted heat

Eq. D of table 9 provides a statistical error of 4.51%, 2.35% and 4.18% for hydrated zeolites at 298.15 K, 400 K and 500 K, respectively. To improve the accuracy of the prediction of the heat capacity of hydrated zeolites, one considers the normalised hydration volume, which is the difference of unit-cell volume (Vu.c.) between the hydrated and dehydrated forms for any

> ( ) *u c hydr u c dehyd Hy Deh H O V V*

.. . .. .

*n* 2

<sup>−</sup> Δ = (45)

 Heat capacity function (298.15- 500 K) Eq. A *Cp H O sil <sup>T</sup>* 2 0.5

**Table 9.** Heat capacity functions for zeolites minerals

equations of [74] and [45].

capacities from Eq. (46).

zeolite per zeolitic water molecule [72]:

*V*

−

Eq. B *Cp H O sil <sup>T</sup> <sup>T</sup>* 2 0.5 <sup>3</sup>

Eq. C *Cp H O sil <sup>T</sup> <sup>T</sup>* 3 5 <sup>2</sup>

Eq. D *Cp H O zeol <sup>T</sup> <sup>T</sup>* 3 5 <sup>2</sup>

**Figure 5.** A – Integral entropy of hydration versus the number of interlayer cation; B – Standard entropy of formation of interlayer water versus the nature of the interlayer cation.

A correlation between the entropy of interlayer cation *<sup>i</sup> S HO* 298 2 ( ) ° with the theoretical hydration entropy, S°hyd Mz+ [73] (Table 8), of the interlayer cation [63] is displayed in Figure 5B with the following equation:

$$\stackrel{\circ}{S}\_{298}^{\circ}(H\_2O)\_i - \stackrel{\circ}{S}\_{298}^{\circ}(H\_2O)\_{\text{(l)}} = 0.06542 \, ^\circ S \stackrel{\circ}{S}\_{hyd} M^{z+} - 0.413 \tag{43}$$

Negative values of the hydration entropy of the hydrated ion lower the hydration entropy of the interlayer water.

#### **5.3. Heat capacity**

**Salts:** An acceptable average heat capacity ranging from 39.0±5.3 J/K-1mol-1 for sulphates and sulphites to 42.0±2.7 J/K-1mol-1for chlorates and chlorides has been proposed [35] and seems to be close to the average value of heat capacity for crystalline water Cp=40.04 J/K-1mol-1 [70]. Only the average heat capacity function of crystalline water was proposed [70]:

$$\Pr(H\_2O)\_{\text{salt}} = 29.75 - 3.448 \, ^\circ 10^{-2} \, ^\circ T \tag{44}$$

**Zeolites:** Three equations of heat capacities ([41], [74] [45] ) were proposed for a set of silicate minerals, including analcime and natrolite, and are displayed in table 9.


**Table 9.** Heat capacity functions for zeolites minerals

to water vapour.

works of [34], adsorption follows desorption, while with our data, desorption follows adsorption. Thus, a maximum entropy difference between hydration and dehydration functions can be depicted and is equal to 18 J/K/mole. These observations show the importance of movement of exchangeable cations from ditrigonal cavities and the rotation of tetrahedrals in the tetrahedral sheets when dry collapsed layers are progressively exposed

**Figure 5.** A – Integral entropy of hydration versus the number of interlayer cation; B – Standard

A correlation between the entropy of interlayer cation *<sup>i</sup> S HO* 298 2 ( ) ° with the theoretical hydration entropy, S°hyd Mz+ [73] (Table 8), of the interlayer cation [63] is displayed in

Negative values of the hydration entropy of the hydrated ion lower the hydration entropy

**Salts:** An acceptable average heat capacity ranging from 39.0±5.3 J/K-1mol-1 for sulphates and sulphites to 42.0±2.7 J/K-1mol-1for chlorates and chlorides has been proposed [35] and seems to be close to the average value of heat capacity for crystalline water Cp=40.04 J/K-1mol-1 [70].

*Cp H O salts <sup>T</sup>* <sup>2</sup>

**Zeolites:** Three equations of heat capacities ([41], [74] [45] ) were proposed for a set of

Only the average heat capacity function of crystalline water was proposed [70]:

silicate minerals, including analcime and natrolite, and are displayed in table 9.

*z*

<sup>2</sup> ( ) 29.75 3.448 \* 10 \* <sup>−</sup> = − (44)

*<sup>i</sup> <sup>l</sup> hyd S H O S H O 6542 S M 413* 298 2 298 2 ( ) ( ) ( ) 0.0 \* 0. ° ° ° + −= − (43)

entropy of formation of interlayer water versus the nature of the interlayer cation.

Figure 5B with the following equation:

of the interlayer water.

**5.3. Heat capacity** 

Heat capacities of the hydration in zeolites have been performed by Neuhoff & Wang [75] on three zeolites (analcime, natrolite and wairakite) and exhibit marked variations in the heat capacity of hydration with temperature. Four zeolites (mordenite, wairakite and two different values from different sources for analcime and natrolite) for which the heat capacity of zeolitic water can be obtained by a difference between the heat capacities measured on both the anhydrous and the hydrated phases [72] are represented in Figure 6A. Also displayed are the heat capacity values of water in all its states (ice, liquid and vapour) for comparison [35]. The heat capacity values of zeolitic water obtained by minimisation techniques (Eq. D of table 9, turquish line) matches with the heat capacity function equations of [74] and [45].

**Figure 6.** Heat capacity of zeolitic water versus temperature: A – Heat capacities of zeolitic water in zeolites, in different prediction models and of water in its various states (l for liquid, c for ice and g for vapour); B – Experimental heat capacities of zeolitic waters in natural zeolites and predicted heat capacities from Eq. (46).

Eq. D of table 9 provides a statistical error of 4.51%, 2.35% and 4.18% for hydrated zeolites at 298.15 K, 400 K and 500 K, respectively. To improve the accuracy of the prediction of the heat capacity of hydrated zeolites, one considers the normalised hydration volume, which is the difference of unit-cell volume (Vu.c.) between the hydrated and dehydrated forms for any zeolite per zeolitic water molecule [72]:

$$
\Delta V\_{Hy-Deh} = \left(V\_{\mu.c.hydr.} - V\_{\mu.c.delyd.}\right) \Bigg/ \tag{45}
$$

$$\begin{aligned} \text{Cp}(H\_2O)\_{z\text{sol.}} &= \left[ 55.85 - 13.073 \, ^\circ \sqrt{\left( \Delta V\_{Hy-D\text{ch}} \right)} \right] \cdot \left[ 0.0062 + 0.0305 \, ^\circ \sqrt{\left( \Delta V\_{Hy-D\text{ch}} \right)} \right] \cdot ^\ast T \\ &+ \left[ 0.008 + 0.0048 \, ^\circ \sqrt{\left( \Delta V\_{Hy-D\text{ch}} \right)} \right] \cdot ^\ast T^{-2} \end{aligned} \tag{46}$$

Thermodynamics of Hydration in Minerals: How to Predict These Entities 361

Δ =Δ − Δ G H T\* S (48)

may be determined from the gravimetric

<sup>=</sup> (50)

and for the partial quantities that

**6. Discussion of the Gibbs free energy of the formation behaviour of the** 

The relationships among the free energy change, enthalpy of reaction and entropy for an

where ΔG stands for free energy change, ΔH for enthalpy of reaction, ΔS for the change in entropy and T for the absolute temperature. This relationship holds true for both the

**Salts**: The free energy of formation of the hydration water of hydrated salts can be calculated from the enthalpy and entropy of hydration given previously. A relationship has been demonstrated by the statistical regressions [35] between the enthalpy and Gibbs free

The predicted values contribute to the variations of solubility products and the cation selectivity constants among the series of salts. An example is given in [78] for two series of

For zeolites and clays, the adsorbent (a single smectite or zeolite phase) is supposed to be inert so that the derived enthalpy and entropy functions formally represent the energy

*<sup>w</sup> n n*

=

where R = 8.314 J.mol-1K-1 is the ideal gas constant, n is the amount of adsorbed water, P is the partial water vapour pressure, and P0 is the saturated water vapour pressure at temperature T (P/P0 = Relative Humidity/100). This value can be determined by plotting R\*T\*Ln(P/Po) versus the water content n (adsorbed or desorbed) and determining the area under the curve from 0 to nw. When nw approaches 0, R\*T\*Ln(P/Po) approaches infinity,

**Zeolites:** Eq. (50) has been used on three clinoptilolites (-Na, -K and -Ca) [22] and show a good correlation between the Gibbs free energies calculated from the adsorption-desorption isotherms and from the contribution of the enthalpy of hydration and the entropy of hydration extracted from the equilibrium at high temperature. The model described by Carey & Bish [62] used a solution solid model between anhydrous and hydrated zeolite, and the variation of the

=

*<sup>P</sup> <sup>Δ</sup>G R\*T\* Ln \*dn*

*n 0 0*

*P*

isotherms of adsorption of the water vapour by the following expression [79]:

*hyd*

() () *<sup>Δ</sup>G HO \* f s* ,298 2 *alts fs <sup>Δ</sup>H HO* ,298 2 *alts* 1.008 59.406 ° ° = + (49)

**hydration water in environmental conditions** 

integral entities that are represented by Δ*G* , Δ*H* and Δ*S*

energy of crystalline water in the following relation:

variations between the adsorbed and the bulk water.

The integral free energy of hydration *ΔGhyd*

making integration somewhat difficult.

are represented by Δ*G* , Δ*H* and Δ*S* .

hydrated salts of MgCl2 and SrCl2.

isothermal process is:

The heat capacities of hydration water given in Figure 6B have been calculated for the four hydrated zeolites (mordenite, analcime, wairakite and natrolite) from the unit-cell volumes of their dehydrated and hydrated forms [72]. It can be observed that the heat capacities of hydration water calculated in this way are much better than those obtained by Eq. D of table 9 (indicated in Fig. 6B as a turquish full line). With this improvement, the errors made on hydrated zeolites decrease to 3.5%, 2.35% and 3.89% at 298.15 K, 400 K and 500 K, respectively.

**Clay minerals:** The heat capacity function comes initially from [67] and is assumed to be independent of the nature of the interlayer cations:

$$\text{Cp}(H\_2O)\_{\text{clays}} = 37.84 + 51.631 \, ^\circ \, 10^{-3} \, ^\circ T - 4.0959 \, ^\circ \, 10^5 \, ^\circ T^{-2} \tag{47}$$

This heat capacity equation has also been used by [68] and by [63] in the behaviour of thermodynamic entities versus temperature.

#### **5.4. Molar volume**

**Salts:** Only [35] provided an average value for different salts: 14.5±1.2 cm3/mole for sulphates and sulphites and 14.1±2.3 cm3/mole for chlorides and chlorates. These average values appear to be similar to the molar volume of crystalline water (V(H2O)(salts)= 13.7 cm3/mole [76]).

**Zeolites:** No study on the variation of the molar volumes of the zeolitic water has been performed. Numerous values of molar volumes of hydrated and anhydrous zeolites have been provided and gathered in such calculations for enthalpy and entropy of hydration [77], [61], [72]. In the study of the relationship between entropy and molar volume within zeolites [72], an average value of the molar volume of hydration water has been calculated and is V(H2O)(zeol)=2.06 cm3/mole, which is not significantly different from those proposed by Helgeson (V(H2O)(zeol)= 8 cm3/mole) estimated from standard molal volumes of zeolites. The weak value of the molar volume of the hydration water is explained by the fact that the anhydrous zeolites exhibit a rigid (Si, Al) framework with different channels or cavities occupied by exchangeable cations. During the hydration process, the water molecules are bound to cations inside the cavities, do not affect the framework and have a weak impact on lattice parameters.

**Clay minerals:** The molar volume of the interlayer water was initially determined by Ransom [67] and also used by [68] and by [63] in thermochemical calculation. Its value is V (H2O)I = 17.22 cm3/mole. This value, greater than those of salts, is explained by the fact that, during hydration, the molar volume of clay minerals increases in one dimension characterised by basal spacing (Miller index 001).

## **6. Discussion of the Gibbs free energy of the formation behaviour of the hydration water in environmental conditions**

360 Thermodynamics – Fundamentals and Its Application in Science

independent of the nature of the interlayer cations:

thermodynamic entities versus temperature.

characterised by basal spacing (Miller index 001).

2 .

respectively.

**5.4. Molar volume** 

lattice parameters.

( ) ( )

− −

(46)

2

<sup>2</sup> ( ) 37.84 51.631 \* 10 \* 4.0959 \* 10 \* − − =+ − (47)

−

( )

0.008 0.0048 \* \*

 ++ Δ 

*Hy Deh*

−

*zeol Hy Deh Hy Deh*

=− Δ + + Δ

*V T*

The heat capacities of hydration water given in Figure 6B have been calculated for the four hydrated zeolites (mordenite, analcime, wairakite and natrolite) from the unit-cell volumes of their dehydrated and hydrated forms [72]. It can be observed that the heat capacities of hydration water calculated in this way are much better than those obtained by Eq. D of table 9 (indicated in Fig. 6B as a turquish full line). With this improvement, the errors made on hydrated zeolites decrease to 3.5%, 2.35% and 3.89% at 298.15 K, 400 K and 500 K,

**Clay minerals:** The heat capacity function comes initially from [67] and is assumed to be

*Cp H O clays <sup>T</sup> <sup>T</sup>* 3 5 <sup>2</sup>

This heat capacity equation has also been used by [68] and by [63] in the behaviour of

**Salts:** Only [35] provided an average value for different salts: 14.5±1.2 cm3/mole for sulphates and sulphites and 14.1±2.3 cm3/mole for chlorides and chlorates. These average values appear

**Zeolites:** No study on the variation of the molar volumes of the zeolitic water has been performed. Numerous values of molar volumes of hydrated and anhydrous zeolites have been provided and gathered in such calculations for enthalpy and entropy of hydration [77], [61], [72]. In the study of the relationship between entropy and molar volume within zeolites [72], an average value of the molar volume of hydration water has been calculated and is V(H2O)(zeol)=2.06 cm3/mole, which is not significantly different from those proposed by Helgeson (V(H2O)(zeol)= 8 cm3/mole) estimated from standard molal volumes of zeolites. The weak value of the molar volume of the hydration water is explained by the fact that the anhydrous zeolites exhibit a rigid (Si, Al) framework with different channels or cavities occupied by exchangeable cations. During the hydration process, the water molecules are bound to cations inside the cavities, do not affect the framework and have a weak impact on

**Clay minerals:** The molar volume of the interlayer water was initially determined by Ransom [67] and also used by [68] and by [63] in thermochemical calculation. Its value is V (H2O)I = 17.22 cm3/mole. This value, greater than those of salts, is explained by the fact that, during hydration, the molar volume of clay minerals increases in one dimension

to be similar to the molar volume of crystalline water (V(H2O)(salts)= 13.7 cm3/mole [76]).

*Cp H O V V T*

( ) 55.85 13.073 \* 0.0062 0.0305 \* \*

The relationships among the free energy change, enthalpy of reaction and entropy for an isothermal process is:

$$
\Delta \mathbf{G} = \Delta \mathbf{H} - \mathbf{T}^\* \Delta \mathbf{S} \tag{48}
$$

where ΔG stands for free energy change, ΔH for enthalpy of reaction, ΔS for the change in entropy and T for the absolute temperature. This relationship holds true for both the integral entities that are represented by Δ*G* , Δ*H* and Δ*S* and for the partial quantities that are represented by Δ*G* , Δ*H* and Δ*S* .

**Salts**: The free energy of formation of the hydration water of hydrated salts can be calculated from the enthalpy and entropy of hydration given previously. A relationship has been demonstrated by the statistical regressions [35] between the enthalpy and Gibbs free energy of crystalline water in the following relation:

$$
\Delta G\_{f,298}^{\circ} \left( H\_2O\_{\text{salt}} \right) = 1.008^\* \Delta H\_{f,298}^{\circ} \left( H\_2O\_{\text{salt}} \right) + 59.406 \tag{49}
$$

The predicted values contribute to the variations of solubility products and the cation selectivity constants among the series of salts. An example is given in [78] for two series of hydrated salts of MgCl2 and SrCl2.

For zeolites and clays, the adsorbent (a single smectite or zeolite phase) is supposed to be inert so that the derived enthalpy and entropy functions formally represent the energy variations between the adsorbed and the bulk water.

The integral free energy of hydration *ΔGhyd* may be determined from the gravimetric isotherms of adsorption of the water vapour by the following expression [79]:

$$
\Delta \tilde{G}\_{hyd} = R'' T^\* \int\_{n=0}^{n=n\_w} L n \left( \frac{P}{P\_0} \right) ^\* d\eta \tag{50}
$$

where R = 8.314 J.mol-1K-1 is the ideal gas constant, n is the amount of adsorbed water, P is the partial water vapour pressure, and P0 is the saturated water vapour pressure at temperature T (P/P0 = Relative Humidity/100). This value can be determined by plotting R\*T\*Ln(P/Po) versus the water content n (adsorbed or desorbed) and determining the area under the curve from 0 to nw. When nw approaches 0, R\*T\*Ln(P/Po) approaches infinity, making integration somewhat difficult.

**Zeolites:** Eq. (50) has been used on three clinoptilolites (-Na, -K and -Ca) [22] and show a good correlation between the Gibbs free energies calculated from the adsorption-desorption isotherms and from the contribution of the enthalpy of hydration and the entropy of hydration extracted from the equilibrium at high temperature. The model described by Carey & Bish [62] used a solution solid model between anhydrous and hydrated zeolite, and the variation of the integral Gibbs free energies of hydration is expressed as a function of a parameter θ characterising the ratio of H2O/(maximum H2O). The main difficulty in these calculations is the determination of the maximal amount of hydration water at the saturation state (p/p°=1).

Thermodynamics of Hydration in Minerals: How to Predict These Entities 363

In each figure, the data of the free energy of hydration (Eq. 50) from the experimental isotherms [15, 16] (blue square for adsorption and red square for desorption) and the free energy of hydration are computed and are reported for comparison of three different works: the Vieillard et al. model [63] (blue, red and black lines for adsorption, desorption and theoretical equilibrium, respectively), the Ransom & Helgeson model [66] (dotted yellow line), and the Vidal & Dubacq model [68] (dotted green line). For the adsorption and desorption reactions, the modelled isotherm calculated by the Vieillard et al. model [63] (blue and red lines, respectively) agrees well with the experimental data acquired by [15, 16]

At low relative humidity, the uncertainty appears very important between the experimental desorption data and the desorption curve from the model in [63] (3 kJ/mol for integral free

The theoretical equilibrium (black line) suggests that the hydration of a Na-SWy-1 begins at only R.H.= 0.15 (Figure 7A), which correlates well with the experimental adsorption but not with the desorption. The integral Gibbs free energy of hydration (figure 7B and 7C) and the Gibbs free energy of the interlayer water (Figure 7D) of the theoretical equilibrium between the water and Na-Smectite (black line) are plotted and correspond to the average values of calculated adsorption-desorption isotherms. The theoretical curve of equilibrium between the water and Na-smectite provided by Ransom & Helgeson [66] displays a small number of moles of the interlayer water close to saturation, while those of Vidal & Dubacq [68] reproduce this stepwise behaviour fairly well because it assumes three solid solutions with 3 end-members (0.7 H2O, 2 H2O, 4.5 H2O). However, it does not match the behaviour of water adsorption versus relative humidity, which is not a step function. The advantage of this approach is that it allows reproducing the observed stepwise evolution of the smectite volume and water content with varying T and relative humidity. At low relative humidity, as the hydration of a Na-SWy-1 begins at only R.H.= 0.15, the computations from [68] and [66] (green and yellow dotted lines) do not match the observed isotherms because there is a strong hydration in the RH= 0.0- 0.18 range. During the interval of R.H.=0.0-0.2, the adsorption of water molecules occurs predominantly on the external surface of tactoids (basal external faces and edges of tactoids). This explains why we observe a very negative integral free energy in the low relative humidity domain for these two models (green dotted line and in a less extent, yellow dotted lines). The necessity to take into account only the hydration water, sensu-stricto, is justified. This procedure has been incompletely performed by Ransom & Helgeson [66] with the external surface of a kaolinite to estimate the amount of surface covering water, while Vidal & Dubacq [68] made no correction on the gravimetric waters. At high relative humidity, Vidal & Dubacq [68] proposed a third step of hydration, with the number of water molecules increasing up to 7 moles per O10(OH)2. All of the experimental isotherms performed by [15, 16] exhibit an increase in the water contents beyond RH = 0.9, identified as osmotic water. The osmotic swelling of the sodium montmorillonite corresponds to an

energy, Fig C) due to the analytical limits and mathematical constraints.

in high relative humidity.

**Clay minerals:** Among the numerous measurements of adsorption-desorption isotherms on various clays, very few works provided a consistent Gibbs free energy of adsorptiondesorption due to the difficulty of measurements in the low relative humidity domains. By choosing the system of water – SWy-1-Na, Figure 7 displays the adsorption-desorption isotherms (Figure 7A), the variation of the integral free energy of hydration versus the number of moles of the interlayer water (Figure 7B) and versus the relative humidity (Figure 7C) and the variation of the free energy of formation of the interlayer water versus the relative humidity (Figure 7D).

**Figure 7.** A - Adsorption-desorption isotherms; B - Integral free energy of hydration versus the number of moles of interlayer water; C - Integral free energy of hydration versus the relative humidity; D - Free energy of formation of the interlayer water versus the relative humidity.

In each figure, the data of the free energy of hydration (Eq. 50) from the experimental isotherms [15, 16] (blue square for adsorption and red square for desorption) and the free energy of hydration are computed and are reported for comparison of three different works: the Vieillard et al. model [63] (blue, red and black lines for adsorption, desorption and theoretical equilibrium, respectively), the Ransom & Helgeson model [66] (dotted yellow line), and the Vidal & Dubacq model [68] (dotted green line). For the adsorption and desorption reactions, the modelled isotherm calculated by the Vieillard et al. model [63] (blue and red lines, respectively) agrees well with the experimental data acquired by [15, 16] in high relative humidity.

362 Thermodynamics – Fundamentals and Its Application in Science

relative humidity (Figure 7D).

integral Gibbs free energies of hydration is expressed as a function of a parameter θ characterising the ratio of H2O/(maximum H2O). The main difficulty in these calculations is the determination of the maximal amount of hydration water at the saturation state (p/p°=1).

**Clay minerals:** Among the numerous measurements of adsorption-desorption isotherms on various clays, very few works provided a consistent Gibbs free energy of adsorptiondesorption due to the difficulty of measurements in the low relative humidity domains. By choosing the system of water – SWy-1-Na, Figure 7 displays the adsorption-desorption isotherms (Figure 7A), the variation of the integral free energy of hydration versus the number of moles of the interlayer water (Figure 7B) and versus the relative humidity (Figure 7C) and the variation of the free energy of formation of the interlayer water versus the

**Figure 7.** A - Adsorption-desorption isotherms; B - Integral free energy of hydration versus the number of moles of interlayer water; C - Integral free energy of hydration versus the relative humidity; D - Free

energy of formation of the interlayer water versus the relative humidity.

At low relative humidity, the uncertainty appears very important between the experimental desorption data and the desorption curve from the model in [63] (3 kJ/mol for integral free energy, Fig C) due to the analytical limits and mathematical constraints.

The theoretical equilibrium (black line) suggests that the hydration of a Na-SWy-1 begins at only R.H.= 0.15 (Figure 7A), which correlates well with the experimental adsorption but not with the desorption. The integral Gibbs free energy of hydration (figure 7B and 7C) and the Gibbs free energy of the interlayer water (Figure 7D) of the theoretical equilibrium between the water and Na-Smectite (black line) are plotted and correspond to the average values of calculated adsorption-desorption isotherms. The theoretical curve of equilibrium between the water and Na-smectite provided by Ransom & Helgeson [66] displays a small number of moles of the interlayer water close to saturation, while those of Vidal & Dubacq [68] reproduce this stepwise behaviour fairly well because it assumes three solid solutions with 3 end-members (0.7 H2O, 2 H2O, 4.5 H2O). However, it does not match the behaviour of water adsorption versus relative humidity, which is not a step function. The advantage of this approach is that it allows reproducing the observed stepwise evolution of the smectite volume and water content with varying T and relative humidity. At low relative humidity, as the hydration of a Na-SWy-1 begins at only R.H.= 0.15, the computations from [68] and [66] (green and yellow dotted lines) do not match the observed isotherms because there is a strong hydration in the RH= 0.0- 0.18 range. During the interval of R.H.=0.0-0.2, the adsorption of water molecules occurs predominantly on the external surface of tactoids (basal external faces and edges of tactoids). This explains why we observe a very negative integral free energy in the low relative humidity domain for these two models (green dotted line and in a less extent, yellow dotted lines). The necessity to take into account only the hydration water, sensu-stricto, is justified. This procedure has been incompletely performed by Ransom & Helgeson [66] with the external surface of a kaolinite to estimate the amount of surface covering water, while Vidal & Dubacq [68] made no correction on the gravimetric waters. At high relative humidity, Vidal & Dubacq [68] proposed a third step of hydration, with the number of water molecules increasing up to 7 moles per O10(OH)2. All of the experimental isotherms performed by [15, 16] exhibit an increase in the water contents beyond RH = 0.9, identified as osmotic water. The osmotic swelling of the sodium montmorillonite corresponds to an

iso-enthalpic effect, with the heat of adsorption being constant. The theoretical equilibrium between water and Na-SWy-1 (black line) matches the experimental dehydration near saturation.

Thermodynamics of Hydration in Minerals: How to Predict These Entities 365

council for scientific Research (CNRS) is gratefully acknowledged. Special thanks to Ph. Blanc, H. Gailhanou and Cl. Fialips for their effective and decisive contribution to this

[1] Reeuwijk van LP (1974) The thermal dehydration of natural zeolites. In: H. V, B.V. Z,

[2] Tardy Y, Mercury L, Roquin C, Vieillard P (1999) Le concept d'eau *ice-like*: hydratationdeshydratation des sels, hydroxydes, zéolites, argiles et matières organiques vivantes

[3] Mercury L, Vieillard P, Tardy Y (2001) "Thermodynamics of Icepolymorphs and "Ice-

[4] Johnson GK, Flotow HE, O'Hare PAG, Wise WS (1982) Thermodynamic studies of

[5] Johnson GK, Flotow HE, O'Hare PAG (1983) Thermodynamic studies of

[6] Johnson GK, Tasker IR, Flotow HE, O'Hare PAG (1992) Thermodynamic studies of

[7] Gailhanou H, Blanc P, Rogez J, Mikaelian G, Kawaji H, Olives J, et al. (2012) Thermodynamic properties of illite IMt-2, smectite MX-80 and beidellite SBld-1 by calorimetric methods: Enthalpies of formation, heat capacities , entropies and Gibbs free

[8] Kiseleva I, Navrotsky A, Belitskii IA, Fursenko BA (1996) Thermochemistry and phase

[9] Sun P, Deore S, Navrotsky A (2007) Enthalpy of formation and dehydration of lithium

[10] Sun P, Navrotsky A (2008) Enthalpy of formation and dehydration of alkaline earth

[11] Yang SY, Navrotsky A (2000) Energetics of formation and hydration of ion-exchanged

[12] Shim SH, Navrotsky A, Gaffney TR, Macdougall JE (1999) Chabazite: Energetics of hydration, enthalpy of formation, and effect of cations on stability. Amer. Mineral. 84:

[13] Gailhanou H, van Miltenburg JC, Rogez J, Olives J, Amouric M, Gaucher EC, et al. (2007) Thermodynamic properties of anhydrous smectite MX-80, illite IMt-2 and mixed-layer illite-smectite ISCz-1 as determined by calorimetric methods. Part I: Heat capacities, heat contents and entropies. Geochim. Cosmochim. Acta 71:

cation exchanged zeolite beta. Microporous Mesoporous Mat. 109: 147-155.

ou inertes. Compt. Rendus Académie Science Paris 329: 377-388.

like" water in hydrates and hydroxides. Appl. Geochem. 16: 161-181.

zeolites : analcime and dehydrated analcime. Amer. Mineral. 67: 736-748.

zeolites:natrolite mesolite and scolecite. Amer. Mineral. 68: 1137-1145.

mordenite, dehydrated mordenite and gibbsite. Am. Mineral. 77: 85-93.

energies of formation. . Am. Mineral. (submitted).

equilibria in calcium zeolites. Am. Mineral. 81: 658-667.

zeolite Y. Microporous Mesoporous Materials 37: 175-186.

and sodium zeolite beta. Microporous Mesoporous Mat. 98: 29-40.

document.

**8. References** 

1870-1882.

5463-5473.

editors. Wageningen 1-88.

The cationic exchange and hydration reactions are intimately related, and the cationic exchange reactions not only modify the interlayer composition but also imply changes in the interlayer water content. These changes can be quite significant and can be evaluated using the Vieillard et al. model [63], which provides a theoretical number of moles of interlayer water and integral thermodynamic properties (H, S and G) for the SWy-1 montmorillonites. Thus, the cationic exchange reactions may depend on an external parameter, such as the relative humidity.

## **7. Conclusions and perspectives**

This chapter demonstrated some useful relationships in the prediction of thermodynamic entities. Some important areas are still lacking and deserve further study.

Among the anhydrous and hydrated salts, the relationships between the entropy and molar volume may exist and will improve the prediction of the entropy of hydration. The function of the average heat capacities for the hydration water is lacking and will help in characterising the behaviour of the crystalline water versus temperature.

In the zeolite field, the relationships among the molar volume, water content and relative humidity may exist within zeolites families with a reversibility in the hydration-dehydration processes. These relationships will contribute to a better understanding of the behaviour of partial molar enthalpy versus the degree of saturation or versus the relative humidity.

In the clay Minerals, three following recommendations are needed: 1) to dissociate the hydration water from the gravimetric water; 2) to measure the heat of hydration among the dioctahedral and trioctahedral clays; and 3) to find a relationship between the hysteresis and the entropy of the hydration-dehydration by the molecular modellings. With these three aims, the knowledge of mechanisms of the hydration within clay minerals will advance our understanding of cationic exchange in soils under various physical and chemical constraints.

## **Author details**

Philippe Vieillard *CNRS IC2MP, UMR 7285, Hydrasa, France* 

## **Acknowledgement**

Financial support from the French National Radioactive Waste Management Agency (ANDRA), from the French Geological Survey (BRGM) and from the national council for scientific Research (CNRS) is gratefully acknowledged. Special thanks to Ph. Blanc, H. Gailhanou and Cl. Fialips for their effective and decisive contribution to this document.

## **8. References**

364 Thermodynamics – Fundamentals and Its Application in Science

dehydration near saturation.

**7. Conclusions and perspectives** 

relative humidity.

constraints.

**Author details** 

Philippe Vieillard

**Acknowledgement** 

*CNRS IC2MP, UMR 7285, Hydrasa, France* 

iso-enthalpic effect, with the heat of adsorption being constant. The theoretical equilibrium between water and Na-SWy-1 (black line) matches the experimental

The cationic exchange and hydration reactions are intimately related, and the cationic exchange reactions not only modify the interlayer composition but also imply changes in the interlayer water content. These changes can be quite significant and can be evaluated using the Vieillard et al. model [63], which provides a theoretical number of moles of interlayer water and integral thermodynamic properties (H, S and G) for the SWy-1 montmorillonites. Thus, the cationic exchange reactions may depend on an external parameter, such as the

This chapter demonstrated some useful relationships in the prediction of thermodynamic

Among the anhydrous and hydrated salts, the relationships between the entropy and molar volume may exist and will improve the prediction of the entropy of hydration. The function of the average heat capacities for the hydration water is lacking and will help in

In the zeolite field, the relationships among the molar volume, water content and relative humidity may exist within zeolites families with a reversibility in the hydration-dehydration processes. These relationships will contribute to a better understanding of the behaviour of partial molar enthalpy versus the degree of saturation or versus the relative humidity.

In the clay Minerals, three following recommendations are needed: 1) to dissociate the hydration water from the gravimetric water; 2) to measure the heat of hydration among the dioctahedral and trioctahedral clays; and 3) to find a relationship between the hysteresis and the entropy of the hydration-dehydration by the molecular modellings. With these three aims, the knowledge of mechanisms of the hydration within clay minerals will advance our understanding of cationic exchange in soils under various physical and chemical

Financial support from the French National Radioactive Waste Management Agency (ANDRA), from the French Geological Survey (BRGM) and from the national

entities. Some important areas are still lacking and deserve further study.

characterising the behaviour of the crystalline water versus temperature.


[14] Barrer RM, Cram PJ (1971) Heats of immersion of outgassed and ion-exchanged zeolites. In: L.B. FEMS, editor. Molecular sieves zeolites. Washington D.C.: American Chemical Society. pp. 105-131.

Thermodynamics of Hydration in Minerals: How to Predict These Entities 367

[30] Keren R, Shainberg I (1980) Water vapor isotherms and heat of immersion of Na/Ca-Montmorillonite systems - 3. Thermodynamics. Clays Clay Minerals 28: 204-210. [31] Hall PL, Astill DM (1989) Adsorption of water by homoionic exchange forms of

[32] Valueva GP, Goryainov SV (1992) Chabazite during dehydration: thermochemical and

[33] Xu W, Johnston CT, Parker P, Agnew SF (2000) Infrared study of water sorption on Na-, Li-, Ca- and Mg-exchanged (SWy-1 AND SAz-1) montmorillonite. Clays Clay Minerals

[34] Fu MH, Zhang ZZ, Low PF (1990) Changes in the properties of a montmorillonite-water system during the adsorption and desorption of water: hysteresis. Clays Clay Minerals

[35] Mercury L, Vieillard P, Tardy Y (2001) "Thermodynamics of Icepolymorphs and "Ice-

[36] Cox JD, Wagman DD, Medvedev VA (1989) CODATA Key values for thermodynamics.

[37] Tardy Y, Mercury L, Roquin C, Vieillard P (1999) Le concept d'eau *ice-like*: hydratationdeshydratation des sels, hydroxydes, zéolites, argiles et matières organiques vivantes

[38] Chipera SJ, Apps JA (2001) Geochemical stability of natural zeolites. In: Bish DLM, D.W., editor. Natural Zeolites:Occurence, properties Applications. Blacksburg (Va)117-

[39] Glasser L, Jenkins HDB (2007) The thermodynamic solvate difference rule: Solvation parameters and their use in interpretation of the role of bound solvent in condensed-

[41] Robinson GR, Haas JL (1983) Heat capacity, relative enthalpy and calorimetric entropy of silicate minerals : an empirical method of prediction. Amer. Mineral. 68: 541-

[42] Moiseev GK, Sestak J (1995) Some calculations methods for estimation of thermodynamical and thermochemical properties of inorganic compounds. Prog Cryst

[43] Burylev BP (1988) Thermodynamics of crystal hydrates of rare-earth metal chlorides.

[44] Burylev BP, Kritskaya EB, Moisov LP, Kostenko NB (1995) Calculation of the standard

[45] Van Hinsberg VJ, Vriend SP, Schumacher JC (2005) A new method to calculate end-member thermodynamic properties of minerals from their constituent polyhedra II: heat capacity, compressibility and thermal expansion. J. Metamorph Geol. 23: 681-

enthalpies of formation of hydrous sodium salts. Inorg. Mater. 31: 727-729.

[40] Sidorov YI (1988) Estimating silicate standard entropies. Geochem. Intern. 80-85.

Wyoming montmorillonite (SWy-1). Clays Clay Minerals 37: 355-363.

Raman spectroscopy study. Russian Geologica Geophys. 33: 68-75.

like" water in hydrates and hydroxides. Appl. Geochem. 16: 161-181.

ou inertes. Compt. Rendus Académie Science Paris 329: 377-388.

New York: Hemisphere Publishing Corporation. 271 p.

phase solvates. Inorg. Chem. 46: 9768-9778.

48: 120-131.

38: 485-492.

161.

443.

693.

Growth Charact. 30:23-81.

Sov. Radiochem. 29): 653-655.


[30] Keren R, Shainberg I (1980) Water vapor isotherms and heat of immersion of Na/Ca-Montmorillonite systems - 3. Thermodynamics. Clays Clay Minerals 28: 204-210.

366 Thermodynamics – Fundamentals and Its Application in Science

sodium-exchanged form. Langmuir 8: 2730-2739.

Chemical Society. pp. 105-131.

200.

151.

286.

[14] Barrer RM, Cram PJ (1971) Heats of immersion of outgassed and ion-exchanged zeolites. In: L.B. FEMS, editor. Molecular sieves zeolites. Washington D.C.: American

[15] Cases JM, Bérend I, Besson G, François M, Uriot JP, Thomas F, et al. (1992) Mechanism of adsorptionand desorption of water vapor by homoionic montmorillonite. I.The

[16] Bérend I, Cases JM, François M, Uriot JP, Michot L, Masion A, et al. (1995) Mechanism of adsorption and desorptionof water vaporby homoionic montmorillonites:2. The Li+,

[17] Cases JM, Bérend I, François M, Uriot JP, Michot L, Thomas F (1997) Mechanism of adsorption and desorption of water vapor by homoionic montmorillonite: 3. the Mg2+,

[18] Keren R, Shainberg I (1975) Water vapor isotherms and heat of immersion of Na/Ca-Montmorillonite systems. I Homoionic clay. Clays Clay Minerals 23: 193-

[19] Keren R, Shainberg I (1979) Water vapor isotherms and heat of immersion of Na/Ca-Montmorillonite systems - 2. Mixed systems. Clays Clay Minerals 27: 145-

[20] Dios Cancela G, Huertas FJ, Romero Taboada E, Sánchez-Rasero F, Hernández Laguna A (1997) Adsorption of Water Vapor by Homoionic Montmorillonites. Heats of

[21] van Olphen H (1965) Thermodynamics of interlayer adsorption of water in clays. I.--

[22] Carey JW, Bish DL (1996) Equilibrium in the clinoptilolite - H2O system. Amer.

[23] Fialips CI, Carey JW, Bish DL (2005) Hydration-dehydration behavior and

[24] Koster van Groos AF, Guggenheim S (1984) The effect of pressure on the dehydration reaction of interlayer water in Na-montmorillonite (SWy-1). Amer. Mineral. 69: 872-879. [25] Wu T-C, Bassett WA, Huang W-L, Guggenheim S, Koster van Groos AF (1997) Montmorillonite under high H 2 O pressures; stability of hydrate phases, rehydration

[26] Koster v, Groos, A.F., Guggenheim S (1986) Dehydration of K exchanged montmorillonite at elevated temperatures and pressures. Clays Clay Minerals 34: 281-

[27] Petrova N, Mizota T, Fujiwara K (2001) Hydration heats of zeolites for evaluation of

[28] Hill TL (1949) Statistical mechanics of adsorption V. Thermodynamics and heat of

[29] Hill TL (1950) Statistical Mechanics of Adsorption. IX. Adsorption Thermodynamics

thermodynamics of chabazite. Geochim. Cosmochim. Acta 69: 2293-2308.

hysteresis, and the effect of interlayer cations. Amer. Mineral. 82: 69-78.

heat exchangers. J. Therm. Anal. Calorim. 64: 157-166.

and Solution Thermodynamics. J. Chem. Phys.s 18: 246-256.

adsorption. J. Chem. Phys. 17: 520-535.

Na+, K+, Rb+ and Cs+ exchanged forms. Clays Clay Minerals 43: 324-336.

Ca2+,Sr2+ and Ba2+ exchanged forms. Clays Clay Minerals 45: 8-22.

Adsorption and Desorption. J. Colloid Interface Sci. 185: 343-354.

Sodium vermiculite. J. Coll. Sci. 20: 822-837.

Mineral. 81: 952-962.


[46] Van Hinsberg VJ, Vriend SP, Schumacher JC (2005) A new method to calculate end-member thermodynamic properties of minerals from their constituent polyhedra I: enthalpy, entropy and molar volume. J. Metamorph. Geol. 23: 165-179.

Thermodynamics of Hydration in Minerals: How to Predict These Entities 369

[62] Carey JW, Bish DL (1997) Calorimetric measurement of the enthalpy of hydration of

[63] Vieillard P, Blanc P, Fialips CI, Gailhanou H, Gaboreau S (2011) Hydration thermodynamics of the SWy-1 montmorillonite saturated with alkali and alkaline-earth

[64] Fripiat J, Cases J, Francois M, Letellier M (1982) Thermodynamic and microdynamic behavior of water in clay suspensions and gels. J. Colloid Interface Sci.

[65] Koster van Groos AF, Guggenheim S (1987) Dehydration of a Ca- and a Mg-exchanged

[66] Ransom B, Helgeson HC (1994) A Chemical and Thermodynamic Model of Aluminous Dioctahedral-2/1 Layer Clay Minerals in Diagenetic Processes - Regular Solution Representation of Interlayer Dehydration in Smectite. Amer. J. Sci. 294: 449-

[67] Ransom BH, H.C. (1994) Estimation of the standard molal heat capacities, entropies, and volumes of 2:1 clay minerals. Geochim. Cosmochim. Acta 58: 4537-4547. [68] Vidal O, Dubacq B (2009) Thermodynamic modelling of clay dehydration, stability and compositional evolution with temperature, pressure and H2O activity. Geochim.

[69] Dubacq B, Vidal O, De Andrade V (2010) Dehydration of dioctahedral aluminous phyllosilicates: thermodynamic modelling and implications for thermobarometric

[70] Helgeson HC, Delany JM, Nesbitt HW, Bird DK (1978) Summary and critique of the thermodynamic properties of the rock-forming minerals. Am. J. Sci. 278-A: 1-229. [71] Johnson GK, Tasker IR, Jurgens R, O'Hare PAG (1991) Thermodynamic studies of

[72] Vieillard P (2010) A predictive model for the entropies and heat capacities of zeolites.

[73] Marcus Y (1994) A simple empirical model describing the thermodynamics of hydration of ions of widely varying charges, sizes, and shapes. Biophys. Chem. 51: 111-

[74] Berman RG, Brown TH (1985) Heat capacity of minerals in the system Na2O-K2O-CaO-MgO-FeO-Fe2O3-Al2O3-SiO2-TiO2-H2O-CO2, representation, estimation and high

[75] Neuhoff PS, Wang J (2007) Heat capacity of hydration. Amer. Mineral. 92: 1358-1367. [76] Fyfe WS, Turner FJ, Verhoogen J. Metamorphic reactions and metamorphic facies. In:

[77] Cruciani G (2006) Zeolites uponheating: Factors governing their thermal stability and

[78] Lassin A, Duplay J, Touret O, Tardy Y (2000) Thermodynamic treatment of hydrates.

montmorillonite (SWy-1) at elevated pressures. Amer. Mineral. 72: 292-298.

cations: a predictive model. Geochim. Cosmochim. Acta 75: 5664-5685.

clinoptilolite. Clays Clay Minerals 45: 826-833.

89: 378-400.

484.

127.

Cosmochim. Acta 73: 6544-6564.

Eur. J. Mineral. 22:823-836.

Memoirs GSoA, editor.1958. p. 253 p.

estimates. Contri. Miner. Petrol. 159: 159.

Zeolites-Clinoptilolite. J.Chem. Thermodyn. 23: 475-484.

temperature extrapolation. Contrib. Miner. Petrol. 89: 168-183.

structural changes. J. Phys. Chem. Solids 67: 1973-1994.

Comptes Rendus Acad. Sci. Paris SER II A 330: 533-540.


[62] Carey JW, Bish DL (1997) Calorimetric measurement of the enthalpy of hydration of clinoptilolite. Clays Clay Minerals 45: 826-833.

368 Thermodynamics – Fundamentals and Its Application in Science

of monovalent cations. J. Chem. Res.: 444-445.

temperatures. U. S. Geol. Surv. Bull., 1259, 256 p.

zeolites. . Amer. Mineral. 94: 565-577.

Molecular Sieve Zeolites 155-170.

York: Harcurt Brace Jovanovich. 1-360 p.

zeolites. . Microporous Mesoporous Mat. 132: 335-351.

Edn. Zeolites 17: 1-230.

Petrol. 119: 155-165.

of America. pp. 403-452.

260-264.

Mineral. 74: 1023-1031.

664.

Res.: 446-447.

[46] Van Hinsberg VJ, Vriend SP, Schumacher JC (2005) A new method to calculate end-member thermodynamic properties of minerals from their constituent polyhedra I:

[47] Chermak JA, Rimstidt JD (1989) Estimating the thermodynamic properties (ΔG°f and ΔH°f) of silicate minerals at 298 K from the sum of polyhedral contributions. Amer.

[48] Li J, Li B, Gao S (2000) Calculation of thermodynamic properties of hydrated borates by

[49] Chen FR, Ewing RC, Clark SB (1999) The Gibbs free energies and enthalpies of formation of U6+ phases: An empirical method of prediction. Amer. Mineral. 84: 650-

[50] Vieillard P, Jenkins HDB (1986) Empirical relationships for estimation of enthalpies of formation of simples hydrates. Part I Hydrates of alkali- metal cations, of hydrogen and

[51] Vieillard P, Jenkins HDB (1986) Empirical relationships for estimation of enthalpies of formation of simples hydrates. Part 2. Hydrates of alkaline-earth metal cations. J. Chem.

[52] Vieillard P, Jenkins HDB (1986) Empirical relationships for estimation of enthalpies of formation of simples hydrates. Part 3. Hydrates of transition metal cations

[53] Robie RA, Waldbaum DR (1968) Thermodynamic properties of minerals and related substances at 298.15°K (25°C) and one atmosphere (1.O13 bars) pressure and at higher

[54] Vieillard P, Mathieu R (2009) A predictive model for the enthalpies of hydration of

[55] Meier WM, Olson DH, Baerlocher C (1996) Atlas of Zeolite Structure types: 4th revised

[56] Meier WM, Olson DH (1971) Zeolites Frameworks In: Flanigen EM, Sand LB, editors.

[57] Barrer RM (1982) Hydrothermal chemistry of Zeolites. Press A, editor. London New

[58] Carey JW (1995) A thermodynamic formulation of hydrous cordierite. Contrib. Miner.

[59] Armbruster T (1993) Dehydration mechanism of clinoptilolite and heulandite: Singlecrystal X-ray-study of Na-poor and K-Mg-rich clinoptilolite at 100K. Amer. Mineral. 78:

[60] Bish DL, Carey JW (2001) Thermal behavior of natural zeolites. In: D.W. BDLM, editor. Natural zeolites:occurrence, properties Applications. Blacksburg: Mineralogical Society

[61] Mathieu R, Vieillard P (2010) A predictive model for the enthalpies of formation of

(Cr2+,Fe2+,Co2+,Ni2+,Cu2+,Zn2+,Cd2+and UO22+). J. Chem. Res.: 448-449.

group contribution method. Physics and Chemistry of Minerals 27: 342-346.

enthalpy, entropy and molar volume. J. Metamorph. Geol. 23: 165-179.

	- [79] Barshad I (1960) Thermodynamics of water adsorption and desorption on montmorillonite. Clays Clay Minerals: 84-101.

**Chapter 14** 

© 2012 Hwang et al., licensee InTech. This is an open access chapter distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

© 2012 The Author(s). Licensee InTech. This chapter is distributed under the terms of the Creative Commons Attribution License http://creativecommons.org/licenses/by/3.0), which permits unrestricted use, distribution,

**Thermodynamics and Kinetics in the Synthesis of** 

The preparation of monodisperse nanoparticles with uniform size, shape and composition has been intensively pursed because of their scientific and technological interests [1-4]. The major advantage of monodisperse particles may be attributed to the uniform properties of individual particles, which makes the property of whole particles strictly controllable. They have been widely used in industries such as pharmacy, catalysts, sensors, film precursors, and information storage. The property of nanoparticles is much more sensitive to their size than that of micro-particles. For example, the florescence of monodisperse CdSe/ZnS core/shell nanoparticles depends strongly on their size [5]. The superparamagnetism also depends strongly on the size of nanoparticles [6]. The properties of these particular size

In general, to prepare the monodisperse nanoparticles in solution, the size had to be selected after somewhat polydisperse nanoparticles were produced. The separation procedures are very laborious and expensive because the size of nanoparticles is too small to be sorted. Moreover, the production yield of monodisperse nanoparticles decreases markedly due to the loss of nanoparticles during the separation procedures. Therefore, many efforts have been made to synthesize directly monodisperse nanoparticles without size selection procedures.

Recently, several methods have been developed successfully to synthesize gram quantity monodisperse nanoparticles directly without size selection procedures. One of them is a slow heating method developed by Hyeon et al. [6,7], which utilizes the burst nucleation followed by sustained growth of particles. Although nanoparticles have some size distribution in the nucleation stage, they became gradually monodisperse during growth. In this case, the principle for the direct synthesis of monodisperse nanoparticles could be approached by sustained growth of nuclei formed by burst nucleation with some initial size

and reproduction in any medium, provided the original work is properly cited.

nanoparticles show great potentials in the field of bio-medicals and electronics.

**Monodisperse Nanoparticles** 

Additional information is available at the end of the chapter

http://dx.doi.org/10.5772/50324

**1. Introduction** 

distribution [8].

Nong-Moon Hwang, Jae-Soo Jung and Dong-Kwon Lee
