A molecular-dynamics simulation study of dielectric relaxation in a 1,4-polybutadiene melt

Understanding relaxation times in supercooled liquids and glasses

We formulate models of segmental relaxation in glasses and supercooled liquids in terms of Poisson processes and stochastic resetting and analyze their properties. This mathematical language allows us to set forth clear, consistent definitions for the various terms used throughout the literature on glassy dynamics. We provide useful algorithms for stochastic simulations of such models and show how they may be properly parameterized from time autocorrelation functions (ACFs), obtained by either simulation or experiment. Interestingly, we find that the time derivative of an ACF provides considerable insight into the distribution of relaxation times or rates and is, therefore, the primary object of analysis. These results allow for physically and mathematically robust modeling of segmental dynamics in molecular and polymeric glasses.

Structure and dynamics of a 1,4-polybutadiene melt in an alumina nanopore: A molecular dynamics simulation

We present results of Molecular Dynamics simulations of a chemically realistic model of 1,4-polybutadiene (PBD)confined in a cylindrical alumina nanopore of diameter 10 nm. The simulations are done at three different temperaturesabove the glass transition temperature T g . We investigate the density layering across the nanopore as well as theorientational ordering in the polymer melt, brought about by the confinement, on both the segmental and chain scales.For the chain scale ordering, the magnitude and orientation of the axes of the gyration tensor ellipsoid of single chainsare studied and are found to prefer to align parallel to the pore axis. Even though double bonds near the wall arepreferentially oriented along the pore walls, studying the nematic order parameter indicates that there is no nematicordering at the melt-wall interface. As for the dynamics in the melt, we focus here on the mean-square-displacement ofthe monomers for several layers across the nanopore as well as the movement of the chain center of mass which bothdisplay a slowing down of the dynamics in the layer at the wall. We also show the strong adsorption of the monomersto the pore wall at lower temperatures.

Role of torsional potential in chain conformation, thermodynamics, and glass formation of simulated polybutadiene melts

For polymer chains, the torsional potential is an important intramolecular energy influencing chain flexibility and segmental dynamics. Through molecular dynamics simulations of an atomistic model for melts of cis-trans-1,4-polybutadiene (PBD), we explore the effect of the torsions on conformational properties (bond vector correlations and mean-square internal distances), fundamental thermodynamic quantities (density, compressibility, internal energy, and specific heat), and glass transition temperature T_g. This is achieved by systematically reducing the strength of the torsional potential, starting from the chemically realistic chain (CRC) model with the full potential toward the freely rotating chain (FRC) model without the torsional potential. For the equilibrium liquid, we find that the effect of the torsions on polymer conformations is very weak. Still weaker is the influence on the monomer density ρ and isothermal compressibility κ_T of the polymer liquid, both of which can be considered as independent of the torsional potential. We show that a van der Waals-like model proposed by Long and Lequeux [Eur. Phys. J. E 4, 371 (2001)] allows us to describe very well the temperature (T) dependence of ρ and κ_T. We also find that our data obey the linear relation between 1/k_BTρκ_T and 1/T (with the Boltzmann constant k_B) that has recently been predicted and verified on the experiment by Mirigian and Schweizer [J. Chem. Phys. 140, 194507 (2014)]. For the equilibrium liquid, simulations result in a specific heat, at constant pressure and at constant volume, which increases on cooling. This T dependence is opposite to the one found experimentally for many polymer liquids, including PBD. We suggest that this difference between simulation and experiment may be attributed to quantum effects due to hydrogen atoms and backbone vibrations, which, by construction, are not included in the classical united-atom model employed here. Finally, we also determine T_g from the density-temperature curve monitored in a finite-rate cooling process. While the influence of the torsional potential on ρ(T) is vanishingly small in the equilibrium liquid, the effect of the torsions on T_g is large. We find that T_g decreases by about 150 K when going from the CRC to the FRC model.

Gradient of Segmental Dynamics in Stereoregular Poly(Methyl Methacrylate) Melts Confined between Pristine or Oxidized Graphene Sheets

Segmental dynamics in unentangled isotactic, syndiotactic, and atactic poly(methyl methacrylate) (i-, a-, and s-PMMA) melts confined between pristine graphene, reduced graphene oxide, RGO, or graphene oxide, GO, sheets is studied at various temperatures, well above glass transition temperature, via atomistic molecular dynamics simulations. The model RGO and GO sheets have different degrees of oxidization. The segmental dynamics is studied through the analysis of backbone torsional motions. In the vicinity of the model nanosheets (distances less than ≈2 nm), the dynamics slows down; the effect becomes significantly stronger with increasing the concentration of the surface functional groups, and hence increasing polymer/surface specific interactions. Upon decreasing temperature, the ratios of the interfacial segmental relaxation times to the respective bulk relaxation times increase, revealing the stronger temperature dependence of the interfacial segmental dynamics relative to the bulk dynamics. This heterogeneity in temperature dependence leads to the shortcoming of the time-temperature superposition principle for describing the segmental dynamics of the model confined melts. The alteration of the segmental dynamics at different distances, d, from the surfaces is described by a temperature shift, ΔTseg(d) (roughly speaking, shift of a characteristic temperature). Next, to a given nanosheet, i-PMMA has a larger value of ΔTseg than a-PMMA and s-PMMA. This trend correlates with the better interfacial packing and longer trains of i-PMMA chains. The backbone torsional autocorrelation functions are shown in the frequency domain and are qualitatively compared to the experimental dielectric loss spectra for the segmental α-relaxation in polymer nanocomposites. The εT″(f) (analogous of dielectric loss, ε″(f), for torsional motion) curves of the model confined melts are broader (toward lower frequencies) and have lower amplitudes relative to the corresponding bulk curves; however, the peak frequencies of the εT″(f) curves are only slightly affected.

Universal localization transition accompanying glass formation: insights from efficient molecular dynamics simulations of diverse supercooled liquids

The origin of the precipitous dynamic arrest known as the glass transition is a grand open question of soft condensed matter physics. It has long been suspected that this transition is driven by an onset of particle localization and associated emergence of a glassy modulus. However, progress towards an accepted understanding of glass formation has been impeded by an inability to obtain data sufficient in chemical diversity, relaxation timescales, and spatial and temporal resolution to validate or falsify proposed theories for its physics. Here we first describe a strategy enabling facile high-throughput simulation of glass-forming liquids to nearly unprecedented relaxation times. We then perform simulations of 51 glass-forming liquids, spanning polymers, small organic molecules, inorganics, and metallic glass-formers, with longest relaxation times exceeding one microsecond. Results identify a universal particle-localization transition accompanying glass formation across all classes of glass-forming liquid. The onset temperature of non-Arrhenius dynamics is found to serve as a normalizing condition leading to a master collapse of localization data. This transition exhibits a non-universal relationship with dynamic arrest, suggesting that the nonuniversality of supercooled liquid dynamics enters via the dependence of relaxation times on local cage scale. These results suggest that a universal particle-localization transition may underpin the glass transition, and they emphasize the potential for recent theoretical developments connecting relaxation to localization and emergent elasticity to finally explain the origin of this phenomenon. More broadly, the capacity for high-throughput prediction of glass formation behavior may open the door to computational inverse design of glass-forming materials.

Relaxation processes and glass transition of confined polymer melts: A molecular dynamics simulation of 1,4-polybutadiene between graphite walls

Molecular dynamics simulations of a chemically realistic model for 1,4-polybutadiene in a thin film geometry confined by two graphite walls are presented. Previous work on melts in the bulk has shown that the model faithfully reproduces static and dynamic properties of the real material over a wide temperature range. The present work studies how these properties change due to nano-confinement. The focus is on orientational correlations observable in nuclear magnetic resonance experiments and on the local intermediate incoherent neutron scattering function, F_s(q_z, z, t), for distances z from the graphite walls in the range of a few nanometers. Temperatures from about 2T_g down to about 1.15T_g, where T_g is the glass transition temperature in the bulk, are studied. It is shown that weakly attractive forces between the wall atoms and the monomers suffice to effectively bind a polymer coil that is near the wall. For a wide regime of temperatures, the Arrhenius-like adsorption/desorption kinetics of the monomers is the slowest process, while very close to T_g the Vogel-Fulcher-Tammann-like α-relaxation takes over. The α-process is modified only for z≤1.2 nm due to the density changes near the walls, less than expected from studies of coarse-grained (bead-spring-type) models. The weakness of the surface effects on the glass transition in this case is attributed to the interplay of density changes near the wall with the torsional potential. A brief discussion of pertinent experiments is given.

Chain relaxation in thin polymer films: turning a dielectric type-B polymer into a type-A' one

A molecular dynamics simulation study of chain relaxation in a thin polymer film is presented, studying the dielectric response of a random copolymer of cis and trans 1,4-polybutadiene, a type B polymer without net chain dipole moment, confined between graphite walls. We stress the orientational effect of the attractive walls, inducing polarization in the vicinity of the walls, while the center of the film stays bulk-like. This polarization leads to a net dipole moment of the adsorbed chains, which is perpendicular to their end-to-end vector, which we termed as type A' behavior. In this situation, the dipole moment relaxes only upon desorption of the chains from the wall, a dynamic process which occurs on timescales much longer than the bulk relaxation time of the polymer.

Dielectric α-relaxation of 1,4-polybutadiene confined between graphite walls : Molecular dynamics investigations through numerical simulations of polymer molecules relaxation

We present results of Molecular Dynamics (MD) simulations of a chemically realistic model of 1,4-polybutadiene confined by crystalline graphite walls. The simulations cover a large range of temperatures from T ≈ 2T g to T ≈ 1.15T g, where relevant time scales are accessible using such computational methods. We investigate the dielectric relaxation close to the walls in comparison to the one in the center of the film, and study the latter as a function of the film thickness from the walls. The segmental dynamics in the film is slowed down close to the walls, in comparison to the bulk. In addition to the α-process, the relaxation exhibits an additional long time decay, the so-called wall desorption process. We focus here on the α-process and find no significant shift of the dielectric T g as a function of layer thickness, in agreement with recent dielectric experiments. These findings can be correlated with the importance of the dihedral dynamics for all relaxation processes in polymers, which is unaltered except for the first nanometer next to the walls.

On the analysis of conformational dynamics in polymers with several rotational isomers

The ability of different correlation functions to shed some light onto the conformational dynamics of an amorphous polymer has been analyzed. The study has been performed on a polyethylene model polymer, which has been simulated at decreasing temperatures towards its glass transition, via the molecular dynamics technique. Three rotational isomers are allowed by the considered torsional potential. The correlation times associated with the evaluated transition rates have shown to be Arrhenius in nature, with activation energies resulting basically from internal rotation barriers. Overall torsional autocorrelation functions have been calculated. We have observed that they are dominated by slow events. Alternatively, a set of torsional autocorrelation functions associated with every isomeric state has been evaluated. Stretched exponential fits lead to correlation times that display Vogel-Fulcher temperature dependence.

Atomistic molecular dynamics simulation of the temperature and pressure dependences of local and terminal relaxations in cis-1,4-polybutadiene

The dynamics of cis-1,4-polybutadiene (cis-1,4-PB) over a wide range of temperature and pressure conditions is explored by conducting atomistic molecular dynamics (MD) simulations with a united atom model on a 32-chain C128 cis-1,4-PB system. The local or segmental dynamics is analyzed in terms of the dipole moment time autocorrelation function (DACF) of the simulated polymer and its temperature and pressure variations, for temperatures as low as 195 K and pressures as high as 3 kbars. By Fourier transforming the DACF, the dielectric spectrum, epsilon* = epsilon' + i epsilon" = epsilon*omega, is computed and the normalized epsilon"/epsilon(max)" vs omega/omega(max) plot is analyzed on the basis of the time-temperature and time-pressure superposition principles. The relative contribution of thermal energy and volume to the segmental and chain relaxation processes are also calculated and evaluated in terms of the ratio of the activation energy at constant volume to the activation energy at constant pressure, Q(V)/Q(P). Additional results for the temperature and pressure dependences of the Rouse times describing terminal relaxation in the two polymers show that, in the regime of the temperature and pressure conditions covered here, segmental and chain relaxations are influenced similarly by the pressure and temperature variations. This is in contrast to what is measured experimentally [see, e.g., G. Floudas and T. Reisinger, J. Chem. Phys. 111, 5201 (1999); C. M. Roland et al.,J. Polym. Sci. Part B, 41, 3047 (2003)] for other, chemically more complex polymers that pressure has a stronger influence on the dynamics of segmental mode than on the dynamics of the longest normal mode, at least for the regime of temperature and pressure conditions covered in the present MD simulations.

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.

Structural relaxation and dynamic heterogeneity in a polymer melt at attractive surfaces

Molecular dynamics simulations of polymer melts at flat and structured surfaces reveal that, for the former, slow dynamics and increased dynamic heterogeneity for an adsorbed polymer is due to densification of the polymer in a surface layer, while, for the latter, the energy topography of the surface plays the dominant role in determining dynamics of interfacial polymer. The dramatic increase in structural relaxation time for polymer melts at the attractive structured surface is largely the result of dynamic heterogeneity induced by the surface and does not resemble dynamics of a bulk melt approaching T(g).