Search NASA⌕ Search

SEARCH · Search NASA

Results for “molecular dynamics of polymers”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
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 OSTI2022

Structure/property relationship of semicrystalline polymers during tensile deformation: a molecular dynamics approach

Molecular dynamics (MD) simulation with an ab initio numerical model is conducted to study the structure–property relationship of semicrystalline polymers during uniaxial deformation. The effects of chain length, temperature, and strain rate on mechanical properties are discussed. The influences of microstructural evolution such as bond length, entanglement density, and chain orientation are also studied quantitatively. The temperature will greatly affect the chain conformation in amorphous domains, and the results revealed that the interaction of amorphous and crystalline domains played a crucial role during stretching. Accordingly, the yielding of semicrystalline polymers follows different mechanisms at temperatures above and below T g . A melt-recrystallization scheme is observed during yielding at higher temperatures, while destruction of crystal structures is observed at lower temperatures at the yield point. The correlated effects of different temperatures and strain rates on mechanical properties are examined. This research is part of our efforts to develop a digital twin of a real experiment for efficient optimization of polymer material properties.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Dimension-free path-integral molecular dynamics without preconditioning

Convergence with respect to imaginary-time discretization (i.e., the number of ring-polymer beads) is an essential part of any path-integral-based molecular dynamics (MD) calculation. However, an unfortunate property of existing non-preconditioned numerical integration schemes for path-integral molecular dynamics—including essentially all existing ring-polymer molecular dynamics (RPMD) and thermostatted RPMD (T-RPMD) methods—is that for a given MD time step, the overlap between the exact ring-polymer Boltzmann–Gibbs distribution and that sampled using MD becomes zero in the infinite-bead limit. This has clear implications for hybrid Metropolis Monte Carlo/MD sampling schemes, and it also causes the divergence with bead number of the primitive path-integral kinetic-energy expectation value when using standard RPMD or T-RPMD. We show that these and other problems can be avoided through the introduction of “dimension-free” numerical integration schemes for which the sampled ring-polymer position distribution has non-zero overlap with the exact distribution in the infinite-bead limit for the case of a harmonic potential. Most notably, we introduce the BCOCB integration scheme, which achieves dimension freedom via a particular symmetric splitting of the integration time step and a novel implementation of the Cayley modification [R. Korol et al., J. Chem. Phys. 151, 124103 (2019)] for the free ring-polymer half-steps. More generally, we show that dimension freedom can be achieved via mollification of the forces from the external physical potential. The dimension-free path-integral numerical integration schemes introduced here yield finite error bounds for a given MD time step, even as the number of beads is taken to infinity; these conclusions are proven for the case of a harmonic potential and borne out numerically for anharmonic systems that include liquid water. The numerical results for BCOCB are particularly striking, allowing for nearly three-fold increases in the stable time step for liquid water with respect to the Bussi–Parrinello (OBABO) and Leimkuhler (BAOAB) integrators, while introducing negligible errors in the calculated statistical properties and absorption spectrum. Importantly, the dimension-free, non-preconditioned integration schemes introduced here preserve ergodicity and global second-order accuracy, and they remain simple, black-box methods that avoid additional computational costs, tunable parameters, or system-specific implementations.

Korol, Roman (ORCID:0000000193076351)↗
DOE OSTI2026

Structure and Flow-Viscosity of Filled-Polymer-Based 3D Printing Ink: Exploration through Coarse-Grained Molecular Dynamics

The addition of nanofiller particles to a polymer matrix has long been known to enhance or modify the composite’s mechanical and rheological properties. However, quantitatively capturing such changes with molecular level simulations remains computationally challenging. Toward that goal, we performed coarse-grained molecular dynamics of a nanocomposite system at a fixed (25 vol %) filler loading under nonspecific, weak polymer–filler interactions representative of a broad class of technologically important materials. We report several interesting results, including: (1) the equilibrium chain-configuration remains Gaussian-like as in an unfilled melt; (2) smaller filler particles display a stronger tendency to cluster; (3) larger fillers act as plasticizers by reducing the entanglement density and accelerating the chain mobility; and (4) fillers enhance the tensile response modulus, with the effect being stronger for larger particles. We also simulate cluster breakup, yielding, and elongational flow under an applied time-linear tensile strain and study the flow viscosity as a function of filler-size and chain-length.

Materials science↗
DOE OSTI2024

Coupling and Decoupling between Stickers and Backbones in Associating Polymers with Terminally Functionalized Side Chains

Associating polymers have recently attracted tremendous attention as self-healing polymers due to the dynamic reversibility of sticker association–dissociation processes. These sticker dynamics inevitably slow significantly as segmental dynamics decelerate below the glass transition temperature that limits the operational temperatures of associating polymers. Using broadband dielectric spectroscopy on a new set of associating polymers composed side chains of various lengths and terminal carboxylic-acid sticker functionality, we discovered an unexpected decoupling near the glass transition (T g ) between the slower backbone and fast sticker motions. These two motions have distinct activation energies for the motion of the sulfur linkage adjacent to the backbone (~22 to 26 kJ/mol) and the carboxylic acid groups (~5 to 8 kJ/mol). A second decoupling appears in the polymer with the longest side chain length at T g – 20 K, reflecting the dynamics of various carboxylic acid binding configurations. Furthermore, our findings open the door to design associating polymers with the capability of sticker association/dissociation at low temperatures.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Understanding Interfacial Block Copolymer Structure and Dynamics

Here, block copolymer (BCP) structure and dynamics were studied using small-angle neutron scattering (SANS), neutron spin echo (NSE) spectroscopy, and molecular dynamics (MD) simulations to obtain a fundamental understanding of the impact of an interfacial block on chain dynamics. A glassy block acted as the interface, and the dynamics of a rubbery block was studied. The rubbery block was protonated near the interface in one sample and near the chain end in another sample to observe the interfacial effect on the rubbery polymer. Analysis of the structure and dynamics revealed that the interfacial rubbery block was confined in layered morphologies and exhibited much slower dynamics than the chain-end rubbery block that was dispersed in the rubbery matrix. The interfacial rubbery block showed weaker dynamical relaxation than that at the chain end, and it also had critically important length scale dependence. Dynamical slowing was only observed at length scales significantly larger than the characteristic segmental length, and the disparity between interfacial and chain-end dynamics increased with increasing length.

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

Segmental Dynamics Measured by Quasi-Elastic Neutron Scattering and Ion Transport in Chemically Distinct Polymer Electrolytes

We investigate the segmental dynamics and ion transport in two chemically distinct polymer electrolytes, poly(2-cyanoethyl acrylate) (PCEA) and poly(ethylene oxide) (PEO), and their mixtures with lithium bis(trifluoromethane)sulfonimide (LiTFSI) salt. Quasi-elastic neutron scattering experiments reveal slow dynamics in PCEA/LiTFSI relative to that in PEO/LiTFSI, translating to monomeric friction coefficients that are orders of magnitude different. In spite of the enhanced salt dissociation in PCEA due to the presence of polar groups, ion transport is largely dominated by the effect of increased monomeric friction in the pure polymer. Conductivity in these systems is quantified through a simple expression combining salt dissociation, the monomeric friction in the pure polymer, and the effect of added salt on the monomeric friction.

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

Relaxation dynamics of water in the vicinity of cellulose nanocrystals

Water molecules near cellulose nanocrystals (CNCs; produced via the sulfuric acid-catalyzed hydrolysis of wood pulp) are believed to relax slower than those in the bulk liquid, which may result in unique properties of CNC aqueous dispersions. Herein this study analyzed the polarization behavior of water molecules in CNC aqueous dispersions and other reference samples using a dielectric relaxation spectroscopy (DRS) technique in the microwave frequency range (0.2–20 GHz). As the CNC concentration increases, two slow relaxation components become prominent. The comparison with DRS data of aqueous dispersions of nanoporous silica, polyvinyl alcohol (PVA), and hairy CNC (HCNC) with amorphous chains protruding from both ends suggested that these slow relaxation modes of water near CNC surfaces cannot be attributed to direct hydrogen bonding interactions with the hydroxyl (OH) groups exposed and immobilized at the solid surface. Instead, they are similar to the water molecules interacting with OH groups attached to flexible polymer chains. Molecular dynamics (MD) simulations of the polarization behavior of water near the (110) facet of cellulose Iβ crystals confirmed that the interactions of water molecules with the cellulose crystal surface do not cause slower relaxations in the frequency range studied via the DRS. These results indicated that the CNC surface cannot be depicted with the crystallographic facets of cellulose Iβ; instead, it resembles a polymer-brush surface on which the short glucan residues or fragments of the strong acid-catalyzed hydrolysis process are swollen and extended into the aqueous phase.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Quantum Simulations of Radiation Damage in a Molecular Polyethylene Analog

Abstract An atomic‐level understanding of radiation‐induced damage in simple polymers like polyethylene is essential for determining how these chemical changes can alter the physical and mechanical properties of important technological materials such as plastics. Ensembles of quantum simulations of radiation damage in a polyethylene analog are performed using the Density Functional Tight Binding method to help bind its radiolysis and subsequent degradation as a function of radiation dose. Chemical degradation products are categorized with a graph theory approach, and occurrence rates of unsaturated carbon bond formation, crosslinking, cycle formation, chain scission reactions, and out‐gassing products are computed. Statistical correlations between product pairs show significant correlations between chain scission reactions, unsaturated carbon bond formation, and out‐gassing products, though these correlations decrease with increasing atom recoil energy. The results present relatively simple chemical descriptors as possible indications of network rearrangements in the middle range of excitation energies. Ultimately, the work provides a computational framework for determining the coupling between nonequilibrium chemistry in polymers and potential changes to macro‐scale properties that can aid in the interpretation of future radiation damage experiments on plastic materials.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2026

Molecular Dynamics Investigation of Oil Wetting on Synthetic Polymer Substrates

Understanding the interaction of polymer surfaces with nonpolar, low surface tension liquids, or whether a substrate is oleophobic versus oleophilic, is critical for applications ranging from antifingerprint coatings to oil–water separation membranes and oil spill remediation. Despite its technological importance, the molecular mechanisms underpinning polymer–oil wetting are not well characterized. Here, we employ molecular dynamics simulations to investigate the behavior of n-hexadecane in contact with chemically distinct polymer surfaces, spanning eight constitutional unit chemistries as well as amorphous and crystalline morphologies. This permits a critical examination of both thermodynamic and dynamic descriptors of oleophobicity, including oil contact angle, dewetting free energy, interfacial diffusivity, and a proposed “ghost probe energy,” which does not require explicit simulation of oil–polymer interactions. We find that the oil contact angle does not reliably distinguish oleophobic behavior across polymer chemistries, whereas other metrics provide clearer and more consistent differentiation. Analysis of the results reveals that polymer–oil wetting behavior is primarily governed by interfacial van der Waals interactions and modulated by surface flexibility and morphology. Collectively, this work establishes a computational framework for characterizing oil wetting and provides additional insight into what molecular-level factors dictate trends in oleophobicity.

lipids↗
DOE OSTI2026

Associative polymers with controlled sticker placement: How reversible bond distribution and density govern polymer dynamics

Associative polymers with precisely arranged stickers offer opportunities to program material properties with molecular precision. Yet, it remains unclear how the placement and fraction of stickers dictate structure, dynamics, and macroscopic properties. By developing a model unentangled polymer system with hydrogen-bonding stickers, we show that randomly distributed stickers neither form clusters nor change flow properties, whereas stickers placed at chain ends drive nanocluster formation even at low concentrations. Adding more end stickers produces a rubbery plateau spanning eight decades in frequency with two distinct relaxation timescales, in contrast to the single plateau predicted by the classic sticky Rouse model. These results demonstrate that sticker distribution dictates whether associative polymers undergo nanocluster formation or microphase separation, while substantial alterations in dynamics and viscoelasticity require both sticker aggregation and thermomechanical stability of associated domains. Our findings resolve a longstanding debate on associative polymer dynamics and provide molecular design rules for programmable soft materials.

36 MATERIALS SCIENCE↗
DOE OSTI2024

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↗
DOE OSTI2022

Microscopic Dynamics of Inverse Wormlike Micelles Probed Using X-ray Photon Correlation Spectroscopy

Wormlike micelles (WLMs) are ubiquitous viscoelastic modifiers that share properties with polymer solutions. While 5 their macroscopic rheology is well-understood, their microscopic dynamics are less studied because they span a large range of time 6 and length scales. In this work, we demonstrate the use of X-ray photon correlation spectroscopy to interrogate the segmental 7 dynamics of inverse WLM solutions swollen with a rubidium chloride solution. We observe a diffusive scaling of the dynamics and 8 extract a temperature-dependent diffusion coefficient, which we associate with the thermal interactions of the slow segmental 9 dynamics near entanglement points. Further, we probe this relaxation process across the unbranched to branched topological transition and 10 find no microstructural evidence of branch formation in the slow mode. Instead, we observe that the dynamics become more 11 homogeneous and prominent as the temperature is reduced and water content increases.

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

Molecular dynamics simulations of energy accommodation between gases and polymers for ultra-low thermal conductivity insulation

Determining the energy accommodation between gases and solids is essential to developing porous thermal insulation materials with ultra-low effective thermal conductivity that reduce energy use, greenhouse gas emissions, and fossil fuel consumption. The energy accommodation coefficients of most gases, however, have been rarely studied, especially with respect to solids that have relatively high thermal resistivity, e.g., polymers. Here, by using all-atom nonequilibrium molecular dynamics simulations, we reveal the accommodation coefficients of He, Ar, N 2 , and O 2 with polymers, mainly polystyrene. We find that their values are around 0.51, 0.72, 0.79, and 0.90, respectively, suggesting a critical reexamination of the commonly used theoretical maximum value of 1. We have also conducted experiments and validated the value for air, which is about 0.81. Such a change in accommodation coefficients can lead to a reduction of about 70%, 50%, 35%, and 20% in the thermal conductivity of He, Ar, N 2 , and O 2 gases in nano pores (below 100 nm) or at low pressures (below 1 millibar). With these new accommodation coefficients, we find that in a 10 nm pore with ambient pressure at 300 K, the gas thermal conductivity of He, Ar, N 2 , and O 2 in porous polystyrene can be as low as 9.7 × 10 –4 , 3.4 × 10 –4 , 7.3 × 10 –4 , and 8.5 × 10 –4 W·m –1 ·K –1 , respectively, which are two to three orders of magnitude lower than their bulk values, promising higher thermal resistivity of insulation materials. This work reveals the fundamental energy exchange between gases and polymers, providing important guidance for designing high-performance thermal insulation materials for various applications.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Fingerprinting Brownian Motions of Polymers under Flow

Here, we present a quantitative approach to the self-dynamics of polymers under steady flow by employing a set of complementary reference frames and extending the spherical harmonic expansion technique to dynamic density correlations. Application of this method to nonequilibrium molecular dynamics simulations of polymer melts reveals a number of universal features. For both unentangled and entangled melts, the center-of-mass motions in the flow frame are described by superdiffusive, anisotropic Gaussian distributions, whereas the isotropic component of monomer self-dynamics in the center-of-mass frame is strongly suppressed. Spatial correlation analysis shows that the heterogeneity of monomer self-dynamics increases significantly under flow.

36 MATERIALS SCIENCE↗
DOE OSTI2026

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↗
DOE OSTI2022

Molecular insights into the structure-property relationships of 3D printed polyamide reverse-osmosis membrane for desalination

3D-printing is an emerging method for manufacturing polyamide (PA) reserve osmosis (RO) membranes for water treatment and desalination, which can precisely control membrane structural properties, such as thickness, roughness, and resolution. However, the synthesis-structure (i.e., degree of cross-linking (DC), m-phenylenediamine/trimesoyl chloride (MPD/TMC) ratio, and membrane thickness) to property (permeability and water-salt selectivity) relationships for these membranes has not been well understood. At the same time, a microscopic understanding of the physical mechanism of water and salt transport is needed to guide the design of high-performance 3D-printed membranes and improve the printing efficiency. Furthermore, the atomic-scale transport features and energetics of water and salt ions are studied at high pressure for the 3D-printed PA RO membranes with the different DCs and MPD/TMC ratios through non-equilibrium molecular dynamics (NEMD) simulations. Factoring in membrane structure properties, rejection ratio of salt ions and pressure-dependent water flux, 3D-printed PA membranes having an MPD/TMC ratio of 3.0:2.0 and a DC between 80%~90% attains ideal performance: high water flux, high rejection of salt ions, and excellent structural integrity. Mechanistically, water permeability for highly cross-linked PA RO membranes depends on the temporary on-and-off channels that allow water molecules to jump from one cavity to another at high pressure. In addition, higher pressures cause rapid compaction of PA membranes’ free volume and membrane thickness. Membrane failure at high pressure is determined by the DC and MPD/TMC ratios-dependent compressive yield strength. In short, these findings provide physical insights for optimizing existing PA membranes and designing next-generation desalination membranes at the molecular level.

3D-printed PA RO membrane↗