Next Article in Journal
Gum Tragacanth (GT): A Versatile Biocompatible Material beyond Borders
Next Article in Special Issue
Fractional Dynamics of HIV with Source Term for the Supply of New CD4+ T-Cells Depending on the Viral Load via Caputo–Fabrizio Derivative
Previous Article in Journal
Zeolites at the Molecular Level: What Can Be Learned from Molecular Modeling
Previous Article in Special Issue
Removal of Cu(II) Contamination from Aqueous Solution by Ethylenediamine@β-Zeolite Composite
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Finite Pore Volume GAB Adsorption Isotherm Model as a Simple Tool to Estimate a Diameter of Cylindrical Nanopores

by
Sylwester Furmaniak
1,*,
Piotr A. Gauden
2,
Maria Leżańska
3,
Radosław Miśkiewicz
4,
Anna Błajet-Kosicka
1 and
Piotr Kowalczyk
5
1
Stanisław Staszic State University of Applied Sciences in Piła, Podchorążych Street 10, 64-920 Piła, Poland
2
Carbon Materials Application in Electrochemistry and Environmental Protection Research Group, Faculty of Chemistry, Nicolaus Copernicus University in Toruń, Gagarin Street 7, 87-100 Toruń, Poland
3
Research Group for Modeling and Synthesis of Novel Materials, Faculty of Chemistry, Nicolaus Copernicus University in Toruń, Gagarin Street 7, 87-100 Toruń, Poland
4
Faculty of Economics and Management, University of Szczecin, Aleja Papieża Jana Pawła II 22A, 70-453 Szczecin, Poland
5
School of Engineering and Information Technology, Murdoch University, Murdoch, WA 6150, Australia
*
Author to whom correspondence should be addressed.
Molecules 2021, 26(6), 1509; https://doi.org/10.3390/molecules26061509
Submission received: 24 January 2021 / Revised: 2 March 2021 / Accepted: 3 March 2021 / Published: 10 March 2021

Abstract

:
The finite pore volume Guggenheim–Anderson–de Boer (fpv-GAB) adsorption isotherm model has been considered as a simple tool which not only enables us to analyze the shape of isotherms theoretically, but also provides information about pore diameter. The proposed methodology is based on the geometrical considerations and the division of the adsorption space into two parts: the monolayer and the multilayer space. The ratio of the volumes of these two spaces is unambiguously related to the pore diameter. This ratio can be simply determined from the N2 adsorption isotherm by its fitting with the use of fpv-GAB model. The volume ratio is equal to the ratio of the adsorption capacities in the monolayer and the multilayer—two of the best-fit parameters. The suggested approach has been verified using a series of isotherms simulated inside ideal carbon nanotubes. The adsorption data for some real adsorbents has also been used during tests. The studies performed have proven that diameters estimated with the use of the proposed method are comparable with the geometrical sizes or diameters published by others and based on the application of more sophisticated methods. For pores wider than 3 nm, the relative error does not exceed a few percent. The approach based on the fpv-GAB model reflects well the differences in pore sizes for the series of materials. Therefore, it can be treated as a convenient tool to compare various samples.

1. Introduction

Nanoporous materials are still of wide interest. This interest results from their many applications, which include, among others, water filters [1], energy and gas storage [2,3], medical implants and drug delivery [4,5,6,7], and catalysis [8]. Nanoporous media are also often used in the chemical industry and in broadly understood nanotechnology [9,10]. The main factor determining properties of such materials is their porous structure, i.e., the geometry and the size of pores and the distribution of these parameters, since the porosity is very often heterogeneous. Among the materials of various pore geometry—which may also be irregular—one can distinguish those with cylindrical pores. This inhomogeneous group includes carbon nanotubes [11,12], some ordered carbons [13] and mineral adsorbents, like MCM-41 [14], AlPO4-5 [15], SBA-15 [16], and alumina [17], as other nanoporous materials, also those ones with cylindrical pores, have various applications—see for example [12,18,19,20,21,22,23,24,25,26,27].
The direct determination of a nanopore’s size and its distribution is usually a very difficult task, especially in the case of disordered and amorphous materials with heterogeneous pores. Therefore, the porosity is usually characterized indirectly. The most popular methods used for this purpose are based on the measurements of the low-temperature adsorption isotherms of inert gases—mainly N2 [28,29]. Such approaches allow us to determine the pore size distribution (PSD) curve with the use of the global adsorption isotherm equation and the base of the local isotherms for pores of well-defined geometry differing in size [30,31]. The local isotherms are generated using more or less sophisticated mathematical models or from the results of molecular simulations. The fitting of an experimental isotherm by the theoretical one, leading to the PSD, also requires the use of a complicated mathematical and/or computational methodology. However, in practice, it is often very convenient to use a single parameter (single numerical value) instead of the full PSD curve to compare or to characterize the materials. Such a parameter is usually the average size or the size of the dominant pores. These values can often be determined using simpler methods.
Previously, we proposed a mathematical model to describe adsorption isotherms in cylindrical pores that leads to the estimation of the pore diameter during fitting of the data [32]. The formalism of this approach, the so-called nanotube polymolecular Fowler–Guggenheim adsorption isotherm (NT-PFG), is complicated since the adsorption space is divided into layers of molecular thickness. As a consequence, it is required to model the adsorption amounts separately for each subsequent layer. Moreover, from the adsorption mechanism standpoint, the exact distinction of the subsequent layers, with the exception of the monolayer, and several following layers, is problematic. This is because the borders between them become more and more blurred as the distance from the pore wall increases [33,34]. However, the reliability of the predictions of this model were successfully verified for experimental and simulated data [32]. In the current study, we have suggested a simplified approach based on the division of the cylindrical space into two parts only: the monolayer and the multilayer. The ratio of the volumes of these two spaces is unambiguously related to the pore diameter and can be used to calculate it. This value can be estimated from the adsorption isotherm as the ratio of adsorption capacities in the monolayer and the multilayer. We have proposed the use of a theoretical adsorption isotherm model to determine these adsorption capacities. We have applied the finite pore volume Guggenheim–Anderson–de Boer model (fpv-GAB) [35] as one of the few models that allow us to fit type IV adsorption isotherms according to the IUPAC classification and are based on the assumption of the realistic mechanism of the adsorption process. Such an approach allows the use of an isotherm description by the theoretical model to estimate the pore diameter. This is quite convenient since the adsorption data are often analyzed by fitting with the use of different theoretical equations of isotherms. The reliability of the proposed methodology has been verified using the series of N2 adsorption isotherms simulated inside carbon nanotubes (CNTs) of different diameters and the experimentally measured data for selected real materials of the cylindrical pores.

2. Materials and Methods

2.1. Mathematical Formalism of Adsorption Modeling

We have based our analysis on the assumption that the space available for adsorption inside a cylindrical pore can be unambiguously divided into two parts: that occupied by the monolayer (contact layer) of the adsorbed fluid, and the central part, where the multilayer adsorption is observed (see Figure 1a). The geometric ratio of volumes of these parts (x) is equal to:
x = V m o n o l a y e r V m u l t i l a y e r = D e f f 2 ( D e f f 2 λ ) 2 ( D e f f 2 λ ) 2
where Deff is the effective pore diameter and λ is the monolayer thickness. For N2 used as the adsorbate in the present work we have assumed λ = 0.35 nm [36]. Since the ratio x is unambiguously related to the pore size (see Figure 1b), knowing the x value, the effective pore size Deff can be deduced from Equation (1):
D e f f = 2 λ ( 1 + 1 + x + 1 x )
The division of the adsorption space into two parts is justified by the adsorption mechanism. The completion of the monolayer is usually easily identifiable on the adsorption isotherm as a change of the isotherm slope. The filling of the remaining space is indicated by the final plateau on the isotherm. This plateau occurs only if the size of the pores allows them to fill completely. Typically, such a mechanism is reflected by the type IV of the isotherm according to the IUPAC classification. Further considerations are limited to such cases as the complete pore filling by the adsorbed N2 and the type IV of its adsorption isotherm.
If we assume that the density of the adsorbate filling the pore is the same in each point of a cylindrical pore, the geometric ratio of the volumes in Equation (2) can be replaced by the ratio of adsorption capacities of the contact layer (a0) and the remaining pore space (asec,s):
x = a 0 a s e c , s
Both adsorption capacities can be deduced from the adsorption isotherm. The above notations, i.e., a0 and asec,s, have been introduced to ensure compatibility with further mathematical considerations.
Unfortunately, most of the equations used to describe type IV isotherms have the dual-mode form. Such models are just the combination of two equations [37]. Examples of such approaches are Dubinin–Izotova adsorption isotherm (bimodal Dubinin–Astakhov equation) [38] or bimodal CMMS equation proposed by Rutherford [39]. These models treat independently both stages, i.e., the monolayer formation and the filling of the remained space. The adsorption isotherm equations are the sums of two independent parts. Hence, these and other similar approaches do not reflect the mechanism of the adsorption process. The molecules forming the multilayer are bound on the upper surface of the monolayer, i.e., on the molecules adsorbed in the contact layer. The quantity of fluid adsorbed in the multilayer should be dependent on adsorption in the monolayer. One of the few theoretical models able to generate the type IV isotherm and correctly reflect the mechanism of the process is the recently proposed fpv-GAB model [35]. As in the formalism of the original GAB [40,41,42] or Brunauer–Emmett–Teller (BET) [43] equations the fpv-GAB approach assumes the formation of vertical complexes of the adsorbed molecules, but their size is restricted. This restriction is equivalent to the finite volume of adsorption space and allows the generation of, among others, type IV isotherms, while the basic GAB equation can describe only type I, II, and III isotherms.
In the fpv-GAB approach the quantity of fluid adsorbed in the contact layer (on the primary adsorption sites) is described by Langmuir equation:
a p r i m = a 0 K h 1 + K h
where a0 is the concentration of primary centers equivalent to the monolayer capacity in the studied case, K is Langmuir constant reflecting the strength of solid–fluid interactions, and h = p/ps denotes the relative pressure, i.e., the fluid pressure in the gaseous phase (p) normalized with respect to the saturated vapor pressure ps at the given temperature. The formation of further layers is modeled using generalized Barton (GB) equation [35]:
a s e c = c ( a p r i m + a s e c ) ( 1 k a s e c n ) h
where asec is the amount adsorbed in the multilayers next to the contact layer, i.e., on the secondary sites, aprim is the pressure-dependent adsorption amount in the monolayer calculated from Equation (4), c is a constant related to the energetics of adsorption in this space, i.e., fluid–fluid interactions, and n is the fitting parameter that ensures the best fit of experimental isotherm by Equation (5) in the high-pressure region. The constant k is related to the above-mentioned adsorption capacities of the remaining pore space (asec,s), i.e., the adsorption amount in the multilayer at the saturation pressure. Since Equation (5) should lead to asec = asec,s for h = 1, the constant k can be expressed as [35]:
k = c ( a p r i m , s + a s e c , s ) a s e c , s c a s e c , s n ( a p r i m , s + a s e c , s )
where aprim,s is the adsorption amount on primary sites for h = 1:
a p r i m , s = a 0 K 1 + K
The total amount adsorbed is the sum of monolayer and multilayer contributions:
a = a p r i m + a s e c
The fpv-GAB model is defined by the combination of Equations (4)–(8).
A shortcoming of the fpv-GAB approach can be the lack of lateral interactions between molecules adsorbed in the contact layer. This component of the interaction is particularly important when the adsorption occurs on homogeneous surfaces. Our previous studies [32,44,45] showed that the fit of the adsorption isotherms in the monolayer region for such systems is unsatisfactory without considering the effects of horizontal interactions in this layer. One of the simplest ways to include such interactions in the modeling of adsorption in the contact layers is to replace the Langmuir formula (Equation (4)) by Fowler–Guggenheim (FG) equation [46]:
h = Θ 1 K ( 1 Θ 1 ) exp [ A Θ 1 ]
where K is the constant related to the strength of solid–fluid interactions, analogically to the same denoted constant in Equation (4), A is the parameter related to the strength of fluid–fluid lateral interactions, and Θ1 is the relative adsorption amount in the contact layer:
Θ 1 = a p r i m a 0
The model, obtained by the modification of the original fpv-GAB approach (Equations (4)–(8)) and replacing Langmuir equation by the FG model, can be denoted as fpv-GAB-li, i.e., the fpv-GAB model with lateral interactions in the first layer. It is defined by the combination of Equations (4)–(6) and (8)–(10). Here, aprim,s in Equation (6) is determined from Equations (9) and (10) for h = 1. For A = 0, when the effects connected with lateral interactions are neglected, the fpv-GAB-li approach simplifies to the basic fpv-GAB one. It should be noted that the fpv-GAB-li approach has been proposed for the first time in the current work.
The mathematical formalism of both models, i.e., fpv-GAB and fpv-GAB-li, causes that calculation of the adsorption amount for the given value of relative pressure requires the use of numerical methods to determine the contribution from the multilayer (Equation (5)) in the case of both models and also the contribution from the monolayer in the case of fpv-GAB-li (Equations (9) and (10)). In the current study the bisection method has been used.
Appearing in Equation (3), the constants a0 and asec,s as the best-fit parameters represent adsorption capacities in monolayer and multilayer. In the current study, the values calculated according to Equation (3) have been used to determine the pore diameter from Equation (2).

2.2. Monte Carlo Simulations of N2 Adsorption Isotherms inside CNTs

To verify the proposed methodology, the N2 adsorption has been simulated inside ideal, infinite, single walled carbon nanotubes (CNTs) of different effective diameters. We selected 30 zigzag type CNTs (see the first column of Tables S1 or S2 in the Supplementary Materials), with the diameters ranging from 1.06 nm to 12.10 nm. The diameter of the narrowest nanotube, Deff = 1.06 nm, is about three times greater than the monolayer thickness. This is a minimum pore size that enables us to observe two-step adsorption—first the formation of a monolayer close to the nanotube wall, followed by the filling of the remaining pore space. The adsorption isotherms have been simulated at T = 77.3 K (N2 boiling point) and relative pressure (p/ps) in the range 1.0 × 10−8–1.0 using hyper parallel tempering Monte Carlo (HPTMC) technique [47]. The values of the activity have been calculated as for an ideal gas. The cylindrical simulation boxes have consisted in fragments of CNT with the length of 4.23 nm. Periodic boundary conditions have been applied in the direction parallel to the tube axes. The formalism of the interaction energy calculations has been analogical to that described in [44]. The interaction parameters have been taken from [48,49].

2.3. Experimental Data

Nine experimental isotherms of N2 adsorption in cylindrical pores, measured at N2 boiling point, have been taken from the literature. They report the N2 adsorption inside carbon nanohorns (CNHs) [50] and in two series of MCM-41 materials [51,52]. Despite the fact that two parts, i.e., the conical and cylindrical ones, can be divided in the structure of CNHs [53], the internal pores of these materials are often modeled as tubular, i.e., CNHs are approximated by CNTs [54]. The presence of the conical part can be neglected because its volume is negligible. Previously, it was shown by us, via simulation, that there are no significant differences in the shapes of the low-temperature N2 isotherms inside infinite CNH and CNT of the same diameter, when the length of CNH is close to that observed for the real CNHs [55]. The MCM-41 type materials have cylindrical geometry of pores and a monomodal distribution of the pore diameters confirmed by XRD, TEM, and adsorption measurements [14,56]. The first series of the analyzed MCM-41 consists of three aluminosilicates samples modified with Al [51]. They were denoted according to the Si/Al molar ratios of the gels used for their preparation (60, 30, and 15) as Al-MCM-41(60), Al-MCM-41(30) and Al-MCM-41(15), respectively, although the final Si/Al ratios in the obtained samples are equal to 36.8, 19.9 and 12.2, respectively [51]. The second series consists of five pure-silica MCM-41 materials prepared using alkyltrimethylammonium halides with different length of alkyl group (C12, C16 and C18) [52]. The samples prepared using a surfactant containing hexadecyl group (C16) as a template were synthesized by three different methods [52]. Therefore, the samples in the second series have been denoted as MCM-41-12, MCM-41-16A, MCM-41-16B, MCM-41-16C and MCM-41-18. In all the experimental isotherms we have neglected the experimental points for the relative pressure above 0.9. The significant increase of the adsorption amount observed in this range results mainly from the adsorption in macropores and on the external surfaces of the real samples, not analyzed here.

2.4. Fitting of Adsorption Isotherms

The genetic algorithm proposed by Storn and Price [57] has been applied to fit simulated and experimental isotherms by the fpv-GAB (Equations (4)–(8)) and fpv-GAB-li (Equations (5), (6) and (8)–(10)) models. The best-fit parameters have been: a0, K, c, n, and asec,s (fpv-GAB) and a0, K, A, c, n, and asec,s (fpv-GAB-li). The fit quality has been evaluated using the determination coefficient (DC):
D C = 1 i ( a t h e o , i a s i m / e x p , i ) 2 i ( a s i m / e x p , i a ¯ s i m / e x p ) 2
where atheo,i and asim/exp,i are, respectively, the adsorbed amounts predicted by the model and the observed amounts (simulated or experimental) for at the i-th point of the isotherm, and āsim/exp is the value of the simulated/experimental adsorption averaged over the whole isotherm. DC is equal to 1 for the perfect fit and decreases when the quality of the fit decreases.

3. Results and Discussion

We begin the discussion of the accuracy of the proposed methodology by the analysis of simulated isotherms of N2 adsorption inside CNTs, for which the effective pore diameters are perfectly known. Figure 2 shows selected simulated isotherms and their fits by fpv-GAB and fpv-GAB-li models. Tables S1 and S2 in Supplementary Materials present the obtained values of the best-fit parameters for all the model CNTs studied here. Both models allow good fit of the simulated data with the determination coefficient DC > 0.99. Quite systematic (close to monotonic) changes in the values of the best-fit parameters reflect the evolution of the shape of the isotherms as the diameter of the CNT increases. When the nanotubes become wider, the decreasing surface curvature causes the reduction of the effective energy of the solid–fluid interactions and the pressure related to monolayer formation increases. In consequence the values of K parameter in both models decrease, as well as the A value in the fpv-GAB-li approach. Similarly, the pressure of the filling of the central part of the pore also increases. The shape of the theoretical isotherm in this region is determined mainly by the values of the c and n parameters. The changes of both values are opposite, i.e., c decreases and n increases. It should be noted that the effects of both parameters can be formally compensating to some extent. The most important differences are observed in the region of the multilayer adsorption. Since the pore volume increases considerably as the diameter increases, the values of asec,s also rise systematically. Surprisingly, the adsorption capacities of the monolayer (the a0 parameter) also increase. Such behavior was analyzed and explained by Salmas and Androutsopoulos [58]: due to geometric restrictions, the uptake of molecules on the highly curved cylindrical surface is much less efficient than on the flat surface. This effect reduces as the pore diameter increases.
The careful analysis of the theoretical curves in Figure 2 indicates that the fpv-GAB-li model fits the simulated isotherms better than fpv-GAB model (the values of DC are significantly higher). The main differences are observed in the low-pressure range, where the contact layer is formed. As the fpv-GAB-li model includes the lateral interactions between adsorbate molecules in this layer, it reproduces the steplike character of the monolayer formation better than the fpv-GAB model.
Figure 3a shows the diameters of the model CNTs calculated using Equation (2) and the fitted monolayer and multilayer capacities (Tables S1 and S2 in Supplementary Materials). In general, both isotherm models lead to a realistic estimation of the pore diameters. The absolute errors (with the exception of the widest CNTs) do not exceed ca. 0.4 nm and 0.6 nm when using the best-fit parameters from the fpv-GAB and fpv-GAB-li models, respectively (Figure 3b).
The higher absolute errors for the widest nanotubes presumably result from imprecise estimation of the adsorption capacity in the multilayer. The isotherms for these systems contain only a few points in the final plateau. The relatively high errors for the narrow CNTs can be explained in two ways. On the one hand, the previously mentioned geometric restrictions [58] cause the noneffective packing of molecules in the contact layer. In such a case, the assumption of the uniform adsorbate density in the whole pore volume is not fulfilled and, according to Equation (2), the underestimation of the monolayer volume leads to the overestimation of the pore diameter. On the other hand, the density of the contact layer may slightly increase after its completion as the gas pressure approaches the saturation value [59,60,61,62]. The applied mathematical models do not account for such an effect. As is shown in Figure S1 in the Supplementary Materials, the monolayer capacity, i.e., the value of a0 parameter, corresponds with the adsorption amount at quite low value of relative pressure, while the shape of the isotherm suggests that some further increase in the adsorption amount for this layer may occur.
Surprisingly, although the fpv-GAB-li model fits better the adsorption isotherms, that is the fpv-GAB model that gives higher accuracy of the pore diameter determination. In fact, the evaluated contact layer capacity a0 is slightly lower in fpv-GAB-li model (Tables S1 and S2 in Supplementary Materials), which in turn leads to a slightly higher overestimation of the pore diameter than in fpv-GAB model. Nevertheless, the relative errors (Figure 3c) of Deff assessment are small and do not exceed 10 % for the nanotubes of diameter higher than ca. 3.5 nm. In the range 4.5–10 nm the errors are even less than 5%. Calculations of these quantities (Δ) include the division by the large geometric diameters. Thus, the relative errors are relatively small for the wider CNTs.
The verification of our model against the simulated isotherms justifies its application to analyze the experimental data. Figure 4 shows the results of the fitting of the experimental N2 adsorption isotherms for the selected studied systems. The results for the remaining systems and all fitting parameters are given in Supplementary Materials (Figures S2 and S3, Tables S3 and S4). As in all cases, the optimized value of A in fpv-GAB-li model is practically equal to 0; this approach simplifies to the basic fpv-GAB one (the values of the other best-fit parameters in both models become the same). As a consequence, only one theoretical curve is shown in each figure. We suggest that the apparent lack of the influence of lateral interaction on the adjustment of the experimental isotherms is a consequence of the heterogeneity of the pore walls in the real samples. If the energetic heterogeneity of the solid–fluid potential is large, the adsorption process begins on the sites of the highest energy, regardless the (weak) lateral fluid–fluid interactions. As the gas pressure increases, less adsorbing sites begin to be filled. The adsorbed amount increases progressively and the sharp step corresponding to the monolayer completion is not observed. Such behavior was observed for many the real adsorbents—see for example [32,44,45]. Nevertheless, the quality of the fits for the isotherms studied here is quite good, as indicated by the high DC values collected in Tables S3 and S4 in Supplementary Materials. The theoretical model reflects well the shape of the isotherms. The only exception is the medium-pressure range of the isotherms of the MCM-41-12 to MCM-41-18 materials. Here, the theoretical model does not precisely reproduce the experimental points at the beginning of the multilayer formation pressure range.
Table 1 gives the pore diameters calculated using Equation (2) and the adsorption capacities of the monolayer and the multilayer predicted by the applied models. For comparison, the diameters estimated using more sophisticated methods [32,50,51,52,63] are also shown. In general, the pore sizes calculated using our methodology are very similar to those reported in the literature. For example, in the case of CNHs, the estimated diameter agrees with the results of the GCMC simulations of Ohba et al. [50] and with the pore size estimated previously by us using the layer-by-layer NT-PFG adsorption isotherm model [32]. This justifies our statement that the proposed methodology is a simpler alternative to the previous complex NT-PFG model. As in [32], the calculated diameter is also significantly lower than the one estimated using the IDBdB method [63]. However, some results, such as observations with the use of high-resolution transmission electron microscopy [50], suggest that the real diameter is much lower in comparison to the predictions of the last method. In case of the MCM-41 materials, the diameters calculated with Equation (2) are in good agreement with those reported in the literature [51,52]. The differences do not exceed 0.3 nm. With the exception of the samples with the widest pores, i.e., MCM-41-16B and MCM-41-16C, the proposed approach overestimates the pore size similarly to in the case of the model CNTs—Figure 3. However, the trend of the diameter changes in both considered series of MCM-41 materials is reflected quite well. Therefore, the proposed methodology can be used as the simple tool for the rapid estimation of cylindrical pore diameters and their comparison between different samples.

4. Conclusions

We proposed a simple approach which allows the estimation of the diameter of cylindrical pores with the use of a N2 adsorption isotherm. It was based on the geometrical considerations and the division of adsorption space into two parts: the monolayer and the multilayer space. The ratio of the volumes of these two spaces is unambiguously related to the pore diameter. We have suggested to use the fpv-GAB adsorption isotherm model to fit the adsorption data. The ratio of the adsorption capacities in the monolayer and the multilayer obtained as the best-fit parameters is equivalent to the volume ratio. We have verified the proposed methodology with the use of the series of isotherms simulated inside CNTs and experimentally measured for different adsorbents. It has been proven that this approach not only describes adsorption isotherms well, but also it provides an estimation of the pore diameter that is comparable with the geometrical values or the results of the application of more sophisticated methods. With the exception of narrower pores (Deff < 3 nm), the relative error does not exceed a few percent. Moreover, our method reflects the trend of the changes in the sizes of the considered series of materials and can be used to compare different samples. Therefore, the fpv-GAB isotherm model can be treated as a simple tool which not only enables analysis of the shape of the isotherm theoretically, but also provides information about the pore diameter.

Supplementary Materials

The following are available online at Table S1: The values of best-fit parameters obtained using fpv-GAB model for the description of all the simulated N2 adsorption isotherms inside CNTs; Table S2: As in Table S2 but for fpv-GAB-li model; Table S3: The values of best-fit parameters obtained using fpv-GAB model for the description of all the studied experimental N2 adsorption isotherms; Table S4: As in Table S3 but for fpv-GAB-li model; Figure S1: The comparison of the simulated N2 adsorption isotherm inside (28,0) CNT and its shape predicted by the fpv-GAB and fpv-GAB-li models and the contributions from the monolayer predicted by the models and the obtained values of adsorption capacities in this layer (the value of a0 parameter); Figure S2: The results of the fitting of experimental N2 adsorption isotherms by fpv-GAB and fpv-GAB-li models for the samples Al-MCM-41(60) and Al-MCM-41(15); Figure S3: As in Figure S2 but for samples MCM-41-12, MCM-41-16A, MCM-41-16C and MCM-41-18.

Author Contributions

Conceptualization, S.F.; methodology, S.F., P.A.G. and R.M.; software, S.F.; validation, S.F., P.A.G. and R.M.; formal analysis, S.F., P.A.G., M.L., A.B.-K. and P.K.; investigation, S.F. and P.A.G.; data curation, S.F., P.A.G., M.L.; writing—original draft preparation, S.F., P.A.G., M.L. and P.K.; writing—review and editing, P.A.G., M.L. and A.B.-K.; visualization, S.F. and P.A.G.; supervision, S.F. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

The authors extent their sincere appreciation and thanks to Lucyna Firlej (Université de Montpellier, France) for proof-reading the manuscript, insightful discussion and valuable comments. S.F., P.A.G. and A.B.-K. acknowledge the use of the computer cluster at Poznań Supercomputing and Networking Centre (Poznań, Poland).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Thomas, S.; Pasquini, D.; Leu, S.-Y.; Gopakumar, D.A. (Eds.) Nanoscale Materials in Water Purification; Elsevier: Amsterdam, The Netherlands, 2019. [Google Scholar]
  2. Kaneko, K.; Rodríguez-Reinoso, F. (Eds.) Nanoporous Materials for Gas. Storage; Springer: Singapore, 2019. [Google Scholar]
  3. Li, F.; Bashir, S.; Liu, J.L. (Eds.) Nanostructured Materials for Next-Generation Energy Storage and Conversion. Fuel Cells; Springer: Berlin, Germany, 2018. [Google Scholar]
  4. Bar-Cohen, Y. (Ed.) Advances in Manufacturing and Processing of Materials and Structures; CRC Press: Boca Raton, FL, USA, 2019. [Google Scholar]
  5. Grumezescu, A.M. (Ed.) Nanomaterials for Drug Delivery and Therapy; Elsevier: Amsterdam, The Netherlands, 2019. [Google Scholar]
  6. Wang, S. Ordered mesoporous materials for drug delivery. Microporous Mesoporous Mater. 2009, 117, 1–9. [Google Scholar] [CrossRef]
  7. Panczyk, T.; Warzocha, T.P. Monte Carlo study of the properties of a carbon nanotube functionalized by magnetic nanoparticles. J. Phys. Chem. C 2009, 113, 19155–19160. [Google Scholar] [CrossRef]
  8. Xiao, F.-S.; Wang, L. (Eds.) Nanoporous Catalysts for Biomass Conversion; Wiley: New York, NY, USA, 2017. [Google Scholar]
  9. Bhushan, N. (Ed.) Springer Handbook of Nanotechnology; Springer: Berlin, Germany, 2010. [Google Scholar]
  10. Vajtai, R. (Ed.) Springer Handbook of Nanomaterials; Springer: Berlin, Germmany, 2013. [Google Scholar]
  11. Iijima, S. Helical microtubules of graphitic carbon. Nature 1991, 354, 56–58. [Google Scholar] [CrossRef]
  12. Harris, P.J.F. Carbon Nanotube Science: Synthesis, Properties and Applications; Cambridge University Press: Cambridge, UK, 2009. [Google Scholar]
  13. Sakintuna, B.; Yürüm, Y. Templated porous carbons:  a review article. Ind. Eng. Chem. Res. 2005, 44, 2893–2902. [Google Scholar] [CrossRef]
  14. Kresge, C.T.; Leonowicz, M.E.; Roth, W.J.; Vartuli, J.C.; Beck, J.S. Ordered mesoporous molecular sieves synthesized by a liquid-crystal template mechanism. Nature 1992, 359, 710–712. [Google Scholar] [CrossRef]
  15. Bennett, J.M.; Cohen, J.P.; Flanigen, E.M.; Pluth, J.J.; Smith, J.V. Crystal-structure of tetrapropylammonium hydroxide-aluminum phosphate number-5. ACS Symp. Ser. 1983, 218, 109–118. [Google Scholar]
  16. Zhao, D.; Huo, Q.; Feng, J.; Chmelka, B.F.; Stucky, G.D. Nonionic triblock and star diblock copolymer and oligomeric surfactant syntheses of highly ordered, hydrothermally stable, mesoporous silica structures. J. Am. Chem. Soc. 1998, 120, 6024–6036. [Google Scholar] [CrossRef]
  17. Losic, D.; Santos, A. (Eds.) Nanoporous Alumina. Fabrication, Structure, Properties and Applications; Springer: Heidenberg, Germany, 2015. [Google Scholar]
  18. Guldi, D.M.; Martín, N. (Eds.) Carbon Nanotubes and Related Structures. Synthesis, Characterization, Functionalization, and Applications; Wiley: Weinheim, Germany, 2010. [Google Scholar]
  19. Zhao, D.; Wan, Y.; Zhou, W. Ordered Mesoporous Materials; Wiley: Weinheim, Germany, 2013. [Google Scholar]
  20. Bhattacharyya, S.; Lelong, G.; Saboungi, M.-L. Recent progress in the synthesis and selected applications of MCM-41: A short review. J. Exp. Nanosci. 2006, 1, 375–395. [Google Scholar] [CrossRef]
  21. Kalbasi, R.J.; Izadi, E. Hydrothermal synthesis of pure AlPO4-5 without fluoride medium: Synthesis, characterization and application as a support. J. Porous Mater. 2013, 20, 547–556. [Google Scholar] [CrossRef]
  22. Singh, S.; Kumar, R.; Setiabudi, H.D.; Nanda, S.; Vo, D.-V.N. Advanced synthesis strategies of mesoporous SBA-15 supported catalysts for catalytic reforming applications: A state-of-the-art review. Appl. Catal. A 2018, 559, 57–74. [Google Scholar] [CrossRef]
  23. Srinivasan, N.R.; Bandyopadhyaya, R. Highly accessible SnO2 nanoparticle embedded SBA-15 mesoporous silica as a superior photocatalyst. Microporous Mesoporous Mater. 2012, 149, 166–171. [Google Scholar] [CrossRef]
  24. Karnjanakom, S.; Guan, G.; Asep, B.; Hao, X.; Kongparakul, S.; Samart, C.; Abudula, A. Catalytic upgrading of bio-oil over Cu/MCM-41 and Cu/KIT-6 prepared by β-cyclodextrin-assisted coimpregnation method. J. Phys. Chem. C 2016, 120, 3396–3407. [Google Scholar] [CrossRef]
  25. Nagpure, S.; Browning, J.F.; Rankin, S.E. Incorporating poly(3-hexyl thiophene) into orthogonally aligned cylindrical nanopores of titania for optoelectronics. Microporous Mesoporous Mater. 2017, 240, 65–72. [Google Scholar] [CrossRef] [Green Version]
  26. Wu, M.; Ma, T. Recent progress of counter electrode catalysts in dye-sensitized solar cells. J. Phys. Chem. C 2014, 118, 16727–16742. [Google Scholar] [CrossRef]
  27. De Souza, L.V.; Tkachenko, O.; Cardoso, B.N.; Pizzolato, T.M.; Dias, S.L.P.; Vasconcellos, M.A.Z.; Arenas, L.T.; Costa, T.M.H.; Moro, C.C.; Benvenutti, E.V. Strategy to control the amount of titania dispersed on SBA-15 surface preserving its porosity, aiming to develop a sensor for electrochemical evaluation of antibiotics. Microporous Mesoporous Mater. 2019, 287, 203–210. [Google Scholar] [CrossRef]
  28. Sing, K.S.W.; Everett, D.H.; Haul, R.A.W.; Moscou, L.; Pieroti, R.A.; Rouquérol, J.; Siemieniewska, T. Reporting physisorption data for gas/solid systems with special reference to the determination of surface area and porosity. Pure Appl. Chem. 1985, 57, 603–619. [Google Scholar] [CrossRef]
  29. Thommes, M.; Kaneko, K.; Neimark, A.V.; Olivier, J.P.; Rodriguez-Reinoso, F.; Rouquérol, J.; Sing, K.S.W. Physisorption of gases, with special reference to the evaluation of surface area and pore size distribution (IUPAC Technical Report). Pure Appl. Chem. 2015, 87, 1051–1069. [Google Scholar] [CrossRef] [Green Version]
  30. Rouquérol, J.; Rouquérol, F.; Sing, K.S.W.; Llewellyn, P.; Maurin, G. Adsorption by Powders and Porous Solids. Principles, Methodology and Applications; Academic Press: Amsterdam, The Netherlands, 2014. [Google Scholar]
  31. Bandosz, T.J.; Biggs, M.J.; Gubbins, K.E.; Hattori, Y.; Iiyama, T.; Kaneko, K.; Pikunic, J.; Thomson, K.T. Molecular Models of Porous Carbons. In Chemistry and Physics of Carbon; Radovic, L.R., Ed.; Marcel Dekker Inc.: New York, NY, USA, 2003; Volume 28, pp. 41–228. [Google Scholar]
  32. Furmaniak, S.; Terzyk, A.P.; Gauden, P.A.; Rychlicki, G. Simple model of adsorption in nanotubes. J. Colloid Interface Sci. 2006, 295, 310–317. [Google Scholar] [CrossRef]
  33. Gu, C.; Gao, G.-H.; Yu, Y.-X. Density functional study of hydrogen adsorption at low temperatures. J. Chem. Phys. 2003, 119, 448–495. [Google Scholar] [CrossRef] [Green Version]
  34. Wu, X.-J.; Ning, Z.-F.; Cheng, Z.-L.; Wang, Q.; Qi, R.-R.; Huang, L.; Zhang, W.-T. Simplified local density model for gas adsorption in cylindrical carbon pores. Appl. Surf. Sci. 2019, 491, 335–349. [Google Scholar] [CrossRef]
  35. Furmaniak, S. Influence of activated carbon porosity and surface oxygen functionalities’ presence on adsorption of acetonitrile as a simple polar volatile organic compound. Environ. Technol. 2015, 36, 1984–1999. [Google Scholar] [CrossRef]
  36. Hagymassy, J., Jr.; Brunauer, S.; Mikhail, R.S. Pore structure analysis by water vapor adsorption: I. t-curves for water vapour. J. Colloid Interface Sci. 1969, 29, 485–491. [Google Scholar] [CrossRef]
  37. Inglezakis, V.J.; Poulopoulos, S.G.; Kazemian, H. Insights into the S-shaped sorption isotherms and their dimensionless forms. Microporous Mesoporous Mater. 2018, 272, 166–176. [Google Scholar] [CrossRef]
  38. Izotova, T.I.; Dubinin, M.M. Microporous structure of activated carbons. Zhurnal Fizicheskoi Khimii 1965, 39, 2769–2803. [Google Scholar]
  39. Rutherford, S.W. Modeling Water Adsorption in Carbon Micropores:  Study of Water in Carbon Molecular Sieves. Langmuir 2006, 22, 702–708. [Google Scholar] [CrossRef] [PubMed]
  40. Guggenheim, E.A. Application of Statistical Mechanics; Clarendon Press: Oxford, UK, 1966. [Google Scholar]
  41. Anderson, R.B. Modification of the Brunauer, Emmett and Teller equations. J. Am. Chem. Soc. 1946, 68, 686–691. [Google Scholar] [CrossRef]
  42. De Boer, J.M. The Dynamic Character of Adsorption; Clarendon Press: Oxford, UK, 1953. [Google Scholar]
  43. Brunauer, S.; Emmett, P.H.; Teller, E. Adsorption of gases in multimolecular layers. J. Am. Chem. Soc. 1938, 60, 309–319. [Google Scholar] [CrossRef]
  44. Furmaniak, S.; Terzyk, A.P.; Gauden, P.A.; Harris, P.J.F.; Wiśniewski, M.; Kowalczyk, P. Simple model of adsorption on external surface of carbon nanotubes—A new analytical approach basing on molecular simulation data. Adsorption 2010, 16, 197–213. [Google Scholar] [CrossRef] [Green Version]
  45. Furmaniak, S.; Terzyk, A.P.; Rychlicki, G.; Wiśniewski, M.; Gauden, P.A.; Kowalczyk, P.; Werengowska, K.M.; Dulska, K. The system: Carbon tetrachloride—closed carbon nanotubes analysed by a combination of molecular simulations, analytical modelling and adsorption calorimetry. J. Colloid Interface Sci. 2010, 349, 321–330. [Google Scholar] [CrossRef]
  46. Fowler, R.H.; Guggenheim, E.A. Statistical Thermodynamics. The Theory of the Properties of Matter in Equilibrium; Cambridge University Press: Cambridge, UK, 1965. [Google Scholar]
  47. Yan, Q.; de Pablo, J.J. Hyper-parallel tempering Monte Carlo: Application to the Lennard-Jones fluid and the restricted primitive model. J. Chem. Phys. 1999, 111, 9509–9516. [Google Scholar] [CrossRef] [Green Version]
  48. Potoff, J.J.; Siepmann, L.I. Vapor-liquid equilibria of mixtures containing alkanes, carbon dioxide, and nitrogen. AIChE J. 2001, 47, 1676–1682. [Google Scholar] [CrossRef]
  49. Steele, W.A. The Interaction of Gases with Solid Surfaces; Pergamon: New York, NY, USA, 1974. [Google Scholar]
  50. Ohba, T.; Murata, K.; Kaneko, K.; Steele, W.A.; Kokai, F.; Takahashi, K.; Kasuya, D.; Yudasaka, M.; Iijima, S. N2 adsorption in an internal nanopore space of single-walled carbon nanohorn:  GCMC simulation and experiment. Nano Lett. 2001, 1, 371–373. [Google Scholar] [CrossRef]
  51. Rozwadowski, M.; Leżańska, M.; Erdmann, K. Mesoporous molecular sieves modified with carbonaceous deposits. Adsorption 2005, 11, 363–377. [Google Scholar] [CrossRef]
  52. Kruk, M.; Jaroniec, M.; Kim, J.M.; Ryoo, R. Characterization of highly ordered MCM-41 silicas using X-ray diffraction and nitrogen adsorption. Langmuir 1999, 15, 5279–5284. [Google Scholar] [CrossRef]
  53. Ajima, K.; Yudasaka, M.; Murakami, T.; Maigné, A.; Shiba, K.; Iijima, S. Carbon nanohorns as anticancer drug carriers. Molec. Pharm. 2005, 2, 475–480. [Google Scholar] [CrossRef]
  54. Tanaka, H.; Kanoh, H.; Yudasaka, M.; Iijima, S.; Kaneko, K. Quantum effects on hydrogen isotope adsorption on single-wall carbon nanohorns. J. Am. Chem. Soc. 2005, 127, 7511–7516. [Google Scholar] [CrossRef]
  55. Furmaniak, S.; Terzyk, A.P.; Kaneko, K.; Gauden, P.A.; Kowalczyk, P.; Itoh, T. The first atomistic modelling-aided reproduction of morphologically defective single walled carbon nanohorns. Phys. Chem. Chem. Phys. 2013, 15, 1232–1240. [Google Scholar] [CrossRef]
  56. Beck, J.S.; Vartuli, J.C.; Roth, W.J.; Leonowicz, M.T.; Kresge, C.; Schmitt, K.; Chu, C.; Olson, O.; Sheppard, E.; McCullen, S.; et al. A new family of mesoporous molecular sieves prepared with liquid crystal templates. J. Am. Chem. Soc. 1992, 114, 10834–10843. [Google Scholar] [CrossRef]
  57. Storn, R.; Price, K. Differential evolution—A simple and efficient heuristic for global optimization over continuous spaces. J. Global Optim. 1997, 11, 341–359. [Google Scholar] [CrossRef]
  58. Salmas, C.E.; Androutsopoulos, G.P. Rigid sphere molecular model enables an assessment of the pore curvature effect upon realistic evaluations of surface areas of mesoporous and microporous materials. Langmuir 2005, 21, 11146–11160. [Google Scholar] [CrossRef]
  59. Furmaniak, S.; Terzyk, A.P.; Gauden, P.A.; Wesołowski, R.P.; Kowalczyk, P. Ar, CCl4, and C6H6 adsorption outside and inside of the bundles of multi-walled carbon nanotubes—Simulation study. Phys. Chem. Chem. Phys. 2009, 11, 4982–4995. [Google Scholar] [CrossRef] [PubMed]
  60. Gauden, P.A.; Furmaniak, S.; Włoch, J.; Terzyk, A.P.; Zieliński, W.; Kowalczyk, P.; Kurzawa, J. The influence of geometric heterogeneity of closed carbon nanotube bundles on benzene adsorption from the gaseous phase—Monte Carlo simulations. Adsorption 2016, 22, 639–651. [Google Scholar] [CrossRef] [Green Version]
  61. Zeng, Y.; Horio, K.; Horikawa, T.; Nakai, K.; Do, D.D.; Nicholson, D. On the evolution of the heat spike in the isosteric heat versus loading for argon adsorption on graphite—A new molecular model for graphite & reconciliation between experiment and computer simulation. Carbon 2017, 122, 622–634. [Google Scholar]
  62. Xu, H.; Zeng, Y.; Do, D.D.; Nicholson, D. On the nonwetting/wetting behavior of carbon dioxide on graphite. J. Phys. Chem. C 2019, 123, 9112–9120. [Google Scholar] [CrossRef]
  63. Kowalczyk, P.; Hołyst, R.; Tanaka, H.; Kaneko, K. Distribution of carbon nanotube sizes from adsorption measurements and computer simulation. J. Phys. Chem. B 2005, 109, 14659–14666. [Google Scholar] [CrossRef]
Figure 1. (a) The scheme of the division of cylindrical adsorption space into the monolayer and the multilayer. Deff is the effective diameter and λ is the monolayer thickness; (b) the ratio of monolayer and multilayer volumes plotted as the function of the effective diameter. The diameters are reduced by the monolayer thickness.
Figure 1. (a) The scheme of the division of cylindrical adsorption space into the monolayer and the multilayer. Deff is the effective diameter and λ is the monolayer thickness; (b) the ratio of monolayer and multilayer volumes plotted as the function of the effective diameter. The diameters are reduced by the monolayer thickness.
Molecules 26 01509 g001
Figure 2. Simulated isotherms of N2 adsorption inside selected carbon nanotubes (CNTs) (points) and their fits by the models: fpv-GAB (Equations (4)–(8))—blue lines and fpv-GAB-li (Equations (5), (6) and (8)–(10))—red lines. Only selected simulated points are shown for clarity. The inserts in lower panels show the low-pressure part of isotherms, related to the formation of N2 monolayer on the pore wall.
Figure 2. Simulated isotherms of N2 adsorption inside selected carbon nanotubes (CNTs) (points) and their fits by the models: fpv-GAB (Equations (4)–(8))—blue lines and fpv-GAB-li (Equations (5), (6) and (8)–(10))—red lines. Only selected simulated points are shown for clarity. The inserts in lower panels show the low-pressure part of isotherms, related to the formation of N2 monolayer on the pore wall.
Molecules 26 01509 g002
Figure 3. (a) Effective pore diameters Deff,model calculated using Equation (2) and fpv-GAB (blue circles) or fpv-GAB-li (red triangles) fits of simulated isotherms of N2 adsorption in CNTs. The diameters are plotted as the function of the geometrical diameter of CNTs (Deff,geo). The dashed line represents Deff,model = Deff,geo; (b,c) the absolute (b) and relative (c) differences between Deff,model and Deff,geo. The horizontal dashed lines correspond to Deff,model = Deff,geo.
Figure 3. (a) Effective pore diameters Deff,model calculated using Equation (2) and fpv-GAB (blue circles) or fpv-GAB-li (red triangles) fits of simulated isotherms of N2 adsorption in CNTs. The diameters are plotted as the function of the geometrical diameter of CNTs (Deff,geo). The dashed line represents Deff,model = Deff,geo; (b,c) the absolute (b) and relative (c) differences between Deff,model and Deff,geo. The horizontal dashed lines correspond to Deff,model = Deff,geo.
Molecules 26 01509 g003
Figure 4. The results of the fitting of experimental N2 adsorption isotherms by fpv-GAB (Equations (4)–(8)) and fpv-GAB-li (Equations (4)–(6) and (8)–(10)) models for the samples CNHs, Al-MCM-41(30) and MCM-41-16B. The insets show the same data in logarithmic scale of the relative pressure. The points represent experimental data and lines reflect the predictions of the models. Since fpv-GAB-li equation is simplified to fpv-GAB (A ≈ 0) only one theoretical line is plotted for each system. In the case of CNHs and MCM-41-16B only selected experimental points are shown for clarity.
Figure 4. The results of the fitting of experimental N2 adsorption isotherms by fpv-GAB (Equations (4)–(8)) and fpv-GAB-li (Equations (4)–(6) and (8)–(10)) models for the samples CNHs, Al-MCM-41(30) and MCM-41-16B. The insets show the same data in logarithmic scale of the relative pressure. The points represent experimental data and lines reflect the predictions of the models. Since fpv-GAB-li equation is simplified to fpv-GAB (A ≈ 0) only one theoretical line is plotted for each system. In the case of CNHs and MCM-41-16B only selected experimental points are shown for clarity.
Molecules 26 01509 g004
Table 1. The comparison of effective pore diameters for all the studied samples calculated from Equation (2) and the best-fit parameters obtained from the fitting of N2 adsorption isotherms by theoretical models and the values obtained from more sophisticated methods published by others.
Table 1. The comparison of effective pore diameters for all the studied samples calculated from Equation (2) and the best-fit parameters obtained from the fitting of N2 adsorption isotherms by theoretical models and the values obtained from more sophisticated methods published by others.
SampleDeff [nm]Method 1References
The Models
(Current Work)
Literature Data
CNHs2.832.9
3.38 ± 0.48
2.96
GCMC
IDBdB
NT-PFG
[50]
[63]
[32]
Al-MCM-41(60)3.223.10BJH[51]
Al-MCM-41(30)3.563.54BJH[51]
Al-MCM-41(15)3.743.56BJH[51]
MCM-41-123.373.07
3.08
BJH
GM
[52]
MCM-41-16A4.023.85
3.88
BJH
GM
[52]
MCM-41-16B4.474.66
4.55
BJH
GM
[52]
MCM-41-16C5.075.12
5.11
BJH
GM
[52]
MCM-41-184.364.16
4.22
BJH
GM
[52]
1 GCMC—the grand canonical Monte Carlo simulations, IDBdB—the improved Derjaguin–Broekhoff–de Boer method, NT-PFG—the nanotube polymolecular Fowler–Guggenheim adsorption isotherm, BJH—the Barrett–Joyner–Halenda method, GM—the geometrical considerations assuming hexagonal arrays of uniform pores of circular geometry.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Furmaniak, S.; Gauden, P.A.; Leżańska, M.; Miśkiewicz, R.; Błajet-Kosicka, A.; Kowalczyk, P. The Finite Pore Volume GAB Adsorption Isotherm Model as a Simple Tool to Estimate a Diameter of Cylindrical Nanopores. Molecules 2021, 26, 1509. https://doi.org/10.3390/molecules26061509

AMA Style

Furmaniak S, Gauden PA, Leżańska M, Miśkiewicz R, Błajet-Kosicka A, Kowalczyk P. The Finite Pore Volume GAB Adsorption Isotherm Model as a Simple Tool to Estimate a Diameter of Cylindrical Nanopores. Molecules. 2021; 26(6):1509. https://doi.org/10.3390/molecules26061509

Chicago/Turabian Style

Furmaniak, Sylwester, Piotr A. Gauden, Maria Leżańska, Radosław Miśkiewicz, Anna Błajet-Kosicka, and Piotr Kowalczyk. 2021. "The Finite Pore Volume GAB Adsorption Isotherm Model as a Simple Tool to Estimate a Diameter of Cylindrical Nanopores" Molecules 26, no. 6: 1509. https://doi.org/10.3390/molecules26061509

APA Style

Furmaniak, S., Gauden, P. A., Leżańska, M., Miśkiewicz, R., Błajet-Kosicka, A., & Kowalczyk, P. (2021). The Finite Pore Volume GAB Adsorption Isotherm Model as a Simple Tool to Estimate a Diameter of Cylindrical Nanopores. Molecules, 26(6), 1509. https://doi.org/10.3390/molecules26061509

Article Metrics

Back to TopTop