# Towards a Comprehensive Characterization of the Low-Temperature Autoxidation of Di-n-Butyl Ether

Nesrine Belhadj, Maxence Lailliau, Roland Benoit, Philippe Dagaut

► To cite this version:

Nesrine Belhadj, Maxence Lailliau, Roland Benoit, Philippe Dagaut. Towards a Comprehensive Characterization of the Low-Temperature Autoxidation of Di-n-Butyl Ether. Molecules, 2021, 26 (23), pp.7174. 10.3390/molecules26237174. hal-03500230

HAL Id: hal-03500230

https://hal.science/hal-03500230

Submitted on 22 Dec 2021

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

# Towards a Comprehensive Characterization of the Low-Temperature Autoxidation of Di-n-Butyl Ether

Nesrine Belhadj <sup>1,2 ID</sup>, Maxence Lailliau <sup>1,2 ID</sup>, Roland Benoit <sup>1 ID</sup> and Philippe Dagaut <sup>1,* ID</sup>

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

check for
updates

**Citation:** Belhadj, N.; Lailliau, M.; Benoit, R.; Dagaut, P. Towards a Comprehensive Characterization of the Low-Temperature Autoxidation of Di-n-Butyl Ether. *Molecules* **2021**, *26*, 7174. https://doi.org/10.3390/molecules26237174

Academic Editor: Igor Jerković

Received: 5 November 2021

Accepted: 24 November 2021

Published: 26 November 2021

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

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

**Copyright:** © 2021 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/).

**Abstract:** In the present study, we investigated the oxidation of 2500 ppm of di-n-butyl ether under fuel-rich conditions ($\varphi = 2$) at low temperatures (460–780 K), a residence time of 1 s, and 10 atm. The experiments were carried out in a fused silica jet-stirred reactor. Oxidation products were identified and quantified in gas samples by gas chromatography and Fourier transform infrared spectrometry. Samples were also trapped through bubbling in cool acetonitrile for high-pressure liquid chromatography (HPLC) analyses. 2,4-dinitro-phenylhydrazine was used to derivatize carbonyl products and distinguish them from other isomers. HPLC coupled to high resolution mass spectrometry (Orbitrap Q-Exactive®) allowed for the detection of oxygenated species never observed before, i.e., low-temperature oxidation products ($C_8H_{12}O_{4,6}$, $C_8H_{16}O_{3,5,7}$, and $C_8H_{18}O_{2,5}$) and species that are more specific products of atmospheric oxidation, i.e., $C_{16}H_{34}O_4$, $C_{11}H_{24}O_3$, $C_{11}H_{22}O_3$, and $C_{10}H_{22}O_3$. Flow injection analyses indicated the presence of high molecular weight oxygenated products ($m/z > 550$). These results highlight the strong similitude in terms of classes of oxidation products of combustion and atmospheric oxidation, and through autoxidation processes. A kinetic modeling of the present experiments indicated some discrepancies with the present data.

**Keywords:** di-n-butyl ether; jet-stirred reactor; low-temperature oxidation; atmospheric chemistry; high resolution mass spectrometry; liquid chromatography; kinetic modeling

## 1. Introduction

The use of fossil fuels contributes to anthropogenic emissions of CO₂ and the formation of particles in the atmosphere, which are harmful to the ecosystem and cause health problems [1]. In addition, the increase in global energy needs and the depletion of fossil resources create socio-economic problems [2] which lead to the exploration of other energy sources such as nuclear, hydro, wind, solar, and renewable fuels from biomass which do not compete with food resources [3]. The chemical composition of a fuel induces specific physico-chemical characteristics, e.g., viscosity, density, boiling point, etc., as well as combustion and oxidation specificities. These characters influence the performance of internal combustion engines and the formation of pollutants and particles [4]. The development of future oxygenated lignocellulosic biofuels has put di-n-butyl ether (DBE) at the ranking of the most promising lignocellulosic biofuel and fuel additive [5]. Various experimental studies have highlighted the combustion characteristics of di-n-butyl ether and its effect on the performance of engines and injection systems [6–8], particle emissions [9], soot formation [10,11], and flame characteristics [12]. The combustion of di-n-butyl ether has been explored in different experimental systems and over various temperature and pressure ranges. Fan et al. [13] investigated its pyrolysis kinetics using a flow reactor experiment at low and atmospheric pressures. They identified and quantified several intermediates and products using photoionization molecular-beam mass spectrometry

<footnote n="1">¹ Centre National de la Recherche Scientifique–Institut National des Sciences de l’Ingénierie et des Systèmes, Institut de Combustion, Aérothermique, Réactivité, Environnement, 1C Avenue de la Recherche Scientifique, CEDEX 2, 45071 Orléans, France; nesrine.belhadj@cnrs-orleans.fr (N.B.); maxence.lailliau@cnrs-orleans.fr (M.L.); roland.benoit@cnrs-orleans.fr (R.B.)</footnote>

<footnote n="2">² Collegium Sciences et Techniques, Rue de Chartres, Université d'Orléans, 45100 Orléans, France</footnote>

<footnote>* Correspondence: philippe.dagaut@cnrs-orleans.fr</footnote>
(PI-MBMS). Thion et al. [14] studied the oxidation of 1000 ppm of di-n-butyl ether at low to high temperatures (470 to 1250 K) in a fused silica jet-stirred reactor (JSR) at 1 and 10 atm. Oxidized species have been identified by gas chromatography-mass spectrometry (GC-MS) and Fourier transform infrared spectrometry (FTIR). GC analyses have also been performed using mass spectrometry, flame ionization, and thermal conductivity detectors. The authors observed two negative temperature coefficient regions (NTC). They also proposed a detailed chemical kinetic reaction mechanism to simulate their JSR experiments and model literature data. Tran et al. [15] studied the oxidation of di-n-butyl ether at 400–1100 K and the nearly atmospheric pressure using different experimental systems, i.e., a plug flow reactor (PFR) combined with electron ionization molecular-beam mass spectrometry (EI-MBMS) as well as two different jet-stirred reactors with either online gas chromatography or synchrotron vacuum ultraviolet photoionization molecular-beam mass spectrometry (SVUV-PI-MBMS). They confirmed two negative temperature coefficient zones around 500–550 K and 650–750 K, as reported earlier by Thion et al. [14]. About 40 species have been reported, among which 20 low-temperature products were not previously observed, such as hydroperoxides (CH₄O₂, C₃H₈O₂, and C₈H₁₈O₃), ketohydroperoxides (C₈H₁₆O₄), and some highly oxygenated molecules (C₈H₁₄O₅ and C₈H₁₆O₆,₈) [15]. They also proposed a detailed chemical kinetic reaction mechanism to simulate their experiments [15].

More recently, an experimental study has been carried out in our research team [16]. A jet-stirred reactor has been used to study the cool flame oxidation of 5000 ppm of di-n-butyl ether (i.e., 520 K and 480 to 670 K at 1 and 10 atm, respectively). A rapid compression machine (RCM) has also been used to perform ignition experiments of DBE/an air mixture at 550–630 K and 5 bar. High resolution mass spectrometry (HRMS) analyses of JSR and RCM chemical samples have revealed the formation of hydroperoxides (C₈H₁₈O₃), ketohydroperoxides (C₈H₁₆O₄), and a range of highly oxygenated molecules (HOMs), namely C₈H₁₄O₅,₇ and C₈H₁₆O₆,₈. However, that study focused on what can be called combustion chemistry, i.e., fuel + X• (•OH, H•, O•, HO₂•, O₂...) → R• + XH; R• + O₂ → ROO• → •QOOH (H-shift); •QOOH + O₂ → •OOQOOH → HOOQ'=O + •OH, and further H-shifts and O₂ additions. The comparison of the JSR and RCM results indicated strong similitude in terms of low-temperature products detected using these two experimental devices. However, that study, as well as other previous studies, overlooked the formation of other oxygenated products via oxidation channels considered more representative of atmospheric chemistry, i.e., 2 ROO• → ROOOOR → 2 RO• + O₂; RO• → •QOH (H-shift); •QOH + O₂ → •OOQOH → OQ'OH + •OH, and further H-shifts and O₂ additions; ROO• + R'OO• → ROOR' + O₂; ROO• + R'OO• → ROH + R'(-H)=O + O₂. Thus, a more comprehensive study is necessary to clarify the oxidation processes occurring during the cool flame oxidation of di-n-butyl ether and to assess the occurrence of pathways other than extended oxidation routes observed earlier [17–19]. Moreover, whereas several studies have shown a similarity between combustion processes and those involved in atmospheric chemistry, this has not been addressed for DBE.

The aim of the present work is to better characterize the low-temperature oxidation chemistry of DBE and to probe the possible formation of products of importance in an atmospheric chemistry such as highly oxidized molecules (HOMs), products of RO• isomerization and oxidation, and accretion products. To this end, the oxidation of 2500 ppm of di-n-butyl ether was performed at low temperature (460–780 K), a residence time of 1 s, and high pressure (10 atm) in a fused silica jet-stirred reactor. Products of DBE oxidation through both ‘combustion’ and ‘atmospheric’ oxidation routes were tracked using gas chromatography, Fourier transform infrared spectrometry, and ultra-high pressure/ high-pressure liquid chromatography (UHPLC/HPLC) coupled to HRMS. An existing kinetic reaction mechanism [14] was used to model the present experiments.

## 2. Results and Discussion

Di-n-butyl ether low-temperature (460–780 K) and high-pressure (10 atm) oxidation leads to the formation of several molecular intermediates and products, including low mass
chemical species, i.e., H₂, CO, CO₂, H₂O, alkanes (methane, ethane, and propane), olefines (ethylene, propene, 1,3-butadiene, and 1-butene), carbonyls (formaldehyde, acetaldehyde, propionaldehyde, butyraldehyde, acetone, acrolein, and 2-butanone), carboxylic acids (formic acid, acetic acid, and butyric acid), alcohols (methanol, ethanol, and 1-butanol), and low mass cyclic ethers (oxirane and 2,3-DHF). All these compounds were identified by FTIR, gas chromatography (GC), or liquid chromatography (LC-HRMS) using standards. Their profiles were plotted and compared to modeling results using a literature kinetic reaction scheme [14] (see Supplementary Materials, Figures S1–S3).

The present study focuses on characterizing low-temperature intermediates of oxidation such as DBE-yl-hydroperoxides (C₈H₁₈O₃, 5), unsaturated DBE-yl-hydroperoxides (C₈H₁₆O₃, 5, 7), ketohydroperoxides (C₃H₆O₃, C₄H₈O₃, C₅H₁₀O₄, C₆H₁₂O₄, and C₈H₁₆O₄), diones (C₈H₁₄O₃), cyclic ethers and butoxy-carbonyls (C₈H₁₆O₂), and highly oxygenated molecules (C₈H₁₂O₄, 6, C₈H₁₄O₅, 7, and C₈H₁₆O₆, 8).

Figure 1 presents the variation of the mole fraction of di-n-butyl ether as a function of temperature. As can be seen from that figure, two NTCs were observed under the present conditions. This is in line with what was reported earlier [14,15]. We observed maximal fuel conversion at ~540 and ~660 K. The kinetic model tends to represent the GC data fairly well, whereas those from Orbitrap analyses were scaled to the GC mole fraction at the lowest oxidation temperature common to the two types of experiments (480 K), i.e., GC, FTIR, and HPLC-MS.

<figure><img src="image_4.png" /><figcaption>**Figure 1.** DBE mole fraction profile during its low-temperature (460–780 K) and high-pressure (10 atm) oxidation in a JSR (2500 ppm of DBE, φ = 2, and residence time of 1 s). The line corresponds to modeling results and symbols to both GC and UHPLC experimental results (UHPLC-positive APCI-HRMS, C₈H₁₉O⁺, m/z 131.1430). Simulations were performed with the model of Thion et al. [14].</figcaption></figure>

The normalized rate of consumption analysis was performed with PSR [20] at 540 K and 660 K. It indicated that DBE is essentially consumed by H-atom abstraction by hydroxyl radicals at both temperatures. However, the formation of hydroxyl radicals proceeds via different routes when temperature increases from 540 to 660 K. In the first cool flame, the production of hydroxyl radicals is predominantly due to the decomposition of C₈·OOQOOH radicals deriving from the fuel (54.3%) and, to some extent, to C₈ KHPs' decomposition (29.1%). Minor contributions of the decomposition of C₃H₇OOH (1.6%) and C₄·OOQOOH (3.7%) were predicted. At 660K, 25% of the hydroxyl radical production was due to the decomposition of C₁–C₃ ROOH: C₃H₇OOH (12%), C₂H₅OOH (7.1%), and CH₃OOH (5.8%). At 660 K, the decomposition of C₈·OOQOOH radicals and C₈ KHPs represents only 13.1% and 11.7%, respectively. A small contribution to ·OH formation also comes from the successive decomposition of C₄ and C₃ intermediates: C₄·OOQOOH → ·OH + KHP (3%) → ·OH + carbonyls (2.6%) and C₃ KHPs → carbonyls + ·OH (4.4%).
This behavior is responsible for the observed two-stage cool flame oxidation of DBE as initially reported by Thion et al. [14], notably similar to what was reported recently for the cool flame oxidation of di-n-propyl ether in a JSR at 700 Torr [21].

## 2.1. Oxygenated Intermediates Formed in Combustion via the RO₂• Route

### 2.1.1. Hydroperoxides (C₈H₁₈O₃), Dihydroperoxides (C₈H₁₈O₅), Unsaturated Hydroperoxides (C₈H₁₆O₃), Unsaturated Dihydroperoxides (C₈H₁₆O₅), and Unsaturated Trihydroperoxides (C₈H₁₆O₇)

The peroxidation of the fuel's radicals (Reaction 1) and H-atom abstraction by RO₂• (Reactions 2 and 3) can yield hydroperoxides (ROOH), i.e., C₈H₁₈O₃ in the case of di-n-butyl ether.

$$R^{\bullet} + O_2 \rightleftharpoons RO_2^{\bullet}, \text{ i.e., } \cdot C_8H_{17}O + O_2 \rightleftharpoons \cdot OOC_8H_{17}O \quad (1)$$

$$RO_2^{\bullet} + R'H \rightarrow ROOH + R'\cdot, \text{ i.e., } \cdot OOC_8H_{17}O + R'H \rightarrow C_8H_{18}O_3 + R'\cdot \quad (2)$$

$$RO_2^{\bullet} + HO_2^{\bullet} \rightarrow ROOH + O_2, \text{ i.e., } \cdot OOC_8H_{17}O + HO_2^{\bullet} \rightarrow C_8H_{18}O_3 + O_2 \quad (3)$$

Further reactions can yield dihydroperoxides (di-ROOH), i.e., C₈H₁₈O₅, as shown in Scheme 1.

<figure><img src="image_5.png" /><figcaption>**Scheme 1.** Proposed formation pathways for hydroperoxides (C₈H₁₈O₃) and dihydroperoxides (C₈H₁₈O₅) during the low-temperature oxidation of di-n-butyl ether.</figcaption></figure>

Due to the use of high-resolution mass spectrometry, signals corresponding to C₈H₁₈O₃ and C₈H₁₈O₅ were detected using a C₁₈ UHPLC column for products' separation and ionization by positive and negative APCI (C₈H₁₉O₃⁺, *m/z* 163.1327 and C₈H₁₇O₅⁻, *m/z* 193.1081). Profiles are plotted in Figure 2. One can note that both C₈H₁₈O₃ and C₈H₁₈O₅ hydroperoxides reach a maximum at 540 K, which corresponds to the maximum amplitude of the first DBE cool flame. A second moderate increase of concentration (signal) was observed at ~640 K, which is close to the maximum amplitude of the second DBE cool flame (~660 K).

C₈H₁₈O₃ hydroperoxides isomers can undergo H-atom abstraction followed by O₂ addition, H-shift, and loss of a HO₂ radical to form unsaturated hydroperoxides (C₈H₁₆O₃; Scheme 2a-e). The degradation of these unsaturated hydroperoxides (C₈H₁₆O₃) can produce carbonyl species (e.g., 2-butenal, acrolein, acetaldehyde, and formaldehyde). These products were identified here. Mole fraction profiles and simulations for acrolein, acetaldehyde, and formaldehyde are presented in the Supplementary Materials (Figures S2 and S3). Additionally, 2-butenal was identified in di-n-butyl ether oxidation samples through the formation and detection of the DNPH derivative (crotonaldehyde + DNPH) standard.
<figure><img src="image_6.png" /><figcaption>**Figure 2.** Experimental UHPLC HRMS signal profiles for ROOH, C<sub>8</sub>H<sub>18</sub>O<sub>3</sub> (positive APCI, C<sub>8</sub>H<sub>19</sub>O<sub>3</sub><sup>+</sup>), and di-ROOH, C<sub>8</sub>H<sub>18</sub>O<sub>5</sub> (negative APCI, C<sub>8</sub>H<sub>17</sub>O<sub>5</sub><sup>-</sup>), formed during the low-temperature oxidation of di-n-butyl ether.</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>**Scheme 2.** Proposed formation pathways for unsaturated hydroperoxides (C<sub>8</sub>H<sub>16</sub>O<sub>3</sub>) during the low-temperature oxidation of di-n-butyl ether (a–e).</figcaption></figure>

Unsaturated di-hydroperoxides (C<sub>8</sub>H<sub>16</sub>O<sub>5</sub>) and unsaturated tri-hydroperoxides (C<sub>8</sub>H<sub>16</sub>O<sub>7</sub>) can also be formed during the low-temperature oxidation of di-n-butyl ether. Different formation routes can be proposed, as illustrated in Scheme 3a–c.
<figure><img src="image_8.png" /><figcaption>**Scheme 3.** Proposed formation routes for (a,b) unsaturated di-hydroperoxides ($C_8H_{16}O_5$) and (c) unsaturated tri-hydroperoxides ($C_8H_{16}O_7$) during the low-temperature oxidation of di-n-butyl ether.</figcaption></figure>

$C_{18}$-UHPLC HRMS analyses revealed the presence of $C_8H_{16}O_3$, $C_8H_{16}O_5$, and $C_8H_{16}O_7$ species, which can correspond to unsaturated hydroperoxides ($C_8H_{16}O_3$), unsaturated di-hydroperoxides ($C_8H_{16}O_5$), and unsaturated tri-hydroperoxides ($C_8H_{16}O_7$). However, these chemical formulae can correspond to other molecules formed via atmospheric oxidation routes. Thus, we will discuss these experimental results later in the article (Section 2.2).

### 2.1.2. Formation of Ketohydroperoxides ($C_3H_6O_3$, $C_4H_8O_3$, $C_5H_{10}O_4$, $C_6H_{12}O_4$, and $C_8H_{16}O_4$), Diketones ($C_8H_{14}O_3$), Cyclic Ethers ($C_8H_{16}O_2$), and Highly Oxygenated Molecules ($C_8H_{12}O_{4.6}$, $C_8H_{14}O_{5.7}$, and $C_8H_{16}O_{6.8}$)

As mentioned in the introduction, the combustion of fuels is based on the reaction sequences initiated by the abstraction of a hydrogen atom from the fuel. In the case of di-n-butyl ether, the oxidation proceeds via the initial formation of a •$C_8H_{17}O$ radical: $C_8H_{18}O + X^• (•OH, H^•, O^•, HO_2^•, ... ) → •C_8H_{17}O + HX$, followed by a molecular oxygen addition and H-shift → •OOC$C_8H_{17}O → HOO^•C_8H_{16}O$, as well as a second $O_2$ addition and H-shift → HOOC$C_8H_{16}O(OO^•)→ HOO^•C_8H_{15}O(OOH)$ which can lose •OH and form $O=C_8H_{15}O(OOH)$ corresponding to keto-hydroperoxides ($C_8H_{16}O_4$). If the H-shift involves a non-oxygenated CH group, a third $O_2$ addition followed by a H-shift to form (HOO)$_2$•$C_8H_{14}O(OOH)$ can occur, forming the first highly oxygenated molecule $O=C_8H_{14}O(OOH)_2$ known as keto-di-hydroperoxides ($C_8H_{16}O_6$), or the process can undergo a fourth $O_2$ addition followed by a H-atom transfer to form (HOO)$_3$•$C_8H_{13}O(OOH)$, which can decompose and yield •OH + keto-tri-hydroperoxides ($C_8H_{16}O_8$): (O=)•$C_8H_{13}O(OOH)_3$. The multi-step autoxidation products (keto-hydroperoxides $C_8H_{16}O_4$, keto-dihydroperoxides $C_8H_{16}O_6$, and keto-tri-hydroperoxides $C_8H_{16}O_8$) can undergo H-abstraction, lose OH from the OOH group, and form species with one or more carbonyl groups [22]:
$C_8H_{16}O_4$ (keto-hydroperoxides) + X → HX + OH + $C_8H_{14}O_3$ (diketones)

$C_8H_{16}O_6$ (keto-di-hydroperoxides) + X → HX + OH + $C_8H_{14}O_5$ (di-keto-hydroperoxides)

$C_8H_{14}O_5$ (di-keto-hydroperoxides) + X → HX + OH + $C_8H_{12}O_4$ (tri-ketones)

$C_8H_{16}O_8$ (keto-tri-hydroperoxides) + X → HX + OH + $C_8H_{14}O_7$ (di-keto-di-hydroperoxides)

$C_8H_{14}O_7$ (di-keto-di-hydroperoxides) + X → HX + OH + $C_8H_{12}O_6$ (tri-keto-hydroperoxides)

Proposed pathways leading to the formation of $C_3, C_4, C_5$, and $C_6$ KHPs are presented in Scheme 4. These ketohydroperoxides ($C_3H_6O_3, C_4H_8O_3, C_5H_{10}O_4$, and $C_6H_{12}O_4$) can be formed from different low-temperature oxidation intermediates, i.e., from the decomposition of $C_8$-ketohydroperoxides (Scheme 4a-c) and $C_8$-peroxy radicals' reactions (Scheme 4d). Their decomposition produces hydroxyl radicals contributing to the observed second cool flame (620–720 K), as shown in Figure 1.

<figure><img src="image_9.png" /><figcaption>**Scheme 4.** Proposed formation routes for $C_{3,4,5,6}$ ketohydroperoxides (a–d).</figcaption></figure>

Due to the use of liquid chromatography (HPLC and UHPLC) coupled to high resolution mass spectrometry and positive or negative APCI modes, signals corresponding to $C_{3,4,5,6,8}$ ketohydroperoxides $C_3H_6O_3$ ($C_3H_5O_3^-$, m/z 89.0243), $C_4H_8O_3$ ($C_4H_9O_3^+$, m/z 105.0543), $C_5H_{10}O_4$ ($C_5H_{11}O_4^+$, m/z 135.0651), $C_6H_{12}O_4$ ($C_6H_{13}O_4^+$, m/z 149.0807), and $C_8H_{16}O_4$ ($C_8H_{17}O_4^+$, m/z 177.1121) were detected (chromatograms are given in Supplementary Materials, Figure S4). Their signal profiles were plotted and compared to the modeling results (Figure 3). One should note that only the formation of $C_3H_6O_3$ and $C_8H_{16}O_4$ was considered in the model of Thion et al. [14]. For convenience, LC-HRMS data for $C_3$ and $C_8$ KHPs were scaled to their maximum simulated mole fractions. For the other KHPs, we present the variation of the signal of their positive ions as a function of temperature. Considering the low signals obtained for $C_5$ and $C_6$ KHPs, the data were multiplied by 100. As can be seen from Figure 3, the kinetic model predicts a maximum mole fraction of $C_8H_{16}O_4$ at about 540 K, which is in agreement with the experimental data. The model predicts the maximum concentration of $C_3H_6O_3$ at 580 K, which is close to what was observed experimentally. For the other KHPs, the maximum signal was recorded at ~580 K ($C_6H_{12}O_4$) and 540 K ($C_5H_{10}O_4$ and $C_4H_8O_3$).
<figure><img src="image_10.png" /><figcaption>**Figure 3.** Comparison of the mole fraction profiles for the C₃H₆O₃, C₄H₈O₃, C₅H₁₀O₄, C₆H₁₂O₄, and C₈H₁₆O₄ KHPs formed during the low-temperature oxidation of 2500 ppm of DBE. The continuous line represents the computed C₈H₁₆O₄ mole fractions, while the dashed line represents the computed C₃H₆O₃ mole fractions. Experimental data are represented by symbols.</figcaption></figure>

As previously shown [22], diones (here, di-keto-ethers) can be produced from ketohydroperoxides, i.e., C<sub>8</sub>H<sub>16</sub>O<sub>4</sub> (keto-hydroperoxides) + X &rarr; HX + OH + C<sub>8</sub>H<sub>14</sub>O<sub>3</sub> (diketones). A signal corresponding to di-keto-ethers (C<sub>8</sub>H<sub>14</sub>O<sub>3</sub>) was detected by C<sub>18</sub>-UHPLC HRMS as well as positive APCI (C<sub>8</sub>H<sub>15</sub>O<sub>3</sub><sup>+</sup>, *m/z* 159.1015). The chromatograms revealed the presence of several C<sub>8</sub>H<sub>14</sub>O<sub>3</sub> isomers (found in the Supplementary Materials, Figure S5, presenting a C<sub>8</sub>H<sub>15</sub>O<sub>3</sub><sup>+</sup> UHPLC-HRMS chromatogram). Only butyric anhydride was identified. The unavailability of diketo-ethers standards made it difficult to identify them. The UHPLC-HRMS experimental profiles for butyric anhydride and the sum of the other diones are presented in Figure 4a.

<figure><img src="image_11.png" /><figcaption>**Figure 4.** Comparison between mole fractions/signal profiles of diones (C<sub>8</sub>H<sub>14</sub>O<sub>3</sub>), cyclic ethers (C<sub>8</sub>H<sub>16</sub>O<sub>2</sub>), and other C<sub>8</sub>H<sub>16</sub>O<sub>2</sub> isomers formed during the low-temperature oxidation of di-n-butyl ether. (a): Comparison between UHPLC-HRMS experimental results for butyric anhydride and the sum of diones signal. (b): Comparison between experimental GC and UHPLC-HRMS data for cyclic ethers (C<sub>8</sub>H<sub>16</sub>O<sub>2</sub>) and simulations (line). (c): Comparison between GC and UHPLC-HRMS experimental mole fractions of butyl butyrate (C<sub>8</sub>H<sub>16</sub>O<sub>2</sub>). (d) UHPLC-HRMS-positive APCI signal profile of R'(-H)=O species corresponding to C<sub>8</sub>H<sub>16</sub>O<sub>2</sub> butoxy-carbonyl isomers.</figcaption></figure>
Cyclic ethers ($C_8H_{16}O_2$) can be formed during di-n-butyl ether oxidation [14,15]. They are produced by the decomposition of $QO_2H$ radicals: •QOOH → cyclic ethers + •OH, although in the case of DBE, one can write $HOO•C_8H_{16}O → C_8H_{16}O_2 + •OH$. These cyclic ethers have the same molecular formula as esters or keto-ethers ($C_8H_{16}O_2$), such as butyl butyrate (n-butyl butanoate), 4-butoxybutanal, 4-butoxybutan-2-one, and 1-butoxybutan-2-one. n-Butyl butanoate can be formed from DBE-yl-hydroperoxides ($C_8H_{18}O_3$) and DBE-yl-hydroperoxy radicals ($C_8H_{17}O_3•$) via H-abstractions from ROOH or via the recombination/disproportionation of $RO_2•$ [15]. 4-Butoxybutanal can be formed through ROO• self-reaction: $ROO• + R'OO• → ROH + R'(-H)=O + O_2$ (cf. Reaction 12 in Section 2.2). Unfortunately, the unavailability of 4-butoxybutanal, 4-butoxybutan-2-one, and 1-butoxybutan-2-one standards made their identification impossible.

$C_8H_{16}O_2$ species formed during the low-temperature oxidation of di-n-butyl ether were detected through $C_{18}$-UHPLC-positive APCI HRMS analyses. The chromatograms revealed the presence of several $C_8H_{16}O_2$ isomers ($C_8H_{17}O_2^+$, m/z 145.1222; see Supplementary Materials, Figure S5). 2,4-DNPH derivatization was applied to distinguish cyclic ethers and butoxy-carbonyl $C_8H_{16}O_2$ isomers (4-butoxybutanal, 4-butoxybutan-2-one, and 1-butoxybutan-2-one). Cyclic ethers do not react with 2,4-DNPH, whereas carbonyls do. Therefore, the intensity of their chromatographic peaks decreases. Computed and experimental (GC, UHPLC) mole fraction profiles of the sum of cyclic ethers and n-butyl butanoate, and the sum of butoxy-carbonyls isomers are presented in Figure 4b–d.

One should note that unlike the modeling and UHPLC-HRMS data, which reveal the formation of several cyclic ethers (ten species in the kinetic model, i.e., 4-ethyl-2-propyl-1,3-dioxolane, 4-methyl-2-propyl-1,3-dioxane, 2-butoxytetrahydrofuran, etc., and at least five in the UHPLC-HRMS analyses), GC analyses allowed for the detection of only two cyclic ethers. Their maxima of concentration are not reached at the same temperature (640 K in GC, 680 K in UHPLC, and 720 K in the simulation), which could be explained by the absence of the detection of some cyclic ethers in the GC analyses. The kinetic model underestimates the formation of cyclic ethers by about 50%.

Concerning butyl butyrate profiles (Figure 4c), we noticed a very good qualitative agreement between the GC and UHPLC results. In addition, it is noticeable that butyl butyrate reached maxima at 520 K and 640 K. This could be due to its formation via different intermediate species.

Figure 4d shows that the butoxy-carbonyl profile (sum of 4-butoxybutanal, 4-butoxybutan-2-one, and 1-butoxybutan-2-one) peaks at 540 K.

One should note that KHPs are isomeric forms of hydroperoxyl-cyclic ethers. We have considered that the chemical formulae detected in this work for these isomers correspond mostly to KHPs. This assumption is based on results of simulations performed for n-hexane under conditions close to those of the present study [23] using a literature mechanism [18], including the formation of functionalized cyclic ethers (hydroperoxy and di-hydroperoxy-cyclic ethers). The model predicted 16 and 160 times less hydroperoxy and di-hydroperoxy-cyclic ethers, respectively, than KHPs. This result seems to validate our assumption. Scheme 5 gives a simplified view of the reaction pathways forming these products.

As presented earlier [16] highly oxygenated molecules ($C_8H_{16}O_6$ and $C_8H_{16}O_8$) resulting from a third and fourth $O_2$ addition on fuel radicals were detected using flow injection (FIA-HRMS) analyses. However, in the present study, we developed a $C_{18}$-UHPLC HRMS analytical method. Chromatographic analyses revealed the presence of several $C_8H_{16}O_6$ and $C_8H_{16}O_8$ peaks (see Supplementary Materials, Figure S6). In addition, the other multi-oxygenated intermediates ($C_8H_{12}O_{4,6}$ and $C_8H_{14}O_{5,7}$), resulting from the decomposition of $C_8H_{16}O_6$, and $C_8H_{16}O_8$, were also detected using UHPLC HRMS with positive and negative APCI. $C_8H_{14}O_5$ (DBE-di-ketohydroperoxides) was observed by Tran et al. [15] and in our earlier study [16], and $C_8H_{14}O_7$ (DBE-di-keto-dihydroperoxides) was also observed in our previous study [16], but $C_8H_{12}O_{4,6}$ has never been observed previously. Figure 5a shows UHPLC-HRMS experimental profiles of DBE-keto-dihydroperoxides ($C_8H_{16}O_6$), DBE-diketo-hydroperoxides ($C_8H_{14}O_5$), and DBE-triketones ($C_8H_{12}O_4$). Figure 5b shows
UHPLC-HRMS experimental profiles of DBE-keto-trihydroperoxides (C₈H₁₆O₈), DBE-di-keto-dihydroperoxides (C₈H₁₄O₇), and DBE-tri-keto-hydroperoxides (C₈H₁₂O₆). The C₈H₁₆O₆, C₈H₁₄O₅, C₈H₁₂O₄, and C₈H₁₆O₈ experimental results were obtained using negative APCI. Since these species are not efficiently ionized in the same APCI mode, the C₈H₁₂O₆ and C₈H₁₄O₇ data were obtained with positive APCI.

<figure><img src="image_12.png" /><figcaption>**Scheme 5.** Schematic representation of the low-temperature oxidation pathways including the formation of KHPs, di and tri-hydroperoxides, and functionalized cyclic ethers isomers. (1) to (4) refer to first to fourth O₂ addition.</figcaption></figure>

<figure><img src="image_13.png" /><figcaption>**Figure 5.** Experimental UHPLC-HRMS signal profiles for multi-oxygenated molecules formed during the low-temperature oxidation of di-n-butyl ether: (a) C₈H₁₆O₆ (negative APCI, C₈H₁₅O₆⁻, m/z 207.0873), C₈H₁₄O₅ (negative APCI, C₈H₁₃O₅⁻, m/z 189.0768), and C₈H₁₂O₄ (negative APCI, C₈H₁₁O₄⁻, m/z 171.0662); and (b) C₈H₁₆O₈ (negative APCI, C₈H₁₅O₈⁻, m/z 239.0771), C₈H₁₄O₇ (positive APCI, C₈H₁₅O₇⁺, m/z 223.0812), and C₈H₁₂O₆ (positive APCI, C₈H₁₃O₆⁺, m/z 205.0706).</figcaption></figure>

## 2.2. Oxygenated Products Formed via Atmospheric Oxidation Routes (RO₂• and RO••)

Additional oxidation channels not taken into account in the combustion were considered in the atmospheric chemistry, apart from the usual oxidation routes (Reactions (4) and (5)).

$$R + X^{\bullet} \rightarrow R^{\bullet} + XH \quad (4)$$

$$R^{\bullet} + O_2 \rightarrow ROO^{\bullet} \quad (5)$$

Among them, one can find reactions which could also occur during combustion but are generally ignored (Reactions (7)–(12)):

$$2 ROO^{\bullet} \rightarrow ROOOOR \rightarrow 2 RO^{\bullet} + O_2 \quad (6)$$

$$RO^{\bullet} \rightarrow \bullet QOH \ (H-shift) \quad (7)$$

$$\bullet QOH + O_2 \rightarrow \bullet OOQOH \rightarrow OQ^{\bullet}OH + \bullet OH \quad (8)$$
$$ \bullet OOQOH \rightarrow HOO\bullet Q'OH \ (\text{H-shift}) \qquad (9) $$

$$ HOO\bullet Q'OH + O_2 \rightarrow \bullet OO(HOO)Q'OH \rightarrow O(HOO)Q''OH + \bullet OH \qquad (10) $$

$$ ROO\bullet + R'OO\bullet \rightarrow ROOR' + O_2 \qquad (11) $$

$$ ROO\bullet + R'OO\bullet \rightarrow ROH + R'_{(-H)} = O + O_2 \qquad (12) $$

Besides, ROH can also be formed through H-atom abstraction by an alkoxy radical:

$$ RO\bullet + XH \rightarrow ROH + X\bullet \qquad (13) $$

Through Reactions (4)–(10), peroxy radicals ($ROO\bullet$ with $R = C_8H_{17}O$) can form $RO\bullet$ radicals (here, $\bullet OC_8H_{17}O$), which undergo successive reactions: H-shift isomerization to form $HO\bullet C_8H_{16}O$ followed by $O_2$ addition and H-shift yielding ($HOO)(HO)\bullet C_8H_{15}O$. That radical can lose OH by $O-O$ scission in $OOH$ to form $C_8H_{16}O_3$ corresponding to $O=C_8H_{15}O(HO)$. Alternatively, it can undergo $O_2$ addition and H-atom transfer to form ($HOO)_2(HO)\bullet C_8H_{14}O$. That radical can either lose an OH radical and form $C_8H_{16}O_5$ corresponding to ($HOO)(HO)(O=)C_8H_{14}O$ or undergo $O_2$ addition and an H-shift to give ($HOO)_3(HO)\bullet C_8H_{13}O$, which can lose $\bullet OH$ and form $C_8H_{16}O_7$ corresponding to ($HOO)_2(HO)(O=)C_8H_{13}O$. One should note that oxygenated molecules ($C_8H_{16}O_{3,5,7}$) formed through atmospheric chemistry pathways present the same chemical formulas as some molecules formed in the fuel's combustion mechanism. These include $C_8H_{16}O_3$, unsaturated DBE-yl-hydroperoxides (Scheme 2a–e), $C_8H_{16}O_5$, unsaturated DBE-yl-di-hydroperoxides, and $C_8H_{16}O_7$, which can be unsaturated DBE-yl-tri-hydroperoxides (Scheme 3a–c). In order to distinguish these two chemical groups, 2,4-DNPH derivatization was performed (see details in Section 3.2). RP-C$_{18}$ UHPLC-HRMS analyses showed the presence of $C_8H_{16}O_{3,5,7}$-dinitrophenylhydrazone derivatives ($C_{14}H_{20}N_4O_{6,8,10}$, respectively). In addition, chromatograms of the $C_8H_{16}O_{3,5,7}$ species obtained before and after DNPH derivatization revealed a decrease of the intensity of some chromatographic peaks, which indicates the presence of carbonyl compounds. The other peaks showing no decreased intensity should be attributed to species having no carbonyl function, e.g., unsaturated hydroperoxides. The chromatograms are given in the Supplementary Materials, Figure S7.

The UHPLC-HRMS analyses of $C_8H_{16}O_{3,5,7}$ (without DNPH derivatization) species were performed using APCI (+/-) modes. However, in Figure 6, only the strongest signals are presented for $C_8H_{16}O_{3,5}$ (APCI+, $C_8H_{17}O_3^+$, m/z 161.1172, and $C_8H_{17}O_5^+$, m/z 193.1070) and $C_8H_{16}O_7$ (APCI-, $C_8H_{15}O_7^-$, m/z 223.0823) isomers. One should note that $C_8H_{16}O_7$ isomers were ionized only in the negative APCI mode. Figure 6 shows UHPLC-HRMS signal profiles of $C_8H_{16}O_{3,5,7}$ species formed in combustion reactions. Based on chromatographic analyses with and without the addition of 2,4-DNPH (see Supplementary Materials, Figure S7), the reported data for $C_8H_{16}O_7$ (Figure 6a) correspond to different co-eluted compounds, at least one of them containing a carbonyl function which reacted with 2,4-DNPH. UHPLC-HRMS data for $C_8H_{16}O_{3,5}$ are presented in Figure 6b. These products can be formed via atmospheric oxidation routes starting with the isomerization of the initial $RO\bullet$ radical into $\bullet QOH$ (H-shift), followed by an $O_2$ addition (Reactions (7)–(10)), i.e., $RO\bullet \rightarrow \bullet QOH$ (H-shift), $\bullet QOH + O_2 \rightarrow \bullet OOQOH \rightarrow OQ'OH + \bullet OH$, $\bullet OOQOH \rightarrow HOO\bullet QOH$ (H-shift), and $HOO\bullet QOH + O_2 \rightarrow \bullet OO(HOO)QOH \rightarrow O(HOO)Q''OH + \bullet OH$.

Let us now focus on the species formed through Reactions (11)–(13), yielding $ROOR'$, ROH, and $R'_{(-H)}=O + ROH$, respectively.

UHPLC-HRMS analyses with positive APCI ionization allowed for the detection of molecules corresponding to $ROOR'$ species. Among them, one could detect $C_{16}H_{34}O_4$ (with $R=R'=C_8H_{17}O$), $C_{11}H_{24}O_3$ (with $R=C_8H_{17}O$ and $R'=C_3H_7$), $C_{11}H_{22}O_3$ (with $R=C_8H_{15}O$ and $R'=C_3H_7$), and $C_{10}H_{22}O_3$ (with $R=C_8H_{17}O$ and $R'=C_2H_5$). The evolution of $ROOR'$ experimental profiles are given in Figure 7. One can see that all $ROOR'$ ($C_{16}H_{34}O_4$, $C_{11}H_{24}O_3$, $C_{11}H_{22}O_3$, and $C_{10}H_{22}O_3$) reach their maximum intensity at 540 K, which corresponds to the maximum intensity of the first DBE cool flame. The relative importance of their peak intensity varies as follows: $C_{16}H_{34}O_4 > C_{11}H_{22}O_3 > C_{11}H_{24}O_3 > C_{10}H_{22}O_3$.
<figure><img src="image_14.png" /><figcaption>**Figure 6.** Experimental UHPLC-HRMS signal profiles of C<sub>8</sub>H<sub>16</sub>O<sub>3,5,7</sub> formed during the low-temperature oxidation of DBE. (a) Species formed in combustion: unsaturated DBE-yl-hydroperoxides (C<sub>8</sub>H<sub>16</sub>O<sub>3</sub>), unsaturated DBE-yl-di-hydroperoxides (C<sub>8</sub>H<sub>16</sub>O<sub>5</sub>), and unsaturated DBE-yl-tri-hydroperoxides (C<sub>8</sub>H<sub>16</sub>O<sub>7</sub>) and carbonyls, i.e., (HOO)<sub>2</sub>(HO)(O=C)<sub>8</sub>H<sub>13</sub>O. (b) Species formed in the atmospheric RO• oxidation route: O=C<sub>8</sub>H<sub>15</sub>O(HO) (C<sub>8</sub>H<sub>16</sub>O<sub>3</sub>) and (HOO)(HO)(O=C)<sub>8</sub>H<sub>14</sub>O (C<sub>8</sub>H<sub>16</sub>O<sub>5</sub>).</figcaption></figure>

<figure><img src="image_15.png" /><figcaption>**Figure 7.** Experimental signal for ROOR' species formed during the low-temperature oxidation of di-n-butyl ether: C<sub>16</sub>H<sub>34</sub>O<sub>4</sub> (C<sub>16</sub>H<sub>35</sub>O<sub>4</sub><sup>+</sup>, m/z 291.2528), C<sub>11</sub>H<sub>24</sub>O<sub>3</sub> (C<sub>11</sub>H<sub>25</sub>O<sub>3</sub><sup>+</sup>, m/z 205.1797), C<sub>11</sub>H<sub>22</sub>O<sub>3</sub> (C<sub>11</sub>H<sub>23</sub>O<sub>3</sub><sup>+</sup>, m/z 191.1641), and C<sub>10</sub>H<sub>22</sub>O<sub>3</sub> (C<sub>10</sub>H<sub>23</sub>O<sub>3</sub><sup>+</sup>, m/z 203.1640). The analyses were performed using UHPLC HRMS and positive APCI.</figcaption></figure>

Both Reactions (12) and (13) can form ROH molecules. In the case of di-n-butyl ether, butoxy-butanol species (C<sub>8</sub>H<sub>18</sub>O<sub>2</sub>) were detected using UHPLC HRMS analyses (positive APCI, C<sub>8</sub>H<sub>19</sub>O<sub>2</sub><sup>+</sup>, m/z 147.1378; Figure 8a). We can note that the C<sub>8</sub>H<sub>18</sub>O<sub>2</sub> signal reached a maximum at 540 K. Based on simulated RO• and ROO• mole fractions (Figure 8b) showing RO• radicals dominating, one can expect that the production of ROH is mainly due to H-atom abstraction reactions by RO• radicals (Reaction (13)). The formation of ROH and R'<sub>-H</sub>=O, involving ROO• radicals, should be a minor channel under the present conditions as it seems to be based on the observed minor formation of butoxy-carbonyls (Figure 8a).
<figure><img src="image_16.png" /><figcaption>**Figure 8.** (a) Experimental profile of ROH (C<sub>8</sub>H<sub>18</sub>O<sub>2</sub>) and the sum of R'(-H)=O (butoxy-ethers) formed during the low-temperature oxidation of di-n-butyl ether. The analyses were performed using UHPLC HRMS and positive APCI. (b) Simulated RO• and ROO• radical mole fractions vs. temperature using the model of Thion et al. [14].</figcaption></figure>

Concerning the formation of an aldehyde R'(-H)=O via Reaction 12, producing 4-butoxybutanal (C<sub>8</sub>H<sub>16</sub>O<sub>2</sub>), DNPH derivatization was used to distinguish between C<sub>8</sub>H<sub>16</sub>O<sub>2</sub> isomers (cyclic ethers and carbonyl-ethers). Experimental UHPLC-HRMS results are presented in Section 2.1.2 (Figure 4d) and Section 2.2 (Figure 8a).

## 2.3. Further Characterization of Oxygenated Products

To further characterize the products of the oxidation of DBE, additional analyses and data processing were performed. A DBE oxidation sample collected at 540 K was analyzed using the flow injection analysis technique (FIA). The analyses were performed using the negative APCI-HRMS parameters given in Table S1 of the Supplementary Materials. As can be seen from the mass spectrum shown in Figure 9, the mass range extended to m/z > 550. This observation is somewhat unexpected for the oxidation of DBE in a cool flame, whereas it would likely be less surprising if DBE had been oxidized under simulated atmospheric conditions where accretion products have been reported for a range of organics [24].

<figure><img src="image_17.png" /><figcaption>**Figure 9.** FIA-APCI(-)-HRMS data obtained with a DBE oxidation sample collected at 540 K.</figcaption></figure>

Unfortunately, no data are available for the atmospheric oxidation of DBE to be compared to the present results. However, in a previous study of ours concerning limonene oxidation in a JSR at atmospheric pressure and 590 K [25], we showed that products of combustion and atmospheric oxidation presented strong similitude in terms of the chemical formulae detected and the variation of the oxidation state of carbon atoms as a function of the number of carbon atoms in the detected formulae. The oxidation state of carbon atoms (OSc ≈ 2O/C-H/C) allows for evaluating the degree of oxidation of a large range of organic species (alcohols, carbonyls, carboxylic acids, esters, and ethers, but not peroxides) [26].
We computed the variation of OSc as a function of the number of carbon atoms, namely nC, for presently detected molecular formulae obtained by HRMS (Figure 10).

<figure><img src="image_18.png" /><figcaption>**Figure 10.** Variation of the OSc as a function of the number of carbon atoms in molecular formulae detected by negative APCI FIA-HRMS of a DBE oxidation sample collected at 540 K.</figcaption></figure>

One can clearly distinguish several regions corresponding to low-volatility oxygenated organic aerosols (LV-OOA), semi-volatile oxygenated organic aerosols (SV-OOA), and water-soluble organic aerosols (WSOC), as defined by Kroll et al. [26]. Since di-n-butyl ether is composed of eight carbon atoms ($C_8H_{18}O$), for nC < 8, the observed products result from oxidative fragmentation. For 8 < nC < 16, we observed products of addition or condensation, and for nC > 16, we observed products of oligomerization. A comparison of Figure 10 with results published earlier for the oxidation of limonene in a JSR and with a range of atmospheric oxidation studies [25] indicates unexpected similar trends regarding the presence of LV-OOA, SV-OOA, and WSOC for the two fuels, namely DBE and limonene (see Supplementary Material, Figure S8).

# 3. Materials and Methods

## 3.1. JSR Experiments

The oxidation of DBE was performed in a jet-stirred reactor, which has 4 injectors (nozzles of 1 mm I.D.) providing stirring [27–29]. Di-n-butyl ether (>99% purity, Sigma-Aldrich®, Saint-Louis, MO, USA) was atomized by a nitrogen ($N_2$) gas flow and vaporized in a heated chamber. Diluted oxygen and a fuel–nitrogen mixture were sent separately to the injectors to avoid oxidation before reaching the JSR. The liquid fuel was pumped using a Shimadzu® LC10 AD-VP HPLC pump equipped with a Shimadzu® DGU-20A3 online degasser. $N_2$ and $O_2$ flow rates were controlled by mass flow controllers (Brooks Instruments, Hatfield, PA, USA). A Pt-Pt/Rh—10% (0.1 mm O.D., located inside a thin wall jacket) was used to verify good thermal homogeneity (gradients of <1 K/cm) along the vertical axis of the JSR. The experiments were performed twice. First, we collected gas samples for online FTIR and gas phase analysis (GC). Secondly, we collected samples by bubbling gas from the sampling sonic probe into cooled acetonitrile (0 °C). The liquid samples were stored in a freezer (−15 °C) and analyzed by both HPLC and UHPLC coupled to high resolution mass spectrometry (Orbitrap Q-Exactive®, Thermo Fisher Scientific, Waltham, MA, USA). This analytical protocol has been used in previous works of ours where its effectiveness was demonstrated [16,23]. The JSR experimental conditions
are shown in Table 1. Excellent stability was observed for the fuel, oxygen, and nitrogen flowrates, temperature, and pressure in the JSR during the experiments.

<table><caption>**Table 1.** Experimental conditions for the DBE oxidation in a JSR.</caption><thead><tr><th>Parameters</th><th>Values</th></tr></thead><tbody><tr><td>Reactive mixture</td><td>0.25% of DBE, 1.5% of O<sub>2</sub>, 98.25% of N<sub>2</sub></td></tr><tr><td>Equivalence ratio (φ)</td><td>2</td></tr><tr><td>Residence time (s)</td><td>1</td></tr><tr><td>Pressure (atm)</td><td>10</td></tr><tr><td>Temperature (K)</td><td>460–780</td></tr><tr><td>Volume of acetonitrile (mL)</td><td>20</td></tr><tr><td>Bubbling time (min)</td><td>60</td></tr></tbody></table>

## 3.2. Analytical Procedures

Different complementary analytical methods were developed to analyze di-n-butyl ether oxidation products. Fourier transform infrared spectroscopy was used to detect formic and acetic acids, namely H₂O, CO₂, CO, and CH₂O. Gas chromatography (GC) with a flame ionization detector (FID) to quantify hydrocarbons and oxygenates up to C₈ was employed. A GC with a thermal conductivity detector (TCD) was employed to quantify H₂ and O₂, and a GC-quadrupole mass spectrometer (q-MS) was used for ≤C₈ products' identification.

Liquid samples were analyzed by high-pressure and ultra-high-pressure liquid chromatography (HPLC/UHPLC) coupled to an Orbitrap Q-Exactive® mass spectrometer, which has very high resolution (140,000), high mass accuracy (<1 ppm), and a large *m/z* range (50–1800). However, due to the limitation of this mass spectrometer to detect molecules with *m/z* < 50, e.g., CO₂, CH₄, C₂H₄, C₃H₆, C₂H₄O, CH₂O, C₂H₄O, etc., these molecules were characterized either by FTIR or GC, or after the chemical derivatization of carbonyls with 2,4-dinitrophenyl hydrazine (moistened with water; >99% purity from Sigma-Aldrich®) Then, *m/z* increased and derivation products could be detected (carbonyl + 2,4-DNPH = H₂O + 2,4-DNPHydrazone). To this end, 50 µL of 2,4-DNPH saturated solution (dissolution of 2 g in 100 mL of ACN) and 20 µL of diluted phosphoric acid from Sigma-Aldrich® (100 µL of H₃PO₄ in 1 mL of ACN) were added to 1 mL of di-n-butyl ether oxidation samples. The 2,4-DNPH derivatization reaction was carried out at ambient temperatures for 2 h.

Liquid chromatography analyses were mostly performed in the reversed-phase (RP-C₁₈ Luna Omega column, Phenomenex®, Torrance, CA, USA). Additional high-performance liquid chromatography (HPLC)–HRMS analyses were performed with a Supelco® Ascentis Silica column. To this end, a Vanquish™ system from Thermo Fisher Scientific®, Waltham, MA, USA was used. Soft ionization using atmospheric pressure chemical ionization (APCI) in positive ([M+H]⁺) and negative ([M-H]⁻) modes was used to ionize products after chromatographic separation. Direct infusion was also used (see Supplementary Materials, Table S1). Then, the sample was injected into the ionization chamber of the mass spectrometer using a syringe pump. FTIR, GC, and HPLC/UHPLC analytical settings are detailed in the Supplementary Materials, Table S1. The HPLC-Orbitrap data are only qualitative because no calibration could be done in this work. The other data obtained by GC and FTIR are quantitative because the chemical species (reactants and products) were calibrated using standards. Uncertainties on the mole fractions are of the order of 10% in GC and 15% in FTIR. HOMs are minor products (in the ppm-range); the carbon balance could not be impacted by the fact they are not quantified in this work.

## 3.3. Kinetic Modeling

The present experiments were simulated using the PSR computer code [20]. The detailed kinetic reaction mechanism of Thion et al. [14] was used. It included 2768 reactions (most of them reversible) and 467 chemical species describing both low and high-
temperature oxidation of DBE. This mechanism is limited to two additions of O₂ to fuel radicals (•C₈H₁₇O), which yield ketohydroperoxides (C₈H₁₆O₄).

**4. Conclusions**

The low-temperature (460–780 K) and high-pressure (10 atm) oxidation of 2500 ppm of di-n-butyl ether was studied in a jet-stirred reactor. Due to the use of high-resolution mass spectrometry (Orbitrap Q-Exactive) coupled to liquid chromatography (HPLC and UHPLC) and atmospheric pressure chemical ionization (APCI), a wide range of previously observed and newly observed low-temperature oxidation products were detected. They include, besides light ubiquitous hydrocarbons’ oxidation products, species resulting from the one to four O₂ additions on the fuel’s radicals, i.e., hydroperoxides and di-hydroperoxides (C₈H₁₈O₃–5), unsaturated hydroperoxides, di and tri-hydroperoxides (C₈H₁₆O₃, 5, 7), ketohydroperoxides (C₃H₆O₃, C₄H₈O₃, C₅H₁₀O₄, C₆H₁₂O₄, and C₈H₁₆O₄), and highly oxygenated molecules (C₈H₁₆O₆, 8, C₈H₁₄O₅, 7, and C₈H₁₂O₄, 6). Oxidation products more specific to atmospheric oxidation chemistry, i.e., C₈H₁₆O₃, 5, 7 with carbonyl and hydroxyl groups, and one or more hydroperoxyl groups, i.e., C₈H₁₈O₂ corresponding to butoxy-butanol species, and ROOR’ compounds, including C₁₆H₃₄O₄, C₁₁H₂₄O₃, C₁₁H₂₂O₃, and C₁₀H₂₂O₃, were also detected. The present results allowed for highlighting the finding that common reaction pathways pertaining to combustion and atmospheric autoxidation chemistry occur under the present JSR experimental conditions. Further studies involving complementary techniques, e.g., SVUV-PI-MBMS, would be useful for confirming the proposed identification of isomeric products such as cyclic ethers vs. carbonyls. FIA indicated the formation of high molecular weight oxygenated products (m/z > 550) under the present experimental conditions. Accretion products are of high importance for the formation of atmospheric particulates. Then, it would be interesting to investigate the formation of the newly detected species in this study using a smog chamber or flow reactors to further clarify similarities in terms of oxidation pathways in combustion and atmospheric oxidation conditions.

**Supplementary Materials:** Figure S1: Evolution of O₂, H₂O, CO, CO₂, and H₂, and formic acid experimental and computed mole fraction profiles obtained for the low-temperature oxidation of 2500 ppm of di-n-butyl ether at 10 atm, φ = 2, and a residence time of 1 s; Figure S2: Comparison between experimental mole fractions of organic species and simulations for the oxidation of 2500 ppm of di-n-butyl ether in JSR (pressure of 10 atm, low-temperature 460–780 K, φ = 2, and a residence time of 1 s); Figure S3: Mole fraction profiles of aldehydes (formaldehyde, acetaldehyde, propionaldehyde, and butyraldehyde), ketones (acetone and 2-butanone), carboxylic acids (acetic acid and butyric acid), and 2,3-DHF formed during the low-temperature oxidation of 2500 ppm of di-n-butyl ether at 10 atm; Figure S4: Chromatograms showing C₃, 4, 5, 6, 8 ketohydroperoxides formed during the oxidation of 2500 ppm of di-n-butyl ether in a JSR (pressure: 10 atm; temperature: 460–780 K; φ = 2; and residence time = 1 s); Figure S5: Top: chromatogram showing the formation of C₈H₁₄O₃ diones during the oxidation of di-n-butyl ether. Butyric anhydride was identified at 7.18 min. Analyses were performed using a C18 UHPLC column with positive APCI HRMS (C₈H₁₅O₃⁺ with m/z 159.1015). Bottom: chromatogram showing the formation of cyclic ethers, butyl butyrate, and butoxy-carbonyls (C₈H₁₆O₂) during the oxidation of di-n-butyl ether; Figure S6: Chromatograms showing HOMs formed during the low-temperature oxidation of 2500 ppm of di-n-butyl ether at 10 atm; Figure S7: Chromatograms showing C₈H₁₆O₃, 5, 7 obtained before and after DNPH derivatization; Figure S8: Variation of the OSc as a function of the number of carbon atoms in molecular formulae detected by the negative FIA-HRMS of a limonene oxidation sample collected at 590 K and 1 atm; and Table S1: Conditions for GC, normal, and reverse-phase liquid chromatography and FTIR analyses.

**Author Contributions:** Conceptualization, P.D.; methodology, P.D., N.B., M.L. and R.B.; validation, P.D.; formal analysis, N.B., M.L., R.B. and P.D.; investigation, N.B. and M.L.; data curation, N.B., M.L. and R.B.; writing—original draft preparation, N.B. and P.D.; writing—review and editing, P.D., N.B. and M.L.; supervision, P.D. All authors have read and agreed to the published version of the manuscript.
**Funding:** This research study was funded by the Labex CAPRYSSES (grant number ANR-11-LABX-0006-01) from the Labex VOLTAIRE (grant number ANR-10-LABX-100-01) and from the Région Centre-Val-de-Loire, EFRD, and CPER (projects PROMESTOCK and APROPOR-E).

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** Data are available from the corresponding author.

**Conflicts of Interest:** The authors declare no conflict of interest.

## References

1. Lelieveld, J.; Klingmüller, K.; Pozzer, A.; Burnett, R.T.; Haines, A.; Ramanathan, V. Effects of fossil fuel and total anthropogenic emission removal on public health and climate. *Proc. Natl. Acad. Sci. USA* **2019**, *116*, 7192–7197. [CrossRef] [PubMed]

2. Olson, C.; Lenzmann, F. The social and economic consequences of the fossil fuel supply chain. *MRS Energy Sustain.* **2016**, *3*, 6. [CrossRef]

3. Muscat, A.; de Olde, E.M.; de Boer, I.J.M.; Ripoll-Bosch, R. The battle for biomass: A systematic review of food-feed-fuel competition. *Glob. Food Secur.* **2020**, *25*, 100330. [CrossRef]

4. Rotavera, B.; Taatjes, C.A. Influence of functional groups on low-temperature combustion chemistry of biofuels. *Prog. Energy Combust. Sci.* **2021**, *86*, 100925. [CrossRef]

5. Martinez, M.; Duret, X.; Pham Minh, D.; Nzihou, A.; Lavoie, J.-M. Conversion of lignocellulosic biomass in biobutanol by a NOVEL thermal process. *Int. J. Energy Prod. Manag.* **2019**, *4*, 298–310. [CrossRef]

6. Guan, L.; Tang, C.; Yang, K.; Jun, M. Effect of di-n-butyl ether blending with soybean-biodiesel on spray and atomization characteristics in a common-rail fuel injection system. *Fuel* **2015**, *140*, 116–125. [CrossRef]

7. Guan, L.; Tang, C.; Yang, K.; Mo, J.; Huang, Z. Experimental and Kinetic Study on Ignition Delay Times of Di-n-butyl Ether at High Temperatures. *Energy Fuels* **2014**, *28*, 5489–5496. [CrossRef]

8. Yuanhang, G.; Liu, W.; Han, D. Comparative Study on Spray Auto-Ignition of Di-n-Butyl Ether and Diesel Blends at Engine-Like Conditions. *J. Energy Resour. Technol.* **2021**, *143*, 042302.

9. Garcia, A.; Monsalve-Serrano, J.; Villalta, D.; Zubel, M.; Pischinger, S. Potential of 1-octanol and di-n-butyl ether (DNBE) to improve the performance and reduce the emissions of a direct injected compression ignition diesel engine. *Energy Convers. Manag.* **2018**, *177*, 563–571. [CrossRef]

10. Gao, Z.; Zhu, L.; Zou, X.; Liu, C.; Tian, B.; Huang, Z. Soot reduction effects of dibutyl ether (DBE) addition to a biodiesel surrogate in laminar coflow diffusion flames. *Proc. Combust. Inst.* **2019**, *37*, 1265–1272. [CrossRef]

11. Ruiz-Rodriguez, I.; Cracknell, R.; Parkes, M.; Megaritis, T.; Ganippa, L. Experimental study of the effect of C8 oxygenates on sooting processes in high pressure spray flames. *Combust. Flame* **2020**, *220*, 235–246. [CrossRef]

12. Wullenkord, J.; Tran, L.-S.; Böttchers, J.; Graf, I.; Kohse-Höinghaus, K. A laminar flame study on di-n-butyl ether as a potential biofuel candidate. *Combust. Flame* **2018**, *190*, 36–49. [CrossRef]

13. Fan, X.; Liu, Z.; Yang, J.; Yang, B. Pyrolysis of Lignocellulosic Biofuel Di-n-butyl Ether (DBE): Flow Reactor Experiments and Kinetic Modeling. *Energy Fuels* **2021**, *35*, 14077–14086. [CrossRef]

14. Thion, S.; Togbe, C.; Serinyel, Z.; Dayma, G.; Dagaut, P. A chemical kinetic study of the oxidation of dibutyl-ether in a jet-stirred reactor. *Combust. Flame* **2017**, *185*, 4–15. [CrossRef]

15. Tran, L.S.; Wullenkord, J.; Li, Y.Y.; Herbinet, O.; Zeng, M.R.; Qi, F.; Kohse-Hoinghaus, K.; Battin-Leclerc, F. Probing the low-temperature chemistry of di-n-butyl ether: Detection of previously unobserved intermediates. *Combust. Flame* **2019**, *210*, 9–24. [CrossRef]

16. Belhadj, N.; Benoit, R.; Dagaut, P.; Lailliau, M.; Serinyel, Z.; Dayma, G.; Khaled, F.; Moreau, B.; Foucher, F. Oxidation of di-n-butyl ether: Experimental characterization of low-temperature products in JSR and RCM. *Combust. Flame* **2020**, *222*, 133–144. [CrossRef]

17. Wang, Z.; Popolan-Vaida, D.M.; Chen, B.; Moshammer, K.; Mohamed, S.Y.; Wang, H.; Sioud, S.; Raji, M.A.; Kohse-Höinghaus, K.; Hansen, N.; et al. Unraveling the structure and chemical mechanisms of highly oxygenated intermediates in oxidation of organic compounds. *Proc. Natl. Acad. Sci. USA* **2017**, *114*, 13102–13107. [CrossRef] [PubMed]

18. Wang, Z.D.; Sarathy, S.M. Third O-2 addition reactions promote the low-temperature auto-ignition of n-alkanes. *Combust. Flame* **2016**, *165*, 364–372. [CrossRef]

19. Wang, Z.; Zhang, L.; Moshammer, K.; Popolan-Vaida, D.M.; Shankar, V.S.B.; Lucassen, A.; Hemken, C.; Taatjes, C.A.; Leone, S.R.; Kohse-Hoeinghaus, K.; et al. Additional chain-branching pathways in the low-temperature oxidation of branched alkanes. *Combust. Flame* **2016**, *164*, 386–396. [CrossRef]

20. Glarborg, P.; Kee, R.J.; Grcar, J.F.; Miller, J.A. PSR: A Fortran Program for Modeling Well-Stirred Reactors; SAND86-8209; Sandia National Laboratories: Livermore, CA, USA, 1986.

21. Fan, X.F.; Sun, W.Y.; Gao, Y.; Hansen, N.; Chen, B.J.; Pitsch, H.; Yang, B. Chemical insights into the multi-regime low-temperature oxidation of di-n-propyl ether: Jet-stirred reactor experiments and kinetic modeling. *Combust. Flame* **2021**, *233*, 11592. [CrossRef]

22. Blin-Simiand, N.; Jorand, F.; Sahetchian, K.; Brun, M.; Kerhoas, L.; Malosse, C.; Einhorn, J. Hydroperoxides with zero, one, two or more carbonyl groups formed during the oxidation of n-dodecane. *Combust. Flame* **2001**, *126*, 1524–1532. [CrossRef]
23. Belhadj, N.; Lailliau, M.; Benoit, R.; Dagaut, P. Experimental and kinetic modeling study of n-hexane oxidation. Detection of complex low-temperature products using high-resolution mass spectrometry. *Combust. Flame* **2021**, *233*, 111581. [CrossRef]

24. Bianchi, F.; Kurtén, T.; Riva, M.; Mohr, C.; Rissanen, M.P.; Roldin, P.; Berndt, T.; Crounse, J.D.; Wennberg, P.O.; Mentel, T.F.; et al. Highly Oxygenated Organic Molecules (HOM) from Gas-Phase Autoxidation Involving Peroxy Radicals: A Key Contributor to Atmospheric Aerosol. *Chem. Rev.* **2019**, *119*, 3472–3509. [CrossRef] [PubMed]

25. Benoit, R.; Belhadj, N.; Lailliau, M.; Dagaut, P. On the similarities and differences between the products of oxidation of hydrocarbons under simulated atmospheric conditions and cool flames. *Atmos. Chem. Phys.* **2021**, *21*, 7845–7862. [CrossRef]

26. Kroll, J.H.; Donahue, N.M.; Jimenez, J.L.; Kessler, S.H.; Canagaratna, M.R.; Wilson, K.R.; Altieri, K.E.; Mazzoleni, L.R.; Wozniak, A.S.; Bluhm, H.; et al. Carbon oxidation state as a metric for describing the chemistry of atmospheric organic aerosol. *Nat. Chem.* **2011**, *3*, 133–139. [CrossRef] [PubMed]

27. Dagaut, P.; Cathonnet, M.; Boettner, J.C.; Gaillard, F. Kinetic modeling of ethylene oxidation. *Combust. Flame* **1988**, *71*, 295–312. [CrossRef]

28. Dagaut, P.; Cathonnet, M.; Boettner, J.C.; Gaillard, F. Kinetic Modeling of Propane Oxidation. *Combust. Sci. Technol.* **1987**, *56*, 23–63. [CrossRef]

29. Dagaut, P.; Cathonnet, M.; Rouan, J.P.; Foulatier, R.; Quilgars, A.; Boettner, J.C.; Gaillard, F.; James, H. A Jet-Stirred Reactor for Kinetic-Studies of Homogeneous Gas-Phase Reactions at Pressures up to 10-Atmospheres (~1 MPa). *J. Phys. E-Sci. Instrum.* **1986**, *19*, 207–209. [CrossRef]