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.     

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.     

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

Our theory of dynamic barriers, slow relaxation, and the glass transition of polymers melts is numerically applied using parameters relevant to real materials. The numerical results are found to be in qualitative agreement with all the approximate analytic expressions previously derived with quantitative differences on the order of approximately 20-30% or much less. The analytic prediction of a universal temperature dependence of the alpha relaxation time, and its intimate connection with the idea of a nearly universal crossover time, is established. Inter-relations between the breadth of the deeply supercooled regime, two definitions of the dynamic fragility, and the magnitude of the fast local Arrhenius process at the glass transition temperature are demonstrated and system-specific limitations identified. A quantitative application to segmental relaxation over 16 orders of magnitude in a polyvinylacetate melt yields encouraging results regarding the accuracy of the theory. The theoretical relaxation time results are well fit by multiple empirical forms (generally containing an assumed singular aspect) using parameters consistent with experimental studies. No physical significance is ascribed to this finding, but it does provide additional support for the temperature dependence of the alpha relaxation process predicted by the theory.     

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.     

Theory of nanoparticle diffusion in unentangled and entangled polymer melts

We propose a statistical dynamical theory for the violation of the hydrodynamic Stokes-Einstein (SE) diffusion law for a spherical nanoparticle in entangled and unentangled polymer melts based on a combination of mode coupling, Brownian motion, and polymer physics ideas. The non-hydrodynamic friction coefficient is related to microscopic equilibrium structure and the length-scale-dependent polymer melt collective density fluctuation relaxation time. When local packing correlations are neglected, analytic scaling laws (with numerical prefactors) in various regimes are derived for the non-hydrodynamic diffusivity as a function of particle size, polymer radius-of-gyration, tube diameter, degree of entanglement, melt density, and temperature. Entanglement effects are the origin of large SE violations (orders of magnitude mobility enhancement) which smoothly increase as the ratio of particle radius to tube diameter decreases. Various crossover conditions for the recovery of the SE law are derived, which are qualitatively distinct for unentangled and entangled melts. The dynamical influence of packing correlations due to both repulsive and interfacial attractive forces is investigated. A central finding is that melt packing fraction, temperature, and interfacial attraction strength all influence the SE violation in qualitatively different directions depending on whether the polymers are entangled or not. Entangled systems exhibit seemingly anomalous trends as a function of these variables as a consequence of the non-diffusive nature of collective density fluctuation relaxation and the different response of polymer-particle structural correlations to adsorption on the mesoscopic entanglement length scale. The theory is in surprisingly good agreement with recent melt experiments, and new parametric studies are suggested.     

Fragility of glass-forming polymer liquids

The fragility of polymeric glass-forming liquids is calculated as a function of molecular structural parameters from a generalized entropy theory of polymer glass-formation that combines the Adam-Gibbs (AG) model for the rate of structural relaxation with the lattice cluster theory (LCT) for polymer melt thermodynamics. Our generalized entropy theory predicts the existence of distinct high and low temperature regimes of glass-formation that are separated by a thermodynamically well-defined crossover temperature T(I) at which the product of the configurational entropy and the temperature has an inflection point. Since the predicted temperature dependence of the configurational entropy and structural relaxation time are quite different in these temperature regimes, we introduce separate definitions of fragility for each regime. Experimentally established trends in the fragility of polymer melts with respect to variations in polymer microstructure and pressure are interpreted within our theory in terms of the accompanying changes in the chain packing efficiency.     

Moderately and strongly supercooled liquids: a temperature-derivative study of the primary relaxation time scale

The primary relaxation time scale tau(T) derived from the glass forming supercooled liquids (SCLs) is discussed within ergodic-cluster Gaussian statistics, theoretically justified near and above the glass-transformation temperature T(g). An analysis is given for the temperature-derivative data by Stickel et al. on the steepness and the curvature of tau(T). Near the mode-coupling-theory (MCT) crossover T(c), these derivatives separate by a kink and a jump, respectively, the moderately and strongly SCL states. After accounting for the kink and the jump, the steepness remains a piecewise conitnuous function, a material-independent equation for the three fundamental characteristic temperatures, T(g), T(c), and the Vogel-Fulcher-Tamman (VFT) T(0), is found. Both states are described within the heterostructured model of solidlike clusters parametrized in a self-consistent manner by a minimum set of observable parameters: the fragility index, the MCT slowing-down exponent, and the chemical excess potential of Adam and Gibbs model (AGM). Below the Arrhenius temperature, the dynamically and thermodynamically stabilized clusters emerge with a size of around of seven to nine and two to three molecules above and close to T(g) and T(c), respectively. On cooling, the main transformation of the moderately into the strongly supercooled state is due to rebuilding of the cluster structure, and is attributed to its rigidity, introduced through the cluster compressibility. It is shown that the validity of the dynamic AGM (dynamically equivalent to the standard VFT form) is limited by the strongly supercooled state (T(g) < T < T(c)) where the superrigid cooperative rearranging regions are shown to be well-chosen parametrized solidlike clusters. Extension of the basic parameter set by the observable kinetic and diffusive exponents results in prediction of a subdiffusion relaxation regime in SCLs that is distinct from that established for amorphous polymers.     

Slow dynamics of the high density Gaussian core model

We numerically study crystal nucleation and glassy slow dynamics of the one-component Gaussian core model (GCM) at high densities. The nucleation rate at a fixed supercooling is found to decrease as the density increases. At very high densities, the nucleation is not observed at all in the time window accessed by long molecular dynamics (MD) simulation. Concomitantly, the system exhibits typical slow dynamics of the supercooled fluids near the glass transition point. We compare the simulation results of the supercooled GCM with the predictions of mode-coupling theory (MCT) and find that the agreement between them is better than any other model glassformers studied numerically in the past. Furthermore, we find that a violation of the Stokes-Einstein relation is weaker and the non-Gaussian parameter is smaller than canonical glassformers. Analysis of the probability distribution of the particle displacement clearly reveals that the hopping effect is strongly suppressed in the high density GCM. We conclude from these observations that the GCM is more amenable to the mean-field picture of the glass transition than other models. This is attributed to the long-ranged nature of the interaction potential of the GCM in the high density regime. Finally, the intermediate scattering function at small wavevectors is found to decay much faster than its self part, indicating that dynamics of the large-scale density fluctuations decouples with the shorter-ranged caging motion.     

Nuclear magnetic resonance investigation of dynamics in poly(ethylene oxide)-based lithium polyether-ester-sulfonate ionomers

Nuclear magnetic resonance spectroscopy has been utilized to investigate the dynamics of poly(ethylene oxide)-based lithium sulfonate ionomer samples that have low glass transition temperatures. (1)H and (7)Li spin-lattice relaxation times (T(1)) of the bulk polymer and lithium ions, respectively, were measured and analyzed in samples with a range of ion contents. The temperature dependence of T(1) values along with the presence of minima in T(1) as a function of temperature enabled correlation times and activation energies to be obtained for both the segmental motion of the polymer backbone and the hopping motion of lithium cations. Similar activation energies for motion of both the polymer and lithium ions in the samples with lower ion content indicate that the polymer segmental motion and lithium ion hopping motion are correlated in these samples, even though lithium hopping is about ten times slower than the segmental motion. A divergent trend is observed for correlation times and activation energies of the highest ion content sample with 100% lithium sulfonation due to the presence of ionic aggregation. Details of the polymer and cation dynamics on the nanosecond timescale are discussed and complement the findings of X-ray scattering and quasi-elastic neutron scattering experiments.     

Positron annihilation and relaxation dynamics from dielectric spectroscopy and nuclear magnetic resonance: cis-trans-1,4-poly(butadiene)

We report a joint analysis of positron annihilation lifetime spectroscopy (PALS), dielectric spectroscopy (BDS), and nuclear magnetic resonance (NMR) on cis-trans-1,4-poly(butadiene) (c-t-1,4-PBD). Phenomenological analysis of the orthopositronium lifetime τ(3)-T dependence by linear fitting reveals four characteristic PALS temperatures: T(b1)(G)=0.63T(g)(PALS), T(g)(PALS), T(b1)(L)=1.22T(g)(PALS), and T(b2)(L)=1.52T(g)(PALS). Slight bend effects in the glassy and supercooled liquid states are related to the fast or slow secondary β process, from neutron scattering, respectively, the latter being connected with the trans-isomers. In addition, the first bend effect in the supercooled liquid coincides with a deviation of the slow effective secondary β(eff) relaxation related to the cis-isomers from low-T Arrhenius behavior to non-Arrhenius one and correlates with the onset of the primary α process from BDS. The second plateau effect in the liquid state occurs when τ(3) becomes commensurable with the structural relaxation time τ(α)(T(b2)). It is also approximately related to its crossover from non-Arrhenius to Arrhenius regime in the combined BDS and NMR data. Finally, the combined BDS and NMR structural relaxation data, when analyzed in terms of the two-order parameter (TOP) model, suggest the influence of solidlike domains on both the annihilation behavior and the local and segmental chain mobility in the supercooled liquid. All these findings indicate the influence of the dynamic heterogeneity in both the primary and secondary relaxations due to the cis-trans isomerism in c-t-1,4-PBD and their impact into the PALS response.     

Fluorescence lifetime fluctuations of single molecules probe the local environment of oligomers around the glass transition temperature

Single molecule fluorescence experiments have been performed on a BODIPY-based dye embedded in oligo(styrene) matrices to probe the density fluctuations and the relaxation dynamics of chain segments surrounding the dye molecules. The time-dependent fluorescence lifetime of the BODIPY probe was recorded as an observable for the local density fluctuations. At room temperature, the mean fraction of holes surrounding the probes is shown to be unaffected by the molecular weight in the glassy state. In contrast, the free volume increases significantly in the supercooled regime. These observations are discussed in the framework of the entropic theories of the glass transition.     

Decoupling charge transport from the structural dynamics in room temperature ionic liquids

Light scattering and dielectric spectroscopy measurements were performed on the room temperature ionic liquid (RTIL) [C4mim][NTf2] in a broad temperature and frequency range. Ionic conductivity was used to estimate self-diffusion of ions, while light scattering was used to study structural relaxation. We demonstrate that the ionic diffusion decouples from the structural relaxation process as the temperature of the sample decreases toward T(g). The strength of the decoupling appears to be significantly lower than that expected for a supercooled liquid of similar fragility. The structural relaxation process in the RTIL follows well the high-temperature mode coupling theory (MCT) scenario. Using the MCT analysis we estimated the dynamic crossover temperature in [C4mim][NTf2] to be T(c) ~ 225 ± 5 K. However, our analysis reveals no sign of the dynamic crossover in the ionic diffusion process.     

Time dependent stretching of aging dynamics in a generalized hydrodynamic model for supercooled liquids

The nonequilibrium dynamics and aging behavior of a supercooled liquid is investigated from an analysis of the correlation of density fluctuations at two different times. The dynamic correlation functions are computed by solving numerically the equations of nonlinear fluctuating hydrodynamics. The aging time dependence follows a modified stretched exponential form with a relaxation time which is dependent on the aging time. This is similar to the behavior seen in the aging data of dielectric response functions of a typical glass forming liquid.     

Glassy dynamics and kinetic vitrification of isotropic suspensions of hard rods

A nonlinear Langevin equation (NLE) theory for the translational center-of-mass dynamics of hard nonspherical objects has been applied to isotropic fluids of rigid rods. The ideal kinetic glass transition volume fraction is predicted to be a monotonically decreasing function beyond an aspect ratio of two. The functional form of the decrease is weaker than the inverse aspect ratio. Vitrification occurs at lower volume fractions for corrugated tangent bead rods compared to their smooth spherocylinder analogs. The ideal glass transition signals a crossover to activated dynamics, which is estimated to be observable before the nematic phase boundary is encountered if the aspect ratio is less than roughly 25. Calculations of the glassy elastic shear modulus and absolute yield stress reveal a roughly exponential growth with volume fraction. The dependence of entropic barriers and mean barrier hopping times on concentration for rods of variable aspect ratios can be collapsed quite well based on a difference volume fraction variable that quantifies the distance from the ideal glass boundary. Full numerical solution of the NLE theory via stochastic trajectory simulation was performed for tangent bead rods, and the results were compared to their hard sphere analogs. With increasing shape anisotropy the characteristic length scales of the nonequilibrium free energy increase and the magnitude of the localization well and entropic barrier curvatures decreases. These changes result in a significant aspect ratio dependence of dynamical properties and time correlation functions including weaker intermediate time subdiffusive transport, stronger two-step decay of the incoherent dynamic structure factor, longer mean alpha relaxation time, and stronger wavevector-dependent decoupling of relaxation times and the self-diffusion constant. The theoretical results are potentially testable via computer simulation, confocal microscopy, and dynamic light scattering.     

Dynamic yielding, shear thinning, and stress rheology of polymer-particle suspensions and gels

The nonlinear rheological version of our barrier hopping theory for particle-polymer suspensions and gels has been employed to study the effect of steady shear and constant stress on the alpha relaxation time, yielding process, viscosity, and non-Newtonian flow curves. The role of particle volume fraction, polymer-particle size asymmetry ratio, and polymer concentration have been systematically explored. The dynamic yield stress decreases in a polymer-concentration- and volume-fraction-dependent manner that can be described as apparent power laws with effective exponents that monotonically increase with observation time. Stress- or shear-induced thinning of the viscosity becomes more abrupt with increasing magnitude of the quiescent viscosity. Flow curves show an intermediate shear rate dependence of an effective power-law form, becoming more solidlike with increasing depletion attraction. The influence of polymer concentration, particle volume fraction, and polymer-particle size asymmetry ratio on all properties is controlled to a first approximation by how far the system is from the gelation boundary of ideal mode-coupling theory (MCT). This emphasizes the importance of the MCT nonergodicity transition despite its ultimate destruction by activated barrier hopping processes. Comparison of the theoretical results with limited experimental studies is encouraging.     

Examination of dynamic facilitation in molecular dynamics simulations of glass-forming liquids

Using data from molecular dynamics computer simulations of the one-component Dzugutov liquid and of BKS silica in metastable equilibrium supercooled states, we examine ideas introduced by Garrahan and Chandler (GC) in their dynamic facilitation (DF) model of the glass transition. Utilizing a recently introduced measure of DF, we find that DF is important for particle motion in both the supercooled Dzugutov liquid and in the BKS silica melt, that mobility propagates continuously, and that this effect becomes increasingly pronounced with decreasing T. We show that, in both systems, dynamic facilitation is strongest on the time scale of the late-beta relaxation, where clusters of highly mobile neighboring particles escaping from their cages are largest and, except for the silicon atoms in BKS silica, stringlike motion is most prominent. By comparing the two systems, we show that the temperature dependence of one measure of DF as the mode-coupling temperature is approached from high temperature is similar, once the temperature dependence of the structural relaxation time in each system is scaled out.     

Superdipole liquid scenario for the dielectric primary relaxation in supercooled polar liquids

We propose a dynamic structure of coupled dynamic molecular strings for supercooled small polar molecule liquids and accordingly we obtain the Hamiltonian of the rotational degrees of freedom of the system. From the Hamiltonian, the strongly correlated supercooled polar liquid state is renormalized to a normal superdipole liquid state. This scenario describes the following main features of the primary or alpha-relaxation dynamics in supercooled polar liquids: (1) the average relaxation time evolves from a high temperature Arrhenius to a low temperature non-Arrhenius or super-Arrhenius behavior; (2) the relaxation function crosses over from the high temperature exponential to low temperature nonexponential form; and (3) the temperature dependence of the relaxation strength shows non-Curie features. According to the present model, the crossover phenomena of the first two characteristics arise from the transition between the superdipole gas and the superdipole liquid. The model predictions are quantitatively compared with the experimental results of glycerol, a typical glass former.     

Intermolecular forces and the glass transition

Random first-order transition theory is used to determine the role of attractive and repulsive interactions in the dynamics of supercooled liquids. Self-consistent phonon theory, an approximate mean field treatment consistent with random first-order transition theory, is used to treat individual glassy configurations, whereas the liquid phase is treated using common liquid-state approximations. Free energies are calculated using liquid-state perturbation theory. The transition temperature, T*A, the temperature where the onset of activated behavior is predicted by mean field theory; the lower crossover temperature, T*C, where barrierless motions actually occur through fractal or stringy motions (corresponding to the phenomenological mode coupling transition temperature); and T*K, the Kauzmann temperature (corresponding to an extrapolated entropy crisis), are calculated in addition to T*g, the glass transition temperature that corresponds to laboratory cooling rates. Relationships between these quantities agree well with existing experimental and simulation data on van der Waals liquids. Both the isobaric and isochoric behavior in the supercooled regime are studied, providing results for DeltaCV and DeltaCp that can be used to calculate the fragility as a function of density and pressure, respectively. The predicted variations in the alpha-relaxation time with temperature and density conform to the empirical density-temperature scaling relations found by Casalini and Roland. We thereby demonstrate the microscopic origin of their observations. Finally, the relationship first suggested by Sastry between the spinodal temperature and the Kauzmann temperatures, as a function of density, is examined. The present microscopic calculations support the existence of an intersection of these two temperatures at sufficiently low temperatures.     

Monte Carlo simulation studies of the interfaces between polymeric and other solids as models for fiber-matrix interactions in advanced composite materials

As a coarse-grained model for dense amorphous polymer systems interacting with solid walls (i.e., the fiber surface in a composite), the bond fluctuation model of flexible polymer chains confined between two repulsive surfaces is studied by extensive Monte Carlo simulations. Choosing a potential for the length of an effective bond that favors rather long bonds, the full temperature region from ordinary polymer melts down to the glass transition is accessible. It is shown that in the supercooled state near the glass transition an “interphase” forms near the walls, where the structure of the melt is influenced by the surface. This “interphase” already shows up in static properties, but also has an effect on monomer mobilities and the corresponding relaxation behavior of the polymer matrix. The thickness of the interphase is extracted from monomer density oscillations near the walls and is found to be strongly temperature dependent. It is ultimately larger than the gyration radius of the polymer chains. Effects of shear deformation on this model are simulated by choosing asymmetric jump rates near the moving wall (large jump rate in the direction of motion, and a small rate against it). It is studied how this dynamic perturbation propagates into the bulk of the polymer matrix.     

In search of temporal power laws in the orientational relaxation near isotropic-nematic phase transition in model nematogens

Recent Kerr relaxation experiments by Gottke et al. have revealed the existence of a pronounced temporal power law decay in the orientational relaxation near the isotropic-nematic phase transition (INPT) of nematogens of rather small aspect ratio, kappa (kappa approximately 3-4). We have carried out very long (50 ns) molecular dynamics simulations of model (Gay-Berne) prolate ellipsoids with aspect ratio 3 in order to investigate the origin of this power law. The model chosen is known to undergo an isotropic to nematic phase transition for a range of density and temperature. The distance dependence of the calculated angular pair correlation function correctly shows the emergence of a long range correlation as the INPT is approached along the density axis. In the vicinity of INPT, the single particle second rank orientational time correlation function exhibits power law decay, (t(-alpha)) with exponent alpha approximately 2/3. More importantly, we find the sudden appearance of a pronounced power-law decay in the collective part of the second rank orientational time correlation function at short times when the density is very close to the transition density. The power law has an exponent close to unity, that is, the correlation function decays almost linearly with time. At long times, the decay is exponential-like, as predicted by Landau-de Gennes mean field theory. Since Kerr relaxation experiments measure the time derivative of the collective second rank orientational pair correlation function, the simulations recover the near independence of the signal on time observed in experiments. In order to capture the microscopic essence of the dynamics of pseudonematic domains inside the isotropic phase, we introduce and calculate a dynamic orientational pair correlation function (DOPCF) obtained from the coefficients in the expansion of the distinct part of orientational van Hove time correlation function in terms of spherical harmonics. The DOPCF exhibits power law relaxation when the pair separation length is below certain critical length. The orientational relaxation of a local director, defined in terms of the sum of unit vectors of all the ellipsoidal molecules, is also found to show slow power law relaxation over a long time scale. These results have been interpreted in terms of a newly developed mode coupling theory of orientational dynamics near the INPT. In the present case, the difference between the single particle and the collective orientational relaxation is huge which can be explained by the frequency dependence of the memory kernel, calculated from the mode coupling theory. The relationship of this power law with the one observed in a supercooled liquid near its glass transition temperature is explored.