Comparative Analysis of Statistical Thickness Models for the Determination of the External Specific Surface and the Surface of the Micropores of Materials: The Case of a Clay Concrete Stabilized Using Sugar Cane Molasses ()
1. Introduction
When characterizing a porous material for use in industry, specific surface area is a key property. However, textural properties are crucial for the very wide range of applications for which the material is used [1] . This is the case for clay soils stabilized with cane molasses for use in road construction. These specific elements are required for the physical adsorption of nitrogen and provide information on textural properties, in particular the total accessible surface area and porosity distribution.
For the record, the results given in the form of curves relating the quantities of nitrogen adsorbed as a function of the ratio between partial pressure and saturation vapour pressure for nitrogen, are called isotherms. Adsorption isotherms can also be used to calculate the external surface area and the microporous surface area using the t-plot method.
The external specific surface (or external area) represents the extent of a non-porous material as well as the surface developed by pores large enough to allow the formation of a multimolecular layer whose thickness t increases regularly with (P/P0). Similarly, the BET method is generally used to characterize the specific surface area of porous materials. It can be used to obtain the external specific surface area σext by including the presence of micropores whose surface is not accessible.
The t-plot method can be used to determine the external specific surface σext by excluding the micropores, and gives a more representative surface in terms of accessibility, by considering the thickness of the multimolecular layer of adsorbed nitrogen [1] .
Nevertheless, several statistical thickness models are given in the literature. Studies have already been carried out by NGOUALLAT et al. (2022), within the framework of microstructure analysis and the determination of isotherms and specific surfaces [2] . These studies show that the quantity of molasses in the materials does not modify the type IV nitrogen adsorption isotherm, which remains and a type H4 hysteresis loop in all the samples, which justifies the monolayer and multilayer absorption mechanism; these are therefore mesoporous materials [3] .
From a microstructural point of view, we have also observed the presence of inter-aggregate pores (mesopores) in the various samples analyzed, which suggests an evolution of the open soil structure towards a dense granular matrix [4] . However, the aim of this study is to compare four statistical thickness models: the Carbon Black model, the Harkins and Jura model, the Hasley model and the Jaroniec et al. model, in the case of clayey concrete stabilized with sugarcane molasses, in order to find the optimum model that best simulates the phenomenon.
2. Materials and Methods
The t-plot statistical thickness models most commonly used to determine the textural properties of materials are: the Harkins and Jura model, the Carbon Black model, the Halsey model and the Jaroniec et al. model [5] [6] . Other models exist, and we present those used in this article.
2.1. Statistical Thickness: Harkins and Jura Model
The empirical value of the statistical thickness t is expressed in the form of the equation of Harkins et al. [1] [7] [8] [9] [10] , valid for as long as (P/P0) < 0.8:
(1)
2.2. Statistical Thickness: Carbon Black Model
The Magee proposes a calculation of the statistical thickness t for carbon black valid for 0.2 < (P/P0) < 0.5 [1] [6] [11] :
(2)
2.3. Statistical Thickness: Halsey Model
Halsey’s statistical thickness is represented by the following equation [11] [12] :
(3)
2.4. Statistical Thickness: Jaroniec et al. Model
Jaroniec’s statistical thickness is given by the following equation [8] :
(4)
2.5. External Specific Surface of the Material
The external specific surface
is obtained by determining the slope (standard linear regression) from the graph of the quantity of nitrogen adsorbed per gram of sample (Qa) as a function of the statistical thickness t. The value of
is calculated using equation Halsey statistical thickness is represented by the following equation [11] Halsey statistical thickness is represented by the following equation [1] [13] :
(5)
With
P, the value of the slope of the curve va = f (t);
15.47: a constant related to the conversion of the volume of nitrogen and the units in m2/g.
If the statistical thickness is given in nm, a conversion factor is added, giving:
(6)
2.6. Microporous Surface of the Material
The difference between the specific surface area value determined by the BET method and the external specific surface area obtained by the t method is used to determine the microporous surface area according to equation [1] [13] :
(7)
2.7. Model Selection Criteria
Four criteria are used to compare and select models: the adjusted R2 criterion, the Akaike information criterion (AIC), the Bayesian information criterion (BIC) and the Chi-square. The best model is the one with the lowest AIC, the lowest BIC, the highest R2adj and the lowest Chi-square [14] . The values of the criteria were obtained using Origin pro software.
2.8. Physical Adsorption of Nitrogen
Nitrogen adsorption experiments on clay concretes were carried out using a Micro Active for ASAP 2460 Version 2.01 apparatus.
2.9. Geotechnical Characterization of Clay Soil
Table 1 gives the geotechnical characteristics of the clay soil.
2.10. Sugar Cane Molasses
The sugar cane molasses used comes from the “Société Agricole de Raffinage Industriel du Sucre (SARIS-Congo)”, a sugar industry organized in the town of Nkayi, Republic of Congo. The molasses used has the following characteristics: the Brix value is 82.85%, which represents a sugar content of 82.85%. The corresponding polarity is 29.07% and the purity 35.09% [2] .
2.11. Composition of Stabilized Clay Concretes
Table 2 shows the composition of stabilized and non-stabilized clay concretes.
![]()
Table 1. Geotechnical characterization of the clay soil used [15] [16] .
![]()
Table 2. Quantity of materials used [2] [15] .
![]()
Table 3. Experimental nitrogen adsorption ratio and statistical thicknesses for clay concrete without molasses.
3. Results and Discussion
3.1. Samples without Cane Molasses
The experimental ratio of nitrogen absorption on clayey concrete, in particular the relative pressure (P/P0) and the quantity of nitrogen adsorbed, is given in Table 3. The statistical thickness values according to the models obtained using: formula 1 for the Harkins and Jura model, formula 2 for the Carbon Black model, formula 3 for the Halsey model, formula 4 for the Jaroniec et al. model, are also presented in Table 3.
The shape of the curves in Figure 1 reveals capillary condensation in the pores that begins as the slope starts to increase, which justifies the presence of condensed water molecules in the mesopores, between the particles and grains in the clay concrete, which is obvious because clay concrete is water-based with an initial water content of 63% (Table 2).
Figure 2 shows the linear interpolation of the curves for the amount of nitrogen adsorbed per gram of sample.
Linear interpolation of the curves showing the quantities of nitrogen adsorbed as a function of statistical thickness (Figure 2), enabled us to obtain linear regression lines of the form Qads = Pt + Q0 (y = Qads, x = t (statistical thickness)). The slope P is used to calculate the external surface of the material (formula 6). The ordinate at the origin Q0 is used to obtain the volume of the micropores. However, in this study we are only interested in the surface area of the micropores. Table 4 gives the values of the model selection criteria clay concrete.
![]()
Figure 1. Quantity of nitrogen adsorbed per gram of sample for clay concrete without molasses as a function of statistical thickness t (carbon black, Halsey, Harkins and Jura, Jaroniec).
![]()
Figure 2. Linear interpolation of the curves for the quantity of nitrogen adsorbed per gram of sample for clay concrete without molasses as a function of the statistical thickness t (carbon black, Halsey, Harkins and Jura, Jaroniec).
![]()
Table 4. Comparison of models using the AIC, BIC, chi-square and adjusted coefficient of determination criteria for clay concrete.
![]()
Table 5. Experimental ratio and statistical thicknesses for clay concrete stabilised at 8%.
The Hasley model has the highest
value, the lowest
value, the lowest BIC value and the lowest AIC value (Table 4), indicating that it is the best performing of the four models, followed by the carbon black model, then the Jaroniec model and finally the Harkins and Jura model.
3.2. Samples Stabilized with 8% Cane Molasses
Table 5 gives the ratio of nitrogen absorption on clay concrete stabilized at 8%, i.e. the relative pressure (P/P0) and the corresponding quantity of adsorbed nitrogen. It also gives the statistical thickness values calculated according to the Carbon Black model, the Harkins and Jura model, the Halsey model and the Jaroniec model.
The shape of the curves (Figure 3) reveals multilayer adsorption on a surface with low porosity, which can be explained by the occupation of the pores by sugarcane molasses. We observe a straightening of the curves (linearization) compared with the curves obtained for non-stabilized clay concrete. This corresponds to the absence of capillary condensation in the mesopores between the particles and the grains. In fact, the presence of molasses, mainly made up of sucrose, glucose and fructose [2] , which are molecules rich in hydroxyl groups (OH) in the mesopores of the material, favours the establishment of hydrogen bonds with the water molecules. The latter are then integrated into the molasses macromolecules in the material’s mesopores. The linearization of these curves will become clearer as the molasses content of the clay concrete increases.
The linear interpolation of the curves for the quantity of nitrogen adsorbed per gram of sample (Figure 3) is given in Figure 4.
![]()
Figure 3. Amount of nitrogen adsorbed per gram of sample for clay concrete stabilised at 8% as a function of statistical thickness t (carbon black, Halsey, Harkins and Jura, Jaroniec).
![]()
Figure 4. Linear interpolation of the curves for the quantity of nitrogen adsorbed per gram of sample for clay concrete stabilized with 8% molasses as a function of the statistical thickness t (carbon black, Halsey, Harkins and Jura, Jaroniec).
![]()
Table 6. Comparison of the models using the AIC, BIC, chi-square and adjusted coefficient of determination criteria for clay concrete stabilized at 8%.
![]()
Table 7. Experimental ratio and statistical thicknesses for 12% stabilized clay concrete.
Table 6 gives the values of the model selection criteria for clay concrete stabilized with 8% molasses.
The Harkins and Jura model has the highest R2adj value, the lowest χ2 value, the lowest BIC value and the lowest AIC value (Table 6), indicating that it is the best performing of the four models, followed by the Jaroniec model, then the Halsey model and finally the Carbon black model.
3.3. Samples Stabilized with 12% Cane Molasses
Table 7 gives the experimental ratio of nitrogen absorption on clay concrete stabilized at 12%: the relative pressure (P/P0) and the quantity of nitrogen adsorbed. The statistical thickness values according to the models are calculated using formulae 1, 2, 3 and 4 (page 3).
The shape of the curves (Figure 5) shows a linearization with respect to the curve at 0% and 8%, which corresponds to multilayer adsorption on a surface with low porosity, which can be justified by the occupation of the pores by sugarcane molasses.
The linear interpolation of the curves for the quantity of nitrogen adsorbed per gram of sample (Figure 5) is shown in Figure 6.
![]()
Figure 5. Amount of nitrogen adsorbed per gram of sample for clay concrete stabilized at 12% as a function of statistical thickness t (carbon black, Halsey, Harkins and Jura, Jaroniec).
![]()
Figure 6. Linear interpolation of the curves for the quantity of nitrogen adsorbed per gram of sample for clay concrete stabilized with 12% molasses as a function of the statistical thickness t (carbon black, Halsey, Harkins and Jura, Jaroniec).
Table 8 gives the values of the model selection criteria for 12% stabilized concrete.
The Harkins and Jura model has the highest R2adj value, the lowest χ2 value, the lowest BIC value and the lowest AIC value (Table 8), indicating that it is the best performing of the four models, followed by the Jaroniec model, then the Halsey model and finally the Carbon black model.
3.4. Samples Stabilized with 16% Cane Molasses
The experimental ratio of nitrogen absorption on clay concrete stabilized at 16% and the statistical thickness values calculated using the model formulae (page 3) are shown in Table 9.
The curves in Figure 7 are more linear than at 8% and 12% and at 0%, which corresponds to multilayer adsorption on a surface with low porosity, which may be justified by the prior occupation of the pores by sugarcane molasses. The adjusted correlation coefficients obtained for the Harkins and Jura model: 0.99836 (Table 6) for clay concrete stabilized at 8%, 0.99825 (Table 8) for clay concrete stabilized at 12%, 0.99936 (Table 10) for clay concrete stabilized at 16%, show an increase in linearization (Figure 8) with molasses content, which can be explained by the absence of condensation in the pores, as the molasses molecules occupy these pores.
![]()
Table 8. Comparison of the models using the AIC, BIC, chi-square criteria and the adjusted coefficient of determination for clay concrete stabilized at 12%.
![]()
Table 9. Experimental ratio and statistical thicknesses for clay concrete stabilized at 16%.
The linear interpolation of the curves for the quantity of nitrogen adsorbed per gram of sample is shown in Figure 8.
![]()
Figure 7. Amount of nitrogen adsorbed per gram of sample for clay concrete stabilized at 16% as a function of statistical thickness t (carbon black, Halsey, Harkins and Jura, Jaroniec).
![]()
Figure 8. Linear interpolation of the curves for the quantity of nitrogen adsorbed per gram of sample for clay concrete stabilized with 16% molasses as a function of the statistical thickness t (carbon black, Halsey, Harkins and Jura, Jaroniec).
Table 10 gives the values of the model selection criteria for 16% stabilized concrete.
The Harkins and Jura model has the highest R2adj value, the lowest χ2 value, the lowest BIC value and the lowest AIC value (Table 10), indicating that it is the best performing of the four models, followed by the Jaroniec model, then the Halsey model and finally the Carbon black model.
Table 11 gives the values of slope and external specific surface as a function of the statistical thickness models.
Table 12 gives the BET specific surface area and micropore surface area values for stabilized and non-stabilized concretes according to the models used.
The results obtained show that the Halsey model performs best for non-stabilized clay soil. However, the external specific surface area values obtained by this model are higher than the BET specific surface area values, so it is not possible to obtain the micropore surface area by this model. Thus, the Carbon Black model becomes the most relevant for non-stabilized clayey concrete, followed by the Jaroniec model, and finally by the Harkins and Jura model.
![]()
Table 10. Comparison of models using the AIC, BIC, chi-square and adjusted coefficient of determination criteria for clay stabilized at 16%.
![]()
Table 11. Slope and external specific surface values as a function of statistical thickness models.
![]()
Table 12. BET specific surface area and micropore surface area according to the models used.
The results obtained for stabilized clay concretes show that the Harkins and Jura model is the most relevant, followed by the Jaroniec model and, finally, the Carbon Black model. The Halsey model does not provide micropore surfaces for the same reason as mentioned above.
This change compared to unsterilized clayey concrete leads us to believe that stabilization using molasses brings about modifications to the external surface of clayey concrete by occupying the pores [4] . Ngouallat Mfoutou (2020) reports that 53.108 molasses molecules occupy the accessible surface of clay concrete [3] .
The change in the shape of the curve from the curved form for clay concrete without molasses to the more linear form for stabilized clay concretes also expresses the change in texture of clay concrete from a porous material to a less and less porous material.
According to Tchemmou and Gherbi (2018), the t-plot method can be used to determine the external surface of microporous materials of the zeolite type [7] . For Magee (1995) cited Moulin (2018) [1] , the Carbon black model is designed for the specific case of carbon blacks. Nevertheless, the results obtained show that the Carbon black model is also relevant for the case of clay concretes. Ngouallat (2022) used the Carbon Black statistical thickness model to obtain the textural properties of clay concretes stabilized with molasses [2] .
According to Jeffrey Kevin (2008), the Carbon black statistical thickness model can also be used to determine the external surface of cetyltrimethyl ammonium bromide [3] .
The carbon black and Halsey equations (Eq. 2 and 3) are applicable to carbon black and related materials [15] .
According to Yijing Zheng (2008), the t-plot method (Halsey model, Harkin and Jura model, Carbon black model) used for the physical characterization of common adsorbents may not be applicable to carbon nanotubes [5] .
The Harkins and Jura equation has been developed for well-selected alumina samples; however, it is applicable to other materials, including graphitized carbon blacks [17] .
The most commonly applied thickness curves are those of Harkin-Jura and Halsey (Webb and Orr, 1997) [6] . Utpalendu Kuila and Manika Prasad (2011), have used Halsey thickness curves to estimate micropore volume and “open area”, mesopore area, macropore area and external area (i.e. total area excluding micropore area) in natural clay minerals and shales [18] .
4. Conclusion
The gist of this study has been to compare statistical thickness models for the determination of textural properties for the case of clay concrete and for stabilized clay concretes. The results showed that:
- Halsey’s model can be used to obtain external surfaces. However, it does not allow the surface of the micropores to be obtained. It is not suitable for clay concrete and for clay concretes stabilized with molasses.
- The Carbon Black model, the Jaroniec model and the Harkins and Jura model can be used for clay concrete and for stabilized clay concrete.
- The Carbon Black model is the most relevant for clay concrete.
- The Harkins and Jura model is the most appropriate for molasses-stabilized clay concrete.