ORIGINAL RESEARCH article
Mesoporous Silica vs. Organosilica Composites to Desulfurize Diesel
- 1LAQV-REQUIMTE, Departamento de Química e Bioquímica, Faculdade de Ciências, Universidade Do Porto, Porto, Portugal
- 2CENIMAT/I3N, Departamento de Ciência dos Materiais, Faculdade de Ciências e Tecnologia, Universidade Nova de Lisboa, Caparica, Portugal
- 3Faculdade de Ciências, Centro de Química e Bioquímica and CQE, Universidade de Lisboa, Lisbon, Portugal
- 4Grupo de Energía y Química Sostenibles (EQS), Instituto de Catálisis y Petroleoquímica, CSIC, Madrid, Spain
The monolacunary Keggin-type [PW11O39]7− (PW11) heteropolyanion was immobilized on porous framework of mesoporous silicas, namely SBA-15 and an ethylene-bridged periodic mesoporous organosilica (PMOE). The supports were functionalized with a cationic group (N-trimethoxysilypropyl-N, N, N-trimethylammonium, TMA) for the successful anchoring of the anionic polyoxometalate. The PW11@TMA-SBA-15 and PW11@TMA-PMOE composites were evaluated as heterogeneous catalysts in the oxidative desulfurization of a model diesel. The PW11@TMA-SBA-15 catalyst showed a remarkable desulfurization performance by reaching ultra-low sulfur levels (<10 ppm) after only 60 min using either a biphasic extractive and catalytic oxidative desulfurization (ECODS) system (1:1 MeCN/diesel) or a solvent-free catalytic oxidative desulfurization (CODS) system. Furthermore, the mesoporous silica composite was able to be recycled for six consecutive cycles without any apparent loss of activity. The promising results have led to the application of the catalyst in the desulfurization of an untreated real diesel supplied by CEPSA (1,335 ppm S) using the biphasic system. The system has proved to be a highly efficient process by reaching desulfurization values higher than 90% for real diesel during three consecutive cycles.
Over the last years, ordered mesoporous silicas (OMS) have attracted researchers' attention in catalysis, due to its long range order, high surface areas and well-defined pore size (2–50 nm) (Esquivel et al., 2014; Wang et al., 2014). Moreover, the surface of these materials can be easily modified through reaction with organosilanes, incorporating proper functional groups on their surfaces that will guarantee the effective immobilization of appropriate compounds, such as active catalytic centers that will form novel heterogeneous catalysts.
Different mesoporous silica families (SBA-n, PMOs, etc.) have been used to create active heterogeneous catalysts, some of these having Keggin-type polyoxometalates (POMs) as active centers (Zhu et al., 2009). The morphology of these silica based materials and the functional groups present at their surface are crucial to the structural robustness of these catalysts.
Periodic mesoporous organosilicas (PMOs) are a recent class of ordered organic-inorganic hybrid mesoporous materials. PMOs allow the tuning of the superficial and physical properties besides being structurally robust (Esquivel et al., 2014; Park et al., 2014). Usually, the preparation of PMOs is conducted, in the presence of a structure-directing agent, by hydrolysis and condensation reactions of bridged silsesquioxane precursors with general formula (R′O)3-Si–R–Si–(OR′)3. The synthetic procedure for PMOs is similar to the preparation of mesoporous silicas, such as SBA-15 (Van Der Voort et al., 2013; Granadeiro et al., 2016). The functional organic moieties in PMOs are present at the surface but also in the channel walls, which are responsible for the structural features of the material (rigidity/flexibility) (Van Der Voort et al., 2013). The unique properties of PMOs make them suitable candidates for catalytic applications; however, the number of reports found in the literature is rather scarce and, in particular, its application in oxidative desulfurization is still inexistent (Karimi et al., 2012).
The demand for ultra-low sulfur fuels has motivated the development of new desulfurization methods, such as oxidative desulfurization (ODS), that brings technological, environmental, and economic advantages (Ribeiro et al., 2017; Yang et al., 2017; Liu et al., 2019). In ODS, the presence of an active catalyst with proper textural and surface chemical properties, as well as, an appropriate oxidant is the key to the success of this process. The monolacunary [PW11O39]7−Keggin polyoxometalate has shown remarkable performance in oxidative reactions (Casuscelli et al., 2004; Granadeiro et al., 2013; Ribeiro et al., 2013a; Singh et al., 2015; Coronel and da Silva, 2018), which has motivated its application in the oxidative desulfurization of fuels (Abdalla et al., 2009; Abdalla and Li, 2012; Wu et al., 2014; Mirante et al., 2018; Ribeiro et al., 2018b). In this work, novel composites have been prepared through the impregnation of the [PW11O39]7− (PW11) heteropolyanion in silica materials functionalized with N-trimethoxysilypropyl-N, N, N-trimethylammonium (TMA). Two different mesoporous silica supports were selected to prepare novel composites: the ordered mesoporous silica SBA-15 and an ethylene-bridged (PMOE) periodic mesoporous organosilica. Both materials were tested in the oxidative desulfurization of simulant diesel under biphasic and solvent-free systems. In this work, a POM-based PMO (PW11@TMA-PMOE) was used for the first time in the oxidative desulfurization of model fuel. Moreover, the composite exhibiting better desulfurization performance (PW11@TMA-SBA-15) with the model fuel was also tested in the sulfur removal of an untreated diesel supplied by CEPSA, under biphasic and solvent-free systems.
Materials and Methods
The reagents used for the synthesis of the materials and the ODS experiments were purchased to chemical suppliers and used without further purification (Ribeiro et al., 2018b). The materials prepared were characterized as previously described by our group for other mesoporous silica-based composites (Ribeiro et al., 2017). The ODS studies with model diesel were followed by GC-FID in a previously described setup while for the experiments with real diesel the sulfur content was determined by ultraviolet fluorescence at Galp (Julião et al., 2017).
Synthesis of PW11@TMA-SBA-15 Composite
The tetra-n-butylammonium [TBA, (C4H9)4N] salt of monolacunary phosphotungstate [PW11O39]7− (PW11) was prepared according to the literature (Ribeiro et al., 2018b). The mesoporous silica SBA-15 functionalized with N,N,N-trimethylammonium groups (TMA-SBA-15) was initially prepared followed by the impregnation of PW11 according to previously reported methods (Mirante et al., 2017; Ribeiro et al., 2018a).
TMA-SBA-15: Anal. Found (%): N, 1.4; C, 7.6; H, 2.2; loading of TMA: 0.098 mmol per g; selected FT-IR (cm−1): 3,436 (vs), 2,360 (w), 1,868 (w), 1,635 (m), 1,489 (m), 1,479 (m), 1,419 (w), 1,086 (vs), 951 (m), 806 (s), 692 (w), 678 (w), 553 (w), 460 (s); selected FT-Raman (cm−1): 3,028 (vs), 2,981 (sh), 2,972 (vs), 2,931 (m), 2,925 (sh), 2,823 (m), 1,450 (s), 912 (m).
PW11@TMA-SBA-15: Anal. Found (%) W, 21.3; Si, 4.7; loading of POM: 0.099 mmol per g, Si/W (molar) = 1.46; selected FT-IR (cm−1): 3,444 (vs), 2,360(m), 2,341 (m), 1,644 (m), 1,447 (m), 1,195 (sh), 1,085 (vs), 948 (m), 914 (w), 856 (w), 809 (m), 458 (s); selected FT-Raman (cm−1): 3,035 (s), 2,979 (s), 2,922 (s), 2,902 (m), 1,444 (m), 976 (vs), 970 (vs).
Synthesis of PW11@TMA-PMOE Composite
For the preparation of TMA-PMOE (Zhang et al., 2008; Granadeiro et al., 2016), Pluronic P123 (0.096 mmol) and KCl (47 mmol) were first dissolved in aqueous HCl (1.6 M; 29 mmol) at 40°C during 4 h. BTEE (2.84 mmol) and TMA (0.71 mmol) were added followed by stirring for 20 h at 40°C. The mixture was placed in an autoclave and heated at 100°C for 24 h. The solid was recovered by filtration and dried in a desiccator under silica gel. Finally, the solid was refluxed for 24 h in an acidic ethanol solution for the extraction of the co-polymer.
The PW11@TMA-PMOE composite was prepared by adding an acetonitrile solution of PW11 (1 g in 20 mL) to the TMA-PMOE (0.5 g) and allowing the mixture to stir at room temperature for 72 h. The solid was recovered by filtration, washed with acetonitrile and dried in a desiccator over silica gel.
TMA-PMOE: Anal. Found (%): N, 0.85; C, 19.79; H, 4.78; loading of TMA: 0.61 mmol per g; selected FT-IR (cm−1): 3,434 (vs), 2,886 (w), 2,360 (w), 1,637 (m), 1,488 (w), 1,413 (m), 1,270 (m), 1,159 (s), 1,093 (sh), 1,031 (vs), 912 (s), 763 (m), 696 (m), 441 (s); selected FT-Raman (cm−1): 2,972 (sh), 2,891 (vs), 2,802 (w), 1,441 (sh), 1,412 (m), 1,271 (m), 995 (m), 902 (w), 770 (m), 511 (s).
PW11@TMA-PMOE: Anal. Found (%) W, 10.24; Si, 1.55; loading of POM: 0.050 mmol per g, Si/W (molar) = 1.01; selected FT-IR (cm−1): 3,444 (vs), 2,975 (w), 2,898 (m), 2,360 (m), 2,343 (w) 1,652 (m), 1,488 (w), 1,417 (m), 1,270 (m), 1,159 (s), 1,079 (sh), 1,033 (vs), 908 (s), 815 (w), 767 (m), 696 (m), 501 (w), 441 (s); selected FT-Raman (cm−1): 3,031 (w), 2,968 (w), 2,891 (vs), 2,806 (w), 1,450 (m), 1,412 (s), 1,271 (m), 980 (s), 968 (sh), 928 (w), 760 (m), 515 (s).
Oxidative Desulfurization Processes Using Model Diesel
The ODS experiments were performed using a multicomponent model diesel composed by 1-BT, DBT, 4-MDBT, and 4,6-DMDBT (500 ppm each) in n-octane. The ODS studies were performed using either a biphasic extractive and catalytic oxidative desulfurization (ECODS) or a solvent-free catalytic oxidative desulfurization (CODS) system. For the biphasic system, a 1:1 mixture of model diesel/MeCN (1.5 mL) was added to the heterogeneous composite (the amount corresponding to 3 μmol of POM), stirred for 10 min and an aliquot of the upper phase oil was taken. Afterwards, aqueous hydrogen peroxide 30% (40 μL; H2O2/S molar = 8) was added initiating the catalytic step. For the solvent-free studies, model diesel (750 μL), 3 μmol of active catalyst and H2O2 oxidant (H2O2/S molar = 4) were used. A final liquid-liquid extraction, was performed to remove the oxidized sulfur compounds, using an extraction solvent such as MeCN. The reactions were monitored by GC analysis using tetradecane as a standard. At the end of oxidation, centrifugation was carried out to separate the solid catalyst, which was washed with ethanol and dried in a desiccator over silica gel. At the end of each cycle, the catalyst was recovered and reused in a new ODS cycle under the same reactional conditions.
Oxidative Desulfurization Process Using Untreated Diesel
An untreated diesel sample supplied by CEPSA (containing about 1,335 ppm of sulfur) was also desulfurized using the PW11@TMA-SBA-15 catalyst. The untreated diesel was mixed with the heterogeneous composite (an equivalent amount containing 3 μmol of PW11) in acetonitrile and with a H2O2/Sulfur ratio equal to 8. The mixture was heated at 70°C for 2 h, after which, the diesel was separated, washed with acetonitrile and separated by decantation. The CEPSA diesel was also used under solvent-free conditions by adding the catalyst and oxidant to diesel and heating at 70°C during 2 h. The mixture was centrifuged and the diesel washed with acetonitrile during 10 min. The catalyst was recycled for three cycles by recovering the catalyst, washing with ethanol and drying in a desiccator over silica gel overnight to be used in another ODS cycle.
Results and Discussion
The synthesis of the TMA-functionalized silica supports and the preparation of the PW11-based composites is represented in Scheme 1. The synthesis of ethylene-bridged PMO has been performed by co-condensation of the bridged bis-silane [1,2-bis (triethoxysilyl) ethane; BTEE] and the terminal silane (TMA) in the presence of the micelles of the surfactant (Pluronic P123) in acid medium. In the end, an ethanol extraction is performed to remove the surfactant. The formation of SBA-15 is also based in the same surfactant using TEOS as silica precursor in acid medium. The surfactant is removed by calcination after an aging period. The functionalization of SBA-15 was performed via post-grafting with the same terminal silane (TMA; Wahab et al., 2005; Van Der Voort et al., 2013).
An impregnation method was used to immobilize the PW11 in TMA-functionalized silicas by electrostatic interactions. Several characterization techniques were used to assess the successful preparation of the materials.
The FT-IR spectra of the prepared TMA-functionalized supports and the resulting PW11 composites are displayed in Figure 1. The spectrum of PW11@TMA-SBA-15 presents a similar profile to that of TMA-SBA-15. The typical bands assigned to the siliceous support located at 1,100-400 cm−1 range namely the νas(Si–O–Si), νs(Si–O–Si), and δ(O–Si–O) vibrational modes, respectively (Pires et al., 2016; Mirante et al., 2017; Ribeiro et al., 2017, 2018a), mask the bands that could be assigned to the POM incorporation and no extra bands can be recognized. Nevertheless, the intense bands observed in the 1,010–860 cm−1 range of the FT-Raman spectrum (Figure 2A) confirm the presence of the POM in the composite. The bands associated with the vibrational modes of the amine groups can also be observed, in particular, in the 3,035–2,902 and 1,450–1,412 cm−1 ranges for ν(C-H) and δ(CH2), respectively. The presence of POM in this composite material was also confirmed by the presence of W in the EDS spectrum (Figure 3A) and elemental analysis with a loading of PW11 of 0.099 mmol per g of material.
Figure 1. FT-IR spectra of (A) TMA-SBA-15 and (B) TMA-PMOE composites with PW11 before and after catalysis.
Figure 3. SEM images and EDS spectra of the PW11 composites: (A) PW11@TMA-SBA-15 and (B) PW11@TMA-PMOE.
The FT-IR spectrum of the TMA-PMOE material (Figure 1B) presents the typical bands of silsesquioxane frameworks namely the intense band in at 1,093 cm−1 assigned to νas(Si–O–Si), as well as the bands located at 912, 763, and 441 cm−1 associated with the νas(Si–OH), νs(Si–O–Si), and δ(Si–O–Si), respectively (Zhu et al., 2002; Zhang et al., 2005, 2008; Granadeiro et al., 2016). The band at 2,898 cm−1 (stretching) and at 1,413 cm−1 (bending) can be associated to the C–H vibrational modes of bridging-ethylene and the alkyl groups of TMA (Zhang et al., 2005, 2008; Granadeiro et al., 2016). The absence of bands located around 1,340 and 1,380 cm−1 suggest that the Pluronic P123 surfactant has been successfully removed during the extraction process (Li et al., 2007). The FT-IR spectrum of the PW11@TMA-PMOE composite suggests that the structure of the support, previously described, has been maintained. The presence of PW11 in the composite material is not evident in the FT-IR spectrum, due to the presence of intense bands arising from TMA-PMOE. As previously discussed, FT-Raman allows a better identification of the POM vibrational bands due to the lower intensity of the bands from the siliceous supports (Zhang and Yang, 2008; Abdalla and Li, 2012; Julião et al., 2015).
The FT-Raman spectrum of PW11@TMA-PMOE (Figure 2B) is mainly dominated by the bands from PMOE, in particular, the bands in the 3,000–2,800 cm−1 range ascribed to ν(C-H), at 1,412 and 1,271 cm−1 assigned to the twisting and wagging modes of CH2, respectively, and the band at 515 cm−1 assigned to the vibration of the ethylene unit against the siliceous framework together with δ(Si–O–Si) vibrations (Hoffmann et al., 2007; Granadeiro et al., 2016). The presence of PW11 is suggested by the band at 980 cm−1 which can be attributed to the νas(W-Od) vibrations (Ribeiro et al., 2018b). Further suggest the presence of PW11 in the composite material was confirmed by EDS (Figure 3B) with the presence of W in the EDS spectrum, as well as by elemental analysis of W with a loading of 0.057 mmol of POM per 1 g of material.
The SEM image of the TMA-SBA-15 support (Figure S1) presents the typical hexagonal elongated particles of the mesoporous SBA-15 framework. The SEM image of its analog composite (Figure 3A) reveals that the morphology of TMA-SBA-15 was retained after the POM incorporation. The SEM image of PW11@TMA-PMOE (Figure 3B) still exhibits the same morphology as the support (Figure S1) composed by ropelike structures with several micrometers in length (Inagaki et al., 1999; Granadeiro et al., 2016). The chemical composition of the composites was evaluated by EDS. The results reveal, besides silicon as the main element, the presence of tungsten which is consistent with the incorporation of PW11 on the final composites (Figure 3). TEM analysis of both PW11-composites were also performed for a more detailed observation of the morphology (Figure 4). The images show highly ordered and uniform pores organized in long channels as well as the typical honeycomb-shaped hexagonal mesopores. The results obtained from electronic microscopy strongly suggest the structural preservation of the supports in the final composites.
Figure 4. TEM images of the (A,B) PW11@TMA-SBA-15 and (C,D) PW11@TMA-PMOE composites at different magnifications.
The powder XRD patterns of the TMA-functionalized supports and the resulting PW11 composites are presented in Figure 4. The powder XRD of the PW11@TMA-SBA-15 composite (Figure 5A) exhibit the typical low-angle three peaks of SBA-15 materials, with a shift to higher 2θ and with lower intensity for the (110) and (200) reflections, as previously reported in other POM-incorporated SBA-15 composites (Zhao et al., 1998; Kruk et al., 2000; Ribeiro et al., 2013b). The absence of peaks from the PW11 points out that the POM has been successfully incorporated. In the case of the mesoporous organosilicas support (TMA-PMOE), the pattern also exhibits the same three peaks as seen in the SBA-15 materials. The highly ordered mesostructure of the prepared PMO was accomplished by the addition of a KCl as additive which improved the interaction between the organosilica oligomers and the surfactant (Zhang et al., 2008; Granadeiro et al., 2016). The powder XRD pattern of the PW11@TMA-PMOE composite is similar to the TMA-PMOE support, indicating structural preservation of the support.
Figure 5. Powder XRD patterns of (A) TMA-SBA-15 and (B) TMA-PMOE composites with PW11 (ac stands for after catalysis).
The textural properties of SBA-15 and ethylene-bridge PMO materials were evaluated by N2 adsorption experiments. Both materials exhibited type IV isotherms with H1 hysteresis loops typical of mesoporous materials (Figure 6). The results obtained for the starting supports and the PW11-composites are summarized in Table 1. The SBET and Vp values obtained for both composites are smaller when compared with the corresponding support, which confirms the incorporation of POM inside the pore channels (Granadeiro et al., 2016; Mirante et al., 2017; Ribeiro et al., 2018b).
Figure 6. N2 adsorption-desorption isotherms of the TMA-SBA-15 and PW11@TMA-SBA-15 composite (left); TMA-PMOE and PW11@TMA-PMOE (right).
The composites were also studied by 31P MAS-NMR (Figure 7) in order to check the structure of the POM inside the silica materials. The spectrum of the PW11@TMA-SBA-15 composite presents a broad peak centered at −10.41 ppm with a prominent shoulder at −12.75 ppm. The shoulder corresponds to free PW11 anion while the broad peak may be resultant from the interaction of the POM with the siliceous matrix. In fact, a downfield shift in the 31P NMR signal of POMs has been reported as a result of the interaction with Si-OH2+ groups of the silica support (Morales et al., 2018). This result indicates that the PW11 structure was retained after its incorporation in the TMA-SBA-15 support. The spectrum of PW11@TMA-PMOE presents three different peaks with approximately similar intensities at −10.64, −12.79, and −15.16 ppm. The two peaks located at −10.64 and −12.79 ppm should correspond to the free PW11 and PW11 interacting with the support, respectively, as previously discussed. The peak at −15.16 ppm should correspond to the PW11 anion with occupied lacuna which is known to promote an upfield shift of the 31P signal (Guo et al., 2002; Granadeiro et al., 2013).
Figure 7. 31P MAS NMR spectra of PW11, and the corresponding PW11@TMA-SBA-15 and PW11@TMA-PMOE composites.
The TMA-functionalized SBA-15 and ethylene-bridged PMO supports as well as the resulting PW11 composites were analyzed by 13C CP MAS NMR spectroscopy (Figure 8). The spectrum of the TMA-SBA-15 support exhibited four peaks located at 70.33, 55.19, 18.59, and 10.94 ppm (Figure 8A). The peak located at 55.19 ppm correspond to the methyl group and the others can be assigned to the C3 (70.33), C2 (18.59), and C1(10.94) carbon atoms of the TMA group, respectively, Si-1CH2-2CH2-3CH2-N+(CH3)3. The nonexistence of 13C signals in the 67–77 ppm range gives a positive indication that the surfactant has been efficiently removed (Hoffmann et al., 2007). The spectrum of the PW11@TMA-SBA-15 also presents chemical shifts similar to those of the support material, namely at 70.74, 55.16, 18.61, and 10.53 ppm.
Figure 8. 13C CP MAS NMR spectra of the TMA-functionalized supports and the resulting PW11 composites.
The spectrum of TMA-PMOE exhibits an intense peak at 6.63 ppm corresponding to the –CH2-CH2- groups (Figure 8B), and three peaks at 10.75, 18.85, and 70.50 ppm assigned to the C1, C2, and C3 of TMA, respectively. The spectrum, also presents a peak at 61.37 ppm ascribable to CH3-CH2-O groups most likely due to incomplete hydrolysis of 1,2-bis (triethoxysilyl) ethane or during the ethanolic extraction (Zhang et al., 2005, 2008; Granadeiro et al., 2016).
The 29Si MAS NMR solid state spectrum of the TMA-SBA-15 support presents a broad band at δ ≈ −112 ppm corresponding to Q4 species as well as two shoulders at δ ≈ −105 ppm and −93 ppm assigned to Q3 and Q2 species, where Qn = Si(OSi)4−n(OH)n, n = 2–4 (Figure 9A) (Zhang and Yang, 2008; Mauder et al., 2009; Ribeiro et al., 2018b). An identical profile could be observed in the spectrum of PW11@TMA-SBA-15 suggesting that the siliceous framework was preserved in the final composite. The 29Si MAS spectrum of the ethylene-bridge PMO (Figure 9B) presents the characteristic Tn signals attributed to [C–Si(OSi)2(OH)] (T2 at −60.8 ppm) and [C–Si(OSi)3] (T3 at −67.1 ppm) and some Qn signals (Si sites attached to four oxygen atom) between −100 and −115 ppm. This is indicative of some cleavage of the Si-C bond during the synthesis and surfactant extraction process (Granadeiro et al., 2016; Huang et al., 2017). The incorporation of PW11 in the TMA-PMOE support did not result in significant changes in the 29Si MAS spectrum of the composite when compared to the starting support material, which also indicates the preservation of the siliceous structure.
Oxidative Desulfurization Processes Using Model Diesel
A model diesel with ~2,000 ppm S was used for the desulfurization studies which were carried out at 70°C. The PW11@TMA-SBA-15 composite and the PW11@TMA-PMOE composite were tested as heterogeneous catalysts. The oxidative desulfurization of the model diesel was accomplished using ether a biphasic ECODS system or a solvent-free CODS system. For both systems, the equivalent amount of composite containing 3 μmol of PW11 were used.
An initial liquid-liquid extraction was performed (10 min of stirring at 70°C) for the biphasic system in the presence of the catalyst. During this step, some sulfur-containing compounds were removed from the model oil to the solvent phase, until the transfer equilibrium is reached. Afterwards, the oxidative catalytic stage was initiated by adding the oxidant (ratio H2O2/S = 8, at 70°C). In this step, the sulfur compounds were simultaneously oxidized and extracted to the MeCN phase. The solvent-free system, begins with the catalytic stage, in the absence of extraction solvent (ratio H2O2/S = 4, at 70°C), followed by a final extraction step with MeCN or water to remove the oxidized sulfur compounds.
Figure 10 displays the results obtained for the biphasic system using both heterogeneous catalysts: PW11@TMA-SBA-15 and PW11@TMA-PMOE. It can be observed that during the initial extraction step (10 min stirring) between 55 and 60% of the sulfur-containing compounds are transferred into the MeCN phase. This transference follows the order 1-BT > DBT > 4-MDBT > 4,6-DMDBT, with 1-BT being the most easily extracted, which can be justified by the size and geometry of each compound (Ribeiro et al., 2013a, 2016). In the oxidative catalytic step, the non-oxidized sulfur compounds, mostly present in the solvent phase, are oxidized and simultaneously more sulfur compounds are transferred to the solvent phase (Julião et al., 2015). The PW11@TMA-SBA-15 catalyst was able to achieve complete desulfurization for DBT, 4-MDBT and 4,6-DMDBT and 93.9% for 1-BT, after 30 min of the catalytic step initiation. After 60 min, only 2 ppm of 1-BT remained in the model diesel. The oxidative reactivity follows the expected order, with 1-BT being the most difficult compound to be oxidized, as a result of steric hindrance and the different electronic densities of the sulfur atoms (Ribeiro et al., 2013a, 2018b; Julião et al., 2015; Mirante et al., 2017). Regarding the PW11@TMA-PMOE catalyst, the desulfurization efficiency was slightly lower reaching after 60 min of catalytic oxidation, 92.8% for 1-BT, 98.2% for DBT, 99.0% for 4-MDBT and 99.3% for 4,6-DMDBT, resulting in a total desulfurization of 96.9%.
Figure 10. Individual desulfurization percentages (left) and total desulfurization (right) using the biphasic system (1:1 model diesel/MeCN; H2O2/S = 8; 70°C), using the PW11@TMA-SBA-15 and PW11@TMA-PMOE catalysts.
The desulfurization results for the solvent-free system using both composites are presented in Figure 11. During the initial 10 min of catalytic oxidation, the PW11@TMA-PMOE catalyst is faster in the oxidation of sulfur compounds, achieving 72.0% of total oxidation, while the PW11@TMA-SBA-15 reached 45.3%. This might be related to the higher hydrophobicity of the PW11@TMA-PMOE composite that possesses more affinity with the diesel phase than PW11@TMA-SBA-15. Nevertheless, after 30 min of the process, the performance of both catalysts is similar, achieving total conversion for DBT, 4-MDBT and 4,6-DMDBT, and 84.2 and 86.8% conversion for 1-BT, using PW11@TMA-PMOE and PW11@TMA-SBA-15, respectively. At the end of 60 min of desulfurization, the SBA-15 composite reached ultra-low levels of sulfur with only 2 ppm of 1-BT (99.6% conversion) remaining in the model diesel, while with the PMOE composite only 89.6% of 1-BT has been converted (52 ppm remaining). In fact, even continuing the reaction up to 120 min, the conversion of 1-BT is still 93.8% using the PW11@TMA-PMOE. It seems that, despite the initial fast reaction rate (first 30 min), desulfurization gradually becomes slower and total conversion could not be reached even after 120 min. A similar behavior was also observed with the biphasic system since total desulfurization only reached 98.1 and 98.3%, after 120 and 240 min reaction, respectively.
Figure 11. Individual desulfurization percentages (left) and total oxidative desulfurization (right), using the solvent-free system (ratio H2O2/S = 4 at 70°C) and PW11@TMA-SBA-15 and PW11@TMA-PMOE as catalysts (containing 3 μmol of PW11 active center).
Increasing the amount of oxidant in the solvent-free system to a ratio of H2O2/S = 8, ultra-low levels of sulfur could be achieved (7 ppm) with PW11@TMA-PMOE at the end of 60 min. Contrastingly, increasing the amount of oxidant (H2O2/S = 4 to H2O2/S = 8) the desulfurization efficiency of PW11@TMA-SBA-15 shows a slight decrease (Figure 12). In summary, the PW11@TMA-PMOE needs a higher amount of oxidant than the SBA-15 composites, which exhibits better desulfurization efficiencies with a H2O2/S = 4 ratio than H2O2/S = 8 (Ribeiro et al., 2018b).
Figure 12. Total conversion for sulfur oxidation presented in the model diesel, using the solvent-free system at 70°C and PW11@TMA-SBA-15 catalyst in the presence of two different H2O2/S ratios.
Recyclability of PW11@TMA-SBA-15
The recycling ability of PW11@TMA-SBA-15 was evaluated for several consecutive cycles using both desulfurization systems (biphasic and solvent-free). The desulfurization percentages obtained for eight cycles with the PW11@TMA-SBA-15 catalyst after 60 min of oxidation using the biphasic or the solvent-free systems are displayed in Figure 13. The results show that both systems maintain the desulfurization performance along six consecutive cycles. After the sixth cycle, some loss of catalytic activity was detected probably due to the active site deactivation by the presence of sulfones strongly adsorbed in the catalyst surface (Mirante et al., 2017, 2018).
Figure 13. Recycling studies with PW11@TMA-SBA-15 composite after 60 min of the oxidative step under the solvent-free (H2O2/S = 4) and biphasic (H2O2/S = 8) systems at 70°C.
The recovered PW11 composites were studied by several characterization techniques to access their robustness after catalytic use. The ICP-OES of the PW11@TMA-SBA-15 composite suggest that pratically no leaching of PW11 occurs using the solvent-free system, since the Si/W (molar) ratios before (1.46) and after catalysis (1.44 after eight desulfurization cycles) are very similar. The analysis performed after one desulfurization cycle, using the biphasic system, detected some leaching of the PW11 from the support material shown by the increase in the Si/W (molar) ratio from 1.46 to 1.79 after catalytic use. The ICP results for PW11@TMA-PMOE revealed that, after its use in the solvent-free system, the Si/W ratio was maintained (1.01 before and 1.00 after catalytic use).
The powder XRD patterns of PW11@TMA-SBA-15-ac and PW11@TMA-PMOE-ac after use in the biphasic system, exhibit similar profiles to the patterns of the as-prepared materials (Figure 5), which shows that both composites are robust and stable catalysts maintaining the structure of the siliceous support.
Regarding FT-IR spectra, the main vibrational bands of the composites after catalytic use remain practically unchanged, which is consistent with its structural retention (Figure 1). The morphology of the composites also seems to have been preserved after the desulfurization process as observed by SEM and TEM (Figures S2, S3), while the corresponding EDS spectra of each composite confirm the presence of PW11 by still exhibiting tungsten. Moreover, the EDS spectra of the PW11@TMA-SBA-15-ac also reveals the presence of S, resulting from oxidized sulfur-containing compounds from the model diesel that precipitated in the presence of the composite (Figure S2).
The 31P MAS NMR spectra of the composites before and after catalysis, using both desulfurization systems, have also been compared (Figure S4). The spectrum of PW11@TMA-SBA-15 after a biphasic cycle presents a broad peak at −12.68 ppm with a shoulder at −10.92 ppm. The results means that, after the biphasic cycle, the predominant species in the material is now the PW11 interacting with the support rather than the free PW11 as in the as-prepared composite. After eight cycles under the solvent-free system, the spectrum displays a main peak at −15.19 ppm that can be assigned to the PW11 with a saturated lacuna, probably by sulfone or peroxo interactions.
The PW11@TMA-PMOE after one biphasic desulfurization cycle exhibits a broad single peak centered at −10.72 ppm, while after one solvent-free cycle, an additional peak is also observed at −12.27 ppm. These peaks correspond to free PW11 (−10.64 ppm) and PW11 interacting with the silica support (−12.27 ppm), showing the structural retention of these species from the as-prepared material.
Oxidative Desulfurization Process Using Untreated Diesel
The PW11@TMA-SBA-15 was selected for application in the desulfurization of an untreated diesel supplied by CEPSA (1335 ppm S). Figure S5 presents the GC-FID/SCD analysis of this diesel. These results reveal that dibenzothiophenes and dibenzothiophene derivatives are the main sulfur compounds present in this untreated diesel. These studies were performed using the biphasic and also the solvent-free conditions at 70°C with H2O2/S = 8 during 2 h of oxidation. After the oxidative catalytic process, a liquid-liquid extraction with 1:1 MeCN/diesel at room temperature was performed for all treated diesel samples, during 10 min with stirring. At the end of the first cycle, the best performance was obtained with the biphasic system (93.1%), performing an initial extraction (1:1 diesel/MeCN) before oxidation step and also a final extraction after oxidation (Figure 14). On the other hand, under the solvent-free system, a lower desulfurization efficiency was reached (75%).
Figure 14. Desulfurization results of a real untreated diesel obtained after 2 h, catalyzed by PW11@TMA-SBA-15 at 70°C, using the solvent-free biphasic systems with H2O2/S = 8.
The solid PW11@TMA-SBA-15 catalyst was further recycled for two more consecutive cycles following the previously described procedure for simulant diesel. The biphasic system showed higher recycling ability than the solvent-free system, since catalyst efficiency is maintained for three consecutive cycles. On the other hand, an increase of desulfurization efficiency is observed from the second to the third cycle in the solvent-free system. This must be related to the formation of catalytic active intermediates during the previous desulfurization cycles. These active intermediates can be attributed to active peroxo compounds.
In summary, the PW11@TMA-SBA-15 using the biphasic system revealed to be a promising process for the desulfurization of untreated diesel reaching a desulfurization of 93.1% from CEPSA diesel. Moreover, the catalyst maintained its catalytic activity over three consecutive desulfurization cycles.
In this work, two different POM-based silica composites were prepared via impregnation of PW11 on TMA-functionalized mesoporous silica (SBA-15 and ethylene-bridged PMO). The cationic functional group promotes the immobilization of the anionic PW11 by electrostatic interaction. The surface modification of SBA-15 was accomplished by post-synthetic grafting, while the introduction of functional groups in PMO was achieved in situ (“co-condensation”).
The desulfurization performance of PW11@TMA-SBA-15 and PW11@TMA-PMOE composites was evaluated using a multicomponent model diesel. The PW11@TMA-PMOE composite achieved 96.9% of desulfurization after 60 min of oxidative reaction, while the PW11@TMA-SBA-15 allowed to reach ultra-low levels of sulfur (<10 ppm) under the same period of time, using the biphasic system.
Complete conversions of DBT, 4-MDBT, and 4,6-DMDBT could be achieved for both catalysts after only 30 min under the solvent-free system. The PW11@TMA-PMOE catalyst revealed to be slightly less active than the PW11@TMA-SBA-15 catalyst, since the complete desulfurization of 1-BT could not be reached, after 60 min of the process. In comparison, the PW11@TMA-SBA-15 reached ultra-low levels of sulfur (2 ppm) for the same period of time. Moreover, the PW11@TMA-SBA-15 has shown a remarkable recycling ability, in both desulfurization systems, by maintaining its catalytic efficiency for six consecutive cycles.
The robustness of the composites was confirmed by characterization studies of the recovered solid catalysts suggesting their structural and chemical preservation after catalytic use.
The remarkable performance of PW11@TMA-SBA-15 with simulant diesel has led to its application with an untreated diesel, under biphasic and solvent-free systems. Furthermore, recycling tests were also performed using both systems for three consecutive cycles. The best result was obtained using the biphasic system removing 93.1% of sulfur compounds from the diesel after only 2 h and maintaining its remarkable desulfurization efficiency for two additional cycles. These successful recycling studies using a real untreated diesel, strongly indicates that the PW11@TMA-SBA-15 catalyst under biphasic conditions can be a promising system for reaching sulfur-free fuels.
Data Availability Statement
All datasets generated for this study are included in the article/Supplementary Material.
SR (supervised by BC and SB) performed all the experimental work related with the preparation of the materials and the oxidative desulfurization reactions. CG contributed with the discussion of the characterization data and manuscript preparation. MC conducted the MAS NMR experiments. JP performed the textural characterization of the samples. JC-M co-supervised the desulfurization experiments with real diesel.
This work was partly funded through the projects REQUIMTE-LAQV [FCT (Fundação para a Ciência e a Tecnologia) Ref. LAQV, REQUIMTE (POCI-01-0145-FEDER-007265, UID/QUI/50006/2019)], and CENIMAT/I3N (POCI-01-0145-FEDER-007688, UID/CTM/50025/2013), financed by national funds through the FCT/MEC and when appropriate co-financed by the Fundo Europeu de Desenvolvimento Regional (FEDER) under the PT2020 Partnership Agreement. This work was funded by national funds (OE), through FCT—Fundação para a Ciência e a Tecnologia, I.P., in the scope of the framework contract foreseen in the numbers 4, 5, and 6 of the article 23, of the Decree-Law 57/2016, of August 29, changed by Law 57/2017, of July 19. The FCT and the European Union are acknowledged for the post-graduation grant SFRH/BD/95571/2013 (to SR) co-funded by Ministério da Ciência, Tecnologia e Ensino Superior (MCTES), and the European Social Fund through the program POPH of QREN. The authors also acknowledge the Portuguese Nuclear Magnetic Resonance Network (PTNMR—ROTEIRO/0031/2013, PINFRA/22161/2016) and the research group of Prof. Isabel Gonçalves from CICECO Laboratory, University of Aveiro, Portugal.
Conflict of Interest
The 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.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fchem.2019.00756/full#supplementary-material
Abdalla, Z. E. A., and Li, B. (2012). Preparation of MCM-41 supported (Bu4N)4H3(PW11O39) catalyst and its performance in oxidative desulfurization. Chem. Eng. J. 200–202, 113–121. doi: 10.1016/j.cej2012.06.004
Abdalla, Z. E. A., Li, B., and Tufail, A. (2009). Direct synthesis of mesoporous (C19H42N)4H3(PW11O39)/SiO2 and its catalytic performance in oxidative desulfurization. Colloids Surf. A 341, 86–92. doi: 10.1016/j.colsurfa.2009.03.042
Casuscelli, S. G., Crivello, M. E., Perez, C. F., Ghione, G., Herrero, E. R., Pizzio, L. R., et al. (2004). Effect of reaction conditions on limonene epoxidation with H2O2 catalyzed by supported Keggin heteropolycompounds. Appl. Catal. A 274, 115–122. doi: 10.1016/j.apcata.2004.05.043
Esquivel, D., Van Der Voort, P., and Romero-Salguero, F. J. (2014). Designing advanced functional periodic mesoporous organosilicas for biomedical applications. Aims Mater. Sci. 1, 70–86. doi: 10.3934/matersci.2014.1.70
Granadeiro, C. M., Barbosa, A. D. S., Silva, P., Paz, F. A. A., Saini, V. K., Pires, J., et al. (2013). Monovacant polyoxometalates incorporated into MIL-101(Cr): novel heterogeneous catalysts for liquid phase oxidation. Appl. Catal. A 453, 316–326. doi: 10.1016/j.apcata.2012.12.039
Granadeiro, C. M., Ribeiro, S. O., Kaczmarek, A. M., Cunha-Silva, L., Almeida, P. L., Gago, S., et al. (2016). A novel red emitting material based on polyoxometalate@periodic mesoporous organosilica. Micropor. Mesopor. Mater. 234, 248–256. doi: 10.1016/j.micromeso.2016.07.023
Guo, Y., Yang, Y., Hu, C., Guo, C., Wang, E., Zou, Y., et al. (2002). Preparation, characterization and photochemical properties of ordered macroporous hybrid silica materials based on monovacant Keggin-type polyoxometalates. J. Mater. Chem. 12, 3046–3052. doi: 10.1039/b110873p
Hoffmann, F., Güngerich, M., Klar, P. J., and Fröba, M. (2007). Vibrational Spectroscopy of Periodic Mesoporous Organosilicas (PMOs) and their precursors: a closer look. J. Phys. Chem. C 111, 5648–5660. doi: 10.1021/jp0668596
Huang, X., Li, W., Wang, M., Tan, X., Wang, Q., Wang, C., et al. (2017). A facile template route to periodic mesoporous organosilicas nanospheres with tubular structure by using compressed CO2. Sci. Rep. 7:45055. doi: 10.1038/srep45055
Inagaki, S., Guan, S., Fukushima, Y., Ohsuna, T., and Terasaki, O. (1999). Novel mesoporous materials with a uniform distribution of organic groups and inorganic oxide in their frameworks. J. Am. Chem. Soc. 121, 9611–9614. doi: 10.1021/ja9916658
Julião, D., Gomes, A. C., Pillinger, M., Cunha-Silva, L., de Castro, B., Gonçalves, I. S., et al. (2015). Desulfurization of model diesel by extraction/oxidation using a zinc-substituted polyoxometalate as catalyst under homogeneous and heterogeneous (MIL-101(Cr) encapsulated) conditions. Fuel Process. Technol. 131, 78–86. doi: 10.1016/j.fuproc.2014.10.030
Julião, D., Valença, R., Ribeiro, J. C., de Castro, B., and Balula, S. S. (2017). Efficient eco-sustainable ionic liquid-polyoxometalate desulfurization processes for model and real diesel. Appl. Catal. A 537, 93–99. doi: 10.1016/j.apcata.2017.02.021
Karimi, B., Mirzaei, H. M., and Mobaraki, A. (2012). Periodic mesoporous organosilica functionalized sulfonic acids as highly efficient and recyclable catalysts in biodiesel production. Catal. Sci. Technol. 2, 828–834. doi: 10.1039/c2cy00444e
Li, C., Liu, J., Shi, X., Yang, J., and Yang, Q. (2007). Periodic mesoporous organosilicas with 1,4-diethylenebenzene in the mesoporous wall: synthesis, characterization, and bioadsorption properties. J. Phys. Chem. C 111, 10948–10954. doi: 10.1021/jp071093a
Liu, Y.-Y., Leus, K., Sun, Z., Li, X., Depauw, H., Wang, A., et al. (2019). Catalytic oxidative desulfurization of model and real diesel over a molybdenum anchored metal-organic framework. Micropor. Mesopor. Mater. 277, 245–252. doi: 10.1016/j.micromeso.2018.11.004
Mauder, D., Akcakayiran, D., Lesnichin, S. B., Findenegg, G. H., and Shenderovich, I. G. (2009). Acidity of sulfonic and phosphonic acid-functionalized SBA-15 under almost water-free conditions. J. Phys. Chem. C 113, 19185–19192. doi: 10.1021/jp907058y
Mirante, F., Dias, L., Silva, M., Ribeiro, S. O., Corvo, M. C., de Castro, B., et al. (2018). Efficient heterogeneous polyoxometalate-hybrid catalysts for the oxidative desulfurization of fuels. Catal. Commun. 104, 1–8. doi: 10.1016/j.catcom.2017.10.006
Mirante, F., Ribeiro, S. O., de Castro, B., Granadeiro, C. M., and Balula, S. S. (2017). Sustainable desulfurization processes catalyzed by Titanium-Polyoxometalate@TM-SBA-15. Top. Catal. 60, 1140–1150. doi: 10.1007/s11244-017-0801-5
Morales, D. M., Frenzel, R. A., Romanelli, G. P., and Pizzio, L. R. (2018). Synthesis and characterization of nanoparticulate silica with organized multimodal porous structure impregnated with 12-phosphotungstic acid for its use in heterogeneous catalysis. Mol. Catal. doi: 10.1016/j.mcat.2018.10.005. [Epub ahead of print].
Pires, J., Borges, S., Carvalho, A., Pereira, C., Pereira, A. M., Fernandes, C., et al. (2016). Magnetically recyclable mesoporous iron oxide–silica materials for the degradation of acetaminophen in water under mild conditions. Polyhedron 106, 125–131. doi: 10.1016/j.poly.2016.01.007
Ribeiro, S., Granadeiro, C. M., Silva, P., Almeida Paz, F. A., de Biani, F. F., Cunha-Silva, L., et al. (2013b). An efficient oxidative desulfurization process using terbium-polyoxometalate@MIL-101(Cr). Catal. Sci. Technol. 3, 2404–2414. doi: 10.1039/c3cy00287j
Ribeiro, S., Barbosa, A. D. S., Gomes, A. C., Pillinger, M., Gonçalves, I. S., Balula, S. S., et al. (2013a). Catalytic oxidative desulfurization systems based on Keggin phosphotungstate and metal-organic framework MIL-101. Fuel Process. Technol. 116, 350–357. doi: 10.1016/j.fuproc.2013.07.011
Ribeiro, S. O., Duarte, B., de Castro, B., Granadeiro, C. M., and Balula, S. S. (2018a). Improving the catalytic performance of Keggin [PW12O40]3− for oxidative desulfurization: ionic liquids versus SBA-15 composite. Materials 11:1196. doi: 10.3390/ma11071196
Ribeiro, S. O., Granadeiro, C. M., Almeida, P. L., Pires, J., Capel-Sanchez, M. C., Campos-Martin, J. M., et al. (2018b). Oxidative desulfurization strategies using Keggin-type polyoxometalate catalysts: biphasic versus solvent-free systems. Catal. Today 333, 226–236. doi: 10.1016/j.cattod.2018.10.046
Ribeiro, S. O., Juliao, D., Cunha-Silva, L., Domingues, V. F., Valenca, R., Ribeiro, J. C., et al. (2016). Catalytic oxidative/extractive desulfurization of model and untreated diesel using hybrid based zinc-substituted polyoxometalates. Fuel 166, 268–275. doi: 10.1016/j.fuel.2015.10.095
Ribeiro, S. O., Nogueira, L. S., Gago, S., Almeida, P. L., Corvo, M. C., de Castro, B., et al. (2017). Desulfurization process conciliating heterogeneous oxidation and liquid extraction: organic solvent or centrifugation/water? Appl. Catal. A 542, 359–367. doi: 10.1016/j.apcata.2017.05.032
Singh, S., Patel, A., and Prakashan, P. (2015). One pot oxidative esterification of aldehyde over recyclable cesium salt of nickel substituted phosphotungstate. Appl. Catal. A 505, 131–140. doi: 10.1016/j.apcata.2015.07.032
Van Der Voort, P., Esquivel, D., De Canck, E., Goethals, F., Van Driessche, I., and Romero-Salguero, F. J. (2013). Periodic Mesoporous Organosilicas: from simple to complex bridges; a comprehensive overview of functions, morphologies and applications. Chem. Soc. Rev. 42, 3913–3955. doi: 10.1039/C2CS35222B
Wahab, M. A., Imae, I., Kawakami, Y., and Ha, C.-S. (2005). Periodic mesoporous organosilica materials incorporating various organic functional groups: synthesis, structural characterization, and morphology. Chem. Mater. 17, 2165–2174. doi: 10.1021/cm0480059
Wang, D., Liu, N., Zhang, J., Zhao, X., Zhang, W., and Zhang, M. (2014). Oxidative desulfurization using ordered mesoporous silicas as catalysts. J. Mol. Catal. A 393, 47–55. doi: 10.1016/j.molcata.2014.05.026
Wu, N., Li, B., Liu, Z., and Han, C. (2014). Synthesis of Keggin-type lacunary 11-tungstophosphates encapsulated into mesoporous silica pillared in clay interlayer galleries and their catalytic performance in oxidative desulfurization. Catal. Commun. 46, 156–160. doi: 10.1016/j.catcom2013.12.013
Yang, H., Jiang, B., Sun, Y., Zhang, L., Huang, Z., Sun, Z., et al. (2017). Heterogeneous oxidative desulfurization of diesel fuel catalyzed by mesoporous polyoxometallate-based polymeric hybrid. J. Hazard. Mater. 333, 63–72. doi: 10.1016/j.jhazmat.2017.03.017
Zhang, L., Liu, J., Yang, J., Yang, Q., and Li, C. (2008). Direct synthesis of highly ordered amine-functionalized mesoporous ethane-silicas. Micropor. Mesopor. Mater. 109, 172–183. doi: 10.1016/j.micromeso.2007.04.050
Zhang, L., Yang, Q., Zhang, W.-H., Li, Y., Yang, J., Jiang, D., et al. (2005). Highly ordered periodic mesoporous ethanesilica synthesized under neutral conditions. J. Mater. Chem. 15, 2562–2568. doi: 10.1039/b504030b
Zhao, D., Feng, J., Huo, Q., Melosh, N., Fredrickson, G. H., Chmelka, B. F., et al. (1998). Triblock copolymer syntheses of mesoporous silica with periodic 50 to 300 angstrom pores. Science 279, 548–552. doi: 10.1126/science.279.5350.548
Zhu, H., Jones, D. J., Zajac, J., Dutartre, R., Rhomari, M., and Rozière, J. (2002). Synthesis of periodic large mesoporous organosilicas and functionalization by incorporation of ligands into the framework wall. Chem. Mater. 14, 4886–4894. doi: 10.1021/cm011742
Keywords: oxidative desulfurization, heterogeneous catalysis, mesoporous silica, polyoxometalate (POM), fuels
Citation: Ribeiro SO, Granadeiro CM, Corvo MC, Pires J, Campos-Martin JM, de Castro B and Balula SS (2019) Mesoporous Silica vs. Organosilica Composites to Desulfurize Diesel. Front. Chem. 7:756. doi: 10.3389/fchem.2019.00756
Received: 03 July 2019; Accepted: 22 October 2019;
Published: 14 November 2019.
Edited by:Pascal Granger, Lille University of Science and Technology, France
Reviewed by:Guohua Liu, Shanghai Normal University, China
Pavel Nikulshin, All-Russia Research Institute of Oil Refining, Russia
Copyright © 2019 Ribeiro, Granadeiro, Corvo, Pires, Campos-Martin, de Castro and Balula. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Salete S. Balula, firstname.lastname@example.org