## Tetrahedrality dictates dynamics in hard spheres

Susana Mar´ ın-Aguilar, Henricus H. Wensink, Giuseppe Foffi, and Frank Smallenburg Laboratoire de Physique des Solides, CNRS, Universit´ Paris-Sud, Universit´ Paris-Saclay, 91405 Orsay, France e e

Glasses are ubiquitous amorphous solids that remain one of the big mysteries in condensed matter. Despite the vast body of literature on glasses, a unifying approach to link the structure and dynamics of glasses is still missing[1, 2]. A growing set of evidence, however, indicates the microscopic local geometry as a key ingredient[3, 4]. This originated from the seminal work of Frank [5], who conjectured that glasses may be the result of the local tendency of liquids to form icosahedral structures, which are not capable of globally filling space regularly. Here, we show that, for a fundamental glass model, dynamics can be fully understood by simply counting the number of tetrahedra. Both globally and locally, these local structures directly predict dynamical slowdown. After more than 60 years of Frank's Conjecture, it might not be the icosahedra that matter for glasses, but rather the tetrahedra inside them.

for glass theory [16-21]. As such, hard-sphere mixtures are an ideal starting point for testing new ways to link structure and glassy dynamics. Here, we quantify local order via the tetrahedrality of the local structure (TLS), a simple idea which consists of counting the number of tetrahedral clusters around each particle. Various order parameters based on local (poly)tetrahedral order have previously been shown to correlate with glassy dynamics in soft spheres [22], granular systems [23], and hard spheres [24, 25]. Our results show that TLS not only has an impressively strong correlation with local dynamics, but is also able to quantitatively predict the diffusivity of a vast range of dense hard-sphere mixtures. As perfect icosahedra can be broken down into tetrahedra, TLS is a natural extension of Frank's idea, but focuses on the most elementary building block of the three-dimensional fluid. Hence, we argue that tetrahedra are the slow structures that are ultimately responsible for dynamic arrest.

The elusiveness of the glass transition arises from the lack of an accepted theoretical framework. In general, two main approaches have emerged [4]. On a global level, different theoretical approaches have been unified under the random first order phase transition (RFOT) paradigm [6, 7]. This is essentially built on the original idea of Adam and Gibbs [8] who hypothesized the existence of large cooperative rearranging regions that grow in size when the system gets deeper into the glassy regime. On a local geometrical level, Frank's conjecture [5] has been a great source of inspiration and different observables have been proposed to try to pin down the relation between local geometry and the extremely long timescales observed in glasses. In particular, several authors have looked at the role of locally favored clusters with perfect or modified icosahedral structures as well as, more recently, at the local packing capability [4]. For comprehensive comparisons of most of these methods, we suggest the recent reviews in Refs. 3 and Ref. 4. In general, the emerging picture is that icosahedra and closely related clusters are long-lived structures that correlate strongly with the slow dynamics in glasses.

In this Letter, we examine the link between local structure and dynamics using simulations of hard-sphere mixtures. Hard spheres are a prototypical model glass former [9-13], and in mixtures their dynamics can be tuned not only via the packing fraction, but also by varying the size ratio and composition [14, 15]. Moreover, hard spheres can be perfectly reproduced in colloidal experiments, which have become a standard testing ground

We begin our study by exploring the structure and dynamics of dense binary hard-sphere mixtures, varying the size ratio q = σ S /σ L (with σ S L ( ) the diameter of the small (large) spheres), and the composition x L , which denotes the fraction of large spheres. For each choice of q and x L , we prepare the system at a packing fraction of η = 0 575, allow it to equilibrate and measure the dif-. fusion time τ D , defined as the time it takes a particle to diffuse over a distance σ L . Note that we only consider systems which avoid crystallization; as a result, we cannot probe compositions close to 0 or 1, or size ratios close to q = 1.

The dynamics of binary hard-sphere mixtures are known to vary significantly upon changing the composition x L and size ratio q [14, 15]. Our systematic investigation confirms this observation, as shown in Fig. 1a where we plot the diffusion time τ D as a function of x L for different size ratios q . Intriguingly, even though the packing fractions of all systems are the same, τ D varies by more than an order of magnitude. Moreover, the dependence of τ D on x L is non-monotonous, and shows qualitatively different behavior for different size ratios. For mixtures with a small size ratio q &lt; ∼ 0 75, . the diffusion time is a convex function of x L , showing a clear single minimum. In contrast, for higher size ratios, a maximum in the diffusion times appears, corresponding to some of the slowest systems we investigated.

The composition and size ratio of a hard-sphere mixture control the geometry of local packings in the fluid. Hence, the strong variation in diffusivity in Fig. 1a hints at strong variations in the local structure of these different systems as well. As local icosahedral order is known

FIG. 1: a) Diffusion time τ D as a function of composition x L for binary hard-sphere mixtures with various size ratios q as indicated. τ is the time unit of the simulation. The star indicates the point corresponding to the analysis in Fig. 3 b) Fraction of particles inside an icosahedral cluster for the same systems. The inset shows a typical icosahedral cluster. c) Average number of tetrahedra per particle. The inset shows a typical tetrahedral cluster.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00425v1-figure-884b4f23-f535-4c72-a55b-e7d4ab5fbe0b.jpg)

to slow down dynamics [3, 5, 26], we measure the fraction of particles involved in at least one local icosahedral cluster (as determined via the Topological Cluster Classification algorithm [27]). As shown in Fig. 1b, the fraction of particles in an icosahedral environment indeed shows behavior similar to that of the diffusion time. In particular, the slowest systems in our study (those with q glyph[similarequal] 0 8 and . x L glyph[similarequal] 0 3, marked by a star in Fig. . 1a) have the highest degree of icosahedral order. However, Fig. 1b clearly does not reproduce all of the behavior of the diffusion time in Fig. 1a, and hence it is unlikely that icosahedra can be seen as the only slow structural motif in the system. Indeed, recent work has shown that a

FIG. 2: Diffusion time for all investigated hard-sphere mixtures as a function of tetrahedrality. Different colors indicate different packing fractions. Within each packing fraction, open symbols correspond to binary mixtures with different size ratios and compositions. Closed symbols are polydisperse systems with different packing fractions. The dashed lines are exponential fits to the binary data for each packing fraction.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00425v1-figure-9dea05c7-110d-453b-92ae-3d7e47673e81.jpg)

variety of complex clusters in hard-sphere mixtures can have long lifetimes [28], and searching for these motifs in our systems reveals a whole family of clusters that commonly appear together with perfect icosahedral clusters (see SI). Intuitively, it is therefore likely that each of these contributes, to some degree, to the slowdown of the system.

Interestingly, one feature shared by all of these local structural motifs is their incorporation of a large number of smaller tetrahedral clusters: groups of four particles in which each pair are considered nearest neighbors, based on a modified Voronoi construction (see Methods). Hence, a natural question to ask is whether the overall tetrahedrality of the local structure (TLS) can predict dynamical behavior. As essentially all particles are involved in multiple tetrahedral clusters, we quantify tetrahedrality by measuring the average number n tet of tetrahedra a particle is involved in (see Methods), and plot the results in Fig. 1c for each binary mixture. This simple structural order parameter captures the behavior of the diffusion time almost perfectly, reproducing both the convexity of τ D for low q and its maximum at high q .

The strong correlation between tetrahedrality and diffusivity becomes even more clear when we plot the relationship between these two quantities directly. In Fig. 2,

we plot τ D vs 〈 n tet 〉 for all investigated size ratios and compositions, and for a range of packing fractions. For each packing fraction, n tet provides an excellent predictor of the diffusion time, revealing an approximately exponential relationship between τ D and 〈 n tet 〉 (dashed lines). Note that just the set of blue triangles in Fig. 2 covers all 60 different systems in Fig. 1, with vastly different size ratios and compositions. Moreover, this data collapse is not restricted to binary systems: polydisperse systems fall on the same lines for all investigated polydispersities, ranging from 1% to 20% (closed symbols in Fig. 2). This strongly suggests that this behavior is universal for all mixtures of hard spheres, with the possible exception of extreme size ratios where depletion effects could start to play a role [29].

Thus far we have examined the relationship between globally averaged TLS and diffusivity. We now turn our attention to the impact of tetrahedral clusters on dynamics at a local level. To this end, we explore the relationship between the number of tetrahedra n tet ( ) i a given particle i is involved in, and the absolute distance δr i over which it moves in a given time interval δt . In Fig. 3a and b, we show a typical snapshot of one of our slowest systems, with particles colored by either their tetrahedrality (a) or their displacement after a time interval δt = 500 τ in a typical trajectory (b). There is clear evidence of correlation between the two quantities, with regions of low tetrahedrality matching regions of high mobility (red), and vice versa (blue). However, examining the displacement in a specific trajectory provides only a limited view of particle mobility. After all, in a given trajectory, the ability of a particle to move depends not only on its environment, but also on the initial velocities of all particles. To average out this thermal noise, we measure the dynamic propensity D i of a particle: its average absolute displacement, taken over an ensemble of simulations starting from the same initial configuration [22, 30]. In Fig. 3c, we color each particle according to D i (again taken at δt = 500 τ ), and indeed reveal a striking correlation between n tet and D i . Intriguingly, both quantities show heterogeneity over approximately the same length scale. Clearly, local tetrahedrality is an excellent predictor for the dynamics of a particle in the near future.

This correlation can be quantified explicitly by calculating the Spearman's rank-order correlation between n tet ( ) i and D i . In Fig. 3b, we plot this correlation as a function of the time interval δt for both the small and large particles. At very short timescales ( δt &lt; ∼ 0 1 . τ ), before any particles escape their cages, there is little correlation between local TLS and dynamic propensity. At intermediate timescales, we find a strong negative correlation between n tet ( ) i and D i , confirming that particles involved in fewer tetrahedra have higher mobility. The level of correlation attained by TLS outperforms most of the local observables investigated so far [4]. Finally, for

FIG. 3: a,b,c) Snapshot of a glassy system at packing fraction η = 0 575, size ratio . q = 0 8 and composition . x L = 0 3 (de-. noted by a star in Fig. 1), with particles colored according to different criteria. a) According to the number of tetrahedra n tet ( ) i a particle is involved in, with red particles involved in fewer tetrahedra, and blue particles in more. b) According to the absolute displacement δr i after a time interval δt = 500 τ in a random trajectory, with red indicating fast particles and blue indicating slow ones. c) According to the dynamic propensity D i over the same time interval. In all snapshots, the color gray indicates the average for both large and small particles. c) Spearman's rank correlation between the number of tetrahedral clusters a particle is involved in n tet and either its displacement δr i (open symbols) or its dynamic propensity D i (closed symbols), measured over a time interval δt .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00425v1-figure-7d9493c9-b4ab-456b-b067-cbc307a4e96e.jpg)

timescales longer than the diffusion time τ D ≈ 10 , the 4 memory of the initial configuration is lost, and the correlations start decaying back to zero. An examination of these correlations in other hard-sphere mixtures provides similar results, with weaker correlations for systems with higher diffusivity (see SI). In faster systems, dynamics are less heterogeneous, and hence less predictable based on local structure.

The results in this Letter paint an elegant and comprehensive picture of the relationship between local structure and dynamics in hard-sphere mixtures, both on the global and local level. On the local level, examining tetra-

hedrality reveals large-scale regions of particles involved in a large number of tetrahedra, corresponding directly to areas of low mobility. Globally, the tetrahedrality of a hard-sphere mixture directly predicts its diffusivity, resulting in a data collapse of a vast variety of hard-sphere mixtures onto the same exponential curve. These predictions can be directly tested in experimental realizations of colloidal hard-sphere mixtures, and are strongly reminiscent of theoretical descriptions of glasses in terms of activation energies for collective rearrangement [8], and random first-order transition theory [1, 6, 7]. Most importantly, our results demonstrate that the dynamics of hard-sphere mixtures can be understood purely by looking for the most fundamental three-dimensional building block for the fluid: the simple tetrahedron.

## METHODS

We perform standard EDMD simulations [31] of mixtures of hard spheres in the microcanonical ensemble (i.e. at constant number of particles N , volume V , and energy E ). All particles have identical mass m . The time unit of the simulation is given by τ = √ βmσ 2 , with σ the diameter of the large spheres (for binary mixtures), or the average sphere size (for polydisperse mixtures), and β = 1 /k B T with k B Boltzmann's constant and T the temperature.

For the binary mixtures, we vary the size ratio q = σ S /σ L (with σ S L ( ) the diameter of the small (large) spheres), as well as the fraction of large spheres x L and the packing fraction η . In particular, we explore size ratios 0 6 . ≤ q ≤ 0 85, . compositions 0 2 . ≤ x L ≤ 0 65, and packing fractions 0 5 . . ≤ η ≤ 0 58. . Systems which crystallized were excluded from all analysis. All systems contained at least N = 700 particles, and measurements of local correlations in the slower systems were done using N = 2000 particles. For the polydisperse systems, we used mixtures of spheres of 15 different equally spaced diameters, with the number of spheres of each diameter chosen from a Gaussian distribution with standard deviation 0 08 . ≤ s/σ ≤ 0 20. .

In order to characterize the structure of the system, we find all tetrahedral clusters, defined as groups of four particles which are all considered nearest neighbors of each other. Nearest neighbors are determined via the same modified Voronoi construction used in the Topological Cluster Classification (TCC) algorithm [27], which differs from standard Voronoi constructions in two ways. First, two particles are only considered nearest neighbors if the line connecting them passes through the facet of the Voronoi cells shared by the two particles. Second, a cutoff parameter f c = 0 82 is used to . improve detection of four-membered rings of particles. Both modifications result in fewer connections between particles that are further apart than the typical first neighbor shell. For more details, see Ref. [27]. Tetrahedra are then detected directly from the set of nearest-neighbor bonds, and each particle i is labeled with the number n tet ( ) of tetrahedra it is involved in. i The global number of tetrahedra per particle is then defined as

$$\langle n _ { \text{tet} } \rangle = \frac { 1 } { N } \sum _ { i } n _ { \text{tet} } ( i ). \quad \quad ( 1 ) \quad [ 1 4$$

In addition to locating tetrahedral clusters, we also measure the fraction of particles involved in icosahedral clusters in our system. For this, we use the TCC algorithm [27] directly, again with a cutoff parameter f c = 0 82. .

Global dynamics are characterized via the diffusion time τ D . For binary mixtures, τ D = σ 2 L /D bin , with D bin the long-time selfdiffusion coefficient averaged over all particles. For the polydisperse systems, τ D = σ /D 2 poly , with D poly the long-time self-diffusion coefficient of particles with a diameter equal to the average diameter

- σ . Note that averaging over the diffusion of all species results in slightly lower diffusion times, due to the strong contribution from the highly mobile smallest species.

To characterize the dynamics on a per-particle basis, we measure the absolute displacement δr i ( δt ) for each particle i after a time interval δt . As this quantity fluctuates heavily between different runs starting from the same configuration, we also make use of the dynamic propensity D δt i ( ), which is defined as the average value of δr i over many ( ∼ 200) runs starting from the same configuration but with randomized particle velocities. To quantify the correlation between n tet ( ) and either i δr i or D i , we use Spearman's rank correlation coefficient, which measures how well the relation between two sets of data can be described by a monotonically increasing or decreasing function.

- [1] Berthier, L. &amp; Biroli, G. Theoretical perspective on the glass transition and amorphous materials. Rev. Mod. Phys. 83 , 587 (2011).
- [2] Stillinger, F. H. &amp; Debenedetti, P. G. Glass transition thermodynamics and kinetics. Annu. Rev. Condens. Matter Phys. 4 , 263-285 (2013).
- [3] Royall, C. P. &amp; Williams, S. R. The role of local structure in dynamical arrest. Phys. Rep. 560 , 1-75 (2015).
- [4] Tanaka, H., Tong, H., Shi, R. &amp; Russo, J. Revealing key structural features hidden in liquids and glasses. Nat. Rev. Phys. 1 , 333-348 (2019).
- [5] Frank, F. C. Supercooling of liquids. Proc. Royal Soc. A 215 , 43-46 (1952).
- [6] Kirkpatrick, T. R., Thirumalai, D. &amp; Wolynes, P. G. Scaling concepts for the dynamics of viscous liquids near an ideal glassy state. Phys. Rev. A 40 , 1045 (1989).
- [7] Bouchaud, J.-P. &amp; Biroli, G. On the adam-gibbskirkpatrick-thirumalai-wolynes scenario for the viscosity increase in glasses. J. Chem. Phys. 121 , 7347-7354 (2004).
- [8] Adam, G. &amp; Gibbs, J. H. On the temperature dependence of cooperative relaxation properties in glassforming liquids. J. Chem. Phys. 43 , 139-146 (1965).
- [9] Bengtzelius, U., Gotze, W. &amp; Sjolander, A. Dynamics of supercooled liquids and the glass transition. J. Phys. C 17 , 5915 (1984).
- [10] Barrat, J., Gotze, W. &amp; Latz, A. The liquid-glass transition of the hard-sphere system. J. Phys. Condens. Matter 1 , 7163 (1989).
- [11] Pusey, P. et al. Hard spheres: crystallization and glass formation. Phil. Trans. R. Soc. A 367 , 4993-5011 (2009).
- [12] Parisi, G. &amp; Zamponi, F. Mean-field theory of hard sphere glasses and jamming. Rev. Mod. Phys. 82 , 789 (2010).
- [13] Sanz, E. et al. Avalanches mediate crystallization in a hard-sphere glass. Proc. Natl. Acad. Sci. USA 111 , 7580 (2014).
- [14] G¨ otze, W. &amp; Voigtmann, T. Effect of composition changes on the structural relaxation of a binary mixture. Phys. Rev. E 67 , 021502 (2003).
- [15] Foffi, G., G¨tze, W., Sciortino, F., Tartaglia, P. &amp; Voigto mann, T. Mixing effects for the structural relaxation in binary hard-sphere liquids. Phys. Rev. Lett. 91 , 085701 (2003).
- [16] Pusey, P. N. &amp; van Megen, W. Observation of a glass transition in suspensions of spherical colloidal particles. Phys. Rev. Lett. 59 , 2083 (1987).

- [17] Van Megen, W. &amp; Underwood, S. Glass transition in colloidal hard spheres: Mode-coupling theory analysis. Phys. Rev. Lett. 70 , 2766 (1993).
- [18] Weeks, E. R., Crocker, J. C., Levitt, A. C., Schofield, A. &amp; Weitz, D. A. Three-dimensional direct imaging of structural relaxation near the colloidal glass transition. Science 287 , 627-631 (2000).
- [19] Berthier, L. et al. Direct experimental evidence of a growing length scale accompanying the glass transition. Science 310 , 1797-1800 (2005).
- [20] Brambilla, G. et al. Probing the equilibrium dynamics of colloidal hard spheres above the mode-coupling glass transition. Phys. Rev. Lett. 102 , 085703 (2009).
- [21] Kegel, W. K. &amp; van Blaaderen, A. Direct observation of dynamical heterogeneities in colloidal hard-sphere suspensions. Science 287 , 290-293 (2000).
- [22] Tong, H. &amp; Tanaka, H. Revealing hidden structural order controlling both fast and slow glassy dynamics in supercooled liquids. Phys. Rev. X 8 , 011041 (2018).
- [23] Xia, C. et al. The structural origin of the hard-sphere glass transition in granular packing. Nat. Commun. 6 , 8409 (2015).
- [24] Anikeenko, A. V. &amp; Medvedev, N. N. Polytetrahedral nature of the dense disordered packings of hard spheres. Phys. Rev. Lett. 98 , 235504 (2007).
- [25] Charbonneau, B., Charbonneau, P. &amp; Tarjus, G. Geometrical frustration and static correlations in a simple glass former. Phys. Rev. Lett. 108 , 035701 (2012).
- [26] Mar´ ın-Aguilar, S., Wensink, H. H., Foffi, G. &amp; Smallenburg, F. Designing slower glasses by manipulating their local structure. arXiv preprint arXiv:1812.00764 (2018).
- [27] Malins, A., Williams, S. R., Eggers, J. &amp; Royall, C. P. Identification of structure in condensed matter with the topological cluster classification. J. Chem. Phys. 139 , 234506 (2013).
- [28] Malins, A., Eggers, J., Royall, C. P., Williams, S. R. &amp; Tanaka, H. Identification of long-lived clusters and their link to slow dynamics in a model glass former. J. Chem. Phys. 138 , 12A535 (2013).
- [29] Anderson, V. J. &amp; Lekkerkerker, H. N. Insights into phase transition kinetics from colloid science. Nature 416 , 811 (2002).
- [30] Widmer-Cooper, A., Harrowell, P. &amp; Fynewever, H. How reproducible are dynamic heterogeneities in a supercooled liquid? Phys. Rev. Lett. 93 , 135701 (2004).
- [31] Rapaport, D. C. The art of molecular dynamics simulation (Cambridge university press, 2004).

## Acknowledgements

The authors thank Robert Botet, Paddy Royall, and Laura Filion for useful discussions. S.M.A. acknowledges CONACyT for funding.
Susana Marín-Aguilar Henricus H. Wensink Giuseppe Foffi Frank Smallenburg
Tetrahedrality dictates dynamics in hard spheres
Glasses are ubiquitous amorphous solids that remain one of the big mysteries
in condensed matter. Despite the vast body of literature on glasses, a unifying
approach to link the structure and dynamics of glasses is still missing. A
growing set of evidence, however, indicates the microscopic local geometry as a
key ingredient. This originated from the seminal work of Frank, who conjectured
that glasses may be the result of the local tendency of liquids to form
icosahedral structures, which are not capable of globally filling space
regularly. Here, we show that, for a fundamental glass model, dynamics can be
fully understood by simply counting the number of tetrahedra. Both globally and
locally, these local structures directly predict dynamical slowdown. After more
than 60 years of Frank's Conjecture, it might not be the icosahedra that matter
for glasses, but rather the tetrahedra inside them.