50th Anniversary Perspective: Challenges in the Dynamics and Kinetics of Glass-Forming Polymers

How to Distinguish Nonexponentiality and Nonlinearity in Isothermal Structural Relaxation of Glass-Forming Materials

The nonexponentiality and nonlinearity are two essential features of the structural relaxation in any glass-forming material, which seem to be inextricably bound together by the material time. It is shown that the temperature down-jump and up-jump experiments of the same magnitude ΔT = T0 - T to the same temperature T provide a clue for their separation. The isothermal structural relaxation can be quantified using the stabilization period on the logarithmic time scale log(tm/t0). It is described as the sum of the nonexponentiality term 1.181/ß and the nonlinearity term (σ/2.303)ΔT for the temperature down-jump, and as their difference for the temperature up-jump. The material parameter σ = -(∂lnτ/∂Tf)i quantifies variation of the relaxation time with structural changes at the inflection point of the relaxation curve and is formulated for the most widely used phenomenological models. The asymmetry of approach to equilibrium after the temperature down-jump and up-jump was first described by Kovacs in 1963. A detailed analysis of this asymmetry is provided, and a simple method for the estimation of the parameters characterizing the nonexponentiality (ß) and nonlinearity (σ) is proposed. The applicability of this method is tested using previously reported isothermal experimental data as well as calculated data for aging of polymers and other glass-forming materials. This concept illuminates differences in structural relaxation kinetics in a simple and consistent way that can be useful in the design of novel materials and the evaluation of their physical aging treatment.

Are Critical Fluctuations Responsible for Glass Formation?

The dynamic heterogeneities occurring just before the transition to the glassy phase have been named as the cause of amorphization in supercooled systems. Numerous studies conducted so far have confirmed this hypothesis, and based on it, a widely accepted solution to the puzzle of glass transition has been developed. This report focuses on verifying the existence of a strong pretransitional anomaly near the glass transition Tg. For this purpose, supercooled liquid-crystalline systems with a strong rod-like structure were selected. Based on the obtained experimental data, we demonstrate in this article that the previously postulated dynamic heterogeneities exhibit a critical characteristic, meaning a strong pretransitional anomaly can be observed with the described critical exponent α=0.5. Due to this property, it can be concluded that these heterogeneities are critical fluctuations, and consequently, the transition to the glassy state can be described based on the theory of critical phenomena. To measure the pretransitional anomaly near Tg in supercooled liquid-crystalline systems, broadband dielectric spectroscopy (BDS) and nonlinear dielectric effect (NDE) methods were applied. The exponent α provides insight into the nature and intensity of critical fluctuations in the system. A value of α=0.5 suggests that the fluctuations become increasingly intense as the system approaches the critical point, contributing to the divergence in specific heat. Understanding the role of critical fluctuations in the glass transition is crucial for innovating and improving a wide range of materials for energy storage, materials design, biomedical applications, food preservation, and environmental sustainability. These advancements can lead to materials with superior properties, optimized manufacturing processes, and applications that meet the demands of modern technology and sustainability challenges.

Nonlinear susceptibilities and higher-order responses related to physical aging: Wiener-Volterra approach and extended Tool-Narayanaswamy-Moynihan models

Large-amplitude thermal excursions imposed on deeply supercooled liquids modulate the nonlinear time evolution of their structural rearrangements. The consequent aftereffects are treated within a Wiener-Volterra expansion in laboratory time that allows one to calculate the associated physical-aging and thermal response functions. These responses and the corresponding higher-harmonic susceptibilities are illustrated using calculations based on the Tool-Narayanaswamy-Moynihan (TNM) model. The conversion from laboratory to material time is thoroughly discussed. Similarities and differences to field-induced higher-harmonic susceptibilities are illustrated using Lissajous and Cole-Cole plots and discussed in terms of aging nonlinearity parameters. For the Lissajous plots, banana-type shapes emerge, while the Cole-Cole plots display cardioidic and other visually appealing patterns. For application beyond the regime in which conventional single-parameter aging concepts work, the Wiener-Volterra material-time-series is introduced as the central tool. Calculations and analyses within this general framework in conjunction with suitable choices of higher-order memory kernels and employing correspondingly extended TNM models yield at least qualitative agreement with recent large-perturbation physical aging experiments. Implications for differential scanning calorimetry and related methods are discussed. The introduced concepts and analyses provide a solid foundation for a generalized description of nonlinear thermal out-of-equilibrium dynamics of glass forming materials, differing from the nonlinear responses known from rheology and dielectric spectroscopy.

Surface mobility gradient and emergent facilitation in glassy films

Confining glassy polymers into films can substantially modify their local and film-averaged properties. We present a lattice model of film geometry with void-mediated facilitation behaviors but free from any elasticity effect. We analyze the spatially varying viscosity to delineate the transport properties of glassy films. The film mobility measurements reported by Yang et al., Science, 2010, 328, 1676 are successfully reproduced. The flow exhibits a crossover from a simple viscous flow to a surface-dominated regime as the temperature decreases. The propagation of a highly mobile front induced by the free surface is visualized in real space. Our approach provides a microscopic treatment of the observed glassy phenomena.

The glass transition and enthalpy recovery of polystyrene nanorods using Flash differential scanning calorimetry

The glass transition (Tg) behavior and enthalpy recovery of polystyrene nanorods within an anodic aluminum oxide (AAO) template (supported nanorods) and after removal from AAO (unsupported nanorods) is studied using Flash differential scanning calorimetry. Tg is found to be depressed relative to the bulk by 20 ± 2 K for 20 nm-diameter unsupported polystyrene (PS) nanorods at the slowest cooling rate and by 9 ± 1 K for 55 nm-diameter rods. On the other hand, bulk-like behavior is observed in the case of unsupported 350 nm-diameter nanorods and for all supported rods in AAO. The size-dependent Tg behavior of the PS unsupported nanorods compares well with results for ultrathin films when scaled using the volume/surface ratio. Enthalpy recovery was also studied for the 20 and 350 nm unsupported nanorods with evolution toward equilibrium found to be linear with logarithmic time. The rate of enthalpy recovery for the 350 nm rods was similar to that for the bulk, whereas the rate of recovery was enhanced for the 20 nm rods for down-jump sizes larger than 17 K. A relaxation map summarizes the behavior of the nanorods relative to the bulk and relative to that for the 20 nm-thick ultrathin film. Interestingly, the fragility of the 20 nm-diameter nanorod and the 20 nm ultrathin film are identical within the error of measurements, and when plotted vs departure from Tg (i.e., T - Tg), the relaxation maps of the two samples are identical in spite of the fact that the Tg is depressed 8 K more in the nanorod sample.

Solid-that-Flows Picture of Glass-Forming Liquids

This perspective article reviews arguments that glass-forming liquids are different from those of standard liquid-state theory, which typically have a viscosity in the mPa·s range and relaxation times on the order of picoseconds. These numbers grow dramatically and become 1012 - 1015 times larger for liquids cooled toward the glass transition. This translates into a qualitative difference, and below the "solidity length" which is roughly one micron at the glass transition, a glass-forming liquid behaves much like a solid. Recent numerical evidence for the solidity of ultraviscous liquids is reviewed, and experimental consequences are discussed in relation to dynamic heterogeneity, frequency-dependent linear-response functions, and the temperature dependence of the average relaxation time.

Strain-driven Kovacs-like memory effect in glasses

Studying complex relaxation behaviors is of critical importance for understanding the nature of glasses. Here we report a Kovacs-like memory effect in glasses, manifested by non-monotonic stress relaxation during two-step high-to-low strains stimulations. During the stress relaxation process, if the strain jumps from a higher state to a lower state, the stress does not continue to decrease, but increases first and then decreases. The memory effect becomes stronger when the atomic motions become highly collective with a large activation energy, e.g. the strain in the first stage is larger, the temperature is higher, and the stimulation is longer. The physical origin of the stress memory effect is studied based on the relaxation kinetics and the in-situ synchrotron X-ray experiments. The stress memory effect is probably a universal phenomenon in different types of glasses.

Highly stable petroleum pitches provide access to the deep glassy state

Differential scanning calorimetry (DSC) was used to study the fast aging behavior of two petroleum pitch materials despite being only three to five years old. We observe that these highly aromatic pitches with broad distributions of both molecular weight and aromaticity exhibit large enthalpic relaxation endotherms in initial DSC heating scans, and 20-32 °C reductions in the fictive temperature and 0.35-0.87 of θK, which are indicative of aged glasses similar to ultrastable glasses and 20 MA aged amber. Quantifying the degree of thermodynamic stability relative to the Kauzmann temperature vs. the aging time demonstrates that these materials age just as quickly as low fragility metallic glasses. Additionally, we observe that pitches age faster than polymers reported in the literature when compared using down-jump experiments. We hypothesize that the fraction of higher aromaticity of pitch molecules plays a crucial role in faster dynamics. The unique aging behavior and the ability to produce pitches in bulk quantities using pilot-scale equipment, while being possible to tailor their molecular composition, make them a useful material for studying complex aging dynamics in the deep glassy state.

How the Presence of Crystalline Phase Affects Structural Relaxation in Molecular Liquids: The Case of Amorphous Indomethacin

The influence of partial crystallinity on the structural relaxation behavior of low-molecular organic glasses is, contrary to, e.g., polymeric materials, a largely unexplored territory. In the present study, differential scanning calorimetry was used to prepare a series of amorphous indomethacin powders crystallized to various extents. The preparations stemmed from the two distinct particle size fractions: 50-125 µm and 300-500 µm. The structural relaxation data from the cyclic calorimetric measurements were described in terms of the phenomenological Tool-Narayanaswamy-Moynihan model. For the 300-500 µm powder, the crystalline phase forming dominantly on the surface led to a monotonous decrease in the glass transition by ~6 °C in the 0-70% crystallinity range. The activation energy of the relaxation motions and the degree of heterogeneity within the relaxing matrix were not influenced by the increasing crystallinity, while the interconnectivity slightly increased. This behavior was attributed to the release of the quenched-in stresses and to the consequent slight increase in the structural interconnectivity. For the 50-125 µm powder, distinctly different relaxation dynamics were observed. This leads to a conclusion that the crystalline phase grows throughout the bulk glassy matrix along the internal micro-cracks. At higher crystallinity, a sharp increase in Tg, an increase in interconnectivity, and an increase in the variability of structural units engaged in the relaxation motions were observed.

Structure-Property Relationship, Glass Transition, and Crystallization Behaviors of Conjugated Polymers

Conjugated polymers have gained considerable interest due to their unique structures and promising applications in areas such as optoelectronics, photovoltaics, and flexible electronics. This review focuses on the structure-property relationship, glass transition, and crystallization behaviors of conjugated polymers. Understanding the relationship between the molecular structure of conjugated polymers and their properties is essential for optimizing their performance. The glass transition temperature (Tg) plays a key role in determining the processability and application of conjugated polymers. We discuss the mechanisms underlying the glass transition phenomenon and explore how side-chain interaction affects Tg. The crystallization behavior of conjugated polymers significantly impacts their mechanical and electrical properties. We investigate the nucleation and growth processes, as well as the factors that influence the crystallization process. The development of the three generations of conjugated polymers in controlling the crystalline structure and enhancing polymer ordering is also discussed. This review highlights advanced characterization techniques such as X-ray diffraction, atomic force microscopy, and thermal analysis, which provide insights into molecular ordering and polymer-crystal interfaces. This review provides an insight of the structure-property relationship, glass transition, and crystallization behaviors of conjugated polymers. It serves as a foundation for further research and development of conjugated polymer-based materials with enhanced properties and performance.

Bond Orientation-Determined Enthalpic Stress in Polymer Glasses upon Deformation

The mechanical properties of polymer glass are determined by both intermolecular local packing structures and aligned intrachain configurations. These configurations involve multiple space scales, and the underlying mechanism is not well understood yet. By applying mechanical stimulation to cold-drawn polymer glasses, the present simulation work shows a one-to-one correspondence between arising retractive stress and the segment orientation parameter on the length scale of the intrachain connecting bond. Such retractive stress is a newly produced enthalpic stress when segment orientation on the length scale of bonds and particle mobility coexist. This reveals a potential mechanism of how the intrachain orientation on the length scale of bonds influences the mechanical behaviors of polymer glasses.

Why Volume and Dynamics Decouple in Nanocomposite Matrices: Space that Cannot Be Accessed is Not Free

Polymer nanocomposites have important material applications and are an ongoing focus of many molecular level investigations, however, puzzling experimental results exist. For example, specific volumes for some polymer nanocomposite matrices are 2% to 4% higher than for the neat polymer; in a pure polymer melt this would correspond to a pressure change of 40 to 100 MPa, and a decrease in isothermal segmental relaxation times of 3 to 5 orders of magnitude. However, the nanocomposite segmental dynamics do not show any speed up. We can explain this apparent uncoupling of dynamics from specific volume, and the key is to consider the system expansivity, i.e., the temperature dependence of the volumetric data, together with the concept of limiting volume at close liquid packing. Using pressure, volume, temperature data as a path to both, we are able to predict the effect of nanoadditives on the accessible, i.e., free, space in the material, which is critical for facilitating molecular rearrangements in dense systems. Our analysis explains why an increase in specific volume in a material may not always lead to faster segmental dynamics.

Dynamic Heterogeneity and Kovacs' Memory Effects in Supercooled Water

Understanding the properties of supercooled water is important for developing a comprehensive theory for liquid water and amorphous ices. Because of rapid crystallization for deeply supercooled water, experiments on it are typically carried out under conditions in which the temperature and/or pressure are rapidly changing. As a result, information on the structural relaxation kinetics of supercooled water as it approaches (metastable) equilibrium is useful for interpreting results obtained in this experimentally challenging region of phase space. We used infrared spectroscopy and the fast time resolution obtained by transiently heating nanoscale water films to investigate relaxation kinetics (aging) in supercooled water. When the structural relaxation of the water films was followed using a temperature jump protocol analogous to the classic experiments of Kovacs, similar memory effects were observed. In particular, after suitable aging at one temperature, water's structure displayed an extremum versus the number of heat pulses upon changing to a second temperature before eventually relaxing to a steady-state structure characteristic of that temperature. A random double well model based on the idea of dynamic heterogeneity in supercooled water accounts for the observations.

From Subaging to Hyperaging in Structural Glasses

We demonstrate nonequilibrium scaling laws for the aging and equilibration dynamics in glass formers that emerge from combining a relaxation equation for the static structure with the equilibrium scaling laws of glassy dynamics. Different scaling regimes are predicted for the evolution of the structural relaxation time τ with age (waiting time t_{w}), depending on the depth of the quench from the liquid into the glass: "simple" aging (τ∼t_{w}) applies for quenches close to the critical point of mode-coupling theory (MCT) and implies "subaging" (τ≈t_{w}^{δ} with δ<1) as a broad equilibration crossover for quenches to nearly arrested equilibrium states; "hyperaging" (or superaging, τ∼t_{w}^{δ^{'}} with δ^{'}>1) emerges for quenches deep into the glass. The latter is cut off by non-mean-field fluctuations that we account for within a recent extension of MCT, the stochastic β-relaxation theory (SBR). We exemplify the scaling laws with a schematic model that quantitatively fits simulation data.

Combined description of pressure-volume-temperature and dielectric relaxation of several polymeric and low-molecular-weight organic glass-formers using SL-TS2 approach

We apply our recently-developed mean-field "SL-TS2" (two-state Sanchez-Lacombe) model to simultaneously describe dielectric α-relaxation time, τ_α, and pressure-volume-temperature (PVT) data in four polymers (polystyrene, poly(methylmethacrylate), poly(vinyl acetate) and poly(cyclohexane methyl acrylate)) and four organic molecular glass formers (ortho-terphenyl, glycerol, PCB-62, and PDE). Previously, it has been shown that for all eight materials, the Casalini-Roland thermodynamical scaling, τ_α = f(Tvγsp) (where T is temperature and v_sp is specific volume) is satisfied (R. Casalini and C. M. Roland, Phys. Rev. E, 2004, 69(6), 62501). It has also been previously shown that the same scaling emerges naturally (for sufficiently low pressures) within the "SL-TS2" framework (V. V. Ginzburg, Soft Matter, 2021, 17, 9094-9106). Here, we fit the ambient pressure curves for the relaxation time and the specific volume as functions of temperature for the eight materials and observe a good agreement between theory and experiment. We then use the Casalini-Roland scaling to convert those results into "master curves", thus enabling predictions of relaxation times and specific volumes at elevated pressures. The proposed approach can be used to describe other glass-forming materials, both low-molecular-weight and polymeric.

Distance-as-time in physical aging

Although it has been known for half a century that the physical aging of glasses in experiments is described well by a linear thermal-history convolution integral over the so-called material time, the microscopic definition and interpretation of the material time remains a mystery. We propose that the material-time increase over a given time interval reflects the distance traveled by the system's particles. Different possible distance measures are discussed, starting from the standard mean-square displacement and its inherent-state version that excludes the vibrational contribution. The viewpoint adopted, which is inspired by and closely related to pioneering works of Cugliandolo and Kurchan from the 1990s, implies a "geometric reversibility" and a "unique-triangle property" characterizing the system's path in configuration space during aging. Both of these properties are inherited from equilibrium, and they are here confirmed by computer simulations of an aging binary Lennard-Jones system. Our simulations moreover show that the slow particles control the material time. This motivates a "dynamic-rigidity-percolation" picture of physical aging. The numerical data show that the material time is dominated by the slowest particles' inherent mean-square displacement, which is conveniently quantified by the inherent harmonic mean-square displacement. This distance measure collapses data for potential-energy aging well in the sense that the normalized relaxation functions following different temperature jumps are almost the same function of the material time. Finally, the standard Tool-Narayanaswamy linear material-time convolution-integral description of physical aging is derived from the assumption that when time is replaced by distance in the above sense, an aging system is described by the same expression as that of linear-response theory.

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.

Relaxation Dynamics of Biomass-Derived Copolymers With Promising Gas-Barrier Properties

This article presents an experimental study on the relaxation dynamics of a series of random copolymers based on bio-friendly comonomers with interesting gas barrier properties. We analyze the relaxation response in the glassy and ultraviscous regime of poly (trimethylene furanoate/sebacate) random copolymers via dielectric spectroscopy. We report lower values of dynamic fragility [a dimensionless index introduced in 1985 (Angell, Relaxations in Complex Systems, 1985)] in comparison to popular polyesters widely used in industry, such as poly (ethylene terephthalate), suggesting that the amorphous phase of these furanoate-based polyesters adopt an efficient chain packing. This is consistent with their low permeability to gases. We also discuss on different equations (phenomenological and theory-based approaches) for fitting the temperature-evolution of the alpha relaxation time.

"Inner clocks" of glass-forming liquids

Providing a physically sound explanation of aging phenomena in non-equilibrium amorphous materialsis a challenging problem in modern statistical thermodynamics. The slow evolution of physical propertiesafter quenches of control parameters is empirically well interpreted via the concept of material time (orinternal clock), based on the Tool-Narayanaswamy-Moynihan (TNM) model. Yet, the fundamental reasonsof its striking success remain unclear. We propose a microscopic rationale behind the material time onthe basis of the linear laws of irreversible thermodynamics and its extension that treats the correspondingkinetic coefficients as state functions of a slowly evolving material state. Our interpretation is based onthe recognition that the same mathematical structure governs both the Tool model and the recently devel-oped non-equilibrium extension of the self-consistent generalized Langevin equation theory (NE-SCGLE),guided by the universal principles of Onsager's theory of irreversible processes. This identification opensthe way for a generalization of the material-time concept to aging systems where several relaxation modeswith very different equilibration processes must be considered, and partially frozen glasses manifest theappearance of partial ergodicity breaking, and hence materials with multiple very distinct inner clocks.

Temperature dependence of aging dynamics in highly non-equilibrium model polymer glasses

A central feature of the non-equilibrium glassy “state” is its tendency to age toward equilibrium, obeying signatures identified by Kovacs over 50 years ago. The origin of these signatures, their fate far from equilibrium and at high temperatures, and the underlying nature of the glassy “state” far from equilibrium remain unsettled. Here, we simulate physical aging of polymeric glasses, driven much farther from equilibrium and at much higher temperatures than possible in experimental melt-quenched glasses. While these glasses exhibit Kovacs’ signatures of glassy aging at sufficiently low temperatures, these signatures disappear above the onset TA of non-Arrhenius equilibrium dynamics, suggesting that TA demarcates an upper bound to genuinely glassy states. Aging times in glasses after temperature up-jumps are found to obey an Arrhenius law interpolating between equilibrium dynamics at TA and at the start of the temperature up-jump, providing a zero-parameter rule predicting their aging behavior and identifying another unrecognized centrality of TA to aging behavior. This differs qualitatively from behavior of our glasses produced by temperature down-jumps, which exhibit a fractional power law decoupling relation with equilibrium dynamics. While the Tool–Narayanaswamy–Moynihan model can predict the qualitative single-temperature behavior of these systems, we find that it fails to predict the disappearance of Kovacs signatures above TA and the temperature dependence of aging after large temperature up-jumps. These findings highlight a need for new theoretical insights into the aging behavior of glasses at ultra-high fictive temperatures and far from equilibrium.

Predicting nonlinear physical aging of glasses from equilibrium relaxation via the material time

The noncrystalline glassy state of matter plays a role in virtually all fields of materials science and offers complementary properties to those of the crystalline counterpart. The caveat of the glassy state is that it is out of equilibrium and therefore exhibits physical aging, i.e., material properties change over time. For half a century, the physical aging of glasses has been known to be described well by the material-time concept, although the existence of a material time has never been directly validated. We do this here by successfully predicting the aging of the molecular glass 4-vinyl-1,3-dioxolan-2-one from its linear relaxation behavior. This establishes the defining property of the material time. Via the fluctuation-dissipation theorem, our results imply that physical aging can be predicted from thermal-equilibrium fluctuation data, which is confirmed by computer simulations of a binary liquid mixture.

Isotope effects on the structural transformation and relaxation of deeply supercooled water

We have examined the structure of supercooled liquid D2O as a function of temperature between 185 and 255 K using pulsed laser heating to rapidly heat and cool the sample on a nanosecond timescale. The liquid structure can be represented as a linear combination of two structural motifs, with a transition between them described by a logistic function centered at 218 K with a width of 10 K. The relaxation to a metastable state, which occurred prior to crystallization, exhibited nonexponential kinetics with a rate that was dependent on the initial structural configuration. When the temperature is scaled by the temperature of maximum density, which is an isostructural point of the isotopologues, the structural transition and the non-equilibrium relaxation kinetics of D2O agree remarkably well with those for H2O.

The dynamics of freestanding films: predictions for poly(2-chlorostyrene) based on bulk pressure dependence and thoughtful sample averaging

In this paper we model the segmental relaxation in poly(2-chlorostyrene) 18 nm freestanding films, using only data on bulk samples to characterize the system, and predict film relaxation times (τ) as a function of temperature that are in semi-quantitative agreement with film data. The ability to translate bulk characterization into film predictions is a direct result of our previous work connecting the effects of free surfaces in films with those of changing pressure in the bulk. Our approach combines the Locally Correlated Lattice (LCL) equation of state for prediction of free volume values (V_free) at any given density (ρ), which are then used in the Cooperative Free Volume (CFV) rate model to predict τ(T, V_free). A key feature of this work is that we calculate the locally averaged density profile as a function of distance from the surface, ρ_av(z), using the CFV-predicted lengthscale, L_coop(z), over which rearranging molecular segments cooperate. As we have shown in the past, ρ_av(z) is significantly broader than the localized profile, ρ(z), which translates into a relaxation profile, τ(z), exhibiting a breadth that mirrors experimental and simulated results. In addition, we discuss the importance of averaging the log of position dependent relaxation times across a film sample (〈log τ(z)〉), as opposed to averaging the relaxation times, themselves, in order to best approximate a whole sample-averaged value that can be directly compared to experiment.

Combined description of polymer PVT and relaxation data using a dynamic "SL-TS2" mean-field lattice model

We develop a combined model to describe the pressure-volume-temperature (PVT) thermodynamics and the α- and β-relaxation time dynamics in glass-forming amorphous materials. The PVT results are described using a two-state modification of the Sanchez-Lacombe equation of state (SL-EoS). The minimization of the Sanchez-Lacombe free energy expression allows one to calculate the density (or specific volume) and the "solid fraction" as a function of temperature and pressure. The solid fraction, ψ(T,P), is then substituted into the TS2 mean-field model (V. V. Ginzburg, Soft Matter, 2020, 16, 810), and the equilibrium relaxation times and the glass transition temperature are calculated. Finally, we use a dynamic model of a Tool-Narayanaswami-Moynihan (TNM)-type to describe the density relaxation as a function of the cooling rate. We applied the new framework to describe experimental data for polystyrene and poly(methylmethacrylate) and found a good qualitative and quantitative agreement.

The effect of pyrolysis on the chemical, thermal and rheological properties of pitch

Pitch-based carbon fibers are of considerable interest as high-performance materials. There are reports over the last several decades detailing (i) methods of improving pitch-based carbon fiber performance, and (ii) reducing the cost of production via novel processing techniques. However, there remain considerable challenges in producing high-performance pitch-based carbon fibers consistently on an industrial scale. This is arguably due to the difficulty of scaling the melt-spinning process to compensate for variability in pitch feedstock quality and a lack of understanding of processing-structure-performance relationships. This work focuses on the early stages of heat treatment (pyrolysis) of isotropic pitch and its effect on the chemical, thermal, and rheological properties of the pitch, which help determine its processability. More specifically, we quantify significant changes in chemical structure, M_w, T_g, T_s, and shear and extensional rheology as a function of pyrolysis time at 400 °C. The extensional rheology, in particular, shows that the 'stretchability' of the pitch samples strongly depends on pyrolysis severity, and is important for characterizing 'drawability'. Using a novel analysis of the uniaxial stretching kinematics, we show an isothermal 'drawability window' that allows for the largest axial and radial Hencky strains at constant rate. We hypothesize that this extensional drawability window could facilitate the successful processing of pitch into high quality fiber, minimizing the trial-and-error approach currently used in the field.

Tuning the dynamic fragility of acrylic polymers by small molecules: the interplay of molecular structures

This report studied changes in the dynamic fragility (m) of poly(butyl methacrylate) (PBMA) by introducing guest hindered phenols capable of forming two or three intermolecular hydrogen bonds (inter-HBs) per molecule with the host polymer. The small molecules effectively decrease the m value, even if they apparently increase the glass transition temperature (T_g) of mixtures. The reduction in m was confirmed by enthalpy relaxation in two aspects: adding the guest molecule leads to a stronger cooling rate dependence of the limiting fictive temperature together with an apparent increase in aging rate of PBMA hybrids at low concentrations. By varying the molecule size and steric hindrance of the hydroxyl group on the hindered phenols, we clarified that m is primarily governed by the strength of inter-HB interactions, while the T_g value of mixtures depends on a combined effect of additive bulkiness and HB interaction. The anomalous dynamics was further rationalized not only by the HB-induced flexibility balance between side groups and backbone, but also by the reduction of cooperative rearranging sizes and alleviation of long-chain connectivity in such HB-driven hybrids.

Time-resolved thermal lens investigation of glassy dynamics in supercooled liquids: Theory and experiments

This work reports results on the simultaneous spectroscopy of the specific heat and thermal expansivity of glycerol by making use of a wideband time-resolved thermal lens (TL) technique. An analytical model is presented which describes TL transients in a relaxing system subjected to impulsive laser heating. Experimentally, a set of TL waveforms, from 1 ns to 20 ms, has been recorded for a glycerol sample upon supercooling, from 300 to 200 K. The satisfactory fitting of the TL signals to the model allows the assessment of relaxation strength and relaxation frequency of the two quantities up to sub-100 MHz, extending the specific heat and thermal expansion spectroscopy by nearly three and eight decades, respectively. Fragility values, extracted from the relaxation behavior of the specific heat and the thermal expansion coefficient, are found to be similar, despite a substantial difference in relaxation strength.

Supercooled melt structure and dynamics of single-chain nanoparticles: A computer simulation study

By using coarse-grained molecular dynamics simulations, we have investigated the structure and dynamics of supercooled single-chain cross-linked nanoparticle (SCNP) melts having a range of cross-linking degrees ϕ. We find a nearly linear increase in glass-transition temperature (T_g) with increasing ϕ. Correspondingly, we have also experimentally synthesized a series of polystyrene-based SCNPs and have found that the measured T_g estimated from differential scanning calorimetry is qualitatively consistent with the trend predicted by our simulation estimates. Experimentally, an increase in T_g as large as ΔT_g = 61 K for ϕ = 0.36 is found compared with their linear chain counterparts, indicating that the changes in dynamics with cross-links are quite appreciable. We attribute the increase in T_g to the enlarged effective hard-core volume and the corresponding reduction in the free volume of the polymer segments. Topological constraints evidently frustrate the local packing. In addition, the introduction of intra-molecular cross-linking bonds slows down the structural relaxation and simultaneously enhances the local coupling motion on the length scales within SCNPs. Consequently, a more pronounced dynamical heterogeneity (DH) is observed for larger ϕ, as quantified by measuring the dynamical correlation length through the four-point susceptibility parameter, χ₄. The increase in DH is directly related to the enhanced local cooperative motion derived from intra-molecular cross-linking bonds and structural heterogeneity derived from the cross-linking process. These results shed new light on the influence of intra-molecular topological constraints on the segmental dynamics of polymer melts.

On a fundamental description of the Kovacs' kinetic signatures in glass-forming systems

The time-evolution equation for the time-dependent static structure factor of the non-equilibrium self-consistent generalized Langevin equation (NE-SCGLE) theory was used to investigate the kinetics of glass-forming systems under isochoric conditions. The kinetics are studied within the framework of the fictive temperature (T_F) of the glassy structure. We solve for the kinetics of T_F(t) and the time-dependent structure factor and find that they are different but closely related by a function that depends only on temperature. Furthermore, we are able to solve for the evolution of T_F(t) in a set of temperature-jump histories referred to as the Kovacs' signatures. We demonstrate that the NE-SCGLE theory reproduces all the Kovacs' signatures, namely, intrinsic isotherm, asymmetry of approach, and memory effect. In addition, we extend the theory into largely unexplored, deep glassy state, regions that are below the notionally "ideal" glass temperature.

Future directions of chemical theory and computation

Theoretical and computational chemistry aims to develop chemical theory and to apply numerical computation and simulation to reveal the mechanism behind complex chemical phenomena via quantum theory and statistical mechanics. Computation is the third pillar of scientific research together with theory and experiment. Computation enables scientists to test, discover, and build models/theories of the corresponding chemical phenomena. Theoretical and computational chemistry has been advanced to a new era due to the development of high-performance computational facilities and artificial intelligence approaches. The tendency to merge electronic structural theory with quantum chemical dynamics and statistical mechanics is of increasing interest because of the rapid development of on-the-fly dynamic simulations for complex systems plus low-scaling electronic structural theory. Another challenging issue lies in the transition from order to disorder, from thermodynamics to dynamics, and from equilibrium to non-equilibrium. Despite an increasingly rapid emergence of advances in computational power, detailed criteria for databases, effective data sharing strategies, and deep learning workflows have yet to be developed. Here, we outline some challenges and limitations of the current artificial intelligence approaches with an outlook on the potential future directions for chemistry in the big data era.

Isothermal Crystallization and Time-Temperature Transformation of Amorphous Nifedipine: A Case of Polymorphism Formation and Conversion

Crystallization of active pharmaceutical ingredients (APIs) from the supercooled liquid state is an important issue in determining the stability of amorphous pharmaceutical dispersions. In the present study, the isothermal crystallization from the supercooled liquid state of the pharmaceutical compound nifedipine was investigated by both rheological and differential scanning calorimetry (DSC) measurements, and the crystallization kinetics was fitted to the Johnson-Mehl-Avrami (JMA) equation. Both the crystallization induction time and completion time from the two methods were used to construct the time-temperature-transformation (TTT) diagram for nifedipine. A model based on a modification of classical homogeneous nucleation and crystal growth theory was employed to fit the induction and completion time curves. Both DSC and rheological methods give similar results for the crystallization kinetics of the nifedipine. From the crystallization kinetics modeling, the solid-liquid interfacial surface tension σ_SL of nifedipine was estimated and the value was found to be consistent with prior results obtained from melting point depression measurements as a function of crystal size. Evidence is shown that for temperatures below 110 °C, at the early stage of nucleation, NIF first nucleates into the metastable β'-form and later converts into the stable α-form during the isothermal crystallization. We are also able to report the heat of fusion of the γ'-NIF based on the calorimetric experiments.

Waiting-time dependent non-equilibrium phase diagram of simple glass- and gel-forming liquids

Under numerous circumstances, many soft and hard materials are present in a puzzling wealth of non-equilibrium amorphous states, whose properties are not stationary and depend on preparation. They are often summarized in unconventional "phase diagrams" that exhibit new "phases" and/or "transitions" in which time, however, is an essential variable. This work proposes a solution to the problem of theoretically defining and predicting these non-equilibrium phases and their time-evolving phase diagrams, given the underlying molecular interactions. We demonstrate that these non-equilibrium phases and the corresponding non-stationary (i.e., aging) phase diagrams can indeed be defined and predicted using the kinetic perspective of a novel non-equilibrium statistical mechanical theory of irreversible processes. This is illustrated with the theoretical description of the transient process of dynamic arrest into non-equilibrium amorphous solid phases of an instantaneously quenched simple model fluid involving repulsive hard-sphere plus attractive square well pair interactions.

Experimental test of a predicted dynamics-structure-thermodynamics connection in molecularly complex glass-forming liquids

Understanding in a unified manner the generic and chemically specific aspects of activated dynamics in diverse glass-forming liquids over 14 or more decades in time is a grand challenge in condensed matter physics, physical chemistry, and materials science and engineering. Large families of conceptually distinct models have postulated a causal connection with qualitatively different "order parameters" including various measures of structure, free volume, thermodynamic properties, short or intermediate time dynamics, and mechanical properties. Construction of a predictive theory that covers both the noncooperative and cooperative activated relaxation regimes remains elusive. Here, we test using solely experimental data a recent microscopic dynamical theory prediction that although activated relaxation is a spatially coupled local-nonlocal event with barriers quantified by local pair structure, it can also be understood based on the dimensionless compressibility via an equilibrium statistical mechanics connection between thermodynamics and structure. This prediction is found to be consistent with observations on diverse fragile molecular liquids under isobaric and isochoric conditions and provides a different conceptual view of the global relaxation map. As a corollary, a theoretical basis is established for the structural relaxation time scale growing exponentially with inverse temperature to a high power, consistent with experiments in the deeply supercooled regime. A criterion for the irrelevance of collective elasticity effects is deduced and shown to be consistent with viscous flow in low-fragility inorganic network-forming melts. Finally, implications for relaxation in the equilibrated deep glass state are briefly considered.

A Simple New Way To Account for Free Volume in Glassy Dynamics: Model-Free Estimation of the Close-Packed Volume from PVT Data

In this article we focus on the important role of well-defined free volume (V_free) in dictating the structural relaxation times, τ, of glass-forming liquids and polymer melts. Our definition of V_free = V - V_hc, where V is the total system volume, means the use of V_free depends on determination of V_hc, the system's volume in the limiting closely packed state. Rejecting the historically compromised use of V_free as a dynamics-dependent fitting function, we have successfully applied a clear thermodynamics-based route to V_hc using the locally correlated lattice (LCL) model equation of state (EOS). However, in this work we go further and show that V_hc can be defined without the use of an equation of state by direct linear extrapolation of a V(T) high-pressure isobar down to zero temperature (T). The results from this route, tested on a dozen experimental systems, yield ln τ vs 1/V_free isotherms that are linear with T-dependent slopes, consistent with the general ln τ ∼ f(T) × (1/V_free) form of behavior we have previously described. This functional form also results by implementing a simple mechanistic explanation via the cooperative free volume (CFV) rate model, which assumes that dynamic relaxation is both thermally activated and that it requires molecular segmental cooperativity. With the degree of the latter, and thus the activation energy, being determined by the availability of free volume, the new route we demonstrate here for determination of V_free expands the potential for understanding and predicting local dynamic relaxation in glass-forming materials.

Single-parameter aging in a binary Lennard-Jones system

This paper studies physical aging by computer simulations of a 2:1 Kob-Andersen binary Lennard-Jones mixture, a system that is less prone to crystallization than the standard 4:1 composition. Starting from thermal-equilibrium states, the time evolution of the following four quantities is monitored by following up and down jumps in temperature: potential energy, virial, average squared force, and the Laplacian of the potential energy. Despite the fact that significantly larger temperature jumps are studied here than in typical similar experiments, to a good approximation, all four quantities conform to the single-parameter-aging scenario derived and validated for small jumps in experiments [T. Hecksher, N. B. Olsen, and J. C. Dyre, J. Chem. Phys. 142, 241103 (2015)]. As a further confirmation of single-parameter aging with a common material time for the four different quantities monitored, their relaxing parts are found to be almost identical for all temperature jumps.