Electronic structure, NBO, and radical stability analyses of thymol and carvacrol along metabolic transformation pathways: a DFT study
Abstract
Understanding the proposed reaction mechanisms and radical stabilization processes leading to the formation of the two bioactive monoterpenes, thymol and carvacrol, remains challenging because of complex reaction networks, tautomeric interconversions, and electron delocalization effects. In this study, density functional theory (DFT) calculations were performed at the B3LYP/6-311 + G(d) level using Gaussian 09 software. Ten proposed reaction pathways and twenty transition states were identified and analyzed, including key tautomeric transformations and intermediate rearrangements. Activation energies and Gibbs free energies were evaluated to assess the feasibility of the proposed reaction routes. To gain deeper mechanistic insight, molecular electrostatic potential (MEP), Natural Bond Orbital (NBO), and spin density analyses were performed for selected intermediates, transition states, and radical species. MEP analysis revealed notable changes in charge distribution during tautomeric interconversions, while NBO analysis indicated strong donor–acceptor interactions and extensive π-electron delocalization within the aromatic framework. Spin density analysis confirmed resonance-assisted delocalization of the unpaired electron following hydrogen abstraction, supporting the stabilization of the corresponding phenoxy radicals. Comparative analysis revealed slightly greater electronic delocalization in thymol than in carvacrol. Overall, the present computational results provide molecular-level insight into the proposed reaction mechanisms, electronic structure evolution, and radical stabilization behavior associated with the formation of thymol and carvacrol.
Related
- Quaternion Dirac--Coulomb--Breit Integral Transformation for Relativistic Four-Component Correlated Electronic Structure Theory
- Thermodynamic stability and electronic structure of thorium boron carbides
- First-principles and cluster expansion investigation of Ti-V-Mn AB₂ alloys: stability, electronic structure, and mechanical properties.
- Manipulating the stability of crystallographic and magnetic sub-lattices: A first-order magnetoelastic transformation in transition metal based Laves phase
- VIBE-QC: an open-source, general-purpose electronic-structure code for molecules and solids developed with generative AI
BibTeX
@article{faal2026electronic,
title = {{Electronic structure, NBO, and radical stability analyses of thymol and carvacrol along metabolic transformation pathways: a DFT study}},
author = {Masoud Faal and Ali Mehrafarin and Hassanali Naghdi Badi and Sadegh Asadi},
journal = {Scientific Reports},
year = {2026},
month = oct,
doi = {10.1038/s41598-026-74640-y},
url = {https://www.nature.com/articles/s41598-026-74640-y},
note = {CC BY 4.0},
}
RIS
TY - JOUR
TI - Electronic structure, NBO, and radical stability analyses of thymol and carvacrol along metabolic transformation pathways: a DFT study
AU - Masoud Faal
AU - Ali Mehrafarin
AU - Hassanali Naghdi Badi
AU - Sadegh Asadi
PY - 2026
DA - 2026/10/05
JO - Scientific Reports
DO - 10.1038/s41598-026-74640-y
UR - https://www.nature.com/articles/s41598-026-74640-y
AB - Understanding the proposed reaction mechanisms and radical stabilization processes leading to the formation of the two bioactive monoterpenes, thymol and carvacrol, remains challenging because of complex reaction networks, tautomeric interconversions, and electron delocalization effects. In this study, density functional theory (DFT) calculations were performed at the B3LYP/6-311 + G(d) level using Gaussian 09 software. Ten proposed reaction pathways and twenty transition states were identified and analyzed, including key tautomeric transformations and intermediate rearrangements. Activation energies and Gibbs free energies were evaluated to assess the feasibility of the proposed reaction routes. To gain deeper mechanistic insight, molecular electrostatic potential (MEP), Natural Bond Orbital (NBO), and spin density analyses were performed for selected intermediates, transition states, and radical species. MEP analysis revealed notable changes in charge distribution during tautomeric interconversions, while NBO analysis indicated strong donor–acceptor interactions and extensive π-electron delocalization within the aromatic framework. Spin density analysis confirmed resonance-assisted delocalization of the unpaired electron following hydrogen abstraction, supporting the stabilization of the corresponding phenoxy radicals. Comparative analysis revealed slightly greater electronic delocalization in thymol than in carvacrol. Overall, the present computational results provide molecular-level insight into the proposed reaction mechanisms, electronic structure evolution, and radical stabilization behavior associated with the formation of thymol and carvacrol.
SN - 2045-2322
ER -