Universal scaling, dynamic fragility, segmental relaxation, and vitrification in polymer melts

Temperature Dependence of Structural Relaxation in Glass-Forming Liquids and Polymers

Understanding the microscopic mechanism of the transition of glass remains one of the most challenging topics in Condensed Matter Physics. What controls the sharp slowing down of molecular motion upon approaching the glass transition temperature Tg, whether there is an underlying thermodynamic transition at some finite temperature below Tg, what the role of cooperativity and heterogeneity are, and many other questions continue to be topics of active discussions. This review focuses on the mechanisms that control the steepness of the temperature dependence of structural relaxation (fragility) in glass-forming liquids. We present a brief overview of the basic theoretical models and their experimental tests, analyzing their predictions for fragility and emphasizing the successes and failures of the models. Special attention is focused on the connection of fast dynamics on picosecond time scales to the behavior of structural relaxation on much longer time scales. A separate section discusses the specific case of polymeric glass-forming liquids, which usually have extremely high fragility. We emphasize the apparent difference between the glass transitions in polymers and small molecules. We also discuss the possible role of quantum effects in the glass transition of light molecules and highlight the recent discovery of the unusually low fragility of water. At the end, we formulate the major challenges and questions remaining in this field.

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.

From boiling point to glass transition temperature: transport coefficients in molecular liquids follow three-parameter scaling

The phenomenon of the glass transition is an unresolved problem in condensed matter physics. Its prominent feature, the super-Arrhenius temperature dependence of the transport coefficients, remains a challenge to be described over the full temperature range. For a series of molecular glass formers, we combined τ(T) collected from dielectric spectroscopy and dynamic light scattering covering a range 10(-12) s < τ(T) < 10(2) s. Describing the dynamics in terms of an activation energy E(T), we distinguish a high-temperature regime characterized by an Arrhenius law with a constant activation energy E(∞) and a low-temperature regime for which E(coop)(T) ≡ E(T)-E(∞) increases exponentially while cooling. A scaling is introduced, specifically E(coop)(T)/E(∞) [proportionality] exp[-λ(T/T(A)-1)], where λ is a fragility parameter and T(A) a reference temperature proportional to E(∞). In order to describe τ(T) still the attempt time τ(∞) has to be specified. Thus, a single interaction parameter E(∞) describing the high-temperature regime together with λ controls the temperature dependence of low-temperature cooperative dynamics.

Theory of aging, rejuvenation, and the nonequilibrium steady state in deformed polymer glasses

The nonlinear Langevin equation theory of segmental relaxation, elasticity, and mechanical response of polymer glasses is extended to describe the coupled effects of physical aging, mechanical rejuvenation, and thermal history. The key structural variable is the amplitude of density fluctuations, and segmental dynamics proceeds via stress-modified activated barrier hopping on a dynamic free-energy profile. Mechanically generated disorder (rejuvenation) is quantified by a dissipative work argument and increases the amplitude of density fluctuations, thereby speeding up relaxation beyond that induced by the landscape tilting mechanism. The theory makes testable predictions for the time evolution and nonequilibrium steady state of the alpha relaxation time, density fluctuation amplitude, elastic modulus, and other properties. Model calculations reveal a rich dependence of these quantities on preaging time, applied stress, and temperature that reflects the highly nonlinear competition between physical aging and mechanical disordering. Thermal history is "erased" in the long-time limit, although the nonequilibrium steady state is not the literal "fully rejuvenated" freshly quenched glass. The present work provides the conceptual foundation for a quantitative treatment of the nonlinear mechanical response of polymer glasses under a variety of deformation protocols.

Segmental dynamics in polymers: from cold melts to ageing and stressed glasses

Recent progress in developing statistical mechanical theories of supercooled polymer melts and glasses is reviewed. The focus is on those approaches that are either explicitly formulated for polymers, or are applications of more generic theories to interpret polymeric phenomena. These include two configurational entropy theories, a percolated free volume distribution model, and the activated barrier hopping nonlinear Langevin theory. Both chemically-specific and universal aspects are discussed. After a brief summary of classic phenomenological approaches, a discussion of the relevant length scales and key experimental phenomena in both the supercooled liquid and glassy solid state is presented including ageing and nonlinear mechanical response. The central concepts that underlie the theories in the molten state are then summarized and key predictions discussed, including the glass transition in oriented polymer liquids and deformed rubber networks. Physical ageing occurs in the nonequilibrium glass, and theories for its consequences on the alpha relaxation are discussed. Very recent progress in developing a segment scale theory for the dramatic effects of external stress on polymer glasses, including acceleration of relaxation, yielding, plastic flow and strain hardening, is summarized. The article concludes with a discussion of outstanding theoretical challenges.

Theory of physical aging in polymer glasses

A statistical segment scale theory for the physical aging of polymer glasses is proposed and applied. The approach is based on a nonlinear stochastic Langevin equation of motion and the concept of an effective free energy which quantifies temporary localization, collective barriers, and the activated segment hopping process. The key collective structural variable that plays the role of the dynamic order parameter for aging is the experimentally measurable nanometer and longer wavelength amplitude of density fluctuations, S0 . The degree of local cooperativity, and the bare activation energy of the high-temperature Arrhenius process, are determined in the molten state by utilizing experimental measurements of the glass temperature and dynamic crossover time, respectively. A first-order kinetic equation with a time varying rate is proposed for the temporal evolution of S0 which is self-consistently and nonlinearly coupled with the mean segmental relaxation time. The theory has been applied to study physical aging of the alpha relaxation time, shear relaxation modulus, amplitude of density fluctuations, cohesive energy, absolute yield stress, and fictive temperature of polymethylmethacrylate and other glasses over a range of temperatures. Temperature-dependent logarithmic and effective power-law aging is predicted at intermediate times. Time-aging time superposition is found for the mechanical relaxation function. A strongly asymmetric aging response is predicted for up and down temperature jump experiments. Comparison of the approach with the classic phenomenological model is presented.

Theory of relaxation and elasticity in polymer glasses

The recently developed activated barrier hopping theory of deeply supercooled polymer melts [K. S. Schweizer and E. J. Saltzman, J. Chem. Phys. 121, 1984 (2004)] is extended to the nonequilibrium glass state. Below the kinetic glass temperature T(g), the exact statistical mechanical relation between the dimensionless amplitude of long wavelength density fluctuations, S(0), and the thermodynamic compressibility breaks down. Proper extension of the theory requires knowledge of the nonequilibrium S(0) which x-ray scattering experiments find to consist of a material specific and temperature-independent quenched disorder contribution plus a vibrational contribution which varies roughly linearly with temperature. Motivated by these experiments and general landscape concepts, a simple model is proposed for S(0)(T). Deep in the glass state the form of the temperature dependence of the segmental relaxation time is found to depend sensitively on the magnitude of frozen in density fluctuations. At the (modest) sub-T(g) temperatures typically probed in experiment, an effective Arrhenius behavior is generically predicted which is of nonequilibrium origin. The change in apparent activation energy across the glass transition is determined by the amplitude of frozen density fluctuations. For values of the latter consistent with experiment, the theory predicts a ratio of effective activation energies in the range of 3-6, in agreement with multiple measurements. Calculations of the shear modulus for atactic polymethylmethacrylate above and below the glass transition temperature have also been performed. The present work provides a foundation for the formulation of predictive theories of physical aging, the influence of deformation on the alpha relaxation process, and rate-dependent nonlinear mechanical properties of thermoplastics.

Ideal glass transitions, shear modulus, activated dynamics, and yielding in fluids of nonspherical objects

An extension of naive ideal mode coupling theory (MCT) and its generalization to treat activated barrier hopping and glassy dynamics in fluids and suspensions composed of nonspherical hard core objects is proposed. An effective center-of-mass description is adopted. It corresponds to a specific type of pre-averaging of the dynamical consequences of orientational degrees of freedom. The simplest case of particles composed of symmetry-equivalent interaction sites is considered. The theory is implemented for a homonuclear diatomic shape of variable bond length. The naive MCT glass transition boundary is predicted to be a nonmonotonic function of the length-to-width or aspect ratio and occurs at a nearly unique value of the dimensionless compressibility. The latter quantifies the amplitude of long wavelength thermal density fluctuations, thereby (empirically) suggesting a tight connection between the onset of localization and thermodynamics. Localization lengths and elastic shear moduli for different aspect ratio and volume fraction systems approximately collapse onto master curves based on a reduced volume fraction variable that quantifies the distance from the ideal glass transition. Calculations of the entropic barrier height and hopping time, maximum restoring force, and absolute yield stress and strain as a function of diatomic aspect ratio and volume fraction have been performed. Strong correlations of these properties with the dimensionless compressibility are also found, and nearly universal dependences have been numerically identified based on property-specific nondimensionalizations. Generalization of the approach to rigid rods, disks, and variable shaped molecules is possible, including oriented liquid crystalline phases.

Non-Gaussian effects, space-time decoupling, and mobility bifurcation in glassy hard-sphere fluids and suspensions

Brownian trajectory simulation methods are employed to fully establish the non-Gaussian fluctuation effects predicted by our nonlinear Langevin equation theory of single particle activated dynamics in glassy hard-sphere fluids. The consequences of stochastic mobility fluctuations associated with the space-time complexities of the transient localization and barrier hopping processes have been determined. The incoherent dynamic structure factor was computed for a range of wave vectors and becomes of an increasingly non-Gaussian form for volume fractions beyond the (naive) ideal mode coupling theory (MCT) transition. The non-Gaussian parameter (NGP) amplitude increases markedly with volume fraction and is well described by a power law in the maximum restoring force of the nonequilibrium free energy profile. The time scale associated with the NGP peak becomes much smaller than the alpha relaxation time for systems characterized by significant entropic barriers. An alternate non-Gaussian parameter that probes the long time alpha relaxation process displays a different shape, peak intensity, and time scale of its maximum. However, a strong correspondence between the classic and alternate NGP amplitudes is predicted which suggests a deep connection between the early and final stages of cage escape. Strong space-time decoupling emerges at high volume fractions as indicated by a nondiffusive wave vector dependence of the relaxation time and growth of the translation-relaxation decoupling parameter. Displacement distributions exhibit non-Gaussian behavior at intermediate times, evolving into a strongly bimodal form with slow and fast subpopulations at high volume fractions. Qualitative and semiquantitative comparisons of the theoretical results with colloid experiments, ideal MCT, and multiple simulation studies are presented.

Activated hopping and dynamical fluctuation effects in hard sphere suspensions and fluids

Single particle Brownian dynamics simulation methods are employed to establish the full trajectory level predictions of our nonlinear stochastic Langevin equation theory of activated hopping dynamics in glassy hard sphere suspensions and fluids. The consequences of thermal noise driven mobility fluctuations associated with the barrier hopping process are determined for various ensemble-averaged properties and their distributions. The predicted mean square displacements show classic signatures of transient trapping and anomalous diffusion on intermediate time and length scales. A crossover to a stronger volume fraction dependence of the apparent nondiffusive exponent occurs when the entropic barrier is of order the thermal energy. The volume fraction dependences of various mean relaxation times and rates can be fitted by empirical critical power laws with parameters consistent with ideal mode-coupling theory. However, the results of our divergence-free theory are largely a consequence of activated dynamics. The experimentally measurable alpha relaxation time is found to be very similar to the theoretically defined mean reaction time for escape from the barrier-dominated regime. Various measures of decoupling have been studied. For fluid states with small or nonexistent barriers, relaxation times obey a simple log-normal distribution, while for high volume fractions the relaxation time distributions become Poissonian. The product of the self-diffusion constant and mean alpha relaxation time increases roughly as a logarithmic function of the alpha relaxation time. The cage scale incoherent dynamic structure factor exhibits nonexponential decay with a modest degree of stretching. A nearly universal collapse of the different volume fraction results occurs if time is scaled by the mean alpha relaxation time. Hence, time-volume fraction superposition holds quite well, despite the presence of stretching and volume fraction dependent decoupling associated with the stochastic barrier hopping process. The relevance of other origins of dynamic heterogeneity (e.g., mesoscopic domains), and comparison of our results with experiments, simulations, and alternative theories, is discussed.

Entropy theory of polymer glass formation revisited. I. General formulation

A generalized entropy theory of glass formation is developed by merging the lattice cluster theory for the thermodynamics of semiflexible polymer melts at constant pressure with the Adam-Gibbs relation between the structural relaxation time and the configurational entropy. Since experimental studies have suggested that the relative rigidity of the chain backbone and the side groups is an essential parameter governing the nature of glass formation in polymers, we incorporate this rigidity disparity parameter, along with monomer structure, into our new theoretical description of the polymer fluid thermodynamics. Our entropy theory is compared with alternative theories that describe the rate of structural relaxation in glass-forming liquids in terms of an activated rate process.

Derivation of a microscopic theory of barriers and activated hopping transport in glassy liquids and suspensions

A recently proposed microscopic activated barrier hopping theory [K. S. Schweizer and E. J. Saltzman, J. Chem. Phys. 119, 1181 (2003)] of slow single-particle dynamics in glassy liquids, suspensions, and gels is derived using nonequilibrium statistical mechanics. Fundamental elements underlying the stochastic nonlinear Langevin equation description include an inhomogeneous liquid or locally solid-state perspective, dynamic density-functional theory (DDFT), a local equilibrium closure, and a coarse-grained free-energy functional. A dynamic Gaussian approximation is not adopted which is the key for avoiding a kinetic ideal glass transition. The relevant excess free energy is of a nonequilibrium origin and is related to dynamic force correlations in the fluid. The simplicity of the approach allows external perturbations to be rather easily incorporated. Dynamic heterogeneity enters naturally via mobility fluctuations associated with the stochastic barrier-hopping process. The derivation both identifies the limitations of the theory and suggests new avenues for its systematic improvement. Comparisons with ideal mode-coupling theory, alternative DDFT approaches and a field theoretic path-integral formulation are presented.

Theory of glassy dynamics in conformationally anisotropic polymer systems

A mode coupling theory for the ideal glass transition temperature, or crossover temperature to highly activated dynamics in the deeply supercooled regime, T(c), has been developed for anisotropic polymer liquids. A generalization of a simplified mode coupling approach at the coarse-grained segment level is employed which utilizes structural and thermodynamic information from the anisotropic polymer reference interaction site model theory. Conformational alignment or/and coil deformation modifies equilibrium properties and constraining interchain forces thereby inducing anisotropic segmental dynamics. For liquid-crystalline polymers a small suppression of T(c) with increasing nematic or discotic orientational order is predicted. The underlying mechanism is reduction of the degree of coil interpenetration and intermolecular repulsive contacts due to segmental alignment. For rubber networks chain deformation results in an enhanced bulk modulus and a modest elevation of T(c) is predicted. The theory can also be qualitatively applied to systems that undergo nonuniversal local deformation and alignment, such as polymer thin films and grafted brush layers, and large elevations or depressions of T(c) are possible. Extension to treat directionally dependent collective barrier formation and activated hopping is possible.