© 2026 Aqora Quantum S.A.S.

Terms Privacy Legal Notice

The Journal of Physical Chemistry Letters

Elucidating Dynamical Behaviors in Kinetically Constrained Models via Energy-Activity Double-Biased Matrix Product State Analysis

Yoong Hee Lee, Jay-Hak Lee, YounJoon Jung

Published Oct 18, 2024• 1 citation

Abstract

We investigate dynamical phase transitions in two representative kinetically constrained models: the 1D Fredrickson–Andersen and East Models. A recently developed energy-activity double-bias approach utilizing both s and g fields, conjugated with dynamical activity and trajectory energy, is combined with matrix product state (MPS) methods. It is demonstrated that MPS methods facilitate the numerical approximation of large-deviation statistics of dynamics by determining the eigenvalues of tilted dynamical generators under the influence of double-biasing fields. Specifically, by focusing on the g -field, a nearly “half-filled” state at moderate negative g values is identified, indicating the potential existence of an anomalous phase. Additionally, dynamical quantities under various s, g, and T conditions are obtained via tensor networks, showing good qualitative consistency with mean-field results and our previous extensive numerical simulation results obtained by the path sampling method. Our study introduces novel methodologies for examining decoupled dynamical behaviors that, although energetically active, remain dynamically inactive within the system. This approach offers a fresh perspective on the theoretical framework and computational strategies for studying dynamical phase transitions in kinetically constrained systems.

DOI: 10.1021/acs.jpclett.4c02197

Related papers

Using Matrix Product States to Study the Dynamical Large Deviations of Kinetically Constrained Models

Mari Carmen Bañuls, Juan P. Garrahan

Here we demonstrate that tensor network techniques-originally devised for the analysis of quantum many-body problems-are well suited for the detailed study of rare event statistics in kinetically constrained models (KCMs). As concrete examples, we consider the Fredrickson-Andersen and East models, two paradigmatic KCMs relevant to the modeling of glasses. We show how variational matrix product states allow us to numerically approximate-systematically and with high accuracy-the leading eigenstates of the tilted dynamical generators, which encode the large deviation statistics of the dynamics. Via this approach, we can study system sizes beyond what is possible with other methods, allowing us to characterize in detail the finite size scaling of the trajectory-space phase transition of these models, the behavior of spectral gaps, and the spatial structure and "entanglement" properties of dynamical phases. We discuss the broader implications of our results.

arXiv:1903.01570•Physical Review Letters•Nov 12, 2019•70 citations

Finite Time Large Deviations via Matrix Product States

Luke Causer, Mari Carmen Bañuls, Juan P. Garrahan

Recent work has shown the effectiveness of tensor network methods for computing large deviation functions in constrained stochastic models in the infinite time limit. Here we show that these methods can also be used to study the statistics of dynamical observables at arbitrary finite time. This is a harder problem because, in contrast to the infinite time case, where only the extremal eigenstate of a tilted Markov generator is relevant, for finite time the whole spectrum plays a role. We show that finite time dynamical partition sums can be computed efficiently and accurately in one dimension using matrix product states and describe how to use such results to generate rare event trajectories on demand. We apply our methods to the Fredrickson-Andersen and East kinetically constrained models and to the symmetric simple exclusion process, unveiling dynamical phase diagrams in terms of counting field and trajectory time. We also discuss extensions of this method to higher dimensions.

arXiv:2108.11418•Physical Review Letters•Mar 4, 2022•29 citations

Dynamics and large deviation transitions of the XOR-Fredrickson-Andersen kinetically constrained model

Luke Causer, Igor Lesanovsky, Mari Carmen Bañuls +1 more

We study a one-dimensional classical stochastic kinetically constrained model (KCM) inspired by Rydberg atoms in their "facilitated" regime, where sites can flip only if a single of their nearest neighbors is excited. We call this model "XOR-FA" to distinguish it from the standard Fredrickson-Andersen (FA) model. We describe the dynamics of the XOR-FA model, including its relation to simple exclusion processes in its domain wall representation. The interesting relaxation dynamics of the XOR-FA is related to the prominence of large dynamical fluctuations that lead to phase transitions between active and inactive dynamical phases as in other KCMs. By means of numerical tensor network methods we study in detail such transitions in the dynamical large deviation regime.

arXiv:2006.08673•Physical review. E•Nov 25, 2020•32 citations

Dynamical Large Deviations of Two-Dimensional Kinetically Constrained Models Using a Neural-Network State Ansatz

Corneel Casert, Tom Vieijra, Stephen Whitelam +1 more

We use a neural-network ansatz originally designed for the variational optimization of quantum systems to study dynamical large deviations in classical ones. We use recurrent neural networks to describe the large deviations of the dynamical activity of model glasses, kinetically constrained models in two dimensions. We present the first finite size-scaling analysis of the large-deviation functions of the two-dimensional Fredrickson-Andersen model, and explore the spatial structure of the high-activity sector of the South-or-East model. These results provide a new route to the study of dynamical large-deviation functions, and highlight the broad applicability of the neural-network state ansatz across domains in physics.

arXiv:2011.08657•Physical Review Letters•Sep 16, 2021•20 citations

Dynamical heterogeneity and large deviations in the open quantum East glass model from tensor networks

Luke Causer, Mari Carmen Bañuls, Juan P. Garrahan

We study the nonequilibrium dynamics of the dissipative quantum East model via numerical tensor networks. We use matrix product states to represent evolution under quantum-jump unravellings for sizes beyond those accessible to exact diagonalization. This allows us to demonstrate that dynamical heterogeneity accompanies slow relaxation, in analogy with what is seen in classical glassy systems. Furthermore, using variational matrix product operators we (i) compute the spectral gap of the Lindbladian and show that glassiness is enhanced in the presence of weak quantum fluctuations compared to the pure classical case, and (ii) obtain the dynamical large deviations by calculating the leading eigenvector of the tilted Lindbladian and find clear evidence for a first-order active-inactive dynamical phase transition. We also show how to directly sample the rare quantum trajectories associated to the large deviations.

arXiv:2404.03750•Physical review. B./Physical review. B•Feb 13, 2025•8 citations

Optimal sampling of dynamical large deviations via matrix product states

Luke Causer, Mari Carmen Bañuls, Juan P. Garrahan

The large deviation statistics of dynamical observables is encoded in the spectral properties of deformed Markov generators. Recent works have shown that tensor network methods are well suited to compute accurately the relevant leading eigenvalues and eigenvectors. However, the efficient generation of the corresponding rare trajectories is a harder task. Here, we show how to exploit the matrix product state approximation of the dominant eigenvector to implement an efficient sampling scheme which closely resembles the optimal (so-called "Doob") dynamics that realizes the rare events. We demonstrate our approach on three well-studied lattice models, the Fredrickson-Andersen and East kinetically constrained models, and the symmetric simple exclusion process. We discuss how to generalize our approach to higher dimensions.

arXiv:2103.01265•Physical review. E•Jun 28, 2021•27 citations

We investigate dynamical phase transitions in two representative kinetically constrained models: the 1D Fredrickson–Andersen and East Models. A recently developed energy-activity double-bias approach utilizing both s and g fields, conjugated with dynamical activity and trajectory energy, is combined with matrix product state (MPS) methods. It is demonstrated that MPS methods facilitate the numerical approximation of large-deviation statistics of dynamics by determining the eigenvalues of tilted dynamical generators under the influence of double-biasing fields. Specifically, by focusing on the g -field, a nearly “half-filled” state at moderate negative g values is identified, indicating the potential existence of an anomalous phase. Additionally, dynamical quantities under various s, g, and T conditions are obtained via tensor networks, showing good qualitative consistency with mean-field results and our previous extensive numerical simulation results obtained by the path sampling method. Our study introduces novel methodologies for examining decoupled dynamical behaviors that, although energetically active, remain dynamically inactive within the system. This approach offers a fresh perspective on the theoretical framework and computational strategies for studying dynamical phase transitions in kinetically constrained systems.

Here we demonstrate that tensor network techniques-originally devised for the analysis of quantum many-body problems-are well suited for the detailed study of rare event statistics in kinetically constrained models (KCMs). As concrete examples, we consider the Fredrickson-Andersen and East models, two paradigmatic KCMs relevant to the modeling of glasses. We show how variational matrix product states allow us to numerically approximate-systematically and with high accuracy-the leading eigenstates of the tilted dynamical generators, which encode the large deviation statistics of the dynamics. Via this approach, we can study system sizes beyond what is possible with other methods, allowing us to characterize in detail the finite size scaling of the trajectory-space phase transition of these models, the behavior of spectral gaps, and the spatial structure and "entanglement" properties of dynamical phases. We discuss the broader implications of our results.

Recent work has shown the effectiveness of tensor network methods for computing large deviation functions in constrained stochastic models in the infinite time limit. Here we show that these methods can also be used to study the statistics of dynamical observables at arbitrary finite time. This is a harder problem because, in contrast to the infinite time case, where only the extremal eigenstate of a tilted Markov generator is relevant, for finite time the whole spectrum plays a role. We show that finite time dynamical partition sums can be computed efficiently and accurately in one dimension using matrix product states and describe how to use such results to generate rare event trajectories on demand. We apply our methods to the Fredrickson-Andersen and East kinetically constrained models and to the symmetric simple exclusion process, unveiling dynamical phase diagrams in terms of counting field and trajectory time. We also discuss extensions of this method to higher dimensions.

We study a one-dimensional classical stochastic kinetically constrained model (KCM) inspired by Rydberg atoms in their "facilitated" regime, where sites can flip only if a single of their nearest neighbors is excited. We call this model "XOR-FA" to distinguish it from the standard Fredrickson-Andersen (FA) model. We describe the dynamics of the XOR-FA model, including its relation to simple exclusion processes in its domain wall representation. The interesting relaxation dynamics of the XOR-FA is related to the prominence of large dynamical fluctuations that lead to phase transitions between active and inactive dynamical phases as in other KCMs. By means of numerical tensor network methods we study in detail such transitions in the dynamical large deviation regime.

We use a neural-network ansatz originally designed for the variational optimization of quantum systems to study dynamical large deviations in classical ones. We use recurrent neural networks to describe the large deviations of the dynamical activity of model glasses, kinetically constrained models in two dimensions. We present the first finite size-scaling analysis of the large-deviation functions of the two-dimensional Fredrickson-Andersen model, and explore the spatial structure of the high-activity sector of the South-or-East model. These results provide a new route to the study of dynamical large-deviation functions, and highlight the broad applicability of the neural-network state ansatz across domains in physics.

We study the nonequilibrium dynamics of the dissipative quantum East model via numerical tensor networks. We use matrix product states to represent evolution under quantum-jump unravellings for sizes beyond those accessible to exact diagonalization. This allows us to demonstrate that dynamical heterogeneity accompanies slow relaxation, in analogy with what is seen in classical glassy systems. Furthermore, using variational matrix product operators we (i) compute the spectral gap of the Lindbladian and show that glassiness is enhanced in the presence of weak quantum fluctuations compared to the pure classical case, and (ii) obtain the dynamical large deviations by calculating the leading eigenvector of the tilted Lindbladian and find clear evidence for a first-order active-inactive dynamical phase transition. We also show how to directly sample the rare quantum trajectories associated to the large deviations.

The large deviation statistics of dynamical observables is encoded in the spectral properties of deformed Markov generators. Recent works have shown that tensor network methods are well suited to compute accurately the relevant leading eigenvalues and eigenvectors. However, the efficient generation of the corresponding rare trajectories is a harder task. Here, we show how to exploit the matrix product state approximation of the dominant eigenvector to implement an efficient sampling scheme which closely resembles the optimal (so-called "Doob") dynamics that realizes the rare events. We demonstrate our approach on three well-studied lattice models, the Fredrickson-Andersen and East kinetically constrained models, and the symmetric simple exclusion process. We discuss how to generalize our approach to higher dimensions.