Molecular Dynamics Simulation of Di ﬀ usion Behavior of CH 4 , CO 2 , and N 2 in Mid-Rank Coal Vitrinite

: The di ﬀ usion characteristics of CH 4 , CO 2 , and N 2 in coal are important for the study of CO 2 -enhanced coalbed methane (CO 2 -ECBM) recovery, which has become the most potential method for carbon sequestration and natural gas recovery. However, quantitative research on the di ﬀ usion characteristics of CH 4 and the invasive gases (CO 2 and N 2 ) in coal, especially those in micropores, still faces enormous challenges. In this paper, the self-, Maxwell’s, and transport di ﬀ usions of CO 2 , CH 4 , and N 2 in mid-rank coal vitrinite (MRCV) macromolecules were simulated based on the molecular dynamics method. The e ﬀ ects of the gas concentration, temperature, and pressure on the di ﬀ usion coe ﬃ cients were examined via the comparison of various ranks. The results indicated that the di ﬀ usion coe ﬃ cients have the order of D (N 2 ) > D (CO 2 ) > D (CH 4 ) in their saturated adsorption states. However, when MRCV adsorbed the same amounts of CH 4 , CO 2 , and N 2 , the self- and transport di ﬀ usion coe ﬃ cients followed the order of D S (N 2 ) > D S (CO 2 ) > D S (CH 4 ) and D t (CO 2 ) > D t (N 2 ) > D t (CH 4 ), respectively. Independent of the gas species, all these di ﬀ usion coe ﬃ cients decreased with increasing gas concentration and increased with increasing temperature. In the saturated adsorption state, the di ﬀ usion activation energies of CH 4 , CO 2 , and N 2 were ordered as CH 4 (27.388 kJ / mol) > CO 2 (11.832 kJ / mol) > N 2 (10.396 kJ / mol), indicating that the di ﬀ usion processes of CO 2 and N 2 occur more easily than CH 4 . The increase of temperature was more conducive to the swelling equilibrium of coal. For the pressure dependence, the di ﬀ usion coe ﬃ cients ﬁrst increased until the peak pressure (3 MPa) and then decreased with increasing pressure. In contrast, the di ﬀ usion activation energy ﬁrst decreased and then increased with increasing pressure, in which the peak pressure was also 3 MPa. The swelling rate changed more obviously in high-pressure conditions.


Introduction
As a clean energy source, which occurs in the adsorption state, coalbed methane (CBM) is an important, unconventional natural gas resource in the United States, Australia, China, and other countries [1][2][3], primarily consisting of CH 4 and a small amount of heavy hydrocarbons, CO 2 , N 2 , and H 2 O [2]. Compared with CH 4 , coal seams can adsorb CO 2 more easily. Injecting CO 2 and N 2 into the coal seams can promote CH 4 to escape from the coal seams [4][5][6][7][8]. Therefore, the unmineable coal seams have the potential to store CO 2 and enhanced coalbed methane recovery. During the production process of CBM, CH 4 in the adsorption state is first desorbed from the optimal sites in the micropores, then diffused into the mesopores, macropores, and cleats of the coal seam through the pore network, and finally enters the production well through a continuous seepage flow [9]. Therefore, Table 1. Proximate analysis and element analysis of the TL-1 coal sample [28].

Macromolecular Structure Construction
Before construction of the 3D macromolecular structure, the size of the aromatics was determined by 13 C NMR spectroscopy, the size of the aliphatic carbon and heteroatoms was determined by infrared spectroscopy, and the 2D chemical structure of the TL-1 vitrinite was characterized. A large number of studies have shown that 13 C NMR can effectively characterize the 2D chemical molecular structure [29]. The 13 C chemical shift of the constructed vitrinite molecular structure was calculated using ACD/C NMR prediction software. After thousands of structural unit adjustments, the 13 C NMR spectrum of the constructed macromolecular structure and the 13 C NMR spectrum of the TL-1 vitrinite sample were matched (Figure 1a).
Geometric optimization and energy annealing of the 2D molecular structure (Figure 1b) were conducted based on the Materials Studio (Accelrys) software. Considering that the TL-1 sample's mid-rank coal structure order had not been completely formed [30], we used the AC module in the software to form a 3D unit box of 14 2D unit structures (Figure 1c). Geometric optimization was The experimental peak of 13 C NMR spectroscopy showed that there were three distinct peaks between 10-75 ppm, 100-160 ppm, and 175-225 ppm. Among these peaks, 10-75 ppm is classified as methyl and methylene in the aliphatic carbon region, and 100-160 ppm belongs to the aromatic carbon region, including protonated aromatic carbon and bridged aromatic carbon, and 175-225 ppm belongs to the carbonyl group containing oxygen. Since the sample was a mid-rank coal vitrinite sample, some of the aliphatic chains had fallen off, and the proportion of aromatic carbon was larger than that of aliphatic carbon, while the aromatic carbon mainly existed in the form of a benzene ring and a naphthalene ring. In the 3D molecular structure (Figure 1c), the relative atomic mass was 3358, and the unit cell size was 3.878 nm × 3.878 nm × 3.878 nm.

Monte Carlo Simulation
Using the sorption module in the Materials Studio (Accelrys) software, the adsorption behavior of CH4/CO2/N2 was studied by the grand canonical ensemble Monte Carlo (GCMC) simulation method. The interactions in the simulations included van der Waals interactions and Coulomb interactions, and the van der Waals interactions are described by the Lennard-Jones (LJ 9-6) potential energy model. The potential energy model describing the combined effects of the van der Waals forces and Coulomb forces is described in Equation 1 [31,32]: where Eij represents the interaction energy, kJ/mol; q i and q j are the atomic charges in the system, C, r ij is the distance between atom i and atom j, nm; is the dielectric constant, 8.854 × 10 -12 F/m; and the van der Waals force parameters, r o ij and ε ij , can be calculated based on r o (the atomic radius) and ε (the potential depth) (Equations 2 and 3, respectively) [31]: The experimental peak of 13 C NMR spectroscopy showed that there were three distinct peaks between 10-75 ppm, 100-160 ppm, and 175-225 ppm. Among these peaks, 10-75 ppm is classified as methyl and methylene in the aliphatic carbon region, and 100-160 ppm belongs to the aromatic carbon region, including protonated aromatic carbon and bridged aromatic carbon, and 175-225 ppm belongs to the carbonyl group containing oxygen. Since the sample was a mid-rank coal vitrinite sample, some of the aliphatic chains had fallen off, and the proportion of aromatic carbon was larger than that of aliphatic carbon, while the aromatic carbon mainly existed in the form of a benzene ring and a naphthalene ring. In the 3D molecular structure (Figure 1c), the relative atomic mass was 3358, and the unit cell size was 3.878 nm × 3.878 nm × 3.878 nm.

Monte Carlo Simulation
Using the sorption module in the Materials Studio (Accelrys) software, the adsorption behavior of CH 4 /CO 2 /N 2 was studied by the grand canonical ensemble Monte Carlo (GCMC) simulation method. The interactions in the simulations included van der Waals interactions and Coulomb interactions, and the van der Waals interactions are described by the Lennard-Jones (LJ 9-6) potential energy model. The potential energy model describing the combined effects of the van der Waals forces and Coulomb forces is described in Equation (1) [31,32]: where E ij represents the interaction energy, kJ/mol; q i and q j are the atomic charges in the system, C, r ij is the distance between atom i and atom j, nm; ε 0 is the dielectric constant, 8.854 × 10 −12 F/m; and the van der Waals force parameters, r o ij and ε ij , can be calculated based on r o (the atomic radius) and ε (the potential depth) (Equations (2) and (3), respectively) [31]: Considering the geological characteristics of the coal reservoir pressure, five temperature points of 283 K, 298 K, 318 K, 338 K, and 358 K were selected, respectively. The pressure range was between 0.1 and 10 MPa. To accurately calculate, the constant temperature and the constant pressure, the point-by-point calculation method was adopted, and the calculation was divided into eight pressure points. In this paper, the Compass II force field [33] was used to calculate the intermolecular force. The key energy equations are as follows (Equations (4)-(6)) [24]: where E total is the total energy, E Val and E Non are the valence electron energy and the nonbond energy, respectively, and E Van , E ele , and E hy are the van der Waals forces, electrostatic forces, and hydrogen bond energies, respectively. Based on this point of view, Gibbs defined the concept of excess adsorption, which is considers that not all gas in the adsorbed phase is on the surface of the adsorbent. The absolute adsorbed amount contains the adsorbed molecules in the form of the adsorbed phases, and in the form of the gas phase. The adsorbed molecules are distributed in the adsorption region, and the excess adsorption capacity only includes the adsorbed molecules in the form of the adsorbed phases ( Figure 2a) [31,34]. When the adsorption pressure is low, the excess adsorption capacity and absolute adsorption capacity are basically the same, due to the small difference between the bulk gas phase density and the adsorption phase density. However, when the experimental pressure increases, the larger difference between the two density values results in a larger difference between the excess adsorption capacity and the absolute adsorption capacity. According to the Gibbs concept, the excess adsorption capacity can be obtained by Equation (7) [35]: where N (excess) and N (absolute) are the excess adsorption and absolute adsorption capacities, respectively, mol/g; V P is the adsorbed phase volume, cm 3 ; ρ g is the density of the gas phase, which can be calculated by the SRK equation, g/cm 3 .
Considering the geological characteristics of the coal reservoir pressure, five temperature points of 283 K, 298 K, 318 K, 338 K, and 358 K were selected, respectively. The pressure range was between 0.1 and 10 MPa. To accurately calculate, the constant temperature and the constant pressure, the point-by-point calculation method was adopted, and the calculation was divided into eight pressure points. In this paper, the CompassⅡforce field [33] was used to calculate the intermolecular force.
The key energy equations are as follows (Equations 4-6) [24]: ENon=EVan+Eele+Ehy, EVal =EBo+EAn+ETor+EIn, where Etotal is the total energy, EVal and ENon are the valence electron energy and the nonbond energy, respectively, and EVan, Eele, and Ehy are the van der Waals forces, electrostatic forces, and hydrogen bond energies, respectively. Based on this point of view, Gibbs defined the concept of excess adsorption, which is considers that not all gas in the adsorbed phase is on the surface of the adsorbent. The absolute adsorbed amount contains the adsorbed molecules in the form of the adsorbed phases, and in the form of the gas phase. The adsorbed molecules are distributed in the adsorption region, and the excess adsorption capacity only includes the adsorbed molecules in the form of the adsorbed phases ( Figure  2a) [31,34]. When the adsorption pressure is low, the excess adsorption capacity and absolute adsorption capacity are basically the same, due to the small difference between the bulk gas phase density and the adsorption phase density. However, when the experimental pressure increases, the larger difference between the two density values results in a larger difference between the excess adsorption capacity and the absolute adsorption capacity. According to the Gibbs concept, the excess adsorption capacity can be obtained by Equation 7 [35]: where N excess and N absolute are the excess adsorption and absolute adsorption capacities, respectively, mol/g; V P is the adsorbed phase volume, cm 3 ; ρ g is the density of the gas phase, which can be calculated by the SRK equation, g/cm 3 . Figure 2. The excess adsorption capacity and the absolute adsorption capacity modified from Zhang et al [34]. Figure 2. The excess adsorption capacity and the absolute adsorption capacity modified from Zhang et al. [34].
In Figure 2, the area of Region a is the amount of adsorption that is completely present in the form of the adsorbed phase, that is, the amount of excess adsorption. Region b represents the amount of adsorption present in the adsorption zone in the form of the gas phase. Region c is the free gas present in the form of a gas phase. The sum of Regions a and b is the absolute adsorption capacity. The sum of Region of a, b, and c is the total amount of gas in the system, N, and the obtained result is the total gas content. Therefore, the excess adsorption capacity can be expressed as [35]: where N is the total gas content in the simulation, mol/g; V g is the gas phase volume, cm 3 ; V P is the free volume, cm 3 ; which can be obtained by selecting the corresponding gas molecules as probes using the Connolly Surface software tool, cm 3 .

Molecular Dynamics Simulation
A molecular dynamics simulation was carried out by using the dynamics tasks in the Focite module of the Material Studio software. In the molecular dynamics simulation, the geometrical optimization and energy optimization of the dominant adsorption configuration were carried out. The initial rate obeyed the Boltzmann random distribution, and the control of the temperature and pressure was performed by the Berendsen method and the Andersen thermal bath temperature control method, respectively. First, a 500 ps NVE ensemble was used for crystal relaxation to ensure the optimal structure; then, a 500 ps NVT ensemble was used to balance the system, and finally, a 500 ps NPT ensemble was used for the molecular dynamics simulations to ensure sufficient time to form a stable linear relationship between MSD and time.

Self-, Maxwell's, and Transport Diffusions
The diffusion coefficients were divided into the self-diffusion coefficient (D s ), the Maxwell's diffusion coefficient (D c ), and the transport diffusion coefficient (D t ), which are used to characterize the gas molecular mobility and transport in porous materials [36]. Self-diffusion is the Brownian motion of particles in thermodynamic equilibrium without the driving force, while the particle motion driven by a chemical potential or concentration gradient is transport diffusion. In molecular dynamics, the Einstein equation is used to calculate the self-diffusion coefficient [21,37]: where N is the number of diffused particles in the system, r i (t) is the position vector of a particle at time t, r i (0) is the position vector of the particle at the initial time, and t is the simulation time. According to Maxwell's theory, the Maxwell's diffusion coefficient can be obtained as follows [38]: The relationship between the self-diffusion and Maxwell's diffusion is as follows: where c is the gas concentration, mol/m 3 ; and V is the volume, m 3 .
where f is the degree of fugacity, MPa; according to the P-R equation [39]: where p is the pressure, MPa, z is the compression coefficient, a is a gravitational parameter, and b is the van der Waals covolume. The diffusion behavior of gas in coal molecules is also an activation process, that is, the relationship between the transport diffusion process and activation energy is where Ea is the diffusion activation energy, J·mol −1 ; D 0 is the diffusion prefactor, m 2 /s; R is the universal gas constant, J·mol −1 ·k −1 ; and T is the thermodynamic temperature, K. When gas molecules interact with coal molecules, strain deformation occurs along with the volume of the coal molecules, and the volumetric strain expansion rate is defined as [12] Swelling ratio where V i and V 0 are the volume before and after the volume strain, A 3 .

Pores in Coal Macromolecule
Previous studies have shown that the micropore volume and specific surface area account for a large proportion of the total pores, which is important for gas adsorption and diffusion [1]. Therefore, it is necessary to define the pores in the macromolecular structure. The same pores were measured using probe molecules with different radii, and the results were different. In gas adsorption studies, it is particularly important that the gas probe molecules detect the free pore volume. In this paper, each gas molecule operating as a probe molecule scanned the pores one by one ( Figure 3). Different pore sizes of probe molecules could detect different pore volumes. Some spaces inside the pores were very small, which large-diameter probe molecules could not detect, and so part of the pore volume was ignored by the probe molecules. However, this part of the space system could be detected when a probe molecule with a smaller radius was used. The surface of the van der Waals surface detected in the macromolecular structure was the micropore surface, and the micropore volume and surface area information could be obtained by detecting the van der Waals surface in the macromolecular structure.
using probe molecules with different radii, and the results were different. In gas adsorption studies, it is particularly important that the gas probe molecules detect the free pore volume. In this paper, each gas molecule operating as a probe molecule scanned the pores one by one ( Figure 3). Different pore sizes of probe molecules could detect different pore volumes. Some spaces inside the pores were very small, which large-diameter probe molecules could not detect, and so part of the pore volume was ignored by the probe molecules. However, this part of the space system could be detected when a probe molecule with a smaller radius was used. The surface of the van der Waals surface detected in the macromolecular structure was the micropore surface, and the micropore volume and surface area information could be obtained by detecting the van der Waals surface in the macromolecular structure. Taking the CO 2 (R = 0.165 nm) probe molecule as an example, the sample pores were calculated ( Figure 4). The pores in the macromolecular structure were formed by vacancies in the absence of atoms ( Figure 4b). In general, the pore wall atoms were mainly composed of aliphatic side chains and a small amount of oxygen-containing functional groups. Part of the side chain of the sample had fallen off in this mid-rank coal, resulting in the formation of some oxygen-containing functional groups in the pore structure. Figure 4c shows that there were a large number of pores in the model with diverse and irregular shapes. Most were elliptic, and some pores were connected by several small pores to form a pore network, which provided an important channel for gas transport and desorption.  Taking the CO2 (R = 0.165 nm) probe molecule as an example, the sample pores were calculated ( Figure 4). The pores in the macromolecular structure were formed by vacancies in the absence of atoms ( Figure 4b). In general, the pore wall atoms were mainly composed of aliphatic side chains and a small amount of oxygen-containing functional groups. Part of the side chain of the sample had fallen off in this mid-rank coal, resulting in the formation of some oxygen-containing functional groups in the pore structure. Figure 4c shows that there were a large number of pores in the model with diverse and irregular shapes. Most were elliptic, and some pores were connected by several small pores to form a pore network, which provided an important channel for gas transport and desorption. Different sizes of gas probe molecules can detect the pore surface area and volume. In the same pore space, the larger the probe molecule radius, the smaller the measured surface area is. If the pore inlet size is too small, the molecular probe with a diameter that is larger than the pore diameter cannot access the pore; however, the molecular probe with a diameter that is smaller than the pore diameter probe molecules can enter the pore ( Figure 5). This results in differences in the detected pore volume. Therefore, according to different research purposes, appropriately sized probe molecules were selected. There are two types of pores in the macromolecular structure. The first type is the accessible pores in Figure 5b, which means that gas can enter the pore from the outside; the other type is the Different sizes of gas probe molecules can detect the pore surface area and volume. In the same pore space, the larger the probe molecule radius, the smaller the measured surface area is. If the pore inlet size is too small, the molecular probe with a diameter that is larger than the pore diameter cannot access the pore; however, the molecular probe with a diameter that is smaller than the pore diameter probe molecules can enter the pore ( Figure 5). This results in differences in the detected pore volume. Therefore, according to different research purposes, appropriately sized probe molecules were selected. Taking the CO2 (R = 0.165 nm) probe molecule as an example, the sample pores were calculated ( Figure 4). The pores in the macromolecular structure were formed by vacancies in the absence of atoms ( Figure 4b). In general, the pore wall atoms were mainly composed of aliphatic side chains and a small amount of oxygen-containing functional groups. Part of the side chain of the sample had fallen off in this mid-rank coal, resulting in the formation of some oxygen-containing functional groups in the pore structure. Figure 4c shows that there were a large number of pores in the model with diverse and irregular shapes. Most were elliptic, and some pores were connected by several small pores to form a pore network, which provided an important channel for gas transport and desorption. Different sizes of gas probe molecules can detect the pore surface area and volume. In the same pore space, the larger the probe molecule radius, the smaller the measured surface area is. If the pore inlet size is too small, the molecular probe with a diameter that is larger than the pore diameter cannot access the pore; however, the molecular probe with a diameter that is smaller than the pore diameter probe molecules can enter the pore ( Figure 5). This results in differences in the detected pore volume. Therefore, according to different research purposes, appropriately sized probe molecules were selected. There are two types of pores in the macromolecular structure. The first type is the accessible pores in Figure 5b, which means that gas can enter the pore from the outside; the other type is the inaccessible pores (the gray part of Figure 5a) which is mainly located inside the molecular structure. For an inaccessible pore, the atomic contact with the external environment does not occur, and the There are two types of pores in the macromolecular structure. The first type is the accessible pores in Figure 5b, which means that gas can enter the pore from the outside; the other type is the inaccessible pores (the gray part of Figure 5a) which is mainly located inside the molecular structure. For an inaccessible pore, the atomic contact with the external environment does not occur, and the external gas cannot enter the pore. It is also difficult for gas molecules to enter such pores during an experiment.
The pores in the model were measured using He, N 2 , CO 2 , and CH 4 probe molecules with radii of 0.13 nm, 0.152 nm, 0.165 nm, and 0.19 nm, respectively. The measured pore volumes were 0.053 mL/g, 0.039 mL/g, 0.03426 mL/g, and 0.026 mL/g, respectively. Figure 6 shows that the pore volumes measured by the different sizes of probe molecules vary greatly. As the molecular radius of helium gas is smaller than those of N 2 , CO 2 , and CH 4 , the pores that can be detected by helium gas may not be accessible to N 2 , CO 2 , and CH 4 , which means that helium gas has a stronger pore-measuring capacity. The pore volume measured by helium was approximately twice the pore volume of the pores measured by the CH 4 . As shown in Figure 6b,c, the pore morphology measured by N 2 was roughly similar to that measured by CO 2 , except for some differences in the small pores. This phenomenon was because the CO 2 -detected pores were mostly detectable by N 2 , but that the CO 2 molecules could not enter the smaller spaces at the edges. In the methane adsorption experiment, the helium gas volume was measured first, resulting in a difference in the calculation of the adsorption capacity. Therefore, selecting the appropriate size of the probe molecule in the simulations is an important prerequisite for accurately calculating the excess adsorption capacity.
Energies 2019, 12, x FOR PEER REVIEW 9 of 21 external gas cannot enter the pore. It is also difficult for gas molecules to enter such pores during an experiment.
The pores in the model were measured using He, N2, CO2, and CH4 probe molecules with radii of 0.13 nm, 0.152 nm, 0.165 nm, and 0.19 nm, respectively. The measured pore volumes were 0.053 ml/g, 0.039 ml/g, 0.03426 ml/g, and 0.026 ml/g, respectively. Figure 6 shows that the pore volumes measured by the different sizes of probe molecules vary greatly. As the molecular radius of helium gas is smaller than those of N2, CO2, and CH4, the pores that can be detected by helium gas may not be accessible to N2, CO2, and CH4, which means that helium gas has a stronger pore-measuring capacity. The pore volume measured by helium was approximately twice the pore volume of the pores measured by the CH4. As shown in Figure 6 b and c, the pore morphology measured by N2 was roughly similar to that measured by CO2, except for some differences in the small pores. This phenomenon was because the CO2-detected pores were mostly detectable by N2, but that the CO2 molecules could not enter the smaller spaces at the edges. In the methane adsorption experiment, the helium gas volume was measured first, resulting in a difference in the calculation of the adsorption capacity. Therefore, selecting the appropriate size of the probe molecule in the simulations is an important prerequisite for accurately calculating the excess adsorption capacity.

Adsorption Results
The energy changes and adsorption configurations for different amounts of CH4, CO2, and N2 are shown in Figure 7a. The Eele basically did not change with the increasing number of adsorption, while the interaction energy between the EVan and the Ebon decreased with the increasing number of the adsorbates, indicating that the vitrinite of coal attracted the adsorbates in this adsorption process. Then the Eele decreased, indicating that the vitrinite no longer attracted the adsorbates and began to exclude them. When Ebon reached the lowest value, the amount of adsorbate reached the maximum. The adsorption of vitrinite reaches the saturation state after adsorbing 28CH4, 33CO2, and 17N2 (Figure 7b), and the Ebon energies of the three gases reaching the maximum adsorption sites were 237.876, -287.417 kcal/mol, and -96.057 kcal/mol, respectively, indicating that the order of the adsorption capacity of the macromolecular structures for the three gases was CO2 > CH4 > N2. In addition, as shown in Figure 7b, the dominant adsorption sites occupied by CO2 covered the CH4/N2 sites, indicating that CO2 can compete with CH4 and N2 for adsorption by preferentially occupying the adsorption sites. The gas diffusion characteristics were studied by the molecular dynamics method, using the dominant adsorption sites for the three gases as the initial configuration.

Adsorption Results
The energy changes and adsorption configurations for different amounts of CH 4 , CO 2 , and N 2 are shown in Figure 7a. The E ele basically did not change with the increasing number of adsorption, while the interaction energy between the E Van and the E bon decreased with the increasing number of the adsorbates, indicating that the vitrinite of coal attracted the adsorbates in this adsorption process. Then the E ele decreased, indicating that the vitrinite no longer attracted the adsorbates and began to exclude them. When E bon reached the lowest value, the amount of adsorbate reached the maximum. The adsorption of vitrinite reaches the saturation state after adsorbing 28CH 4 , 33CO 2 , and 17N 2 (Figure 7b), and the E bon energies of the three gases reaching the maximum adsorption sites were 237.876, −287.417 kcal/mol, and −96.057 kcal/mol, respectively, indicating that the order of the adsorption capacity of the macromolecular structures for the three gases was CO 2 > CH 4 > N 2 . In addition, as shown in Figure 7b, the dominant adsorption sites occupied by CO 2 covered the CH 4 /N 2 sites, indicating that CO 2 can compete with CH 4 and N 2 for adsorption by preferentially occupying the adsorption sites. The gas diffusion characteristics were studied by the molecular dynamics method, using the dominant adsorption sites for the three gases as the initial configuration.  Figure 8 shows the adsorption isotherm curves of the average loading per cell and excess adsorption capacity of CO2, CH4, and N2 under different temperature and pressure conditions. The average loading per cell refers to the number of gas molecules adsorbed by the 3D coal molecular structure unit constructed above. The adsorption capacity of CH4, CO2, and N2 decreased with increasing temperature, which is due to the adsorption behavior. During the exothermic process, the temperature rise increased the kinetic energy of the gas molecules, which increased the rate of movement, making it easier for gas molecules to escape the nucleus of the vitrinite, and thus the van der Waals function could be reduced. The adsorption capacity of CH4 and N2 increased with increasing pressure; specifically, the adsorption rate of CH4 and N2 increased rapidly when P < 5MPa at the low-pressure stage. In the middle stage of the curve, the adsorption capacity increased slowly for P > 6 MPa until P < 8 MPa. When P > 8 MPa at the high-pressure stage, the adsorption curves of the two groups tended to be stable and reached the saturated adsorption state at this time, regardless of whether the average loading per unit or the excess adsorption capacity was larger for CH4 and N2 at the same temperature or pressure. Compared with CH4 and N2, CO2 has more preferential adsorption characteristics, which is consistent with the above simulation results for the saturated adsorption sites. This finding is closely related to the difference in the interaction between each gas and the coal molecular structure and the molecular dynamics diameter and polarity. Compared with the excess adsorption of CH4 and N2, the change of excess adsorption capacity of CO2 was different. Under the low pressure, the excess adsorption capacity showed an increasing trend. When the pressure was larger than 5 MPa, the excess adsorption capacity reached its maximum and began to decrease. This is due to the excess adsorption capacity of CO2 being related to the pore volume. Under high pressure, the volume phase was still compressed, the pores are filled up, and the density of the adsorbed phase leveled off. Once the bulk phase density was higher than the adsorbed phase density, the excess adsorption capacity showed a decreasing trend.  Figure 8 shows the adsorption isotherm curves of the average loading per cell and excess adsorption capacity of CO 2 , CH 4 , and N 2 under different temperature and pressure conditions. The average loading per cell refers to the number of gas molecules adsorbed by the 3D coal molecular structure unit constructed above. The adsorption capacity of CH 4 , CO 2 , and N 2 decreased with increasing temperature, which is due to the adsorption behavior. During the exothermic process, the temperature rise increased the kinetic energy of the gas molecules, which increased the rate of movement, making it easier for gas molecules to escape the nucleus of the vitrinite, and thus the van der Waals function could be reduced. The adsorption capacity of CH 4 and N 2 increased with increasing pressure; specifically, the adsorption rate of CH 4 and N 2 increased rapidly when P < 5 MPa at the low-pressure stage. In the middle stage of the curve, the adsorption capacity increased slowly for P > 6 MPa until P < 8 MPa. When P > 8 MPa at the high-pressure stage, the adsorption curves of the two groups tended to be stable and reached the saturated adsorption state at this time, regardless of whether the average loading per unit or the excess adsorption capacity was larger for CH 4 and N 2 at the same temperature or pressure. Compared with CH 4 and N 2 , CO 2 has more preferential adsorption characteristics, which is consistent with the above simulation results for the saturated adsorption sites. This finding is closely related to the difference in the interaction between each gas and the coal molecular structure and the molecular dynamics diameter and polarity. Compared with the excess adsorption of CH 4 and N 2 , the change of excess adsorption capacity of CO 2 was different. Under the low pressure, the excess adsorption capacity showed an increasing trend. When the pressure was larger than 5 MPa, the excess adsorption capacity reached its maximum and began to decrease. This is due to the excess adsorption capacity of CO 2 being related to the pore volume. Under high pressure, the volume phase was still compressed, the pores are filled up, and the density of the adsorbed phase leveled off. Once the bulk phase density was higher than the adsorbed phase density, the excess adsorption capacity showed a decreasing trend. A comparison of the CH4 adsorption simulation at 298K and experimental data [28] at 303.15 K is shown in Figure 9. Overall, the simulation results are consistent with the general trend of the experimental data changes with increasing pressure, and the adsorption capacity reached a maximum when the pressure reached the saturation pressure value. The maximum value of the simulated adsorption was slightly smaller than the experimental value. When P < 5 MPa, the simulated values of the adsorption capacity and adsorption rate were significantly higher than the experimental values. When P > 5 MPa, the relative differences of the two variables gradually decreased, and the difference was the smallest when the maximum adsorption capacity was reached. The simulated maximum adsorption capacity was 13.75 cm 3 /g, and the experimental value was 14.78 cm 3 /g. There are many reasons for the difference at the low-pressure stage. Firstly, the experimental sample may have contained ash and other mineral components inside the components of the handselected vitrinite group. Compared with the pure organic matter vitrinite, the macromolecular structure model constructed in this paper had a relatively low adsorption capacity. The presence of inorganic mineral components may have hindered the adsorption process. The smaller difference under high-pressure conditions may be because that the adsorption of minerals and ash is more   A comparison of the CH 4 adsorption simulation at 298K and experimental data [28] at 303.15 K is shown in Figure 9. Overall, the simulation results are consistent with the general trend of the experimental data changes with increasing pressure, and the adsorption capacity reached a maximum when the pressure reached the saturation pressure value. The maximum value of the simulated adsorption was slightly smaller than the experimental value. When P < 5 MPa, the simulated values of the adsorption capacity and adsorption rate were significantly higher than the experimental values. When P > 5 MPa, the relative differences of the two variables gradually decreased, and the difference was the smallest when the maximum adsorption capacity was reached. The simulated maximum adsorption capacity was 13.75 cm 3 /g, and the experimental value was 14.78 cm 3 /g. There are many reasons for the difference at the low-pressure stage. Firstly, the experimental sample may have contained ash and other mineral components inside the components of the hand-selected vitrinite group. Compared with the pure organic matter vitrinite, the macromolecular structure model constructed in this paper had a relatively low adsorption capacity. The presence of inorganic mineral components may have hindered the adsorption process. The smaller difference under high-pressure conditions may be because that the adsorption of minerals and ash is more obvious under low-pressure conditions, and the adsorption is not exerted under high-pressure conditions, while the adsorption of vitrinite plays a dominant role. Secondly, the model mainly simulates the micropore structure in the vitrinite macromolecular structure, and does not consider the structure of mesopores and macropores while the experimental sample contained both micropores, mesopores, and macropores structures. Since the adsorption mainly occurs in the micropores structure under low-pressure conditions, the micropore filling effect occurs in a short time, so the simulated value relative to the experimental adsorption capacity was more notable under low-pressure conditions. Finally, in this paper, in order to calculate more realistically and effectively and obtain the excess adsorption capacity, we selected the appropriate gas as the probe molecule to measure the free pore volume. For example, the free pore volume of CH 4 was measured by a CH 4 (r = 0.13 nm) molecular probe. However, in the experimental process, helium gas (r = 0.19 nm) was selected as the probe molecule to measure the free pore volume, which meant that the free pore volume measured relative to the gas molecules was small, resulting in relatively small differences in the experimental and calculated values. obvious under low-pressure conditions, and the adsorption is not exerted under high-pressure conditions, while the adsorption of vitrinite plays a dominant role. Secondly, the model mainly simulates the micropore structure in the vitrinite macromolecular structure, and does not consider the structure of mesopores and macropores while the experimental sample contained both micropores, mesopores, and macropores structures. Since the adsorption mainly occurs in the micropores structure under low-pressure conditions, the micropore filling effect occurs in a short time, so the simulated value relative to the experimental adsorption capacity was more notable under low-pressure conditions. Finally, in this paper, in order to calculate more realistically and effectively and obtain the excess adsorption capacity, we selected the appropriate gas as the probe molecule to measure the free pore volume. For example, the free pore volume of CH4 was measured by a CH4 (r = 0.13 nm) molecular probe. However, in the experimental process, helium gas (r = 0.19 nm) was selected as the probe molecule to measure the free pore volume, which meant that the free pore volume measured relative to the gas molecules was small, resulting in relatively small differences in the experimental and calculated values.

Self-, Corrected, and Transport Diffusion Coefficients
Transport diffusion means that the attenuation of the mass transport and density fluctuations in the system occurs more macroscopically than that of self-diffusion, and transport diffusion is more closely related to the application of modern technology and engineering. An increasing number of scholars have paid attention to transport diffusion [22]. To investigate the influence of different concentrations of adsorbed gas on the diffusion behavior under saturated conditions, CH4 and CO2 were taken as the control group with the same number of dominant saturated adsorption sites of N2. The gas diffusion coefficients were obtained by molecular dynamics simulations and compared with the diffusion experimental values using the same coal rank, temperature, and pressure in the paper of Saghafi. The result indicates that CH4 and CO2 were relatively close to the experimental values in the literature.

Self-, Corrected, and Transport Diffusion Coefficients
Transport diffusion means that the attenuation of the mass transport and density fluctuations in the system occurs more macroscopically than that of self-diffusion, and transport diffusion is more closely related to the application of modern technology and engineering. An increasing number of scholars have paid attention to transport diffusion [22]. To investigate the influence of different concentrations of adsorbed gas on the diffusion behavior under saturated conditions, CH 4 and CO 2 were taken as the control group with the same number of dominant saturated adsorption sites of N 2 . The gas diffusion coefficients were obtained by molecular dynamics simulations and compared with the diffusion experimental values using the same coal rank, temperature, and pressure in the paper of Saghafi. The result indicates that CH 4 and CO 2 were relatively close to the experimental values in the literature.
The diffusion coefficients of CO 2 , CH 4 , and N 2 simulated by Song et al. [11] in low-rank coal vitrinite molecules were 27.83-79.56 × 10 −11 , 16.23-38.45 × 10 −11 , and 42.37-86.45 × 10 −11 m 2 /s, respectively. Zhao et al. [12] obtained the diffusion coefficients of CO 2 and CH 4 based on the Wiser coal molecular simulations, which were 77.3 × 10 −11 m 2 /s and 36.3 × 10 −11 m 2 /s, respectively. However, these values are much smaller than those found by Hu et al. [36] based on the Wiser coal molecular model; they obtained CO 2 and CH 4 diffusion coefficients of 1.2 × 10 −9 and 1.0 × 10 −9 m 2 /s, respectively. The reasons for the difference mainly include two aspects. On the one hand, based on the different coal molecular models, there was a large difference in the number ratio of the oxygen-containing functional groups and aliphatic side chains. Previous studies have shown that different types of chemical groups have complex effects on the gas adsorption capacity, which is mainly related to the interaction between the groups and gases and the distribution of the groups in space. A study by Yu et al. [1], based on the slite pore model, found that the order of the adsorption capacity of the groups was aromatic ring > aliphatic side chain > carboxyl group > hydroxyl group when adsorbing a single methane molecule, but this provided more gas when the carboxyl group and the hydroxyl group formed the pore wall atoms. At this point, the adsorption capacity was greater than those of the aromatic ring and methyl group. Therefore, the difference in the number of chemical groups in each model is the main reason for the difference in the simulated diffusion results. On the other hand, the different force field selections and parameter settings in the simulations will also lead to large differences in the simulation results, as reported by Song et al. [11], who based their work on the Dreiding force field, and by Hu et al. [36], who used the Compass II force field. In molecular dynamics, the self-diffusion requires a long simulation time to ensure a stable linear relationship between MSD and time. Therefore, the length of the simulation will also cause a difference in the results. In this paper, the simulation time was 500 ps for accurate calculation results.
As summarized in Table 2, the diffusion coefficients of CO 2 , CH 4 , and N 2 were ordered as D S < D c < D t , indicating that the diffusion capacity increases with an increasing observation scale from the microscale to the macroscale, in which the diffusion coefficient decreases with the increase in adsorption capacity, indicating that the increase in gas concentration is not conducive to diffusion. The increase in concentration corresponds to a high pressure and high gas density, and the intermolecular force is enhanced to hinder the tendency of gas to diffuse out of the pore. In the saturated adsorption state, the self-diffusion and transport diffusion were both ordered as N 2 > CO 2 > CH 4 . However, when there were 17 adsorption sites, the order was D S (N 2 ) > D S CO 2 ) > D S (CH 4 ) and D t (CO 2 ) > D t (N 2 ) > D t (CH 4 ), respectively, which are opposite to the order of the gas molecular dynamics diameter, namely σ CH 4 (0.38 nm) > σ N 2 (0.36 nm) > σ CO 2 (0.33 nm). The molecular dynamics diameter is a sensitive measure of the ability to move in a highly restricted environment [40]. CH 4 cannot enter pores with pore diameters smaller than 0.38 nm, while CO 2 can enter and exit. The smaller the kinetic diameter is, the wider the range of the pore diameter is. In addition to the difference in diameter, this phenomenon is related to the tensile morphology of gas molecules. CO 2 is smaller and more linear than CH 4 and is more accessible to the restricted pore space volume. Figure 10 shows the difference in the adsorption sites of the five CO 2 , CH 4 , and N 2 adsorbed molecules configurations in the model. As shown in Figure 10a-c, when the three gases existed individually, CO 2 was distributed in the narrow pore throats according to the linear dominance, and the CH 4 tetrahedral structure was distributed in the interior of the large voids; N 2 may have been distributed in the interior of the small pores due to its triple bond form. As shown in Figure 10d-f, when CH 4 and CO 2 were in the same pore, CO 2 was closer to the pore throat position, and CH 4 was inside the pore. When N 2 and CO 2 were in the same pore, CO 2 was closer to the pore throat position, and N 2 was distributed in the pores. The molecular dynamics diameter and morphology were important factors affecting adsorption and diffusion. In addition, Larsen [41] proposed another theory whereby the CO 2 adsorption isotope is greater than that of nonpolar hydrocarbons, such as CH 4 , which is responsible for the faster diffusion of CO 2 in coal. the CO2 adsorption isotope is greater than that of nonpolar hydrocarbons, such as CH4, which is responsible for the faster diffusion of CO2 in coal. Einstein's diffusion law can be used to determine the migration path of molecules in the different micropores by the diffusion trajectory of a single gas molecule, which can clearly discriminate the diffusion path and displacement of gas molecules from one pore to another. The gas molecules are constantly moving inside the coal molecules. In most cases, small-amplitude vibrations occur in the pores. When the gas molecule vibrates to a certain position in the pores, the molecular structure of the vitrinite is twisted or moved, and the gas molecule is directed to the other pores. Jumping provides a channel and again results in a small amplitude vibration in the other aperture. The diffusion trajectories and displacements of the five stable configurations of CO2, CH4, and N2 are shown in Figure 11. As shown by the comparison of Figure 11a-c, when the capacity of the adsorbed gas was constant, the displacement change intensity was ordered as N2 > CO2 > CH4, and the simulated jump time between 0 and 50 ps was because the system energy that had not been completely stabilized. When N2 shifted within 100 ps, a cross-pore jump occurred, and the CO2 and CH4 cross-pore movement points appeared at 200 and 450 ps, respectively, which are consistent with the calculation results listed in Table 2 above. As shown in Figure 11d, the contrast between e and a, b, and c increased with the increase in the concentration of the adsorbed gas, and the vibration amplitude of the gas decreased significantly. The time of the large-jump point was shortened. It should be noted that the large-jump point at this time was not the cross-pore jump. Instead, relatively large-amplitude vibrations occurred in the same pore. Einstein's diffusion law can be used to determine the migration path of molecules in the different micropores by the diffusion trajectory of a single gas molecule, which can clearly discriminate the diffusion path and displacement of gas molecules from one pore to another. The gas molecules are constantly moving inside the coal molecules. In most cases, small-amplitude vibrations occur in the pores. When the gas molecule vibrates to a certain position in the pores, the molecular structure of the vitrinite is twisted or moved, and the gas molecule is directed to the other pores. Jumping provides a channel and again results in a small amplitude vibration in the other aperture. The diffusion trajectories and displacements of the five stable configurations of CO 2 , CH 4 , and N 2 are shown in Figure 11. As shown by the comparison of Figure 11a-c, when the capacity of the adsorbed gas was constant, the displacement change intensity was ordered as N 2 > CO 2 > CH 4 , and the simulated jump time between 0 and 50 ps was because the system energy that had not been completely stabilized. When N 2 shifted within 100 ps, a cross-pore jump occurred, and the CO 2 and CH 4 cross-pore movement points appeared at 200 and 450 ps, respectively, which are consistent with the calculation results listed in Table 2 above. As shown in Figure 11d, the contrast between e and a, b, and c increased with the increase in the concentration of the adsorbed gas, and the vibration amplitude of the gas decreased significantly. The time of the large-jump point was shortened. It should be noted that the large-jump point at this time was not the cross-pore jump. Instead, relatively large-amplitude vibrations occurred in the same pore. the CO2 adsorption isotope is greater than that of nonpolar hydrocarbons, such as CH4, which is responsible for the faster diffusion of CO2 in coal. Einstein's diffusion law can be used to determine the migration path of molecules in the different micropores by the diffusion trajectory of a single gas molecule, which can clearly discriminate the diffusion path and displacement of gas molecules from one pore to another. The gas molecules are constantly moving inside the coal molecules. In most cases, small-amplitude vibrations occur in the pores. When the gas molecule vibrates to a certain position in the pores, the molecular structure of the vitrinite is twisted or moved, and the gas molecule is directed to the other pores. Jumping provides a channel and again results in a small amplitude vibration in the other aperture. The diffusion trajectories and displacements of the five stable configurations of CO2, CH4, and N2 are shown in Figure 11. As shown by the comparison of Figure 11a-c, when the capacity of the adsorbed gas was constant, the displacement change intensity was ordered as N2 > CO2 > CH4, and the simulated jump time between 0 and 50 ps was because the system energy that had not been completely stabilized. When N2 shifted within 100 ps, a cross-pore jump occurred, and the CO2 and CH4 cross-pore movement points appeared at 200 and 450 ps, respectively, which are consistent with the calculation results listed in Table 2 above. As shown in Figure 11d, the contrast between e and a, b, and c increased with the increase in the concentration of the adsorbed gas, and the vibration amplitude of the gas decreased significantly. The time of the large-jump point was shortened. It should be noted that the large-jump point at this time was not the cross-pore jump. Instead, relatively large-amplitude vibrations occurred in the same pore.

Effect of the Temperature on Gas Diffusion
The laws of the self-diffusion and transport diffusion coefficients of CH4, CO2, and N2 with the temperature are shown in Figure 12 (the fixed pressure was 1 MPa). The three self-diffusion and transport diffusion coefficients increased with temperature, which was the same growth trend as the Charrière et al. [23] unique model simulation results. As shown in Figure 12a-c, the influence of the temperature on the self-diffusion coefficient can be divided into three stages-rapid rise, slow rise, and rapid rise again. The initial rapid rise is due to the increase in the molecular vibration amplitude of a single gas molecule in a certain pore with increasing temperature. When the molecule moves to the position of the orifice, the resistance of the pore wall increases, and the diffusion rate is decreases, Finally, the self-diffusion coefficient rapidly rises again due to the molecule from a pore jumping to another pore and continuing to make a large vibration. The effect of the temperature on transport diffusion is manifested in the interaction force and adsorption affinity between gas and coal. Compared with the simulation results of Charrière et al. [23], which were based on the unique model with a uniform particle size of 0.3 nm, the results in Figure 12d,e are several orders of magnitude larger since Charrière is based on the unique model with a uniform particle size of 0.3 nm. On one hand, the vitrinite molecular model is based on the atomic level and combined with intermolecular interaction. On the other hand, the unique model is a graphite slit pore composed of an aromatic ring,

Effect of the Temperature on Gas Diffusion
The laws of the self-diffusion and transport diffusion coefficients of CH 4 , CO 2 , and N 2 with the temperature are shown in Figure 12 (the fixed pressure was 1 MPa). The three self-diffusion and transport diffusion coefficients increased with temperature, which was the same growth trend as the Charrière et al. [23] unique model simulation results. As shown in Figure 12a-c, the influence of the temperature on the self-diffusion coefficient can be divided into three stages-rapid rise, slow rise, and rapid rise again. The initial rapid rise is due to the increase in the molecular vibration amplitude of a single gas molecule in a certain pore with increasing temperature. When the molecule moves to the position of the orifice, the resistance of the pore wall increases, and the diffusion rate is decreases, Finally, the self-diffusion coefficient rapidly rises again due to the molecule from a pore jumping to another pore and continuing to make a large vibration. The effect of the temperature on transport diffusion is manifested in the interaction force and adsorption affinity between gas and coal. Compared with the simulation results of Charrière et al. [23], which were based on the unique model with a uniform particle size of 0.3 nm, the results in Figure 12d,e are several orders of magnitude larger since Charrière is based on the unique model with a uniform particle size of 0.3 nm. On one hand, the vitrinite molecular model is based on the atomic level and combined with intermolecular interaction. On the other hand, the unique model is a graphite slit pore composed of an aromatic ring, and the vitrinite molecular model contains different types of chemical groups. The adsorption order of the groups mentioned above was aromatic ring > aliphatic side chain > carboxyl group > hydroxyl group. Therefore, the diffusion simulation results in this paper were relatively large. Compared with the unique model, the vitrinite model can better reflect the heterogeneity of the pore structure.
Energies 2019, 12, x FOR PEER REVIEW 16 of 21 and the vitrinite molecular model contains different types of chemical groups. The adsorption order of the groups mentioned above was aromatic ring > aliphatic side chain > carboxyl group > hydroxyl group. Therefore, the diffusion simulation results in this paper were relatively large. Compared with the unique model, the vitrinite model can better reflect the heterogeneity of the pore structure. The migration and diffusion of gas in coal molecules is an activation process in accordance with the Arrhenius law [40]. As shown in Figure 13, ln(D t ) has a good linear relationship with 1/T (R 2 = 0.69-0.913), which also conforms to Equation 18. From the slope, a set of gas diffusion activation energies at different temperatures can be calculated. The diffusion activation energy also indicates the dependence of diffusion on the temperature, which can characterize the difficulty of migration and diffusion. The lower the activation energy is, the smaller the migration barrier is that needs to be overcome, and the faster the diffusion process is. The order of the activation energy of CH4, CO2, and N2 in the adsorption saturated state, as shown in Figure 13, was CH4 (27.388 kJ/mol) > CO2 (11.832 kJ/mol) > N2 (10.396 kJ/mol), which indicates that CO2 and N2 are more prone to diffusion behavior during diffusion than CH4.
Studies [20,23,29,41] have shown that the diffusion process of different concentrations of gas molecules in coal molecules causes coal volume swelling. As shown in Figure 13d-f, the volume swelling ratio increased with increasing temperature, independent of the gas concentration, indicating that the temperature promotes a swelling balance. On the one hand, the relationship between the swelling ratio and temperature is determined by the decrease in the Gibbs free energy required for the chemical thermodynamic processes (∆G < 0 and ∆Gm = ∆Hm -T ∆Sm). As the adsorption entropy increases during the diffusion process (∆S > 0), the temperature rise is favorable for ∆G < 0, resulting in an increase in the swelling ratio. On the other hand, some studies [42] have The migration and diffusion of gas in coal molecules is an activation process in accordance with the Arrhenius law [42]. As shown in Figure 13, ln(D t ) has a good linear relationship with 1/T (R 2 = 0.69-0.913), which also conforms to Equation (18). From the slope, a set of gas diffusion activation energies at different temperatures can be calculated. The diffusion activation energy also indicates the dependence of diffusion on the temperature, which can characterize the difficulty of migration and diffusion. The lower the activation energy is, the smaller the migration barrier is that needs to be overcome, and the faster the diffusion process is. The order of the activation energy of CH 4 , CO 2 , and N 2 in the adsorption saturated state, as shown in Figure 13, was CH 4 (27.388 kJ/mol) > CO 2 (11.832 kJ/mol) > N 2 (10.396 kJ/mol), which indicates that CO 2 and N 2 are more prone to diffusion behavior during diffusion than CH 4 .
Studies [20,23,29,43] have shown that the diffusion process of different concentrations of gas molecules in coal molecules causes coal volume swelling. As shown in Figure 13d-f, the volume swelling ratio increased with increasing temperature, independent of the gas concentration, indicating that the temperature promotes a swelling balance. On the one hand, the relationship between the swelling ratio and temperature is determined by the decrease in the Gibbs free energy required for the chemical thermodynamic processes (∆G < 0 and ∆Gm = ∆Hm -T ∆Sm). As the adsorption entropy increases during the diffusion process (∆S > 0), the temperature rise is favorable for ∆G < 0, resulting in an increase in the swelling ratio. On the other hand, some studies [41] have shown that the partial dissolution of CO 2 in the adsorption process leads to the rearrangement effect of the coal molecular structure, which is reflected in the macroscopic volume swelling, and that, with a rise in temperature, the solubility of CO 2 increases and the volume changes more notably. Compared with CO 2 , CH 4 is mainly present in the adsorbed state, and so the volume swelling ratio is lower than for CO 2 . The higher swelling ratio of N 2 at the same adsorption capacity of CH 4 and CO 2 may be related to the strength of the interaction between the triple bonds and the coal molecules. In addition, Figure 13 shows that the volume swelling ratio increased as the capacity of the gas increased. shown that the partial dissolution of CO2 in the adsorption process leads to the rearrangement effect of the coal molecular structure, which is reflected in the macroscopic volume swelling, and that, with a rise in temperature, the solubility of CO2 increases and the volume changes more notably. Compared with CO2, CH4 is mainly present in the adsorbed state, and so the volume swelling ratio is lower than for CO2. The higher swelling ratio of N2 at the same adsorption capacity of CH4 and CO2 may be related to the strength of the interaction between the triple bonds and the coal molecules. In addition, Figure 13 shows that the volume swelling ratio increased as the capacity of the gas increased.

Effect of the Pressure on Gas Diffusion
At present, the influence of the pressure on gas diffusion is controversial. Based on the different models, there are notable differences in the diffusion law with the pressure. Based on the unique model, Charrière et al. [23] found that the diffusion coefficient increased with increasing pressure. The results of Cui et al. [4] based on the bidisperse model, show that the diffusion rate of gas in coal is significantly reduced with the increase of pressure due to the high shrinkage of the micropores. In this paper, based on the vitrinite microporous model, the molecular dynamics simulation results of the diffusion coefficient with the pressure changes are shown in Figure 14. The results show that the self-diffusion and transport diffusion first increased and then decreased with increasing pressure, independent of the gas type and concentration. Additionally, the peak pressures of CH4, CO2, and N2 were all fixed values of 3 MPa. Before the peak pressure, the increase in pressure aggravated the collision between the gas molecules or between the gas molecules and the vitreous group molecules, leading to an increase in the diffusion coefficient. The decrease in the diffusion coefficient after the peak pressure was caused by the strong collision between the gas molecules and the molecules and the interaction with coal adsorption swelling. Figure 14d and e compare the transport diffusion coefficients calculated with those of Zhao et al. [12], that were based on the Wiser model. Generally, the trend of the change was consistently rising at first and then decreasing, but the peak pressures of   Figure 13. Relationship between ln(D t ) and 1/T for CH 4 (a), CO 2 (b), and N 2 (c) diffusions in vitrinite macromolecules, and the swelling ratios for CH 4 (d), CO 2 (e), and N 2 (f).

Effect of the Pressure on Gas Diffusion
At present, the influence of the pressure on gas diffusion is controversial. Based on the different models, there are notable differences in the diffusion law with the pressure. Based on the unique model, Charrière et al. [23] found that the diffusion coefficient increased with increasing pressure. The results of Cui et al. [4] based on the bidisperse model, show that the diffusion rate of gas in coal is significantly reduced with the increase of pressure due to the high shrinkage of the micropores. In this paper, based on the vitrinite microporous model, the molecular dynamics simulation results of the diffusion coefficient with the pressure changes are shown in Figure 14. The results show that the self-diffusion and transport diffusion first increased and then decreased with increasing pressure, independent of the gas type and concentration. Additionally, the peak pressures of CH 4 , CO 2 , and N 2 were all fixed values of 3 MPa. Before the peak pressure, the increase in pressure aggravated the collision between the gas molecules or between the gas molecules and the vitreous group molecules, leading to an increase in the diffusion coefficient. The decrease in the diffusion coefficient after the peak pressure was caused by the strong collision between the gas molecules and the molecules and the interaction with coal adsorption swelling. Figure 14d,e compare the transport diffusion coefficients calculated with those of Zhao et al. [12], that were based on the Wiser model. Generally, the trend of the change was consistently rising at first and then decreasing, but the peak pressures of CH 4 and CO 2 were lower than in this paper and were not fixed values of 1.5 and 1.25 MPa, respectively, which may be related to the differences in coal rank and pore structure of the model. CH4 and CO2 were lower than in this paper and were not fixed values of 1.5 and 1.25 MPa, respectively, which may be related to the differences in coal rank and pore structure of the model. In contrast to the law of the diffusion coefficient with the pressure, the law of diffusion activation energy decreased with the pressure, as shown in Figure 15a-c, and the peak pressure was is 3 MPa. This indicates that the lower the energy barrier is, the easier it is for the gas to overcome the energy, threshold, and therefore the higher the diffusion coefficient will be. Conversely, when the peak pressure is reached, the energy barrier increases, and the gas hardly overcomes the barrier, and the diffusion coefficient decreases. The relationship between the volume swelling rate of the vitrinite and the pressure is shown in Figure 15d-f. Overall, it was consistent with the law in temperature change and increased with the increase in pressure. Under low-pressure conditions (1-3 MPa), the rate of change was low, and the rate of change was more notable under high-pressure conditions, indicating that the volume swelling effect is more notable under high-pressure conditions, which is consistent with Busch's [29] research results. Volume swelling can lead to the narrowing of the micropore orifice, which is in contrast to the gas under high-pressure conditions, and the potential energy of the pore wall to the gas molecules is enhanced, and thus the diffusion coefficient is reduced. This condition corresponds to a decrease in the diffusion coefficient under the above high-pressure conditions. In contrast to the law of the diffusion coefficient with the pressure, the law of diffusion activation energy decreased with the pressure, as shown in Figure 15a-c, and the peak pressure was is 3 MPa. This indicates that the lower the energy barrier is, the easier it is for the gas to overcome the energy, threshold, and therefore the higher the diffusion coefficient will be. Conversely, when the peak pressure is reached, the energy barrier increases, and the gas hardly overcomes the barrier, and the diffusion coefficient decreases. The relationship between the volume swelling rate of the vitrinite and the pressure is shown in Figure 15d-f. Overall, it was consistent with the law in temperature change and increased with the increase in pressure. Under low-pressure conditions (1-3 MPa), the rate of change was low, and the rate of change was more notable under high-pressure conditions, indicating that the volume swelling effect is more notable under high-pressure conditions, which is consistent with Busch's [29] research results. Volume swelling can lead to the narrowing of the micropore orifice, which is in contrast to the gas under high-pressure conditions, and the potential energy of the pore wall to the gas molecules is enhanced, and thus the diffusion coefficient is reduced. This condition corresponds to a decrease in the diffusion coefficient under the above high-pressure conditions. Based on the molecular structure model of TL-1 vitrinite, the self-diffusion and transport diffusion coefficients of CO2, CH4, and N2 were obtained by the Monte Carlo method and molecular dynamics simulation. The effects of the gas concentration, temperature, and pressure were examined. The main conclusions are as follows: (1) Based on the simulation of the molecular pore structure characteristics of coal vitrinite, it has been shown that the pore volume of the models measured based on the method of using a probe with different radius scales varied greatly, so selecting the appropriate probe to measure the free volume is an important prerequisite for the accurate calculation of the adsorption capacity.
(2) The GCMC results show that the order of the gas adsorption capacity of the three gases was CO2 > CH4 > N2 under the same temperature and pressure conditions. Increasing the temperature is not conducive to gas adsorption. Among these gases, the excess CO2 adsorption capacity reached the saturation pressure value, and the saturation pressure value increased when the temperature was higher than 298 K. The simulated excess adsorption capacity of CH4 in coal was relatively consistent with the experimental results.
(3) The self-diffusion and transport diffusion coefficients in the saturated adsorption state of CO2, CH4, and N2 were ranked as D S (N2) > D S (CO2) > D S (CH4) and D t (N2) > D t (CO2) > D t (CH4), respectively. When the number of adsorption sites was the same, the orders were D S (N2) > D S (CO2) > D S (CH4) and D t (CO2) > D t (N2) > D t (CH4), respectively, which are closely related to the gas molecular dynamics diameter, tensile morphology, and allosteric adsorption. As the gas concentration increased, the three gases self-and transport diffusion coefficients decreased. Based on the five diffusion configurations and the self-diffusion trajectory, the order of the three gases vibration amplitudes was N2 > CO2 > CH4, and the jump time of N2 was the shortest.
(4) The diffusion coefficients of the three gases in coal increased with increasing temperature, and first increasing and then decreasing with increasing pressure. The peak pressure was 3 MPa. The order of the diffusion activation energies of CH4, CO2, and N2 under the saturated adsorption

Conclusions
Based on the molecular structure model of TL-1 vitrinite, the self-diffusion and transport diffusion coefficients of CO 2 , CH 4 , and N 2 were obtained by the Monte Carlo method and molecular dynamics simulation. The effects of the gas concentration, temperature, and pressure were examined. The main conclusions are as follows: (1) Based on the simulation of the molecular pore structure characteristics of coal vitrinite, it has been shown that the pore volume of the models measured based on the method of using a probe with different radius scales varied greatly, so selecting the appropriate probe to measure the free volume is an important prerequisite for the accurate calculation of the adsorption capacity.
(2) The GCMC results show that the order of the gas adsorption capacity of the three gases was CO 2 > CH 4 > N 2 under the same temperature and pressure conditions. Increasing the temperature is not conducive to gas adsorption. Among these gases, the excess CO 2 adsorption capacity reached the saturation pressure value, and the saturation pressure value increased when the temperature was higher than 298 K. The simulated excess adsorption capacity of CH 4 in coal was relatively consistent with the experimental results.
(3) The self-diffusion and transport diffusion coefficients in the saturated adsorption state of CO 2 , CH 4 , and N 2 were ranked as D S (N 2 ) > D S (CO 2 ) > D S (CH 4 ) and D t (N 2 ) > D t (CO 2 ) > D t (CH 4 ), respectively. When the number of adsorption sites was the same, the orders were D S (N 2 ) > D S (CO 2 ) > D S (CH 4 ) and D t (CO 2 ) > D t (N 2 ) > D t (CH 4 ), respectively, which are closely related to the gas molecular dynamics diameter, tensile morphology, and allosteric adsorption. As the gas concentration increased, the three gases self-and transport diffusion coefficients decreased. Based on the five diffusion configurations and the self-diffusion trajectory, the order of the three gases vibration amplitudes was N 2 > CO 2 > CH 4 , and the jump time of N 2 was the shortest.
(4) The diffusion coefficients of the three gases in coal increased with increasing temperature, and first increasing and then decreasing with increasing pressure. The peak pressure was 3 MPa. The order of the diffusion activation energies of CH 4 , CO 2 , and N 2 under the saturated adsorption condition were CH 4 (27.388 kJ/mol) > CO 2 (11.832 kJ/mol) > N 2 (10.396 kJ/mol). The diffusion activation energies first decreased and then increased with the pressure, and the peak values was also 3 MPa. At a high pressure, the effect of volume swelling was more notable.
Author Contributions: J.L. and Y.W. designed the project; J.L. and S.L. performed the simulation and analyzed the data; J.L. wrote the paper and Y.W. corrected it. S.L. modified the formats.