# Characterization of Pyrolysis Oils Using a Combination of GCxGC/TOFMS and GC/HRMS Analysis: The Impact of Data Processing Parameters

Xiangdong Chen, Carlos Rincon, Benoît Gadenne, José Dugay, Michel Sablier,

Jérôme Vial

► To cite this version:

Xiangdong Chen, Carlos Rincon, Benoît Gadenne, José Dugay, Michel Sablier, et al.. Characterization of Pyrolysis Oils Using a Combination of GCxGC/TOFMS and GC/HRMS Analysis: The Impact of Data Processing Parameters. Separations, 2025, 12 (9), pp.239. 10.3390/separations12090239. hal-05325049

HAL Id: hal-05325049

https://hal.science/hal-05325049v1

Submitted on 21 Oct 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Article

# Characterization of Pyrolysis Oils Using a Combination of GC×GC/TOFMS and GC/HRMS Analysis: The Impact of Data Processing Parameters

Xiangdong Chen<sup>1,2,*</sup> (D), Carlos Rincon<sup>2</sup>, Benoît Gadenne<sup>2</sup>, José Dugay<sup>1</sup>, Michel Sablier<sup>1,D</sup> and Jérôme Vial<sup>1,D</sup>

## Abstract

Human population growth and increasing transportation demands have led to rising global tire consumption and associated waste. In response, various material and energy recovery strategies, such as pyrolysis, have been developed to produce high-value-added products such as pyrolysis oils, which can be reused as materials or fuels. However, these oils often contain heteroatom-containing compounds (e.g., nitrogen, oxygen, sulfur) that can hinder their valorization and must therefore be identified and removed. To characterize heteroatomic compounds present in distillation fractions of pyrolysis oils, GC×GC/TOFMS and GC/HRMS were employed. For non-target analysis, data processing parameters were optimized using a Central Composite Design (CCD). The most influential parameters for GC×GC/TOFMS were the minimum number of mass-to-charge ratio (m/z) signals kept in the deconvoluted spectra (minimum stick count) and peak signal-to-noise ratio (S/N), while for GC/HRMS, optimization focused on the m/z S/N threshold, peak S/N, and total ion current (TIC). Under optimal conditions, 129 and 92 heteroatomic compounds were identified via GC×GC/TOFMS and GC/HRMS, respectively, within a single distillation fraction, with 57 compounds identified using both techniques. Notably, GC×GC/TOFMS exclusively identified 72 compounds, while there were only 5 unique to GC/HRMS. These results highlight the effectiveness of GC×GC/TOFMS in characterizing heteroatomic compounds in complex mixtures, while also underlining the complementary value of GC/HRMS.

**Keywords:** central composite design; tire recycling; data processing; pyrolysis

## 1. Introduction

Currently, different types of pollution, having a negative impact on the environment, are generated by transport (whether public or private). They do not only consist of air pollution (the emission of CO₂ or SO₂ by using fuels for instance), but also of industrial waste, especially tire waste. With the growth of the human population, and consequently the growth of transport use, the demand for tires has increased around the world. In fact, the production rate of all types of tires is continuing to grow as their demand in developed countries is rapidly increasing. In 2023, the global consumption of tires and tire products

<figure><img src="image_2.png" /></figure>

check for
updates

Academic Editor: Victoria Samanidou

Received: 23 July 2025

Revised: 1 September 2025

Accepted: 2 September 2025

Published: 4 September 2025

**Citation:** Chen, X.; Rincon, C.; Gadenne, B.; Dugay, J.; Sablier, M.; Vial, J. Characterization of Pyrolysis Oils Using a Combination of GC×GC/TOFMS and GC/HRMS

**Analysis:** The Impact of Data Processing Parameters. *Separations* **2025**, *12*, 239. https://doi.org/10.3390/separations12090239

**Copyright:** © 2025 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

<footnote n="1">¹ UMR CNRS CBI, Laboratoire des Sciences Analytiques, Bioanalytiques et Miniaturisation, ESPCI Paris, PSL Research University, 10 Rue Vauquelin, CEDEX 05, 75231 Paris, France; jose.dugay@espci.fr (J.D.); michel.sablier@espci.fr (M.S.); jerome.vial@espci.fr (J.V.)</footnote>

<footnote n="2">² Manufacture Françaises des Pneumatiques Michelin, Place des Carmes Deschaux, CEDEX 09, 63040 Clermont-Ferrand, France; carlos-andres.rincon-mejia@michelin.com (C.R.); benoit.gadenne@michelin.com (B.G.)</footnote>

<footnote>* Correspondence: xiangdong.chen@espci.fr; Tel.: +33-652-194-740</footnote>
was 17 billion tons, an increase of 1% from 2022 [1]. According to the literature, the annual rates of tire accumulation are considerable: around 1.5 billion waste tires are produced per year after their end of life, which will increase to 5 billion by the end of 2030 [2,3]. Different methods to manage tire waste have been considered in the past few decades, with the most common being landfilling, direct reutilization, retreading, and material and energy recovery [2,4]. With restrictions on waste tire disposal in landfills, 95% of wasted tires were collected and treated for material recycling and energy recovery in Europe in 2019 [5].

Pyrolysis constitutes an alternative tire waste treatment method that can obtain high-value-added products [6,7]: solid residues (such as recovered carbon black and metals), gases (such as methane and hydrogen, which can be recycled as a syngas) and the product of interest in this study, pyrolysis oils, which may be reused as a source of materials and fuel [7–9]. Pyrolysis oils are a complex matrix containing high amounts of aliphatic, aromatic and heterocyclic compounds (usually more than thousands of compounds, some of which have similar boiling points) [7,8,10]. For this reason, the distillation of pyrolysis oils is commonly performed to simplify the matrix and to obtain different fractions for targeted valorization processes.

Several analytical methods have been considered to characterize the chemical composition of these samples. Gas chromatography coupled with low-resolution mass spectrometry (GC/MS) [6,11,12] and liquid chromatography coupled with low-resolution mass spectrometry (LC/MS) were applied to understand the composition of the distillation fractions of pyrolysis oils and evaluate their applications [13,14]. However, the analysis of these samples with these conventional methods is complicated. In order to obtain a detailed characterization of the pyrolysis oils of waste tires, high-resolution MS and/or high chromatographic resolution is necessary. Zhang et al. (2023) have analyzed light fractions of waste tire pyrolysis oil (WTPO) using GC/MS and a GC/nitrogen chemiluminescence detector (NCD) in order to characterize the nitrogen compounds present in WTPO [15]. Due to the high complexity of heavy distillation fractions, Fourier-Transform Ion Cyclotron Resonance Mass Spectrometry (FT-ICR MS) was employed to detect nitrogen compounds instead of GC/MS [15]. Li et al. (2023) demonstrated the advantages of Comprehensive two-Dimensional Gas Chromatography combined with Time of Flight Mass Spectrometry (GC×GC/TOFMS) over GC/MS for analyzing heavy lubricant oil due to its high separation capacity [16]. Furthermore, the authors compared the characterization results of GC×GC/TOFMS and Atmospheric Pressure Chemical Ionization (APCI)-Orbitrap MS in the same study [16].

More specifically, the detailed characterization of the volatile fraction of pyrolysis oils (plastic waste [17] and crude tire [18]) has previously been conducted using GC×GC/TOFMS and Gas Chromatography combined with High-Resolution Mass Spectrometry (GC/HRMS) [19]. GC×GC increases peak capacity by utilizing the ideal orthogonality of two columns; for instance, a non-polar column in the first dimension to separate compounds according to their volatility and a polar column in the second dimension to separate them according to polarity [20,21]. Reversed configuration was also considered to complete the normal configuration. The choice of column configuration should be made with care: Mahé et al. (2011) compared four normal configurations with four reverse configurations to determine which combination was the most well adapted for heavy petroleum cuts [22]. Therefore, compounds that coeluted in one-dimensional gas chromatography (GC-1D) can often be resolved in GC×GC, making it particularly useful for complex mixtures. Moreover, due to the cryogenic modulation process, compounds are collected and focused into narrower and more concentrated modulated peaks that are reinjected in the second column, hence increasing the sensitivity of the analysis compared with
GC-1D [23]. Therefore, trace compounds can be more easily detected. HRMS allows the determination of the *m/z* of ions with high resolving power and mass accuracy: ions with very close *m/z* values can be distinguished, making it efficient for the identification of unknown compounds, isotopic patterns and elemental compositions [24,25]. HRMS is usually coupled with GC for the targeted and untargeted analysis of volatile and semi-volatile compounds, including those containing heteroatoms such as nitrogen, oxygen and sulfur [26,27]. According to previous studies, heteroatomic compounds present in pyrolysis oils hinder their valorization. For example, sulfur and nitrogen compounds are sources of toxic and corrosive gas such as SOₓ and NOₓ [7, 11, 15, 28, 29].

In this study, GC×GC/TOFMS and GC/HRMS analytical protocols were used for the untargeted analysis of pyrolysis oils from car waste tires, with a particular focus on the characterization of heteroatom-containing compounds. Two distillation fractions of pyrolysis oil were characterized by this analytical strategy. After chromatographic acquisition, the data processing software, ChromaTOF for GC×GC/TOFMS and TraceFinder for GC/HRMS, were used to deconvolute coeluted peaks. Here, different software parameters were studied and optimized in an innovative manner as they can have a critical impact on the detection of trace-level compounds. Until now, the determination of these parameters has rarely been discussed in detail in the literature, with several articles roughly reporting the applied values [21,30]. Therefore, in the present study, the key parameters for each data processing software were identified by a CCD that could assess the influence and interactions of software parameters and model the responses in view of optimization [31–33]. The number of detected peaks and true positives were considered as responses. The optimization criteria were to minimize the number of detected peaks and maximize the true positives for identification. Using the identified optimal settings, compound identification of selected fractions was performed with each software. The resulting peak attributions were then compared to evaluate which approach, GC×GC/TOFMS or GC/HRMS, was more suitable for the untargeted analysis of heteroatom-containing compounds in the pyrolysis oils.

## 2. Materials and Methods

### 2.1. Sample Preparation

Tire pyrolysis oil was purchased from Enviro Systems (Gothenburg, Sweden). This oil was subjected to distillation, yielding eight fractions with boiling points ranging from 95 °C to 549 °C. Fractions 369 °C to 509 °C and 509 °C to 549 °C could not be easily analyzed by GC×GC/TOFMS: the nitrogen cryogenic modulator can trap the heaviest molecules and these two heavy fractions risked being partially eluted. The six remaining fractions were considered for this study. Among them, two light fractions named fraction A (95–142 °C) and B (142–169 °C) were chosen. Fraction A and B could be considered representative of pyrolysis oil fractions in terms of their complexity and were analyzed using GC×GC/TOFMS and GC/HRMS. The distillation fractions were diluted tenfold with dichloromethane before use. A C₈–C₄₀ alkane mixture (Sigma-Aldrich, Saint-Quentin-Fallavier, France) was used to determine the linear retention indices of the identified compounds. As the original concentration of the standard (1000 ppm) was too high, it was diluted 20-fold in dichloromethane before use.

### 2.2. GC×GC/TOFMS Analysis

Once prepared, the distillation fractions of the pyrolysis oils were analyzed using a GC×GC/TOFMS system (Pegasus BT4D, LECO, Villepinte, France). The Pegasus BT4D combines a TOFMS detector and an Agilent 7890 GC device equipped with a cryogenic modulator and a secondary oven. A normal configuration was chosen. A non-polar
column, Rxi-5ms column (5% diphenyl/95% dimethylpolysiloxane; 30 m × 0.25 mm × 0.25 µm; Restek, Lisses, France) was used in the first dimension (1D). In the second dimension (2D), a medium polar DB-1701 (14% cyanopropyl-phenyl-methylpolysiloxane; 52 cm × 0.18 mm × 0.18 µm; Agilent, Les Ulis, France) was installed in the secondary oven. This column configuration allowed us to separate the compounds according to their volatility. Heteroatomic compounds, the compounds of interest, were located at the top of the 2D chromatogram. Other column configurations could also have been considered, especially the reverse configuration, but this was beyond the scope of the present study. Helium was used as the carrier gas at a constant flow rate of 1 mL/min. The Split/Splitless liquid (SSL) injector was maintained at 300 °C with a split ratio of 100:1. The oven temperature program began at 35 °C (held for 2 min) and increased at 3 °C/min to 250 °C (held for 5 min). The secondary oven followed the same temperature program with a 5 °C offset. Modulation was performed using a QuadJet™ nitrogen cryogenic modulator, operating with a 4 s modulation period composed of two 2 s stages: 1.2 s hot jet and 0.8 s cold jet. After compound separation, analytes were transferred to the TOFMS via a deactivated silica transfer line (31 cm × 0.18 mm) heated to 250 °C. Mass spectra were recorded over a mass range of 45–500 *m/z* with an acquisition rate of 200 spectra/s. To calculate the retention indices, a solution of C₈–C₄₀ alkanes (50 ppm) was injected under the same analytical conditions [34].

## 2.3. GC/HRMS Analysis

For the GC/HRMS analysis of the distillation fractions, a Q-Exactive GC Orbitrap system coupled with a Trace 1310 GC (Thermo Fisher Scientific, Waltham, MA, USA) was employed. A non-polar HP-5ms column (5% diphenyl/95% dimethylpolysiloxane; 30 m × 0.25 mm × 0.25 µm; Agilent) was installed. This column was nearly the same one as the column used in the first dimension of GC×GC/TOFMS in order to provide a similar separation. Helium was used as the carrier gas at a constant flow rate of 1 mL/min. The SSL injector was maintained at 300 °C with a split ratio of 20:1. The oven temperature program started at 40 °C (held for 3 min), increased at a rate of 5 °C/min to 300 °C (held for 5 min). These temperature program settings are quite standard and are routinely used in this technique to allow a reasonable analysis time. After chromatographic separation, analytes were transferred to the HRMS through the non-polar column, of which the last centimeter also served as a heated transfer line maintained at 300 °C. Mass spectra were acquired over a mass range of 50–500 *m/z* with a mass resolution set at 60,000 (corresponding to an acquisition rate of 7 spectra/s). A solution of C₈–C₄₀ alkanes (50 ppm) was injected under the same analytical conditions as the real samples, enabling the calculation of retention indices for compound identification.

## 2.4. Data Processing Software

Once the distillation fractions of the pyrolysis oils were analyzed by GC×GC/TOFMS and GC/HRMS, the resulting chromatograms were processed using dedicated software. Due to the complex nature of pyrolysis oil matrices, many peaks in both one-dimensional and two-dimensional chromatograms were coeluted. Therefore, deconvolution was essential for accurate peak identification. ChromaTOF (v5.55.41, Leco Corporation, St. Joseph, MI, USA) was used for processing the GC×GC chromatograms, while TraceFinder (v4.1) equipped with the deconvolution plugin (v1.7, ThermoFisher Scientific instruments, Waltham, MA, USA) was used for the GC/Orbitrap data. After deconvolution and data processing, untargeted analysis was performed by comparing the obtained mass spectra with entries in the National Institute of Standards and Technology (NIST) mass spectral library (NIST 20). Both software platforms provided preliminary lists of identified compounds,
which were subsequently verified manually using retention indices to ensure accuracy of the identification.

## 2.5. Level of Identification Confidence

In this study, the detection of heteroatom-containing compounds arises from proposed identifications automatically generated through the data processing method. However, their identities must be manually verified using their calculated retention index in both GC×GC/TOFMS and GC/HRMS and/or their position in the 2D chromatogram in the case of GC×GC/TOFMS. A detected compound was considered as a true positive if its mass spectrum and retention index (semi-standard non-polar column) match those of the proposed compound identified by the software. During GC×GC/TOFMS analysis, the second dimension of the 2D chromatograms provides additional information on the polarity of detected compounds, enabling the rejection of several false positives, defined as detected compounds which were in reality not present in the distillation fraction. In contrast, GC/HRMS offers structural information through the exact molecular masses of fragment ions.

To improve the confidence level of the identification, further investigations should be conducted using chemical standards of the identified heteroatomic compounds. This approach was not undertaken in the present study, as it falls beyond the scope of this paper.

## 2.6. Selection of Experimental Design

A CCD was applied to optimize the key parameters of each data processing software. The CCD enabled the quadratic modeling of the experimental responses as functions of the k parameters shown in Equation (1):

$$y = \beta_0 + \sum_{i=1}^{k} \beta_i x_i + \sum_{i=1}^{k} \beta_{ii} x_i^2 + \sum_{i=1}^{k} \sum_{j=i+1}^{k} \beta_{ij} x_i x_j + \varepsilon \quad (1)$$

The CCD design included three types of points: 2k factorial points, 2k axial or star points at a distance α from the center (to estimate curvature), and center points (to estimate experimental error). In our case, α was set to 1 for the sake of simplicity. The responses considered in the CCD were the number of detected heteroatomic compounds and the number of true positives (heteroatomic compounds identified by mass spectrum and retention index comparisons, Section 2.5). The overall objective was to minimize the number of detected peaks (to reduce the number of peaks requiring manual verification) while maximizing the number of true positives (to ensure reliable identification). JMP 18.2.0 was used to build experimental designs and perform statistical analyses.

For GC×GC/TOFMS data, a CCD with 2 factors was established: the minimum signal-to-noise ratio threshold (MSN) and the minimum stick count (MSC). These two parameters were selected based on a preliminary study investigating the influence of all the parameters: MSN and MSC.

* MSN specifies the lowest signal-to-noise ratio required for a signal to be recognized as a true peak. If the MSN is too high, trace compounds cannot be detected. Therefore, an appropriate S/N threshold value should be carefully selected. Boudard et al. (2024) chose 1000 as a minimum S/N threshold for body odor samples [21] while Stephanuto et al. (2017) used an S/N threshold of 100 for profiling volatile aromas of beer upon Thermal desorption(TD)-GC×GC/TOFMS analysis [35]. In the context of TD injection, ambient air contributes significantly to background noise, which justifies the use of a relatively high minimum S/N to balance between the number of detected peaks and the reliability of true positive identifications. However, in the case of liquid injections, the background noise is significantly lower than in TD injection. Consequently, a
lower minimum S/N threshold is more suitable to enable the untargeted detection of heteroatom-containing compounds at the trace level. Bean et al. (2015) used an S/N of 10 to analyze metabolites for biomarker discovery [36], which is comparable to the value used in our study for liquid injection.

* MSC refers to the minimum number of ions in a mass spectrum that must be present for a detected peak to be considered as a valid compound, thereby reducing the risk of false positives due to artefacts. Thus, increasing the MSC leads to a decrease in the number of detected peaks, due to the stricter criteria imposed for peak formation. Bean et al. (2015) used an MSC of 2 for data processing (called apexing masses) [36].

Consequently, MSN and MSC appear to have the most significant impact. Within the frame of the experimental design, 9 experiments were performed to examine the influence of these two parameters on the two responses mentioned above. Table 1 indicates the experiments and factor levels used for GC×GC/TOFMS.

<table><caption>**Table 1.** Experimental table of experiments for the GG × GC/TOFMS parameters and factor levels.</caption><thead><tr><th>Experiment</th><th>MSC</th><th>MSN</th></tr></thead><tbody><tr><td>1</td><td>0</td><td>0</td></tr><tr><td>2</td><td>-1</td><td>1</td></tr><tr><td>3</td><td>1</td><td>-1</td></tr><tr><td>4</td><td>1</td><td>1</td></tr><tr><td>5</td><td>-1</td><td>0</td></tr><tr><td>6</td><td>0</td><td>1</td></tr><tr><td>7</td><td>0</td><td>-1</td></tr><tr><td>8</td><td>-1</td><td>-1</td></tr><tr><td>9</td><td>1</td><td>0</td></tr></tbody><thead><tr><th>Factor level</th><th>MSC</th><th>MSN</th></tr></thead><tbody><tr><td>-1</td><td>3</td><td>5</td></tr><tr><td>0</td><td>5</td><td>55</td></tr><tr><td>1</td><td>7</td><td>105</td></tr></tbody></table>

For GC/HRMS, a CCD with 3 factors was established: the signal-to-noise ratio threshold for mass-to-charge ratio *m/z* (*m/z* S/N), the peak S/N threshold (peak S/N) and the TIC. These three parameters were selected based on a preliminary study investigating the influence of all data processing parameters on TraceFinder.

* The *m/z* S/N threshold corresponds to the minimum ratio between an ion's signal height and the baseline noise for the ion to be included in the mass spectrum. Increasing the *m/z* S/N threshold led to a decrease in the number of detected peaks, due to stricter filtering of low-intensity ions.

* The peak S/N threshold specifies the minimum ratio between the peak height and the baseline noise for a signal to be classified as a valid chromatographic peak.

* The TIC intensity threshold sets the minimum total ion current (TIC) intensity required for a peak to be retained.

The effect of these three factors was studied in 15 experiments. Table 2 shows the experiments and factor levels used for GC/HRMS.
<table><caption>**Table 2.** Experimental table for the GC/HRMS parameters and factor levels.</caption><thead><tr><th>Number of Experiments</th><th>m/z S/N</th><th>TIC</th><th>S/N</th></tr></thead><tbody><tr><td>1</td><td>-1</td><td>1</td><td>-1</td></tr><tr><td>2</td><td>1</td><td>-1</td><td>-1</td></tr><tr><td>3</td><td>0</td><td>0</td><td>1</td></tr><tr><td>4</td><td>0</td><td>1</td><td>0</td></tr><tr><td>5</td><td>0</td><td>0</td><td>0</td></tr><tr><td>6</td><td>-1</td><td>0</td><td>0</td></tr><tr><td>7</td><td>0</td><td>-1</td><td>0</td></tr><tr><td>8</td><td>0</td><td>0</td><td>-1</td></tr><tr><td>9</td><td>1</td><td>-1</td><td>1</td></tr><tr><td>10</td><td>1</td><td>0</td><td>0</td></tr><tr><td>11</td><td>-1</td><td>1</td><td>1</td></tr><tr><td>12</td><td>-1</td><td>-1</td><td>-1</td></tr><tr><td>13</td><td>1</td><td>1</td><td>1</td></tr><tr><td>14</td><td>-1</td><td>-1</td><td>1</td></tr><tr><td>15</td><td>1</td><td>1</td><td>-1</td></tr></tbody><thead><tr><th>Factor level</th><th>m/z S/N</th><th>TIC</th><th>S/N</th></tr></thead><tbody><tr><td>-1</td><td>3</td><td>10<sup>4</sup></td><td>3</td></tr><tr><td>0</td><td>51.5</td><td>5.0005 × 10<sup>7</sup></td><td>6.5</td></tr><tr><td>1</td><td>100</td><td>10<sup>8</sup></td><td>10</td></tr></tbody></table>

To achieve both a minimum number of detected peaks and a maximum number of identified true positives, a desirability approach was applied using Derringer functions to transform each response of interest into a scale from 0 (fully undesirable) to 1 (fully desirable) [37,38]. Optimal parameter values were determined when desirability became maximum. The optimized values were applied to compare the results obtained from GC×GC/TOFMS and GC/HRMS, in order to assess and contrast the analytical performance of both methods on two chosen fractions. It must be noted that in the present case, each assay corresponded to the application of a reprocessing algorithm with different settings but on the same raw data. Replicates would therefore lead to the same mathematical value for the responses and would be of no interest.

## 3. Results and Discussion

### 3.1. GC×GC/TOFMS and GC/HRMS Chromatograms of Fraction A

Figure 1 presents the chromatograms obtained by GC × GC/TOFMS (a) and GC/HRMS (b). The organization of the chemical families was observed in the 2D chromatogram: alkane/alkene compounds, monoaromatic compounds and heteroatom-containing compounds. Examples of heteroatomic compounds are circled in red in Figure 1a. Three main classes of heteroatom-containing compounds were identified: oxygen-, nitrogen-, and sulfur-containing species. Heteroatomic compounds were usually located above hydrocarbon compounds as a result of their higher polarity, which facilitated their identification. However, with GC/HRMS, heteroatomic compounds were usually co-eluted with hydrocarbon compounds and were thus difficult to identify. Only major compounds, such as for example toluene, p-xylene and o-xylene, could still be easily visualized and identified as shown in Figure 1b.
<figure><img src="image_3.png" /></figure>

<figure><img src="image_4.png" /></figure>

**Figure 1.** Chromatograms of pyrolysis oil fraction A (95–142 °C) obtained by (a) GC×GC/TOFMS and (b) GC/HRMS. Three heteroatomic compounds are circled in red in (a): from left to right, red circles identify, respectively, 4-Methyl-2-pentanone, cyclopentanone and benzonitrile. Some major compounds in GC/HRMS chromatogram are presented in (b): Toluene, p-xylene and o-xylene.

### 3.2. CCD for GC×GC/TOFMS

#### 3.2.1. Correlation of the Two Responses

To investigate the correlation between the number of detected heteroatomic compounds and the number of true positives, a correlation plot was performed based on the results of CCD. Figure 2 presents the results with ellipse density contours highlighted at a confidence level of 0.95. The more circular the shape of the ellipse, the weaker the correlation between the two responses. In this case, the number of detected heteroatomic compounds and the number of true positives appeared to be weakly correlated, which cor-
responds to a correlation coefficient of 0.54 (which is not significantly different from 0 with a risk of first kind set at 5%). This indicates that they deserve to be studied independently.

<figure><img src="image_5.png" /><figcaption>**Figure 2.** The correlation of the number of detected heteroatomic compounds and the number of true positives based on the results of CCD of GC×GC on pyrolysis oil fraction A (95–142 °C) using the multivariate correlation of JMP. Ellipse density at confidence level of 0.95 was highlighted in red.</figcaption></figure>

### 3.2.2. Significance of Coefficients for Each Response

By using a CCD, it becomes possible to assess the extent to how each factor influences the two responses in terms of their linear contribution, interaction and curvature. Figure 3 illustrates the significance of these two factors and their interaction on two responses: the number of detected heteroatom-containing peaks and the number of true positives in fraction A. In this case, [MSN] and [MSC] had a significant negative effect and the quadratic term [MSN]² and [MSC] × [MSN] a significant positive effect on the number of detected heteroatomic compounds. In addition, [MSC], [MSN], [MSN]² and [MSC]² had a significant negative effect on the number of true positives, whereas the interaction term [MSC] × [MSN] had a positive effect. When applying the same CCD approach on fraction B, the results obtained were similar to those described for fraction A.
<figure><img src="image_6.png" /><figcaption>**Figure 3.** GC×GC/TOFMS CCD results: standardized values of the coefficients of the model for (a) the number of detected peaks and (b) the number of true positives. The effects of MSN, MSC, their interaction, and their quadratic terms on both responses were evaluated by the t-ratio using the scaled estimates. The red lines represent the 5% significance threshold. When a coefficient exceeds the red line, the corresponding coefficient is considered statistically significant.</figcaption></figure>

### 3.2.3. Optimization with Desirability for Fraction A

To achieve the best compromise between minimizing the number of detected heteroatom-containing compounds and maximizing the number of true positives, a desirability approach using Derringer functions was employed. Desirability was defined to range linearly from 0 (least desirable outcome) to 1 (most desirable), with a value of 1 corresponding to a minimum number of detected heteroatomic peaks and a maximum number of true positives. Figure 4 shows the optimal point that satisfies both objectives simultaneously. The maximum overall desirability was reached at an MSC of 3.936 and an MSN of 60. Given that the MSC and MSN values must be integers in practice, the MSC value was rounded up to 4 and the MSN value to 60 to identify the heteroatomic compounds. These parameters values were also used as the basis for comparison with results obtained from GC/HRMS.
<figure><img src="image_7.png" /><figcaption>**Figure 4.** Variation in the predicted number of detected heteroatomic compounds and the predicted number of true positives as a function of MSN and MSC by modeled responses, using prediction profiler of JMP based on the results of CCD of GC×GC. The last row indicates the desirability as a function of MSC and MSN, with optimum values in red (MSC = 3.936, MSN = 60), and the corresponding 95% confidence interval. The last column represents the individual desirability values for each response. The red points represent the optimal parameter settings that yielded the maximum desirability. The four blue crosses represent the values of experimental responses at the optimum values of MSC and MSN: experimental values were all included in the 95% confidence interval.</figcaption></figure>

### 3.2.4. CCD Relevance and Experimental Validation for Fraction A

To illustrate the relevance of using a CCD—and thus a quadratic model—response surfaces were generated for each response. Figure 5 displays the response surfaces for the number of detected heteroatom-containing compounds (Figure 5a) and for the number of true positives (Figure 5b). In both cases, the surfaces deviated significantly from a flat plane and exhibited curvature characteristic of quadratic variation, indicating that the response behavior was not linear. This result demonstrates the relevance of using CCD as an experimental design strategy. In addition, the R square value for the number of detected heteroatomic compounds was 0.9763, and for the number of true positives was 0.9851, which showed the high performance of predicted model. Residual plots are provided in Figure S2: the residuals are randomly distributed around the line y = 0.
<figure><img src="image_8.png" /><figcaption>**Figure 5.** Response surface plots of (a) the number of heteroatomic compounds detected as a function of MSC and MSN and (b) the number of true positives as a function of MSC and MSN using the surface profiler in JMP based on the results of CCD of GC×GC. The color scale ranges from purple to red, where purple indicates the minimum response and red the maximum.</figcaption></figure>

The optimized parameter values predicted by the model were then validated experimentally. As shown in Figure 4, the predicted number of detected heteroatomic compounds was 3082 and of true positives was 129.6 at an MSC of 3.936 (rounded up to 4 for the subsequent study) and an MSN of 60, with 95% confidence intervals of [2412.6, 3751.6] and [127.5, 131.6], respectively. When applying these settings, 2848 heteroatom-containing compounds were detected, and 129 true positives corresponding to heteroatom-containing compounds were identified. Figure 4 highlights this result, with the blue crosses corresponding to the experimental validation at the optimum set of values for MSC and MSN. This result closely matched the predicted value, thereby validating the chosen optimization strategy and optimum.
### 3.3. CCD for GC/HRMS Data of Fraction A

The same strategy was applied to data obtained by GC/HRMS in order to optimize software parameters. Detailed results are shown in Supplementary Material. Figure S1 illustrates the significance of these three factors and their interaction on two responses. However, several results presented a behavior different from those previously observed for GC×GC/TOFMS data. As in Section 3.2.1, a multivariate analysis was conducted, and the results are represented by a scatterplot matrix. Figure 6 shows the correlation plot of the two responses based on the results of CCD of GC/HRMS, with the highlighted ellipse density at a confidence level of 0.95. Compared to Figure 2, the two responses seem to correlate more in the case of GC/HRMS, as the shape of the ellipse was flattened. Therefore, the behavior of these two responses should be rather similar in this case.

<figure><img src="image_9.png" /><figcaption>**Figure 6.** The correlation of the number of detected heteroatomic compounds and the number of true positives based on the results of CCD of GC/HRMS on pyrolysis oil fraction A (95–142 °C) using the multivariate correlation of JMP. Ellipse density at confidence level of 0.95 is highlighted in red. The black dots represent the experimental data points selected in JMP.</figcaption></figure>

As a consequence, the two responses decreased in a correlated manner as a function of the parameters. According to the desirability curve (Figure 7), it was not possible to identify a set of conditions that would simultaneously maximize the number of true positives and minimize the number of detected heteroatomic compounds. As a result, the most optimal compromise was found at the following parameter settings: m/z S/N = 3, TIC intensity threshold = 104, and peak S/N = 3, which provided the maximum true positives (predicted value of 89 with a confidence interval of [79.5, 98.8]) and the maximum detected heteroatomic compounds to inspect manually (predicted value of 1230 with a 95% confidence interval of [989.2, 1470.9]). Experimentally, the set of values (m/z S/N = 3, TIC = 104, S/N = 3) provided 1306 detected heteroatomic compounds and 92 true positives. Figure 6
highlights this result with six blue crosses at the optimum set of values for MSC and MSN.
Both values were in the 95% confidence interval, which validated the chosen model.

<figure><img src="image_10.png" /></figure>

**Figure 7.** Variation in the number of detected heteroatomic compounds and the number of true positives as a function of *m/z* S/N, TIC and S/N predicted by modeled responses, using the prediction profiler of JMP based on the results of CCD of GC/HRMS. The last row indicates the desirability as a function of *m/z* S/N, TIC intensity and S/N represented at optimum values (m/z S/N = 3, TIC = 10⁴ and S/N = 3) and the corresponding 95% confidence level. The last column represents the individual desirability values for each response. The six blue crosses represent the values of experimental responses at the optimum values (m/z S/N = 3, TIC = 10⁴ and S/N = 3): experimental values were all included in the 95% confidence interval.

3.4. *Confirmation of These Optimal Values Through Another Sample: Fraction B (142–169 °C)*

Using the optimized parameter set of MSC = 4 and MSN = 60, 2928 heteroatomic compounds were detected in fraction B (142–169 °C) upon GC×GC/TOFMS analysis, of which 159 were true positives—covering 99.3% of all true positives. In contrast, the MSC = 3 and MSN = 5 setting ((−1, −1) point in CCD) led to the detection of 7020 compounds to recover all 160 true positives. The only compound not detected with (MSC = 4, MSN = 60) was a low-intensity trace compound. This highlights the power of the CCD experimental design, as well as the relevance of the (MSC = 4, MSN = 60) setting, which achieved near-complete identification while reducing the number of detected heteroatom-containing compounds by more than half (58.3%). Table S1 shows the identification of the 160 true positives detected.

For GC/HRMS, the set of values (m/z S/N = 3, TIC = 10⁴, S/N = 3) provided 1390 detected heteroatomic compounds and 154 true positives. With this set, the maximum true positives and the maximum detected heteroatomic compounds were obtained. Consequently, no experimental validation was necessary, as this set of values was a point in the CCD.
3.5. *Comparison of GC×GC/TOFMS and GC/HRMS Results*

Using the optimized parameter sets defined in Sections 3.2.3 and 3.3, a comparison of compounds identified between GC×GC/TOFMS and GC/HRMS was carried out on fraction A. The results for fraction B were similar to those of fraction A. As previously mentioned, 129 and 92 heteroatom-containing compounds were identified using GC×GC/TOFMS and GC/HRMS, respectively, under optimum conditions. The lists of compounds identified using both methods were compared: a total of 57 compounds were commonly identified using both techniques. Table S2 highlights these 57 compounds. When comparing the results obtained for GC×GC/TOFMS and GC/HRMS data processing, false positives and false negatives could be defined. In fact, the false positives corresponded to compounds identified and proposed by the software, but which did not actually exist. These false positives mainly appeared in GC/HRMS, as the proposed heteroatomic compounds were based solely on the measured mass spectrum rather than being controlled by the retention index value and, consequently, the presence of these compounds was not confirmed in the 2D chromatogram. For instance, the oxygenated heteroatomic compound, furan 2-ethyl-5-methyl, was detected twice by GC/HRMS, but only once in the 2D chromatogram of fraction A (with an associated retention time that is in accordance with its retention index). False negatives are compounds that were present but not identified. They are typically overlooked either due to the parameter values applied during data processing, or simply because of limited separation capacity caused by the complexity of the samples and the chromatographic configurations. False positives and false negatives were not considered as responses in the CCD, as their identification necessitated the comparison between GC×GC and GC/HRMS results.

No false positives were observed with GC×GC/TOFMS, while 26 false positives were detected with GC/HRMS. This observation suggested that the enhanced chromatographic separation provided by GC×GC/TOFMS was more effective in eliminating false identifications. However, false negatives were observed in both cases: 5 compounds were missed by GC×GC/TOFMS and 72 by GC/HRMS. Tables S3 and S4 show these compounds. The large number of false negatives in GC/HRMS can be attributed to the absence of a second dimension of separation by polarity, which limits the resolution of coeluted species. Conversely, the 5 compounds detected by GC/HRMS but not by GC×GC/TOFMS exhibited very low intensity. These compounds were not initially detected by the GC×GC/TOFMS software but after the manual inspection of the 2D chromatogram, they were detected and identified.

These results demonstrate that the two analytical methods are complementary, allowing cross-verification of compound identifications. In particular, the results obtained with GC×GC/TOFMS could be considered to confirm the presence or absence of specific compounds detected in GC/HRMS.

For example, hexanal, 3-methyl- was detected in fraction A by GC/HRMS. Surprisingly, this compound was not observed in the upper region of the GC×GC chromatogram, a zone where polar compounds are usually eluted in the second dimension under our conditions of GC×GC analysis. This suggests that this chemical species was probably a false positive when detected by GC/HRMS and processed with the software Trace Finder. On the other hand, GC/HRMS detected quinoline 2,4-dimethyl, a compound that is present in very small amounts and is polar, but ChromaTOF's data processing did not detect it after GC×GC/TOFMS analysis. Following a manual inspection of processed data, the presence of quinoline 2,4-dimethyl was indeed confirmed in the 2D chromatogram. These two examples highlight how GC×GC/TOFMS can help eliminate false positives detected during the GC/HRMS analysis and how the combined use of both analytical approaches enhanced the performance of heteroatomic compound identification in pyrolysis oil fractions.
**4. Conclusions and Perspectives**

The GC×GC/TOFMS and GC/HRMS data processing methods have been developed to analyze distillation fractions of pyrolysis oil from waste tires. In particular, heteroatom-containing compounds were investigated by combining these protocols. Optimization of parameters was first carried out with an experimental design (CCD), yielding a value set for both MSC and MSN parameters of 4 and 60, respectively, to obtain the maximum true positives and the minimum detected heteroatomic peaks during GC×GC/TOFMS data processing. In the case of GC/HRMS data, *m/z* S/N, TIC intensity, and S/N were optimized, with values of 3, 10⁴, and 3, respectively, chosen for *m/z* S/N, TIC intensity, and S/N in the case of GC/HRMS. The curvature of the response surfaces justified the use of the CCD experimental design.

GC×GC/TOFMS proved to be more effective for the detailed identification of heteroatom-containing compounds in light distillation fractions, owing to the second dimension separation based on polarity. In total, 129 heteroatomic compounds were identified in distillation fraction A and 160 in fraction B using GC×GC/TOFMS. However, the results could still be complemented by GC/HRMS, which enabled the additional detection of five trace compounds in fraction A based on the high mass accuracy (within 5 ppm) provided by the Orbitrap analyzer. The optimized parameter set was validated on another distillation fraction (boiling range: 142–169 °C), where 159 true positives were recovered out of 160 confirmed heteroatomic compounds (confirmed with the similarity of the mass spectrum and retention index), with only 2968 detected compounds—compared to 7020 detected compounds using the initial parameter set (MSC = 3, MSN = 5, level (−1, −1) in CCD).

In this study, only two light distillation fractions of pyrolysis oil from waste tires were considered. GC×GC/TOFMS allows for a detailed characterization of heteroatom-containing compounds on these two fractions. For further study, other types of distillation fractions from different types of waste tires should be considered to confirm the optimal values obtained from the CCD. However, the results obtained with these two fractions show that applying a CCD to rationalize key parameters is appropriate for optimizing the experimental parameters associated with specific analytical techniques, enabling full utilization of their potential. Furthermore, it would be interesting to consider GC×GC/HRMS for the analysis of distillation fractions of pyrolysis oils, given that the LECO HRMS is 10 to 100 times less sensitive than the LECO BT-4D TOFMS.

**Supplementary Materials:** The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/separations12090239/s1. Figure S1. GC/HRMS CCD results: standardized values of the coefficients of the model for the number of detected peaks (a) and the number of true positives (b); Figure S2. Residual plots of the number of detected heteroatomic compounds (a) and the number of true positives (b) in the case of GC×GC/TOFMS; Table S1. Identification of 160 true positives in fraction B by GC×GC/TOFMS; Table S2. 57 common heteroatomic compounds detected by GCXGC/TOFMS and GC/HRMStitle; Table S3. 5 compounds detected by GC/HRMS but not by GCXGC/TOFMS; Table S4. 72 compounds detected by GCXGC/TOFMS but not by GC/HRMS.

**Author Contributions:** Conceptualization, X.C., M.S., C.R., J.D. and J.V.; methodology, X.C. and J.V.; software, X.C.; validation, X.C., M.S., C.R., J.D. and J.V.; formal analysis, X.C.; investigation, X.C.; resources, J.V., C.R. and M.S.; data curation, X.C.; writing—original draft preparation, X.C.; writing—review and editing, X.C., M.S., C.R., B.G. and J.V.; visualization, X.C.; supervision, M.S., J.D., C.R. and J.V.; project administration, M.S., C.R., J.D. and J.V.; funding acquisition, B.G. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research received no external funding.

**Data Availability Statement:** Data is contained within the article or Supplementary Materials.
**Conflicts of Interest:** Authors Carlos Rincon and Benoît Gadenne were employed by the company Manufacture Françaises des Pneumatiques Michelin. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

## References

1. Tyre Production Worldwide | Tyre Industry of JAPAN | Japan Automobile Tyre Manufacturers Association, Inc. Available online: https://www.jatma.or.jp/english/tyre_industry/tyreproduction.html (accessed on 14 March 2024).

2. Moasas, A.M.; Amin, M.N.; Khan, K.; Ahmad, W.; Al-Hashem, M.N.A.; Deifalla, A.F.; Ahmad, A. A Worldwide Development in the Accumulation of Waste Tires and Its Utilization in Concrete as a Sustainable Construction Material: A Review. *Case Stud. Constr. Mater.* **2022**, *17*, e01677. [CrossRef]

3. Thomas, B.S.; Gupta, R.C. A Comprehensive Review on the Applications of Waste Tire Rubber in Cement Concrete. *Renew. Sustain. Energy Rev.* **2016**, *54*, 1323–1333. [CrossRef]

4. Arabiourrutia, M.; Lopez, G.; Artetxe, M.; Alvarez, J.; Bilbao, J.; Olazar, M. Waste Tyre Valorization by Catalytic Pyrolysis—A Review. *Renew. Sustain. Energy Rev.* **2020**, *129*, 109932. [CrossRef]

5. Available online: www.etrma.org/wp-content/uploads/2021/05/20210520_ETRMA PRES-RELEASE_ELT-2019.pdf (accessed on 9 June 2025).

6. Kumar Singh, R.; Ruj, B.; Jana, A.; Mondal, S.; Jana, B.; Kumar Sadhukhan, A.; Gupta, P. Pyrolysis of Three Different Categories of Automotive Tyre Wastes: Product Yield Analysis and Characterization. *J. Anal. Appl. Pyrolysis* **2018**, *135*, 379–389. [CrossRef]

7. Williams, P.T. Pyrolysis of Waste Tyres: A Review. *Waste Manag.* **2013**, *33*, 1714–1728. [CrossRef]

8. Antoniou, N.; Zabaniotou, A. Features of an Efficient and Environmentally Attractive Used Tyres Pyrolysis with Energy and Material Recovery. *Renew. Sustain. Energy Rev.* **2013**, *20*, 539–558. [CrossRef]

9. Rofiqulislam, M.; Haniu, H.; Rafiqulalambeg, M. Liquid Fuels and Chemicals from Pyrolysis of Motorcycle Tire Waste: Product Yields, Compositions and Related Properties. *Fuel* **2008**, *87*, 3112–3122. [CrossRef]

10. Martínez, J.D.; Puy, N.; Murillo, R.; García, T.; Navarro, M.V.; Mastral, A.M. Waste Tyre Pyrolysis—A Review. *Renew. Sustain. Energy Rev.* **2013**, *23*, 179–213. [CrossRef]

11. Campuzano, F.; Abdul Jameel, A.G.; Zhang, W.; Emwas, A.-H.; Agudelo, A.F.; Martínez, J.D.; Sarathy, S.M. On the Distillation of Waste Tire Pyrolysis Oil: A Structural Characterization of the Derived Fractions. *Fuel* **2021**, *290*, 120041. [CrossRef]

12. Laresgoiti, M.F.; Caballero, B.M.; de Marco, I.; Torres, A.; Cabrero, M.A.; Chomón, M.J. Characterization of the Liquid Products Obtained in Tyre Pyrolysis. *J. Anal. Appl. Pyrolysis* **2004**, *71*, 917–934. [CrossRef]

13. Amarasekara, A.S.; Reyes, C.D.G. Acidic Ionic Liquid Catalyzed Liquefaction of Corn Cobs and Switchgrass in Acetone: Analysis of Bio-Oils Using LC-MS and GC-MS. *J. Anal. Appl. Pyrolysis* **2020**, *145*, 104752. [CrossRef]

14. Joseph, J.; Rasmussen, M.J.; Fecteau, J.P.; Kim, S.; Lee, H.; Tracy, K.A.; Jensen, B.L.; Frederick, B.G.; Stemmler, E.A. Compositional Changes to Low Water Content Bio-Oils during Aging: An NMR, GC/MS, and LC/MS Study. *Energy Fuels* **2016**, *30*, 4825–4840. [CrossRef]

15. Zhang, Y.; Li, S.; Zhang, Q.; Zhao, Y.; Liu, M.; Zhang, D.; Cai, X.; Wang, N.; Wang, W. Structural Characterization and Transformation of Nitrogen Compounds in Waste Tire Pyrolysis Oils. *J. Chromatogr. A* **2023**, *1702*, 464093. [CrossRef] [PubMed]

16. Li, W.; Manheim, J.M.; Fu, Y.; Laaksonen, T.; Kilaz, G.; Kenttämaa, H.I. Comparison of APCI Orbitrap MS and GCxGC/EI TOF MS for the Hydrocarbon Analysis of Heavy Base Oils. *Fuel* **2023**, *343*, 127993. [CrossRef]

17. Dao Thi, H.; Djokic, M.R.; Van Geem, K.M. Detailed Group-Type Characterization of Plastic-Waste Pyrolysis Oils: By Comprehensive Two-Dimensional Gas Chromatography Including Linear, Branched, and Di-Olefins. *Separations* **2021**, *8*, 103. [CrossRef]

18. Mohan, A.; Dutta, S.; Balusamy, S.; Madav, V. Liquid Fuel from Waste Tires: Novel Refining, Advanced Characterization and Utilization in Engines with Ethyl Levulinate as an Additive. *RSC Adv.* **2021**, *11*, 9807–9826. [CrossRef]

19. Burdová, H.; Pilnaj, D.; Kuráň, P. Application of Low-Energy-Capable Electron Ionization with High-Resolution Mass Spectrometer for Characterization of Pyrolysis Oils from Plastics. *J. Chromatogr. A* **2023**, *1711*, 464445. [CrossRef]

20. Klein, R.; Dugay, J.; Vial, J.; Thiébaut, D.; Colombet, G.; Barreteau, D.; Gruntz, G. Hyphenation of Thermodesorption into GC GC-TOFMS for Odorous Molecule Detection in Car Materials: Column Sets and Adaptation of Second Column Dimensions to TD Pressure Constraints. *Separations* **2024**, *11*, 162. [CrossRef]

21. Boudard, E.; Moumane, N.; Dugay, J.; Vial, J.; Thiébaut, D. Body Volatilome Study Strategy for COVID-19 Biomarker Identification Considering Exogenous Parameters. *Separations* **2024**, *11*, 336. [CrossRef]

22. Mahé, L.; Dutriez, T.; Courtiade, M.; Thiébaut, D.; Dulot, H.; Bertoncini, F. Global Approach for the Selection of High Temperature Comprehensive Two-Dimensional Gas Chromatography Experimental Conditions and Quantitative Analysis in Regards to Sulfur-Containing Compounds in Heavy Petroleum Cuts. *J. Chromatogr. A* **2011**, *1218*, 534–544. [CrossRef]

23. Adahchour, M.; Beens, J.; Vreuls, R.J.J.; Brinkman, U.A.T. Recent Developments in Comprehensive Two-Dimensional Gas Chromatography (GC×GC). *TrAC Trends Anal. Chem.* **2006**, *25*, 438–454. [CrossRef]
24. Liu, Y.; D'Agostino, L.A.; Qu, G.; Jiang, G.; Martin, J.W. High-Resolution Mass Spectrometry (HRMS) Methods for Nontarget Discovery and Characterization of Poly- and per-Fluoroalkyl Substances (PFASs) in Environmental and Human Samples. *TrAC Trends Anal. Chem. 2019*, *121*, 115420. [CrossRef]

25. Staš, M.; Chudoba, J.; Kubička, D.; Pospíšil, M. Chemical Characterization of Pyrolysis Bio-Oil: Application of Orbitrap Mass Spectrometry. *Energy Fuels* **2015**, *29*, 3233–3240. [CrossRef]

26. Legeard, T.; Tisse, S.; Vaccaro, M.; Moufarrej, L.; Mignot, M.; Castilla, C.; Schmitz, I.; Portet-Koltalo, F.; Méausoone, C.; Monteil, C.; et al. Nontargeted Screening of Air Samples Using TD-GC-HRMS to Identify Volatile Compounds as Markers of an Industrial Plant Fire in Rouen, France. *Atmos. Pollut. Res.* **2025**, *16*, 102328. [CrossRef]

27. Guillemant, J.; Lacoue-Nègre, M.; Berlioz-Barbier, A.; Albrieux, F.; De Oliveira, L.P.; Joly, J.-F.; Duponchel, L. Towards a New Pseudo-Quantitative Approach to Evaluate the Ionization Response of Nitrogen Compounds in Complex Matrices. *Sci. Rep.* **2021**, *11*, 6417. [CrossRef]

28. Williams, P.T.; Bottrill, P. Sulfur-polycyclic aromatic hydrocarbons in tyre pyrolysis oil. *Fuel* **1995**, *74*, 736–742. [CrossRef]

29. Kucinska-Lipka, J.; Janik, H.; Balas, A. Progress in Used Tyres Management in the European Union: A Review. *Waste Manag.* **2012**, *32*, 1742–1751. [CrossRef]

30. Cerasa, M.; Balducci, C.; Giannelli Moneta, B.; Santoro, S.; Perilli, M.; Nikiforov, V. Critical Insights into Untargeted GC-HRMS Analysis: Exploring Volatile Organic Compounds in Italian Ambient Air. *Separations* **2025**, *12*, 35. [CrossRef]

31. Vial, J.; Jardy, A. Utilisation des plans d'expériences pour évaluer la robustesse d'une méthode d'analyse quantitative par Chromatographie enPhase Liquide. *Analisis 1998*, *26*, 15-24. [CrossRef]

32. Sarazin, C.; Delaunay, N.; Costanza, C.; Eudes, V.; Gareil, P.; Vial, J. On the Use of Response Surface Strategy to Elucidate the Electrophoretic Migration of Carbohydrates and Optimize Their Separation. *J. Sep. Sci.* **2012**, *35*, 1351–1358. [CrossRef]

33. Alinat, E.; Delaunay, N.; Archer, X.; Vial, J.; Gareil, P. Multivariate Optimization of the Denitration Reaction of Nitrocelluloses for Safer Determination of Their Nitrogen Content. *Forensic Sci. Int.* **2015**, *250*, 68–76. [CrossRef]

34. Van Den Dool, H.A.; Kratz, P.D. A Generalization of the Retention Index System Including Linear Temperature Programmed Gas—Liquid Partition Chromatography. *J. Chromatogr. A* **1963**, *11*, 463–471. [CrossRef] [PubMed]

35. Stefanuto, P.-H.; Perrault, K.A.; Dubois, L.M.; L'Homme, B.; Allen, C.; Loughnane, C.; Ochiai, N.; Focant, J.-F. Advanced Method Optimization for Volatile Aroma Profiling of Beer Using Two-Dimensional Gas Chromatography Time-of-Flight Mass Spectrometry. *J. Chromatogr. A* **2017**, *1507*, 45–52. [CrossRef] [PubMed]

36. Bean, H.D.; Hill, J.E.; Dimandja, J.-M.D. Improving the Quality of Biomarker Candidates in Untargeted Metabolomics via Peak Table-Based Alignment of Comprehensive Two-Dimensional Gas Chromatography–Mass Spectrometry Data. *J. Chromatogr. A* **2015**, *1394*, 111–117. [CrossRef] [PubMed]

37. Furno, L.; Combès, A.; Thiébaut, D.; Méré, A.; Passade-Boupat, N.; Vial, J. Liquid Chromatography Column Screening for the Analysis of Corrosion Inhibitor Molecules Using Derringer Desirability Functions. *J. Sep. Sci.* **2024**, *47*, e70046. [CrossRef]

38. Cuzuel, V.; Sizun, A.; Cognon, G.; Rivals, I.; Heulard, F.; Thiébaut, D.; Vial, J. Human Odor and Forensics. Optimization of a Comprehensive Two-Dimensional Gas Chromatography Method Based on Orthogonality: How Not to Choose between Criteria. *J. Chromatogr. A* **2018**, *1536*, 58–66. [CrossRef] [PubMed]

**Disclaimer/Publisher's Note:** The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.