Removing Simultaneously Sulfur and Nitrogen from Fuel under a Sustainable Oxidative Catalytic System

: An effective process to remove nitrogen-based compounds from fossil fuels without harming the process of sulfur removal is an actual gap in reﬁneries. A success combination of desulfurization and denitrogenation processes capable of completely removing the most environmental contaminates in diesel under sustainable conditions was achieved in this work, applying polyoxometalates as catalysts, hydrogen peroxide as oxidant, and an immiscible ionic liquid as an extraction solvent. The developed process based in simultaneous oxidative desulfurization (ODS) and oxidative denitrogenation (ODN) involved initial extraction of sulfur and nitrogen compounds followed by catalytic oxidation. Keggin-type polyoxomolybdates revealed much higher reusing capacity than the related polyoxotungstate. Effectively, the ﬁrst catalysts practically allowed complete sulfur and nitrogen removal only in 1 h of reaction and for ten consecutive cycles, maintaining the original catalyst and ionic liquid samples.


Introduction
Long-term energy outlooks agree that fossil fuels will remain the dominant energy for at least more two decades, since the demand for energy is increasing with the growth of the world's population, economy, and increasing living standards. At present, fossil fuels cover~88% of the energy requirement in the world [1]. The inevitable drawbacks of the current fossil fuel consumption are related to the environment degradation and public health. In the transportation sector, fossil fuels will be the main energy source for the near future. Motor vehicles powered on fossil fuels are among the major sources of air pollution in urban areas [2,3]. Fuels such as fuel oil, jet fuel, gasoline, diesel, etc., contain various pollutants, with the sulfur-and nitrogen-based organic compounds (SCs and NCs, respectively) as the majority of the contaminants [4]. The utilization of fossil fuels without purification poses a devastating effect on the environment and humans, mainly because of the emission of sulfur oxides (SOx) and nitrogen oxides (NOx) from SCs and NCs, respectively [5,6]. The global impact of SOx emissions resulting from the burning of fossil fuels has led to the imposition of environmental restrictions for the sulfur level in transportation fuels (10 ppm in EU) [7,8].
The standard industrial methods for the removal of sulfur and nitrogen-based compounds from crude oils are hydrodesulfurization (HDS) and hydrodenitrogenation (HDN) [9,10]. HDS is an effective process that needs high temperature (>350 • C) and pressure (up to 6 MPa), and also requires high amounts of hydrogen consumption [11][12][13]. SCs and NCs in fuels are converted to hydrogen sulfide and ammonia, respectively, via Sustain. Chem. 2021, 2 383 the catalytic reduction with H 2 [14] In the HDS process, NCs severely affect the catalyst and cause corrosion of the refinery equipment due to their high reactivity [14][15][16]. Therefore, it is essential to remove NCs before the desulfurization treatment of fuel, or preferably, to remove NCs and SCs simultaneously by an alternative process, capable of depth desulfurization and denitrogenation, reducing the impact of fuels on the environment by decreasing NOx and SOx emissions. Oxidative desulfurization (ODS) and oxidative denitrogenation (ODN) processes have been the focus of promising research in recent years due to their efficiency for removal of the most refractory S and N-content present in fuels [10,17,18]. This technology conciliates oxidative catalysis and liquid-liquid extraction steps [15,16]. Consequently, the efficiency of combined ODS/ODN processes depends of the catalyst performance and the ability of the extraction solvent. The sustainability and the cost-effectiveness of the process can be assured by the use of recovered and recyclable catalysts and the application of environmental-friendly solvents, such as ionic liquids [19,20]. The former are able to behave as an extraction solvent and also as an immobilization medium for the homogeneous catalyst, containing catalyst and subtract in the same liquid phase [21,22]. Hydrogen peroxide is generally the preferred oxidant to use in these processes due to its high active oxygen content and the fact that water is the sole by-product [23,24]. Our research group has been dedicated to the design of efficient catalysts for ODS based on polyoxometalates (POMs) [18,[25][26][27][28][29] POMs are polyanions incorporating various transition metals and are well-known by their structural diversity and potential in several areas of application [30]. In general, POMs are powerful catalysts for a great variety of oxidative reactions. In particular, Keggin-type POMs ([X n+ M 12 O 40 ] (8-n)-(X: block p or d heteroatom) have been successfully used as efficient catalysts in ODS processes [31][32][33][34][35]. However, the investigation combining ODS and ODN is in a primitive stage and only one publication could be found in the literature demonstrating the ability of POMs as catalysts in this mixed systems [36].

Characterization Methods and Instrumentation
Fourier-transformed Infrared (FT-IR) spectra were acquired between 400 and 4000 cm −1 on a Bruker Tensor 27 Spectrometer using KBr pellets, and all the representations are shown in arbitrary unities of transmittance. Solution 31 P nuclear magnetic resonance (NMR) spectra were recorded in CD 3 CN or D 2 O at 162 MHz with a Bruker Avance III 400 spectrometer, and the chemical shifts are given with respect to external 85% H 3 PO 4 . Catalytic reactions were periodically monitored by GC-FID analysis carried out in a Bruker 430-GC-FID chromatograph (Germany). Hydrogen was used as carrier gas (55 cm.s −1 ) and fused silica Supelco capillary columns SPB-5 (30 m × 0.25 mm i.d.; 25 µm film thickness) were used.

[BPy] 3 Salts Preparation
[BPy] 3 12 ) were prepared following an experimental procedure adapted from previously reported methods [34,35]. Briefly, [BPy]Br (2.5 mmol) was dissolved in 2.5 mL of deionized water, and an aqueous solution of phosphomolybdic or phosphotungstic acid (0.5 mmol) in deionized water (2.5 mL) was added dropwise under magnetic stirring at room temperature. The mixture was left stirring for 30 min, after which the precipitate was filtrated, washed thoroughly with deionized water, and dried in a desiccator overnight.   Figure S1 in Supporting Information), and it was found to be approximately 5% to ODS and 4% to ODN.

Catalysts Synthesis and Characterization
The catalysts [BPy]PMo 12 and [BPy]PW 12 are POM salts previously reported by our group [34,35] Several weak/medium bands in the range 1580-1100 cm −1 are observed due to aromatic ν(C-N), ν(C-C), and δ(C-H) vibrations. A strong band at 1485 cm −1 is attributable to the conjugation of ν(C=C) and ν(C=N) bonds, typical of pyridinium salts. The strong one at 682 cm −1 corresponds to the aromatic out-of-plane hydrogen deformation [40]. 31 P NMR analysis reveals the typical chemical shifts for each structure [34,35]. The spectrum recorded for [Na]PMo 12 registers two single peaks attributed to a different content of Na cations in the POM structure.
Sustain. Chem. 2021, 2, FOR PEER REVIEW 4 infrared spectra of both [BPy]3 salts contain weak bands at 3126, 3085, and 3065 cm −1 corresponding to ν(C-H) aromatic bond vibrations, as well as at 2964, 2930, and 2874 cm −1 attributed to the ν(C-H) of the aliphatic chain, a medium intensity band at 1632 cm −1 corresponding to ν(C=N) vibrations that is typical of the quaternary nitrogen atom in a heterocyclic ring. Several weak/medium bands in the range 1580-1100 cm −1 are observed due to aromatic ν(C-N), ν(C-C), and δ(C-H) vibrations. A strong band at 1485 cm −1 is attributable to the conjugation of ν(C=C) and ν(C=N) bonds, typical of pyridinium salts. The strong one at 682 cm −1 corresponds to the aromatic out-of-plane hydrogen deformation [40]. 31 P NMR analysis reveals the typical chemical shifts for each structure [34,35]. The spectrum recorded for [Na]PMo12 registers two single peaks attributed to a different content of Na cations in the POM structure.

ODS/ODN Catalytic Studies
All ODS/ODN studies were performed at 70 °C in a biphasic liquid

ODS/ODN Catalytic Studies
All ODS/ODN studies were performed at 70 • C in a biphasic liquid-liquid system based on [BMIM]PF 6 as extraction solvent and an immiscible multicomponent S/N Model Diesel containing the most representative refractory sulfur and nitrogen content present in fuels, namely benzothiophene (BT), dibenzothiophene (DBT), 4-methyldibenzothiophene (MDBT), and 4,6-dimethyldibenzothiophene (DMDBT) from SCs, and indole (IND) and quinoline (QUI) from NCs. Catalytic systems based on [Na] 3  These results indicate that the nature of the cation in the POM catalyst does not seem to have an important influence in its catalytic performance. On the other hand, a distinct combined ODS/ODN profile was found using polyoxotungstate and polyoxomolybdate, where the first showed to be the most effective, guaranteeing total desulfurization and denitrogena-tion after just 1 h of reaction. The removal rates of each component of the S/N Model Diesel are the following: IND > QUI > DBT > MDBT > DMDBT~BT. The removal of NCs is largely achieved during the initial extraction step. The difference in removal rates between IND and QUI relates to the higher basicity of the proton-donor present in IND [41,42]. During the oxidation catalytic step (after the 10 min), the difference of reactivity between the SCs can be explained by the decrease in electron density on the sulfur atom, which hinders oxidation [18]. The nature of the obtained oxidized products was identified by the analysis of the extraction phase, since, in the model diesel phase, no oxidized products were found. The oxidized products of sulfur compounds were the corresponding sulfones. Only vestigial amounts of sulfoxide from BT were found. From the QUI oxidation, the quinoline N-oxide was identified.
Sustain. Chem. 2021, 2, FOR PEER REVIEW 5 [BPy]PMo12 catalyst. These results indicate that the nature of the cation in the POM catalyst does not seem to have an important influence in its catalytic performance. On the other hand, a distinct combined ODS/ODN profile was found using polyoxotungstate and polyoxomolybdate, where the first showed to be the most effective, guaranteeing total desulfurization and denitrogenation after just 1 h of reaction. The removal rates of each component of the S/N Model Diesel are the following: IND > QUI > DBT > MDBT > DMDBT ~ BT. The removal of NCs is largely achieved during the initial extraction step.
The difference in removal rates between IND and QUI relates to the higher basicity of the proton-donor present in IND [41,42]. During the oxidation catalytic step (after the 10 min), the difference of reactivity between the SCs can be explained by the decrease in electron density on the sulfur atom, which hinders oxidation [18]. The nature of the obtained oxidized products was identified by the analysis of the extraction phase, since, in the model diesel phase, no oxidized products were found. The oxidized products of sulfur compounds were the corresponding sulfones. Only vestigial amounts of sulfoxide from BT were found. From the QUI oxidation, the quinoline N-oxide was identified.  Figure 3 displays the denitrogenation and the desulfurization profiles of the combined ODS/ODS processes. These results demonstrated that complete denitrogenation is achieved faster than desulfurization, since higher extraction of NCs occurred during the first 10 min of the process, which is associated to the higher oxidative facility of NCs than SCs. Figure 3 also compares the desulfurization profiles that occurred in a combined ODS/ODN process and a single ODS process, i.e., using a Model Diesel containing only SCs. These experiments were performed using the same reactional conditions. Without duties, the desulfurization profiles are similar in combined S/N and single S processes. These results indicate that the conciliation of extraction and oxidative catalytic steps did not result in a competitive S and N removal process, and these may occur in parallel, sharing the oxidant and the catalyst action. This is a remarkable advantage compared to the combined hydrodesulfuriation/hydrodenitrogenation processes [14].  Figure 3 displays the denitrogenation and the desulfurization profiles of the combined ODS/ODS processes. These results demonstrated that complete denitrogenation is achieved faster than desulfurization, since higher extraction of NCs occurred during the first 10 min of the process, which is associated to the higher oxidative facility of NCs than SCs. Figure 3 also compares the desulfurization profiles that occurred in a combined ODS/ODN process and a single ODS process, i.e., using a Model Diesel containing only SCs. These experiments were performed using the same reactional conditions. Without duties, the desulfurization profiles are similar in combined S/N and single S processes. These results indicate that the conciliation of extraction and oxidative catalytic steps did not result in a competitive S and N removal process, and these may occur in parallel, sharing the oxidant and the catalyst action. This is a remarkable advantage compared to the combined hydrodesulfuriation/hydrodenitrogenation processes [14]. Sustain  The recycling ability of each POM/[BMIM]PF6 catalytic system was accessed by consecutive reutilization cycles that consisted in the preservation of the POM/[BMIM]PF6 phase for various consecutive cycles. This reusing process is the desired sustainable procedure of using the same portion of catalyst and ionic liquid extraction solvent for consecutive combined ODS/ODN cycles. After each cycle, the treated S/N Model Diesel is removed from the system and replaced by a novel portion of S/N Model Diesel, as well the H2O2 oxidant. During the reutilization cycles, an interlayer aqueous phase increased in volume between the S/N Model Diesel and [BMIM]PF6 phase due to the successive addition of aqueous oxidizing agent (Figure 4a). Figure 4b display the results obtained from reusing the three different POM/[BMIM]PF6 catalytic systems. These results were achieved after 1 h. The [BPy]PW12/[BMIM]PF6 system presented a decrease of catalytic efficiency, mainly after the 3 rd cycle, since its combined ODS/ODN efficiency decreased from 98% from the 3 rd cycle to 63% for the 4 th cycle. In fact, the efficiency found for the 4 th cycle was only attributed to the initial extraction step, indicating the absence of the oxidative catalytic performance of [BPy]PW12 catalyst. In the 5 th combined desulfurization/denitrogenation cycle, the catalyst did not promote any oxidative catalysis, and the initial extraction of the S/N compounds decreased to 46%. This decrease of S/N extraction is a consequence of absence of S/N oxidation, i.e., without decreasing non-oxidized SCs and NCs in the [BMIM]PF6 phase, the extraction of more S/N from Diesel did not occur.   6 system presented a decrease of catalytic efficiency, mainly after the 3 rd cycle, since its combined ODS/ODN efficiency decreased from 98% from the 3 rd cycle to 63% for the 4 th cycle. In fact, the efficiency found for the 4 th cycle was only attributed to the initial extraction step, indicating the absence of the oxidative catalytic performance of [BPy]PW 12 catalyst. In the 5 th combined desulfurization/denitrogenation cycle, the catalyst did not promote any oxidative catalysis, and the initial extraction of the S/N compounds decreased to 46%. This decrease of S/N extraction is a consequence of absence of S/N oxidation, i.e., without decreasing nonoxidized SCs and NCs in the [BMIM]PF 6 phase, the extraction of more S/N from Diesel did not occur. Reusing experiments using the [BPy]PW 12 /[BMIM]PF 6 system were not performed after the 5 th cycle. On the other hand, using the polyoxomolybdate catalytic systems, i.e., [Na]PMo 12 and [BPy]PMo 12 , the loss of the catalytic activity was not verified for ten consecutive combined ODS/ODN cycles and near complete removal of S/N from Model Diesel was achieved after 1 h. A small decrease in activity is only observed during the 10 th cycle (93 and 98%, for [Na]PMo 12 and [BPy]PMo 12 , respectively).   After the reusing cycles, the stability of the homogeneous POM catalysts was investigated by 31 P NMR. The immiscible aqueous and [BMIM]PF6 phases (Figure 4b) were separated, analyzed, and the results are displayed in Figure 5. In general, all studied catalysts suffered structural transformations during the various reusing cycles. Furthermore, the structural transformations may alter solubility properties, since POM derivative peaks could be found in both the Ionic Liquid and aqueous phase. The structural modification of Na[PMo12] and [BPy]PMo12 was not coincident, and different active POM derivative fragments were found in aqueous and [BMIM]PF6 phases. On the other hand, the analysis of both phases obtained after the 5th cycle using the [BPy]PW12 catalyst revealed similar distribution of POM transformed species. It ends up being implied that the structural transformations associated with [BPy]3PW12O40 catalytic activity ultimately led to its deactivation. In a future work, the stability and the facility of the recovery of these homogeneous catalysts will be treated by their immobilization in appropriate supporting materials.
(a) (b) After the reusing cycles, the stability of the homogeneous POM catalysts was investigated by 31 P NMR. The immiscible aqueous and [BMIM]PF 6 phases ( Figure 4b) were separated, analyzed, and the results are displayed in Figure 5. In general, all studied catalysts suffered structural transformations during the various reusing cycles. Furthermore, the structural transformations may alter solubility properties, since POM derivative peaks could be found in both the Ionic Liquid and aqueous phase. The structural modification of Na[PMo 12 ] and [BPy]PMo 12 was not coincident, and different active POM derivative fragments were found in aqueous and [BMIM]PF 6 phases. On the other hand, the analysis of both phases obtained after the 5 th cycle using the [BPy]PW 12 catalyst revealed similar distribution of POM transformed species. It ends up being implied that the structural transformations associated with [BPy] 3 PW 12 O 40 catalytic activity ultimately led to its deactivation. In a future work, the stability and the facility of the recovery of these homogeneous catalysts will be treated by their immobilization in appropriate supporting materials. After the reusing cycles, the stability of the homogeneous POM catalysts was investigated by 31 P NMR. The immiscible aqueous and [BMIM]PF6 phases (Figure 4b) were separated, analyzed, and the results are displayed in Figure 5. In general, all studied catalysts suffered structural transformations during the various reusing cycles. Furthermore, the structural transformations may alter solubility properties, since POM derivative peaks could be found in both the Ionic Liquid and aqueous phase. The structural modification of Na[PMo12] and [BPy]PMo12 was not coincident, and different active POM derivative fragments were found in aqueous and [BMIM]PF6 phases. On the other hand, the analysis of both phases obtained after the 5th cycle using the [BPy]PW12 catalyst revealed similar distribution of POM transformed species. It ends up being implied that the structural transformations associated with [BPy]3PW12O40 catalytic activity ultimately led to its deactivation. In a future work, the stability and the facility of the recovery of these homogeneous catalysts will be treated by their immobilization in appropriate supporting materials.

Conclusions
This work reported a novel efficient strategy capable of combining complete desulfurization and denitrogenation to treat a multicomponent S/N Model Diesel containing the most refractory elements present in Diesel. Keggin-type POMs showed to be active catalysts to perform parallel oxidative desulfurization (ODS) and oxidative denitrogenation (ODS) under sustainable conditions, i.e., using hydrogen peroxide as oxidant and an ionic liquid as extraction solvent ([BMIM]PF6). Complete removal of sulfur and nitrogen compounds was achieved only after 1 h. Furthermore, a reusing test was performed by using the homogeneous POM catalyst and [BMIM]PF6 for consecutive ODS/ODN cycles. The polyoxomolybdates [Na]PMo12 and [BPy]PMo12 showed high reusing capacity for ten consecutive cycles without any appreciable reduction of catalytic performance. The kinetic similarity obtained between these two polyoxomolybdates-based compounds indicated that the nature of the cation does not seem to influence its activity. The reusing capacity was drastically decreased when the polyoxotungstate [BPy]PW12 was used and, in this case, only three consecutive ODS/ODN cycles were successfully performed. Lastly, the stability analysis of the POM catalysts after reusing cycles indicates that the polyoxomolybdate is in fact a precursor for the active catalytic species, since its Keggin structures are modified in other active compounds. On the other hand, the polyoxotungstate also suffered structural modification but into inactive species. Therefore, future work is already planned to prepare heterogeneous POM catalysts to avoid their fast structural modification and also to facilitate catalyst recovery.

Conclusions
This work reported a novel efficient strategy capable of combining complete desulfurization and denitrogenation to treat a multicomponent S/N Model Diesel containing the most refractory elements present in Diesel. Keggin-type POMs showed to be active catalysts to perform parallel oxidative desulfurization (ODS) and oxidative denitrogenation (ODS) under sustainable conditions, i.e., using hydrogen peroxide as oxidant and an ionic liquid as extraction solvent ([BMIM]PF 6 ). Complete removal of sulfur and nitrogen compounds was achieved only after 1 h. Furthermore, a reusing test was performed by using the homogeneous POM catalyst and [BMIM]PF 6 for consecutive ODS/ODN cycles. The polyoxomolybdates [Na]PMo 12 and [BPy]PMo 12 showed high reusing capacity for ten consecutive cycles without any appreciable reduction of catalytic performance. The kinetic similarity obtained between these two polyoxomolybdates-based compounds indicated that the nature of the cation does not seem to influence its activity. The reusing capacity was drastically decreased when the polyoxotungstate [BPy]PW 12 was used and, in this case, only three consecutive ODS/ODN cycles were successfully performed. Lastly, the stability analysis of the POM catalysts after reusing cycles indicates that the polyoxomolybdate is in fact a precursor for the active catalytic species, since its Keggin structures are modified in other active compounds. On the other hand, the polyoxotungstate also suffered structural modification but into inactive species. Therefore, future work is already planned to prepare heterogeneous POM catalysts to avoid their fast structural modification and also to facilitate catalyst recovery.