Detailed Group-Type Characterization of Plastic-Waste Pyrolysis Oils: By Comprehensive Two-Dimensional Gas Chromatography Including Linear, Branched, and Di-Oleﬁns

: Plastic-waste pyrolysis oils contain large amounts of linear, branched, and di-oleﬁnic compounds. This makes it not obvious to determine the detailed group-type composition in particular to the presence of substantial amounts of N-, S-, and O-containing heteroatomic compounds. The thorough evaluation of different column combinations for two-dimensional gas chromatography (GC × GC), i.e., non-polar × polar and polar × non-polar , revealed that the second combination had the best performance, as indicated by the bi-dimensional resolution of the selected key compounds. By coupling the GC × GC to multiple detectors, such as the ﬂame ionization detector (FID), a sulfur chemiluminescence detector (SCD), a nitrogen chemiluminescence detector (NCD), and a mass spectrometer (MS), the identiﬁcation and quantiﬁcation were possible of hydrocarbon, oxygen-, sulfur-, and nitrogen-containing compounds in both naphtha (C5–C11) and diesel fractions (C7–C23) originating from plastic-waste pyrolysis oils. Group-type quantiﬁcation showed that large amounts of α -oleﬁns (36.39 wt%, 35.08 wt%), iso-oleﬁns (8.77 wt%, 9.06 wt%), and dioleﬁns (4.21 wt%, 4.20 wt%) were present. Furthermore, oxygen-containing compounds (alcohols, ketones, and ethers) could be distinguished from abundant hydrocarbon matrix, by employing Stabilwax as the ﬁrst column and Rxi-5ms as the second column. Ppm levels of sulﬁdes, thiophenes, and pyridines could also be quantiﬁed by the use of selective SCD and NCD detectors.


Introduction
Disposal and upgrading of plastic waste is one of the main challenges of the 21st century [1,2]. Chemical recycling as an alternative to landfill through (catalytic) pyrolysis [3][4][5] is a very promising conversion technology to produce liquid energy carriers, waxes, and chemicals while minimizing negative impacts on the environment [6][7][8].
Catalytic and non-catalytic pyrolysis are processes in which organic material (i.e., plastic waste) is heated at high temperatures (300-900 • C), in the absence of oxygen, to break down its macromolecular constituents into smaller molecules with or without the aid of a catalyst. The product consists of non-condensables, a liquid fraction (aka the plasticwaste pyrolysis oil), and char [9][10][11][12]. For example, it has been shown that pyrolysis oil obtained by pyrolysis of post-consumer High-Density Polyethylene (HDPE) bears fuel properties similar to crude oil [13,14]. In general, these pyrolysis oils are complex mixtures of n-paraffins, iso-paraffins, olefins, diolefins, iso-olefins, naphthenes, and aromatics with a wide carbon number distribution (C9-C25). It also contains a relatively low abundance of heteroatom compounds such as oxygen-, nitrogen-, and sulfur-containing compounds, as well as halogens and metals [14][15][16][17][18]. However, the compositions thereof vary widely depending on the plastic source used and the operating conditions applied (e.g., temperature, pressure, and residence time) [19][20][21].

Sample Preparation
Two fractions of plastic-waste pyrolysis oils were produced by using a pilot plant unit for chemical recycling in the Laboratory of Chemical Technology at Ghent University. The first one is a light fraction with a boiling range similar to naphtha and is called "plastic waste naphtha". The second one is a heavier product with a wider boiling-point range and is referred to as "plastic waste diesel". Four different samples of each fraction were prepared for the separate analysis on GC × GC-FID, GC × GC-SCD, GC × GC-NCD, and GC × GC-MS. An internal standard was added to the sample for the GC × GC-FID/SCD/NCD analysis. The selected internal standard should be easily detectable on the selective detector, inert, and well-separated from the hydrocarbon matrix. Therefore, 3-chlorothiophene was chosen as an internal standard for the FID and SCD analyses, and 2-chloropyridine was chosen for the NCD analysis in line with previous work [16]. Furthermore, the amount of internal standard added to the sample should result in a similar peak height as the main compounds that need to be quantified. This results in 2.12 wt% and 2.14 wt% of 3-chlorothiophene added to the plastic-waste naphtha and diesel fraction for the FID analysis. For the SCD analysis, 90 ppm of 3-chlorothiophene was added, while for the NCD analysis, 8 ppm 2-chloropyridine was added.

Elemental Analyzer Setup
The elemental composition of the diesel and naphtha fraction was determined by using a Thermo Scientific™ FLASH 2000 Series Elemental Analyzer (EA) (Interscience, Louvain-la-Neuve, Belgium) equipped with both a Thermal Conductivity Detector (TCD) and a Flame Photometric Detector (FPD). The EA instrument consists of two different configurations, namely combustion and pyrolysis modes. The first one is used to determine the content of carbon, hydrogen, nitrogen, and sulfur (CHNS), while the latter is used to calculate elemental oxygen content (O). The information of the elemental content (CHNSO) could be obtained by the same setup in two separate operation modes. Because the samples are quite volatile, a liquid injection technique was applied to avoid interference of oxygen and nitrogen from the air.

GC × GC-FID/SCD/NCD/MS Setup
Two fractions were analyzed by using Thermo Scientific TRACE GC × GC setups (Interscience, Louvain-la-Neuve, Belgium). The setups are composed of a dual-state cryogenic modulator (liquid CO 2 ) and a detector FID/SCD/NCD or MS. The different column sets and modulator were placed in the same oven. Various GC × GC column configurations, i.e., non-polar × polar (NP), polar × non-polar (PN), and mid-polar × nonpolar (MPN), were tested in this study. An overview of the tested column combinations is given in Table 1.
The first GC × GC device was equipped with an FID and an SCD detector (Thermo Scientific, Interscience Belgium), and the second device was equipped with an NCD and an MS detector (Thermo Scientific, Interscience Belgium). For the FID, flow rates of H 2 , air, and N 2 (make-up gas) were set at 35, 350, and 35 mL min −1 , respectively. The FID temperature was set at 300 • C, and the data-acquisition rate was 100 Hz. For the SCD, flow rates of H 2 and air were set at 45 and 5 mL min −1 , respectively, while the burner temperature was set at 800 • C. For the NCD, the H 2 and O 2 flow rates were set at 5 and 10 mL min −1 , respectively, while the burner temperature was set at 920 • C. The data-acquisition rate of the SCD and NCD was set at 100 Hz. For the MS, the data-acquisition rate was set at 30 spectra s −1 , with the scanning range set from 30 to 400 amu. The GC-MS interface (transfer line) temperature was set at 280 • C, and the ion source temperature was set at 300 • C. The MS detector used electron ionization (70 eV). Helium was used as a carrier gas at a constant flow rate. The flow rates of carrier gas depended on the column combination and detector used. The modulation period was optimized to obtain a maximal resolution in the first dimension without causing a wrap-around. Therefore, different modulation regimes were applied depending on the column combination used, as well as the type of fraction. The GC system was operated in programmed temperature conditions and the applied programs for each column combination and sample type. A summary of GC × GC parameters applied in this work is given in Table 2.

Data Acquisition and Quantification
Thermo Scientific's Chrom-Card data system was used for data acquisition and processing of the FID, SCD, and NCD, while Thermo Scientific's XCalibur software was applied for data acquired with MS. The raw data of GC × GC-FID/SCD/NCD were ex- ported to a cdf file, subsequently processed by GC Image (Zoex Corporation, Houston, TX, USA) for quantification. The tentative identification of the peaks is based on the combined information from the literature, Kovats retention indices, and the group-tile chromatogram obtained by the orthogonal separation of GC × GC. The identification of hydrocarbon compounds was performed by the use of GC × GC-MS, which could not be applied for sulfur-and nitrogen-containing compounds because of the interference of the hydrocarbon matrix. To assist in the identification of the sulfur-containing compounds, a standard mixture consisting of carbon disulfide, thiophene, dimethyl disulfide, 2-methyl thiophene, and 3-propyl thiophene was analyzed by the GC × GC-SCD setup to determine retention times thereof in both dimensions of GC × GC-SCD chromatogram (see Figure S1, Supporting Materials). Furthermore, a standard mixture consisting of pyridine, 2-methylpyridine, 2-ethylpyridine, and 5-ethyl-2-methylpyridine ( Figure S2, Supporting Materials) was analyzed to support the identification of the nitrogen-containing compounds in the respective GC × GC-NCD chromatogram.
The peak names and peak volumes were exported as a csv file which was subsequently processed by using an in-house macro file. Thanks to the stability of the FID response, the quantification of the present compound groups was conducted from the FID chromatogram [55], based on the internal standard calibration method and the effective carbon number developed by Djokic et al. [56] and Dijkmans et al. [23,57]. In particular, a weight fraction of each compound is calculated based on the known amount of internal standard, using the following equation: where wt% i is the weight fraction of compound I; wt% IS is the weight fraction of internal standard; f i is the relative response factor of compound i, which is calculated based on effective carbon number approach [58]; f IS is the relative response factor of the internal standard (3-chlorothiophene), which is determined as 1.94 based on the calibration; A i is the peak area of compound i; and A IS is peak area of the internal standard. Furthermore, the quantification of sulfur-and nitrogen-containing compounds was based on SCD and NCD chromatograms applying the same methodology for FID. However, the response of the SCD and the NCD detectors is linear and equimolar to the amount of sulfur and nitrogen present in the sample. As a result, the following equation was applied, which was reported elsewhere [16].
where Ni is the number of sulfur or nitrogen atoms in compound i, and M IS is the molecular mass of the sulfur or nitrogen standard (3-chlorothiophene and 2-chloropyridine, respectively). The mass fractions of all the quantified compounds in both samples are in a range of 95 and 105 wt%. Subsequently, all of these mass fractions are normalized to 100 wt%.

Column Combination Selection and Evaluation
The performance and suitability of different column combinations were evaluated by comparing the chromatographic resolution obtained for each analysis by calculating the bi-dimensional resolution of the internal standard, each detected compound, and the intraclass bi-dimensional resolution [27,54]. The bi-dimensional resolution was introduced by Giddings [59] and used by Dutriez [60], according to the following formula: While R 1 is the first-dimension resolution, R 2 is the second-dimension resolution; ∆t r1 is the difference of retention time between compounds A and B in the first dimension; ∆t r2 is the difference of retention time between compounds A and B in the second dimension; ω A1 is the peak width of compound A in the first dimension; ω B1 is the peak width of compound B in the first dimension; ω A2 is the peak width of compound A in the second dimension; and ω B2 is the peak width of compound B in the second dimension. The average bidimensional resolutions in the FID chromatograms of the most abundant chemical groups have been compared with different column combinations. The separation between the peaks is acceptable if the value of resolution R 2D ≥ 1, and the separation between the peaks is good if the value of resolution R 2D ≥ 1.5.

Plastic-Waste Diesel Pyrolysis Oil
Initially, the diesel fraction was analyzed by using a normal-phase GC × GC column combination. In particular, a non-polar RTX-1 PONA column (50 m × 0.25 mm × 0.25 µm) in the first dimension was configured with a mid-polar BPX-50 column (2 m × 0.18 mm × 0.18 µm) in the second dimension. This combination gives the most abundant chemical groups present in the diesel sample as illustrated on the GC × GC-FID color plot in Figure 1. It should be noted that the ramp temperature plays an important role in the separation of iso-paraffins and iso-olefins. Namely, the ramp temperature at a rate of 2 • C/min instead of 3 • C/min considerably improved the separation between these adjacent groups. Furthermore, the chosen internal standard 3-chlorothiophene is nicely separated from the chemical components of the sample, with a calculated bi-dimensional resolution of 11.93.
The performance and suitability of different column combinations were evaluated by comparing the chromatographic resolution obtained for each analysis by calculating the bi-dimensional resolution of the internal standard, each detected compound, and the intraclass bi-dimensional resolution [27,54]. The bi-dimensional resolution was introduced by Giddings [59] and used by Dutriez [60], according to the following formula: While R1 is the first-dimension resolution, R2 is the second-dimension resolution; Δtr1 is the difference of retention time between compounds A and B in the first dimension; Δtr2 is the difference of retention time between compounds A and B in the second dimension; ωA1 is the peak width of compound A in the first dimension; ωB1 is the peak width of compound B in the first dimension; ωA2 is the peak width of compound A in the second dimension; and ωB2 is the peak width of compound B in the second dimension. The average bi-dimensional resolutions in the FID chromatograms of the most abundant chemical groups have been compared with different column combinations. The separation between the peaks is acceptable if the value of resolution R2D ≥ 1, and the separation between the peaks is good if the value of resolution R2D ≥ 1.5.

Plastic-Waste Diesel Pyrolysis Oil
Initially, the diesel fraction was analyzed by using a normal-phase GC × GC column combination. In particular, a non-polar RTX-1 PONA column (50 m × 0.25 mm × 0.25 μm) in the first dimension was configured with a mid-polar BPX-50 column (2 m × 0.18 mm × 0.18 μm) in the second dimension. This combination gives the most abundant chemical groups present in the diesel sample as illustrated on the GC × GC-FID color plot in Figure  1. It should be noted that the ramp temperature plays an important role in the separation of iso-paraffins and iso-olefins. Namely, the ramp temperature at a rate of 2 °C/min instead of 3 °C/min considerably improved the separation between these adjacent groups. Furthermore, the chosen internal standard 3-chlorothiophene is nicely separated from the chemical components of the sample, with a calculated bi-dimensional resolution of 11.93.  The tentative identification of resulting peaks from GC × GC-FID was achieved by a reproducible analysis using GC × GC-MS with the identical column combination and an optimized carrier gas flow [61,62]. The result showed that the abundant chemical groups found in the diesel fraction are n-paraffins, α-olefins, iso-paraffins, iso-olefins, diolefins, naphthenes, and aromatics, which eluted and separated in the group-tile structure on the chromatogram. However, the preliminary result of the elemental analysis revealed that a small amount of oxygen-containing compounds is present in the sample, which could not be detected by using the normal-phase GC × GC combination. This could be explained by the co-elution of oxygenates and the abundant hydrocarbon matrix [27]. The biggest hurdle is the identification of diolefins due to the presence of unsaturated naphthenes, which have the same molecular mass and very similar mass fragmentation patterns on GC × GC-MS chromatogram with diolefins. Nevertheless, these diolefins are not completely separated from iso-olefins, and also neighboring peaks of the high amounts of both paraffins and olefins, which could cause an error in terms of the quantification. As aforementioned, the content of olefins and diolefins results in negative effects on quality control, as well as the technical process of fuels and feedstocks due to the gum-formation property through chemical secondary reaction (e.g., Diels-Alder reaction) [32,63]. Moreover, the heat value of plastic pyrolysis oils is reduced due to the presence of this unsaturated groups [64].
Taking advantages of reversed-phase column configuration, a Stabilwax (60 m × 0.25 mm × 0.25 µm) column in the first dimension coupled to a Rxi-5 ms (2 m × 0.1 mm × 0.1 µm) column in the second dimension was investigated [27,30,31,65]. The result is shown in Figure 2. It should be noted that modulation of 18 s is the best choice to avoid wrap-around phenomena and to obtain a satisfactory chromatogram.
As shown in Figure 2, the chosen internal standard 3-chlorothiophene does not overlap with any components present in the sample, which is confirmed by a high bi-dimensional resolution of 11.26. As such, the quantitative and qualitative performances on FID and MS of the sample are more accurate. A comparison of average bi-dimensional resolutions of the abundant components in the diesel fraction is given in Table 3. The results reveal that the group-type separation is improved remarkably not only for diolefins, iso-olefins, and α-olefins but also for n-paraffins, iso-paraffins, and monoaromatics, whereas the separation of naphthenes is deteriorated but still acceptable compared to the normal-phase combination. This finding is consistent with the result reported by Ristic et al. in an effort to improve the separation of unsaturated hydrocarbons and oxygen-containing compounds in a pyrolysis shale oil sample in which a polar × non-polar column set was applied [27]. Although both normal-phase and reversed-phase combinations resulted in high intraclass resolutions, while the latter combination could separate the oxygen-containing compounds from the hydrocarbon matrix, which were confirmed by the elemental analysis. The oxygenates did not arise in groups (or isomers), but rather sporadically, e.g., 2-ethyl-1hexanol, 2-butyltetrahydrofuran, pentanal, and hexanal with a weight percent of 0.03 wt% corresponding to 0.01 wt% of elemental oxygen. These compounds were co-eluted with the high hydrocarbon content, using the normal-phase mode. The identification of these individual compounds is based on GC × GC-MS analysis combined with the NIST library. The tentative identification of resulting peaks from GC × GC-FID was achieved by a reproducible analysis using GC × GC-MS with the identical column combination and an optimized carrier gas flow [61,62]. The result showed that the abundant chemical groups found in the diesel fraction are n-paraffins, α-olefins, iso-paraffins, iso-olefins, diolefins, naphthenes, and aromatics, which eluted and separated in the group-tile structure on the chromatogram. However, the preliminary result of the elemental analysis revealed that a small amount of oxygen-containing compounds is present in the sample, which could not be detected by using the normal-phase GC × GC combination. This could be explained by the co-elution of oxygenates and the abundant hydrocarbon matrix [27]. The biggest hurdle is the identification of diolefins due to the presence of unsaturated naphthenes, which have the same molecular mass and very similar mass fragmentation patterns on GC × GC-MS chromatogram with diolefins. Nevertheless, these diolefins are not completely separated from iso-olefins, and also neighboring peaks of the high amounts of both paraffins and olefins, which could cause an error in terms of the quantification. As aforementioned, the content of olefins and diolefins results in negative effects on quality control, as well as the technical process of fuels and feedstocks due to the gum-formation property through chemical secondary reaction (e.g., Diels-Alder reaction) [32,63]. Moreover, the heat value of plastic pyrolysis oils is reduced due to the presence of this unsaturated groups [64].
Taking advantages of reversed-phase column configuration, a Stabilwax (60 m × 0.25 mm × 0.25 μm) column in the first dimension coupled to a Rxi-5 ms (2 m × 0.1 mm × 0.1 μm) column in the second dimension was investigated [27,30,31,65]. The result is shown in Figure 2. It should be noted that modulation of 18 s is the best choice to avoid wraparound phenomena and to obtain a satisfactory chromatogram. As shown in Figure 2, the chosen internal standard 3-chlorothiophene does not overlap with any components present in the sample, which is confirmed by a high bidimensional resolution of 11.26. As such, the quantitative and qualitative performances on FID and MS of the sample are more accurate. A comparison of average bi-dimensional resolutions of the abundant components in the diesel fraction is given in Table 3. The results reveal that the group-type separation is improved remarkably not only for diolefins, iso-olefins, and α-olefins but also for n-paraffins, iso-paraffins, and monoaromatics, whereas the separation of naphthenes is deteriorated but still acceptable compared to the normal-phase combination. This finding is consistent with the result It can be concluded that the reversed-phase set (Stabilwax and Rxi-5 ms) considerably improved the separation of the hydrocarbons (n-paraffins, α-olefins, diolefins, iso-olefins, iso-paraffins, and monoaromatics), as illustrated by an improvement in average bi-dimensional resolutions, and the set could identify oxygen-containing compounds compared to the normal-phase configuration (RTX-1 PONA and BPX-50). Therefore, the quantification of chemical components in the diesel fraction described in Figure 3 was obtained by using the GC × GC separation employing polar × non-polar column set. The diesel fraction consists of n-paraffins (28.36 wt%), α-olefins (36.39 wt%), and mononaphthenes + unsaturated mononaphthenes (14.62 wt%) as the most prominent chemical groups. Furthermore, it is composed of significant amounts of iso-paraffins (4.05 wt%), iso-olefins (8.77 wt%), diolefins (4.21 wt%), and dinaphthenes (2.30 wt%). In addition, small amounts of monoaromatics (0.72 wt%) and naphthenoaromatics (0.53 wt%) were detected in this fraction. Compared to hydrocarbons, only a small content of oxygenates was spotted with 0.03 wt%, corresponding to 0.01 wt% of total elemental oxygen in the EA result. The composition of the diesel sample is shown in Figure 3. The detailed chemical groups in the diesel fraction can be found in Supplementary Materials Table S1.

Plastic-Waste Naphtha Pyrolysis Oil
Analogously, the methodology applied for the diesel fraction was also exploited on the naphtha fraction. More specifically, a combination of a Stabilwax (60 m × 0.25 mm × 0.25 μm) column in the first dimension coupled to a Rxi-5 ms (2 m × 0.1 mm × 0.1 μm) column in the second dimension was examined on the naphtha fraction. The obtained result is shown in Figure 4.

Plastic-Waste Naphtha Pyrolysis Oil
Analogously, the methodology applied for the diesel fraction was also exploited on the naphtha fraction. More specifically, a combination of a Stabilwax (60 m × 0.25 mm × 0.25 µm) column in the first dimension coupled to a Rxi-5 ms (2 m × 0.1 mm × 0.1 µm) column in the second dimension was examined on the naphtha fraction. The obtained result is shown in Figure 4.

Plastic-Waste Naphtha Pyrolysis Oil
Analogously, the methodology applied for the diesel fraction was also exploited on the naphtha fraction. More specifically, a combination of a Stabilwax (60 m × 0.25 mm × 0.25 μm) column in the first dimension coupled to a Rxi-5 ms (2 m × 0.1 mm × 0.1 μm) column in the second dimension was examined on the naphtha fraction. The obtained result is shown in Figure 4. This combination revealed an acceptable separation for diolefins from iso-olefins and unsaturated mononaphthenes with the same carbon number, which was confirmed by the bi-dimensional value of 1.87 (Table 4). Furthermore, the analysis could also detect the oxygenates existing in the naphtha fraction from the hydrocarbon matrix, which was disclosed from the elemental analysis result. Besides, the chosen internal standard (3chlorothiophene) did not overlap with any components in the sample, whose bidimensional resolution is equal to 16.48. However, the presence of light hydrocarbons with carbon number distribution from C5 to C11 in the fraction, which eluted quite fast This combination revealed an acceptable separation for diolefins from iso-olefins and unsaturated mononaphthenes with the same carbon number, which was confirmed by the bi-dimensional value of 1.87 (Table 4). Furthermore, the analysis could also detect the oxygenates existing in the naphtha fraction from the hydrocarbon matrix, which was disclosed from the elemental analysis result. Besides, the chosen internal standard (3-chlorothiophene) did not overlap with any components in the sample, whose bi-dimensional resolution is equal to 16.48. However, the presence of light hydrocarbons with carbon number distribution from C5 to C11 in the fraction, which eluted quite fast on the reversed-phase combination, gave a less effective separation of hydrocarbon classes, such as n-paraffins and α-olefins. As a result, the peak identification in the early elution region of the resulting chromatogram corresponding to the major distribution of n-paraffins, α-olefins, iso-olefins, and iso-paraffins groups could be wrongly assigned. To cope with this chromatographic obstacle, many operating conditions were adjusted, such as performing the analysis with a lower ramp temperature or setting the carrier gas (He) at a slower rate. Although many efforts were paid, the resulting chromatograms have not been improved. The outcome above has shown that the use of a non-polar RTX-1 PONA column (50 m × 0.25 mm × 0.25 µm) in the first dimension with a mid-polar BPX-50 column (2 m × 0.18 mm × 0.18 µm) in the second dimension could provide a reasonable separation of the abundant hydrocarbon groups in the diesel fraction; this column set was, thus, tested on this naphtha fraction [27]. The result is presented in Figure 5. The outcome above has shown that the use of a non-polar RTX-1 PONA column (50 m × 0.25 mm × 0.25 μm) in the first dimension with a mid-polar BPX-50 column (2 m × 0.18 mm × 0.18 μm) in the second dimension could provide a reasonable separation of the abundant hydrocarbon groups in the diesel fraction; this column set was, thus, tested on this naphtha fraction [27]. The result is presented in Figure 5. The result showed a considerable improvement in the separation of chemical groups. The separability of two-column combinations for the naphtha fraction was evaluated for the most abundant chemical groups given in Table 4. In addition, a high value of bidimensional resolution of 10.39 was found for the internal standard (3-chlorothiophene).
The use of a normal-phase column set reveals a considerable separation improvement on n-paraffins, α-olefins, and monoaromatics. Furthermore, the reversed- The result showed a considerable improvement in the separation of chemical groups. The separability of two-column combinations for the naphtha fraction was evaluated for the most abundant chemical groups given in Table 4. In addition, a high value of bi-dimensional resolution of 10.39 was found for the internal standard (3-chlorothiophene).
The use of a normal-phase column set reveals a considerable separation improvement on n-paraffins, α-olefins, and monoaromatics. Furthermore, the reversed-phase column could separate oxygen-containing compounds from the abundant hydrocarbon matrix, while only benzofuran and butoxylbenzene were detected by using the normal-phase column set. It should be underlined that many other column combinations, as well as different analysis operations, were tested to obtain a satisfactory separation of all chemical groups present in the light naphtha fraction. For example, the use of mid-polar combinations of IL-59 × Rxi-5ms or BPX-50 × Rxi-5ms, or the use of longer BPX-50 column (4 m in length instead of 2 m) in the second dimension in normal-phase mode (RTX-1 PONA × BPX-50) was performed at various analytical conditions. However, these attempts could not give better results.
In a conclusion, the use of the reversed-phase column combination provides a significant separation capacity for diolefins and oxygenates in both diesel and naphtha fractions. However, this combination bears some limitations for the lighter naphtha fraction. Therefore, the quantification of the naphtha fraction is more accurate when using GC × GC-FID in the normal-phase combination (RTX-1 PONA × BPX-50). The naphtha fraction contains n-paraffins (15.22 wt%), α-olefins (35.08 wt%), and mononaphthenes (19.87 wt%) as the abundant chemical classes. The following contributions are of monoaromatics (9.10 wt%), iso-olefins (9.06 wt%), dinaphthenes (5.52 wt%), diolefins (4.20 wt%), and iso-paraffins (1.80 wt%). Moreover, a small number of naphthenoaromatics and oxygen-containing compounds (benzofuran and butoxylbenzene) were detected with values of 0.08 wt% and 0.06 wt%, respectively. The chemical groups in the naphtha fraction are demonstrated in Figure 6 with the carbon region between C5 and C11. The detailed chemical groups in the naphtha fraction can be found in Supplementary Materials Table S2. monoaromatics (9.10 wt%), iso-olefins (9.06 wt%), dinaphthenes (5.52 wt%), diolefins (4.20 wt%), and iso-paraffins (1.80 wt%). Moreover, a small number of naphthenoaromatics and oxygen-containing compounds (benzofuran and butoxylbenzene) were detected with values of 0.08 wt% and 0.06 wt%, respectively. The chemical groups in the naphtha fraction are demonstrated in Figure 6 with the carbon region between C5 and C11. The detailed chemical groups in the naphtha fraction can be found in Supplementary Materials Table S2. As aforementioned, the implementation of the reversed-phase column combination could distinguish the oxygenates from the abundant hydrocarbon groups. As a result, the content determination thereof in the reversed-phase column combination is better compared to the normal-phase column set. The naphtha fraction contains an amount of alcohol, ketone, and ether chemical classes, which occupy 0.84 wt%, corresponding to 0.14 wt% of elemental oxygen. Because of the limitations of both the column combinations in the comprehensive quantification of chemical classes in the naphtha fraction, two-column sets could be thus employed to complement each other depending on the analytical purpose. As aforementioned, the implementation of the reversed-phase column combination could distinguish the oxygenates from the abundant hydrocarbon groups. As a result, the content determination thereof in the reversed-phase column combination is better compared to the normal-phase column set. The naphtha fraction contains an amount of alcohol, ketone, and ether chemical classes, which occupy 0.84 wt%, corresponding to 0.14 wt% of elemental oxygen. Because of the limitations of both the column combinations in the comprehensive quantification of chemical classes in the naphtha fraction, two-column sets could be thus employed to complement each other depending on the analytical purpose.

GC × GC-SCD
GC × GC-SCD was used to quantify sulfur-containing compounds. Figure 7 shows the GC × GC-SCD color plot of the diesel and naphtha fractions by the use of RTX-1 PONA column (50 m × 0.25 mm × 0.25 µm) in the first dimension, combined with BPX-50 column (2 m × 0.18 mm × 0.18 µm) in the second dimension. The GC × GC parameters are given in Table 2. Figure 7a indicates a small amount of carbon disulfide, sulfides, and C1 thiophenes found in the diesel fraction with the sulfur content of 6.83 ppm, corresponding to 0.002 wt%. On the other hand, the carbon disulfide, sulfides, C1 thiophenes, C2 thiophenes, C3 thiophenes, and C4 thiophenes were detected in the naphtha fraction with the higher sulfur content of 210.74 ppm corresponding to 0.04 wt% in Figure 7b. In Figure 7b, carbon disulfide lost the modulation on the second dimension; it could be explained by its low boiling point and its presence in a higher concentration compared to other sulfur-containing compounds. However, the quantification thereof would not be an error, because it is wellseparated from others. It should be noted that the use of GC × GC-SCD could reveal the presence of a small amount of sulfur-containing compounds, which overlap with the abundant hydrocarbon matrix, using GC × GC-FID, and also could not be detected by the elemental analysis. The detailed information on sulfur-containing compounds in the two diesel and naphtha fractions is given in Table 5.

GC × GC-SCD
GC × GC-SCD was used to quantify sulfur-containing compounds. Figure 7 shows the GC × GC-SCD color plot of the diesel and naphtha fractions by the use of RTX-1 PONA column (50 m × 0.25 mm × 0.25 μm) in the first dimension, combined with BPX-50 column (2 m × 0.18 mm × 0.18 μm) in the second dimension. The GC × GC parameters are given in Table 2.  Figure 7a indicates a small amount of carbon disulfide, sulfides, and C1 thiophenes found in the diesel fraction with the sulfur content of 6.83 ppm, corresponding to 0.002 wt%. On the other hand, the carbon disulfide, sulfides, C1 thiophenes, C2 thiophenes, C3 thiophenes, and C4 thiophenes were detected in the naphtha fraction with the higher sulfur content of 210.74 ppm corresponding to 0.04 wt% in Figure 7b. In Figure 7b, carbon disulfide lost the modulation on the second dimension; it could be explained by its low boiling point and its presence in a higher concentration compared to other sulfurcontaining compounds. However, the quantification thereof would not be an error, because it is well-separated from others. It should be noted that the use of GC × GC-SCD could reveal the presence of a small amount of sulfur-containing compounds, which overlap with the abundant hydrocarbon matrix, using GC × GC-FID, and also could not be detected by the elemental analysis. The detailed information on sulfur-containing compounds in the two diesel and naphtha fractions is given in Table 5.

GC × GC-NCD
In addition to oxygen-and sulfur-containing compounds, nitrogen compounds are also commonly found in plastic-waste pyrolysis oils [16,24]. Analogous to sulfur-containing compounds, the nitrogen-containing compounds often overlap with the abundant hydrocarbon matrix when applying GC × GC-FID due to their low concentration. Therefore, an elemental NCD detector coupled with GC × GC was implemented for the detection of nitrogen-containing compounds in both diesel and naphtha fractions. A summary of GC × GC parameters and column combinations are presented in Table 2. For the diesel fraction, no nitrogen-containing compound was detected in GC × GC-NCD chromatogram. Figure 8 shows the GC × GC-NCD color plot of the naphtha fraction by the use of RTX-1 PONA column (50 m × 0.25 mm × 0.25 µm) in the first dimension, combined with BPX-50 column (2 m × 0.18 mm × 0.18 µm) in the second dimension.
In addition to oxygen-and sulfur-containing compounds, nitrogen compounds are also commonly found in plastic-waste pyrolysis oils [16,24]. Analogous to sulfurcontaining compounds, the nitrogen-containing compounds often overlap with the abundant hydrocarbon matrix when applying GC × GC-FID due to their low concentration. Therefore, an elemental NCD detector coupled with GC × GC was implemented for the detection of nitrogen-containing compounds in both diesel and naphtha fractions. A summary of GC × GC parameters and column combinations are presented in Table 2. For the diesel fraction, no nitrogen-containing compound was detected in GC × GC-NCD chromatogram. Figure 8 shows the GC × GC-NCD color plot of the naphtha fraction by the use of RTX-1 PONA column (50 m × 0.25 mm × 0.25 μm) in the first dimension, combined with BPX-50 column (2 m × 0.18 mm × 0.18 μm) in the second dimension. It can be added that the confirmation of retention time of nitrogen compounds in the naphtha fraction was performed through the injection of a standard mixture containing pyridine, 2-methylpyridine, 2-ethylpyridine, and 5-ethyl-2-methylpyridine. The contribution of nitrogen-containing compounds to the naphtha fraction caused by pyridine, C1 pyridines, C2 pyridines, and C3 pyridines with a total amount of 30.46 ppm, corresponding to 0.05 wt%, is seen in Table 6. It can be added that the confirmation of retention time of nitrogen compounds in the naphtha fraction was performed through the injection of a standard mixture containing pyridine, 2-methylpyridine, 2-ethylpyridine, and 5-ethyl-2-methylpyridine. The contribution of nitrogen-containing compounds to the naphtha fraction caused by pyridine, C1 pyridines, C2 pyridines, and C3 pyridines with a total amount of 30.46 ppm, corresponding to 0.05 wt%, is seen in Table 6.

Elemental Analysis
The CHNSO contents from EA analysis are compared with results obtained from the GC × GC-FID/SCD/NCD analysis in Table 7. It was shown that the hydrogen-to-carbon ratio (H/C) of the EA analysis is consistent with the result obtained from using GC × GC-FID. In particular, the H/C ratio of the diesel fraction is 2.02 for GC × GC and 2.03 for EA. The corresponding H/C values of the naphtha fraction are identical (1.93) in both analytical methods. Although the fact that the naphtha fraction of a fossil origin has the H/C value higher than a diesel fraction of the same origin, that is not in the case when these fractions come from plastic pyrolysis oils [66,67]. Furthermore, the oxygen elemental composition determined by EA is higher in comparison with the obtained results of GC × GC-FID in both studied samples. The smaller amount of oxygen components in the GC × GC results could be explained by the interference of the prominent hydrocarbon matrix in both naphtha and diesel fractions. Moreover, the presence of sulfur and nitrogen compounds in both samples could not be detected by EA, because their concentrations are below the instrumental detection limit. However, the use of GC × GC coupled to the SCD and NCD elemental detectors could allow for the determination of sulfur and nitrogen components at ppm levels.
The results have shown that the saturated property of hydrocarbons, as well as the concentration of heteroatom-containing compounds, is significantly improved in these diesel and naphtha fractions in comparison with plastic-waste pyrolysis oil without further upgrading the process [16].

Conclusions
A new method for the quantification of plastic-waste pyrolysis oils was developed using a two-dimensional gas chromatography GC × GC coupled to four different detectors (FID, SCD, NCD, and MS). With the StabilWax × Rxi-5ms/RTX-1 PONA × BPX-50 column combination, it was possible to accurately characterize the chemical compositions of two plastic-waste pyrolysis oils with different boiling points. One of the strengths of the developed method is the separation of different types of olefins by optimizing the column configuration and operating conditions. A combination of a polar × non-polar column configuration (StabilWax × Rxi-5ms) was shown as the best choice for the diesel fraction in terms of separability of diolefins from iso-olefins and mononaphthenes and oxygen-containing compounds from a complicated hydrocarbon matrix. This did not affect the separation of the other abundantly present classes of compounds, e.g., n-paraffins, iso-paraffins, and monoaromatics. Therefore, this column combination can be used for quantitation, even for samples containing small amounts of oxygenates.
Alternatively, using the non-polar × polar (RTX-1 PONA × BPX-50) combination was superior for the group-type separation of the naphtha fraction. With the selective SCD and NCS detectors, it was possible to detect trace amounts of sulfur and nitrogen compounds.