skip to main content

Explore Research Products in the NSF-PAR

Title: LigninGraphs: lignin structure determination with multiscale graph modeling
Abstract

Lignin is an aromatic biopolymer found in ubiquitous sources of woody biomass. Designing and optimizing lignin valorization processes requires a fundamental understanding of lignin structures. Experimental characterization techniques, such as 2D-heteronuclear single quantum coherence (HSQC) nuclear magnetic resonance (NMR) spectra, could elucidate the global properties of the polymer molecules. Computer models could extend the resolution of experiments by representing structures at the molecular and atomistic scales. We introduce a graph-based multiscale modeling framework for lignin structure generation and visualization. The framework employs accelerated rejection-free polymerization and hierarchical Metropolis Monte Carlo optimization algorithms. We obtain structure libraries for various lignin feedstocks based on literature and new experimental NMR data for poplar wood, pinewood, and herbaceous lignin. The framework could guide researchers towards feasible lignin structures, efficient space exploration, and future kinetics modeling. Its software implementation in Python, LigninGraphs, is open-source and available on GitHub.

Graphical abstract

  more » « less
NSF-PAR ID:
10368553
Author(s) / Creator(s):
; ; ; ;
Publisher / Repository:
Springer Science + Business Media
Date Published:
Journal Name:
Journal of Cheminformatics
Volume:
14
Issue:
1
ISSN:
1758-2946
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ; ; ; ; ( , Energy & Fuels)
    Hongwei Wu (Ed.)
    Pyrolytic lignin is a fraction of pyrolysis oil that contains a wide range of phenolic compounds that can be used as intermediates to produce fuels and chemicals. However, the characteristics of the raw lignin structure make it difficult to establish a pyrolysis mechanism and determine pyrolytic lignin structures. This study proposes dimer, trimer, and tetramer structures based on their relative thermodynamic stability for a hardwood lignin model in pyrolysis. Different configurations of oligomers were evaluated by varying the positions of the guaiacyl (G) and syringyl (S) units and the bonds βO4 and β5 in the hardwood model lignin through electronic structure calculations. The homolytic cleavage of βO4 bonds is assumed to occur and generate two free radical fragments. These can stabilize by taking hydrogen radicals that may be in solution during the intermediate liquid (pathway 1) formation before the thermal ejection. An alternative pathway (pathway 2) could occur when the radicals use intramolecular hydrogen, turning themselves into stable products. Subsequently, a demethylation reaction can take place, thus generating a methane molecule and new oligomeric lignin-derived molecules. The most probable resulting structures were studied. We used FTIR and NMR spectra of selected model compounds to evaluate our calculation approach. Thermophysical properties were calculated using group contribution methods. The results give insights into the lignin oligomer structures and how these molecules are formed. They also provide helpful information for the design of pyrolysis oil separation and upgrading equipment. 
    more » « less
  2. ; ( , ChemPhysChem)
    Abstract

    Potential energy surface (PES) analyses at the SMD[MP2/6–311++G(d,p)] level and higher‐level energies up to MP4(fc,SDTQ) are reported for the fluorinated tertiary carbamateN‐ethyl‐N‐(2,2,2‐trifluoroethyl) methyl carbamate (VII) and its parent systemN,N‐dimethyl methyl carbamate (VI). Emphasis is placed on the analysis of the rotational barrier about the CN carbamate bond and its interplay with the hybridization of theN‐lone pair (NLP). All rotational transition state (TS) structures were found by computation of 1D relaxed rotational profiles but only 2D PES scans revealed the rotation‐inversion paths in a compelling fashion. We found four unique chiral minima ofVII, one pair each ofE‐andZ‐rotamers, and we determined theeightunique rotational TS structures associated with every possibleE/Z‐isomerization path. It is a significant finding that all TS structures featureN‐pyramidalization whereas the minima essentially contain sp2‐hybridized nitrogen. We will show that the TS stabilities are affected by the synergetic interplay between NLP/CO2repulsion minimization, NLP→σ*(CO) negative hyperconjugation, and two modes of intramolecular through‐space electrostatic stabilization. We demonstrate how Boltzmann statistics must be applied to determine the predicted experimental rotational barrier based on the energetics of all eight rotamerization pathways. The computed barrier forVIIis in complete agreement with the experimentally measured barrier of the very similar fluorinated carbamateN‐Boc‐N‐(2,2,2‐trifluoroethyl)‐4‐aminobutan‐1‐olII. NMR properties ofVIIwere calculated with a variety of density functional/basis set combinations and Boltzmann averaging over theE‐andZ‐rotamers at our best theoretical level results in good agreement with experimental chemical shifts δ(13C) andJ(13C,19F) coupling constants ofII(within 6 %).

     
    more » « less
  3. ; ; ; ; ; ; ( , Journal of the American Ceramic Society)
    Abstract

    Borates and borosilicates are potential candidates for the design and development of glass formulations with important industrial and technological applications. A major challenge that retards the pace of development of borate/borosilicate based glasses using predictive modeling is the lack of reliable computational models to predict the structure‐property relationships in these glasses over a wide compositional space. A major hindrance in this pursuit has been the complexity of boron‐oxygen bonding due to which it has been difficult to develop adequate B–O interatomic potentials. In this article, we have evaluated the performance of three B–O interatomic potential models recently developed by Bauchy et al [J.Non‐Cryst. Solids, 2018, 498, 294–304], Du et al [J. Am. Ceram. Soc.https://doi.org/10.1111/jace.16082] and Edèn et al [Phys. Chem. Chem. Phys., 2018, 20, 8192–8209] aiming to reproduce the short‐to‐medium range structures of sodium borosilicate glasses in the system 25 Na2OxB2O3(75 − x) SiO2(x = 0‐75 mol%). To evaluate the different force fields, we have computed at the density functional theory level the NMR parameters of11B,23Na, and29Si of the models generated with the three potentials and the simulated MAS NMR spectra compared with the experimental counterparts. It was observed that the rigid ionic models proposed by Bauchy and Du can both reliably reproduce the partitioning between BO3and BO4species of the investigated glasses, along with the local environment around sodium in the glass structure. However, they do not accurately reproduce the second coordination sphere of silicon ions and the Si–O–T (T = Si, B) and B‐O‐T distribution angles in the investigated compositional space which strongly affect the NMR parameters and final spectral shape. On the other hand, the core‐shell parameterization model proposed by Edén underestimates the fraction of BO4species of the glass with composition 25Na2O 18.4B2O356.6SiO2but can accurately reproduce the shape of the11B and29Si MAS‐NMR spectra of the glasses investigations due to the narrower B–O–T and Si‐O‐T bond angle distributions. Finally, the effect of the number of boron atoms (also distinguishing the BO3and BO4units) in the second coordination sphere of the network former cations on the NMR parameters have been evaluated.

     
    more » « less
  4. ; ; ; ; ; ; ; ( , ChemSusChem)
    Abstract

    This study aimed to use a persulfate together with transition metal ions as the reagent to effectively depolymerize lignin into monophenolic compounds under mild conditions (ambient pressure, temperature <100 °C). The Box‐Behnken experimental design in combination with the response surface methodology was applied to obtain optimized reaction conditions. The results showed that this reagent could depolymerize up to 99 % of lignin dimers to mainly veratraldehyde. This reaction also successfully depolymerized industrial lignins with a high yield of phenolic oils and monophenolic compounds. Quantum chemistry calculations using the density functional theory level indicated that the persulfate free radical attacks Cβto break the β‐O‐4 bond of lignin through a five‐membered ring mechanism. This mechanism using persulfate free radicals has a lower activation barrier than that using hydroxyl radicals. Gel permeation chromatography and 2D‐NMR spectroscopy demonstrated the effective cleavage of the β‐O‐4 bonds of lignin after depolymerization.

     
    more » « less
  5. ; ; ; ; ; ; ; ; ; ; et al ( , ChemSusChem)
    Abstract

    Integrating multidisciplinary research in plant genetic engineering and renewable deep eutectic solvents (DESs) can facilitate a sustainable and economic biorefinery. Herein, we leveraged a plant genetic engineering approach to specifically incorporate C6C1monomers into the lignin structure. By expressing the bacterialubiCgene in sorghum,p‐hydroxybenzoic acid (PB)‐rich lignin was incorporated into the plant cell wall while this monomer was completely absent in the lignin of the wild‐type (WT) biomass. A DES was synthesized with choline chloride (ChCl) and PB and applied to the pretreatment of the PB‐rich mutant biomass for a sustainable biorefinery. The release of fermentable sugars was significantly enhanced (∼190 % increase) compared to untreated biomass by the DES pretreatment. In particular, the glucose released from the pretreated mutant biomass was up to 12 % higher than that from the pretreated WT biomass. Lignin was effectively removed from the biomass with the preservation of more than half of the β‐Ο‐4 linkages without condensed aromatic structures. Hydrogenolysis of the fractionated lignin was conducted to demonstrate the potential of phenolic compound production. In addition, a simple hydrothermal treatment could selectively extract PB from the same engineered lignin, showing a possible circular biorefinery. These results suggest that the combination of PB‐based DES and engineered PB‐rich biomass is a promising strategy to achieve a sustainable closed‐loop biorefinery.

     
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email