SEARCH · Search NASA

Results for “Dynamic Structural Analysis”

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 235 records · Page 13

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Entropy Analysis of FPGA Interconnect and Switch Matrices for Physical Unclonable Functions

Random variations in microelectronic circuit structures represent the source of entropy for physical unclonable functions (PUFs). In this paper, we investigate delay variations that occur through the routing network and switch matrices of a field-programmable gate array (FPGA). The delay variations are isolated from other components of the programmable logic, e.g., look-up tables (LUTs), flip-flops (FFs), etc., using a feature of Xilinx FPGAs called dynamic partial reconfiguration (DPR). A set of partial designs is created to fix the placement of a time-to-digital converter (TDC) and supporting infrastructure to enable the path delays through the target interconnect and switch matrices to be extracted by subtracting out common-mode delay components. Delay variations are analyzed in the different levels of routing resources available within FPGAs, i.e., local routing and across-chip routing. Data are collected from a set of Xilinx Zynq 7010 devices, and a statistical analysis of within-die variations in delay through a set of the randomly-generated and hand-crafted interconnects is presented.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2026

Object-Based Evaluation of Dynamical and Statistical Downscaled Precipitation Products over CONUS

High-resolution precipitation data, generated through dynamical downscaling (DD) or statistical downscaling (SD) of global climate model output, provide critical information for regional climate assessment and adaptation planning. Most downscaling development and validation have focused on accurate gridscale precipitation construction and ignored the spatial structure of precipitation across model grids and at the event scale. However, many applications, e.g., hydrologic modeling and the analysis using the downscaled precipitation, require a reasonable representation of the spatial structure of precipitation within watersheds. Therefore, a set of standard metrics to evaluate the representation of the spatial structure of individual storms across diverse downscaled precipitation products is desired. To address this need, we conducted an object-based evaluation of precipitation in decades-long DD and SD products over the contiguous United States (CONUS). Specifically, we evaluate their ability to reproduce various features of precipitation objects in the observations: total volume, precipitation area, peak intensity, and spatial structure. Multiple metrics (bias, Perkins score, and nonparametric statistical tests) are used to quantify model performance. Our evaluation reveals notable variations in performance among individual products across different climate zones and seasons, as well as between extreme and nonextreme events. In general, most DD products exhibit balanced performance across the four precipitation object features, while SD products vary more significantly in their performance across products. Based on this comprehensive evaluation, we provide guidance on choosing downscaled products for specific regions, seasons, and precipitation object features. These findings and recommendations can inform precipitation-relevant modeling and analysis over CONUS, guide future downscaling technique developments, and provide actionable information for climate impact assessment and adaptation.

Downscaling↗

DOE OSTI2024

Formation and coarsening of epitaxially-supported metal nanoclusters

This mini-review describes developments over the last ~30 years in characterizing the nucleation & growth of epitaxially-supported metal nanoclusters (NCs) or islands during vapor deposition, as well as their post-deposition coarsening. A beyond-mean-field treatment for homogeneous nucleation & growth corrects the deficiencies of traditional treatments in describing, e.g., the island size distribution, but also necessitates consideration of the spatial distribution of islands and their capture zones. We discuss advances in modeling capabilities, including those based upon on an ab-initio level treatment of periphery diffusion kinetics, for description of the non-equilibrium growth shapes of these NCs, focusing on 2D NCs. For post-deposition coarsening of arrays of NCs, there is generally a competition between Ostwald Ripening (OR) and Smoluchowski Ripening (SR). SR is also known as Particle Migration & Coalescence. For 2D NCs in homoepitaxial systems, conventional OR is observed on pristine fcc(111) surfaces, dramatically enhanced OR in the presence of even trace amounts of chalcogens for Cu(111) and Ag(111), and anomalous OR on anisotropic fcc(110) surfaces. The unexpected discovery of SR for fcc(100) homoepitaxial systems prompted extensive analysis of the underlying diffusivities of 2D NCs as a function of size, as well as of NC coalescence dynamics. A comprehensive understanding of these processes is now available. Self-assembly of 3D NCs during deposition, issues related to heterogeneous nucleation, directed assembly, NC growth structure selection, and coarsening are addressed. Finally, for SR of 3D epitaxial NCs, recent insights into the size-dependence of diffusivity are described.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Shear-Induced Conformations of Salt-Free Polyelectrolytes in Semidilute Solutions

Here, we present coarse-grained molecular dynamics simulations of salt-free polyelectrolyte chains in semidilute solutions under simple shear flow, with full hydrodynamic interactions and explicit dipolar solvent. At equilibrium, chain orientation statistics follow a pseudo-Voigt distribution, and the structural correlation length and chain end-to-end vector autocorrelation function exhibit scaling behavior consistent with theoretical predictions for polyelectrolytes. Under shear, chains transition from coiled to stretched states and the end-to-end vector autocorrelation function reveals oscillatory dynamics at high Weissenberg numbers. Analysis of the gyration tensor and shear strain distributions identified three distinct chain populations with directional alignment along and against the shear gradient. Compared with their neutral polymer counterparts, polyelectrolytes exhibit stronger shear thinning and enhanced chain alignment under the same Weissenberg number, which is attributed to electrostatic interactions and shear-induced counterion release. These findings provide molecular insight into the distinct flow response of charged polymers and have implications for tailoring the rheological properties of polyelectrolyte-based materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Insights into distorted lamellar phases with small-angle scattering and machine learning

Lamellar phases are essential in various soft matter systems, with topological defects significantly influencing their mechanical properties. In this report, we present a machine-learning approach for quantitatively analyzing the structure and dynamics of distorted lamellar phases using scattering techniques. By leveraging the mathematical framework of Kolmogorov–Arnold networks, we demonstrate that the conformations of these distorted phases – expressed as superpositions of complex waves – can be reconstructed from small-angle scattering intensities. Through the contour analysis of wave field phase singularities, we obtain the statistics of the spatial distribution of topological defects. Furthermore, we establish that the temporal evolution of these defects can be derived from the time-dependent traveling wave field, informed by the dispersion relation of spectral components. This method opens new avenues for investigating the dynamics of distorted lamellar phases using various dynamic scattering techniques such as neutron spin echo and X-ray photon correlation spectroscopy. These findings enhance our microscopic understanding of how defects influence the physical properties of lamellar materials, with implications for both equilibrium and non-equilibrium states in general lamellar systems.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Strain-Driven Evolution of Structural and Carrier Recombination Dynamics in Germanium

Strain engineering in group-IV semiconductors, in particular germanium (Ge), offers a promising route for advancing next-generation optoelectronic device performance. High minority carrier lifetime serves as a direct indicator of superior material quality and device efficiency. Several factors govern the minority carrier lifetime in both strained Ge and strained GeSn, such as strain level, epilayer thickness, defects, surface roughness, and Sn alloying. In this study, the impact of biaxial strain on effective minority carrier lifetimes of tensile-strained Ge (0% ..epsilon..-Ge/Al?As, 0.78% ..epsilon..-Ge/In0.115?Ga0.885?As, and 0.93% ..epsilon..-Ge/In0.16?Ga0.84?As) and compressively strained GeSn (-0.53% c-Ge0.97?Sn0.03/Al?As) epilayers, grown by molecular beam epitaxy, was experimentally investigated via the microwave-reflection photoconductive decay (..mu..-PCD) technique as a simultaneous function of thickness and strain. The structural and strain-relaxation properties were analyzed by high-resolution x-ray diffraction and Raman spectroscopy, while surface morphology analysis using atomic force microscopy provided insights into the surface roughness of strained Ge-based epilayers. Defect characterization by cross-section transmission electron microscopy revealed superior crystalline quality of ..epsilon..-Ge in one epilayer, while twin boundary formation was observed in the other. We investigated the carrier recombination dynamics in epitaxially strained Ge by extracting the bulk recombination lifetime and surface recombination velocity using the thickness-dependent minority carrier lifetime. High effective minority carrier lifetimes of approximately 60-340 ns were achieved at high strain levels with minimal dependence on the thickness of the strained Ge-based epilayers by means of the ..mu..-PCD technique. Minority carrier lifetimes in Ge-based epilayers exhibit a quadratic dependence on the biaxial strain, ranging from -0.88% (compressive) to +1.6% (tensile). The analysis of the carrier recombination dynamics in strained Ge-based heterostructures offers a pathway to designing group-IV based optoelectronic devices.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Barriers and Opportunities for Energy Technology Adoption in Juneau, Alaska

This report presents findings from a qualitative study examining barriers and opportunities for air source heat pump (ASHP) and electric vehicle (EV) adoption in Juneau, Alaska, with a particular focus on manufactured and multifamily housing. The analysis draws on community insights from end users and middle actors to better understand how technology adoption unfolds in contexts with distinct logistical, infrastructural, and housing constraints. The report is organized according to key barriers and opportunities identified through stakeholder input, providing a structured understanding of adoption dynamics across technologies and housing types. These insights are intended to inform program design and support more effective electrification strategies tailored to local conditions. The study team employed qualitative methods to capture both in-depth and high-level perspectives on technology adoption. Data collection included: 1) two 2-hour focus groups with a total of five end users and seven middle actors, enabling detailed and structured discussion and 2) ten semistructured interviews with manufactured home owners, multifamily landlords, and one tenant, providing complementary insights across housing contexts. Focus groups captured accounts of shared challenges and opportunities while interviews offered more concise reflections on individual experiences. Together, these methods enabled a more comprehensive understanding of both systemic barriers and lived experiences with ASHPs and EVs. The findings reveal that adoption of electrification technologies is shaped by a combination of economic, logistical, and informational factors that vary across housing types, technology characteristics, user groups, and other demographic factors. Addressing these factors requires tailored strategies that reflect local conditions and user experiences. The insights in this report can provide a foundation for organizations such as AEL&P to refine program design, support more effective outreach, and anticipate shifts in energy demand associated with increased electrification. More broadly, the study highlights the importance of incorporating community perspectives when developing electrification initiatives to ensure they are both practical and responsive to real-world constraints.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2026

Structural Aspects of Neutron Survival Probabilities

The neutron survival probability (and related quantities including probabilities of extinction and initiation) is a central element of the broader stochastic theory of neutron populations and finds application in fields including reactor start-up, analysis of reactor power bursts and criticality accidents, and safeguards. In a full neutron transport formulation, the equation governing the single-neutron survival probability is a backward or adjoint-like integro-partial differential equation with the added complexity of being highly nonlinear. Analogous formulations of this equation exist in the context of many approximate theories of neutron transport, with the point kinetics formulation having received significant theoretical attention since the 1940s. This work continues this tradition by providing a novel analysis of the single-neutron survival probability equation using the tools of boundary layer theory. The analysis reveals that the “fully dynamic” solution of the single-neutron survival probability equation—and some key probability distributions derived from it—may be cast as a singular perturbation around the underlying quasi-static single-neutron probability of initiation. In this perturbation solution, the expansion parameter is the ratio of the neutron generation time to a macroscopic time scale characterizing the overall system evolution; this interpretation illuminates some of the fundamental structural aspects of neutron survival phenomena.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

Electrochemical Corrosion and Catalysis Dynamics of Tin Oxide during Water Oxidation

Metal oxide corrosion severely limits anodic electrocatalysis, particularly at high potentials in acidic environments, where degradation pathways remain poorly defined. This study establishes explicit connections between corrosion and electrocatalysis on tin oxide during water oxidation by examining the roles of lattice defects, reactive oxygen species, interfacial pH variations, and speciation of corroded tin in acid. We first demonstrate the presence of structural defects such as oxygen vacancies and substoichiometric Sn(II) species by integrating electron paramagnetic resonance spectroscopy, ultraviolet photoelectron spectroscopy, and Mott–Schottky analysis. Kohn–Sham density functional theory calculations reveal that explicit water structures thermodynamically stabilize reaction intermediates and lower reaction overpotentials. Moreover, we propose that water dissociation leads to hydrogen-bonding networks formed by H* and OH* intermediates, which may span the entire catalyst surface and decrease the interfacial pH to drive corrosion. In contrast, the electrochemical generation of reactive oxygen species is shown to play a minor role in catalyst corrosion during water oxidation using inductively coupled plasma mass spectrometry coupled with selective chemical scavengers. Square-wave voltammetry combined with rotating ring-disk electrodes is used to reveal that under open-circuit conditions, only Sn(IV) cations chemically dissolve from tin oxide, while both Sn(IV) and Sn(II) species electrochemically corrode during water oxidation. Our results unveil a dynamic and complicated interplay between corrosive and catalytic pathways on metal oxide electrocatalysts: a decrease in interfacial pH due to water oxidation exacerbates Sn(II)/Sn(IV) corrosion. Subsequently, the electrochemical corrosion of Sn(II)/Sn(IV) facilitates product formation from lattice oxygen, while the redeposition of corroded Sn(II) as Sn(IV) can enable oxygen exchange with water. By elucidating the roles of defects and interfacial chemistry, this work provides a roadmap for engineering improved electrocatalysts that balance activity and stability, a critical step toward scalable and durable energy technologies.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Robust Self-Healing Adhesives Based on Dynamic Urethane Exchange Reactions

Thermoset polyurethanes (PUs) have been successfully reprocessed as covalent adaptable networks (CANs) by catalyzing carbamate exchange. Here we extend bond exchange beyond the internal network cross-links to create a dynamic urethane adhesive. Interfacing PU CANs to substrates with nucleophilic functional groups creates adhesives capable of reversible transcarbamoylation with the substrate, which has not been demonstrated previously by CAN adhesives. Two types of thermoset PU films were synthesized, one containing the green carbamate exchange catalyst Zr(tmhd)4 and the other containing no catalyst. Although otherwise identical in chemical and network properties, as indicated by FT-IR spectroscopy and dynamic mechanical thermal analysis (DMTA), the film containing catalyst showed dynamic bond exchange behavior through stress relaxation analysis. When evaluated as an adhesive, the CAN film exhibited self-healing properties and retained its adhesive strength for five cycles, which is attributed to reversible covalent bonding to the glass substrate. Furthermore, this work expands industrially relevant CANs to structural adhesives and demonstrates their potential value in an application that presently employs PUs as single-use materials.

adhesive↗

DOE OSTI2025

Visualizing degradation mechanisms in a gas-fed CO 2 reduction cell via operando X-ray tomography

We utilize operando X-ray computed tomography, coupled with real-time electrochemical analysis, to reveal the underlying failure mechanisms of membrane electrode assemblies (MEAs) for electrochemical CO 2 reduction (eCO 2 R). Through operando imaging, we can obtain unprecedented insights into the dynamic behavior of the MEA under different operating conditions, revealing critical changes in interface interactions, phase distribution, and structural integrity over time. Our findings identify phenomena giving rise to the transition from CO 2 R to the hydrogen evolution reaction (HER), as evidenced by shifts in cathode potential and CO 2 R selectivity. The formation of inhomogeneous precipitates at the gas diffusion electrode disrupts the CO 2 supply and reduces the active sites for eCO 2 R, resulting in a shift toward H2 production during low current density operation. Additionally, under high current density conditions, rapid water crossover up to the microporous layer/gas diffusion layer promotes the transition from CO 2 R to HER, further shifting cell potential toward anodic direction. Oscillating voltage conditions reveal the dissolution and regrowth of precipitates, providing direct visualization of the competing selectivity of CO 2 R and HER. This work offers new insight into the degradation mechanisms of MEAs, with implications for the design of more durable CO 2 R systems.

Lee, Sol A [California Institute of Technology (Ca↗

DOE OSTI2025

TMD Handbook

This handbook provides a comprehensive review of transverse-momentum-dependent parton distribution functions and fragmentation functions, commonly referred to as transverse momentum distributions (TMDs). TMDs describe the distribution of partons inside the proton and other hadrons with respect to both their longitudinal and transverse momenta. They provide unique insight into the internal momentum and spin structure of hadrons, and are a key ingredient in the description of many collider physics cross sections. Understanding TMDs requires a combination of theoretical techniques from quantum field theory, nonperturbative calculations using lattice QCD, and phenomenological analysis of experimental data. The handbook covers a wide range of topics, from theoretical foundations to experimental analyses, as well as recent developments and future directions. It is intended to provide an essential reference for researchers and graduate students interested in understanding the structure of hadrons and the dynamics of partons in high energy collisions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2026

Insights into Rotational and Translational Dynamics in Mixtures of Ethylene Glycol and Choline Chloride Using Nuclear Magnetic Resonance Techniques

This work examines molecular dynamics and interactions in ethylene glycol–choline chloride (EG–ChCl) mixtures across 0–33 mol % ChCl, spanning the true eutectic region near 17–20 mol % and the commonly used 1:2 formulation. We combine pulsed-field-gradient (PFG) diffusion, fast-field-cycling (FFC) relaxometry, temperature-dependent 13 C T 1 , and nuclear Overhauser effect spectroscopy (NOESY) to disentangle local from macroscopic dynamics. PFG and FFC show that both translational and average rotational motions largely track the strong increase in viscosity with ChCl content, with ethylene glycol consistently diffusing faster than the choline cation and no global dynamical anomaly at the eutectic composition. More subtle, site-specific composition effects nevertheless emerge. The ratio of the diffusion coefficient of the hydroxyl group of choline to the diffusion coefficient of the methyl group of choline displays a shallow minimum in the 17–25 mol % region, indicating a modest change in how the hydroxyl-bearing end of choline samples the underlying translational motion relative to the methyl groups. 13 C T 1 analysis shows that rotational correlation times at 25 °C generally increase with ChCl, reflecting viscosity-coupled slowing, while the CH 2 –N α site exhibits a small but reproducible deviation from this monotonic trend near the eutectic. NOESY spectra at similar compositions reveal enhanced cross-relaxation between EG and choline protons, consistent with increased headgroup–solvent contact density rather than a wholesale structural rearrangement. Overall, our multitechnique study demonstrates that EG–ChCl dynamics are predominantly viscosity-dominated, with the eutectic region acting as a subtle dynamical crossover where specific choline segments become maximally coupled to the hydrogen-bond network. These insights refine the structure–dynamics picture of choline-chloride DESs and provide practical guidance for tuning composition in electrochemical, separation, and catalytic applications.

diffusion↗

DOE OSTI2026

Performing Numerical Analysis of Cybersecurity Options Using Dynamic Risk Analysis Tool EMRALD

Cyberattacks can have many different attack paths, durations, and goals. There are also many different mitigation options involving hardware, software, and/or humans. Considering a cyber threat should involve defense-in-depth methods and a quantitative or numerical evaluation of overall effectiveness against dynamic, time-dependent attacks to make cost and risk-informed decisions. Typical cyberattack modeling methods only provide a qualitative evaluation. The main areas of cybersecurity are confidentiality, integrity, and availability. For companies with cyber-physical systems such as advanced nuclear reactors, cyber-related safety is a requirement set by North American Electric Reliability and the U.S. Nuclear Regulatory Commission. They are also concerned about availability or reliability as a business case. As cyber threats are evolving to a business-for-hire structure, more attacks may focus on disrupting business success and reliability, causing financial and economic stability risk. Companies want to know business reliability and recovery from those threats, and that requires modeling physical behavior of the targets. Dynamic-state-based and Markov-based modeling provides a method for better cyber scenario modeling with different tools having issues such as state-base explosion. Dynamic modeling enables time and conditional features not found in other numerical evaluation methods. EMRALD (Event Modeling Risk Assessment using Lined Diagrams) is a dynamic risk analysis modeling and simulation tool and has features that reduce modeling issues. It has been used to model different time-dependent events including plant behavior and operator procedures. As a general modeling tool, EMRALD can also be used to model cyberattack scenarios with varying mitigation options and quantify effectiveness, producing numerical data for risk-informed decisions. This paper uses EMRALD to demonstrate that dynamic numerical risk analysis can be used for cyber threat modeling to provide insights for design decision-making and optimize defense strategies. Keywords: cyber modeling; cyber-physical systems; numerical cyber modeling

97 - MATHEMATICS AND COMPUTING↗

DOE OSTI2024

SAM: A Modern System Code for Advanced Non-LWR Safety Analysis

The System Analysis Module (SAM), developed at Argonne National Laboratory and by collaborators at other organizations, is for advanced non–light water reactor safety analysis. SAM aims to provide fast-running, modest-fidelity, whole-plant transient analysis capabilities that are essential for fast-turnaround design scoping and engineering analyses of advanced reactor concepts. To facilitate code development, SAM utilizes the MOOSE object-oriented application framework, its underlying finite element library, and linear and nonlinear solvers to leverage modern advanced software environments and numerical methods. SAM aims to solve tightly coupled physical phenomena, including fission reaction, heat transfer, fluid dynamics, and thermal-mechanical responses in advanced reactor structures, systems, and components with high accuracy and efficiency. Finally, this paper gives an overview of the SAM code development, including goals and functional requirements, physical models, current capabilities, verification and validation, software quality assurance, and examples of simulations for advanced nuclear reactor applications.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2025

Effective field theory factorization for diffraction

We derive a factorization formula for coherent and incoherent ep diffraction using the soft collinear effective theory, utilizing multiple power expansion parameters to handle different kinematic regions. This goes beyond the known hard-collinear diffractive factorization to address the small-x Regge dynamics and Pomeron exchange from first principles. The effective field theory analysis also uncovers and factorizes an important irreducible incoherent background generated by color-nonsinglet exchange, dubbed “quasi-diffraction”, for which we calculate the associated Sudakov suppression. For unpolarized scattering we show that there are four diffractive structure functions at leading power, and point out the importance of studying ${F}_{3,4}^D$ through asymmetries, in addition to ${F}_{2,L}^D$. For the quasi-diffractive background, we make model independent predictions for ratios of the corresponding structure functions in a perturbative kinematic region. Our analysis also makes predictions for six leading-power spin-dependent structure functions. Finally, we provide connections to diffractive parton distributions, and assess the Ingelman-Schlein model. Our work lays a path for further QCD-based studies of diffraction.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

A protocol and analysis of year-long simulations of global storm-resolving models and beyond

We propose a protocol to evaluate and analyze year-long simulations of global storm-resolving models (GSRMs). The proposed protocol complements an earlier 40-day simulation protocol under the DYAMOND (DYnamics of the Atmospheric general circulation Modeled On Non-hydrostatic Domains) project to allow the analysis of the seasonal cycle and associated climatic relevant phenomena. This intercomparison aims to reveal how GSRMs, which can simulate mesoscale convective systems (MCSs) in the global domain, reproduce atmospheric large-scale structures related to convection beyond month-long simulations. The intercomparison for one-year simulations is conducted by either atmosphere-only models or atmosphere–ocean coupled models with atmospheric horizontal mesh sizes less than 5 km. We recommend the continuous four seasons from March 2020 to February 2021 as a target period for the intercomparison but with options for many groups to join more flexibly. The output variables are collected at 0.25° resolution, and archives of a small set of native grid variables are encouraged to analyze tropical cyclones and MCSs. Through the proposed global storm-resolving simulation, we will evaluate the climatological distributions of the atmospheric large-scale circulations, such as the Intertropical Convergence Zone (ITCZ), monsoon, midlatitude jets, their time evolution, and the upscale impacts on them. We present sample analyses from a one-year simulation using the 3.5 km mesh Nonhydrostatic Icosahedral Atmospheric Model (NICAM), revealing the realistic zonal contrast of tropical precipitation, no double ITCZ structure, the reasonable midlatitude jet position and intensity but a weak bias of storm track activities, and a warm bias over the Eurasia during boreal winter. We also clarify the cross-scale interaction, such as the effects of cold pools on mean precipitation over the Maritime Continent through the precipitation diurnal cycle and the effects of resolved gravity waves on midlatitude mean flows. The proposed one-year simulation protocol is referred to as the “Sendai Protocol.” This protocol is not unique or definite for evaluating GSRMs; we prospect a hierarchical set of experiments from short-term to multi-year simulations as GSRM intercomparisons.

54 ENVIRONMENTAL SCIENCES↗

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↗