Structural and conformational dynamics of supercooled polymer melts: insights from first-principles theory and simulations

Dynamics of Acrylamide Hydrogels, Polymers, and Monomers in Water Measured with Optical Heterodyne-Detected Optical Kerr Effect Spectroscopy

The ultrafast dynamics of acrylamide monomers (AAm), polyacrylamide (PAAm), and polyacrylamide hydrogels (PAAm-HG) in water were studied using optical heterodyne-detected optical Kerr effect (OHD-OKE) spectroscopy. Previous ultrafast infrared (IR) measurements of the water dynamics showed that at the same concentration of the acrylamide moiety, AAm, PAAm, and PAAm-HG exhibited identical water dynamics and that these dynamics slowed with increasing concentration. In contrast to the IR measurements, OHD-OKE experiments measure the dynamics of both the water and the acrylamide species, which occur on different time scales. In this study, the dynamics of all the acrylamide systems slowed with increasing concentration. We found that AAm exhibits tetraexponential decays, the longest component of which followed Debye-Stokes-Einstein behavior except for the highest concentration, 40% (w/v). Low concentrations of PAAm followed a single power law decay, while high concentrations of PAAm and all concentrations of PAAm-HG decayed with two power laws. The highest concentrations, 25% and 40%, of PAAm and PAAm-HG showed nearly identical dynamics. We interpreted this result as reflecting a similar extent of chain-chain interactions. At low concentrations, PAAm displays non-Markovian, single-chain dynamics (single power law), but PAAm displays entangled chain-chain interactions at high concentrations (two power laws). PAAm-HG has chain-chain interactions at all concentrations that arise from the cross-linking. At high concentrations, the dynamics of the entangled of PAAm become identical within error as those of the cross-linked PAAm-HG.

Multi-component generalized mode-coupling theory: predicting dynamics from structure in glassy mixtures

The emergence of glassy dynamics and the glass transition in dense disordered systems is still not fully understood theoretically. Mode-coupling theory (MCT) has shown to be effective in describing some of the non-trivial features of glass formation, but it cannot explain the full glassy phenomenology due to the strong approximations on which it is based. Generalized mode-coupling theory (GMCT) is a hierarchical extension of the theory, which is able to outclass MCT by carefully describing the dynamics of higher-order correlations in its generalized framework. Unfortunately, the theory has so far only been developed for single-component systems and as a result works poorly for highly polydisperse materials. In this paper, we solve this problem by developing GMCT for multi-component systems. We use it to predict the glassy dynamics of the binary Kob-Andersen Lennard-Jones mixture, as well as its purely repulsive Weeks-Chandler-Andersen analogue. Our results show that each additional level of the GMCT hierarchy gradually improves the predictive power of GMCT beyond its previous limit. This implies that our theory is able to harvest more information from the static correlations, thus being able to better understand the role of attraction in supercooled liquids from a first-principles perspective.

Glassy dynamics of a binary Voronoi fluid: a mode-coupling analysis

The binary Voronoi mixture is a fluid model whose interactions are derived from the Voronoi-Laguerre tessellation of the configurations of the system. The resulting interactions are local and many-body. Here we perform molecular-dynamics (MD) simulations of an equimolar mixture that is weakly polydisperse and additive. For the first time we study the structural relaxation of this mixture in the supercooled-liquid regime. From the simulations we determine the time- and temperature-dependent coherent and incoherent scattering functions for a large range of wave vectors, as well as the mean-square displacements of both particle species. We perform a detailed analysis of the dynamics by comparing the MD results with the first-principles-based idealized mode-coupling theory (MCT). To this end, we employ two approaches: fits to the asymptotic predictions of the theory, and fit-parameter-free binary MCT calculations based on static-structure-factor input from the simulations. We find that many-body interactions of the Voronoi mixture do not lead to strong qualitative differences relative to similar analyses carried out for simple liquids with pair-wise interactions. For instance, the fits give an exponent parameter λ ≈ 0.746 comparable to typical values found for simple liquids, the wavevector dependence of the Kohlrausch relaxation time is in good qualitative agreement with literature results for polydisperse hard spheres, and the MCT calculations based on static input overestimate the critical temperature, albeit only by a factor of about 1.2. This overestimation appears to be weak relative to other well-studied supercooled-liquid models such as the binary Kob-Andersen Lennard-Jones mixture. Overall, the agreement between MCT and simulation suggests that it is possible to predict several microscopic dynamic properties with qualitative, and in some cases near-quantitative, accuracy based solely on static two-point structural correlations, even though the system itself is inherently governed by many-body interactions.

Multiple character of non-monotonic size-dependence for relaxation dynamics in polymer-particle and binary mixtures

Adding plasticizers is a well-known procedure to reduce the glass transition temperature in polymers. It has been recently shown that this effect shows a non-monotonic dependence on the size of additive molecules (2019 J. Chem. Phys. 150 024903). In this work, we demonstrate that, as the size of the additive molecules is changed at fixed concentration, multiple extrema emerge in the dependence of the system's relaxation time on the size ratio. The effect occurs on all relevant length scales including single monomer dynamics, decay of Rouse modes and relaxation of the chain's end-to-end vector. A qualitatively similar trend is found within mode-coupling theoretical results for a binary hard-sphere mixture. An interpretation of the effect in terms of local packing efficiency and coupling between the dynamics of minority and majority species is provided.

Understanding, predicting, and tuning the fragility of vitrimeric polymers

Fragility is an empirical property that describes how abruptly a glass-forming material solidifies upon supercooling. The degree of fragility carries important implications for the functionality and processability of a material, as well as for our fundamental understanding of the glass transition. However, the microstructural properties underlying fragility still remain poorly understood. Here, we explain the microstructure-fragility link in vitrimeric networks, a novel type of high-performance polymers with unique bond-swapping functionality and unusual glass-forming behavior. Our results are gained from coarse-grained computer simulations and first-principles mode-coupling theory (MCT) of star-polymer vitrimers. We first demonstrate that the vitrimer fragility can be tuned over an unprecedentedly broad range, from fragile to strong and even superstrong behavior, by decreasing the bulk density. Remarkably, this entire phenomenology can be reproduced by microscopic MCT, thus challenging the conventional belief that existing first-principles theories cannot account for nonfragile behaviors. Our MCT analysis allows us to rationally identify the microstructural origin of the fragile-to-superstrong crossover, which is rooted in the sensitivity of the static structure factor to temperature variations. On the molecular scale, this behavior stems from a change in dominant length scales, switching from repulsive excluded-volume interactions to intrachain attractions as the vitrimer density decreases. Finally, we develop a simplified schematic MCT model which corroborates our microscopically founded conclusions and which unites our findings with earlier MCT studies. Our work sheds additional light on the elusive structure-fragility link in glass-forming matter and provides a first-principles-based platform for designing amorphous materials with an on-demand dynamic response.

The unique role of bond length in the glassy dynamics of colloidal polymers

Bond length is generally not considered as a controllable variable for molecular polymers. Hence, no experimental, simulation or theoretical research, to our knowledge, has examined the influence of bond length on the glassy dynamics of polymers. Recently, a new class of assembling materials called "colloidal polymers" has been synthesized. These colloidal polymers have advantages over molecular polymers in the visibility and flexibility of tuning, for example, the size and shape of the "monomers", the interaction, and the bond length. Dense suspension of colloidal polymers will become a very promising ideal model system for exploring the fundamental problems in the glass transition of chain "molecules". Here, we study the static structure and activated dynamics of hard-sphere colloidal polymers by generalizing the colloidal nonlinear Langevin equation theory to colloidal polymers. Surprisingly, we find that the bond length plays a critical and unique role in many aspects. For instance, the universal relations of the characteristic local lengths and the activated barrier versus the "degree of supercooling", and the structural relaxation versus local vibrational motion are found to be dependent on bond length and independent of chain length and rigidity. We hope that our findings inspire future experimental and simulation research studies on the glassy dynamics of colloidal polymers.

Are polymers standard glass-forming systems? The role of intramolecular barriers on the glass-transition phenomena of glass-forming polymers

Traditionally, polymer melts have been considered archetypal glass-formers. This has been mainly due to the fact that these systems can easily be obtained as glasses by cooling from the melt, even at low cooling rates. However, the macromolecules, i.e. the structural units of polymer systems in general, are rather different from the standard molecules. They are long objects ('chains') made by repetition of a given chemical motif (monomer) and have intra-macromolecular barriers that limit their flexibility. The influence of these properties on, for instance, the glass-transition temperature of polymers, is a topic that has been widely studied by the polymer community almost from the early times of polymer science. However, in the framework of the glass-community, the relevant influence of intra-macromolecular barriers and chain connectivity on glass-transition phenomena of polymers has started to be recognized only recently. The aim of this review is to give an overview and to critically revise the results reported on this topic over the last years. From these results, it seems to be evident that there are two different mechanisms involved in the dynamic arrest in glass-forming polymers: (i) the intermolecular packing effects, which dominate the dynamic arrest of low molecular weight glass-forming systems; and (ii) the effect of intra-macromolecular barriers combined with chain connectivity. It has also been shown that the mode coupling theory (MCT) is a suitable theoretical framework to discuss these questions. The values found for polymers for the central MCT parameter--the so-called λ-exponent--are of the order of 0.9, clearly higher than the standard values (λ ≈ 0.7) found in systems where the dynamic arrest is mainly driven by packing effects ('standard' glass-formers). Within the MCT, this is a signature of the presence of two competing mechanisms of dynamic arrest, as it has been observed in short-ranged attractive colloids or two component mixtures with dynamic asymmetry. Moreover, recent MD-simulations of a 'bead-spring' polymer model, but including intra-macromolecular potential of different strengths, confirm that the high λ-values found in polymers are due to the effect of intra-macromolecular barriers. Although there are still open questions, these results allow to conclude that there is a fundamental difference between the nature of the glass transition in polymers and in simple (standard) glass-formers.

Simulated glass-forming polymer melts: dynamic scattering functions, chain length effects, and mode-coupling theory analysis

We present molecular-dynamics simulations for a fully flexible model of polymer melts with different chain length N ranging from short oligomers (N = 4) to values near the entanglement length (N = 64). For these systems we explore the structural relaxation of the supercooled melt near the critical temperature T c of mode-coupling theory (MCT). Coherent and incoherent scattering functions are analyzed in terms of the idealized MCT. For temperatures T > T c we provide evidence for the space-time factorization property of the β relaxation and for the time-temperature superposition principle (TTSP) of the α relaxation, and we also discuss deviations from these predictions for T ≈ T c. For T larger than the smallest temperature where the TTSP holds we perform a quantitative analysis of the dynamics with the asymptotic MCT predictions for the late β regime. Within MCT a key quantity, in addition to T c, is the exponent parameter λ. For the fully flexible polymer models studied we find that λ is independent of N and has a value (λ = 0.735 ) typical of simple glass-forming liquids. On the other hand, the critical temperature increases with chain length toward an asymptotic value T c (∞) . This increase can be described by T c (∞) - T c(N) ∼ 1/N and may be interpreted in terms of the N dependence of the monomer density ρ, if we assume that the MCT glass transition is ruled by a soft-sphere-like constant coupling parameter Γ c = ρ c T c (-1/4), where ρ c is the monomer density at T c. In addition, we also estimate T c from a Hansen-Verlet-like criterion and MCT calculations based on structural input from the simulation. For our polymer model both the Hansen-Verlet criterion and the MCT calculations suggest T c to decrease with increasing chain length, in contrast to the direct analysis of the simulation data.

Continuous-time random-walk approach to supercooled liquids. II. Mean-square displacements in polymer melts

The continuous-time random walk (CTRW) describes the single-particle dynamics as a series of jumps separated by random waiting times. This description is applied to analyze trajectories from molecular dynamics (MD) simulations of a supercooled polymer melt. Based on the algorithm presented by Helfferich et al. [Phys. Rev. E 89, 042603 (2014)], we detect jump events of the monomers. As a function of temperature and chain length, we examine key distributions of the CTRW: the jump-length distribution (JLD), the waiting-time distribution (WTD), and the persistence-time distribution (PTD), i.e., the distribution of waiting times for the first jump. For the equilibrium (polymer) liquid under consideration, we verify that the PTD is determined by the WTD. For the mean-square displacement (MSD) of a monomer, the results for the CTRW model are compared with the underlying MD data. The MD data exhibit two regimes of subdiffusive behavior, one for the early α process and another at later times due to chain connectivity. By contrast, the analytical solution of the CTRW yields diffusive behavior for the MSD at all times. Empirically, we can account for the effect of chain connectivity in Monte Carlo simulations of the CTRW. The results of these simulations are then in good agreement with the MD data in the connectivity-dominated regime, but not in the early α regime where they systematically underestimate the MSD from the MD.

Applicability of mode-coupling theory to polyisobutylene: a molecular dynamics simulation study

The applicability of Mode Coupling Theory (MCT) to the glass-forming polymer polyisobutylene (PIB) has been explored by using fully atomistic molecular dynamics simulations. MCT predictions for the so-called asymptotic regime have been successfully tested on the dynamic structure factor and the self-correlation function of PIB main-chain carbons calculated from the simulated cell. The factorization theorem and the time-temperature superposition principle are satisfied. A consistent fitting procedure of the simulation data to the MCT asymptotic power-laws predicted for the α-relaxation regime has delivered the dynamic exponents of the theory-in particular, the exponent parameter λ-the critical non-ergodicity parameters, and the critical temperature T(c). The obtained values of λ and T(c) agree, within the uncertainties involved in both studies, with those deduced from depolarized light scattering experiments [A. Kisliuk et al., J. Polym. Sci. Part B: Polym. Phys. 38, 2785 (2000)]. Both, λ and T(c)/T(g) values found for PIB are unusually large with respect to those commonly obtained in low molecular weight systems. Moreover, the high T(c)/T(g) value is compatible with a certain correlation of this parameter with the fragility in Angell's classification. Conversely, the value of λ is close to that reported for real polymers, simulated "realistic" polymers and simple polymer models with intramolecular barriers. In the framework of the MCT, such finding should be the signature of two different mechanisms for the glass-transition in real polymers: intermolecular packing and intramolecular barriers combined with chain connectivity.

Hydrodynamic and viscoelastic effects in polymer diffusion

We develop a fluctuating hydrodynamics approach to study the impact of the hydrodynamic and viscoelastic interactions on the motion of the center of mass of a polymer as well as on the relaxation of Rouse modes, either in a Θ solvent or in a melt of identical unentangled chains. We show that this method allows us to describe the effect of hydrodynamic interactions beyond the Zimm (for a single chain in a Θ solvent) or the Rouse models (for an unentangled melt). In the latter case, we recover the same important effect of the viscoelastic hydrodynamic interactions on the center-of-mass diffusion, first described in Farago et al (2011 Phys. Rev. Lett. 107 178301).

Mode-coupling approach to polymer diffusion in an unentangled melt. I. The effect of density fluctuations

We quantitatively assess the effect of density fluctuation modes on the dynamics of a tagged polymer in an unentangled melt. To this end, we develop a density-based mode-coupling theory (dMCT) using the Mori-Zwanzig approach and projecting the fluctuating force onto pair-density fluctuation modes. The effect of dynamical density fluctuations on the center-of-mass (c.m.) dynamics is also analyzed based on a perturbative approach and we show that dMCT and perturbation techniques yield identical results. The c.m. velocity autocorrelation function (c.m. VAF) exhibits a slow power law relaxation in the time range between the monomer time t_{1} and the Rouse relaxation time t_{N}. We obtain an analytical expression for the c.m. VAF in terms of molecular parameters. In particular, the c.m. VAF scales as -N^{-1}t^{-5/4} (where N is the number of monomer units per chain) in the relevant time regime. The results are qualitatively accounted for by the dynamical correlation hole effect. The predicted -t^{-5/4} dependence of the c.m. VAF is supported by data of non-momentum-conserving computer simulations. However, the comparison shows that the theory significantly underestimates the amplitude of the effect. This issue is discussed and an alternative approach is addressed in the second part of this series [Farago et al., Phys. Rev. E 85, 051807 (2012), the following paper].

Anomalous diffusion of a polymer chain in an unentangled melt

Contrary to common belief, hydrodynamic interactions in polymer melts are not screened beyond the monomer length and are important in transient regimes. We show that viscoelastic hydrodynamic interactions (VHIs) lead to anomalous dynamics of a tagged chain in an unentangled melt at t<t(N), the Rouse time. The center-of-mass (c.m.) mean-square displacement is enhanced by a large factor increasing with chain length. We develop a theory of VHI-controlled chain dynamics yielding a negative c.m. velocity autocorrelation function which agrees with our molecular dynamics simulations without any fitting parameter. It is also shown that Langevin friction strongly affects the short-t c.m. dynamics, also captured by our theory. The transient VHI effects thus provide the dominant contribution to the subdiffusive c.m. motion universally observed in simulations and experiments.

First-principle approach to rescale the dynamics of simulated coarse-grained macromolecular liquids

We present a detailed derivation and testing of our approach to rescale the dynamics of mesoscale simulations of coarse-grained polymer melts (I. Y. Lyubimov, J. McCarty, A. Clark, and M. G. Guenza, J. Chem. Phys. 132, 224903 (2010)). Starting from the first-principle Liouville equation and applying the Mori-Zwanzig projection operator technique, we derive the generalized Langevin equations (GLEs) for the coarse-grained representations of the liquid. The chosen slow variables in the projection operators define the length scale of coarse graining. Each polymer is represented at two levels of coarse graining: monomeric as a bead-and-spring model and molecular as a soft colloid. In the long-time regime where the center-of-mass follows Brownian motion and the internal dynamics is completely relaxed, the two descriptions must be equivalent. By enforcing this formal relation we derive from the GLEs the analytical rescaling factors to be applied to dynamical data in the coarse-grained representation to recover the monomeric description. Change in entropy and change in friction are the two corrections to be accounted for to compensate the effects of coarse graining on the polymer dynamics. The solution of the memory functions in the coarse-grained representations provides the dynamical rescaling of the friction coefficient. The calculation of the internal degrees of freedom provides the correction of the change in entropy due to coarse graining. The resulting rescaling formalism is a function of the coarse-grained model and thermodynamic parameters of the system simulated. The rescaled dynamics obtained from mesoscale simulations of polyethylene, represented as soft-colloidal particles, by applying our rescaling approach shows a good agreement with data of translational diffusion measured experimentally and from simulations. The proposed method is used to predict self-diffusion coefficients of new polyethylene samples.

Structural relaxation of polydisperse hard spheres: comparison of the mode-coupling theory to a Langevin dynamics simulation

We analyze the slow glassy structural relaxation as measured through collective and tagged-particle density correlation functions obtained from Brownian dynamics simulations for a polydisperse system of quasi-hard spheres in the framework of the mode-coupling theory (MCT) of the glass transition. Asymptotic analyses show good agreement for the collective dynamics when polydispersity effects are taken into account in a multicomponent calculation, but qualitative disagreement at small q when the system is treated as effectively monodisperse. The origin of the different small-q behavior is attributed to the interplay between interdiffusion processes and structural relaxation. Numerical solutions of the MCT equations are obtained taking properly binned partial static structure factors from the simulations as input. Accounting for a shift in the critical density, the collective density correlation functions are well described by the theory at all densities investigated in the simulations, with quantitative agreement best around the maxima of the static structure factor and worst around its minima. A parameter-free comparison of the tagged-particle dynamics however reveals large quantitative errors for small wave numbers that are connected to the well-known decoupling of self-diffusion from structural relaxation and to dynamical heterogeneities. While deviations from MCT behavior are clearly seen in the tagged-particle quantities for densities close to and on the liquid side of the MCT glass transition, no such deviations are seen in the collective dynamics.

MD simulation of concentrated polymer solutions: structural relaxation near the glass transition

We examine by molecular dynamics simulations the relaxation of polymer-solvent mixtures close to the glass transition. The simulations employ a coarse-grained model in which polymers are represented by bead-spring chains and solvent particles by monomers. The interaction parameters between polymer and solvent are adjusted such that mixing is favored. We find that the mixtures have one glass transition temperature T(g) or critical temperature T(c) of mode-coupling theory (MCT). Both T(g) and T(c) (> T(g)) decrease with increasing solvent concentration φ(S). The decrease is linear for the concentrations studied (up to φ(S) = 25%). Above T(c) we explore the structure and relaxation of the polymer-solvent mixtures on cooling. We find that, if the polymer solution is compared to the pure polymer melt at the same T, local spatial correlations on the length scale of the first peak of the static structure factor S(q) are reduced. This difference between melt and solution is largely removed when comparing the S(q) of both systems at similar distance to the respective T(c). Near T(c) we investigate dynamic correlation functions, such as the incoherent intermediate scattering function φ(q)(s)(t), mean-square displacements of the monomers and solvent particles, two non-Gaussian parameters, and the probability distribution P(ln r; t) of the logarithm of single-particle displacements. In accordance with MCT we find, for instance, that φ(q)(s)(t) obeys the time-temperature superposition principle and has α relaxation times τ(q)(s) which are compatible with a power law increase close (but not too close) to T(c). In divergence to MCT, however, the increase of τ(q)(s) depends on the wavelength q, small q values having weaker increase than large ones. This decoupling of local and large-length scale relaxation could be related to the emergence of dynamic heterogeneity at low T. In the time window of the α relaxation an analysis of P(ln r; t) reveals a double-peak structure close to T(c). The first peak corresponds to "slow" particles (monomer or solvent) which have not moved much farther than 10% of their diameter in time t, whereas the second occurs at distances of the order of the particle diameter. These "fast" particles have succeeded in leaving their nearest-neighbor cage in time t. The simulation thus demonstrates that large fluctuations in particle mobility accompany the final structural relaxation of the cold polymer solution in the vicinity of the extrapolated T(c).

Subdiffusion and lateral diffusion coefficient of lipid atoms and molecules in phospholipid bilayers

We use a long, all-atom molecular-dynamics (MD) simulation combined with theoretical modeling to investigate the dynamics of selected lipid atoms and lipid molecules in a hydrated diyristoyl-phosphatidylcholine lipid bilayer. From the analysis of a 0.1 micros MD trajectory, we find that the time evolution of the mean-square displacement, <[deltar(t)]2>, of lipid atoms and molecules exhibits three well-separated dynamical regions: (i) ballistic, with <[deltar(t)]2> approximately t2 for t<or approximately 10 fs; (ii) subdiffusive, with <[deltar(t)]2> approximately tbeta with beta<1 for 10 ps<or approximately t<or approximately 10 ns; and (iii) Fickian diffusion, with <[deltar(t)]2> approximately t for t>or approximately 30 ns. We propose a memory-function approach for calculating <[deltar(t)]2> over the entire time range extending from the ballistic to the Fickian diffusion regimes. The results are in very good agreement with the ones from the MD simulations. We also examine the implications of the presence of the subdiffusive dynamics of lipids on the self-intermediate scattering function and the incoherent dynamic structure factor measured in neutron-scattering experiments.

Dynamic arrest in polymer melts: competition between packing and intramolecular barriers

We present molecular dynamics simulations of a simple model for polymer melts with intramolecular barriers. We investigate structural relaxation as a function of the barrier strength. Dynamic correlators can be consistently analyzed within the framework of the mode coupling theory of the glass transition. Control parameters are tuned in order to induce a competition between general packing effects and polymer-specific intramolecular barriers as mechanisms for dynamic arrest. This competition yields unusually large values of the so-called mode coupling theory exponent parameter and rationalizes qualitatively different observations for simple bead-spring and realistic polymers. The systematic study of the effect of intramolecular barriers presented here also establishes a fundamental difference between the nature of the glass transition in polymers and in simple glass formers.