A comparison of neutron scattering studies and computer simulations of polymer melts

Complex Monte Carlo Light-Driven Dynamics of Monomers in Functionalized Bond Fluctuation Model Polymer Chains

We study Monte Carlo dynamics of the monomers and center of mass of a model polymer chain functionalized with azobenzene molecules in the presence of an inhomogeneous linearly polarized laser light. The simulations use a generalized Bond Fluctuation Model. The mean squared displacements of the monomers and the center of mass are analyzed in a period of Monte Carlo time typical for a build-up of Surface Relief Grating. Approximate scaling laws for mean squared displacements are found and interpreted in terms of sub- and superdiffusive dynamics for the monomers and center of mass. A counterintuitive effect is observed, where the monomers perform subdiffusive motion but the resulting motion of the center of mass is superdiffusive. This result disparages theoretical approaches based on an assumption that the dynamics of single monomers in a chain can be characterized in terms of independent identically distributed random variables.

Anomalous Dynamics in Macromolecular Liquids

Macromolecular liquids display short-time anomalous behaviors in disagreement with conventional single-molecule mean-field theories. In this study, we analyze the behavior of the simplest but most realistic macromolecular system that displays anomalous dynamics, i.e., a melt of short homopolymer chains, starting from molecular dynamics simulation trajectories. Our study sheds some light on the microscopic molecular mechanisms responsible for the observed anomalous behavior. The relevance of the correlation hole, a unique property of polymer liquids, in relation to the observed subdiffusive dynamics, naturally emerges from the analysis of the van Hove distribution functions and other properties.

Comparison of molecular dynamics simulations of water with neutron and X-ray scattering experiments

The atomistic structure and dynamics obtained from molecular dynamics (MD) simulations with the example of TIP3P (rigid and flexible) and TIP4P/2005 (rigid) water is compared to neutron and X-ray scattering data at ambient conditions. Neutron and X-ray diffractograms are calculated from the simulations for four isotopic substitutions as well as the incoherent intermediate scattering function for neutrons. The resulting curves are compared to each other and to published experimental data. Differences between simulated and measured intermediate scattering functions are quantified by fitting an analytic model to the computed values. The sensitivity of the scattering curves to the parameters of the MD simulations is demonstrated on the example of two parameters, bond length and angle.

Insight into the Structure and Dynamics of Polymers by Neutron Scattering Combined with Atomistic Molecular Dynamics Simulations

Combining neutron scattering and fully atomistic molecular dynamics simulations allows unraveling structural and dynamical features of polymer melts at different length scales, mainly in the intermolecular and monomeric range. Here we present the methodology developed by us and the results of its application during the last years in a variety of polymers. This methodology is based on two pillars: (i) both techniques cover approximately the same length and time scales and (ii) the classical van Hove formalism allows easily calculating the magnitudes measured by neutron scattering from the simulated atomic trajectories. By direct comparison with experimental results, the simulated cell is validated. Thereafter, the information of the simulations can be exploited, calculating magnitudes that are experimentally inaccessible or extending the parameters range beyond the experimental capabilities. We show how detailed microscopic insight on structural features and dynamical processes of various kinds has been gained in polymeric systems with different degrees of complexity, and how intriguing questions as the collective behavior at intermediate length scales have been faced.

A practical method to account for random phase approximation effects on the dynamic scattering of multi-component polymer systems

Investigations of polymer systems that rely on the interpretation of dynamical scattering results as, e.g., the structure factor S(Q, t) of single chains or chain sections may require the inclusion of effects, as described within the framework of the random phase approximation (RPA) for polymers. To do this in practice for the dynamic part of S(Q, t) beyond the initial slope is a challenge. Here, we present a method (and software) that allows a straightforward assessment of dynamical RPA effects and inclusion of these in the process/procedures of model fitting. Examples of applications to the interpretation of neutron spin-echo data multi-component polymer melts are shown.

Chain end mobilities in polymer melts--a computational study

The Rouse model can be regarded as the standard model to describe the dynamics of a short polymer chain under melt conditions. In this contribution, we explicitly check one of the fundamental assumptions of this model, namely, that of a uniform friction coefficient for all monomers, on the basis of MD simulation data of a poly(ethylene oxide) (PEO) melt. This question immediately arises from the fact that in a real polymer melt, the terminal monomers have on average more intermolecular neighbors than the central monomers, and one would expect that exactly these details affect the precise value of the friction coefficient. The mobilities are determined by our recently developed statistical method, which provides detailed insights into the local polymer dynamics. Moreover, it yields complementary information to that obtained from the mean square displacement (MSD) or the Rouse mode analysis. It turns out that the Rouse assumption of a uniform mobility is fulfilled to a good approximation for the PEO melt. However, a more detailed analysis reveals that the underlying microscopic dynamics are highly affected by different contributions from intra- and intermolecular excluded volume interactions, which cannot be taken into account by a modified friction coefficient. Minor deviations occur only for the terminal monomers on larger time scales, which can be attributed to the presence of two different escape mechanisms from their first coordination sphere. These effects remain elusive when studying the dynamics with the MSD only.

Generation of Well-Relaxed All-Atom Models of Large Molecular Weight Polymer Melts: A Hybrid Particle-Continuum Approach Based on Particle-Field Molecular Dynamics Simulations

A procedure based on Molecular Dynamics (MD) simulations employing soft potentials derived from self-consistent field (SCF) theory (named MD-SCF) able to generate well-relaxed all-atom structures of polymer melts is proposed. All-atom structures having structural correlations indistinguishable from ones obtained by long MD relaxations have been obtained for poly(methyl methacrylate) (PMMA) and poly(ethylene oxide) (PEO) melts. The proposed procedure leads to computational costs mainly related on system size rather than to the chain length. Several advantages of the proposed procedure over current coarse-graining/reverse mapping strategies are apparent. No parametrization is needed to generate relaxed structures of different polymers at different scales or resolutions. There is no need for special algorithms or back-mapping schemes to change the resolution of the models. This characteristic makes the procedure general and its extension to other polymer architectures straightforward. A similar procedure can be easily extended to the generation of all-atom structures of block copolymer melts and polymer nanocomposites.

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].

Mode-coupling approach to polymer diffusion in an unentangled melt. II. The effect of viscoelastic hydrodynamic interactions

A mode-coupling theory (MCT) version (called hMCT thereafter) of a recently presented theory [Farago, Meyer, and Semenov, Phys. Rev. Lett. 107, 178301 (2011)] is developed to describe the diffusional properties of a tagged polymer in a melt. The hMCT accounts for the effect of viscoelastic hydrodynamic interactions (VHIs), that is, a physical mechanism distinct from the density-based MCT (dMCT) described in the first paper of this series. The two versions of the MCT yield two different contributions to the asymptotic behavior of the center-of-mass velocity autocorrelation function (c.m. VAF). We show that in most cases the VHI mechanism is dominant; for long chains and prediffusive times it yields a negative tail ∝-N^{-1/2}t^{-3/2} for the c.m. VAF. The case of non-momentum-conserving dynamics (Langevin or Monte Carlo) is discussed as well. It generally displays a distinctive behavior with two successive relaxation stages: first -N^{-1}t^{-5/4} (as in the dMCT approach), then -N^{-1/2}t^{-3/2}. Both the amplitude and the duration of the first t^{-5/4} stage crucially depend on the Langevin friction parameter γ. All results are also relevant for the early time regime of entangled melts. These slow relaxations of the c.m. VAF, thus account for the anomalous subdiffusive regime of the c.m. mean square displacement widely observed in numerical and experimental works.

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.

Characterization of the "simple-liquid" state in a polymeric system: coherent and incoherent scattering functions

By means of time-of-flight neutron scattering, we have characterized the dynamic structure factor and the hydrogen motions of polyethylene above its melting point. As signatures of simple-liquid dynamics, we observe that (i) in the explored dynamic window, the intermediate scattering functions display a single-step decay, without reminiscence of cage effects, and (ii) the structural relaxation as observed at the intermolecular structure factor peak shows a weak temperature dependence, indicating a low degree of interchain cooperativity. However, stretched functional forms are observed for all length scales and temperatures investigated. An apparent direct crossover from the microscopic regime to Rouse-like dynamics suggests an essential role of connectivity in the observed stretching. Finally, while at momentum-transfer values above the structure factor peak the relation between coherent and incoherent characteristic times is reasonably described by the de Gennes narrowing, at larger length scales it is not reproduced by any existing approach.

Cooperative dynamics in homopolymer melts: a comparison of theoretical predictions with neutron spin echo experiments

We present a comparison between theoretical predictions of the generalized Langevin equation for cooperative dynamics (CDGLE) and neutron spin echo data of dynamic structure factors for polyethylene melts. Experiments cover an extended range of length and time scales, providing a compelling test for the theoretical approach. Samples investigated include chains with increasing molecular weights undergoing dynamics across the unentangled to entangled transition. Measured center-of-mass (com) mean-square displacements display a crossover from subdiffusive to diffusive dynamics. The generalized Langevin equation for cooperative dynamics relates this anomalous diffusion to the presence of the interpolymer potential, which correlates the dynamics of a group of slowly diffusing molecules in a dynamically heterogeneous liquid. Theoretical predictions of the subdiffusive behavior, of its crossover to free diffusion, and of the number of macromolecules undergoing cooperative motion are in quantitative agreement with experiments.

Glass transition in 1,4-polybutadiene: Mode-coupling theory analysis of molecular dynamics simulations using a chemically realistic model

We present molecular dynamics simulations of the glass transition in a chemically realistic model of 1,4-polybutadiene (PBD). Around 40 K above the calorimetric glass transition of this polymer the simulations reveal a well-developed two-stage relaxation of all correlation functions. We have analyzed the time-scale separation between vibrational degrees of freedom (subpicosecond dynamics) and the alpha relaxation behavior (nanosecond to microsecond dynamics) using the predictions of mode-coupling theory (MCT). Our value for the mode-coupling critical temperature Tc agrees perfectly with prior experimental estimates for PBD. The predictions of MCT for the scaling behavior of the so-called beta relaxation, the plateau regime separating vibrational dynamics and the alpha relaxation, are well fulfilled. Furthermore, we are able to derive a consistent set of MCT exponents, completely characterizing the scaling behavior of relaxation processes in the vicinity of Tc. For the temperature dependence of the alpha relaxation we find deviations from MCT predictions which we trace to the influence of torsional barriers on the atomic motions.

Molecular dynamics simulation study of the alpha- and beta-relaxation processes in a realistic model polymer

Molecular dynamics simulations of a melt of freely rotating chains of 1,4-polybutadiene (FRC-PBD) have been performed over a wide range of temperature. Removal of the dihedral barriers in FRC-PBD allows for complete resolution of the Johari-Goldstein beta process from the primary alpha process in the simulation time window. We find that relaxation in the beta regime occurs as the result of large-angle excursions of all backbone dihedrals. While during the beta process each dihedral visits all available states, the complete relaxation of torsional autocorrelation function (alpha process) occurs only when the polymer matrix shows significant motion and each dihedral populates all states with equilibrium probability.

A molecular dynamics simulation study of the α-relaxation in a 1,4-polybutadiene melt as probed by the coherent dynamic structure factor

The dynamic coherent structure factor Scoh(q,t) for a 1,4-polybutadiene (PBD) melt has been investigated using atomistic molecular dynamics simulations. The relaxation of Scoh(q,t) at q = 1.44 angstroms(-1) and q = 2.72 angstroms(-1), corresponding to the first and second peaks in the static structure factor for PBD, was studied in detail over a wide range of temperature. It was found that time-temperature superposition holds for the alpha-relaxation for both q values over a wide temperature range and that the alpha-relaxation can be well described by a stretched (Kohlrauch-William-Watts) exponential with temperature independent but q dependent amplitude and stretching exponent. The alpha-relaxation times for both q values were found to exhibit the same non-Arrhenius temperature dependence, indicating that the same physical processes are responsible for relaxation on both length scales. The alpha-relaxation time was found to depend strongly upon the dynamical range of data utilized in determining the relaxation time, accounting for qualitative discrepancies between alpha-relaxation times reported here and those extracted for PBD from experimentally measured Scoh(q,t).

Anomalous pressure dependence of the structure factor in 1,4-polybutadiene melts: a molecular dynamics simulation study

Neutron scattering has shown the first diffraction peak in the structure factor of a 1,4-polybutadiene melt under compression to move to larger q values as expected but to decrease significantly in intensity. Simulations reveal that this behavior does not result from loss of structure in the polymer melt upon compression but rather from the generic effects of differences in the pressure dependence of the intermolecular and intramolecular contributions to the melt structure factor and differences in the pressure dependence of the partial structure factors for carbon-carbon and carbon-deuterium intermolecular correlations. This anomalous pressure dependence is not seen for protonated melts.

Intramolecular caging in polybutadiene due to rotational barriers

We present molecular dynamics simulations of a chemically realistic model of 1,4-polybutadiene and a freely rotating chain model derived from the first model by neglecting all dihedral potentials. We show that the presence of energy barriers hindering dihedral rotation leads to an intermediate plateau regime in the tagged particle mean-squared displacement reminiscent of the cage effect underlying the mode-coupling description of the liquid-glass transition. This intramolecular caging, however, occurs already at temperatures well above the glass transition regime. Because of its different physical origin, it also does not comply with the theoretical predictions of the mode-coupling theory. Consequences for the applicability of the mode-coupling theory to the glass transition in polymer melts are discussed.