Study on the CO Formation Mechanism during Coal Ambient Temperature Oxidation

: The CO formation rules of coal were analyzed by a self-developed testing device under ambient temperature. The changes of functional groups caused by oxidation were obtained using Fourier-transform infrared spectroscopy (FTIR). The experimental results showed that CO was generated during the ambient temperature oxidation. The highest concentration level of CO could be 389 ppm. The methylene and aldehyde groups on the side chains were involved in the reaction. For the quantum mechanical approach, we employed the density functional theory with the 6–31 G (d, p) basis set. Density functional theory–based computations interpreted the possible reaction sites on a coal molecule by electronic static potential analysis. The rationality of the predicted reactions was also evaluated by transition state analysis and energy analysis. This research theoretically proved that coal could be oxidized to carbon monoxide under ambient temperatures and gave the possible reaction paths.


Introduction
CO is often used as an essential index gas for early prediction of coal spontaneous combustion for its high sensitivity and easy detection [1][2][3][4]. According to China's Coal Mine Safety Regulation, underground CO concentration should not be higher than 24 ppm. Otherwise, measures must be taken to prevent coal spontaneous combustion. However, many field practices show that, in the mine area, the CO concentration often exceeds the standard without apparent signs of spontaneous combustion [5]. Sometimes, CO can be detected in the mining face or air return way without apparent signs of spontaneous combustion, such as tarry odor, coal wall sweating, and smoking. Therefore, it is of considerable significance for mine safety production to study the mechanism of the spontaneous combustion of coal, especially the mechanism of oxidation to form CO at low temperatures.
Many researchers have studied the mechanism of coal oxidation at low temperatures [6][7][8]. By using the method of density functional theory (DFT), Wang and Deng [9] analyzed the change rules of different functional groups in the molecular structure of coal and simulated the physical adsorption, chemical adsorption, and chemical reaction processes of different functional groups of coal with oxygen. Wang Deming et al. [10] described the oxidation mechanism of coal spontaneous combustion by the DFT method. They proposed a series of cyclic chain reactions. Deng Jun [11] found eight index gases of coal spontaneous combustion and corresponding characteristic temperatures by thermal analysis and large-scale experiments on coal spontaneous combustion.
The oxidation of coal at low temperatures has been well studied. However, there is very little research on the oxidation of coal, especially the formation mechanism of CO under ambient temperature. Therefore, a study of the CO formation mechanism during the coal ambient temperature oxidation is necessary. This kind of study can give us a deeper understanding of the initial stages of coal spontaneous combustion. Our previous research has already done some work on the calculation of the CO production rate of coal under ambient temperatures and at critical temperatures [12]. To follow up on the previous research, the demand for this research is to theoretically prove that coal can be oxidized to carbon monoxide under ambient temperature and give the possible reaction paths. In this paper, we investigate the oxidation of Australia White (AW) coal to CO by self-made equipment for coal oxidation under ambient temperature. The changes of functional groups in coal during oxidation under ambient temperature were analyzed by infrared spectroscopy. Density functional theory (DFT)-based computations interpreted the reactive sites and electronic static potential (ESP) on the coal molecule. The rationality of the predicted reaction was also evaluated by the transition state analysis and energy analysis. This study may be helpful to understand the formation mechanism of CO under ambient temperature.

Coal Samples
The coal samples were from an AW coal mine that had shown self-ignition tendencies. Fresh coal samples weighing up to 50 kg were cut from the working face. The coal samples were then moved to the laboratory. Before the test, the oxide layer on the coal sample surface was stripped. The remaining unoxidized coal samples were crushed and sieved with different particle diameters (2.00~3.35 mm). The proximate analysis and elemental analysis of the coal samples are shown in Table 1. Note: R0max is the maximum reflectance of the vitrinite in coal, Mad is the mass fraction of moisture in the sample on dry air basis, Aad is the mass fraction of ash in the sample on dry air basis, and Vdaf is the mass fraction of volatile matter in the sample on the dry ash-free basis. w(O*) =1-w(C)-w(H)-w(N)-w(S).

Ambient Temperature Oxidation
The experiment of coal oxidation under ambient temperature was carried out in a hermetically sealed reactor (see Figure 1). The sensors, including carbon monoxide, oxygen, and temperature, were set up and connected to the computer. The testing device consists of a cylindrical chamber with a volume of 1.718 × 10 −2 m 3 (the height is 350 mm, and the bottom radius is 125 mm). The measurement range of the CO sensor is 0-1000 ppm, and the resolution is 1 ppm. For the O 2 sensor, the measurement range is 0-25 VOL%, and the resolution is 0.1 VOL%. All the data of the sensors were recorded by the computer automatically. In particular, there was no heating part in the reactor. Also, no air can flow in or out of though the device during the reaction. The coal samples were crushed and screened into samples of 2.00~3.35 mm particle size, and then the coal samples (3.5 kg) were sealed into the reaction chamber. The experiment stopped when the concentration of CO stopped changing. The room temperature was around 15 • C during the experiment.

FTIR Study
The experiments were also carried out using the Fourier-transform infrared spectroscopy (FTIR) (MIRacle 10 SHIMADZU; JPN). The FTIR study can help to obtain a better understanding of coal oxidization to carbon monoxide under ambient temperature. The samples before and after the ambient temperature oxidation were tested respectively for the comparison. Five fresh coal samples (10 g each) were selected for the FTIR test. Another five coal samples (10 g each) were ground into fine particles and placed with the air under ambient temperature. The FTIR test for the second group was conducted after the completion of the ambient temperature oxidation experiment. Functional groups of coal samples before and after the oxidation were obtained, respectively. The test parameters were set as follows: the scan range was 400-4000 cm −1 , the resolution was 2 cm −1 , and the scan times were 64.

DFT Calculation
By analyzing the FTIR test results and combining previous studies on the coal oxidation reaction process and mechanism, the possible reaction of AW coal samples during the ambient temperature oxidation process was predicted. The geometries of reactants, products, intermediates, and transition state molecules in the reaction process were optimized using the Gaussian 09 program package with DFT at the level of B3 LYP/6-31 G (d, p). The potential distribution on the surface of the reactant was calculated. Combined with the FTIR experiment results, the reaction locations that may be involved in the oxidation reaction were analyzed. By calculating the energies of reactants, products, intermediates, and transition state molecules, the possibility of the predicted reaction paths was judged from the energy perspective.

Ambient Temperature Oxidation Results
The ambient temperature oxidation experiment consists of two parts: the first oxidation and second oxidation. The primary oxidation lasted for 118 h, and then open the reactor, and refresh the coal samples by the air. The secondary oxidation lasted for 38 h. The experiment results can be seen in Figure 2 and Figure 3.

FTIR Study
The experiments were also carried out using the Fourier-transform infrared spectroscopy (FTIR) (MIRacle 10 SHIMADZU; JPN). The FTIR study can help to obtain a better understanding of coal oxidization to carbon monoxide under ambient temperature. The samples before and after the ambient temperature oxidation were tested respectively for the comparison. Five fresh coal samples (10 g each) were selected for the FTIR test. Another five coal samples (10 g each) were ground into fine particles and placed with the air under ambient temperature. The FTIR test for the second group was conducted after the completion of the ambient temperature oxidation experiment. Functional groups of coal samples before and after the oxidation were obtained, respectively. The test parameters were set as follows: the scan range was 400-4000 cm −1 , the resolution was 2 cm −1 , and the scan times were 64.

DFT Calculation
By analyzing the FTIR test results and combining previous studies on the coal oxidation reaction process and mechanism, the possible reaction of AW coal samples during the ambient temperature oxidation process was predicted. The geometries of reactants, products, intermediates, and transition state molecules in the reaction process were optimized using the Gaussian 09 program package with DFT at the level of B3 LYP/6-31 G (d, p). The potential distribution on the surface of the reactant was calculated. Combined with the FTIR experiment results, the reaction locations that may be involved in the oxidation reaction were analyzed. By calculating the energies of reactants, products, intermediates, and transition state molecules, the possibility of the predicted reaction paths was judged from the energy perspective.

Ambient Temperature Oxidation Results
The ambient temperature oxidation experiment consists of two parts: the first oxidation and second oxidation. The primary oxidation lasted for 118 h, and then open the reactor, and refresh the coal samples by the air. The secondary oxidation lasted for 38 h. The experiment results can be seen in   During the primary oxidation, the CO appeared after 1.5 h. The maximum of CO concentration level can be 389 ppm after 79 h. CO production leads to O2 consumption. After 66 h, the O2 concentration level reached a minimum of 15.3%. During the secondary oxidation, the CO appeared after 4 h. The maximum CO concentration level can be 80 ppm after 26 h. After 29 h, the reaction stopped and the O2 concentration level reached a minimum of 17.9%.
We can conclude that the coal samples can react with O2 and produce CO under ambient temperature. The CO during the primary oxidation appears earlier and have more production than the secondary oxidation.

FTIR Analysis
In the initial stage of coal oxidation, only a small number of active groups participate in the reactions. During the burning stage, a large number of secondary active groups are generated and a chain reaction occurs rapidly, and the reaction strength increases rapidly [13]. The FTIR analysis was further carried out to get more information on the difference in chemical composition during the oxidation processes. By comparing the FTIR test results of coal samples, the functional groups that might be involved in the reaction were inferred.
As shown in Figure 4, several peaks of oxidized coal have changed in the spectrum compared with that of fresh coal. Peaks decreased after oxidation at 3300 cm -1 , 2860-2920 cm -1 , and 1610 cm -1 , which could be assigned to the "hydroxyl group (-OH)," "stretching vibration of methyl (-CH3) asymmetric stretching vibration of methylene (-CH2)," and "carbonyl group(C = O)." The peak of carboxyl appears around 2321-2400 cm -1 and was strengthened after the oxidation.  During the primary oxidation, the CO appeared after 1.5 h. The maximum of CO concentration level can be 389 ppm after 79 h. CO production leads to O2 consumption. After 66 h, the O2 concentration level reached a minimum of 15.3%. During the secondary oxidation, the CO appeared after 4 h. The maximum CO concentration level can be 80 ppm after 26 h. After 29 h, the reaction stopped and the O2 concentration level reached a minimum of 17.9%.
We can conclude that the coal samples can react with O2 and produce CO under ambient temperature. The CO during the primary oxidation appears earlier and have more production than the secondary oxidation.

FTIR Analysis
In the initial stage of coal oxidation, only a small number of active groups participate in the reactions. During the burning stage, a large number of secondary active groups are generated and a chain reaction occurs rapidly, and the reaction strength increases rapidly [13]. The FTIR analysis was further carried out to get more information on the difference in chemical composition during the oxidation processes. By comparing the FTIR test results of coal samples, the functional groups that might be involved in the reaction were inferred.
As shown in Figure 4, several peaks of oxidized coal have changed in the spectrum compared with that of fresh coal. Peaks decreased after oxidation at 3300 cm -1 , 2860-2920 cm -1 , and 1610 cm -1 , which could be assigned to the "hydroxyl group (-OH)," "stretching vibration of methyl (-CH3) asymmetric stretching vibration of methylene (-CH2)," and "carbonyl group(C = O)." The peak of carboxyl appears around 2321-2400 cm -1 and was strengthened after the oxidation. During the primary oxidation, the CO appeared after 1.5 h. The maximum of CO concentration level can be 389 ppm after 79 h. CO production leads to O 2 consumption. After 66 h, the O 2 concentration level reached a minimum of 15.3%. During the secondary oxidation, the CO appeared after 4 h. The maximum CO concentration level can be 80 ppm after 26 h. After 29 h, the reaction stopped and the O 2 concentration level reached a minimum of 17.9%.
We can conclude that the coal samples can react with O 2 and produce CO under ambient temperature. The CO during the primary oxidation appears earlier and have more production than the secondary oxidation.

FTIR Analysis
In the initial stage of coal oxidation, only a small number of active groups participate in the reactions. During the burning stage, a large number of secondary active groups are generated and a chain reaction occurs rapidly, and the reaction strength increases rapidly [13]. The FTIR analysis was further carried out to get more information on the difference in chemical composition during the oxidation processes. By comparing the FTIR test results of coal samples, the functional groups that might be involved in the reaction were inferred.
As shown in Figure 4, several peaks of oxidized coal have changed in the spectrum compared with that of fresh coal. Peaks decreased after oxidation at 3300 cm −1 , 2860-2920 cm −1 , and 1610 cm −1 , which could be assigned to the "hydroxyl group (-OH)," "stretching vibration of methyl (-CH 3 ) asymmetric stretching vibration of methylene (-CH 2 )," and "carbonyl group (C = O)." The peak of carboxyl appears around 2321-2400 cm −1 and was strengthened after the oxidation.  These results indicate that methyl or methylene and carbonyl group in XM coal may react with oxygen under ambient temperature. Evaporation of water from coal also results in the reduction of hydroxyl functional groups. There was no more evidence that the hydroxyl group was involved in oxidation reactions or not. As an important product, the carboxyl group appears in the oxidation reaction of coal under ambient temperature. The carboxyl group plays a significant role during the reaction.
According to the Transition-state theory, a chemical reaction system can be divided into the following three regions: a) the reactant region, b) the transition state region, and c) the product region [14]. A lot of work has been done on the mechanism of oxidation of the active groups in coal molecules. The active groups of coal molecules involved in the oxidation process are divided into bridged bonds and side chains. According to the previous study on the molecular structure of coal [15], it is inferred that the methyl-ether, the alkylidene in bridged bonds, and the methoxyl aldehyde in the side chain are involved in the oxidation progress. Based on the results of previous studies and combine with the FTIR results of this study, the possible reactions to carbon monoxide under ambient temperature are inferred as follows: Reaction path a: Step1. (1) In this step, the chemisorptions take place between O2 and reactant 1 (R1) and produce the intermediate (M1). Step2. ( The intermediate (M1) creates H2 O and reactant 2 (R2) with the carboxy group. Step3.
(3) These results indicate that methyl or methylene and carbonyl group in XM coal may react with oxygen under ambient temperature. Evaporation of water from coal also results in the reduction of hydroxyl functional groups. There was no more evidence that the hydroxyl group was involved in oxidation reactions or not. As an important product, the carboxyl group appears in the oxidation reaction of coal under ambient temperature. The carboxyl group plays a significant role during the reaction.
According to the Transition-state theory, a chemical reaction system can be divided into the following three regions: (a) the reactant region, (b) the transition state region, and (c) the product region [14]. A lot of work has been done on the mechanism of oxidation of the active groups in coal molecules. The active groups of coal molecules involved in the oxidation process are divided into bridged bonds and side chains. According to the previous study on the molecular structure of coal [15], it is inferred that the methyl-ether, the alkylidene in bridged bonds, and the methoxyl aldehyde in the side chain are involved in the oxidation progress. Based on the results of previous studies and combine with the FTIR results of this study, the possible reactions to carbon monoxide under ambient temperature are inferred as follows: Reaction path a: Step 1.
Energies 2020, 13, x FOR PEER REVIEW 5 of 12 These results indicate that methyl or methylene and carbonyl group in XM coal may react with oxygen under ambient temperature. Evaporation of water from coal also results in the reduction of hydroxyl functional groups. There was no more evidence that the hydroxyl group was involved in oxidation reactions or not. As an important product, the carboxyl group appears in the oxidation reaction of coal under ambient temperature. The carboxyl group plays a significant role during the reaction.
According to the Transition-state theory, a chemical reaction system can be divided into the following three regions: a) the reactant region, b) the transition state region, and c) the product region [14]. A lot of work has been done on the mechanism of oxidation of the active groups in coal molecules. The active groups of coal molecules involved in the oxidation process are divided into bridged bonds and side chains. According to the previous study on the molecular structure of coal [15], it is inferred that the methyl-ether, the alkylidene in bridged bonds, and the methoxyl aldehyde in the side chain are involved in the oxidation progress. Based on the results of previous studies and combine with the FTIR results of this study, the possible reactions to carbon monoxide under ambient temperature are inferred as follows: Reaction path a: Step1. (1) In this step, the chemisorptions take place between O2 and reactant 1 (R1) and produce the intermediate (M1). Step2.
In this step, the chemisorptions take place between O 2 and reactant 1 (R1) and produce the intermediate (M 1 ).
Energies 2020, 13, x FOR PEER REVIEW 5 of 12 These results indicate that methyl or methylene and carbonyl group in XM coal may react with oxygen under ambient temperature. Evaporation of water from coal also results in the reduction of hydroxyl functional groups. There was no more evidence that the hydroxyl group was involved in oxidation reactions or not. As an important product, the carboxyl group appears in the oxidation reaction of coal under ambient temperature. The carboxyl group plays a significant role during the reaction.
According to the Transition-state theory, a chemical reaction system can be divided into the following three regions: a) the reactant region, b) the transition state region, and c) the product region [14]. A lot of work has been done on the mechanism of oxidation of the active groups in coal molecules. The active groups of coal molecules involved in the oxidation process are divided into bridged bonds and side chains. According to the previous study on the molecular structure of coal [15], it is inferred that the methyl-ether, the alkylidene in bridged bonds, and the methoxyl aldehyde in the side chain are involved in the oxidation progress. Based on the results of previous studies and combine with the FTIR results of this study, the possible reactions to carbon monoxide under ambient temperature are inferred as follows: Reaction path a: Step1. (1) In this step, the chemisorptions take place between O2 and reactant 1 (R1) and produce the intermediate (M1). Step2.
Energies 2020, 13, x FOR PEER REVIEW 5 of 12 These results indicate that methyl or methylene and carbonyl group in XM coal may react with oxygen under ambient temperature. Evaporation of water from coal also results in the reduction of hydroxyl functional groups. There was no more evidence that the hydroxyl group was involved in oxidation reactions or not. As an important product, the carboxyl group appears in the oxidation reaction of coal under ambient temperature. The carboxyl group plays a significant role during the reaction.
According to the Transition-state theory, a chemical reaction system can be divided into the following three regions: a) the reactant region, b) the transition state region, and c) the product region [14]. A lot of work has been done on the mechanism of oxidation of the active groups in coal molecules. The active groups of coal molecules involved in the oxidation process are divided into bridged bonds and side chains. According to the previous study on the molecular structure of coal [15], it is inferred that the methyl-ether, the alkylidene in bridged bonds, and the methoxyl aldehyde in the side chain are involved in the oxidation progress. Based on the results of previous studies and combine with the FTIR results of this study, the possible reactions to carbon monoxide under ambient temperature are inferred as follows: Reaction path a: Step1. (1) In this step, the chemisorptions take place between O2 and reactant 1 (R1) and produce the intermediate (M1). Step2. ( The intermediate (M1) creates H2 O and reactant 2 (R2) with the carboxy group. Step3. (3) The reactant 2 (R 2 ) decomposes into product 1 (P 1 ) and CO.

Calculation Model
Coal is a complex organic macromolecule with a heterogeneous chemical composition [16]. The Wiser model of coal is a reasonable and accurate chemical structure model in the field of coal chemistry. It reflects the new progress in the study of the chemical structure of coal. It can explain many chemical reactions of coal [17]. When investigating the reaction involving the local structure of molecules, some relatively small aromatic molecules, such as toluene and dibenzyl benzene, can often be used to represent the structure of coal [18,19]. Based on the Wise coal molecular model, and combined with the infrared analysis results, the feature part was abstracted to form the local microstructure model (see Figure 5). Then, quantum chemical calculations of the local microstructure model can reveal the relationship between the structure of the coal and the specific reaction from a microscopic perspective.

ESP Analysis
The reactivity of coal depends on its structure, especially on its microscopic electronic structure. The active site of the reactants can be inferred from the ESP and electron density on the surface of each molecule [20]. Based on this experimental fact, the study of ESP of molecules has been successfully applied to the researches of the relationships between structures and properties [21][22][23].
The REDOX reaction between coal molecules and oxygen molecules leads to the spontaneous combustion of the coal. The oxygen molecules, as oxidants, gain electrons in the reaction, and some certain functional groups of coal molecules, as reductants, lose electrons in the reaction. Therefore, it is meaningful to study the ESP of coal molecules. First, the simplified molecular structure of coal was optimized. The optimization path is shown in Figure 6. The reactant (R 1 ) with aldehyde group decomposes into the product (P 1 ) and CO.

Calculation Model
Coal is a complex organic macromolecule with a heterogeneous chemical composition [16]. The Wiser model of coal is a reasonable and accurate chemical structure model in the field of coal chemistry. It reflects the new progress in the study of the chemical structure of coal. It can explain many chemical reactions of coal [17]. When investigating the reaction involving the local structure of molecules, some relatively small aromatic molecules, such as toluene and dibenzyl benzene, can often be used to represent the structure of coal [18,19]. Based on the Wise coal molecular model, and combined with the infrared analysis results, the feature part was abstracted to form the local microstructure model (see Figure 5). Then, quantum chemical calculations of the local microstructure model can reveal the relationship between the structure of the coal and the specific reaction from a microscopic perspective.

Calculation Model
Coal is a complex organic macromolecule with a heterogeneous chemical composition [16]. The Wiser model of coal is a reasonable and accurate chemical structure model in the field of coal chemistry. It reflects the new progress in the study of the chemical structure of coal. It can explain many chemical reactions of coal [17]. When investigating the reaction involving the local structure of molecules, some relatively small aromatic molecules, such as toluene and dibenzyl benzene, can often be used to represent the structure of coal [18,19]. Based on the Wise coal molecular model, and combined with the infrared analysis results, the feature part was abstracted to form the local microstructure model (see Figure 5). Then, quantum chemical calculations of the local microstructure model can reveal the relationship between the structure of the coal and the specific reaction from a microscopic perspective.

ESP Analysis
The reactivity of coal depends on its structure, especially on its microscopic electronic structure. The active site of the reactants can be inferred from the ESP and electron density on the surface of each molecule [20]. Based on this experimental fact, the study of ESP of molecules has been successfully applied to the researches of the relationships between structures and properties [21][22][23].
The REDOX reaction between coal molecules and oxygen molecules leads to the spontaneous combustion of the coal. The oxygen molecules, as oxidants, gain electrons in the reaction, and some certain functional groups of coal molecules, as reductants, lose electrons in the reaction. Therefore, it is meaningful to study the ESP of coal molecules. First, the simplified molecular structure of coal was optimized. The optimization path is shown in Figure 6.

ESP Analysis
The reactivity of coal depends on its structure, especially on its microscopic electronic structure. The active site of the reactants can be inferred from the ESP and electron density on the surface of each molecule [20]. Based on this experimental fact, the study of ESP of molecules has been successfully applied to the researches of the relationships between structures and properties [21][22][23].
The REDOX reaction between coal molecules and oxygen molecules leads to the spontaneous combustion of the coal. The oxygen molecules, as oxidants, gain electrons in the reaction, and some certain functional groups of coal molecules, as reductants, lose electrons in the reaction. Therefore, it is meaningful to study the ESP of coal molecules. First, the simplified molecular structure of coal was optimized. The optimization path is shown in Figure 6.
It can be seen from Figure 6 that, after 25 steps of optimization, the total energy of the structure gradually tends to be stable and reaches the minimum value. The energy curve shows that the final structure is the optimal structure. The ESP distribution can be obtained by frequency and energy calculations based on the optimal structure. ESP mapped on the isosurface (0.001 a.u.) of the simplified coal molecular structure was plotted in Figure 7. Moreover, the negative potential distribution on the surface of the coal molecules can be seen in Table 2. It can be seen from Figure 6 that, after 25 steps of optimization, the total energy of the structure gradually tends to be stable and reaches the minimum value. The energy curve shows that the final structure is the optimal structure. The ESP distribution can be obtained by frequency and energy calculations based on the optimal structure. ESP mapped on the isosurface (0.001 a.u.) of the simplified coal molecular structure was plotted in Figure 7. Moreover, the negative potential distribution on the surface of the coal molecules can be seen in Table 2.  It can be seen from Table 2 that the negative potential sites mainly appear near the oxygen atoms of oxygen-containing functional groups and some benzene rings. There are large negative potential values near O81 and O88 (-0.06726 a.u. and -0.06628 a.u., respectively). It indicates that these two functional groups are more likely to react with oxygen. The ESP analysis results are consistent with the predicted react path a and b.  It can be seen from Figure 6 that, after 25 steps of optimization, the total energy of the structure gradually tends to be stable and reaches the minimum value. The energy curve shows that the final structure is the optimal structure. The ESP distribution can be obtained by frequency and energy calculations based on the optimal structure. ESP mapped on the isosurface (0.001 a.u.) of the simplified coal molecular structure was plotted in Figure 7. Moreover, the negative potential distribution on the surface of the coal molecules can be seen in Table 2.  It can be seen from Table 2 that the negative potential sites mainly appear near the oxygen atoms of oxygen-containing functional groups and some benzene rings. There are large negative potential values near O81 and O88 (-0.06726 a.u. and -0.06628 a.u., respectively). It indicates that these two functional groups are more likely to react with oxygen. The ESP analysis results are consistent with the predicted react path a and b.  Table 2. Negative potential distribution on the surface of coal molecules.

No.
Site Potential (a.u.) It can be seen from Table 2 that the negative potential sites mainly appear near the oxygen atoms of oxygen-containing functional groups and some benzene rings. There are large negative potential values near O81 and O88 (−0.06726 a.u. and −0.06628 a.u., respectively). It indicates that these two functional groups are more likely to react with oxygen. The ESP analysis results are consistent with the predicted react path a and b.
The same conclusion can be drawn from Figure 7. We can see that there are positive and negative potential distributions on the surface of coal molecules. The side chains with hydroxyl groups attached to α-carbon atoms have a negative electrostatic potential. The maximum is −0.06726 a.u. Due to the electron conjugation, an unshared pair of electrons in the -OH group tends to become the electron donor. Moreover, the oxygen atoms in the hydroxyl group are more electronegative than the carbon atoms, leading to the inductive effect of electron absorption. These two effects weaken the C-H bond and make it easy to form stable free radicals.
Similarly, there is a negative electrostatic potential near the aldehyde group at the edge of the benzene ring. The maximum is −0.06628 a.u. The carbonyl group in the aldehyde group produces π-π conjugation with the C-C bond attached to it. Moreover, the aromatic ring produces π-π conjugation with the same C-C bond. Therefore, the positive charge of the carbonyl group disperses into the aromatic ring due to the delocalization.
According to the calculation results, the extremum of the electrostatic potential occurs near the oxygen and hydrogen atoms of the oxygen-containing functional groups (particularly near the side chains with hydroxyl groups attached to α-carbon atoms and the aldehyde group at the edge of the benzene ring). This indicates that, compared with a benzene ring, methyl group, and other structures, oxygen-containing functional groups are more likely to react with oxygen molecules. Low metamorphic coal contains more oxygen-containing functional groups, which brings more extreme points to the surface electrostatic potential of coal molecules, thus providing more reaction sites. This conclusion is consistent with the previous studies.

Energy Analysis
The transition state structure refers to the highest point of energy in the reaction path. The reaction rate can be calculated by the height of the energy barrier. The difficulty of the reaction can be evaluated by the transition state of the reaction. The default calculation settings for temperature and pressure are 298.15 k and 1 atm, respectively. The enthalpy and Gibbs free energy of the molecules can be obtained by optimizing the molecular structure to the lowest energy. By optimizing the transition state to the energy saddle point (QST3), the optimized transition state, the corresponding enthalpy, and Gibbs free energy can be obtained.
Moreover, the reaction rate constant was calculated based on the transition state theory [24].
where k is the reaction rate constant (s −1 ); σ is reaction path degeneracy (the value is 4); k B is The Boltzmann constant; T is the temperature of the coal sample (K); h is the Planck constant; R is Gas constant; P is the reaction pressure (kPa); ∆Ea is the activation energies (kcal/mol) where ∆Ea is the activation energies (kcal/mol); H TS is the enthalpy of the transition state (kcal/mol); H Reactant is the enthalpy of the reactant (kcal/mol). The standard thermodynamic parameters, including enthalpy (H), Gibbs free energy (G), enthalpy changes (∆H), Gibbs free energy changes (∆G), activation energies (Ea), and reaction rate constant (k), are summarized in Table 3. For reaction path a, the ∆H and ∆G of step 1 and step 2 are −61.16 kcal/mol, −56.30 kcal/mol, −78.03 kcal/mol, and −80.38 kcal/mol, respectively. This means that these two reactions are exothermic and spontaneous.
Step 3 is an endothermic reaction because the ∆H is 24.02 kcal/mol. The ∆G of step 3 is 19.08 kcal/mol which means that the reactions cannot take place spontaneously. The Ea and k of step 3 show that the reaction will take place with the energy generated by step 1 and step 2. For reaction path b, we can conclude that it is an exothermic and spontaneous reaction by ∆H and ∆G. By comparing the reaction rates, it can be seen that reaction path b is more difficult than reaction path a. From Figures 8 and 9, we can see that the Ea of the three steps in the reaction path a is 16.24 kcal/mol, 25.68 kcal/mol and 34.55 kcal/mol. This means the chemisorptions in step 1 will take place easily and widely. Some of the intermediate (M 1 ) can decompose into reactant 2 (R 2 ) with carboxy group (we have detected the increase in the carboxyl group by FTIR), and less reactant 2 (R 2 ) continue to decompose into product 1 (P 1 ) and CO; The Ea of reaction path b are 84.76 kcal/mol. They are higher than the Ea in reaction path a. But it will not be too hard to make the reaction take place with the energy generated by reaction path a.
Energies 2020, 13, x FOR PEER REVIEW 10 of 12 For reaction path a, the △H and △G of step 1 and step 2 are -61.16 kcal/mol, -56.30 kcal/mol, -78.03 kcal/mol, and -80.38 kcal/mol, respectively. This means that these two reactions are exothermic and spontaneous.
Step 3 is an endothermic reaction because the △H is 24.02 kcal/mol. The △G of step 3 is 19.08 kcal/mol which means that the reactions cannot take place spontaneously. The Ea and k of step 3 show that the reaction will take place with the energy generated by step 1 and step 2. For reaction path b, we can conclude that it is an exothermic and spontaneous reaction by △H and △G. By comparing the reaction rates, it can be seen that reaction path b is more difficult than reaction path a.
From Figure 8 and Figure 9, we can see that the Ea of the three steps in the reaction path a is 16.24 kcal/mol, 25.68 kcal/mol and 34.55 kcal/mol. This means the chemisorptions in step 1 will take place easily and widely. Some of the intermediate (M1) can decompose into reactant 2 (R2) with carboxy group (we have detected the increase in the carboxyl group by FTIR), and less reactant 2 (R2) continue to decompose into product 1 (P1) and CO; The Ea of reaction path b are 84.76 kcal/mol. They are higher than the Ea in reaction path a. But it will not be too hard to make the reaction take place with the energy generated by reaction path a.

Conclusions
Coal can react with oxygen and generate CO under ambient temperature. The CO formation rules of AW coal under ambient temperature were analyzed by self -developed coal oxidation device. During the experiment, the maximum concentration level of CO can be 389 ppm after 79 h. The FTIR comparison test shows that methylene and aldehyde groups involved in the reaction. The carboxy is For reaction path a, the △H and △G of step 1 and step 2 are -61.16 kcal/mol, -56.30 kcal/mol, -78.03 kcal/mol, and -80.38 kcal/mol, respectively. This means that these two reactions are exothermic and spontaneous.
Step 3 is an endothermic reaction because the △H is 24.02 kcal/mol. The △G of step 3 is 19.08 kcal/mol which means that the reactions cannot take place spontaneously. The Ea and k of step 3 show that the reaction will take place with the energy generated by step 1 and step 2. For reaction path b, we can conclude that it is an exothermic and spontaneous reaction by △H and △G. By comparing the reaction rates, it can be seen that reaction path b is more difficult than reaction path a.
From Figure 8 and Figure 9, we can see that the Ea of the three steps in the reaction path a is 16.24 kcal/mol, 25.68 kcal/mol and 34.55 kcal/mol. This means the chemisorptions in step 1 will take place easily and widely. Some of the intermediate (M1) can decompose into reactant 2 (R2) with carboxy group (we have detected the increase in the carboxyl group by FTIR), and less reactant 2 (R2) continue to decompose into product 1 (P1) and CO; The Ea of reaction path b are 84.76 kcal/mol. They are higher than the Ea in reaction path a. But it will not be too hard to make the reaction take place with the energy generated by reaction path a.

Conclusions
Coal can react with oxygen and generate CO under ambient temperature. The CO formation rules of AW coal under ambient temperature were analyzed by self -developed coal oxidation device. During the experiment, the maximum concentration level of CO can be 389 ppm after 79 h. The FTIR comparison test shows that methylene and aldehyde groups involved in the reaction. The carboxy is a key functional group for coal to produce CO, and it was detected as an intermediate in the oxidation

Conclusions
Coal can react with oxygen and generate CO under ambient temperature. The CO formation rules of AW coal under ambient temperature were analyzed by self -developed coal oxidation device. During the experiment, the maximum concentration level of CO can be 389 ppm after 79 h. The FTIR comparison test shows that methylene and aldehyde groups involved in the reaction. The carboxy is a key functional group for coal to produce CO, and it was detected as an intermediate in the oxidation reaction. There would be some interesting findings if we can associate the generation of carboxyl groups with the spontaneous combustion tendency of coal. This is well worth pursuing further. Based on the FTIR results and combined with previous research, the possible reaction paths were concluded for AW coal oxidation to produce CO under ambient temperature. The ESP calculation showed that the side chains with hydroxyl groups attached to α-carbon atoms and the aldehyde group at the edge of the benzene ring are the possible reaction sites. The standard thermodynamic parameters of the molecules along the predicted reaction paths were calculated, and the results show that the experiment can be carried out spontaneously under ambient temperature. This study may be helpful in understanding the rule of coal oxidation under ambient temperature. This research theoretically proved that coal could be oxidized to carbon monoxide under ambient temperature. Furthermore, under certain conditions, CO can accumulate to a very high concentration. The results of the experiment and analysis show that the presence of a certain concentration of CO under the coal mine may not be caused by a serious fire. The warning value of carbon monoxide needs to be evaluated, respectively. This is of considerable significance to guide coal mine safety management.