Composition Distribution of the Thermal Soluble Organics from Naomaohu Lignite and Structural Characteristics of the Corresponding Insoluble Portions

With cyclohexane (CH), benzene (BE), and ethyl acetate (EA) as solvents, Naomaohu lignite (NL, a typical oil-rich, low-rank coal) from Hami, Xinjiang, was thermally dissolved (TD) to obtain three types of soluble organics (NLCH, NLBE, and NLEA) and the corresponding insoluble portions (NLCH-R, NLBE-R, and NLEA-R). Ultimate analysis, X-ray photoelectron spectroscopy (XPS), Fourier transform infrared spectroscopy (FTIR), thermogravimetric analysis (TG-DTG), and gas chromatography–mass spectrometry (GC/MS) were used to characterize NL and its soluble and insoluble portions. Results showed that, compared with NL, the C element in NLCH-R, NLBE-R, and NLEA-R increased, while the O element decreased significantly, indicating that thermal dissolution is a carbon enrichment process and an effective deoxidation method. The GC/MS results showed that oxygen-containing organic compounds (OCOCs) are dominant in NLCH, NLBE, and NLEA. NLCH is mainly composed of ketones (11.90%) and esters (19.04%), while NLBE and NLEA are composed of alcohols (12.18% and 2.42%, respectively) and esters (66.09% and 84.08%, respectively), with alkyl and aromatic acid esters as the main components. Among them, EA exhibits significant selective destruction for oxygen-containing functional groups in NL. XPS, FTIR, and TG-DTG results showed that thermal dissolution can not only affect the macromolecular network structure of NL, but also improve its pyrolysis reactivity. In short, thermal dissolution can effectively obtain oxygen-containing organic compounds from NL.


Introduction
As petroleum resources rapidly diminish and energy demands increase, the conversion of coal into liquid fuels, chemicals, and carbon materials is garnering more attention.Among these, oil-rich coal (a kind of coal with a tar yield of 7-12% [1]) holds a significant advantage in improving oil and gas conversion efficiency and reducing economic costs [2,3].The "Implementation Plan for Carbon Peak in the Industrial Field of Xinjiang Uygur Autonomous Region", released on 26 July 2023, emphasizes the acceleration of the development of coal, coal-fired power, and coal chemical industry clusters; the utilization of advanced coal capacity and the enhancement of coal-to-oil and gas processes; oil-rich, low-rank coal fractionation and classification; and clean and efficient utilization.It also stresses the vigorous development of the modern coal chemical industry, the expedited construction of Zhundong and Hami national coal-to-oil and gas strategic bases, and the transformation of the coal-to-oil and gas industry towards specialty fuels, high-end chemicals, etc.Furthermore, the plan aims to develop coal-to-olefins, aromatics, and oxygencontaining compounds as basic chemical raw materials, as well as high-end polyolefins, high-performance polyesters, fibers, and other products.
The Hami region in Xinjiang is abundant in coal resources, with large reserves, diverse varieties, and easy extraction.The predicted coal resources amount to 570.8 billion tons, accounting for 12.5% of the national predicted resources and ranking first in Xinjiang.Particularly, the coal in the Naomaohu mining area of Yiwu County exhibits the "three lows and three highs" characteristics (low ash, low sulfur, low phosphorus, high calorific value, high oil content and high volatile matter).The proportion of oil-rich low-rank coal resources exceeds 90%, with oil content exceeding 10% and reaching up to 16.3% at its highest, making it a rare global resource.The oil-rich, low-rank coal in Naomaohu combines coal, oil, and gas properties, and understanding its molecular structural characteristics is crucial for its transition towards specialty fuels, high-end chemicals, and other directions.
Thermal dissolution (TD) is a technique that allows for the separation of organic matter from coal without the use of hydrogen, the addition of catalysts, operating below 400 • C and employing organic solvents [3][4][5][6].It is also considered one of the most effective methods for efficient coal utilization.Thermal dissolution can break weak covalent bonds and mediumstrength covalent bonds in coal, obtaining the soluble fraction of coal and providing insights into the composition and structure of organic matter within coal at a molecular level [7][8][9][10].Typically, polar solvents are preferred for thermal dissolution over non-polar solvents due to their enhanced effectiveness.Masaki et al. [11] demonstrated that a strong polar solvent (N-methyl-2-pyrrolidone) can disrupt hydrogen bonds in low-rank coal, forming new hydrogen bonds with the coal, thereby increasing the thermal dissolution rate.Wang et al. [12] conducted thermal dissolution of Xianfeng lignite using various organic solvents.Among them, tetralin (a hydrogen donor solvent) plays a role in supplying hydrogen during thermal dissolution, resulting in higher dissolution yields under hightemperature conditions.Benzene, as an inert organic solvent, minimally disrupts oxygen bridge bonds in coal during thermal dissolution and does not react with the soluble fraction of coal.Therefore, the soluble fraction dissolved in benzene can reflect the basic situation of free organic compounds in coal.Zhao et al. [6] used benzene and ethanol as solvents to study the thermal dissolution of Shengli lignite (SL) and Xiaolongtan lignite (XL) at 320 • C. The results showed that, compared to XL, SL had a higher content of aromatic hydrocarbons in the soluble organic matter, with the content of aromatic hydrocarbons dissolved in benzene significantly higher than in ethanol.OCOCs were mainly ketones and phenols, with the phenol content dissolved in benzene lower than in XL, while the phenol content dissolved in ethanol was higher than in XL, indicating a lower number of free phenols in SL.Wang et al. [13] used cyclohexane as a solvent for stepwise thermal dissolution of Xianfeng lignite at 200-320 • C. The results showed that cyclohexane, due to its low boiling point, low viscosity, stable chemical properties, and non-reactivity with the coal sample during thermal dissolution, effectively extracted the organic components of Xianfeng lignite, reflecting the composition and structural characteristics of the coal sample through corresponding characterization.Li et al. [10] conducted thermal dissolution of Zhaotong lignite (ZL) in cyclohexane and methanol at 300 • C. The results indicated that cyclohexane could effectively extract the inherent components of ZL without significant cleavage of covalent bonds.Yang et al. [14] studied the thermal dissolution of SL in ethyl acetate (EA) under various conditions, such as temperature, liquid-solid ratio, and time.The results showed that, as a nucleophilic reagent, EA underwent esterification and ester exchange reactions during thermal dissolution, with EA and ethanol (produced from the reaction of EA with water in SL) significantly participating in the depolymerization of SL by attacking the C aryl , C acyl , and C alkyl structures in SL.Hao et al. [15] conducted thermal dissolution of ZL and NL using ethanol, tetrahydrofuran, and ethyl acetate, and separated the thermal dissolution products into three fractions using column chromatography.The results indicated that ethanol exhibited good selectivity for phenols and esters, tetrahydrofuran released highly condensed polycyclic compounds from low-rank coal, and ethyl acetate showed strong selectivity for C=O in low-rank coal, releasing a significant amount of ester compounds during thermal dissolution.
Low-boiling organic solvents have the potential to reduce solvent loss and environmental emissions as well as enhance product yield and quality, and low-energy reaction conditions have the ability to lower production costs, energy consumption, and carbon emissions.Therefore, this study focuses on the thermal dissolution (TD) of Naomaohu lignite (NL) using cyclohexane (CH), benzene (BE), and ethyl acetate (EA) at 300 • C. The composition and structural characteristics of NL were investigated through characterization techniques, including GC/MS, ultimate analysis, FTIR, XPS, and TGA.This research aims to provide a theoretical basis for the transformation of oil-rich coal into specialty fuels and high-end chemicals.

Composition Distribution of NL Thermal Soluble Organics
Figures 1 and S1-S3 present the total ion chromatograms of NL CH , NL BE and NL EA (retention time 22-46 min).Figure 2 and Tables S1-S12 depict the compositional distribution of three thermally soluble organic matters.As shown in Figures S1-S3 and Tables S1-S12, GC/MS detected 33, 25, and 28 compounds in NL CH , NL BE and NL EA , respectively, categorized into alkanes, olefines, aromatics, oxygen-containing organic compounds (OCOCs) and nitrogen-containing organic compounds (NCOCs).Alkanes were further classified into straight-chain, branched-chain, and cycloalkanes, while aromatics were divided into monocyclic and polycyclic aromatics.Additionally, OCOCs were categorized as alcohols, phenols, ethers, carboxylic acids, esters, and ketones, totaling 6 compound classes.
Figure 2 reveals that OCOCs dominate in NL CH (40.46%),NL BE (80.88%), and NL EA (87.11%), indicating that small-molecule compounds with O in NL can be effectively dissociated and removed during the thermal dissolution process, releasing some small molecules (such as alkanes and aromatics).Alkane compounds follow in abundance, with the order of content being NL CH (26.18%) > NL BE (18.27%) > NL EA (7.57%), suggesting that, as organic inert solvents, CH and BE effectively dissolve a large amount of alkane compounds in coal without reacting with the soluble portion of coal.The carbon numbers are mainly distributed in C 11 -C 25 , with the content of C 13 alkane compounds being the highest in both NL CH and NL BE .
As exhibited in Figure 2c, a detailed comparison of OCOCs in NL CH , NL BE , and NL EA reveals that NL CH is predominantly composed of ketones (11.90%) and esters (19.04%), followed by alcohols, phenols, and carboxylic acids, indicating that CH exhibits good selectivity towards ketones and esters.Specifically, ketones mainly include those containing benzene rings, pyridine ketones, and furan ketones.NL BE and NL EA are primarily composed of alcohols (12.18% and 2.42%) and esters (66.09% and 84.08%), respectively.It is evident that ester compounds in NL thermal dissolution products are more widely distributed and diverse in NL EA , possibly due to the easier cleavage of COO in NL by EA.Tables S4, S8 and S11 indicate that ester compounds in NL CH , NL BE , and NL EA are mainly alkyl esters and aromatic acid esters.Notably, the content of aromatic acid esters in NL CH (9.25%) is equal to that of alkyl esters, while NL BE has a higher content of aromatic acid esters (62.61%), and NL EA is dominated by alkyl esters (81.20%), with higher carbon numbers.In summary, because CH and BE mainly extract the organic components present in NL during the thermal dissolution process, OCOCs in NL primarily exist in the form of free alcohols, ketones, and aromatic acid esters.However, OCOCs in NL EA mainly exist in the form of alkyl esters, as EA undergoes ester exchange reactions with aromatic acid esters in NL during thermal dissolution, generating new alkyl esters, as illustrated in Scheme 1.
Molecules 2024, 29, x FOR PEER REVIEW 5 of 14 In summary, because CH and BE mainly extract the organic components present in NL during the thermal dissolution process, OCOCs in NL primarily exist in the form of free alcohols, ketones, and aromatic acid esters.However, OCOCs in NLEA mainly exist in the form of alkyl esters, as EA undergoes ester exchange reactions with aromatic acid esters in NL during thermal dissolution, generating new alkyl esters, as illustrated in Scheme 1.

Scheme 1. Possible pathways for the ester exchange reaction of NL during thermal dissolution (R' represents aryl group, R represents alkyl group).
The aromatic acid esters in NLBE mainly include 1-Butyl-2-isobutyl phthalate (BABE), di-n-butyl phthalate (DP), and Bis(2-ethylhexyl) phthalate (BEP).DP and BEP are used as plasticizers in the plastics industry to enhance flexibility and toughness of plastic products.BABE, a non-competitive α-glucosidase inhibitor with an IC50 value of 38 μM, exhibits blood glucose-lowering effects and is studied for its potential in diabetes treatment.
The thermal dissolution extraction process involves organic solvents penetrating and diffusing through the pores of coal and dissolving (exchanging) the soluble organic matter throughout the coal's macromolecular network structure.As inert organic solvents (primarily dissolving inherent free small molecules in coal), CH and BE yield significantly different thermal dissolution products due to their distinct physical structures and solubility towards soluble organic matter in coal.CH is an alkane, while BE is an aromatic hydrocarbon, and CH's solubility for soluble organic matter in coal (16.7 J 1/2 /cm 3/2 ) is lower than that of BE (18.5 J 1/2 /cm 3/2 ) [16], resulting in a higher alkane content in NLCH than in NLBE, while the content of aromatic acid esters (the main form of free oxygen-containing organic compounds in NL) is much higher in NLBE than in NLCH.

Impact of Thermal Dissolution on Element Distribution in NL
Based on Table 1, after thermal dissolution treatment, the volatile matter yield of the samples significantly increased, while the moisture and ash content decreased.Additionally, the carbon element content followed the order of NLEA-R > NLCH-R > NLBE-R > NL, indicating that coal thermal dissolution is an effective method for carbon enrichment.Compared to NL, the H/C ratios of NLCH-R, NLBE-R and NLEA-R all decreased, with the most significant reductions observed in NLCH-R and NLBE-R, attributed to the large amount of alkane compounds dissolved in CH and BE in the coal samples.The oxygen content in the coal samples showed the order NLEA-R < NLCH-R < NLBE-R < NL, suggesting that thermal dissolution can effectively remove O from NL. Furthermore, compared with NL, the O/C ratio of the insoluble residues significantly decreased, providing further evidence of the deoxygenation and carbon enrichment effects of thermal dissolution treatment.The aromatic acid esters in NL BE mainly include 1-Butyl-2-isobutyl phthalate (BABE), di-n-butyl phthalate (DP), and Bis(2-ethylhexyl) phthalate (BEP).DP and BEP are used as plasticizers in the plastics industry to enhance flexibility and toughness of plastic products.BABE, a non-competitive α-glucosidase inhibitor with an IC50 value of 38 µM, exhibits blood glucose-lowering effects and is studied for its potential in diabetes treatment.
The thermal dissolution extraction process involves organic solvents penetrating and diffusing through the pores of coal and dissolving (exchanging) the soluble organic matter throughout the coal's macromolecular network structure.As inert organic solvents (primarily dissolving inherent free small molecules in coal), CH and BE yield significantly different thermal dissolution products due to their distinct physical structures and solubility towards soluble organic matter in coal.CH is an alkane, while BE is an aromatic hydrocarbon, and CH's solubility for soluble organic matter in coal (16.7 J 1/2 /cm 3/2 ) is lower than that of BE (18.5 J 1/2 /cm 3/2 ) [16], resulting in a higher alkane content in NL CH than in NL BE , while the content of aromatic acid esters (the main form of free oxygen-containing organic compounds in NL) is much higher in NL BE than in NL CH .

Impact of Thermal Dissolution on Element Distribution in NL
Based on Table 1, after thermal dissolution treatment, the volatile matter yield of the samples significantly increased, while the moisture and ash content decreased.Additionally, the carbon element content followed the order of NL EA-R > NL CH-R > NL BE-R > NL, indicating that coal thermal dissolution is an effective method for carbon enrichment.Compared to NL, the H/C ratios of NL CH-R , NL BE-R and NL EA-R all decreased, with the most significant reductions observed in NL CH-R and NL BE-R , attributed to the large amount of alkane compounds dissolved in CH and BE in the coal samples.The oxygen content in the coal samples showed the order NL EA-R < NL CH-R < NL BE-R < NL, suggesting that thermal dissolution can effectively remove O from NL. Furthermore, compared with NL, the O/C ratio of the insoluble residues significantly decreased, providing further evidence of the deoxygenation and carbon enrichment effects of thermal dissolution treatment.

Changes in Functional Groups of Insoluble Portions
The FTIR spectra of NL and its thermally insoluble portions, as shown in Figure 3, reveal that the peak types of NL, NL CH-R , NL BE-R , and NL EA-R are generally similar, but with significant differences in peak intensities, indicating that thermal dissolution treatment affects the macro-molecular network structure of NL.The peak at 3340 cm −1 corresponds to the stretching vibration peak of intermolecular or intramolecular −OH bonds.Compared to NL, the peak intensities of this peak in the three thermally insoluble portions are weaker, attributed to the disruption of hydrogen bonds and other non-covalent bonds by the organic solvents used [17].The absorption peaks near 2926 cm −1 and 2855 cm −1 are related to the stretching vibrations of −CH 2 − groups, while the peak at 1378 cm −1 is attributed to −CH 3 groups.The peak intensities of these peaks in the three thermally insoluble portions are significantly lower than in NL, indicating that a large portion of the aliphatic organic compounds in NL were extracted by CH, BE, and EA during the thermal dissolution process.Moreover, the peak intensities of NL CH-R and NL BE-R are relatively similar, while NL EA-R exhibits stronger intensities, suggesting that CH and BE have a stronger solubilizing effect on aliphatic organic compounds in NL, leading to higher levels of these compounds in NL EA-R compared to NL CH-R and NL BE-R (consistent with the GC/MS results).The peak at 1710 cm −1 corresponds to the stretching vibration of conjugated C=O bonds, mainly from ketones, esters, and carboxylic acids.The strong absorption peak at 1583 cm −1 in all four samples indicates the presence of a significant amount of olefinic and aromatic ring structures in the coal samples.In NL EA-R , the peak at 1271 cm −1 for >CH−OH is stronger, while the intensity near 1110 cm −1 for >CH−O−CH< is weaker, suggesting that >CH−O−CH< bonds in NL were disrupted during the EA thermal dissolution process.

Changes in Functional Groups of Insoluble Portions
The FTIR spectra of NL and its thermally insoluble portions, as shown in Figure 3, reveal that the peak types of NL, NLCH-R, NLBE-R, and NLEA-R are generally similar, but with significant differences in peak intensities, indicating that thermal dissolution treatment affects the macro-molecular network structure of NL.The peak at 3340 cm −1 corresponds to the stretching vibration peak of intermolecular or intramolecular −OH bonds.Compared to NL, the peak intensities of this peak in the three thermally insoluble portions are weaker, attributed to the disruption of hydrogen bonds and other non-covalent bonds by the organic solvents used [17].The absorption peaks near 2926 cm −1 and 2855 cm −1 are related to the stretching vibrations of −CH2− groups, while the peak at 1378 cm −1 is attributed to −CH3 groups.The peak intensities of these peaks in the three thermally insoluble portions are significantly lower than in NL, indicating that a large portion of the aliphatic organic compounds in NL were extracted by CH, BE, and EA during the thermal dissolution process.Moreover, the peak intensities of NLCH-R and NLBE-R are relatively similar, while NLEA-R exhibits stronger intensities, suggesting that CH and BE have a stronger solubilizing effect on aliphatic organic compounds in NL, leading to higher levels of these compounds in NLEA-R compared to NLCH-R and NLBE-R (consistent with the GC/MS results).The peak at 1710 cm −1 corresponds to the stretching vibration of conjugated C=O bonds, mainly from ketones, esters, and carboxylic acids.The strong absorption peak at 1583 cm −1 in all four samples indicates the presence of a significant amount of olefinic and aromatic ring structures in the coal samples.In NLEA-R, the peak at 1271 cm −1 for >CH−OH is stronger, while the intensity near 1110 cm −1 for >CH−O−CH< is weaker, suggesting that >CH−O−CH< bonds in NL were disrupted during the EA thermal dissolution process.According to the Beer-Lambert Law, the intensity (area) of absorption peaks in the infrared spectrum is directly proportional to the concentration of functional groups causing the absorption if KBr pellets contain coal samples of the same concentration.In complex multi-component systems, the infrared bands often overlap, making it difficult to clearly observe the weak covalent bonds in coal.Therefore, curve-fitting analysis is used to determine individual peaks in the coal infrared spectrum [18].This method can identify hidden peaks with fixed frequencies and half-widths but allows for changes in relative intensity.This is an important principle for obtaining accurate and reproducible results.According to the Beer-Lambert Law, the intensity (area) of absorption peaks in the infrared spectrum is directly proportional to the concentration of functional groups causing the absorption if KBr pellets contain coal samples of the same concentration.In complex multi-component systems, the infrared bands often overlap, making it difficult to clearly observe the weak covalent bonds in coal.Therefore, curve-fitting analysis is used to determine individual peaks in the coal infrared spectrum [18].This method can identify hidden peaks with fixed frequencies and half-widths but allows for changes in relative intensity.This is an important principle for obtaining accurate and reproducible results.
As shown in Figure 4 and Table 2, the 3000-3600 cm −1 region represents hydroxyl absorption bands.Compared to NL (39.75%), the self-condensation −OH content (peak area percentage) in NL CH-R (29.29%),NL BE-R (34.79%), and NL EA-R (35.18%) all decreased, indicating that CH, BE, and EA disrupted the weak hydrogen bonds in NL during thermal dissolution.The asymmetric −CH 2 and symmetric −CH 2 absorption peaks at 2930-2900 cm −1 and 2870-2810 cm −1 are relatively high in all samples, with contents above 20% (except for NL BE-R ).
As shown in Figure 4 and Table 2, the 3000-3600 cm −1 region represents hydroxyl absorption bands.Compared to NL (39.75%), the self-condensation −OH content (peak area percentage) in NLCH-R (29.29%),NLBE-R (34.79%), and NLEA-R (35.18%) all decreased, indicating that CH, BE, and EA disrupted the weak hydrogen bonds in NL during thermal dissolution.The asymmetric −CH2 and symmetric −CH2 absorption peaks at 2930-2900 cm −1 and 2870-2810 cm −1 are relatively high in all samples, with contents above 20% (except for NLBE-R). ) 3000 3100 3200 3300 3400 3500 3600 3000 3100 3200 3300 3400 3500 3600 3000 3100 3200 3300 3400 3500 3600 3000 3100 3200 3300 3400 3500 3600    The ratio of the absorption peak areas of −CH 2 to −CH 3 can reflect the relative content of organic structures, denoted as the aliphatic structure parameter S = A(CH 2 )/A(CH 3 ) [22,23].The values for NL (3.71) and NL CH-R (1.53) are higher than those for NL BE-R (0.53) and NL EA-R (0.69), indicating that aliphatic hydrocarbons in NL and NL CH-R mainly exist in long-chain forms with fewer side-chain aliphatic hydrocarbons.Additionally, trisubstituted benzene rings dominate in all four samples, with the relative content ranking as NL EA-R > NL BE-R > NL CH-R > NL > 45%, indicating that hydrogen atoms on the benzene rings of the coal samples are substituted by other atoms or groups.This suggests that aromatic compounds in the coal samples may undergo a series of substitution reactions during the thermal dissolution process, leading to structural changes in the coal samples.

Wavenumber (cm
Figure 5 and Table 3 present the XPS spectra of NL and its insoluble portions, as well as the relative content of various functional groups.In Figure 5, the peak intensity of C 1s increases in the order of NL < NL CH-R < NL BE-R < NL EA-R , while the peak intensity of O 1s is exactly the opposite, aligning with the elemental analysis findings.According to references [4,15,17,24,25], the binding energies of C 1s at 284.0 eV, 285.6 eV, 287.3 eV and 288.7 eV correspond to C−C, C−H, C=O, and COO, respectively.Table 3 reveals that the dominant form of C 1s on the surface of the four samples is C−C, while O 1s is mainly in the form of C=O and C−O, primarily existing as phenols, alcohols, and esters.The results of C 1s indicate that, compared to the original coal, the C=O content of NL decreased from 10.54% to 5.18%, 6.82%, and 5.67% after thermal dissolution with CH, BE, and EA, respectively, while the COO content increased from 4.01% to 4.59%, 4.20%, and 4.61%, respectively.This may be due to the destruction of oxygen bridges (such as C=O) in the coal sample during thermal dissolution in the presence of EA, leading to the formation of esters and carboxylic acid compounds [15,26].Additionally, during the thermal dissolution processes with CH and BE at 300 • C, the capsule structure [27] in NL was disrupted, releasing some small OCOCs, with some C=O-containing compounds dissolving in CH and BE, while some COO-containing compounds were retained in the coal sample structure.

Thermal Decomposition Characteristics of Insoluble Portions
Figure 6 shows the TG-DTG curves of NL and its insoluble portions.As depicted in Figure 6a, the four samples exhibit similar trends of mass loss, with mass losses in the order of NL (47.85%) > NLCH-R (36.43%) > NLBE-R (35.66%) > NLEA-R (35.24%).This is attributed to the partial extraction of soluble components from NL during the thermal   [4,15,17,24,25].From Figure 5 and Table 3, it can be observed that, compared to NL, the C=O content increased in NL CH-R , NL BE-R , and NL EA-R , while the C−O and COO contents relatively decreased.The lowest C−O content was observed in NL CH-R (10.60%), indicating higher content in NL CH (as indicated by GC/MS, primarily existing in the form of aromatic alcohols, alkyl chain alcohols, and phenols).The COO content in NL EA-R decreased from 6.08% to 0.54%, which may be due to the reaction of carboxylic acid compounds in NL with EA during the thermal dissolution process, forming ethyl acetate and acetic acid [15], with possible reaction pathways, as outlined in Scheme 2. In conclusion, the thermal dissolution treatment in the presence of CH, BE, and EA can break down and convert the oxygen-containing functional groups in NL into various forms of oxygen-containing compounds.with possible reaction pathways, as outlined in Scheme 2. In conclusion, the thermal dissolution treatment in the presence of CH, BE, and EA can break down and convert the oxygen-containing functional groups in NL into various forms of oxygen-containing compounds.
Scheme 2. Acidolysis of EA with higher carboxylic acids (R denotes an alkyl group with more than two carbon atoms).

Thermal Decomposition Characteristics of Insoluble Portions
Figure 6 shows the TG-DTG curves of NL and its insoluble portions.As depicted in Figure 6a, the four samples exhibit similar trends of mass loss, with mass losses in the order of NL (47.85%) > NL CH-R (36.43%) > NL BE-R (35.66%) > NL EA-R (35.24%).This is attributed to the partial extraction of soluble components from NL during the thermal dissolution process by organic solvents, leading to a reduction in the mass loss of the insoluble portions.
Molecules 2024, 29, x FOR PEER REVIEW 10 of 14 dissolution process by organic solvents, leading to a reduction in the mass loss of the insoluble portions.Figure 7 and Table 4 present the DTG fitting and attribution of fitting parameters for the four samples, revealing fundamental information on the covalent bonds in NL and its insoluble portions.Combining with Figure 6b, the thermal decomposition process of the coal samples can be divided into four stages.Primarily, all four samples (especially NL) exhibit a significant peak in the mass loss rate at around 80 °C, mainly due to prolonged sample-air contact time (resulting in excessive moisture absorption).Below 400 °C (peaks 1-3), the main processes involve the volatilization of moisture and some low-boiling-point   4 present the DTG fitting and attribution of fitting parameters for the four samples, revealing fundamental information on the covalent bonds in NL and its insoluble portions.Combining with Figure 6b, the thermal decomposition process of the coal samples can be divided into four stages.Primarily, all four samples (especially NL) exhibit a significant peak in the mass loss rate at around 80 • C, mainly due to prolonged sample-air contact time (resulting in excessive moisture absorption).Below 400 • C (peaks 1-3), the main processes involve the volatilization of moisture and some low-boilingpoint substances, as well as the release of bound water and the cleavage of carboxyl groups, accompanied by the breakage of weak bonds, such as >C al −N< and >C al −S−.Comparatively, NL EA-R shows a distinct peak in the mass loss rate between 130 • C and 400 • C, attributed to the thorough swelling and loosening of the sample due to thermal dissolution treatment, promoting the escape of small molecular substances from the coal and enhancing the diffusion of some high-boiling-point organic matter within the sample.Figure 7 and Table 4 present the DTG fitting and attribution of fitting parameters for the four samples, revealing fundamental information on the covalent bonds in NL and its insoluble portions.Combining with Figure 6b, the thermal decomposition process of the coal samples can be divided into four stages.Primarily, all four samples (especially NL) exhibit a significant peak in the mass loss rate at around 80 °C, mainly due to prolonged sample-air contact time (resulting in excessive moisture absorption).Below 400 °C (peaks 1-3), the main processes involve the volatilization of moisture and some low-boiling-point substances, as well as the release of bound water and the cleavage of carboxyl groups, accompanied by the breakage of weak bonds, such as >Cal−N< and >Cal−S−.Comparatively, NLEA-R shows a distinct peak in the mass loss rate between 130 °C and 400 °C, attributed to the thorough swelling and loosening of the sample due to thermal dissolution treatment, promoting the escape of small molecular substances from the coal and enhancing the diffusion of some high-boiling-point organic matter within the sample.4. The temperature range of 380-610 • C represents the main process of organic matter thermal decomposition in the samples, involving the cleavage of various covalent bonds leading to significant mass loss.Between 380 • C and 500 • C, the predominant bond cleavages include >C al −C al <, >C al −H, and >C ar −N, while 440-610 • C mainly involves the dissociation of strong covalent bonds, such as >C ar −C al <, >C ar −O-, and >C ar −S−.At temperatures exceeding 610 • C (peaks 6 and 7), the process primarily involves the condensation of aromatic rings and the decomposition of mineral matter, releasing some small condensed molecules, like H 2 O and CO 2 .Additionally, due to the high reaction temperature, H 2 may also be generated during the condensation reactions of aromatic rings.In conclusion, thermal dissolution can effectively alter the thermal decomposition characteristics of coal samples (with EA exerting the most significant influence).Furthermore, selective cleavage of covalent bonds with different binding energies can be achieved by adjusting the thermal decomposition temperature, providing insights for the targeted extraction of organic matter from the samples.

Raw Materials
NL was collected from the Naomaohu coal mine in Hami, Xinjiang, China.After crushing, it was sifted through a 200-mesh sieve (<74 µm) and dried naturally for 24 h.CH, BE, and EA were purchased from Tianjin Xin Platte Chemical Co., Ltd.(Tianjin, China).All the solvents are analytical reagents and were distilled in a rotary evaporator before use.

Thermal Dissolution
NL (2 g) and CH were added to a 50 mL beaker at a solid-liquid ratio of 1:10 and thoroughly stirred before transferring to a high-pressure reaction vessel equipped with a magnetic stirrer.After purging the vessel with N 2 to displace air, the high-pressure reaction vessel was heated to 300 • C at a rate of 20 • C/min under an initial pressure of 1 MPa and maintained for 2 h.Upon completion of the reaction, the vessel was cooled to room temperature, and the mixture was filtered to separate the filtrate and residue.The residue was washed repeatedly with CH until the filtrate became colorless and transparent, yielding CH-soluble (NL CH ) and -insoluble (NL CH-R ) fractions.The thermal dissolution process of NL is illustrated in Figure 8.A similar procedure was followed to obtain soluble and insoluble fractions of benzene (BE) and ethyl acetate (EA), resulting in samples designated as NL BE , NL BE-R , NL EA, and NL EA-R , respectively.

Analysis Methods
The thermal dissolution products NLCH, NLBE, and NLEA were analyzed using GC/MS (Thermo Fisher Scientific UltiMate 3000 UHPLC-Q Exactive, Waltham, Massachusetts, USA).The instrument was equipped with an HP-5MS chromatographic column, and highpurity helium gas (1 mL/min) was used as the carrier gas.The inlet temperature was set at 250 °C with a temperature program of 60 °C for 3 min, followed by a ramp rate of 5 °C/min up to 300 °C for 10 min.The split ratio was 10:1, and the mass scan range was set from 30 amu to 500 amu.
Ultimate analysis of C, H, N, and S content in NL and the thermally insoluble fraction was conducted using an ultimate analyzer (Elementar UNICUBE, Langenselbold, Germany), with the oxygen content calculated by the difference method.XPS (Thermo Scientific K-Alpha, Waltham, Massachusetts, USA) was employed to determine the surface elemental composition of NL and the thermally insoluble fraction.After preparing pressed pellets of the samples and placing them on the sample holder, the samples were introduced into the instrument chamber.The pressure in the chamber was maintained below 2.0 × 10 −7 mbar, with a spot size of 400 μm, operating at 12 kV and 6 mA filament current.The full-scan pass energy was set at 150 eV, with a step size of 1 eV, while the narrow-scan pass energy was 50 eV, with a step size of 0.1 eV.

Analysis Methods
The thermal dissolution products NL CH , NL BE , and NL EA were analyzed using GC/MS (Thermo Fisher Scientific UltiMate 3000 UHPLC-Q Exactive, Waltham, Massachusetts, USA).The instrument was equipped with an HP-5MS chromatographic column, and high-purity helium gas (1 mL/min) was used as the carrier gas.The inlet temperature was set at 250 • C with a temperature program of 60 • C for 3 min, followed by a ramp rate of 5 • C/min up to 300 • C for 10 min.The split ratio was 10:1, and the mass scan range was set from 30 amu to 500 amu.
Ultimate analysis of C, H, N, and S content in NL and the thermally insoluble fraction was conducted using an ultimate analyzer (Elementar UNICUBE, Langenselbold, Germany), with the oxygen content calculated by the difference method.XPS (Thermo Scientific K-Alpha, Waltham, Massachusetts, USA) was employed to determine the surface elemental composition of NL and the thermally insoluble fraction.After preparing

Figure 3 .
Figure 3. FTIR spectra of NL and its insoluble portions (NLCH−R, NLBE−R, and NLEA−R are thermal insoluble substances of CH, BE, and EA, respectively).

Figure 3 .
Figure 3. FTIR spectra of NL and its insoluble portions (NL CH−R , NL BE−R , and NL EA−R are thermal insoluble substances of CH, BE, and EA, respectively).

Figure 4 .
Figure 4. FTIR fitting results of NL and its insoluble portions (NLCH−R, NLBE−R, and NLEA−R are thermal insoluble substances of CH, BE, and EA, respectively).The solid line represents the raw data, and the dashed line represents the fitted data.

Figure 4 .
Figure 4. FTIR fitting results of NL and its insoluble portions (NL CH−R , NL BE−R , and NL EA−R are thermal insoluble substances of CH, BE, and EA, respectively).The solid line represents the raw data, and the dashed line represents the fitted data.

Figure 5 .
Figure 5. X-ray photoelectron energy spectra and fitting curves of NL and its insoluble portions (NLCH−R, NLBE−R, and NLEA−R are thermal insoluble substances of CH, BE, and EA, respectively).The solid line represents the raw data, and the dashed line represents the fitted data.

Figure 5 .
Figure 5. X-ray photoelectron energy spectra and fitting curves of NL and its insoluble portions (NL CH−R , NL BE−R , and NL EA−R are thermal insoluble substances of CH, BE, and EA, respectively).The solid line represents the raw data, and the dashed line represents the fitted data.

s
2024, 29, x FOR PEER REVIEW 9 of 14

Figure 6 .
Figure 6.TG (a) and DTG (b) profiles of NL and its insoluble portions (NLCH−R, NLBE−R, and NLEA−R are thermal insoluble substances of CH, BE, and EA, respectively).

Figure 6 .
Figure 6.TG and DTG (b) profiles of NL and its insoluble portions (NL CH−R , NL BE−R , and NL EA−R are thermal insoluble substances of CH, BE, and EA, respectively).

Figure 7
Figure7and Table4present the DTG fitting and attribution of fitting parameters for the four samples, revealing fundamental information on the covalent bonds in NL and its

Figure 6 .
Figure 6.TG (a) and DTG (b) profiles of NL and its insoluble portions (NLCH−R, NLBE−R, and NLEA−R are thermal insoluble substances of CH, BE, and EA, respectively).

Figure 7 .Figure 7 .
Figure 7. DTG fitting curve of NL and its insoluble portions (NLCH−R, NLBE−R, and NLEA−R are thermal insoluble substances of CH, BE, and EA, respectively).The solid line represents the raw data, and Figure 7. DTG fitting curve of NL and its insoluble portions (NL CH−R , NL BE−R , and NL EA−R are thermal insoluble substances of CH, BE, and EA, respectively).The solid line represents the raw data, and the dashed line represents the fitted data.The corresponding explanations of the numbers are listed in Table4.

Figure 8 .
Figure 8. Thermal dissolution process of NL (NLCH, NLBE, and NLEA are the thermal soluble substances of CH, BE, and EA, respectively; NLCH−R, NLBE−R and NLEA−R are thermal-insoluble substances of CH, BE, and EA, respectively).

Figure 8 .
Figure 8. Thermal dissolution process of NL (NL CH , NL BE , and NL EA are the thermal soluble substances of CH, BE, and EA, respectively; NL CH−R , NL BE−R and NL EA−R are thermal-insoluble substances of CH, BE, and EA, respectively).

Table 1 .
Proximate and ultimate analyses of NL and its insoluble portions.Possible pathways for the ester exchange reaction of NL during thermal dissolution (R' represents aryl group, R represents alkyl group).
ad: air dry basis; daf: dry and ash-free basis; diff : by difference.Scheme 1.

Table 1 .
Proximate and ultimate analyses of NL and its insoluble portions.
ad : air dry basis; daf : dry and ash-free basis; diff : by difference.
Molecules 2024,29,x FOR PEER REVIEW 9 of 14with possible reaction pathways, as outlined in Scheme 2. In conclusion, the thermal dissolution treatment in the presence of CH, BE, and EA can break down and convert the oxygen-containing functional groups in NL into various forms of oxygen-containing compounds.
Scheme 2. Acidolysis of EA with higher carboxylic acids (R denotes an alkyl group with more than two carbon atoms).

Table 3 .
Distribution of C and O forms in NL and its insoluble portions.