Search NASA⌕ Search

SEARCH · Search NASA

Results for “Coarse-graining”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Coarse-grained resource allocation modeling for decoding and rewiring microbial metabolism

Microbial metabolism is a complex, emergent system driven by the coordinated interplay of intricate and dynamic molecular processes. To elucidate cellular behavior and enable biotechnological applications, quantitative models that address the inherent complexity of metabolism have been developed from a resource allocation perspective. Here, we synthesize recent advances in coarse-grained resource allocation frameworks and their applications in understanding microbial physiology and guiding gene circuit design. Here, these frameworks reveal global regulatory constraints and predict cellular adaptation to nutrient and environmental changes. In addition, they enable the quantification of metabolic costs, the dissection of circuit–host interactions, and the development of strategies for burden mitigation. Collectively, these modeling frameworks provide a powerful platform for uncovering quantitative principles of microbial growth and engineering robust synthetic biological systems.

coarse-grained modeling↗
DOE OSTI2026

Challenges of conventional iterative all-atom and coarse-grained multiscale molecular dynamics

In this work, we evaluate the biomolecular dynamics behaviors when conventionally iterating between all-atom (AA) and coarse-grained (CG) molecular dynamics (MD) simulations over multiple cycles. We implemented the workflow to iterate between AA and CG in OpenMM, namely the iterative multiscale MD (iMMD) simulation workflow. In particular, we aim to identify practical applications for iterating between AA and CG simulations in a conventional manner without any constraints or model modifications. We evaluate the iMMD workflow on four representative systems, spanning folding of two soluble proteins and protein-protein as well as protein-lipid interactions of two membrane proteins. We observe that iteration between AA and CG representations could help the soluble proteins exit undesirable metastable states to fold, resulting from random protein structural distortions due to cycling. Consequently, the most reliable use of iterative AA and CG simulations appears to be to accelerating complex lipid mixing for membrane-bound protein systems rather than sampling protein conformational space. Our work explores the practical usages and limitations for iterative AA and CG simulations using readily available AA and CG force fields. The evaluated iMMD workflow in OpenMM is made available at https://github.com/lanl/iMMD.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Coarse-Grained Simulations of Polyrotaxane Hydrogels under Quiescent and Shear Conditions

Cyclodextrin-based polyrotaxanes (PR) form hydrogels in water through cyclodextrin (CD) aggregation and crystallization. These networks can break under shear flows, making them versatile platforms for extrusion-based 3D printing. To optimize the material properties of 3D-printed PR gels, a microscopic understanding of the structural evolution during 3D printing is necessary. Here, we employ coarse-grained (CG) simulations to reveal the PR assembly process at the molecular level. Our simulations reproduce the experimental crystal morphologies of PR at varying concentrations and chain lengths under quiescent conditions. Using nonequilibrium simulations, we show shear flow ruptures crystalline domain connectivity in PR gels and stacks the lamellae in the gradient direction, allowing the materials to flow during 3D printing. After the cessation of flow, the anisotropic crystal alignment and absence of available dangling PR diminish the intercrystal connectivity. The printed materials are therefore mechanically weaker than the pristine hydrogels, in agreement with experimental results. Nonetheless, by relating the microscopic structural evolution with viscoelastic properties of PR gels and solutions, we elucidate how flow conditions and sample composition affect the 3D printing performance of PR hydrogels.

Smith, Cameron D. [Dartmouth College, Hanover, NH ↗
DOE OSTI2026

Temporal Coarse Graining for Classical Stochastic Noise in Quantum Systems

Simulations of quantum systems with Hamiltonian classical stochastic noise can be challenging when the noise exhibits temporal correlations over a multitude of time scales, such as for 1/f noise in solid-state quantum information processors. Here we present an approach for simulating Hamiltonian classical stochastic noise that performs temporal coarse-graining by effectively integrating out the high-frequency components of the noise. We focus on the case where the stochastic noise can be expressed as a sum of Ornstein-Uhlenbeck processes. Temporal coarse-graining is then achieved by conditioning the stochastic process on a coarse realization of the noise, expressing the conditioned stochastic process in terms of a sum of smooth, deterministic functions and bridge processes with boundaries fixed at zero, and performing the ensemble average over the bridge processes. For Ornstein-Uhlenbeck processes, the deterministic components capture all dependence on the coarse realization, and the stochastic bridge processes are not only independent but taken from the same distribution with correlators that can be expressed analytically, allowing the associated noise propagators to be precomputed once for all simulations. This combination of noise trajectories on a coarse time grid and ensemble averaging over bridge processes has practical advantages, such as a simple concatenation rule, that we highlight with numerical examples.

Albash, Tameem [Sandia National Lab. (SNL-NM), Alb↗
DOE OSTI2026

Topology-Informed Design Rules for Deconstructable Thermoset Copolymer Networks

Existing models of thermoset deconstruction facilitated by incorporating cleavable comonomers rely on a mean-field reverse gel point paradigm, which predicts network dissolution once cleavable bonds reach a critical stoichiometric threshold, but does not account for where those bonds reside within the network architecture. Using reactive coarse-grained molecular dynamics simulations coupled with graph-theoretic analysis, we extend this stoichiometric picture to show that deconstructability is governed by the curing-imprinted network topology rather than stoichiometry alone. This topological organization is hierarchical: at the local scale, the elastic effectiveness of cross-link junctions determines which cross-links constitute the load-bearing scaffold; at the mesoscale, the cross-linking rate kinetically templates that scaffold into topologically modular communities─densely cross-linked clusters connected by sparse bridging strands that sustain network connectivity. Using betweenness centrality to identify nodes that disproportionately lie on intercommunity shortest paths, we demonstrate that effective deconstruction of the network into macromolecular fragments requires cleavable comonomers to intercept these high-centrality bridging strands. We further find that under uniform, disassortative comonomer incorporation, this topological requirement provides a mechanistic basis for extending the reverse gel point to incorporate network topology. We also show that modularity imposes a fundamental limit on fragment uniformity that persists even when the centrality requirement is met. Finally, we demonstrate that chain stiffness provides a nearly independent lever to suppress mechanically redundant cross-links and raise the glass transition temperature without significantly altering the deconstruction outcome. Together, these findings reframe the thermoset design space around network topology and provide actionable guidelines for engineering thermoset copolymers with predictable deconstructability and targeted thermomechanical performance.

coarse-grained molecular dynamics↗
DOE OSTI2026

Can classical DEM simultaneously capture compressibility and flowability of milled biomass?

Accurate prediction of the rheological behavior of biomass is essential for the design and operation of hoppers, feeders, and storage systems in biorefineries. This study examines whether the classical, coarse-grained discrete element method (DEM) formulation can simultaneously reproduce the compressibility and flowability of milled herbaceous biomass, using Miscanthus × giganteus as a representative material. The model represents particles as rigid spheres interacting through Hertz-Mindlin elastic-frictional contacts augmented with an area-dependent cohesion term. Laboratory cyclic compression and wedge-shaped hopper discharge experiments were used as calibration benchmarks. Although the model can independently reproduce each behavior by appropriately tuning particle Young's modulus E and cohesion energy density k, an extensive parametric investigation comprising more than 600 simulations reveals that the optimal parameter regions for compression and hopper flow are distinct and non-overlapping in (E, k) space. Surrogate surface analysis further shows that the corresponding objective-function valleys exhibit similar trends but are approximately parallel and spatially offset, precluding a unified calibration within the explored domain. Sensitivity analysis indicates that compressibility is governed predominantly by stiffness and cohesion, whereas the slope of the mass flow rate-opening relation in hopper discharge is primarily controlled by tangential friction. Extensions incorporating particle size distribution and clumped-sphere representations do not eliminate the incompatibility. These results systematically reveal, for the first time, the structural limitation of simplified DEM formulations in representing biomass rheological behavior, underscoring the necessity for models incorporating additional physical mechanisms, such as particle deformability or enhanced interlocking, to achieve unified predictive capability for biomass handling behavior.

09 BIOMASS FUELS↗
DOE OSTI2026

Retention and surface morphology evaluation of fine-grain dispersion-strengthened tungsten for plasma-facing component applications

This study exposed novel fine-grain dispersion-strengthened tungsten (W) to high fluence, low energy deuterium (D) and helium (He) plasmas to evaluate how material microstructure and composition affect hydrogen retention and surface morphology. Tested materials included fine-grain dispersion-strengthened tungsten (DSW) with 3 wt% zirconium carbide (ZrC) dispersoids, fine-grain dense W without any dispersoids (FGW), and coarse-grained polycrystalline ‘ITER-grade’ W. Samples were exposed to D 2 + and He + plasmas at fusion-relevant fluences (∼10 25 m -2 ) and ion energies (75 eV) over a range of temperatures (200 °C, 300 °C, 450 °C for D, 850 °C for He). Helium ion microscopy was performed on the exposed samples to evaluate surface morphology changes and material integrity. After D plasma exposure, the ZrC dispersoids showed near-surface degradation at exposure temperatures above 300 °C, but no detrimental morphology changes were observed for the adjacent W grains. After He plasma-exposure, nano-structured fuzz formation was observed in the tungsten matrix of all samples. The ZrC dispersoids maintained their integrity despite the surrounding fuzz growth, with clear delineation between the W fuzz and dispersoid regions. Thermal desorption spectroscopy showed that ZrC DSW consistently retained more D than the FGW by about a factor of 2 across all temperatures. At 200 °C and 300 °C, the ITER-W displayed lower D retention than both the DSW and FGW, however at 450 °C ITER-W showed the highest retention, about 50% more than DSW. He retention was comparable across all samples, with the highest retention observed in the fine-grain W, only 26% higher than in ITER-W. These insights on retention behavior will inform further optimization of these novel fine-grained tungsten materials with and without dispersoid additives.

Dispersion-strengthened tungsten↗
DOE OSTI2026

Molecular Modeling of Surfactant Interaction on Phospholipid Bilayers Mimicking Corneal Epithelium

Surfactants found in consumer products can compromise eye corneal membrane integrity upon accidental exposure. Traditional in vitro and in vivo approaches to evaluate membrane–surfactant interaction pose experimental limitations such as species variability, reproducibility, and most often do not provide the overall picture. These limitations motivate the use of in silico models to study phenomena like cellular disruption assays caused by surfactants at the molecular scale. In this work, coarse-grained molecular dynamics simulations have been employed to investigate how nonionic alcohol ethoxylate (AE) and anionic surfactant alcohol ethoxy sulfate (AES) interact with lipid bilayer liposomes that mimic corneal epithelial cell membranes. The spherical liposome is composed of 1,2-dihexadecanoyl-sn-glycero-3-phosphocholine (DPPC), 1,2-di(9Z-octadecenoyl)-sn-glycero-3-phosphoethanolamine (DOPE), 1,2-di(9Z-octadecenoyl)-sn-glycero-3-phospho-l-serine (DOPS), and cholesterol, resembling the composition of the corneal epithelial cells’ membrane bilayer. The simulation consisted of varying degrees of representative surfactant compositions and two initial types of surfactant configurations within or outside the liposome. Our results reveal that both surfactants induce outer leaflet bulging, agreeing with membrane solubilization models. The more highly ethoxylated surfactant, AE, caused more consistent inner leaflet disruption than AES, resulting in significantly more water permeation and membrane thinning. In addition, both surfactants increase the lateral diffusion of lipids within the membrane layers, with higher ethoxylated AE showing a stronger effect than AES. This study demonstrates how surfactant structure and localization influence bilayer membrane integrity, offering mechanistic insights into the irritation potential, thus guiding the rational design of effective surfactant-based formulations.

Lipids↗
DOE OSTI2026

Development of an Amine Oxide Polyzwitterion Brush Martini Model with Polarizable Water and Ions

Abstract In order to accurately simulate the fouling process of proteins onto polyzwitterion brushes, models that accurately capture the hydration properties and chain conformations of such brushes must first be established. We developed a Martini coarse-grained (CG) model for amine oxide polyzwitterion (PNOMA) brushes, a promising class of antifouling materials, in polarizable water and ions by fitting to all-atom bond and angle distributions, monomer hydration free energy, monomer–monomer distance potential of mean force (PMF), and monomer–salt radial distribution functions (RDFs). Martini 2.2P was selected for compatibility with the established polarizable water and ion models. For comparison with PNOMA, we also constructed models for conventional sulfobetaine (PSBMA) and phosphorylcholine (PMPC) polyzwitterions and the polycation PMETAC using established nonbonded bead types from the literature and refitting bond and angle potentials. We simulated each polymer brush chemistry for varying grafting density and chain length, validating brush height scaling relations against experimental data. The CG models captured the relative hydration strengths among different polyzwitterion chemistries, and brush heights extrapolated to higher molecular weights are in agreement with experimental ellipsometry data. We find that chain swelling of the superhydrophilic PNOMA brushes lies between that of the traditional polyzwitterions PSBMA/PMPC and the polycation PMETAC. For PNOMA brushes in NaCl solution, simulated brush height decreases with salt concentration due to the selectively strong interactions between amine oxide and sodium ions.

Walker, Christopher C. [Oak Ridge National Laborat↗
DOE OSTI2026

Molecular Mechanisms Underlying Surfactant-Based Plastics De-Inking

Surfactant-mediated binder removal is critical for de-inking pretreatment in plastics recycling, yet the molecular mechanisms governing surfactant performance remain poorly understood. We used all-atom (AA) and coarse-grained (CG) molecular dynamics (MD) simulations alongside alkaline surfactant washing experiments to investigate interactions between a series of surfactants and a polyether urethane (PEU) binder in solution and on a polyethylene (PE) surface. Experiments reveal a range of de-inking efficiencies, ranging from <25 to >95% depending upon surfactant headgroup charge and tail length. AA simulations reveal that charged surfactants reach stable levels of surfactant coverage, while nonionic surfactants aggregate on the binder. CG umbrella sampling calculations quantify the thermodynamics of binder desorption in water. In ∼0.25 M surfactant solutions, up to a 52% reduction in the free energy barrier is computed, with trends in good agreement (R 2 = 0.92, Pearson’s r = –0.96, Spearman’s ρ = –0.83) with experimental de-inking efficiencies. We find that charged surfactants are more effective than nonionic surfactants for de-inking and propose three regimes of surfactant de-inking processes: good de-inking occurs in surfactants that promote PEU desorption with a low radius of gyration (R g ); moderate de-inking occurs when surfactants stabilize PEU but increase R g ; poor de-inking occurs in surfactants that aggregate on the binder and promote extensive anchoring to the surface. Overall, these molecular-level insights have the potential to guide the design of surfactant formulations for plastics recycling applications.

ink removal↗
DOE OSTI2026

A Robotic High-Throughput Grid-Search Platform for Mapping Phase Behavior in Triblock Copolymer–Homopolymer Blends

We present a high-throughput experimental investigation of the phase behavior in triblock copolymers (PS-b-PB-b-PS and PS-b-PI-b-PS) and polystyrene (PS) homopolymer blends as a function of homopolymer molecular weight (MW) and blend ratio. Using a robotic thin-film processing platform (NOVA) integrated with Grazing Incidence Small-Angle X-ray Scattering (GISAXS) and Atomic Force Microscopy (AFM), we systematically mapped the order–disorder transition (ODT) boundaries and domain spacing evolution across a broad MW range (4.0–101.3 kDa) with varying homopolymer loadings (10% to 90%). The results reveal three distinct regimes: low-MW homopolymers, corresponding to the wet-brush regime produced only gradual domain swelling before disordering at high blend ratios (weight fraction); medium-MW homopolymers, corresponding to thedry-brush regime induced significant domain spacing increase up to 80% followed by earlier disordering, while high-MW homopolymers led to macrophase separation with minimal changes in domain spacing. Additionally, coarse-grained molecular dynamics simulations confirmed our experimental finding that in the low-MW region, the PS homopolymer uniformly distributed in the PS domain. These findings demonstrate that homopolymer molecular weight critically governs both the extent of domain swelling and the onset of disorder in triblock copolymer systems. This high-throughput platform enables the rapid mapping of composition–morphology relationships and can be integrated with AI/ML tools for designing next-generation nanostructured polymers.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Why Seeding Works When Nucleation Barriers Vanish

Crystallization is a process governed by the interplay between nucleation and growth. While crystalline seeds are known to reduce nucleation barriers and accelerate crystallization under nucleation-limited conditions, their influence when nucleation is not the limiting step remains poorly understood. This creates a mechanistic puzzle in systems where nucleation barriers are already negligible, yet seeding still accelerates crystallization. Here, the synthesis of zeolites is a quintessential example of growth-limited crystallization in which seeds accelerate the process despite negligible homogeneous nucleation barriers. Using coarse-grained molecular dynamics simulations─validated across two zeolites and the unrelated case of ice crystallization─we establish that growth-limited crystallization produces many small, misoriented crystallites whose slow coarsening into larger domains controls the emergence of X-ray-detectable crystallinity. Local crystalline order, structural coherence, and X-ray detectability are therefore kinetically decoupled milestones: the first can be reached rapidly while the latter two lag significantly. Seeds resolve this lag by imposing a common orientational registry on nascent crystallites, enabling their coherent coalescence into large seed-bound domains. This early coherence-building step produces a crystallite-size asymmetry that accelerates subsequent coarsening, advancing the onset of X-ray-detectable crystallinity without necessarily increasing the nucleated fraction. We conclude that under growth-limited conditions, the apparent induction period observed in powder X-ray diffraction reflects the time required to build long-range coherence, not the time to form crystalline material. Accordingly, seeds function not primarily by reducing nucleation barriers but by enforcing spatial coherence, thereby shortening the time required to develop long-range order detectable by X-ray diffraction.

Growth-Limited Crystallization↗
DOE OSTI2026

Interfacial Inversion of Stealth Surfactants

Amphiphilic macromolecular surfactants segregate to liquid–liquid interfaces, thereby reducing the interfacial tension and free energy. Here, we investigated “stealth surfactants” in the form of core–shell bottlebrush polymers comprised of pH-responsive diblock copolymer side chains forming a hydrophilic core and a hydrophobic shell, enabling solubility in oil. At liquid–liquid interfaces, these polymers undergo a structural “inversion”, with hydrophilic blocks segregating into the aqueous phase and hydrophobic blocks residing in the oil phase. The reconfiguration kinetics and surfactant properties are influenced by multiple factors, including the molecular weights of the backbone and side chain components, the hydrophilic-to-hydrophobic balance of the side chains, and the pH of the aqueous phase. An observed nonmonotonic dependence of interfacial tension with time is attributed to a progressive structural inversion, where the projected area of the macromolecule onto the interface decreases. To validate this inversion hypothesis, interfacial properties were characterized by sum-frequency generation vibrational spectroscopy, which revealed configurational changes of the core–shell bottlebrush polymers at the fluid interface and revealed a pH-dependent interfacial coverage. Coarse-grained molecular dynamics simulations supported these experimental findings, showing that the pH-responsive core and hydrophobic shell assume a time-averaged configuration with orientations parallel and perpendicular to the plane of the interface, respectively. These findings open routes to design multistimuli-responsive polymeric surfactants and compatibilizers, expanding their potential applications in advanced interfacial systems.

Stealth surfactants↗
DOE OSTI2026

Molecular dynamics based study on the effects of cation size on the local structure and diffusion in polymerized ionic liquids

We have used coarse-grained molecular dynamics simulations to understand the effect of cation size on ion diffusion in polymerized ionic liquids at temperatures well above the glass transition temperature of the polymers. We investigated dependencies of the diffusion constant on the cation radius and static dielectric constant and interpreted these results in terms of underlying structural changes, and decoupling phenomenon. We have found non-monotonic effects of the cation radius on the diffusion constant, with a maximum at an intermediate radius resulting from two different size dependent effects. Changes in the radial distribution functions characterizing spatial distribution of cations with respect to anions result in a monotonic increase in the coordination number. This yields an increase in the effective interaction energy as a function of cation radius, which can explain the decrease in the diffusion constant. This increase is counteracted at small cation radii by a corresponding increase in decoupling of cation and anion dynamics. In addition, the diffusion constant is found to increase on increasing the static dielectric constant for all radii of the cations. These results highlight the significant role of decoupling in designing polymerized ionic liquids with an enhanced diffusion constant of cations.

Gillespie, Colin [ORNL] (ORCID:0000000320287080)↗
DOE OSTI2026

High-throughput small-angle X-ray scattering reveals effective structure factor transitions linked to high-concentration antibody viscosity

High-concentration monoclonal antibody (mAb) formulations are often constrained by elevated viscosity, largely driven by protein–protein interactions, which complicates manufacturing and limits subcutaneous delivery. Early viscosity risk assessment is essential during discovery, yet traditional measurements require large sample volumes, and lack high-throughput capability. Here, we develop a high-throughput small-angle X-ray scattering (SAXS) protocol to detect mAb self-association at dilute concentrations, enabling early predictive insights into high-concentration viscosity. Synchrotron SAXS measurements were conducted for 21 mAbs formulated in a histidine buffer at pH 6.0. An initial subset of 10 mAbs analyzed across 1–150 mg/mL revealed that effective structure factor transitions in the low-q region, indicative of interparticle interactions, consistently emerged below 25 mg/mL. Subsequently, 11 additional mAbs were analyzed at 1–25 mg/mL using automated liquid handling and flow cells to enable high-throughput screening. High-viscosity mAbs exhibited detectable low-q upturns at concentrations ≤10 mg/mL, whereas low-viscosity mAbs showed downturns. A classification criterion based on effective structure factor transitions accurately classified all high- and low-viscosity mAbs at 150 mg/mL, offering a scalable, sample-efficient alternative to conventional methods. These results extend recent findings on the concentration-dependent sensitivity of SAXS to short-range attractions, demonstrating that they can emerge at lower concentrations than previously reported. This study presents the most comprehensive and diverse SAXS dataset for mAbs reported to date within a single formulation, providing a valuable resource for developing and validating coarse-grained models that can more accurately capture intermolecular interactions governing high-concentration solution behavior, thereby enabling rational antibody engineering and improved developability.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Analog and symbolic computation through the Koopman framework

We develop a Koopman operator framework for studying the computational structure of dynamical systems. Specifically, we show that the resolvent of the Koopman operator provides a natural abstraction of halting, yielding a ‘Koopman halting problem’ that is recursively enumerable in general. For symbolic systems, such as those defined on Cantor space, this operator formulation captures reachability between clopen sets, while for equicontinuous systems we prove that the Koopman halting problem is decidable. Our framework demonstrates that absorbing (halting) states in coarse-grained finite automata correspond to Koopman eigenfunctions with eigenvalue one, while cycles in the transition graph impose spectral constraints associated with periodic dynamics. These results provide a unifying perspective on computation in symbolic and analog systems, showing how computational universality is reflected in operator spectra, invariant subspaces, and algebraic structures. Beyond symbolic dynamics, this operator-theoretic lens opens pathways to analyze the computational properties of a broader class of dynamical systems, including polynomial and analog models, and suggests that computational hardness may admit dynamical signatures in terms of Koopman spectral structure.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Seamlessly joining length scales: From atomistic thermal graphs to anisotropic continuum conductivity

Thermal transport in complex solids is governed by local structure, defects, and anisotropy, yet most continuum models still rely on oversimplified and homogenized conductivities. Here, we bridge atomistic and continuum descriptions by building finite element (FE) models directly from the site-projected thermal conductivity (SPTC), an atomic-level decomposition of the Green–Kubo thermal conductivity. We introduce a toolkit, the “Simulator Collection for Atomic-to-Continuum Scales (SCACS)”, which uses a graph neural network to predict SPTC on large atomic structures, coarse-grains these fields into anisotropic conductivity tensors, and embeds them into the heat-flow FE equation with a customized, anisotropy-aware adaptive mesh refinement scheme. Applied to silicon nanostructures, the resulting FE models act as representative volume elements, reproduce bulk conductivities, and capture interfacial and defect-driven anisotropy while maintaining thermodynamic consistency. Additionally, SCACS predicts experimental conductance trends and fields. This work demonstrates a general route for transferring atomistic transport information into device-scale thermal simulations with physics-based approximations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗