Universality of the onset of activated transport in Lennard-Jones liquids with tunable coordination: Implications for the effects of pressure and directional bonding on the crossover to activated transport, configurational entropy, and fragility of glassforming liquids

Temperature-driven narrowing of the insulating gap as a precursor of the insulator-to-metal transition: Implications for the electronic structure of solids

We present a microscopic picture rationalizing the surprisingly steep decrease in the bandgap with temperature in insulators, crystalline or otherwise. The gap narrowing largely results from fluctuations of long-wavelength optical phonons-when the latter are present-or their disordered analogs if the material is amorphous. We elaborate on this notion to show that possibly with the exception of weakly bound solids made of closed-shell electronic configurations, the existence of an insulating gap or pseudogap in a periodic solid implies that optical phonons must be present, too. This means that in an insulating solid, the primitive cell must have at least two atoms and/or that a charge density wave is present, with the possible exception of weakly bonded solids such as rare-gas or ferromagnetic Wigner crystals. As a corollary, a (periodic) elemental solid held together by nonclosed shell interactions and whose primitive unit contains only one atom will ordinarily be a metal, consistent with observation. Consequences of the present picture for Wigner solids are discussed. A simple field theory of the metal-insulator transition is constructed that directly ties long-wavelength optical vibrations with fluctuations of an order parameter for the metal-insulator transition. The order parameter is shown to have at least two components, yet no Goldstone mode arises as a result of the transition.

Classical density functional theory, unconstrained crystallization, and polymorphic behavior

While in principle, classical density functional theory (cDFT) should be a powerful tool for the study of crystallization, in practice this has not so far been the case. Progress has been hampered by technical problems which have plagued the study of the crystalline systems using the most sophisticated fundamental measure theory models. In this paper, the reasons for the difficulties are examined and it is proposed that the tensor functionals currently favored are in fact numerically unstable. By reverting to an older, more heuristic model it is shown that all of the technical difficulties are eliminated. Application to a Lennard-Jones fluid results in a demonstration of power of cDFT to describe crystallization in a highly inhomogeneous system. First, we show that droplets attached to a slightly hydrophobic wall crystallize spontaneously upon being quenched. The resulting crystallites are clearly faceted structures and are predominantly HCP structures. In contrast, droplets in a fully periodic calculational cell remain stable to lower temperatures and eventually show the same spontaneous localization of the density into "atoms" but in an amorphous structure having many of the structural characteristics of a glass. A small change of the protocol leads, at the same temperature, to the formation of crystals, this time with the fcc structure typical of bulk Lennard-Jones solids. The fcc crystals have lower free energy than the amorphous structures which in turn are more stable than the liquid droplets. It is demonstrated that as the temperature is raised, the free energy differences between the structures decrease until the solid clusters become less stable than the liquid droplets and spontaneously melt. The presence of energy barriers separating the various structures is therefore clearly demonstrated.

Aging, Jamming, and the Limits of Stability of Amorphous Solids

Apart from not having crystallized, supercooled liquids can be considered as being properly equilibrated and thus can be described by a few thermodynamic control variables. In contrast, glasses and other amorphous solids can be arbitrarily far away from equilibrium and require a description of the history of the conditions under which they formed. In this paper we describe how the locality of interactions intrinsic to finite-dimensional systems affects the stability of amorphous solids far off equilibrium. Our analysis encompasses both structural glasses formed by cooling and colloidal assemblies formed by compression. A diagram outlining regions of marginal stability can be adduced which bears some resemblance to the quasi-equilibrium replica meanfield theory phase diagram of hard sphere glasses in high dimensions but is distinct from that construct in that the diagram describes not true phase transitions but kinetic transitions that depend on the preparation protocol. The diagram exhibits two distinct sectors. One sector corresponds to amorphous states with relatively open structures, the other to high density, more closely packed ones. The former transform rapidly owing to there being motions with no free energy barriers; these motions are string-like locally. In the dense region, amorphous systems age via compact activated reconfigurations. The two regimes correspond, in equilibrium, to the collisional or uniform liquid and the so-called landscape regime, respectively. These are separated by a spinodal line of dynamical crossovers. Owing to the rigidity of the surrounding matrix in the landscape, high-density part of the diagram, a sufficiently rapid pressure quench adds compressive energy which also leads to an instability toward string-like motions with near vanishing barriers. Conversely, a dilute collection of rigid particles, such as a colloidal suspension leads, when compressed, to a spatially heterogeneous structure with percolated mechanically stable regions. This jamming corresponds to the onset of activation when the spinodal line is traversed from the low density side. We argue that a stable glass made of sufficiently rigid particles can also be viewed as exhibiting sporadic and localized buckling instabilities that result in local jammed structures. The lines of instability we discuss resemble the Gardner transition of meanfield systems but, in contrast, do not result in true criticality owing to being short-circuited by activated events. The locally marginally stable modes of motion in amorphous solids correspond to secondary relaxation processes in structural glasses. Their relevance to the low temperature anomalies in glasses is also discussed.

Dependence of the configurational entropy on amorphous structures of a hard-sphere fluid

The free energy of a hard-sphere fluid for which the average energy is trivial signifies how its entropy changes with packing. The packing η_{f} at which the free energy of the crystalline state becomes lower than that of the disordered fluid state marks the freezing point. For packing fractions η>η_{f} of the hard-sphere fluid, we use the modified weighted density functional approximation to identify metastable free energy minima intermediate between uniform fluid and crystalline states. The distribution of the sharply localized density profiles, i.e., the inhomogeneous density field ρ(x) characterizing the metastable state is primarily described by a pair function g_{s}(η/η_{0}). η_{0} is a structural parameter such that for η=η_{0} the pair function is identical to that for the Bernal random structure. The configurational entropy S_{c} of the metastable hard-sphere fluid is calculated by subtracting the corresponding vibrational entropy from the total entropy. The extrapolated S_{c} vanishes as η→η_{K} and η_{K} is in agreement with other works. The dependence of η_{K} on the structural parameter η_{0} is obtained.

On the mechanism of activated transport in glassy liquids

We explore several potential issues that have been raised over the years regarding the "entropic droplet" scenario of activated transport in liquids, due to Wolynes and co-workers, with the aim of clarifying the status of various approximations of the random first-order transition theory (RFOT) of the structural glass transition. In doing so, we estimate the mismatch penalty between alternative aperiodic structures, above the glass transition; the penalty is equal to the typical magnitude of free energy fluctuations in the liquid. The resulting expressions for the activation barrier and the cooperativity length contain exclusively bulk, static properties; in their simplest form they contain only the bulk modulus and the configurational entropy per unit volume. The expressions are universal in that they do not depend explicitly on the molecular detail. The predicted values for the barrier and cooperativity length and, in particular, the temperature dependence of the barrier are in satisfactory agreement with observation. We thus confirm that the entropic droplet picture is indeed not only internally consistent but is also fully constructive, consistent with the apparent success of its many quantitative predictions. A simple view of a glassy liquid as a locally metastable, degenerate pattern of frozen-in stress emerges in the present description. Finally, we derive testable relationships between the bulk modulus and several characteristics of glassy liquids and peculiarities in low-temperature glasses.

Microscopically based calculations of the free energy barrier and dynamic length scale in supercooled liquids: the comparative role of configurational entropy and elasticity

We compute the temperature-dependent barrier for α-relaxations in several liquids, without adjustable parameters, using experimentally determined elastic, structural, and calorimetric data. We employ the random first order transition (RFOT) theory, in which relaxation occurs via activated reconfigurations between distinct, aperiodic minima of the free energy. Two different approximations for the mismatch penalty between the distinct aperiodic states are compared, one due to Xia and Wolynes (Proc. Natl. Acad. Sci. U. S. A. 2000, 97, 2990), which scales universally with temperature as for hard spheres, and one due to Rabochiy and Lubchenko (J. Chem. Phys. 2013, 138, 12A534), which employs measured elastic and structural data for individual substances. The agreement between the predictions and experiment is satisfactory, given the uncertainty in the measured experimental inputs. The explicitly computed barriers are used to calculate the glass transition temperature for each substance--a kinetic quantity--from the static input data alone. The temperature dependence of both the elastic and structural constants enters the temperature dependence of the barrier over an extended range to a degree that varies from substance to substance. The lowering of the configurational entropy, however, seems to be the dominant contributor to the barrier increase near the laboratory glass transition, consistent with previous experimental tests of the RFOT theory using the XW approximation. In addition, we compute the temperature dependence of the dynamical correlation length, also without using adjustable parameters. These agree well with experimental estimates obtained using the Berthier et al. (Science 2005, 310, 1797) procedure. Finally, we find the temperature dependence of the complexity of a rearranging region is consistent with the picture based on the RFOT theory but is in conflict with the assumptions of the Adam-Gibbs and "shoving" scenarios for the viscous slowing down in supercooled liquids.

Microscopic calculation of the free energy cost for activated transport in glass-forming liquids

Activated transport in liquids--supercooled liquids in particular--occurs via mutual nucleation of alternative, aperiodic minima of the free energy. Xia and Wolynes [Proc. Natl. Acad. Sci. U.S.A. 97, 2990 (2000)] have made a general argument that at temperatures near the ideal glass transition, the surface penalty for this kind of nucleation is largely determined by the temperature and the logarithm of the size of the vibrational fluctuation of rigid molecular units about the local minimum. Here, we independently show how to estimate this surface tension and, hence, the activation barrier for the activated transport for several actual liquids, using their structure factors and knowledge of the finite-frequency elastic constants. In this estimate, the activation free energy, while depending on the configurational entropy, also depends on the elastic modulus as in the "shoving" models. The resulting estimates are however consistent with the estimate provided by Xia and Wolynes' argument near the glass transition and, in addition, reflect the barrier softening effects predicted earlier for fragile substances.