Pith. sign in

REVIEW 4 major objections 5 minor 56 references

A kinetically constrained model exhibiting non-linear diffusion and jamming

T0 review · 4 major / 5 minor · reviewed 2026-08-11 · deepseek-v4-flash

Pith's one-line read The diffusion coefficient of the constrained lattice gas is $D=3(1-\rho)$ up to a jamming transition at density $2/3$, and the two-hole sector is solved exactly with $D=3/8$.

desk verdict A new KCM with exact jammed-entropy and doublon results; the mean-field D=3(1−ρ) formula is elegant but not proven, and the paper should be sent to referees with that distinction made explicit. read the letter →

arxiv 2412.05231 v4 pith:ZKJ62OHE submitted 2024-12-06 cond-mat.stat-mech

classification cond-mat.stat-mech MSC 82C2282C2682C3160K35 PACS 05.40.-a05.60.-k
keywords kineticallyconstrainedmodeltriangularladderdiffusioncoefficientjammingtransitionclassical-quantummappingmean-fieldapproximationdoublonfermionization
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

The paper studies particles on a triangular ladder that can hop only when they are the only particle in their triangle. By recasting the random motion as a quantum Hamiltonian, the authors derive a closed-form diffusion constant $D=3(1-\rho)$ for how density disturbances spread, where $\rho$ is the fraction of occupied sites. They further show that at density $\rho=2/3$ the system jams, with exponentially many configurations becoming frozen, and they solve exactly the two-hole sector with diffusion constant $D=3/8$. A sympathetic reader would care because the model gives an analytically tractable example of a strongly interacting diffusive process with a sharp dynamical transition.

What carries the argument

The load-bearing device is the classical-to-quantum mapping: the rate matrix of the master equation is reinterpreted as a stoquastic Hamiltonian whose ground state is the uniform distribution, so the lowest excitation gap at momentum $Q=2\pi/L$ is $DQ^2$ and yields the diffusion constant. Two calculations carry the result. First, the mean-field replacement $a_i^\dagger a_i\to\rho$ turns the constrained hopping into free fermions with dispersion $\epsilon_k=(3-\cos 2k-2\cos k)(1-\rho)$, whose curvature fixes $D=3(1-\rho)$. Second, the two-hole sector reduces to a $2L$-site decorated chain whose reduced momentum-space matrix gives the exact two-hole diffusion constant $D=3/8$.

What would settle it

Extrapolate the exact lowest gap in the $Q=2\pi/L$ sector at density $\rho=0.6$ for chain lengths up to about $L=100$; if the resulting diffusion constant differs from $D=3(1-\rho)=1.2$ beyond the extrapolation uncertainty, the mean-field formula is only approximate.

Watch

Extended reading notes

Core claim

The paper claims that the diffusion coefficient of this kinetically constrained lattice gas is $D=3(1-\rho)$ throughout the diffusive regime up to the jamming point, with the single-particle value $D_0=3$ recovered at zero density. The derivation maps the classical master equation to a fermionic Hamiltonian via a standard fermionization, then applies a mean-field decoupling that replaces local occupancies by the average density $\rho$. The resulting free-fermion dispersion $\epsilon_k=(3-\cos 2k-2\cos k)(1-\rho)$ has curvature at $k=0$ that gives an effective mass $m=1/[6(1-\rho)]$ and hence $D=3(1-\rho)$. Above $\rho=2/3$, configurations such as $110110110\ldots$ become frozen, and the number of jammed configurations grows exponentially with an entropy the paper computes by saddle-point counting, $s_j(\rho)=(\rho-1)\ln((1-\rho)/(2\rho-1))+(2-3\rho)\ln((3\rho-2)/(2\rho-1))$; numerical spectra follow the mean-field curve closely until this transition.

Load-bearing premise

The load-bearing premise is that every particle moves under the influence of the average occupancy $\rho$ rather than local fluctuations, and the paper provides no small parameter that makes this replacement controlled.

Editorial extensions

If this is right

  • Below density $2/3$, density disturbances spread with diffusion constant $D=3(1-\rho)$, so the transport coefficient is known in closed form for the whole diffusive regime.
  • At $\rho=2/3$ the system acquires a finite density of exactly frozen configurations; states initialized with overlap on them never fully relax.
  • With exactly two holes, the only mobile object is a bound pair of holes, called a doublon, with diffusion constant $D=3/8$ independent of system size.
  • With three or more holes, the diffusion constant vanishes as $c_n/L$, with the first coefficients $c_3=3/8$ and $c_4=3/4$, matching the conjecture $c_n=(3/8)\lfloor n/2 \rfloor$.
  • The quasiparticle mass $m=1/[6(1-\rho)]$ diverges as $\rho\to 1$, making the full-occupancy state dynamically isolated.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • Beyond the paper, if the mean-field curve is exact, the fermionic Hamiltonian likely hides an integrable or frustration-free structure, and a direct construction of excited eigenstates would reveal it.
  • The same mapping on a full two-dimensional triangular lattice should produce a similar linear-in-$(1-\rho)$ diffusion law with a lattice-dependent prefactor; Monte Carlo measurement at $\rho=1/2$ would test whether the mean-field form survives in higher dimensions.
  • The jamming-entropy calculation suggests a clear finite-size signature: above $\rho=2/3$, the fraction of frozen configurations should grow as $\exp[L(s_j(\rho)-s(\rho))]$, which is observable by brute-force enumeration on chains of length about $L=50$.
  • The paper leaves open whether the finite-hole diffusion coefficients follow the simple doublon-counting rule $c_n=(3/8)\lfloor n/2\rfloor$ at larger $n$; a direct spectral calculation for $n=6$ and $n=7$ would be a sharp test.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

4 major / 5 minor

Summary. The paper introduces a classical kinetically constrained lattice gas on a triangular ladder in which a particle may hop only when both nearest-neighbor sites on the ladder are empty. The master operator is mapped to a fermionic Hamiltonian, and a Hartree-Fock treatment in Sec. 3 replaces local density operators by their expectation value ρ, yielding the quasiparticle dispersion ε_k = (3 − cos 2k − 2 cos k)(1 − ρ) and hence the mean-field diffusion constant D_MF = 3(1 − ρ). Section 4 computes the entropy of jammed configurations for ρ ≥ 2/3 in closed form. Section 5 exactly solves the two-hole sector, obtaining the doublon dispersion and D = 3/8 for N = L − 2 in the thermodynamic limit. Section 6 presents exact diagonalization results for finite chains and extrapolations of D(ρ) that are compared with Eq. (28). The central claim is that D(ρ) = 3(1 − ρ) describes the model up to ρ ≈ 2/3 and that the model displays a jamming transition at ρ = 2/3.

Significance. If the mean-field formula were exact, the paper would provide a closed-form, density-dependent diffusion coefficient for a kinetically constrained model together with an exactly solvable jamming entropy and an exact two-hole spectrum. The strongest parts of the paper are rigorous and well supported: the jammed-configuration entropy in Eq. (35), the exact doublon dispersion in Eq. (39), and the numerical verification of D = 3/8 in Fig. 8 are concrete, checkable results. The paper is also honest in labeling the mean-field step as an approximation, and the absence of fitted parameters in Eq. (28) is a virtue. However, the derivation of Eq. (28) is not controlled, and the paper's numerical support weakens precisely in the region where the approximation is most needed; the significance of the central claim therefore depends on whether the mean-field result can be justified or reframed as a conjecture.

major comments (4)
  1. [Sec. 3, Eq. (25)] The Hartree-Fock replacement of the density operators a†_i a_i by the c-number ρ is uncontrolled. The manuscript provides no small parameter, large-dimension limit, or systematic expansion justifying this decoupling, and the resulting quadratic Hamiltonian is the entire origin of Eq. (28). Since this equation is the main quantitative claim, the paper should either supply a controlled approximation scheme with estimates of the neglected terms or explicitly present D_MF = 3(1 − ρ) as a conjecture supported by numerics.
  2. [Sec. 3, Eqs. (26)–(28)] There is a mismatch between the object whose energy is computed and the object that defines diffusion in Sec. 2.3. The quasiparticle state a†_k|Ψ(z)⟩ changes the particle number, whereas H conserves N; the relation [a†_k, H] = ε_k a†_k holds only for the approximate H_HF, not for the exact Hamiltonian. The diffusion constant in Eq. (28) is identified with the fixed-N, Q = 2π/L density-mode gap, but the calculation actually gives the energy of adding a particle to a structureless grand-canonical background. The exact two-hole result D = 3/8 in Eq. (42) versus D_MF → 0 as ρ → 1 makes the nonuniformity concrete, so the agreement at low density is empirical rather than explained.
  3. [Sec. 6, Figs. 1 and 7] The comparison with Eq. (28) is not quantitatively controlled because the extrapolation procedure changes between density regimes: Fig. 7 uses quadratic fits in Q for ρ ≤ 1/2 and linear fits for ρ ≥ 2/3. The caption of Fig. 1 attributes the high-density disagreement to 'progressively less accurate extrapolation,' but no evidence or model of the finite-size corrections is given. At ρ = 3/4 the extrapolated value 0.572 ± 0.020 already differs from D_MF = 0.75 by about 24%, so the claimed validity 'until the jamming point ρ = 2/3' needs a careful discussion of which fit form is justified and what systematic error it introduces.
  4. [Sec. 5, Eq. (42)] The fixed-hole-density regime is not addressed by the mean-field formula. For N = L − n with n = O(1) and L → ∞, the density ρ → 1 but the exact diffusion constant is D = 3/8 for n = 2 and D_n = c_n/L for n ≥ 3. This shows that the limit ρ → 1 is nonuniform in how the thermodynamic limit is taken. The manuscript should state this limitation explicitly when presenting Eq. (28) as a global prediction, because it prevents the formula from being extended to the high-density edge by continuity.
minor comments (5)
  1. [Abstract] The abstract says that at the critical density ρ = 2/3 'exponentially many configurations become jammed,' but the text of Sec. 4 states that at exactly ρ = 2/3 only three jammed configurations appear and the exponential growth occurs for ρ > 2/3; please rephrase to avoid the incorrect implication.
  2. [Sec. 1] There is a typo in 'kynetically constrained process'; it should be 'kinetically constrained process.'
  3. [Sec. 2.1, Eq. (1)] The transition-rate matrix is written for one orientation of the triangle, but the indexing of rows and columns 000,...,111 is not stated explicitly; adding the basis ordering would make the matrix easier to verify.
  4. [Sec. 6.1, Fig. 6] The caption of Fig. 6 lists fits D1P, D2P, etc., but the text does not define these symbols; please define the notation or remove it.
  5. [Sec. 7] The statement that 'we see little deviations from the MF result in the numerics' is vague in light of the deviations already visible at ρ = 3/4 in Fig. 7; please state the quantitative range over which the deviations are small.

Circularity Check

0 steps flagged · score 1.0 of 10

No significant circularity: D=3(1−ρ) is a parameter-free mean-field result, and the exact doublon and jammed-entropy calculations are derived independently.

full rationale

I find no circular step in the paper's derivation chain. The central claim D_MF=3(1−ρ) (Eq. 28) follows from a parameter-free Hartree-Fock treatment (Eq. 25) of the fermionized Hamiltonian: local density operators are replaced by their expectation value ρ, and the coefficient 3 is fixed by the triangular-ladder geometry rather than by any fit to numerical data. The exact doublon diffusion constant D=3/8 (Eq. 42) is obtained by diagonalizing the two-hole sector, and the jammed-configuration entropy s_j(ρ) (Eq. 35) is an independent combinatorial saddle-point calculation; both agree with exact diagonalization. The paper cites the authors' companion work [38] for the model rationale and the diffusion-cascade concept, and these citations are used when explaining why the Q=2π/L sector is relevant, but the mean-field derivation itself does not rest on [38]; no uniqueness theorem or imported ansatz is invoked as a load-bearing premise. Concerns about the uncontrolled character of the Hartree-Fock decoupling and about the differing linear/quadratic extrapolations in Fig. 7 are correctness or numerical-analysis issues, not cases where a prediction is equivalent to its input by construction.

Assumptions & free parameters 1 free parameters · 8 assumptions · 0 invented entities

The paper introduces no new particles or forces. The quasiparticles and doublons are derived excitations, not postulated entities. The main assumptions are the mean-field decoupling (ad hoc to this paper), the transport-limit identification of D from a finite-Q gap, and the standard classical-quantum mapping. The counting of jammed configurations is exact and relies on the precise definition of the kinetic rule. There are no free parameters fitted to data.

free parameters (1)
  • None
    The model parameters p, a, and Δt are set to 1 for convenience and are not fitted. The fugacity z in the coherent state is fixed by the density ρ through the constraint relation |z|² = ρ/(1 − ρ), not by fitting.
assumptions (8)
  • standard math The rate matrix is symmetric (detailed balance), so H has a complete set of eigenvectors and eigenvalues.
    Invoked in Sec. 2.1 and 2.2 to justify the spectral decomposition of the probability evolution.
  • domain assumption The classical-to-quantum mapping via the square root of the probability distribution is valid and yields a stoquastic Hamiltonian.
    Introduced in Sec. 2.2 to reinterpret the rate matrix as a quantum many-body Hamiltonian; standard in the field.
  • standard math The ground state in each particle-number sector is the uniform superposition of all configurations.
    Follows from the doubly stochastic property of the rate matrix; stated in Sec. 2.2.
  • ad hoc to paper The coherent state |Ψ(z)⟩ is a good variational state for the ground state in the large-L limit.
    Assumed in Sec. 3 to justify the mean-field treatment; the saddle-point approximation is standard but not rigorously controlled.
  • ad hoc to paper In the Hartree-Fock approximation, density operators a†_i a_i can be replaced by their expectation value ρ in the interaction terms.
    This is the central approximation in Sec. 3 leading to the free-fermion HHF and the dispersion in Eq. (26); no small parameter is provided.
  • domain assumption The diffusion constant is obtained from the lowest excitation energy at Q = 2π/L via ϵ_Q = D Q², and the diffusion-cascade processes are negligible.
    Discussed in Sec. 2.3; the paper acknowledges ignoring the subtlety of the transport limit and assumes the single-quasiparticle dispersion gives D.
  • domain assumption A configuration is jammed if and only if it contains no particle with both neighboring sites empty, i.e., zeros are isolated and separated by at least two ones.
    Used in Sec. 4 to count jammed configurations; the counting of sequences n_i ≥ 2 follows directly from this characterization.
  • standard math The saddle-point evaluation of the jammed-entropy generating function is valid for large L.
    Used in Sec. 4 to derive the exponential growth rate s_j(ρ).

how reviews work

0 comments
Cite this review

Pith. "Pith review of A kinetically constrained model exhibiting non-linear diffusion and jamming." pith.science (2026). https://pith.science/paper/ZKJ62OHE

@misc{pith2026241205231,
  author       = {Pith},
  title        = {Pith review of: A kinetically constrained model exhibiting non-linear diffusion and jamming},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/ZKJ62OHE}},
  note         = {Machine review of arXiv:2412.05231}
}
read the original abstract

We present a classical kinetically constrained model of interacting particles on a triangular ladder, which displays diffusion and jamming and can be treated by means of a classical-quantum mapping. Interpreted as a theory of interacting fermions, the diffusion coefficient is the inverse of the effective mass of the quasiparticles which can be computed using mean-field theory. At a critical density \r{ho} = 2/3, the model undergoes a dynamical phase transition in which exponentially many configurations become jammed while others remain diffusive. The model can be generalized to two dimensions.

Figures

Figures reproduced from arXiv: 2412.05231 by the authors.

Figure 1
Figure 1. Comparison of mean-field solution (dashed straight-line) against extrapolated numerical results (green diamonds) – see Figs. 7 and 8; Notice that the numerical data follow the mean-field prediction of Eq.(28) quite closely until the jamming point ρ = 2/3. The disagreement at high density may be due to progressively less accurate extrapolation of numerics. structural glasses [25], and many other examples [35]. In thi… view at source ↗
Figure 2
Figure 2. The rules for hopping. In active triangles, where only one particle out of 3 is present, it can hop with the same probability along the green arrows. Where two or more particles are present in a triangle (red cross) such triangle is inactive. On the lower axis, the numeration convention we use. is proportional to the inverse mass of the quasiparticle. By construction, the diffusion coefficient is expected to decreas… view at source ↗
Figure 3
Figure 3. Entropy of jammed configurations (blue) and total entropy of the system (yellow). If the initial probability distribution P(0) has overlap with one of these, then it will never fully relax to the equipartite equilibrium Peq = 1 Z P σ |σ⟩ = 1 Z1/2 |E0⟩, but rather leave behind a ”localized fraction”[18], a jammed fingerprint of the initial state, whose weight will be exponentially small log N /Z = 2(1 − ρ) 2 + 3(1 − … view at source ↗
Figures from the paper (5 more)
Figure 4
Figure 4. Figure 4: Reduced Hamiltonian of the L = 20, N = 18 configuration. It contains 2L = 40 unjammed configurations, which are connected either to 2 or 4 other configurations. Amplitudes are −1/2 over all the off-diagonal matrix elements. 5.1. Mobile doublons The sector N = L − 2 sho…
Figure 5
Figure 5. Figure 5: Energy vs momentum of the doublons eigenstates, in the case N = 15, L = 17. Numerics (dots), and analytics (lines) in Eq.(39). together, in the second they oscillate around a common center of mass): ϵk = 1 4  −2 cos(k) ± √ 2 p cos(2k) + 7 + 6 . (39) At small k we hav…
Figure 6
Figure 6. Figure 6: Diffusion coefficient vs Q = 2π L for ρ → 0 with quadratic fit. L = 50, 60, · · · 990 for 1 particle (blue). L = 50, 55, · · · 395 for 2 particles(red). L = 50, 55, · · · 195 for 3 particles(green). And L = 50, 55, · · · 95 for 4 particles(black). They all extrapolate …
Figure 7
Figure 7. Figure 7: (Left) Diffusion coefficient vs Q = 2π L with a quadratic fit and the values predicted by MFT D = 3(1 − ρ). L = 8, 10, · · · 30 for ρ = 1/2 (blue). L = 9, 12, · · · 33 for ρ = 1/3 (red). L = 12, 16, · · · 36 for ρ = 1/4 (green). L = 15, 20, · · · 40 for ρ = 1/5(orange)…
Figure 8
Figure 8. Figure 8: (Left) Diffusion coefficient vs Q = 2π L for ρ → 1 with quadratic fit. L = 50, 60, · · · 490 for 2 holes (blue). L = 50, 60, · · · 200 for 3 holes(red). L = 50, 55, · · · 115 for 4 holes(green). (Right) For n ≥ 3 holes the coefficient cn defined ad cn = LDN−n,L vs 1/L …

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

56 extracted references · 47 canonical work pages

  1. [1]

    A kinetically constrained model exhibiting non-linear diffusion and jamming

    Introduction In recent years, we have witnessed a shift of attention in the field of statistical physics from the study of thermodynamic, equilibrium properties, which have been the core of the subject since its inception, to the characterization of the approach to equilibrium, or lack thereof. This is due to theoretical and experimental advances, both in...

  2. [2]

    The mean-field solution D = 3(1 − ρ), where ρ is the density of particles, agrees surprisingly well with the numerics in a very large range of densities, see Fig.1, which makes us suspect it is an exact solution. This is more surprising if we consider that the model exhibits the presence of jammed configurations, which appear at a critical density ρ = 2 /...

  3. [3]

    Preliminaries – model, mapping, definitions 2.1. The model and a classical-quantum mapping The random process studied describes random updates of 3 particles configurations where the only mobile configurations are those in which a particle is free on both sides. The rationale behind this choice[38] is to have a simple explicit model of random walkers with...

  4. [4]

    (11) This form is also reminiscent of a PXP spin model on a triangular ladder §

    Fermionization and a mean field approximation In order to proceed with the analysis of the model, we break hi,j,k into Pauli basis‡ and first rewrite the rate matrix as the Hamiltonian of a spin system with sα = σα/2 where σα are Pauli matrices: h1,2,3 = 1 2 − sz 1 ⃗ s2 · ⃗ s3 + 1 2 − sz 2 ⃗ s3 · ⃗ s1 + + 1 2 − sz 3 ⃗ s1 · ⃗ s2 + 1 3 sz 1 + sz 2 + sz 3 − ...

  5. [5]

    conditioned

    (13) The full Hamiltonian is obtained by summing these terms over all the triangular plaquettes. The resulting fermionic Hamiltonian then can be split into a “conditioned” hopping on the triangle in Fig.2: T = 1 2 X ⟨i,j,m⟩ a† i (1 − a† jaj)am + h.c., (14) where the notation ⟨i, j, m⟩ denotes three sites that belong to a triangle, by analogy with standard...

  6. [6]

    They are trivial ground states of the Hamiltonian since H |σ⟩ = 0 for each of these

    Large densities and the appearance of jammed configurations As we increase the density past the critical value ρ = 2 /3 a series of particle configurations appear which are stuck or jammed, namely that cannot be moved by our dynamical rules. They are trivial ground states of the Hamiltonian since H |σ⟩ = 0 for each of these. At exactly ρ = 2 /3 one can se...

  7. [7]

    Dynamics near maximum density: holes and doublons (and more-ons) We now turn to exploiting spectral methods to compute the dynamics of few hole excitations of the fully occupied inert state of N = L particles ( Nh ≡ L − N ). To outline the strategy: (i) recall, that at low hole density ϵ = 1 − ρ the leading contribution to thermodynamics and dynamics is f...

  8. [8]

    National Centre for HPC, Big Data and Quantum Computing (HPC)CN00000013 (CUP D43C22001240001) [MUR Decree n. 341- 15/03/2022] - Cascade Call launched by SPOKE 10 POLIMI: “CQEB

    Acknowledgements A.S. would like to thank the Graduate Center of CUNY for hospitality during a visit in July 2024, when this project got started. The work of A.S. was funded by the European Union - NextGenerationEU under the project NRRP “National Centre for HPC, Big Data and Quantum Computing (HPC)CN00000013 (CUP D43C22001240001) [MUR Decree n. 341- 15/0...

Show all 56 references
  1. [9]

    it contains a single particle, makes a random move of the particle on one of the other two, empty, vertices (if it does not qualify, then nothing happens)

    Exact numerical results The classical model can be simulated by a random process in which on picks a random triangular plaquette and then, if it qualifies, i.e. it contains a single particle, makes a random move of the particle on one of the other two, empty, vertices (if it d...

  2. [10]

    The solution is obtained by a classical- quantum mapping to a model of interacting fermions

    Conclusions and further work We have designed and solved, in the mean-field approximation, a kinetically constrained model of particles hopping on a triangular ladder. The solution is obtained by a classical- quantum mapping to a model of interacting fermions. The diffusion co...

  3. [11]

    J. R. Dorfman and E. G. D. Cohen. Velocity correlation functions in two and three dimensions. Phys. Rev. Lett. , 25:1257–1260, Nov 1970

  4. [12]

    Colloquium: Many-body localization, thermalization, and entanglement

    Dmitry A Abanin, Ehud Altman, Immanuel Bloch, and Maksym Serbyn. Colloquium: Many-body localization, thermalization, and entanglement. Reviews of Modern Physics, 91(2):021001, 2019

  5. [13]

    Metal–insulator transition in a weakly interacting many-electron system with localized single-particle states

    Denis M Basko, Igor L Aleiner, and Boris L Altshuler. Metal–insulator transition in a weakly interacting many-electron system with localized single-particle states. Annals of physics , 321(5):1126–1205, 2006

  6. [14]

    The physics of jamming for granular materials: a review

    Robert P Behringer and Bulbul Chakraborty. The physics of jamming for granular materials: a review. Reports on Progress in Physics , 82(1):012601, 2018

  7. [15]

    Stochastic interacting particle systems out of equilibrium

    Lorenzo Bertini, Alberto De Sole, Davide Gabrielli, Giovanni Jona-Lasinio, and C233569507120473 Landim. Stochastic interacting particle systems out of equilibrium. Journal of Statistical Mechanics: Theory and Experiment , 2007(07):P07014, 2007

  8. [16]

    P. M. Chaikin and T. C. Lubensky. Principles of Condensed Matter Physics. Cambridge University Press, 1995

  9. [17]

    together, in the second they oscillate around a common center of mass): ϵk = 1 4 −2 cos(k) ± √ 2 p cos(2k) + 7 + 6

    Numerics (dots), and analytics (lines) in Eq.(39). together, in the second they oscillate around a common center of mass): ϵk = 1 4 −2 cos(k) ± √ 2 p cos(2k) + 7 + 6 . (39) At small k we have ϵk,1 = 3 8 k2 + O(k4), (40) ϵk,2 = 2 + 1 8 k2 + O(k4). (41) Remembering that the minu...

  10. [18]

    From quantum chaos and eigenstate thermalization to statistical mechanics and thermodynamics

    Luca D’Alessio, Yariv Kafri, Anatoli Polkovnikov, and Marcos Rigol. From quantum chaos and eigenstate thermalization to statistical mechanics and thermodynamics. Advances in Physics , 65(3):239–362, 2016

  11. [19]

    Delacretaz

    Luca V. Delacretaz. Heavy operators and hydrodynamic tails. SciPost Phys. , 9:034, 2020

  12. [20]

    B. Derrida. An exactly soluble non-equilibrium system: The asymmetric simple exclusion process. Physics Reports, 301(1):65–83, 1998

  13. [21]

    Eigenstate thermalization hypothesis

    Joshua M Deutsch. Eigenstate thermalization hypothesis. Reports on Progress in Physics , 81(8):082001, 2018

  14. [22]

    Jamming in hard sphere and disk packings

    Aleksandar Donev, Salvatore Torquato, Frank H Stillinger, and Robert Connelly. Jamming in hard sphere and disk packings. Journal of applied physics , 95(3):989–999, 2004. + We thank Frank Pollman for this suggestion, which opens the way to an unexpected generalization of our w...

  15. [23]

    Theory of spin glasses

    Samuel Frederick Edwards and Phil W Anderson. Theory of spin glasses. Journal of Physics F: Metal Physics , 5(5):965, 1975

  16. [24]

    Long time tails in stationary random media

    MH Ernst, J Machta, JR Dorfman, and H Van Beijeren. Long time tails in stationary random media. i. theory. Journal of statistical physics , 34:477–495, 1984

  17. [25]

    Hydrodynamics of stationary non- equilibrium states for some stochastic lattice gas models

    Gregory Eyink, Joel L Lebowitz, and Herbert Spohn. Hydrodynamics of stationary non- equilibrium states for some stochastic lattice gas models. Communications in mathematical physics, 132(1):253–283, 1990

  18. [26]

    Quantum simulation

    Iulia M Georgescu, Sahel Ashhab, and Franco Nori. Quantum simulation. Reviews of Modern Physics, 86(1):153–185, 2014

  19. [27]

    Quantum simulations with ultracold atoms in optical lattices

    Christian Gross and Immanuel Bloch. Quantum simulations with ultracold atoms in optical lattices. Science, 357(6355):995–1001, 2017

  20. [28]

    Models of interacting bosons with exact ground states: a unified approach

    Zhaoyu Han and Steven A Kivelson. Models of interacting bosons with exact ground states: a unified approach. arXiv preprint arXiv:2408.15319 , 2024

  21. [29]

    Three-dimensional localization of ultracold atoms in an optical disordered potential

    Fred Jendrzejewski, Alain Bernard, Killian Mueller, Patrick Cheinet, Vincent Josse, Marie Piraud, Luca Pezz´ e, Laurent Sanchez-Palencia, Alain Aspect, and Philippe Bouyer. Three-dimensional localization of ultracold atoms in an optical disordered potential. Nature Physics, 8(...

  22. [30]

    T. H. Johnson, S. R. Clark, and D. Jaksch. Dynamical simulations of classical stochastic systems using matrix product states. Phys. Rev. E , 82:036702, Sep 2010

  23. [31]

    Nonequilibrium steady states of stochastic lattice gas models of fast ionic conductors

    Sheldon Katz, Joel L Lebowitz, and Herbert Spohn. Nonequilibrium steady states of stochastic lattice gas models of fast ionic conductors. Journal of statistical physics , 34(3):497–537, 1984

  24. [32]

    Critical behavior in the satisfiability of random boolean expressions

    Scott Kirkpatrick and Bart Selman. Critical behavior in the satisfiability of random boolean expressions. Science, 264(5163):1297–1301, 1994

  25. [33]

    Superconducting qubits: Current state of play

    Morten Kjaergaard, Mollie E Schwartz, Jochen Braum¨ uller, Philip Krantz, Joel I-J Wang, Simon Gustavsson, and William D Oliver. Superconducting qubits: Current state of play. Annual Review of Condensed Matter Physics , 11(1):369–395, 2020

  26. [34]

    Dynamics of repulsion processes

    PL Krapivsky. Dynamics of repulsion processes. Journal of Statistical Mechanics: Theory and Experiment, 2013(06):P06012, 2013

  27. [35]

    Fluctuations of current in nonstationary diffusive lattice gases

    PL Krapivsky and Baruch Meerson. Fluctuations of current in nonstationary diffusive lattice gases. Physical Review E—Statistical, Nonlinear, and Soft Matter Physics , 86(3):031106, 2012

  28. [36]

    Glasses and replicas.Structural Glasses and Supercooled Liquids: Theory, Experiment, and Applications , pages 151–191, 2012

    Marc M´ ezard and Giorgio Parisi. Glasses and replicas.Structural Glasses and Supercooled Liquids: Theory, Experiment, and Applications , pages 151–191, 2012

  29. [37]

    Spin glass theory and beyond: An Introduction to the Replica Method and Its Applications , volume 9

    Marc M´ ezard, Giorgio Parisi, and Miguel Angel Virasoro. Spin glass theory and beyond: An Introduction to the Replica Method and Its Applications , volume 9. World Scientific Publishing Company, 1987

  30. [38]

    Analytic and algorithmic solution of random satisfiability problems

    Marc M´ ezard, Giorgio Parisi, and Riccardo Zecchina. Analytic and algorithmic solution of random satisfiability problems. Science, 297(5582):812–815, 2002

  31. [39]

    Random k-satisfiability problem: From an analytic solution to an efficient algorithm

    Marc M´ ezard and Riccardo Zecchina. Random k-satisfiability problem: From an analytic solution to an efficient algorithm. Physical Review E , 66(5):056126, 2002

  32. [40]

    Michailidis, Dmitry A

    Alexios A. Michailidis, Dmitry A. Abanin, and Luca V. Delacr´ etaz. Corrections to diffusion in interacting quantum systems. Phys. Rev. X , 14:031020, Aug 2024

  33. [41]

    Statistical theory of transport by strongly interacting lattice fermions

    Subroto Mukerjee, Vadim Oganesyan, and David Huse. Statistical theory of transport by strongly interacting lattice fermions. Phys. Rev. B , 73:035113, Jan 2006

  34. [42]

    Many-body localization and thermalization in quantum statistical mechanics

    Rahul Nandkishore and David A Huse. Many-body localization and thermalization in quantum statistical mechanics. Annu. Rev. Condens. Matter Phys. , 6(1):15–38, 2015

  35. [43]

    Localization of interacting fermions at high temperature

    Vadim Oganesyan and David A Huse. Localization of interacting fermions at high temperature. Physical Review B—Condensed Matter and Materials Physics , 75(15):155111, 2007

  36. [44]

    Energy transport in disordered classical spin chains

    Vadim Oganesyan, Arijeet Pal, and David A Huse. Energy transport in disordered classical spin chains. Physical Review B—Condensed Matter and Materials Physics , 80(11):115104, 2009. A kinetically constrained model exhibiting non-linear diffusion and jamming 19

  37. [45]

    Nobel lecture: Multiple equilibria

    Giorgio Parisi. Nobel lecture: Multiple equilibria. Reviews of Modern Physics , 95(3):030501, 2023

  38. [46]

    Mean-field theory of hard sphere glasses and jamming

    Giorgio Parisi and Francesco Zamponi. Mean-field theory of hard sphere glasses and jamming. Reviews of Modern Physics , 82(1):789–845, 2010

  39. [47]

    Exact dynamical equations for kinetically-constrained- models

    Gianmarco Perrupato and Tommaso Rizzo. Exact dynamical equations for kinetically-constrained- models. arXiv preprint arXiv:2212.05132 , 2022

  40. [48]

    Facilitated spin models, mode coupling theory, and ergodic–nonergodic transitions

    Steven J Pitts, Thomas Young, and Hans C Andersen. Facilitated spin models, mode coupling theory, and ergodic–nonergodic transitions. The Journal of Chemical Physics , 113(19):8671– 8679, 2000

  41. [49]

    Diffusion cascade in a model of interacting random walkers

    Abhishek Raj, Paolo Glorioso, Sarang Gopalakrishnan, and Vadim Oganesyan. Diffusion cascade in a model of interacting random walkers. arXiv preprint arXiv:2412.05222 , 2024. https: //arxiv.org/abs/2412.05222

  42. [50]

    Ritort and P

    F. Ritort and P. Sollich. Glassy dynamics of kinetically constrained models. Advances in Physics, 52(4):219–342, 2003

  43. [51]

    An analytical approach to the fredrickson–andersen model in one dimension

    Michael Schulz and Steffen Trimper. An analytical approach to the fredrickson–andersen model in one dimension. International Journal of Modern Physics B , 11(24):2927–2940, 1997

  44. [52]

    Constraint-induced delocalization

    Piotr Sierant, Eduardo Gonzalez Lazo, Marcello Dalmonte, Antonello Scardicchio, and Jakub Zakrzewski. Constraint-induced delocalization. Physical Review Letters, 127(12):126603, 2021

  45. [53]

    Many-body localization in the age of classical computing

    Piotr Sierant, Maciej Lewenstein, Antonello Scardicchio, Lev Vidmar, and Jakub Zakrzewski. Many-body localization in the age of classical computing. arXiv preprint arXiv:2403.07111 , 2024

  46. [54]

    Large scale dynamics of interacting particles

    Herbert Spohn. Large scale dynamics of interacting particles. Springer Science & Business Media, 2012

  47. [55]

    Chaos and quantum thermalization

    Mark Srednicki. Chaos and quantum thermalization. Physical review e , 50(2):888, 1994

  48. [56]

    Quantum field theory of many-body systems: From the origin of sound to an origin of light and electrons

    Xiao-Gang Wen. Quantum field theory of many-body systems: From the origin of sound to an origin of light and electrons . Oxford university press, 2004

Pith tools

Reviewed August 11, 2026 · model on record in the stance chip above.