Oxidative Cleavage of the β-O-4 Linkage of Lignin by Transition Metals: Catalytic Properties and the Performance of Density Functionals

Salicornia dolichostachya organosolv fractionation: towards establishing a halophyte biorefinery

Halophytes are a potential source of lignocellulosic material for biorefinery, as they can be grown in areas unsuitable for the cultivation of crops aimed at food production. To enable the viable use of halophytes in biorefineries, the present study investigated how different organosolv process parameters affected the fractionation of green pressed fibers of Salicornia dolichostachya. We produced pretreated solids characterized by up to 51.3% ± 1.7% cellulose, a significant increase from 25.6% ± 1.3% in untreated fibers. A delignification yield of as high as 60.7%, and hemicellulose removal of as high as 86.1% were also achieved in the current study. The obtained cellulose could be completely converted to glucose via enzymatic hydrolysis within 24 h. The lignin fractions obtained were of high purity, with sugar contamination of only 1.22% w/w and ashes below 1% w/w in most samples. Finally, up to 29.1% ± 0.4% hemicellulose was recovered as a separate product, whose proportion of oligomers to total sugars was 69.9% ± 3.0%. To the best of our knowledge, this is the first report in which Salicornia fibers are shown to be a suitable feedstock for organosolv biomass fractionation. These results expand the portfolio of biomass sources for biorefinery applications.

An organosolv method was developed for the fractionation of fibers of a halophyte plant in a biorefinery approach. Salicornia dolichostachya was used as raw material allowing the production of cellulose, hemicellulose, and lignin fractions.

Concentration effects on optical properties, DFT, crystal characterization and α-glucosidase activity studies: Novel Zn(II) complex

A novel Zn(II) complex of 6-ClpicH and picH was synthesized and its structure was determined by XRD technique. The detailed experimental optical susceptibility and band gap, refractive index, linear polarizability, optical and electrical conductivity parameters in various concentrations were investigated by means of the UV-Vis spectroscopic data. The optical band gap, refractive index (n), linear optical susceptibility (χ⁽¹⁾), third-order nonlinear optical susceptibility (χ⁽³⁾), second- and third-order nonlinear optical (β and γ) parameters were examined by using DFT/M06-L and ωB97XD/6-311++G(d,p) levels. The IC₅₀ value of Zn(II) complex against α-glucosidase was also obtained at 0.44 mM. The experimental band gap of the Zn(II) complex at 13, 33, 44 and 94 µM concentrations in ethanol were found to be 4.38, 4.37, 4.35 and 4.28 eV, respectively. The third-order NLO susceptibility χ⁽³⁾ parameter at 94 µM concentration corresponding to the photon energies of 4.6 and 5.7 eV in the UV-Vis region were observed at 206.6 × 10^-13 and 294.3 × 10^-13 esu, respectively. Besides, the theoretical χ⁽³⁾ values were obtained at 50.58 × 10^-13 and 20.37 × 10^-13 esu by using M06-L level. These results indicate that Zn(II) complex could be an effective third-order NLO candidate material. In brief, the detailed theoretical and experimental structural, spectral and optical properties of the Zn(II) complex were presented comparatively.

Isolating key reaction energetics and thermodynamic properties during hardwood model lignin pyrolysis

Computational studies on the pyrolysis of lignin using electronic structure methods have been largely limited to dimeric or trimeric models. In the current work we have modeled a lignin oligomer consisting of 10 syringyl units linked through 9 β-O-4' bonds. A lignin model of this size is potentially more representative of the polymer in angiosperms; therefore, we used this representative model to examine the behavior of hardwood lignin during the initial steps of pyrolysis. Using this oligomer, the present work aims to determine if and how the reaction enthalpies of bond cleavage vary with positions within the chain. To accomplish this, we utilized a composite method using molecular mechanics based conformational sampling and quantum mechanically based density functional theory (DFT) calculations. Our key results show marked differences in bond dissociation enthalpies (BDE) with the position. In addition, we calculated standard thermodynamic properties, including enthalpy of formation, heat capacity, entropy, and Gibbs free energy for a wide range of temperatures from 25 K to 1000 K. The prediction of these thermodynamic properties and the reaction enthalpies will benefit further computational studies and cross-validation with pyrolysis experiments. Overall, the results demonstrate the utility of a better understanding of lignin pyrolysis for its effective valorization.

DFT Functionals for Modeling of Polyethylene Chains Cross-Linked by Metal Atoms. DLPNO-CCSD(T) Benchmark Calculations

Density functional theory (DFT) functionals for calculations of binding energies (BEs) of the polyethylene (PE) chains cross-linked by selected metal atoms (M) are benchmarked against DLPNO-CCSD(T) and DLPNO-CCSD(T1) data. PEX-M-PEX complexes as compared with plain parallel PEX···PEX chains with X = 3-9 carbon atoms are model species characterized by a cooperative effect of covalent C-M-C bonds and interchain dispersion interactions. The accuracy of DLPNO-CC methods was assessed by a comparison of BEs with the canonical CCSD(T) results for small PE3-M-PE3 complexes. Functionals for PEX···PEX and closed-shell PEX-M-PEX complexes (M = Be, Mg, Zn) were benchmarked against DLPNO-CCSD(T) BEs; open-shell complexes (M = Li, Ag, Au) were benchmarked against the DLPNO-CCSD(T1) method with iterative triples. Three dispersion corrections were combined with 25 DFT functionals for calculations of BEs with respect to PEX-M and PEX fragments employing def2-TZVPP and def2-QZVPP basis sets. Accuracy to within 5% for the closed-shell PEX-M-PEX complexes was achieved with five functionals. Less accurate are functionals for the open-shell PEX-M-PEX complexes; only two functionals deviate by less than 15% from DLPNO-CCSD(T1). Particularly problematic were PEX-Li-PEX complexes. A reasonable overall performance across all complexes in terms of the mean absolute percentage error is found for the range-separated hybrid functionals ωB97X-D3 and CAM-B3LYP/D3(BJ)-ABC.

Understanding Plant Biomass via Computational Modeling

Plant biomass, especially wood, has been used for structural materials since ancient times. It is also showing great potential for new structural materials and it is the major feedstock for the emerging biorefineries for building a sustainable society. The plant cell wall is a hierarchical matrix of mainly cellulose, hemicellulose, and lignin. Herein, the structure, properties, and reactions of cellulose, lignin, and wood cell walls, studied using density functional theory (DFT) and molecular dynamics (MD), which are the widely used computational modeling approaches, are reviewed. Computational modeling, which has played a crucial role in understanding the structure and properties of plant biomass and its nanomaterials, may serve a leading role on developing new hierarchical materials from biomass in the future.

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.

A Genetic Algorithmic Approach to Determine the Structure of Li-Al Layered Double Hydroxides

Layered double hydroxides (LDH) demonstrate significant potential across a range of applications, including as catalysts, delivery vehicles for pharmaceuticals, environmental remediation, and supercapacitors. Explaining the mechanism of LDH action at the atomic scale in these and other applications is challenging, however, due to the difficulty in precisely defining the bulk and surface structure and chemical compositions. Here, we focus on the determination of the structure of lithium-aluminum (Li-Al) LDH, which has shown promise in the catalytic depolymerization of lignin, both directly as the catalyst and as a support for gold nanoparticles. While the relative positions of the Li and Al metals are generally well resolved by X-ray crystallography, it is the structures of the anionic layers, consisting of water and carbonate, that are less well established. Combinatorial analyses of all possible positions and rotations of the water and carbonate in the three-layered Li-AL LDH polytope reveals that the phase space is much too large to examine in any reasonable time frame in a one-by-one structure exploration. To overcome this limitation, we develop and deploy a genetic algorithm (GA) wherein fitness is determined by matching a calculated X-ray diffraction (XRD) pattern for a given structure to the known experimental XRD pattern. The GA approach results in structures of high fitness that portend the bulk Li-Al LDH structure. Importantly, the GA approach offers the potential to determine the structures of other LDH, and more generally layered materials, which are generally difficult to describe given the large chemical and structural space to be explored.

Modified Lennard-Jones potentials for nanoscale atoms

A classical 6-12 Lennard-Jones (LJ) equation has been widely used to model materials and is the potential of choice in studies when the focus is on fundamental issues. Here we report a systematic study comparing the pair interaction potentials within solid-state materials (i.e., [Co₆ Se₈ (PEt₃ )₆ ][C₆₀ ]₂ , [Cr₆ Te₈ (PEt₃ )₆ ][C₆₀ ]₂ , [Ni₉ Te₆ (PEt₃ )₈ ][C₆₀ ]) using density functional theory (DFT) calculations and LJ parametrization. Both classical (6-12 LJ) and modified LJ (mLJ) models were developed. In the mLJ approach, the exponents 6 and 12 are replaced by different integer number n and 2n, respectively, and an additional parameter (α) is introduced to describe intermolecular distance shift arising within the geometric centers' approach (instead of the shortest interatomic distance between particles). A general LJ approach reexamination reveals that in the case of nanoatoms, the attractive term decays with distance as the inverse fourth power, and the dominating at short distances repulsive term decays as the inverse eighth power. The modification of the LJ equation is even more prominent for interaction profiles, where intermolecular distance corresponds to separation between geometric centers of particles. In this approach, the attractive term decays with distance as the inverse 12th power, while the repulsive term decays rapidly (as the inverse 24th power). Thus, the mLJ models (e.g., 4-8 LJ) rather than the 6-12 classical ones seem to be a better choice for the description of binary interactions of nanoatoms. The developed mLJ models and electronic structure characteristics give an insight into the explanation of the unique physicochemical properties of superatomic-based solid-state materials.

Prediction of pKa s of Late Transition-Metal Hydrides via a QM/QM Approach

Three implicit solvation models, the conductor-like polarizable continuum model (C-PCM), the conductor-like screening model (COSMO), and universal implicit solvent model (SMD), combined with a hybrid two layer QM/QM approach (ONIOM), were utilized to calculate the pK_a values, using a direct thermodynamic scheme, of a set of Group 10 transition metal (TM) hydrides in acetonitrile. To obtain the optimal combination of quantum methods for ONIOM calculations with implicit solvation models, the influence of factors, such as the choice of density functional and basis set, the atomic radii used to build a cavity in the solvent, and the size of the model system in an ONIOM scheme, was examined. Additionally, the impact of Grimme's empirical dispersion correction and exact exchange was also investigated. The results were calibrated by experimental data. This investigation provides insight about effective models for the prediction of thermodynamic properties of TM-containing complexes with bulky ligands. © 2019 Wiley Periodicals, Inc.

Thermal behaviors and harmful volatile constituents released from asphalt components at high temperature

Asphalt binder releases lots of heat and harmful volatiles at high temperature. To further understand thermal behaviors, dynamic release and toxic constituents of emitted volatiles during the combustion of asphalt binder, such fractions as saturates, aromatics, resins and asphaltenes (SARA) were first prepared. Thermal behaviors, volatile constituents and combustion residue microstructures of SARA fractions are discussed. Results indicate that polymerization degree of asphalt binder is high and the content of polycyclic aromatic compounds is large. Combustion processes of resins and asphaltenes only show single-stage exothermic reactions, but other two fractions present obvious multi-stage combustion reactions. As the heating rate is raised, the incomplete combustion of SARA fractions is increased, and more volatiles are released. Main volatiles released from SARA fractions are inflammable, toxic, corrosive or explosive compounds, and such common volatiles as acetaldehyde and propane are released from each SARA fraction. More toxic volatiles are released at combustion stage I, but macromolecular volatiles are mainly released at stage II. Volatile release behaviors of saturates and aromatics are more obviously affected by the heating rate. Combustion residues show more intact morphologies from saturates to asphaltenes, and mainly contain C, O and S elements. Asphalt binder is hazardous material at high temperature.

Wavefunction-like Correlation Model for Use in Hybrid Density Functionals

We present Unsöld-W12 (UW12), an approximation to the correlation energy of molecules that is an explicit functional of the single-particle reduced-density matrix. The approximation resembles one part of modern explicitly correlated second-order Møller-Plesset (MP2) theory and is intended as an alternative to MP2 in double-hybrid exchange-correlation functionals. Orbital optimization with UW12 is straightforward, and the UW12 energy is evaluated without a double summation over unoccupied orbitals, leading to a faster basis-set convergence than is seen in double-hybrid functionals. We suggest a one-parameter hybrid exchange-correlation functional XCH-BLYP-UW12. XCH-BLYP-UW12 is similar to double-hybrid functionals, but contains UW12 correlation instead of MP2 correlation. We find that XCH-BLYP-UW12 is more accurate than the existing double-hybrid functional B2-PLYP for small-molecule main-group reaction barrier heights and has roughly the same accuracy as the existing hybrid functional B3LYP for atomization energies.

Molecular Products and Fundamentally Based Reaction Pathways in the Gas-Phase Pyrolysis of the Lignin Model Compound p-Coumaryl Alcohol

The fractional pyrolysis of lignin model compound para-coumaryl alcohol (p-CMA) containing a propanoid side chain and a phenolic OH group was studied using the System for Thermal Diagnostic Studies at temperatures from 200 to 900 °C, in order to gain mechanistic insight into the role of large substituents in high-lignin feedstocks pyrolysis. Phenol and its simple derivatives p-cresol, ethyl-, propenyl-, and propyl-phenols were found to be the major products predominantly formed at low pyrolysis temperatures (<500 °C). A cryogenic trapping technique was employed combined with EPR spectroscopy to identify the open-shell intermediates registered at pyrolysis temperatures above 500 °C. These were characterized as radical mixtures primarily consisting of oxygen-linked conjugated radicals. A comprehensive potential energy surface analysis of p-CMA and p-CMA + H atom systems was performed using various DFT protocols to examine the possible role of concerted molecular eliminations and free-radical mechanisms in the formation of major products. Other significant unimolecular concerted reactions along with formation and decomposition of primary radicals are also described and evaluated. The calculations suggest that a set of the chemically activated secondary radical channels is relevant to the low temperature product formation under fractional pyrolysis conditions.

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.

A look at the density functional theory zoo with the advanced GMTKN55 database for general main group thermochemistry, kinetics and noncovalent interactions

We present the updated and extended GMTKN55 benchmark database for more accurate and extensive energetic evaluation of density functionals and other electronic structure methods with detailed guidelines for method users.