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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Effect of Zwitterionic Additives on Solvation and Transport of Sodium and Potassium Cations in (Ethylene Oxide)10: A Molecular Dynamics Simulation Study

Sodium- (Na + ) and potassium- (K + ) ion batteries are cost-effective alternatives to lithium-ion (Li + ) batteries due to the abundant sodium and potassium resources. Solid polymer electrolytes (SPEs) are essential for safer and more efficient Na + and K + batteries because they often exhibit low ionic conductivity at room temperature. While zwitterionic (ZW) materials enhance Li+ battery conductivity, their potential for Na + and K + transport in batteries remains unexplored. In this study, we investigated the effect of three ZW molecules (ChoPO4, i.e., 2-methacryloyloxyethyl phosphorylcholine, ImSO3, i.e., sulfobetaine ethylimidazole, and ImCO2, i.e., carboxybetaine ethylimidazole) on the dissociation of Na + and K + coordination with ethylene oxide (EO) chains in EO-based electrolytes through molecular dynamics simulations. Our results showed that ChoPO4 possessed the highest cation–EO10 dissociation ability, while ImSO3 exhibited the lowest. Such dissociation ability correlated with the cation–ZW molecule coordination strength: ChoPO4 and ImSO3 showed the strongest and the weakest coordination with cations. However, the cation–ZW molecule coordination could slow the cationic diffusion. The competition of these effects resulted in accelerating or decelerating cationic diffusion. Our simulated results showed that ImCO2 enhanced Na + diffusion by 20%, while ChoPO4 and ImSO3 led to a 10% reduction. For K + , ChoPO4 reduced its diffusion by 40%, while ImCO2 and ImSO3 caused a similar decrease of 15%. These findings suggest that the ZW structure and the cationic size play an important role in the ionic dissociation effect of ZW materials.

25 ENERGY STORAGE↗
DOE OSTI2020

Effects of Ion Size and Dielectric Constant on Ion Transport and Transference Number in Polymer Electrolytes

Salt-doped polymers without added solvents have potential as nonflammable, robust electrolytes for batteries but suffer from low ion conductivity. Using bulky anions with delocalized charge may reduce ion agglomeration and increase conduction. However, size asymmetry between ions may increase preferential solvation of cations versus the larger anions, lowering the transference number t + (fraction of conductivity contributed by the cation, relevant to battery performance). We use coarse-grained molecular dynamics simulations, including a 1/r 4 potential to capture size-dependent solvation effects, to relate polymer and ion chemistry to t + and overall conductivity. At low ion size disparity, increasing dielectric constant improves cation conduction due to the mitigation of ion aggregation and correlated cation-anion motion. However, at high ion size disparity, high dielectric medium results in strong preferential solvation of cations, reducing cation mobility and t + . The trade-off between better cation-anion separation and cation preferential solvation suggests that the strategy of developing high dielectric polymers may enhance performance only at low ion size disparity.

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

Effect of reaction media on hydrogenolysis of polyethylene plastic waste: Polymer-surface interactions in small alkane/polymer blends

The polymer reaction media and its properties can be altered by recycling a fraction of liquid products or adding alkane solvents. Less clear is whether this strategy affects hydrogenolysis. Herein, we investigated the effect of short-chain alkanes C n consisting of n carbons (n=8, 16, and 32) on the upcycling of high-density polyethylene (HDPE) plastic waste to lubricant-range products over Ru/TiO 2 catalysts by multiscale simulations and experiments. First, we trained a force field for polymer/surface interactions on a Ru 22 nanoparticle (NP) supported on TiO 2 . Using replica exchange molecular dynamics simulations, we studied the effect of small hydrocarbons on the adsorption of a surrogate polymer, C 142 , on the catalyst. We found segregation of long chains (C 142 ) at the catalyst surface due to the enthalpy gained by adsorbing more C-C bonds of the long chains, compensating for entropic losses upon adsorption. Short-chain molecules decrease the adsorbed carbons of long chains on the Ru NP due to blocking Ru active sites. Compared to the bulk chains, competitive adsorption results in a broader, heavy-tailed distribution of end-to-end distance of adsorbed chains. Our experiments demonstrated that catalyst activity declines significantly beyond simple dilution due to changes in polymer adsorption, and tuning the reaction media by creating suitable blends impacts hydrogenolysis. Density distributions for a 50:50%wt mixture of PP and PE show that PE chains are segregated at the surface, so they are prone to C-C bond breaking much faster than PP chains. H/D exchange experiments show preferential deuteration of PE, while CH 3 groups of PP remain undeuterated. Furthermore, this may be explained by the preferential sorption of PE over PP, leading to specific distribution in the polymer blend.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2018

Multiscale Modeling of PEEK Using Reactive Molecular Dynamics Modeling and Micromechanics

Polyether ether ketone (PEEK) is a high-performance, semi-crystalline thermoplastic that is used in a wide range of engineering applications, including some structural components of aircraft. The design of new PEEK-based materials requires a precise understanding of the multiscale structure and behavior of semi-crystalline PEEK. Molecular Dynamics (MD) modeling can efficiently predict bulk-level properties of single phase polymers, and micromechanics can be used to homogenize those phases based on the overall polymer microstructure. In this study, MD modeling was used to predict the mechanical properties of the amorphous and crystalline phases of PEEK. The hierarchical microstructure of PEEK, which combines the aforementioned phases, was modeled using a multiscale modeling approach facilitated by NASA's MSGMC. The bulk mechanical properties of semi-crystalline PEEK predicted using MD modeling and MSGMC agree well with vendor data, thus validating the multiscale modeling approach.

Pisani, William A.↗
DOE OSTI2021

Ultrahigh thermal conductivity in three-dimensional covalent organic frameworks

Covalent organic frameworks (COFs) are a new, growing class of crystalline polymers with exceptional characteristics such as low density, mechanical strength, and flexibility. Taking three-dimensional (3D) COF-300 derivatives as an example, we apply molecular dynamics simulations to show that ultrahigh thermal conductivity (>15 W/(mK)) is achieved at a pore size of d = 0.63 nm. This is the first observation of such an outstanding intrinsic thermal conductivity of a 3D polymer in all directions. Further temperature-dependent thermal conductivity calculations and X-ray diffraction (XRD) patterns highlight the importance of maintaining microscopic order for a large thermal conductivity. We also demonstrate that the tunability and anisotropy in COFs provide sufficient room to engineer and obtain desired thermal conductivity. These findings provide a fundamental understanding of structure-thermal conductivity relation in COFs and shine a light on high-thermal-conductivity COFs for thermal management applications.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Machine learning prediction of glass transition temperature of conjugated polymers from chemical structure

Predicting the glass transition temperature (T g ) is of critical importance as it governs the thermomechanical performance of conjugated polymers (CPs). Here, we report a predictive modeling framework to predict T g of CPs through the integration of machine learning (ML), molecular dynamics (MD) simulations, and experiments. With 154 T g data collected, an ML model is developed by taking simplified “geometry” of six chemical building blocks as molecular features, where side-chain fraction, isolated rings, fused rings, and bridged rings features are identified as the dominant ones for T g . MD simulations further unravel the fundamental roles of those chemical building blocks in dynamical heterogeneity and local mobility of CPs at a molecular level. The developed ML model is demonstrated for its capability of predicting T g of several new high-performance solar cell materials to a good approximation. The established predictive framework facilitates the design and prediction of T g of complex CPs, paving the way for addressing device stability issues that have hampered the field from developing stable organic electronics.

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

Evaluating Cryo–TEM Reconstruction Accuracy of Self–Assembled Polymer Nanostructures

Cryogenic transmission electron microscopy (cryo–TEM) combined with single particle analysis (SPA) is an emerging imaging approach for soft materials. However, the accuracy of SPA–reconstructed nanostructures, particularly those formed by synthetic polymers, remains uncertain due to potential packing heterogeneity of the nanostructures. In this study, the combination of molecular dynamics (MD) simulations and image simulations is utilized to validate the accuracy of cryo–TEM 3D reconstructions of self–assembled polypeptoid fibril nanostructures. Using CryoSPARC software, image simulations, 2D classifications, ab initio reconstructions, and homogenous refinements are performed. By comparing the results with atomic models, the recovery of molecular details is assessed, heterogeneous structures are identified, and the influence of extraction location on the reconstructions is evaluated. In conclusion, these findings confirm the fidelity of single particle analysis in accurately resolving complex structural characteristics and heterogeneous structures, exhibiting its potential as a valuable tool for detailed structural analysis of synthetic polymers and soft materials.

36 MATERIALS SCIENCE↗
DOE OSTI2020

An in situ shearing x-ray measurement system for exploring structures and dynamics at the solid–liquid interface

Revealing interfacial structure and dynamics has been one of the essential thematic topics in material science and condensed matter physics. Synchrotron-based x-ray scattering techniques can deliver unique and insightful probing of interfacial structures and dynamics, in particular, in reflection geometries with higher surface and interfacial sensitivity than transmission geometries. We demonstrate the design and implementation of an in situ shearing x-ray measurement system, equipped with both inline parallel-plate and cone-and-plate shearing setups and operated at the advanced photon source at Argonne National Laboratory, to investigate the structures and dynamics of end-tethered polymers at the solid-liquid interface. With a precise lifting motor, a micrometer-scale gap can be produced by aligning two surfaces of a rotating upper shaft and a lower sample substrate. A torsional shear flow forms in the gap and applies tangential shear forces on the sample surface. The technical combination with nanoscale rheology and the utilization of in situ x-ray scattering allow us to gain fundamental insights into the complex dynamics in soft interfaces under shearing. In this work, we demonstrate the technical scope and experimental capability of the in situ shearing x-ray system through the measurements of charged polymers at both flat and curved interfaces upon shearing. Furthermore, through the in situ shearing x-ray scattering experiments integrated with theoretical simulations, we aim to develop a detailed understanding of the short-range molecular structure and mesoscale ionic aggregate morphology, as well as ion transport and dynamics in soft interfaces, thereby providing fundamental insight into a long-standing challenge in ionic polymer brushes with a significant technological impact.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Structure–Thermodynamic Relationship of a Polysaccharide Gel (Alginate) as a Function of Water Content and Counterion Type (Na vs Ca)

Biofilms are the predominant mode of microbial life on Earth, and so a deep understanding of microbial communities-and their impacts on environmental processes-requires a firm understanding of biofilm properties. Because of the importance of biofilms to their microbial inhabitants, microbes have evolved different ways of engineering and reconfiguring the matrix of extracellular polymeric substances (EPS) that constitute the main non-living component of biofilms. This ability makes it difficult to distinguish between the biotic and abiotic origins of biofilm properties. An important route toward establishing this distinction has been the study of simplified models of the EPS matrix. This study builds on such efforts by using atomistic simulations to predict the nanoscale (≤10 nm scale) structure of a model EPS matrix and the sensitivity of this structure to interpolymer interactions and water content. To accomplish this, we use replica exchange molecular dynamics (REMD) simulations to generate all-atom configurations of ten 3.4 kDa alginate polymers at a range of water contents and Ca–Na ratios. Simulated systems are solvated with explicitly modeled water molecules, which allows us to capture the discrete structure of the hydrating water and to examine the thermodynamic stability of water in the gels as they are progressively dehydrated. Our primary findings are that (i) the structure of the hydrogels is highly sensitive to the identity of the charge-compensating cations, (ii) the thermodynamics of water within the gels (specific enthalpy and free energy) are, surprisingly, only weakly sensitive to cation identity, and (iii) predictions of the differential enthalpy and free energy of hydration include a short-ranged enthalpic term that promotes hydration and a longer-ranged (presumably entropic) term that promotes dehydration, where short and long ranges refer to distances shorter or longer than ~0.6 nm between alginate strands.

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

Polyelectrolyte Complex Coacervation across a Broad Range of Charge Densities

Polyelectrolyte complex coacervates of homologous (co)polyelectrolytes with a near-ideally random distribution of a charged and neutral ethylene oxide comonomer were synthesized. The unique platform provided by these building blocks enabled an investigation of the phase behavior across charge fractions 0.10 ≤ f ≤ 1.0. Experimental phase diagrams for f = 0.30–1.0 were obtained from thermogravimetric analysis of complex and supernatant phases and contrasted with molecular dynamics simulations and theoretical scaling laws. At intermediate to high f, a dependence of polymer weight fraction in the salt-free coacervate phase (w P,c ) of w P,c ~ f 0.37±0.01 was extracted; this trend was in good agreement with accompanying simulation predictions. Below f = 0.50, w P,c was found to decrease more dramatically, qualitatively in line with theory and simulations predicting an exponent of 2/3 at f ≤ 0.25. Preferential salt partitioning to either coacervate or supernatant was found to be dictated by the chemistry of the constituent (co)polyelectrolytes.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Exponential Water Uptake in Ionomer Membranes Results from Polymer Plasticization

Water content is the most influential factor in the performance of ion-exchange membranes (IEMs), controlling their mechanical rigidity, ionic conductivity, and degradation rates. Membranes absorb water exponentially at high water activity (a w ), making that region of the sorption isotherm the most influential on membrane properties. Plasticization of the polymer by water has been proposed to cause this exponential uptake at high a w . However, an integrated microscopic picture of the structure, thermodynamic, and mechanical properties of IEMs as a function of a w has remained elusive. In this work, we use large-scale molecular simulations validated with experimental measurements to compute the sorption isotherms, Young’s modulus, polymer dynamics, and structure of IEMs. The simulations unveil that the exponential increase in water uptake (WU) coincides in all cases with the glass to rubber transition of the membrane, as measured through its Young’s modulus and segmental polymer dynamics. Functionalization of the polymer with alkyl groups further contributes to its plasticization, increasing the WU at a given a w and ion-exchange capacity (IEC). The alkyl chains act synergistically with water to plasticize the polymer matrix and allow water penetration in the membrane. The simulations reveal that the width of the water channels depends on the ratio λ of water to ions in the polymer but is independent of its IEC. We conclude that differences in the polymer matrix–not the water channels–are responsible for the distinct uptake response of IEMs to the thermodynamic driving force of water activity.

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

Role of Intradomain Heterogeneity on Ion and Polymer Dynamics in Block Polymer Electrolytes: Investigating Interfacial Mobility and Ion-Specific Dynamics and Transport

Block polymers show promise as solid-state battery electrolytes due to the optimization of conductive and mechanical properties enabled via tuning of block chemistry and length. We investigate a polystyrene-block-poly(oligo-oxyethylene methacrylate) (PS-b-POEM) electrolyte doped with various lithium salts to investigate the role of molecular structure on ion transport properties and on local ion dynamics and associations. Anion charge becomes more delocalized with increasing size, reducing the coupling between salt ions while increasing coupling between ion and polymer chain motions and creating a more mobile overall environment. We observe support for this ion-polymer coupling via 1 H, 7 Li and 19 F NMR spectroscopy, from which we obtain ion-specific mobility transition temperatures that differ from the polymer glass transition temperature. We also note faster transport and weaker local energetic interactions with anion size using temperature-dependent NMR diffusometry. 1 H NMR spectroscopy further elucidates polymer chain dynamics and enables quantification of the temperature-dependent fraction of the conducting block that is immobile near the PS-POEM domain interface. As a result, NMR thus represents a species-specific and timescale-specific platform to quantify phase and interface behavior, and to correlate ion-specific transport with polymer chain dynamics.

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

Learning Latent Representations to Bridge Coarse-Grained and Atomistic Resolutions in Polymer Simulations

We present a machine-learning-based framework for learning reduced-order representations of polymer chain conformations across coarse-grained (CG) and united-atom (UA) fidelities. By employing linear singular value decomposition and nonlinear autoencoders, we compress high-dimensional polymer configurations into latent spaces with minimal loss of structural accuracy. Crucially, we demonstrate a near-perfect linear mapping between CG and UA latent spaces, enabling an efficient super-resolution back-mapping procedure that reconstructs high-fidelity UA configurations from CG simulations. While minor structural inaccuracies occur, they are effectively corrected through a brief molecular dynamics relaxation, forming a practical hybrid machine learning−physics scheme. This approach establishes the key structural prerequisites for accelerated polymer dynamics simulations: a compact and accurate latent encoding of polymer chain conformations and a validated multi-fidelity mapping that permits reconstruction of UA structures from CG configurations. The extension of this framework to explicit time evolution within the latent space, enabling dynamics to be propagated at CG fidelity and decoded to UA resolution only when required, represents a natural and well-motivated direction for future work.

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

How Topological Polymer Loops on the Nanoparticle Surface Control the Mechanical Properties of Nanocomposites

Carbon black (CB) and silica (SiO 2 ) filled elastomers are known to be the most successful polymer nanocomposites (PNCs) in industry, where “bound rubber (BR)” (i.e., polymer chains that are physically or chemically adsorbed on the nanofiller surface) plays a critical role in their reinforcement. Here, we report a molecular-scale mechanism underlying the “BR-induced reinforcement” by integrating neutron scattering experiments and molecular dynamics simulations. Simplified non-cross-linked SiO 2 -filled polybutadiene (PB) and CB-filled PB reveal the critical role of topological polymer loops in the BR for the enhanced mechanical performance. The average loop size on the SiO 2 surface modified with a silane coupling agent is much smaller than that on the CB surface and the loops on the SiO 2 surface are densely formed, preventing interdigitation with the matrix chains. On the other hand, the larger, uncrowded loops formed on the CB surface facilitate the interdigitation with the matrix polymer chains even near the filler surface. In this way, a strong connectivity is established between a matrix and a nanofiller, resulting in an adhesive filler–polymer interface. Furthermore, our findings shed light on rich and complex physics and materials design problems in PNCs, where the topological polymer structure on the nanofiller surface directly controls the macroscopic mechanical properties.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Hydrogen Bonding Inside Anionic Polymeric Brush Layer: Machine Learning-Driven Exploration of the Relative Roles of the Polymer Steric Effect, Charging, and Type of Screening Counterions

This paper employs a combination of all-atom molecular dynamics (MD) simulations and unsupervised machine learning (ML) for studying the water-water hydrogen bonds (HBs) inside the anionic poly-acrylic acid (PAA) brushes modeled using all-atom MD simulations. PAA brush layer with different charge fraction (f), namely f=0, f=0.25, and f=1, is considered. Water-water interactions, both inside and outside the brush layer, are represented through distinct clusters of tupules of variables representing distances associated with the interacting water molecules. While clusters representing the HBs are present for water inside and outside the brushes, several clusters representing the long-range water-water interactions are missing for the water molecules inside the highly charged (f=1) PAA brushes. More importantly, inside highly charged brushes, the edge of the clusters representing the water-water HBs is progressively shortened, as compared to that in the bulk. Both these results stem from the presence of the PAA brushes imparting the steric effect and the charge effect, or the effect associated with enhanced interactions of water molecules with PE charges and counterions, thereby disrupting the water connectivity. This water-charged-species interaction also increases the water-water HB angle, i.e., makes the water-water HBs less stable inside the highly charged PAA brush layer. The narrowing of the clusters representing the HBs and the alteration of the angle characterizing the HBs confirm that the conditions defining the water-water HBs change inside the PAA brush layer as a function of the charges on the PAA brush layer. Furthermore, we show that the use of the generic definition of HBs, as compared to using our simulation-motivated modified definition of water-water HBs, overpredict the number of water-water HBs inside the PAA brush layer. Finally, we employ this all-atom-MD-ML framework to quantify the effect of other types of screening counterions (Li + , Ca 2+ , and Y 3+ ions) in determining the water-water interactions and water-water HB properties inside the PAA brush layer. Furthermore, the findings of the present study, confirming the weakening of water-water HBs inside the PAA brush layer, points to the possibility that the water molecules will be more available for hydrating the brush layer and counterions, thereby leading to a more pronounced wetting of the PAA brush layer.

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

Percolated Ionic Aggregate Morphologies and Decoupled Ion Transport in Precise Sulfonated Polymers Synthesized by Ring-Opening Metathesis Polymerization

Here, we describe a set of precise single-ion conducting polymers that form self-assembled percolated ionic aggregates in glassy polymer matrices and have decoupled transport of metal cations. These precise single-ion conductors (SICs), synthesized by a scalable ring-opening metathesis polymerization, consist of a polyethylene backbone with a sulfonated phenyl group pendant on every fifth carbon and are fully neutralized by a counterion X + (Li + , Na + , or Cs + ). Experimental X-ray scattering measurements and fully atomistic molecular dynamics (MD) simulations are in good agreement. The MD simulations show that the ionic groups nanophase separate from the polymer backbone to form percolating ionic aggregates. Using graph theory, we find that within the Li + - and Na + -neutralized polymers the percolated aggregates exhibit planar and ribbon-like configurations at intermediate length scales, while the percolated aggregates within the Cs + -neutralized polymers are more isotropic. Electrical impedance spectroscopy measurements show that the ionic conductivities exhibit Arrhenius behavior, with conductivities of 10 –7 to 10 –6 S/cm at 180 °C. In the MD simulations, the cations move between sulfonate groups in the percolated aggregates, larger ions travel further, and overall cations travel further than the polymer backbones, indicating a decoupled ion-transport mechanism. Thus, the percolated ionic aggregates in these polymers can serve as pathways to facilitate decoupled ion motion through a glassy polymer matrix.

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

Cation–polymer interactions and local heterogeneity determine the relative order of alkali cation diffusion coefficients in PEGDA hydrogels

Current research efforts are focused on endowing polymer membranes with ion–ion selectivity by incorporating ion–polymer interactions into materials to bias the selective partitioning and or diffusivity of one species over another. However, little is known about the impact of such interactions on the mechanisms of ion transport. In this study, we probe the influence of cation–polymer interactions on cation, anion, and salt diffusivity in a model membrane material, poly(ethylene glycol) diacrylate (PEGDA) by modeling concentrated polyethylene oxide solutions via molecular dynamics simulations. These results are compared to published experimental data for LiCl, NaCl, and KCl diffusion in PEGDA. Experimentally, the order of salt and cation diffusion coefficients for LiCl, NaCl, and KCl deviate from the order in aqueous solutions. Here, simulations identify these deviations to arise from cation–polymer coordination in the membrane. Both the fraction of bound cations and the average binding lifetime increases with decreasing cation hydration free energy (moving down the alkali series), leading to different diffusivity trends in the membrane compared to solution. However, to recover the experimentally observed order of diffusivities cations and salt in our simulations, we needed to incorporate membrane heterogeneity explicitly via a polymer charge scaling procedure. Together, our results indicate that cation–polymer interactions, as well as spatial heterogeneity within the membrane, play a critical role in dictating the observed order of alkali cation and salt diffusion coefficients in membranes.

36 MATERIALS SCIENCE↗