Molecular Glass-forming Liquids Research Articles

Soft pinning: Experimental validation of static correlations in supercooled molecular glass-forming liquids.

Enormous enhancement in the viscosity of a liquid near its glass transition is a hallmark of glass transition. Within a class of theoretical frameworks, it is connected to growing many-body static correlations near the transition, often called "amorphous ordering." At the same time, some theories do not invoke the existence of such a static length scale in the problem. Thus, proving the existence and possible estimation of the static length scales of amorphous order in different glass-forming liquids is very important to validate or falsify the predictions of these theories and unravel the true physics of glass formation. Experiments on molecular glass-forming liquids become pivotal in this scenario as the viscosity grows several folds (), and simulations or colloidal glass experiments fail to access these required long-time scales. Here we design an experiment to extract the static length scales in molecular liquids using dilute amounts of another large molecule as a pinning site. Results from dielectric relaxation experiments on supercooled Glycerol with different pinning concentrations of Sorbitol and Glucose, as well as the simulations on a few model glass-forming liquids with pinning sites, indicate the versatility of the proposed method, opening possible new avenues to study the physics of glass transition in other molecular liquids.

Kinetics of Surface Enrichment of a Polymer in a Glass-Forming Molecular Liquid

X-ray photoelectron spectroscopy has been used to measure the surface concentration and the surface enrichment kinetics of a polymer in a glass-forming molecular liquid. As a model, the bulk-miscible system of maltitol-polyvinylpyrrolidone (PVP) was studied. The PVP concentration is significantly higher at the liquid/vapor interface than in the bulk by up to a factor of 170, and the effect increases with its molecular weight. At a freshly created liquid/vapor interface, the concentration of PVP gradually increases from the bulk value at a rate controlled by bulk diffusion. The polymer diffusion coefficient obtained from the kinetics of surface enrichment agrees with that calculated from viscosity and the Stokes-Einstein equation. Our finding allows prediction of the rate at which the surface composition equilibrates in an amorphous material after milling, fracture, and a change in ambient temperature.

Radical Diffusion Crossover Phenomenon in Glass-Forming Liquids.

We studied the diffusivities of a nitroxide radical at various temperatures in six glass-forming molecular liquids by electron spin resonance. By comparing the radical diffusivities and solvent self-diffusivities, we found that the radical diffusivities are lower than the self-diffusivities at high temperatures and approach them at low temperatures in all liquids. This crossover behavior was considered as evidence that a single-molecule diffusion process transforms into a collective process with temperature lowering. The crossover phenomenon was analyzed by a novel, simple diffusion model, combining collective and single-molecule diffusion processes, and it was compared to the Arrhenius crossover phenomenon. The obtained results suggest that future studies of tracer diffusion could contribute to a better understanding of diffusion mechanisms in glass-forming liquids. The proposed diffusion model could be used to study the crossover phenomena of tracer diffusion measured by other techniques, and it could serve as a base for developing more advanced models.

Comparing Microscopic and Macroscopic Dynamics in a Paradigmatic Model of Glass-Forming Molecular Liquid.

Glass transition is a most intriguing and long-standing open issue in the field of molecular liquids. From a macroscopic perspective, glass-forming systems display a dramatic slowing-down of the dynamics, with the inverse diffusion coefficient and the structural relaxation times increasing by orders of magnitude upon even modest supercooling. At the microscopic level, single-molecule motion becomes strongly intermittent, and can be conveniently described in terms of “cage-jump” events. In this work, we investigate a paradigmatic glass-forming liquid, the Kob–Andersen Lennard–Jones model, by means of Molecular Dynamics simulations, and compare the macroscopic and microscopic descriptions of its dynamics on approaching the glass-transition. We find that clear changes in the relations between macroscopic timescales and cage-jump quantities occur at the crossover temperature where Mode Coupling-like description starts failing. In fact, Continuous Time Random Walk and lattice model predictions based on cage-jump statistics are also violated below the crossover temperature, suggesting the onset of a qualitative change in cage-jump motion. Interestingly, we show that a fully microscopic relation linking cage-jump time- and length-scales instead holds throughout the investigated temperature range.

Long Wavelength Thermal Density Fluctuations in Molecular and Polymer Glass-Forming Liquids: Experimental and Theoretical Analysis under Isobaric Conditions.

We establish via an in-depth analysis of experimental data that the dimensionless compressibility (proportional to the dimensionless amplitude of long wavelength thermal density fluctuations) of one-component normal and supercooled liquids of chemically complex nonpolar and weakly polar molecules and polymers follows extremely well a surprisingly simple and general temperature dependence over an exceptionally wide range of pressures and temperatures. A theoretical basis for this behavior is shown to exist in the venerable van der Waals model and its more modern interpretations. Although associated hydrogen-bonding (and to a lesser degree strongly polar) liquids display modestly more complex behavior, rather simple temperature and pressure dependences are also discovered. A new approach to collapse the temperature- and pressure-dependent dimensionless compressibility data onto a master curve is formulated that differs from the empirical thermodynamic scaling approach. As a practical matter, we also find that the dimensionless compressibility scales well as an inverse power law with temperature with an exponent that is system dependent and decreases with pressure. At very high pressures and low temperatures, the thermal liquid behavior appears to approach (but not reach) a repulsion-dominated random close packing limit. All these findings are relevant to our recent theoretical work on the problem of activated relaxation and vitrification of supercooled molecular and polymeric liquids.

A Dynamically Correlated Network Model for the Collective Dynamics in Glass-Forming Molecular Liquids and Polymers.

The non-Arrhenius behavior of segmental dynamics in glass-forming liquids is one of the most profound mysteries in soft matter physics. In this article, we propose a dynamically correlated network (DCN) model to understand the growing behavior of dynamically correlated regions during cooling, which leads to the viscous slowdown of supercooled liquids. The fundamental concept of the model is that the cooperative region of collective motions has a network structure that consists of string-like parts, and networks of various sizes interpenetrate each other. Each segment undergoes dynamical coupling with its neighboring segments via a finite binding energy. Monte Carlo simulations showed that the fractal dimension of the DCNs generated at different temperatures increased and their size distribution became broader with decreasing temperature. The segmental relaxation time was evaluated based on a power law with four different exponents for the activation energy of rearrangement with respect to the DCN size. The results of the present DCN model are consistent with the experimental results for various materials of molecular and polymeric liquids.

Probing Dynamic Heterogeneity and Amorphous Order Using Rotational Dynamics of Rodlike Particles in Supercooled Liquids

Probing dynamic and static correlation in glass-forming supercooled liquids has been a challenge for decades in spite of extensive research. Dynamic correlation which manifests itself as Dynamic Heterogeneity is ubiquitous in a vast variety of systems starting from molecular glass-forming liquids, dense colloidal systems to collections of cells. On the other hand, the mere concept of static correlation in these dense disordered systems remain somewhat elusive and its existence is still actively debated. We propose a novel method to extract both dynamic and static correlations using rod-like particles as probe. This method can be implemented in molecular glass-forming liquids in experiments as well as in other soft matter systems including biologically relevant systems. We also rationalize the observed log-normal like distribution of rotational decorrelation time of elongated probe molecules in reported experimental studies along with a proposal of a novel methodology to extract dynamic and static correlation lengths in experiments.

Thermodynamic and Kinetic Transitions of Liquids in Nanoconfinement.

Core principles of chemistry are ubiquitously invoked to shed light on the nature of molecular level interactions in nanoconfined fluids, which play a pivotal role in a wide range of processes in geochemistry, biology, and engineering. A detailed understanding of the physicochemical processes involved in the flow, structural transitions, and freezing or melting behavior of fluids confined within nanometer-sized pores of solid materials is thus of enormous importance for both basic research and technological applications.This Account provides a perspective on new insights into the thermodynamic and kinetic transitions of nanoconfined fluids in their stable and metastable forms. After briefly introducing the unique properties of mesoporous silicas from the SBA, MCM, and FDU families that serve as the confinement matrices, combining highly ordered single and bimodal mesopore architectures with tunable pore sizes in the ∼2-15 nm range and narrow size distributions, recent studies on melting/freezing behavior of water confined in these host matrices are reviewed. While differential scanning calorimetry (DSC) reveals a linear relationship between melting point depression and pore size (independent of the pore shape), as predicted by the Gibbs-Thomson relation, variable temperature 2H wide-line nuclear magnetic resonance (NMR) spectroscopy studies confirm the core-shell model of water and give evidence for a layer-by-layer freezing mechanism, which gives rise to an apparent fragile-to-strong transition in the solidification dynamics.In contrast to the freezing/melting behavior of water, the effect of nanoconfinement on the glass transition of supercooled liquids is nonuniversal and the glass transition temperature Tg can either increase or decrease with the dimensionality and extent of confinement. This nonuniversal behavior is exemplified by the two glass-forming molecular liquids, glycerol and ortho-terphenyl (OTP). While glycerol shows an increase in Tg and a pronounced slowdown of the rotational dynamics of the constituent molecules due to a change in the molecular packing between the bulk and the confined liquid, OTP displays a linear and confining-media-dependent depression of Tg with increased confinement that is strongly influenced by the pore-liquid interface characteristics.This Account concludes with a focus on recent experimental evidence of extreme spatial and dynamical heterogeneity in both freezing and glass transition processes. This discovery was enabled by the unique mesoporous structures of SBA-16 and FDU-5, possessing bimodal architectures with two interconnected pore types of different size and shape (spherical and cylindrical). For the very first time, two melting points for water and two glass transitions for supercooled OTP, corresponding to a specific pore type, were observed. Collectively, these observations strongly suggest a close mechanistic connection between the local fluctuations in the structure and dynamics of nanoconfined liquids. While the findings reviewed in this Account provide new insights into thermodynamic and kinetic transitions of fluids, there remain many unanswered questions regarding the effects of nanoconfinement on the fundamental properties of fluids, which offer exciting future opportunities in chemical research.

Determination of Young's Modulus of Active Pharmaceutical Ingredients by Relaxation Dynamics at Elevated Pressures.

A new approach is theoretically proposed to study the glass transition of active pharmaceutical ingredients and a glass-forming anisotropic molecular liquid at high pressures. We describe amorphous materials as a fluid of hard spheres. Effects of nearest neighbor interactions and cooperative motions of particles on glassy dynamics are quantified through a local and collective elastic barrier calculated using the elastically collective nonlinear Langevin equation theory. Inserting two barriers into Kramer's theory gives the structural relaxation time. Then, we formulate a new mapping based on the thermal expansion process under pressure to intercorrelate particle density, temperature, and pressure. This analysis allows us to determine the pressure and temperature dependence of α relaxation. From this, we estimate the effective elastic modulus of amorphous materials and capture the effects of conformation on the relaxation process. Remarkably, our theoretical results agree well with experiments.

Effect of Polymers on Crystallization in Glass-Forming Molecular Liquids: Equal Suppression of Nucleation and Growth and Master Curve for Prediction

Crystal nucleation plays a critical role in the stability of supercooled liquids and glasses and is often controlled through addition of polymers. A dissolved polymer alters both the thermodynamics...

Can glassforming liquids be ‘simple’?

The tendency of various types of liquids to vitrify at moderate cooling rates is discussed. Good glass-forming liquids–covalent melts, polymer and organic liquids–feature high viscosity values of 10−2 – 107 Pa s at temperatures close to the melting point. The glasses obtained by cooling such liquids are nonergodic systems, i.e., their properties are not unambiguously defined by external parameters. At the same time, many glass-forming molecular liquids are usually considered to be ‘simple’ systems described by an effective pair central potential of interaction between particles. In particular, the scaling of thermodynamic and transport characteristics for varying temperature and density and behavior under pressure of the melting temperature and the bulk modulus of molecular liquids are well described by the parameters of a simple effective pair potential. It is shown that the values of the viscosity of glassy molecular liquids in no way correspond to the effective-potential parameters (in contrast to the viscosity of true ‘simple’ fluids, such as liquefied inert gases). Due to the complex structure of the molecules of such liquids, the effective-potential parameters heavily depend on the distance between particles. A complex set of intermolecular and interatomic chemical bonding emerges when such liquids are vitrified. Thus, glass-forming molecular liquids can only be considered ‘simple’ from the point of view of their thermodynamic properties, while the polyatomic structure of their molecules results in a ‘complex’ behavior of transport characteristics and anomalously high viscosity values.

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.

Microscopic origin of the logarithmic relaxation in molecular glass-forming liquids

Logarithmic relaxation is a unique relaxation process exhibited by a few molecular liquids and biomolecules. However, the microscopic origin of logarithmic relaxation is still unclear. To understand the origin of this process, we studied two liquids that exhibit logarithmic relaxation in a dissolved state using quasielastic neutron scattering (QENS) and depolarized dynamic light scattering (DDLS). Although the intermolecular potential of the liquids is drastically different in the dissolved state from the bulk liquids, we observed that the logarithmic relaxation still persists. Our results indicate that the intermolecular potential does not play a role in determining the logarithmic relaxation process. The coupling of rotational and translational relaxation processes could be the origin of the logarithmic relaxation process exhibited by the molecular liquids.

Crystal nucleation rates in glass-forming molecular liquids: D-sorbitol, D-arabitol, D-xylitol, and glycerol

The rate of crystal nucleation has been measured in four glass-forming molecular liquids: D-sorbitol, D-arabitol, D-xylitol, and glycerol. These polyalcohols have similar rates of crystal growth when compared at the same temperature relative to Tg (the glass transition temperature), peaking near 1.4 Tg, while the nucleation rates J are vastly different. In D-sorbitol and D-arabitol, J reaches a maximum of ∼108 m-3 s-1 near 1.1 Tg, whereas J < 10-2 m-3 s-1 in D-xylitol and <1 m-3 s-1 in glycerol based on no nucleation in large samples after long waits. This confirms the fundamentally different mechanisms for nucleation and growth. Near Tg, both nucleation and growth slow down with a similar temperature dependence, suggesting a similar kinetic barrier for the two processes. This temperature dependence is significantly weaker than that of viscosity η, approximately following η-0.75. This indicates that viscosity is a poor representative of the kinetic barrier for nucleation, and a better choice is the crystal growth rate. Under the latter assumption, the classical nucleation theory (CNT) describes our data reasonably well, yielding σ = 0.013 J/m2 for D-sorbitol and 0.026 J/m2 for D-arabitol, where σ is the critical nucleus/liquid interfacial free energy. There is no strong indication that the CNT fails as the length scale for corporative rearrangement exceeds the size of the critical nucleus, as recently suggested for lithium disilicate.

Influence of hydrogen bonding on the elastic properties of molecular glassforming liquids under high pressure

We present ultrasonic study of bulk modulus of liquid glycerol and propylene carbonate (PC) at room temperature under high pressure up to 0.6 and 1 GPa, respectively. Pressure dependence of the isothermal bulk modulus, recalculated from the ultrasonic adiabatic data, is almost linear for the both liquids with relatively similar slopes, in spite of the initial value of bulk modulus being twice higher for glycerol due to the presence of the hydrogen bonds. The values of pressure derivative (7.8 for glycerol and 8.6 for PC) indicate that the intermolecular interaction in the both liquids is close to the Lennard-Jones potential.

Influence of chemistry, interfacial width, and non-isothermal conditions on spatially heterogeneous activated relaxation and elasticity in glass-forming free standing films.

We employ the Elastically Collective Nonlinear Langevin Equation (ECNLE) theory of activated relaxation to study several questions in free standing thin films of glass-forming molecular and polymer liquids. The influence of non-universal chemical aspects on dynamical confinement effects is found to be relatively weak, but with the caveat that for the systems examined, the bulk ECNLE polymer theory does not predict widely varying fragilities. Allowing the film model to have a realistic vapor interfacial width significantly enhances the reduction of the film-averaged glass transition temperature, Tg, in a manner that depends on whether a dynamic or pseudo-thermodynamic averaging of the spatial mobility gradient is adopted. The nature of film thickness effects on the spatial profiles of the alpha relaxation time and elastic modulus is studied under non-isothermal conditions and contrasted with the corresponding isothermal behavior. Modest differences are found if a film-thickness dependent Tg is defined in a dynamical manner. However, adopting a pseudo-thermodynamic measure of Tg leads to a qualitatively new form of the alpha relaxation time gradient where highly mobile layers near the film surface coexist with strongly vitrified regions in the film interior. As a consequence, the film-averaged shear modulus can increase with decreasing film thickness, despite the Tg reduction and presence of a mobile surface layer. Such a behavior stands in qualitative contrast to the predicted mechanical softening under isothermal conditions. Spatial gradients of the elastic modulus are studied as a function of temperature, film thickness, probing frequency, and experimental protocol, and a rich behavior is found.

Comparison of high pressure and nanoscale confinement effects on crystallization of the molecular glass-forming liquid, dimethyl phthalate.

High pressure and nanoscopic confinement are two different strategies commonly employed to modify the physicochemical properties of various materials. Both strategies act mostly by changing the molecular packing. In this work, we performed a comparative study on the effect of compression and confined geometry on crystallization of a molecular liquid. Dielectric spectroscopy was employed to investigate the crystallization of the van der Waals liquid, dimethyl phthalate, in nanoporous alumina of different pore sizes as well as on increased pressure (up to 200 MPa). The analysis of the crystallization kinetics under varying thermodynamic conditions revealed that both strategies affect the crystallization behavior of the sample in very distinct ways. Compression shifts the maximum crystallization rate towards a higher temperature and broadens it. As a result, it is more challenging to avoid crystallization upon cooling the liquid at high pressure. In contrast, when the same material is incorporated into nanopores, crystallization significantly slows down and the maximum rate shifts towards a lower temperature with decreasing pore size. Finally, we show that crystallization in nanoporous alumina is accompanied by pre-crystallization effects upon which a shift of the α-relaxation peak is observed. An equilibration process prior to the initiation of crystallization was detected for the confined material both above and below the glass transition temperature of the interfacial layer, while not in the bulk.

Observation of coupling between vibrational and configurational entropy in a supercooled molecular chalcogenide glass-forming liquid

The temperature dependence of the line widths and frequencies of the vibrational modes of As4S3 molecules is studied in the Ge3As52S45 molecular glass and supercooled liquid at temperatures ranging from ~Tg−189K to up to ~Tg+111K, using Raman spectroscopy. The results reveal the presence of strong dynamical coupling between vibrational and configurational degrees of freedom of the molecules such that both the phonon lifetime and the density of states are affected by the rotational (and possibly translational) molecular dynamics in the supercooled liquid. The physical origin of such coupling is discussed within the framework of the potential energy landscape model of the glass transition.

Isotropic rotation vs. shear relaxation in supercooled liquids with globular cage molecules.

The temperature dependence of the rotational dynamics of P4Se3 molecules in the glass-forming molecular liquid P5Se3 is studied using two-dimensional (31)P nuclear magnetic resonance spectroscopy. Unlike typical molecular glass-forming liquids, the constituent molecules in the P5Se3 liquid perform rapid isotropic rotation without significant translational diffusion in the supercooled regime and this rotational process shows a decoupling in time scale from shear relaxation by nearly six orders of magnitude at the glass transition. This dynamical behavior of liquid-like rotation and localized translation appears to be universal to glass-forming liquids with high-symmetry globular molecules that are characterized by an underlying thermodynamically stable plastic crystal phase.

Secondary relaxation dynamics in rigid glass-forming molecular liquids with related structures.

The dielectric relaxation in three glass-forming molecular liquids, 1-methylindole (1MID), 5H-5-Methyl-6,7-dihydrocyclopentapyrazine (MDCP), and Quinaldine (QN) is studied focusing on the secondary relaxation and its relation to the structural α-relaxation. All three glass-formers are rigid and more or less planar molecules with related chemical structures but have dipoles of different strengths at different locations. A strong and fast secondary relaxation is detected in the dielectric spectra of 1MID, while no resolved β-relaxation is observed in MDCP and QN. If the observed secondary relaxation in 1MID is identified with the Johari-Goldstein (JG) β-relaxation, then apparently the relation between the α- and β-relaxation frequencies of 1MID is not in accord with the Coupling Model (CM). The possibility of the violation of the prediction in 1MID as due to either the formation of hydrogen-bond induced clusters or the involvement of intramolecular degree of freedom is ruled out. The violation is explained by the secondary relaxation originating from the in-plane rotation of the dipole located on the plane of the rigid molecule, contributing to dielectric loss at higher frequencies and more intense than the JG β-relaxation generated by the out-of-plane rotation. MDCP has smaller dipole moment located in the plane of the molecule; however, presence of the change of curvature of dielectric loss, ε″(f), at some frequency on the high-frequency flank of the α-relaxation reveals the JG β-relaxation in MDCP and which is in accord with the CM prediction. QN has as large an in-plane dipole moment as 1MID, and the absence of the resolved secondary relaxation is explained by the smaller coupling parameter than the latter in the framework of the CM.