Abstract
The oxidative instability of biodiesel remains a critical barrier to its widespread adoption despite its advantages as renewable, biodegradable, and low-emission fuel. Antioxidant additives are an established strategy to suppress free radical chain reactions, yet their efficiency is strongly modulated by molecular structure and solvent environment. This is the first comparative density functional theory study of dibrominated dimethoxybenzaldehydes and standard phenolic antioxidants under biodiesel-relevant solvent conditions using the conductor-like polarizable continuum model. Frontier molecular orbitals, Fukui index, ionization potentials, spin density distributions, and natural bond orbital hyperconjugations were systematically analyzed across polar and nonpolar environments. The computational results suggest that bromination is associated with increased electronic softness and electron transfer potential, while also leading to changes in the stability of radical intermediates, especially in ortho-substituted derivatives. Among the dibrominated compounds, IB1 exhibits the most balanced combination of computed properties, whereas IB3, although highly reactive in silico, is predicted to form comparatively less stable radical species. Compared with commercial benchmarks, these halogenated systems constitute a distinct mechanistic class governed by polarization rather than hydroxyl-centered resonance. These computational findings provide guidance for the rational design of next-generation biodiesel stabilizers, pending future experimental validation.
1 Introduction
The transition to sustainable energy sources is a central challenge for global energy security and environmental stewardship (Yusuf et al., 2011; ; Usmani et al., 2023). Increasing the efficiency and operational stability of renewable fuels is vital for accelerating the energy transition and achieving lower carbon emissions. Biodiesel, produced from renewable feedstocks such as vegetable oils, and animal fats, has emerged as a promising alternative to petroleum-derived fuels due to its biodegradability, low toxicity, and reduced greenhouse gas emissions compared with fossil fuels (; ). Integration of biodiesel into existing and future energy infrastructures directly supports clean energy goals and sustainable process design. Despite these advantages, the large-scale deployment of biodiesel remains hindered by its intrinsic oxidative instability. Overcoming these process limitations is essential to enable the broader adoption of sustainable fuels within the context of advanced process engineering and reactor optimization for the energy transition. Oxidation during storage and use accelerates the formation of peroxides, acids, and polymers, leading to increased viscosity, corrosiveness, and sediment deposition, ultimately impairing both fuel quality and engine performance (; Pullen and Saeed, 2012; Saluja et al., 2016; Yang et al., 2013; Pullen and Saeed, 2014; Pantoja et al., 2013; Serqueira et al., 2021). This challenge is further exacerbated in tropical and subtropical climates such as those in Brazil () (Köppen area A4), where elevated temperatures, humidity, and solar irradiation intensify oxidative degradation (Varatharajan and Pushparani, 2018). Agricultural operations are especially impacted: machinery exposed to these conditions and fueled by blends experiences more rapid loss of fuel quality and higher maintenance requirements (; ; Ryu, 2009; Silva et al., 2025).
A well-established strategy to counteract these effects is the incorporation of antioxidant additives, which suppress autoxidation, extend shelf life, and enhance operational reliability (Saxena et al., 2017; Rajamohan et al., 2022; ). Phenolic compounds have been widely employed for this purpose because of their strong capacity to interrupt free radical chain reactions via electron or hydrogen donation (Schober and Mittelbach, 2004; Sundus et al., 2017; ). However, the efficacy of phenolic antioxidants is strongly modulated by their molecular structure and the surrounding medium (). Understanding these structure–activity relationships is therefore essential for rational design of efficient stabilizers. A robust understanding of these structure–activity relationships and their process implications is crucial for advancing energy technologies that align with the goals of cleaner production. Density Functional Theory (DFT) (; ) has become a powerful tool for elucidating the electronic and structural determinants of antioxidant performance. By probing reactivity descriptors, frontier orbital distributions, and ionization potentials, computational methods provide predictive insight into antioxidant activity under diverse conditions. Prior studies of commercial phenolic antioxidants including butylated hydroxyanisole (BHA), butylated hydroxytoluene (BHT), gallic acid (GA), propyl gallate (PG), pyrogallol (PY), and tert butylhydroquinone (TBHQ) have established benchmark profiles of electron donating ability and radical stabilization.
Within this broader landscape of potential stabilizers, phenolic antioxidants remain central because their radical scavenging efficiency has been consistently rationalized in terms of electronic descriptors. The commercial use of this class of compounds has established a robust theoretical benchmark, since their activity can be directly correlated with well-defined structure–activity relationships involving hydrogen donation and resonance stabilization of phenoxyl radicals (Yang et al., 2013; Yaakob et al., 2014; Mittelbach and Schober, 2003). This consolidated body of knowledge provides both a conceptual and computational framework against which new molecules can be evaluated. Dimethoxybenzaldehyde (DMB) derivatives extend this framework by combining methoxy groups, which enhance donation, with carbonyl and halogen substituents that introduce inductive influences capable of modulating frontier orbital energies and reactivity indices (; Wenceslau et al., 2025; ; ; ; ; ; ; ). Bromine substitution, in particular, alters electron density distribution in ways that may enhance radical stabilization but also requires careful consideration regarding fuel compatibility. Positioning DMB derivatives in direct comparison with well-characterized phenolic antioxidants therefore establishes a coherent theoretical basis for probing how substituent effects govern antioxidant performance in biodiesel (; ; Sterpu et al., 2024).
Here, we extend the computational framework to three dibrominated dimethoxybenzaldehyde derivatives 4,5 dibromo 2,3 dimethoxybenzaldehyde (IB1), 2,3 dibromo 5,6 dimethoxybenzaldehyde (IB2), and 4,6 dibromo 2,3 dimethoxybenzaldehyde (IB3) whose antioxidant potential has not been systematically investigated. The antioxidant properties and fuel stabilization potential of dibrominated DMBs derivatives remain largely unknown, as their mechanisms and synergistic effects have not been systematically evaluated in biodiesel contexts. This study addresses this gap by computationally investigating their electronic structure and radical-scavenging efficiency relative to established phenolic antioxidants. Their performance was evaluated through comparison with commercial phenolics, focusing on ionization potential, frontier orbital distributions, Fukui reactivity indices, spin density mapping, and Natural Bond Orbital (NBO) hyperconjugation analysis. By integrating molecular modeling with comparative benchmarking, this study provides a mechanistic basis for assessing brominated derivatives as potential biodiesel stabilizers. The comparative evaluation of DMBs derivatives and commercial phenolics reveals key structure–activity relationships that determine radical scavenging efficiency.
2 Computational procedures
2.1 Electron density calculations
All theoretical calculations were performed using the Gaussian 16 (; ). Molecular geometries were optimized at the hybrid M06-2X (Zhao and Truhlar, 2008), which incorporates long-range exchange corrections and is widely used for thermochemistry, kinetics, and noncovalent interactions. Studies have demonstrated that this functional reliably describes mid-range electronic correlation effects, and non-covalent interactions and is one of the best-performing functionals for modeling the thermodynamics of chemical processes (; ). Furthermore, this functional has shown good results in analyses of the antioxidant activity of compounds (; Thuy and Son, 2022; ). Self-consistent field (SCF) calculations employed tight convergence criteria and an ultrafine integration grid. Geometry optimizations maintained stringent thresholds, and vibrational frequency analyses confirmed that all structures represent true minimum, with no imaginary frequencies. Unless otherwise specified, thermal corrections were obtained at 298.15 K and 1 atm.
Electronic structures were analyzed in terms of Frontier Molecular Orbitals (FMOs) (Zhang and Musgrave, 2007): the Highest Occupied Molecular Orbital (HOMO), associated with electron-donating ability, and the Lowest Unoccupied Molecular Orbital (LUMO), associated with electron-accepting ability (; Rocha et al., 2015; ; NC et al., 2023). Within conceptual DFT, the central descriptors are defined from the total electronic energy as functions of the number of electrons under a fixed external potential . The chemical potential (Equation 1) (Pearson, 1992),the chemical hardness (Equation 2) (Pearson, 1992; Pearson, 2005)and the global electrophilicity index (ω) (Equation 3) (Parr et al., 1999)were calculated to obtain information about the chemical reactivity of the compounds studied in this work. In this context, is the total electronic energy, is the number of electrons, and is the external potential (arising from the nuclei and any fixed fields). The chemical potential () is related to electronegativity () by , the ionization energy () is approximated as , and the electron affinity (A) as (). The chemical hardness () measures the resistance of the electronic cloud to deformation during charge transfer, while the global electrophilicity index () quantifies the stabilization achieved upon electron acquisition.
2.2 Antioxidant analysis
The Fukui index (; ) was applied to identify molecular sites most susceptible to radical attack. Atomic charges were obtained through natural population analysis, and finite-difference approximations were used to derive the reactivity indices. These indices guided the generation of radical species, which were subsequently optimized at the same theoretical level (Weinhold et al., 2016; Reed et al., 1985; ; Yamauchi et al., 2024). There are several methods for determining the antioxidant activity of compounds (Wenceslau et al., 2025), and to determine these property of the dibrominated DMBs derivatives, thermodynamic descriptors were calculated in the gas phase, in the same level of theory. Two main mechanisms can explain the antioxidant activity of compounds. The first, the hydrogen atom transfer (HAT) and the other, proton transfer (; ). Hydrogen transfer processes can occur through homolytic or heterolytic scission of hydrogen atoms. The heterolytic pathway involves either a single electron transfer followed by proton transfer (SET-PT) or sequential proton loss followed by electron transfer (SPLET). The schematic representation of both mechanisms is shown in Scheme 1.
SCHEME 1
Since IB1, IB2, and IB3 have no phenolic groups, the antioxidant activity calculations were conducted only using the ET mechanism, in which one electron was removed from each structure and the ionization potential of the dibrominated DMBs derivatives was calculated using Equation 4,in which , , and are the dibrominated DMBs derivatives, radical cation formed after free radical scavenging, and electron enthalpies. The lower the IP value, the greater the antioxidant activity of the azine. The results were tabulated and compared. Furthermore, spin density distributions were obtained to predict the stability of the radicals formed.
To investigate the influence of the medium on the electron-transfer step, the ionization potentials (IP) were calculated in different solvation environments. For this purpose, the implicit conductor-like polarizable continuum (CPCM) model (Takano and Houk, 2005) was employed, using solvents with distinct dielectric constants to represent hydrophobic and hydrophilic microenvironments relevant to biodiesel-like systems. Benzene (ε = 2.27) and toluene (ε = 2.37) were selected as representatives of nonpolar media because their polarity is similar to the hydrophobic matrix of biodiesel, which is predominantly composed of fatty acid esters with dielectric constants between 3 and 4.5. These solvents therefore allow the low-polarity microenvironment in which antioxidant additives operate to be reasonably mimicked. The inclusion of ethanol (ε = 24.85) was intended to model residual polar components that may remain from the transesterification process, affecting local polarity and the solvation profile of the compounds (Patiño-Camino et al., 2021; Plácido and Capareda, 2016). Water (ε = 78.35) was also considered due to the hygroscopic nature of biodiesel (; ), which enables the incorporation of small amounts of this contaminant and can substantially alter the dielectric properties of the medium. In this manner, it becomes possible to assess how different physicochemical environments influence the ionization process without anticipating any interpretation regarding antioxidant performance.
Radical stabilization was further investigated by spin density distribution and NBO analysis (Weinhold and Landis, 2012; Weinhold and Landis, 2001) through the donor–acceptor hyperconjugation () obtained by the second-order perturbation formula (Equation 5),where or is the Fock matrix element between the natural bond orbitals i and j; is the energy of the antibonding orbital , and is the energy of the bonding orbital ; represents the population occupation of the donor orbital.
To complement these quantum descriptors, a machine learning approach was adopted for estimating radical reaction rate constants (). Calculations leveraged the pySiRC (Sanches-Neto et al., 2021) framework, employing Morgan fingerprints and the XGBoost algorithm to predict for each target molecule under oxidative attack by hydroxyl radicals in model biodiesel environments. Performance benchmarks and method validation followed protocols previously reported (). This allowed direct computational comparison of dibrominated DMBs (IB1, IB2, IB3) and reference commercial antioxidants (BHT, TBHQ, BHA, PG, PY, GA) (Figure 1), supporting a predictive ranking of radical-trapping potentials in biodiesel-relevant conditions.
FIGURE 1
3 Results and discussion
3.1 Molecular modeling analysis
The spatial distributions of the HOMO and LUMO (Figure 2) provide computational insights into the reactivity profiles of the studied compounds (; ). In IB1, the computed HOMO is extensively delocalized over the aromatic backbone and adjacent methoxy and bromine substituents, indicating a pronounced propensity for electron donation from these conjugated regions. The LUMO is primarily localized on the brominated and methoxylated sites, suggestive of a substituent-directed capacity for charge acceptance that may modulate interaction with radical species. IB3, conversely, exhibits both HOMO and LUMO densities that are distinctly more localized, confined to specific aromatic and substituent positions; this spatial restriction is associated with decreased π-delocalization and enhanced electronic softness, factors that favor increased reactivity but potentially at the cost of reduced radical stabilization. IB2 displays intermediate behavior, with modest orbital delocalization and energetic descriptors that fall between IB1 and IB3.
FIGURE 2
Among phenolic benchmarks, BHA and BHT show HOMO density highly concentrated on the hydroxyl moieties and ortho/para positions of the aromatic core. This canonical distribution aligns with classical mechanisms for HAT, where delocalized spin upon oxidation stabilizes the resulting phenoxyl radicals. In contrast, GA, PG, and PY each exhibit broader orbital delocalization in their HOMOs, spanning multiple hydroxyl sites and aromatic centers, which enhances radical stabilization but may slow initial scavenging kinetics (Varatharajan and Pushparani, 2018; ). The complementary LUMO maps reveal diffuse and, in some cases, spatially extended character particularly pronounced in PG, PY, and TBHQ. For these molecules, the computed LUMOs exhibit substantial electron density expansion beyond the aromatic moiety into peripherally located regions, a hallmark indication of Rydberg orbital character. The presence of these Rydberg orbitals endows TBHQ, PY, and BHA with alternative electron-accepting channels and non-classical stabilization pathways for radicals, as observed in both computational theory and experimental radical chemistry.
Quantitative descriptors (Supplementary Table S1) consistently support these qualitative trends. IB1 displays the highest HOMO energy (≈- 8.39 eV) among the DMB derivatives, reflecting its strong electron-donating profile and increased electronic hardness. By contrast, IB3, characterized by the narrowest HOMO–LUMO gap (≈6.53 eV) and lowest ionization energy (≈7.89 eV), is definitively the most electronically “soft” and computationally reactive member within the new inhibitor series. IB2 presents an intermediate electronic configuration. Solvent-dependent analyses further suggest that while IB1 and IB2 maintain their electronic features across environments of varying polarity, IB3’s descriptors, especially its electronic gap and ionization potential, are substantially modulated by medium, reflecting a heightened sensitivity to environmental changes as predicted within the computational protocol. This effect may have implications for radical persistence under operational biodiesel conditions, but experimental corroboration is warranted.
Commercial antioxidants define well-established mechanistic archetypes: GA and PG have the largest calculated HOMO–LUMO gaps (≈7.26–7.37 eV) and highest calculated hardness, which typically correlate with enhanced kinetic stability but lower inherent reactivity. BHA and BHT, with moderate gaps (≈6.80–6.96 eV) and the lowest ionization energy (≈6.99–7.12 eV), are computationally positioned as rapid radical quenchers. TBHQ, sharing a moderate gap and showing pronounced solvent modulation, reflects practical efficiency in nonpolar biodiesel systems.
3.2 Antioxidant potential
Spin density mapping reveals the spatial distribution of unpaired electrons following radical formation. For the dibrominated derivatives (IB1–IB3), the unpaired electron preferentially localizes on carbons adjacent to brominated positions, with partial delocalization into the aromatic π-system (Figure 3). In IB1, spin density is found mainly on C2, C3, and C5; IB2 concentrates spin density at C3 and C6, reflecting a more localized distribution but still partially extended across the aromatic ring. IB3 exhibits localization at C3 and C5, but with reduced overlap into the π-framework, consistent with the structural distortion imposed by ortho bromination.
FIGURE 3
Commercial antioxidants display distinct but complementary stabilization strategies. BHA and BHT show spin densities centered on the phenolic oxygen and adjacent carbons, directly benefiting from strong O–H to π conjugation, which efficiently stabilizes the radical cation. GA and PG exhibit broader spin distribution patterns, associated with their multiple hydroxyl substitutions, while PY and TBHQ exhibit intermediate patterns, with spin density distributed between oxygen substituents and aromatic carbons, reflecting their capacity to quench radicals efficiently in biodiesel-like environments (; ; Rizwanul Fattah et al., 2014; Sui et al., 2021; ; Rodrigues et al., 2020; ). The spin density analysis consolidates the mechanistic framework of the dibrominated series: antioxidant activity is mediated by electron transfer, but the efficiency of radical stabilization is strongly modulated by substitution geometry. In contrast, commercial antioxidants achieve stabilization primarily through well-established O–H centered delocalization pathways, underscoring the advantage of phenolic hydroxyl groups over halogen substitution in ensuring persistent radical quenching.
The NBO analysis (Supplementary Table S2–S4) revealed that the O2 and O3 atoms experience the least stabilization due to the hyperconjugations occurring in the radicals. Notably, IB3 exhibits hyperconjugation from π(C1–C6) π*(O2–C2) with a stabilizing energy of ≈1.07 eV, effectively stabilizing the radical on the O2 atom. When the unpaired electron is located at C2, there are two important hyperconjugations that significantly contribute to stabilizing the IB1 radical: π(C5–C6) π*(C1–C2) with an value of ≈0.569 eV, and (O2) π*(C1–C2) with an value of ≈0.734 eV. In contrast, for the IB2 radical, only one hyperconjugation strongly aids in its stabilization: (O2) π*(C1–C2), with an value of ≈1.15 eV. Finally, in the IB3 radical, the hyperconjugation that contributes most to its stabilization is (O2) σ*(C2–C3). However, the value for this interaction is lower, indicating relatively poor stabilization of the radical in this region of the molecule. When the unpaired electron is located on the C3 atom, there are three significant hyperconjugations involving donor orbitals π(C1–C2), π(C5–C6), and (Br1) of IB1 with the acceptor orbital π*(C3–C4) contribute to its stabilization, with values of ≈0.525, ≈0.703, and ≈0.390 eV, respectively. In the case of IB2 and IB3, only one crucial hyperconjugation each was observed in the stabilization of their radicals: (O2) π*(C1–C2) with an value of ≈1.15 eV for IB2, and (O2) σ*(C2–C3) with an value of ≈0.241 eV for IB3.
Collectively, the NBO results underscore the decisive role of substitution geometry in dictating radical stability. IB1 achieves stabilization through a distributed network of hyperconjugations, IB2 depends on one dominant interaction, and IB3 is structurally penalized by steric congestion, which suppresses orbital overlap. These findings reconcile the paradox observed in global descriptors: although IB3 exhibits the lowest ionization potential and narrowest HOMO–LUMO gap, its radical intermediates lack sufficient hyperconjugative support, reducing its efficacy as an antioxidant. The duality of bromination—enhancing electronic softness while constraining delocalization—emerges here as a fundamental trade-off that limits the utility of this class in biodiesel stabilization.
The antioxidant potential of the dibrominated derivatives was evaluated employing site-specific reactivity indices, radical stability parameters, and ionization potential metrics. Fukui’s function (Figure 4a) consistently identified bromine substituents and the phenolic oxygen as the most electrophilic sites for potential radical attack, with the electronic polarization induced by bromination rendering these sites particularly susceptible to electron abstraction. Notably, the absence of labile hydrogen atoms in DMB molecules indicates HAT is not their operative scavenging pathway; instead, their antioxidant mechanism is primarily governed by electron transfer processes. This mechanism is supported by the computed electronic descriptors and is further illustrated in Scheme 1, where the free radical (R·) captures an electron from the additive, leading to the formation of a radical cation and a neutralized radical species.
FIGURE 4
In this process, R· captures an electron from the additive (DMB), generating a radical cation ([DMB·]+) and a neutralized radical (R–) as schematically illustrated in Scheme:
The relative stability of the radical intermediates, expressed as enthalpy differences (ΔH, Figure 4b), follows the order IB3 > IB2 > IB1. These values are comparable to those of commercial antioxidants such as GA and PG, suggesting that despite structural distortions introduced by bromination, the DMBs can form radicals of competitive stability. In these radicals, the unpaired electron preferentially localizes on C2, C3, and C5 in IB1, C3 and C6 in IB2, and C3 and C5 in IB3 before undergoing free radical scavenging. The IP calculations further refine this mechanistic framework. Gas-phase values are ≈8.48, ≈8.43, and ≈8.17 eV for IB1, IB2, and IB3, respectively, which situate them within the reference range of phenol (Xue et al., 2013a) (≈8.31 eV), glutathione () (≈7.55 eV), phenolic chalcones (Xue et al., 2013b) (mean: ≈7.55 eV), and other chalcones () (≈8.02 eV). Solvent effects significantly modulate these values, with reductions of ≈12% in nonpolar solvents and up to ≈20% in polar solvents. This trend enhances reactivity in the presence of water an inevitable biodiesel contaminant due to its hygroscopic nature thereby ensuring antioxidative activity under practical conditions (; ; Rodriguez et al., 2018).
Commercial antioxidants display similar solvent-dependent behavior (Schober and Mittelbach, 2004; Serqueira et al., 2015; ; ; ). BHA and BHT exhibit pronounced reductions, from ≈7.46–7.55 eV in the gas phase to ≈5.68–5.94 eV in water (Figure 4c), rationalizing their widespread efficiency in biodiesel applications. GA, PG, and PY, by contrast, maintain consistently high IPs with limited solvent sensitivity, conferring robustness but reduced reactivity (; ; ). TBHQ represents an intermediate case, with reactivity enhanced in nonpolar solvents, consistent with its efficiency in biodiesel-like matrices. While IB1 emerges as the most thermodynamically stable, IB3 is the most reactive but least stabilized, and IB2 exhibits intermediate behavior. These mechanistic insights confirm that dibromination can produce competitive descriptors relative to established phenolic antioxidants, although the dual influence of bromine facilitating ET but constraining delocalization—remains a fundamental computationally observed limitation of these structures.
The machine learning-derived reaction rate constants provide predictive (Table 1) insights into the radical-trapping potential of DMB derivatives relative to established commercial antioxidants. The three DMB derivatives exhibit values that situate them within the same theoretical reactivity regime as BHT and GA, but below the highly reactive benchmarks TBHQ, BHA, and particularly PG, which was identified as the most reactive antioxidant among those assessed. This computational ranking aligns qualitatively with literature trends for antioxidant efficiency in biodiesel, where PG, TBHQ, and BHA typically suppress oxidation and extend induction periods more effectively than BHT and GA under Rancimat conditions (Varatharajan and Pushparani, 2018; ; Schober and Mittelbach, 2004; Mittelbach and Schober, 2003; Yamauchi et al., 2024; Singh et al., 2022; ; ). Consequently, the moderate ( values predicted for the DMBs suggest radical scavenging capabilities comparable to BHT, indicating suitability for scenarios requiring controlled stabilization.
TABLE 1
| Compound | Reaction rate constant |
|---|---|
| IB1 | |
| IB2 | |
| IB3 | |
| GA | |
| TBHQ | |
| BHA | |
| PG |
Reaction Rate ( for DMB derivatives and commercial antioxidants.
Furthermore, the similarity between IB1 and IB3 suggests that their structural differences have minimal impact on hydroxyl radical trapping, whereas the lower reactivity of IB2 implies specific electronic or steric constraints. While the present work is strictly theoretical, the derived descriptor–kinetics hierarchy establishes a concrete framework for future validation, serving as a bridge to practical observations. To confirm the predicted performance, the rank ordering established herein—based on IP, softness, NBO hyperconjugation, and ML-derived ( should be directly confronted with a targeted experimental protocol. This validation would necessitate DPPH and ABTS radical-scavenging assays to determine ( values, correlating kinetic rates with computed electronic descriptors, alongside DSC and TGA measurements on treated biodiesel to determine oxidation onset temperatures and mass-loss profiles. Thus, this study provides the theoretical foundation and a concise roadmap for subsequent experimental verification of DMB derivatives as effective biodiesel stabilizers.
4 Final considerations
This study integrated computational descriptors frontier molecular orbitals, Fukui indices, ionization energies, spin density distributions, and NBO hyperconjugation analyses to elucidate the antioxidant mechanisms of dibrominated DMBs derivatives compared to commercial phenolic benchmarks. Our calculations suggest that bromination enhances electronic softness and facilitates electron transfer, positioning DMB additives as computationally promising for radical scavenging in biodiesel stabilization. However, increased softness and planarity deviation, particularly in ortho-substituted derivatives, may restrict conjugation and reduce the persistence of radical intermediates.
Within the DMB series, IB1 is predicted to combine a favorable electron-donating profile with distributed stabilization potential, while IB3 displayed high reactivity but limited radical stability due to local electronic effects. IB2 exhibits intermediate behavior based on our analyses. When benchmarked against commercial antioxidants, distinct mechanistic contrasts are indicated by our computational results. BHA and BHT are characterized in our study as achieve efficiency through low ionization potentials and hydroxyl-centered resonance stabilization; GA and PG rely on wide HOMO–LUMO gaps and extended delocalization across multiple hydroxyl groups; TBHQ benefits from solvent-dependent reactivity optimized for nonpolar biodiesel matrices. The dibrominated derivatives, by contrast, operate through halogen-induced polarization and ET pathways, situating them in a mechanistically distinct class. These results demonstrate that the antioxidant potential of DMB derivatives depends on achieving the right balance between reactivity and radical stability. Optimization of substitution patterns can be guided by computational descriptors established in this study, and may inform the design of future fuel additives, once validated experimentally.
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.
Author contributions
IB: Conceptualization, Data curation, Formal Analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing – review and editing, Writing – original draft. AA: Conceptualization, Data curation, Formal Analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing – original draft. AC: Conceptualization, Data curation, Formal Analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing – review and editing. HN: Conceptualization, Data curation, Formal Analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing – review and editing.
Funding
The author(s) declared that financial support was not received for this work and/or its publication.
Acknowledgments
The authors are grateful to Coordenação de Aperfeiçoamento de Pessoal de Nível Superior, Conselho Nacional de Desenvolvimento Científico e Tecnológico and Fundação de Amparo à Pesquisa de Goiás. Theoretical calculations were performed in the High-Performance Computing Center of the Universidade Estadual de Goiás.
Conflict of interest
The author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declared that generative AI was not used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fceng.2025.1716732/full#supplementary-material
References
1
AguiarA. S. N.Dos SantosV. D.BorgesI. D.NavarreteA.AguirreG.ValverdeC.et al (2022a). Bromine substitution effect on structure, reactivity, and linear and third-order nonlinear optical properties of 2,3-Dimethoxybenzaldehyde. J. Phys. Chem. A126, 7852–7863. 10.1021/acs.jpca.2c04658
2
AguiarA. S. N.BorgesI. D.BorgesL. L.DiasL. D.CamargoA. J.PerjesiP.et al (2022b). New insights on glutathione’s supramolecular arrangement and its in silico analysis as an angiotensin-converting enzyme inhibitor. Molecules27, 7958. 10.3390/molecules27227958
3
AguiarA. S. N.DiasP. G. M.QueirozJ. E.FirminoP. P.CustódioJ. M. F.DiasL. D.et al (2023). Insights on potential photoprotective activity of two butylchalcone derivatives: synthesis, spectroscopic characterization and molecular modeling. Photonics10, 228. 10.3390/photonics10030228
4
AlabuginI. V.GilmoreK. M.PetersonP. W. H. (2011). Wiley Interdisciplinary Reviews: Computational Molecular Science1, 109–141. 10.1002/wcms.6
5
AlmeidaE. S.PortelaF. M.SousaR. M.DanielD.TerronesM. G.RichterE. M.et al (2011). Behaviour of the antioxidant tert-butylhydroquinone on the storage stability and corrosive character of biodiesel. Fuel90, 3480–3484. 10.1016/j.fuel.2011.06.056
6
AlvaresC. A.StapeJ. L.SentelhasP. C.GonçalvesJ. L. de M.SparovekG. (2013). Köppen’s climate classification map for Brazil. Meteorol. Zeitschrift22, 711–728. 10.1127/0941-2948/2013/0507
7
AmranN. A.BelloU.Hazwan RuslanM. S. (2022). The role of antioxidants in improving biodiesel’s oxidative stability, poor cold flow properties, and the effects of the duo on engine performance: a review. Heliyon8, e09846. 10.1016/j.heliyon.2022.e09846
8
ArumugamS.MuthaiyanR.RajendranS. (2023). Assessment of oxidative stability of biodiesel and biodiesel blends. Energy sources.45, 6371–6387. 10.1080/15567036.2023.2217089
9
BarreraN. F.Cabezas-EscaresJ.MuñozF.MurielW. A.GómezT.CalatayudM.et al (2025). Fukui function and Fukui potential for solid-state chemistry: application to surface reactivity. J. Chem. Theory Comput.21, 3187–3203. 10.1021/acs.jctc.5c00086
10
BorgesI.NavarreteA.AguirreG.AguiarA.OliveiraS.CamargoA.et al (2022). Synthesis and molecular modeling study of two bromo-dimethoxybenzaldehydes. J. Braz Chem. Soc.10.21577/0103-5053.20220023
11
BorgesI. D.FariaE. C. M.CustódioJ. F. M.DuarteV. S.FernandesF. S.AlonsoC. G.et al (2022). Insights into chalcone analogues with potential as antioxidant additives in diesel–biodiesel blends. RSC Adv.12, 34746–34759. 10.1039/d2ra07300e
12
ChenP.WangZ.WuJ.XiaH.TaiC.LiR. (2019). Effects of phenolic antioxidants on biodiesel oxidative stability and emission. Environ. Prog. Sustain Energy38, 13203. 10.1002/ep.13203
13
ChristensenE. D.McCormickR. L. (2023). Water contamination impacts on biodiesel antioxidants and storage stability. Energy & Fuels37, 5179–5188. 10.1021/acs.energyfuels.2c03911
14
CuiM.HouX.HanY.ZhangY.LiuZ.LiJ.et al (2025). Real-world emissions and formation mechanism of IVOCs from biodiesel-fueled agricultural machinery. Environ. Sci. Technol.59, 9017–9026. 10.1021/acs.est.4c11004
15
DavisM.AhiduzzamanMd.KumarA. (2018). How will Canada’s greenhouse gas emissions change by 2050? A disaggregated analysis of past and future greenhouse gas emissions using bottom-up energy modelling and sankey diagrams. Appl. Energy220, 754–786. 10.1016/j.apenergy.2018.03.064
16
de AguiarA. S. N.de CarvalhoL. B. R.GomesC. M.CastroM. M.MartinsF. S.BorgesL. L. (2025). Computational insights into the antioxidant activity of luteolin: density functional theory analysis and docking in cytochrome P450 17A1. Pharmaceuticals18, 410. 10.3390/ph18030410
17
de SousaL. S.de MouraC. V. R.de OliveiraJ. E.de MouraE. M. (2014). Use of natural antioxidants in soybean biodiesel. Fuel134, 420–428. 10.1016/j.fuel.2014.06.007
18
DuarteV. S.D. BorgesI.d’OliveiraG. D. C.FariaE. C. M.de AlmeidaL. R.Carvalho-SilvaV. H.et al (2023). Arylsulfonamide chalcones as alternatives for fuel additives: antioxidant activity and machine learning protocol studies. New J. Chem.47, 10003–10015. 10.1039/d3nj00255a
19
DuarteV. S.de PaulaR. L. G.AlmeidaL. R.D’OliveiraG. D. C.PérezC. N.CustódioJ. M. F.et al (2025). Exploring quinolinone–chalcones: synthesis, antioxidant potential and industrial applications in biofuels. Biofuels, Bioprod. Biorefining19, 1765–1783. 10.1002/bbb.2774
20
ErdemirA.LiS.JinY. (2005). Relation of certain quantum chemical parameters to lubrication behavior of solid oxides. Int. J. Mol. Sci.6, 203–218. 10.3390/i6060203
21
EtimA. O.JisieikeC. F.IbrahimT. H.BetikuE. (2022). “Biodiesel and its properties,” in Production of biodiesel from non-edible sources (Elsevier), 39–79. 10.1016/B978-0-12-824295-7.00004-8
22
FernandesD. M.SerqueiraD. S.PortelaF. M.AssunçãoR. M.MunozR. A.TerronesM. G. (2012). Preparation and characterization of methylic and ethylic biodiesel from cottonseed oil and effect of tert-butylhydroquinone on its oxidative stability. Fuel97, 658–661. 10.1016/j.fuel.2012.01.067
23
FrançaF. R. M.dos Santos FreitasL.RamosA. L. D.da SilvaG. F.BrandãoS. T. (2017). Storage and oxidation stability of commercial biodiesel using Moringa oleifera lam as an antioxidant additive. Fuel203, 627–632. 10.1016/j.fuel.2017.03.020
24
FregolenteP. B. L.FregolenteL. V.Wolf MacielM. R. (2012). Water content in biodiesel, diesel, and biodiesel–diesel blends. J. Chem. Eng. Data57, 1817–1821. 10.1021/je300279c
25
FreitasA. V.AlvesG. G. B.PaschoalG. M. A.Lafargue-dit-HauretW.HiornsR. C.BéguéD.et al (2024). A DFT bottom-up approach on non-fullerene acceptors: what makes highly efficient acceptors. J. Mater Sci.59, 10888–10903. 10.1007/s10853-024-09811-1
26
FukuiK. (1982). Role of frontier orbitals in chemical reactions. Science218, 747–754. 10.1126/science.218.4574.747
27
GanievB.MardonovU.KholikovaG. (2023). Molecular structure, HOMO-LUMO, MEP - – analysis of triazine compounds using DFT (B3LYP) calculations. Mater Today Proc.10.1016/j.matpr.2023.09.191
28
GawaiK. R.LokhandeP. D.KodamK. M.SoojhawonI. (2005). Oxidation of carbonyl compounds by whole-cell biocatalyst. World J. Microbiol. Biotechnol.21, 457–461. 10.1007/s11274-004-2467-y
29
GiakoumisE. G.RakopoulosC. D.DimaratosA. M.RakopoulosD. C. (2012). Exhaust emissions of diesel engines operating under transient conditions with biodiesel fuel blends. Prog. Energy Combust. Sci.38, 691–715. 10.1016/j.pecs.2012.05.002
30
hao NiZ.LiF. s.WangH.WangS.GaoS. y.ZhouL. (2020). Antioxidative performance and oil-soluble properties of conventional antioxidants in rubber seed oil biodiesel. Renew. Energy145, 93–98. 10.1016/j.renene.2019.04.045
31
HazratM. A.RasulM. G.KhanM. M. K.MofijurM.AhmedS. F.OngH. C.et al (2021). Techniques to improve the stability of biodiesel: a review. Environ. Chem. Lett.192209–2236. 10.1007/s10311-020-01166-8
32
HeB. B.ThompsonJ. C.RouttD. W.Van GerpenJ. H. (2007). Moisture absorption in biodiesel and its petro-diesel blends. Appl. Eng. Agric.23, 71–76. 10.13031/2013.22320
33
HohenbergP.KohnW. (1964). Inhomogeneous electron gas. Phys. Rev.136, B864–B871. 10.1103/physrev.136.b864
34
Hosseinzadeh-BandbafhaH.KumarD.SinghB.ShahbeigH.LamS. S.AghbashloM.et al (2022). Biodiesel antioxidants and their impact on the behavior of diesel engines: a comprehensive review. Fuel Process. Technol.232, 107264. 10.1016/j.fuproc.2022.107264
35
IleriE.KoçarG. (2014). Experimental investigation of the effect of antioxidant additives on NOx emissions of a diesel engine using biodiesel. Fuel125, 44–49. 10.1016/j.fuel.2014.02.007
36
JenkinM. E.HaymanG. D. (1999). Photochemical ozone creation potentials for oxygenated volatile organic compounds: sensitivity to variations in kinetic and mechanistic parameters. Atmos. Environ.33, 1275–1293. 10.1016/s1352-2310(98)00261-1
37
Kebi̇rogluH.AkF. (2023). Molecular structure, geometry properties, HOMO-LUMO, and MEP analysis of acrylic acid based on DFT calculations. J. Phys. Chem. Funct. Mater.6, 92–100. 10.54565/jphcfum.1343235
38
KimJ. H.ChanK. L.MahoneyN.CampbellB. C. (2011). Antifungal activity of redox-active benzaldehydes that target cellular antioxidation. Ann. Clin. Microbiol. Antimicrob.10, 23. 10.1186/1476-0711-10-23
39
KimB.MaX.ChenC.IeY.CoirE. W.HashemiH.et al (2013). Energy level modulation of HOMO, LUMO, and band‐gap in conjugated polymers for organic photovoltaic applications. Adv. Funct. Mater23, 439–445. 10.1002/adfm.201201385
40
KnotheG. (2007). Some aspects of biodiesel oxidative stability. Fuel Process. Technol.88, 669–677. 10.1016/j.fuproc.2007.01.005
41
KohnW.ShamL. J. (1965). Self-consistent equations including exchange and correlation effects. Phys. Rev.140, A1133–A1138. 10.1103/physrev.140.a1133
42
KreivaitisR.GumbytėM.KazancevK.PadgurskasJ.MakarevičienėV. (2013). A comparison of pure and natural antioxidant modified rapeseed oil storage properties. Ind. Crops Prod.43, 511–516. 10.1016/j.indcrop.2012.07.071
43
LapuertaM.Rodríguez-FernándezJ.RamosA.ÁlvarezB. (2012). Effect of the test temperature and anti-oxidant addition on the oxidation stability of commercial biodiesel fuels. Fuel93, 391–396. 10.1016/j.fuel.2011.09.011
44
LauC. H.LauH. L. N.NgH. K.Thangalazhy-GopakumarS.LeeL. Y.GanS. (2024). Evaluation of synthetic and bio-based additives on the oxidation stability of palm biodiesel: parametric, kinetics and thermodynamics studies. Sustain. Energy Technol. Assessments64, 103738. 10.1016/j.seta.2024.103738
45
LeeC. Y.SharmaA.SemenyaJ.AnamoahC.ChapmanK. N.BaroneV. (2020). Computational study of ortho-substituent effects on antioxidant activities of phenolic dendritic antioxidants. Antioxidants9, 189. 10.3390/antiox9030189
46
LiY.EvansJ. N. S. (1995). The Fukui function: a key concept linking frontier molecular orbital theory and the hard-soft-acid-base principle. J. Am. Chem. Soc.117, 7756–7759. 10.1021/ja00134a021
47
LuqueR.LovettJ. C.DattaB.ClancyJ.CampeloJ. M.RomeroA. A. (2010). Biodiesel as feasible petrol fuel replacement: a multidisciplinary overview. Energy Environ. Sci.3, 1706–1721. 10.1039/c0ee00085j
48
MendesR. A.da MataV. A. S.BrownA.de SouzaG. L. C. (2024). A density functional theory benchmark on antioxidant-related properties of polyphenols. Phys. Chem. Chem. Phys.26, 8613–8622. 10.1039/d3cp04412b
49
MittelbachM.SchoberS. (2003). The influence of antioxidants on the oxidation stability of biodiesel. J. Am. Oil Chem. Soc.80, 817–823. 10.1007/s11746-003-0778-x
50
NcP.KV.KgS.MnR.RA. R.JT.et al (2023). Quantum computations of non-steroidal anti-inflammatory drug molecules using density functional theory. Chem. Phys. Impact7, 100317. 10.1016/j.chphi.2023.100317
51
PantojaS. S.ConceiçãoL. R. V.DaCostaC. E. F.DaZamianJ. R.FilhoG. N. D. R. (2013). Oxidative stability of biodiesels produced from vegetable oils having different degrees of unsaturation. Energy Convers. Manag.74, 293–298. 10.1016/j.enconman.2013.05.025
52
ParrR. G.SzentpályL. V.LiuS. (1999). Electrophilicity index. J. Am. Chem. Soc.121, 1922–1924. 10.1021/ja983494x
53
Patiño-CaminoR.Cova-BonilloA.Rodríguez-FernándezJ.IglesiasT. P.LapuertaM. (2021). Relaxation dynamics of ethanol and N-Butanol in diesel fuel blends from terahertz spectroscopy. J. Infrared Millim. Terahertz Waves42, 772–792. 10.1007/s10762-021-00807-5
54
PearsonR. G. (1992). The electronic chemical potential and chemical hardness. J. Mol. Struct. (Theo&em)255, 261–270. 10.1016/0166-1280(92)85014-c
55
PearsonR. G. (2005). Chemical hardness and density functional theory. J. Chem. Sci.117, 369–377. 10.1007/bf02708340
56
PlácidoJ.CaparedaS. (2016). Conversion of residues and by-products from the biodiesel industry into value-added products. Bioresour. Bioprocess3, 23. 10.1186/s40643-016-0100-1
57
PullenJ.SaeedK. (2012). An overview of biodiesel oxidation stability. Renew. Sustain. Energy Rev.16, 5924–5950. 10.1016/j.rser.2012.06.024
58
PullenJ.SaeedK. (2014). Experimental study of the factors affecting the oxidation stability of biodiesel FAME fuels. Fuel Process. Technol.125, 223–235. 10.1016/j.fuproc.2014.03.032
59
RajamohanS.Hari GopalA.MuralidharanK. R.HuangZ.ParamasivamB.AyyasamyT.et al (2022). Evaluation of oxidation stability and engine behaviors operated by Prosopis juliflora biodiesel/diesel fuel blends with presence of synthetic antioxidant. Sustain. Energy Technol. Assessments52, 102086. 10.1016/j.seta.2022.102086
60
ReedA. E.WeinstockR. B.WeinholdF. (1985). Natural population analysis. J. Chem. Phys.83, 735–746. 10.1063/1.449486
61
Rizwanul FattahI. M.MasjukiH. H.KalamM. A.MofijurM.AbedinM. J. (2014). Effect of antioxidant on the performance and emission characteristics of a diesel engine fueled with palm biodiesel blends. Energy Convers. Manag.79, 265–272. 10.1016/j.enconman.2013.12.024
62
RochaM.Di SantoA.AriasJ. M.GilD. M.AltabefA. B. (2015). Ab-initio and DFT calculations on molecular structure, NBO, HOMO–LUMO study and a new vibrational analysis of 4-(Dimethylamino) benzaldehyde. Spectrochim. Acta A Mol. Biomol. Spectrosc.136, 635–643. 10.1016/j.saa.2014.09.077
63
RodriguesJ. S.do ValleC. P.UchoaA. F. J.RamosD. M.da PonteF. A. F.RiosM. A. d. S.et al (2020). Comparative study of synthetic and natural antioxidants on the oxidative stability of biodiesel from tilapia oil. Renew. Energy156, 1100–1106. 10.1016/j.renene.2020.04.153
64
RodriguezR. del G.ScanlonB. R.KingC. W.ScarpareF. V.XavierA. C.PruskiF. F. (2018). Biofuel-water-land nexus in the last agricultural frontier region of the Brazilian cerrado. Appl. Energy231, 1330–1345. 10.1016/j.apenergy.2018.09.121
65
RyuK. (2009). Effect of antioxidants on the oxidative stability and combustion characteristics of biodiesel fuels in an indirect-injection (IDI) diesel engine. J. Mech. Sci. Technol.23, 3105–3113. 10.1007/s12206-009-0902-6
66
SalujaR. K.KumarV.ShamR. (2016). Stability of biodiesel – a review. Renew. Sustain. Energy Rev.62, 866–881. 10.1016/j.rser.2016.05.001
67
Sanches-NetoF. O.Dias-SilvaJ. R.Keng Queiroz JuniorL. H.Carvalho-SilvaV. H. (2021). Py SiRC: machine learning combined with molecular fingerprints to predict the reaction rate constant of the radical-based oxidation processes of aqueous organic contaminants. Environ. Sci. Technol.55, 12437–12448. 10.1021/acs.est.1c04326
68
SaxenaV.KumarN.SaxenaV. K. (2017). A comprehensive review on combustion and stability aspects of metal nanoparticles and its additive effect on diesel and biodiesel fuelled C.I. engine. Renew. Sustain. Energy Rev.70, 563–588. 10.1016/j.rser.2016.11.067
69
SchoberS.MittelbachM. (2004). The impact of antioxidants on biodiesel oxidation stability. Eur. J. Lipid Sci. Technol.106, 382–389. 10.1002/ejlt.200400954
70
SerqueiraD. S.DornellasR. M.SilvaL. G.de MeloP. G.CastellanA.RuggieroR.et al (2015). Tetrahydrocurcuminoids as potential antioxidants for biodiesels. Fuel160, 490–494. 10.1016/j.fuel.2015.07.104
71
SerqueiraD. S.PereiraJ. F.SquissatoA. L.RodriguesM. A.LimaR. C.FariaA. M.et al (2021). Oxidative stability and corrosivity of biodiesel produced from residual cooking oil exposed to copper and carbon steel under simulated storage conditions: dual effect of antioxidants. Renew. Energy164, 1485–1495. 10.1016/j.renene.2020.10.097
72
SilvaF. R. daBaumgardt da SilvaF. J. L.VandresenF.da Silva LisboaF.SequinelR. (2025). Oxidative stability of biodiesel: challenges and perspectives for the sustainability of a large-scale program in Brazil. Biofuels16, 545–552. 10.1080/17597269.2024.2432154
73
SinghA.PrajapatiP.VyasS.GaurV. K.SindhuR.BinodP.et al (2022). A comprehensive review of feedstocks as sustainable substrates for next-generation biofuels. Bioenergy Res.16, 105–122. 10.1007/s12155-022-10440-2
74
SterpuA. E.SimedreaB. G.ChisT. V.SăpunaruO. V. (2024). Corrosion effect of biodiesel-diesel blend on different metals/alloy as automotive components materials. Fuels5, 17–32. 10.3390/fuels5010002
75
SuiM.ChenY.LiF.WangW.ShenJ. (2021). Study on the mechanism of auto-oxidation of jatropha biodiesel and the oxidative cleavage of C-C bond. Fuel291. 10.1016/j.fuel.2020.120052
76
SundusF.FazalM. A.MasjukiH. H. (2017). Tribology with biodiesel: a study on enhancing biodiesel stability and its fuel properties. Renew. Sustain. Energy Rev.70, 399–412. 10.1016/j.rser.2016.11.217
77
TakanoY.HoukK. N. (2005). Benchmarking the Conductor-like polarizable continuum model (CPCM) for aqueous solvation free energies of neutral and ionic organic molecules. J. Chem. Theory Comput.1, 70–77. 10.1021/ct049977a
78
ThuyP. T.SonN. T. (2022). Thermodynamic and kinetic studies on antioxidant capacity of amentoflavone: a DFT (density functional theory) computational approach. Free Radic. Res.56, 526–535. 10.1080/10715762.2022.2146584
79
UsmaniR. A.MohammadA. S.AnsariS. S. (2023). Comprehensive biofuel policy analysis framework: a novel approach evaluating the policy influences. Renew. Sustain. Energy Rev.183, 113403. 10.1016/j.rser.2023.113403
80
VaratharajanK.PushparaniD. S. (2018). Screening of antioxidant additives for biodiesel fuels. Renew. Sustain. Energy Rev.82, 2017–2028. 10.1016/j.rser.2017.07.020
81
WeinholdF.LandisC. R. (2001). Natural bond orbitals and extensions of localized bonding concepts. Chem. Educ. Res. Pract.2, 91–104. 10.1039/b1rp90011k
82
WeinholdF.LandisC. R. (2012). Discovering chemistry with natural bond orbitals. USA: John Wiley & Sons, Inc.
83
WeinholdF.LandisC. R.GlendeningE. D. (2016). What is NBO analysis and how is it useful?Int. Rev. Phys. Chem.35, 399–440. 10.1080/0144235x.2016.1192262
84
WenceslauP. R. S.AguiarA. S. N.DuarteV. S.de AlmeidaL. R.FrancoC. H. J.de AquinoG. L. B.et al (2025). Comprehensive analysis of pyrazoline analogs: exploring their antioxidant potential as biofuel additives. ACS Omega10, 40843–40856. 10.1021/acsomega.5c00398
85
XueY.ZhangL.LiY.YuD.ZhengY.AnL.et al (2013a). A DFT study on the structure and radical scavenging activity of newly synthesized hydroxychalcones. J. Phys. Org. Chem.26, 240–248. 10.1002/poc.3074
86
XueY.ZhengY.ZhangL.WuW.YuD.LiuY. (2013b). Theoretical study on the antioxidant properties of 2′-hydroxychalcones: H-Atom vs. electron transfer mechanism. J. Mol. Model19, 3851–3862. 10.1007/s00894-013-1921-x
87
YaakobZ.NarayananB. N.PadikkaparambilS.KS. U.PM. A. (2014). A review on the oxidation stability of biodiesel. Renew. Sustain. Energy Rev.35, 136–153. 10.1016/j.rser.2014.03.055
88
YamauchiM.KitamuraY.NaganoH.KawatsuJ.GotohH. (2024). DPPH measurements and structure—activity relationship studies on the antioxidant capacity of phenols. Antioxidants13, 309. 10.3390/antiox13030309
89
YangZ.HolleboneB. P.WangZ.YangC.LandriaultM. (2013). Factors affecting oxidation stability of commercially available biodiesel products. Fuel Process. Technol.106, 366–375. 10.1016/j.fuproc.2012.09.001
90
YusufN. N. A. N.KamarudinS. K.YaakubZ. (2011). Overview on the current trends in biodiesel production. Energy Convers. Manag.52, 2741–2751. 10.1016/j.enconman.2010.12.004
91
ZhangG.MusgraveC. B. (2007). Comparison of DFT methods for molecular orbital eigenvalue calculations. J. Phys. Chem. A111, 1554–1561. 10.1021/jp061633o
92
ZhaoY.TruhlarD. G. (2008). The M06 suite of density functionals for main group thermochemistry, thermochemical kinetics, noncovalent interactions, excited states, and transition elements: two new functionals and systematic testing of four M06-class functionals and 12 other functionals. Theor. Chem. Acc.120, 215–241. 10.1007/s00214-007-0310-x
Summary
Keywords
additives, antioxidant potential, biodiesel, fukui, stability
Citation
Borges ID, Aguiar ASN, Camargo AJ and Napolitano HB (2026) Biodiesel stabilization by dibrominated dimethoxybenzaldehydes: a comprehensive computational perspective. Front. Chem. Eng. 7:1716732. doi: 10.3389/fceng.2025.1716732
Received
30 September 2025
Revised
05 December 2025
Accepted
17 December 2025
Published
09 January 2026
Volume
7 - 2025
Edited by
Andre Luiz Da Silva, University of São Paulo, Brazil
Reviewed by
Ikbal Agah Ince, INSERM U1054 Centre de Biochimie Structurale de Montpellier, France
Mohan Govindasamy, University College of Engineering Villupuram, India
Updates
Copyright
© 2026 Borges, Aguiar, Camargo and Napolitano.
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: Hamilton B. Napolitano, hbnapolitano@gmail.com; Igor D. Borges, dalarmelino@ieee.org
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.