Dynamic fragility in polymers: A comparison in isobaric and isochoric conditions

Amorphous cis-1,4-polybutadiene P-V-T properties from atomistic simulations

CONTEXT

As a result of the diversity of microstructures encountered in cis-1,4-polybutadiene and the variety of measurement methods used, experimental values of variation of glass transition temperature (Tg) with pressure are relatively dispersed. However, atomistic simulations enable access to valuable information for very well-controlled chemistry and structures with a well-defined and systematic acquisition protocol. By varying the temperature and pressure, the specific volume of the melt was computed, yielding results that deviated by only 2% from experimental data. A linear relationship between Tg and pressure was observed, with Tg predicted to be 162 K at zero pressure and a rate of change of Tg with respect to pressure (dTg/dP) of 0.24 K/MPa.

METHOD

The atomistic dilatometry experiments were conducted on a model of amorphous cis-1,4 polybutadiene with an approximate molecular weight of 5400 g/mol using the LAMMPS code and the all-atom forcefield pcff + . The dilatometry process involved cooling and heating at a rate of 9 × 10¹² K/min. The specific volume was calculated by averaging over seven independent configurations for each temperature. The Tait equation was employed to fit the specific volume evolution within the temperature range of 10 to 700 K under different pressures of 0, 60, and 100 MPa.

Perspective: Searching for simplicity rather than universality in glass-forming liquids

This article gives an overview of experimental results on dynamics in bulk glass-forming molecular liquids. Rather than looking for phenomenology that is universal, in the sense that it is seen in all liquids, the focus is on identifying the basic characteristics, or "stylized facts," of the glass transition problem, i.e., the central observations that a theory of the physics of glass formation should aim to explain in a unified manner.

Linking slow dynamics and microscopic connectivity in dense suspensions of charged colloids

The quest to unravel the nature of the glass transition, where the viscosity of a liquid increases by many orders of magnitude, while its static structure remains largely unaffected, remains unresolved. While various structural and dynamical precursors to vitrification have been identified, a predictive and quantitative description of how subtle changes at the microscopic scale give rise to the steep growth in macroscopic viscosity is missing. It was recently proposed that the presence of long-lived bonded structures within the liquid may provide the long-sought connection between local structure and global dynamics. Here we directly observe and quantify the connectivity dynamics in liquids of charged colloids en route to vitrification using three-dimensional confocal microscopy. We determine the dynamic structure from the real-space van Hove correlation function and from the particle trajectories, providing upper and lower bounds on connectivity dynamics. Based on these data, we extend Dyre's model for the glass transition to account for particle-level structural dynamics; this results in a microscopic expression for the slowing down of relaxations in the liquid that is in quantitative agreement with our experiments. These results indicate how vitrification may be understood as a dynamical connectivity transition with features that are strongly reminiscent of rigidity percolation scenarios.

Correlation between Fragility and the Arrhenius Crossover Phenomenon in Metallic, Molecular, and Network Liquids

We report the observation of a distinct correlation between the kinetic fragility index m and the reduced Arrhenius crossover temperature θ_{A}=T_{A}/T_{g} in various glass-forming liquids, identifying three distinguishable groups. In particular, for 11 glass-forming metallic liquids, we universally observe a crossover in the mean diffusion coefficient from high-temperature Arrhenius to low-temperature super-Arrhenius behavior at approximately θ_{A}≈2 which is in the stable liquid phases. In contrast, for fragile molecular liquids, this crossover occurs at much lower θ_{A}≈1.4 and usually in their supercooled states. The θ_{A} values for strong network liquids spans a wide range higher than 2. Intriguingly, the high-temperature activation barrier E_{∞} is universally found to be ∼11k_{B}T_{g} and uncorrelated with the fragility or the reduced crossover temperature θ_{A} for metallic and molecular liquids. These observations provide a way to estimate the low-temperature glassy characteristics (T_{g} and m) from the high-temperature liquid quantities (E_{∞} and θ_{A}).

String model for the dynamics of glass-forming liquids

We test the applicability of a living polymerization theory to describe cooperative string-like particle rearrangement clusters (strings) observed in simulations of a coarse-grained polymer melt. The theory quantitatively describes the interrelation between the average string length L, configurational entropy Sconf, and the order parameter for string assembly Φ without free parameters. Combining this theory with the Adam-Gibbs model allows us to predict the relaxation time τ in a lower temperature T range than accessible by current simulations. In particular, the combined theories suggest a return to Arrhenius behavior near Tg and a low T residual entropy, thus avoiding a Kauzmann "entropy crisis."

Fragility of an Isochorically Confined Polymer Glass

We report the effect of isochoric confinement on the dynamic fragility of a polymeric glass-former, that is, polystyrene (PS). Utilizing silica-capped PS nanospheres as a model system, the fictive temperature (Tf) and the isochoric heat capacity (Cv) are measured as a function of diameter via differential scanning calorimetry. By examining Tf as a function of cooling rate for each sample, the isochoric fragility (mv) is obtained, which decreases significantly as the diameter of the nanospheres is reduced from 260 to 129 nm. Hence, the temperature dependence of structural relaxation near the glass transition is weakened with isochoric confinement.

The physics of the colloidal glass transition

As one increases the concentration of a colloidal suspension, the system exhibits a dramatic increase in viscosity. Beyond a certain concentration, the system is said to be a colloidal glass; structurally, the system resembles a liquid, yet motions within the suspension are slow enough that it can be considered essentially frozen. For several decades, colloids have served as a valuable model system for understanding the glass transition in molecular systems. The spatial and temporal scales involved allow these systems to be studied by a wide variety of experimental techniques. The focus of this review is the current state of understanding of the colloidal glass transition, with an emphasis on experimental observations. A brief introduction is given to important experimental techniques used to study the glass transition in colloids. We describe features of colloidal systems near and in glassy states, including increases in viscosity and relaxation times, dynamical heterogeneity and ageing, among others. We also compare and contrast the glass transition in colloids to that in molecular liquids. Other glassy systems are briefly discussed, as well as recently developed synthesis techniques that will keep these systems rich with interesting physics for years to come.

Glass transition and fragility in the simple molecular glassformer CS(2) from CS(2)-S(2)Cl(2) solution studies

With an interest in finding the fragility for a simple, single component, molecular glassformer, we have determined the dielectric relaxation and glass transition behavior for a series of glasses in the CS(2)-S(2)Cl(2) and CS(2)-toluene systems. Crystallization of CS(2) can be completely avoided down to the composition 20 mol% second component, and the fragility proves almost independent of CS(2) content in each system. Since the glass temperature T(g) obtained from both thermal studies and from dielectric relaxation (using T(g,diel)=T(tau=100 s)) is quite linear over the whole composition range in each system, and since relaxation time data for pure CS(2) fall on the same master plot when scaled by the linearly extrapolated T(g) value, we deduce that pure CS(2) has the same high fragility as the binary solutions. The value is m=86, as for ortho-terphenyl (OTP). Based on observations of independent studies for the vibrational density of states (VDoS) (of inherent structures for OTP and instantaneous, at-temperature structures for CS(2)), we attribute the high fragility to an excess vibrational heat capacity (defined by C(p) (vib, excess)=dS(vib, excess)/d ln T) originating in the behavior of the low frequency modes of the VDoS (the boson peak modes). Both low frequency DoS and anharmonicity increase with increasing temperature, augmenting the configurational entropy drive to the top of the system energy landscape. The surprising implication is that fragility is determined in the vibrational, not configurational, manifold of microstates.

High pressure dynamics of polymer/plasticizer mixtures

Plasticizers are usually added to polymers to give them the desired flexibility and processability by changing the dynamical properties of the polymer chains. It is therefore important to give a quantitative description about how the dynamic behavior of a given polymer is modified by the incorporation of a second component. We analyze in this work, by means of dielectric spectroscopy, the dynamics of poly(vinyl acetate)/diethyl phthalate mixtures, at different concentrations, over a broad range of frequency, pressure, and temperature. The dynamics of these particular mixtures show only one main relaxation process contrarily to what is observed in athermal miscible polymer mixtures. From the dielectric spectra the maximum relaxation time as a function of pressure and temperature was obtained and analyzed. We studied the pressure dependence of the glass transition temperature as well as the fragility of both the neat components and the mixtures at different concentrations (on the rich polymer range). Finally, the experimental data were rationalized within the framework of an Adam-Gibbs (AG) based approach recently developed [G. A. Schwartz et al., J. Chem. Phys. 127, 154907 (2007)]. The model, originally developed for athermal blends, is here modified to take into account the non-negligible interaction between polymer and plasticizer. We found that the temperature-pressure dependence of the alpha-relaxation time is very well described by this AG extended model.

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.

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.

A thermodynamic approach to the fragility of glass-forming polymers

We have connected the dynamic fragility, namely, the steepness of the relaxation-time variation upon temperature reduction, to the excess entropy and heat capacity of a large number of glass-forming polymers. The connection was obtained in a natural way from the Adam-Gibbs equation, relating the structural relaxation time to the configurational entropy. We find a clear correlation for a group of polymers. For another group of polymers, for which this correlation does not work, we emphasize the role of relaxation processes unrelated to the alpha process in affecting macroscopic thermodynamic properties. Once the residual excess entropy at the Vogel temperature is removed from the total excess entropy, the correlation between dynamic fragility and thermodynamic properties is reestablished.

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.

Why liquids are fragile

The fragilities (T(g)-normalized temperature dependence of alpha-relaxation times) of 33 glass-forming liquids and polymers are compared for isobaric, mP, and isochoric, mV, conditions. We find that the two quantities are linearly correlated: mP = (37+/-3) + (0.84+/-0.05)mV. This result has obvious and important consequences, since the ratio mV/mP is a measure of the relative degree to which temperature and density control the dynamics. Moreover, we show that the fragility itself is a consequence of the relative interplay of temperature and density effects near T(g). Specifically, strong behavior reflects a substantial contribution from density (jammed dynamics), while the relaxation of fragile liquids is more thermally activated. Drawing on the scaling law log(tau) = I(T upsilon(gamma)), a physical interpretation of this result in terms of the intermolecular potential is offered.

Effect of chemical structure on the isobaric and isochoric fragility in polychlorinated biphenyls

Pressure-volume-temperature data, along with dielectric relaxation measurements, are reported for a series of polychlorinated biphenyls (PCB), differing in the number of chlorine atoms on their phenyl rings. Analysis of the results reveals that with increasing chlorine content, the relaxation times of the PCB become governed to a greater degree by density rho relative to the effect of temperature T. This result is consistent with the respective magnitudes of the scaling exponent gamma yielding superpositioning of the relaxation times measured at various temperatures and pressures, when plotted versus rho(gamma)/T. While at constant (atmospheric) pressure, fragilities for the various PCB are equivalent, the fragility at constant volume varies inversely with chlorine content. Evidently, the presence of bulkier chlorine atoms on the phenyl rings magnifies the effect which the density has on the relaxation dynamics.

Disentangling density and temperature effects in the viscous slowing down of glassforming liquids

We present a consistent picture of the respective role of density (rho) and temperature (T) in the viscous slowing down of glassforming liquids and polymers. Specifically, based in part upon a new analysis of simulation and experimental data on liquid ortho-terphenyl, we conclude that a zeroth-order description of the approach to the glass transition (in the range of experimentally accessible pressures) should be formulated in terms of a temperature-driven super-Arrhenius activated behavior rather than a density-driven congestion or jamming phenomenon. The density plays a role at a quantitative level, but its effect on the viscosity and the alpha-relaxation time can be simply described via a single parameter, an effective interaction energy that is characteristic of the high-T liquid regime; as a result, rho does not affect the "fragility" of the glassforming system.

Influence of density and temperature on the microscopic structure and the segmental relaxation of polybutadiene

We investigate the influence of temperature and density on the local structure and the dynamics of polybutadiene by controlling both hydrostatic pressure and temperature in polarized neutron diffraction experiments on deuterated polybutadiene and in inelastic incoherent scattering experiments on protonated polybutadiene. We observe that the static structure factor S(Q) does not change along macroscopic isochores. This behavior is contrary to the relaxations observed on the nanosecond and picosecond time scales and viewed by the dynamic incoherent scattering law S(Q,omega), which differ strongly along the same thermodynamic path. We conclude that the static behavior, i.e., S(Q), is dominated by macroscopic density changes, similar to the vibrational excitations in the meV range. However, the relaxation dynamics is more sensitive to thermal energy changes. This is confirmed by the finding that lines of identical relaxation behavior (in time, shape, and Q dependence), isochrones on the 10(-9) sec time scale, clearly cross the constant density lines in the (P,T) plane. Concerning S(Q), we can reasonably relate the variation of the main-peak position to the average neighbor chain distance and deduce crude microscopic thermal expansion and compressibility coefficients. In the low-Q regime, the observed pressure and temperature variation of S(Q) exceeds the compressibility contribution and suggests the existence of additional scattering, which might originate from structural correlations arising at higher temperature and low pressure.

Influence of molecular structure on the dynamics of supercooled van der Waals liquids

Dielectric spectroscopy was carried out on the van der Waals liquid, 1,1(')-di(4-methoxy-5-methylphenyl)cyclohexane (BMMPC) in the supercooled state at pressures up to 218 MPa. The excess wing in this type-A glass former exhibits a response to pressure and temperature changes that is identical to that of the primary structural relaxation peak, indicating that the two processes reflect correlated molecular motions. Under no conditions was a distinct secondary peak observed in BMMPC, unlike the structurally very similar BMPC [1,1(')-bis(p-methoxyphenyl)cyclohexane]. However, the pressure dependences of both the glass temperature and fragility for the two materials are very close. The fragility is a decreasing function of pressure, although there is no concomitant narrowing of the relaxation peak. The pressure dependence of the relaxation times could be described as a simple volume-activated process, with the activation volume at the glass transition having the same magnitude as the molar volume.

Local structure and glass transition of polybutadiene up to 4 GPa

This communication presents a determination of the glass transition of polybutadiene under very high pressure, and raises the problem of the determination of the relative effects of temperature and density on the glass transition, depending on the pressure and temperature conditions. Local structure and slow dynamics were studied, by neutron scattering and calorimetry. To the best of our knowledge in neutron diffraction on soft matter such a high pressure, up to 4 GPa, was achieved.