Huimin
Liu
*a,
Shaoyuan
Sun
a,
Dezheng
Li
a and
Yiming
Lei
*b
aSchool of Chemical and Environmental Engineering, Liaoning University of Technology, Jinzhou, 121001, Liaoning Province, P. R. China. E-mail: liuhuimin08@tsinghua.org.cn
bDepartament de Química (Unitat de Química Inorgànica), Facultat de Ciències, Universitat Autònoma de Barcelona (UAB), Cerdanyola del Valles, 08193, Barcelona, Spain. E-mail: yiming.lei@uab.cat
First published on 24th June 2024
O2-Assisted oxidative dehydrogenation of propane (O2-ODHP) could convert abundant shale gas into propylene as an important chemical raw material, meaning O2-ODHP has practical significance. Thermodynamically, high temperature is beneficial for O2-ODHP; however, high reaction temperature always causes the overoxidation of propylene, leading to a decline in its selectivity. In this regard, it is crucial to achieve low temperatures while maintaining high efficiency and selectivity during O2-ODHP. The use of catalytic technology provides more opportunities for achieving high-efficiency O2-ODHP under mild conditions. Up to now, many kinds of catalytic systems have been elaborately designed, including transition metal oxide catalysts (such as vanadium-based catalysts, molybdenum-based catalysts, etc.), transition metal-based catalysts (such as Pt nanoclusters), rare earth metal oxide catalysts (especially CeO2 related catalysts), and non-metallic catalysts (BN, other B-containing catalysts, and C-based catalysts). In this review, we have summarized the development progress of mainstream catalysts in O2-ODHP, aiming at providing a clear picture to the catalysis community and advancing this research field further.
The emergence of the catalytic energy conversion technique is one of the reliable strategies,9–13 which can achieve the conversion of shale gas into propylene. Propane can undergo thermal cracking to release hydrogen and produce propylene under anaerobic conditions (direct dehydrogenation of propane), eqn (1).14,15 With the assistance of suitable catalysts, at a reaction temperature of 590–630 °C, the one-way conversion rate of propane can reach 33–44%, and the selectivity of propylene is above 80%.14,15 Although catalytic propane conversion is very promising and attractive from a thermodynamic view, the deep dehydrogenation of propane tends to form carbon deposition, which in turn leads to catalyst deactivation.16,17
Owing to the above problem, the O2-assisted oxidative dehydrogenation process using propane as raw material and O2 as a strong oxidant (oxidative dehydrogenation of propane (O2-ODHP), eqn (2)) has attracted significant attention and has been considered as another effective process for producing propylene due to its remarkable cost advantage.18–21 ODHP is a reversible and strongly endothermic reaction. Theoretically, the equilibrium constant of O2-ODHP should increase with the increase of reaction temperature, but practically, its equilibrium constant still remains low at high temperatures. Therefore, from a thermodynamic perspective, O2-ODHP needs to be carried out at high temperatures and low pressure. However, the high reaction temperature conditions intensify the overoxidation of propylene, leading to a decrease in the selectivity of desired propylene. In this case, the fabrication of well-designed catalysts with tunable selectivity and high activity is important to directionally obtain the target product propylene from O2-ODHP at low temperatures.
C3H8 → C3H6 + H2 (+Coke) | (1) |
C3H8 + 1/2O2 → C3H6 + H2O | (2) |
To date, a great number of catalysts have been developed for O2-ODHP (Scheme 1), including transition metal oxide catalysts (such as vanadium-based catalysts, molybdenum-based catalysts, etc.), transition metal-based catalysts (such as Pt nanoclusters), rare earth metal oxide catalysts (especially CeO2 related catalysts) and non-metallic catalysts (BN, other B-containing catalysts and C-based catalysts). Herein, the recent progress of these catalysts in O2-ODHP is summarized and the underlying mechanism is revealed, with the aim to pave pathways for the rational design of more efficient catalysts for ODHP and advance this research field further.
In the conventional ODHP reaction, the conversion of propane into C3H8 is endothermic and volume-increasing, meaning the conditions with high temperatures and low pressures are beneficial for achieving a high propane conversion rate. In practical industrial systems, the ODHP process is always conducted at 550–700 °C to achieve appreciable propylene yields.25 However, the high temperature leads to side or competitive reactions, forming undesired C2H4 and CH4 products due to C–C cracking. Meanwhile, coke is also one of the products from these side reactions. Therefore, it is necessary to consider the ability to dissociate C–H and avoid C–C bond break during the catalyst design steps.26
In order to selectively oxidize C3H8 to C3H6 and H2/H2O (eqn (1) and (2)), O2 as an oxidant has several advantages including low reaction temperatures, resisting coke formation, and extending catalyst longevity.27,28 The main challenge, as mentioned above, is the overoxidation reactions, which oxidize C3H8 to CO and useless CO2. In this review, the representative and recent works on catalyst design for ODHP using O2 as the oxidant will be focused on, discussing the feasible strategies that can improve selectivity and stability under mild conditions.
Although some excellent reviews have summarized the catalytic systems for ODHP, they have focused on the relationship between ODHP performance and catalyst structures.29–31 Other reviews have paid much attention to ODHP with CO2 assistance32 or the comparison of direct and CO2-assisted ODHP (CO2-ODHP).26 Thus, the regulation of O2-ODHP selectivity under mild conditions has not received enough attention. Considering the fact that ODHP technology could convert abundant shale gas into propylene, which is an important chemical raw material, and thus has practical significance, the representative works about utilizing catalyst design to improve the O2-ODHP selectivity should be systematically reviewed to encourage more in-depth research in this area in the future.
The mainstream mechanisms for CO2-assisted ODHP are a direct route with a redox mechanism and an indirect route combining ODHP and RWGS reactions.26 For the former, C3H8 interacts with O species to form C3H6 and water, and generates O vacancies over the catalyst surface. Then, CO2 refills the O vacancies, supplying active O species and producing CO simultaneously. The direct route can be completed at a single active site or a dual site. Regarding the indirect route, ODHP and RWGS reactions could be completed to realize the simultaneous conversion of CO2 and C3H8. During this pathway, the coke should be rapidly removed via the reverse Boudouard reaction.
Many research results about CO2-ODHP have been published. For instance, Furukawa et al. reported a Pt–Co–In ternary nanoalloy on CeO2 with (Pt1−xCox)2In3 pseudo-binary alloy structure.34 The synergistic effect of alloying platinum with indium and cobalt leads to high activity involving C3H6 selectivity and CO2 reduction ability. The existence of Co species offers a high density of states near the Fermi level, reducing the energy barrier of CO2 reduction. The cooperation between the alloy with CO2 activation ability and CeO2 with oxygen releasing ability facilitates coke combustion, thereby enhancing the catalyst stability. Song's group reported an efficient bifunctional CO2-ODHP site capable of activating C–H and CO bonds with inhibited C–C scission capacity.35 In detail, their work proposed a La-modified binuclear Feoxo (La-Feoxo) site stabilized on a silicalite-1 support, achieving 32.7% C3H8 conversion with 83.2% C3H6 selectivity. The improvement should be attributed to the balance in the acid–base property due to FeOx dispersion and La modification. The elaborate active site could stabilize HCOO* intermediates toward enhanced CO2-assisted H removal efficiency.
At present, the main problem during CO2-ODHP is the dry reforming of propane (DRP) competing reaction (C3H8 + 3CO2 → 6CO + 4H2). The DRP reaction could be the dominant reaction under high-temperature conditions, converting too much C3H8 into CO instead of C3H6.36 Therefore, to avoid the undesired DRP process, the C–H/C–C scission ability of C3H8 should be mainly focused on during the development of CO2-ODHP catalysts.37 A few works have summarized and discussed the soft oxidant-assisted ODHP system from three aspects: reaction mechanisms, catalyst composition, and oxidant footprint analysis. Interested readers are recommended to read ref. 29 and 34.
In order to overcome the above disadvantages, chemical looping ODHP (CL-ODHP) combining dehydrogenation of propane and catalyst regeneration is explored and used in industries.43 In a CL-ODHP system, the reaction routes include (i) the conversion of propane into propylene and water via active oxygen species (oxygen carrier) from lattice oxygen (C3H8 + MOx → C3H6 + MOx−1 + H2O); (ii) transfer of the reduced oxygen carrier to the air reactor, where the missing lattice oxygen of the active phase is re-supplied by the air (MOx−1 + 1/2O2 → MOx); and (iii) transfer of the re-oxidized oxygen carriers to ODHP reactors to complete the chemical loop. Therefore, through repeated reduction and oxidation of active species, the oxygen can migrate from the air reactor to the ODHP reactor and take part in the CL-ODHP process. The advantages of CL-ODHP include (i) the decline in energy consumption by H2 oxidation; (ii) reduced investment and operation cost due to the elimination of separation equipment; (iii) avoidance of the potential explosion of C3H8; and (iv) reduced CO2 emissions by indirect flameless combustion of H2.44,45
The development of a suitable catalyst for CL-ODHP is crucial, but compared to the other two ODHP reaction systems, the investigations of catalyst design in the CL-ODHP field are much less. For instance, Gong's group designed multifunctional ferric vanadate–vanadium oxide (FeVO4–VOx) redox catalysts by combining chemical looping-selective H2 combustion and ODHP. VOx provided dehydrogenation sites, and FeVO4 acted as an oxygen carrier for H2 combustion. After 200 chemical looping cycles for the re-oxidation of FeVO4, FeVO4–VOx could maintain an integral performance of 81.3% propylene selectivity at 42.7% propane conversion. More importantly, the authors proposed a hydrogen spillover-mediated coupling mechanism. They believed that the H species generated at VOx sites could transfer to adjacent FeVO4 for combustion, leading to the formation of C3H6. The author's group also reported a MoO3–Fe2O3 catalyst with a similar CL-ODHP mechanism.46 Through introducing Mo species onto Fe2O3, the MoO3–Fe2O3 catalyst achieved very stable catalytic behaviour after 300 cycles with ∼49% propane conversion and ∼90% propylene selectivity. Several works have comprehensively discussed the recent development of the CL-ODHP system,39,47,48 and therefore the subsequent sections will mainly focus on O2-assisted ODHP, which has not been systematically summarized in recent years.
![]() | ||
Fig. 1 (a) Variation of selectivity to propylene (■), ethylene (●), and COx (▲) with vanadium surface densities. Reaction temperature = 600 °C; propane conversion = 40%.52 Copyright 2004 Elsevier. (b) N2 adsorption–desorption isotherms of nV-MCM-41 samples. (c) NH3-TPD profiles of nV-MCM-41 samples.55 Copyright 2018 Elsevier B.V. (d) Propane conversion as a function of temperature on V-MCF catalysts. (Reaction conditions: C3H8![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
For regulating the dispersion of V active sites, Guo's group designed mesoporous nV-MCM-41 molecular sieves via surface engineering.55 Compared to pristine V2O5, at low concentrations in nV-MCM-41, the V species existed in an amorphous state and became incorporated into the silicate framework or dispersed on mesoporous channels. The N2 adsorption–desorption isotherms showed typical type-IV isotherms, suggesting the maintenance of the mesoporous structure after the addition of V species (Fig. 1b). The NH3-TPD characterization was carried out to investigate the acidic characteristics. The results suggested that the isolated tetrahedral-like coordinated VOx species had the properties of Lewis acid sites, which could be adjusted via changing the concentration of V elements (Fig. 1c). The optimal propane conversion rate could reach 58% accompanied by a propene selectivity of ∼92% on 6.8 V-MCM-41. The authors believed that the improvement on catalytic efficiency and selectivity should be ascribed to the large specific surface of MCM-41 with a highly ordered mesoporous structure. But the increasing number of acidic sites reduced propene release from the V active sites, leading to deep oxidation and carbon deposition. This work also confirmed that suitably controlling the number and dispersion of VOx sites on the substrate is a crucial factor for achieving efficient and selective O2-ODHP reaction.
Another example of regulating activity and selectivity via utilizing the morphology and specific absorption sites on supports was reported by Fan et al.54 Fan's group loaded vanadium on MCF at a loading amount of 4.2% and evaluated the catalytic performance in O2-ODHP. Under the reaction temperature of 550 °C and atmospheric pressure conditions, the conversion of propane reached 40.8%, the selectivity of propylene was 68.5%, and the one-way yield was as high as 27.9% (Fig. 1d). The authors suggested that the larger specific surface area and the uniform pore size of the support MCF were beneficial for the rational dispersion of catalytic active sites (VOx, Fig. 1e). Meanwhile, the weak acidity of MCF also promoted the selectivity of propylene.
Considering that the catalytic performance of supported VOx is often influenced by its two- or three-dimensional dispersion, Grant et al. used an SiO2 inert support (im-SiO2) with low reactivity of surface SiOH groups toward the precursor to avoid undesired three-dimensional nanoparticles, which were often formed on conventional supports like Al2O3, TiO2, Nb2O5, or ZrO2.58 Raman spectra suggested that im-SiO2 did not show any sharp signal at 995 cm−1 from 3D V2O5, meaning the enhanced 2D dispersion behavior of V2O5 species (Fig. 1f). Because of the strong structure-sensitivity of O2-ODHP, 2D V2O5 loaded im-SiO2 had faster C3H8 consumption rate and higher TOF. Their results offered new opportunities for SiO2-supported V-based oxides and even other metal oxides applied for O2-ODHP.
Regarding the above vanadium-based catalysts, it is well accepted that their intrinsic activities originate from the strong redox properties of V5+ species. In O2-ODHP, they follow the Mars–van-Krevelen mechanism, where oxygen is in zero order while propane is in first order.60 However, in different catalytic systems, the types of active sites in vanadium-based catalysts might be varied, which include (1) bridged oxygen group (VO); (2) bridged oxygen between vanadium and support (V–O–S); and (3) bridged oxygen between adjacent vanadium atoms (V–O–V). Therefore, making clear the type of active sites in a specific catalytic system and proposing novel reaction mechanisms are the remaining tasks for the future in-depth investigation of vanadium-based catalysts.
In addition to the above controversial active centers, the design of the VOx-based catalyst also faces some challenges. In the O2-ODHP area, extensive studies have been made on VOx-based catalysts.61–63 It has been demonstrated that the catalytic behavior over VOx is related to a series of parameters such as the oxidation state, coordination number, aggregation state, and reducibility of vanadium species. These parameters are not only affected by the loading amount of V but also by the properties of the carrier and the preparation of the catalyst.64 Moreover, through the Mars–van-Krevelen mechanism, hydrocarbon molecules will react with oxygen species in VOx to generate alkene and H2O, and then reduced VOx should be re-oxidized by gaseous oxygen molecules to regenerate the active VOx sites. In this case, the whole process depends on the redox ability of VOx-based catalysts, which is also important for VOx catalyst design.
![]() | ||
Fig. 2 (a) Dependence of the initial propylene formation rate on Cs–Mo/Zr catalysts with different Cs![]() ![]() ![]() ![]() |
Gong et al. used two strategies including the coupling of surface acid catalysis and selective oxidation from lattice oxygen over MoO3–Fe2O3, simultaneously, to promote propylene production.46 They successfully introduced atomically dispersed Mo species on γ-Fe2O3 to create acid sites via the formation of MoO3. The coupling strategies achieved a stable catalytic performance with 49% propane conversion and 90% propylene selectivity for more than 300 cycles. HR-XPS spectra showed the enrichment of Mo, which was beneficial for propane activation, and the increasing density of medium acid sites. NH3-DRIFTS spectra suggested that the addition of Mo species caused both Lewis acid and Brønsted acid sites (Fig. 2c). The enrichment of Mo and modulation of acid properties accelerated the propane activation. Meanwhile, Mo also regulated the lattice oxygen activity, which led to selective oxidative dehydrogenation instead of over-oxidation in γ-Fe2O3 as confirmed by C3H8-D2-TPSR results (Fig. 2d).
Similarly, Xie's group designed the coordination conditions of Mo atoms over single-crystalline Mo2N and MoN monoliths for enhancing the non-oxidative dehydrogenation of propane to propylene.67 Through simulation, the authors predicted the electron structures, suggesting that the high-density Lewis acid sites were from the top-layer Mo with electron deficiency (Fig. 2e). Experimentally, HS-LEISS analysis further confirmed the atomic termination layer of Mo atoms for both Mo2N and MoN. The electron-deficient status was favorable to receive electrons in C–H, accelerating the reactant dissociation and activation (Fig. 2f). The NH3-TPD test clarified that the number of effective active sites (Lewis acid sites) on porous single crystals was higher compared to that on polycrystals. Taking Mo2N as an example, the authors further calculated the energy barrier for each reaction step (Fig. 2g). The results showed that high-density acid sites achieved the low activation energy with a reasonable potential barrier. Therefore, the porous single-crystalline Mo2N and MoN monoliths reached 11% propane conversion and ≈95–97% propylene selectivity under the mild conditions of 500 °C. The above works emphasized two facts: (i) catalysts containing exposed Mo active sites had excellent activity and tunable selectivity in O2-ODHP and (ii) compared to blindly increasing the number of exposed sites, reasonably designing the effective active sites on the catalyst surface could improve the catalytic activity and selectivity considerably.
Despite extensive studies, the performance of molybdenum-based catalysts in O2-ODHP is far from satisfactory, and even not comparable to their vanadium-based counterparts. In addition, molybdenum-based catalysts in O2-ODHP are structure-sensitive. Sometimes, molybdenum-based catalysts obtained by similar preparation methods might exhibit different catalytic performances and even reaction pathways in O2-ODHP.68 All these disadvantages make it a challenging task to study molybdenum-based catalysts in O2-ODHP.
![]() | ||
Fig. 3 (a) Propane conversion and propylene selectivity, 0.3 g catalyst, T = 600 °C, WHSVC3H8 = 0.4 h−1; (b) propylene formation rate tested under similar initial propane conversion; WHSVC3H8 = 4.5, 3.125 and 0.4 h−1 for Co-acac@S-1, Co/S-1 and Co-edta@S-1; (c) the potential energy surface diagram.70 Copyright 2006 Royal Society of Chemistry. (d) XANES spectra of Co-SIM and Co-AIM + NU-1000 before and after activation, along with the data for cobalt reference samples. (e) Reaction coordinate for propane ODH. ΔH503K in kcal mol−1.71 Copyright 2016 American Chemical Society. (f) Schematic representation of the preparation of NU-1000-supported bimetallic catalysts. The promoter ions are anchored via SIM (pink) and Co ions are anchored via AIM (blue). (g) TOF values at 230 °C, as determined by varying the space velocity of propane.72 Copyright 2017 American Chemical Society. |
Cobalt-based catalysts exhibit high activities in O2-ODHP application, but they still suffer from two main shortages. Firstly, the cobalt-based catalysts generally have strong oxidation ability. This property makes it easy for propane overoxidation, which ultimately leads to low selectivity towards propylene. What's more, due to the variability of cobalt valence states, the stability of cobalt-based catalysts during the O2-ODHP process is poor, which remains an important concern.
![]() | ||
Fig. 4 (a) O2-TPD spectra of nano-NiO-sg and meso-NiO-300 samples.76 Copyright 2010 Elsevier. (b) Pore size distribution of mesoporous NiO support and Ni–Mo–O catalysts [(a) Mo/Ni = 0.05, (b) Mo/Ni = 0.1, (c) Mo/Ni = 0.2, (d) Mo/Ni = 0.3, (e) Mo/Ni = 0.4, (f) Mo/Ni = 0.5]. (c) SEM images and particle size distribution (inset) of mesoporous Ni–Mo–O catalysts: Mo/Ni = 0.05. (d) TEM images of mesoporous Ni–Mo–O catalysts: Mo/Ni = 0.5. (e) H2-TPR profiles of mesoporous NiO and Ni-Mo-O catalysts [(a) Mo/Ni = 0.05, (b) Mo/Ni = 0.1, (c) Mo/Ni = 0.2, (d) Mo/Ni = 0.3, (e) Mo/Ni = 0.4, (f) Mo/Ni = 0.5].79 Copyright 2019 Elsevier B.V. (f) Catalytic data obtained at 300 °C for Ti–Ni–O catalysts as a function of TiO2 content.80 Copyright 2005 Elsevier. (g) O2-TPD spectra of pure and Ce/Nb-doped NiO catalysts (200 mg sample). (h) NH3-TPD spectra of pure and Ce/Nb-doped NiO catalysts (200 mg sample).81 Copyright 2010 Springer. (i) Proposed reaction mechanism for the oxidative dehydrogenation of C3H8 to C3H6 by O2 in the presence of HCl over CeO2-based catalysts.82 Copyright 2018 American Chemical Society. |
Liu's group fabricated ordered mesoporous Ni–Mo–O catalysts via a two-step method involving a hard template and impregnation steps.79 The results of pore size distribution displayed the ordered mesoporous characteristic of Ni–Mo–O catalysts with a low ratio of Mo/Ni (Fig. 4b). Similarly, SEM images revealed a mesoporous structure when the concentration of Mo was low (Fig. 4c). TEM images suggested that the elimination of the mesoporous structure with the higher concentration of Mo should be attributed to the blockage of mesoporous channels due to the impregnation of Mo species (Fig. 4d). H2-TPR characterization showed that the addition of Mo could inhibit the reducibility of the Ni species according to the shift of the characteristic peaks toward the high temperature side (Fig. 4e). With an excellent mesoporous morphology and suitable redox ability, the Ni–Mo–O catalyst with a Mo/Ni ratio of 0.4 shows the highest activity with a propylene yield of 15.6%. The authors believed that the Ni–Mo–O catalyst was favorable for the C–C bond cleavage, which was the reason for the high selectivity of propylene. On the other hand, the decreased reducibility also contributed to the high selectivity because pure NiO with high redox ability and a lot of chemisorbed oxygen species would have strong reaction activity during the first step of OHDP.
Besides, multielement-doped NiO catalysts were also developed for O2-ODHP. For instance, Ce3NbNiO catalysts have been prepared by a sol–gel method.81 Compared to pristine NiO, Ce3NbNiO presents a higher propane conversion with higher propene selectivity at a relatively low temperature (250 °C). The enhancement effects of both Ce and Nb were investigated. O2-TPD measurements confirmed that Ce-doping increased the number of electrophilic oxygen species (Fig. 4g), which were mostly responsible for the deep oxidation reactions, so the propane conversion rate was largely improved. For the role of the Nb element, NH3-TPD proved the formation of weak acid sites (Fig. 4h). Since the higher strength and more acid sites over the catalyst surface could decrease the selectivity to olefin. The weak acid sites induced by Nb-doping should contribute to a decline in surface acidity and further improve the propene selectivity of the Ce3NbNiO catalyst.
Besides doping Ni-based catalysts, introducing Ni into other oxide catalysts is another elemental doping idea to improve the activity and selectivity. For instance, Ni doped CeO2 nanorod catalysts for O2-ODHP have been reported.86 Based on the periodic DFT calculations, the authors proposed a potential strategy to improve catalyst performance towards O2-ODHP. In detail, surface hydroxylation enhanced the propylene selectivity via inhibiting the formation of oxygenates and the formation of acetone. Compared to surface oxygen species, hydroxyls were more beneficial for the H abstractors, leading to an elementary mechanistic transformation in propene formation from MvK to LHHW. Thus, the Ni doping and surface hydroxylation (Lewis base addition) could be a reliable strategy for achieving selective O2-ODHP. Another example of Ni as a foreign doping element was reported by Farin et al.87 This work fabricated mesoporous β-NiMoO4 using Ni as the guest ion. Due to the coordination by the copolymer COO− functions, Ni cations were stable in the main β-phase. The presence of Ni species and the large specific surface from the mesoporous structure maintained a high activity level (propane conversion = 14.1%) and also reached high propylene selectivity (propene selectivity = 72%).
Another work about the NiO/CeO2 catalyst was recently reported by Wang's group.82 The authors investigated the relationship between catalytic behavior and exposed facets of CeO2 nanocrystals. The results showed that the [1 1 0] and [1 0 0] facets on the nanorod had the highest activity, and the nanocube enclosed by [1 0 0] facets possessed the most selectivity for propylene evolution. More importantly, the NiO modification achieved a single-pass propylene yield of 55% and outstanding stability with more than 100 h activity. The DFT and mechanistic studies suggested that the determining factors for enhanced activity and selectivity should be the surface oxygen vacancy and the surface chloride coverage (Fig. 4i). In detail, the peroxide species (O22−), which was formed by adsorption of O2 on surface oxygen vacancies, activated chloride and led to a radical-like active chlorine species. Further, the active chlorine species induced the dissociation of the C–H bond in propane, forming propylene as a major product.
![]() | ||
Fig. 5 (a) Raman spectroscopy of MoVOx-1–4.90 Copyright 2019 Wiley. (b) Normalized UV-Raman spectra (385 nm) of bare P25 and vanadia-loaded samples. The signal from the used CaF2 window is marked with an asterisk. (c) Proposed reaction mechanism for the ODH of propane over vanadia-loaded titania (P25) based on operando and transient spectroscopic analyses.93 Copyright 2023, American Chemical Society. (d) Diagram of primary reaction steps.19 Copyright 2009 Springer. (e) Ni K-edge XANES spectra and (f) Ni K-edge FT-EXAFS spectra of Ni1/CA catalysts and reference samples.107 Copyright 2023 Wiley. |
Recently, Hess et al. designed a VOx/TiO2 catalyst.93 Compared to conventional P25, which always caused carbon deposition during the O2-ODHP process, the active sites was located at VOx instead of TiO2. And UV-Raman spectroscopy showed that vanadia-related features at 858, 936, 991, and 1024 cm−1 were attributed to interface V–O–Ti, bridging V–O–V, the terminal vanadyl bond of V2O5, and the terminal vanadyl bond of amorphous vanadia, respectively (Fig. 5b). The O2-ODHP reaction involved preferentially the VO bonds of dimeric species instead of doubly bridged V–O–V bonds. With a high density of V active sites, the oligomeric V with better reducibility avoided the undesired oxidation process. The operando/transient spectroscopic results clarified that TiO2 could transfer H species from propane to V active sites based on Ti–OH groups on anatase (Fig. 5c). Therefore, the determining step of the C–H bond dissociation could be accelerated, improving the whole O2-ODHP efficiency.
As a typical p-type semiconductor, NiO can adsorb various types of oxygen species on its surface even under mild conditions.94 Thus, it has been confirmed that NiO is able to realize O2-ODHP at low temperatures, suggesting its great potential for the practical industry.95,96 However, there are still some challenges in NiO-based catalytic systems. NiO with poor stability is easily reduced to Ni0, and single NiO is hard to use as a stable catalyst with high catalytic performance.97 Although some NiO-based composite oxides, such as Ni–Ce–O, Ni–Nb–O, and Ni–Ti–O, have been developed for improving stability during the O2-ODHP process, the propylene yield rates are still very low.80,89,98 Therefore, how to simultaneously improve stability and reaction efficiency is an urgent problem when designing NiO-based catalysts. Recently, it has been reported that polyoxometalates possess the ability to activate molecular oxygen at moderate temperatures, so this new family has been developed for selective oxidation of light alkanes.99,100 Introducing polyoxometalates into NiO-based catalysts might be a reliable strategy to overcome the current problems.
Besides the above metallic and alloyed transition metal-based catalysts, single-atom catalysts (SACs) are one of the most recent, revolutionary, and rapidly developing research fields in catalytic science.103 SACs inherit the advantages of easy separation and good recyclability of supported metal nanostructured catalysts and also combine the characteristics of homogeneous catalysts with highly homogeneous active centers and an adjustable coordination environment.104 The coordination between single atoms and the support involves strong interaction or charge transfer, conferring unique electronic and geometrical structures to individual metal atoms, differing from that of conventional metal NPs and carrying some charge.105 The high activity of valence electrons, the quantum confinement effect, and the sparse quantum level of metal atoms contribute to the maximum surface free energy of the metal species in SACs, which facilitates chemical interactions between single atoms and the support.105 In the O2-ODHP field, some SAC catalysts were also designed and fabricated. For example, Liu's group achieved the loading of PtSn-SA over commercial nano-Al2O3 with 0.1 wt% Pt.106 By combining in situ characterization and theoretical calculations, the Sn1Pt single-atom alloy structure could be confirmed, showing high efficiency with 47.6 mol gPt−1 h−1 of C3H6 evolution rate and more than 40 h of long-term stability. The results of in situ diffuse reflectance infrared Fourier-transform spectroscopy (DRIFT) suggested that the coke deposition could be inhibited via Sn1Pt single-atom alloy mediated reconstruction, contributing to stability enhancement. Besides the noble-metal SACs, Zhang et al. reduced the size of transition metals and synthesized a Ni-based single-atom catalyst for O2-ODHP.107 X-ray absorption spectroscopy (XAS) confirmed the oxidation state (Ni3+) of Ni species on the CA support (Fig. 5c). EXAFS showed the isolated features of the metal Ni species, meaning the successful construction of the Ni-based SAC catalyst (Fig. 5d). Ni-SA supported on calcium aluminate was highly effective in O2-ODHP (Fig. 5e and f), recording a propylene selectivity of 47.3%, a yield of propylene of 14.2%, and good stability for 20 hours of testing, which was 2–3 times better than that of Ni nanoparticles. XPS characterization elucidated the reason for the improved efficiency and selectivity. The more positive electronic properties of the Ni-SA sites were beneficial for the activation of oxygen molecules (propylene desorption), effectively inhibiting further oxidation of propylene to COx by-products.
Although the above works proved that transition metal-based catalysts are efficient for O2-ODHP, the synthesis process of the clusters and single-atom catalysts is quite complex, requiring accurate and strict control, which limits the industrial application of clusters and single-atom catalysts in O2-ODHP. Meanwhile, in some transition metal-based catalysts containing noble metals, their high cost is also another disadvantageous factor, limiting their large-scale applications and leading to higher requirements in their activity, stability, and recyclability.
The performance of pure CeO2 is still not ideal for industrialization, evidenced by the high propane conversion rate but extremely low propylene selectivity, owing to the occurrence of overoxidation. Modifying CeO2 by NiO,82 Cl,110 F,111 or P112 could improve its performance. For instance, Wan et al. demonstrated that F-modified CeO2 could selectively promote the production of propylene.111 After further modification by an appropriate amount of Cs, the catalyst with the composition 3% Cs2O/CeO2–2CeF3 achieved a propylene yield of more than 30.0% at 500 °C. CeF3 played two functions in 3% Cs2O/CeO2–2CeF3 for O2-ODHP: one is that the surface F replaced some lattice oxygen in CeO2, promoting the migration of oxygen species and the activation of surface oxygen species; and another is that the addition of F resulted in a highly dispersed surface CeO2 active phase. The modification of CeO2 by P could also facilitate the dispersion of the surface active phase CeO2, thus improving propylene selectivity.112 As a result, at the reaction temperature of 550 °C, CePO4 recorded a propylene selectivity of 62.0%, much higher than that of CeO2 (31.0% under the same reaction conditions).
Although there are some studies on the application of rare earth CeO2-based catalysts in O2-ODHP, the enhancement effect from the rare earth metal oxides has not been extensively investigated and utilized in the O2-ODHP area. The reaction mechanism is not clear yet. Thus, more efforts should be devoted to this kind of catalytic system to elucidate the possible theories.
The application of BN in O2-ODHP was initiated in 2016 when the Hermans group reported that commercial hexagonal BN and BN nanotubes exhibited much higher catalytic activity and propylene selectivity in O2-ODHP than other catalysts.118 For example, on a commercial hexagonal BN catalyst, at a reaction temperature of 490 °C, the conversion of propane was 14.0%, the selectivity of propylene could reach 79.0%, and there was almost no complete oxidation product CO2 generated. In addition, BN also exhibited very high stability in O2-ODHP. The rich hydroxyl groups at the edge endowed BN with an initial propane conversion rate of 20.6% and an initial propylene selectivity of 80.2% at a reaction temperature of 530 °C in O2-ODHP.119 Besides, its activity and selectivity remained almost unchanged within 300 hours. The excellent catalytic performance of BN has opened up a new path for the selective cleavage of C–H bonds in alkanes and has become a research hotspot. From then on, several approaches have been adopted to improve the performance of BN in O2-ODHP, which include (1) increasing the specific surface area of BN;120 (2) doping BN with foreign elements;121 (3) preparing BN with a local chemical environment regulated by plasma, and so on.122 Here some typical examples are illustrated.
![]() | ||
Fig. 6 (a) FTIR and (b) B 1s high-resolution spectra of BN-M and uBN-M.124 Copyright 2021 American Chemical Society. (c) The XPS B 1s feature of fresh (dotted) and spent (solid) B-containing catalysts.123 Copyright 2016 Science. (d) The proposed redox reaction cycle in the oxidative dehydrogenation of propane over boron nitride.126 Copyright 2018 Elsevier. (e) TG-MS analysis of Mg-BNSs and their spent forms under an O2-rich atmosphere (20% O2 + 80% N2). (f) FT-IR spectra of Mg-BNSs and the spent Mg-BNSs activated in the reaction atmosphere at different temperatures and durations. (g) B 1s XPS spectrum of Mg-BNSs tested for 100 h over ODHP.18 Copyright 2023 Royal Society of Chemistry. |
In order to study the active sites of the BN catalyst in O2-ODHP, understand the reaction pathway of O2-ODHP on the BN catalyst at the molecular/atomic level, and lay a theoretical foundation for the design of an efficient BN catalyst, researchers analyzed the structure of the BN catalyst and explored the mechanism of BN in O2-ODHP. For example, combining kinetic analysis and theoretical calculations, Hermans et al. proposed that the oxygen terminal site on the hexagonal BN armchair edge B–O–O–B
was the catalytic active site,
B–O–O–B
activated the second hydrogen atom of propane, while breaking the O–O bond to form boron hydroxyl groups and nitrogen oxygen radicals.118 The authors further analyzed the surface elements of the catalyst using XPS and found that in O2-ODHP, BOx species gradually formed on the BN surface and tended to stabilize, suggesting that the BOx species were the active site.125 Hermans et al. used 1H–11B coupling nuclear magnetic resonance and X-ray absorption spectroscopy to investigate the structural changes of the BN catalyst before and after the reaction, determining that the active site of the BN catalyst was a three coordinated boron site (B(OH)xO3−x) with a variable number of hydroxyl groups and bridging oxidation groups, rather than the oxygen terminal site on the armchair edge
B–O–O–B
.123 Lu et al. found that BN exhibited excellent activity and olefin selectivity only after undergoing an activity induction period in O2-ODHP. XPS and in situ XRD characterization confirmed that surface BOx species gradually formed and tended to stabilize during the oxidative dehydrogenation reaction (Fig. 6c), showing a consistent trend with catalyst activity, confirming that BOx species were active sites as proposed by Hermans et al.126,127 However, further combining kinetic analysis, Lu et al. claimed that the B–OH dynamically generated on the catalyst surface was the active site of the catalyst (Fig. 6d).
Besides the above-mentioned several species, boron atoms at the defective sites,128 boron-atom-terminated zigzag,129 and N2O or NOx-type sites130 have also been reported as active sites on BN in O2-ODHP. Therefore, currently, there is still controversy over the active site of BN in O2-ODHP. Researchers need to devote more effects to make clear which one is the intrinsic active site for O2-ODHP over a specific catalytic system.
Different C-based catalysts with 2D or 3D morphologies, including carbon nanotubes (CNTs),135 carbon nanofibers,136 nanodiamonds,137etc., have been applied in O2-ODHP. Generally, these C-based catalysts have been modified with the aim of improving their catalytic performance. For instance, Frank et al. modified carbon nanotubes with boron or phosphorus and investigated the effects of heteroatoms in O2-ODHP.135 It was revealed that the introduction of both boron and phosphorus increased the apparent activation energy of the catalyst and the reaction order of oxygen while affecting little the reaction order of propane. It indicated that boron or phosphorus can inhibit the formation of electrophilic oxygen-active sites, where phosphorus was more effective. The inhibited formation of electrophilic oxygen active sites prohibited the overoxidation of propane and propylene and thus improved the selectivity towards propylene. As a result, the modification of carbon nanotubes with phosphorus enhanced the thermal stability of carbon nanotubes, achieving a propane conversion of 7.0% at a selectivity of propylene of about 50.0% (Fig. 7a and b). Chen et al. doped nitrogen into carbon nanotubes to regulate the types of functional groups on their surface, and investigated their performance in O2-ODHP.138 It was demonstrated that nitrogen doping helped to increase the surface charge density of carbon nanotubes, leading to a strengthened repulsive force between carbon nanotubes and propylene. It was beneficial for reducing the adsorption of propylene on the catalyst, avoiding the overoxidation of propylene, and consequently improving the selectivity towards propylene (Fig. 7c and d).
![]() | ||
Fig. 7 (a) Catalytic performance of carbon nanotube catalysts in ODHP, with VOx–Al2O3 as a reference. (b) Propylene selectivity at 5% propane conversion (■) and reaction rate r (○) as a function of B2O3 loading (WB2O3). Reproduced with permission from ref. 135. Copyright 2009 Wiley. (c) Structural parameters and C3H6 formation rates of spent catalysts. (d) Reaction orders, C3H6 formation rates, and activation energy as a function of graphitic N content. Reproduced with permission.138 Copyright 2013 Royal Society of Chemistry. Deconvolution of (e) S 2p and (f) P 2p XPS signals for PZS@OCNT-800 before and after ODHP reactions (450 °C and 520 °C).139 Copyright 2022 Elsevier. (g) and (h) In situ DRIFT difference spectra at steady state after switching various reaction atmospheres at 500 °C.140 Copyright 2020 Elsevier. (i) UV-DRS of SS-BCNNSs. (j) Dependence of (Ahv)2 on the photon energy. (k) Schematic band structure of SS-BCNNSs. (l) DFT calculated energy profile for ODHP under thermal and photo-thermal conditions.141 Copyright 2023 Elsevier. |
Another strategy was reported by Qi et al., who utilized heteroatom co-doped polymers (PZS) and oxidized CNTs, achieving 63% propylene selectivity at 14.3% propane conversion.139 What's more, the poor stability of CNTs has also been improved, realizing more than 20 h O2-ODHP reactions at 520 °C. The kinetic tests conducted at different propane and oxygen partial pressures demonstrated that the as-prepared CNT-based catalysts have excellent resistance to oxidation and carbon deposition, which is the reason for the enhanced stability. The XPS spectra could detect the P–O and S–C species on the catalyst surface, which contributed to the improvement of propylene selectivity and oxidation resistance (Fig. 7e and f).
In spite of the fact that C-based catalysts exhibited a certain activity in O2-ODHP, due to the complex types of functional groups on the surface of C-based materials (such as hydroxyl group, carbonyl group, anhydride group, lactone group, etc.), the qualitative and quantitative determination of effective functional groups has always been a research difficulty in the field of O2-ODHP.142 In addition, the complex variety of functional groups on the surface of C-based catalysts makes the reaction more complex. The oxidative dehydrogenation process as well as the overoxidation of propane and propylene might occur. The multiple parallel reaction pathways and diverse intermediates make it challenging to use C-based catalysts in O2-ODHP.143
A series of carbon nitride (CN) polymers provided more opportunities for further developing C-based catalysts with high selectivity in the O2-ODHP field. For instance, graphitic CN has been reported as an O2-ODHP catalyst achieving a propane conversion of 12.8% and a selectivity of 74.7% for the desired propylene product.140In situ DRIFTS characterization showed the orderly formation of –CO and –OH species (Fig. 7g and h). Accordingly, the authors speculated the reaction pathway. Specifically, the C
O abstracted H atoms from propane to generate propylene and formed C–OH. Subsequently, the following O2 treatment extracted H atoms from C–OH species and re-generated C
O groups. The switch between C
O and C–OH groups under the reactant atmosphere drove the O2-ODHP process while avoiding the direct interaction between C3H8 and O2. Theoretically, DFT calculations predicted a reasonable O2-ODHP reaction pathway, during which two adjacent C
O groups on the edge of CN can activate two C–H bonds of propane simultaneously, thus avoiding the generation of free C3H7* radicals and leading to the high selectivity of light olefins.
Due to the high light absorption ability of C-based materials, photothermal catalytic technology has emerged as a novel strategy to trigger O2-ODHP reaction under mild conditions. Compared to traditional reaction conditions, the lower reaction temperature during photothermal catalysis benefits the control of propylene selectivity and catalyst stability. One recent example is based on a semiconducting boron CN (BCN) catalyst, which integrated photocatalysis with thermal catalysis.141 The fundamental theory is the cooperation between electrically conductive properties from the graphene phase and insulation from boron nitride (h-BN). Thus, the BCN alloy could possess both excellent catalytic performance of h-BN towards O2-ODHP and the semiconductor behavior of graphene. Ultraviolet-visible diffuse reflectance spectroscopy (UV-DRS) showed the effective visible light response of BCN with a 500 nm absorption edge (Fig. 7i). And the Tauc plot clarified that BCN was a direct band gap material, which had a band gap of 2.83 eV (Fig. 7j). According to ultraviolet photoelectron spectroscopy (UPS) measurements, the authors determined the energy band structure of BCN with 3.92 eV vs. evac conduction band and 6.62 eV vs. evac valence band, indicating the photocatalytic ability for activating O2 (Fig. 7k). At 480 °C, the propane conversion increased from 11.1% to 21.2% via light irradiation, resulting in an increase in the overall yield of olefin from 10.4% to 19.6% and higher propylene selectivity. DFT results clearly showed the advantage of photothermal catalysis in terms of the energy barrier and the rate-determining step (the dehydrogenation of B-OH groups to form BO·) in the BCN catalytic system (Fig. 7l). Indeed, this work confirmed the great potential of integrating photocatalysis with conventional thermal catalysis via semiconductor catalysts to enhance propane conversion and propylene selectivity during the O2-ODHP reaction. The photothermal method will play a key role in achieving carbon neutrality.
Catalyst | Reaction condition | Performance | Ref. |
---|---|---|---|
V2O5/SiO2 | Reactor type: conventional continuous flow quartz microreactor | Selectivity towards propylene is 72.0% | 16 |
Reactant ratio: C3H8![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
|||
Reactant feed: GHSV = 1700 h−1 | |||
Reaction temperature: 450 °C | |||
Reaction pressure: 1 atm | |||
ZrMo2O8/ZrO2 | Reactor type: packed-bed tubular quartz reactor | Selectivity towards propylene is 80.0% | 23 |
Reactants: propane and oxygen at 14.03 and 1.74 kPa | |||
Reaction temperature: 500 °C | |||
Ti–Ni–O | Reactant ratio: C3H8![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
Propane conversion is about 28.4%, selectivity to propylene is 42.5% and propylene yield is 12.1% | 36 |
Reactant feed: 9000 mL h−1 g−1 | |||
Reaction temperature: 300 °C | |||
MoVOx | Reactor type: fixed bed reactor | Propane conversion is about 2.7%, selectivity to propylene is 100% and propylene yield is 2.7% | 45 |
Reactant ratio: C3H8![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
|||
Reaction temperature: 350 °C | |||
Pt–Ce/γ-Al2O3 | Reactor type: freestanding anaerobic-anoxic-oxic flow reactor setup | Propane conversion is about 37.0% and selectivity to propylene is 96.0% | 49 |
Reactant ratio: C3H8![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
|||
Reaction temperature: 576 °C | |||
Ni-CA SACs | Reactor type: fixed bed reactor | Selectivity to propylene is 47.3% and propylene yield is 14.2% | 53 |
Reactant ratio: C3H8![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
|||
Reactant feed: WHSV = 12![]() |
|||
Reaction temperature: 580 °C | |||
Reaction pressure: 1 atm | |||
NiO−CeO2 | Reactor type: fixed-bed flow reactor | Propane conversion is about 69.0%, selectivity to propylene is 80.0% and propylene yield is 55.0% | 57 |
Reactant ratio: C3H8![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
|||
Reactant feed: 48 mL min−1 | |||
Reaction temperature: 500 °C | |||
Reaction pressure: 1 atm | |||
P-CeO2 | Reactor type: fixed-bed quartz tube down-flow reactor | Propane conversion is about 8.7% and selectivity to propylene is 70.1% | 59 |
Reactant ratio: C3H8![]() ![]() ![]() ![]() |
|||
Reactant feed: VHSV = 9000 h−1 | |||
Reaction temperature: 600 °C | |||
h-BN | Reactor type: quartz tube reactor | Propane conversion is about 14.0% and selectivity to propylene is 79.0% | 65 |
Reactants: 0.15 atm O2 and 0.3 atm C3H8 | |||
Reaction temperature: 490 °C. | |||
BNOH | Reactor type: packed-bed quartz microreactor | Propane conversion is about 20.6% and selectivity to propylene is 80.2%. | 66 |
Reactant ratio: C3H8![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
|||
Reactant feed: 192 mL min−1 | |||
Reaction temperature: 530 °C | |||
BN with a larger specific surface area | Reactor type: fixed-bed reactor | Propane conversion is about 50.0% and selectivity to propylene is 70.0% | 67 |
Reactant ratio: C3H8![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
|||
Reactant feed: 20 mL min−1 | |||
Reaction temperature: 525 °C | |||
Reaction pressure: 1 atm | |||
N2 plasma treated BN | Reactor type: packed-bed quartz reactor | Propane conversion is 26.0% and the total selectivity towards propylene and ethene is 89.4% | 69 |
Reactant ratio: C3H8![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
|||
Reactant feed: 48 mL min−1 | |||
Reaction temperature: 520 °C | |||
Reaction pressure: atmospheric pressure | |||
B4C | Reactor type: quartz reactor tube | Propane conversion is about 5.1% and selectivity to propylene is 86.0% | 73 |
Reactant ratio: C3H8![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
|||
Reactant feed: WHSV−1 = 10.5 kg-cat * s mol−1 C3H8 | |||
Reaction temperature: 500 °C | |||
Reaction pressure: 0.30 atm C3H8 and 0.15 atm O2 | |||
Ti2B | Reactor type: quartz reactor tube | Propane conversion is about 5.2% and selectivity to propylene is 85.3% | 73 |
Reactant ratio: C3H8![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
|||
Reactant feed: WHSV−1 = 12.5 kg-cat * s mol−1 C3H8 | |||
Reaction temperature: 500 °C | |||
Reaction pressure: 0.30 atm C3H8 and 0.15 atm O2 | |||
Exfoliated layered MgB2 | Reactor type: fixed-bed reactor | Propane conversion is about 39.8%, selectivity to propylene is 63.5% and propylene yield is 26.3% | 78 |
Reactant ratio: C3H8![]() ![]() ![]() ![]() ![]() ![]() |
|||
Reactant feed: WHSV = 24![]() |
|||
Reaction temperature: 530 °C | |||
Reaction pressure: atmospheric pressure | |||
N doped CNT | Reactor type: immobilized bed quartz reactor | Propane conversion is about 5.0% and selectivity to propylene is 56.5% | 82 |
Reactant ratio: C3H8![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
|||
Reactant feed: GHSV = 13 L g−1 h−1 | |||
Reaction temperature: 400 °C | |||
Reaction pressure: atmospheric pressure |
(1) Despite the variety of catalysts developed for O2-ODHP, to the best of our knowledge, the highest yield towards propylene was achieved over the NiO–CeO2 catalyst, up to 55.0%, with a propane conversion of 69.0% and a selectivity towards propylene of 80.0%. Apparently, their performances are still far away from industrialization. In addition, for most of the catalytic systems, it is still challenging to realize a high propane conversion rate and high selectivity towards propylene simultaneously. Therefore, in the future, it is demanded to develop a novel catalyst, which can achieve high selectivity towards propylene at high propane conversion rates.144,145
(2) Even though characterization and theoretical calculations have been conducted over some of the catalytic systems, the intrinsic mechanism remains unclear. As is known, the catalytic reaction mechanism plays an important role in guiding the development of more advanced and efficient catalysts. Therefore, in the future, more studies should be focused on the investigation of the reaction system, using in situ/operando characterization techniques as well as theoretical simulations, to unravel the intrinsic active site and reaction pathways of the O2-ODHP process.
(3) Considering the two facts that (i) propane as the reactant has to be extracted from other raw materials like shale gas and (ii) propylene is important for downstream chemical reactions to get high value-added products, the coupling between O2-ODHP and other catalytic reactions should be paid more attention to realize the large-scale industrialization and even achieve the green chemistry and carbon neutral aims. For this aim, a cheaper catalyst should be designed and developed to meet the requirements of continuous and large-scale catalytic reactions. What's more, the reaction vessel for conducting different catalytic processes should also be constructed to promote the advanced catalytic systems involving the combination of the shale gas–propane conversion and O2-ODHP, and the coupling between O2-ODHP and reactions converting propylene into value-added products.
This journal is © The Royal Society of Chemistry 2024 |