Argon and Nitrogen Adsorption in Disordered Nanoporous Carbons

Department of Chemical Engineering, North Carolina State University, 113 Riddick ... Raleigh, North Carolina 27695-7905, Laboratoire des Mate´riaux D...
0 downloads 0 Views 1MB Size
Langmuir 2005, 21, 4431-4440

4431

Argon and Nitrogen Adsorption in Disordered Nanoporous Carbons: Simulation and Experiment Jorge Pikunic,†,‡ Philip Llewellyn,§ Roland Pellenq,| and Keith E. Gubbins*,† Department of Chemical Engineering, North Carolina State University, 113 Riddick Labs, Raleigh, North Carolina 27695-7905, Laboratoire des Mate´ riaux Divise´ s, Reveˆ tements, Electroce´ ramiques, MADIREL, UMR 6121 Universite´ de Provence-CNRS, Centre de Saint Je´ roˆ me, 13397 Marseille, Cedex 20, France, and Centre de Recherche en Matie` re Condense´ e et Nanosciences, CNRS-UPR 7251, Campus de Luminy, case 913, 13288 Marseille, Cedex 09, France Received November 18, 2004. In Final Form: February 24, 2005 We report experimental measurements of the isosteric heats of adsorption for argon and nitrogen in two microporous saccharose-based carbons, using a Tian-Calvet microcalorimeter. These data are used to test recently developed molecular models of these carbons, obtained by a constrained reverse Monte Carlo method. grand canonical Monte Carlo simulation is used to calculate the adsorption isotherms and isosteric heats for these systems, and the results for the latter are compared to the experimental data. For argon, excellent quantitative agreement is obtained over the entire range of pore filling. In the case of nitrogen, very good agreement is obtained over the range of coverage 0.25 e Γ/Γ0 e 0.85, but discrepancies are observed at lower and higher coverages. The discrepancy at low coverage may be due to the presence of oxygenated groups on the pore surfaces, which are not taken into account in the model. The differences at high coverage are believed to arise from the presence of a few mesopores, which again are not included in the model. Pair correlation functions (argon-carbon and argon-argon) are determined from the simulations and are discussed as a function of pore filling. Snapshots of the simulations are presented and provide a picture of the pore filling process.

1. Introduction Realistic molecular models of porous carbons are needed to interpret experimental data used to characterize the pore structure of carbons, as well as to predict a wide range of properties (adsorption, heats of adsorption, diffusion rates, etc.) of confined phases within carbons. A variety of models have been proposed [for a review see ref 1]. The simplest and most widely used is the slit pore model. However, this model is oversimplified, since it fails to account for pore connectivity or the highly disordered pore structure of many industrial carbons. Recently, we have proposed a new reconstruction method2 based on a constrained reverse Monte Carlo procedure, for building realistic models of such carbons from experimental X-ray or neutron diffraction data. The method was shown to give good agreement with experimental structure data (small- and wide-angle X-ray diffraction, high resolution transmission electron microscopy) for highly disordered microporous carbons prepared from saccharose. In the work presented here, we address the question of whether such structural models are able to predict the energetics of a fluid adsorbed in the real material. To achieve this, we compare experiment and simulation results for heats of adsorption of argon at 77 K in two carbons produced by pyrolisis of saccharose at two different treatment temperatures. The two materials are termed * To whom correspondence should be addressed. † North Carolina State University. ‡ Current address: Department of Biochemistry, South Parks Road, University of Oxford, Oxford OX1 3QU, United Kingdom. § Universite ´ de Provence-CNRS. | CNRS-UPR 7251. (1) Bandosz, T. J.; Biggs, M. J.; Gubbins, K. E.; Hattori, Y.; Iiyama, T.; Kaneko, K.; Pikunic, J. P.; Thomson, K. T. Chem. Phys. Carbon 2003, 8, 41. (2) Pikunic, J.; Clinard, C.; Cohaut, N.; Gubbins, K. E.; Guet, J.-M.; Pellenq, R. J.-M.; Rannou, I.; Rouzaud, J.-N. Langmuir 2003, 19, 8565.

CS400 (heated at 400 °C) and CS1000 (heated at 1000 °C). Structural models of these materials have been developed and tested against experimental TEM images.2 The simulated images reproduce, at least qualitatively, the main features observed in the experimental images. Therefore, the obvious next step is to test the models by comparing adsorption predictions to experiment. The first step toward such a comparison is to design a suitable experimental system (adsorbate, conditions, and experiment). Since our models are intended to describe the microporosity of disordered porous carbons, it is desirable to use a small molecule that would be able to probe a higher fraction of the micropore volume. In this respect, nitrogen or argon are suitable candidates, and both have been widely used for the characterization of microporous materials using adsorption-based characterization methods.3,4 In our case, however, it is preferable to use argon as an adsorbate to test the carbon models for the following reasons: (i) Argon atoms are spherical and nonpolar, and they can be logically modeled as Lennard-Jones spheres. On the other hand, nitrogen molecules are diatomic and have a quadrupole; although they can be simplistically modeled as Lennard-Jones spheres, it is necessary to include the quadrupole and the anisotropy in shape to obtain quantitative agreement with experimental structure and thermodynamic data, even for bulk systems.5 Thus, argon-argon interactions can be more accurately described with simple potentials than nitrogen-nitrogen interactions. (3) Rouquerol, F.; Rouquerol, J.; Sing, K. In Adsorption by Powders and Porous Solids; Academic Press: London, 1999. (4) Pikunic, J. P.; Lastoskie, C. M.; Gubbins, K. E. In Handbook of Porous Solids; Schuth, F., Sing, K., Weitkamp, J., Eds.; Wiley-VCH: Weinheim, 2003; p 182. (5) Cheung, P. S. Y.; Powles, J. G. Mol. Phys. 1976, 32, 1383.

10.1021/la047165w CCC: $30.25 © 2005 American Chemical Society Published on Web 04/06/2005

4432

Langmuir, Vol. 21, No. 10, 2005

(ii) Although the carbons that we have modeled have only a small amount of heteroatoms (less than 8% of oxygen), these atoms are present as polar chemical groups (e.g., phenolic and carboxylic groups, ether bridges, etc.) that strongly interact with the nitrogen quadrupole. The potential energy resulting from these Coulombic forces may be quite strong compared to carbon atom-adsorbate interactions, especially at low relative pressures (or fractional filling). Consequently, it would be necessary to include the heteroatoms in the adsorbent model in order to make quantitative predictions of the behavior of nitrogen confined in porous carbons. On the other hand, since argon is nonpolar, the interaction of argon with chemical groups is weak, and any induction interactions are small compared to carbon atom-adsorbate dispersive interactions. This has been verified experimentally by comparing gas adsorption microcalorimetry results for argon and nitrogen in the same activated carbon samples.6 Thus, omitting the heteroatoms in the adsorbent models should not have a strong effect in the prediction of adsorption properties of argon. (iii) The carbon-adsorbate interaction potential parameters are usually derived from experiment by fitting the parameters to Henry’s constants (gas-solid virial coefficients) of the adsorbate on a nonporous graphitic surface. However, it has been pointed out that the interaction of nitrogen and other quadrupolar molecules (e.g., carbon dioxide) with carbon segments is highly anisotropic.7,8 The interaction of nitrogen with a carbon segment is largely determined by the position and orientation of the nitrogen molecule relative to the carbon segment. This casts some doubt on the validity of the use of experimentally determined potentials in molecular simulations of nitrogen adsorption in carbons that contain finite graphene segments. Although further investigation is required, it is believed that the main source of anisotropy is the interaction of the nitrogen quadrupole with the π electrons of the carbon segment. This problem is thus much less important in the case of argon, which is nonpolar and spherical. The enhancement in potential energy due to the confinement of a fluid in micropores can be measured experimentally by means of gas adsorption microcalorimetry. In these experiments, the isosteric heat of adsorption is measured as a function of the fractional pore filling. An important feature of gas adsorption microcalorimetry is its high sensitivity in the range of low relative pressures, which is exactly the pressure range at which micropore filling takes place.3 On the other hand, experimental isotherm measurements often lack accuracy at low relative pressures. Therefore, since we need to accurately probe the microporosity of these materials to compare to our model predictions, it is convenient to use gas adsorption microcalorimetry. In this work, we first assess the ability of the structural models of carbons to describe the texture of the real materials and to reproduce the energetics of confined argon. We then report grand canonical Monte Carlo simulations for three nitrogen models. The intention is to determine the importance of accurately describing the adsorbate-adsorbent intermolecular potential in predicting adsorption behavior. Finally, we present argon-argon and argon-carbon pair correlation functions at different fractional fillings (6) Fernandez-Colinas, J.; Denoyel, R.; Grillet, Y.; Rouquerol, F.; Rouquerol, J. Langmuir 1989, 5, 1205. (7) Meyer, E. F. J. Phys. Chem. 1967, 71, 4416. (8) Nicholson, D. Surf. Sci. 1987, 181, L189.

Pikunic et al.

characteristic of the main regions of the isosteric heat of adsorption curves. We use this structural information, along with the isosteric heats of adsorption and adsorption isotherms to describe the adsorption mechanism in disordered nanoporous carbons. 2. Experiment The differential enthalpies of adsorption, or isosteric heats of adsorption, were obtained directly using a Tian-Calvet type microcalorimeter coupled to a manometric apparatus.9 A quasiequilibrium system of gas introduction allows the determination of the high-resolution adsorption isotherms and the corresponding microcalorimetric recording. The microcalorimeter signal, φ, can be related to the differential enthalpy of adsorption, ∆adsh˙ , via the following relationship:

∆adsh˙ )

φ + Vb(∂p/∂t)T,A fa

where Vb is the volume of the microcalorimetric cell and ∂p/∂t is the increase in pressure, p, with time, t. The term Vb(∂p/∂t)T,A is the contribution due to the pressure expansion in the microcalorimetric cell. fa is the molar rate of adsorption. The term Vb(∂p/∂t)T,A becomes relatively negligible during the vertical regions of the isotherm (e.g., during micropore filling). If one also considers that in these regions all the adsorptive entering the system is adsorbed, then the following simple expression for the isosteric heat of adsorption can be obtained:

qst ) ∆adsh˙ ≈ φ/f where f is the constant flow of adsorptive introduced into the sample cell. Equilibrium tests are carried out either by varying the rate of gas introduction or by stopping the gas introduction at various points during the experiment. The saccharose-based carbons used for adsorption and microcalorimetry experiments in this work are the same that were used in diffraction, small angle scattering and TEM experiments earlier reported.2 Around 100 mg of sample was used for these experiments. Prior to each adsorption experiment, the sample was outgassed using sample controlled thermal analysis, SCTA.10,11 To prepare a clean surface with the most reproducibility, the sample was heated using SCTA conditions under a constant residual vacuum pressure of 0.02 mbar up to a final temperature of 150 °C. Under such conditions, 100 mg of sample requires 36 h to be heated to 150 °C. This final temperature was maintained until the residual pressure was less than 5 × 10-3 mbar. The argon and nitrogen used were obtained from Air Liquide (Alphagaz) and were of purity greater than 99.9995%.

3. Simulation Details We used the models of saccharose-based carbons presented in ref 2 to model the structure of CS400 and CS1000. These structural models were developed using a constrained reverse Monte Carlo procedure, and match the carbon-carbon pair correlation function of the real materials as determined from X-ray diffraction experiments. We modeled argon atoms as Lennard-Jones spheres. The potential energy due to the interaction of two argon atoms, as well as that due to the interaction of an argon atom with a carbon atom, is described with the LennardJones 6-12 potential. The interaction potential parameters are summarized in Table 1. Argon-carbon parameters were fitted to experimentally obtained Henry’s law constants for argon on graphitized carbon black.12 (9) Rouquerol, J. Thermochimie; CNRS: Paris, 1971. (10) Rouquerol, J. Thermochim. Acta 1989, 144, 209. (11) Toft Sorensen, O.; Rouquerol, J. Sample Controlled Thermal Analysis; Kluwer Acad.: Dordrecht, 2003. (12) Steele, W. A. J. Phys. Chem. 1978, 82, 817.

Argon and Nitrogen Adsorption in Carbons

Langmuir, Vol. 21, No. 10, 2005 4433

Figure 1. Pore size distribution of the models of CS400 (solid line) and CS1000 (dashed line). The test particle is a LennardJones model of argon. Table 1. Lennard-Jones Potential Parameters Used for Argon in This Work

Figure 2. Adsorbent-adsorbate potential energy distribution for a single argon atom in the models of CS400 (solid line) and CS1000 (dashed line) at 77 K.

We used two different atom-atom potential models for nitrogen, 2LJQ_H and 2LJQ_LB. In both models nitrogen is described with two Lennard-Jones spheres separated by a distance of 1.10 Å, plus an ideal quadrupole placed at the center of the molecule. We used the parameters derived by Cheung and Powles5 (see Table 2). We modeled each carbon atom as a Lennard-Jones sphere. Therefore, the interaction energy between a nitrogen molecule and a carbon atom is given by a sum of two Lennard-Jones interaction energies. The two models differ in the LennardJones potential parameters used for the nitrogen-carbon interaction. For the 2LJQ_H model, we used parameters fitted to the Henry’s law constant of nitrogen on “national” graphite, obtained experimentally.13 For the 2LJQ_LB model, we calculated the Lennard-Jones parameters for the nitrogen-carbon interaction using the LorentzBerthelot rules to combine the parameters for nitrogen with those for carbon obtained by Steele.14 We compare our results to those presented in ref 2 using a model for nitrogen that we term here 1LJ. In the 1LJ model, nitrogen is described as a Lennard-Jones sphere, like argon. The Lennard-Jones parameters for the nitrogen-carbon interactions used with the 1LJ model were obtained by fitting the nitrogen-carbon potential energies to those for a nonspherical model (two Lennard-Jones spheres) in the submonolayer region at 74 K.15 In Figures 1 and 2, we present the pore size distribution and the adsorbent-adsorbate potential energy distribution for a single test particle, respectively, using argon as a test particle. It is important to note that the functions used to characterize the adsorbent models depend on the test particle used to probe the structure. The results shown

in Figure 1 and Figure 2 are very similar to the analogous results for the 1LJ model of nitrogen shown in ref 2. This is not surprising, since both test particles (argon and the 1LJ model of nitrogen) are Lennard-Jones spheres of very similar size. Even though the pore size distributions of both carbon models are very similar, the single-particle adsorbent-adsorbate potential energy distributions are significantly different. The reasons and implications of this observation are presented in ref 2. We performed grand canonical Monte Carlo (GCMC) simulations of argon and models 2LJQ_H and 2LJQ_LB of nitrogen at 77 K and several bulk pressures (i.e., several chemical potentials) in the structural models of CS400 and CS1000. The acceptance criteria and other details of the Monte Carlo simulations in the grand canonical ensemble are given in.16 We applied a cutoff radius of 6σ for the argon-argon and nitrogen-nitrogen interactions, and periodic boundary conditions and minimum image convention in all directions. To speed up the energy calculation, we used an interpolation procedure17 in which the simulation box is divided into a fine cubic energy grid of 200 × 200 × 200 points and a Lennard-Jones sphere is placed at each grid point, as explained in ref 2. For each chemical potential, we calculated the amount adsorbed from the averages in the number of molecules, and used the ideal gas equation of state to calculate the density and pressure of the bulk phase in equilibrium with the adsorbed phase. The vapor pressure, P0, of argon at 77 K, 0.2 atm, was interpolated from molecular simulation results of solid-vapor equilibrium for LennardJones fluids.18 The vapor pressure of nitrogen was estimated using the Antoine equation for Lennard-Jones fluids shown in ref 19 for the 1LJ model (1.2 atm) and from Monte Carlo simulations in the isothermal-isobaric ensemble for the 2LJQ_H and 2LJQ_LB models (1 atm). We calculated the excess amount adsorbed (Γex) at each chemical potential by subtracting the number of molecules of an ideal gas in the volume of the simulation box to the average number of molecules from the corresponding simulation, and dividing by the mass of carbon. The excess amount adsorbed was then normalized with the excess amount adsorbed at P ) P0 to obtain the fractional filling, Γ/Γ0. We also calculated the isosteric heat of adsorption

(13) Cascarini de Torre, L. E.; Flores, E. S.; Llanos, J. L.; Bottani, E. J. Langmuir 1995, 11, 4742. (14) Steele, W. A. In The Interaction of Gases with Solid Surfaces; Pergamon Press: Glasgow, 1974. (15) Bottani, E. J.; Bakaev, V. A. Langmuir 1994, 10, 1550.

(16) Allen, M. P.; Tildesley, D. J. Computer Simulation of Liquids; Clarendon Press: Oxford, 1987. (17) Pellenq, R. J.-M.; Nicholson, D. Langmuir 1995, 11, 1626. (18) Agrawal, R.; Kofke, D. A. Mol. Phys. 1995, 85, 43. (19) Kofke, D. A. Adv. Chem. Phys. 1999, 105, 405.

argon-argon

σ (Å)

/k (K)

3.41

120

argon-carbon

σ (Å)

/k (K)

3.38

58

Table 2. Potential Parameters for Nitrogen-Nitrogen and Nitrogen-Carbon Interactions model

σNN (Å)

NN/k (K)

QNN/NN1/2σNN5/2

σNC (Å)

NC/k (K)

1LJ 2LJQ_LB 2LJQ_H

3.75 3.31 3.31

95.20 35.30 35.30

1.063 1.063

3.36 3.36 3.36

61.4 31.4 34.7

4434

Langmuir, Vol. 21, No. 10, 2005

Figure 3. Argon adsorption isotherms at 77 K in the models of CS400 (circles) and CS1000 (squares). The solid lines are guides to the eye.

from fluctuation theory20 and, for argon, the pair correlation function averaged over 4000 configurations using a standard algorithm that loops over all pairs of particles.16

Pikunic et al.

Figure 4. Isosteric heat of adsorption of argon at 77 K obtained from GCMC simulations in the model of CS400 as a function of fractional filling. Shown are the adsorbate-adsorbate contribution (circles), adsorbent-adsorbate contribution (triangles), and total isosteric heat of adsorption (squares). The solid lines are guides to the eye. Empty symbols indicate points at which pair correlation functions and snapshots are presented.

4. Adsorption Isotherms and Heats of Adsorption 4.1. Argon: A Test of the Structural Models of Porous Carbons. The argon adsorption isotherms at 77 K in CS400 and CS1000 are shown in Figure 3. Both adsorption isotherms are type I, according to the IUPAC classification, which is typical of microporous adsorbents. We note that we model the microporosity of real carbons; type I isotherms are thus expected. The experimental adsorption isotherms in the real materials, however, are not necessarily type I due to the presence of mesoporosity. As shown in Figure 3, the simulated fractional filling is higher in CS1000 than in CS400 for all relative pressures. The bulk pressure at which micropore filling begins in CS1000 is about 2 orders of magnitude lower than that in CS400. As pointed out in ref 2, these observations are consistent with the adsorbent-adsorbate potential energy distribution (Figure 2), but cannot be explained in terms of differences in pore size distributions (Figure 1) between the two models. Knowing the accessible pore regions from the adsorbate-carbon potential grid, one can determine the pore volume and consequently the adsorbate density. The average adsorbate density at P/P0 ) 1 in CS400 and CS1000 is 2.1 g/mL and 2.6 g/mL, respectively. These values are much larger than that for the bulk solid (1.54 g/mL) or liquid (1.39 g/mL). Similar high-density values have been obtained for nitrogen in other complex carbon models21 and for argon and xenon in silicalite zeolite.22-24 This indicates that the Gurvitch rule,3 which states that the confined fluid has the same density as the bulk liquid, is not valid for these materials with pores of the size of a few molecular diameters. We show the isosteric heat of adsorption, as well as the adsorbent-adsorbate and the adsorbate-adsorbate contributions to this quantity, as given by the simulations, in Figures 4 and 5. For both model carbons, the isosteric heat of adsorption is a continuously decreasing function (20) Nicholson, D.; Parsonage, N. G. Computer Simulation and the Statistical Mechanics of Adsorption; Academic Press: London, 1982; p 97. (21) Biggs, M. J.; Buts, A.; Williamson, D. Langmuir 2004, 20, 5786. (22) Douguet, D.; Pellenq, R. J.-M.; Boutin, A.; Fuchs, A. H.; Nicholson, D. Mol. Sim. 1996, 17, 255. (23) Pellenq, R. J.-M.; Coasne, B.; Levitz, P. E. In Adsorption and Transport Phenomena in Nanopores; Nicholson, N., Quirke, N., Eds.; Taylor and Francis: 2004, in press. (24) Pellenq, R. J.-M.; Levitz, P. E. Mol. Simul. 2001, 27, 353.

Figure 5. Same as Figure 4 for the model of CS1000.

Figure 6. Isosteric heat of adsorption of argon at 77 K. Experiment: CS400 (dashed line) and CS1000 (solid line). Simulations: CS400 (circles) and CS1000 (squares).

of fractional filling. This behavior is determined by the adsorbent-adsorbate contribution. In both models, the adsorbate-adsorbate contribution has a maximum at a fractional filling of approximately 0.8, due to the importance of adsorbate-adsorbate repulsive interactions at the higher loadings.17 The heat of adsorption decreases more rapidly in CS1000 than in CS400 at very low and high fractional filling. In Figure 6, we compare our simulation results with experimental gas adsorption microcalorimetry data in the same samples that were used to perform the X-ray diffraction experiments (see section 2). The simulation results are in excellent quantitative agreement with

Argon and Nitrogen Adsorption in Carbons

Langmuir, Vol. 21, No. 10, 2005 4435

Figure 7. Nitrogen adsorption isotherm at 77 K in CS400. Simulation results for models 1LJ (circles, from ref 2), 2LJQ_H (filled squares), and 2LJQ_LB (empty squares). The solid lines are a guide to the eye.

Figure 9. Isosteric heat of adsorption of nitrogen at 77 K in CS400. Experiment (solid line) and simulation of the models 1LJ (circles, from ref 2), 2LJQ_H (filled squares), and 2LJQ_LB (empty squares).

Figure 8. Nitrogen adsorption isotherm at 77 K in CS1000. Simulation results for models 1LJ (circles, from ref 2), 2LJQ_H (filled squares), and 2LJQ_LB (empty squares). The solid lines are a guide to the eye.

Figure 10. Isosteric heat of adsorption of nitrogen at 77 K in CS1000. Experiment (solid line) and simulation of the models 1LJ (circles, from ref 2), 2LJQ_H (filled squares), and 2LJQ_LB (empty squares).

experiments. The small discrepancies at high fractional filling for CS400 may be due to the presence of larger pores or high external surface area in this sample, which is not accounted for in our models. 4.2. Nitrogen: Comparison of Potential Models. In Figures 7 and 8, we present nitrogen adsorption isotherms at 77 K obtained from GCMC simulations in CS400 and CS1000, respectively, for the 2LJQ_H and 2LJQ_LB models (see section 3). All the isotherms are type I, as in the case of argon. We compare the results with those presented in ref 2 for the 1LJ model. At low pressure, the excess amount adsorbed for model 2LJQ_H is similar to that obtained for model 1LJ. This is due to the fact that the potential parameters for the nitrogen-carbon interaction in both models were derived from similar experimental results (fitting to Henry’s law constants of nitrogen in nonporous carbons). The excess amount adsorbed for the model 2LJQ_LB, however, is significantly lower than that for the other two models of nitrogen. This illustrates the deviations caused by deriving adsorbate-adsorbent potential parameters through Lorentz-Berthelot rules. At low pressures, a simple model of the adsorbate (a Lennard-Jones sphere in this case) with parameters derived from experimental results, correctly reproduces the adsorption isotherm obtained with the more detailed model. At high pressure, in contrast, the adsorption isotherm obtained with 2LJQ_H deviates from that of 1LJ and tends to behave more like the isotherm obtained with the 2LJQ_LB model. At high pressure, the adsorbateadsorbate interactions become more important than at low pressure and the repulsive contributions (both ad-

sorbate-adsorbate and adsorbate-adsorbent) to the potential energy are significantly larger. Thus, an adequate atomistic description (molecular geometry) of the adsorbate is more important than accurately obtained adsorbate-adsorbent potential energy parameters in this region. In Figures 9 and 10, we present the isosteric heat of adsorption obtained from simulations of models 2LJQ_H and 2LJQ_LB. We compare our results to those for model 1LJ2 and experiment. As expected, for both CS400 and CS1000, the heats of adsorption for model 2LJQ_H give the best agreement with experimental results. The heats of adsorption calculated with this model are similar to those obtained with 1LJ at low fractional filling, and to those using 2LJQ_LB at high fractional filling. This is consistent with the trends observed in the adsorption isotherms. The simulation results using model 2LJQ_H are in quantitative agreement with experiment in the range 0.25 < Γ/Γo < 0.85. We attribute the disagreement at low filling to the presence of polar groups (containing oxygen) in the real porous material that strongly interact with the nitrogen quadrupole (see section 1). These polar groups are omitted in the simulations and are absent in the carbon material used to fit adsorbate-adsorbent potential parameters. Also, the anisotropy of the nitrogen-carbon interaction is not taken into account in the potential used to describe this interaction (see section 1). In other words, the interaction of nitrogen with carbon atoms with carbon coordination number of 2 and hydrogen coordination number of 1 must be different to that with carbon atoms bonded to three other carbon atoms. Similarly, there is a

4436

Langmuir, Vol. 21, No. 10, 2005

Figure 11. Carbon-argon pair correlation function in CS400 at Γ/Γ0 ) 0.039 (circles); 0.32 (squares); 0.67 (crosses); 1.0 (triangles).

Figure 12. Carbon-argon pair correlation function in CS1000 at Γ/Γ0 ) 0.048 (circles); 0.34 (squares); 0.68 (crosses); 1.0 (triangles).

difference in interaction of nitrogen with carbon atoms in a planar section of graphene and those in curved or defective sections.26 As mentioned in the case of argon, the disagreement at high filling is presumably caused by the presence of larger pores (i.e., larger micropores and mesopores). The structural models for porous carbons presented in ref 2 are intended to describe only the microporosity of the real material. Therefore, agreement in this region should only be expected if the material has no large micropores and mesopores. Moreover, we have estimated the statistical error of the heat of adsorption at fractional filling higher than 0.8 to be between 9 and 15%. 5. Pair Correlation Functions and Adsorption Mechanism In Figures 11 and 12, we present carbon-argon pair correlation functions in the models of CS400 and CS1000 at four different values of fractional filling, indicated with empty symbols in Figures 4 and 5. The first peak in the carbon-argon correlation functions at low fractional filling occurs at an interatomic distance of approximately 3.8 Å, which corresponds to the minimum of the carbon-argon Lennard-Jones potential (21/6σ). This peak, however, shifts to lower distances (ca. 3.6 Å) at higher fractional fillings, showing that under pressure, the adsorbate molecules are pressed against the pore wall and probe the repulsive part of the adsorbate-adsorbent potential energy. The height of the first peak decreases with fractional filling.

Pikunic et al.

Figure 13. Argon-argon pair correlation function in CS400 at Γ/Γ0 ) 0.32 (squares); 0.67 (crosses); 1.0 (triangles).

Figure 14. Argon-argon pair correlation function in CS1000 at Γ/Γ0 ) 0.34 (squares); 0.68 (crosses); 1.0 (triangles).

The decrease is larger at lower fractional filling, especially for CS1000, and almost negligible at higher filling. A decrease in height and shift to lower distances is also observed for the second and third peaks. The area under the first peaks is directly related to the number of carbon atoms in the immediacy of an argon atom. Thus, in average, argon atoms interact with less carbon atoms at higher fractional filling, which is also a reason for the decrease on the heat of adsorption with fractional filling. Comparing the carbon-argon pair correlation functions at the lowest fractional filling, it is obvious that there are longer range correlations in CS1000 than in CS400. This is a consequence of the differences in structural order of the adsorbent (see carbon-carbon pair correlation functions in ref 2). We also present argon-argon pair correlation functions at different values of fractional filling in CS400 (Figure 13) and CS1000 (Figure 14), as well as snapshots of the simulations at the same values of fractional filling (see Figures 15-22). In both models, the first peak of the argon-argon pair correlation function becomes more pronounced as the fractional filling increases, indicating an increase in the number of argon atoms in the first coordination shell. At the two lower values of fractional filling, the first peak appears at a distance of 3.7 Å, which corresponds to the minimum in the argon-argon LennardJones potential. Thus, the increase in the number of atoms in the first coordination shell produces the observed increase in the adsorbate-adsorbate contribution to the isosteric heat of adsorption (Figure 4 and Figure 5). At high fractional filling, however, the first peak shifts to

Argon and Nitrogen Adsorption in Carbons

Langmuir, Vol. 21, No. 10, 2005 4437

Figure 15. (a) Snapshots of GCMC simulations of argon at 77 K in the resulting model of CS400 at Γ/Γ0 ) 0.039. The rods and the spheres represent C-C bonds and argon atoms, respectively. (b) Projection of a 1.25 nm section.

Figure 16. Same as Figure 15 for Γ/Γ0 ) 0.32.

Figure 17. Same as Figure 15 for Γ/Γ0 ) 0.67.

smaller distances. The argon atoms thus probe the repulsive part of the argon-argon interaction potential, causing the observed slight decrease in the contribution to the isosteric heat of adsorption.

At low fractional filling, there is a peak at approximately 8.7 Å in the argon-argon pair correlation functions for CS1000 and 8.3 Å in that for CS400, which correspond to correlations between argon atoms adsorbed at different

4438

Langmuir, Vol. 21, No. 10, 2005

Pikunic et al.

Figure 18. Same as Figure 15 for Γ/Γ0 ) 1.0.

Figure 19. Same as Figure 15 but for the model of CS1000 for Γ/Γ0 ) 0.048.

Figure 20. Same as Figure 19 for Γ/Γ0 ) 0.34.

sites, either in the same pore or in neighboring pores. This distance corresponds to about 2.5 atomic diameters. At higher fractional fillings, these peaks are replaced by ones at lower distances (6.9 and 7.2 Å), corresponding to 2 atomic diameters, due to the complete

filling of the pores and consequent closer packing of the argon atoms. This distance is equal to that at which the second peak of the pair correlation function of the bulk liquid is observed. At intermediate fractional filling, both peaks are present.

Argon and Nitrogen Adsorption in Carbons

Langmuir, Vol. 21, No. 10, 2005 4439

Figure 21. Same as Figure 19 for Γ/Γ0 ) 0.68.

Figure 22. Same as Figure 19 for Γ/Γ0 ) 1.0.

No defined peaks are observed at higher distances. In contrast, a third peak at 3 atomic diameters is observed in the pair correlation function of the bulk liquid. The lack of longer-range correlations in the confined fluid is a consequence of the small pore dimensions. Some structure is observed at low fractional filling, but the peaks do not correspond to those observed in the bulk phase. Instead, at low filling we are probing the pore-pore correlations. More structural features are thus observed in the low filling argon-argon pair correlation function for CS1000 than in that for CS400, due to the higher structural order of the former. It is interesting to note that the more attractive points are not localized in specific regions of the material; instead, they are distributed in different pores (see Figures 15, 16, 19, and 20). 6. Discussion and Conclusions In a previous work,2 Pikunic et al. developed realistic models of two saccharose-based carbons from X-ray diffraction data, density, and composition measurements, using a constrained reverse Monte Carlo method. It was found that simulated TEM images were in qualitative agreement with experimental images; the simulated images capture the most important features observed experimentally and also the change of these features with heat treatment temperature.

In this work, we have tested these models by predicting the isosteric heats of adsorption of argon at 77 K and comparing our predictions to experimental gas adsorption microcalorimetry data. The simulated isosteric heats of adsorption in our models of CS400 and CS1000 are in excellent quantitative agreement with the experiments (Figure 6). This implies that the resulting structural models reproduce the potential energy field “seen” by argon in the real materials and, thus, seem to successfully describe the texture of CS400 and CS1000. The isosteric heats of adsorption decrease as a function of coverage. The shape of the experimental curves, which is strongly related to the energetic heterogeneity of the material, is almost exactly reproduced by our models. The small discrepancies at higher fractional filling might be due to the presence of larger pores or external surfaces that are not accounted for in our model. We found that the decrease in the isosteric heat of adsorption is dominated by the adsorbent-adsorbate contribution, as observed in other microporous adsorbents such as zeolites.17 We also presented adsorption isotherms and heats of adsorption at 77 K for three models of nitrogen in CS400 and CS1000. Isosteric heats of adsorption obtained with the 2LJQ_H model are in very good agreement with experiment in the range 0.25 < Γ/Γ0 < 0.85. We attribute discrepancies at low pressures to the presence of polar

4440

Langmuir, Vol. 21, No. 10, 2005

groups or attractive points that strongly interact with the nitrogen quadrupole. This is consistent with the fact that this discrepancy is not observed in the case of argon, which is spherical and nonpolar. We also note that, if a realistic description of the solid structure is available, it is important to have a relatively accurate description of the adsorbate in order to obtain good agreement with experiment. Adsorbate-adsorbent potential parameters derived from experimental data are important to obtain the correct behavior at low densities, while a good geometrical description of the adsorbate is necessary at higher densities. We calculated the carbon-argon and argon-argon pair correlation functions. These results, aided by simulation snapshots, can help determine the mechanism of the adsorption process. At very low chemical potentials, the first argon atoms are preferentially adsorbed at the points of lowest carbon-argon potential energy (i.e., the most attractive points), which are characterized by a high local density of carbon atoms, as evidenced by the carbonargon pair correlation functions (see Figures 11 and 12). These low-energy points are not grouped, but instead they are dispersed throughout the porous network. Carbon atoms are more closely packed in CS1000,2 producing points of higher local carbon density, and thus lower energy points. It is this difference in the structure of the two carbons, and not a difference in pore size distributions, what causes the higher heats of adsorption and higher fractional fillings at lower pressures observed in CS1000. As we have previously pointed out,2,25 interpreting the adsorption behavior observed in CS400 and CS1000 using the slit pore model, would lead to conclude that the observed differences are due to dissimilarities between the pore size distributions of CS400 and CS1000, which is not necessarily the case (see Figure 1). At slightly higher chemical potentials, argon atoms are adsorbed at higher energy points, which are also uniformly distributed throughout the porous network. The average number of carbon atoms in the immediacy of an argon atom decreases, causing the observed decrease in the isosteric heats of adsorption. The larger energy difference between lower carbon-argon energy points in CS1000 causes the decrease in this material to be steeper than that in CS400. A peak in the argon-argon pair correlation function at approximately 2.5 atomic diameters appears due to correlations of argon atoms adsorbed in neighbor sites. (25) Coasne, B.; Pikunic, J. P.; Pellenq, R. J.-M.; Gubbins, K. E. Presented at the Material Society Fall Meeting, Boston, MA, U.S.A., 2003; P8.5. (26) Klauda, J. B.; Jiang, J.; Sandler, S. I. J. Phys. Chem. B 2004, 108, 9842.

Pikunic et al.

Further increasing the chemical potential causes a gradual progression toward micropore filling. When the pores are completely filled, a peak in the argon-argon pair correlation function appears at 2 atomic diameters, which corresponds to the distance at which the second peak of the bulk liquid pair correlation is observed. No defined peaks, however, are observed at higher distances. This reflects the effect of the small pore dimensions on the structure of the confined fluid, which cannot be considered as being bulk liquid. Our results show that the atomic-level structure of disordered porous carbons is essential in explaining the adsorbent-adsorbate potential energy distribution in these materials, which ultimately determines the adsorption mechanism, and thus macroscopic observables (e.g., adsorption isotherm and isosteric heat of adsorption). In improving adsorption-based characterization methods for porous carbons, we suggest the use of isosteric heats of adsorption instead of (or to complement) adsorption isotherms. The shape of the isosteric heat of adsorption as a function of fractional filling is much more sensitive than that of the adsorption isotherm to the heterogeneity caused by the local structure of the adsorbent. Also, it seems that future developments in adsorption-based characterization methods for porous carbons should focus on obtaining adsorbent-adsorbate potential energy distributions, rather than pore size distributions. If the concept of pore-size distribution is preferred, however, future improvements should include heterogeneity in the density of the pore walls. Finally, many realistic models have been developed and are currently being developed. These models take into account the finite length of carbonaceous segments. It is thus necessary to develop new carbon-adsorbate potentials that take into account the effect of the location of the carbon atoms involved. Acknowledgment. This work was funded by the National Science Foundation, through the Grant CTS0211792; international cooperation was supported by the NSF/CNRS Award No. INT-0089696. Supercomputing access was supported in part by the National Resource Allocation Committee (Grant MCA93S011) through computing resources provided by the National Partnership for Advanced Computational Infrastructure at the San Diego Supercomputer Center and Texas Advanced Computing Center. LA047165W