Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “molecular”

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 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

The use of a digital micromirror array as a temporal gate and spatial-filtering device for laser-induced breakdown spectroscopy and laser ablation molecular isotopic spectrometry

A digital micromirror array (DMMA) is used in a simple and inexpensive approach for spatially resolved temporal gating of atomic emission detection in laser-induced breakdown spectroscopy (LIBS). Selected mirrors in the array are actuated at a time delayed from the primary laser pulse, permitting optical gating with a response time of 160 ns while using a conventional CCD detector. Detector gating is shown to decrease noise from short-lived spectroscopic background (both continuum and N II from air), improving signal-to-background ratios by 22-times for selected samples and emission lines. The utility of rapid temporal gating is also demonstrated for the optical detection of isotopes using the laser-ablation molecular isotope spectrometry (LAMIS) experiment. Additionally, optical temporal gating and spatial filtering are used to isolate regions of the laser-induced plasma, enhancing signal from molecular emission. The potential of the simple technique for use in field-portable LIBS and LAMIS experiments is explored.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

EZFF: Python library for multi-objective parameterization and uncertainty quantification of interatomic forcefields for molecular dynamics

Parameterization of interatomic forcefields is a necessary first step in performing molecular dynamics simulations. This is a non-trivial global optimization problem involving quantification of multiple empirical variables against one or more properties. We present EZFF, a lightweight Python library for parameterization of several types of interatomic forcefields implemented in several molecular dynamics engines against multiple objectives using genetic-algorithm-based global optimization methods. The EZFF scheme provides unique functionality such as the parameterization of hybrid forcefields composed of multiple forcefield interactions as well as built-in quantification of uncertainty in forcefield parameters and can be easily extended to other forcefield functional forms as well as MD engines.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Exploring the interfacial structure and crystallinity for direct growth of Mn 3 Sn(0001) on sapphire (0001) by molecular beam epitaxy

The Kagome antiferromagnet Mn 3 Sn has garnered a lot of attention due to the presence of interesting properties such as anomalous Hall effects and Nernst effects. Until now, few papers have been reported to grow using molecular beam epitaxy with a buffer layer or on other substrates. In this paper, we discuss the synthesis of crystalline Mn 3 Sn layers, prepared on Al 2 O 3 (0001) without a buffer layer using molecular beam epitaxy. The growth is monitored in-situ using reflection high energy electron diffraction and measured ex-situ using X-ray diffraction, Rutherford back-scattering, and cross-sectional scanning transmission electron microscopy. The samples were deposited at 524 ± 5 °C, with an Mn: Sn atomic flux ratio of 3.2:1 for 90 minutes. Orientation relationships between the Mn 3 Sn films and the sapphire substrates are determined from in-plane and out-of-plane measurements. Our analysis indicates that the resulting film is predominantly c−plane oriented. Lastly, the samples prepared in this way were found to be discontiguous, showing a 3-dimensional morphology. According to first-principles calculations, the Mn 3 Sn exhibits a displaced Kagome structure in the very first stages of growth, for 2ML and 4ML growth on Al 2 O 3 (0001). This result is corroborated by calculating the surface formation energies and explains the observed RHEED patterns.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Probing molecular vibrations by monochromated electron microscopy

Chemical bonds fundamentally determine molecular properties and are prevalently characterized by various spectroscopic means such as infrared and Raman spectroscopies. However, the spatial resolution of these conventional approaches is insufficient to reveal nanoscale features. Recently, monochromated electron energy-loss spectroscopy (EELS) in the transmission electron microscope achieved a groundbreaking energy resolution of a few millielectronvolts and enabled direct observation of molecular vibrational spectrum with unmatched spatial resolution. Vibrational EELS is widely applicable to both organic and inorganic matter in the solid state or liquid phase. In this study, we introduce recent advancements and key concepts of this method, compare with other spectroscopic techniques, and discuss future developments for potential applications in research fields centered on catalysts, polymers, and live cells.

Vibrational Spectroscopy↗
DOE OSTI2024

Effect of layer bending on montmorillonite hydration and structure from molecular simulation

Conceptual models of smectite hydration include planar (flat) clay layers that undergo stepwise expansion as successive monolayers of water molecules fill the interlayer regions. However, X-ray diffraction (XRD) studies indicate the presence of interstratified hydration states, suggesting non-uniform interlayer hydration in smectites. Additionally, recent theoretical studies have shown that clay layers can adopt bent configurations over nanometer-scale lateral dimensions with minimal effect on mechanical properties. Therefore, in this study we used molecular simulations to evaluate structural properties and water adsorption isotherms for montmorillonite models composed of bent clay layers in mixed hydration states. Results are compared with models consisting of planar clay layers with interstratified hydration states (e.g. 1W–2W). The small degree of bending in these models (up to 1.5 Å of vertical displacement over a 1.3 nm lateral dimension) had little or no effect on bond lengths and angle distributions within the clay layers. Except for models that included dry states, porosities and simulated water adsorption isotherms were nearly identical for bent or flat clay layers with the same averaged layer spacing. Similar agreement was seen with Na- and Ca-exchanged clays. In conclusion, while the small bent models did not retain their configurations during unconstrained molecular dynamics simulation with flexible clay layers, we show that bent structures are stable at much larger length scales by simulating a 41.6×7.1 nm 2 system that included dehydrated and hydrated regions in the same interlayer.

58 GEOSCIENCES↗
DOE OSTI2023

Sarcocystis cruzi (Hasselmann, 1923) Wenyon, 1926: redescription, molecular characterization and deposition of life cycle stages specimens in the Smithsonian Museum

Currently, 7 named Sarcocystis species infect cattle: Sarcocystis hirsuta, S. cruzi, S. hominis, S. bovifelis, S. heydorni, S. bovini and S. rommeli; other, unnamed species also infect cattle. Of these parasites of cattle, a complete life cycle description is known only for S. cruzi, the most pathogenic species in cattle. The life cycle of S. cruzi was completed experimentally in 1982, before related parasite species were structurally characterized, and before the advent of molecular diagnostics; to our knowledge, no archived frozen tissues from the cattle employed in the original descriptions remain for DNA characterization. Here, we isolated DNA from a paraffin-embedded kidney of a calf experimentally infected with S. cruzi in 1980; we then sequenced portions of 18S rRNA, 28S rRNA, COX1 and Acetyl CoA genes and verified that each shares 99–100% similarity to other available isolates attributed to S. cruzi from naturally infected cattle. We also reevaluated histological sections of tissues of calves experimentally infected with S. cruzi in the original description, exploiting improvements in photographic technology to render clearer morphological detail. Finally, we reviewed all available studies of the life cycle of S. cruzi, noting that S. cruzi was transmitted between bison (Bison bison) and cattle (Bos taurus) and that the strain of parasite derived from bison appeared more pathogenic than the cattle strain. Based on these newfound molecular, morphological and physiological data, we thereby redescribed S. cruzi and deposited reference material in the Smithsonian Museum for posterity.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

N -Arylimide Molecular Balances: A Comprehensive Platform for Studying Aromatic Interactions in Solution

Noncovalent interactions of aromatic surfaces play a key role in many biological processes and in determining the properties and utility of synthetic materials, sensors, and catalysts. However, the study of aromatic interactions has been challenging because these interactions are usually very weak and their trends are modulated by many factors such as structural, electronic, steric, and solvent effects. Recently, N-arylimide molecular balances have emerged as highly versatile and effective platforms for studying aromatic interactions in solution. These molecular balances can accurately measure weak noncovalent interactions in solution via their influence on the folded–unfolded conformational equilibrium. The structure (i.e., size, shape, π-conjugation, and substitution) and nature (i.e., element, charge, and polarity) of the π-surfaces and interacting groups can be readily varied, enabling the study of a wide range of aromatic interactions. These include aromatic stacking, heterocyclic aromatic stacking, and alkyl–π, chalcogen–π, silver–π, halogen–π, substituent–π, and solvent–π interactions. The ability to measure a diverse array of aromatic interactions within a single model system provides a unique perspective and insights as the interaction energies, stability trends, and solvent effects for different types of interactions can be directly compared. Some broad conclusions that have emerged from this comprehensive analysis include: (1) The strongest aromatic interactions involve groups with positive charges such as pyridinium and metal ions which interact with the electrostatically negative π-face of the aromatic surface via cation–π or metal–π interactions. Attractive electrostatic interactions can also form between aromatic surfaces and groups with partial positive charges. (2) Electrostatic interactions involving aromatic surfaces can be switched from repulsive to attractive using electron-withdrawing substituents or heterocycles. These electrostatic trends appear to span many types of aromatic interactions involving a polar group interacting with a π-surface such as halogen–π, chalcogen–π, and carbonyl–π. (3) Nonpolar groups form weak but measurable stabilizing interactions with aromatic surfaces in organic solvents due to favorable dispersion and/or solvophobic effects. Furthermore, a good predictor of the interaction strength is provided by the change in solvent-accessible surface area. (4) Solvent effects modulate the aromatic interactions in the forms of solvophobic effects and competitive solvation, which can be modeled using solvent cohesion density and specific solvent–solute interactions.

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

Molecular Characterization of Water-Soluble Brown Carbon Chromophores in Snowpack from Northern Xinjiang, China

This study reports molecular-level characterization of brown carbon (BrC) attributed to water-soluble organic carbon in six snowpack samples collected from northern Xinjiang, China. Here, the molecular composition and light-absorbing properties of BrC chromophores were unraveled by application of high-performance liquid chromatography (HPLC) coupled to a photodiode array (PDA) detector and high-resolution mass spectrometry. The chromophores were classified into five major types, that is, (1) phenolic/lignin-derived compounds, (2) flavonoids, (3) nitroaromatics, (4) oxygenated aromatics, and (5) other chromophores. Identified chromophores account for ~23–64% of the total light absorption measured by the PDA detector in the wavelength range of 300–370 nm. In the representative samples from urban and remote areas, oxygenated aromatics and nitroaromatics dominate the absorption in the wavelengths below and above 320 nm, respectively. The highly polluted urban sample shows the most complex HPLC-PDA chromatogram, and more other chromophores contribute to the bulk absorption. Phenolic/ lignin-derived compounds are the most light-absorbing species in the soil-influenced sample. Chromophores in two remote samples exhibit ultraviolet–visible features distinct from other samples, which are attributed to flavonoids. Identification of individual chromophores and quantitative analysis of their optical properties are helpful for elucidating the roles of BrC in snow radiative balance and photochemistry.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Molecular Imaging Reveals Two Distinct Mixing States of PM 2.5 Particles Sampled in a Typical Beijing Winter Pollution Case

Mixing states of aerosol particles are crucial for understanding the role of aerosols in influencing air quality and climate. However, a fundamental understanding of the complex mixing states is still lacking because most traditional analysis techniques only reveal bulk chemical and physical properties with limited surface and 3-D information. In this research, 3-D molecular imaging enabled by ToF-SIMS was used to elucidate the mixing states of PM 2.5 samples obtained from a typical Beijing winter haze event. In light pollution cases, a thin organic layer covers separated inorganic particles, while in serious pollution cases, ion exchange and an organic-inorganic mixing surface on large-area particles were observed. Finally, the new results provide key 3-D molecular information of mixing states, which is highly desirable for reducing uncertainty and bias in representing aerosol-cloud interactions in current Earth System Models, and improving the understanding of aerosols on air quality and human health.

3-D molecular imaging↗
DOE OSTI2020

Automated Identification of Molecular Crystals’ Packing Motifs

Packing motifs—patterns in how molecules orient relative to one another in a crystal structure—are an important concept in many subdisciplines of materials science because of correlations observed between specific packing motifs and properties of interest. That said, packing motif data sets have remained small and noisy due to intensive manual labeling processes and insufficient labeling schemes. The most prominent labeling algorithms calculate relative interplanar angles of nearest neighbor molecules to determine the packing motif of a molecular crystal, but this simple approach can fail when neighbors are naively sampled isotropically around the crystal structure. To remedy this issue, here we propose an optimization algorithm, which rotates the molecular crystal structure to find representative molecules that inform the packing motif. We package this algorithm into an automated framework—Autopack—which both optimally rotates the crystal structure and labels the packing motif based on the appropriate neighboring molecules. In this work, we detail the Autopack framework and its performance, which shows improvements compared to previous state-of-the-art labeling methods, providing the first quantitative point of comparison for packing motif labeling algorithms. Furthermore, using Autopack (available at https://ipo.llnl.gov/technologies/software/autopack), we perform the first large-scale study of potential relationships between chemicals’ compositions and packing motifs, which shows that these relationships are more complex than previously hypothesized from studies that used only tens of polycyclic aromatic hydrocarbon molecules. Autopack’s capabilities help pose next steps for crystal engineering research focusing not only on a molecule’s adoption of a specific packing motif but also on new structure–property relationships.

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

First Principles Nonadiabatic Excited-State Molecular Dynamics in NWChem

Computational simulation of non-adiabatic molecular dynamics is an indispensable tool for understanding complex photoinduced processes such as internal conversion, energy transfer, charge separation, and spatial localization of excitons, to name a few. Here, we report an implementation of the fewest-switches surface hopping algorithm in the NWChem computational chemistry program. Here the surface hopping method is combined with linear-response time-dependent density functional theory calculations of adiabatic excited state potential energy surfaces. To treat quantum transitions between arbitrary electronic Born{Oppenheimer states, we have implemented both numerical and analytical differentiation schemes for derivative non-adiabatic couplings. A numerical approach for the time-derivative non-adiabatic couplings together with an analytical method for calculating non-adiabatic coupling vectors is an efficient combination for surface hopping approaches. Additionally, electronic decoherence schemes and a state reassigned unavoided crossings algorithm are also implemented to improve the accuracy of the simulated dynamics and to handle trivial unavoided crossings. We apply our code to study the ultrafast decay of photoexcited benzene, including a detailed analysis of the potential energy surface, population decay time scales, and vibrational coordinates coupled to the excitation dynamics. The development presented in this work is a baseline for future implementations of more sophisticated frameworks for simulating of non-adiabatic molecular dynamics in NWChem.

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

Minimal Experimental Bias on the Hydrogen Bond Greatly Improves Ab Initio Molecular Dynamics Simulations of Water

Experiment Directed Simulations (EDS) is a method within a class of techniques seeking to improve molecular simulations by minimally biasing the system Hamiltonian to reproduce certain experimental observables. In a previous application of EDS to ab initio molecular dynamics (AIMD) simulation based on electronic density functional theory (DFT), the AIMD simulations of water were biased to reproduce its experimentally derived solvation structure. In particular, by solely biasing the O-O pair correlation functions, other structural and dynamical properties that were not biased were improved. In this work, the hypothesis is tested that directly biasing the OH pair correlation (and hence the H-O∙∙∙H hydrogen bonding), will provide an even better improvement of DFT-based water properties in AIMD simulations. The logic behind this hypothesis is that for most electronic DFT descriptions of water the hydrogen bonding is known to be deficient due to anomalous charge transfer and over polarization in the DFT. Using recent advances to the EDS learning algorithm, we thus train a minimal bias on AIMD water that reproduces the O-H radial distribution function derived from the highly accurate MB-pol model of water. Finally, it is then confirmed that biasing the O-H pair correlation alone can lead to improved AIMD water properties, with structural and dynamical properties in even closer to experiment than the previous EDS-AIMD model.

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

Convergent Concordant Mode Approach for Molecular Vibrations: CMA-2

The concordant mode approach (CMA) is a promising new scheme for dramatically increasing the system size and level of theory achievable in quantum chemical computations of molecular vibrational frequencies. Here, we achieve advances in the CMA hierarchy by computations targeting CCSD(T)/cc-pVTZ (coupled cluster singles and doubles with perturbative triples using a correlation-consistent polarized-valence triple-ζ basis set) benchmarks within the G2 molecular test set, executing a statistical analysis for 1501 frequencies from 111 compounds and then separately solving the refractory case of pyridine. First, MP2/cc-pVTZ (second-order Møller–Plesset perturbation theory with the same basis set) proves to be an excellent and preferred choice for generating the underlying (Level B) normal modes of the CMA scheme. Utilizing this Level B within the CMA-0A method reproduces the 1501 benchmark frequencies with a mean absolute error (MAE) of only 0.11 cm –1 and an attendant standard deviation of 0.49 cm –1 . Second, a convergent CMA-2 method is constituted that allows efficient computation of higher level (Level A) frequencies to any reasonable accuracy threshold by using only Hartree–Fock (HF) and MP2 or density functional theory (DFT) data to generate ξ parameters, which select the sparse off-diagonal force field elements for explicit evaluation at Level A. When Level B = MP2/cc-pVTZ, a cutoff of ξ = 0.02 provides an average maximum absolute error per molecule of only 0.17 cm –1 by incurring merely a 33% increase in average cost over CMA-0A. This CMA-2 method also eradicates the 4 problematic CMA-0A outliers of pyridine with even less effort (ξ = 0.04, 22% increase). Finally, the newly developed CMA procedures are shown to be highly successful when applied to 1-(1H-pyrrol-3-yl)ethanol, a new test molecule with diverse types of vibration.

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

Effective Fragment Potential-Based Molecular Dynamics Studies of Diffusion in Acetone and Hexane

To facilitate more reliable descriptions of transport properties in liquids, molecular dynamics (MD) simulations are performed based on the effective fragment potential (EFP) method derived from first-principles quantum mechanics (in contrast to MD based upon empirically fitted potentials). The EFP method describes molecular interactions in terms of Coulomb, polarization/induction, dispersion, exchange-repulsion, and charge-transfer interactions. The EFP MD simulations described in this paper, performed on hexane and acetone, are able to track the mean-square displacement of molecules for sufficient time to reliably extract translational diffusion coefficients. Finally, the results reported here are in reasonable agreement with experiment.

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

Reorganization of Water at Aqueous Aluminum Chloride (AlCl 3 ) Interfaces: Vibrational Sum Frequency Generation and Molecular Dynamics Simulations

AlCl 3 hydration states and complexation are not well understood both in solutions and at the air–aqueous interface despite their potential significance in natural waters and their industrial and energy-related applications. Here, we investigated Al 3+ and Cl – ion behaviors in an AlCl 3 aqueous bulk solution and at the air–aqueous interface using interface-selective vibrational sum frequency generation (SFG), Raman and infrared spectroscopies, molecular dynamics (MD) simulation, as well as molecular-informed reduced modeling. Our reduced modeling reveals relatively long-range effects for Al 3+ as compared to monovalent ions such as Na + indicating that the interfacial depth of trivalent ions can be significantly larger than that of monovalent ions at the air–water interface. MD simulations reveal interfacial stratification and multiple layering of the ions. Compression of the Al 3+ and Cl – distributions with increasing concentrations from 0.5 to 2.5 m is also observed in the subsurface regions. Significant SSP- and PPP-polarized SFG OH spectral intensity increases are observed from 0.5 to 1.5 m and 0.5 to 2.5 m, respectively, indicative of interfacial depth increases and a change in average orientation above 1.5 m. Extensive evaluation of SFG spectra, Fresnel-corrected using several approaches, shows the same trends. The nonmonotonic trend points to a changing structure in surface and subsurface water orientation and hydrogen bonding environment generally consistent with the MD simulation of stratification and water orientation changes. Furthermore, solvent-shared ion pairing is implicated with MD simulation radial distribution analysis and consistent with infrared spectral identification of the hexaaqua aluminum ion in the solution phase. Spectral evidence of a strong Al 3+ hydration shell and the acidic behavior of the Al 3+ ions is obvious in the Raman and infrared spectra of the bulk solution. In conclusion, we show that the MD dipole potential is directly related to the MD second-order susceptibility of the interface, χ SFG–MD (2) , both of which correlate up to ∼35 Å with the spectral observations of increasing and then saturating intensities, suggesting that both ion stratification and interfacial depth determine the water orientations at an air–water interface of 1-3 electrolyte solutions.

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

DFT-Based Calculation of Molecular Hyperpolarizability and SFG Intensity of Symmetric and Asymmetric Stretch Modes of Alkyl Groups

Vibrational sum frequency generation (SFG) spectroscopy has been extensively used for obtaining structural information of molecular functional groups at two-dimensional (2D) interfaces buried in the gas or liquid medium. Although the SFG experiment can be done elegantly, interpreting the measured intensity in terms of molecular orientation with respect to the lab coordinate is quite complicated. One of the main reasons is the difficulty of determining the hyperpolarizability tensors of even simple molecules that govern their SFG responses. The single-bond polarizability derivative (SBPD) model has been proposed to estimate the relative magnitude of SFG-active hyperpolarizability by assuming that the perturbation associated to each vibration is negligible. In this study, density functional theory (DFT) was used to calculate the polarizability and dipole derivative tensors of the CH 3 stretch mode of CH 3 I, CH 3 CH 2 I, CH 3 OH, CH 3 CH 2 OH. Then, the hyperpolarizability tensors of symmetric and asymmetric vibration modes were calculated considering the Boltzmann distribution of representative conformers, which allowed theoretical calculation of their SFG intensities at all polarization combinations as a function of the tilt angle of the CH 3 group with respect to the surface normal direction. Then the ratios of the calculated SFG intensities for the CH 3 peaks used in experimental studies for the CH 3 tilt angle determination were compared. This comparison clearly showed the effect of vibrational coupling among neighboring functional groups. As a result, this study presents new parameters that can be used in determining the average tilt angle of the CH 3 group at the 2D interface with SFG measurements as well as limitations of the method.

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

Accelerated Sequence Design of Star Block Copolymers: An Unbiased Exploration Strategy via Fusion of Molecular Dynamics Simulations and Machine Learning

Star block copolymers (s-BCPs) have potential applications as novel surfactants or amphiphiles for emulsification, compatibilization, chemical transformations, and separations. s-BCPs have chain architectures where three or more linear diblock copolymer arms comprised of two chemically distinct linear polymers, e.g., solvophobic and solvophilic chains, are covalently joined at one point. The chemical composition of each of the subunit polymer chains comprising the arms, their molecular weights, and the number of arms can be varied to tailor the surface and interfacial activity of these architecturally unique molecules. Further, this makes identification of the optimal s-BCP design nontrivial as the total number of plausible s-BCP architectures is experimentally or computationally intractable. In this work, we use molecular dynamics (MD) simulations coupled with a reinforcement learning-based Monte Carlo tree search (MCTS) to identify s-BCP designs that minimize the interfacial tension between polar and nonpolar solvents. We first validate the MCTS approach for the design of small- and medium-sized s-BCPs and then use it to efficiently identify sequences of copolymer blocks for large-sized s-BCPs. The structural origins of interfacial tension in these systems are also identified by using the configurations obtained from MD simulations. Chemical insights into the arrangement of copolymer blocks that promote lower interfacial tension were mined using machine learning (ML) techniques. Overall, this work provides an efficient approach to solve design problems via fusion of simulations and ML and provides important groundwork for future experimental investigation of s-BCPs for various applications.

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

Effect of Solvents on Lignin–Surface Interactions via Molecular Dynamics Simulations

Lignin, an essential building block of lignocellulosic biomass, is a potential abundant source of aromatic monomers for the polymer and chemical industry. Reductive catalytic fractionation (RCF) is one promising process that can produce high yields of phenolic monomers and oligomers from lignin under different catalytic conditions. An important choice in optimizing RCF is the selection of solvent; however, detailed insights into solvent effects on lignin behaviors and interactions remain limited. Here, in this work, we perform all-atom molecular dynamics simulations to study the solvation of lignin, solvent-mediated conformational changes, and the interaction of solvated lignin oligomers with model surfaces. We focus on the behavior of an oligomeric lignin model compound in methanol, ethanol, a binary mixture of ethanol and water, and water at both the RCF reaction temperature (473 K) and room temperature. Analysis of structural features of lignin suggests that these three organic solvent systems favorably solvate lignin, resulting in a more extended conformation suitable for catalytic conversion to valuable chemicals. We further introduce model palladium (Pd) and carbon (C) surfaces to understand how solvent choice impacts adsorption onto a representative catalytic surface and support, and to quantify the competition among the reactant and solvent molecules for the surface. Unbiased simulations suggest strong adsorption of lignin on both Pd and C surfaces at 473 K, with notable solvent-mediated differences in adsorption energies. Additionally, our findings indicate that lignin adsorption is promoted by the entropy change resulting from the displacement of solvent molecules from the surface. This study provides a molecular perspective of adsorption of lignin onto varying surfaces, which is a step towards understanding and optimizing the catalytic conversion of lignin into valuable chemicals.

adsorption↗