ReaxFF study of the oxidation of lignin model compounds for the most common linkages in softwood in view of carbon fiber production

Tracking Thermo-Oxidation Reaction Products and Pathways of Modified Lignin Structures from Reactive Molecular Dynamics Simulations

Thermo-oxidation of biomass is an important process that occurs through a variety of reaction pathways depending on the chemical nature of the molecules and reaction conditions. These processes can be modeled using reactive molecular dynamics to study chemical reactions and the evolution of converted molecules over time. The advantage of this approach is that many molecules can be modeled, but it is challenging to use the large amount of data obtained from such a simulation to determine reaction products and pathways. In this study, we developed a tracking approach to identify the reaction pathways of the dominant reaction products from reactive molecular dynamics simulations. We demonstrated the approach for thermo-oxidation reactions of modified model lignin compounds. For two modified lignin structures, we tracked the evolving chemical species to find the most common reaction products. Subsequently, we monitored specific bonds to determine the individual steps in the reaction process. This combined approach of reactive molecular dynamics and tracking enabled us to identify the most likely thermo-oxidation pathways. The methodology can be used to investigate the thermo-oxidative pathways of a wider range of chemical compounds.

Simulation and experimental evaluation of laser-induced graphene on the cellulose and lignin substrates

In this article, the formation of laser-induced graphene on the two natural polymers, cellulose, and lignin, as precursors was investigated with molecular dynamics simulations and some experiments. These eco-friendly polymers provide significant industrial advantages due to their low cost, biodegradability, and recyclable aspects. It was discovered during the simulation that LIG has numerous defects and a porous structure. Carbon monoxide, H2, and water vapor are gases released by cellulose and lignin substrates. H2O and CO are released when the polymer transforms into an amorphous structure. Later on, as the amorphous structure changes into an ordered graphitic structure, H2 is released continuously. Since cellulose monomer has a higher mass proportion of oxygen (49%) than lignin monomer (29%), it emits more CO. The LIG structure contains many 5- and 7-carbon rings, which cause the structure to have bends and undulations that go out of the plane. In addition, to verify the molecular dynamics simulation results with experimental tests, we used a carbon dioxide laser to transform filter paper, as a cellulose material, and coconut shell, as a lignin material, into graphene. Surprisingly, empirical experiments confirmed the simulation results.

Detailed mechanism study of volatile organic compound decomposition and oxidation removal based on a ReaxFF MD method

Volatile organic compounds (VOCs) are typical air pollutants as well as gaseous wastes that contain energy. Utilization and disposition of VOCs is currently an important research hotspot in the field of atmospheric environment. In this paper, the thermal cracking and oxidation reaction processes of typical VOCs components were modelled and analyzed by combining molecular dynamics and detailed reaction mechanisms, focusing on the effects of temperature, oxygen and other conditions on the conversion of VOCs. The results of molecular dynamics studies show that improving temperature and reaction time benefit the decomposition of VOCs. High temperatures under an inert atmosphere can sufficiently crack the VOCs themselves, but other by-products are generated, which in turn cause secondary pollution. The activation energies derived by ReaxFF-MD calculation are 328 kJ mol-1, 147 kJ mol-1 and 121 kJ mol-1 for toluene, styrene and benzaldehyde respectively, which is consistent with experimental results. Under the oxygen atmosphere, the conversion rate of VOCs is greatly increased and the reaction temperature is significantly reduced. Meanwhile, the oxidation reaction fully converts VOCs into non-polluting products such as CO2 and H2O. Detailed kinetic studies show that initial oxidation of toluene molecules raised by hydrogen abstraction reaction is the dominant step during toluene oxidation, which significantly improved the decomposition efficiency of toluene.

Fiber Melt Spinning and Thermo-Stabilization of Para-Rubber Wood Lignin: An Approach for Fully Biomass Precursor Preparation

Para-rubber wood (PRW) lignin, extracted from agricultural waste, was successfully melt-spun to fibers and thermo-stabilized without employing auxiliary additives. 31P NMR analysis revealed that PRW-lignin contained mainly a syringyl unit of phenolic C5-substituted OH group, which enabled melt flow during fiber spinning, as well as a guaiacyl unit which offered the ability to cross-link during thermo-stabilization. Thermo-stabilized fibers with no fusion were achieved at 250 °C with the heating rate of 0.1 °C/min. Structural changes in the fibers during stabilization were systematically investigated using FTIR and XPS analyses. From the results, changes in the intensities of characteristic bands relating to C-H stretching, aromatic C-H stretching, and C=O stretching indicated structural changes of lignin toward aromaticity via oxidation reactions. XPS analysis of the fibers carbonized at 900, 1000, and 1200 °C revealed an increase in carbon content from 72 to 87 wt %. and a decrease in oxygen content from 28 to 13 wt %. with the increasing carbonization temperature. The weight loss of carbonized fibers was in the range of 73.6 to 88.7%. The high weight loss of fibers carbonized at 1200 °C was explained partly due to the thermal decomposition of disordered carbon. The tensile strength and modulus of carbonized fibers were 163.0 and 275.1 MPa, respectively. This study demonstrates an approach to prepare a fully biomass precursor fiber and contributes to the exploration of the potential use of lignin from biomass waste.

Predicting molecule size distribution in hydrocarbon pyrolysis using random graph theory

Hydrocarbon pyrolysis is a complex process involving large numbers of chemical species and types of chemical reactions. Its quantitative description is important for planetary sciences, in particular, for understanding the processes occurring in the interior of icy planets, such as Uranus and Neptune, where small hydrocarbons are subjected to high temperature and pressure. We propose a computationally cheap methodology based on an originally developed ten-reaction model and the configurational model from random graph theory. This methodology generates accurate predictions for molecule size distributions for a variety of initial chemical compositions and temperatures ranging from 3200 to 5000 K. Specifically, we show that the size distribution of small molecules is particularly well predicted, and the size of the largest molecule can be accurately predicted provided that this molecule is not too large.

Temperature Extrapolation of Molecular Dynamics Simulations of Complex Chemistry to Microsecond Timescales Using Kinetic Models: Applications to Hydrocarbon Pyrolysis

We develop a method to construct temperature-dependent kinetic models of hydrocarbon pyrolysis, based on information from molecular dynamics (MD) simulations of pyrolyzing systems in the high-temperature regime. MD simulations are currently a key tool to understand the mechanism of complex chemical processes such as pyrolysis and to observe their outcomes in different conditions, but these simulations are computationally expensive and typically limited to nanoseconds of simulation time. This limitation is inconsequential at high temperatures, where equilibrium is reached quickly, but at low temperatures, the system may not equilibrate within a tractable simulation timescale. In this work, we develop a method to construct kinetic models of hydrocarbon pyrolysis using the information from the high-temperature high-reactivity regime. We then extrapolate this model to low temperatures, which enables microsecond-long simulations to be performed. We show that this approach accurately predicts the time evolution of small molecules, as well as the size and composition of long carbon chains across a wide range of temperatures and compositions. Further, we show that the range of suitable temperatures for extrapolation can easily be improved by adding more simulations to the training data. Compared to experimental results, our kinetic model leads to similar compositional trends while allowing for more detailed kinetic and mechanistic insights.

A Review on the Effect of the Mechanism of Organic Polymers on Pellet Properties for Iron Ore Beneficiation

Iron ore pellets not only have excellent metallurgical and mechanical properties but are also essential raw materials for improving iron and steel smelting in the context of the increasing global depletion of high-grade iron ore resources. Organic polymers, as important additive components for the production of high-quality pellets, have a significant impact on the formation as well as the properties of pellets. In this review, the mechanisms of organic polymers on the pelletizing properties, bursting temperature, and pellet strength at low and high temperatures, as well as the existing measures and mechanisms to improve the high-temperature strength of the organic binder pellets are systematically summarized. Compared with traditional bentonite additives, the organic polymers greatly improve the pelletizing rate and pellet strength at low temperatures, and significantly reduces metallurgical pollution. However, organic binders often lead to a decrease in pellet bursting temperature and pellet strength at high temperatures, which can be significantly improved by compounding with a small amount of low-cost inorganic minerals, such as bentonite, boron-containing compounds, sodium salts, and copper slag. At the same time, some industrial solid wastes can be rationally used to reduce the cost of pellet binders.

Pyrolysis Mechanism of Wheat Straw Based on ReaxFF Molecular Dynamics Simulations

Biomass has played an increasingly important role in the consumption of energy worldwide because of its renewability and carbon-neutral property. In this work, the pyrolysis mechanism of wheat straw is explored using reactive force field molecular dynamics simulations. A large-scale wheat straw model composed of cellulose, hemicellulose, and lignin is built. After model validation, the temporal evolutions of the main pyrolysis products under different temperatures are analyzed. As the temperature rises, the gas production increases and the tar yield can decrease after peaking. Relatively high temperatures accelerate the generation rates of the main gas and tar species. CO and CO₂ molecules mainly come from the cleavage of CHO₂ radicals, and numerous H₂O molecules are generated on account of dehydration. Moreover, the evolution of six functional groups and pyran and phenyl rings as well as three types of bonds is also presented. It is observed that the phenyl rings reflect improved thermostability. Finally, the pyrolytic kinetics analysis is conducted, and the estimated activation energy of wheat straw pyrolysis is found to be 56.19 kJ/mol. All these observations can help deeply understand the pyrolytic mechanism of wheat straw biomass.

Atoms to fibers: Identifying novel processing methods in the synthesis of pitch-based carbon fibers

Understanding and optimizing the key mechanisms used in the synthesis of pitch-based carbon fibers (CFs) are challenging, because unlike polyacrylonitrile-based CFs, the feedstock for pitch-based CFs is chemically heterogeneous, resulting in complex fabrication leading to inconsistency in the final properties. In this work, we use molecular dynamics simulations to explore the processing and chemical phase space through a framework of CF models to identify their effects on elastic performance. The results are in excellent agreement with experiments. We find that density, followed by alignment, and functionality of the molecular constituents dictate the CF mechanical properties more strongly than their size and shape. Last, we propose a previously unexplored fabrication route for high-modulus CFs. Unlike graphitization, this results in increased sp³ fraction, achieved via generating high-density CFs. In addition, the high sp³ fraction leads to the fabrication of CFs with isometric compressive and tensile moduli, enabling their potential applications for compressive loading.

New Mechanistic Insights into the Lignin β-O-4 Linkage Acidolysis with Ethylene Glycol Stabilization Aided by Multilevel Computational Chemistry

Acidolysis in conjunction with stabilization of reactive intermediates has emerged as one of the most powerful methods of lignin depolymerization that leads to high aromatic monomer yields. In particular, stabilization of reactive aldehydes using ethylene glycol results in the selective formation of the corresponding cyclic acetals (1,3-dioxolane derivatives) from model compounds, lignin, and even from softwood lignocellulose. Given the high practical utility of this method for future biorefineries, a deeper understanding of the method is desired. Here, we aim to elucidate key mechanistic questions utilizing a combination of experimental and multilevel computational approaches. The multiscale computational protocol used, based on ReaxFF molecular dynamics, represents a realistic scenario, where a typical experimental setup can be reproduced confidently given the explicit molecules of the solute, catalyst, and reagent. The nudged elastic band (NEB) approach allowed us to characterize the key intermolecular interactions involved in the reaction paths leading to crucial intermediates and products. The high level of detail obtained clearly revealed for the first time the unique role of sulfuric acid as a proton donor and acceptor in lignin β-O-4 acidolysis as well as the reaction pathways for ethylene glycol stabilization, and the difference in reactivity between compounds with different methoxy substituents.

One-step alcoholysis of lignin into small-molecular aromatics: Influence of temperature, solvent, and catalyst

•

The reactant suspension mode is an effective strategy to deoxy-liquefaction of lignin.

•

The catalyst Cu-C has the optimal catalytic activity and selectivity in methanol.

•

The catalyst Fe-SiC possesses the optimal catalytic deoxygenation in ethanol.

•

The cleavages of C—O ether bonds and C—C bonds directly promote the formation of small-molecular aromatics.

Lignin valorization is a challenge because of its complex structure and high thermal stability. Supercritical alcoholysis of lignin without external hydrogen in a self-made high-pressure reactor is investigated under different temperatures (450–500 °C) and solvents as well as catalysts by using a reactant suspension mode. Small-molecular arenes and mono-phenols (C₇-C₁₂) are generated under short residence time of 30 min. High temperature (500 °C) favors efficient deoxy-liquefaction of lignin (70%) and formation of small-molecular arenes (C₆-C₉). Solvents methanol and ethanol demonstrate much more synergistic effect on efficient deoxy-liquefaction of lignin than propanol. The catalyst Cu-C has the optimal activity and selectivity in methanol (70% of conversion, 83.93% of arenes), whereas Fe-SiC possesses the optimal catalytic deoxygenation in ethanol, resulting in the formation of arenes other than phenols. Further analysis indicates that lignin is converted into arenes by efficient cleavages of C—O ether bonds and C—C bonds under high temperature and pressure.

ReaxFF Simulations of Lignin Fragmentation on a Palladium-Based Heterogeneous Catalyst in Methanol-Water Solution

The interaction of fragments derived from lignin depolymerization with a heterogeneous palladium catalyst in methanol-water solution is studied by means of experimental and theoretical methodologies. Quantum chemistry calculations and molecular dynamics simulations based on the ReaxFF approach are combined effectively to obtain an atomic level characterization of the crucial steps of the adsorption of the molecules on the catalyst, their fragmentation, reactions, and desorption. The main products are identified, and the most important routes to obtain them are explained through extensive computational procedures. The simulation results are in excellent agreement with the experiments and suggest that the mechanisms comprise a fast chemisorption of identified fragments from lignin on the metal interface accompanied by bond breaking, release of some of their hydrogens and oxygens to the support, and eventual desorption depending on the local environment. The strongest connections are those involving the aromatic rings, as confirmed by the binding energies of selected representative structures, estimated at the quantum chemistry level. The satisfactory agreement with the literature, quantum chemistry data, and experiments confirms the reliability of the multilevel computational procedure to study complex reaction mixtures and its potential application in the design of high-performance catalytic devices.

The Rise of Hierarchical Nanostructured Materials from Renewable Sources: Learning from Nature

Mimicking Nature implies the use of bio-inspired hierarchical designs to manufacture nanostructured materials. Such materials should be produced from sustainable sources ( e.g., biomass) and through simple processes that use mild conditions, enabling sustainable solutions. The combination of different types of nanomaterials and the implementation of different features at different length scales can provide synthetic hierarchical nanostructures that mimic natural materials, outperforming the properties of their constitutive building blocks. Taking recent developments in flow-assisted assembly of nanocellulose crystals as a starting point, we review the state of the art and provide future perspectives on the manufacture of hierarchical nanostructured materials from sustainable sources, assembly techniques, and potential applications.

Molecular Dynamics Simulations on the Thermal Decomposition of Meta-Aramid Fibers

The thermal decomposition mechanism of a meta-aramid fiber was simulated at the atomic level using the ReaxFF reactive force field. The simulation results indicated that the main initial decomposition positions of the meta-aramid fiber elements were Caromatic ring⁻N and C=O, which could be used as targets for the modification of meta-aramid fibers. The meta-aramid fiber elements first decomposed into C6⁻C13 and then into smaller segments and micromolecular gases. The temperature was shown to be the key factor affecting the thermal decomposition of the meta-aramid fibers. More complex compositions and stable gases were produced at high temperatures than at lower temperatures. HCN was a decomposition product at high temperature, suggesting that its presence could be used for detecting thermal faults in meta-aramid fibers. Generation path tracing of the thermal decomposition products NH₃ and H₂O was also performed. NH₃ was produced when the NH₂ group captured an H atom adjacent to the system. H₂O was formed after a carbonyl group captured an H atom, became a hydroxyl group, with subsequent intramolecular dehydration or intermolecular hydrogen abstraction.

Persistent free radicals in carbon-based materials on transformation of refractory organic contaminants (ROCs) in water: A critical review

With the increased concentrations and kinds of refractory organic contaminants (ROCs) in aquatic environments, many previous reviews systematically summarized the applications of carbon-based materials in the adsorption and catalytic degradation of ROCs for their economically viable and environmentally friendly behavior. Interestingly, recent studies indicated that carbon-based materials in natural environment can also mediate the transformation of ROCs directly or indirectly due to their abundant persistent free radicals (PFRs). Understanding the formation mechanisms of PFRs in carbo-based materials and their interactions with ROCs is essential to develop their further applications in environment remediation. However, there is no comprehensive review so far about the direct and indirect removal of ROCs mediated by PFRs in amorphous, porous and crystalline carbon-based materials. The review aims to evaluate the formation mechanisms of PFRs in carbon-based materials synthesized through pyrolysis and hydrothermal carbonization processes. The influence of synthesis conditions (temperature and time) and carbon sources on the types as well as the concentrations of PFRs in carbon-based materials are also discussed. In particular, the effects of metals on the concentrations and types of PFRs in carbon-based materials are highlighted because they are considered as the catalysts for the formation of PFRs. The formation mechanisms of reactive species and the further transformation mechanisms of ROCs are briefly summarized, and the surface properties of carbon-based materials including surface area, types and number of functional groups, etc. are found to be the key parameters controlling their activities. However, due to diversity and complexity of carbon-based materials, the exact relationships between the activities of carbon-based materials and PFRs are still uncertain. Finally, the existing problems and current challenges for the ROCs transformation with carbon-based materials are also pointed out.

Depolymerization Pathways for Branching Lignin Spirodienone Units Revealed with ab Initio Steered Molecular Dynamics

Lignocellulosic biomass is an abundant, rich source of aromatic compounds, but direct utilization of raw lignin has been hampered by both the high heterogeneity and variability of linking bonds in this biopolymer. Ab initio steered molecular dynamics (AISMD) has emerged both as a fruitful direct computational screening approach to identify products that occur through mechanical depolymerization (i.e., in sonication or ball-milling) and as a sampling approach. By varying the direction of force and sampling over 750 AISMD trajectories, we identify numerous possible pathways through which lignin depolymerization may occur in pyrolysis or through catalytic depolymerization as well. Here, we present eight unique major depolymerization pathways discovered via AISMD for the recently characterized spirodienone lignin branching linkage that may comprise around 10% weight of all lignin in some softwoods. We extract representative trajectories from AISMD and carry out reaction pathway analysis to identify energetically favorable pathways for lignin depolymerization. Importantly, we identify dynamical effects that could not be observed through more traditional calculations of bond dissociation energies. Such effects include thermodynamically favorable recovery of aromaticity in the dienone ring that leads to near-barrierless subsequent ether cleavage and hydrogen-bonding effects that stabilize newly formed radicals. Some of the most stable spirodienone fragments that reside at most 1 eV above the reactant structure are formed with only 2 eV barriers for C-C bond cleavage, suggesting key targets for catalyst design to drive targeted depolymerization of lignin.

Conformations of Low-Molecular-Weight Lignin Polymers in Water

Low-molecular-weight lignin binds to cellulose during the thermochemical pretreatment of biomass for biofuel production, which prevents the efficient hydrolysis of the cellulose to sugars. The binding properties of lignin are influenced strongly by the conformations it adopts. Here, we use molecular dynamics simulations in aqueous solution to investigate the dependence of the shape of lignin polymers on chain length and temperature. Lignin is found to adopt collapsed conformations in water at 300 and 500 K. However, at 300 K, a discontinuous transition is found in the shape of the polymer as a function of the chain length. Below a critical degree of polymerization, Nc =15, the polymer adopts less spherical conformations than above Nc. The transition disappears at high temperatures (500 K) at which only spherical shapes are adopted. An implication relevant to cellulosic biofuel production is that lignin will self-aggregate even at high pretreatment temperatures.

Ab Initio Screening Approach for the Discovery of Lignin Polymer Breaking Pathways

The directed depolymerization of lignin biopolymers is of utmost relevance for the valorization or commercialization of biomass fuels. We present a computational and theoretical screening approach to identify potential cleavage pathways and resulting fragments that are formed during depolymerization of lignin oligomers containing two to six monomers. We have developed a chemical discovery technique to identify the chemically relevant putative fragments in eight known polymeric linkage types of lignin. Obtaining these structures is a crucial precursor to the development of any further kinetic modeling. We have developed this approach by adapting steered molecular dynamics calculations under constant force and varying the points of applied force in the molecule to diversify the screening approach. Key observations include relationships between abundance and breaking frequency, the relative diversity of potential pathways for a given linkage, and the observation that readily cleaved bonds can destabilize adjacent bonds, causing subsequent automatic cleavage.

Why genetic modification of lignin leads to low-recalcitrance biomass

Molecular dynamics simulations show genetically modified lignin to associate less with hemicellulose than does wild type.

Lignin hydrolysis and phosphorylation mechanism during phosphoric acid-acetone pretreatment: a DFT study

The study focused on the structural sensitivity of lignin during the phosphoric acid–acetone pretreatment process and the resulting hydrolysis and phosphorylation reaction mechanisms using density functional theory calculations. The chemical stabilities of the seven most common linkages (β-O-4, β-β, 4-O-5, β-1, 5-5, α-O-4, and β-5) of lignin in H₃PO₄, CH₃COCH₃, and H₂O solutions were detected, which shows that α-O-4 linkage and β-O-4 linkage tend to break during the phosphoric acid–acetone pretreatment process. Then α-O-4 phosphorylation and β-O-4 phosphorylation follow a two-step reaction mechanism in the acid treatment step, respectively. However, since phosphorylation of α-O-4 is more energetically accessible than phosphorylation of β-O-4 in phosphoric acid, the phosphorylation of α-O-4 could be controllably realized under certain operational conditions, which could tune the electron and hole transfer on the right side of β-O-4 in the H₂PO₄⁻ functionalized lignin. The results provide a fundamental understanding for process-controlled modification of lignin and the potential novel applications in lignin-based imprinted polymers, sensors, and molecular devices.