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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 397 records · Page 22

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Computing inelastic neutron scattering spectra from molecular dynamics trajectories

Abstract Inelastic neutron scattering (INS) provides a weighted density of phonon modes. Currently, INS spectra can only be interpreted for perfectly crystalline materials because of high computational cost for electronic simulations. INS has the potential to provide detailed morphological information if sufficiently large volumes and appropriate structural variety are simulated. Here, we propose a method that allows direct comparison between INS data with molecular dynamics simulations, a simulation method that is frequently used to simulate semicrystalline/amorphous materials. We illustrate the technique by analyzing spectra of a well-studied conjugated polymer, poly(3-hexylthiophene-2,5-diyl) (P3HT) and conclude that our technique provides improved volume and structural variety, but that the classical force field requires improvement before the morphology can be accurately interpreted.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Relationship between microporous structure and light gas transport through glassy polymeric membranes revealed by molecular simulations

Microporous glassy polymers are attractive materials for gas separation membranes, due to their high permeability and tailorable selectivity, provided physical aging can be delayed. The archetypal microporous glassy polymer, PTMSP, can be blended with a hyper-crosslinked isatin–triptycene porous polymer network (PPN) to delay physical aging. However, while PPN is effective at reducing physical aging, it also affects the permeability of light gases through PTMSP. Molecular dynamics simulations were used here to shed fundamental light on the mechanisms responsible for these effects. Atomistic models are developed that satisfactorily reproduce experimental observations such as matrix density and cavity size distributions for neat PTMSP as well as for PTMSP–PPN blends. Analysis of the simulation results suggests that physical aging is delayed because the PPN inclusions slow down PTMSP relaxation while reducing the connectivity between free volume pockets. To understand how PPN inclusions affect light gas permeability, the atomistic models developed are used to probe CO 2 and CH 4 diffusion and sorption. These simulations are conducted for PTMSP matrices exhibiting varying density, towards reproducing experimental permeability data. Interrogating the simulation trajectories, it is discovered that while CH 4 travels preferentially through the free volume cavities, CO 2 preferentially interacts with the available surfaces, especially in the presence of PPN. These differences help interpret experimental observations. The transport of C 2 H 6 and H 2 S through PTMSP matrices was also investigated. These gases, and in particular H 2 S, were found to absorb within PTMSP, yielding very low diffusion coefficients. The ability to predict differences in diffusion pathways for various natural gas components unveils the possibility of engineering membranes to control permeability and selectivity towards large-scale applications.

Bao Le, Tran Thi [Univ. of Oklahoma, Norman, OK (U↗

DOE OSTI2020

In situ multimodal imaging for nanoscale visualization of tribofilm formation

A third of the energy from fuel combustion in passenger car gasoline engines is lost due to friction. Carefully designed engine lubricants can recover some of these losses by reducing friction and wear by forming a nanometer-scale chemico-physico tribofilm between surfaces. Accordingly, attention has focused on developing oil formulations that form low-friction tribofilms. However, analyses of resultant tribofilms are typically conducted after tribo-tests with conventional characterization tools and do not offer insights into tribofilm formation and evolution, precluding information critical to tuning tribofilm properties. In this work, we developed a unique multimodal methodology based on Atomic Force Microscopy (AFM) with local probe heating for in situ tribological studies that activates friction modifiers and simultaneously captures the evolution of friction and surface roughness, with nanometer resolution. As a platform to demonstrate the ability of this methodology to visualize dynamics of tribofilm formation in situ, we apply it to molybdenum-based friction modifiers to distinguish key factors in their functionality and correlate nanoscale AFM and Friction Force Microscopy data to bench tribo-tests used in the industry. To decode the formation mechanisms observed in situ and underlying chemistry of tribofilms, we performed ab initio Molecular Dynamics (AIMD) simulations at comparable conditions. AIMD simulations confirmed both nanoscale and bench tribo-test results and showed deviations in molecular organization in tribofilms that are formed from different molybdenum-based friction modifiers that corroborates with surface functionality. With this innovative methodology, we demonstrate proof-of-principle in situ formation of molybdenum-based tribofilms directly on steel surfaces that could be applied generally to studying tribofilm formation.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Molecular dynamics based study on the effects of cation size on the local structure and diffusion in polymerized ionic liquids

We have used coarse-grained molecular dynamics simulations to understand the effect of cation size on ion diffusion in polymerized ionic liquids at temperatures well above the glass transition temperature of the polymers. We investigated dependencies of the diffusion constant on the cation radius and static dielectric constant and interpreted these results in terms of underlying structural changes, and decoupling phenomenon. We have found non-monotonic effects of the cation radius on the diffusion constant, with a maximum at an intermediate radius resulting from two different size dependent effects. Changes in the radial distribution functions characterizing spatial distribution of cations with respect to anions result in a monotonic increase in the coordination number. This yields an increase in the effective interaction energy as a function of cation radius, which can explain the decrease in the diffusion constant. This increase is counteracted at small cation radii by a corresponding increase in decoupling of cation and anion dynamics. In addition, the diffusion constant is found to increase on increasing the static dielectric constant for all radii of the cations. These results highlight the significant role of decoupling in designing polymerized ionic liquids with an enhanced diffusion constant of cations.

Gillespie, Colin [ORNL] (ORCID:0000000320287080)↗

DOE OSTI2023

Rapid activation of non-oriented mechanophores via shock loading and spallation

Mechanophores (MPs), stimuli-responsive molecules that respond chromatically to mechanochemical reactions, are important for understanding the coupling between mechanics and chemistry as well as in engineering applications. However, the atomic-level understanding of their activation originates from gas phase studies or under simple linear elongation forces directly on molecules or polymer chains containing MPs. The effect of many-body distortions, pervasive in condensed-phase applications, is not understood. Therefore, we performed large-scale molecular dynamics (MD) simulations of a poly(methyl methacrylate)-spiropyran copolymer under dynamic mechanical loading and studied the activation of the MP under various conditions from dynamical compression to tension during unloading. Detailed analysis of the all-atom MD trajectories shows that the MP blocks experience significant many-body intramolecular distortion that can significantly decrease the activation barrier as compared with when deformation rates are slow relative to molecular relaxation time scales. Here we find that the reactivity of MPs under material compression states is governed by many-body effects of intramolecular torsions, whereas under tension, the reactions are governed by tensile stresses.

36 MATERIALS SCIENCE↗

NASA NTRS1972

Investigations to improve carbon dioxide control with amine and molecular sieve type sorbers

The optimization trends and operating parameters of an integral molecular sieve bed heat exchanger were investigated. The optimum combination of substrate and coating for the HS-B porous polymer was determined based on the CO2 dynamic capacity in the presence of water vapor. Full size HS-B canister performance was evaluated. An Amine CO2 Concentrator utilizing IR-45 sorber material and available Manned Orbiting Laboratory hardware was designed, fabricated and tested for use as an experiment in the NASA 90-day space simulator test of 1970. It supported four men in the simulator for 71 days out of the 90-day test duration.

Bertrand, J. F.↗

DOE OSTI2020

Polymer Chain Diffusion in All-Polymer Nanocomposites: Confinement vs Chain Acceleration

All-polymer nanocomposites, in which soft, polymer-based nanoparticles are dispersed in the polymer matrix, have received great interest lately due to their potential use in a range of applications, including drug delivery and self-healing materials. However, the impact of this new class of nanoparticles on the dynamics of a linear polymer matrix in an all-polymer nanocomposite is still largely unknown. In this work, we report that the addition of polystyrene soft nanoparticles accelerates the diffusion of high molecular weight linear PS chains over a range of nanoparticle loadings. Our results show that at nanoparticle loadings below 1%, the diffusion of the linear matrix increases, presumably via a constraint release mechanism. At loadings above 1%, the increase in diffusion is mitigated by confinement effects of the nanoparticles. Thus, the response of these all-polymer nanocomposites is dominated by the balance of entropic confinement of the chain, which slows diffusion and a constraint release mechanism that speeds up the diffusion. However, the diffusion of the linear chain in the all-polymer nanocomposite is faster than that of the same chain in the melt at all loadings, diverging from the behavior of most nanocomposites with hard, impenetrable nanoparticles. Thus, the mechanism that accelerates the chain diffusion dominates in these systems. This behavior is unusual and fundamentally different than what has been reported for nanocomposites with hard inorganic nanoparticles, indicating that new perspectives are needed for these materials where the control of loading can either accelerate or decelerate the dynamics of the matrix in a distinctive manner.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Exploring the Origins of Anti-Ambipolarity in BBL Polymer: Links to Redox Chemistry, Electronic Structure, and Structural Dynamics

We examine the intrinsic physical-chemical properties of the conjugated ladder-type polymer poly(benzimidazobenzophenanthroline) (BBL) in response to electron transfer. We aim at explaining the origin of the anti-ambipolar behavior behind the observed BBL nonlinear response associated with specific device architectures. To elucidate this point, we use theory and computation based on first principles, including density functional theory optimizations, ab initio molecular dynamics, time-dependent DFT, and Marcus-theory analysis. Our results reveal that this redox response is not simply monotonic but follows an alternating odd/even pattern in which gap narrowing and reopening occur sequentially before near-gapless behavior emerges at high charging. Converging theoretical evidence in this work demonstrates that bell shaped conductivity in BBL originates in its fundamental electronic structure and supramolecular organization.

FOS: Physical sciences↗

DOE OSTI2024

Source of Processable Vitrimer Viscosities: Swap Frequencies and Steric Factors

Vitrimers exhibit high, processable viscosities, where other polymers do not, and are among the most promising polymers for closed-loop material circularity. We sought to investigate the underlying chemical kinetic factors that result in high viscosities for vitrimers, which are crucial to designing vitrimers with tunable viscosity. To interrogate these factors, we achieved the first simulated predictions of real vitrimer viscosities, using a novel kinetic Monte Carlo molecular dynamics method, overcoming the time and length scale gaps to predict experimental bulk viscosities. The vitrimer architecture investigated is based on poly(dimethylsiloxane) chains and vinylogous urethane bond swaps. We probed the effects of the extent of free swapping groups, %F, the activation energy, E A , and the steric factor, ρ. The steric factor is related to the intrinsic reaction probability for molecules with sufficient energy. All three factors were found to be significant, but the role of ρ was found to be the biggest and also the most underappreciated. The results show that the inclusion of accurate ρ is of critical importance for viscosity predictions, with the evidence suggesting that the typical assumption of ρ = 1 is not valid for vitrimers and that, indeed, very low steric factors are present in bond-swap vitrimers such that values of ρ < 10 –10 may be typical. This greatly influences the bond exchange rates and, ultimately, the viscosities. Recognition of this result is necessary for the prediction of vitrimer viscosities from molecular simulations and to make vitrimers by design from molecular dynamics. We also investigated the effects that E A , ρ, and the number of free swapping groups have upon vitreous range temperatures, TV, with respect to achieving a specific viscosity (η V = 1 × 10 8 Pa·s), as well as for a commonly reported higher viscosity extrapolation (η V = 1 × 10 12 Pa·s). The evidence suggests that vitrimers may follow universal curves for E A vs T V , as a function of ρ. Finally, this study achieves the first of these comparisons of molecular simulations to experiments and reveals critical insights toward creating vitrimers by design, while providing a route for the prediction of T V from kinetic Monte Carlo molecular dynamics simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Rational Enzyme Design for Controlled Functionalization of Acetylated Xylan for Cell-Free Polymer Biosynthesis

Xylan O-acetyltransferase 1 (XOAT1) is involved in O-acetylating the backbone of hemicellulose xylan. Recent structural analysis of XOAT1 showed two unequal lobes forming a cleft that is predicted to accommodate and position xylan acceptors into proximity with the catalytic triad. Here, we used docking and molecular dynamics simulations to investigate the optimal orientation of xylan in the binding cleft of XOAT1 and identify putative key residues (Gln445 and Arg444 on Minor lobe & Asn312, Met311 and Asp403 on Major lobe) involved in substrate interactions. Site-directed mutagenesis coupled with biochemical analyses revealed the major lobe of XOAT1 is important for xylan binding. Mutation of single key residues yielded XOAT1 variants with various enzymatic efficiencies that are applicable to one-pot synthesis of xylan polymers with different degrees of O-acetylation. Taken together, our results demonstrate the effectiveness of computational modeling in guiding enzyme engineering aimed at modulating xylan and redesigning plant cell walls.

93 SC - Biological and Environmental Research (BER↗

DOE OSTI2021

Topological Frustration as a New Parameter to Tune Morphology Revealed through Exploring the Continuum between A-B-C 3-Arm Star and Linear Triblock Polymers

Block polymers assemble into a variety of phase-separated morphologies based on volume fraction (φ) and interactions (χ) of the respective blocks. The arrangement of three different polymer blocks could either be a 3-arm star, with each block having one terminus attached to a common junction point or a linear A-B-C architecture. A versatile strategy is reported to synthesize a series of well-defined graft polymers that lie along the unexplored continuum between a 3-arm star and an A-B-C linear triblock polymer architecture. Using the technique of single-molecule insertion, precise control over the position of graft arm C along the B chain was achieved. A series of discrete graft polymers (PMMA-b-PS-g-PEO) with fixed φ and prescribed ω values that lie on the continuum between a 3-arm star (ω = 0) and linear triblock polymer (ω = 1) were synthesized. Morphological studies using small-angle X-ray scattering and conventional and energy-filtered transmission electron microscopy reveal the transition between lamellae, perforated lamellae, and cylindrical morphologies with systematic variation in the ω values, a trend attributed to the topological frustration and the associated χ values between the three blocks. Molecular dynamic simulations of coarse-grained models were found to predict phase diagrams that are consistent with the experimentally observed morphologies. Furthermore, our results suggest that changes in ω lead to topological frustration which is an important additional new design parameter that can be used to tune the morphology of multiblock polymers in addition to φ and χ.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

DL_POLY Quantum 2.0: A modular general-purpose software for advanced path integral simulations

DL_POLY Quantum 2.0, a vastly expanded software based on DL_POLY Classic 1.10, is a highly parallelized computational suite written in FORTRAN77 with a modular structure for incorporating nuclear quantum effects into large-scale/long-time molecular dynamics simulations. This is achieved by presenting users with a wide selection of state-of-the-art dynamics methods that utilize the isomorphism between a classical ring polymer and Feynman’s path integral formalism of quantum mechanics. Here, the flexible and user-friendly input/output handling system allows the control of methodology, integration schemes, and thermostatting. DL_POLY Quantum is equipped with a module specifically assigned for calculating correlation functions and printing out the values for sought-after quantities, such as dipole moments and center-of-mass velocities, with packaged tools for calculating infrared absorption spectra and diffusion coefficients.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Deep learning potential of mean force between polymer grafted nanoparticles

Grafting polymer chains on the surfaces of nanoparticles is a well-known route to control their self-assembly and distribution in a polymer matrix. A wide variety of self-assembled structures are achieved by changing the grafting patterns on the surface of an individual nanoparticle. However, an accurate estimation of the effective potential of mean force between a pair of grafted nanoparticles that determines their assembly and distribution in a polymer matrix is an outstanding challenge in nanoscience. We address this problem via deep learning. As a proof of concept, here we report a deep learning framework that learns the interaction between a pair of single-chain grafted spherical nanoparticles from their molecular dynamics trajectory. Subsequently, we carry out the deep learning potential of mean force-based molecular simulation that predicts the self-assembly of a large number of single-chain grafted nanoparticles into various anisotropic superstructures, including percolating networks and bilayers depending on the nanoparticle concentration in three-dimensions. The deep learning potential of mean force-predicted self-assembled superstructures are consistent with the actual superstructures of single-chain polymer grafted spherical nanoparticles. This deep learning framework is very generic and extensible to more complex systems including multiple-chain grafted nanoparticles. Finally, we expect that this deep learning approach will accelerate the characterization and prediction of the self-assembly and phase behaviour of polymer-grafted and unfunctionalized nanoparticles in free space or a polymer matrix.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Metal–Ligand Based Mechanophores Enhance Both Mechanical Robustness and Electronic Performance of Polymer Semiconductors

The backbone of diketopyrrolopyrrole-thiophene-vinylene-thiophene-based polymer semiconductors (PSCs) is modified with pyridine (Py) or bipyridine ligands to complex Fe(II) metal centers, allowing the metal–ligand complexes to act as mechanophores and dynamically crosslink the polymer chains. Here, mono- and bi-dentate ligands are observed to exhibit different degrees of bond strengths, which subsequently affect the mechanical properties of these Wolf-type-II metallopolymers. The counter ion also plays a crucial role, as it is observed that Py-Fe mechanophores with non-coordinating BPh 4 – counter ions (Py-FeB) exhibit better thin film ductility with lower elastic modulus, as compared to the coordinating chloro ligands (Py-FeC). Interestingly, besides mechanical robustness, the electrical charge carrier mobility can also be enhanced concurrently when incorporating Py-FeB mechanophores in PSCs. This is a unique observation among stretchable PSCs, especially that most reports to date describe a decreased mobility when the stretchability is improved. Next, it is determined that improvements to both mobility and stretchability are correlated to the solid-state molecular ordering and dynamics of coordination bonds under strain, as elucidated via techniques of grazing-incidence X-ray diffraction and X-ray absorption spectroscopy techniques, respectively. This study provides a viable approach to enhance both the mechanical and the electronic performance of polymer-based soft devices.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Simulations of Glass Transition and Mechanical Behavior of Off-Stoichiometric Crosslinked Polymers

Here, this work explores the influence of blend composition, network architecture, and hydrogen bonding on the material properties of crosslinked epoxy networks, focusing on the glass transition temperature (T g ) and Young’s modulus (Y). We used coarse-grained molecular dynamics simulations to simulate varying compositions of stiff and flexible components in epoxy monomer blends with varying excess of curative. We find that, without hydrogen bonding, networks of any composition show a monotonically increasing T g with decreasing excess curative, consistent with theory. In contrast, we find that when hydrogen bonding is introduced, the binary blend networks show significant enhancement in T g for lightly crosslinked systems. This result contributes to an explanation of the anomalous T g behavior observed experimentally in these systems. We further find that Y is generally enhanced by hydrogen bonds, especially below T g , demonstrating that hydrogen bonding has a significant influence on mechanical properties and can allow access to other desirable dynamic behavior, especially self-healing.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRSDate not supplied

Multiscale Modeling of Thermoplastics Using Atomistic-informed Micromechanics

A multiscale model was developed for predicting the thermoelastic behavior of semi-crystalline thermoplastic materials for composite aerospace applications. At the highest scale containing the semi-crystalline spherulite in an amorphous matrix, the generalized method of cells, or high fidelity method of cells, was used to perform the homogenization calculations to obtain the effective properties. Models were developed assuming a cubic, or spherical shape, for the spherulite to understand if the morphology of the spherulite affects the effective thermoelastic properties. The generalized method of cells was used to model at the repeating unit cells at the subscales of the microstructure including the lamellae stacks and granular crystal blocks. The scales are integrated using the multiscale micromechanics method in the NASA Multiscale Analysis Tool. Data from molecular dynamics simulations were used as inputs for the amorphous and crystalline constituents. Convergence studies were performed to determine the best level of discretization for the repeating unit cell at the highest scale. Effective Young’s modulus, shear modulus, Poisson’s ratio, coefficient of thermal expansion, and thermal conductivity were predicted for polyether ether ketone and polyether ketone ketone, and very good agreement between the model utilizing the cubic spherulite and the experimental data, where available, was observed for polyether ketone ketone. Normalization of the data for the bulk polyether ketone ketone, against amorphous data, improved the predictions as compared to experimental data. Overall, the high fidelity method of cells predicted a stiffer response then the generalized method of cells as the crystallinity was increased. The shape of the spherulite had a minimal effect on the predicted bulk properties of the polymers.

thermoplastics↗

DOE OSTI2020

Temperature dependence of nuclear quantum effects on liquid water via artificial neural network model based on SCAN meta-GGA functional

We investigate the temperature dependence of nuclear quantum effects (NQEs) on structural and dynamic properties of liquid water by training a neural network force field using first-principles molecular dynamics (FPMD) based on the strongly constrained and appropriately normed meta-generalized gradient approximation exchange-correlation approximation. The FPMD simulation based on density functional theory has become a powerful computational approach for studying a wide range of condensed phase systems. However, its large computational cost makes it difficult to incorporate NQEs in the simulation and investigate temperature dependence of various properties. To circumvent this difficulty, we use an artificial neural network model and employ the thermostatted ring polymer MD approach for studying the temperature dependence of NQEs on various properties. The NQEs generally bring the radial distribution functions closer to the experimental measurements. Translational diffusivity and rotational dynamics of water molecules are both slowed down by the NQEs. The competing inter-molecular and intra-molecular quantum effects on hydrogen bonds, as discussed by Habershon, Markland, and Manolopoulos [J. Chem. Phys. 131(2), 024501 (2019)], can explain the observed temperature dependence of the NQEs on the dynamical properties in our simulation.

Yao, Yi (ORCID:0000000299576979)↗

DOE OSTI2022

Chain Redistribution Stabilizes Coexistence Phases in Block Copolymer Blends

The nanoscale morphologies of block copolymer (BCP) thin films are determined by chain architecture. Experimental studies of thin film blends of different BCP chain types have demonstrated that blending can stabilize new motifs, such as coexistence phases. Here, we deploy coarsegrained molecular dynamics (MD) simulations in order to better understand the self-assembly behavior of BCP blend thin films. Here, we consider blends of lamella- and cylinder-forming BCP chains, studying their morphological makeup, the chain distribution within the morphology, and the underlying polymer chain conformations. Our simulations show that there are local concentration deviations at the scale of the morphological objects that dictate the local structure, and that BCP chains redistribute within the morphology so as to stabilize the structure. Underlying these effects are measurable distortions in the BCP chain conformations. The conformational freedom afforded by BCP blending stabilizes defects and allows coexistence phases to appear, while also leading to kinetic trapping effects. These results highlight the power of blending in designing the morphology that forms.

36 MATERIALS SCIENCE↗