Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “hydrogen bonds”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Energy localization efficiency in 1,3,5-trinitro-2,4,6-triaminobenzene pore collapse mechanisms

Atomistic and continuum scale modeling efforts have shown that the shock-induced collapse of porosity can occur via a wide range of mechanisms dependent on pore morphology, the shockwave pressure, and material properties. The mechanisms that occur under weaker shocks tend to be more efficient at localizing thermal energy but do not result in high, absolute temperatures or spatially large localizations compared to mechanisms found under strong shock conditions. However, the energetic material 1,3,5-trinitro-2,4,6-triaminobenzene (TATB) undergoes a wide range of collapse mechanisms that are not typical of similar materials, leaving the collapse mechanisms and the resultant energy localization from the collapse, i.e., hotspots, relatively uncharacterized. Therefore, we present the pore collapse simulations of cylindrical pores in TATB for a wide range of pore sizes and shock strengths that trigger viscoplastic collapses that occur almost entirely perpendicular to the shock direction for weak shocks and hydrodynamic-like collapses for strong shocks that do not break the strong hydrogen bonds of the TATB basal planes. The resulting hotspot temperature fields from these mechanisms follow trends that differ considerably from other energetic materials; hence, we compare them under normalized temperature values to assess the relative efficiency of each mechanism to localize energy. The local intra-molecular strain energy of the hotspots is also assessed to better understand the physical mechanisms behind the phenomena that lead to a latent potential energy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Electronic structure theory on modeling short-range noncovalent interactions between amino acids

While short-range noncovalent interactions (NCIs) are proving to be of importance in many chemical and biological systems, these atypical bindings happen within the so-called van der Waals envelope and pose an enormous challenge for current computational methods. We introduce SNCIAA, a database of 723 benchmark interaction energies of short-range noncovalent interactions between neutral/charged amino acids originated from protein x-ray crystal structures at the “gold standard” coupled-cluster with singles, doubles, and perturbative triples/complete basis set [CCSD(T)/CBS] level of theory with a mean absolute binding uncertainty less than 0.1 kcal/mol. Subsequently, a systematic assessment of commonly used computational methods, such as the second-order Møller−Plesset theory (MP2), density functional theory (DFT), symmetry-adapted perturbation theory (SAPT), composite electronic-structure methods, semiempirical approaches, and the physical-based potentials with machine learning (IPML) on SNCIAA is carried out. It is shown that the inclusion of dispersion corrections is essential even though these dimers are dominated by electrostatics, such as hydrogen bondings and salt bridges. Overall, MP2, ωB97M-V, and B3LYP+D4 turned out to be the most reliable methods for the description of short-range NCIs even in strongly attractive/repulsive complexes. SAPT is also recommended in describing short-range NCIs only if the δMP2 correction has been included. The good performance of IPML for dimers at close-equilibrium and long-range conditions is not transferable to the short-range. We expect that SNCIAA will assist the development/improvement/validation of computational methods, such as DFT, force-fields, and ML models, in describing NCIs across entire potential energy surfaces (short-, intermediate-, and long-range NCIs) on the same footing.

Chemistry↗
DOE OSTI2024

Monovalent ion–graphene oxide interactions are controlled by carboxylic acid groups: Sum frequency generation spectroscopy studies

Graphene oxide (GO) is a two-dimensional, mechanically strong, and chemically tunable material for separations. Elucidating GO–ion–water interactions at the molecular scale is highly important for predictive understanding of separation systems. However, direct observations of the nanometer region by GO surfaces under operando conditions are not trivial. Therefore, thin films of GO at the air/water interface can be used as model systems. With this approach, we study the effects of alkali metal ions on water organization near graphene oxide films at the air/water interface using vibrational sum frequency generation (SFG) spectroscopy. Here, we also use an arachidic acid Langmuir monolayer as a benchmark for a pure carboxylic acid surface. Theoretical modeling of the concentration-dependent sum frequency signal from graphene oxide and arachidic acid surfaces reveals that the adsorption of monovalent ions is mainly controlled by the carboxylic acid groups on graphene oxide. An in-depth analysis of sum frequency spectra reveals at least three distinct water populations with different hydrogen bonding strengths. The origin of each population can be identified from concentration dependent variations of their SFG signal. Interestingly, an interfacial water structure seemed mostly insensitive to the character of the alkali cation, in contrast to similar studies conducted at the silica/water interface. However, we observed an ion-specific effect with lithium, whose strong hydration prevented direct interactions with the graphene oxide film.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Deprotonated sulfamic acid and its homodimers: Does sulfamic acid adopt zwitterion during cluster growth?

Here we present a joint experimental and computational study on the geometric and electronic structures of deprotonated sulfamic acid (SA) clusters [(SA) n –H] – (n = 1, 2) employing negative ion photoelectron spectroscopy and high-level ab initio calculations. The photoelectron spectra provide the vertical/adiabatic detachment energy (VDE/ADE) of the sulfamate anion (SM – ) H 2 N•SO 3 – at 4.85 ± 0.05 and 4.58 ± 0.08 eV, respectively, and the VDE and ADE of the SM – •SA dimer at 6.41 ± 0.05 and 5.87 ± 0.08 eV, respectively. The significantly increased electron binding energies of the dimer confirm the enhanced electronic stability upon the addition of one SA molecule. The CCSD(T)-predicted VDEs/ADEs agree excellently with the experimental data, confirming the identified structures as the most stable ones. Two types of dimer isomers possessing different hydrogen bonding (HB) motifs are identified, corresponding to SM – binding to a zwitterionic SA (SM – •SA z ) and a canonical SA (SM – •SA c ), respectively. Two N–H$\cdots$O HBs and one superior O–H$\cdots$O HB are formed in the lowest-lying SM – •SA c , while SM – •SA z has three moderate N–H$\cdots$O HBs, with the former being 4.71 kcal/mol more stable. Further theoretical analyses reveal that the binding strength advantage of SM – •SA c over SM – •SA z arises from its significant contributions of orbital interactions between fragments, illustrating that sulfamate strongly interacts with its parent SA acid and preferably chooses the canonical SA in the subsequent cluster formations. Given the prominent presence of SA, this study provides the first evidence that the canonical dimer model of sulfamic acid should exist as a superior configuration during cluster growth.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

An experimental and computational view of the photoionization of diol–water clusters

In the interstellar medium, diols and other prebiotic molecules adsorb onto icy mantles surrounding dust grains. Water in the ice may affect the reactivity and photoionization of these diols. Ethylene glycol (EG), 1,2-propylene glycol, and 1,3-propylene glycol clusters with water clusters were used as a proxy to study these interactions. The diol–water clusters were generated in a continuous supersonic molecular beam, photoionized by synchrotron-based vacuum ultraviolet light from the Advanced Light Source, and subsequently detected by reflectron time-of-flight mass spectrometry. The appearance energies for the detected clusters were determined from the mass spectra, collected at increasing photon energy. Clusters of both diol fragments and unfragmented diols with water were detected. The lowest energy geometry optimized conformers for the observed EG–water clusters and EG fragment–water clusters have been visualized using density functional theory (DFT), providing insight into hydrogen bonding networks and how these affect fragmentation and appearance energy. As the number of water molecules clustered around EG fragments (m/z 31 and 32) increased, the appearance energy for the cluster decreased, indicating a stabilization by water. This trend was supported by DFT calculations. Fragment clusters from 1,2-propylene glycol exhibited a similar trend, but with a smaller energy decrease, and no trend was observed from 1,3-propylene glycol. In conclusion, we discuss and suggest that the reactivity and photoionization of diols in the presence of water depend on the size of the diol, the location of the hydroxyl group, and the number of waters clustered around the diol.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Prediction of hydration energies of adsorbates at Pt(111) and liquid water interfaces using machine learning

Aqueous phase heterogeneous catalysis is important to various industrial processes, including biomass conversion, Fischer–Tropsch synthesis, and electrocatalysis. Accurate calculation of solvation thermodynamic properties is essential for modeling the performance of catalysts for these processes. Explicit solvation methods employing multiscale modeling, e.g., involving density functional theory and molecular dynamics have emerged for this purpose. Although accurate, these methods are computationally intensive. This study introduces machine learning (ML) models to predict solvation thermodynamics for adsorbates on a Pt(111) surface, aiming to enhance computational efficiency without compromising accuracy. In particular, ML models are developed using a combination of molecular descriptors and fingerprints and trained on previously published water–adsorbate interaction energies, energies of solvation, and free energies of solvation of adsorbates bound to Pt(111). These models achieve root mean square error values of 0.09 eV for interaction energies, 0.04 eV for energies of solvation, and 0.06 eV for free energies of solvation, demonstrating accuracy within the standard error of multiscale modeling. Feature importance analysis reveals that hydrogen bonding, van der Waals interactions, and solvent density, together with the properties of the adsorbate, are critical factors influencing solvation thermodynamics. Furthermore, these findings suggest that ML models can provide rapid and reliable predictions of solvation properties. This approach not only reduces computational costs but also offers insights into the solvation characteristics of adsorbates at Pt(111)–water interfaces.

Adsorption↗
DOE OSTI2025

Measurement of electrical characteristics of Ge-based diodes with thermal GeO2 films exposed to controlled humidity conditions

We investigated the impact of the interaction between gas-phase water molecules and thermally oxidized GeO2 on Ge on the performance of metal–oxide–semiconductor (MOS) structures. A vacuum-integrated setup was developed to form MOS diodes with a GeO2/Ge structure that had been exposed to controlled humidity conditions and to subsequently measure their electrical characteristics in situ. The capacitance–voltage (C–V) curves exhibited a significant negative shift, indicating the generation of positive charges at the GeO2/Ge interface when the GeO2 surface was exposed to humidity levels above approximately 1%. According to a previous study using electron spectroscopy, this threshold corresponds to the humidity level at which molecular water begins to grow on a GeO2/Ge structure. It is likely that gas-phase water molecules infiltrating the GeO2 film bind to local OH sites via hydrogen bonding, leading to the formation of positive fixed charges at the GeO2/Ge interface.

Sano, Shuto↗
DOE OSTI2025

Microscopic insights into the solvation of polyethylene glycol chains in water: A machine learning potential approach

Polyethylene glycol (PEG) is a structurally simple, nontoxic, and water-soluble polymer widely utilized in medical and pharmaceutical applications. Notably, when a PEG chain is immersed in water, the surrounding water molecules play a key role in driving conformational changes of this macromolecule. In this study, we explore the solvation behavior of PEG under mechanical strain using molecular dynamics simulations, with an interatomic potential obtained from machine learning. Our focus is on the transition from the favored coil-like conformation to an extended one under external force. Through analyses of radial distribution functions, hydrogen bonding, and solvation dynamics, we uncover how mechanical stretching influences the local hydration environment. Furthermore, we disentangle the enthalpic and entropic contributions to the conformational stability of PEG in water. Surprisingly, our neural network potential model identifies dewetting of PEG C-atoms, and not water H-bonding with PEG O-atoms, as the main enthalpic driving force for the coiling of PEG in water.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Machine learning-accelerated path integral molecular dynamics simulations of reactive organic electrolytes

Hydrogen bonded electrolytes that exhibit accelerated proton transport via sequential reactive hops have drawn interest for their promise in clean energy applications. Molecular dynamics simulations of these electrolytes offer the opportunity to uncover microscopic mechanistic details that could be used to design and tune the properties of candidate electrolyte technologies. However, accurately modeling the proton transfer reactions and transport properties that give rise to high charge conductivites in these electrolytes proves computationally challenging because of the need to perform lengthy condensed phase simulations, treating both the electronic and nuclear degrees of freedom quantum mechanically. In this paper, we demonstrate that such a modeling task can be efficiently achieved with the use of density functional theory (DFT)-trained machine learning potentials (MLP) to accelerate path integral molecular dynamics (PIMD) simulations. We highlight the practical utility of this approach by using it to benchmark how closely PIMD simulations employing different DFT exchange–correlation functionals reproduce the composition-dependent densities, diffusion coefficients, and electrical conductivities of mixtures consisting of imidazole and levulinic acid. Even with the speedup afforded by our MLPs, PIMD simulations remain quite expensive. Furthermore, in order to render PIMD more computationally tractable, we introduce and benchmark the accuracy of a ring polymer contraction approach that leverages a computationally efficient short-range MLP to accelerate our PIMD simulations by an additional factor of four.

Chemical bonding↗
DOE OSTI2019

Structure of the cytochrome aa 3 -600 heme-copper menaquinol oxidase bound to inhibitor HQNO shows TM0 is part of the quinol binding site

Virtually all proton-pumping terminal respiratory oxygen reductases are members of the heme-copper oxidoreductase superfamily. Most of these enzymes use reduced cytochrome c as a source of electrons, but a group of enzymes have evolved to directly oxidize membrane-bound quinols, usually menaquinol or ubiquinol. All of the quinol oxidases have an additional transmembrane helix (TM0) in subunit I that is not present in the related cytochrome c oxidases. The current work reports the 3.6-Å-resolution X-ray structure of the cytochrome aa 3 -600 menaquinol oxidase from Bacillus subtilis containing 1 equivalent of menaquinone. The structure shows that TM0 forms part of a cleft to accommodate the menaquinol-7 substrate. Crystals which have been soaked with the quinol-analog inhibitor HQNO (N-oxo-2-heptyl-4-hydroxyquinoline) or 3-iodo-HQNO reveal a single binding site where the inhibitor forms hydrogen bonds to amino acid residues shown previously by spectroscopic methods to interact with the semiquinone state of menaquinone, a catalytic intermediate.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

A cocrystal structure of dengue capsid protein in complex of inhibitor

Dengue virus (DENV) was listed as a top 10 public health threat by the World Health Organization in 2019. No clinically approved anti-DENV drug is currently available. Here we report the high-resolution cocrystal structure (1.5 Å) of the DENV-2 capsid protein in complex with an inhibitor that potently suppresses DENV-2 but not other DENV serotypes. The inhibitor induces a “kissing” interaction between two capsid dimers. The inhibitor-bound capsid tetramers are assembled inside virions, resulting in defective uncoating of nucleocapsid when infecting new cells. Resistant DENV-2 emerges through one mutation that abolishes hydrogen bonds in the capsid structure, leading to a loss of compound binding. Structure-based analysis has defined the amino acids responsible for the inhibitor’s inefficacy against other DENV serotypes. The results have uncovered an antiviral mechanism through inhibitor-induced tetramerization of the viral capsid and provided essential structural and functional knowledge for rational design of panserotype DENV capsid inhibitors.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Structural basis for differential recognition of phosphohistidine-containing peptides by 1-pHis and 3-pHis monoclonal antibodies

In 2015, monoclonal antibodies (mAbs) that selectively recognize the 1-pHis or 3-pHis isoforms of phosphohistidine were developed by immunizing rabbits with degenerate Ala/Gly peptides containing the nonhydrolyzable phosphohistidine (pHis) analog- phosphotriazolylalanine (pTza). Here, we report structures of five rabbit mAbs bound to cognate pTza peptides: SC1-1 and SC50-3 that recognize 1-pHis, and their 3-pHis–specific counterparts, SC39-4, SC44-8, and SC56-2. These cocrystal structures provide insights into the binding modes of the pTza phosphate group that are distinct for the 1- and 3-pHis mAbs with the selectivity arising from specific contacts with the phosphate group and triazolyl ring. The mode of phosphate recognition in the 3-pHis mAbs recapitulates the Walker A motif, as present in kinases. The complementarity-determining regions (CDRs) of four of the Fabs interact with the peptide backbone rather than peptide side chains, thus conferring sequence independence, whereas SC44-8 shows a proclivity for binding a GpHAGA motif mediated by a sterically complementary CDRL3 loop. Specific hydrogen bonding with the triazolyl ring precludes recognition of pTyr and other phosphoamino acids by these mAbs. Kinetic binding experiments reveal that the affinity of pHis mAbs for pHis and pTza peptides is submicromolar. Bound pHis mAbs also shield the pHis peptides from rapid dephosphorylation. The epitope–paratope interactions illustrate how these anti-pHis antibodies are useful for a wide range of research techniques and this structural information can be utilized to improve the specificity and affinity of these antibodies toward a variety of pHis substrates to understand the role of histidine phosphorylation in healthy and diseased states.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Computational prediction of the effect of amino acid changes on the binding affinity between SARS-CoV-2 spike RBD and human ACE2

The association of the receptor binding domain (RBD) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike protein with human angiotensin-converting enzyme 2 (hACE2) represents the first required step for cellular entry. SARS-CoV-2 has continued to evolve with the emergence of several novel variants, and amino acid changes in the RBD have been implicated with increased fitness and potential for immune evasion. Reliably predicting the effect of amino acid changes on the ability of the RBD to interact more strongly with the hACE2 can help assess the implications for public health and the potential for spillover and adaptation into other animals. Here, we introduce a two-step framework that first relies on 48 independent 4-ns molecular dynamics (MD) trajectories of RBD-hACE2 variants to collect binding energy terms decomposed into Coulombic, covalent, van der Waals, lipophilic, generalized Born solvation, hydrogen bonding, π-π packing, and self-contact correction terms. The second step implements a neural network to classify and quantitatively predict binding affinity changes using the decomposed energy terms as descriptors. The computational base achieves a validation accuracy of 82.8% for classifying single–amino acid substitution variants of the RBD as worsening or improving binding affinity for hACE2 and a correlation coefficient of 0.73 between predicted and experimentally calculated changes in binding affinities. Both metrics are calculated using a fivefold cross-validation test. Our method thus sets up a framework for screening binding affinity changes caused by unknown single– and multiple–amino acid changes offering a valuable tool to predict host adaptation of SARS-CoV-2 variants toward tighter hACE2 binding.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Imaging active site chemistry and protonation states: NMR crystallography of the tryptophan synthase α-aminoacrylate intermediate

NMR-assisted crystallography-the integrated application of solid-state NMR, X-ray crystallography, and first-principles computational chemistry-holds significant promise for mechanistic enzymology: by providing atomic-resolution characterization of stable intermediates in enzyme active sites, including hydrogen atom locations and tautomeric equilibria, NMR crystallography offers insight into both structure and chemical dynamics. Here, this integrated approach is used to characterize the tryptophan synthase α-aminoacrylate intermediate, a defining species for pyridoxal-5'-phosphate-dependent enzymes that catalyze β-elimination and replacement reactions. For this intermediate, NMR-assisted crystallography is able to identify the protonation states of the ionizable sites on the cofactor, substrate, and catalytic side chains as well as the location and orientation of crystallographic waters within the active site. Most notable is the water molecule immediately adjacent to the substrate β-carbon, which serves as a hydrogen bond donor to the ε-amino group of the acid-base catalytic residue βLys87. From this analysis, a detailed three-dimensional picture of structure and reactivity emerges, highlighting the fate of the L-serine hydroxyl leaving group and the reaction pathway back to the preceding transition state. Reaction of the α-aminoacrylate intermediate with benzimidazole, an isostere of the natural substrate indole, shows benzimidazole bound in the active site and poised for, but unable to initiate, the subsequent bond formation step. When modeled into the benzimidazole position, indole is positioned with C3 in contact with the α-aminoacrylate C $\bar{β}$ and aligned for nucleophilic attack. Here, the chemically detailed, three-dimensional structure from NMR-assisted crystallography is key to understanding why benzimidazole does not react, while indole does.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Lithiating magneto-ionics in a rechargeable battery

Magneto-ionics, real-time ionic control of magnetism in solid-state materials, promise ultralow-power memory, computing, and ultralow-field sensor technologies. The real-time ion intercalation is also the key state-of-charge feature in rechargeable batteries. Here, we report that the reversible lithiation/delithiation in molecular magneto-ionic material, the cathode in a rechargeable lithium-ion battery, accurately monitors its real-time state of charge through a dynamic tunability of magnetic ordering. The electrochemical and magnetic studies confirm that the structural vacancy and hydrogen-bonding networks enable reversible lithiation and delithiation in the magnetic cathode. Coupling with microwave-excited spin wave at a low frequency (0.35 GHz) and a magnetic field of 100 Oe, we reveal a fast and reliable built-in magneto-ionic sensor monitoring state of charge in rechargeable batteries. The findings shown herein promise an integration of molecular magneto-ionic cathode and rechargeable batteries for real-time monitoring of state of charge.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Cryo-EM structures of light-harvesting 2 complexes from Rhodopseudomonas palustris reveal the molecular origin of absorption tuning

The genomes of some purple photosynthetic bacteria contain a multigene puc family encoding a series of α- and β-polypeptides that together form a heterogeneous antenna of light-harvesting 2 (LH2) complexes. To unravel this complexity, we generated four sets of puc deletion mutants in Rhodopseudomonas palustris , each encoding a single type of pucBA gene pair and enabling the purification of complexes designated as PucA-LH2, PucB-LH2, PucD-LH2, and PucE-LH2. The structures of all four purified LH2 complexes were determined by cryogenic electron microscopy (cryo-EM) at resolutions ranging from 2.7 to 3.6 Å. Uniquely, each of these complexes contains a hitherto unknown polypeptide, γ, that forms an extended undulating ribbon that lies in the plane of the membrane and that encloses six of the nine LH2 αβ-subunits. The γ-subunit, which is located near to the cytoplasmic side of the complex, breaks the C9 symmetry of the LH2 complex and binds six extra bacteriochlorophylls (BChls) that enhance the 800-nm absorption of each complex. The structures show that all four complexes have two complete rings of BChls, conferring absorption bands centered at 800 and 850 nm on the PucA-LH2, PucB-LH2, and PucE-LH2 complexes, but, unusually, the PucD-LH2 antenna has only a single strong near-infared (NIR) absorption peak at 803 nm. Comparison of the cryo-EM structures of these LH2 complexes reveals altered patterns of hydrogen bonds between LH2 αβ-side chains and the bacteriochlorin rings, further emphasizing the major role that H bonds play in spectral tuning of bacterial antenna complexes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

The birth and evolution of solvated electrons in the water

The photo-induced radiolysis of water is an elementary reaction in biology and chemistry, forming solvated electrons, OH radicals, and hydronium cations on fast time scales. Here, we use an optical-pump terahertz-probe spectroscopy setup to trigger the photoionization of water molecules with optical laser pulses at ~400 nm and then time-resolve the transient solvent response with broadband terahertz (THz) fields with a ~90 fs time resolution. We observe three distinct spectral responses. The first is a positive broadband mode that can be attributed to an initial diffuse, delocalized electron with a radius of (22 ± 1) Å, which is short lived (<200 fs) because the absorption is blue-shifting outside of the THz range. The second emerging spectroscopic signature with a lifetime of about 150 ps is attributed to an intermolecular mode associated with a mass rearrangement of solvent molecules due to charge separation of radicals and hydronium cations. After 0.2 ps, we observe a long-lasting THz signature with depleted intensity at 110 cm -1 that is well reproduced by ab initio molecular dynamics. We interpret this negative band at 110 cm -1 as the solvent cage characterized by a weakening of the hydrogen bond network in the first and second hydration shells of the cavity occupied by the localized electron.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Lignin with controlled structural properties by N-heterocycle-based deep eutectic solvent extraction

The complex and heterogeneous nature of the lignin macromolecule has presented a lasting barrier to its utilization. To achieve high lignin yield, the technical lignin extraction process usually severely modifies and condenses the native structure of lignin, which is a critical drawback for its utilization in conversion processes. In addition, there is no method capable of separating lignin from plant biomass with controlled structural properties. Here, we developed an N-heterocycle-based deep eutectic solvent formed between lactic acid and pyrazole (La–Py DES) with a binary hydrogen bonding functionality resulting in a high affinity toward lignin. In this study, up to 93.7% of lignin was extracted from wheat straw biomass at varying conditions from 90 °C to 145 °C. Through careful selection of treatment conditions as well as lactic acid to pyrazole ratios, lignin with controlled levels of ether linkage content, hydroxyl group content, and average molecular weight can be generated. Under mild extraction conditions (90 °C to 120 °C), light-colored native-like lignin can be produced with up to 80% yield, whereas ether linkage–free lignin with low polydispersity can be obtained at 145 °C. Overall, this study offers a new strategy for native lignin extraction and generating lignin with controlled structural properties.

59 BASIC BIOLOGICAL SCIENCES↗