Close packed transitions in slit-shaped pores: Density functional theory study of methane adsorption capacity in carbon

An important feature of improving lattice gas models and classical isotherms is the incorporation of a pore size dependent capacity, which has hitherto been overlooked. In this paper, we develop a model for predicting the temperature dependent variation in capacity with pore size. The model is based on the analysis of a lattice gas model using a density functional theory approach at the close packed limit. Fluid-fluid and solid-fluid interactions are modeled by the Lennard-Jones 12-6 potential and Steele's 10-4-3, potential respectively. The capacity of methane in a slit-shaped carbon pore is calculated from the characteristic parameters of the unit cell, which are extracted by minimizing the grand potential of the unit cell. The capacities predicted by the proposed model are in good agreement with those obtained from grand canonical Monte Carlo simulation, for pores that can accommodate up to three adsorbed layers. Single particle and pair distributions exhibit characteristic features that correspond to the sequence of buckling and rhombic transitions that occur as the slit pore width is increased. The model provides a useful tool to model continuous variation in the microstructure of an adsorbed phase, namely buckling and rhombic transitions, with increasing pore width. (C) 2002 American Institute of Physics.

Adsorption of argon from sub- to supercritical conditions on graphitized thermal carbon black and in graphitic slit pores: A grand canonical Monte Carlo simulation study

In this paper we consider the adsorption of argon on the surface of graphitized thermal carbon black and in slit pores at temperatures ranging from subcritical to supercritical conditions by the method of grand canonical Monte Carlo simulation. Attention is paid to the variation of the adsorbed density when the temperature crosses the critical point. The behavior of the adsorbed density versus pressure (bulk density) shows interesting behavior at temperatures in the vicinity of and those above the critical point and also at extremely high pressures. Isotherms at temperatures greater than the critical temperature exhibit a clear maximum, and near the critical temperature this maximum is a very sharp spike. Under the supercritical conditions and very high pressure the excess of adsorbed density decreases towards zero value for a graphite surface, while for slit pores negative excess density is possible at extremely high pressures. For imperfect pores (defined as pores that cannot accommodate an integral number of parallel layers under moderate conditions) the pressure at which the excess pore density becomes negative is less than that for perfect pores, and this is due to the packing effect in those imperfect pores. However, at extremely high pressure molecules can be packed in parallel layers once chemical potential is great enough to overcome the repulsions among adsorbed molecules. (c) 2005 American Institute of Physics.

Simulation study of ammonia adsorption on graphitized carbon black

GCMC simulations are applied to the adsorption of sub-critical ammonia on graphitized carbon black at 240 K. The carbon black was modelled both with and without carbonyl functional groups. Large differences are seen between the amount adsorbed for different carbonyl configurations at low pressure (P < 10kPa). Once a single layer is formed on the carbon black, the adsorption behaviour is similar between the model surfaces with and without functional groups. Simulation isotherms are qualitatively similar to the few experimental isotherms available in the literature for ammonia on highly graphitized carbon black. The mode of adsorption up to monolayer coverage is exhaustively shown to be two-dimensional clustering using various techniques. A comparison between experiment and simulation isosteric heats shows that a surface without functional groups cannot reproduce the experimental isosteric heats of adsorption, even comparing with the experimental results of carbon black heat treated at 3373 K. The addition of carbonyls produces isosteric heats with similar features to those in the literature if the separation between the carbonyls is small.

Effect of pore constriction on adsorption behaviour in nanoporous carbon slit pore: a computer simulation

A Grand Canonical Monte Carlo simulation (GCMC) method is used to study the effects of pore constriction on the adsorption of argon at 87.3 K in carbon slit pores of infinite and finite lengths. It is shown that the pore constriction affects the pattern of adsorption isotherm. First, the isotherm of the composite pore is greater than that of the uniform pore having the same width as the larger cavity of the composite pore. Secondly, the hysteresis loop of the composite pore is smaller than and falls between those of uniform pores. Two types of hysteresis loops have been observed, irrespective of the absence or presence of constriction and their presence depend on pore width. One hysteresis loop is associated with the compression of adsorbed particles and this phenomenon occurs after pore has been filled with particles. The second hysteresis loop is the classical condensation-evaporation loop. The hysteresis loop of a composite pore depends on the sizes of the larger cavity and the constriction. Generally, it is found that the pore blocking effect is not manifested in composite slit pores, and this result does not support the traditional irkbottle pore hypothesis.

The adsorption of water in finite carbon pores

Grand canonical Monte Carlo simulations were applied to the adsorption of SPCE model water in finite graphitic pores with different configurations of carbonyl functional groups on only one surface and several pore sizes. It was found that almost all finite pores studied exhibit capillary condensation behaviour preceded by adsorption around the functional groups. Desorption showed the reverse transitions from a filled to a near empty pore resulting in a clear hysteresis loop in all pores except for some of the configurations of the 1.0nm pore. Carbonyl configurations had a strong effect on the filling pressure of all pores except, in some cases, in 1.0nm pores. A decrease in carbonyl neighbour density would result in a higher filling pressure. The emptying pressure was negligibly affected by the configuration of functional groups. Both the filling and emptying pressures increased with increasing pore size but the effect on the emptying pressure was much less. At pressures lower than the pore filling pressure, the adsorption of water was shown to have an extremely strong dependence on the neighbour density with adsorption changing from Type IV to Type III to linear as the neighbour density decreased. The isosteric heat was also calculated for these configurations to reveal its strong dependence on the neighbour density. These results were compared with literature experimental results for water and carbon black and found to qualitatively agree.

A wavelet-based adaptive technique for adsorption problems involving steep gradients

A general, fast wavelet-based adaptive collocation method is formulated for heat and mass transfer problems involving a steep moving profile of the dependent variable. The technique of grid adaptation is based on sparse point representation (SPR). The method is applied and tested for the case of a gas–solid non-catalytic reaction in a porous solid at high Thiele modulus. Accurate and convergent steep profiles are obtained for Thiele modulus as large as 100 for the case of slab and found to match the analytical solution.

Computationally efficient solution techniques for adsorption problems involving steep gradients in bidisperse particles

A piecewise uniform fitted mesh method turns out to be sufficient for the solution of a surprisingly wide variety of singularly perturbed problems involving steep gradients. The technique is applied to a model of adsorption in bidisperse solids for which two fitted mesh techniques, a fitted-mesh finite difference method (FMFDM) and fitted mesh collocation method (FMCM) are presented. A combination (FMCMD) of FMCM and the DASSL integration package is found to be most effective in solving the problems. Numerical solutions (FMFDM and FMCMD) were found to match the analytical solution when the adsorption isotherm is linear, even under conditions involving steep gradients for which global collocation fails. In particular, FMCMD is highly efficient for macropore diffusion control or micropore diffusion control. These techniques are simple and there is no limit on the range of the parameters. The techniques can be applied to a variety of adsorption and desorption problems in bidisperse solids with non-linear isotherm and for arbitrary particle geometry.

Adsorption in mesopores: A molecular-continuum model with application to MCM-41

The classical model of surface layering followed by capillary condensation during adsorption in mesopores, is modified here by consideration of the adsorbate solid interaction potential. The new theory accurately predicts the capillary coexistence curve as well as pore criticality, matching that predicted by density functional theory. The model also satisfactorily predicts the isotherm for nitrogen adsorption at 77.4 K on MCM-41 material of various pore sizes, synthesized and characterized in our laboratory, including the multilayer region, using only data on the variation of condensation pressures with pore diameter. The results indicate a minimum mesopore diameter for the surface layering model to hold as 14.1 Å, below which size micropore filling must occur, and a minimum pore diameter for mechanical stability of the hemispherical meniscus during desorption as 34.2 Å. For pores in-between these two sizes reversible condensation is predicted to occur, in accord with the experimental data for nitrogen adsorption on MCM-41 at 77.4 K.

Atomic Latin squares of order eleven

A Latin square is pan-Hamiltonian if the permutation which defines row i relative to row j consists of a single cycle for every i j. A Latin square is atomic if all of its conjugates are pan-Hamiltonian. We give a complete enumeration of atomic squares for order 11, the smallest order for which there are examples distinct from the cyclic group. We find that there are seven main classes, including the three that were previously known. A perfect 1-factorization of a graph is a decomposition of that graph into matchings such that the union of any two matchings is a Hamiltonian cycle. Each pan-Hamiltonian Latin square of order n describes a perfect 1-factorization of Kn,n, and vice versa. Perfect 1-factorizations of Kn,n can be constructed from a perfect 1-factorization of Kn+1. Six of the seven main classes of atomic squares of order 11 can be obtained in this way. For each atomic square of order 11, we find the largest set of Mutually Orthogonal Latin Squares (MOLS) involving that square. We discuss algorithms for counting orthogonal mates, and discover the number of orthogonal mates possessed by the cyclic squares of orders up to 11 and by Parker's famous turn-square. We find that the number of atomic orthogonal mates possessed by a Latin square is not a main class invariant. We also define a new sort of Latin square, called a pairing square, which is mapped to its transpose by an involution acting on the symbols. We show that pairing squares are often orthogonal mates for symmetric Latin squares. Finally, we discover connections between our atomic squares and Franklin's diagonally cyclic self-orthogonal squares, and we correct a theorem of Longyear which uses tactical representations to identify self-orthogonal Latin squares in the same main class as a given Latin square.

Superchemistry: dynamics of coupled atomic and molecular Bose condensates

We analyze the dynamics of a dilute, trapped Bose-condensed atomic gas coupled to a diatomic molecular Bose gas by coherent Raman transitions. This system is shown to result in a new type of “superchemistry,” in which giant collective oscillations between the atomic and the molecular gas can occur. The phenomenon is caused by stimulated emission of bosonic atoms or molecules into their condensate phases.

Quantum dynamics of three coupled atomic Bose-Einstein condensates

The simplest model of three coupled Bose-Einstein condensates is investigated using a group theoretical method. The stationary solutions are determined using the SU(3) group under the mean-field approximation. This semiclassical analysis, using system symmetries, shows a transition in the dynamics of the system from self trapping to delocalization at a critical value for the coupling between the condensates. The global dynamics are investigated by examination of the stable points, and our analysis shows that the structure of the stable points depends on the ratio of the condensate coupling to the particle-particle interaction, and undergoes bifurcations as this ratio is varied. This semiclassical model is compared to a full quantum treatment, which also displays a dynamical transition. The quantum case has collapse and revival sequences superimposed on the semiclassical dynamics, reflecting the underlying discreteness of the spectrum. Nonzero circular current states are also demonstrated as one of the higher-dimensional effects displayed in this system.

Stimulated Raman adiabatic passage from an atomic to a molecular Bose-Einstein condensate

The process of stimulated Raman adiabatic passage (STIRAP) provides a possible route for the generation of a coherent molecular Bose-Einstein condensate (BEC) from an atomic BEC. We analyze this process in a three-dimensional mean-field theory, including atom-atom interactions and nonresonant intermediate levels. We find that the process is feasible, but at larger Rabi frequencies than anticipated from a crude single-mode lossless analysis, due to two-photon dephasing caused by the atomic interactions. We then identify optimal strategies in STIRAP allowing one to maintain high conversion efficiencies with smaller Rabi frequencies and under experimentally less demanding conditions.

Adsorption characteristics of carbonized bark for phenol and pentachlorophenol

The potential of using carbonized slash pine bark as a substitute for activated carbon was examined in this study. The bark was carbonized by slow heating in nitrogen for 6.5 h to 672 degrees C. The BET-N-2 surface area, average micropore and mesopore diameter, and micropore volume were 332 m(2) g(-1) 21.7 Angstrom, and 0.125 cm(3) g(-1), respectively. The adsorption capacities for phenol and pentachlorophenol (PCP) at pH 2 and pH 8 were evaluated. The Langmuir equation provided a slightly better fit than the Freundlich equation to two sets of phenol data. The calculated Freundlich constants, K = 0.41 - 0.58 mmol/g/(mmol dm(-3))(1/n) and 1/n = 0.30 - 0.41, were lower and higher, respectively, than literature values for activated carbons. The adsorption capacity of the carbonized bark was much lower for PCP than for phenol. The protonated and anionic PCP isotherms were Type II or III, respectively, in the Brunauer classification. The BET equation provided the best fit to protonated PCP isotherm data. The anionic PCP data were fitted to both the BET model and an equation used in the literature to represent phosphate adsorption on activated carbons. Nonlinear regression of the data for both phenol and PCP adsorption with the Freundlich, Langmuir and BET equations generally gave more accurate parameters, compared with the use of linearized equations to obtain the parameters. (C) 1998 SCI.

Modification of MCM-41 by surface silylation with trimethylchlorosilane and adsorption study

Siliceous MCM-41 samples were modified by silylation using trimethylchlorosilane (TMCS). The surface coverage of functional groups was studied systematically in this work. The role of surface silanol groups during modification was evaluated using techniques of FTIR and Si-29 CP/MAS NMR. Adsorption of water and benzene on samples of various hydrophobicities was measured and compared. It was found that the maximum degree of surface attachments of trimethylsilyl (TMS) groups was about 85%, corresponding to the density of TMS groups of 1.9 per nm(2). The degree of silylation is found to linearly increase with increasing pre-outgassing temperature prior to silylation. A few types of silanol groups exist on MCM-41 surfaces, among which both free and geminal ones are responsible for active silylation. Results of water adsorption show that aluminosilicate MCM-41 materials are more or less hydrophilic, giving a type IV isotherm, similar to that of nitrogen adsorption, whereas siliceous MCM-41 are hydrophobic, exhibiting a type V adsorption isotherm. The fully silylated Si-MCM-41 samples are more hydrophobic giving a type III adsorption isotherm. Benzene adsorption on all MCM-41 samples shows type IV isotherms regardless of the surface chemistry. Capillary condensation occurs at a higher relative pressure for the silylated MCM-41 than that for the unsilylated sample, though the pore diameter was found reduced markedly by silylation. This is thought attributed to the diffusion constriction posed by the attached TMS groups. The results show that the surface chemistry plays an important role in water adsorption, whereas benzene adsorption is predominantly determined by the pore geometry of MCM-41.

Experimental and theoretical investigations of adsorption hysteresis and criticality in MCM-41: Studies with O-2, Ar, and CO2

MCM-41 materials of six different pore diameters were prepared and characterized using X-ray diffraction, transmission electron microscopy, helium pycnometry, small-angle neutron scattering, and gas adsorption (argon at 77.4 and 87.4 K, nitrogen and oxygen at 77.4 K, and carbon dioxide at 194.6 K). A recent molecular continuum model of the authors, previously used for adsorption of nitrogen at 77.4 K, was applied here for adsorption of argon, oxygen, and carbon dioxide. While model predictions of single-pore adsorption isotherms for argon and oxygen are in satisfactory agreement with experimental data, significant deviation was found for carbon dioxide, most likely due to its high quadrupole moment. Predictions of critical pore diameter, below which reversible condensation occurs: were possible by the model and found to be consistent with experimental estimates, for the adsorption of the various gases. On the other hand, existing models such as the Barrett-Joyner-Halenda (BJH), Saito-Foley, and Dubinin-Astakhov models were found to be inadequate, either predicting an incorrect pore diameter or not correlating the isotherms adequately. The wall structure of MCM-41 appears to be close to that of amorphous silica, as inferred from our skeletal density measurements.