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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Robust chain aggregation of low-entropy rigid ladder polymers in solution

Conjugated polymers have been widely investigated where ladder-type conjugated polymers receive more attention due to their rigid backbones and extraordinary properties. However, the understanding of how the rigid conformation of ladder polymers translates to material properties is still limited. Here, we systematically investigated the solution aggregation properties of a carbazole-derived conjugated ladder polymer (LP) and its analogous non-ladder control polymer (CP) via light scattering, neutron scattering, and UV-vis absorption spectroscopy characterization techniques, revealing a highly robust, temperature-insensitive aggregation behavior of the LP. The experimental findings were further validated by computational molecular dynamics simulations. Here we found that the peak positions and intensities of the UV spectra of the LP remained constant between 20 °C and 120 °C in chlorobenzene solution. The polymer also showed a stable hydrodynamic radius measured by dynamic light scattering from 20 °C to 70 °C in the chlorobenzene solution. Using small-angle neutron scattering, no Guinier region was reached in the measured q range down to 0.008 Å –1 , even at elevated temperature. In contrast, the non-ladder control polymer CP was fully soluble in the chlorobenzene solvent without the observation of any notable aggregates. The Brownian dynamics simulation showed that during polymer aggregation, the entropy change of the LP was significantly less negative than that of the non-ladder control polymer. These findings revealed the low entropy nature of rigid conjugated ladder polymers and the low entropy penalty for their aggregation, which is promising for highly robust intermolecular interactions at high temperatures. Such a unique thermodynamic feature of rigid ladder polymers can be leveraged in the design and application of next-generation electronic and optoelectronic devices that function under unconventional high temperature conditions.

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

Building and Breaking Carbon Composites with REACTER

Carbon-based composites have become indispensable materials in aerospace and other high-performance applications, yet obtaining a detailed, nanoscale understanding of their morphology and failure mechanisms using only experimental methods remains a difficult challenge. REACTER is a versatile computational modeling tool for atomistic molecular dynamics simulations designed to model chemical reactions at the speed and length scales of classical force fields. In this work, several recent features of REACTER were applied to the creation and subsequent mechanical testing of two classes of carbon composites: carbon nanotube (CNT) composites and carbon fiber reinforced polymers (CFRP). A network of CNTs was grown dynamically using the new ‘create atoms’ feature of REACTER. The CNT filler was embedded into a polyarylacetylene (PAA) matrix by simulated in situ polymerization to obtain the final composite model. To generate the second class of carbon composite, fully carbonized (graphitic) carbon fiber morphologies were created by the method of Desai et al. [1], but using the advanced reaction constraints framework of REACTER. Two fiber models were created, representing a circular carbon fiber core and a flat surface, and similarly infiltrated with resin to obtain the final CFRP structure. Failure mechanisms were elucidated by simulating mechanically induced bond breaking, as characterized by third order DFT-based tight-binding simulations, via a reaction constraint on the total potential energy of the involved atoms.

Polymer↗
DOE OSTI2024

Uncovering Backbone Conformation for Rigid DPP-Based Donor–Acceptor Conjugated Polymer Using Deuterium Labeling and Neutron Scattering

The conjugated polymer’s backbone conformation dictates the delocalization of electrons, ultimately affecting its optoelectronic properties. Most conjugated polymers can be viewed as semirigid rods with their backbone embedded among long alkyl side chains. Thus, it is challenging to experimentally quantify the conformation of a conjugated backbone. Here, we performed contrast variation neutron scattering on rigid conjugated donor–acceptor (D–A) diketopyrrolopyrrole (DPP) polymers with selectively deuterated side chains to measure the conjugated backbone conformation. We first synthesized DPPbased polymers with deuterated side chains, confirmed by NMR and FTIR. Using contrast variation neutron scattering, we found that the DPP-based conjugated polymers are much more rigid than poly(3-alkylthiophenes), with persistence length (L p ) at 16–18 nm versus 2–3 nm. More importantly, in contrast to the relatively flexible poly(3-alkylthiophenes) whose backbone is more flexible than the whole polymer, we found that the backbone of DPP-based polymers has the same L p value compared to the whole polymer chain. This indicates that side chain interference on backbone conformation is not present for the semirigid polymer, which is further confirmed by coarse-grained molecular dynamics (CG-MD) simulations. Our work provides a novel protocol to probe polymer’s backbone conformation and paradigm-shifting understanding of the backbone conformation of semirigid conjugated polymers.

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

Large-scale computational polymer solubility predictions and applications to dissolution-based plastic recycling

Dissolution-based plastic recycling is a promising approach to separate and recover high quality pure polymer resins from multicomponent plastic waste by exploiting differences in polymer solubility. The design of a dissolution-based polymer recycling process requires the selection of appropriate solvent systems and operating temperatures to dissolve only target polymers. Determining these parameters experimentally is challenging due to the wide range of solvents and temperatures possible for a given set of target polymers. In this work, we report a computational scheme that employs molecular dynamics simulations and the Conductor-like Screening Model for Realistic Solvents to predict polymer solubilities. Using this scheme, we established a computational solubility database for 8 common polymers and 1007 solvents at multiple temperatures and measured selected solubilities experimentally to validate computational predictions. Analysis of functional groups within this large database then provides chemical heuristics relating the molecular structures of good and non-solvents for selected polymers. We further developed a tool that automates the selection of solvents for all possible sequences in which target polymers can be selectively dissolved to guide the design of dissolution-based plastic recycling processes. Here, we demonstrate the application of these methods via multiple experimental case studies of representative dissolution-based polymer recycling processes in which pure polymer resins were successfully recovered from physical mixtures of polymers.

42 ENGINEERING↗
DOE OSTI2021

Spatiotemporal mapping of mesoscopic liquid dynamics

The study of liquid dynamics at mesoscopic scales is still strewn with difficulty due to limitations in theory and experiment. Historically, significant attention has been given to the analysis of space-time correlation functions and their frequency-Fourier transforms at a few discrete wave numbers. The massive computing power afforded by modern high performance computing clusters and the advent of a wide-angle neutron spin-echo spectrometer, however, have unlocked a more intuitive and fruitful approach to this problem. Using molecular dynamics simulations, here we demonstrate the benefits of spatiotemporally mapping intermediate scattering functions on a dense grid of correlation times and wave numbers. Four model systems are investigated: a Lennard-Jones liquid, a coarse-grained bead-spring polymer, a molten sodium chloride, and a poly(ethylene oxide) melt. We show that the spatiotemporal mapping approach is particularly useful for elucidating the mesoscopic dynamics in these liquids, where several underlying mechanisms, such as molecular relaxations, hydrodynamic modes, and nonhydrodynamic excitations, are potentially at play. Therefore, compared to the traditional method, direct visualization of density space-time correlation functions on two-dimensional color maps permits appraisals of complicated dynamical behavior at mesoscales in a global manner. For example, the scaling relations between space and time for different types of molecular motions can be straightforwardly identified on these plots, without any model-dependent analysis. Additionally, we show how theoretical ideas regarding collective mesoscopic dynamics, such as the classical hydrodynamic theory, the convolution approximation, and a recently proposed phenomenological model, can be discussed in terms of the global features of spatiotemporal maps of intermediate scattering functions. The new perspective offered by the spatiotemporal mapping method should prove useful for the study of liquid dynamics in general.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Local conformations and heterogeneities in structures and dynamics of isotactic polypropylene adsorbed onto carbon fiber

Carbon fiber (CF) reinforced polymers (CFRPs) have experienced widespread use in various industries. One of the most important parameters that controls the macroscopic property of CFRPs is the interface between a polymer matrix and CF. There is growing evidence to suggest the formation of a bound polymer layer (BPL), i.e., polymer chains that physically adsorb on a filler surface. However, this interface is always in contact with the thicker part of a polymer matrix, rendering its understanding a difficult task. Therein, we use CF-reinforced isotactic polypropylene (iPP) as a rational CFRP. To characterize the BPL on the CF surface, we extracted it from the CFRP using solvent-rinsing with p-xylene. The physical and thermal properties of the BPL were characterized by differential scanning calorimetry and thermogravimetric analysis, while its microscopic structures and dynamics were probed by small-angle neutron scattering and quasi-elastic neutron scattering (QENS) techniques. Subsequently, we employed atomistic molecular dynamics (MD) simulations to complement the QENS results above the bulk melting temperature and reveal details that were experimentally inaccessible. In this work, we observed that the degree of crystallinity of the BPL was quite lower than the bulk, while the melting temperature of the BPL remained the same as the bulks. Within the given length and time scales probed by QENS and MD, we also observed that most of the bound chains were mobile, with the formation of a high-density region (less than 1 nm in thickness) near the CF surface. The segmental dynamics of the bound chains probed by both QENS and MD were also much faster than those of the free chains, possibly due to the presence of a free surface region at the topmost surface of the BPL. Furthermore, the MD results demonstrated that the backbone chains and side groups lie nearly flat on the CF surface, which is the driving force for the flattening process of the iPP BPL to overcome the conformational entropy loss in the total free energy.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Effect of Solvents on Lignin–Surface Interactions via Molecular Dynamics Simulations

Lignin, an essential building block of lignocellulosic biomass, is a potential abundant source of aromatic monomers for the polymer and chemical industry. Reductive catalytic fractionation (RCF) is one promising process that can produce high yields of phenolic monomers and oligomers from lignin under different catalytic conditions. An important choice in optimizing RCF is the selection of solvent; however, detailed insights into solvent effects on lignin behaviors and interactions remain limited. Here, in this work, we perform all-atom molecular dynamics simulations to study the solvation of lignin, solvent-mediated conformational changes, and the interaction of solvated lignin oligomers with model surfaces. We focus on the behavior of an oligomeric lignin model compound in methanol, ethanol, a binary mixture of ethanol and water, and water at both the RCF reaction temperature (473 K) and room temperature. Analysis of structural features of lignin suggests that these three organic solvent systems favorably solvate lignin, resulting in a more extended conformation suitable for catalytic conversion to valuable chemicals. We further introduce model palladium (Pd) and carbon (C) surfaces to understand how solvent choice impacts adsorption onto a representative catalytic surface and support, and to quantify the competition among the reactant and solvent molecules for the surface. Unbiased simulations suggest strong adsorption of lignin on both Pd and C surfaces at 473 K, with notable solvent-mediated differences in adsorption energies. Additionally, our findings indicate that lignin adsorption is promoted by the entropy change resulting from the displacement of solvent molecules from the surface. This study provides a molecular perspective of adsorption of lignin onto varying surfaces, which is a step towards understanding and optimizing the catalytic conversion of lignin into valuable chemicals.

adsorption↗
DOE OSTI2021

Composite entanglement topology and extensional rheology of symmetric ring-linear polymer blends

Extensive molecular simulations are applied to characterize the equilibrium dynamics, entanglement topology, and nonlinear extensional rheology of symmetric ring-linear polymer blends with systematically varied ring fraction ϕR. Chains with degree of entanglement Z≈14 are mixed to produce 10 well-entangled systems with ϕR varying from neat linear to neat ring melts. Z is large enough that except for very large ϕR, the rings are threaded by multiple linear chains in equilibrium. Primitive path analysis is used to visualize and quantify the structure of the composite ring-linear entanglement network. We measure the quantity of ring-linear threading and linear-linear entanglement as a function of ϕR and identify with simple arguments a ring fraction ϕR≈0.4 where the topological constraints of the composite entanglement network are maximized. These topological analyses are used to rationalize the ϕR-dependence of ring and linear chain dynamics, conformations, and rheology. Simulations of startup uniaxial elongation flows demonstrate the extensional stress overshoot observed in recent filament stretching experiments and characterize how it depends on the blend composition and entanglement topology. The overshoot is driven by an overstretching and recoil of ring polymers due to the convective unthreading of rings from linear chains.

O’Connor, Thomas C.↗
DOE OSTI2025

Shock-induced chemistry and high strain-rate viscoelastic behavior of a phenolic polymer

We use impact experiments and a finite element model (up to 1.2 GPa), and molecular dynamics simulations (up to 60 GPa), to examine the behavior of a phenolic polymer under shock compression, spanning both nonreactive and reactive regimes. In the nonreactive regime, relaxation following compression at strain rates of ∼105 s−1 can be explained by viscoelasticity observed at ordinary laboratory rates (≲1 s−1) by accounting for the temperature dependence of the phenolic β-transition. Reasonable agreement is found between the measured shock Hugoniot up to 1.2 GPa and molecular dynamics simulation for cross-linked structures of comparable density. We also observed a first-order mechanical transition near 0.36 GPa shock stress and estimated a spall strength of 0.102 GPa and Hugoniot elastic limit of 1–2 GPa. The shock stress is found to vary up to 24% among phenolics made with different resin and/or cure processes. Finally, molecular dynamics simulations are used to identify a reactive regime at shock pressures ≳20 GPa that is characterized by chemically driven, rate-dependent relaxation processes, including dehydrogenation and dehydration reactions that promote the formation of a dense, highly cross-linked carbonaceous solid and the release of light volatiles.

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

Role of solvation site segmental dynamics on ion transport in ethylene-oxide based side-chain polymer electrolytes

Ion conducting capability is often imparted to polymeric materials through short polyether side-chains, and yet the impact of this graft polymer architecture on ion solvation and conduction has not been fully explored. Here, we use a combination of impedance spectroscopy, vibrational spectroscopy, and atomistic molecular dynamics (MD) to compare the conductivity, ionic interactions, and polymer dynamics in a series of graft polyether electrolytes. We find that in poly[(oligo ethylene oxide)methyl ether methacrylate] (POEM), a widely used graft polymer electrolyte, the ionic conductivity drops more than an order of magnitude as the side-chain length is decreased from nine ethylene oxide (EO) units to three. This difference in conductivity is unexplained by differences in the calorimetric glass transition temperature (T g ), which varies only slightly with side-chain length. Furthermore, through vibrational spectroscopy and MD simulations we demonstrate that both linear and graft polyethers solvate Li⁺ ions effectively and dissociate them from large counterions, irrespective of side-chain length. Li⁺ ions do, however, show preferential solvation by EO units far from the methacrylate backbone. Similarly, EO units far from the backbone show enhanced segmental dynamics, while those near the immobile methacrylate group move substantially more slowly, as quantified by bond vector autocorrelation relaxation times. This heterogeneity in both ion solvation and local segmental relaxation explains variation in ion conductivity where material-averaged properties such as T g and number of free ions fail to do so. Importantly, the ionic conductivity is dictated primarily by the segmental mobility of the EO units which form effective solvation sites, rather than system-wide dynamics.

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

Predicting the Mechanical Response of Polyhydroxyalkanoate Biopolymers Using Molecular Dynamics Simulations

Polyhydroxyalkanoates (PHAs) have emerged as a promising class of biosynthesizable, biocompatible, and biodegradable polymers to replace petroleum-based plastics for addressing the global plastic pollution problem. Although PHAs offer a wide range of chemical diversity, the structure–property relationships in this class of polymers remain poorly established. In particular, the available experimental data on the mechanical properties is scarce. In this contribution, we have used molecular dynamics simulations employing a recently developed forcefield to predict chemical trends in mechanical properties of PHAs. Specifically, we make predictions for Young’s modulus, and yield stress for a wide range of PHAs that exhibit varying lengths of backbone and side chains as well as different side chain functional groups. Deformation simulations were performed at six different strain rates and six different temperatures to elucidate their influence on the mechanical properties. Our results indicate that Young’s modulus and yield stress decrease systematically with increase in the number of carbon atoms in the side chain as well as in the polymer backbone. In addition, we find that the mechanical properties were strongly correlated with the chemical nature of the functional group. The functional groups that enhance the interchain interactions lead to an enhancement in both the Young’s modulus and yield stress. Finally, we applied the developed methodology to study composition-dependence of the mechanical properties for a selected set of binary and ternary copolymers. Overall, our work not only provides insights into rational design rules for tailoring mechanical properties in PHAs, but also opens up avenues for future high throughput atomistic simulation studies geared towards identifying functional PHA polymer candidates for targeted applications.

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

Molecular Alignment of a Meta-Aramid on Carbon Nanotubes by In Situ Interfacial Polymerization

Molecularly organized nanocomposites of polymers and carbon nanotubes (CNTs) have great promise as high-performance materials; in particular, conformal deposition of polymers can control interfacial properties for mechanical load transfer, electrical or thermal transport, or electro/chemical transduction. However, controllability of polymer-CNT interaction remains a challenge with common processing methods that combine CNTs and polymers in melt or in solution, often leading to nonuniform polymer distribution and CNT aggregation. In this work, we demonstrate CNTs within net-shape sheets can be controllably coated with a conformal coating of meta-aramid by simultaneous capillary infiltration and interfacial polymerization. We determine that π-interaction between the polymer and CNTs results in chain alignment parallel to the CNT outer wall. Subsequent nucleation and growth of the precipitated aramid forms a smooth continuous layered sheath around the CNTs. These findings motivate future investigation of mechanical properties of the resulting composites, and adaptation of the in situ polymerization method to other substrates.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
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 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 OSTI2022

Understanding curing dynamics of arylacetylene and phthalonitrile thermoset blends

Tetrakis(phenylethynyl)benzene (TPEB), a high char yield arylacetylene resin, was blended with a resorcinol-based PEEK™-like oligomeric phthalonitrile resin (Res). Here, we observed the influence of each polymer precursor on the structural changes and material dynamics of different mixtures of these resins during thermal curing to produce thermosets with high char retention. We provide insight to help understand fundamental processes that control and tune this complex thermally driven curing mechanism, which, to date, has not been well understood. We show that an equivalent blend of Res and TPEB mitigates the exothermic curing process and exhibits a low viscosity, which may facilitate processing of these thermosetting blends for composite manufacturing. To obtain a more comprehensive overview of the different dynamics of these polymers during curing, we correlated thermodynamic changes with rheology, materials characterization, and molecular dynamics derived from neutron scattering. Our efforts shed more light on the thermally driven chemical changes, physical transformation, and gelation dynamics that drive or accompany softening and curing of these two distinct crosslinking resin groups.

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

Critical Role of the Interfacial Layer in Associating Polymers with Microphase Separation

Polymers with dynamic (transient) bonds, often called associating polymers, have been attracting significant attention in recent years due to their unique viscoelastic properties, self-healing ability, and recyclability. Nevertheless, understanding the mechanisms and the factors controlling their macroscopic properties remains limited due to the higher complexity introduced by the dynamic bonds. Here in this study, small-angle X-ray scattering (SAXS), broadband dielectric spectroscopy (BDS), and rheology were applied to unravel the structure and dynamics of telechelic associating polymers with different molecular weights. SAXS measurements revealed phase separation of the functional end groups with an average cluster size of ~2–3 nm and the distance between clusters controlled by the chain length. Borrowing the interfacial layer model analysis of BDS data from the polymer nanocomposite field, we demonstrated the presence of an interfacial polymer layer with a thickness of ~0.7–0.9 nm surrounding these clusters. Rheological measurements showed quantitatively that the presence of the interfacial layer significantly alters the viscoelastic behavior of these materials, indicating the crucial role of the interfacial layer in defining the macroscopic mechanical properties of the studied telechelic materials. The presented results emphasize that phase separation of the functional groups in associating polymers leads to very significant changes of the viscoelastic properties, opening a promising avenue in the design of novel functional materials.

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

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↗