Computational Prediction of α/β Selectivities in the Pyrolysis of Oxygen-Substituted Phenethyl Phenyl Ethers

Initial Reactivity of Linkages and Monomer Rings in Lignin Pyrolysis Revealed by ReaxFF Molecular Dynamics

The initial conversion pathways of linkages and their linked monomer units in lignin pyrolysis were investigated comprehensively by ReaxFF MD simulations facilitated by the unique VARxMD for reaction analysis. The simulated molecular model contains 15 920 atoms and was constructed on the basis of Adler's softwood lignin model. The simulations uncover the initial conversion ratio of various linkages and their linked aryl monomers. For linkages and their linked monomer aryl rings of α-O-4, β-O-4 and α-O-4 & β-5, the C_α/C_β ether bond cracking dominates the initial pathway accounting for at least up to 80% of their consumption. For the linkage of β-β & γ-O-α, both the C_α-O ether bond cracking and its linked monomer aryl ring opening are equally important. Ring-opening reactions dominate the initial consumption of other 4-O-5, 5-5, β-1, β-2, and β-5 linkages and their linked monomers. The ether bond cracking of C_α-O and C_β-O occurs at low temperature, and the aryl ring-opening reactions take place at relatively high temperature. The important intermediates leading to the stable aryl ring opening are the phenoxy radicals, the bridged five-membered and three-membered rings and the bridged six-membered and three-membered rings. In addition, the reactivity of a linkage and its monomer aryl ring may be affected by other linkages. The ether bond cracking of α-O-4 and β-O-4 linkages can activate its neighboring linkage or monomer ring through the formed phenoxy radicals as intermediates. The important intermediates revealed in this article should be of help in deepening the understanding of the controlling mechanism for producing aromatic chemicals from lignin pyrolysis.

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.

Coupling and Reactions of 5-Hydroxyconiferyl Alcohol in Lignin Formation

The catechol alcohols, caffeyl and 5-hydroxyconiferyl alcohol, may be incorporated into lignin either naturally or through genetic manipulation. Due to the presence of o-OH groups, these compounds form benzodioxanes, a departure from the interunit connections found in lignins derived from the cinnamyl alcohols. In nature, lignins composed of caffeyl and 5-hydroxyconiferyl alcohol are linear homopolymers and, as such, may have properties that make them amenable for use in value-added products, such as lignin-based carbon fibers. In the current work, results from density functional theory calculations for the reactions of 5-hydroxyconiferyl alcohol, taking stereochemistry into account, are reported. Dehydrogenation and quinone methide formation are found to be thermodynamically favored for 5-hydroxyconiferyl alcohol, over coniferyl alcohol. The comparative energetics of the rearomatization reactions suggest that the formation of the benzodioxane linkage is under kinetic control. Ring-opening reactions of the benzodioxane groups show that the bond dissociation enthalpy of the α-O cleavage reaction is lower than that of the β-O reaction. The catechol lignins represent a novel form of the polymer that may offer new opportunities for bioproducts and genetic targets.

Multicomponent kinetic analysis and theoretical studies on the phenolic intermediates in the oxidation of eugenol and isoeugenol catalyzed by laccase

Laccase catalyzes the oxidation of natural phenols and thereby is believed to initialize reactions in lignification and delignification. Numerous phenolic mediators have also been applied in laccase-mediator systems. However, reaction details after the primary O-H rupture of phenols remain obscure. In this work two types of isomeric phenols, EUG (eugenol) and ISO (trans-/cis-isoeugenol), were used as chemical probes to explore the enzymatic reaction pathways, with the combined methods of time-resolved UV-Vis absorption spectra, MCR-ALS, HPLC-MS, and quantum mechanical (QM) calculations. It has been found that the EUG-consuming rate is linear to its concentration, while the ISO not. Besides, an o-methoxy quinone methide intermediate, (E/Z)-4-allylidene-2-methoxycyclohexa-2,5-dienone, was evidenced in the case of EUG with the UV-Vis measurement, mass spectra and TD-DFT calculations; in contrast, an ISO-generating phenoxyl radical, a (E/Z)-2-methoxy-4-(prop-1-en-1-yl) phenoxyl radical, was identified in the case of ISO. Furthermore, QM calculations indicated that the EUG-generating phenoxyl radical (an O-centered radical) can easily transform into an allylic radical (a C-centered radical) by hydrogen atom transfer (HAT) with a calculated activation enthalpy of 5.3 kcal mol(-1) and then be fast oxidized to the observed eugenol quinone methide, rather than an O-radical alkene addition with barriers above 12.8 kcal mol(-1). In contrast, the ISO-generating phenoxyl radical directly undergoes a radical coupling (RC) process, with a barrier of 4.8 kcal mol(-1), while the HAT isomerization between O- and C-centered radicals has a higher reaction barrier of 8.0 kcal mol(-1). The electronic conjugation of the benzyl-type radical and the aromatic allylic radical leads to differentiation of the two pathways. These results imply that competitive reaction pathways exist for the nascent reactive intermediates generated in the laccase-catalyzed oxidation of natural phenols, which is important for understanding the lignin polymerization and may shed some light on the development of efficient laccase-mediator systems.

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.

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.

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

Lignin is an underused but major component of biomass. One possible area of utilization is the production of carbon fiber. A necessary processing step is the stabilization of lignin fiber (typically in an oxygen environment) before high temperature treatment. We investigate oxidative, thermal conversion of lignin using computational methods. Dilignol model compounds for the most common (seven) linkages in softwood are chosen to represent the diverse structure of lignin. We perform molecular dynamics simulation where the potential energy surface is described by a reactive force field (ReaxFF). We calculate overall activation energies for model conversion and reveal initial mechanisms of formaldehyde formation. We record fragmentation patterns and average carbon oxidation numbers at various temperatures. Most importantly, we identify mechanisms for stabilizing reactions that result in cyclic and rigid connections in softwood lignin fibers that are necessary for further processing into carbon fibers.

Pyrolysis of phenethyl phenyl ether tethered in mesoporous silica. Effects of confinement and surface spacer molecules on product selectivity

There has been expanding interest in exploring porous metal oxides as a confining environment for organic molecules resulting in altered chemical and physical properties including chemical transformations. In this paper, we examine the pyrolysis behavior of phenethyl phenyl ether (PPE) confined in mesoporous silica by covalent tethers to the pore walls as a function of tether density and the presence of cotethered surface spacer molecules of varying structure (biphenyl, naphthyl, octyl, and hexadecyl). The PPE pyrolysis product selectivity, which is determined by two competitive free-radical pathways cycling through the two aliphatic radical intermediates (PhCH·CH(2)OPh and PhCH(2)CH·OPh), is shown to be significantly different from that measured in the liquid phase as well as for PPE tethered to the exterior surface of nonporous silica nanoparticles. Tailoring the pore surface with spacer molecules further alters the selectivity such that the PPE reaction channel involving a molecular rearrangement (O-C phenyl shift in PhCH(2)CH·OPh), which accounts for 25% of the products in the liquid phase, can be virtually eliminated under pore confinement conditions. The origin of this change in selectivity is discussed in the context of steric constraints on the rearrangement path inside the pores, surface and pore confinement effects, pore surface curvature, and hydrogen bonding of PPE with residual surface silanols supplemented by nitrogen physisorption data and molecular dynamics simulations.