<figure><img src="image_1.png" /></figure>

HAL
open science

# Zero-temperature glass transition in two dimensions

Ludovic Berthier, Patrick Charbonneau, Andrea Saverio Ninarello, Misaki Ozawa, Sho Yaida

<figure><img src="image_2.png" /></figure>

## ► To cite this version:

Ludovic Berthier, Patrick Charbonneau, Andrea Saverio Ninarello, Misaki Ozawa, Sho Yaida. Zero-temperature glass transition in two dimensions. Nature Communications, 2019, 10, pp.1508. 10.1038/s41467-019-09512-3. hal-02101358

HAL Id: hal-02101358

https://hal.science/hal-02101358

Submitted on 10 Jun 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_3.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
ARTICLE

https://doi.org/10.1038/s41467-019-09512-3

OPEN

# Zero-temperature glass transition in two dimensions

Ludovic Berthier<fnref n="1" />, Patrick Charbonneau²,³, Andrea Ninarello⁴, Misaki Ozawa<fnref n="1" /> & Sho Yaida<fnref n="5" />

Liquids cooled towards the glass transition temperature transform into amorphous solids that have a wide range of applications. While the nature of this transformation is understood rigorously in the mean-field limit of infinite spatial dimensions, the problem remains wide open in physical dimensions. Nontrivial finite-dimensional fluctuations are hard to control analytically, and experiments fail to provide conclusive evidence regarding the nature of the glass transition. Here, we develop Monte Carlo methods for two-dimensional glass-forming liquids that allow us to access equilibrium states at sufficiently low temperatures to directly probe the glass transition in a regime inaccessible to experiments. We find that the liquid state terminates at a thermodynamic glass transition which occurs at zero temperature and is associated with an entropy crisis and a diverging static correlation length. Our results thus demonstrate that a thermodynamic glass transition can occur in finite dimensional glass-formers.

<footnote n="1">¹Laboratoire Charles Coulomb (L2C), University of Montpellier, CNRS, F-34095 Montpellier, France. ²Department of Chemistry, Duke University, Durham, NC 27708, USA. ³Department of Physics, Duke University, Durham, NC 27708, USA. ⁴CNR-ISC, UOS Sapienza, Piazzale A. Moro 2, IT-00185 Roma, Italy.</footnote>

<footnote n="5">⁵Facebook AI Research, Facebook Inc., Menlo Park, CA 94025, USA. Correspondence and requests for materials should be addressed to L.B. (email: ludovic.berthier@umontpellier.fr)</footnote>
Difficult scientific problems can drastically simplify in some unphysical limits. For instance, very large dimensions ($d \to \infty$), where *d* is the spatial dimensions) give relevant fluctuations a simple mean-field character¹, and one-dimensional ($d=1$) models can often be treated exactly. Yet these two solvable limits are crude idealizations of our three-dimensional reality. The rich theoretical arsenal developed to interpolate between them has revealed the highly nontrivial role of spatial fluctuations in all areas of science. In particular, as the number of spatial dimensions decreases, a phase transition may change nature or even disappear. Dimensionality thus provides a key tool for understanding the essence of many natural phenomena.

The glass transition from a viscous liquid to an amorphous solid is no exception². Its mean-field description, which becomes mathematically exact as $d \to \infty$, explains the dramatic slowdown of glass-forming liquids through the rarefaction of the number of glassy metastable states upon approaching a critical temperature, $T_K^{3,4}$. The configurational entropy, $s_{conf}$, which is the logarithm of the number of such states, becomes subextensive when $T \le T_K$. The equilibrium glass transition thus corresponds to an entropy crisis, a hypothesis first suggested by Kauzmann in his visionary analysis of experimental data⁵ and initially formalized by Gibbs and DiMarzio⁶ in the context of a lattice polymer model.

The broad discussion that has since ensued² has notably tried to describe the role of finite-*d* fluctuations beyond the mean-field framework$^{7-12}$, relating in particular the vanishing of $s_{conf}$ to a diverging point-to-set correlation length, the key quantity for characterizing nonperturbative fluctuations in glass formers$^{13}$. These fluctuations, however, make it difficult to examine finite- dimensional glass formers analytically, even for simple models composed of point-particles such as those we study here. Exploring a broader diversity of models, from polymer$^{14}$ to anisotropic pat- chy$^{15}$ models, may yet provide additional theoretical insight.

Meanwhile, Kauzmann's intuition has been repeatedly vali- dated by experiments$^{16,17}$, but the conceptual and technical limits of his results have not been lifted. Current experiments access essentially the same restricted temperature range as his 70-year old work. Theory and experiments thus currently fail to assess the status of the Kauzmann transition in finite *d*, or whether new mechanisms qualitatively change the underlying physics$^{18,19}$. Experimentally, it thus remains controversial whether the trend discovered by Kauzmann survives at much lower temperatures; entropy could go smoothly to zero$^{20,21}$, or to a finite residual value as temperature vanishes$^{15,22,23}$.

In this context, computer simulations are especially valuable. They allow direct measurements of both the configurational entropy and the point-to-set correlation length for realistic models of finite-dimensional glass formers². The recent devel- opment of the swap Monte Carlo algorithm (SWAP) further allows the exploration of a temperature regime that experiments cannot easily access²⁴, even using ultrastable glassy materials²⁵. This has consolidated and extended Kauzmann's experimental findings for three-dimensional glass formers²⁶. Here, we report that SWAP is so efficient in *d* = 2 that it provides access to a temperature regime equivalent to experimental timescales 10¹⁸ larger than the age of the universe. This remarkable advance gives very strong evidence of a thermodynamic glass transition at $T_K = 0$ for *d* = 2, accompanied by an entropy crisis and the divergence of the point-to-set correlation length. Our results thus illuminate the low-dimensional fate of the glass transition and shed light on the nature of glassy dynamics in *d* = 2$^{27-30}$.

## Results

**Model and macroscopic behavior.** We study a two-dimensional mixture of soft particles interacting with a $1/r^{12}$ purely repulsive

power-law pair potential and a size polydispersity chosen to minimize demixing, fractionation, and crystallization (see Methods). The average particle diameter is used as unit length, and the strength of the interaction potential as unit temperature. SWAP is implemented following the methodology recently validated for $d = 3^{24}$. Systems ranging from $N = 300$ to $N = 20,000$ particles within a periodic box are used to carefully track finite-size effects in both dynamics and thermodynamics. We mainly present results of $N = 1000$. Whereas experimental systems are typically composed of more complex particles (such as large molecules or polymers), the exact mean-field theory has thus far only been developed for the same type of point particles as we simulate here. In addition, such models have become a standard to study fundamental aspects of the glass transition, and are good representations of colloidal glasses.

Figure 1a shows that the static structure factor $S(k)$ evolves smoothly over a broad temperature range, from the onset temperature $T_{onset} = 0.250$ down to $T = 0.026$, which is the lowest temperature for which our strict equilibrium criteria are met. The typical low-temperature configuration depicted in Fig. 1b shows that particles of different sizes are well mixed, and that local ordering is extremely weak. In fact, no crystallization event was ever observed in our simulations, and the correlation lengths extracted from the pair correlation function

<figure><img src="image_4.png" /></figure>

**Fig. 1** Statics and dynamics of the *d* = 2 glass former. **a** The smooth evolution of the static structure factor from $T_{onset}$ down to the lowest studied temperature *T* = 0.026 indicates that the system remains fully amorphous at all *T*. **b** Snapshot of an equilibrium configuration at *T* = 0.026. **c** Arrhenius representation of the structural relaxation time $\tau_\alpha$ using SWAP and normal Monte Carlo dynamics, rescaled by the relaxation time at the onset temperature. The mode-coupling temperature, $T_{MCT}$ (gray dashed line), and the estimated range of experimental glass temperature, $T_g$ (navy strip), are indicated. The Arrhenius fit to the low-*T* data provides a lower bound for the growth of $\tau_\alpha$. SWAP can equilibrate systems down to $T \approx 0.3T_g$, where the Arrhenius fit gives $\tau_\alpha^{\text{normal}}/\tau_0 \sim 10^{46}$
for translational and bond-orientational orders evolve modestly with *T* (see Supplementary Note 1). In other words, the model is an excellent glass former.

The bulk dynamics and equilibration are captured by the bond-orientational order time correlation, $C_ψ(t)$, which is not affected by long-time tails observed in simple two-dimensional fluids$^{27,31}$. The $1/e$ decay of $C_ψ(t)$ robustly defines bulk relaxation timescales $\tau_\alpha$ both for SWAP ($\tau_\alpha^{\text{SWAP}}$) and normal ($\tau_\alpha^{\text{normal}}$) Monte Carlo dynamics (Fig. 1c). We normalize these timescales by $\tau_0 \equiv \tau_\alpha^{\text{normal}}(T_{\text{onset}})$. In agreement with earlier works$^{27}$, we find that translational correlation functions suffer from large finite-size effects, but that subtracting long-range Mermin–Wagner translational fluctuations results in system-size independent measurements$^{28–30}$ consistent with bond-orientational dynamics (see Supplementary Fig. 1). The normal dynamics exhibits a well-known super-Arrhenius growth of $\tau_\alpha$. Fitting its temperature evolution to a power-law divergence situates the mode-coupling crossover at $T_{\text{MCT}} = 0.123$, which is roughly the lowest temperature accessible with this dynamics. Earlier work showed that thermalization can be achieved below $T_{\text{MCT}}$ using Monte Carlo simulations$^{32}$. Following ref. $^{24}$, we estimate the narrow range within which the experimental glass temperature takes place as $T_g \in [0.0738, 0.0907]$. (Henceforth we set $T_g = 0.082$.) The lower end of this interval stems from an Arrhenius fit which provides a lower bound to the true $\tau_\alpha$. By all estimates, SWAP dynamics is clearly much faster than the normal one. The speedup is about 5 orders of magnitude at $T_{\text{MCT}}$, 10 at $T_g$, and the Arrhenius lower bound suggests a formidable 42 order-of-magnitude speedup at $T = 0.026$. Using an atomistic value, $\tau_0 = 10^{-10}$ s, converts this estimate to $\tau_\alpha = 10^{36}$ s, or approximately $10^{18}$ times the age of the universe. Such a “cosmological” speedup leaves no doubt that the SWAP equilibration algorithm largely bypasses the slowdown associated with the glass transition in $d=2$.

**Configurational entropy.** This computational advance permits the study of the $d=2$ configurational entropy and its relationship to the putative entropy crisis far beyond the previous work$^{33}$. Extending earlier work on $d=3$ systems$^{26}$, we obtain independent estimates of $s_{\text{conf}}$ using state-of-the-art methodologies, see Fig. 2a. Technical details are described in Supplementary Note 2. The first estimate stems from subtracting the vibrational contribution, measured by minimizing the potential energy of the system to an inherent structure and obtaining its vibrational spectrum, from the total liquid entropy$^{34}$. This potential energy landscape (PEL) approach needs to be complemented, for polydisperse systems, with an independent measure of the mixing entropy$^{35}$. Because minor but systematic additional adjustments are then required, two sets of PEL estimates are reported in Fig. 2a. The two are quantitatively close and similarly decrease with $T$, which confirms that methodological details do not affect our results in any essential way. This approach extends $s_{\text{conf}}$ measurements from $1.5T_g$ in earlier $d=2$ simulations$^{33}$ down to a temperature five times smaller, $0.3T_g$.

Our second estimate directly measures the glass entropy by performing a thermodynamic integration from the well-controlled harmonic solid limit. This approach, which is inspired by the Frenkel–Ladd method for crystals$^{36}$, was recently adapted to polydisperse amorphous solids$^{37}$. Because it does not count the number of inherent structures but measures instead the entropy of constrained glassy states, it is also very close in spirit (although not equivalent$^{37}$) to the free-energy measurement that makes use of the Franz–Parisie potential$^{38}$. The Frenkel–Ladd estimate is smaller than the PEL ones, as expected, but exhibits a similar temperature dependence.

From the data in Fig. 2a, $s_{\text{conf}}$ seemingly vanishes close to $T_K = 0$. This behavior sharply contrasts with that of three-dimensional

glass formers, for which evidence suggests that $T_K > 0$,$^{5,16,17,26}$. The impending entropy crisis is expected to give rise to large-scale fluctuations with a growing point-to-set correlation length$^{13}$. We use the computational tools developed in refs. $^{26,39,40}$ to analyze the thermodynamic properties of liquids confined within spherical cavities of radius $R$ drawn from a reference equilibrium configuration (see Supplementary Note 3). The distribution $P(Q)$ of the core cavity overlap $Q$ among the confined equilibrium glassy configurations is then analyzed. The point-to-set correlation length, $\xi_{\text{PTS}}$, is determined from the decay with $R$ of the average overlap. This length is then transformed into a third estimate, $s_{\text{conf}} \propto \xi_{\text{PTS}}^{-(d-\theta)}$ with $\theta = 1$. In $d=2$, this choice of $\theta$ is natural because it both saturates the bound $\theta \le d - 1$ and satisfies the wetting relation $\theta = d/2$.$^{3}$ The resulting $s_{\text{conf}}(T) = \xi_{\text{PTS}}(T_g)/\xi_{\text{PTS}}(T)$ in Fig. 2a again has a similar temperature evolution as other estimates.

Figure 2b shows that rescaling all configurational entropies by their value at $T_g$ collapses the entire set of measurements. This robustness is non-trivial because all four estimates make different types of approximations. The agreement of their temperature dependence may thus resolve earlier discrepancies and debates regarding conflicting estimates of the configurational entropy$^{41,42}$.

One expects $s_{\text{conf}}$ to vanish linearly, $s_{\text{conf}} \propto (T - T_K)$, but this scaling arguably has a quadratic correction at higher temperatures. We thus perform a quadratic fit to the low-temperature regime, $T < T_g$. This fitting yields $|T_K| \le 0.003$ for all cases. These estimates of $T_K$ are 10 times smaller than our lowest temperature, $T = 0.026$, and 30 times smaller than $T_g$. The scaling behavior implied by this observation is presented in Supplementary Note 4. Known alternatives to an entropy crisis invoke a change in the concavity of $s_{\text{conf}}^{14,22,43}$ and should be accompanied by a maximum in the specific heat $c_V$,$^{18,20,21}$; we observe neither the convexity (Fig. 2a) nor the specific heat maximum (Fig. 2c). As $T \to T_K$, $c_V$ instead monotonically increases towards a finite value that is larger than the Dulong–Petit law. These observations thus strongly support the occurrence of a non-trivial entropy crisis at $T_K = 0$. The only alternative left is a change of behavior occurring at temperatures even lower than those we can study directly.

**Point-to-set length scale.** The thermodynamic glass transition at $T_K = 0$ also coincides with a divergence of the point-to-set correlation length. We illustrate the physical meaning of this length scale in Fig. 3a in the form of a $(T, 1/R)$ diagram reminiscent of both the Franz–Parisie thermodynamic construction$^{38}$ and of the random pinning approach$^{44,45}$. Upon decreasing the cavity size at a given temperature, the system crosses over from a low-Q regime at large $R$ to a high-Q regime at small $R$, as illustrated by the snapshots in Fig. 3a. For any $T > 0$, this crossover around $R \approx \xi_{\text{PTS}}$ corresponds to a finite-size version of the random first-order glass transition with a rarefaction of the number of locally available states as $R$ decreases$^{46}$. The evolution of $P(Q)$ in Fig. 3b indeed exhibits features reminiscent of phase coexistence near an incipient random first-order transition. The crossover also sharpens as $T$ decreases, suggesting that the growing correlation length transforms it into a genuine thermodynamic phase transition as $T \to T_K = 0$. In absolute values, $\xi_{\text{PTS}} \approx 6.5$ at $T = 0.028$, which represents a very large static correlation length for glassy models$^{26,39,40}$. It implies that large clusters comprising about 120 particles are statically correlated, and should thus move collectively to restructure the liquid. These results are consistent with the sharp decay of the configurational entropy in Fig. 2 and the expected dramatic increase of the relaxation time in Fig. 1.
<figure><img src="image_5.png" /><figcaption>**Fig. 2** Zero-temperature Kauzmann transition. **a** Decrease of the configurational entropy with temperature using the potential energy landscape (PEL), Frenkel-Ladd (FL), and point-to-set (PTS) length estimates. The error bars for FL correspond to the ambiguity of defining the plateau regime in the mean squared displacement of the FL construction. **b** Once rescaled by their value at $T_g$, all estimates evolve nearly identically, and the collapsed data are well fitted by a quadratic function of $T$ for $T < T_g$ (dashed blue line: $s_{conf}(T)/s_{conf}(T_g) = 0.01 + 1.48(T/T_g) - 0.49(T/T_g)^2$ indicates the quadratic fit for the point-to-set estimate). All results are consistent with a linearly vanishing $s_{conf}$ at $T_K = 0$. **c** The specific heat, $c_V$, obtained from the derivative of the potential energy increases monotonically above the Dulong-Petit law for $d = 2$ (dashed horizontal line), which is also consistent with a thermodynamic transition at $T_K = 0$</figcaption></figure>

## Discussion

The problem of the glass transition has two fundamental facets: thermodynamics and dynamics. While the current study focused on the thermodynamics of $d = 2$ glass formers, its dynamical counterpart, which involves obtaining a detailed functional form of the structural relaxation time, remains for now out of reach of computational work. Our results nonetheless suggest that in $d = 2$ the divergence of the relaxation time must take place at zero (rather than at finite) temperature. By identifying the thermodynamic properties that underlie the nature of glassy dynamics in $d = 2^{27-30}$, our results provide additional evidence that a thermodynamic transition can occur in finite-dimensional systems, and that the lower critical dimension for the long-range amorphous order is $d_L = 2$ (see Supplementary Note 5). This finding lends indirect support to previous observations in $d = 3^{26}$, and will surely guide future analytical work.

## Methods

**Model.** The glass-forming model we consider consists of particles with purely repulsive soft-sphere interactions, and a continuous size polydispersity. Particle

diameters, $\sigma_i$, are randomly drawn from a distribution of the form: $f(\sigma) = A\sigma^{-3}$, for $\sigma \in [\sigma_{min}, \sigma_{max}]$, where A is a normalization constant. The size polydispersity is quantified by $\delta = \sqrt{\sigma_i^2 - \sigma_{max}^2/\sigma}$, where $\cdot \cdot \equiv \int d\sigma f(\sigma)(\cdot \cdot)$, and is here set to $\delta = 0.23$ by imposing $\sigma_{min}/\sigma_{max} = 0.45$. The average diameter, $\bar{\sigma}$, sets the unit of length. The soft-sphere interactions are pairwise and described by an inverse power-law potential

$$v_{ij}(r) = v_0 \left(\frac{\sigma_{ij}}{r}\right)^{12} + c_0 + c_1 \left(\frac{r}{\sigma_{ij}}\right)^2 + c_2 \left(\frac{r}{\sigma_{ij}}\right)^4, \quad (1)$$

$$\sigma_{ij} = \frac{(\sigma_i + \sigma_j)}{2} (1 - \varepsilon |\sigma_i - \sigma_j|), \quad (2)$$

where $v_0$ sets the unit of energy (and temperature with Boltzmann constant $k_B = 1$), and $\varepsilon = 0.2$ quantifies the degree of non-additivity of particle diameters. We introduce $\varepsilon > 0$ to the model in order to suppress fractionation and thus enhance glass form ability$^{24,47}$. The constants, $c_0$, $c_1$, and $c_2$, enforce a vanishing potential and the continuity of the first and second derivatives of the potential at the cut-off distance $r_{cut} = 1.25\sigma_{ij}$. We simulate a system with N particles within a square cell of area V under periodic boundary conditions, at number density $\rho = N/V = 1.01$. Most simulations have N = 1000, but systems with N = 300, 3000, 8000, and 20,000 are also studied.
<figure><img src="image_6.png" /><figcaption>**Fig. 3** Approaching the random first-order transition. **a** Phase diagram showing the low-Q region for large cavities and high-Q region for small cavities, separated by the boundary determined by the point-to-set correlation length, $\xi_{PTS}$. The dashed blue line is the same quadratic fit (after unit conversion) as in Fig. 2b. Inset: Representative configurations with overlap field for $T = 0.035$ at $R = 6.6$ (low Q, white) and 4.8 (high Q, dark). **b** Evolution of the probability distribution function of overlap $P(Q)$ at $T = 0.035$ from $R = 4.8$ to $R = 6.6$. Bimodality signals a first-order-like phase coexistence</figcaption></figure>

**Observables.** We monitor the system structure with two common liquid state quantities: the pair-distribution function $g(r)$, and the structure factor $S(k) = \langle p_{-k}p_k \rangle / N$, where $p_k = \sum_i e^{ik}r_i$ is the Fourier-space density. Orientational correlations are also considered, and are quantified using the six-fold bond-orientational order parameter [?]48

$$ \psi_6 = \frac{1}{N} \sum_{j=1}^{N} \psi_6^j \quad \text{where} \quad \psi_6^j = \frac{1}{n_j} \sum_{k=1}^{n_j} \exp(i6\theta_{jk}), \quad (3) $$

where the sum is performed over the $n_j$ first neighbors of the j-particle. These neighbors are defined as particles with $r_{ij}/\sigma_{ij} < 1.33$, which is the location of the distance of the first minimum in the rescaled radial distribution function $g(r/\sigma_{ij})$. The angle $\theta_{jk}$ then measures the orientation of the axis between the two particles with respect to the x-axis. Because these correlations are orientationally invariant the choice of x-axis is made without loss of generality. Orientational correlations are then monitored through the two-point bond-orientational correlation function

$$ g_6(r) = \langle \psi_6(r) \psi_6^*(0) \rangle, \quad (4) $$

where $\psi_6(r) = \sum_{i=1}^{N} \delta(|r - r_i|)\psi_6^i$. The radial decay of the hexatic order correlation function, $g_6(r)/g(r)^{48}$, provides an hexatic correlation length $\xi_6$, as presented in Supplementary Note 1.

Translational dynamics is characterized by first measuring the intermediate scattering function

$$ F_s(k, t) = \frac{1}{N} \left\langle \sum_{j=1}^{N} \exp \left[ ik \cdot (r_j(t) - r_j(0)) \right] \right\rangle \quad (5) $$

at the wave number $k$ corresponding to the first peak of $S(k)$. The relaxation time of the density fluctuations, $\tau_{\alpha}^{TR}$, is then extracted from the exponential decay of the scattering function, i.e., $F_s(k, \tau_{\alpha}^{TR}) = e^{-1}$. Orientational dynamics is characterized similarly, replacing the Fourier-space density by the bond-orientational correlation function in Eq. (3) defined by

$$ C_{\psi_6}(t) = \frac{1}{N} \left\langle \sum_{i=0}^{N} \psi_6^i(t) [\psi_6^i(0)]^* \right\rangle \quad (6) $$

In order to extract the bond-orientational relaxation time $\tau_{\alpha}$, we use $C_{\psi_6}(\tau_{\alpha}) = e^{-1}$.

**Equilibration and the glass ceiling.** Normal Monte-Carlo (MC) simulations allow only local particle displacements, drawing a random displacement vector on the $(x, y)$ axis in the interval $[-\Delta r_{max}, \Delta r_{max}]$ with $\Delta r_{max} = 0.6$ and moving a randomly chosen particle following a Metropolis acceptance criterion. Compounding $N$ such displacement attempts defines a MC step, which is used as unit of time in this work. To ensure equilibration, we monitor both static and dynamical observables. Starting from a high-temperature liquid configuration, we quench the system at the final temperature and wait for the potential energy of the system to stop aging on a time window of ~$10^6$ MC steps. We first estimate $\tau_{\alpha}$ on simulations long enough to allow few decorrelations of $C_{\psi_6}(t)$, and then perform simulations for 220$\tau_{\alpha}$. The system is left to equilibrate during the first $20\tau_{\alpha}$, static and dynamical observables are computed over the following $200\tau_{\alpha}$. Swap MC simulations include attempts at exchanging random pairs of particle diameters, which replace particle displacements with probability $p_{swap} = 0.2$. This algorithm defines the SWAP dynamics. The same equilibration and measuring protocol as for normal MC is then followed. Static observables monitor ordering and phase separation in the system, as discussed in Supplementary Note 1, whereas dynamical observables quantify the relaxation and equilibration timescales.

In Supplementary Fig. 1, we report orientational $\tau_{\alpha}$ and translational $\tau_{\alpha}^{TR}$ relaxation times for both normal and SWAP dynamics. Because the relaxation of local orientational degrees of freedom is slower, the associated timescale is used as reference. We perform three different fits to the $\tau_{\alpha}$ results for the physical dynamics, in order to extract the temperatures relevant to the dynamical slowing down. First, we fit $\tau_{\alpha}$ to a power-law function, as is predicted in the context of the mode-coupling theory49,

$$ \tau_{\alpha} \propto (T - T_{\text{MCT}})^{-\gamma}, \quad (7) $$

over the interval $\tau_{\alpha} \in (\tau_0, 10^3\tau_0)$. The resulting $T_{MCT} = 0.123$ roughly corresponds to the lowest temperature at which normal dynamics can reach equilibrium in simulations of reasonable duration24.

Next, we estimate the laboratory glass transition temperature, $T_g$, at which experiments with atomic and molecular glass formers cannot be equilibrated anymore. At $T_g$ relaxation times have increased by 12 orders of magnitude with respect to their value at the onset of the supercooled dynamics50. We thus fit the relaxation times both to a Vogel-Fulcher-Tallman (VFT) law

$$ \tau_{\alpha} \propto \exp \left( \frac{A}{T - T^{\text{VFT}}} \right), \quad (8) $$

and to an Arrhenius law

$$ \tau_{\alpha} \propto \exp \left( \frac{B}{T} \right), \quad (9) $$

where $A$ and $B$ are fitting constants. These two expressions respectively overestimate and underestimate the increase of relaxation times in experimental glass-formers51,52. We fit Eq. (8) using the whole temperature range $T < T_{onset}$ whereas we fit Eq. (9) only to $T < 0.16$ to ensure that the result serves as a proper lower bound on the relaxation time. Extrapolating up to the temperature at which $\log_{10}(\tau_{\alpha}/\tau_0) \simeq 12$ gives $T_g^{VFT} = 0.0907$ and $T_g^{Arr} = 0.0738$. These two temperatures are, by construction, upper and lower bounds for $T_g$, and thus define an experimental glass-ceiling regime (blue shaded region)26 in Figs. 1, 2 and 3 as well as Supplementary Fig. 1. In all cases, SWAP dynamics equilibrates well beyond this experimentally limited regime, reaching $T = 0.026$. Supplementary Fig. 1 also shows the fitting curves to the dynamics. The mode-coupling power-law prediction describes the growth of the relaxation times only within the first three decades of the glassy regime, but at lower temperatures it overestimates the results by many orders of magnitude. Whereas Eq. (8) adequately describes these same results over more than four decades, an Arrhenius law captures barely two decades.

## Data availability

The data necessary to reproduce the figures in this paper are publicly available through the Duke University Libraries Digital Repository (https://doi.org/10.7924/r46w9b248)53.
Received: 5 December 2018 Accepted: 10 March 2019

Published online: 03 April 2019

## References

1. Chaikin, P. M. & Lubensky, T. C. *Principles of Condensed Matter Physics* (Cambridge University Press, Cambridge, UK 2000).

2. Berthier, L. & Biroli, G. Theoretical perspective on the glass transition and amorphous materials. *Rev. Mod. Phys.* **83**, 587 (2011).

3. Lubchenko, V. & Wolynes, P. G. Theory of structural glasses and supercooled liquids. *Annu. Rev. Phys. Chem.* **58**, 235 (2007).

4. Charbonneau, P., Kurchan, J., Parisi, G., Urbani, P. & Zamponi, F. Glass and jamming transitions: from exact results to finite-dimensional descriptions. *Annu. Rev. Condens. Matter Phys.* **8**, 265 (2017).

5. Kauzmann, W. The nature of the glassy state and the behavior of liquids at low temperatures. *Chem. Rev.* **43**, 219 (1948).

6. Gibbs, J. H. & DiMarzio, E. A. Nature of the glass transition and the glassy state. *J. Chem. Phys.* **28**, 373–383 (1958).

7. Dzero, M., Schmalian, J. & Wolynes, P. G. Activated events in glasses: the structure of entropic droplets. *Phys. Rev. B* **72**, 100201 (2005).

8. Franz, S. First steps of a nucleation theory in disordered systems. *J. Stat. Mech. Theory Exp.* **2005**, P04001 (2005).

9. Angelini, M. C. & Biroli, G. Real space Migdal–Kadanoff renormalisation of glassy systems: recent results and a critical assessment. *J. Stat. Phys.* **167**, 476 (2017).

10. Rulquin, C., Urbani, P., Biroli, G., Tarjus, G. & Tarzia, M. Nonperturbative fluctuations and metastability in a simple model: from observables to microscopic theory and back. *J. Stat. Mech. Theory Exp.* **2016**, 023209 (2016).

11. Biroli, G., Cammarota, C., Tarjus, G. & Tarzia, M. Random-field ising-like effective theory of the glass transition. I. Mean-field models. *Phys. Rev. B* **98**, 174205 (2018).

12. Biroli, G., Cammarota, C., Tarjus, G. & Tarzia, M. Random field ising-like effective theory of the glass transition. II. Finite-dimensional models. *Phys. Rev. B* **98**, 174206 (2018).

13. Bouchaud, J.-P. & Biroli, G. On the Adam–Gibbs–Kirkpatrick–Thirumalai– Wolynes scenario for the viscosity increase in glasses. *J. Chem. Phys.* **121**, 7347 (2004).

14. Dudowicz, J., Freed, K. F. & Douglas, J. F. Generalized entropy theory of polymer glass formation. *Adv. Chem. Phys.* **137**, 125–222 (2008).

15. Smallenburg, F. & Sciortino, F. Liquids more stable than crystals in particles with limited valence and flexible bonds. *Nat. Phys.* **9**, 554 (2013).

16. Richert, R. & Angell, C. A. Dynamics of glass-forming liquids. V. On the link between molecular dynamics and configurational entropy. *J. Chem. Phys.* **108**, 9016 (1998).

17. Tatsumi, S., Aso, S. & Yamamuro, O. Thermodynamic study of simple molecular glasses: universal features in their heat capacity and the size of the cooperatively rearranging regions. *Phys. Rev. Lett.* **109**, 045701 (2012).

18. Tarjus, G., Kivelson, S. A., Nussinov, Z. & Viot, P. The frustration-based approach of supercooled liquids and the glass transition: a review and critical assessment. *J. Phys. Condens. Matter* **17**, R1143 (2005).

19. Chandler, D. & Garrahan, J. P. Dynamics on the way to forming glass: bubbles in space-time. *Annu. Rev. Phys. Chem.* **61**, 191 (2010).

20. Kivelson, D. & Tarjus, G. The Kauzmann paradox interpreted via the theory of frustration-limited-domains. *J. Chem. Phys.* **109**, 5481–5486 (1998).

21. Debenedetti, P. G., Stillinger, F. H. & Shell, M. S. Model energy landscapes. *J. Phys. Chem. B* **107**, 14434 (2003).

22. Xu, W.-S., Douglas, J. F. & Freed, K. F. Generalized entropy theory of glass-formation in fully flexible polymer melts. *J. Chem. Phys.* **145**, 234509 (2016).

23. Donev, A., Stillinger, F. H. & Torquato, S. Do binary hard disks exhibit an ideal glass transition? *Phys. Rev. Lett.* **96**, 225502 (2006).

24. Ninarello, A., Berthier, L. & Coslovich, D. Models and algorithms for the next generation of glass transition studies. *Phys. Rev. X* **7**, 021039 (2017).

25. Ediger, M. D. Perspective: highly stable vapor-deposited glasses. *J. Chem. Phys.* **147**, 210901 (2017).

26. Berthier, L. et al. Configurational entropy measurements in extremely supercooled liquids that break the glass ceiling. *Proc. Natl Acad. Sci. USA* **114**, 11356 (2017).

27. Flenner, E. & Szamel, G. Fundamental differences between glassy dynamics in two and three dimensions. *Nat. Commun.* **6**, 7392 (2015).

28. Shiba, H., Yamada, Y., Kawasaki, T. & Kim, K. Unveiling dimensionality dependence of glassy dynamics: 2D infinite fluctuation eclipses inherent structural relaxation. *Phys. Rev. Lett.* **117**, 245701 (2016).

29. Vivek, S., Kelleher, C. P., Chaikin, P. M. & Weeks, E. R. Long-wavelength fluctuations and the glass transition in two dimensions and three dimensions. *Proc. Natl Acad. Sci. USA* **114**, 1850 (2017).

30. Illing, B. et al. Mermin–Wagner fluctuations in 2D amorphous solids. *Proc. Natl Acad. Sci. USA* **114**, 1856 (2017).

31. Isobe, M. & Alder, B. J. Generalized bond order parameters to characterize transient crystals. *J. Chem. Phys.* **137**, 194501 (2012).

32. Santen, L. & Krauth, W. Absence of thermodynamic phase transition in a model glass former. *Nature* **405**, 550 (2000).

33. Sengupta, S., Karmakar, S., Dasgupta, C. & Sastry, S. Adam–Gibbs relation for glass-forming liquids in two, three, and four dimensions. *Phys. Rev. Lett.* **109**, 095705 (2012).

34. Sciortino, F., Kob, W. & Tartaglia, P. Inherent structure entropy of supercooled liquids. *Phys. Rev. Lett.* **83**, 3214 (1999).

35. Ozawa, M. & Berthier, L. Does the configurational entropy of polydisperse particles exist? *J. Chem. Phys.* **146**, 014502 (2017).

36. Frenkel, D. & Ladd, A. J. C. New Monte Carlo method to compute the free energy of arbitrary solids. Application to the fcc and hcp phases of hard spheres. *J. Chem. Phys.* **81**, 3188 (1984).

37. Ozawa, M., Parisi, G. & Berthier, L. Configurational entropy of polydisperse supercooled liquids. *J. Chem. Phys.* **149**, 154501 (2018).

38. Franz, S. & Parisi, G. Phase diagram of coupled glassy systems: a mean-field study. *Phys. Rev. Lett.* **79**, 2486–2489 (1997).

39. Biroli, G., Bouchaud, J.-P., Cavagna, A., Grigera, T. S. & Verrocchio, P. Thermodynamic signature of growing amorphous order in glass-forming liquids. *Nat. Phys.* **4**, 771 (2008).

40. Berthier, L., Charbonneau, P. & Yaida, S. Efficient measurement of point-to-set correlations and overlap fluctuations in glass-forming liquids. *J. Chem. Phys.* **144**, 024501 (2016).

41. Stillinger, F. H. Supercooled liquids, glass transitions, and the Kauzmann paradox. *J. Chem. Phys.* **88**, 7818 (1988).

42. Biroli, G. & Monasson, R. From inherent structures to pure states: some simple remarks and examples. *Europhys. Lett.* **50**, 155 (2000).

43. Saika-Voivod, I., Sciortino, F. & Poole, P. H. Free energy and configurational entropy of liquid silica: fragile-to-strong crossover and polyamorphism. *Phys. Rev. E* **69**, 041503 (2004).

44. Kob, W. & Berthier, L. Probing a liquid to glass transition in equilibrium. *Phys. Rev. Lett.* **110**, 245702 (2013).

45. Cammarota, C. & Biroli, G. Ideal glass transitions by random pinning. *Proc. Natl Acad. Sci. USA* **109**, 8850–8855 (2012).

46. Biroli, G. & Bouchaud, J.-P. The random first-order transition theory of glasses: a critical assessment. In *Structural Glasses and Supercooled Liquids* Ch. 2, 31–113 (John Wiley & Sons, Ltd., Hoboken, New Jersey, USA 2012).

47. Berthier, L., Coslovich, D., Ninarello, A. & Ozawa, M. Equilibrium sampling of hard spheres up to the jamming density and beyond. *Phys. Rev. Lett.* **116**, 238002 (2016).

48. Russo, J. & Tanaka, H. Assessing the role of static length scales behind glassy dynamics in polydisperse hard disks. *Proc. Natl Acad. Sci. USA* **112**, 6920 (2015).

49. Göëtze, W. *Complex dynamics of glass-forming liquids: A mode-coupling theory*, Vol. **145**, (Oxford University Press, Oxford 2008).

50. Ediger, M. D., Angell, C. A. & Nagel, S. R. Supercooled Liquids and Glasses. *J. Phys. Chem.* **100**, 13200–13212 (1996).

51. Elmatad, Y. S., Chandler, D. & Garrahan, J. P. Corresponding states of structural glass formers. *J. Phys. Chem. B* **113**, 5563–5567 (2009).

52. Hecksher, T., Nielsen, A. I., Olsen, N. B. & Dyre, J. C. Little evidence for dynamic divergences in ultraviscous molecular liquids. *Nat. Phys.* **4**, 737–741 (2008).

53. Berthier, L., Charbonneau, P., Ninarello, A., Ozawa, M. & Yaida, S. Data and scripts from: Zero-temperature glass transition in two dimensions. *Duke Digital Repository*. https://doi.org/10.7924/r46w9b248 (2019).

## Acknowledgements

The authors thank G. Tarjus for stimulating discussions. This research was supported by a grant from the Simons Foundation (#454933, L. Berthier and #454937, P. Charbonneau). Part of the computations was carried out through the Duke Compute Cluster.

## Author contributions

L.B., P.C., A.N., M.O., and S.Y. designed research, performed research, analyzed data, and wrote the paper.
**Additional information**

**Supplementary Information** accompanies this paper at https://doi.org/10.1038/s41467-019-09512-3.

**Competing interests:** The authors declare no competing interests.

**Reprints and permission** information is available online at http://npg.nature.com/reprintsandpermissions/

**Journal peer review information:** *Nature Communications* thanks Pablo Debenedetti and the other anonymous reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_7.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2019