#
Finite-Size Scaling of a First-Order Dynamical Phase Transition:

Adaptive Population Dynamics and an Effective Model

###### Abstract

We analyze large deviations of the time-averaged activity in the one dimensional Fredrickson-Andersen model, both numerically and analytically. The model exhibits a dynamical phase transition, which appears as a singularity in the large deviation function. We analyze the finite-size scaling of this phase transition numerically, by generalizing an existing cloning algorithm to include a multi-canonical feedback control: this significantly improves the computational efficiency. Motivated by these numerical results, we formulate an effective theory for the model in the vicinity of the phase transition, which accounts quantitatively for the observed behavior. We discuss potential applications of the numerical method and the effective theory in a range of more general contexts.

###### pacs:

05.40.-a, 05.10.-a, 05.70.LnIntroduction – Systems far from equilibrium display a wide spectrum of complex behavior dorfman_introduction_1999 (); kampen_stochastic_2007 (). For example, thermodynamic phase transitions are usually forbidden in one-dimensional systems, but a variety of dynamical phase transitions are still observed hinrichsen2000 (); evans2005 (); bertini_current_2005 (); bodineau_distribution_2005 (); garrahanjacklecomtepitardvanduijvendijkvanwijland (). Such transitions can appear in far-from-equilibrium states that are defined by restricting (or conditioning) trajectories so that time-averaged observables take non-typical values Touchette20091 (). They can be related to physical properties of systems where metastability is important, especially glassy systems merolle_spacetime_2005 (); jack_space-time_2006 (); garrahanjacklecomtepitardvanduijvendijkvanwijland (); garrahan_first-order_2009 (); kurchan_six_2009 (); hedges_dynamic_2009 (); garrahan_kinetically_2010 (); chandler_dynamics_2010 (); limmer2014 ().

In some cases these transitions can be studied analytically bodineau_cumulants_2007 (); bertini_current_2005 (); bertini_non_2006 (); bodineau_distribution_2005 (); garrahanjacklecomtepitardvanduijvendijkvanwijland (), but in practical applications one must often resort to numerical methods: These access the relevant far-from-equilibrium states by rare-event sampling algorithms giardina_direct_2006 (); hedges_dynamic_2009 (); pitard_dynamic_2011 (); speck2012 (); malins2012 (), employing for instance population dynamics or path sampling. Such methods tend to perform poorly in the vicinity of dynamical phase transitions, just as conventional sampling methods tend to fail close to equilibrium phase transitions. For the equilibrium case, advanced methods exist that solve this problem, including finite-size scaling analysis privmanfisher (); borgskotecky () and multi-canonical sampling PhysRevLett.68.9 (); PhysRevLett.86.2050 (); bruce2003 (). For dynamical phase transitions, some progress has been made in this direction speck2012 (); limmer2014 (); gingrich2015 () but accurate calculations are numerically expensive and suffer from significant finite-size effects.

Here, we analyze a dynamical phase transition garrahanjacklecomtepitardvanduijvendijkvanwijland () in the Fredrickson–Andersen (FA) model fredrickson_kinetic_1984 (). We combine a state-of-the-art numerical approach NemotoBouchetLecomteJack () with a theoretical analysis. We show that numerical results and theoretical predictions for finite-size scaling near the phase transition agree quantitatively. By combining these ingredients we obtain a full description of the transition, at a modest computational cost. The phase transition is a prototype for transitions in a range of systems garrahanjacklecomtepitardvanduijvendijkvanwijland (); elmatad2010 (); hedges_dynamic_2009 (); limmer2014 (), so we argue that these new methods and insights have broad potential application in this field.

Model – The one-dimensional (1d) FA model fredrickson_kinetic_1984 () is a kinetically constrained model (KCM) that consists of spins on a periodic lattice. The th spin takes values (down) or (up) and the configuration is . We define an operator that flips the state of spin , so that . The kinetic constraint of the model is that spin can flip only if at least one of its neighbors is up. The transition rates between configurations reflect this constraint, they are:

(1) |

where enforces the kinetic constraint and is a parameter that depends on the temperature in the model fredrickson_kinetic_1984 (). The rates obey detailed balance and the model’s equilibrium distribution follows a Bernoulli law, . Despite this trivial distribution, the kinetic constraint in the model leads to rich behavior, related to dynamical heterogeneity in glassy systems doi:10.1080/0001873031000093582 (); garrahan_kinetically_2010 (); hedges_dynamic_2009 ().

Dynamical phase transitions – We define the dynamical activity where is the total number of spin flips during the time interval . The phase transitions that we consider take place in ensembles of trajectories that are restricted to a given value of the activity. In the limit , the activity converges to its equilibrium value : to estimate the probability of rare trajectories with , we consider the cumulant generating function (CGF)

(2) |

where denotes an ensemble average. Dynamical phase transitions are associated with singularities in : they are analogous to thermodynamic phase transitions, with corresponding to the thermodynamic free energy garrahan_first-order_2009 () and corresponding to an intensive thermodynamic field that is used to drive the system through its phase transition. We also define

(3) |

which specifies the dependence of the mean activity on the field , analogous to the dependence of the order parameter on its conjugate field in thermodynamics.

The dynamical phase transition that occurs in the FA model separates a high-activity state [with ] from an inactive (glass) state [with . It is defined in a joint limit of large time and large system size . In this work, we first take and then take . The phase transition is first-order, and the order parameter exhibits a discontinuous jump at garrahanjacklecomtepitardvanduijvendijkvanwijland (). However, the large- limit is not accessible numerically, and for finite the activity is a smooth function of , whose representative behavior is shown in Fig. 1. The crossover sharpens as increases: to analyze the transition, one must consider the finite-size scaling of .

Cloning algorithm with feedback – To perform this finite-size scaling, we require a numerical method that provides accurate results for a range of system sizes. To this end, we generalize a recently-proposed adaptive method NemotoBouchetLecomteJack () to Markov jump processes. The method is based on a cloning algorithm giardina_direct_2006 (); giardina_simulating_2011 () which uses a population of clones (or copies) of the system. We fix a time interval and the dynamics of the model are propagated over intervals of length , such that the total time is . For each interval, one calculates a weighting factor for clone :

(4) |

where is the number of spin flips for clone , evaluated over the whole time interval ; also one defines . After each time interval, clones are duplicated or removed, to enforce the conditioning on the activity. In this step, each clone generates a number of offspring proportional to its weight . The mean activity can then be obtained as the average of over the final population NemotoBouchetLecomteJack ().

This algorithm provides accurate results when is sufficiently large NemotoGuevaraLecomte (), but in practice this may require a very large number of clones, which is computationally expensive. To avoid this issue, we combine the existing cloning algorithm giardina_direct_2006 (); giardina_simulating_2011 (); hedges_dynamic_2009 (); pitard_dynamic_2011 (); speck2012 (); malins2012 (); NemotoGuevaraLecomte () with a modification of the dynamics jack_large_2010 (); 1742-5468-2010-10-P10007 (); PhysRevLett.111.120601 (); PhysRevLett.112.090602 (), following NemotoBouchetLecomteJack (). In order to aid sampling of trajectories with non-typical activity, we modify the transition rates of the model as

(5) |

where is an effective potential or control potential jack2015_eff (). The weight factors are also modified, by replacing in (4) with

(6) |

where is the escape rate from configuration , and is obtained in the same way but using the modified rates (5).

In the limit of large , the results of the algorithm are independent of the choice of . However, an appropriate choice can dramatically improve the accuracy of results obtained with finite populations. In particular, there exists an optimal control potential for which a population of is already sufficient for convergence. This optimal potential is given (up to an arbitrary constant) by where is the probability of observing configuration within the steady state of the cloning algorithm NemotoBouchetLecomteJack (). Calculating this optimal control directly is not feasible in practice: instead we restrict to control potentials that involve interactions between each spin and its nearest neighbors at distance , so that the change in the effective potential on flipping spin is

(7) |

for some function . To obtain the most suitable values for these potentials, we use a feedback scheme: we run the cloning algorithm, estimate the probabilities of particular local arrangements of the spins, and update the effective potential based on that choice. By repeating this procedure, one can optimize the choice of the control potential (see NemotoBouchetLecomteJack () and SM1 () for details).

Results – Fig. 1 shows the performance of the algorithm, close to the dynamical phase transition. We plot as a function of , as obtained from the cloning algorithm, using the feedback method to determine suitable effective interactions. The interaction range is , and we take between 50 and 800. We set , which is sufficiently large to converge to the large- limit. As increases, the estimates of converge to a smooth curve, indicating that these clone populations are large enough to achieve accurate results. By contrast, the inset to Fig. 1 shows that the original cloning method (with ) deviates significantly from the correct result, compared with the feedback method for the same number of copies. This tendency is observed throughout the whole range of , irrespective of the presence of the dynamical coexistence SM2 ().

The results of Fig. 1 show the expected crossover from high to low activity, consistent with the existence of a dynamical phase transition near . To analyze the finite-size scaling of this transition, we define as the maximal susceptibility

(8) |

Let the finite-size transition point be the value of at which this maximum occurs. For large systems, we expect and garrahanjacklecomtepitardvanduijvendijkvanwijland (); bodineautoninelli (). Fig. 2 shows the dependence of on the system size : the results are consistent with an exponential divergence of as , in contrast to traditional finite-size scaling at thermodynamic transitions, where scales as a power of borgskotecky (); elmatad2010 (). To gain insight into these phase transitions and explain the numerical results in Figs. 1 and 2, we now present some theoretical arguments.

Analogy with a two-dimensional thermodynamic system on a cylinder – These dynamical phase transitions in one dimension can be mapped to thermodynamic transitions in two dimensions (2d) jack_large_2010 (); elmatad2010 (). Recalling that we have taken the limit before taking , the relevant geometry for the 2d system is a long cylinder, with the system size in the dynamical system corresponding to the perimeter of the cylinder. For equilibrium systems in such geometries, the behavior near phase coexistence is sketched in Fig. 3 privman_finite-size_1983 (); borgs_crossover_1992 (): the two phases form domains arranged along the cylinder. The typical domain length scales exponentially in : the reason is that these 2d domains are separated by domain walls of length which run around the cylinder, and the associated interfacial free-energy cost scales as , so the density of domain walls is of order . From (8), is analogous to a susceptibility in the thermodynamic transition; it is also equal to a time-integral of the autocorrelation function of garrahan_first-order_2009 (). Hence scales with the relaxation time of the system. Identifying this relaxation time with the domain size in Fig. 3, the susceptibility therefore diverges as for large . The numerical data in Fig. 2 are consistent with such a divergence, indicating that the thermodynamic analogy can predict properties of the dynamical transition.

Effective interfacial model of the dynamical phase transition – We now introduce a simplified effective model for the domains in Fig. 2. Following Section 4.1 of bodineaulecomtetoninelli (), we assume that typical configurations in the model include a single active domain of size (Fig. 3). Within this domain, the system is close to its active (equilibrium) state; in the remainder of the system, the system is inactive and there are no up spins. The system contains at least one up spin so . We will show that this simplified model makes quantitatively accurate predictions for the dynamical phase transition. (For thermodynamic transitions, similar results may be available via spectral properties of the transfer matrix privman_finite-size_1983 (); borgs_crossover_1992 (): our analysis here is different, and is based on the dynamical nature of the phase transition.)

The dynamical rules of the FA model mean that the value of increases with rate and decreases with rate bodineaulecomtetoninelli (). We take reflecting boundary conditions at . (For systems that are predominately active, we interpret as the size of the largest inactive domain in the system, which has a typical value of order .) The activity in this effective model is obtained by assuming that the spins in the active domain flip with typical rate , but there are no spin flips outside this region, due to the kinetic constraint. Hence the analog of is .

The result is an effective model where undergoes a random walk whose hop rates are biased to the right (positive ), but conditioned on a relatively small time-averaged position. In the limit of large , the model can be solved exactly, as shown in supplemental material SM3 (). We summarize the main results: There is a dynamical phase transition at a field , where solves

(9) |

Here, and . Also, the susceptibility diverges as

(10) |

Finally, the scaling form of the activity near the phase transition is

(11) |

where and are specified in SM3 (). This last result is similar to that obtained in a mean-field FA model 1742-5468-2014-10-P10001 () and that for 2d equilibrium phase coexistence of ferromagnets on a cylinder geometry borgs_crossover_1992 (). These similarities indicate that the scaling function (11) might be a general property of first-order phase transitions with exponentially diverging susceptibilities. We also remark that our effective model yields the interfacial free-energy cost (9,10), which is not available from the 2d equilibrium approach of Refs. privman_finite-size_1983 (); borgs_crossover_1992 ().

The theoretical predictions (9-11) are shown in Figs. 1 and 2, together with the numerical results. The agreement is excellent, despite the simplicity of the model. The conclusion is that the finite-size scaling of the phase transition is dominated by the dynamical properties of the interface between the active and inactive regions, and this interface is accurately described by the effective model. Moreover, in the analogy with the classical phase transition on a cylinder, we can interpret the parameter in terms of an interfacial tension between the active and inactive domains shown in Fig. 2.

Discussion – There are two key outcomes of this work. First, we have shown that the cloning-with-feedback algorithm used here allows accurate characterization of dynamical phase transitions for a range of system sizes, with much greater computational efficiency than the original cloning scheme. Second, we have shown how the finite-size scaling of first-order dynamical phase transitions can be understood qualitatively by mapping them to classical phase transitions in cylindrical geometries; it can also be analyzed quantitatively by mapping to the effective interfacial model.

We expect both the numerical and theoretical methods to apply generally for dynamical phase transitions of this type: for example, application to other KCMs garrahanjacklecomtepitardvanduijvendijkvanwijland (); elmatad2010 () should be straightforward. We also anticipate application to atomistic systems that support similar phase transitions hedges_dynamic_2009 (); pitard_dynamic_2011 (); speck2012 (); malins2012 (); limmer2014 (). Moreover, the transitions considered here are directly related to quantum phase transitions in spin chains, for which results similar to (9) have been derived heyl_dynamical_2013 (). The effective interfacial model presented here provides a clear physical interpretation of such results, whose implications for quantum systems remain to be explored.

We also highlight several useful features of the numerical algorithm used here. The computational cost of the cloning algorithm scales linearly in the time . This allows the large- limit to be converged numerically. Hence, the only parameter in the finite size scaling is , which allows direct comparison with the theory presented here. This analysis is significantly simpler than finite-size scaling via path sampling, where both and must be varied together hedges_dynamic_2009 (); speck2012 (). The cloning algorithm can also be applied in systems where detailed balance is broken, where path-sampling methods are not directly applicable. Cloning methods are also related to Diffusion Quantum Monte Carlo Anderson1975_63_4 (): it would be interesting to investigate how the cloning-with-feedback method might be applied in that context casula_beyond_2006 ().

Another advantage of this method is that the control potential determined numerically provides physical insight into these dynamical phase transitions. In the active phase close to the transition, the control potential acts to suppress the number of up spins (reducing the activity), but one also finds an effective attraction between up spins SM4 (). This attraction is weak but decays slowly in space, which acts to stabilize the large spatial domains shown in Fig. 3 jack2014_east (); jack2015_hyper (). Based on the effective model, we expect that the optimal control potential should depend primarily on these domain sizes, so it naturally includes long-range interactions. It would be interesting to obtain a better understanding of optimal control potentials close to dynamical phase transitions, especially since incorporating such information into numerical methods now has the potential to significantly improve their performance. For example, one might consider transitions in other spin models garrahan_first-order_2009 () as well as exclusion processes bodineau_distribution_2005 (); jack2015_hyper (); 1751-8121-45-17-175001 (). In any case, we stress that while the ansatz (7) is much simpler than the optimal control, it still results in a significant improvement of the computational efficiency.

In conclusion, we have shown how a combination of numerical and theoretical methods provide a detailed insight into the dynamical phase transition in the FA model. The basic ideas of the method are quite general, such as the modification of the cloning algorithm with a feedback procedure to determine the optimal force, or the interfacial model as a coarse-grained description of systems near coexistence. For first-order dynamical transitions, we believe that effective interfacial models should apply rather generally. The numerical method has even broader potential application, although the choice of a suitable control potential will depend on the problem of interest – this remains to be explored.

###### Acknowledgements.

T. N. gratefully acknowledges the support of Fondation Sciences Mathématiques de Paris – EOTP NEMOT15RPO, PEPS LABS and LAABS Inphyniti CNRS project. V. L. acknowledges support by the ANR-15-CE40-0020-03 Grant LSD and by the ERC Starting Grant 680275 MALIG. T. N. and V. L. are grateful to G. Semerjian for discussions.## References

- (1) J. Robert Dorfman. An introduction to chaos in nonequilibrium statistical mechanics. Number 14 in Cambridge lecture notes in physics. Cambridge University Press Cambridge ; New York 1999.
- (2) N. G. van Kampen. Stochastic processes in physics and chemistry. North-Holland personal library. Elsevier Amsterdam ; Boston 3rd ed edition 2007.
- (3) H Hinrichsen. Non-equilibrium critical phenomena and phase transitions into absorbing states. Adv. Phys. 49, 815 (2000). Advances In Physics.
- (4) M R Evans and T Hanney. Nonequilibrium statistical mechanics of the zero-range process and related models. Journal of Physics A: Mathematical and General 38, R195 (2005).
- (5) Lorenzo Bertini, Andrea De Sole, Davide Gabrielli, Giovanni Jona-Lasinio, and Claudio Landim. Current Fluctuations in Stochastic Lattice Gases. Phys. Rev. Lett. 94, 030601 (2005).
- (6) Thierry Bodineau and Bodineau Derrida. Distribution of current in nonequilibrium diffusive systems and phase transitions. Phys. Rev. E 72, 066110 (2005).
- (7) Juan P. Garrahan, Robert L. Jack, Vivien Lecomte, Estelle Pitard, Kristina van Duijvendijk, and Frédéric van Wijland. Dynamical First-Order Phase Transition in Kinetically Constrained Models of Glasses. Phys. Rev. Lett. 98, 195702 (2007).
- (8) Hugo Touchette. The large deviation approach to statistical mechanics. Physics Reports 478, 1 (2009).
- (9) Mauro Merolle, Juan P. Garrahan, and David Chandler. Space-time thermodynamics of the glass transition. PNAS 102, 10837 (2005).
- (10) Robert L. Jack, Juan P. Garrahan, and David Chandler. Space-time thermodynamics and subsystem observables in a kinetically constrained model of glassy materials. J. Chem. Phys. 125, 184509 (2006).
- (11) Juan P. Garrahan, Robert L. Jack, Vivien Lecomte, Estelle Pitard, Kristina van Duijvendijk, and Frédéric van Wijland. First-order dynamical phase transition in models of glasses: an approach based on ensembles of histories. J. Phys. A 42, 075007 (2009).
- (12) Jorge Kurchan. Six out of equilibrium lectures. arXiv:0901.1271 [cond-mat] (2009). arXiv: 0901.1271.
- (13) Lester O. Hedges, Robert L. Jack, Juan P. Garrahan, and David Chandler. Dynamic Order-Disorder in Atomistic Models of Structural Glass Formers. Science 323, 1309 (2009).
- (14) Juan P. Garrahan, Peter Sollich, and Cristina Toninelli. Kinetically Constrained Models. arXiv:1009.6113 and Chapter of “Dynamical heterogeneities in glasses, colloids, and granular media”, Eds.: L. Berthier, G. Biroli, J-P Bouchaud, L. Cipelletti and W. van Saarloos (Oxford University Press, 2011) (2010).
- (15) David Chandler and Juan P. Garrahan. Dynamics on the Way to Forming Glass: Bubbles in Space-Time. Ann. Rev. Phys. Chem. 61, 191 (2010).
- (16) David T. Limmer and David Chandler. Theory of amorphous ices. Proceedings of the National Academy of Sciences 111, 9413 (2014).
- (17) Thierry Bodineau and Bernard Derrida. Cumulants and large deviations of the current through non-equilibrium steady states. Comptes Rendus Physique 8, 540 (2007).
- (18) L. Bertini, A. De Sole, D. Gabrielli, G. Jona-Lasinio, and C. Landim. Non Equilibrium Current Fluctuations in Stochastic Lattice Gases. Journal of Statistical Physics 123, 237 (2006).
- (19) Cristian Giardinà, Jorge Kurchan, and Luca Peliti. Direct Evaluation of Large-Deviation Functions. Phys. Rev. Lett. 96, 120603 (2006).
- (20) Estelle Pitard, Vivien Lecomte, and Frédéric Van Wijland. Dynamic transition in an atomic glass former: a molecular dynamics evidence. Europhys. Lett. 96, 56002 (2011).
- (21) Thomas Speck and David Chandler. Constrained dynamics of localized excitations causes a non-equilibrium phase transition in an atomistic model of glass formers. The Journal of Chemical Physics 136, 184509 (2012).
- (22) Thomas Speck, Alex Malins, and C. Patrick Royall. First-Order Phase Transition in a Model Glass Former: Coupling of Local Structure and Dynamics. Phys. Rev. Lett. 109, 195703 (2012).
- (23) Vladimir Privman and Michael E. Fisher. Finite-size effects at first-order transitions. Journal of Statistical Physics 33, 385 (1983).
- (24) Christian Borgs and Roman Kotecký. A rigorous theory of finite-size scaling at first-order phase transitions. Journal of Statistical Physics 61, 79 (1990).
- (25) Bernd A. Berg and Thomas Neuhaus. Multicanonical ensemble: A new approach to simulate first-order phase transitions. Phys. Rev. Lett. 68, 9 (1992).
- (26) Fugao Wang and D. P. Landau. Efficient, Multiple-Range Random Walk Algorithm to Calculate the Density of States. Phys. Rev. Lett. 86, 2050 (2001).
- (27) Alastair D. Bruce and Nigel B. Wilding. Computational strategies for mapping equilibrium phase diagrams. Adv. Chem. Phys. 127, 1 (2003).
- (28) Todd R. Gingrich and Phillip L. Geissler. Preserving correlations between trajectories for efficient path sampling. The Journal of Chemical Physics 142 (2015).
- (29) Glenn H. Fredrickson and Hans C. Andersen. Kinetic Ising Model of the Glass Transition. Phys. Rev. Lett. 53, 1244 (1984).
- (30) Takahiro Nemoto, Freddy Bouchet, Robert L. Jack, and Vivien Lecomte. Population-dynamics method with a multicanonical feedback control. Phys. Rev. E 93, 062123 (2016).
- (31) Yael S. Elmatad, Robert L. Jack, David Chandler, and Juan P. Garrahan. Finite-temperature critical point of a glass transition. Proceedings of the National Academy of Sciences 107, 12793 (2010).
- (32) Felix Ritort and Peter Sollich. Glassy dynamics of kinetically constrained models. Advances in Physics 52, 219 (2003).
- (33) Cristian Giardinà, Jorge Kurchan, Vivien Lecomte, and Julien Tailleur. Simulating rare events in dynamical processes. J. Stat. Phys. 145, 787 (2011).
- (34) Takahiro Nemoto, Esteban Guevara Hidalgo, and Vivien Lecomte. Finite-time and finite-size scalings in the evaluation of large-deviation functions: Analytical study using a birth-death process. Phys. Rev. E 95, 012102 (2017).
- (35) Robert L. Jack and Peter Sollich. Large Deviations and Ensembles of Trajectories in Stochastic Models. Prog. Theor. Phys. Supplement 184, 304 (2010).
- (36) Vladislav Popkov, Gunter M Schütz, and Damien Simon. ASEP on a ring conditioned on enhanced flux. Journal of Statistical Mechanics: Theory and Experiment 2010, P10007 (2010).
- (37) R Chetrite and Hugo Touchette. Nonequilibrium Microcanonical and Canonical Ensembles and Their Equivalence. Phys. Rev. Lett. 111, 120601 (2013).
- (38) Takahiro Nemoto and Shin-ichi Sasa. Computation of Large Deviation Statistics via Iterative Measurement-and-Feedback Procedure. Phys. Rev. Lett. 112, 090602 (2014).
- (39) R. L. Jack and P. Sollich. Effective interactions and large deviations in stochastic processes. Eur. Phys. J. Spec. Top. 224, 2351 (2015).
- (40) See Supplemental Material for the details of the feedback algorithm, first section.
- (41) See Supplemental Material for numerical evidence supporting the feedback algorithm compared to the standard one, second section, which includes Ref. SecondPart ().
- (42) Thierry Bodineau and Cristina Toninelli. Activity Phase Transition for Constrained Dynamics. Communications in Mathematical Physics 311, 357 (2012).
- (43) Vladimir Privman and Michael E. Fisher. Finite-size effects at first-order transitions. J Stat Phys 33, 385 (1983).
- (44) Christian Borgs and John Z. Imbrie. Crossover finite-size scaling at first-order transitions. J Stat Phys 69, 487 (1992).
- (45) Thierry Bodineau, Vivien Lecomte, and Cristina Toninelli. Finite Size Scaling of the Dynamical Free-Energy in a Kinetically Constrained Model. Journal of Statistical Physics 147, 1 (2012).
- (46) See Supplemental Material for the derivation of Eqs. (9), (10) and (11), third section.
- (47) T. Nemoto, V. Lecomte, S. Sasa, and F. van Wijland. Finite-size effects in a mean-field kinetically constrained model: dynamical glassiness and quantum criticality. Journal of Statistical Mechanics: Theory and Experiment 2014, P10001 (2014).
- (48) M. Heyl, A. Polkovnikov, and S. Kehrein. Dynamical Quantum Phase Transitions in the Transverse-Field Ising Model. Phys. Rev. Lett. 110, 135704 (2013).
- (49) James B. Anderson. A randomâwalk simulation of the Schrödinger equation: . The Journal of Chemical Physics 63, 1499 (1975).
- (50) Michele Casula. Beyond the locality approximation in the standard diffusion Monte Carlo method. Phys. Rev. B 74, 161102 (2006).
- (51) See Supplemental Material for a study of the effective interactions around the dynamical phase transition, fourth section.
- (52) Robert L Jack and Peter Sollich. Large deviations of the dynamical activity in the East model: analysing structure in biased trajectories. Journal of Physics A: Mathematical and Theoretical 47, 015003 (2014).
- (53) Robert L. Jack, Ian R. Thompson, and Peter Sollich. Hyperuniformity and Phase Separation in Biased Ensembles of Trajectories for Diffusive Systems. Phys. Rev. Lett. 114, 060601 (2015).
- (54) Vivien Lecomte, Juan P Garrahan, and Frédéric van Wijland. Inactive dynamical phase of a symmetric exclusion process on a ring. Journal of Physics A: Mathematical and Theoretical 45, 175001 (2012).
- (55) Esteban Guevara Hidalgo, Takahiro Nemoto, and Vivien Lecomte. Finite-Time and -Size Scalings in the Evaluation of Large Deviation Functions - Part II: Numerical Approach in Continuous Time. arXiv:1607.08804 [cond-mat] (2016).

See pages ,1,,2,,3,,4,,5,,6,,7,,8,,9,,10,,11,,12,,13,,14 of Supplemental_Material.pdf