Misaki Ozawa

Atsushi Ikeda

Kunimasa Miyazaki

Walter Kob
email: walter.kob@umontpellier.fr

Ideal glass states are not purely vibrational: Insight from randomly pinned glasses

Keywords: numbers: 64.70.kj, 63.50.Lm, 64.70.Q-

We use computer simulations to probe the thermodynamic and dynamic properties of a glassformer that undergoes an ideal glass-transition because of the presence of randomly pinned particles. We find that even deep in the equilibrium glass state the system relaxes to some extent because of the presence of localized excitations that allow the system to access different inherent structures, giving thus rise to a non-trivial contribution to the entropy. By calculating with high accuracy the vibrational part of the entropy, we show that also in the equilibrium glass state thermodynamics and dynamics give a coherent picture and that glasses should not be seen as a disordered solid in which the particles undergo just vibrational motion but instead as a system with a highly nonlinear internal dynamics.

The nature of the glass transition is one of the most challenging research topics in condensed matter physics and has therefore been in the focus of a multitude of studies [1][2][3][4][5]. Many aspects of the slow and complex dynamics of supercooled liquids and the resulting glass transition have been successfully explained in terms of the potential energy landscape (PEL) [6][7][8]10]. In this framework the configurational space of the system is partitioned into basins of attractions of the local energy minima (the inherent structures, IS) of the potential energy and the dynamics at low temperatures is characterized as the motion through the complex pathway that connects neighboring basins. The conventional description of glasses is that the glass state corresponds to a vibrational motion around the IS, which in real space means that the particles are trapped by the cages formed by their neighbors and vibrate around a fixed amorphous configuration. However, recent experiments challenge this view since they seem to suggest that certain types of relaxation processes are present even in the glass, implying that the motion of the atoms is more complex than pure vibrations [11][12][13]. However, it is difficult to decide whether or not these relaxation processes are indeed an equilibrium property of the sample or just related to aging. The existence of equilibrium relaxation processes in the glass state implies that even at low temperatures the system explores a complex landscape, i.e., it can access in a finite time many different local minima. This in turn has the consequence that the conventional definition of the configurational entropy s conf as the difference between the total entropy s tot of the system and its purely vibrational part s vib should be questioned. Thus the study of (equilibrium!) relaxation processes in the equilibrium glass state will allow to advance our understanding of the meaning of the glass state on the microscopic level.

Advancing on this question has so far been hampered by the fact that it was impossible to generate equilib-rium glasses (also called "ideal glasses"), i.e., equilibrium structures at very low temperatures. This situation has recently changed since it has been realized that if one pins (immobilizes) randomly a finite fraction c of the particles [14], the fluid particles, i.e., non-pinned ones, undergo an equilibrium glass transition if c is increased beyond a certain threshold [15,[START_REF] Berthier | [END_REF]. Numerical simulations of simple glass-formers have confirmed that this pinning approach does indeed allow to observe an equilibrium glass transition at which the entropy shows a marked bend and hence to access the equilibrium glass state [9,10]. Thus these results have opened the door to study the properties of equilibrium glasses and in the following we will demonstrate that even in the ideal glass state relaxation processes are present, implying that in the glass the entropy is not just given by the vibrational contribution.

We simulate a binary mixture of N Lennard-Jones (LJ) particles in three dimensions [6] of which a fraction c are permanently pinned. N is 300 or 1200 and we use the standard LJ units for the length, energy, and temperature, setting the Boltzmann constant k B = 1. First we have equilibrated the system without pinned particles at a given temperature T , then the positions of cN particles are frozen and we study the static and dynamical properties of the remaining (1 -c)N mobile particles of the system. To study the static properties, we use the parallel tempering (PT) molecular dynamics method with 24 replicas [20,21] (see Ref. [9,10] for details on the pinning procedure and the PT). For the dynamical properties we use the standard Monte Carlo (MC) dynamics simulation [22], where an elementary move is a random displacement of a randomly chosen particle within a linear box size δ = 0.15. One MC step, which is our unit of time, consists of (1 -c)N such attempts. Despite the huge acceleration of the sampling due to the PT, we found that runs up to 2 • 10 10 steps were needed to get statistically significant results. This, and the necessary averaging over the independent realizations of the pinned configurations (typically 25 for N = 300 and 10 for N = 1200).

To probe the relaxation dynamics of the system we focus on T = 0.45 and increase c, thus crossing at around c = 0.1 the boundary between fluid and glass state (see inset of Fig. 1(b)). Figure 1(a) shows the mean squared displacement (MSD) for different values of c:

∆(t) = (1-c)N i=1 [ |r i (t) -r i (0)| 2 ]/(1 -c)N . Here • • • and [• • • ]
are the thermal and disorder averages, respectively. Note that we have averaged the MSD over both species of particles. At intermediate times ∆(t) has a marked plateau that is related to the usual cage effect observed in glassy systems [1] and the horizontal dashed line indicates the plateau height of ∆(t) at c = 0.2. Surprisingly we find that even in the equilibrium glass phase, i.e., at this temperature for c ≥ 0.10 (see SI), ∆(t) shows at long times a marked increase above this plateau, indicating that the particles can leave their cage. To understand this behavior, we analyze the van-Hove correlation function (see SI) and find that the particles undergo an exchange motions, i.e., the particles tend to be replaced by the same kind of particles (see SI). This seems to explain the upturn of the mean squared displacement. More surprising is the results regarding the collective overlap

Q(t) = i,j [ θ(a -|r i (t) -r j (0)|) ]
, where θ(x) is the Heaviside function and a = 0.3. This function probes the collective relaxation of the system and it is not affected by the particle exchange motions. Figure 1(b) shows that Q(t) decays slightly after the plateau at intermediate times, before it decays to its long time limit (horizontal dashed lines

Q(t → ∞) = Q (static)
), a quantity that can be calculated with high accuracy directly from the PT simulations. The very high value of Q (static) > 0.85 and the strong c-dependence of this quantity, see Ref. [10], demonstrates that the explored state points are indeed glass states and hence we can conclude that the system shows subtle and non-trivial relaxation dynamics even in the ideal glass. In the SI we show that this is not an out-of-equilibrium effect. The presence of this relaxation dynamics is at odds with the results of Ref. [10] that in the glass phase the configurational entropy s conf seems to be zero (see red circles in Fig. 3(b)), since s conf = 0 implies that there are no states into which the system can move to relax. In that work s conf was estimated from s tot -s harm , where s harm is the entropy of the strictly harmonic solid, i.e., a quantity that can be obtained directly from the vibrational density of states. Our present results rise thus the question whether s conf can indeed be approximated reliably by this difference [23][START_REF] Berthier | Proc. Natl. Acad. Sci[END_REF][START_REF] Johari | [END_REF].

To advance on this point, we have accurately determined s vib by taking into account the effect of anharmonic vibrations, using two independent approaches. The first one makes use of the idea of Frenkel and Ladd [15,[26][27][28][START_REF] Williams | [END_REF] of introducing a series of systems which parametrically interpolates between the original system and an Einstein solid, and then to carry out a thermodynamic integration to calculate the vibrational entropy of the original system. In practice we have introduced the Hamiltonian

βH(α) = βH(0) + α (1-c)N i=1 |r i -r 0i | 2
, where H(0) is the original Hamiltonian, β = 1/T , α is a spring constant, and r 0i is the equilibrium configuration of the particle i in the original system. We then measure the

MSD ∆ α = (1-c)N i=1 [ |r i -r 0i | 2 ]/(1 -c)N from which we obtain the entropy s FL = s E (α max ) + αmax 0 dα∆ α , (1)
where s E (α max ) is the entropy of the Einstein solid. Figure 2 shows ∆ α for 10 -1 ≤ α ≤ 10 7 . For α > 10 4 , this function follows very closely 3/(2α), the behavior of the Einstein solid, and hence we can replace ∆ α by this expression if α is large. The entropy for this Einstein solid is then given by s

E (α max ) = 3 2 -3 ln Λ -3 2 ln αmax π
, where Λ is the de Broglie thermal wavelength. In practice we have set α max = 10 7 .

At small α, the dependence of ∆ α on α becomes weak suggesting that in this parameter range the particles are vibrating in a cage created by the original Hamiltonian and not by the one of the Einstein solid. However, at the smallest α, ∆ α is not completely flat, since at long times the particles are able to escape slowly from their cages. The height of the plateau in ∆ α allows to estimate the amplitude of the vibrations in the real system, i.e., to determine the harmonic and anharmonic component of the motion. Indeed, the height of the plateau in ∆ α is consistent with that of ∆(t) from Fig. 1(a). To estimate the entropy that is due to the vibrational motion we can evaluate the integral given by Eq. ( 1) by replacing the contribution to the integral for α < α min = 2 by α min ∆ αmin , thus removing in this manner the contribution of the relaxational part to the entropy. The so obtained vibrational entropy s FL is shown in Fig. 3(a) (green diamonds). As expected, the c-dependence of s FL is smooth and shows no apparent singularity.

The second approach to estimate the anharmonic contribution to the vibrational entropy is to determine the difference between the potential energy of the system and its inherent structure energy and then to subtract the harmonic part, giving thus the anharmonic part of the vibrational energy [7,8]. One finds that the so obtained energy is negative which makes that the resulting value for the vibrational entropy, s harm + s anh , is smaller than s harm (see blue triangles in Fig. 3(a)), a result that agrees with previous studies [7,8]. We see that this quantity agrees very well with our estimate for the vibrational entropy as obtained from the Frenkel-Ladd procedure, indicating that we have determined it with good precision.

Figure 3(b) shows the c-dependence of ∆s = s tots vib . We see that the improved estimate for s vib makes that now ∆s no longer goes to zero even at large c. Instead it shows a kink at the concentration at which the order parameter had a jump [10], indicating that at this point the thermodynamic properties of the system have a singular behavior, i.e., that the phase diagram determined in Ref. [10] is not altered by this improved estimate of s vib . Note that we show data for the two system sizes and within the accuracy of the data we see no finite size effects.

To determine the origin of the finite value of ∆s in the ideal glass phase we probe the potential energy landscape of the system [6][7][8]10]. Since at low T the configuration space can be decomposed into the basins of attraction of the inherent structures and vibrations around the IS, an investigation of the PEL should help to understand the nature of the motion of the system. Figure 4(a) shows the time dependence of the IS energy e IS (t), [32,33], for 10 4 steps of MC dynamics at c = 0.2, a timescale that corresponds to the vibrations inside the cage (see Fig. 1(a)). We see that e IS (t) remains basically constant but shows some spikes indicating that the system accesses for a short period an excited state before it falls back to the original IS, implying that on this time scale only a single IS is relevant [10]. Figure 4(d) shows the superposition of the configurations for these IS's and we see that , we find that e IS (t) starts to make larger excursions and that these are no longer reversible, indicating that the system explores new IS. This is the reason why on this time scale the correlation function decay slightly, see Fig. 1(b). Figure 4(e) shows that now there are indeed many different IS's, but since the corresponding positions of the particles form small clusters we can conclude that the configurations are quite similar. Finally we show in Fig. 4(c) the evolution of e IS (t) as obtained from a PT run. We see that in this case the value of e IS (t) fluctuates quite significantly, but that these fluctuations are still smaller than the sample to sample fluctuations and therefore we can infer that there are indeed many different IS's even in the equilibrium glass state. This conclusion is corroborated by the real space image of the configurations, Fig. 4(f), which now shows quite a few small clusters, i.e., the positions of the particles in the different IS's differ by a small amount. The presence of these different IS's is thus the reason for the finite value of ∆s. We emphasize that the clusters seen in Fig. 4(f) do not depend on the length of the PT run.

Our results show that the dynamic and thermodynamic properties of a system with pinned particles give a coherent picture of the equilibrium glass state. Although quantities like the self intermediate scattering function [3,34,35] or the self van Hove function (see SI) do decay to zero, the collective functions do not, once the critical pinning concentration is surpassed. At this transition point the order parameter shows a jump [10] and the entropy a kink. Our finding that even in the equilibrium glass state particles are able to explore more than one IS indicates that care has to be taken in the definition of the configurational entropy: s conf should not be identified as the difference between the total entropy and the purely vibrational part, since such a definition gives rise to a non-zero s conf even in the ideal glass state. Instead s conf is related to the number of local free energy minima [START_REF] Biroli | [END_REF] and for the case of the pinned system such a minimum is a collection of IS that are geometrically close together. The motion of the particles inside this local free energy minimum is the reason for the partial relaxation of the collective correlation functions. Although the present study concerns a pinned system, it can be expected that bulk systems have a similar behavior since local inhomogeneities in the structure will make that certain regions in the very deeply supercooled liquid are already frozen in, thus leading to a very high value of the viscosity, whereas other regions are still mobile. So the evoked problem with the correct definition of the configurational entropy is likely to exist also in the case of bulk systems.

In Fig. 1(a) of the main text we show that the mean squared displacement (MSD) ∆(t) of a tagged particle shows at intermediate times the usual plateau found in glassy systems [1] but that, surprisingly, the MSD turns at long times upward even when the system is in the ideal glass phase. In order to demonstrate that this upturn is not related to any kind of non-equilibrium dynamics, see e.g. Ref [2], we show in Fig. 5(a) the MSD as obtained for different origins of time. In practice we have taken the trajectory as obtained from the parallel tempering (PT) run and selected five different configurations from that trajectory. For each one of these configurations we have started a standard MC run and have obtained in this manner five different MSDs. The figure shows that these curves are identical with high precision, i.e., they do not depend on the starting configurations. This is hence strong evidence that the upturn of ∆(t) seen at large times is not related to aging dynamics since in that case one finds that the time for the upturn shifts to increasingly longer times if the age of the system, i.e., the time at which the initial configuration for the MC run was selected, is increased [2]. The same conclusion is reached if one considers the time dependence of the overlap, again calculated starting from different starting configurations, see Fig 5(b).

Dynamics as characterized by the van Hove function

In Fig. 1 of the main text we show the time-dependence of the self dynamics (i.e., the mean squared displacement) as well as a collective correlation function, i.e., the collective overlap function Q(t). In Ref. [3] we have argued that because of the presence of the pinned particles the self correlation functions behave very differently from the collective ones in that the latter will decay to a plateau with finite value whereas the former will go to zero as in a normal liquid. To understand better the motion of the particles in real space it is useful to look at the self and distinct parts of the van Hove correlation functions [4].

The self part of the van Hove function is defined by

G α s (r, t) = 1 (1 -c)N α (1-c)Nα i=1 [ δ(r -|r i (0) -r i (t)|) ] , (2)
where N α is the number of particles of species α. Here α, β ∈ {A, B}. The distinct part of the van Hove function is defined by

G αα d (r, t) = N (1 -c)N 2 α (1-c)Nα i=1 (1-c)N β j( =i) [ δ(r -|r i (0) -r j (t)|) ] , (3)
G αβ d (r, t) = N (1 -c)N α N β (1-c)Nα i=1 (1-c)N β j [ δ(r -|r i (0) -r j (t)|) ] (α = β). (4)
The r-dependence of G s (r, t) is shown in Fig. 6 for temperature T = 0.45 and different pinning concentrations c. For c = 0.0167, i.e., when the system is very similar to a glassy bulk system, the G s (r, t) shows at very short times the Gaussian behavior expected for a harmonic system since on this time scale the particle explores its local cage. With increasing time the particle leaves that cage and accesses larger distances before becoming diffusive. If the concentration is increased to c = 0.1, Fig. 6(b), we see that at large times the distribution function has a marked peak at r ≈ 1, i.e., the particle hops by a nearest neighbor distance. Small peaks are seen at larger distances at r ≈ 2 and 2.7, which correspond to the second and third nearest neighbor distances [5]. If c is increased even more, Fig. 6(c), the form of G s (r, t) is not altered appreciably, although the probability at large distances decreases by more than one order of magnitude. Thus we can conclude that the presence of the pinned particles change the nature of the motion from a flow-like motion at low c to a hopping motion at intermediate and high c.

The self part of the van Hove functions discussed so far were for the A particles which are larger and thus less mobile than the B particles. In Figs. 6 (d)-(f) we show the same quantity for the B particles. In agreement with the results for the bulk system [6], we see that already for very low concentration there is a small peak at r ≈ 1, i.e., the particles show some tendency for a hopping dynamics. If c is increased the peak becomes more pronounced but remains smaller than the corresponding peak for the A particles.

The r-dependence of the distinct part of the van Hove function is shown in Fig. 7. By definition, G αβ d (r, t = 0) is the radial distribution function, which is zero at r 1 and with several distinct peaks corresponding to the nearest or next nearest distances. All panels show that the hole of G αβ d (r, t = 0) at r 1 is filled up by other particles at longer times. For the A-A component, panels (a)-(c), one observes that, for a fixed t, the peak developing at r = 0 becomes sharper and the height of the peak becomes lower as the c is increased. We find a very similar behavior for the B-B component, panels (d)-(f). The peaks are slightly higher than those of the A-A component because the B particles can undergo a hopping motion more easily since they are smaller than the A particle. On the contrary, for the A-B component, panels (g)-(i), the growth of the peak is somewhat weaker. These results imply that, in the hopping motion, the particles tend to be replaced by the same kind of particles.

Note that, although the exchange of two particles of the same type allows that the particles undergo a diffusive motion at large times, the underlying structure does not change (i.e., it is like the diffusion of defects in a crystal). Hence this type of motion will not give rise to a change in the overlap function. This is different for the exchange of an A particle with a B particle since an exchange of unlike particles will lead to a non-equivalent configurations. As mentioned above, we recognize that in this case the peak at the origin is now much smaller which is explained by the fact that the size of the A particle is significantly larger than the one of the B particles and hence it is difficult to swap the position of the unlike particles. Nevertheless, at sufficiently long times a small peak can be observed even in this case and hence the system is indeed able to access new configurations, in agreement with the findings on the entropy and PEL analysis presented in the main text.

Thermodynamics Anharmonic contribution for the vibrational entropy

In this section we show how to evaluate the anharmonic contribution for s vib . Details of this method are described in Refs. [7,8] for the bulk system. The extension of the method to the pinned system is straightforward.

The anharmonic contribution of the total potential energy U anh (c, T ) is given by

U anh (c, T ) = U (T ) -U IS (c, T ) - 3 2 (1 -c)N T, (5)
where U (T ) and U IS (c, T ) are, respectively, the equilibrium bulk potential energy and the inherent structure energy of the pinned system. The last term in Eq. ( 5) is the energy of the harmonic vibrations of (1 -c)N mobile particles. Making a T -expansion of U anh (c, T ) around T = 0 gives

U anh (c, T ) = k=2 C k (c)T k , (6)
where C k (c) are T -independent coefficients. Note that the sum starts at k = 2 since the system is completely harmonic in the low temperature limit. The anharmonic contribution for the vibrational entropy S anh (c, T ) is given by

S anh (c, T ) = T 0 dT ′ 1 T ′ ∂U anh (c, T ′ ) ∂T ′ . (7)
Note that we set S anh (c, T = 0) = 0.

Substituting Eq. ( 6) to Eq. ( 7), we obtain

S anh (c, T ) = k=2 k k -1 C k (c)T k-1 . (8)
To evaluate S anh we have used the following numerical procedure: First we use the simulations to obtain U anh (c, T ) and then the coefficients C k (c)'s are obtained by fitting U anh (c, T ) to a polynomial function of T . In practice we use the first two terms, C 2 (c) and C 3 (c), for this fitting [7]. Finally S anh (c, T ) is evaluated using Eq. ( 8). It is clear that the signs of C k (c)'s determine the sign of S anh (c, T ). In Fig. 8 We find that s anh (c, T ) is small and negative for all c's and T 's considered, in agreement with the results from Refs. [7,8].

Overlap and Finite Size Effects

In Fig. 9(a) we show the distribution P (Q) of the overlap Q αβ between two configurations α and β which is defined by

Q αβ = 1 N N i,j θ(a -|r α i -r β j |), (9)
where θ(x) is the Heaviside function, r α i is the position of the i-th particle in the configuration α, and the lengthscale a is set to 0.3 [9]. For small c we see a single peak at small Q which shows that for weak pinning the overlap between two typical equilibrium configurations is small. For intermediate c we find a pronounced double peak structure indicating the coexistence of states that have a low overlap and states that have a high overlap with each other. If c is increased further one finds only a single peak at large Q, i.e., most configurations are very similar. These results are for N = 1200 particles and are qualitatively very similar to the one for N = 300 presented in Ref. [10]. To show the N -independence more quantitatively, we present in Fig. 9(b) the c-dependence of the mean overlap Q (static) which is given by the first moment of P (Q). We see that within the accuracy of the data the curves for N = 300 coincide with the one for N = 1200. We also observe that the N -dependence of s tot -s vib at the lowest temperature is small (see Fig. 9(c)). This result is somewhat surprising since one expects noticeable finite size effect in the vicinity of the ideal glass transition due to some growing static correlation lengths [11,12] and, thus of the N -dependence of the order-parameter Q (static) . One possibility to rationalize the absence of such finite size effects is that the finite size exponents is very small in the temperature range which we explore [13]. In order to verify this, one must study lower temperatures with more efficient algorithms or consider other model systems [14,15].

On this point, we note that our results are not very sensitive to the temperature, see e.g. Fig. 9(c) for the difference s tot -s vib , and therefore we expect them to be relevant also at lower temperatures.

Potential energy landscape

Snapshot

In Fig. 4 of the main text we show the superposition of snapshots of IS configurations for N = 300 particles. In Fig. 10 we present similar data for N = 1200. Panel (a) shows the snapshots for c = 0.0167, i.e., the case where the system is still in the liquid-like regime. We see that the position of the particles basically fill out all the accessible space, with the exception of the region around the pinned particles. For c = 0.08, i.e., slightly below the transition, some empty regions can be spotted, indicating that now the particles can no longer access the whole space, panel (b). If the concentration is increased to c = 0.167, i.e., a state that is deep in the glass state, most of the IS are very similar and hence the available configuration space has become very small, panel (c). However, even in this case we find that there are small clusters of particle positions that indicate that a glass state is composed of many different IS, which is thus the reason for having a non-trivial entropy.

FIG. 1 .

 1 FIG. 1. (a) Time dependence of the mean squared displacement ∆(t) for different concentrations of pinned particles, c, at T = 0.45 for N = 1200. The horizontal dashed line is the height of the plateau for c = 0.2. (b) The dynamic overlap function Q(t) at T = 0.45 for N = 1200 and concentrations c at which the system is in the equilibrium glass state. The horizontal dashed arrows are the static overlap functions Q (static) evaluated from the parallel tempering simulations. The inset shows the state points for the presented Q(t) in the phase diagram from Ref. [10].

FIG. 2 .

 2 FIG. 2. The mean-squared displacement ∆α at T = 0.45 for N = 1200 used for the Frenkel-Ladd thermodynamic integration. The horizontal dashed line is the plateau height for c = 0.2 shown in Fig. 1(a). The solid line corresponds to the behavior of the Einstein solid, ∆α = 3/(2α).

FLFIG. 3 .

 3 FIG. 3. (a) c-dependence of the total entropy stot, as well as different estimates of the vibrational entropy. (b) c-dependence of stots vib . The open and filled symbols are for N = 300 and N = 1200, respectively. The horizontal arrows locate TK(c) where the skewness of the overlap distribution becomes zero [10].

FIG. 4 .

 4 FIG. 4. (a-c): eIS(t) at T = 0.45 and c = 0.2 for N = 300 samples with (a) 10 4 and (b) 10 8 MC time step, and (c) 2 × 10 9 PT time step. Different lines correspond to different samples. (d-f): Superposition of the corresponding snapshots (sample 1) of the IS configurations with (d) 10 4 and (e) 10 8 MC time step, and (f) 2 × 10 9 PT time step. Pinned particles are shown in red, and mobile A and B particles are shown in blue and gray, respectively.

FIG. 5 .

 5 FIG. 5. (a): Time dependence of the mean squared displacement ∆(t) of a single sample at T = 0.45 and c = 0.167 for N = 1200. The starting configurations correspond to different times tw in the PT trajectory and are given in the legend. (b) Same as in panel (a) but now for the collective overlap function Q(t) defined in the main text.

FIG. 6 .

 6 FIG. 6. Self part of the van Hove correlation function at T = 0.45 and different values of c for N = 1200. (a)-(c): Species A, (d)-(f): Species B.

FIG. 7 .

 7 FIG. 7. Distinct part of the van Hove correlation function at T = 0.45 and different values of c for N = 1200. (a)-(c): A-A correlation, (d)-(f): B-B correlation, (g)-(i): A-B correlation.

  (a), the simulation data of U anh (c, T ) are plotted for several values of c. The solid lines are fits down to T = 0. The resulting s anh (c, T ) = S anh (c, T )/(1 -c)N are shown in Fig. 8(b).

FIG. 8 .

 8 FIG. 8. (a) Temperature dependence of the anharmonic potential energy U anh (c, T ) for several values of c's. (b) Temperature dependence of the anharmonic contribution of the vibrational entropy s anh (c, T ) = S anh (c, T )/(1 -c)N per a particle for several values of c. The data in (a) (b) are from N = 1200.

FIG. 9 .

 9 FIG. 9. (a) Distribution function P (Q) of the overlap for different concentration of pinned particles. T = 0.45 and N = 1200. (b) c-dependence of the average overlap for N = 300 and N = 1200 at T = 0.45. (c) T -dependence of the configurational entropy along c = 0.2. The open and filled symbols are for N = 300 and N = 1200, respectively.

FIG. 10 .

 10 FIG. 10. Snapshot of IS for N = 1200, T = 0.45. The size of the box is 10.0, pinned particles are shown in red, and mobile A and B particles are shown in blue and gray, respectively. (a) c = 0.0167 (fluid), (b) c = 0.08 (fluid close to the ideal glass transition), and (c) c = 0.1667 (glass).

We thank G. Biroli, C. Dasgupta, K. Hukushima, Y. Jin, S. Karmakar, and K. Kim for helpful discussions. The numerical simulations were partially performed at Research Center of Computational Science (RCCS), Okazaki, Japan. W. K. acknowledges ANR-15-

Supplemental Information

Dynamics Absence of aging