Thermodynamic consequences of the kinetic nature of the glass transition

Role of anisotropy in understanding the molecular grounds for density scaling in dynamics of glass-forming liquids

Molecular Dynamics (MD) simulations of glass-forming liquids play a pivotal role in uncovering the molecular nature of the liquid vitrification process. In particular, much focus was given to elucidating the interplay between the character of intermolecular potential and molecular dynamics behaviour. This has been tried to achieve by simulating the spherical particles interacting via isotropic potential. However, when simulation and experimental data are analysed in the same way by using the density scaling approaches, serious inconsistency is revealed between them. Similar scaling exponent values are determined by analysing the relaxation times and pVT data obtained from computer simulations. In contrast, these values differ significantly when the same analysis is carried out in the case of experimental data. As discussed thoroughly herein, the coherence between results of simulation and experiment can be achieved if anisotropy of intermolecular interactions is introduced to MD simulations. In practice, it has been realized in two different ways: (1) by using the anisotropic potential of the Gay-Berne type or (2) by replacing the spherical particles with quasi-real polyatomic anisotropic molecules interacting through isotropic Lenard-Jones potential. In particular, the last strategy has the potential to be used to explore the relationship between molecular architecture and molecular dynamics behaviour. Finally, we hope that the results presented in this review will also encourage others to explore how 'anisotropy' affects remaining aspects related to liquid-glass transition, like heterogeneity, glass transition temperature, glass forming ability, etc.

Water Promotes Melting of a Metal-Organic Framework

Water is one of the most reactive and abundant molecules on Earth, and it is thus crucial to understand its reactivity with various material families. One of the big unknown questions is how water in liquid and vapor forms impact the fast-emerging class of metal-organic frameworks (MOFs). Here, we discover that high-pressure water vapor drastically modifies the structure and hence the dynamic, thermodynamic, and mechanical properties of MOF glasses. In detail, we find that an archetypical MOF (ZIF-62) is extremely sensitive to heat treatments performed at 460 °C and water vapor pressures up to ∼110 bar. Both the melting and glass transition temperatures decrease remarkably (by >100 °C), and simultaneously, hardness and Young's modulus increase by up to 100% under very mild treatment conditions (<20 bar of hydrothermal pressure). Structural analyses suggest water to partially coordinate to Zn in the form of a hydroxide ion by replacing a bridging imidazolate-based linker. The work provides insight into the role of hot-compressed water in influencing the structure and properties of MOF glasses and opens a new route for systematically changing the thermodynamics and kinetics of MOF liquids and thus altering the thermal and mechanical properties of the resulting MOF glasses.

Density Scaling of Translational and Rotational Molecular Dynamics in a Simple Ellipsoidal Model near the Glass Transition

In this paper, we show that a simple anisotropic model of supercooled liquid properly reflects some density scaling properties observed for experimental data, contrary to many previous results obtained from isotropic models. We employ a well-known Gay–Berne model earlier parametrized to achieve a supercooling and glass transition at zero pressure to find the point of glass transition and explore volumetric and dynamic properties in the supercooled liquid state at elevated pressure. We focus on dynamic scaling properties of the anisotropic model of supercooled liquid to gain a better insight into the grounds for the density scaling idea that bears hallmarks of universality, as follows from plenty of experimental data collected near the glass transition for different dynamic quantities. As a result, the most appropriate values of the scaling exponent γ are established as invariants for a given anisotropy aspect ratio to successfully scale both the translational and rotational relaxation times considered as single variable functions of density^γ/temperature. These scaling exponent values are determined based on the density scaling criterion and differ from those obtained in other ways, such as the virial–potential energy correlation and the equation of state derived from the effective short-range intermolecular potential, which is qualitatively in accordance with the results yielded from experimental data analyses. Our findings strongly suggest that there is a deep need to employ anisotropic models in the study of glass transition and supercooled liquids instead of the isotropic ones very commonly exploited in molecular dynamics simulations of supercooled liquids over the last decades.

Coupling between structural relaxation and diffusion in glass-forming liquids under pressure variation

We theoretically investigate structural relaxation and activated diffusion of glass-forming liquids at different pressures using both Elastically Collective Nonlinear Langevin Equation (ECNLE) theory and molecular dynamics (MD) simulations.

The cooperative free volume rate model for segmental dynamics: Application to glass-forming liquids and connections with the density scaling approach⋆

In this paper, we apply the cooperative free volume (CFV) rate model for pressure-dependent dynamics of glass-forming liquids and polymer melts. We analyze segmental relaxation times, [Formula: see text] , as a function of temperature (T and free volume ( [Formula: see text] , and make substantive comparisons with the density scaling approach. [Formula: see text] , the difference between the total volume (V and the volume at close-packing, is predicted independently of the dynamics for any temperature and pressure using the locally correlated lattice (LCL) equation-of-state (EOS) analysis of characteristic thermodynamic data. We discuss the underlying physical motivation in the CFV and density scaling models, and show that their key, respective, material parameters are connected, where the CFV b parameter and the density scaling [Formula: see text] parameter each characterize the relative sensitivity of dynamics to changes in T , vs. changes in V . We find [Formula: see text] , where [Formula: see text] is the value predicted by the LCL EOS at the ambient [Formula: see text] . In comparing the predictive power of the two models we highlight the CFV advantage in yielding a universal linear collapse of relaxation data using a minimal set of parameters, compared to the same parameter space yielding a changing functional form in the density scaling approach. Further, we demonstrate that in the low data limit, where there is not enough data to characterize the density scaling model, the CFV model may still be successfully applied, and we even use it to predict the correct [Formula: see text] parameter.

Communication: Glass transition and melting lines of an ionic liquid

The phase diagram of the ionic liquid 1-butyl-1-methylpyrrolidinium bis(trifluoromethanesufonyl)imide, [Pyrr_1,4][NTf₂], was explored by synchroton X-ray diffraction and Raman scattering measurements as a function of temperature and pressure. Glass transition T_g(p) and melting T_m(p) temperatures were obtained from atmospheric pressure up to ca. 2.0 GPa. We found that both the T_g(p) and T_m(p) curves follow essentially the same pressure dependence. The similarity of pressure coefficients, dT_g/dp ≈ dT_m/dp, is explained within the non-equilibrium thermodynamics approach for the glass transition by assuming that one of the Ehrenfest equations is appropriated for T_g(p), whereas T_m(p) follows the Clausius-Clapeyron equation valid for the first-order transitions. The results highlight that ionic liquids are excellent model systems to address fundamental questions related to the glass transition.

Cooling-rate induced softening in a colloidal glass

Contrary to crystalline solids, amorphous solids always become softer when vitrifying the melts under higher cooling rates. Understanding this phenomenon is of utmost importance in providing a basis for the mechanical-performance control of amorphous solids. However, the underlying mechanisms leading to this cooling-rate-induced softening of amorphous solids have remained elusive, especially the dynamic reasons are neglected. Here, we use a colloidal glass as the model system to directly study this issue. Shear modulus is used as the representative parameter to monitor the stress-bearing properties of colloidal glass. The space-spanning immobile particles, whose population is sensitive to the cooling rate, are found to make the dominant contribution to the shear modulus. The rapid solidification induced softening of colloidal glass is observed to originate from fewer immobile particles formed at higher cooling rates.

In search of invariants for viscous liquids in the density scaling regime: investigations of dynamic and thermodynamic moduli

In this paper, we report the nontrivial results of our investigations of dynamic and thermodynamic moduli in search of invariants for viscous liquids in the density scaling regime by using selected supercooled van der Waals liquids as representative materials. Previously, the dynamic modulus M_p-T (defined in the pressure-temperature representation by the ratio of isobaric activation energy and activation volume) as well as the ratio B_T/M_p-T (where B_T is the thermodynamic modulus defined as the inverse isothermal compressibility) have been suggested as some kinds of material constants. We have established that they are not valid in the explored wide range of temperatures T over a dozen decades of structural relaxation times τ. The temperature dependences of M_p-T and B_T/M_p-T have been elucidated by comparison with the well-known measure of the relative contribution of temperature and density fluctuations to molecular dynamics near the glass transition, i.e., the ratio of isochoric and isobaric activation energies. Then, we have implemented an idea to transform the definition of the dynamic modulus M_p-T from the p-T representation to the V-T one. This idea relied on the disentanglement of combined temperature and density fluctuations involved in isobaric parameters and has resulted in finding an invariant for viscous liquids in the density scaling regime, which is the ratio of thermodynamic and dynamic moduli, B_T/M_V-T. In this way, we have constituted a characteristic of thermodynamics and molecular dynamics, which remains unchanged in the supercooled liquid state for a given material, the molecular dynamics of which obeys the power density scaling law.

Activation volume of selected liquid crystals in the density scaling regime

In this paper, we demonstrate and thoroughly analyze the activation volumetric properties of selected liquid crystals in the nematic and crystalline E phases in comparison with those reported for glass-forming liquids. In the analysis, we have employed and evaluated two entropic models (based on either total or configurational entropies) to describe the longitudinal relaxation times of the liquid crystals in the density scaling regime. In this study, we have also exploited two equations of state: volumetric and activation volumetric ones. As a result, we have established that the activation volumetric properties of the selected liquid crystals are quite opposite to such typical properties of glass-forming materials, i.e., the activation volume decreases and the isothermal bulk modulus increases when a liquid crystal is isothermally compressed. Using the model based on the configurational entropy, we suggest that the increasing pressure dependences of the activation volume in isothermal conditions and the negative curvature of the pressure dependences of isothermal longitudinal relaxation times can be related to the formation of antiparallel doublets in the examined liquid crystals. A similar pressure effect on relaxation dynamics may be also observed for other material groups in case of systems, the molecules of which form some supramolecular structures.

Predicting Nanoscale Dynamics of a Glass-Forming Liquid from Its Macroscopic Bulk Behavior and Vice Versa

The properties of a molecular liquid confined at the nanometer length scale can be very distinct from the bulk. For that reason, the macro- and the nanoscopic behaviors of glass-forming liquids are regarded as two nonconnected realms, governed by their own rules. Here, we show that the glassy dynamics in molecular liquids confined to nanometer pores might obey the density scaling relation, ρ^γ/T, just like in bulk fluids. Even more surprisingly, the same value of the scaling exponent γ superposes the α-relaxation time measured at different state points in nanoscale confinement and upon increased pressure. We report this remarkable finding for van der Waals liquids tetramethyl-tetraphenyl-trisiloxane (DC704) and polyphenyl ether (5PPE), considered as simple, single-parameter liquids. Demonstrating that the density scaling idea can be fulfilled in both environments opens an exciting possibility to predict the dynamic features of the nanoconfined system close to its glass-transition temperature from the high-pressure studies of the bulk liquid. Likewise, we can describe the viscous liquid dynamics at any given combination of temperature and pressure based on analysis of its behavior in nanopores.

Ultrastable glasses portray similar behaviour to ordinary glasses at high pressure

Pressure experiments provide a unique opportunity to unravel new insights into glass-forming liquids by exploring its effect on the dynamics of viscous liquids and on the evolution of the glass transition temperature. Here we compare the pressure dependence of the onset of devitrification, T_on, between two molecular glasses prepared from the same material but with extremely different ambient-pressure kinetic and thermodynamic stabilities. Our data clearly reveal that, while both glasses exhibit different dT_on/dP values at low pressures, they evolve towards closer calorimetric devitrification temperature and pressure dependence as pressure increases. We tentatively interpret these results from the different densities of the starting materials at room temperature and pressure. Our data shows that at the probed pressures, the relaxation time of the glass into the supercooled liquid is determined by temperature and pressure similarly to the behaviour of liquids, but using stability-dependent parameters.

The implications of various molecular interactions on the dielectric behavior of cimetidine and cimetidine hydrochloride

We employed broadband dielectric spectroscopy to characterize the molecular dynamics of cimetidine base and cimetidine hydrochloride, materials with similar structural skeletons but involving different molecular interactions (ionic vs. non-ionic).