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

Full‐dimensional coupled‐channel statistical approach to atom‐triatom systems and applications to H/D + O 3 reaction

Abstract The statistical quantum model (SQM), which assumes that the reactivity is controlled by entrance/exit channel quantum capture probabilities, is well suited for chemical reactions with a long‐lived intermediate complex. In this work, a time‐independent coupled‐channel implementation of the SQM approach is developed for atom‐triatom systems in full dimensionality. As SQM treats the capture dynamics quantum mechanically, it is capable of handling quantum effects such as tunneling. A detailed study of the H/D + O 3 capture dynamics was performed by applying the newly developed SQM method on an accurate global potential energy surface. Agreement with previous ring polymer molecular dynamics (RPMD) results on the same potential energy surface is excellent except for very low temperatures. The SQM results are also in reasonably good agreement with available experimental rate coefficients. The strong H/D kinetic isotope effect underscores the dominant role of quantum tunneling under an entrance channel barrier at low temperatures.

Yang, Dongzheng↗
DOE OSTI2022

Importance of Nuclear Quantum Effects on Aqueous Electrolyte Transport under Confinement in Ti 3 C 2 MXenes

Protons display a high chemical activity and strongly affect the charge storage capability in confined interlayer spaces of two-dimensional (2D) materials. As such, an accurate representation of proton dynamics under confinement is important for understanding and predicting charge storage dynamics in these materials. While often ignored in atomistic-scale simulations, nuclear quantum effects (NQEs), e.g., tunneling, can be significant under confinement even at room temperature. Here, using the thermostatted ring polymer molecular dynamics implementation of path integral molecular dynamics (PIMD) in conjunction with the ReaxFF force field, density functional tight binding (DFTB), and NequIP neural network potential simulations, we investigate the role of NQEs on proton and water transport in bulk water and aqueous electrolytes under confinement in Ti 3 C 2 MXenes. Although overall NQEs are relatively small, especially in bulk, we find that they can alter both quantitative values and qualitative trends on both proton transport and water self-diffusion under confinement relative to classical MD predictions. Therefore, our results suggest the need for NQEs to be considered to simulate aqueous systems under confinement for both qualitative and quantitative accuracy.

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

Microstructural characterization of a star-linear polymer blend under shear flow by using rheo-SANS

We present an investigation into the dynamic relaxation mechanisms of a polybutadiene blend composed of a four-arm star (10 wt. %) and a linear polymer matrix in the presence of an applied shear flow. Our focus was the response of the star polymer, which cannot be unambiguously assessed via linear viscoelastic measurements since the signature of the star polymer can barely be detected due to the dominant contribution of the linear matrix. By utilizing small-angle neutron scattering (SANS) coupled with a Couette shear device and a deuterated matrix polymer, we investigated the dynamics of the minority star component of the blend. Our results confirm that the stars deform anisotropically with increasing shear rate. We have compared the SANS data with predictions from the well-established scattering adaptation of the state-of-the-art tube model for entangled linear polymer melts undergoing shear, i.e., Graham, Likhtman, Milner, and McLeish (GLaMM) approach, appropriately modified following earlier studies in order to apply to the star. This modified model, GLaMM-R, includes the physics necessary to understand stress relaxation in both the linear and nonlinear flow regimes, i.e., contour length fluctuations, constraint release, convective constraint release, and chain retraction. The full scattering signal is due to the minority star component and, although the contribution of the linear chains is hidden from the neutron scattering, they still influence the star polymer molecular dynamics, with the applied shear rate ranging from approximately 8 to 24 s –1 , below the inverse relaxation time of the linear component. This study provides another confirmation that the combination of rheology and neutron scattering is an indispensable tool for investigating the nonlinear dynamics of complex polymeric systems.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Two-temperature warm dense hydrogen as a test of quantum protons driven by orbital-free density functional theory electronic forces

We consider a steady-state (but transient) situation in which a warm dense aggregate is a two-temperature system with equilibrium electrons at temperature Te, ions at Ti, and Te ≠ Ti. Such states are achievable by pump–probe experiments. For warm dense hydrogen in such a two-temperature situation, we investigate nuclear quantum effects (NQEs) on structure and thermodynamic properties, thereby delineating the limitations of ordinary ab initio molecular dynamics. We use path integral molecular dynamics (PIMD) simulations driven by orbital-free density functional theory (OFDFT) calculations with state-of-the-art noninteracting free-energy and exchange-correlation functionals for the explicit temperature dependence. We calibrate the OFDFT calculations against conventional (explicit orbitals) Kohn–Sham DFT. We find that when the ratio of the ionic thermal de Broglie wavelength to the mean interionic distance is larger than about 0.30, the ionic radial distribution function is meaningfully affected by the inclusion of NQEs. Moreover, NQEs induce a substantial increase in both the ionic and electronic pressures. This confirms the importance of NQEs for highly accurate equation-of-state data on highly driven hydrogen. For Te > 20 kK, increasing Te in the warm dense hydrogen has slight effects on the ionic radial distribution function and equation of state in the range of densities considered. In addition, we confirm that compared with thermostatted ring-polymer molecular dynamics, the primitive PIMD algorithm overestimates electronic pressures, a consequence of the overly localized ionic description from the primitive scheme.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Quantum Simulations of Vibrational Strong Coupling via Path Integrals

A quantum simulation of vibrational strong coupling (VSC) in the collective regime via thermostated ring-polymer molecular dynamics (TRPMD) is reported. For a collection of liquid-phase water molecules resonantly coupled to a single lossless cavity mode, the simulation shows that as compared with a fully classical calculation, the inclusion of nuclear and photonic quantum effects does not lead to a change in the Rabi splitting but does broaden polaritonic line widths roughly by a factor of 2. Moreover, under thermal equilibrium, both quantum and classical simulations predict that the static dielectric constant of liquid water is largely unchanged inside vs outside the cavity. Furthermore, this result disagrees with a recent experiment demonstrating that the static dielectric constant of liquid water can be resonantly enhanced under VSC, suggesting either limitations of our approach or perhaps other experimental factors that have not yet been explored.

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

Zero-cost corrections to influence functional coefficients from bath response functions

Recent work has shown that it is possible to circumvent the calculation of the spectral density and directly generate the coefficients of the discretized influence functionals using data from classical trajectory simulations. However, the accuracy of this procedure depends on the validity of the high temperature approximation. In this work, an alternative derivation based on the Kubo formalism is provided. This enables the calculation of additional correction terms that increases the range of applicability of the procedure to lower temperatures. Because it is based on the Kubo-transformed correlation function, this approach allows the direct use of correlation functions obtained from methods such as ring-polymer molecular dynamics and centroid molecular dynamics in determining the influence functional coefficients for subsequent system-solvent simulations. The accuracy of the original procedure and the corrected procedure is investigated across a range of parameters. It is interesting that the correction term comes at zero additional cost. Furthermore, it is possible to improve upon the correction using zero-cost physical intuition and heuristics making the method even more accurate.

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

Dynamics of poly[n]catenane melts

Inspired by advances in the chemical synthesis of interlocking polymer architectures, extensive molecular dynamics simulations have been conducted to study the dynamical properties of poly[n]catenanes-polymers composed entirely of interlocking rings-in the melt state. Both the degree of polymerization (number of links) and the number of beads per ring are systematically varied, and the results are compared to linear and ring polymers. Herein, a simple Rouse-like model is presented, and its analytical solution suggests a decomposition of the dynamics into "ring-like" and "linear-like" regimes at short and long times, respectively. In agreement with this picture, multiple sub-diffusive regimes are observed in the monomer mean-squared-displacements even though interchain entanglement is not prevalent in the system. However, the Rouse-type model does not account for the topological effects of the mechanical bonds, which significantly alter the dynamics at intermediate length scales both within the rings and at the chain segment scales. The stress relaxation in the system is extremely rapid and may be conveniently separated into ring-like and linear-like contributions, again in agreement with the Rouse picture. However, the viscosity has a non-monotonic dependence on the ring size for long chains, which disagrees strongly with theoretical predictions. This unexpected observation cannot be explained in terms of chain disentanglement and is inconsistent with other measures of polymer relaxation. Possible mechanisms for this behavior are proposed and implications for materials design are discussed.

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

Energy Transfer Mechanisms in Large Low-Bandgap Polymers from Time-Resolved Experiments and Nonadiabatic Molecular Dynamics Calculations

Conjugated polymers offer unprecedented chemical tunability for modulating energy transfer in a multitude of infrared light applications. In this work, we use a combination of time-resolved spectroscopic experiments and nonadiabatic molecular dynamics calculations to probe the photochemistry and nonradiative transitions in a recently synthesized narrow bandgap donor–acceptor conjugated polymer based on alternating cyclopentadithiophene and electronegative benzothiadiazole heterocycles. Using large-scale semi-empirical nonadiabatic molecular dynamics, which can treat a large 260-atom hexamer, we calculate an S 5 → S 1 lifetime of 34.75 fs, which is consistent with our time-resolved spectroscopic data. Our simulations suggest that vibronic motions of the central carbons in the cyclopentadithiophene functional groups are predominantly involved in the nonradiative transitions, and the excitation becomes more localized on a monomer fragment over time. The combined use of time-resolved experiments and nonadiabatic molecular dynamics calculations in this work provides mechanistic insight into chemical functionalities that can be tuned to enhance energy transfer in other prospective low-bandgap polymer materials.

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

Accelerating amorphous polymer electrolyte screening by learning to reduce errors in molecular dynamics simulated properties

Polymer electrolytes are promising candidates for the next generation lithium-ion battery technology. Large scale screening of polymer electrolytes is hindered by the significant cost of molecular dynamics (MD) simulation in amorphous systems: the amorphous structure of polymers requires multiple, repeated sampling to reduce noise and the slow relaxation requires long simulation time for convergence. Here, we accelerate the screening with a multi-task graph neural network that learns from a large amount of noisy, unconverged, short MD data and a small number of converged, long MD data. We achieve accurate predictions of 4 different converged properties and screen a space of 6247 polymers that is orders of magnitude larger than previous computational studies. Further, we extract several design principles for polymer electrolytes and provide an open dataset for the community. Our approach could be applicable to a broad class of material discovery problems that involve the simulation of complex, amorphous materials.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Hoobas: A highly object-oriented builder for molecular dynamics

Polydispersity and random sequences are ubiquitous features of polymers, and molecular dynamics simulations can help elucidate the impact of disorder in polymer systems. However, currently available packages for building polymer topologies do not enable the user to include randomness in a straightforward fashion. Here, we introduce Hoobas, a molecular builder package that easily handles polydispersity using a prototype-builder design pattern. This enables fast and easy building of systems comprised of thousands of distinct objects. It is written in the Python programming language, which ensures compatibility with a wide range of molecular dynamics packages and tools, as well as easy integration into most workflows.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

The Work of Mechanical Degradation in Elongating Polymer Melts

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

36 MATERIALS SCIENCE↗
DOE OSTI2021

Role of Low Molecular Weight Polymers on the Dynamics of Silicon Anodes During Casting

This work probes the slurry architecture of a high silicon content electrode slurry with and without low molecular weight polymeric dispersants as a function of shear rate to mimic electrode casting conditions for poly(acrylic acid) (PAA) and lithium neutralized poly(acrylic acid) (LiPAA) based electrodes. Rheology coupled ultra-small angle neutron scattering (rheo-USANS) was used to examine the aggregation and agglomeration behavior of each slurry as well as the overall shape of the aggregates. The addition of dispersant has opposing effects on slurries made with PAA or LiPAA binder. With a dispersant, there are fewer aggregates and agglomerates in the PAA based silicon slurries, while LiPAA based silicon slurries become orders of magnitude more aggregated and agglomerated at all shear rates. The reorganization of the PAA and LiPAA binder in the presence of dispersant leads to a more homogeneous slurry and a more heterogeneous slurry, respectively. This reorganization ripples through to the cast electrode architecture and is reflected in the electrochemical cycling of these electrodes.

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

Molecular Dynamics Investigation of Nanoscale Hydrophobicity of Polymer Surfaces: What Makes Water Wet?

The wettability of a polymer surface–related to its hydrophobicity or tendency to repel water–can be crucial for determining its utility, such as for a coating or a purification membrane. While wettability is commonly associated with the macroscopic measurement of a contact angle between surface, water, and air, the molecular physics that underlie these macroscopic observations are not fully known, and anticipating relative behavior of different polymers is challenging. To address this gap in molecular-level understanding, we use molecular dynamics simulations to investigate and contrast interactions of water with six chemically distinct polymers: polytetrafluoroethylene, polyethylene, polyvinyl chloride, poly(methyl methacrylate), Nylon 66, and polyvinyl alcohol. We show that several prospective quantitative metrics for hydrophobicity agree well with experimental contact angles. Moreover, the behavior of water in proximity to these polymer surfaces can be distinguished with analysis of interfacial water dynamics, extent of hydrogen bonding, and molecular orientation–even when macroscopic measures of hydrophobicity are similar. The predominant factor dictating wettability is found to be the extent of hydrogen bonding between polymer and water, but the precise manifestation of hydrogen bonding and its impact on surface-water structure varies. In the absence of hydrogen bonding, other molecular interactions and polymer mechanics control hydrophobic ordering. Furthermore, these results provide new insights into how polymer chemistry specifically impacts water-polymer interactions and translates to surface hydrophobicity. Such factors may facilitate the design or processing of polymer surfaces to achieve targeted wetting behavior, and presented analyses can be useful in studying the interfacial physics of other systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗