Search NASA⌕ Search

SEARCH · Search NASA

Results for “polymer simulations, 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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Exploration of Tertiary Structure in Sequence-Defined Polymers Using Molecular Dynamics Simulations

Peptoids are a class of sequence-defined biomimetic polymers with peptide-like backbones and side chains located on backbone nitrogens rather than alpha carbons. These materials demonstrate a strong ability for precise control of single-chain structure, multiunit self-assembly, and macromolecular assembly through careful tuning of sequence due to the diversity of available side chains, although the driving forces behind these assemblies are often not understood. Prior experimental work has shown that linked 15mer peptoids can mimic the protein helical hairpin structure by leveraging the chirality-inducing nature of bulky side chains and hydrophobicity, but there are still gaps in our understanding of the relationship between sequence, stability, and particular secondary or tertiary structure. Here, we present a molecular dynamics (MD) study on the folding behavior of these polymers into hairpins, discussing the differences in structure from sequences with various characteristics in water and acetonitrile, and then compare the handedness preference of common helical motifs between solvents.

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

The Work of Mechanical Degradation in Elongating Polymer Melts

Molecular dynamics simulations are used to study the mechanical degradation of well-entangled polymer melts during uniaxial extensional flow. Simulations measure the transient rise in extensional stresses and relate them to the molecular alignment and scission of chain backbones. Intermolecular entanglements couple chain scission in space and time, making degradation sensitive to deformation history and strain rate in ways not displayed by dilute polymer solutions. The rate of chain scission is nonmonotonic and peaks at strains corresponding to the maximum extensibility of entanglement segments but prior to the full extension of chain backbones. We measure a specific work per scission event w* and decompose it into separate contributions associated with chain alignment, chemical bond breaking, and scission-induced plasticity. We find chain scission in melts requires activating plastic dissipation that is multiple orders of magnitude larger than the chemical work required to break a covalent backbone bond. Our findings underscore the critical need to consider bulk polymer mechanics and rheology in designing efficient mechanical degradation and mechanochemical processes.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Solubilization of PET in binary mixtures of HFIP and DCM

The dissolution of polyethylene terephthalate (PET) is a critical step for a solvent-based process, yet it typically requires highly corrosive or toxic solvents. Here, we investigate the solubilization and conformational behavior of PET in binary mixtures of hexafluoro-2-propanol (HFIP) and dichloromethane (DCM) as a strategy to reduce HFIP usage while maintaining effective dissolution. Small-angle neutron scattering (SANS) measurements reveal that PET remains molecularly dissolved in HFIP/DCM mixtures up to 50 vol% DCM. Analysis of PET chain conformations shows a transition from Gaussian behavior at low HFIP fractions to more swollen chains at intermediate compositions, accompanied by a counter-intuitive minimum in the radius of gyration at 50% HFIP. Complementary SANS measurements of the binary solvents demonstrate that compositional heterogeneity is maximized at this same solvent composition, suggesting a direct coupling between solvent microstructure and polymer dimensions. Molecular dynamics simulations corroborate the experimental findings, revealing solvent domain formation, preferential solvation of PET by HFIP, and a “caging” effect arising from solvent heterogeneity that leads to polymer coil compaction. Together, these results provide molecular-level insight into polymer behavior in mixed solvent systems and establish HFIP/DCM mixtures as a promising, more sustainable solvent platform for the PET post-process.

Arifuzzaman, Arif [ORNL]↗
DOE OSTI2025

Investigating Morphology and Diffusion in Simulations of Precise Anion-Conducting Polymers

Using atomistic molecular dynamics simulations, we investigate the morphology and transport properties of a new class of polymers which are functionalized with quaternary ammonium groups for use as anion exchange membranes. The polymers are precision polyolefins with either a trimethylammonium (p5CNMe3) or a dimethyl-hexyl ammonium (p5CNMe2Hx) pendant group at every fifth carbon along a polyethylene backbone. Simulations are performed at hydration levels of 5, 10, 15, and 20 water molecules per ammonium group. The hydrated polymers form nanoscale, percolated hydrophilic domains (water channels) in the hydrophobic polymer matrix that become wider with increasing water content. Water and hydroxide anion diffusion coefficients also increase with increasing water content. The morphology of the water domains is similar in both polymers, while the diffusion coefficients are somewhat lower in p5CNMe2Hx at fixed water content. Furthermore, the diffusion coefficients in both polymers fall on the same curve as a function of the fractal dimension of the percolated water channels, which appears to be a useful scalar measure of the effects of the nanoscale morphology on water and hydroxide anion transport.

Anion exchange membrane↗
DOE OSTI2025

Flow-Driven Stretch Fluctuations Govern the Nonlinear Viscoelasticity of Elongating Associative Polymer Networks

We use nonequilibrium molecular dynamics simulations to verify recent tube-model predictions that associative polymer networks exhibit broad stretch fluctuations during elongational flow. Simulations further show that these fluctuating dynamics give rise to the rate-dependent extensional viscosity 𝜂 𝐸 measured in filament stretching experiments on H-bonding networks. Simulations model bivalent associative networks with a reactive bead-spring model for varying association strength and extensional strain rate. We observe that stretch fluctuations are driven by a new form of chain tumbling, where chains continually collapse and elongate as their associations break and reform within the advecting network. This produces a broad, nearly uniform distribution of chain stretch over a wide range of strain rates, manifesting as a rate-independent plateau in the extensional stress. Our results show that the nonlinear viscoelasticity of associative networks is dominated by large fluctuations in molecular response, which cannot be captured by current mean-field models.

biomimetic & bio-inspired materials↗
DOE OSTI2026

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),↗
DOE OSTI2026

Computational methods in solution-based plastics purification

Plastic waste can be recycled into resins with near-virgin properties by solution-based purification processes that selectively dissolve polymers, remove contaminants, or detach printing residues. Here, in this review, we examine computational methods for predicting the behavior governing solution-based plastic purification, motivated by the vast polymer–solvent–contaminant compositional space. We discuss thermodynamic and machine learning methods for predicting polymer–solvent and polymer–contaminant interaction and review physics-based molecular dynamics simulations that resolve molecular-scale phenomena within polymer matrices inaccessible to screening methods. We highlight how these methods have informed experimental design for dissolution-based recycling and solvent-based contaminant removal. Finally, we discuss the prospective role of agentic AI in integrating these computational tools with real-time sorting data to adapt purification conditions to the compositional variability of real post-consumer feedstocks. This review charts a path toward computationally guided solution-based purification workflows that can respond to the complexity inherent in plastic waste streams.

Altamimi, Ali [Univ. of Wisconsin, Madison, WI (Un↗
DOE OSTI2025

Patchy nanoparticles by atomic stencilling

Stencilling, in which patterns are created by painting over masks, has ubiquitous applications in art, architecture and manufacturing. Modern, top-down microfabrication methods have succeeded in reducing mask sizes to under 10 nm, enabling ever smaller microdevices as today’s fastest computer chips. Meanwhile, bottom-up masking using chemical bonds or physical interactions has remained largely unexplored, despite its advantages of low cost, solution-processability, scalability and high compatibility with complex, curved and three-dimensional (3D) surfaces. Here we report atomic stencilling to make patchy nanoparticles (NPs), using surface-adsorbed iodide submonolayers to create the mask and ligand-mediated grafted polymers onto unmasked regions as ‘paint’. We use this approach to synthesize more than 20 different types of NP coated with polymer patches in high yield. Polymer scaling theory and molecular dynamics (MD) simulation show that stencilling, along with the interplay of enthalpic and entropic effects of polymers, generates patchy particle morphologies not reported previously. These polymer-patched NPs self-assemble into extended crystals owing to highly uniform patches, including different non-closely packed superlattices. We propose that atomic stencilling opens new avenues in patterning NPs and other substrates at the nanometre length scale, leading to precise control of their chemistry, reactivity and interactions for a wide range of applications, such as targeted delivery, catalysis, microelectronics, integrated metamaterials and tissue engineering.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Shear Response of Ionizable Polymer Melts at the Crossover from Ionomers to Polyelectrolytes

Nonlinear shear response of polymers is affected by inherent barriers for diffusive motion, including entanglements and topology. In ionizable polymers, ionic clusters further constrain the intrinsic dynamics of the polymers, significantly enhancing their viscosity. Here, using fully atomistic molecular dynamics simulations, the nonlinear shear response of ionizable polymers is presented, across the transition from the ionomer regime where distinctive clusters dominate the structure to the polyelectrolyte regime where clusters percolate, in polystyrene randomly sulfonated with fractions of SO 3 − groups of f = 0.20 and 0.35, in pristine and tetrahydrofuran (THF) swollen polyelectrolyte melts. For f = 0.20, the ionic clusters first fracture into smaller clusters followed by splitting into individual ionic groups and eventually reform. At higher f, the clusters morph in shape but do not break under high shear. At very high shear rates, all of the chains stretch and recoil rapidly. As the shear rate is reduced, some chains stretch and recoil, while others remain largely unaffected by the shear. Macroscopically, for all systems, the shear viscosity displays initially an elastic response, followed by nonlinear shear stress overshoot and, eventually, a steady state. The evolution of viscosity with time and shear reflects that of the ionizable domains.

cluster chemistry↗
DOE OSTI2024

Effect of stoichiometry and hydration level on water domain size and transport in poly(aryl piperidinium) alkaline anion-exchange membranes

Alkaline water electrolysis holds promise in decarbonizing the global economy by enabling renewable hydrogen production with non-precious group metal catalysts. Anion exchange membranes are an important component of alkaline water electrolyzers and would ideally be durable while allowing for high hydroxide conductivity. The poly(aryl piperidinium) (PAP) class of polymers has attracted recent interest due to their good mechanical robustness and high ionic conductivity. Here, in this work, we perform atomistic molecular dynamics (MD) simulations of several PAP polymers at experimentally relevant hydration levels and polymer ion exchange capacities (IECs) to gain nanoscale insight into their properties and to help elucidate the trade-offs that result from tuning the IECs through the polymer stoichiometry. Our MD-predicted macroscopic polymer properties were found to be in good agreement with experimentally available polymer swelling ratios, water-occupied volumes, X-ray scattering, and ionic conductivities. The models show that for hydration levels greater than 8H 2 O per cation a single water cluster will form that percolates through the system. The growth in water cluster size results in large polymer swelling, the creation of larger channels with widths of 7 Å or larger, and nanophase separation between the hydrophilic domains and the polymer with characteristic length scales of approximately 20–30 Å. The experimentally observed lack of a strong X-ray scattering peak at low wavevectors can be explained by a cancellation between the polymer-polymer/water-water and polymer-water correlations and not a loss in nanophase separation. The overlap in coordination environments of the hydroxide oxygen and polymer nitrogen atoms implies that vehicular diffusion between cationic groups could play a role in hydroxide transport. The polymers' hydroxide and water diffusion constants increase by approximately an order of magnitude between hydration levels of 8 and 20H 2 O per cation. However, there are diminishing returns in hydroxide diffusion constant once the IEC exceeds 2.4 meq/g.

08 HYDROGEN↗
DOE OSTI2026

Understanding the Influence of Chain Architecture on the Transport Quantities of Polymer Electrolytes with Covalently Bonded Anions

Here, we use a combination of experiments and coarse-grained molecular dynamics simulations to elucidate the structure–property relationships in polymer electrolytes obtained by the copolymerization of poly(vinyl ethylene carbonate─lithium styrene bis(trifluoromethanesulfonyl)imide) or p(VEC-LiSTFSI). Experiments show that the conductivity reduces with increasing anion (i.e., STFSI) fraction on the chain, and the cation transference number (t + ) is found to be dependent on the anion fraction. Furthermore, a significant fraction of unpolymerized VEC monomers are observed. Since it is inherently difficult to experimentally control the chain architecture and the amount of unpolymerized VEC in these systems, we perform coarse-grained molecular dynamics simulations on model polymer systems with different chain architectures to mimic the plausible experimental systems. Specifically, we look at the differences in transference numbers arising from (i) a random copolymer of VEC and STFSI monomers; (ii) a blend of VEC-STFSI copolymer with VEC monomers; and (iii) a ternary blend of the VEC homopolymer, STFSI homopolymers, and VEC monomers. The ternary blend model demonstrates the closest resemblance with the experimental transference numbers and diffusivities. The lithium diffusivity obtained from the coarse-grained models with VEC monomers (plasticizers) is about 1.5 times that of the model without VEC monomers, showing that the plasticizing effect of VEC monomers is modest. We rationalize the experimental observations based on aggregate and cluster analyses obtained from molecular simulations. This work reveals that polymer electrolyte chain architecture and plasticizers can critically influence the transport properties, and these parameters should be considered when designing single ion conducting polymeric electrolytes.

cluster distribution↗
DOE OSTI2025

Probing Nanorod Assembly and Dynamics in Polymer Nanocomposites in Equilibrium and Shear

Coarse-grained molecular dynamics simulations are used to examine the structure and dynamics of nanorod assemblies in polymer melts under equilibrium and simple shear. We show that as the concentration of nanorods increases, there is a transition from an isotropic phase to a two-phase region in which the nanorods phase separate into a dilute phase and dense bundles of hexagonally packed nanorods. The onset of the two-phase region is below that predicted by Onsager theory, which we attribute to an effective increase in the diameter of nanorods due to a layer of polymer bound to the rod surfaces. Equilibrium simulations show that increasing polymer chain length N enhances nanorod bundling at fixed nanorod concentration. Simulations of systems undergoing simple shear show that flow enhances nanorod alignment and bundling relative to those of equilibrium systems with similar properties. Finally, simulations reveal that increasing nanorod length enhances nanorod alignment under shear at equivalent shear rates but not at equivalent Péclet numbers. Overall, our simulations highlight that polymer matrix-nanorod attraction (i.e., chemistry), shear rate, and matrix chain length are desirable design variables to control the structure of nanorod-containing soft materials under simple shear.

diffusion↗
DOE OSTI2024

A transferable classical force field to describe glyme based lithium solvate ionic liquids

A non-polarizable force field for lithium (Li + ) and bis(trifluoromethanesulfonyl)imide (TFSI – ) ions solvated in diglyme at around 0.2 mol fraction salt concentration was developed based on ab initio molecular dynamics (AIMD) simulations and a modified polymer consistent force field model. A force–torque matching based scheme, in conjunction with a genetic algorithm, was used to determine the Lennard-Jones (LJ) parameters of the ion–ion and ion–solvent interactions. This force field includes a partial charge scaling factor and a scaling factor for the 1–4 interactions. The resulting force field successfully reproduces the radial distribution function of the AIMD simulations and shows better agreement compared to the unmodified force field. The new force field was then used to simulate salt solutions with glymes of increasing chain lengths and different salt concentrations. The comparison of the MD simulations, using the new force field, with experimental data at different salt concentrations and AIMD simulations on equimolar concentrations of the triglyme system demonstrates the transferability of the force field parameters to longer glymes and higher salt concentrations. Furthermore, the force field appears to reproduce the features of the experimental x-ray structure factors, suggesting accuracy beyond the first solvation shell, for equimolar salt solutions using both triglyme and tetraglyme as the solvent. Altogether, the new force field was found to accurately reproduce the molecular descriptions of LiTFSI-glyme systems not only at various salt concentrations but also with glymes of different chain lengths. Thus, the new force field provides a useful and accurate tool to perform in silico studies of this family of systems at the atomistic level.

25 ENERGY STORAGE↗
DOE OSTI2026

Molecular insights into CO 2 -to-bicarbonate transformation in functionalized anion exchange ionomers for electrochemical separations

Bipolar membrane (BPM) electrochemical processes are a promising platform for carbon dioxide (CO 2 ) separations, but the molecular level thermodynamic and kinetic understanding of CO 2 -to-bicarbonate (HCO 3 − ) transformation remain poorly understood. This study employs a multiscale computational approach to systematically explore the adsorption and reactive transformation of CO 2 in five anion exchange ionomer systems. Classical molecular dynamics (MD) simulation results demonstrate that polymers with imidazolium groups significantly reduce CO 2 diffusion and enhance (OH − )–CO 2 interactions due to stronger electrostatic and π-interactions. Compared to the commonly used quaternary ammonium ionomers, imidazolium-functionalized ionomers show improved CO 2 proximity and interaction strength. Ab initio MD and density functional theory (DFT) calculations reveal that the benzyl-substituted imidazolium (IM-Ben) substantially reduces the energy barrier for HCO 3 − formation (∼72 meV lower) compared to the alkyl-substituted IM-nBu, while also mitigating imidazolium deprotonation under moderate hydration conditions. Transition state analysis shows IM-Ben forms more extensive hydrogen-bonding networks, which stabilize the transition state structure and contribute to a lower energy barrier for bicarbonate formation. These findings highlight the advantage of the adjacent benzyl moiety in enabling efficient CO 2 -to-bicarbonate transformation via hydrated hydroxide ion counterions, offering mechanistic insights and clear molecular design principles for optimizing anion exchange ionomers at bipolar membrane interfaces for electrochemical CO 2 separation applications.

Bipolar membranes, Reactive transformation of CO2,↗
DOE OSTI2025

High-voltage water-scarce hydrogel electrolytes enable mechanically safe stretchable Li-ion batteries

Soft Li-ion batteries, based on conventional organic electrolytes, face performance degradation challenges due to moisture penetration and safety concerns due to possible leakage of toxic fluorine compounds and flammable solvents under mechanical damage. We design a water-scarce hydrogel electrolyte with fluorine-free lithium salt to achieve wide electrochemical stability window (up to 3.11 volts) in ambient air without hermetic packaging while balancing high stretchability (1348%), ion conductivity (41 millisiemens per centimeter), and self-healing capabilities for mechanically and chemically safe stretchable Li-ion batteries. Molecular synergy between hydrophilicity and lithiophilicity of zwitterionic polymer backbone is revealed by molecular dynamics simulations. The battery exhibits capacity retention under harsh mechanical stresses—enduring stretching, twisting, folding, and multiple through-punctures by a needle—while self-healing from repeated through cuts by a razor blade. Stable ambient operation for 1 month over 500 charge-discharge cycles (average coulomb efficiency, 95%) is achieved. A prototype self-healing electronic system with embedded soft batteries demonstrates practical application as a durable embodied energy source.

Science & Technology - Other Topics↗
DOE OSTI2026

Specific ion effects on ion transport in charged polymer membranes

Selective separation of like-charged ions is a central challenge in applications such as critical mineral recovery. Electrochemical membrane–based separations offer promising pathways to address this need, but a limited fundamental understanding of ion transport in charged polymer membranes hampers the development of highly selective materials. This study elucidates the role of specific ion effects (SIEs) in ion transport through charged polymer membranes by integrating experimental measurements of ion mobility and in situ ion/ion and ion/water interactions in model charged polymer membranes with results from molecular dynamics simulations on analogous systems. We demonstrate that solvent-mediated ion interactions drive pronounced SIEs within the studied membranes, with the ion softness, i.e., the malleability of ion hydration shells, from hard/soft acid/base (HSAB) theory emerging as a key predictor of transport properties. HSAB theory also explains the observed mechanism of solvent-mediated ion interactions. Our findings offer a mechanistic framework for designing membranes with tailored ion selectivity, potentially enabling efficient separations of chemically similar ions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗