The significance of intra-molecular hydrogen bonding in the ß-O-4 linkage of lignin

OH-Initiated Reactions of para-Coumaryl Alcohol Relevant to the Lignin Pyrolysis. Part III. Kinetics of H-Abstraction by H, OH, and CH3 Radicals

Lignin is the most complex component of biomass, and development of a detailed chemical kinetic model for biomass pyrolysis mainly relies on the understanding of the lignin decomposition kinetics. para-Coumaryl alcohol (p-CMA, HOPh-CH═CH-CH₂OH), the focus of our analysis, is the simplest of the lignin monomers (monolignols) containing a typical side-chain double bond and both alkyl- and phenolic-type OH-groups. In parts I and II of our work (Asatryan, R. J. Phys. Chem. A 2019, 123, 2570-2585; Hudzik, J. M. J. Phys. Chem. A 2020, current issue), we created a detailed potential energy surface (PES) and performed a kinetic analysis of chemically activated, unimolecular, and bimolecular reactions pathways for p-CMA + OH. Reaction pathways analyzed include dissociation, intramolecular abstraction, group transfer, and elimination processes. The α- and β-carbon addition reactions generate 1,3- (RA1) and 1,2-diol (RB1) adduct radicals, respectively. Well depths are approximately 29 and 41 kcal/mol below the p-CMA + OH entrance level. Kinetic analysis aides in determining the major pathways for our conventional and fractional pyrolysis experiments. The current paper focuses on the H-abstraction reactions via H, OH, and CH₃ light ("pool") radicals from p-CMA. The thermochemical properties of all stable, radical, and transition-state species were determined using the ωB97XD density functional theory (DFT) and higher-level CBS-QB3 composite methods. Barrier heights from the prereaction complexes, for OH-radical abstractions, to the transition states for the propanoid side chain are compared to the model H-abstraction reactions of allyl alcohol (AA) with OH and p-CMA with H and CH₃ radicals. The lowest-energy, most stable, p-CMA radical formed is at the C9 allylic position (p-CMA-C9j) with exothermicity of 26.63, 41.32, and 27.34 kcal/mol for H, OH, and CH₃, respectively. For OH-radical abstraction at this position, our findings are consistent with corresponding data on AA + OH at 37.44 kcal/mol and similar to that of RB1. A similar stable radical with an exothermicity of 34.95 kcal/mol occurs for the phenol hydroxyl group, generating the p-CMA-O4j radical. H-abstraction pathways are considered in relation to other major pathways previously considered for p-CMA + OH reactions including H-atom shifts, dehydration, and β-scission reactions. Derived rate coefficients for substituted phenols can be utilized in detailed kinetic models for lignin/biomass pyrolysis.

The Self-Assembly of Lignin and Its Application in Nanoparticle Synthesis: A Short Review

Lignin serves as a significant contributor to the natural stock of non-fossilized carbon, second only to cellulose in the biosphere. In this review article, we focus on the self-assembly properties of lignin and their contribution to its effective utilization and valorization. Traditionally, investigations on self-assembly properties of lignin have aimed at understanding the lignification process of the cell wall and using it for efficient delignification for commercial purposes. In recent years (mainly the last three years), an increased number of attempts and reports of technical-lignin nanostructure synthesis with controlled particle size and morphology have been published. This has renewed the interests in the self-assembly properties of technical lignins and their possible applications. Based on the sources and processing methods of lignin, there are significant differences between its structure and properties, which is the primary obstacle in the generalized understanding of the lignin structure and the lignification process occurring within cell walls. The reported studies are also specific to source and processing methods. This work has been divided into two parts. In the first part, the aggregation propensity of lignin based on type, source and extraction method, temperature, and pH of solution is discussed. This is followed by a critical overview of non-covalent interactions and their contribution to the self-associative properties of lignin. The role of self-assembly towards the understanding of xylogenesis and nanoparticle synthesis is also discussed. A particular emphasis is placed on the interaction and forces involved that are used to explain the self-association of lignin.

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.

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.

Spectroscopic investigation, computed IR intensity, Raman activity and vibrational frequency analysis on 3-bromoanisole using HF and DFT (LSDA/MPW1PW91) calculations

In this work, the experimental and theoretical spectra of 3-bromoanisole (3-BA) are studied. FT-IR and FT-Raman spectra of title molecule in the liquid phase have been recorded in the region 4000-100 cm(-1). The structural and spectroscopic data of the molecule in the ground state have been calculated by using Hartree-Fock and density functional method (LSDA and MPW1PW91) with the 6-31G (d, p) and 6-311G (d, p) basis sets. The vibrational frequencies are calculated and scaled values have been compared with the experimental FT-IR and FT-Raman spectra. The observed and calculated frequencies are found in good agreement. The DFT-LSDA/6-311G (d, p) calculations have been found are more reliable than the ab initio HF/6-31G (d, p) calculations for the vibrational study of 3-BA. The optimized geometric parameters (bond lengths and bond angles) are compared with experimental values of the molecule. The alteration of vibrational bands due to the substitutions in the base molecule is also investigated from their characteristic region of linked spectrum.

Ab Initio Study of Hydrogen-Bond Formation between Cyclic Ethers and Selected Amino Acid Side Chains

Binding energies for hydrogen-bonded complexes of six cyclic ethers with five hydrogen-bond donor molecules that mimic selected amino acid side chains have been calculated at the MP2/6-31G*, MP2/6-31+G*, MP2/6-311++G**(single point), and MP2/aug-cc-pvtz levels, using geometries obtained with or without counterpoise corrections throughout the geometry optimization. The calculated basis set superposition error (BSSE) amounts to 10-20% and 5-10% of the uncorrected binding energies for the neutral and ionic species, respectively, at the MP2/aug-cc-pvtz level. The authors conclude that the O...H distances in the hydrogen bonds and binding energies for the studied systems may be determined with uncertainties of up to 0.08 A and 1-2 kcal/mol, respectively, in comparison with the MP2/aug-cc-pvtz values at a reasonable computational cost by performing standard geometry optimization at the MP2/6-31+G* level. Hydrogen-bond formation energies are more negative for cyclic ethers compared to their counterparts with a C=C double bond in the ring next to the oxygen atom. The less negative hydrogen-bonding energy and the increased O...H separation have been attributed to the reduced basicity of the ether oxygen when the lone pairs can enter conjugation with the pi-electrons of the Calpha=Cbeta double bond. The present study is the first step toward the development of an affordable computational level for estimating the binding energies of natural product, fused ring ether systems to the human estrogen receptor.

Molecular dynamics simulations of a guaiacyl beta-O-4 lignin model compound: examination of intramolecular hydrogen bonding and conformational flexibility

The dynamical conformational behavior of a guaiacyl beta-O-4 lignin model compound has been investigated by molecular simulations. The potential energy surface of the molecule in vacuum has been examined by means of an adiabatic map, showing a large accessible conformational space with multiple energy minima separated by low barriers. Molecular dynamics simulations have been performed in vacuum and with explicit solvent molecules for 10 and 2.1 ns, respectively. Molecular dynamics trajectories recorded in vacuum have shown the molecule to be flexible and to visit a large number of conformations. Many intramolecular H-bonds have been observed, existing for more than 90% of the total simulation time. The presence of explicit solvent molecules induces a significant broadening of some regions of the accessible conformational space and also largely reduces the statistical significance of intramolecular H-bonding. Intramolecular H-bonds observed in vacuum do not persist significantly and are preferentially exchanged with intermolecular H-bonds to the surrounding solvent molecules. The theoretical results are in good agreement with experimental NMR data that do not support the existence of strong and persistent intramolecular H-bonds in solution but instead indicate that H-bonds to solvent predominate. Finally, both molecular modeling and NMR approaches predict the guaiacyl beta-O-4 structure to be flexible and indicate that intramolecular H-bonds are not strong and persistent enough to confer rigidity to the molecule in solution.