Home Journals Articles Updates About
CN
image
收藏切换
Reaction network database for fast pyrolysis of vanillyl alcohol based on molecular wavefunction descriptors
收藏切换
PDF
Jingyuan SHI, Chenchen LI, Xiaoxue CHENG, Qifan LING, Mao MU, Shuang WANG*, Ding JIANG*
Journal of Fuel Chemistry and Technology | 2026, 54(7) : 20250400-1458 - 20250400-1465
Less
收藏切换
Journal of Fuel Chemistry and Technology | 2026, 54(7): 20250400-1458-20250400-1465
• DATA PAPER •
Reaction network database for fast pyrolysis of vanillyl alcohol based on molecular wavefunction descriptors
Full
Jingyuan SHI, Chenchen LI, Xiaoxue CHENG, Qifan LING, Mao MU, Shuang WANG*, Ding JIANG*
Affiliations
  • School of Energy and Power Engineering, Jiangsu University, Zhenjiang 212013, China
Published: 2026-07-10 doi: 10.1016/S1872-5813(26)60710-4
Outline
收藏切换

Lignin pyrolysis is a promising avenue for the sustainable production of high-value phenolic chemicals. Nevertheless, unraveling the intricate radical reaction network remains a major bottleneck to optimizing product selectivity. This work constructs a standardized DFT computational database, which systematically describes the fast pyrolysis of vanillyl alcohol at 823.15 K. The database features out three key components, including primary reaction pathways, thermodynamic energy barriers, and atomic-level electronic fingerprints. Ultimately, the construction of this database not only provides an important reference for data-driven catalyst development, but also lays a theoretical foundation for the precise regulation of lignin depolymerization.

lignin  /  vanillyl alcohol  /  pyrolysis  /  phenols  /  DFT  /  木质素  /  香草醇  /  热解  /  酚类  /  DFT
史景元, 李晨琛, 程晓雪, 凌祺凡, 穆茂, 王爽, 蒋鼎. 基于分子波函数描述符的香草醇快速热解反应网络数据库. 燃料化学学报（中英文）, 2026 , 54 (7) : 20250400-1458 -20250400-1465 . DOI: 10.1016/S1872-5813(26)60710-4
Jingyuan SHI, Chenchen LI, Xiaoxue CHENG, Qifan LING, Mao MU, Shuang WANG, Ding JIANG. Reaction network database for fast pyrolysis of vanillyl alcohol based on molecular wavefunction descriptors[J]. Journal of Fuel Chemistry and Technology, 2026 , 54 (7) : 20250400-1458 -20250400-1465 . DOI: 10.1016/S1872-5813(26)60710-4
Full Text
Less
Lignin, the most abundant renewable aromatic polymer in nature, is an ideal feedstock for producing high-value phenolic compounds via pyrolysis[1−2]. However, the pyrolysis process is hindered by extremely low phenolic yields, typically only 10%−20%, thereby limiting its industrial-scale application[3]. This bottleneck is primarily attributed to the extremely complex radical reaction network during pyrolysis, in which complicated radicals inevitably trigger secondary reactions such as repolymerization[4]. Although the final product distribution has been extensively studied, there remain significant knowledge gaps regarding the evolution of key radical intermediates and their competitive reaction pathways[5].
Density Functional Theory (DFT) is a powerful tool for elucidating such complex reaction mechanisms[6]. However, constructing the potential energy surface of radical network that involves the major reaction pathways poses significant computational challenges. More importantly, single energy analysis (i.e., energy barriers) is often difficult to reveal the intrinsic driving forces for reactions. And electron-density-based reactivity descriptors, which are crucial for understanding chemical reactivity, remain relatively scarce in the existing databases. At present, there is a lack of a publicly available and standardized database that can systematically document the detailed reaction networks, energetic parameters, and key electronic structure descriptors of lignin model compounds under different reaction conditions of fast pyrolysis.
To address this gap, we constructed a multi-dimensional reaction network database for the fast pyrolysis of vanillyl alcohol[7], a representative lignin model compound, based on DFT calculations. This dataset covers the reaction pathways of vanillyl alcohol fast pyrolysis at 823.15 K, including the primary reaction pathways, the secondary rearrangements, and the both global and local reactivity indices of key intermediates. The comprehensive dataset provided herein not only elucidates the intrinsic limitations of the pyrolysis reactions, but also supplies high-quality descriptors for feature engineering, thereby facilitating the development of machine-learning-based reaction prediction models and the inverse design of novel pyrolysis catalysts.
1　Data summary
Less
1.1　Data profile
The dataset contains raw density functional theory (DFT) calculation data for the fast pyrolysis reaction network of vanillyl alcohol at 823.15 K (Table 1). It covers optimized geometries, vibrational frequencies, and electronic energies for all stationary points involved in the reaction network, including reactants, intermediates, transition states, and products. In addition, the dataset provides key secondary reaction pathways, initial bond dissociation energies, global reactivity descriptors, atom-resolved electronic fingerprints, spin population analyses of radical species, and molecular wavefunction-based descriptors for characterizing surface reactivity, such as Average Local Ionization Energy (ALIE), Local Electron Affinity Energy (LEAE), and Local Electron Attachment Energy (LEA).
1.2　Description of data generation methods and procedures
Database construction consists of two steps: raw data generation from DFT calculations, data extraction and data structuring. The detail workflow is as follows. First, the structures of reactants, key intermediates, transition states, and products involved in the fast pyrolysis of vanillyl alcohol (823.15 K) were built. Subsequently, DFT calculations were performed using Gaussian 16. Finally, the calculation output files were subjected to data extraction and preprocessing to generate the raw data package and the structured data. The operating details are described as follows:
Computational settings: All calculations were performed using the Gaussian 16 package with the B3LYP functional and the 6-311G (d, p) basis set.
All relevant species along the target reaction pathways, including reactants, intermediates, products, and those involved in secondary reaction pathways, were fully geometry-optimized. All optimized structures were subjected to vibrational frequency analysis to obtain Gibbs free energies at 823.15 K. All stable structures are presented as local minima on the potential energy surface. Each obtained transition state has only one imaginary frequency, indicating that it corresponds to a first-order saddle point. Furthermore, intrinsic reaction coordinate (IRC) calculations were performed to verify that each transition state correctly connects the corresponding reactant and the product on the one-dimensional potential energy profile[8−9].
Energy and barrier calculations: Free energies were obtained from frequency calculations with zero-point energy corrections. Reaction barriers were determined as the free energy differences between the transition states and the corresponding reactants.
Bond energy calculations: Key bond energies of the initial vanillyl alcohol species were calculated out based on the electronic energy differences.
Extraction of electronic structures and reactivity descriptors: These parameters were obtained based on the wavefunction analysis using the Multiwfn program. Notably, to keep the scaling consistent with that widely used in the literature, the nucleophilicity index was computed by referencing the energy of tetracyanoethylene (TCE) at the B3LYP/6-31G* level as implemented in Multiwfn program[10−12].
1.3　Description of data examples
This dataset provides multi-dimensional reactivity data spanning from macroscopic thermodynamic potential energy surfaces to microscopic electronic structures, in order to characterize the fast pyrolysis behavior of vanillyl alcohol. To facilitate rapid indexing and reusing, this section first describes the file organization of the dataset and the workflow of data generation. Then, the detailed data using key intermediates as examples are presented.
1.3.1　Dataset generation workflow and organizational structure
Computational workflow: This database uses DFT as the underlying data-generation approach and follows the computational procedure as shown in Figure 1, which mainly consists of the following three core steps.
Structure search and optimization: Based on the Py-GC/MS data of the fast pyrolysis of vanillyl alcohol, initial structures of key species involved in the reaction network were systematically constructed (as shown in Figure 3(a) and (c)). Subsequently, geometry optimizations and frequency calculations were carried out at the B3LYP/6-311G (d, p) level.
Multi-electronic-state wavefunction generation: Based on the optimized equilibrium geometries, single-point energy calculations were performed for the neutral (N), cationic (N−1) and anionic (N+1) systems, and the corresponding wavefunction files were output.
Descriptor analysis: The wavefunction files were batch processed using Multiwfn to systematically compute (i) the descriptors covering global reactivity and atom-level local electronic-structure, (ii) the spin population analysis for radicals, and (iii) the average local ionization energy (ALIE), local electron affinity energy (LEAE), and local electron attachment energy (LEA) for characterizing molecular surface activity.
Data file organizational structure: This database adopts a standardized hierarchical directory structure for storage to ensure data readability and traceability. As shown in Figure 2, the dataset is organized according to a hierarchical logic of “reaction pathway-species node”. In this framework, each structure (e.g., “SR-1”) is assigned an independent storage unit. The unit specifically includes the following information:
Gaussian calculation records: These files include output log files (.log) and checkpoint files (.chk), which contain complete molecular geometries and thermodynamic data.
Wavefunction analysis files: These files include wavefunction files (.wfn) for the neutral (N), cationic (N−1) and anionic (N+1) states, which are used for subsequent detailed analysis of electronic properties.
1.3.2　Secondary cracking reaction network
The final product distribution of vanillyl alcohol pyrolysis largely depends on the secondary evolution of the primary products. To elucidate the complex competition mechanism among radical reaction pathways, this section highlights representative secondary cracking pathways. As shown in Figure 3, this database not only covers the formation pathways of the primary products (Figure 3(a) and (b)) but also emphases on the secondary cracking processes of key intermediates. Figure 3(c) constructs a reaction network centered on the consumption of hydrogen radicals, whereas Figure 3(d) provides the corresponding data of potential energy surface. These results systematically document the thermodynamic evolution from the intermediates to the final small-molecule products.
1.3.3　Electronic-structure characteristics of key species
To facilitate an in-depth investigation of the intrinsic factors governing reaction kinetics, this database provides global electronic fingerprints for all key species along the pathways as shown in Figure 3. Tables 2 and 3 present an example for the pathway in Figure 3(c). These features (including energy, chemical hardness, electrophilicity and nucleophilicity indices, average local ionization energy, local electron affinity energy, and local electron attachment energy, etc.) constitute the basis for understanding molecular reactivity. Based on these standardized data, users can not only evaluate the thermodynamic stability of intermediates but also quantitatively interpret the electronic origin for the competitive advantages among different reaction pathways (e.g., why SR-1 tends to follow a specific rearrangement pathway).
1.3.4　Analysis of the local active sites
To facilitate mechanistic analysis of radical reactions at the microscopic level, this database further provides atom-resolved electronic descriptors. Table 4 takes the intermediate SR-1 as an example. The atomic numbering is shown on the left of Figure 4, and the table presents a multi-dimensional local feature matrix that includes Fukui indices and atomic charges. Table 5 reports the spin populations for intermediate SR-2, whose atomic numbering is displayed on the right of Figure 4. These data enable effective identification of intramolecular active sites and help to explain the regioselectivity rules of radical attacks. In addition, these standardized features also provide essential training data for constructing high-accuracy machine-learning potentials.
1.4　Data quality control statement
The quality control of this dataset is ensured through the following measures:
Standardized computational workflow: All calculations in this database strictly follow unified computational parameters and procedures to ensure the comparability of data across all the species. Specifically, all calculations were performed using the Gaussian 16 package at the B3LYP/6-311G (d, p) level of theory. All thermodynamic parameters were obtained at 823.15 K.
Structure and convergence verification: All the computational tasks were confirmed to be converged. All stationary points were validated by frequency calculations. Stable structures were confirmed as minima, and transition states were confirmed as first-order saddle points.
Based on the IRC calculations verification, each transition state correctly connected the corresponding reactant and product on the one-dimensional potential energy profile, thereby ensuring the accuracy of the computational results.
Data consistency and scaling validation: All key data in this dataset were extracted from the original computational output files, either automatically or manually, and were cross-checked to ensure consistency between the summarized data and the raw data, which effectively avoided transcription errors. For electronic-structure descriptors, a specific clarification was provided. To maintain consistency with the scaling widely adopted in the literature, the TCE at the B3LYP/6-31G* level was used as the reference when calculating the nucleophilicity index.
2　Value and potential applications of the dataset
Less
This dataset establishes the first multi-dimensional database for the fast pyrolysis of lignin model compounds that integrates thermodynamic barriers and electronic-structure descriptors, addressing the lack of a systematic “reaction network-electronic property” correlation map in existing resources. The core value and specific application scenarios of this database can be summarized in the following four aspects:
A theoretical benchmark for fast pyrolysis mechanistic studies: This dataset systematically provides thermodynamic data for vanillyl alcohol under fast pyrolysis conditions. It can serve as a benchmark dataset for mechanistic studies of lignin model compound pyrolysis, enabling the comparisons of reaction pathways across different catalytic systems. It also allows quantitative evaluation of the catalytic efficiencies for specific elementary reactions, thereby assisting in an in-depth elucidation of catalytic mechanisms.
Reactivity analysis based on electronic effects: Unlike traditional databases that contain only energetic parameters, this dataset integrates key electronic descriptors such as chemical hardness and nucleophilicity/electrophilicity indices, which provides a new perspective for understanding the competition mechanism among radical pathways.
In practical applications, researchers can use these data to interpret regioselectivity rules, and such electron-density-based analysis can reveal out the kinetic phenomena that cannot be clarified by conventional thermodynamic analysis. This knowledge provides a microscopic basis for tuning reaction pathways.
Data-driven machine-learning model development: With a standardized data structure and atom-level feature labels, this dataset is an ideal resource for training high-accuracy machine-learning models. First, this dataset can be used to develop reactivity prediction models, enabling predictions for lignin oligomer systems with more complex structure. Second, based on the provided high-fidelity configurations and energy data, machine-learning potentials can be trained, thereby enabling reactive molecular dynamics simulations at the nanosecond timescale. Finally, by integrating the multi-dimensional electronic descriptors in this dataset, researchers can build quantitative models to reveal out the relationship between structure and property of “electronic microenvironment-product selectivity” , which provides key theoretical criteria for the rational screening and inverse design of efficient lignin depolymerization catalysts.
Catalyst design: The electrophilicity/nucleophilicity indices in this database can be used to guide the design of catalyst active sites for the pyrolysis of lignin model compounds. Taking the vanillyl alcohol pyrolysis system as an example, it is highly dependent on a hydrogen-radical-mediated chain reaction network, and an ideal catalyst design should promote specific pathways while minimizing disturbances to the hydrogen radical concentration in the system or the excessive consumption of key intermediates. Using the electrophilicity/nucleophilicity indices provided by the database, researchers can predict the adsorption tendencies of intermediates on the active sites of catalyst surface. For example, for radical intermediates with high nucleophilicity indices, catalyst materials featuring strongly electrophilic surface active sites should be avoided in the design process to prevent radical quenching or non-selective deep cracking caused by overly strong adsorption.
3　Conclusion
Less
Based on DFT, this work establishes a multi-dimensional reaction network database for the fast pyrolysis of the lignin model compound vanillyl alcohol at 823.15 K. This database not only involves complex radical reaction pathways and thermodynamic energy barriers, but also innovatively integrates electronic-structure fingerprints of key intermediates, which fills the gap in “reaction network-electronic property” correlation data in this field. Standardized computational workflows and rigorous quality control were employed to ensure the accuracy, consistency, and reproducibility of the data. The public release of this dataset will not only provide a reliable theoretical benchmark for mechanistic studies of lignin pyrolysis, but also offer valuable foundational data support for the rational design of new catalysts and the refinement of reaction kinetic models.
References
Less
1
OKONKWO C E, HUSSAIN S Z, ONYEAKA H, et al. Lignin polyphenol: From biomass to innovative food applications, and influence on gut microflora[J]. Ind Crops Prod, 2023, 206: 117696.
2
GUI C, WANG L, LIU G, et al. The catalytic valorization of lignin from biomass for the production of liquid fuels[J]. Energies, 2025, 18(6): 1478.
3
JUNG K A, WOO S H, LIM S -R, et al. Pyrolytic production of phenolic compounds from the lignin residues of bioethanol processes[J]. Chem Eng J, 2015, 259: 107−116.
4
BAI X, KIM K H, BROWN R C, et al. Formation of phenolic oligomers during fast pyrolysis of lignin[J]. Fuel, 2014, 128: 170−179.
5
HOSOYA T, KAWAMOTO H, SAKA S. Secondary reactions of lignin-derived primary tar components[J]. J Anal Appl Pyrolysis, 2008, 83(1): 78−87.
6
CUI D, ZHANG M, YIN H, et al. Study of the mechanism of substituent effects in lignin pyrolysis using density functional theory[J]. J Energy Inst, 2025, 119: 102009.
7
LIU C, DENG Y, WU S, et al. Study on the pyrolysis mechanism of three guaiacyl-type lignin monomeric model compounds[J]. J Anal Appl Pyrolysis, 2016, 118: 123−129.
8
FUKUI K. Formulation of the reaction coordinate[J]. J Phy Chem, 1970, 74(23): 4161−4163.
9
MAEDA S, HARABUCHI Y, ONO Y, et al. Intrinsic reaction coordinate: Calculation, bifurcation, and automated search[J]. Int J Quantum Chem, 2015, 115(5): 258−269.
10
LU T, CHEN F. Multiwfn: A multifunctional wavefunction analyzer[J]. J Comput Chem, 2012, 33(5): 580−592.
11
LU T. A comprehensive electron wavefunction analysis toolbox for chemists, Multiwfn[J]. J Chem Phys, 2024, 161(8).
12
LU T, CHEN Q. Realization of conceptual density functional theory and information‐theoretic approach in Multiwfn Program[M]. Conceptual Density Functional Theory: Towards a New Chemical Reactivity Theory (Vol. 2), Weinhein: Wiley-VCH, 2022.
13
JIANG D, LI H, CHENG X, et al. Catalytic role of biochar surface groups in phenol formation from lignin pyrolysis: Experimental and DFT insights[J]. Chem Eng J, 2025, 521: 166613.
Appendix
Less
Year 2026 volume 54 Issue 7
PDF
133
56
Cite this Article
BibTeX
Article Info
doi: 10.1016/S1872-5813(26)60710-4
  • Receive Date：2025-12-10
  • Online Date：2026-09-23
  • Published：2026-07-10
Article Data
Affiliations
History
  • Received：2025-12-10
  • Revised：2026-01-23
  • Accepted：2026-02-25
Affiliations
    School of Energy and Power Engineering, Jiangsu University, Zhenjiang 212013, China

Corresponding:

Tel: 0511-88780215, E-mail:
References
Share
https://castjournals.cast.org.cn/joweb/rlhxxb/EN/10.1016/S1872-5813(26)60710-4
Share to
QR

Scan QR to access full text

Cite this article
BibTeX
Citations
表12种不同金属材料的力学参数
科
Family		 属数
Number of
genus		 种数
Number of
species		 占总种数比例
Percentage of
total species (%)		 属
Genus		 种数
Number of
species		 占总种数比例
Percentage of total
species (%)
鹅膏菌科Amanitaceae 2 11 5.26 鹅膏菌属 Amanita 10 4.78
小菇科 Mycenaceae 2 12 5.74 丝盖伞属 Inocybe 5 2.39
多孔菌科 Polyporaceae 8 14 6.70 蜡蘑属 Laccaria 5 2.39
红菇科 Russulaceae 3 23 11.00 小皮伞属 Marasmius 6 2.87
小菇属 Mycena 11 5.26
光柄菇属 Pluteus 5 2.39
红菇属 Russula 17 8.13
栓菌属 Trametes 5 2.39
关闭全屏
  • BibTeX
  • EndNote
  • RefWorks
  • TxT
Links
Articles
Latest Articles
Most Read
Collections
Updates
Events
News
Multimedia
About
About Us
Contact
No. 86 Xueyuan South Road, Haidian District, Beijing
100081
010-62662699(Operation and Technical Support)
qkjq@cast.org.cn

Copyright © 2025 China Association for Science and Technology. All rights reserved. For all open access content, the relevant licensing terms apply.
Sponsored by the Office of the Leading Group for Cybersecurity and Informatization of CAST, and supported by Science and Technology Review Publishing House