Computational prediction of adsorption of small molecules in porous materials has great impact on the basic and applied research in chemical engineering and material sciences. In this work, we report an approach based on grand canonical ensemble Monte Carlo (GCMC) simulations and ab initio force fields. We calculated the adsorption curves of ammonia in ZSM-5 zeolite and hydrogen in MOF-5 (a metal-organic-framework material). The predictions agree well with experimental data. Because the predictions are based on the first principle force fields, this approach can be used for the adsorption prediction of new molecules or materials without experimental data as guidance.Porous materials, adsorption, ab initio force field, Grand Canonical ensemble Monte Carlo Adsorption of gaseous or liquid small molecules in porous materials has great potentials in wide-spread applications. For example, zeolites, a group of crystalline silicon/alumina oxides with nano-scale porous structures, have been widely used in the fine chemical and oil industries [1][2][3] for their capabilities of ion-exchange, catalysis, and adsorption; Metal organic frameworks (MOFs), a new type of materials with extremely large surfaces and controllable pore sizes, have been explored for the hydrogen storage purpose with great interest [4][5][6][7] . The adsorption of small molecules in porous materials is extremely important for the application of these materials. In zeolites, the pore structures, capabilities of adsorption of reactants and acidities play significant roles in determining their catalytic properties. For MOFs, the ability of adsorbing hydrogen molecules on their surfaces directly reflects their capacity of hydrogen storage. Therefore, to study the adsorption of small molecules in porous materials is of great interest.Significant progress has been made in studying adsorption of small molecules in porous materials using experimental methods. Researchers have gained substantial knowledge of the chemical and physical properties of porous materials using techniques such as the BET isotherms measurement and Temperature-Programmed Desorption (TPD). However, due to inherent limitations of experimental measurements, it is difficult to reveal the mechanism of small molecules on the surfaces at the microscopic level, where the theoretical or computational approaches can be very valuable. In recent years, with the fast development of quantum chemistry and molecular modeling methods, computational chemistry methods have been widely used in studying the adsorption of small molecules in porous materials. The advantage of the computational approach is that it can provide more precise description of the behavior of guest molecules in the porous materials, explain the experimental observed phenomena and perhaps predict the adsorption under different conditions. Among various computational approaches, the grand canonical ensemble Monte Carlo (GCMC) method is capable for predicting the phase equilibrium of gas molecules and the adsorbent directly.Underlying th...