Search NASA⌕ Search

SEARCH · Search NASA

Results for “ring polymer molecular dynamics”

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

Two-dimensional ring polymer molecular dynamics determination of the MnO + + H 2 /D 2 reaction rates on a Δ-machine learned potential energy surface

In this work, we investigate the impact of nuclear quantum effects in the kinetics of the MnO + + H 2 reaction, a prototypical system for gas-phase H 2 activation by transition metal oxide ions. The DFT based potential energy surfaces (PESs) for the lowest-lying quintet and septet spin states reported in our previous work are improved by 2953 newly calculated CCSD(T)/AVDZ points using a delta-machine learning (Δ-ML) method. To examine nuclear quantum effects, the rate coefficients are computed using ring-polymer molecular dynamics (RPMD). Due to the mechanistic complexity of the reaction, two reaction coordinates are necessary to map out the free-energy surface and an extended RPMD rate theory is developed. The calculated RPMD rate coefficients on the new PES are in better agreement with experimental data. The calculated kinetic isotope effects (KIEs) range from 1.6 to 1.8, also in good agreement with the experimental values. Additionally, the converged RPMD rate coefficients are 2.3 to 3.1 times higher than their classical counterparts across the studied temperature range, suggesting the presence of moderate nuclear quantum effects in this reaction.

Liu, Yang [University of New Mexico, Albuquerque, ↗

Mean-Field Ring Polymer Rates Using a Population Dividing Surface

Mean-field Ring Polymer Molecular Dynamics (MF-RPMD) offers a computationally efficient method for the simulation of reaction rates in multi-level systems. Previous work has established that, to model a nonadiabatic state-to-state reaction accurately, the dividing surface must be chosen to explicitly sample kinked ring polymer configurations where at least one bead is in a different electronic state than the others. Building on this, we introduce a population difference coordinate and a kink-constrained dividing surface, and we test the accuracy of the resulting mean-field rate theory on a series of bilinear coupling model systems. Here, we demonstrate that this new MF-RPMD rate approach is efficient to implement and quantitatively accurate for models over a wide range of driving forces, coupling strengths, and temperatures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Path integral molecular dynamics: A high-fidelity approach to quantum dynamics of electrons

We investigate electron transport in the uniform electron gas using ring-polymer molecular dynamics (RPMD). Working in the weakly coupled, non-degenerate regime, we use RPMD to probe how the onset of quantum diffraction effects at high temperature reshapes electron–electron collisions and leads to a classical-to-quantum crossover in macroscopic transport properties. Static thermodynamics obtained with RPMD are consistent with the weak-coupling equation of state, confirming correct quantum Boltzmann sampling. Real-time transport extracted from mean square displacements exhibits the expected ballistic-to-diffusive transition and a systematic reduction of the electronic self-diffusivity as quantum effects strengthen, due to quantum diffraction modifying electron–electron collisions. Direct ring-polymer scattering simulations reveal diffractive “softening” of binary deflections, providing a micro-to-macro link between collision physics and diffusion. The present study establishes RPMD as a quantitative, trajectory-based tool for electron transport across the classical–quantum crossover and furnishes benchmarks for improving Coulomb-log interpolation models. We outline extensions to multi-component plasmas and a path to incorporate Fermi–Dirac statistics within path-integral dynamics.

Electronic transport↗

Full‐dimensional coupled‐channel statistical approach to atom‐triatom systems and applications to H/D + O 3 reaction

Abstract The statistical quantum model (SQM), which assumes that the reactivity is controlled by entrance/exit channel quantum capture probabilities, is well suited for chemical reactions with a long‐lived intermediate complex. In this work, a time‐independent coupled‐channel implementation of the SQM approach is developed for atom‐triatom systems in full dimensionality. As SQM treats the capture dynamics quantum mechanically, it is capable of handling quantum effects such as tunneling. A detailed study of the H/D + O 3 capture dynamics was performed by applying the newly developed SQM method on an accurate global potential energy surface. Agreement with previous ring polymer molecular dynamics (RPMD) results on the same potential energy surface is excellent except for very low temperatures. The SQM results are also in reasonably good agreement with available experimental rate coefficients. The strong H/D kinetic isotope effect underscores the dominant role of quantum tunneling under an entrance channel barrier at low temperatures.

Yang, Dongzheng↗

Heterogeneous flexibility can contribute to chromatin segregation in the cell nucleus

The highly and slightly condensed forms of chromatin, heterochromatin and euchromatin, respectively, segregate in the cell nucleus. Heterochromatin is more abundant in the nucleus periphery. Here, in this work, we study the mechanism of heterochromatin segregation by modeling interphase chromosomes as diblock ring copolymers confined in a rigid spherical shell using molecular dynamics simulations. In our model, heterochromatin and euchromatin are distinguished by their bending stiffnesses only, while an interaction potential between the spherical shell and chromatin is used to model lamin-associated proteins. Our simulations indicate that in the absence of attractive interactions between the nuclear shell and the chromatin, most heterochromatin segregates towards the nuclear interior due to the depletion of less flexible heterochromatin segments from the nuclear periphery. This inverted chromatin distribution, which is opposite to the conventional case with heterochromatin dominating at the periphery, is in accord with experimental observations in rod cells. This “inversion” is also found to be independent of the heterochromatin concentration and chromosome number. The chromatin distribution at the periphery found in vivo can be recovered by further increasing the bending stiffness of heterochromatin segments or by turning on attractive interactions between the nuclear shell and heterochromatin. Our results indicate that the bending stiffness of chromatin could be a contributor to chromosome organization along with differential effects of HP ⁢1⁢𝛼-driven phase segregation and of loop extruders and interactions with the nuclear envelope and topological constraints.

biomolecular & subcellular processes↗

Enhancing Thermal Transport in Polymeric Composites Via Engineered Noncovalent Filler–Polymer Interactions

Understanding thermal transport mechanisms in polymeric composites allows us to expand the boundaries of thermal conductivity in them, either increasing it for more efficient heat dissipation or decreasing it for better thermal insulation. But, these mechanisms are not fully understood. Systematic experimental investigations remain limited. Practical strategies to tune the interfacial thermal resistance (ITR) between fillers and polymers and the thermal conductivity of composites remain elusive. Here, we studied the thermal transport in representative polymer composites, using polyethylene (PE) or polyaniline (PANI) as matrices and graphite as fillers. PANI, with aromatic rings in its backbone, interacts with graphite through strong noncovalent π–π stacking interactions, whereas PE lacks such interactions. We can then quantify how π–π stacking interactions between graphite and polymers enhance thermal transport in composites. PE/graphite and PANI/graphite composites with the same 1.5% filler volume fractions show a ∼22.82% and ∼34.85% enhancement in thermal conductivity compared to pure polymers, respectively. Calculated ITRs in PE/graphite and PANI/graphite are ~6 x 10 -8 m 2 KW -1 and ~1 x 10 -8 m 2 KW -1 , respectively, highlighting how π–π stacking interactions reduce ITR. Molecular dynamics (MD) simulations suggest that π–π stacking interactions between PANI chains and graphite surfaces enhance alignment of PANI's aromatic rings with graphite surfaces. This allows more carbon atoms from PANI chains to interact with graphite surfaces at a shorter distance compared to PE chains. Finally, our work indicates that tuning the π–π stacking interactions between polymers and fillers is an effective approach to reduce the ITR and enhance the thermal conductivity of composites.

atoms↗

Molecular-Level Insights into the Influence of Ionic Liquids on the Structure and Dynamics of Neutral and Charged Polyimides

We use all-atom molecular dynamics simulations to investigate how increasing ionic functionalization influences polyimide (PI) behavior in ionic liquid (IL) solvation environments. Here, we examine three polymer systems with varying charge densities: a neutral polyimide [N−PI] containing imidazole rings and progressively introduce charge through quaternization to create singly [C−PI] + and doubly charged [C−2PI] 2+ variants across a wide range of IL concentration (0−90 wt %). Through comprehensive structural, mechanical, and electrostatic analyses, we reveal that polymer charge density plays a central role in shaping IL organization and interaction with the polymer matrix. At low IL content, charged systems exhibit strong electrostatic complexation, leading to chain compaction, localized ESP environments, and elevated dielectric constants. As IL concentration increases, the effect of the different polymer charge states becomes less significant. Notably, an intermediate composition regime at approximately 50 wt % IL is associated with changes in IL-rich domain connectivity and overall system behavior across all three PI systems.

36 MATERIALS SCIENCE↗

ToPolyAgent: AI agents for coarse-grained bead-spring topological polymer simulations

We introduce ToPolyAgent, a multi-agent AI framework for performing coarse-grained molecular dynamics (MD) simulations of topological polymers through natural language instructions. By integrating large language models (LLMs) with domain-specific computational tools, ToPolyAgent supports both interactive and autonomous simulation workflows across diverse polymer architectures, including linear, ring, brush, and star polymers, as well as dendrimers. The system consists of four LLM-powered agents: a Config Agent for generating initial polymer–solvent configurations, a Simulation Agent for executing LAMMPS-based MD simulations and conformational analyses, a Report Agent for compiling markdown reports, and a Workflow Agent for streamlined autonomous operations. Interactive mode incorporates user feedback loops for iterative refinements, while autonomous mode enables end-to-end task execution from detailed prompts. We demonstrate ToPolyAgent's versatility through case studies involving diverse polymer architectures under varying solvent conditions, thermostats, and simulation lengths. Furthermore, we highlight its potential as a research assistant by directing it to investigate the effect of interaction parameters on the linear polymer conformation, and the influence of grafting density on the persistence length of the brush polymer. By coupling natural language interfaces with rigorous simulation tools, ToPolyAgent lowers barriers to complex computational workflows and advances AI-driven materials discovery in polymer science. It lays the foundation for autonomous and extensible multi-agent scientific research ecosystems.

Ding, Lijie [Oak Ridge National Laboratory (ORNL),↗

Effects of high-pressure hydrogen exposure on filler-elastomer adhesion

Elastomers are known to gain enhanced mechanical properties through compounding with nanosized filler particles such as silica or carbon black. Filler dispersion and filler-polymer interfacial strength are key contributing factors to this improvement. The interfacial strength is critical to part lifetime in pressurized gas sealing applications such as O-rings, where weak binding between the filler particle and polymer matrix can lead to internal void structures. With the aim to build a fundamental understanding of precursors to pressurized hydrogen-induced failure in elastomers, we use all-atom molecular dynamics simulations to study the impact of hydrogen oversaturation on filler-polymer interaction strength. We systematically study the interface between a commonly used elastomer, ethylene-propylene-diene monomer (EPDM) and silica by varying gas concentration, crosslink density, and surface chemistry. Our simulations predict that decompression leads to a localization of excess gas near the interface. In conclusion, we demonstrate that this localized gas can weaken interfacial adhesion and quantify the interaction using thermodynamic approaches.

EPDM↗

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

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

Chemical bonding↗

Evidence of short chains in liquid sulfur

High energy x-ray pair distribution function measurements show the average coordination number of the first shell in liquid sulfur is 1.86 ± 0.04 across the λ-transition, not precisely 2.0 as widely accepted. This indicates that upon melting, liquid sulfur does not comprise solely of S 8 rings but also possesses a significant number of short chains. Intensities of the pre-peak and first diffraction peak of the x-ray structure factor and third peak height of the pair distribution function all show deviations at the λ-transition temperature T λ , associated with the break-up of S 8 rings and the start of oligomer polymerization. A significant number of non-bonded or loosely bonded “interstitial atoms,” with an average coordination number of 0.20 ± 0.005, are also observed in the so-called “forbidden zone” between the first and second shells upon melting. The number of interstitial atoms is found to decrease to a minimum at the λ-transition, but the majority persist into the high temperature polymerized liquid. Furthermore, the existence of short chains and nearby interstitial atoms represent the two main factors required to initiate the S 8 -ring to chain transition, as proposed by recent molecular dynamics simulations.

Chemical bonding↗

Origins of Enhanced Ion Transport in Nanostructured Anion-Conducting Polyelectrolytes

Ion-conducting polymer chemistry and microstructure profoundly impact membrane water uptake and ionic conductivity. Water uptake strongly impacts ionic conductivity; yet excess water uptake compromises ion-exchange membrane mechanical properties. Although nanophase separation has been proposed to overcome this trade-off, it is unclear how polymer backbone architecture governs ionic nanostructure and its subsequent impact on water uptake and conductivity. Here, we integrate experiments and molecular dynamics simulations to elucidate the role of backbone chemistry in governing ionic nanostructure, hydration behavior, and ion transport in anion-conducting polyelectrolytes (ACPs). We systematically investigate hydrocarbon polynorbornene (PNB)-based ACPs with three distinct backbone architectures: vinyl-addition polymerization (VAP), ring-opening metathesis polymerization (ROMP), and hydrogenated ROMP. While maintaining comparable ion exchange capacities (IECs) and identical side-chain chemistry, we isolate the effects of backbone structure. We show that nanophase-separated ionic nanostructures originate in the dry state and evolve upon hydration through heterogeneous water uptake, with water preferentially partitioning into ion-rich domains. This nanophase separation arises from a delicate interplay between ionic segregation propensity and the entropic barrier imposed by backbone stiffness. Specifically, flexible backbones intensify attractive ion–ion interactions by reducing the entropic penalty for backbone deformation, promoting nanophase separation, while rigid backbones suppress ionic nanostructure formation. Nanophase-separated ion domains locally concentrate water upon hydration, which in turn enables the connectivity required for fast transport at lower water concentration values. Furthermore, these findings demonstrate that backbone chemistry can be tuned as a design lever to promote nanophase separation and enhance ion transport without excessive water uptake.

Anions↗

Measuring Topological Constraint Relaxation in Ring-Linear Polymer Blends

Polymers are an effective test bed for studying topological constraints in condensed matter due to a wide array of synthetically available chain topologies. When linear and ring polymers are blended together, emergent rheological properties are observed as the blend can be more viscous than either of the individual components. This emergent behavior arises since ring-linear blends can form long-lived topological constraints as the linear polymers thread the ring polymers. Here, in this work, we demonstrate how the Gauss linking integral can be used to efficiently evaluate the relaxation of topological constraints in ring-linear polymer blends. For majority-linear blends, the relaxation rate of topological constraints depends primarily on reptation of the linear polymers, resulting in the diffusive time τ d,R for rings of length N R blended with linear chains of length N l to scale as τ d,R ~ N$^2_R$N$^{3.4}_L$.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Sub‐5 Ångstrom Porosity Tuning in Calixarene‐Derived Porous Liquids via Supramolecular Complexation Construction

Abstract Sub‐Ångstrom‐level porosity engineering, which is appealing in gas separations, has been demonstrated in solid carbon, polymer, and framework materials but rarely achieved in the liquid phase. In this work, a gas molecular sieving effect in the liquid phase at sub‐5 Ångstrom scale is created via sophisticated porosity tuning in calixarene‐derived porous liquids (PLs). Type II PLs are constructed via supramolecular complexation between the sodium salts of calixarene derivatives and crown ether solvents. The chemical structure variation and assembly behavior of the porous host upon PL construction are monitored by spectroscopy‐, X‐ray‐, and neutron‐scattering techniques. The presence of permanent porosity in calixarene‐derived PLs is verified by pressure swing gas uptake, altered CO 2 physisorption behavior, and molecular simulations. Sub‐5 Ångstrom porosity tuning within the PL phase is achieved by introducing bulky substituted groups on the benzene ring of the calixarene host, which then greatly affects the dynamic motion and transport behavior of CO 2 molecules and the Xe uptake performance. The approach being demonstrated in this work represents a promising pathway to tune and leverage the porosity effect for enhanced gas uptake capacity and selectivity in liquid sorbents.

Li, Errui [Department of Chemistry University of T↗

Sub-5 Ångstrom Porosity Tuning in Calixarene-Derived Porous Liquids via Supramolecular Complexation Construction

Sub-Ångstrom-level porosity engineering, which is appealing in gas separations, has been demonstrated in solid carbon, polymer, and framework materials but rarely achieved in the liquid phase. In this work, a gas molecular sieving effect in the liquid phase at sub-5 Ångstrom scale is created via sophisticated porosity tuning in calixarene-derived porous liquids (PLs). Type II PLs are constructed via supramolecular complexation between the sodium salts of calixarene derivatives and crown ether solvents. The chemical structure variation and assembly behavior of the porous host upon PL construction are monitored by spectroscopy-, X-ray-, and neutron-scattering techniques. The presence of permanent porosity in calixarene-derived PLs is verified by pressure swing gas uptake, altered CO 2 physisorption behavior, and molecular simulations. Sub-5 Ångstrom porosity tuning within the PL phase is achieved by introducing bulky substituted groups on the benzene ring of the calixarene host, which then greatly affects the dynamic motion and transport behavior of CO 2 molecules and the Xe uptake performance. Further, the approach being demonstrated in this work represents a promising pathway to tune and leverage the porosity effect for enhanced gas uptake capacity and selectivity in liquid sorbents.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Polymer Macrocycles: A novel topology to control dynamics of rubbery materials (Final Report)

One of the current challenges in polymer physics and engineering is to understand the behavior of circular macromolecules. Though such molecules have been made by dilute solution ring closure methods in the past, such procedures have been limited to relatively small, unentangled rings as well as to small sample quantities. Initial studies were also hampered by issues surrounding purity of the rings, though this problem has been somewhat overcome by the use of liquid chromatography at critical conditions (LCCC) methods that separate rings and linear chains. Furthermore, in spite of the continued progress a full understanding of ring dynamics remains elusive. The present work has addressed multiple aspects of the behavior of ring molecules by using a novel synthetic route: reversible radical recombination redox polymerization (R3P) to produce poly(3,6-dioxa-1,8-octanedithiol) (PolyDODT) molecules. The synthesis conditions were chosen to produce linear (LDODT), cyclic (RDODT) and linear-cyclic mixtures (LRDODT) in order to investigate dilute solution and rheological responses of the materials. Importantly, the R3P synthesis makes larger quantities and higher molecular weights than do the dilute solution ring closure methods and also has the potential to be scaled to industrial quantities, thus becoming commercializable. We emphasize two important outcomes of the work. The first is that the rheological investigation shows that very large molecular weight the PolyDODT rings permitted the investigation of ring dynamics from glassy towards terminal flow for molecular sizes corresponding to approximately 300 entanglements of the linear counterpart. At the same time, these very large sized rings could be diluted in a solvent in order to reduce the molecular entanglement density to less than unity. The results demonstrate that ring viscosity follows a Rouse-like linear dependence on molecular weight or entanglement number to approximately 15 entanglements. This range of "Rouse-like" response is greater than what was achieved in all prior works where synthesis conditions have limited the molecular sizes achievable. After the entanglement number $Z_w$ increases beyond approximately 15, we find that the viscosity-molecular weight scaling undergoes a rapid change from the Rouse-like $η\sim{Z}^1_w$ to a very strong power-law scaling that is approximately $η\sim{Z}^{5.8}_w$. This aspect of the work strongly suggests that prior investigations on highly pure rings have invariably been performed on unentangled systems simply because ring molecules have very high entanglement thresholds. The results are further supported by plateau modulus measurements that show that in the "Rouse-like" viscosity regime, for the RDODT samples there is a corresponding lack of a rubbery entanglement plateau and once entanglement coupling is seen, the rubbery plateau also appears in the dynamic modulus data. Furthermore, upon appropriate scaling, we show that the unentangled polyDODT rings (i.e., the diluted high molecular weight rings) show virtually identical dynamics (upon appropriate scaling to account for the chain dilutions) to the dynamics of similarly unentangled LCCC fractionated polystyrene rings.

36 MATERIALS SCIENCE↗

Intermediate time sub-diffusion and stress relaxation in ring polymer melts

The slow dynamics of non-concatenated ring melts remains a frontier problem in polymer science with implications for many soft material environments including cellular biophysics. Here, in this work, we report large-scale simulations of model ring melts that analyze the monomer and center-of-mass (CM) mean square displacements (MSD) and stress relaxation function on intermediate time and length scales. The degree of dynamical slowing down is characterized by the maximally sub-diffusive fractional time scaling exponents. The data span an exceptionally wide range of ring degrees of polymerization and stiffnesses and are not successfully organized based on the classic measure linear chain entanglement, N/N e . Rather, we find that the crossover degree of polymerization, N D , based on ring macromolecular caging that successfully allows master curves to be constructed for the long-time CM self-diffusion constant also collapses these temporal dynamic scaling exponents. Different properties display different exponents and exhibit one or two regimes of linear variation with the logarithm of N D / N . A distinct crossover of the CM-MSD and stress relaxation exponents emerges at sufficiently large N or stiffness that is not found for the monomer MSD, indicating a novel form of dynamic decoupling. This crossover aligns with the predicted critical degree of polymerization for transitioning from a weak to strong caging regime, indicative of activated transport. The latter may reflect the emergence of an intermolecular collective contribution to stress in analogy with dense soft colloidal matter. Suggestions are made for future theoretical work to address the rich patterns of behavior discovered.

Anomalous diffusion↗

Simplifying the Chemical Design of Nonfused-Ring Electron Acceptors─Lessons Learned from Thienothiophene and Benzodithiophene Cores

Nonfused-ring electron acceptors (NFREA)s have been proposed as alternatives to fused-ring acceptors in organic photovoltaics (OPV)s. The simpler molecular design of NFREAs results in a reduced synthetic complexity, potentially lowering manufacturing costs, and easier access to desired optoelectronic properties by structural modifications. However, they have not seen the same rise in popularity as fused nonfullerene acceptors (NFA)s due to their lower performances in devices. In this work, we explore structure–property relationships of three NFREA molecules based on thieno[3,2-b]thiophene (TT) and benzo[1,2-b:4,5-b′]dithiophene (BDT) cores and malononitrile functionalized-isatin end group acceptors. We describe their synthetic routes, computational analysis, and photophysical and electronic properties. Here, we investigated the effect of changing the core from TT to BDT on electronic properties of the molecules─a result of raising/lowering of the frontier molecular orbital (FMO) levels. We also explored the effect of branched vs nonbranched alkyl groups on the morphological and blend-capabilities in bulk heterojunction (BHJ) OPV devices. These effects were then reflected in the overall device performances, primarily due to favorable FMO energy alignment between the NFREA and polymer. To investigate the behavior of the excited states of the NFREAs within the active layer of the devices, we performed time-resolved microwave conductivity and transient absorption experiments, which confirmed that the TT molecules exhibited very low conductance while the BDT molecule showed moderate photoconductance. Our combined theoretical, synthetic, device, and charge carrier dynamics studies of these materials reveal the central role that structure–property relationships play in future molecular designs of NFREAs that can perform as well as their fused counterparts.

14 SOLAR ENERGY↗