Experimental Confirmation of Different Mechanisms of Evaporation

Nov 1, 2002 - Adsorption Hysteresis of Nitrogen and Argon in Pore Networks and Characterization of Novel Micro- and Mesoporous Silicas. Matthias ...
0 downloads 9 Views 170KB Size
9830

Langmuir 2002, 18, 9830-9837

Experimental Confirmation of Different Mechanisms of Evaporation from Ink-Bottle Type Pores: Equilibrium, Pore Blocking, and Cavitation Peter I. Ravikovitch and Alexander V. Neimark* The Center for Modeling and Characterization of Nanoporous Materials, TRI/Princeton, 601 Prospect Avenue, Princeton, New Jersey 08540-0625 Received June 27, 2002. In Final Form: September 24, 2002 We present a new classification of physical mechanisms of adsorption hysteresis in ink-bottle type pores confirmed by experimental studies of capillary condensation of N2, Ar, and Kr at 77.4 and 87.4 K in ordered 3D cage-like structures of FDU-1 and SBA-16 nanoporous templated silicas with narrow pore size distribution. An analysis of the hysteresis loops and scanning desorption isotherms on FDU-1 with ∼15 nm spherical pores reveals three mechanisms of evaporation: (1) evaporation from the blocked cavities controlled by the size of connecting pores (classical ink-bottle or pore blocking effect); (2) spontaneous evaporation caused by cavitation of the stretched metastable liquid; and, for the first time, (3) nearequilibrium evaporation in the region of hysteresis from unblocked cavities that have access to the vapor phase. Studies of the temperature dependence of the hysteresis loop showed a transition between the cavitation and the pore blocking regimes of evaporation, which opens up the possibility to tune the experimental conditions in order to characterize the pore neck size distribution. A similar transition was also observed for Kr in SBA-16 silica with ∼8.5 nm pores. The pressure of cavitation was found to depend on the pore geometry. This conclusion questions the conventional assumption that the pressure of cavitation, which determines the lower closure point of the hysteresis loop, is a unique function of the adsorbate and the temperature. The experimental results agree with the nonlocal density functional theory (NLDFT) of capillary condensation hysteresis in spherical cavities reported earlier (Ravikovitch, P. I.; Neimark, A. V. Langmuir 2002, 18, 1550).

1. Introduction Following the classical works of Kraemer,1 McBain,2 and Cohan,3 the phenomenon of capillary hysteresis in pores of varying cross section and pore networks is usually explained by the pore blocking effect.4,5 The pore blocking effect is most pronounced in ink-bottle type pores, the simplest image of which is a spheroidal cavity fitted with a narrow cylindrical neck. A qualitative explanation of the capillary hysteresis in the ink-bottle pore is as follows. The pore filling in the process of condensation follows the formation of a liquid film on the cavity wall and thus is controlled by the radius of curvature of the cavity, rc. The pressure of condensation pc is given by the Kelvin-Laplace equation, pc/p0 ) exp(-2γvl/RTrc), where p0 is the saturation vapor pressure at the given temperature T, γ is the liquid-vapor surface tension, vl is the liquid molar volume, R is the universal gas constant. The evaporation in the process of desorption occurs after the formation of a hemispherical meniscus in the pore neck and thus is controlled by the radius of the neck, rn. The pressure of desorption is determined by the Kelvin-Laplace equation for the equilibrium meniscus in the cylindrical neck pd/p0 ) exp(-2γvl/RTrn). This concept has been elaborated later by Broekhoff and de Boor.6 Similarly to the mechanism of filling of * Corresponding author. E-mail: [email protected]. (1) Kraemer, E. O. Treatise on Physical Chemistry; Taylor, H. S., Ed.; Van Nostrand: New York, 1931; p 1661. (2) McBain, J. W. J. Am. Chem. Soc. 1935, 57, 699. (3) Cohan, L. H. J. Am. Chem. Soc. 1938, 60, 433. (4) Everett, D. H. The Solid-Gas Interface; Flood, E. A., Ed.; Dekker: New York, 1967; Vol. 2. (5) Gregg, S. J.; Sing, K. S. W. Adsorption, Surface Area and Porosity; Academic Press: New York, 1982. (6) Broekhoff, J. C. P.; de Boer, J. H. J. Catal. 1968, 10, 153.

cylindrical pores,7 they associated the condition of condensation in the cavity with the achievement of the limit of stability of metastable adsorption films on the spherical wall. Broekhoff and de Boor used the Derjaguin principle of the additive action of capillary and surface forces8 to determine the condition of stability of adsorbed films and also the condition of equilibrium between the adsorbed film and the condensed liquid for cylindrical and spherical pores. The equations they derived for the pressure of spontaneous condensation and the pressure of equilibrium are referred to as the Derjaguin-Broekhoff-de Boer (DBdB) equations.6 For many years, the DBdB equations for the limit of stability of adsorption films on the spherical wall and for the equilibrium meniscus in the cylindrical neck represented the best quantitative description of the capillary condensation hysteresis in ink-bottle type pores. The concept of pore blocking in ink-bottle type pores was successfully employed to describe the triangular hysteresis loops of type H2 by IUPAC classification typical for porous glasses.4 The pore structure of these materials is disordered and represents a network of channels with prominent enlargements and constrictions as shown by high-resolution microscopy and scattering data9a and by an off-lattice reconstruction method.9b Gelb and Gubbins10 mimicked this structure by molecular simulations of the process of formation of porous glasses. To account for networking effects, in the early 80s, Wall and Brown,11a (7) Broekhoff, J. C. P.; de Boer, J. H. J. Catal. 1967, 9, 8. (8) Derjaguin, B. V. Acta Phys. Chem. URSS 1940,12, 181. (9) (a) Levitz, P. E.; Ehret, G.; Sinha, S. K.; Drake, J. M. J. Chem. Phys. 1991, 95, 6151. (b) Pellenq, R.-J. M.; Rousseau, B.; Levitz, P. E. Phys. Chem. Chem. Phys. 2001, 3, 1207. (10) Gelb, L. D.; Gubbins, K. E. Langmuir 1998, 14, 2097. (11) (a) Wall, G. C.; Brown, R. J. C. J. Colloid Interface Sci. 1981, 82, 141. (b) Mason, G. J. Colloid Interface Sci. 1982, 88, 36.

10.1021/la026140z CCC: $22.00 © 2002 American Chemical Society Published on Web 11/01/2002

Mechanisms of Evaporation from Ink-Bottle Type Pores

Mason,11b and Neimark12 and later Seaton13 developed the percolation theory of capillary condensation and desorption in the network of cavities with multiple necks, assuming that the condensation in a cavity is controlled by its size while the desorption from a cavity is controlled by the smallest neck in a chain of necks connecting the cavity with the vapor phase. In particular, Neimark12 showed that the behavior of scanning isotherm of Xe on porous glass4 follows the predictions of the percolation theory and developed a method for calculations of the neck size distribution. Recently, the problem of adsorption hysteresis in inkbottle type pores and the role of pore network connectivity attracted renewed attention.9b,14,15 Using molecular dynamics simulations Sarkisov and Monson15a have questioned the existance of pore blocking effects in pores of simple geometry. Monte Carlo simulations of adsorption in model disordered systems have shown that pore blocking effects9b and network connectivity14 have little influence on the shape of the hysteresis loop. Kierlik et al.16 and Woo et al.15b developed lattice DFT models of sorption in disordered media and constructed, without invoking any networking effects, the hysteresis loops similar in shape to the experimental adsorption isotherms on porous glasses. Advances in experimental investigations were related to the development of new nanoporous materials with an ordered structure of pores, most of which were synthesized using regular surfactant and micellar mesophases as templates.17-20 The pores in these materials resemble with the greatest accuracy the simple geometrical models of cylindrical channels17 or spherical cavities18-20 which can be investigated theoretically in great detail. That is, the ordered nanoporous materials provide a unique opportunity to compare theoretical predictions and experimental data without invoking adjustable parameters. In the preceding paper,21 we developed a nonlocal density functional theory (NLDFT) model for capillary (12) (a) Neimark, A. V. Rep. Acad. Sci. USSR (Doklady Akademii Nauk SSSR) 1983, 273, 384. (b) Neimark, A. V. Russian J. Phys. Chem. (Z. Fiz. Khim.) 1986, 60, 1745. (c) Stud. Surf. Sci. Catal. 1991, 62, 67. (13) Seaton, N. A. Chem. Eng. Sci. 1991, 46, 1895. (14) Gelb, L. D.; Gubbins, K. E. Fundamentals of Adsorption 7; Kaneko, K., Kanoh, H., Hanzawa, Y., Eds.; IK International: ChibaCity, 2002; p 333. (15) (a) Sarkisov, L. D.; Monson, P. A. Langmuir 2001, 17, 7600. (b) Woo, H.-J.; Sarkisov, L. D.; Monson, P. A. Langmuir 2001, 17, 7472. (16) (a) Kierlik, E.; Rosinberg, M. L.; Tarjus, G.; Viot, P. Phys. Chem. Chem. Phys. 2001, 3, 1201-1206. (b) Kierlik, E.; Monson, P. A.; Rosinberg, M. L.; Sarkisov, L.; Tarjus, G. Phys. Rev. Lett. 2001, 87, 055701. (17) (a) Kresge, C. T.; Leonowicz, M. E.; Roth, W. J.; Vartuli, J. C.; Beck, J. S. Nature 1992, 359, 710. (b) Zhao, D. Y.; Feng, J. L.; Huo, Q. S.; Melosh, N.; Fredrickson, G. H.; Chmelka, B. F.; Stucky, G. D. Science 1998, 279, 548. (c) Inagaki, S.; Guan, S.; Ohsuna, T.; Terasaki, O. Nature 2002, 416, 34. (d) Van Der Voort, P.; Ravikovitch, P. I.; De Jong, K. P.; Benjelloun, M.; Van Bavel, E.; Janssen, A. H.; Neimark, A. V.; Weckhuysen, B. M.; Vansant, E. F. J. Phys. Chem. B 2002, 106, 5873. (18) (a) Yu, C.; Yu, Y.; Zhao, D. Chem. Commun. 2000, 575. (b) Yu, C.; Yu, Y.; Miao, L.; Zhao, D. Microporous Mesoporous Mater. 2001, 44-45, 65. (c) Matos, J. R.; Mercuri, L. P.; Kruk, M.; Jaroniec, M. Langmuir 2002, 18, 884. (d) Matos, J. R.; Kruk, M.; Mercuri, L. P.; Jaroniec, M.; Asefa, T.; Coombs, N.; Ozin, G. A.; Kamiyama, T.; Terasaki, O. Chem. Mater. 2002, 14, 1903. (19) (a) Guan, S.; Inagaki, S.; Ohsuna, T.; Terasaki, O. J. Am. Chem. Soc. 2000, 122, 5660. (b) Sakamoto, Y.; Dı´az, I.; Terasaki, O.; Zhao, D.; Pe´rez-Pariente, J.; Kim, J. M.; Stucky, G. D. J. Phys. Chem. B 2002, 106, 3118. (c) Lukens, W. W.; Yang, P.; Stucky, G. D. Chem. Mater. 2001, 13, 28. (d) Lukens, W. W.; Stucky, G. D. Chem. Mater. 2002, 14, 1665. (e) Kim, J. M.; Stucky, G. D. Chem. Commun. 2000, 1159. (20) (a) Zhao, D. Y.; Huo, Q. S.; Feng, J. L.; Chmelka, B. F.; Stucky, G. D. J. Am. Chem. Soc. 1998, 120, 6024. (b) Sakamoto, Y. H.; Kaneda, M.; Terasaki, O.; Zhao, D. Y.; Kim, J. M.; Stucky, G.; Shim, H. J.; Ryoo, R. Nature 2000, 408, 449. (c) Van Der Voort, P.; Benjelloun, M.; Vansant, E. F. J. Phys. Chem. B 2002, 106, 9027. (21) Ravikovitch, P. I.; Neimark, A. V. Langmuir 2002, 18, 1550.

Langmuir, Vol. 18, No. 25, 2002 9831

Figure 1. NLDFT N2 adsorption-desorption isotherm at 77 K in 15.5 nm spherical cavity (thick solid line). The theoretical pressures of liquidlike spinodal, equilibrium, and vaporlike spinodal are denoted, respectively, as psl, pe, and psv. Experimentally observed pressures of spontaneous capillary condensation and cavitation are denoted as pc and pcav, respectively (pc ≈ psv and psl < pcav < pe). Three possible regimes of evaporation from the cavity are classified with respect to the pressure pd of desorption from the neighboring pores or from the pore neck (see text in Discussion section).

condensation in spherical cavities. The theory allows us to construct stable and metastable equilibrium states, which possess a spherical symmetry compliant with the cavity geometry. The isotherm of these states forms a hysteresis loop comprising two branches, the adsorption branch of low density, or vaporlike states, and the desorption branch of high density, or liquidlike states (Figure 1). The pressure of the equilibrium transition pe separates the stable and metastable states. At p < pe, the adsorption branch corresponds to the stable adsorption films while the desorption branch corresponds to the metastable stretched liquid. The desorption branch terminates at the liquidlike spinodal, psl, the limit of mechanical stability of stretched condensed liquid. At p > pe, the desorption branch corresponds to the stable condensed liquid while the adsorption branch corresponds to the metastable adsorption film. The adsorption branch terminates at the vaporlike spinodal, psv, which represents the limit of mechanical stability of the metastable adsorption film. In general, the pressures of condensation and desorption are controlled by the nucleation barriers at given experimental conditions22 and cannot be determined from the NLDFT calculations of thermodynamic equilibrium states. Thus, the theoretical spinodal pressures give the upper and lower bounds for the experimental hysteresis loop. A comparison of the theoretical results with the reference experiments allowed us to make several conclusions about the mechanisms of hysteresis in inkbottle type pores. The results of NLDFT calculations were found to be in quantitative agreement with experimental data on adsorption of nitrogen at 77.4 K in templated nanoporous silica with ordered structure of cage-like pores.21 In sufficiently large pores, which exhibit wide hysteresis loops, the condensation was shown to occur near the vaporlike spinodal, in accord with the classical scenario of spontaneous condensation (condensation pressure is (22) Vishnyakov, A.; Neimark, A. V. J. Phys. Chem. B 2001, 105, 7009.

9832

Langmuir, Vol. 18, No. 25, 2002

denoted as pc in Figure 1). In the further discussion, this condensation mechanism is referred to as the nearspinodal condensation. The NLDFT predictions for the equilibrium and vaporlike spinodal pressures asymptotically merge with the DBdB equations6 for cavity diameters exceeding ca. 10 nm. At the same time, we argued that there are two sorption mechanisms which are not accounted for in the DBdB theory. First, the DBdB theory rules out the possibility of equilibrium capillary condensation and/or desorption in spherical cavities. We have shown that in sufficiently small cavities (for example for materials reported in refs 19a,19b) the condensation occurs reversibly without a hysteresis and, most importantly, the pressure of condensation corresponds to the equilibrium pressure predicted by the NLDFT.21 With respect to sorption of nitrogen at 77.4 K, this mechanism of equilibrium condensation and desorption is expected in cavities with the diameter between ∼3 and ∼ 6 nm. In such strong confinements, the energy barriers for nucleation of new phases are so small that at the standard conditions of adsorption experiments the metastable states are short lived and both the condensation and the desorption occur near equilibrium. Second, the classical picture of hysteresis in ink-bottle pores implies that the desorption process is controlled exclusively by the diameter of the pore neck. We have shown that the pore blocking effect takes place only when the neck diameter is greater than a certain characteristic size (ca. 4 nm for nitrogen).21 Desorption from a cavity fitted with narrower necks occurs at the conditions of cavitation of the metastable stretched fluid and does not appear to depend on the neck diameter. The pressure of cavitation, pcav, is determined by the onset of nucleation at given experimental conditions and may appreciably exceed the liquidlike spinodal pressure, psl. (Figure 1). In this case, the whole adsorption-desorption cycle in the cavity is controlled by its geometry, and the isotherm would be unchanged even if all the necks were closed yet still permeable to vapor. Probably, the absence of pore blocking effects was observed in molecular simulations14,15a at these particular conditions. In the further discussion, the mechanism in which the desorption is controlled by the size of pore necks is referred to as the pore blocking controlled desorption, in accord with the conventional picture of pore blocking effects. The desorption mechanism of spontaneous evaporation of the stretched metastable liquid is referred to as cavitation. Starting from the original idea put forward by Schofield23 and elaborated in the works of Kadlec and Dubinin24 and Burgess and Everett,25 the latter mechanism was discussed in the literature based on the limiting tensile strength hypothesis to explain the lower closure point of adsorption hysteresis loops, which is known to be insensitive to the nature of the adsorbent. For sorption of nitrogen at 77.4 K, it is well documented that the lower closure point of adsorption hysteresis varies in a narrow range of relative pressures, from ∼0.42 to ∼0.5.5 In this paper we present an elaborated picture of possible pathways of adsorption hysteresis in ink-bottle type pores along with the experimental evidence for three different mechanisms of evaporation, (Figure 1). These mechanisms are: (i) desorption caused by cavitation of the metastable adsorbate, (ii) pore blocking controlled desorption, and also, for the first time, (iii) near-equilibrium desorption in the region of hysteresis from the cavities that have (23) Schofield, R. K. Discuss. Faraday Soc. 1948, 3, 105. (24) Kadlec, O.; Dubinin, M. M. J. Coll. Interface Sci. 1969, 31, 479. (25) Burgess, C. G. V.; Everett, D. H. J. Colloid Interface Sci. 1970, 33, 611.

Ravikovitch and Neimark

immediate access to the vapor phase and thus are effectively unblocked. We show that the temperature decrease may lead to transition from the cavitation to the pore blocking regime of desorption. This observation opens up an opportunity to tune the experimental conditions in order to characterize the pore neck size distribution. 2. Experimental Results The isotherms of capillary condensation and desorption were measured using Autosorb-1C volumetric instrument (Quantachrome) on well-characterized samples of ordered nanoporous materials with cagelike pore structures. Nitrogen and argon isotherms at 77.4 and 87.4 K were measured on a sample of FDU-1.18a,18b FDU-1 possesses a regular periodic structure with cagelike pores of ∼15 nm in diameter, as calculated by the NLDFT method.21 The argon isotherm at 77.4 K and krypton isotherm at 87.4 K were measured on a sample of SBA-16 silica20c (∼8.5 nm cagelike pores by the NLDFT method). SBA-16 has an ordered Im3m structure recently confirmed by electron crystallography.20b FDU-1 and SBA-16 represent the best currently available examples of ordered ink-bottle pore systems. Figures 2a-c show adsorption isotherms of N2 and Ar on FDU-1. Repeated measurements are denoted by different symbols. High stability of the material was confirmed by the repeated adsorption measurements performed after the sample was stored at ambient conditions for a period of two years. No change in the isotherm was detected except for the ∼6% reduction in the pore volume and surface area due to aging. Ar isotherms at 77 K are plotted using the saturation pressure of the supercooled liquid Ar (∼230 Torr) rather than the experimentally observed saturation pressure of bulk solid Ar (∼205 Torr). This choice of the saturation pressure gives a consistent comparison between the Ar isotherms at 77 and at 87 K and was shown to give consistent results when the experimental data at 77 K are compared with the simulated and NLDFT isotherms.26 For all three experimental conditions, the adsorption isotherms exhibit the following characteristic features (see Figure 2). The initial portion of the isotherm corresponds to adsorption in micropores and formation of a polymolecular film on the pore walls. FDU-1 contains intrawall micropores (their amount Niw in shown in Figure 2a) that are typical for the materials prepared using polymeric surfactants.27,28 The sharp capillary condensation step between points A and B indicates a relatively narrow pore size distribution. At point B all mesopores are filled and the uptake at higher pressures corresponds to adsorption in macropores, presumably spaces between particles. Desorption is considerably delayed and essentially does not occur until the point F (for N2 at 77 K and Ar at 87 K) or point E (for Ar at 77 K). N2 and Ar isotherms at 87 K behave alike and exhibit a sharp evaporation step FG at the relative pressures of ∼0.5 and ∼0.44, respectively. However, starting from point E, the Ar isotherm at 77 K is qualitatively different from the others. It exhibits a two-stage evaporation regime: a sharp step FG is preceded by a more gradual part EF. An important feature of the desorption branches is the presence of a small substep CD, which is more clearly seen on Ar isotherms. This substep transforms into a well(26) Neimark, A. V.; Ravikovitch, P. I.; Grun, M.; Schuth, F.; Unger, K. K. J. Colloid Interface Sci. 1998, 207, 159. (27) Kruk, M.; Jaroniec, M.; Ko, C. H.; Ryoo, R. Chem. Mater. 2000, 12, 1961. (28) Ravikovitch, P. I.; Neimark, A. V. J. Phys. Chem. B 2001, 105, 6817.

Mechanisms of Evaporation from Ink-Bottle Type Pores

Langmuir, Vol. 18, No. 25, 2002 9833

Figure 2. Adsorption isotherms on FDU-1 silica: (a) N2 at 77.4 K, thick solid line denotes NLDFT isotherm in 15.5 nm spherical pore shifted upward by Niw to account for adsorption in intrawall pores; (b) Ar at 87.3 K; (c) Ar at 77.4 K; (d) adsorption isotherms on SBA-16 [Ar at 77.4 K (open squares) and Kr at 87.3 K (solid circles)] plotted using the saturation pressures of the bulk solid Kr. For (a) and (b), vertical dotted lines indicate positions of the equilibrium (left) and adsorption spinodal transitions (right). Symbols denote different experiments.

pronounced step displayed on scanning desorption isotherms, thus indicating a distinctive mechanism of desorption which will be discussed in detail below. The scanning desorption isotherms were measured starting from the partially saturated sample (indicated as point B1 for one of the scanning curves; the others are similar). The position of the first step C1D1 on the scanning desorption curves B1C1D1F1G (B1C1D1E1F1G for Ar at 77 K) coincides with the position of the substep CD on the main desorption branch. As the pressure is reduced further, the scanning isotherms repeat all the features of the main desorption curves, namely the region of gradual desorption is followed by the one-stage evaporation F1G for N2 and Ar at 87 K, and the two-stage evaporation for Ar at 77 K (with steps, E1F1 and F1G). Figure 2d presents Ar and Kr adsorption-desorption isotherms on SBA-16 silica. The Ar isotherm at 77 K exhibits a single evaporation step, much like N2 and Ar (87 K) isotherms on FDU-1. However, the Kr isotherm on SBA-16 at 87 K exhibits a two-stage evaporation similar to the Ar isotherm on FDU-1 at 77 K. The substep on the main desorption branch has not been detected for SBA16.

3. Discussion We observe several new features of adsorption-desorption hysteresis loops, which are well-defined in ordered cagelike structures and are smeared out in disordered materials: (1) there is a substep on the initial plateau region of the main desorption branch; (2) the scanning desorption isotherms have two distinct steps and are qualitatively different from the smooth scanning curves reported earlier for H1 and H2 type hysteresis loops;4,5 (3) as the temperature decreases from 87 K to 77 K, the abrupt step on the desorption isotherm of Ar is split into two steps. These characteristic features of adsorption-desorption isotherms can be explained on the basis of the nonlocal density functional theory NLDFT of capillary condensation in spherical cavities.21 The comparison of the experimental isotherm with the model NLDFT isotherm calculated for a spherical cavity of 15.5 nm in diameter (Figure 2a) shows that the capillary condensation proceeds in accord with theoretical predictions. The calculated NLDFT isotherm, which was shifted upward by Niw to account for adsorption in intrawall pores (∼0.3 cm3/g) is in good agreement with the experimental

9834

Langmuir, Vol. 18, No. 25, 2002

data. As it is generally accepted for such wide mesopores, the polymolecular adsorption film becomes metastable6,7,29 and condensation is thus delayed until the free energy barrier for nucleation of the liquidlike phase is overcome by fluctuations.22 The capillary condensation occurs practically at the limit of stability of polymolecular adsorption film (vaporlike spinodal point) predicted theoretically,21,30-31 where the energy barrier for nucleation vanishes.22 Thus, the step AB (Figures 2a-c) corresponds to the near-spinodal condensation. Using NLDFT adsorption isotherms in spherical pores, we calculated the size distribution of mesopores by the method described earlier.21 The distribution was centered at 15 nm with the standard deviation of ∼1 nm. Both for N2 and Ar in 15.5 nm spherical pore, the calculated relative pressure of the thermodynamic equilibrium between the vaporlike and liquidlike phases (vertical dotted lines in Figures 2a and 2b) matches closely the position of the substep CD on the main desorption isotherm and the step C1D1 on the scanning desorption isotherms (NLDFT calculations for Ar were performed with the parameters reported earlier for siliceous materials of MCM-41 type26). We therefore conjecture that the steps CD and C1D1 correspond to the equilibrium evaporation from the cavities. Until recently, the possibility that the liquid condensed in spheroidal cavities may evaporate at the equilibrium conditions has not been discussed in the literature. In the classical description of adsorption in ink-bottle pores6 it has been explicitly assumed that the evaporation is controlled by the size of the pore necks. As related to the equilibrium desorption, the authors assumed that “there is no means of establishing an equilibrium path of desorption in the case of spheroidal cavities”.6 As was noted in the Introduction, in the preceding paper21 we have shown that in highly ordered materials with sufficiently small cagelike pores (∼3-6 nm in diameter) the reversible capillary condensation step corresponds to the equilibrium transition, as calculated by the NLDFT for spherical pores. The validity of this interpretation has been confirmed by the fact that the pore diameters predicted by NLDFT were in excellent agreement with the pore sizes obtained from geometrical models of these structures.21 In sufficiently narrow pores the free energy barrier for nucleation is relatively low while the fluctuations are strong, which leads to the regime of reversible capillary condensation/ evaporation.32 Wider spherical pores can also exhibit near-equilibrium evaporation. In the blocked cavities larger than ∼6 nm, the barrier for the liquid-vapor nucleation (bubble formation) is very high. Therefore, desorption does not occur until either the liquid evaporates from the neighboring pores, which block access to the vapor phase (classical pore blocking effect), or spontaneous cavitation of the metastable liquid trapped in the cavity occurs (see possible scenarios in Figure 3). The equilibrium evaporation, however, can occur when the cavity has a direct contact with the vapor phase, so that the liquid-gas interface can be formed, and thus the cavity is effectively unblocked. This can happen when the cavity has a short neck or a wide window opening. We conjecture that the substep CD corresponds to the evaporation from unblocked (29) Saam, W. F.; Cole, M. W. Phys. Rev. B 1975, 11, 1088. (30) Neimark, A. V.; Ravikovitch, P. I.; Vishnyakov, A. Phys. Rev. E 2000, 62, R1493. (31) Neimark, A. V.; Ravikovitch, P. I. Microporous Mesoporous Mater. 2001, 44-45, 697. (32) Ravikovitch, P. I.; Domhnaill, S. C.O Ä .; Neimark, A. V.; Schu¨th, F.; Unger K. K. Langmuir 1995, 11, 4765.

Ravikovitch and Neimark

Figure 3. Schematic illustration of three different desorption mechanisms from ink-bottle type pores (see text).

cavities which have immediate access to the outer surface of the primary particles of FDU-1. The number of such cavities is small, and thus we observe just a small substep on the desorption branch. When the particle is partially saturated, as in the case of the scanning desorption curves, a larger number of cavities remains unblocked and has direct access to the vapor phase. Evaporation from these pores gives rise to the step C1D1 on the scanning desorption curves. The simplest situation of two adjacent overlapping cavities of slightly different sizes (R2 > R1) is shown in Figure 3. In this case the cavity R1 is filled while the cavity R2 has only a polymolecular adsorption film on its interior surface. Upon decrease in pressure, the cavity R1 will empty at the pressure that is close to the equilibrium conditions for a spherical pore of radius R1. Once all cavities that have access to the vapor phase empty, the desorption transits into the plateau region, which corresponds to the evaporation from the polymolecular films in the emptied pores and to the reduction of density of condensed liquid in the blocked cavities. From the scanning isotherm we can determine the fraction of cavities that were blocked in the beginning of desorption at point B1 and remained filled after the evaporation from unblocked pores had occurred. Indeed, in the plateau region, the scanning isotherm is determined by the fraction of filled cavities, R, and is given by

Nscan) (1 - R)Nads + RNdes In Figure 4, we demonstrate the validity of this relation for the scanning isotherm of Ar at 87 K. The fractions of filled pores were determined from the best fit to experimental data in the region indicated by arrows. Evaporation from the blocked pores along the main and scanning desorption isotherms occurs abruptly at the same pressure, which lies between the equilibrium pressure, pe, and the liquidlike spinodal pressure, psl (Figure 2). While in the classical models of independent ink-bottled pores1,2,6 and percolation models of pore networks,11-13 the onset of evaporation from the blocked pores is assumed to be controlled entirely by the sizes of pore necks, in the case of FDU-1 we observe spontaneous cavitation, which occurs upon achievement of the conditions of instability of condensed liquid. The cavitation pressure, pcav, is determined by the nucleation conditions and may exceed the liquidlike spinodal pressure, psl (Figure 1). Abrupt desorption steps at the same relative pressure (pcav/p0 ∼ 0.5 for N2 at 77 K) were observed for other cagelike structures with different pore sizes.18,20a A similar behavior was recently reported for materials consisting of open and blocked cylindrical pores.17d The spontaneous desorption occurs in the cavities with the necks narrower than the neck from which the evaporation occurs at the pressure of cavitation pcav/p0. The limiting size of necks that cannot be opened at pressures larger than the pressure of

Mechanisms of Evaporation from Ink-Bottle Type Pores

Langmuir, Vol. 18, No. 25, 2002 9835

Figure 4. Fit of Ar scanning desorption isotherms on FDU-1 at 87 and 77 K (lines) in the pressure interval indicated by arrows. At 87 K the fraction of blocked pores on the upper and lower scanning curves is 0.52 and 0.32, respectively. At 77 K, the fraction of blocked pores is 0.44.

cavitation pcav/p0 are estimated by the NLDFT model for the infinitely long cylindrical pores as ∼5 nm for N2 at 77 K and for Ar at 87 K, and ∼4 nm for Ar at 77 K. Since the neck sizes in the FDU-1 sample were smaller than 5 nm, the desorption of N2 at 77 K and Ar at 87 K occurred in the regime of cavitation. Figures 1 and 3 summarize three possible evaporation regimes from ink-bottle type pores. The conditions at which one or the other mechanism is operative depend on the pressure pd at which the neighboring pores empty (note that the pore neck itself is considered as the neighboring pore). Each cavity is characterized by two characteristic pressures-the pressure of equilibrium transition, pe, and the pressure of spontaneous cavitation, pcav. When the neighboring pores empty at the pressure lower than pcav, evaporation occurs via cavitation. When the neighboring pores empty at the intermediate pressure, between pcav and the equilibrium pressure pe, we deal with the classical pore blocking effect and such a cavity-neck assembly can be truly called an ink-bottle pore. When the neighboring pores empty at the pressure that exceeds pe, the cavity becomes effectively unblocked and will empty at the equilibrium pressure. Two possible examples, a cavity with a wide cylindrical neck and slightly overlapping spherical pores, are shown in Figure 3. To elucidate the two-stage behavior of the desorption isotherm of Ar at 77 K, one can compare this isotherm (Figure 2c) with the Ar isotherm on a sample of SBA-16 material20c (Figure 2d), which contains ∼8.5 nm wide cavities with the standard deviation of ∼0.5 nm, connected by necks 6-7 nm) ink-bottle pores with ∼4 nm neck (N2 at 77 K). NLDFT bulk spinodal line estimates the limit of stability of the stretched metastable liquid. The line of equilibrium evaporation from 12.5 σff cylindrical pore approximately describes evaporation from the pore neck.

by the different temperature dependencies of the conditions of cavitation and the conditions of evaporation from necks. The temperature dependencies of the chemical potentials of cavitation and evaporation are shown in Figure 5, which illustrates schematically two different mechanisms of evaporation from a large ink-bottle type pore. The NLDFT calculations were performed for N2. The chemical potential of cavitation was estimated from the spinodal of bulk liquid. The equilibrium chemical potential in the ∼4 nm infinitely long cylindrical pore was used to approximate the conditions of evaporation from the neck. Two scenarios of the process of desorption can be suggested as the chemical potential decreases from zero (saturation). At a relatively high (but subcritical) temperature, the liquid condensed in the cavity becomes unstable and evaporates before the conditions for the equilibrium evaporation from the necks are met. This cavitation regime is observed for N2 at 77 K and Ar at 87 K on FDU-1 and Ar at 77 K on SBA-16. Note that the pressure of cavitation does not appear to depend significantly on the size of cavities at a given temperature; the cavity diameter itself is irrelevant as long as it is sufficiently larger than the

9836

Langmuir, Vol. 18, No. 25, 2002

pore neck. Below the crossover temperature, which depends on the neck size, evaporation from the cavity occurs after the pore neck opens. Thus, the portions EF and E1F1 of Ar isotherms at 77 K (Figure 2c) correspond to the desorption controlled by pore blocking. The slope of these portions reflects the size distribution of necks (from ∼5 nm and below) and the pore network connectivity. Due to the distribution of pore necks below ∼4 nm, the desorption process is not fully completed upon the achievement of the conditions of cavitation at points F and F1 (which corresponds to the crossover point in Figure 5) and further decreases in the chemical potential leads to spontaneous evaporation reflected by sharp steps FG and F1G. A decrease of the pore neck diameter has an effect similar to an increase of the temperature: it leads to the transition from the pore blocking to the cavitation regime of evaporation. Thus, to characterize the neck size distribution, the temperature should be chosen low enough to provide the pore blocking regime of desorption. For Kr at 87 K, which is below its bulk triple point by 28.5 K, we still are able to observe the pore filling of SBA-16 (Figure 2d). On the desorption path, one can see a transition from the gradual pore blocking regime of desorption to the abrupt, spontaneous evaporation, which is qualitatively similar to that observed for Ar on FDU-1 at 77 K. Note that for SBA-16 we do not observe the substep related to the evaporation from the unblocked pores facing the outer surface of particles. This is probably due to the larger size of the primary particles in SBA-16, which makes this amount too small to be detected. It should be noted that Ar at 77.4 K and Kr at 87.3 K are below their bulk triple points by 6.5 and 28.5 K, respectively. Generally, solidification/melting transition for a fluid confined in silica pores is suppressed.33-35 In particular, Ar at 77 K behaves as a liquid in ∼3-5 nm cylindrical pores of MCM-41, and condensation/evaporation pressures are predicted rather accurately by the NLDFT, provided that the experimental pressures are reduced by the saturation pressure of the supercooled liquid Ar.26,36 What is remarkable is that the condensation and evaporation behavior for Ar at 77 K in ∼15 nm pores of FDU-1 and, especially, Kr at 87 K in ∼8.5 nm SBA-16 is similar to that for fluids above the triple point. In particular, when the experimental condensation and equilibrium pressures of Ar in FDU-1 are reduced by the saturation pressure of supercooled liquid Ar,26 the results are in agreement with the NLDFT calculations in spherical pores. However, we cannot totally exclude the possibility of solidification of Ar in FDU-1 or Kr in SBA-16 because the pores in these materials are larger than in MCM-41. Adsorption behavior of Ar in FDU-1 is consistent with the studies of Ar sorption both above and below the triple point in ∼10 nm (cylindrical model) Vycor and Gelsil glasses,35 which exhibit a capillary condensation/solidification step at about the same pressure. These studies concluded that the solid Ar condensate undergoes a firstorder transition very much like the liquid.35 It should also be noted that below the triple point there exists an upper (33) (a) Warnock, J.; Awschalom, D. D.; Shafer, M. W. Phys. Rev. Lett. 1986, 57, 1752. (b) Molz, E.; Wong, A. P. Y.; Chan, M. H. W.; Beamish, J. R. Phys. Rev. B 1993, 48, 5741. (34) Gelb, L. D.; Gubbins, K. E.; Radhakrishnan, R.; SliwinskaBartkowiak, M. Rep. Prog. Phys. 1999, 62, 1573. (35) (a) Huber, P.; Knorr, K. Phys. Rev. B 1999, 60, 12657. (b) Wallacher, D.; Knorr, K. Phys. Rev. B 2001, 63, 104202. (36) Ravikovitch, P. I.; Haller, G. L.; Neimark, A. V. Adv. Colloid Interface Sci. 1998, 76-77, 203.

Ravikovitch and Neimark

Figure 6. Adsorption isotherms of Ar at 77.4 K on ordered materials with cagelike and cylindrical pores: (a) 3.6 nm cylindrical MCM-41 (black circles); (b) 8.5 nm cages,