Open access peer-reviewed chapter

# Kinetic Approach to Multilayer Sorption: Equations of Isotherm and Applications

By Francisco S. Pantuso, María L. Gómez Castro, Claudia C. Larregain, Ethel Coscarello and Roberto J. Aguerre

Submitted: August 11th 2018Reviewed: November 22nd 2018Published: March 20th 2019

DOI: 10.5772/intechopen.82669

• isotherm
• roughness
• multilayer
• fractal
• free energy

## 1. Introduction

It is well known that the knowledge and understanding of water adsorption isotherms is of great importance in food technology. This knowledge is highly important for the design and optimization of drying equipment, packaging of foods, prediction of quality, and stability during storage.

In order to describe the overall sorption over the whole region of relative pressures of water, an isotherm for multilayer sorption must be used.

In 1938, Brunauer, Emmett, and Teller (BET) [1] extended Langmuir’s monolayer theory [2, 3, 4, 5] to multilayer adsorption. The BET equation derived was applied to a wide variety of gases on surfaces as well as to the sorption of water vapor by food materials [6, 7, 8].

But, the simple BET equation gives a good agreement with experimental data only at relative pressures lower than 0.35 of adsorbate. A great number of researchers have been analyzed this worrying fact, and numerous modifications have been proposed to the BET model to amend this problem [9, 10, 11].

Among these, the three parameters of GAB equation [12, 13, 14, 15], introduce a modification to the BET sorption model. The GAB model is basically similar to BET ones in its assumptions. These authors propose that the state of the adsorbed molecules beyond the first layer is the same but different from that in the liquid state. This equation describes satisfactorily the sorption of water vapor in foods up to water activities of 0.8–0.9 [16, 17, 18, 19]. The main advantage of the GAB equation is that its parameters have physical meaning. This equation has been adopted by West European Food Researchers [20].

For water activities higher than 0.8–0.9, most of the food materials show values of moisture content larger than that predicted by the GAB model. This flaw indicates that state of the adsorbed molecules beyond the first layer introduced by the GAB model is limited to a certain number of sorption layers. Then turn up as plausible to assume a third stage for the water molecules in the outer zone with true liquid-like properties, as postulated by the original BET model.

A three-zone model for the structure of water near water/solid interfaces was proposed by Drost-Hansen [21]; in this model, beyond the monolayer, a zone of ordered molecular structures of water is expected to exist adjacent to a surface, the ordering extending into the bulk liquid. This is a transition region over which one structure decays into another. At sufficiently large distances from the surface, bulk water structure exists.

The BET model and its modifications were developed for an energetically homogeneous flat surface without lateral interaction and are not suitable for highly rough surfaces [22].

This roughness plays a significant role in the determination of the adsorption characteristics [23, 24, 25], since the shape of the adsorbent surface influences the accessibility of the adsorbate to the active adsorption sites. In this chapter, their fractal dimension will characterize the roughness of the adsorbing surfaces. In addition, taking into account the model of the three zones, the derivation of an equation is presented for BET type multilayer isotherms on rough surfaces. This equation takes into account the influence of the adsorbate-adsorbent interaction of all the adsorbed layers.

It is shown that under certain conditions, this equation is reduced to the known classical forms. The capacity of the different isothermal equations to adjust the equilibrium moisture in the food is analyzed.

## 2. Mathematical model

Brunauer, Emmett, and Teller proposed an adsorption surface divided into n segments, having 1, 2, 3, …, i number of layers of adsorbed molecules. According to this model, adsorption and desorption occur at the top of these segments. So, the equilibrium between the uncovered surface so and the first layer s1is:

a1PPs0=b1s1expE1RTE1

where a1 and b1 are adsorption and desorption coefficients, the same meaning as in BET theory, E1 is the heat of adsorption of the first layer, R is the gas constant, T is the temperature, and P is the vapor pressure of adsorbate. Between any successive layers, the equilibrium can be expressed as:

aiPPsi1=bisiexpEiRTE2

Being si1and sithe surfaces at the top of the respective i−1 and i layers. Considering that Rlnbi/aiis the sorption entropy of the i-layer, Eq. (2) can be written in a more convenient form:

PP=sisi1expΔGiRTE3

where ΔGiis the sorption free energy of the i-layer. This development differs from the classical BET model in that ΔGifor all layers above the first is not considered equal to free energy of bulk liquid adsorbate, ΔGL. Assuming that the free energy of sorption for the i-layer differs from the free energy of bulk liquid by a certain amount, it can be written in general that:

ΔGi=ΔGL+ΔGieE4

where ΔGiedifferentiates the state of the adsorbed molecules from that of the molecules in the pure liquid. Substituting Eq. (4) in Eq. (3), it results:

si=ωisi1E5

where

ωi=PPexpΔGLRTexpΔGieRTE6

Defining

hi=expΔGieRTE7

and given that P0=PexpΔGL/RT, it results:

ωi=PP0hi=xhiE8

being x=P/P0. The fraction of surface occupied by 1st, 2nd, …, ith layer follows the relation:

si=s1j=2iωjE9

Combining Eqs. (8) and (9), we have

si=s1j=2ixhjE10

As s1=ω1s0, it results:

si=h1s0xij=2ihj=Cs0xij=2ihjE11

where C=h1=expΔG1e/RTis the constant C of BET theory.

Given that the adsorbate molecules, considered as spheres, when adsorbed osculate the surface, for a fractal surface the relationship between the surfaces at the top, si, and the bottom, si, for a molecular stack of i-layers is [26]:

si=si2i12DE12

where D is the fractal dimension; when Eq. (12) is substituted in Eq. (11), it gives:

si=Cs0xi2i1D2j=2ihjE13

According to the BET theory, the monolayer capacity, Nm, is:

Nm=1σi=0si=s0σ1+Cx+i=2Cxi2i1D2j=2ihjE14

where σ is the cross-sectional area of water molecule. The total amount of adsorbent in a given layer n is:

Nnt=2n12Dσi=nsi=Cs0σ2n12Di=n2i1D2xij=2ihjE15

The total number of molecules, N, that form the adsorbed film is:

N=1σs1+i=2sik=1i2k12D=Cs0σx+i=2xi2i1D2j=2ihjk=1i2k12DE16

but, from Eq. (15), it is also:

N=i=1Nit=Cs0σx+i=22i12Dj=ixj2j1D2k=2jhkE17

Finally, combining Eqs. (14) and (16), the following general equation for sorption isotherms is found:

NNm=Cx+i=2xi2i1D2j=2ihik=1i2k12D1+Cx+i=2Cxi2i1D2j=2ihjE18

But combining Eqs. (14) and (17), other equivalent form of the equation for sorption isotherms is reached:

NNm=Cx+i=22i12Dj=ixj2j1D2k=2jhk1+Cx+i=2Cxi2i1D2j=2ihjE19

Eq. (18) is therefore the isotherm equation for multilayer adsorption on fractal surfaces that takes into account the variation of the free energy of adsorption with successive layers.

### 2.1 Applications of the model to smooth surfaces

For a nonfractal surface (D = 2), Eqs. (18) and (19) reduce, respectively, to:

NNm=Cx+i=2ixij=2ihi1+Cx+Ci=2xij=2ihjE20
NNm=Cx+i=2j=ixjk=2jhk1+Cx+Ci=2xij=2ihjE21

It is interesting to comment that for hi = 1 and i ≥ 2 (free energy of the multilayer equal to the free energy of bulk water), Eq. (20) reduces to BET equation (see Eq. (22) in Table 1).

hiEquationReferences
1NNm=Cx1x1x+Cx(22)[1]
kNNm=Ckx1kx1kx+Ckx(23)[27]
ii1NNm=Cx1+x1x1x2+Cx(24)[28]
i1iNNm=Cx1x1Cln1x(25)[28]

### Table 1.

Equations derived from Eq. (20).

Even more, if for the second and higher layers the free energy of the adsorbate differs from that of pure liquid in a constant amount (h = k), Eq. (18) reduces to GAB equation (see Eq. (23) in Table 1).

Assuming that the adsorbate properties approach to the pure liquid as i increases (limiΔGie=0), h2 > h3 > … > 1 or h2 < h3 < … < 1, depending on arrangement of the adsorbate in the multilayer region (see Eqs. (24) and (25) in Table 1).

Eqs. (22) and (23) are frequently used in bibliography and tested for different adsorbate/adsorbent systems [29, 30, 31].

The ability of Eq. (24) to fit experimental data of water sorption on different food products is presented in Table 2.

MaterialTemperatureCNm%, d.b.E%*References
Milk products
Edam cheese (a)25°C27.62.14.7[32]
Emmental cheese (a)25°C77.71.910.3
Yoghurt (a)25°C25.22.82.2
Fruits
Apple (d)20°C60.34.67.1[33]
Apricots (a)30°C1.73.910.4
Banana (a)25°C1.54.818.1[32]
Figs (a)30°C1.74.95.8[33]
Pear (a)25°C6.15.26.2[32]
Pineapple (a)25°C4.74.725
Plums (a)30°C1.94.910.8[33]
Sultana raisins (a)30°C48.44.58.3
Vegetables
Carrots (s)37°C2.94.55.8[34]
Onion, tender (a)25°C64.74.15[32]
Spinach (s)37°C64.33.21.2[34]
Sugar beet (d)25°C14.34.36.4[35]

### Table 2.

Food sorption isotherms fitted with Eq. (24).

E%=100nNpNe/Npp:predictede:experimental

(a), adsorption; (d), desorption, (s), sorption.

In Table 3, the fitting test corresponding to Eq. (25) can be seen.

MaterialTemperatureCNm%, d.b.E%*References
Starchy foods
Barley (d)25°C16.68.47.5[36]
Corn shelled (d)25°C21.68.57.5[37]
Flour (a)30°C57.76.89.1[38]
Hard red winter (d)25°C23.19.35.4[39]
Native manioc starch (d)25°C22.49.56.8[37]
Native potato starch (a)25°C11.911.53.1
Oats (d)25°C16.288.2[37]
Rough rice (d)60°C8.25.73.2[40]
Rye (d)25°C26.48.37[37]
Sorghum (a)38°C19.08.46.1[41]
Starch (a)30°C9.69.41.6[38]
Tapioc starch (d)25°C19.910.88.2[37]
Wheat, durum (d)25°C22.19.15[39]
Nuts and oilseeds
Moroccan sweet almonds (a)25°C6.53.72.2[42]
Para nut (a)25°C9.42.74.7[32]
Peas, dried (a)25°C75.18.33.3[42]
Pecan nut (a)25°C6.82.94.7[32]
Rapeseed, guile (d)25°C10.74.63.2[43]
Rapeseed, tower (a)25°C8.452.3
Soybean seed (s)15°C5.685.4[44]
Meats
Beef, minced (a)30°C4.49.25.4[45]
Beef, raw/minced (s)10 °C2.314.38.4[46]
Cod, freeze dried/unsalted (a)25°C5.612.25.8[47]
Fish flour (a)25°C4.26.85.3[48]
Mullet roe, unsalted (a)25°C5.34.81.4[49]
Mullet, white muscle (a)25°C6.110.13.7

### Table 3.

Food sorption isotherms fitted with Eq. (25).

E%=100nNpNe/Npp:predictede:experimental

(a), adsorption; (d), desorption, (s), sorption.

Eq. (25) gives a good agreement using data of amilaceous materials, nuts, and meats, whereas Eq. (24) shows a good fitting with fruits, some vegetables, and milk products.

#### 2.1.1 Limited sorption

If the number of the adsorbed layers cannot exceed some finite number n, then for D = 2, hi = 1, and i ≥ 2, from Eq. (18)

NNm=Cx1x1n+1xn+nxn+11x+CxCxn+1E28

also obtained by Brunauer et al. [1].

But from Eq. (19), D = 2, hi = 1, and i ≥ 2

NNm=Cx1x1xn1x+CxE29

known as Pickett [10] or Rounsley [50] isotherm equation.

### 2.2 The fractal isotherm

If the surface of the adsorbent behaves like a fractal with 2 < D < 3 and assuming that the free energy distribution in the adsorbed film is the same like in BET theory (hi = 1), the following fractal isotherm equations can be obtained from Eqs. (18) and (19), respectively:

NNm=Cx+i=2xi2i1D2j=1i2j12D1+Cx+Ci=2xi2i1D2E30
NNm=Cx+i=22i12Dj=ixj2j1D21+Cx+Ci=2xi2i1D2E31

To illustrate the effect of roughness on the shape of the isotherms, in Figure 1, the influence of D values, for C = 20, can be seen.

Table 4 illustrates the usefulness of Eq. (30) on starchy materials.

MaterialTemperatureCNm%, d.b.DE%*Source
Starches
Wheat (a)30°C9.18.12.84.3[38]
Corn (d)25°C12.59.82.91.0[37]
Potato (a)25°C4.09.52.83.3[51]
Tapioca (d)25°C14.310.22.92.4[37]
Native manioc (a)25°C15.89.32.81.7[37]
Cereals
Rough rice (d)25°C10.79.02.92.6[52]
Rough rice (d)40°C8.38.92.91.3[53]
Sorghum (d)20°C16.49.62.90.9
Sorghum (d)50°C6.98.82.91.1
Wheat durum (d)25°C12.68.72.80.9[39]
Corn, continental (d)20°C12.59.52.91.0[53]
Corn, continental (d)50°C4.07.92.81.9

### Table 4.

Food sorption isotherms fitted with Eq. (30).

E%=100nNpNe/Npp:predictede:experimental

(a), adsorption; (d), desorption; (s) sorption.

In all the modeled materials, it is found that the value of the parameter D is approximately 2.8. This indicates that the tested products show a high roughness.

### 2.3 A four parameters equation

Assuming that the total surface area available for sorption is formed by two types of surfaces or regions: (a) a region representing a fraction, α, of the total adsorbing surface that only adsorbs a limited number of adsorbate layers, that is, internal surface such as pores; (b) the remainder fraction, (1 − α), of the total adsorbing surface, where unlimited sorption may occur, that is, external surface and macropores where large values of n are required to fill them. From Eq. (29), it can be written [54]:

NNm=Cx1αxn1x1x+CxE32

In Figure 2, data of champignon mushroom and its fit using Eqs. (32) and (29) for different temperatures are presented.

In Table 5, the fitting of moisture sorption on different products using Eq. (32) can be seen.

MaterialT°CCNm%, d.b.αnE%*References
Champignon mushroom (a)
(Agaricus bisporus)
302.4220.280.84781.410.19[55]
402.3618.520.85561.371.49
502.6115.590.85881.501.91
601.7917.800.90371.261.78
701.6415.680.89451.302.71
Casein (a)2512.045.740.98643.720.97[51]
Casein (d)259.727.960.89232.340.77
Coffee (a)202.213.230.878414.943.62
Dextrin (a)1012.0213.010.87520.990.98
Potato starch (a)208.578.030.92303.471.21
Anis (a)2515.034.320.99009.421.15[32]
Banana (a)250.4118.650.64323.879.06
Cardamom (a)2525.306.020.67342.640.67
Celery (a)255.397.060.84069.093.43
Chamomile (a)2516.786.070.50104.670.94
Emmenthal (a)259.683.420.960413.251.66
Cinnamon (a)2520.186.260.99993.630.19
Clove (a)2529.724.250.66383.081.45
Coriander (a)2510.955.870.78142.210.84
Eggplant (a)256.567.760.43967.964.00
Fennel (a)250.3370.060.97770.227.60
Forelle (a)456.635.120.918512.794.86
Ginger (a)2515.237.300.80492.291.11
Huhn (a)456.895.390.925310.931.99
Joghurt (a)255.185.150.942119.538.07
Laurel (a)2517.584.380.81255.952.37
Lentils (a)2517.676.950.87013.630.76
Marjoram (a)2520.244.940.40062.912.37
Mint (a)2512.697.420.64931.900.88
Nutmeg (a)2526.934.570.81702.820.86
Para nut (a)2526.851.810.71004.100.90
Pears (a)251.6412.310.92007.962.65
Pecan nut (a)2513.601.950.98694.901.38
Pineapple (a)250.4626.240.96882.633.95
Savory (a)2528.366.610.84973.441.97
Thyme (a)2523.394.770.67294.431.44
Rice, rough (d)403.1211.300.91481.450.56[56]
502.6011.020.92641.341.68
601.9411.200.95201.391.67
701.5611.720.95081.282.10
801.4010.940.95441.452.17
Meat, raw minced (s)1011.416.480.980710.712.42[46]
307.825.990.971110.204.91
5013.674.780.974612.801.54
Lard (s)2513.960.360.66264.101.07
Mullet roe, unsalted (a)253.147.370.79570.650.73[49]
Mullet, white muscle (a)259.467.120.96814.861.73
Cod, unsalted (a)2513.027.780.63363.674.79[47]

### Table 5.

Food sorption isotherms fitted with Eq. (32).

E%=100nNpNe/Npp:predictede:experimental

(a), adsorption; (d), desorption; (s), sorption.

Eq. (32) gives a good agreement with experimental data. The inclusion of the alpha parameter allows modifying the amplitude of the isotherm plateau. The modeling of the sigmoid isotherms is facilitated, typical form found in the adsorption of water in food products.

## 3. Conclusions

In the framework of the BET model, a general isotherm equation was obtained that includes the roughness of the adsorbent surface and characterizes the transition region between the monolayer and the outer zone where the adsorbate has the properties of the bulk liquid through a free energy excess that differentiates the adsorbed phase from the bulk liquid.

This general equation, depending on the simplifications assumed, gives the classical BET and GAB equations. But, taking into account an asymptotic reduction of the free energy excess, for a flat surface, two different equations were obtained. One of them appears useful to model starchy materials, extending the isotherm plateau, and the other, successfully model fruit isotherms, reducing the isotherm plateau.

Considering only the roughness, and assigning bulk liquid properties at all layers beyond the first, an equation that includes the fractal dimension is obtained.

This fractal dimension can vary from 2 to 3. The rising of its value result in an outspread the isotherm plateau. Particularly, for highly rough surfaces, the multilayer growing is limited by geometrical restrictions. In this case, the magnitude of the interactions practically has no effect on the shape of the isotherm.

It results from the present analysis that modifications of the BET model based only on the three-zone model or geometric considerations conduct to similar results.

So, Eqs. (25), (30), and (32) predict lower sorption capacity with the increment of water activity, giving better agreements with experimental isotherms than the classical BET equation.

This fact forewarns that the fractal dimension in the model could be affected from unsuitable accounting for the adsorbate-adsorbent interactions.

## Acknowledgments

The authors gratefully acknowledge the Morón University for providing financial support.

## More

© 2019 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

## How to cite and reference

### Cite this chapter Copy to clipboard

Francisco S. Pantuso, María L. Gómez Castro, Claudia C. Larregain, Ethel Coscarello and Roberto J. Aguerre (March 20th 2019). Kinetic Approach to Multilayer Sorption: Equations of Isotherm and Applications, Food Engineering, Teodora Emilia Coldea, IntechOpen, DOI: 10.5772/intechopen.82669. Available from:

### Related Content

Next chapter

#### Risk Management of Egg and Egg Products: Advanced Methods Applied

By Poliana Mendes de Souza, Regiane de Melo, Miriam Aparecida de Aguilar Santos, Fabiana Regina Lima and Kássia Héllen Vieira

#### Structure and Function of Food Engineering

Edited by Ayman Amer Eissa

First chapter

#### Antioxidant, Anticancer Activity, and Other Health Effects of a Nutritional Supplement (Galaxy(r))

By Gary M. Booth, Tory L. Parker and Christopher M. Lee

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

View all Books