Transition from Arrhenius to non-Arrhenius temperature dependence of structural relaxation time in glass-forming liquids: continuous versus discontinuous scenario

Normal-to-Supercooled Liquid Transition in Molecular Glass-Formers: A Hidden Structural Transformation Fuelled by Conformational Interconversion

Molecular dynamics and transport coefficients change significantly around the so-called Arrhenius crossover in glass-forming systems. In this article, we revisit the dynamic processes occurring in a glass-forming macrocyclic crown thiaether MeBzS2O above its glass transition, revealing two crossover temperatures: TB at 309 and TA at 333 K. We identify the second one as the Arrhenius crossover that is closely related to the normal-to-supercooled liquid transition in this compound. We show that the transformation occurring at this point goes far beyond molecular dynamics (where the temperature dependence of structural relaxation times changes its character from activation-like to super-Arrhenius), being reflected also in the internal structure and diffraction pattern. In this respect, we found a twofold local organization of the nearest-neighbor molecules via weak van der Waals forces, without the formation of any medium-range order or mesophases. The nearest surrounding of each molecule evolves structurally in time due to the ongoing fast conformational changes. We identify several conformers of MeBzS2O, demonstrating that its lowest-energy conformation is preferred mainly at lower temperatures, i.e., in the supercooled liquid state. Its increased prevalence modifies locally the short-range intermolecular order and promotes vitrification. Consequently, we indicate that the Arrhenius transition is fuelled rather by conformational changes in this glass-forming macrocyclic crown thiaether, which is a different scenario from the so-far existing concepts. Our studies combine broadband dielectric spectroscopy (BDS), X-ray diffraction, Fourier transform infrared (FTIR) spectroscopy, molecular dynamics (MD) simulations, and density functional theory (DFT) calculations.

Revisiting Dynamic Processes and Relaxation Mechanisms in a Heterocyclic Glass-Former: Direct Observation of a Transient State

Despite decades of studies, a clear understanding of near-Tg phenomena remains challenging for glass-forming systems. This review delves into the intricate molecular dynamics of the small, heterocyclic thioether, 6-methyl-2,3-dihydro-1,4-benzodithiine (MeBzS2), with a particular focus on its near-Tg cold crystallization and relaxation mechanisms. Investigating isothermal crystallization kinetics at various temperatures reveals a significant interplay between its molecular dynamics and recrystallization from a supercooled liquid. We also identify two independent interconversion paths between energetically privileged conformers, characterized by strained transition states. We demonstrate that these spatial transformations induce substantial alterations in the dipole moment orientation and magnitude. Our investigation also extends to the complex salt PdCl2(MeBzS2), where we observe the transient conformers directly, revealing a direct relationship between their abundance and the local or macroscopic electric field. The initially energetically privileged isomers in an undisturbed system become less favored in the presence of an external electric field or ions, resulting even in an unexpected inversion of states. Consequently, we confirm the intramolecular character of secondary relaxation in MeBzS2 and its mechanism related to conformational changes within the heterocyclic ring. The research is based on the combination of broadband dielectric spectroscopy, X-ray diffraction, and quantum density functional theory calculations.

Arrhenius Crossover Temperature of Glass-Forming Liquids Predicted by an Artificial Neural Network

The Arrhenius crossover temperature, TA, corresponds to a thermodynamic state wherein the atomistic dynamics of a liquid becomes heterogeneous and cooperative; and the activation barrier of diffusion dynamics becomes temperature-dependent at temperatures below TA. The theoretical estimation of this temperature is difficult for some types of materials, especially silicates and borates. In these materials, self-diffusion as a function of the temperature T is reproduced by the Arrhenius law, where the activation barrier practically independent on the temperature T. The purpose of the present work was to establish the relationship between the Arrhenius crossover temperature TA and the physical properties of liquids directly related to their glass-forming ability. Using a machine learning model, the crossover temperature TA was calculated for silicates, borates, organic compounds and metal melts of various compositions. The empirical values of the glass transition temperature Tg, the melting temperature Tm, the ratio of these temperatures Tg/Tm and the fragility index m were applied as input parameters. It has been established that the temperatures Tg and Tm are significant parameters, whereas their ratio Tg/Tm and the fragility index m do not correlate much with the temperature TA. An important result of the present work is the analytical equation relating the temperatures Tg, Tm and TA, and that, from the algebraic point of view, is the equation for a second-order curved surface. It was shown that this equation allows one to correctly estimate the temperature TA for a large class of materials, regardless of their compositions and glass-forming abilities.

Gigahertz elastic modulus and OH stretching frequency correlate with Jones-Dole's B-coefficient in aqueous solutions of the Hofmeister series

The ability of salts to change the macroscopic viscosity of their aqueous solutions is described by the Jones-Dole equation with B-coefficient for the linear concentration term. The sign and value of this coefficient are often considered as a measure of the salt's structure-making/breaking ability, while the validity of this assignment is still under discussion. Here, by applying Raman and Brillouin scattering spectroscopy to various salts from the Hofmeister series, we studied a possible relation between macroscopic Jones-Dole's B-coefficient and the microscopic dynamic response. Raman spectroscopy provides information about molecular vibrations and Brillouin spectroscopy about acoustic phonons with wavelengths of hundreds of nanometers. It has been found that Jones-Dole's B-coefficient correlates linearly with the coefficients, describing the concentration dependences of the average OH stretching frequency, real and imaginary parts of gigahertz elastic modulus. These relationships have been interpreted to mean that the OH stretching frequency is a measure of the ion-induced changes in the water network that cause changes in both viscosity and gigahertz relaxation. Depolarized inelastic light scattering revealed that the addition of structure-making ions not only changes the frequency of the relaxation peak but also increases the low-frequency part of the relaxation susceptibility. It was shown that the ion-induced increase in the gigahertz elastic modulus can be described by changes in the relaxational susceptibility without a noticeable change in the instantaneous elastic modulus. The isotropic Raman contribution associated with the tetrahedral-like environment of H₂O molecule does not correlate with Jones-Dole's B-coefficient, suggesting a minor influence of these tetrahedral-like configurations on viscosity.

On the Mutual Relationships between Molecular Probe Mobility and Free Volume and Polymer Dynamics in Organic Glass Formers: cis-1,4-poly(isoprene)

We report on the reorientation dynamics of small spin probe 2,2,6,6-tetramethylpiperidinyl-1-oxyl (TEMPO) in cis-1,4-poly(isoprene) (cis-1,4-PIP10k) from electron spin resonance (ESR) and the free volume of cis-1,4-PIP10k from positron annihilation lifetime spectroscopy (PALS) in relation to the high-frequency relaxations of cis-1,4-PIP10k using light scattering (LS) as well as to the slow and fast processes from broadband dielectric spectroscopy (BDS) and neutron scattering (NS). The hyperfine coupling constant, 2A_zz '(T), and the correlation times, τ _c(T), of cis-1,4-PIP10k/TEMPO system as a function of temperature exhibit several regions of the distinct spin probe TEMPO dynamics over a wide temperature range from 100 K up to 350 K. The characteristic ESR temperatures of changes in the spin probe dynamics in cis-1,4-PIP10k/TEMPO system are closely related to the characteristic PALS ones reflecting changes in the free volume expansion from PALS measurement. Finally, the time scales of the slow and fast dynamics of TEMPO in cis-1,4-PIP10k are compared with all of the six known slow and fast relaxation modes from BDS, LS and NS techniques with the aim to discuss the controlling factors of the spin probe reorientation mobility in polymer, oligomer and small molecular organic glass-formers.

Universal correlations between the fragility and interparticle repulsion of glass-forming liquids

A recently published analytical model describing and predicting elasticity, viscosity, and fragility of metallic melts is applied for the analysis of about 30 nonmetallic glassy systems, ranging from oxide network glasses to alcohols, low-molecular-weight liquids, polymers, plastic crystals, and even ionic glass formers. The model is based on the power-law exponent λ representing the steepness parameter of the repulsive part of the inter-atomic or inter-molecular potential and the thermal-expansion parameter α_T determined by the attractive anharmonic part of the effective interaction. It allows fitting the typical super-Arrhenius temperature variation of the viscosity or dielectric relaxation time for various classes of glass-forming matter, over many decades. We discuss the relation of the model parameters found for all these different glass-forming systems to the fragility parameter m and detect a correlation of λ and m for the non-metallic glass formers, in accord with the model predictions. Within the framework of this model, the fragility of glass formers can be traced back to microscopic model parameters characterizing the intermolecular interactions.

On Viscous Flow in Glass-Forming Organic Liquids

The two-exponential Sheffield equation of viscosity η(T) = A1·T·[1 + A2·exp(Hm/RT)]·[1 + C·exp(Hd/RT)], where A1, A2, Hm, C, and Hm are material-specific constants, is used to analyze the viscous flows of two glass-forming organic materials-salol and α-phenyl-o-cresol. It is demonstrated that the viscosity equation can be simplified to a four-parameter version: η(T) = A·T·exp(Hm/RT)]·[1 + C·exp(Hd/RT)]. The Sheffield model gives a correct description of viscosity, with two exact Arrhenius-type asymptotes below and above the glass transition temperature, whereas near the Tg it gives practically the same results as well-known and widely used viscosity equations. It is revealed that the constants of the Sheffield equation are not universal for all temperature ranges and may need to be updated for very high temperatures, where changes occur in melt properties leading to modifications of A and Hm for both salol and α-phenyl-o-cresol.

Second-order-derivative analysis of structural relaxation time in the elastic model of glass-forming liquids

Recently it was shown [V. N. Novikov and A. P. Sokolov, Phys. Rev. E 92, 062304 (2015)10.1103/PhysRevE.92.062304] that the second derivative with respect to inverse temperature of the structural relaxation time in some supercooled molecular liquids has a sharp maximum. It marks the point at which the apparent activation energy begins to saturate with decreasing temperature. The elastic model of glass-forming liquids expresses the temperature dependence of the structural relaxation time through that of the shear modulus. In this paper, we test whether this model is able to predict the maximum of the second derivative. We confirm its presence in the elastic model by analyzing the temperature dependence of the Brillouin light scattering in salol. This is a very subtle feature of the temperature dependence, which is greatly enhanced when taking derivatives. Its presence in the Brillouin data provides strong support to the elastic model of glass-forming liquids.

Composition- and temperature-dependent liquid structures in Al-Cu alloys: an ab initio molecular dynamics and x-ray diffraction study

The composition- and temperature-dependent liquid structures in eight Al_rich-Cu binary alloys (from hypoeutectic Al₉₃Cu₇ to hypereutectic Al₇₀Cu₃₀) have been experimentally and computationally studied by x-ray diffraction (XRD) experiments and ab initio molecular dynamics (AIMD) simulations. The remarkable agreements of structure factors for all liquid Al_rich-Cu alloys obtained from high-temperature high-energy XRD measurements and AIMD simulations have been achieved, which consolidates the analyses of structural evolutions in Al_rich-Cu liquids during the cooling processing by AIMD simulations. The heat capacity of liquid Al_rich-Cu alloys continuously increases and presents no abnormal peak when reducing the temperature, which differs from the reported prediction for 55-atom Al_rich-Cu nanoliquids. The diffusivities of Al and Cu undergo an increasing deviation from Arrhenius behavior by tuning Cu concentration from 7 to 30 atomic percentages, correlated to the local ordering in these liquids by means of coordination number, bond-angle distribution, Honeycutt-Andersen index, bond-orientational order and Voronoi tessellation analyses. Upon cooling, the microstructure of the liquid Al_rich-Cu alloys inclines to form Al₂Cu crystal-like local atomic ordering, especially in the hypereutectic liquids. The favorable short-range ordering between Cu and Al atoms could cause the non-Arrhenius diffusion behavior.

Qualitative change in structural dynamics of some glass-forming systems

Analysis of the temperature dependence of the structural relaxation time τ(α)(T) in supercooled liquids revealed a qualitatively distinct feature-a sharp, cusplike maximum in the second derivative of logτ(α)(T)at some T(max). It suggests that the super-Arrhenius temperature dependence of τ(α)(T) in glass-forming liquids eventually crosses over to an Arrhenius behavior at T<T(max), and there is no divergence of τ(α)(T) at nonzero T. T(max) can be above or below T(g), depending on the sensitivity of τ(T) to a change in the liquid's density quantified by the exponent γ in the scaling τ(α)(T)∼exp(A/Tρ(-γ)). These results might turn the discussion of the glass transition in a different direction-toward the origin of the limiting activation energy for structural relaxation at low T.

Computer simulations of glasses: the potential energy landscape

We review the current state of research on glasses, discussing the theoretical background and computational models employed to describe them. This article focuses on the use of the potential energy landscape (PEL) paradigm to account for the phenomenology of glassy systems, and the way in which it can be applied in simulations and the interpretation of their results. This article provides a broad overview of the rich phenomenology of glasses, followed by a summary of the theoretical frameworks developed to describe this phenomonology. We discuss the background of the PEL in detail, the onerous task of how to generate computer models of glasses, various methods of analysing numerical simulations, and the literature on the most commonly used model systems. Finally, we tackle the problem of how to distinguish a good glass former from a good crystal former from an analysis of the PEL. In summarising the state of the potential energy landscape picture, we develop the foundations for new theoretical methods that allow the ab initio prediction of the glass-forming ability of new materials by analysis of the PEL.