Hierarchical modeling N2 adsorption on the surface of and within a C60 crystal: from quantum mechanics to molecular simulation

The adsorption of N(2) on the surface of, and within, a C(60) face-centered cubic crystal has been studied using a hierarchical approach. First, an ab initio potential between N(2) and C(60) is obtained from a recently developed quantum mechanical hybrid method, and then the adsorption behavior is predicted using Monte Carlo simulation. On the crystal surface, N(2) adsorption isotherm at 77.3 K is of type II. The adsorption simulated with the ab initio potential is slightly greater than that with the empirical Steele potential derived from experimental N(2) adsorption on planar graphite, and both are in fairly good agreement with measured results. With increasing pressure, N(2) molecules are found to sequentially occupy three favorable sites: the octahedral sites, the tetrahedral sites, and the top of C(60) molecules. Finally multiple layers form and wetting occurs as the bulk N(2) saturation pressure is reached. The isosteric heat of adsorption exhibits two maxima and finally approaches the enthalpy of vaporization of bulk N(2). Within the crystal, the N(2) adsorption isotherm at 77.3 K is of type I, and the use of ab initio potential leads to significantly greater adsorption than the Steele potential. N(2) molecules are observed to intercalate only the octahedral sites, and the isosteric heat of adsorption is nearly a constant. As in our previous work of N(2) and O(2) adsorption in the C(168) schwarzite (Jiang et al. J. Phys. Chem. B 2004, 108, 9852), this work demonstrates the importance of an accurate adsorbate-adsorbent interaction potential in the determination of gas adsorption behavior.