Research Article

Equilibrium phase diagram of a randomly pinned glass-former

Misaki Ozawa, Walter Kob, Atsushi Ikeda, and Kunimasa Miyazaki
PNAS June 2, 2015 112 (22) 6914-6919; first published May 14, 2015; https://doi.org/10.1073/pnas.1500730112

  1. Edited by Pablo G. Debenedetti, Princeton University, Princeton, NJ, and approved April 16, 2015 (received for review January 13, 2015)

This article has a Letter. Please see:

See related content:

Significance

Confirming by experiments or simulations whether or not an ideal glass transition really exists is a daunting task, because at this point the equilibration time becomes astronomically large. Recently it has been proposed that this difficulty can be bypassed by pinning a fraction of the particles in the glass-forming system. Here we study numerically a liquid with such random pinned particles and identify the ideal glass transition point TK at which the configurational entropy vanishes, thus realizing for the first time, to our knowledge, a glass with zero entropy. We find that as the fraction of pinned particles increases, the TK line crosses the dynamical transition line, implying the existence of an end point at which theory predicts a new type of criticality.

Abstract

We use computer simulations to study the thermodynamic properties of a glass-former in which a fraction c of the particles has been permanently frozen. By thermodynamic integration, we determine the Kauzmann, or ideal glass transition, temperature TK(c) at which the configurational entropy vanishes. This is done without resorting to any kind of extrapolation, i.e., TK(c) is indeed an equilibrium property of the system. We also measure the distribution function of the overlap, i.e., the order parameter that signals the glass state. We find that the transition line obtained from the overlap coincides with that obtained from the thermodynamic integration, thus showing that the two approaches give the same transition line. Finally, we determine the geometrical properties of the potential energy landscape, notably the T- and c dependence of the saddle index, and use these properties to obtain the dynamic transition temperature Td(c). The two temperatures TK(c) and Td(c) cross at a finite value of c and indicate the point at which the glass transition line ends. These findings are qualitatively consistent with the scenario proposed by the random first-order transition theory.

Upon cooling, glass-forming liquids show a dramatic increase of their viscosities and relaxation times before they eventually fall out of equilibrium at low temperatures (1, 2). This laboratory glass transition is a purely kinetic effect because it occurs at the temperature at which the relaxation time of the system crosses the time scale imposed by the experiment, e.g., via the cooling rate. Despite the intensive theoretical, numerical, and experimental studies of the last five decades, the mechanism responsible for the slowing down and thus for the (kinetic) glass transition is still under debate and hence a topic of intense research. From a fundamental point of view the ultimate goal of these studies is to find an answer to the big question in the field: Does a finite temperature exist at which the dynamics truly freezes and, if so, is this ideal glass transition associated with a thermodynamic singularity or is it of kinetic origin (3⇓⇓–6)?

Support for the existence of a kinetic transition comes from certain lattice gas models with a “facilitated dynamics” (6). In these models, the dynamics is due to the presence of “defects” and hence for such systems the freezing is not related to any thermodynamic singularity. However, the first evidence that there does indeed exist a thermodynamic singularity goes back to Kauzmann, who found that the residual entropy (the difference of the entropy of the liquid state from that of the crystalline state) vanishes at a finite temperature TK if it is extrapolated to temperatures below the laboratory glass transition (7). Subsequently many theoretical scenarios that invoke the presence of a thermodynamic transition have been proposed (8⇓–10). One of these is the so-called “random first-order transition” (RFOT) theory which, inspired by the exact solution of a mean-field spin glass, predicts that at TK the glass-former does indeed undergo a thermodynamic transition at which the residual, or configurational entropy Sc (the logarithm of the number of the states which are available to the system) vanishes and concomitantly breaks the replica symmetry (10, 11). A further appealing feature of RFOT is that it seems to reconcile in a natural way the (free) energy-landscape scenario and mode-coupling theory (MCT), a highly successful theory that describes the relaxation dynamics at intermediate temperatures (12).

Despite all these advances, the arguments put forward in the various papers must be considered as phenomenological because compelling and undisputed experimental or numerical evidence to prove or disprove any of these theories and scenarios is still lacking. The only exception is hard spheres in infinite dimensions, for which mean-field theory should become exact (13), but even in this case some unexpected problems are present (see ref. 14). This lack of understanding is mainly due to the steep increase of the relaxation times which hampers the access to the transition point of thermally equilibrated systems, and hence most of the efforts to identifying the transition point, if it exists, resort to unreliable extrapolation.

Randomly Pinned Systems

Recently a novel idea to bypass this difficulty has been proposed (15⇓⇓–18). By freezing, or pinning, a fraction of the degrees of freedom of the system, the ideal glass transition temperature has been predicted to rise to a point at which experiments and simulations in equilibrium are feasible, thus allowing one to probe the nature of this transition. In ref. 15 the authors have studied the effect of pinning for the case of a mean-field spin-glass model which is known to exhibit a dynamical MCT transition at a temperature Td and a thermodynamic transition at a lower temperature TK. It was demonstrated that by pinning a fraction c of the degree of freedoms of the spins (selected at random) in the equilibrated system, both TK(c) and Td(c) increase with c. Thus, by equilibrating the nonpinned system at an intermediate temperature and subsequently increasing c, one can access and probe the ideal glass state in thermal equilibrium [Note that changing c does not perturb this equilibrium (19).] It was found that at sufficiently large c the two lines TK(c) and Td(c) merge and terminate at a critical point with a universality class of the random-field Ising model (15, 20).

Very recently, it has been tested whether this approach to detect TK in mean-field models can also be used in realistic glass-formers in finite dimensions (21). It was found that pinned systems do indeed show a behavior that agrees qualitatively with the theoretical predictions, thus giving encouraging evidence that the nature of the ideal glass transition can be studied in equilibrium. In these numerical studies, the overlap q and its distribution [P(q)] have been used to identify the amorphous-fluid phase diagram in the T–c plane for pinned systems (21).

Despite these results, it is not clear if the so-obtained amorphous state is the bona fide ideal glass for which the configurational entropy Sc vanishes and, in view of the conceptual importance of Sc, this is a very disturbing situation. We recall that Sc is related to the number of available states of the system (10, 11) and is also a key quantity that controls the slow dynamics in the activation regime in which the relaxation time τα is related to Sc via the Adam–Gibbs relation, ln τα∝1/TSc (10, 22). Finally, Sc is intimately related to the Landau-like free energy associated with the overlap of two coupled replicas (23).

Although from the simulations reported in ref. 21 the existence of TK has been inferred from the behavior of the overlap distribution [P(q)] (discussed in more detail below), such an analysis of the overlap q alone may not be conclusive to demonstrate the existence of the arrested phase, because one cannot exclude the possibility that other scenarios, such as the purely kinetic one, can also explain the observed features of [P(q)]. Thus, evaluating Sc directly and identifying the temperature at which it vanishes is crucial to disentangle conflicting theoretical scenarios.

In the present work, we use computer simulation to determine the ideal glass transition temperature TK(c) as a point at which Sc vanishes for a canonical glass-former with pinned particles. For the first time, to our knowledge, this is done without invoking any kind of extrapolation. We also calculate the overlap distribution [P(q)] and find that TK obtained from Sc(TK)=0 and from [P(q)] agree very well up to a finite value of c. Furthermore, we analyze the geometrical properties of the potential energy landscape (PEL) and use it to evaluate the dynamic transition temperature Td(c) as a point at which the saddles of the energy landscape vanish. We find that Td(c) merges with TK(c) exactly at the point at which the two mentioned TK(c) depart from each other, strongly indicating the existence of a critical point which is predicted by the mean-field analysis (15⇓–17).

Results

We study a standard glass-forming model: A 3D binary Lennard-Jones mixture (24). The number of particles is N=150 and 300, but most of the results are for N=300 (see Materials and Methods for details). Throughout the present study, the system has been prepared at each temperature by randomly choosing a fraction c of particles from the thermally equilibrated samples and quenching their positions (see Materials and Methods for details).

Entropy and Configurational Entropy.

To obtain the entropy of the pinned system S, we used thermodynamic integration to determine the entropy of a given configuration of pinned particles and subsequently calculated S by averaging over the realizations of pinned particles (see Materials and Methods for details).

Fig. 1A shows the entropy per (unpinned) particle s≡S/N(1−c). as a function of the fraction of pinned particles at several temperatures T, and we recognize that with increasing c the entropy decreases rapidly. For all temperatures this decrease is linear at small c but then the curves bend at intermediate c and follow a weaker c dependence. For T ≲ 0.5 this bent becomes sharp, strongly indicating that a thermodynamic glass transition takes place.

Fig. 1.
Fig. 1.

(A) Entropy of the system, s, as evaluated from the thermodynamic integration, as a function of c (symbols). The entropies of the disordered solid states svib, obtained using the harmonic approximation, are drawn as solid lines. The dashed lines are a linear extrapolation from the low-c sides. (B) The configurational entropy sc=s−svib. The error bars have been estimated from the sample-to-sample fluctuations.

This becomes more evident by evaluating the configurational entropy obtained by subtracting from S the vibrational entropy Svib. To estimate Svib, we have determined the inherent structures (25) and calculated the eigenfrequencies ωa. Using the harmonic approximation, one can then approximate Svib=∑a{1−log(βℏωa)}. svib≡Svib/N(1−c) is shown in Fig. 1A as well (solid lines) and we see that it shows basically a linear decrease with c, a trend which is due to the suppression of the low-frequency modes in the density of states.

We can now estimate the configurational entropy Sc as the difference Sc=S−Svib (26) and in Fig. 1B we show the c dependence of sc=Sc/N(1−c) for various temperatures. This figure shows that, for T≲0.5, sc quickly decreases with increasing c and becomes basically zero at a finite value of c, indicating that the system has entered the ideal glass state in which the entropy is basically due to harmonic vibrations. [The fact that the difference s−svib is slightly negative is related to the harmonic approximation used to estimate svib. We have estimated the anharmonic correction to svib by considering the nonlinear contribution of the potential energy at the inherent structures (27) and found that the corrections are very small and negative. Evaluation of the exact values of this correction is, however, difficult and beyond the scope of the present work. In principle, one can avoid this difficulty by evaluating the configurational entropy directly (28). However, at present such an approach is computationally way too expensive for systems close to the ideal glass transition.] For T≳0.55, the approach of sc to zero is milder and the bent is less sharp, indicating that the transition becomes a cross-over.

We define the ideal glass transition point cK(T), or TK(c), as the point at which sc becomes zero. As the temperature is lowered, cK(T) decreases and we show the resulting phase diagram, the details of which will be discussed below in Fig. 4. Finally, we mention that the presented results are for N=300. However, we have also simulated systems with N=150 and found that the results do not depend significantly on N (see SI Text for details).

Overlap Approach.

An alternative method to locate and characterize the thermodynamic transition is to study the overlap qαβ between two configurations α and β: qαβ=N−1∑i,jθ(a−|riα−rjβ|), where θ is the Heaviside function, riα is the position of particle i in configuration α, and the length scale a is 0.3 (21). RFOT predicts that at the glass transition the average value [〈q〉] will increase quickly from a small value in the fluid phase to a large value in the glass phase (16). Here 〈⋯〉 and [⋯] stand for the thermal and disorder averages, respectively. We have computed the overlap distribution P(q) using replica exchange molecular dynamics (21) (SI Text) and in Fig. 2 we present [P(q)] for T= 0.7 and 0.45. For T=0.7, [P(q)] remains single-peaked for all c, and the peak position shifts continuously toward larger q as c increases (Fig. 2A). A qualitatively different behavior is observed at T=0.45 (Fig. 2B): [P(q)] is single-peaked at low and high c, but has a double-peak structure at intermediate c, thus signaling the coexistence of the fluid and glass phase, which indicates that the transition from the fluid phase to the glass phase is first-order–like (21).

Fig. 2.
Fig. 2.

Distribution of the overlap [P(q)] for T=0.7 (A) and T=0.45 (B). (C) The average overlap [〈q〉] obtained from [P(q)] as a function of c for several temperatures. The error bars have been calculated using the jackknife method.

The c dependence of the average overlap [〈q〉] is shown in Fig. 2C. For high temperatures, T≳0.6, [〈q〉] smoothly increases with c, reflecting the continuous shift of the single peak of [P(q)] as shown in Fig. 2A. For T≲0.5, [〈q〉] shows a quick increase at intermediate values of c, in agreement with the presence of the double-peak structure seen in [P(q)] at these T. It suggests that a first-order–like transition, rounded by finite size effect, takes place, in qualitative agreement with the result for a system of harmonic spheres (21). Note that, within the accuracy of our data, we see that the amplitude of the (smeared out) jump in [〈q〉] seems to vanish around 0.55≲T≲0.6, thus indicating that at around that temperature the line of first-order transition ends in a critical point that is second-order–like.

From this approach with the overlap, we can define the ideal glass transition temperature TK(q)(c) as the temperature at which the skewness of [P(q)] vanishes (SI Text), and in Fig. 4 we have included the resulting TK(q)(c) as well. For small and intermediate c, we find a very good agreement between TK(c) and TK(q), thus showing that the two very different approaches do give the same ideal glass transition temperature. This is thus very strong evidence that at this temperature the system does indeed undergo a thermodynamic phase transition from the fluid to an ideal glass state. For temperatures above the second-order–like critical point the curves TK(c) and TK(q)(c) differ. We will rationalize this in Discussion below. Finally, we mention that the curves TK(c) and TK(q)(c) seem to extrapolate in a smooth manner to the Kauzmann temperature of the bulk which has been estimated from extrapolation (26). This shows that the present measurement of the bulk TK is compatible with the one from previous estimates.

PEL and MCT.

In the past, it has been found that the slow dynamics of glass-forming systems is closely related to the features of the PEL (27) and in the following we will use these relations to characterize the relaxation dynamics of the pinned system.

Fig. 3A shows the T dependence of the average inherent structure energy [〈eIS〉]. For the bulk system, c=0, [〈eIS〉] is basically constant at high temperatures but then steeply decreases below a cross-over temperature T≈1, a temperature which signals that the relaxation dynamics becomes strongly influenced by the PEL (29, 30). As c increases, the value of [〈eIS〉] at high T moves steadily upward, which is reasonable because the energy of the system is literally pinned at the higher energy levels due to the presence of the pinned particles. Concomitantly the cross-over temperature increases with c and the cross-over becomes smeared out, completely disappearing at the highest c. The vanishing of this cross-over with growing c indicates thus that the pinning qualitatively affects the nature of the PEL and of the relaxation dynamics. For instance, it is found that with increasing c the fragility of the system decreases and shows at high c an Arrhenius dependence (21, 31).

Fig. 3.
Fig. 3.

(A) T dependence of the averaged inherent structure [〈eIS〉] for several values of c. (Inset) Same quantity as a function of 1/T. (B) The average normalized saddle index k as a function of its energy eSP for several values of c. The threshold energy eth is defined by a linear extrapolation from the high-energy side.

Fig. 4.
Fig. 4.

Phase diagram of the randomly pinned system. The filled circles show the ideal glass line TK at which the configurational entropy vanishes. The diamonds show TK(q)(c) determined from the skewness of [P(q)]. The squares are dynamic transition points Td(c). The open circles show the ideal glass line TK′ determined by the linear extrapolation of Sc(c) to vanish from Fig. 1B. The iso-relaxation times are drawn by +, ×, and *. The star denoted at c=0 is a putative ideal glass transition point TK≈0.3 for the bulk reported in ref. 26.

In Fig. 3A, Inset, the low-temperature behavior of [〈eIS〉] is shown. It clearly demonstrates that [〈eIS〉] is inversely proportional to T for all c. According to the energy landscape scenario, this is an indicator that the distribution of eIS is Gaussian (32).

Another important quantity that connects the glassy dynamics of a system with its PEL is the saddle index K, i.e., the number of negative eigenvalues of the Hessian matrix at a stationary point of the PEL. For bulk systems it has been found that K shows a linear dependence on eSP, the bare energy of a saddle point (33, 34). Because the value of eSP at which K goes to zero (this value is often denoted as “threshold energy” eth) corresponds to the average energy of the inherent structures 〈eIS〉 at the critical temperature of MCT, one can extract from the geometrical properties of the PEL the value of Td without having to do any fit to dynamical data (33, 34).

We use a standard method to determine numerically the energy and index of saddles for the pinned system (see SI Text for details), and in Fig. 3B we plot the average normalized saddle index k=K/3N(1−c) as a function of its corresponding eSP. In agreement with previous studies of the PEL, we find that k decreases linearly as a function of eSP and hence we can obtain eth(c) from a linear fit (included in the figure as well). We see immediately that eth increases with c and, together with the data of Fig. 3A and eth(c)=[〈eIS〉](Td(c)), we can conclude that Td(c) increases with c. The resulting c dependence of Td is included in Fig. 4 as well and will be discussed in the next section.

We have also evaluated Td by calculating the relaxation time τα from the time-dependent density correlation function (Material and Methods) and by fitting τα with the MCT power law ταα|T−Td (fit)|−γ. We find that the so-obtained values of Td (fit) agree well with those obtained from the PEL (see SI Text for details). Note that one needs several fit parameters to obtain Td (fit) from the dynamic data, whereas basically no fit parameters are required for Td from the PEL. In Fig. 4, the iso-τα lines are also plotted. The graph shows that τα quickly increases with c and that the lines asymptotically approach the TK line from the high-T side.

Discussion

In Fig. 4 we summarize the results of the previous sections in the form of a phase diagram in the c–T plane. The ideal glass transition lines TK(c) determined from Sc=0 and TK(q)(c) obtained from [P(q)] are plotted as filled circles and diamonds, respectively. At low T and c the two temperatures basically coincide but around (Tc,cc)≈(0.55,0.16) they start to depart from each other in that TK(q)(c) increases continuously, whereas TK(c) bends and its c dependence becomes weaker. Note that this separation occurs at the same point at which [P(q)] changes from the bimodal to the single-peaked shape (Fig. 2), thus indicating that in this region of parameter space the system has a second-order–like critical point.

The theoretical calculations for a mean-field spin-glass model show that TK(c) should terminate at a finite c and that this end point is a critical point of the universality class of the random field Ising model (15⇓–17, 20). Furthermore, the theory predicts that at this end point the coexistence line TK(c) and the dynamical line Td(c) merge. Fig. 4 shows that this prediction is indeed compatible with our data in that the two lines do cross near (Tc,cc).

From the figure we also recognize that, beyond the end point, TK(c) obtained from Sc=0 almost matches with Td(c). This result is reasonable because in this range of T and c the particles are strongly confined by the labyrinthine structure imposed by the pinned particles, i.e., the system resides mainly at the bottom of a free-energy minimum where both the configurational entropy and the number of the saddles vanish simultaneously. In contrast with TK(c), TK(q)(c) raises continuously even beyond the end point. This can be understood by recalling that this line is defined by the points at which the skewness becomes zero, i.e., the point at which [P(q)] changes from being left-skewed to right-skewed. Because this change of sign in the skewness occurs also in the region of the phase diagram beyond the critical point, the line TK(q) will extend into that region, similar to the Widom line present in a standard liquid–gas transition (35, 36).

As it is evident from Fig. 1A, TK, defined as the point at which Sc=S−Svib=0, becomes somewhat ill-defined beyond the end point (T≳0.55 in Fig. 1A), because the cross-over from the fluid to the glass phase becomes broad. To estimate the effect of this ambiguity, we have also included in Fig. 4 TK′ defined by the linear extrapolation Sc of Fig. 1B from the low-c side (open circles). We see that below the end point T ≲ 0.5, TK′ is indistinguishable from TK, whereas beyond the end point, T≳0.6, they bifurcate and TK′ becomes comparable with TK(q).

To the best of our knowledge, the present study is the first report of a system in finite dimensions that shows the existence of an ideal glass state in equilibrium, i.e., a state in which the configurational entropy is zero at a finite T. The Kauzmann temperatures reported in the past have all relied on somewhat questionable extrapolation procedures, leaving thus room for debate over the very existence of a thermodynamic transition (6, 8, 37).

Our findings are inconsistent with recent simulation studies in which the T- and c dependence of the relaxation dynamics has been studied (38). In ref. 38, the structural relaxation time τα was fitted with the Vogel–Fulcher relation τα∼exp[A/(T−T0)], with T0 as fit parameter. This relation can be directly derived from the Adam–Gibbs as well as the RFOT theory, both of which assert that ln⁡τα∝1/TSc, and thus the Vogel–Fulcher temperature T0 is predicted to be identical to TK. Furthermore, the authors of ref. 38 fitted their data also with the MCT power law ταα|T−Td|−γ to determine the c dependence of Td. It was found that whereas Td increases moderately with c, T0 remains constant. We have plotted the relaxation time τα as a function of Sc for finite c and found that the Adam–Gibbs relation is violated (see SI Text for details). Thus, we conclude that in the case of pinned systems one cannot deduce the Vogel–Fulcher law from the Adam–Gibbs relation.

Because the results presented here are all obtained in thermodynamic equilibrium without referring to any kind of extrapolation, we are confident that the phase diagram presented in Fig. 4 does indeed reflect the properties of the system and is not an artifact of the analysis. Further evidence that the simulated system is really in equilibrium is the observation that the entropy obtained by thermodynamic integration from the high-temperature limit matches with that obtained from the low-temperature side (via harmonic approximation) in the glass phase. It is also reassuring that all three methods, the thermodynamic integration (vanishing entropy), the overlap distribution (the discontinuous jump of q), and the geometric change of the PEL, consistently point to the same end point, thus giving strong evidence that this point really exists. Also suggestive is that each combination of pairs among the three methods is compatible beyond the end point, which is reminiscent of the Widom line in the standard gas–liquid phase transition.

At this stage we can conclude that the phase diagram as predicted by the RFOT theory is confirmed at least qualitatively. What remains to be done is to probe the relaxation dynamics in the vicinity of the critical end point because one can expect that this dynamics is rather unusual (39) and to establish its universality class (20, 40). Furthermore, it will also be important to see whether the predicted phase diagram can also be observed in real experiments. Although this will be not easy, for certain systems such as colloids or granular media it should be possible.

Materials and Methods

Model.

The system we use is a binary mixture of Lennard-Jones particles (24). Both species A and B have the same mass and the composition ratio is NA:NB=80:20. The interaction potential between two particles is given by vαβ(r)=4ϵαβ{(r/σαβ)12−(r/σαβ)6}, where α,β∈{A,B}. We set ϵAA=1.0,ϵAB=1.5,ϵBB=0.5,σAA=1.0,σAB=0.8, and σBB=0.88. vαβ(r) is truncated and shifted at r=2.5σαβ. We show energy in units of ϵAA, with the Boltzmann constant kB=1, and length in units of σAA. Time units are defined by Monte Carlo sweeps (see below). Simulations are performed at constant density ρ≈1.2. The number of particles is N=150 and 300, and most of the results in the present study are for N=300.

Making Pinned Configurations.

The configuration of the pinned particles is generated by making first a replica exchange run for the bulk system, i.e., c=0, using eight replicas (41). This allows us to generate relatively quickly many equilibrium configurations that are completely independent, i.e., between consecutive configurations the mean-squared displacement of a tagged particle is more than 100. Next we use a “template” to identify the cN particles that will be permanently pinned. Details on how to create the template can be found in ref. 21. This approach has the advantage that the pinned particles cover the space in a relatively uniform manner, thus avoiding the creation of dense regions or large empty regions and hence suppressing strong sample-to-sample fluctuations of the thermodynamic properties of the system.

Simulation Methods.

Thermodynamics.

To sample thermodynamic properties efficiently at low-T and large-c region, we use the replica exchange method (41). The maximum number of replicas is 24. More detail is presented in SI Text and in ref. 21. The total CPU time to obtain the presented results is about 580 y of single core time.

Dynamics.

We use the Monte Carlo (MC) dynamics simulation to calculate dynamical observables (42). The rule of the MC dynamics is the following: In an elementary move, one of the (1−c)N unpinned particles is chosen at random. Then the particle is displaced at random within a cubic box of linear size δ=0.15 and the standard Metropolis rule is used to decide whether or not the move is accepted. One MC step consists of (1−c)N such attempts and we set this as a unit of time scale. The relaxation time τα is determined by FsA(k,τα)=e−1, where FsA(k,t) is the self part of the intermediate scattering function of the free particles of species A for the wave vector k at the peak of the corresponding structure factor. We have averaged over 20–30 different realizations of pinned particles to calculate FsA(k,t).

Entropy.

To calculate the entropy S(c,T) of the pinned system, we have first determined the entropy of the system for a given configuration of pinned particles, and then taken the average over the realization of the configuration of pinned particles. For this we have calculated the potential energy at temperatures ranging from the target temperature up to the ideal gas limit at T=∞, while keeping the pinning configuration fixed. We have evaluated the entropy of the system with that pinning configuration using thermodynamic integration. For this integration, we have used a grid in the inverse temperature β=1/T of width Δβ that ranged between 0.01 and 0.1, depending on the temperature, and integrated the potential energy as a function of β. Special care was taken in the very high temperature regime, to accurately and rapidly achieve the convergence to high-temperature ideal gas limit. The high-temperature expansion of the potential energy of Lennard-Jones fluid can be written as U=Aβ−1/4+Bβ−2/4+Cβ−3/4 (43). We have used simulations at very high temperatures to determine the coefficients A, B, and C, and then carried out the thermodynamic integration analytically.

Analysis of the Saddles.

To locate the saddles of the PEL of the system, we have made a minimization of the squared gradient potential W=1/2|∇U|2 (33, 34). Minimization of W is performed by the Broyden–Fletcher–Goldfarb–Shanno algorithm (44). Similar to the minimization of U used to calculate the inherent structures, W includes the position of the pinned as well as unpinned particles, but only the position of the latter is optimized. After having located a saddle with energy eSP, the Hessian matrix was diagonalized and we counted the fraction of negative eigenvalues k(eSP). The raw data are shown in SI Text and in Fig. 3B we present the average index as a function of eSP. The threshold energy eth is defined by k(eth)=0.

Acknowledgments

We thank G. Biroli, C. Cammarota, D. Coslovich, and K. Kim for helpful discussions. M.O. acknowledges the financial support by Grant-in-Aid for Japan Society for the Promotion of Science Fellows (26.1878). W.K. acknowledges the Institut Universitaire de France. A.I. acknowledges JSPS KAKENHI 26887021. K.M. and M.O. acknowledge KAKENHI 24340098, 25103005, 25000002, and the JSPS Core-to-Core Program. The simulations have been done in Research Center for Computational Science, Okazaki, Japan, at the HPC@LR, and the CINES (Grant c2014097308).

Footnotes

  • ↵1To whom correspondence should be addressed. Email: walter.kob{at}univ-montp2.fr.
View Abstract
Back to top
Article Alerts
Email Article
Citation Tools
Request Permissions
Equilibrium pinned glass-former
Misaki Ozawa, Walter Kob, Atsushi Ikeda, Kunimasa Miyazaki
Proceedings of the National Academy of Sciences Jun 2015, 112 (22) 6914-6919; DOI: 10.1073/pnas.1500730112
Digg logo Reddit logo Twitter logo Facebook logo Google logo Mendeley logo
Proceedings of the National Academy of Sciences: 112 (22)
Table of Contents

Submit

Article Classifications

Jump to section

You May Also be Interested in

Similar Articles