Self-Diffusion of Supercooled Tris-naphthylbenzene

Importance of Experiments That Can Test Theories Critically

General dynamic and thermodynamic properties of complex materials, including amorphous polymers and molecular glass-formers, have been established from the wealth of experimental data accumulated over the years. Naturally, these general properties attract researchers to construct theories and models to address and explain them. Often more than one theory with contrasting or even conflicting theoretical bases can equally explain a general property rather well. The correct explanation becomes unclear, and progress is stopped. The resolution of the problem comes when an innovative experiment is performed with insightful results that can critically test the premise and assumptions of each theory. This important role played by experimentalists is exemplified by the contributions of Mark Ediger in several general properties considered in this paper: (1) dynamics of the components in binary polymer blends; (2) breakdown of the Stokes-Einstein and the Debye-Stokes-Einstein relations; (3) enhancement of surface mobility and in relation to formation of ultrastable glasses; and (4) the Johari-Goldstein β-relaxation in ultrastable glasses. Different theories proposed to explain these properties are discussed, including the Coupling Model of the author.

Anion dynamics and motional decoupling in a glycerol-choline chloride deep eutectic solvent studied by one- and two-dimensional 35Cl NMR

Chlorine-35 is among the few nuclides that provide an experimental handle on the anion dynamics in choline based deep eutectic solvents. By combining several nuclear magnetic resonance (NMR) techniques, the present work examines the Cl- motions within glyceline, a glycerol : choline chloride 2 : 1 solution, in a large temperature range down to the glass transition temperature Tg. The applied methods include spin relaxometry, second-order line shape analysis, as well as two-dimensional central-transition exchange and stimulated-echo spectroscopy. The finding of unstructured central-transition NMR spectra characterized by a relatively small average quadrupolar coupling attests to a highly disordered, essentially nondirectional anionic coordination in glyceline. For temperatures larger than about 1.3Tg the chlorine motions are well coupled to those of the glycerol and the choline moieties. At lower temperatures the local translational anion dynamics become Arrhenian and increasingly faster than the motion of glyceline's matrix molecules. Upon further cooling, the overall ionic conductivity continues to display a super-Arrhenius behavior, implying that the choline cations rather than the Cl anions dominate the long-range charge transport also near Tg.

Molecular Dynamics and Diffusion in Amorphous Solid Dispersions Containing Imidacloprid

The main goal of this study is to develop an experimental toolbox to estimate the self-diffusion coefficient of active ingredients (AI) in single-phase amorphous solid dispersions (ASD) close to the glass transition of the mixture using dielectric spectroscopy (DS) and oscillatory rheology. The proposed methodology is tested for a model system containing the insecticide imidacloprid (IMI) and the copolymer copovidone (PVP/VA) prepared via hot-melt extrusion. For this purpose, reorientational and the viscoelastic structural (α-)relaxation time constants of hot-melt-extruded ASDs were obtained via DS and shear rheology, respectively. These were then utilized to extract the viscosity as well as the fragility index of the dispersions as input parameters to the fractional Stokes-Einstein (F-SE) relation. Furthermore, a modified version of Almond-West (AW) formalism, originally developed to describe charge diffusion in ionic conductors, was exercised on the present model system for the estimation of the AI diffusion coefficients based on shear modulus relaxation times. Our results revealed that, at the calorimetric glass-transition temperature (T_g), the self-diffusion coefficients of the AI in the compositional range from infinite dilution up to 60 wt % IMI content lied in the narrow range of 10^-18-10^-20 m² s^-1, while the viscosity values of the dispersions at T_g varied between 10⁸ Pa s and 10¹⁰ Pa s. In addition, the phase diagram of the IMI-PVP/VA system was determined using the melting point depression method via differential scanning calorimetry (DSC), while mid-infrared (IR) spectroscopy was employed to investigate the intermolecular interactions within the solid dispersions. In this respect, the findings of a modest variation in melting point at different compositions stayed in agreement with the observations of weak hydrogen bonding interactions between the AI and the polymer. Moreover, IR spectroscopy showed the intermolecular IMI-IMI hydrogen bonding to have been considerably suppressed, as a result of the spatial separation of the AI molecules within the ASDs. In summary, this study provides experimental approaches to study diffusivity in ASDs using DS and oscillatory rheology, in addition to contributing to an enhanced understanding of the interactions and phase behavior in these systems.

Upper bound of fragility from spatial fluctuations of shear modulus and boson peak in glasses

It is shown that the normalized rms fluctuation of the shear modulus on the medium-range order scale in glasses correlates with fragility: the higher fragility, the smaller the fluctuation amplitude. The latter is calculated within the heterogeneous elasticity theory using the data on the boson peak in glasses. On a smaller scale corresponding to cooperative structural relaxation, the normalized rms fluctuation of the infinite-frequency shear modulus was estimated using the data on the decoupling of viscosity and diffusion in supercooled liquids. These fluctuations are much smaller in amplitude, and, in contrast, they increase with increasing fragility. Extrapolation predicts intersection of both rms fluctuations and disappearing of the boson peak at the upper limit to fragility ≈180.

Polymorphic selectivity in crystal nucleation

Crystal nucleation rates have been measured in the supercooled melts of two richly polymorphic glass-forming liquids: ROY and nifedipine (NIF). ROY or 5-methyl-2-[(2-nitrophenyl)amino]-3-thiophenecarbonitrile is known for its crystals of red, orange, and yellow colors and many polymorphs of solved structures (12). Of the many polymorphs, ON (orange needles) nucleates the fastest with the runner up (Y04) trailing by a factor of 103 when compared under the same mobility-limited condition, while the other unobserved polymorphs are slower yet by at least 5 orders of magnitude. Similarly, of the six polymorphs of NIF, [Formula: see text]′ nucleates the fastest, [Formula: see text]′ is slower by a factor of 10, and the rest are slower yet by at least 5 decades. In both systems, the faster-nucleating polymorphs are not built from the lowest-energy conformers, while they tend to have higher energies and lower densities and thus greater similarity to the liquid phase by these measures. The temperature ranges of this study covered the glass transition temperature Tg of each system, and we find no evidence that the nucleation rate is sensitive to the passage of Tg. At the lowest temperatures investigated, the rates of nucleation and growth are proportional to each other, indicating that a similar kinetic barrier controls both processes. The classical nucleation theory provides an accurate description of the observed nucleation rates if the crystal growth rate is used to describe the kinetic barrier for nucleation. The quantitative rates of both nucleation and growth for the competing polymorphs enable prediction of the overall rate of crystallization and its polymorphic outcome.

Surface diffusion of a glassy discotic organic semiconductor and the surface mobility gradient of molecular glasses

Surface diffusion has been measured in the glass of an organic semiconductor, MTDATA, using the method of surface grating decay. The decay rate was measured as a function of temperature and grating wavelength, and the results indicate that the decay mechanism is viscous flow at high temperatures and surface diffusion at low temperatures. Surface diffusion in MTDATA is enhanced by 4 orders of magnitude relative to bulk diffusion when compared at the glass transition temperature T_g. The result on MTDATA has been analyzed along with the results on other molecular glasses without extensive hydrogen bonds. In total, these systems cover a wide range of molecular geometries from rod-like to quasi-spherical to discotic and their surface diffusion coefficients vary by 9 orders of magnitude. We find that the variation is well explained by the existence of a steep surface mobility gradient and the anchoring of surface molecules at different depths. Quantitative analysis of these results supports a recently proposed double-exponential form for the mobility gradient: log D(T, z) = log D_v(T) + [log D₀ - log D_v(T)]exp(-z/ξ), where D(T, z) is the depth-dependent diffusion coefficient, D_v(T) is the bulk diffusion coefficient, D₀ ≈ 10^-8 m²/s, and ξ ≈ 1.5 nm. Assuming representative bulk diffusion coefficients for these fragile glass formers, the model reproduces the presently known surface diffusion rates within 0.6 decade. Our result provides a general way to predict the surface diffusion rates in molecular glasses.

Surface Diffusion Is Controlled by Bulk Fragility across All Glass Types

Surface diffusion is vastly faster than bulk diffusion in some glasses, but only moderately enhanced in others. We show that this variation is closely linked to bulk fragility, a common measure of how quickly dynamics is excited when a glass is heated to become a liquid. In fragile molecular glasses, surface diffusion can be a factor of 10^{8} faster than bulk diffusion at the glass transition temperature, while in the strong system SiO_{2}, the enhancement is a factor of 10. Between these two extremes lie systems of intermediate fragility, including metallic glasses and amorphous selenium and silicon. This indicates that stronger liquids have greater resistance to dynamic excitation from bulk to surface and enables prediction of surface diffusion, surface crystallization, and formation of stable glasses by vapor deposition.

Glass Dynamics Deep in the Energy Landscape

When a liquid is cooled, progress down the energy landscape is arrested near the glass transition temperature T_g. In principle, lower energy states can be accessed by waiting for further equilibration, but the rough energy landscape of glasses quickly leads to kinetics on geologically slow time scales below T_g. Over the past decade, progress has been made probing deeper into the energy landscape via several techniques. By looking at bulk and surface diffusion, using layered deposition that promotes equilibration, imaging glass surfaces with faster dynamics below T_g, and optically exciting glasses, experiments have moved into a regime of ultrastable, low energy glasses that was difficult to access in the past. At the same time, both simulations and energy landscape theory based on a random first order transition (RFOT) have tackled systems that include surfaces, optical excitation, and interfacial dynamics. Here we review some of the recent experimental work, and how energy landscape theory illuminates glassy dynamics well below the glass transition temperature by making direct connections between configurational entropy, energy landscape barriers, and the resulting dynamics.

Rheology based estimates of self- and collective diffusivities in viscous liquids

The self-diffusion coefficient of viscous liquids is estimated on the basis of a simple analysis of their rheological shear spectra. To this end, the Almond-West approach, previously employed to access single-particle diffusivities in ionic conductors, is generalized for application to molecular dynamics in supercooled liquids. Rheology based estimates, presented for indomethacin, ortho-terphenyl, and trinaphthylbenzene, reveal relatively small, yet systematic differences when compared with diffusivity data directly measured for these highly viscous liquids. These deviations are discussed in terms of mechanical Haven ratios, introduced to quantify the magnitude of collective translational effects that have an impact on the viscous flow.

A mechanism for ageing in a deeply supercooled molecular glass

Measurements of the decay of electric fields, formed spontaneously within vapour-deposited films of cis-methyl formate, provide the first direct assessment of the energy barrier to secondary relaxation in a molecular glass. At temperatures far below the glass transition temperature, the mechanism of relaxation is shown to be through hindered molecular rotation. Magnetically-polarised neutron scattering experiments exclude diffusion, which is demonstrated to take place only close to the glass transition temperature.

Diffusion-controlled and `diffusionless' crystal growth: relation between liquid dynamics and growth kinetics of griseofulvin

Understanding how liquid dynamics govern crystallization is critical for maintaining the physical stability of amorphous pharmaceutical formulations. In the present study, griseofulvin (GSF), a classic antifungal drug, was used as the model system to investigate the correlations between crystal growth kinetics and liquid dynamics. The temperature dependence of the kinetic part of the bulk crystal growth in a supercooled liquid of GSF was weaker than that of the structural relaxation time τα and scaled as τα −0.69. In the glassy state, GSF exhibited the glass-to-crystal (GC) growth behavior, whose growth rate was too fast to be under the control of the α-relaxation process. Moreover, from the perspective of τα, the GC growth of GSF also satisfied the general condition for GC growth to exist: D/u < 7 pm, where D is the diffusion coefficient and u the speed of crystal growth. Also compared were the fast surface crystal growth rates u s and surface relaxation times τsurface predicted by the random first-order transition theory. Here, the surface crystal growth rate u s of GSF exhibited a power-law dependence upon the surface structural relaxation time: u s ∝ τsurface −0.71, which was similar to that of the bulk growth rate and τα. These findings are important for understanding and predicting the crystallization of amorphous pharmaceutical solids both in the bulk and at the surface.

Physical Bond Breaking in Associating Copolymer Liquids

We combine ideas from polymer and glassy liquid physics to construct a new model for the bond-breaking time scale of attractive sticker groups in associating copolymer liquids that form transient networks. The activated event is argued to be a two-step process, involving first the release of the nonsticker dynamic caging constraints that defines the primary alpha relaxation, followed by attractive stickers surmounting an association free-energy barrier subject to a local frictional resistance which can be strongly affected by relaxation-diffusion decoupling. The ideas embedded in the model produce a consistent and good description of the bond-breaking time scale for diverse polymer chemistries and architectures as a function of temperature and fraction of sticky groups. Chemically sensible values for association free energies are deduced. In strong contrast, the existing phenomenological models are shown to incur qualitative failures.

Surface diffusion in glasses of rod-like molecules posaconazole and itraconazole: effect of interfacial molecular alignment and bulk penetration

The method of surface grating decay has been used to measure surface diffusion in the glasses of two rod-like molecules posaconazole (POS) and itraconazole (ITZ). Although structurally similar antifungal medicines, ITZ forms liquid-crystalline phases while POS does not. Surface diffusion in these systems is significantly slower than in the glasses of quasi-spherical molecules of similar volume when compared at the glass transition temperature T_g. Between the two systems, ITZ has slower surface diffusion. These results are explained on the basis of the near-vertical orientation of the rod-like molecules at the surface and their deep penetration into the bulk where mobility is low. For molecular glasses without extensive hydrogen bonds, we find that the surface diffusion coefficient at T_g decreases smoothly with the penetration depth of surface molecules and the trend has the double-exponential form for the surface mobility gradient observed in simulations. This supports the view that these molecular glasses have a similar mobility vs. depth profile and their different surface diffusion rates arise simply from the different depths at which molecules are anchored. Our results also provide support for a previously observed correlation between the rate of surface diffusion and the fragility of the bulk liquid.

Pressure densified 1,3,5-tri(1-naphthyl)benzene glass. I. Volume recovery and physical aging

The effects of pressure densification on 1,3,5-tri(1-naphthyl)benzene (TNB) are assessed from volumetric and calorimetric measurements. The pressure densified glass (PDG) has higher density than conventional glass (CG), but unlike ultrastable TNB glass prepared using vapor deposition which also has elevated density, TNB PDG exhibits higher enthalpy and lower thermal stability than when formed at ambient pressure. PDG also exhibits anomalous physical aging. Rather than evolving monotonically toward the equilibrium density, there is an overshoot to a lower density state. Only when the density of the PDG becomes equivalent to the corresponding CG does the density begin a slow approach toward equilibrium.

Does a Growing Static Length Scale Control the Glass Transition?

Several theories of the glass transition propose that the structural relaxation time τ_{α} is controlled by a growing static length scale ξ that is determined by the free energy landscape but not by the local dynamic rules governing its exploration. We argue, based on recent simulations using particle-radius-swap dynamics, that only a modest factor in the increase in τ_{α} on approach to the glass transition may stem from the growth of a static length, with a vastly larger contribution attributable, instead, to a slowdown of local dynamics. This reinforces arguments that we base on the observed strong coupling of particle diffusion and density fluctuations in real glasses.

Study of the upper-critical dimension of the East model through the breakdown of the Stokes-Einstein relation

We investigate the dimensional dependence of dynamical fluctuations related to dynamic heterogeneity in supercooled liquid systems using kinetically constrained models. The d-dimensional spin-facilitated East model with embedded probe particles is used as a representative super-Arrhenius glass forming system. We examine the existence of an upper critical dimension in this model by considering decoupling of transport rates through an effective fractional Stokes-Einstein relation, D∼τ^-1+ω, with D and τ the diffusion constant of the probe particle and the relaxation time of the model liquid, respectively, and where ω>0 encodes the breakdown of the standard Stokes-Einstein relation. To the extent that decoupling indicates non-mean-field behavior, our simulations suggest that the East model has an upper critical dimension at least above d = 10 and argue that it may actually be infinite. This result is due to the existence of hierarchical dynamics in the East model in any finite dimension. We discuss the relevance of these results for studies of decoupling in high dimensional atomistic models.

Influence of Hydrogen Bonding on the Surface Diffusion of Molecular Glasses: Comparison of Three Triazines

Surface grating decay measurements have been performed on three closely related molecular glasses to study the effect of intermolecular hydrogen bonds on surface diffusion. The three molecules are derivatives of bis(3,5-dimethyl-phenylamino)-1,3,5-triazine and differ only in the functional group R at the 2-position, with R being C₂H₅, OCH₃, and NHCH₃, and referred to as "Et", "OMe", and "NHMe", respectively. Of the three molecules, NHMe forms more extensive intermolecular hydrogen bonds than Et and OMe and was found to have slower surface diffusion. For Et and OMe, surface diffusion is so fast that it replaces viscous flow as the mechanism of surface grating decay as temperature is lowered. In contrast, no such transition was observed for NHMe under the same conditions, indicating significantly slower surface diffusion. This result is consistent with the previous finding that extensive intermolecular hydrogen bonds slow down surface diffusion in molecular glasses and is attributed to the persistence of hydrogen bonds even in the surface environment. This result is also consistent with the lower stability of the vapor-deposited glass of NHMe relative to those of Et and OMe and supports the view that surface mobility controls the stability of vapor-deposited glasses.

Growth rate of crystalline ice and the diffusivity of supercooled water from 126 to 262 K

Understanding deeply supercooled water is key to unraveling many of water's anomalous properties. However, developing this understanding has proven difficult due to rapid and uncontrolled crystallization. Using a pulsed-laser-heating technique, we measure the growth rate of crystalline ice, G(T), for 180 K < T < 262 K, that is, deep within water's "no man's land" in ultrahigh-vacuum conditions. Isothermal measurements of G(T) are also made for 126 K ≤ T ≤ 151 K. The self-diffusion of supercooled liquid water, D(T), is obtained from G(T) using the Wilson-Frenkel model of crystal growth. For T > 237 K and P ∼ 10-8 Pa, G(T) and D(T) have super-Arrhenius ("fragile") temperature dependences, but both cross over to Arrhenius ("strong") behavior with a large activation energy in no man's land. The fact that G(T) and D(T) are smoothly varying rules out the hypothesis that liquid water's properties have a singularity at or near 228 K at ambient pressures. However, the results are consistent with a previous prediction for D(T) that assumed no thermodynamic transitions occur in no man's land.

Dielectric and specific heat relaxations in vapor deposited glycerol

Recently [S. Capponi, S. Napolitano, and M. Wübbenhorst, Nat. Commun. 3, 1233 (2012)], vapor deposited glasses of glycerol have been found to recover their super-cooled liquid state via a metastable, ordered liquid (MROL) state characterized by a tremendously enhanced dielectric strength along with a slow-down of the relaxation rate of the structural relaxation. To study the calorimetric signature of this phenomenon, we have implemented a chip-based, differential AC calorimeter in an organic molecular beam deposition setup, which allows the simultaneous measurement of dielectric relaxations via interdigitated comb electrodes and specific heat relaxation spectra during deposition and as function of the temperature. Heating of the as-deposited glass just above the bulk Tg and subsequent cooling/reheating revealed a step-wise increase in cp by in total 9%, indicating unambiguously that glycerol, through slow vapour deposition, forms a thermodynamically stable glass, which has a specific heat as low as that of crystalline glycerol. Moreover, these glasses were found to show excellent kinetic stability as well as evidenced by both a high onset-temperature and quasi-isothermal recovery measurements at -75 °C. The second goal of the study was to elucidate the impact of the MROL state on the specific heat and its relaxation to the super-cooled state. Conversion of "MROL glycerol" to its "normal" (ordinary liquid, OL) state revealed a second, small (∼2%) increase of the glassy cp, a little gain (<10%) in the relaxed specific heat, and no signs of deviations of τcal from that of normal "bulk" glycerol. These findings altogether suggest that the MROL state in glycerol comprises largely bulk-type glycerol that coexist with a minor volume fraction (<10%) of PVD-induced structural anomalies with a crystal-like calorimetric signature. Based on the new calorimetric findings, we have proposed a new physical picture that assumes the existence of rigid polar clusters (RPCs) and conclusively explains the extraordinary high kinetic stability of the MROL state, its specific calorimetric signature, the enhanced strength, and apparent slow-down of the dielectric α-relaxation. In this new picture, the incredibly slow and strengthened dielectric response is ascribed to driven rotational diffusion of whole RPCs, a mechanism that perfectly couples to the relaxation time of the "normal" glycerol fraction. First considerations based on the strength and the retardation of the dielectric RPCs' response yield independently a size estimate for the RPCs in the order of 4-5 nm. Finally, we have discussed possible crystallisation and reorganisation effects, which give rise to pronounced out-of phase components of the specific heat at higher temperatures.

Enhanced diffusion and mobile fronts in a simple lattice model of glass-forming liquids

The diffusion of mobility in bulk and thin film fluids near their glass transition is examined with a kinetic lattice model, and compared to recent experiments on bulk liquids and vapor-deposited thin film glasses. The "limited mobility" (LM) lattice model exhibits dynamic heterogeneity of mobility when the fluid is near its kinetic arrest transition; a finite-parameter second-order critical point in the LM model bearing strong resemblance to the glass transition in real fluids. The spatial heterogeneity of mobility near kinetic arrest leads to dynamics that violate the Stokes-Einstein relation. To make connections with experiment, LM model simulations of self-diffusion constants in fluids near kinetic arrest are compared to those in two organic glass-formers. In addition, simulations of mobility in films that have been temperature-jumped above kinetic arrest (starting from an arrested state) are carried out. The films develop a "front" of mobility at their free surface that progresses into the film interior at a constant rate, thereby mobilising the entire film to fluidity. The velocity of the front scales with the self-diffusion constant for analogous bulk systems-an observation consistent with experiments on vapor-deposited molecular thin films.

Synthesis and high-throughput characterization of structural analogues of molecular glassformers: 1,3,5-trisarylbenzenes

We report the synthesis and characterization of an analogous series of small organic molecules derived from a well-known glass former, 1,3-bis(1-naphthyl)-5-(2-naphthyl)benzene (α,α,β-TNB). Synthesized molecules include α,α,β-TNB, 3,5-di(naphthalen-1-yl)-1-phenylbenzene (α,α-P), 9-(3,5-di(naphthalen-1-yl)phenyl)anthracene (α,α-A), 9,9'-(5-(naphthalen-2-yl)-1,3-phenylene)dianthracene (β-AA) and 3,3',5,5'-tetra(naphthalen-1-yl)-1,1'-biphenyl (α,α,α,α-TNBP). The design of molecules was based on increasing molecular weight with varied π-π interactions in one or more substituents. The synthesis is based on Suzuki cross-coupling of 1-bromo-3-chloro-5-iodobenzene with arylboronic acids, which allows attachment of various substituents to tailor the chemical structure. The bulk compounds were characterized using NMR spectroscopy and differential scanning calorimetry (DSC). Thin films of these compounds were produced using physical vapor deposition and were subsequently annealed above the glass transition temperatures (Tg). For each molecular glass, cooling rate-dependent glass transition temperature measurements (CR-Tg) were performed using ellipsometry as a high-throughput method to characterize thin film properties. CR-Tg allows rapid characterization of glassy properties, such as Tg, apparent thermal expansion coefficients, apparent activation energy at Tg and fragility. DSC measurements confirmed the general trend that increasing molecular weight leads to increasing melting point (Tm) and Tg. Furthermore, CR-Tg provided evidence that the introduction of stronger π-interacting substituents in the chosen set of structural analogues increases fragility and decreases the ability to form glasses, such that β-AA has the largest fragility and highest tendency to crystallize among all the compounds. These strong interactions also significantly elevate Tg and promote more harmonic intermolecular potentials, as observed by decreasing value of the apparent thermal expansion coefficient.

Fast surface diffusion of amorphous o-terphenyl and its competition with viscous flow in surface evolution

Surface self-diffusion coefficients have been measured for the model molecular glass o-terphenyl (OTP) through surface-grating decay driven by capillarity. The decay mechanism transitions from viscous flow at high temperatures to surface diffusion at low temperatures; for 1000 nm wavelength gratings, the transition occurs at Tg + 11 K. The surface diffusion of OTP is 10(8) times faster than bulk diffusion at Tg and even faster at lower temperatures because of its weaker temperature dependence. At Tg, OTP has approximately the same bulk diffusivity as the previously studied molecular liquid indomethacin, but its surface diffusion is 100 times faster. While the molecular glass-formers exhibit transitions from viscous flow to surface diffusion as the mechanism of capillarity-driven surface flattening, polystyrenes and silicates show no such transition under comparable conditions, suggesting slower surface diffusion on these materials and a general dependence of surface diffusion on intermolecular forces. The velocity of surface crystal growth on molecular glasses is proportional to surface diffusivity, indicating a common kinetic barrier for both processes for temperatures below Tg.

Termination of Solid-State Crystal Growth in Molecular Glasses by Fluidity

Fast crystal growth can abruptly emerge as molecular liquids are cooled to become glasses, a phenomenon presently unknown for other materials. This glass-to-crystal (GC) mode can cause crystallization rates orders of magnitude faster than those predicted by standard models. While GC growth is known for 12 systems, its features vary greatly with growth rates spanning a factor of 10(4). We report that the general condition for GC growth to exist is that liquid diffusion be slow relative to crystal growth according to D/u < 7 pm. This condition holds for all liquids exhibiting GC growth and suggests that the phenomenon is a solid-state process terminated by fluidity. GC growth must solidify several molecular layers before rearrangement by diffusion. We propose that GC growth propagates by a nonequilibrium crystal/liquid interface 3 nm wide that solidifies by local mobility. These results explain the prevalence of GC growth among organic liquids and guide its discovery in other materials.

Structural relaxation of vapor-deposited molecular glasses and supercooled liquids

The properties of vapor-deposited molecular glasses largely depend on deposition conditions, and stable and/or dense glasses are formed with several compounds.

Simple picture of supercooled liquid dynamics: dynamic scaling and phenomenology based on clusters

Although it is now well established that in glassy liquids, slow structural relaxation accompanies a correlated structural rearrangement, the role of such a correlation in the transport anomaly, and thus in the slow dynamics, remains unclear. In this paper, we argue from a hydrodynamic viewpoint that a correlated structure (cluster) with a characteristic size ξ sustains the long-lived stress and dynamically couples with the hydrodynamic fluctuations; therefore, the dynamics of this cluster is the origin of the mesoscopic nature of anomalous hydrodynamic transport. Based on this argument, we derive a dynamic scaling law for τ(α) (or η, where η is the macroscopic shear viscosity) as a function of ξ: τ(α)([proportionality]η)[proportionality]ξ(4). We provide a simple explanation for basic features of anomalous transport, such as the breakdown of the Stokes-Einstein relation and the length-scale-dependent decoupling between viscosity and diffusion. The present study further suggests a different physical picture: Through the coarse graining of smaller-scale fluctuations (</~ξ), the supercooled liquid dynamics can be regarded as the dynamics of normal (cluster) liquids composed of units with a typical size of ξ. Although the correlation length of hydrodynamic transport ξ and the dynamic heterogeneity size ξ(DH), which is determined by the usual four-point correlation function, reflect some aspects of the cooperative effects, the correspondence between ξ and ξ(DH) is not one to one. We highlight the possibility that ξ(DH) overestimates the actual collective transport range at a low degree of supercooling.

Coherent neutron scattering and collective dynamics on mesoscale

By combining, and modestly extending, a variety of theoretical concepts for the dynamics of liquids in the supercooled regime, we formulate a simple analytic model for the temperature and wavevector dependent collective density fluctuation relaxation time that is measurable using coherent dynamic neutron scattering. Comparison with experiments on the ionic glass-forming liquid Ca-K-NO3 in the lightly supercooled regime suggests the model captures the key physics in both the local cage and mesoscopic regimes, including the unusual wavevector dependence of the collective structural relaxation time. The model is consistent with the idea that the decoupling between diffusion and viscosity is reflected in a different temperature dependence of the collective relaxation time at intermediate wavevectors and near the main (cage) peak of the static structure factor. More generally, our analysis provides support for the ideas that decoupling information and growing dynamic length scales can be at least qualitatively deduced by analyzing the collective relaxation time as a function of temperature and wavevector, and that there is a strong link between dynamic heterogeneity phenomena at the single and many particle level. Though very simple, the model can be applied to other systems, such as molecular liquids.

An explanation of the differences in diffusivity of the components of the metallic glass Pd43Cu27Ni10P20

Bartsch et al. [Phys. Rev. Lett. 104, 195901 (2010)] reported measurements of the diffusivities of different components of the multi-component bulk metallic glass Pd43Cu27Ni10P20. The diffusion of the largest Pd and the smallest P was found to be drastically different. The Stokes-Einstein relation breaks down when considering the P constituent atom, while the relation is obeyed by the Pd atom over 14 orders of magnitude of change in Pd diffusivity. This difference in behavior of Pd and P poses a problem challenging for explanation. With the assist of a recent finding in metallic glasses that the β-relaxation and the diffusion of the smallest component are closely related processes by Yu et al. [Phys. Rev. Lett. 109, 095508 (2012)], we use the Coupling Model to explain the observed difference between P and Pd quantitatively. The same model also explains the correlation between property of the β-relaxation with fragility found in the family of (CexLa1-x)68Al10Cu20Co2 with 0 ≤ x ≤ 1.

Molecular dynamics simulation of amorphous indomethacin

Molecular dynamics (MD) simulations have been conducted using an assembly consisting of 105 indomethacin (IMC) molecules and 12 water molecules to investigate the underlying dynamic (e.g., rotational and translational diffusivities and conformation relaxation rates) and structural properties (e.g., conformation, hydrogen-bonding distributions, and interactions of water with IMC) of amorphous IMC. These properties may be important in predicting physical stability of this metastable material. The IMC model was constructed using X-ray diffraction data with the force-field parameters mostly assigned by analogy with similar groups in Amber-ff03 and atomic charges calculated with the B3LYP/ccpVTZ30, IEFPCM, and RESP models. The assemblies were initially equilibrated in their molten state and cooled through the glass transition temperature to form amorphous solids. Constant temperature dynamic runs were then carried out above and below the T(g) (i.e., at 600 K (10 ns), 400 K (350 ns), and 298 K (240 ns)). The density (1.312 ± 0.003 g/cm(3)) of the simulated amorphous solid at 298 K was close to the experimental value (1.32 g/cm(3)) while the estimated T(g) (384 K) was ~64 degrees higher than the experimental value (320 K) due to the faster cooling rate. Due to the hindered rotation of its amide bond, IMC can exist in different diastereomeric states. Different IMC conformations were sufficiently sampled in the IMC melt or vapor, but transitions occurred rarely in the glass. The hydrogen-bonding patterns in amorphous IMC are more complex in the amorphous state than in the crystalline polymorphs. Carboxylic dimers that are dominant in α- and γ-crystals were found to occur at a much lower probability in the simulated IMC glasses while hydrogen-bonded IMC chains were more easily identified patterns in the simulated amorphous solids. To determine molecular diffusivity, a novel analytical method is proposed to deal with the non-Einsteinian behavior, in which the temporal evolution of the apparent diffusivity D(t) is described by a relaxation model such as the KWW function and extrapolated to infinite time. The diffusion coefficient found for water diffusing in amorphous indomethacin at 298 K (2.7 × 10(-9) cm(2)/s) compares favorably to results obtained in experimental IMC glasses (0.9-2.0 × 10(-9) cm(2)/s) and is mechanistically associated with β-relaxation processes that are dominant in sub-T(g) glasses.

Dynamic scaling for anomalous transport in supercooled liquids

The anomalous mesoscopic transport in supercooled liquids was investigated using three-dimensional molecular dynamics simulation. We show that the dynamic correlation length, ξ, can be identified as a crossover length between the microscopic and macroscopic transports. We also find that in highly supercooled states, where a clear stress plateau is observed in the stress autocorrelation, cooperative transport, in both the (longitudinal) density diffusion and the (transverse) viscous relaxation, can be characterized by the single length scale, ξ. The present identification of the dynamic correlation length has an advantage over other characterization methods in that it directly interprets the anomalous hydrodynamic transport in terms of the growing length scale. In the context of the present study, we provide a simple explanation for the long-standing problem of the breakdown of the Stokes-Einstein relation.

Translation-rotation decoupling and nonexponentiality in room temperature ionic liquids

Using a combination of light scattering techniques and broadband dielectric spectroscopy, we have measured the temperature dependence of structural relaxation time and self diffusion in three imidazolium-based room temperature ionic liquids: [bmim][NTf(2)], [bmim][PF(6)], and [bmim][TFA]. A detailed analysis of the results demonstrates that self diffusion decouples from structural relaxation in these systems as the temperature is decreased toward T(g). The degree to which the dynamics are decoupled, however, is shown to be surprisingly weak when compared to other supercooled liquids of similar fragility. In addition to the weak decoupling, we demonstrate that the temperature dependence of the structural relaxation time in all three liquids can be well described by a single Vogel-Fulcher-Tamann function over 13 decades in time from 10(-11) s up to 10(2) s. Furthermore, the stretching of the structural relaxation is shown to be temperature independent over the same range of time scales, i.e., time temperature superposition is valid for these ionic liquids from far above the melting point down to the glass transition temperature. We suggest that these phenomena are interconnected and all result from the same underlying mechanism--strong and directional intermolecular interactions.

On the dynamics of liquids in their viscous regime approaching the glass transition

Recently, Mallamace et al. (Eur. Phys. J. E 34, 94 (2011)) proposed a crossover temperature, T(×), and claimed that the dynamics of many supercooled liquids follow an Arrhenius-type temperature dependence between T(×) and the glass transition temperature T(g). The opposite, namely super-Arrhenius behavior in this viscous regime, has been demonstrated repeatedly for molecular glass-former, for polymers, and for the majority of the exhaustively studied inorganic glasses of technological interest. Therefore, we subject the molecular systems of the Mallamace et al. study to a "residuals" analysis and include not only viscosity data but also the more precise data available from dielectric relaxation experiments over the same temperature range. Although many viscosity data sets are inconclusive due to their noise level, we find that Arrhenius behavior is not a general feature of viscosity in the T(g) to T(×) range. Moreover, the residuals of dielectric relaxation times with respect to an Arrhenius law clearly reveal systematic curvature consistent with super-Arrhenius behavior being an endemic feature of transport properties in this viscous regime. We also observe a common pattern of how dielectric relaxation times decouple slightly from viscosity.

Space-time correlated two-particle hopping in glassy fluids: structural relaxation, irreversibility, decoupling, and facilitation

The microscopic nonlinear Langevin equation (NLE) theory of correlated two-particle dynamics in dense fluids of spherical particles is extended to construct a predictive model of multiple correlated hopping and recaging events of a pair of tagged particles as a function of their initial separation. Modest coarse graining over the liquid structural disorder allows contact to be made with various definitions of irreversible particle motion within the context of a multistate Markov model. The correlated space-time hopping process that underlies structural relaxation can also be analyzed in the context of kinetically constrained models. The dependence of microscopically defined mean persistence and exchange times, their distributions, and relaxation-diffusion decoupling on hard-sphere fluid volume fraction is derived from a model in which irreversible jumps serve as the nucleating persistence event. For a subset of questions, the predictions of the two-particle theory are compared with results from the earlier single-particle NLE approach.

Shear-transformation-zone theory of viscosity, diffusion, and stretched exponential relaxation in amorphous solids

The shear-transformation-zone (STZ) theory has been remarkably successful in accounting for broadly peaked, frequency-dependent, viscoelastic responses of amorphous systems near their glass temperatures T_{g}. This success is based on the theory's first-principles prediction of a wide range of internal STZ transition rates. Here, I show that the STZ rate distribution causes the Newtonian viscosity to be strongly temperature dependent; and I propose that it is this temperature dependence, rather than any heterogeneity-induced enhancement of diffusion, that is at least in part responsible for Stokes-Einstein violations near T_{g}. I also show that stretched-exponential relaxation of density fluctuations emerges naturally from the same distribution of STZ transition rates that predicts the viscoelastic behavior. To be consistent with observations of Fickian diffusion near T_{g}, however, an STZ-based diffusion theory somehow must include the cascades of correlated displacement events that are seen in low-temperature numerical simulations.

Breaking Through the Glass Ceiling: Recent Experimental Approaches to Probe the Properties of Supercooled Liquids near the Glass Transition

Experimental measurements of the properties of supercooled liquids at temperatures near their glass transition temperatures, Tg, are requisite for understanding the behavior of glasses and amorphous solids. Unfortunately, many supercooled molecular liquids rapidly crystallize at temperatures far above their Tg, making such measurements difficult to nearly impossible. In this Perspective, we discuss some recent alternative approaches to obtain experimental data in the temperature regime near Tg. These new approaches may yield the additional experimental data necessary to test current theoretical models of the dynamical slowdown that occurs in supercooled liquids approaching the glass transition.