Next Article in Journal
Additively Manufactured Multifunctional Composite Parts with the Help of Coextrusion Continuous Carbon Fiber: Study of Feasibility to Print Self-Sensing without Doped Raw Material
Next Article in Special Issue
A Comparison of Three Simulation Techniques for Modeling the Fan Blade–Composite Abradable Rub Strip Interaction in Turbofan Engines
Previous Article in Journal
Powder Metallurgy Preparation and Characterization of Titanium-Titanium Diboride Composite Targeted for Dental Implant
Previous Article in Special Issue
Mechanical Properties of Uncured Thermoset Tow Prepreg: Experiment and Finite Element Analysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Review on Characterization of Biochar Derived from Biomass Pyrolysis via Reactive Molecular Dynamics Simulations

1
Mechanical Engineering Department, South Dakota State University, Brookings, SD 57007, USA
2
Agricultural and Biosystems Engineering Department, South Dakota State University, Brookings, SD 57007, USA
*
Author to whom correspondence should be addressed.
J. Compos. Sci. 2023, 7(9), 354; https://doi.org/10.3390/jcs7090354
Submission received: 31 July 2023 / Revised: 20 August 2023 / Accepted: 23 August 2023 / Published: 25 August 2023
(This article belongs to the Special Issue Characterization and Modelling of Composites, Volume III)

Abstract

:
Biochar is a carbon-rich solid produced during the thermochemical processes of various biomass feedstocks. As a low-cost and environmentally friendly material, biochar has multiple significant advantages and potentials, and it can replace more expensive synthetic carbon materials for many applications in nanocomposites, energy storage, sensors, and biosensors. Due to biomass feedstock species, reactor types, operating conditions, and the interaction between different factors, the compositions, structure and function, and physicochemical properties of the biochar may vary greatly, traditional trial-and-error experimental approaches are time consuming, expensive, and sometimes impossible. Computer simulations, such as molecular dynamics (MD) simulations, are an alternative and powerful method for characterizing materials. Biomass pyrolysis is one of the most common processes to produce biochar. Since pyrolysis of decomposing biomass into biochar is based on the bond-order chemical reactions (the breakage and formation of bonds during carbonization reactions), an advanced reactive force field (ReaxFF)-based MD method is especially effective in simulating and/or analyzing the biomass pyrolysis process. This paper reviewed the fundamentals of the ReaxFF method and previous research on the characterization of biochar physicochemical properties and the biomass pyrolysis process via MD simulations based on ReaxFF. ReaxFF implicitly describes chemical bonds without requiring quantum mechanics calculations to disclose the complex reaction mechanisms at the nano/micro scale, thereby gaining insight into the carbonization reactions during the biomass pyrolysis process. The biomass pyrolysis and its carbonization reactions, including the reactivity of the major components of biomass, such as cellulose, lignin, and hemicellulose, were discussed. Potential applications of ReaxFF MD were also briefly discussed. MD simulations based on ReaxFF can be an effective method to understand the mechanisms of chemical reactions and to predict and/or improve the structure, functionality, and physicochemical properties of the products.

1. Introduction

Biochar is a carbon-rich solid material that can be produced by oxygen-limited thermochemical processes (e.g., torrefaction, pyrolysis, gasification, combustion, etc.) [1,2,3,4,5,6]. Many types of organic materials, including agricultural residues and forest solid waste, can be used as feedstock for biochar production. This biochar usually contains carbon (C) (40–90%), oxygen (O) (1–40%), hydrogen (H) (1–10%), nitrogen (N) (0–5%), sulfur (S) (<1%), and other trace elements [7,8,9]. Depending on biomass feedstock species, reactor types, and operating conditions, such as particle size, heating rate, residence time, carrier gas, reaction temperature, pressure, etc., the proportions and physicochemical properties of different compounds of biochar may vary greatly. Biochar has various significant advantages and potentials due to its remarkable properties [10,11,12,13]. For thousands of years, it has been used as a soil amendment to improve agricultural productivity, but recently, the use of biochar as a precursor, which is activated or functionalized into carbon materials, such as graphite, graphene, carbon black, or synthetic carbon fibers (CFs), carbon nanofibers (CNFs), carbon nanotubes (CNTs), carbon nanospheres (CNSs), etc., has become more and more attractive for many applications in nanocomposites, energy storage, sensors, biosensors, etc. Carbon materials produced from non-renewable sources such as coal, oil, and natural gas may produce excessive greenhouse gas (GHG) emissions, accelerate the consumption of fossil fuel resources, and are not economically viable. In contrast, biomass resources are abundant, available, renewable, low-cost, harmless, high in carbon content, and low in sulfur and ash [10,11,12,13,14,15,16,17]. In addition, applications of biochar can promote carbon sequestration and reduce CO2 emissions into the atmosphere, improving resilience to global climate change [10,15]. Cancer-targeted drug delivery systems based on carbon nanostructures are another promising application due to their ability to selectively recognize specific receptors overexpressed in cancer cells [18,19]. Therefore, various methods, such as laser ablation, gamma irradiation, chemical vapor deposition (CVD), arc-discharge, catalytic chemical vapor deposition (CCVD), hydrothermal carbonization [18,19], plasma technique [8,20], and plasma-enhanced chemical vapor deposition (PECVD), have been developed to synthesize different carbon nanomaterials (e.g., CNTs) [14]. Thus, the biochar market is expected to reach USD 3.14 billion by 2025 [14,15,21,22].
Biochar is a complex carbonaceous material with many physical and chemical parameters that control its reactivity towards inorganic and organic substances in aqueous solution [23]. Some chemical properties of biochar vary with the type of parent biomass feedstock species, with the most pronounced and consistent changes in volume and surface chemistry occurring with production conditions, biochar surface charge and ion exchange capacity, environmental influences, etc. [24]. Biochar has also attracted considerable attention as a filler material in polymer matrices and as a reinforcement. The fabrication of bio-composites employs a variety of methods. The effect of adding biochar to the overall composites shows great promise in improving the performance of the overall composites [25]. Various modification methods were proposed to enhance certain functions of biochar. However, these modifications may also lead to structural uncertainty, additional energy consumption, secondary pollution, and/or additional costs. The most commonly used biochar modification methods can be classified according to the purposes of modification, such as surface area enlargement, persistent free radical manipulation, magnetism introduction, and redox potential enhancement. More importantly, the balance considerations of biochar designs, such as the balance between effectiveness and stability, function and risk, effectiveness and cost, etc., were systematically analyzed and discussed [26,27]. To determine the synergistic effects of biomass feedstock species, reactor types, operating conditions, and post-treatments on biochar activation and functionalization, the surface chemistry, physiochemistry, structural, and molecular characterization of the produced biochar, as well as the relationships between process parameters and biochar surface morphology and internal microstructure, have been intensively studied since the last few decades [28,29,30]. The surface chemistry, structural, and molecular characterization of the produced biochar have been investigated using various advanced characterization techniques. For example, the chemical characterization of biochar can be tested using elemental analysis, Fourier transform infrared spectroscopy (FTIR), X-ray photoelectron spectroscopy (XPS), and scanning electron microscope—energy dispersive spectrometer (SEM-EDS), as well as statistical analysis of experimental data [31]. Regarding surface and structural characterization, biochar and activated carbons can have varying degrees of aromaticity, which generally refers to the level of defects and the size of the aromatic rings in the bulk material. Since micropores formed by stacking faults between layers of aromatic ring clusters in porous carbons are responsible for the high effective surface area in these materials, it is crucial to characterize these phases formed of sp2 hybridized carbons. Raman spectroscopy (Raman), scanning electron microscopy (SEM), specific surface area and pore size analyzer, high-resolution transmission electron microscopy (HR-TEM), X-ray diffraction (XRD), XPS, etc., provide simple and non-destructive means to achieve this goal, and diffuse reflectance infrared Fourier transform spectroscopy (DRIFTS) can be used semi-quantitatively [32,33]. To characterize the structure of biochar at the molecular level, HR-TEM, Raman, XPS, etc., can be used [30,34]. The 3D morphology and distribution of biochar pores can be studied by using a density analyzer, X-ray computed tomography (XCT), and SEM [35,36,37]. In summary, many modern characterization techniques have been reported for the characterization of biochar, such as SEM, FTIR, XRD, thermogravimetric analysis (TGA), nuclear magnetic resonance spectroscopy (NMR), Brunauer–Emmett–Teller (BET), Raman, density analyzer, proximate and ultimate analysis, etc. [38,39].
Biochar can be produced from biomass using various methods. Biochar physicochemical properties, including chemical composition, carbon content, ash content, surface area and morphology, pore size and distribution, functionality, etc., may be affected by the biomass feedstock species nature, pretreatment, pyrolysis conditions (such as reaction temperature, heating rate, retention time), functionalization strategies (e.g., magnetic biochar, plasticized biochar, and co-composed biochar), and so on. The physical chemistry of biochar functioning includes the thermodynamics of biochar adsorption/desorption, the kinetics of biochar adsorption/desorption, the meaning of reaction order and Langmuir isotherm from kinetic considerations, the dynamics of water and nutrients in the biochar pore system, and the trapping and decomposition mechanisms of pollutants in biochar [40,41,42,43]. However, understanding biochar seems more experimental and empirical than relying on a well-structured theoretical framework. From previous studies, the properties and functionalities of biochar depend on many factors, including the composition, morphology, size, internal structure, pore size and distribution of biochar clusters, processing methods, etc. However, there is still a lack of quantitative analysis and a comprehensive understanding of the effects of these factors and their interactions on the physicochemical properties and functionalities of biochar, which are critical for producing high-quality biochar as a precursor to effectively developing innovative carbon materials for new applications. The effects of these factors are usually not in a monotonically changing pattern (monotonically increasing or monotonically decreasing). There is usually an inflection point (threshold) where the material performance is at a maximum or minimum point, which can be fully considered when conducting biochar production and effective new/novel carbon material designs.
Depending on the problem and the spatial and time scales of interest, various computer modeling-based approaches to materials design/research exist, ranging from quantum mechanics (QM) to continuum simulations. Molecular dynamics (MD) or first-principles simulations are ideal for studying the properties of materials at the nanoscale. MD is an atomistic-scale simulation that describes the interactions between atoms through interatomic potentials. In the MD method, electronic effects are averaged, and the time evolution of atomic positions and velocities is calculated according to Newton’s equations of motion. The electron correlation approximation is based on the Born–Oppenheimer theory that the MD time step used to describe the atomic motion is sufficient for electrons to achieve their ground stable state compared to the nuclei due to the difference in mass. Interatomic potentials (force fields) are established from the first principles or experimentally to describe the interactions between the atoms, including the effect of electrons, in terms of reproducible forces. The reliability of the interatomic potential determines the accuracy of the MD simulations and is also related to the ability to bridge the mesoscale methods [44,45,46,47,48,49]. With the development of MD in the past few decades, MD has developed into a ubiquitous, versatile, and powerful computational method for basic scientific research in biology, chemistry, biomedicine, physics, etc. [50]. Driven by the rapid development of supercomputing technology in recent decades, MD has entered the field of engineering as a first-principles predictive method for material properties, physicochemical processes, and even as a design tool. These developments have far-reaching implications and are discussed in recent papers focusing on MD of combustion and energy systems, gas/liquid/solid fuel oxidation, pyrolysis, catalytic combustion, electrochemistry, nanoparticle synthesis, heat transfer, phase change, and fluid mechanics. Due to the practical availability of MD simulations of large-scale reactive chemical systems, the reactive force field (ReaxFF) system was developed and successfully applied to biomass systems [29,30,50,51,52,53,54,55,56,57,58,59]. ReaxFF uses a general relationship between chemical bond distance and bond order to separate atoms. The other valence terms (angular and torsional) present in the force field are defined in terms of the same bond orders so that all these terms smoothly go to zero when the bonds break. In addition, ReaxFF has Coulomb and Morse (van der Waals) potentials to describe non-bond interactions between all atoms. These non-bond interactions are shielded at short distances, so Coulomb and van der Waals interactions become constant as the distances approach zero. The parameters of the ReaxFF were derived from quantum chemistry (QC) calculations on bond dissociation and reactions of small molecules, as well as heat of formation and geometric data for some stable hydrocarbons [51].
The current theoretical framework of MD methodology covers both classical and reactive MD. The ReaxFF MD can simulate chemical reactions such as carbonization in biomass pyrolysis within the MD framework using QC calculations and/or force field representations of experimental data. Therefore, this paper reviews the fundamentals of the ReaxFF methods and MD simulations based on ReaxFF MD to characterize the physicochemical properties and biomass pyrolysis process. The simulation of the carbonization reaction in biomass pyrolysis, including the reactivity of the three major components: cellulose, lignin, and hemicellulose, was discussed. The potential applications of ReaxFF MD were also briefly discussed.

2. Fundamentals of ReaxFF MD Simulations

Modeling atomic and molecular systems requires computationally intensive QM methods such as, but not limited to, QM-based ab initio calculations [60] or density functional theory (DFT) [46,48,58,61,62,63]. These methods have successfully predicted various properties of atomic-scale chemical systems [53]. Due to the inherent non-locality of QM, the scalability of these methods ranges from O(N3) to O(N7), depending on the method used and the approximations involved. This greatly limits the size of simulated systems to a few thousand atoms, even on massively parallel platforms. Furthermore, the scaling of the computational cost as a function of the system size (~O(N3) for the DFT methods and ~O(N) for the ReaxFF method) still limits the possibility of running QM or ReaxFF simulations long enough to produce large carbonization structures [64]. On the other hand, classical approximations of quantum systems, while computationally (relatively) easy to implement, produce simpler models that lack fundamental chemical properties such as reactivity and charge transfer. The recent work by van Duin et al. [51,52,53,54,55,56,57,58,59,60,61,62,63,64] carefully incorporated limited non-locality (to simulate quantum behavior) via empirical bond order potentials, overcoming the limitations of non-reactive classical MD approximations while largely preserving the computational simplicity of classical MD. The reactive system force field, ReaxFF, was developed. ReaxFF uses a general relationship between bond distances and bond orders on the one hand, and bond orders and bond energies on the other, resulting in the proper dissociation of bonds into separated atoms. The other valence terms (angular and torsional) present in the force field are defined in terms of the same bond orders so that all these terms smoothly go to zero when the bonds break. In addition, ReaxFF has Coulomb and Morse (van der Waals) potentials to describe all non-bonded interactions between all atoms.

2.1. General Force Field

The ReaxFF interatomic potentials are powerful computational tools for exploring, developing, and optimizing material properties. Although methods based on QM principles provide valuable theoretical guidance at the electronic level, they are often too computationally intensive for simulations that consider the complete dynamic evolution of a system. Like empirical non-reactive force fields, reactive force fields divide the system’s energy into various partial energy contributions. Russo and van Duin [53] described the ReaxFF potential forms in detail. Therefore, only a brief overview of the central concepts of the method is presented here. ReaxFF uses a bond-order form combined with a description of polarizable charges between atoms. This enables ReaxFF to accurately model the covalent and electrostatic interactions of various materials, as shown in Figure 1.
The energy contributions to the ReaxFF potentials are summarized below [49,51,55,57,58]:
E s y s t e m   =   E b o n d + E o v e r + E a n g l e + E t o r s + E v d W a a l s + E C o u l o m b + E s p e c i f i c
where the contributions to the total system energy, Esystem, are the atomic bonds (continuous functions of interatomic distances, describing the energies associated with forming bonds between atoms), atomic over-coordination penalties (energy penalties that preventing atoms from over-coordination, based on valence rules, e.g., stiff energy penalties are applied if carbon atoms form more than four bonds), atomic valence angle and torsion (the energies related to three-body valence angle strain and four-body torsional angle strain), atomic non-bonding Coulombic and van der Waals energies (electrostatic and dispersive contributions calculated between all atoms, regardless of connectivity and bond-order), and the atomic specific term (system specific term that is generally not included unless required to capture properties particular to the system of interest, such as atomic lone-pairs, atomic conjugation, hydrogen binding, and C2 corrections).

2.2. Bond Order and Bond Energy

A basic assumption of ReaxFF is that the bond order BOij between atomic pairs can be obtained directly from the potential energy, which is divided into bond-order-dependent and bond-order-independent contributions. The bond order is calculated directly from the interatomic distances using an empirical formula:
B O i j = B O i j π + B O i j π π = e x p p b o , 1 r i j r o σ p b o , 2 + e x p p b o , 3 r i j r o π p b o , 4 + e x p p b o , 5 r i j r o π π p b o , 6
where BOij is the bond order between atoms i and j, rij is the interatomic distance, ro is the equilibrium bond length, and Pbo is the empirical parameter. In Equation (2), the first, second, and third exponential terms determine the bond order contributions of single bond (σ-bond), double bond (π-bond), and triple bond (double π-bond). Each bonding term, p, and each bonding equilibrium distance, ro, are parameterized to yield bond strengths and distances consistent with quantum mechanical predictions for species separated by distance rij [51,53,65,66,67], as shown in Figure 2, where two carbon atoms are considered. The equilibrium distance, ro, is about 0.154 nm. This results in a maximum bond order level of 3 for C–C. For C–H and H–H bonds, only the σ-bond contribution is considered, resulting in a maximum bond order level of 1. All bond-order functions fall smoothly to zero.
The classical MD simulations use force field potentials to describe the time evolution of any particle (atoms, molecules, or radicals) for a set of initial conditions based on Newton’s second law (Newton’s equations of motion). For atoms and/or molecular systems, the force field potentials (as uncorrected bond order in Figure 2) are calculated for each pair of atoms in the simulation as functions of the distance between atoms only. Solving Newton’s equations of motion for a system constitutes the central task of classical MD, just as solving the Schrödinger equation is the central task of QM. However, in classical MD, subatomic electronic structures and dynamics are not calculated, thus excluding intrinsic QM events such as chemical reactions (i.e., chemical bond breaking and formation) and charge transfer. Whereas reactive MD seeks to balance power (chemical reactivity), accuracy, and computational cost. The reactive force field, particularly the ReaxFF method, expresses the energy terms (Equation (1)) as functions of bond order (or bond distance), such as σ-bond, π-bond, double π-bond, and the summation of the three bond orders—corrected bond order, as shown in Figure 2). The energy functions create a continuous surface connecting the reactants and produce and are generated through all intermediate states, which makes ReaxFF useful for simulating chemical reactions [50]. Based on the comparison of the schematic representations of the bond orders by the classical MD (uncorrected in Figure 2) and the ReaxFF MD (corrected in Figure 2), the difference, especially around the equilibrium bond distance, has been clearly shown.

2.3. Atomic Over-Coordination

According to the bonding valence theory, the total bond order of C should not exceed 4, and that of H should not exceed 1, except in the case of hyper valence. However, even after correcting the original bond orders BOij, some degree of over-coordination may still remain in the molecule. To fix this, an over-coordination penalty term has been added to the force field.
For over-coordinated atoms (the degree to which the sum of the uncorrected bond orders around the atomic center deviates from its valence Vali, Δi > 0), Equation (3) imposes an energy penalty on the system. The form of Equation (3) ensures that Eover quickly vanishes to zero for under-coordinated systems (Δi < 0). Δi is calculated using Equation (4), using the corrected bond orders BOij instead of the uncorrected bond orders.
E o v e r = p o v e r · i · 1 1 + e x p λ 1 · i
i = j = 1 n b o n d B O i j V a l i

2.4. Valence Angle and Torsional Angle Interaction Terms

One of the disadvantages of the non-reactive force fields is their strict description of the angular and torsional interactions between atoms within the simulation. These types of interactions are often described as simple harmonic relationships, and the same harmonic potentials apply no matter how strong or weak the bond is. However, in ReaxFF, these angular and torsional interactions are also bond order dependent. This means that when an atom breaks a bond and leaves the molecule, the forces exerted on it due to the angle and torsion relative to the rest of the molecule diminish smoothly with the bond order. This is the expression in Equation (5):
E a n g l e = 1 exp λ · B O 1 3 1 exp λ · B O 2 3 · k a k b · exp k b · ϕ ϕ o 2
where BO1 and BO2 are the bond orders of each of the two bonds connecting the three atoms within the angle, λ is an angle parameter set to obtain agreement with quantum values, ka and kb are the harmonic force constants determining the depth and width of the angular potential, respectively, ϕ is the angle, and ϕo is the equilibrium angle.

2.5. Coulombic and van der Waals Interaction Terms

The ReaxFF is also able to calculate the polarization of intramolecular charges. This is achieved by using the electronegativity and hardness parameters of each element in the system. These values have also been optimized using QM data. Equation (6) shows how to calculate the polarization:
E q n = χ n + 2 · q n · η n + C · j = 1 n q j r n , j 3 + 1 / γ n , j 3 1 / 3 , i = 1 n q i = 0
where χn is the electronegativity of element n, ηn is the hardness of element n, and γn,j is the shielding parameter between atoms n and j [51,53,65,66,67]. This method is based on the electronegativity equalization method (EEM) [53,68,69,70]. These charge values are determined for each time step of the simulation and depend on the geometry of the system.
Due to the rigid connectivity associated with non-reactive force fields, Coulomb and van der Waals forces are usually only calculated between pairs of atoms that do not share bonds or valence angles with each other. However, in the ReaxFF, the Coulomb and van der Waals forces are calculated between all pairs of atoms regardless of their connectivity [51,53,65,66,67]. To avoid excessively repulsive or attractive non-bonded interactions at short distances, Coulomb and van der Waals interactions are shielded in ReaxFF. This is achieved by using a shielding term, γ, as shown in the Coulombic calculation of Equation (7):
E C o u l o m b = C · q i · q j r i j 3 + 1 / γ i j 3 1 / 3
where qi and qj are the charges of the two atoms, γij is the shielding parameter, and C is the coefficient. Atomic charges are calculated using the EEM approach [53,68,69,70]. The EEM charge derivation method is similar to the charge-equilibrium (QEq) scheme [70], with the only difference, aside from the parameter definition, that EEM does not use an iterative scheme for hydrogen charges (such as QEq), while QEq uses a stricter Slater orbital approach to account for charge overlap. However, the γij in Equation (7) can be optimized to reproduce the QEq orbital overlap correction [51].
In addition to overlapping-dependent valence interactions, there are also repulsive interactions at short interatomic distances due to orthogonalization of the Pauli principle and attractive energies at long distances due to dispersion. These interactions, consisting of van der Waals and Coulomb forces, are included in all atomic pairs, thus avoiding awkward changes in the energy description during bond dissociation. In this respect, ReaxFF is similar in nature to the central valence force fields used earlier in vibrational spectroscopy [37]. To account for the van der Waals interactions, a distance-corrected Morse potential is used (Equations (8) and (9)). Excessively high repulsive forces between bonded atoms (1–2 interactions) and valence-sharing atoms (1–3 interactions) can be avoided by including shielding interactions (Equation (9)) [51,67].
E v d W a a l s = T a p · D i j · exp α i j · 1 f 13 r i j r v d W 2 · exp 1 2 · α i j · 1 f 13 r i j r v d W
f 13 r i j = r i j p v d W 1 + 1 γ v d w p v d W 1 1 / p v d W 1
where
T a p = 20 R c u t 7 · r i j 7 70 R c u t 6 · r i j 6 + 84 R c u t 5 · r i j 5 35 R c u t 4 · r i j 4 + 1
where the terms in Equation (10) are used to avoid energy discontinuities outside the non-bonded cut-off radius (Rcut = 1.0 nm) of the ReaxFF. The terms in this polynomial are chosen to ensure that the first, second, and third derivatives of the non-bonded interactions with respect to the distance are all continuous and zero at the cut-off boundary. Therefore, this force field can adequately describe the weak van der Waals interactions even in the long range up to 1.0 nm.

2.6. Conjugated System Term

In theoretical chemistry, a conjugated system is a system in which p-orbitals of delocalized electrons are linked in a molecule, usually reducing the total energy of the molecule and increasing its stability. It is usually represented as having alternated single and multiple bonds. Lone pairs, radicals, or carbonium ions may be part of the system, which may be cyclic, linear, or mixed [71]. A conjugated system term is the overlap of one p-orbital with another p-orbital at an adjacent σ bond [72]. Molecules that contain conjugated systems of orbitals and electrons are called conjugated molecules, which have overlapping p-orbitals on three or more atoms. Some simple organic conjugated molecules are 1,3, butadiene, benzene, and allylic carbocations [73,74]. The largest conjugated systems are found in graphene, graphite, conductive polymers, and carbon nanotubes.
Equations (11) and (12) describe the contribution of conjugation effects to molecular energy. The conjugation energy contribution is greatest when the bond-order value of consecutive bonds is 1.5 (as in benzene and other aromatic compounds) [51].
E c o n j = f 12 B O i j , B O j k , B O k l · λ 26 · 1 + cos 2 ω i j k l 1 · sin Θ i j k · sin Θ j k l
f 12 B O i j , B O j k , B O k l = exp λ 27 · B O i j 1 1 2 2 · exp λ 27 · B O j k 1 1 2 2 · exp λ 27 · B O k l 1 1 2 2

2.7. Force Field Optimization Procedure

The force fields of reactive MD are optimized using a sequential single-parameter search technique as described by van Duin et al. [75]. In general, the search technique aims to reproduce the heat of formation to within 4.0 kcal/mol. Bond lengths are within 0.001 nm, and bond angles are within 2° of their literature values. In order to use the QC data during the force field optimization, structures related to these data are added to the force field training set. During the force field optimization, the energies of all molecules used to form the heat and geometric data comparisons are continuously minimized, while the structures related to the QC data are kept fixed or optimized with appropriate bond length or torsional angle constraints.

3. Results and Discussion

For all the above reasons, computational MD methods for biochar offer an alternative approach to understanding the mechanisms of biomass pyrolysis processes and carbonization reactions, predicting the outcomes and compositions of biochar after the processes, and characterizing the performance and functionalities of the biochar without costing a lot and taking a lot of time to produce, characterize, and experimentally study the insight into the final biochar and carbon materials. Furthermore, computational methods can be successfully employed to screen better biomass feedstocks or materials with suitable properties early in the development of various processes, which can provide opportunities for further research and evaluation and provide complementary information to experimentally collected data.

3.1. The Compositions and Physicochemical Properties of Biomass Feedstocks and the Produced Biochar

As mentioned above, biochar can be produced through thermochemical processes (e.g., torrefaction, pyrolysis, gasification, combustion, etc.). Biomass pyrolysis is generally defined as the thermal decomposition of biomass organic matrix in a non-oxidizing atmosphere to produce liquid bio-oil, solid biochar, and non-condensable gas products (also named syngas) [76]. Therefore, pyrolysis is demonstrated as one of the most promising technologies for converting biomass into biochar [77]. Biomass, like forest wastes and agricultural residues, is the most common feedstock for biochar production, which mainly consists of cellulose, hemicellulose, and lignin. The contents of cellulose, hemicellulose, and lignin vary significantly depending on the feedstock species, parts (leaf, root, stem, etc.) of plants, and growth conditions. The cellulose content may range from 40–60%, the hemicellulose content is between 15 and 30%, and the lignin content may vary from 10–25%. In addition to these three major components, a small fraction of extracts and inorganic ash are also present in biomass feedstocks as non-structural components that do not constitute the cell walls or cell layers [78]. Agricultural residues and forest wastes usually contain carbon (C) (40–55%), oxygen (O) (30–45%), hydrogen (H) (5–10%), nitrogen (N) (0.5–3%), sulfur (S) (<1%), and other trace elements [7,8,9], which may result in significant variations in the properties and compositions of produced biochar. For example, biochar is produced from cassava rhizomes, cassava stalks, and corncobs using a patented kiln by farmers. It was found that the biochar derived from cassava stems yielded the highest BET surface area among the biochar products, while the biochar produced from corncobs yielded the highest C (81.35%), and highest H content (2.42%). This was because the corncobs contain the highest C content (41.66%), the highest H content (6.84%), the lowest N content (0.74%), and the lowest O content (50.76%) [79]. In contrast, the cassava rhizomes contained the lowest C content (37.60%), the lowest H content (6.04%), the highest N content (1.27%), and the highest O content (55.3%). The study also showed that the biochar produced by slow pyrolysis with a longer reaction time at a lower heating rate and temperature was of high quality, stable C, and had a low H/C ratio. The high BET surface area and total pore volume of this biochar make it suitable for soil amendment, helping to reduce soil density, improve soil moisture and aeration, and reduce the leaching of plant nutrients from the rhizosphere.
Pyrolysis models usually pre-consider the reactions of each individual component of biomass without considering the potential interactions between the components. The behavior of biomass pyrolysis can be considered a comprehensive performance of the thermal cracking behavior of the three main components (cellulose, hemicellulose, and lignin). Table 1 lists the physicochemical properties, and Table 2 shows the major elemental compositions of different biochar samples from the example above [79,80,81,82,83,84]. All these data were used for the MD simulation model generation.

3.2. Carbonization Reactions in Biomass Pyrolysis Processes

MD simulations based on ReaxFF can be a powerful tool to analyze the biomass pyrolysis process. By adopting the bond-order formalism in a classical approach, ReaxFF implicitly describes the chemical bonds without costly QM calculations, giving insight into the biomass pyrolysis process and its carbonization. It can describe the bond breaking and formation during the chemical reactions, thereby exploring the complex carbonization reaction mechanisms at a nano/microscale, including different pyrolysis processes of various feedstocks as well as other chemical process mechanisms for producing biochar or other carbonaceous materials [29,30,52,55,59,64,76,77,78,85,86,87,88,89,90,91,92,93,94,95,96,97,98,99,100], such as thermochemical reactions in combustion and energy systems [50,53,57,101], energetic and dissociative water properties under various conditions [56], properties of carbon nano-rings, carbon nanotube bundles, and crosslinked epoxy resins [66,67,102], inclusion of geometry-dependent charge calculations [75], to name a few. Biomass can produce biochar through thermochemical processes such as pyrolysis, gasification, and combustion [77]. Chen et al. [57] constructed a simplified biomass model containing all three main components of cellulose, lignin, and hemicellulose and applied it to the simulations of pyrolysis and combustion processes under various oxidative and humid environments. The pyrolysis models usually presuppose the reactions of each individual component without considering the potential interactions between the major components.

3.2.1. Reactivity of Cellulose in Pyrolysis

Cellulose is the most abundant organic compound on earth. Cellulose, with the chemical formula (C6H10O5)n, is a polysaccharide composed of linear chains of hundreds to tens of thousands of 1,4-β-D-glucopyranose units. The schematic diagram of the long linear chain molecule structure of cellulose is shown in Figure 3.
By weight, cellulose accounts for 40–50% of biomass. It is difficult to understand the complex reaction process and detailed reaction mechanism of biomass pyrolysis only through experimental methods. Studying pyrolysis mechanisms and the product pathway of cellulose is of great significance for producing new carbon materials and exploring innovative applications. When using pyrolysis–gas chromatography–mass spectrometry (Py-GC/MS) to study the pyrolysis mechanism of cellulose, since Py-GC/MS cannot capture the reactions related to free radicals at the molecular level in a short time, it cannot provide detailed information on the pyrolysis mechanism [78,103]. ReaxFF MD simulations are effective methods to reveal the internal reaction mechanisms from a microscopic perspective. The decomposition of cellulose was classified into three categories using ReaxFF MD [104]: (1) Depolymerization reactions; (2) other chain scission reactions; and (3) release of low molecular weight products such as glycolaldehyde, water, formaldehyde, and formic acid. Since the simulated temperature is higher than the experimental temperature and the chemical and electrostatic ambient temperature of crystalline (or amorphous) cellulose does not exist, levoglucosan (LGA) is not observed in the simulation, and the kinetic parameters do not depend on molecular weight or initial conformation. Then, post-processing tools were used to parse the bond information. Specific algorithms were developed to search for LGA among decomposition products, taking into account both individual molecules and LGA end groups. Similar studies have been performed on glycolaldehyde, formaldehyde, formic acid, and hydroxymethyl radical. Some smaller products, such as water, carbon monoxide, carbon dioxide, and hydroxyl radicals, can be identified by their chemical compositions.
It should be pointed out that the experimental techniques, including product-specific Py-GC/MS experiments, cannot directly monitor the temporal evolution trend of specific products during cellulose pyrolysis because the radial reaction time is many orders of magnitude shorter than the experimental techniques allow, resulting in a lower concentration of pyrolyzed products. However, the simplicity of observing the evolution of different products as a function of time and temperature through ReaxFF MD simulations provides a feasible method for computationally probing the evolution of the pyrolysis products for experimental or industrial applications [97].
Zhang, et al. [97] combined large-scale models with GPU-based ReaxFF MD simulations using a canonical ensemble (conservation of substance quantity (N), volume (V), and temperature (T) of species, also known as the NVT ensemble) in a periodic cubic box and a unique cheminformatics-based reaction analysis tool (VARxMD) and studied the cellulose pyrolysis process and revealed the evolution and reaction mechanism of cellulose pyrolysis products. The three major products are hydroxyl-acetone, propyl aldehyde, and glycolaldehyde, in addition to levoglucosan, CO2, CO, H2O, etc., which are closely related to the experimental literature. Both the overall spectral product evolution and the underlying detailed chemical reactions of cellulose pyrolysis have been revealed. For example, the reaction pathway of hydroxyl-acetone is to undergo a series of depolymerization reactions such as homolysis, elimination, and ring opening at 800K to form the unstable compound C6H10O6. The unstable C6H10O6 then undergoes ring formation, rearrangement of unsaturated bonds, and hydroxyl radicals falling off to generate fragment C6H9O5. C6H9O5 decomposes into two species by bond dissociation, namely 2-hydroxyl-malonaldehyde (C3H4O3) and 2-hydroxyl-propyl-aldehyde radical (C3H5O2). 2-hydroxyl-propyl-aldehyde radial rearranges and reacts with the hydrogen radical to generate hydroxyl-aceton C3H6O2. In the reaction pathway of glycolaldehyde, C2H4O3, at 800 K, the generation of C2H4O3 comes from the bond dissociation of different compounds, such as C6H10O6 and C11H19O9 released by depolymerization of cellulose [97]. Since C6H10O6 is the initial reactant for the formation of hydroxyl-acetone and 2-hydroxy-propionaldehyde, the same initial fragment may undergo different reactions in the cellulose pyrolysis system to produce different compounds, and the three main products compete with each other for formation at high temperatures. The weight percentage of the pyrolysis products was found to be a function of temperature. When the temperature is low, the decomposition process predominates. As the temperature increases, the rate of thermal degradation of cellulose is accelerated, accompanied by the appearance of fragments and inorganic gas molecules. The inorganic gas production also increases with temperature between 500 and 1400 K, which is consistent with Py-GC/MS experiments. The pyrolysis products can be grouped roughly into syngas, tar (or bio-oil), and biochar [104]. Syngas and its gas compounds have carbon atoms equal to or less than 4, and bio-oil is a collection of compounds of C5–C9, C10–C19, C20–C29, and C30–C39. The rest (C40+) can be treated as biochar. The gases in syngas are released rapidly and increase monotonically with increasing temperature. When the temperature is increased, biochar is rapidly reduced, so a lower temperature favors a higher yield of biochar. The above results are in good agreement with the experimental data [105]. Analysis of the molecular products formed by cellulose pyrolysis in ReaxFF MD simulations at different time steps and temperatures ranging from 800K to 1400 K revealed that two compounds, namely C2H4O2 and C6H10O5, dominated the system. With the help of VARxMD, it was found that most of C2H4O2 is glycolaldehyde and C6H10O5 is levoglucosan or the precursor of levoglucosan, which is consistent with the experimental data of Py-GC/MS. At 700 K, the precursor of levoglucosan constitutes about 80% of the C6H10O5 compounds, which drops to only 15% at 1400 K. Since levoglucosan is a large molecule for the cellulose model, both levoglucosan and its precursors are considered levoglucosan products. The evolution trend of the large model simulation is closer to that of the Py-GC/MS experiments than the small model [97].

3.2.2. Reactivity of Lignin in Pyrolysis

Lignin is a polyphenolic polymer. Three types of phenylpropanoid units are generally considered the main precursors for lignin biosynthesis: coniferyl, sinapyl, and p-coumaryl alcohol (see Figure 4), which structurally give rise to guaiacyl (G), syringyl (S), and p-hydroxyphenyl (H) units, respectively, and are linked by different C–C bonds and ether bonds, such as α-O-4 bonds and β-O-4 bonds [106]. Lignin is the most recalcitrant of the three components of lignocellulosic biomass. Lignin, the second most abundant natural polymer after cellulose, is pyrolyzed differently from cellulose and hemicellulose due to the differences in chemical structure and characteristics [76]. Generally speaking, the lignin in softwoods is mainly composed of guaiacyl units and contains a small number of p-hydroxyphenyl units; in contrast, the lignin in hardwoods is mainly composed of guaiacyl units and syringyl units and contains a small number of p-hyfroxyphenyl units. The lignin in grasses typically contains all three types of monlignol units, with peripheral groups (i.e., hydroxycinnamic acids) incorporated into its core structure. The lignin macromolecules are mainly linked by C–C and C–O bonds between their phenylpropanoid structural units, among which aryl ether bonds (β-O-4) are the most common and important interunit linkages (see Figure 5).
It can be seen from Figure 5 that lignin contains a variety of oxygen-containing functional groups, including methoxyl, hydroxyl, carboxyl, and carbonyl groups, etc., which significantly affect the reactivity of lignin. The content of methoxyl groups in lignin is related to the formation of lignin pyrolysis char, i.e., lignin with high methoxyl group contents produces less char during the pyrolysis process.
Many lignin models have been proposed for ReaxFF MD simulations in the last century, including the softwood lignin models by Alder, Freudenberg, Brukraft, Forss, and Glasser; the hardwood lignin model by Nimz; and the pine kraft lignin model by Marton [76,107]. The corresponding details of the four lignin models are listed in Table 3 [76,107].
The simulation results show that the Alder, Frendenberg, and Nimz lignin models have three pyrolysis stages. The first stage refers to the main process, starting with the formation of initial pyrolysis products and ending with the complete consumption of source lignin molecules. The second stage starts with the cracking of primary pyrolysis products into secondary pyrolysis products and ends with a maximum. The third stage begins with the reduction of secondary pyrolysis products. However, the Marton lignin model only has the first and second stages. Consumptions of all α-O-4 linkages were similar in the softwood and hardwood lignin models. The differences in pyrolysis product evolution and linkage behavior between hardwood, softwood, and kraft lignin can be attributed to the different reactions of linkages and their linked monomers induced by different oxygen-containing substituents [76]. The ring structure evolutions of the four lignin models are almost identical. The similarities and differences in the pyrolysis of the different lignin models suggested that the simulation work can provide new insights into the high-value utilization of lignin.
The study of the ReaxFF MD simulated pyrolysis of lignin model will be helpful for better understanding the reaction behavior and pyrolysis characteristics in the pyrolysis of lignin, so as to optimize the pyrolysis process of lignin, and even the whole biomass feedstock. Furthermore, lignin is heterogeneous and does not have a well-defined primary structure. The phenylpropane units are organized into a 3D amorphous polymeric network with varying degrees of aggregation. Abundant forms of biomass are a promising alternative to fossil fuels [108]. Lignin depolymerization is a difficult process that requires specific reaction conditions [109,110] and selected catalysts [108] to increase the yield of high aromatic monomers and produce value-added fuels and chemicals.

3.2.3. Reactivity of Hemicellulose in Pyrolysis

Hemicellulose has a heteropolymeric structure (lower molecular weight than cellulose) that is composed of a variety of sugar monomers, including glucose, galactose, mannose, xylose, arabinose, 4-O-methyl glucuronic acid, and galacturounic acid residues. The composition of hemicellulose from different biomass materials varies, and its structure is more complex than that of cellulose. Xylan is often used as a typical hemicellulose model [111,112], as shown in Figure 6. Previous work has provided important insights into the relationship between the distribution of pyrolysis products and the structural features of xylosyl hemicelluloses [113,114,115,116]. The main differences between xylan and cellulose pyrolysis are as follows: (1) xylan melts and generates bubbles during pyrolysis, producing many unidentified didehydrated pentose, while levoglucosan is the main compound in the dellulose bio-oil; (2) xylan and pyranose molecules formed during the pyrolysis of xylan tend to form char through multi-step dehydration reactions [113,114,115,116,117,118,119,120].
Furthermore, ReaxFF MD simulations showed that cellulose and hemicellulose were the main sources of CO and CO2 production, although there were slightly more CO2 molecules in hemicellulose due to the potential presence of carboxyl and carbonyl groups. Due to the high O content, cellulose and hemicellulose are the main sources of C2H2O2 and CH3CHO. Most of the C2H4O2 molecules come from the degeneration of the cellulose. In addition, CH3OH is mainly derived from cellulose due to the cleavage of the hydroxyl-rich pyran rings. The dissociation of -OCH3 radicals in the lignin also provides the precursors for the formation of CH3OH [59]. Table 4 lists the percentage of each gas produced during the pyrolysis process of the three main components of the biomass (wheat straw) at a temperature of 2000 K [59].
Functional groups play an important role in the pyrolysis process and determine the chemical properties of organic compounds. Cellulose, hemicellulose, and lignin in biomass have different kinds and numbers of functional groups, which lead to their different properties. In addition, phenyl rings in lignin and pyran rings in cellulose and hemicellulose also have potential effects on the thermal stability of organics. Therefore, the study of its evolution and behavior is beneficial to understanding the mechanism of biomass pyrolysis [59]. The functional groups are generally divided into six categories: ether groups (R–O–R), hydroxyl groups (R–OH), aldehyde groups (R–CHO), ester groups (R–COO–R), carboxyl groups (R–COOH), and carbonyl groups (R–CO–R). The ether and hydroxyl groups are the most abundant, while the other four functional groups are much less abundant. The hydroxyl groups are concentrated in the cellulose and hemicellulose in the form of alcoholic hydroxyl groups [121]. The ether groups are mainly derived from the pyran rings and glycosidic bonds in the cellulose and hemicellulose, as well as the methoxy or other carbon structures connected to the phenyl rings in the lignin [122]. The aldehyde and carbonyl groups are concentrated in lignin, while hemicellulose provides carboxyl groups. Furthermore, the ester groups are distributed in hemicellulose and lignin.

3.3. Additional Notes

While computational methods are inexpensive, they can sometimes provide insight information where experimental methods are impossible. However, computational methods have their own limitations. In general, any model is limited by the size of the structure being modeled, the period of reaction time it can simulate, and ultimately the processing power available. QM-based ab initio and DFT can provide more accurate representations of materials by calculating the potentials acting on all atoms from first principles, but they are computationally resource intensive, limiting their usefulness for predicting the behavior of bulk materials.
Classical MD approximates the interatomic potentials using force fields that have been parameterized using ab initio calculations and empirical data to approximate the behavior of the system as accurately as possible, with well-defined bonds, angles, dihedrals, and impropers explicitly.
ReaxFF is a reactive implementation of MD where bonding is implicit and the bond orders between atomic pairs are functions of interatomic distances. By setting the parameters of hydrocarbon reactions or carbonization processes, reactive MD has been applied to the combustion, oxidation, and pyrolysis of coal [57,90,91,92,93], the pyrolysis reaction processes and mechanisms of polyethylene [94], and polycarbonate [95,96], to name a few. Among other things, modeling the thermal decomposition of cured epoxy resins and their mechanical responses elucidated unique insights into decomposition steps and failure modes, which were difficult to determine experimentally [98]. Furthermore, studies of the pyrolysis of polymeric ablative materials for thermal protection systems, a popular choice of ablator using ReaxFF, are another application area [102]. The potential applications of reactive MD simulations can be extended to processes involving chemical reactions, such as the design and evaluation of ablative materials, the design and development of energetic materials, the design of hydrogen storage materials, the study of petrochemical manufacturing methods, etc.
Many key factors, such as the biomass feedstock species, reactor types and operating conditions, post-treatment of biochar activation, physiochemistry, morphology, performance, and functionalities of the biomass, synergistically affect the pyrolysis and carbonization processes and the production of inorganic gas, organic gas, biochar, bio-oil, etc. [29,30,44,57,111,113]. The reactive processes and mechanisms that can be simulated and analyzed quantitatively and accurately remain challenging.

4. Conclusions

Each biochar produced is unique due to the biomass feedstock species, the reactor types and operating conditions, and the synergistic effects of different factors. The biochar production processes are complex and involve chemical reactions that determine the physicochemical properties of the final biochar produced. Classical MD simulations cannot effectively analyze the mechanisms of biomass pyrolysis processes and carbonization reactions This paper has reviewed the fundamentals of the ReaxFF method and the corresponding MD simulation findings based on ReaxFF to characterize the biomass pyrolysis processes and the physicochemical properties of the produced biochar, as well as the carbonization reactions in pyrolysis. The possible composition of the biomass and the physicochemical properties of the produced biochar have been summarized. The models and reactivities of the three major components of biomass: Cellulose, lignin, and hemicellulose in pyrolysis have been discussed. The potential applications of ReaxFF MD were briefly discussed. The MD simulations based on ReaxFF can be an effective method for understanding various processes involving chemical reactions, such as the design and evaluation of ablative materials, the design and development of energetic materials, the design of hydrogen storage materials, the investigation of petrochemical manufacturing methods, in addition to the carbonization mechanisms of the biomass pyrolysis process. The information derived from the simulations will be helpful for further optimization of the processes.

Author Contributions

Conceptualization, Z.H. and L.W.; methodology, Z.H. and L.W.; formal analysis, Z.H.; investigation, Z.H.; resources, Z.H. and L.W.; data curation, Z.H.; writing—original draft preparation, Z.H.; writing—review and editing, Z.H. and L.W.; visualization, Z.H.; supervision, Z.H.; project administration, Z.H.; funding acquisition, Z.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Not applicable.

Acknowledgments

This work was supported by the Department of Mechanical Engineering in the J. J. Lohr College of Engineering and the Department of Agricultural & Biosystems Engineering in the College of Agriculture, Food & Environmental Sciences at South Dakota State University and are gratefully acknowledged.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Wei, F.; Cao, J.-P.; Zhao, X.-Y.; Ren, J.; Wang, J.-X.; Fan, X.; Wei, X.-Y. Nitrogen evolution during fast pyrolysis of sewage sludge under inert and reductive atmospheres. Energy Fuels 2017, 31, 7191–7196. [Google Scholar] [CrossRef]
  2. Sun, H.; Wu, C. Autothermal CaO looping biomass gasification for renewable syngas production. Environ. Sci. Technol. 2019, 53, 9298–9305. [Google Scholar] [CrossRef] [PubMed]
  3. Ren, J.; Liu, Y.-L.; Zhao, X.-Y.; Cao, J.-P. Biomass thermochemical conversion: A review on tar elimination from biomass catalytic gasification. J. Energy Inst. 2020, 93, 1083–1098. [Google Scholar] [CrossRef]
  4. Li, H.; Xu, B.; Jin, H.; Luo, K.; Fan, J. Molecular dynamics investigation on the lignin gasification in supercritical water. Fuel Process. Technol. 2019, 192, 203–209. [Google Scholar] [CrossRef]
  5. Marangwanda, G.T.; Madyira, D.M.; Babarinde, T.O. Combustion models for biomass: A review. Energy Rep. 2020, 6, 664–672. [Google Scholar] [CrossRef]
  6. Kan, T.; Strezov, V.; Evans, T.J. Lignocellulosic biomass pyrolysis: A review of product properties and effects of pyrolysis parameters. Renew. Sustain. Energy Rev. 2016, 57, 1126–1140. [Google Scholar] [CrossRef]
  7. Garcia, B.; Alves, O.; Rijo, B.; Lourinho, G.; Nobre, C. Biochar: Production, applications, and market prospects in Portugal. Environments 2022, 9, 95. [Google Scholar] [CrossRef]
  8. Gupta, R.; Dubey, M.; Kharel, P.; Gu, Z.; Fan, Q. Biochar activated by oxygen plasma for supercapacitors. J. Power Sources 2015, 274, 1300–1305. [Google Scholar] [CrossRef]
  9. Clough, T.J.; Condron, L.M.; Kammann, C.; Müller, C. A review of biochar and soil nitrogen dynamics. Agronomy 2013, 3, 275–293. [Google Scholar] [CrossRef]
  10. Steiner, C.; McLaughlin, H.; Harley, A.; Flora, G.; Larson, R.; Reed, A. United States Biochar Initiative 2010: U.S.-Focused Biochar Report–Assessment of Biochar’s Benefits for the United States of America; United States Center for Energy and Environmental Security; United States Biochar Initiative: Morgantown, WV, USA, 2010; p. 84. [Google Scholar]
  11. Kumar, A.; Singh, J.; Churasia, A. Production of biodiesel from soybean oil biomass as renewable energy source. J. Environ. Biol. 2016, 37, 1303–1307. [Google Scholar]
  12. Paletto, A.; Bernardi, S.; Pieratti, E.; Teston, F.; Romagnoli, M. Assessment of environmental impact of biomass power plants to increase the social acceptance of renewable energy technologies. Heliyon 2019, 5, e02070. [Google Scholar] [CrossRef]
  13. Hodges, D.G.; Chapagain, B.; Watcharaanantapong, P.; Poudyal, N.C.; Kline, K.L.; Dale, V.H. Opportunities and attitudes of private forest landowners in supplying woody biomass for renewable energy. Renew. Sustain. Energy Rev. 2019, 113, 109205. [Google Scholar] [CrossRef]
  14. Omoriyekomwan, J.E.; Tahmasebi, A.; Zhang, J.; Yu, J. Synthesis of Super-Long Carbon Nanotubes from Cellulosic Biomass under Microwave Radiation. Nanomaterials 2022, 12, 737. [Google Scholar] [CrossRef]
  15. Hersh, B.; Mirkouei, A.; Sessions, J.; Rezaie, B.; You, Y. A review and future directions on enhancing sustainability benefits across food-energy-water systems: The potential role of biochar-derived products. AIMS Environ. Sci. 2019, 6, 379–416. [Google Scholar] [CrossRef]
  16. Bartoli, M.; Giorcelli, M.; Jagdale, P.; Rovere, M.; Tagliaferro, A. A review of non-soil biochar applications. Materials 2020, 13, 261. [Google Scholar] [CrossRef] [PubMed]
  17. Allohverdi, T.; Mohanty, A.K.; Roy, P.; Misra, M. A review on current status of biochar uses in agriculture. Molecules 2021, 26, 5584. [Google Scholar] [CrossRef] [PubMed]
  18. Iannazzo, D.; Celesti, C.; Espro, C.; Ferlazzo, A.; Giofrè, S.V.; Scuderi, M.; Scalese, S.; Gabriele, B.; Mancuso, R.; Ziccarelli, I.; et al. Orange-peel-derived nanobiochar for targeted cancer therapy. Pharmaceutics 2022, 14, 2249. [Google Scholar] [CrossRef] [PubMed]
  19. Fernandez, M.E.; Ledesma, B.; Román, S.; Bonelli, P.R.; Cukierman, A.L. Development and characterization of activated hydrochars from orange peels as potential adsorbents for emerging organic contaminants. Bioresour. Technol. 2015, 183, 221–228. [Google Scholar] [CrossRef]
  20. Wang, T.; Liu, J.; Zhang, Y.; Zhang, H.; Chen, W.-Y.; Norris, P.; Pan, W.-P. Use of a non-thermal plasma technique to increase the number of chlorine active sites on biochar for improved mercury removal. Chem. Eng. J. 2018, 331, 536–544. [Google Scholar] [CrossRef]
  21. Wang, J.; Wang, S. Preparation, Modification and Environmental Application of Biochar: A Review. J. Clean. Prod. 2019, 227, 1002–1022. [Google Scholar] [CrossRef]
  22. Adhamash, E.; Pathak, R.; Qiao, Q.; Zhou, Y.; McTaggart, R. Gamma-radiated biochar carbon for improved supercapacitor performance. RSC Adv. 2020, 10, 29910–29917. [Google Scholar] [CrossRef]
  23. Alam, M.S.; Alessi, D.S. Chapter 4—Modeling the Surface Chemistry of Biochars. In Biochar from Biomass and Waste—Fundamentals and Applications; Ok, Y.S., Tsang, D.C.W., Bolan, N., Novak, J.M., Eds.; Elsevier: Amsterdam, The Netherlands, 2019; pp. 59–72. [Google Scholar]
  24. Mukherjee, A.; Zimmerman, A.R.; Harris, W. Surface chemistry variations among a series of laboratory-produced biochars. Geodema 2011, 163, 247–255. [Google Scholar] [CrossRef]
  25. Das, C.; Tamrakar, S.; Kiziltas, A.; Xie, X. Incorporation of biochar to improve mechancal, thermal and electrical properties of polymer composites. Polymers 2021, 13, 2663. [Google Scholar] [CrossRef] [PubMed]
  26. Zhang, P.; Duan, W.; Peng, H.; Pan, B.; Xing, B. Functional biochar and its balanced design. ACS Anvironmental 2022, 2, 115–127. [Google Scholar] [CrossRef] [PubMed]
  27. Chen, C.-K.; Nguyen, N.-T.; Le, T.-T.; Duong, C.-C.; Duong, T.-T. Specifically designed amine functional group doped sludge biochar for inorganic and organic arsenic removal. Sustain. Environ. Res. 2021, 31, 28. [Google Scholar] [CrossRef]
  28. Bamdad, H.; Papari, S.; MacQuarrie, S.; Hawboldt, K. Study of surface heterogeneity and nitrogen functinalization of biochars: Molecular modeling approach. Carbon 2021, 171, 161–170. [Google Scholar] [CrossRef]
  29. Hong, D.; Gao, P.; Wang, C. A comprehensive understanding of the synergistic effect during co-pyrolysis of polyvinyl chloride (PVC) and coal. Energy 2022, 239, 122258. [Google Scholar] [CrossRef]
  30. Wang, M.; Gao, J.; Xu, J.; Du, Q.; Zhang, Y. Effect of H2O on the transformation of sulfur during demineralized coal pyrolysis: Molecular dynamics simulation using ReaxFF. Energy Fuels 2021, 35, 2379–2390. [Google Scholar] [CrossRef]
  31. Cui, L.; Fan, Q.; Sun, J.; Quan, G.; Yan, J.; Hina, K.; Wang, H.; Zhang, Z.; Hussain, Q. Change in surface characteristics and adsorption properties of 2, 4, 6-trichlorophenol following Fenton-like aging of biochar. Sci. Rep. 2021, 11, 4293. [Google Scholar] [CrossRef]
  32. Muretta, J.E. Investigation of the Structure, Chemistry, and Functional Performance of Novel Biochars. Ph.D. Dissertation, Montana Technological University, Butte, MT, USA, 2022. [Google Scholar]
  33. Feng, D.; Guo, D.; Zhang, Y.; Sun, S.; Zhao, Y.; Shang, Q.; Sun, H.; Wu, J.; Tan, H. Functionalized construction of biochar with hierarchical pore strucyures and surface O-/N-containing groups for phenol adsorption. Chem. Eng. J. 2021, 410, 127707. [Google Scholar] [CrossRef]
  34. Xiao, X.; Chen, B. A direct observation of the fine aromatic clusters and molecular structures of biochars. Environ. Sci. Technol. 2017, 61, 5473–5482. [Google Scholar] [CrossRef]
  35. Jimenez-Cordero, D.; Heras, F.; Alonso-Morales, N.; Gilarranz, M.A.; Rodriguez, J.J. Porous structure and morphology of granular chars from flash and conventional pyrolysis of grape seeds. Biomass Bioenergy 2013, 54, 123–132. [Google Scholar] [CrossRef]
  36. Brewer, C.E.; Chuang, V.J.; Masiella, C.A.; Gonnermann, H.; Gao, X.; Dugan, B.; Driver, L.E.; Panzacchi, P.; Zygourakis, K.; Davies, C.A. New approaches to measuring biochar density and porosity. Biomass Bioenergy 2014, 66, 176–185. [Google Scholar] [CrossRef]
  37. Vidal-Vidal, Á.; Faza, O.N.; López, C.S. CO2 complexes with five-membered heterocycles: Structure, topology, and spectroscopic characterization. J. Phys. Chem. A 2017, 121, 9118–9130. [Google Scholar] [CrossRef] [PubMed]
  38. Yaashikaa, P.R.; Kumar, P.S.; Varijani, S.; Saravanan, A. A critical review on the biochar production techniques, characterizaiton, stability and applications for circular bioeconomy. Biotechnol. Rep. 2020, 28, e00570. [Google Scholar] [CrossRef]
  39. Keiluwrit, M.; Nico, P.S.; Johnson, M.G.; Kleber, M. Dynamic molecular structure of plant biomass-derived black carbon (biochar). Environ. Sci. Technol. 2010, 44, 1247–1253. [Google Scholar] [CrossRef]
  40. Kaal, J.; Cortizas, A.M.; Reyes, O.; Soliño, M. Molecular characterization of Ulex europaeus biochar obtained from laboratory heat treatment experiments—A pyrolysis—GC/MS study. J. Anal. Appl. Pyrolysis 2012, 95, 205–212. [Google Scholar] [CrossRef]
  41. Conte, P.; Bertani, R.; Sgarbossa, P.; Bambina, P.; Schmidt, H.-P.; Raga, R.; Papa, G.L.; Martino, D.F.C.; Meo, P.L. Recent developments in understanding biochar’s pysical-chemistry. Agronomy 2021, 11, 615. [Google Scholar] [CrossRef]
  42. Zuhara, S.; Mackey, H.R.; Al-Ansari, T.; Mckay, G. A review of propsects and current scenarios of biomass co-pyrolysis for water treatment. Biomass Convers. Biorefinery 2022. [Google Scholar] [CrossRef]
  43. Mehdi, R.; Khoja, A.H.; Naqvi, S.R.; Gao, N.; Amin, N.A.S. A review on production and surface modifications of biochar materials via biomass pyrolysis process for supercapacitor applications. Catalysts 2022, 12, 798. [Google Scholar] [CrossRef]
  44. Najmi, L.; Hu, Z. Review on molecular dynamics simulations of effects of carbon nanotubes (CNTs) on electrical and thermal conductivities of CNT-modified polymeric composites. J. Compos. Sci. 2023, 7, 165. [Google Scholar] [CrossRef]
  45. Zhao, N.; Lv, Y.; Yang, X. A new 3D conceptual structure modeling of biochars by molecular mechanic and molecular dynamic simulation. J. Solids Sediments 2017, 17, 641–655. [Google Scholar] [CrossRef]
  46. Dash, B. Carbon dioxide capture by nitrogen containing organic materials—A density functional theory investigation. Comput. Theor. Chem. 2018, 1128, 1–14. [Google Scholar] [CrossRef]
  47. Mrozik, W.; Minofar, B.; Thongsamer, T.; Wiriyaphong, N.; Khawkomol, S.; Plaimart, J.; Vakros, J.; Karapanagioti, H.; Vinitnantharat, S.; Werner, D. Valorisation of agricultural waste derived biochars in aquaculture to remove organic micropollutants from water—Experimental study and molecular dynamics simulations. J. Environ. Manag. 2021, 300, 113717. [Google Scholar] [CrossRef] [PubMed]
  48. Ouachtak, H.; Guerdaoui, A.E.; Haouti, R.E.; Haounati, R.; Ighnih, H.; Toubi, Y.; Alakhras, F.; Rehman, R.; Hafid, N.; Addi, A.A.; et al. Combined molecular dynamics simulations and experimental studies of the removal of cationic dyes on the eco-friendly adsorbent of activated carbon decorated montmorillonite Mt@AC. RSC Adv. 2023, 13, 5027–5044. [Google Scholar] [CrossRef]
  49. Gooneie, A.; Schuschnigg, S.; Holzer, C. A review of multiscale computational methods in polymeric materials. Polymers 2017, 9, 16. [Google Scholar] [CrossRef]
  50. Mao, Q.; Feng, M.; Jiang, X.Z.; Ren, Y.; Luo, K.H.; van Duin, A.C.T. Classical and reactive molecular dyanmics: Principles and applications in combustion and energy systems. Prog. Energy Combust. Sci. 2023, 97, 101084. [Google Scholar] [CrossRef]
  51. van Duin, A.C.T.; Dasgupta, S.; Lorant, F.; Goddard III, W.A. ReaxFF: A reactive force filed for hydrocarbons. J. Physcal Chem. A 2001, 105, 9396–9409. [Google Scholar] [CrossRef]
  52. Chenoweth, K.; van Duin, A.C.T.; Goddard III, W.A. ReaxFF force field for molecular dynamics simulations of hydrocarbon oxidation. J. Phys. Chem. A 2008, 112, 1040–1053. [Google Scholar] [CrossRef]
  53. Russo, M.F., Jr.; van Duin, A.C.T. Atomistic-scale simulations of chemical reactions: Bridging from quantum chemistry to engineering. Nucl. Instrum. Methods Phys. Res. B 2011, 269, 1549–1554. [Google Scholar] [CrossRef]
  54. Aktulga, H.M.; Pandit, S.A.; van Duin, A.C.T.; Grama, A.Y. Reactive molecular dynamics: Numerical methods and algorithmic techniques. SIAM J. Sci. Comput. 2012, 34, C1–C23. [Google Scholar] [CrossRef]
  55. Beste, A. ReaxFF study of the oxidation of softwood lignin in view of carbon fiber production. Energy Fuels 2014, 28, 7007–7013. [Google Scholar] [CrossRef]
  56. Senftle, T.P.; Hong, S.; Islam, M.M.; Kylasa, S.B.; Zheng, Y.; Shin, Y.K.; Junkermeier, C.; Engel-Herbert, R.; Janik, M.J.; Aktulga, H.M.; et al. The ReaxFF reactive force-field: Development, applications and future directions. NPJ Comput. Mater. 2016, 2, 15011. [Google Scholar] [CrossRef]
  57. Chen, C.; Zhao, L.; Wang, J.; Lin, S. Reactive molecular dynamics simulations of biomass pyrolysis and combustion under various oxidative and humidity environments. Ind. Eng. Chem. Res. 2017, 56, 12278–12288. [Google Scholar] [CrossRef]
  58. Leven, I.; Hao, H.; Das, A.K.; Head-Gordon, T. A reactive force field with coarse-grained electrons for liquid water. J. Phys. Chem. Lett. 2020, 11, 9240–9247. [Google Scholar] [CrossRef]
  59. Liu, Z.; Ku, X.; Jin, H. Pyrolysis mechanism of wheat straw based on ReaxFF molecular dynamics simulations. ACS Omega 2022, 7, 21075–21085. [Google Scholar] [CrossRef]
  60. Friesner, R.A. Ab initio quantum chemistry: Methodology and applications. Proc. Natl. Acad. Sci. USA 2005, 102, 6648–6653. [Google Scholar] [CrossRef]
  61. Lee, H.M.; Youn, I.S.; Saleh, M.; Lee, J.W.; Kim, K.S. Interactions of CO2 with various functional molecules. Phys. Chem. Chem. Phys. 2015, 17, 10925. [Google Scholar] [CrossRef]
  62. Zhao, N.; Zhao, C.; Liu, K.; Zhang, W.; Tsang, D.C.W.; Yang, Z.; Yang, X.; Yan, B.; Morel, J.L.; Qiu, R. Experimental and DFT investigation on N-functionalized biochars for enhanced removel of Cr(VI). Environ. Pollut. 2021, 291, 118244. [Google Scholar] [CrossRef]
  63. Car, R.; Parrinello, M. Unified approach for molecular dynamics and density-functional theory. Phys. Rev. Lett. 1985, 55, 2471–2474. [Google Scholar] [CrossRef]
  64. Purse, M.; Holmes, B.; Sacchi, M.; Howlin, B. Simulation the complete pyrolysis and charring process of phenol-formaldehyde resins using reactive molecular dynamics. J. Mater. Sci. 2022, 57, 7600–7620. [Google Scholar] [CrossRef]
  65. van Duin, A.C.T.; Damsté, J.S.S. Computational chemical investigation into isorenieratene cyclisation. Org. Geochem. 2003, 34, 515–526. [Google Scholar] [CrossRef]
  66. Chen, N.; Lusk, M.T.; van Duin, A.C.T.; Goddard III, W.A. Mechanical properties of connected carbon nanorings via molecular dynamics simulation. Phys. Rev. B 2005, 72, 085416. [Google Scholar] [CrossRef]
  67. Han, S.S.; Kang, J.K.; Lee, H.M.; van Duin, A.C.T.; Goddard III, W.A. Liquefaction of H2 molecules upon exterior surfaces of carbon nanotube bundles. Appl. Phys. Lett. 2005, 86, 203108. [Google Scholar] [CrossRef]
  68. Mortier, W.J.; Ghosh, S.K.; Shankar, S. Electronegativity-equalization method for the calculation of atomic charges in molecules. J. Am. Chem. Soc. 1986, 108, 4315–4320. [Google Scholar] [CrossRef]
  69. Janssens, G.O.A.; Baekelandt, B.G.; Toufar, H.; Mortier, W.J.; Schoonheydt, R.A. Comparison of cluster and infinite crystal calculations on zeolites with the electronegativity equalization method (EEM). J. Phys. Chem. 1995, 99, 3251–3258. [Google Scholar] [CrossRef]
  70. Rappe, A.K.; Goddard, W.A., III. Charge equilibration for molecular dynamics simulations. J. Phys. Chem. 1991, 95, 3358–3363. [Google Scholar] [CrossRef]
  71. Johannes, T. Zur Kenntnis der ungesättigten verbindungen. Justus Liebig’s Ann. Der Chem. 1899, 306, 87–142. [Google Scholar]
  72. UPAC. Compendium of Chemical Terminology, 2nd ed.; McNaught, A.D., Wilkinson, A., Eds.; The “Gold Book”; Blackwell Scientific Publications: Oxford, UK, 1997; ISBN 0-9678550-9-8. [Google Scholar] [CrossRef]
  73. March, J. Advanced Organic Chemistry Reactions, Mechanisms and Structure, 3rd ed.; John Wiley & Sons, Inc.: New York, NY, USA, 1985; ISBN 0-471-85472-7. [Google Scholar]
  74. Smith, J.G. Chapter 16: Conjugation, Resonance, and Dienes. In Organic Chemistry, 5th ed.; Iswar Chandra Vidyasagar College: Belonia, Tripura, India, 2022; pp. 604–640. ISBN 9780078021558. [Google Scholar]
  75. van Duin, A.C.T.; Baas, J.M.A.; van de Graaf, B.J. Delft molecular mechanics: A new approach to hydrocarbon force fields. Inclusion of a geometry-dependent charge calculation. J. Chem. Soc. Faraday Trans. 1994, 90, 2881–2895. [Google Scholar] [CrossRef]
  76. Li, G.; Zheng, F.; Huang, Q.; Wang, J.; Niu, B.; Zhang, Y.; Long, D. Molecular insight into pyrolysis processes via reactive force field molecular dynamics: A state-of-the-art review. J. Anal. Appl. Pyrolysis 2022, 166, 105620. [Google Scholar] [CrossRef]
  77. Bridgwater, A.V. Review of fast pyrolysis of biomass and product upgrading. Biomass Bioenergy 2012, 38, 68–94. [Google Scholar] [CrossRef]
  78. Wang, S.; Dai, G.; Yang, H.; Luo, Z. Lignocellulosic biomass pyrolysis mechanism: A state-of-the-art review. Prog. Energy Combust. Sci. 2017, 62, 33–86. [Google Scholar] [CrossRef]
  79. Wijitkosum, S.; Jiwnok, P. Elemental composition of biochar obtained from agricultural waste for soil amendment and carbon sequestration. Appl. Sci. 2019, 9, 3980. [Google Scholar] [CrossRef]
  80. Ippolito, J.A.; Cui, L.; Kammann, C.; Wrage-Mönnig, N.; Estavillo, J.M.; Fuertes-Mendizabal, T.; Cayuela, M.L.; Sigua, G.; Novak, J.; Spokas, K.; et al. Feedstock choice, pyrolysis temperature and type influence biochar characteristics: A comprehensive meta-data analysis review. Biochar 2020, 2, 421–438. [Google Scholar] [CrossRef]
  81. Libra, J.A.; Ro, K.S.; Kammann, C.; Funke, A.; Berge, N.D.; Neubauer, Y.; Titirici, M.-M.; Uhner, C.F.; Bens, O.; Kern, J.; et al. Hydrothermal carbonizationof biomass residuals: A comparative review of the chemistry, processes and applications of wet and dry pyrolysis. Biofuels 2011, 2, 71–106. [Google Scholar] [CrossRef]
  82. Brownsort, P.A. Biomass Pyrolysis Processes: Performance Parameters and Their Influence on Biochar System benefits. Master’s Dissertation, The University of Edinburgh, Edinburgh, Scotland, UK, 2009. [Google Scholar]
  83. Antal, M.J.; Grønli, M. The art, science, and technology of charcoal production. Ind. Eng. Chem. Res. 2003, 42, 1619–1640. [Google Scholar] [CrossRef]
  84. Lucian, M.; Fiori, L. Hydrothermal carbonization of waste biomass: Process design, modeling, energy efficiency and cost analysis. Energies 2017, 10, 211. [Google Scholar] [CrossRef]
  85. Arvelos, S.; Abrahão, O., Jr.; Hori, C.E. ReaxFF molucular dynamics study on the pyrolysis process of cyclohexanone. J. Analtical Appl. Pyrolysis 2019, 141, 104620. [Google Scholar] [CrossRef]
  86. Xu, F.; Liu, H.; Wang, Q.; Pan, S.; Zhao, D.; Liu, Q.; Liu, Y. ReaxFF-based molecular dyanamics simulation of the intial pyrolysis mechanism of lignite. Fuel Process. Technol. 2019, 195, 106147. [Google Scholar] [CrossRef]
  87. Xu, F.; Liu, H.; Wang, Q.; Pan, S.; Zhao, D.; Liu, Y. Study of non-isothermal pyrolysis mechanism of lignite using ReaxFF molecular dynamics simulations. Fuel 2019, 256, 115884. [Google Scholar] [CrossRef]
  88. Lele, A.; Kwon, H.; Ganeshan, K.; Xuan, Y.; van Duin, A.C.T. ReaxFF molecular dynamics study on pyrolysis of bicyclic compunds for aviation fuel. Fuel 2021, 297, 120724. [Google Scholar] [CrossRef]
  89. Guo, J.; Zhang, Y. Reactive molecular dynamics simulation on degradation of tetracycline antibiotics treated by cold atmospheric plasmas. Molecules 2023, 28, 3850. [Google Scholar] [CrossRef] [PubMed]
  90. Gao, M.; Li, X.; Guo, L. Pyrolysis simulations of Fugu coal by large-scale ReaxFF molecular dynamics. Fuel Process. Technol. 2018, 178, 197–205. [Google Scholar] [CrossRef]
  91. Zheng, M.; Li, X.; Guo, L. Investigation of N behavior during coal pyrolysis and oxidation using ReaxFF molecular dynamics. Fuel 2018, 233, 867–876. [Google Scholar] [CrossRef]
  92. Bhoi, S.; Banerjee, T.; Mohanty, K. Molecular dynamic simulation of spontaneous combustion and pyrolysis of brown coal using ReaxFF. Fuel 2014, 136, 326–333. [Google Scholar] [CrossRef]
  93. Liu, X.; Li, X.; Liu, J.; Wang, Z.; Kong, B.; Gong, X.; Yang, X.; Lin, W.; Guo, L. Study of high density polyethylene (HDPE) pyrolysis with reactive molecular dynamics. Polym. Degrad. Stab. 2014, 104, 62–70. [Google Scholar] [CrossRef]
  94. Zhao, T.; Li, T.; Xin, Z.; Zou, L.; Zhang, L. A ReaxFF-based molecular dynamics simulation of the pyrolysis mechanism for polycarbonate. Energy Fuels 2018, 32, 2156–2162. [Google Scholar] [CrossRef]
  95. Liu, Q.; Liu, S.; Lv, Y.; Hu, P.; Huang, Y.; Kong, M.; Li, G. Atomic-scale insight into the pyrolysis of polycarbonate by ReaxFF-based reactive molecular dynamics simulation. Fuel 2021, 287, 119484. [Google Scholar] [CrossRef]
  96. Zheng, M.; Wang, Z.; Li, X.; Qiao, X.; Song, W.; Guo, L. Initial reaction mechanisms of cellulose pyrolysis revealed by ReaxFF molecular dynamics. Fuel 2016, 177, 130–141. [Google Scholar] [CrossRef]
  97. Diao, Z.; Zhao, Y.; Chen, B.; Duan, C.; Song, S. ReaxFF reactive force field for molecular dynamics simulations of epoxy resin thermal decomposition with model compound. J. Anal. Appl. Pyrolysis 2013, 104, 618–624. [Google Scholar] [CrossRef]
  98. Zhang, X.; Wu, Y.; Chen, X.; Wen, H.; Xiao, S. Theoretical study on decomposition mechanism of insulating epoxy resin cured by anhydride. Polymers 2017, 9, 341. [Google Scholar] [CrossRef] [PubMed]
  99. Li, G.; Hu, P.; Luo, W.; Zhang, J.; Yu, H.; Chen, F.; Zhang, F. Simulation of pyrolysis of crosslinked epoxy resin using ReaxFF molecular dynamics. Comput. Theor. Chem. 2021, 1200, 113240. [Google Scholar] [CrossRef]
  100. Natali, M.; Kenny, J.M.; Torre, L. Science and technology of polymeric ablative materials for thermal protection systems and propulsion devices: A review. Prog. Mater. Sci. 2016, 84, 192–275. [Google Scholar] [CrossRef]
  101. Castro-Marcano, F.; Kamat, A.M.; Russo, M.F., Jr.; van Duin, A.C.T.; Mathews, J.P. Combustion of an Illinois No. 6 coal char simulated using an atomistic char representation and the ReaxFF reactive force field. Combust. Flame 2012, 159, 1272–1285. [Google Scholar] [CrossRef]
  102. Odegard, G.M.; Jensen, B.D.; Gowtham, S.; Wu, J.; He, J.; Zhang, Z. Predicting mechanical response of crosslinked epoxy using ReaxFF. Chem. Phys. Lett. 2014, 591, 175–178. [Google Scholar] [CrossRef]
  103. Wang, S.; Guo, X.; Liang, T.; Zhou, Y.; Luo, Z. Mechanism research on cellulose pyrolysis by Py-GC/MS and subsequent density functional theory studies. Bioresour. Technol. 2012, 104, 722–728. [Google Scholar] [CrossRef]
  104. Paajamen, A.; Vaari, J. High-temperature decomposition of the cellulose molecule: A stochastic molecular dynamics study. Cellulose 2017, 24, 2713–2725. [Google Scholar] [CrossRef]
  105. Luo, Z.; Wang, S.; Liao, Y.; Cen, K. Mechanism study of cellulose rapid pyrolysis. Ind. Eng. Chem. Res. 2004, 43, 5605–5610. [Google Scholar] [CrossRef]
  106. Davin, L.B.; Lewis, N.G. Lignin primary structures and dirigent sites. Curr. Opin. Biotechnol. 2005, 16, 407–415. [Google Scholar] [CrossRef]
  107. Zhang, T.; Li, X.; Guo, L.; Guo, X. Reaction mechanisms in pyrolysis of hardwood, softwood, and kraft lignin revealed by ReaxFF MD simulations. Energy Fuels 2019, 33, 11210–11225. [Google Scholar] [CrossRef]
  108. Zhang, T.; Li, X.; Guo, L. Initial reactivity of linkages and monomer rings in lignin pyrolysis revealed by ReaxFF molecular dynamics. Langmuir 2017, 33, 11646–11657. [Google Scholar] [CrossRef] [PubMed]
  109. Santi, A.D.; Monti, S.; Barcaro, G.; Zhang, Z.; Barta, K. New mechanistic insights into the lignin β-O-4 linkage acidolysis with ethylene glycol stabilization aided by multilevel computational chemistry. ACS Sustain. Chem. Eng. 2021, 9, 2388–2399. [Google Scholar] [CrossRef] [PubMed]
  110. Gupta, S. Understanding Biomass Chemistry Using Multiscale Molecular Modeling Approach. In Catalysis for Clean Energy and Environmental Sustainability; Pant, K.K., Gupta, S.K., Ahmed, E., Eds.; Springer: Basel, Switzerland, 2021. [Google Scholar]
  111. Guidicianni, P.; Gargiulo, V.; Alfè, M.; Ragucci, R.; Ferreiro, A.I.; Rabaҫal, M.; Costa, M. Slow pyrolysis of xylan as pentose model compound for hardwood hemicellulose: A study of the catalytic effect of Na ions. J. Anal. Appl. Pyrolysis 2019, 137, 266–275. [Google Scholar] [CrossRef]
  112. Gargiulo, V.; Giudicianni, P.; Alfè, M.; Ragucci, R. About the influence of doping approach on the alkali metal catalyzed sloe pyrolysis of xylan. J. Chem. 2019, 2019, 9392571. [Google Scholar] [CrossRef]
  113. Gargiulo, V.; Ferreiro, A.I.; Giudicianni, P.; Tomaselli, S.; Costa, M.; Ragucci, R.; Alfe, M. Insights about the effect of composition, branching and molecular weight on the slow pyrolysis of xylose-based polysaccharides. J. Anal. Appl. Pyrolysis 2022, 161, 105369. [Google Scholar] [CrossRef]
  114. Shen, D.K.; Gu, S.; Bridgwater, A.V. Study on the pyrolytic behavior of xylan-based hemicellulose using TG-FTIR and Py-GC-FTIR. J. Anal. Appl. Pyrolysis 2010, 87, 199–206. [Google Scholar] [CrossRef]
  115. Chen, D.; Can, K.; Zhuang, X.; Gan, Z.; Zhou, J. Insight into biomass pyrolysis mechanism based on cellulose, hemicellulose, and lignin: Evolution of volatiles and kinetics, elucidation of reaction pathways, and characterization of gas, biochar and bio-oil. Combust. Flame 2022, 242, 112142. [Google Scholar] [CrossRef]
  116. Qu, T.; Guo, W.; Shen, L.; Xiao, J.; Zhao, K. Experimental study of biomass pyrolysis based on three major components: Hemicellulose, cellulose, and lignin. Ind. Eng. Chem. Res. 2011, 50, 10424–10433. [Google Scholar] [CrossRef]
  117. Zhou, X.; Li, W.; Mabon, R.; Broadbelt, L.J. A critical review on hemicellulose pyrolysis. Energy Technol. 2017, 5, 52–79. [Google Scholar] [CrossRef]
  118. Usino, D.O.; Supriyanto, P.Y.; Pettersson, A.; Richards, T. Influence of temperature and time on initial pyrolysis of cellulose and xylan. J. Anal. Appl. Pyrolysis 2020, 147, 104782. [Google Scholar] [CrossRef]
  119. Gao, Z.; Li, N.; Wang, Y.; Niu, W.; Yi, W. Pyrolysis behavior of xylan-based hemi-cellulose in a fixed bed reactor. J. Anal. Appl. Pyrolysis 2020, 146, 104772. [Google Scholar] [CrossRef]
  120. Patwardhan, P.R.; Brown, R.C.; Shanks, B.H. Product distribution from the fast pyrolysis of hemicellulose. ChemSusChem 2011, 4, 636–643. [Google Scholar] [CrossRef] [PubMed]
  121. Liu, Z.; Ku, X.; Jin, H.; Yang, S. Research on the microscopic reaction mechanism of cellulose pyrolysis using the molecular dynamics simulation. J. Anal. Appl. Pyrolysis 2021, 159, 105333. [Google Scholar] [CrossRef]
  122. Nimz, H. Beech lignin—Proposal of a constitutional scheme. Angew. Chem. Int. Ed. 1974, 13, 313–321. [Google Scholar] [CrossRef]
Figure 1. Some bond orders considered in ReaxFF potential forms.
Figure 1. Some bond orders considered in ReaxFF potential forms.
Jcs 07 00354 g001
Figure 2. Schematic representations of uncorrected and corrected bond order in terms of interatomic distance vs. carbon–carbon bond order to be considered.
Figure 2. Schematic representations of uncorrected and corrected bond order in terms of interatomic distance vs. carbon–carbon bond order to be considered.
Jcs 07 00354 g002
Figure 3. Schematic structure of cellulose (C6H10O5)n.
Figure 3. Schematic structure of cellulose (C6H10O5)n.
Jcs 07 00354 g003
Figure 4. Schematic basic units of lignin polymer: (a) coniferyl, (b) sinapyl, and (c) p-coumaryl alcohol structures.
Figure 4. Schematic basic units of lignin polymer: (a) coniferyl, (b) sinapyl, and (c) p-coumaryl alcohol structures.
Jcs 07 00354 g004
Figure 5. Schematic structure of lignin, where the purple wavy lines represent the continuing of the molecule.
Figure 5. Schematic structure of lignin, where the purple wavy lines represent the continuing of the molecule.
Jcs 07 00354 g005
Figure 6. Schematic structure of xylan, where the purple wavy lines represent the continuing of the molecule.
Figure 6. Schematic structure of xylan, where the purple wavy lines represent the continuing of the molecule.
Jcs 07 00354 g006
Table 1. Major physicochemical properties of biochar [80].
Table 1. Major physicochemical properties of biochar [80].
Pyrolysis TypeSSA (m2/g)CEC (cmol/kg)AEC (cmol/kg)CCE (%)PV (m3/t)APS (nm)Ash (%)pHEC (dS/m)
Fast183 ± 17.344.9 ± 3.624.90 ± 3.456.10 ± 1.122.04 ± 0.8152.3 ± 40.219.2 ± 19.28.7 ± 0.14.43 ± 0.50
Slow98.6 ± 3.5348.1 ± 3.125.33 ± 1.5111.2 ± 0.983.66 ± 1.271190 ± 56522.0 ± 0.518.7 ± 0.05.85 ± 1.58
Feedstock sourceSSA (m2/g)CEC (cmol/kg)AEC (cmol/kg)CCE (%)PV (m3/t)APS (nm)Ash (%)pHEC
(dS/m)
Wood based184 ± 11.423.9 ± 1.875.65 ± 1.809.04 ± 1.177.01 ± 3.0774.6 ± 44.410.2 ± 0.438.3 ± 0.16.20 ± 2.85
Crop wastes98.2 ± 5.4556.3 ± 3.924.51 ± 1.966.12 ± 0.972.05 ± 0.912320 ± 115021.1 ± 0.548.9 ± 0.15.72 ± 0.67
Other grasses63.4 ± 8.8463.316.42.05 ± 1.053.36 ± 3.30268 ± 12518.0 ± 1.018.9 ± 0.15.20 ± 0.93
Manures/biosolids52.2 ± 4.2366.1 ± 8.007.77 ± 7.5214.2 ± 1.560.82 ± 0.3027.3 ± 12.544.6 ± 0.978.9 ± 0.13.98 ± 0.41
SSA—specific surface area; CEC—cation exchange capacity; AEC—anion exchange capacity; CCE—calcium carbonate equivalent; PV—total pore volume; APS—average particle size; EC—electrical conductivity.
Table 2. The major elemental compositions of biochar [79,80].
Table 2. The major elemental compositions of biochar [79,80].
Pyrolysis TypeC (wt.%)H (wt.%)O (wt.%)N (wt.%)S (wt.%)
Fast60.6 ± 0.473.37 ± 0.0819.1 ± 0.381.63 ± 0.060.085 ± 0.009
Slow60.8 ± 0.343.36 ± 0.0918.4 ± 0.291.63 ± 0.040.055 ± 0.004
Feedstock sourceC (wt.%)H (wt.%)O (wt.%)N (wt.%)S (wt.%)
Wood based70.5 ± 0.393.38 ± 0.0817.7 ± 0.350.95 ± 0.030.044 ± 0.007
Crop wastes61.4 ± 0.413.28 ± 0.1018.1 ± 0.381.54 ± 0.060.039 ± 0.006
Other grasses63.6 ± 0.725.11 ± 0.5020.9 ± 0.741.80 ± 0.140.051 ± 0.021
Manures/biosolids41.6 ± 0.682.73 ± 0.1016.5 ± 0.702.42 ± 0.060.089 ± 0.006
Corncobs/cassava rhizomes/cassava stems62.95–81.352.24–2.7315.23–33.441.22–1.65
Table 3. Constituent details of four lignin models [76,107].
Table 3. Constituent details of four lignin models [76,107].
Lignin ModelSpeciesConstituentsNumber of C6H3 Units
HGS
AlderSoftwoodC6400H7200O23204056040
FreudenbergSoftwoodC6980H7640O228020048040
NimzHardwoodC10854H11940O406242612396
MartonKraft ligninC5080H5040O1640S401204400
Table 4. Percentage of each gas produced resulting from the three main components of the biomass (wheat straw) during the pyrolysis process at 2000 K.
Table 4. Percentage of each gas produced resulting from the three main components of the biomass (wheat straw) during the pyrolysis process at 2000 K.
CO (%)CO2 (%)CH4 (%)H2O (%)
Cellulose48.2947.0453.3383.96
Lignin7.324.6840.005.30
Hemicellulose44.3948.286.6710.74
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Hu, Z.; Wei, L. Review on Characterization of Biochar Derived from Biomass Pyrolysis via Reactive Molecular Dynamics Simulations. J. Compos. Sci. 2023, 7, 354. https://doi.org/10.3390/jcs7090354

AMA Style

Hu Z, Wei L. Review on Characterization of Biochar Derived from Biomass Pyrolysis via Reactive Molecular Dynamics Simulations. Journal of Composites Science. 2023; 7(9):354. https://doi.org/10.3390/jcs7090354

Chicago/Turabian Style

Hu, Zhong, and Lin Wei. 2023. "Review on Characterization of Biochar Derived from Biomass Pyrolysis via Reactive Molecular Dynamics Simulations" Journal of Composites Science 7, no. 9: 354. https://doi.org/10.3390/jcs7090354

APA Style

Hu, Z., & Wei, L. (2023). Review on Characterization of Biochar Derived from Biomass Pyrolysis via Reactive Molecular Dynamics Simulations. Journal of Composites Science, 7(9), 354. https://doi.org/10.3390/jcs7090354

Article Metrics

Back to TopTop