Macromolecule simulation and CH4 adsorption mechanism of coal vitrinite

Abstract The microscopic mechanism of interactions between CH4 and coal macromolecules is of significant practical and theoretical importance in CBM development and methane storage. Under periodic boundary conditions, the optimal energy configuration of coal vitrinite, which has a higher torsion degree and tighter arrangement, can be determined by the calculation of molecular mechanics (MM) and molecular dynamics (MD), and annealing kinetics simulation based on ultimate analysis, 13C NMR, FT IR and HRTEM. Macromolecular stabilization is primarily due to the van der Waals energy and covalent bond energy, mainly consisting of bond torsion energy and bond angle energy. Using the optimal configuration as the adsorbent, GCMC simulation of vitrinite adsorption of CH4 is conducted. A saturated state is reached after absorbing 17 CH4s per coal vitrinite molecule. CH4 is preferentially adsorbed on the edge, and inclined to gathering around the branched chains of the inner vitrinite sites. Finally, the adsorption parameters are calculated through first principle DFT. The adsorbability order is as follows: aromatic structure> heteroatom rings > oxygen functional groups. The adsorption energy order is as follows: Top  C O] and [C O C]. CH4 distributed in the distance of 0.99–16 A to functional groups in the type of monolayer adsorption and the average distance order manifest as [ C O] (1.64 A)  O C] (1.89 A)  COOH] (3.78 A)  C O] and [C O-C] whereas is rather dispersed about [-COOH] and [CH3]. Simulation and experiment data are both in strong agreement with the Langmuir and D-A isothermal adsorption model and the D-A model fit better than Langmuir model. Preferential adsorption sites and orientations in vitrinite are identical to those of graphite/graphene. However, the energy of the most preferential location is much lower than that of graphite/graphene. CH4 is more easily absorbed on the surface of vitrinite. Adsorbability varies considerably at different adsorption locations and sites on the surface of vitrinite. Crystal parameter of vitrinite is a = b = c = 15.8 A and majority of its micropores are blow 15.8 A, indicating that the vitrinite have the optimum adsorption aperture. It can explain its higher observed adsorption capacities for CH4 compared with graphite/graphene.