SEARCH · Engineering Papers

Results for “Noncovalent interactions”

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.

SourceResearch connectionFrom yearThrough year

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

At least 109 records · Page 6

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Overall protein structure quality assessment using hydrogen-bonding parameters

Atomic model refinement at low resolution is often a challenging task. This is mostly because the experimental data are not sufficiently detailed to be described by atomic models. To make refinement practical and ensure that a refined atomic model is geometrically meaningful, additional information needs to be used such as restraints on Ramachandran plot distributions or residue side-chain rotameric states. However, using Ramachandran plots or rotameric states as refinement targets diminishes the validating power of these tools. Therefore, finding additional model-validation criteria that are not used or are difficult to use as refinement goals is desirable. Hydrogen bonds are one of the important noncovalent interactions that shape and maintain protein structure. These interactions can be characterized by a specific geometry of hydrogen donor and acceptor atoms. Systematic analysis of these geometries performed for quality-filtered high-resolution models of proteins from the Protein Data Bank shows that they have a distinct and a conserved distribution. Here, it is demonstrated how this information can be used for atomic model validation.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2024

Crystal Engineering of Hydrogen Bonding for Direct Air Capture of CO2: A Quantum Crystallography Perspective

Rising atmospheric CO2 levels demand efficient and sustainable carbon capture solutions. Direct air capture (DAC) via crystallizing hydrogen-bonded frameworks such as carbonate salts has emerged as a promising approach. This review explores the potential of crystal engineering, in tandem with advanced quantum crystallography techniques and computational modeling, to unlock the full potential of DAC materials. We examine the critical role of hydrogen bonding and other noncovalent interactions within a family of bis-guanidines that governs the formation of carbonate salts with high CO2 capture capacity and low regeneration energies for utilization. Quantum crystallography and charge density analysis prove instrumental in elucidating these interactions. A case study of a highly insoluble carbonate salt of a 2,6-pyridine-bis-(iminoguanidine) exemplifies the effectiveness of these approaches. However, challenges remain in the systematic and precise determination of hydrogen atom positions and atomic displacement parameters within DAC materials using quantum crystallography, and limitations persist in the accuracy of current energy estimation models for hydrogen bonding interactions. Future directions lie in exploring diverse functional groups, designing advanced hydrogen-bonded frameworks, and seamlessly integrating experimental and computational modeling with machine learning. This synergistic approach promises to propel the design and optimization of DAC materials, paving the way for a more sustainable future.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Improved Predictive Tools for Structural Properties of Metal–Organic Frameworks

The accurate determination of structural parameters is necessary to understand the electronic and magnetic properties of metal–organic frameworks (MOFs) and is a first step toward accurate calculations of electronic structure and function for separations and catalysis. Theoretical structural determination of metal-organic frameworks is particularly challenging because they involve ionic, covalent, and noncovalent interactions, which must be treated in a balanced fashion. Here, we apply a diverse group of local exchange-correlation functionals (PBE, PBEsol, PBE-D2, PBE-D3, vdW-DF2, SOGGA, MN15-L, revM06-L, SCAN, and revTPSS) to a broad test set of MOFs to seek the most accurate functionals to study various structural aspects of porous solids, in particular to study lattice constants, unit cell volume, two types of pore size characteristics, bond lengths, bond angles, and torsional angles). The recently developed meta functionals revM06-L and SCAN, without adding any molecular mechanics terms, are able to predict more accurate structures than previously recommended functionals, both those without molecular mechanics terms (PBE, PBEsol, vdW-DF2, and revTPSS) and those with them (PBE-D2 and PBE-D3). To provide a broader test, these two functionals are also tested for lattice constants and band gaps of unary, binary, and ternary semiconductors.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Non-Covalent Interactions and Helical Packing in Thiophene-Phenylene Copolymers: Tuning Solid-State Ordering and Charge Transport for Organic Field-Effect Transistors

In this study, we introduce two thiophene-phenylene-thiophene (TPT) polymers designed to leverage noncovalent intramolecular interactions to regulate main-chain conformation and enhance solid-state ordering. By incorporating unsubstituted thiophene (T) or bithiophene (2T) units, we reveal striking divergence in the thermal, morphological, and optoelectronic properties of the resulting films, facilitated by these noncovalent interactions. Using a combination of computational and experimental approaches, we show that annealing yields remarkably different polymer conformations and, consequently, charge transport properties. TPT-T undergoes a significant structural transformation, adopting a more planar backbone conformation and a highly crystalline, edge-on molecular orientation. In contrast, the introduction of a single additional thiophene unit in TPT-2T leads to a more isotropic molecular orientation with a slight preference for face-on alignment, resulting in a heterogeneous film structure that hinders charge transport despite achieving tighter molecular packing. Remarkably, despite being composed of achiral components, TPT-2T develops chirality upon annealing, indicating the formation of a helical conformation. Organic field-effect transistor measurements reveal that the well-ordered alignment in annealed TPT-T films results in higher charge carrier mobility and a narrower distribution of mobility values than in TPT-2T. These findings provide critical insights into the structure−property relationships of conjugated polymers, offering guidance for optimizing molecular design and processing strategies for highperformance organic electronic materials.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Systematic Evaluation of Counterpoise Correction in Density Functional Theory

A widespread belief persists that the Boys–Bernardi function counterpoise (CP) procedure “overcorrects” supramolecular interaction energies for the effects of basis-set superposition error. To the extent that this is true for correlated wave function methods, it is usually an artifact of low-quality basis sets. The question has not been considered systematically in the context of density functional theory, however, where basis-set convergence is generally less problematic. We present a systematic assessment of the CP procedure for a representative set of functionals and basis sets, considering both benchmark data sets of small dimers and larger supramolecular complexes. The latter include layered composite polymers with ~150 atoms and ligand–protein models with ~300 atoms. Provided that CP correction is used, we find that intermolecular interaction energies of nearly complete-basis quality can be obtained using only double-ζ basis sets. Furthermore, this is less expensive as compared to triple-ζ basis sets without CP correction. CP-corrected interaction energies are less sensitive to the presence of diffuse basis functions as compared to uncorrected energies, which is important because diffuse functions are expensive and often numerically problematic for large systems. Our results upend the conventional wisdom that CP “overcorrects” for basis-set incompleteness. In small basis sets, CP correction is mandatory in order to demonstrate that the results do not rest on error cancellation.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2025

Intermolecular Interactions in Direct Air Capture Materials: Insights from Charge Density Analysis

Direct air capture (DAC) materials enable the removal of CO 2 from the atmosphere, but improving their efficiency requires a detailed understanding of the intermolecular interactions that govern CO 2 sorption and release. Here, we present an experimental electron density study of methylglyoxal-bis(iminoguanidine) (MGBIG), a promising DAC material, using high-resolution X-ray and neutron diffraction data combined with quantum crystallographic analysis. This approach bridges theoretical and experimental data by quantifying electron density distributions and revealing how hydrogen bonds stabilize CO 2 -derived carbonate phases and may influence the desorption behavior. We identify distinct hydrogen-bonding environments in two crystalline carbonate phases: P1, a transient kinetic product, and P3, a thermodynamically stable phase. Multipolar refinement and electrostatic potential and multipole moment calculations precisely map electron density distributions, revealing key hydrogen bonds involved in CO 2 capture. Topological analysis of electron density highlights a cooperative hydrogen-bonding network in the thermodynamically favored P3 phase, where enhanced electron density delocalization and water-mediated interactions contribute to a more stable lattice. Energetic analyses confirm that stronger hydrogen bonding networks enhance the stability of P3 with a binding energy of −607.0 kJ/mol and greater lattice stability (−847.3 kJ/mol) compared to P1 (−302.5 and −571.0 kJ/mol, respectively). Electrostatic potential maps further illustrate polarization patterns that may influence the stability of the binding of CO 2 and release conditions. These findings establish a direct experimental framework for linking electron density distributions to intermolecular interactions in DAC materials, providing a rational design strategy for optimizing sorbents with improved CO 2 capture efficiency and reduced energy demands.

Electron density↗

DOE OSTI2021

Fluorescence Microscopic Investigations of Molecular Dynamics in Self–Assembled Nanostructures

Many analytical methods employ self-assembled nanostructured materials as chemical recognition media. Molecular permeation through these materials exhibits unique selectivity owing to nanoconfinement-induced enhancement of permeant–nanostructure interactions. This Personal Account introduces our efforts to investigate the detailed dynamics of single or a small number of molecules in nanostructured materials. Here, we developed new experimental and analysis approaches built upon laser-based fluorescence microscopy to measure the detailed translational and orientational dynamics of molecules diffusing in horizontally-oriented, cylindrical nanostructures, including surfactant micelles, silica mesopores, block copolymer microdomains, and bolaamphiphile-based organic nanotubes. Our studies clarified nanoscale details on the structural/chemical heterogeneity of the nanostructures, and their impacts on molecular mass transport dynamics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

The Role of Cooperative Interactions Among Surfaces, Solvents, and Reactive Intermediates on Catalysis at Liquid–Solid Interfaces (Final Report DE-SC0020224)

This project established quantitative links between inner‑sphere chemistry (active metal identity, coordination, and zeolite topology) and outer‑sphere organization (solvent identity, hydrogen‑bond networks, and pore condensation) that govern rates, activation barriers, and selectivities for alkene epoxidation and epoxide ring‑opening at solid–liquid and quasi‑liquid–solid interfaces. We deconvoluted contributions from covalent interactions at active sites and noncovalent, solvent‑mediated interactions within pores by pairing well‑defined metal substituted zeolites with controlled solvent environments. We then mapped those contributions onto measurable kinetics (ΔH‡, ΔS‡), adsorption thermodynamics (ITC), and in situ spectroscopy. The transferrable outcomes include a set of design rules that include the following understandings. First, tune silanol ((SiOH)x) density and pore topology to organize solvent networks that selectively stabilize transition states. Second, exploit activity‑coefficient‑normalized rates and adsorption– barrier correlations to diagnose when solvent reorganization rather than surface chemistry limits performance. Third, use partial pore condensation (e.g., acetonitrile, water but also generalizable to other solvents) to elicit liquid‑like stabilization effects even in nominally vapor‑phase reactors. Collectively, these results provide strategies to increase epoxidation rates, improve oxidant utilization (i.e., selectivities), and steer regioselectivity in zeolite‑based catalytic processes relevant to sustainable oxidation chemistry. These outcomes should be transferable to other classes of reactions that proceed in microporous materials and under confinement provided by organized solvents (e.g., electrochemical double layers).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Mixed Molecular and Dissociative Water Adsorption on Hydroxylated TiO 2 (110): An Infrared Spectroscopy and Ab Initio Molecular Dynamics Study

For this work, we have investigated the structure and dynamics of water (D 2 O) adsorbed on TiO 2 (110) for coverages between 0 and 1 monolayer with infrared reflection absorption spectroscopy (IRAS) and ab initio molecular dynamics (AIMD) simulations. For D 2 O coverages as low as 0.4 monolayers (ML) on a hydroxylated surface, IR spectra typical of hydrogen-bonded chains of water molecules are observed. However, for D 2 O coverages ≥ 0.3 ML, a sharp, high frequency peak is also observed in the p-polarized spectra that is red-shifted relative to the bridging hydroxyl peak. This new peak is not observed for water adsorbed on an oxidized surface. Based on the AIMD simulations and comparisons with previous IR spectra for TiO 2 nanoparticles, the new peak is assigned to terminal hydroxyl groups produced by dissociative adsorption of some of the water on TiO 2 (110). The simulations indicate that water dissociation is related to the presence of defect electrons in the system, but not due to direct interactions between adsorbed water and bridging hydroxyls.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Investigating the Variation in Solvation Interactions of Choline Chloride-Based Deep Eutectic Solvents Formed Using Different Hydrogen Bond Donors

Choline chloride-based deep eutectic solvents (DESs) have gained widespread recognition as green solvents in catalysis and separation science as their designer properties permit solvation properties to be modulated by the choice of the hydrogen bond donor (HBD) as well as its molar ratio. Despite being one of the most popular classes of DESs, very little is currently understood regarding how their solvation characteristics differ among themselves as well as from other classes of DESs. Previous studies have shown that the catalytic activity, extraction efficiency, and solubilities of natural compounds can be influenced by the type of choline salt and HBD, but no study has systematically related these parameters with its solvation properties. In this study, inverse gas chromatography was employed for the first time to study the individual solvation characteristics of a diverse range of choline chloride and acetate-based DESs composed of amides, diols, organic acids, and carbohydrates. Furthermore, hydrogen bond acidity was found to be a dominant interaction for all DESs, especially those composed of organic acid HBDs. DESs composed of isomers of butane diol and hexane diol as HBDs exhibited higher hydrogen bond basicity and dispersive-type interactions compared to those composed of urea, acetamide, and organic acids with the position of the hydroxyl functional group influencing their solvation interactions. Choline acetate-based DESs exhibited lower dipolarity and stronger hydrogen bond basicity and dispersive-type interactions compared to choline chloride DESs. The solvation models developed in this study were used to interpret and explain DES behavior in previously reported studies involving catalysis, organic synthesis, and extractions demonstrating that the measured solute–solvent interactions can serve to predict the performance of choline chloride-based DESs when used in various applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Enhancing Substrate–Metal Catalyst Affinity via Hydrogen Bonding: Pd(II)-Catalyzed β-C(sp 3 )–H Bromination of Free Carboxylic Acids

Achieving sufficient substrate-metal catalyst affinity is a fundamental challenge for developing synthetically useful C–H activation reactions of weakly coordinating native substrates. While hydrogen bonding has been harnessed to bias site selectivity in existing C(sp 2 )–H activation reactions, the potential for designing catalysts with hydrogen bond donors (HBDs) to enhance catalyst substrate affinity, and thereby facilitate otherwise unreactive C(sp 3 )–H activation, remains to be demonstrated. Herein, we report the discovery of a ligand scaffold containing a remote amide motif that can form a favorable meta-macrocyclic hydrogen bonding interaction with the aliphatic acid substrate. The utility of this ligand scaffold is demonstrated through the development of an unprecedented C(sp 3 )–H bromination of α-tertiary and α-quaternary free carboxylic acids, which proceeds in exceedingly high mono-selectivity. Furthermore, the geometric relationship of the NHAc hydrogen bond donor and the coordinating quinoline ligand is crucial for forming the meta-macrocyclophane-like hydrogen bonding interaction, providing a guideline for the future design of catalysts employing secondary interactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Molecular Dynamics Investigation of Nanoscale Hydrophobicity of Polymer Surfaces: What Makes Water Wet?

The wettability of a polymer surface–related to its hydrophobicity or tendency to repel water–can be crucial for determining its utility, such as for a coating or a purification membrane. While wettability is commonly associated with the macroscopic measurement of a contact angle between surface, water, and air, the molecular physics that underlie these macroscopic observations are not fully known, and anticipating relative behavior of different polymers is challenging. To address this gap in molecular-level understanding, we use molecular dynamics simulations to investigate and contrast interactions of water with six chemically distinct polymers: polytetrafluoroethylene, polyethylene, polyvinyl chloride, poly(methyl methacrylate), Nylon 66, and polyvinyl alcohol. We show that several prospective quantitative metrics for hydrophobicity agree well with experimental contact angles. Moreover, the behavior of water in proximity to these polymer surfaces can be distinguished with analysis of interfacial water dynamics, extent of hydrogen bonding, and molecular orientation–even when macroscopic measures of hydrophobicity are similar. The predominant factor dictating wettability is found to be the extent of hydrogen bonding between polymer and water, but the precise manifestation of hydrogen bonding and its impact on surface-water structure varies. In the absence of hydrogen bonding, other molecular interactions and polymer mechanics control hydrophobic ordering. Furthermore, these results provide new insights into how polymer chemistry specifically impacts water-polymer interactions and translates to surface hydrophobicity. Such factors may facilitate the design or processing of polymer surfaces to achieve targeted wetting behavior, and presented analyses can be useful in studying the interfacial physics of other systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Insights into the Solvation Characteristics of Zwitterionic Deep Eutectic Solvents Using Multiple Polarity Scales

Deep eutectic solvents (DESs) formed by mixing zwitterionic hydrogen bond acceptors (HBAs) with diol- and organic acid-based hydrogen bond donors (HBDs) are popular sustainable solvents. In addition to possessing an amphiphilic character, zwitterionic DESs offer unique solvation properties, unlike those composed of halide-containing HBAs. While all other popular classes of DESs have been investigated, no study has attempted to characterize the solvation characteristics of zwitterionic DESs. In this study, the solvation properties of betaine [Bet]-, sulfobetaine-, and l-carnitine [l-Car]-based DESs were thoroughly investigated for the first time using solvatochromic dye methods and inverse gas chromatography (IGC). Herein, the Nile red polarity scale revealed that DESs composed of organic acids were generally more polar than diol-based DESs. Kamlet–Taft solvent parameters revealed that glycerol-based DESs offered exceptionally high hydrogen bond donating ability compared to those composed of 1,4-butane diol and triethylene glycol. The Abraham solvation parameter model revealed hydrogen bond basicity to be the most dominant interaction, and it was generally higher for diol-based HBDs compared to organic acid-based HBDs. Additionally, sulfobetaine HBAs possessing longer alkyl chains generally offered higher dispersive-type interactions and hydrogen bond basicity compared to [Bet] and [l-Car] salts. This is the first study to combine and compare results from solvatochromic dye methods and IGC in an attempt to comprehensively explain and predict the solvation characteristics of zwitterionic DESs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Water Structure and Dynamics Near the Surfaces of Silicalite-1

Zeolites are crystalline microporous aluminosilicates that are commonly used as industrial sorbents and membranes. Solvent structuring and dynamics near zeolite crystal surfaces are thought to influence interfacial transport and molecular recognition processes, as well as fundamental aspects of their crystallization from solution. Here, in this work, we use molecular dynamics (MD) simulations to investigate the behavior of interfacial water near the exposed (010), (100), and (101) crystal faces of silicalite-1, one of the most widely studied zeolites. The MD simulations reveal that water’s translational and orientational order is strongly influenced by the distinct corrugations, pore apertures, and functional group distributions presented on each surface. Specifically, we observe two distinct hydration layers near each surface. Water molecules in the contact layer are hydrogen bonded to two or three exposed surface silanol groups. The relative populations of water molecules that are doubly and triply hydrogen bonded to the surface strongly depend on the distribution of the exposed silanols. The interactions with exposed surface silanols also influence water’s local orientational order and distribution across each surface. We also show that the structuring of the solvent near the different faces of silicalite-1 strongly impacts the interfacial dynamics. The translational and orientational relaxation dynamics of water are slowest in the contact layers and correlate with structural ordering near each face.

42 ENGINEERING↗

DOE OSTI2020

Effects of Non-Electrostatic Intermolecular Interactions on the Phase Behavior of pH-Sensitive Polyelectrolyte Complexes

Polyelectrolyte complexes (PECs) offer enormous material tunability and desirable functionalities, and consequently have found broad utility in biomedical and material industries. While poly(acrylic acid) (PAA) and poly( allylamine hydrochloride) (PAH) are a commonly used pairing, various aspects of the phase behavior of PAA-PAH complexes have not been sufficiently quantified. In this work, we present a comprehensive experimental study depicting the binodal phase boundaries for the PAA-PAH complexes prepared under acidic, neutral, and basic conditions using thermogravimetric analysis, turbidimetry, and optical microscopy. Under neutral and basic conditions, phase behaviors of the complexes were largely similar to one another and followed general expectations of PEC phase behavior, except for unusually high resistance to disruption of the complex with added salt. Stable complexes are observed up to 4 M NaCl concentrations. Under acidic conditions, strikingly different phase behaviors of the PAA-PAH complexes were observed. The polymer content in the complex phase increased initially, followed by an expected decrease as salt was added to the complexes. This behavior may result from a combination of associative phase separation of PAA and PAH chains, influenced by electrostatic interactions, and segregative phase separation, which can be ascribed to the influence of a combination of the hydrophobic interactions of the aliphatic polymer backbone and the interpolymer hydrogen bonding of unionized acrylic monomer units. Our systematic investigations over a range of pH detailing these discrepancies in the PAA-PAH phase behavior are expected to clarify the inconsistencies among the reports in the literature and provide the material design strategies for practical use of the PAA-PAH complexes and multilayer assemblies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Adjusting the Energy Profile for CH–O Interactions Leads to Improved Stability of RNA Stem-Loop Structures in MD Simulations

The role of ribonucleic acid (RNA) in biology continues to grow, but insight into important aspects of RNA behavior is lacking, such as dynamic structural ensembles in different environments, how flexibility is coupled to function, and how function might be modulated by small molecule binding. In the case of proteins, much progress in these areas has been made by complementing experiments with atomistic simulations, but RNA simulation methods and force fields are less mature. It remains challenging to generate stable RNA simulations, even for small systems where well-defined, thermostable structures have been established by experiments. Further many different aspects of RNA energetics have been adjusted in force fields, seeking improvements that are transferable across a variety of RNA structural motifs. In this work, the role of weak CH···O interactions is explored, which are ubiquitous in RNA structure but have received less attention in RNA force field development. By comparing data extracted from high-resolution RNA crystal structures to energy profiles from quantum mechanics and force field calculations, it is shown that CH···O interactions are overly repulsive in the widely used Amber RNA force fields. A simple, targeted adjustment of CH···O repulsion that leaves the remainder of the force field unchanged was developed. Then, the standard and modified force fields were tested using molecular dynamics (MD) simulations with explicit water and salt, amassing over 300 μs of data for multiple RNA systems containing important features such as the presence of loops, base stacking interactions as well as canonical and noncanonical base pairing. In this work and others, standard force fields lead to reproducible unfolding of the NMR-based structures. Including a targeted CH···O adjustment in an otherwise identical protocol dramatically improves the outcome, leading to stable simulations for all RNA systems tested.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Covalent Drug Binding in Live Cells Monitored by Mid-Infrared Quantum Cascade Laser Spectroscopy: Photoactive Yellow Protein as a Model System

The detection of drug-target interactions in live cells enables analysis of therapeutic compounds in a native cellular environment. Recent advances in spectroscopy and molecular biology have facilitated the development of genetically encoded vibrational probes like nitriles that can sensitively report on molecular interactions. Nitriles are powerful tools for measuring electrostatic environments within condensed media like proteins, but such measurements in live cells have been hindered by low signal-to-noise ratios. In this study, we design a spectrometer based on a double-beam quantum cascade laser (QCL)-based transmission infrared (IR) source with balanced detection that can significantly enhance sensitivity to nitrile vibrational probes embedded in proteins within cells compared to a conventional FTIR spectrometer. Here, using this approach, we detect small-molecule binding in Escherichia coli, with particular focus on the interaction between para-Coumaric acid (pCA) and nitrile-incorporated photoactive yellow protein (PYP). This system effectively serves as a model for investigating covalent drug binding in a cellular environment. Notably, we observe large spectral shifts of up to 15 cm –1 for nitriles embedded in PYP between the unbound and drug-bound states directly within bacteria, in agreement with observations for purified proteins. Such large spectral shifts are ascribed to the changes in the hydrogen-bonding environment around the local environment of nitriles, accurately modeled through high-level molecular dynamics simulations using the AMOEBA force field. Our findings underscore the QCL spectrometer’s ability to enhance sensitivity for monitoring drug–protein interactions, offering new opportunities for advanced methodologies in drug development and biochemical research.

chromophores↗

DOE OSTI2022

Data-Centric Development of Lignin Structure–Solubility Relationships in Deep Eutectic Solvents Using Molecular Simulations

Lignin is a natural source of aromatic chemicals with significant potential as an abundant, renewable feedstock for value-added products. Deep eutectic solvents (DES)–solvents composed of a hydrogen bond donor (HBD) and acceptor (HBA) in varying ratios–have emerged as a highly tunable class of solvents for lignin solubilization. However, the variety of possible DES compositions and limited molecular-scale understanding of lignin solubility makes solvent selection a challenge without laborious trial-and-error experimentation. To address these challenges, we use classical molecular dynamics (MD) simulations to study the interactions of lignin model compounds with various DES–water systems. Quantitative parameters (descriptors) were calculated by postprocessing the MD results and used to train a regression model that predicts experimentally determined solubilities of lignin model compounds. This approach revealed that the most important descriptors of solubility are the system temperature, solute hydrophilicity, and metrics quantifying hydrogen bonding. Maximizing the interactions between solute–HBD (hydrophobic group), water–HBD (hydrophilic group), and water–HBA molecules led to the highest model compound solubility. Our results support a hydrotropic mechanism in which extensive DES–water hydrogen bonding and favorable HBD interactions with the solute promote high solubility. We applied the regression model derived using model compounds to predict the solubility of representative lignin oligomers. The model predicted lignin oligomers’ solubilities in good agreement with experiments, indicating that the simulations of model compounds can be extended to predict the solubility of larger lignin compounds across a range of solvent compositions and temperatures. Furthermore, these findings provide new molecular-scale insight into lignin solubilization mechanisms and a new method for computationally screening potential solvent systems for lignin valorization.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗