SLIDE 1 Absorbing-state phase transitions
Leonardo T. Rolla
Argentinian National Research Council at the University of Buenos Aires NYU-ECNU Institute of Mathematical Sciences at NYU-Shanghai XXIII Brazilian School on Probability – July, 2019
SLIDE 2
Background
SLIDE 3
Context
Non-equilibrium Statistical Mechanics
SLIDE 4
Context
Non-equilibrium Statistical Mechanics Critical phenomena
self-similar shapes large fluctuations long-range correlations avalanches
SLIDE 5
Critical phenomena
Many physical systems behave as: β < βc → well behaved, short-range correlations
SLIDE 6
Critical phenomena
Many physical systems behave as: β < βc → well behaved, short-range correlations β = βc → long range, power laws, intricate behavior
SLIDE 7
Critical phenomena
Many physical systems behave as: β < βc → well behaved, short-range correlations β = βc → long range, power laws, intricate behavior β > βc → well behaved, short-range correlations
SLIDE 8
Critical phenomena
Many physical systems behave as: β < βc → well behaved, short-range correlations β = βc → long range, power laws, intricate behavior β > βc → well behaved, short-range correlations But why do we observe critical behavior outside a controlled environment?
SLIDE 9
Critical phenomena (cont)
Physics literature: Late 80’s – Self-organized criticality
A system that finds a critical state all by itself
SLIDE 10
Critical phenomena (cont)
Physics literature: Late 80’s – Self-organized criticality
A system that finds a critical state all by itself
Late 90’s – Relate it to ordinary phase transitions
Absorbing-state phase transitions
SLIDE 11
Driven-dissipative dynamics
SLIDE 12
Driven-dissipative dynamics
SLIDE 13
Driven-dissipative dynamics
SLIDE 14
Driven-dissipative dynamics
SLIDE 15
Driven-dissipative dynamics
SLIDE 16
Driven-dissipative dynamics
Infinite-volume conservative system
SLIDE 17
Driven-dissipative dynamics
Infinite-volume conservative system
System goes to an absorbing state if the density of particles is below ζc, and remains unstable if the density is above ζc
SLIDE 18
Model and predictions
SLIDE 19
Activated Random Walks
SLIDE 20
Assumptions
◮ Jumps to nearest-neighbors
SLIDE 21
Assumptions
◮ Jumps to nearest-neighbors ◮ Particles start active
SLIDE 22
Assumptions
◮ Jumps to nearest-neighbors ◮ Particles start active ◮ At t = 0, i.i.d. Poisson(ζ) particles
SLIDE 23
Assumptions
◮ Jumps to nearest-neighbors ◮ Particles start active ◮ At t = 0, i.i.d. Poisson(ζ) particles ◮ 0 < λ ∞
SLIDE 24
Fixation vs activity
Fixation: each site is eventually stable Activity: each site is visited infinitely many times
SLIDE 25
Fixation vs activity
Fixation: each site is eventually stable Activity: each site is visited infinitely many times Dichotomy: either fixation a.s. or activity a.s.
SLIDE 26
Fixation vs activity
Fixation: each site is eventually stable Activity: each site is visited infinitely many times Dichotomy: either fixation a.s. or activity a.s. Monotonicity:
SLIDE 27
Predictions
ζ λ 1 ∞ ∞
Activity Fixation
SLIDE 28 Predictions
ζ λ 1 ∞ ∞
Activity Fixation
ζc(λ) should not depend
SLIDE 29 Predictions
ζ λ 1 ∞ ∞
Activity Fixation
ζc(λ) should not depend
At and near ζ = ζc: no fixation power laws rich scaling limits bursts of activity
SLIDE 30
Difficulties
Lack of attractiveness
Overcome by using constructions other than Harris’
SLIDE 31
Difficulties
Lack of attractiveness
Overcome by using constructions other than Harris’
Conservation of particles
Rules out “energy vs. entropy” approaches
SLIDE 32
Phase transition results
SLIDE 33 d = 1 directed d = 1 biased d = 1 unbiased d = 2 unbiased d ≥ 3 unbiased d ≥ 2 biased ζ λ 1 ∞ ζ λ 1 ∞ ζ λ 1 ∞ ζ λ 1 ∞ ζ λ 1 ∞ ζ λ 1 ∞
SLIDE 34 d = 1 directed d = 1 biased d = 1 unbiased d = 2 unbiased d ≥ 3 unbiased d ≥ 2 biased ζ λ 1 ∞ ζ λ 1 ∞ ζ λ 1 ∞ ζ λ 1 ∞ ζ λ 1 ∞ ζ λ 1 ∞
Hoffman, Sidoravicius. Unpublished (2004) | Cabezas, R, Sidoravicius. J Stat Phys (2014)
SLIDE 35 d = 1 directed d = 1 biased d = 1 unbiased d = 2 unbiased d ≥ 3 unbiased d ≥ 2 biased ζ λ 1 ∞ ζ λ 1 ∞ ζ λ 1 ∞ ζ λ 1 ∞ ζ λ 1 ∞ ζ λ 1 ∞
R, Sidoravicius. Invent Math (2012)
SLIDE 36 d = 1 directed d = 1 biased d = 1 unbiased d = 2 unbiased d ≥ 3 unbiased d ≥ 2 biased ζ λ 1 ∞ ζ λ 1 ∞ ζ λ 1 ∞ ζ λ 1 ∞ ζ λ 1 ∞ ζ λ 1 ∞
- Shellef. ALEA (2010) | Amir, Gurel-Gurevich. Electron Commun Probab (2010)
SLIDE 37 d = 1 directed d = 1 biased d = 1 unbiased d = 2 unbiased d ≥ 3 unbiased d ≥ 2 biased scaling limit ζ λ 1 ∞ ζ λ 1 ∞ ζ λ 1 ∞ ζ λ 1 ∞ ζ λ 1 ∞ ζ λ 1 ∞
Cabezas, R, Sidoravicius. J Stat Phys (2014), Probab Theory Relat Fields (2018)
SLIDE 38 d = 1 directed d = 1 biased d = 1 unbiased d = 2 unbiased d ≥ 3 unbiased d ≥ 2 biased scaling limit ζ λ 1 ∞ ζ λ 1 ∞ ζ λ 1 ∞ ζ λ 1 ∞ ζ λ 1 ∞ ζ λ 1 ∞
- Taggi. Electron J Probab (2016)
SLIDE 39 d = 1 directed d = 1 biased d = 1 unbiased d = 2 unbiased d ≥ 3 unbiased d ≥ 2 biased scaling limit ζ λ 1 ∞ ζ λ 1 ∞ ζ λ 1 ∞ ζ λ 1 ∞ ζ λ 1 ∞ ζ λ 1 ∞
Sidoravicius, Teixeira. Electron J Proba (2017)
SLIDE 40 d = 1 directed d = 1 biased d = 1 unbiased d = 2 unbiased d ≥ 3 unbiased d ≥ 2 biased scaling limit ζ λ 1 ∞ ζ λ 1 ∞ ζ λ 1 ∞ ζ λ 1 ∞ ζ λ 1 ∞ ζ λ 1 ∞
R, Tournier. Ann Inst H Poincar´ e Probab Statist (2018)
SLIDE 41 d = 1 directed d = 1 biased d = 1 unbiased d = 2 unbiased d ≥ 3 unbiased d ≥ 2 biased scaling limit ζ λ 1 ∞ ζ λ 1 ∞ ζ λ 1 ∞ ζ λ 1 ∞ ζ λ 1 ∞ ζ λ 1 ∞
Basu, Ganguly, Hoffman. Comm Math Phys (2018)
SLIDE 42 d = 1 directed d = 1 biased d = 1 unbiased d = 2 unbiased d ≥ 3 unbiased d ≥ 2 biased scaling limit ζ λ 1 ∞ ζ λ 1 ∞ ζ λ 1 ∞ ζ λ 1 ∞ ζ λ 1 ∞ ζ λ 1 ∞
Stauffer, Taggi. Ann Probab (2018)
SLIDE 43 slow d = 1 directed d = 1 biased d = 1 unbiased d = 2 unbiased d ≥ 3 unbiased d ≥ 2 biased fast scaling limit ζ λ 1 ∞ ζ λ 1 ∞ ζ λ 1 ∞ ζ λ 1 ∞ ζ λ 1 ∞ ζ λ 1 ∞
Basu, Ganguly, Hoffman, Richey. Ann Inst H Poincar´ e Probab Statist (2019+)
SLIDE 44 slow d = 1 directed d = 1 biased d = 1 unbiased d = 2 unbiased d ≥ 3 unbiased d ≥ 2 biased fast scaling limit ζ λ 1 ∞ ζ λ 1 ∞ ζ λ 1 ∞ ζ λ 1 ∞ ζ λ 1 ∞ ζ λ 1 ∞
- Taggi. Ann Inst H Poincar´
e Probab Statist (2019+)
SLIDE 45 slow d = 1 directed d = 1 biased d = 1 unbiased d = 2 unbiased d ≥ 3 unbiased d ≥ 2 biased fast scaling limit ζ λ 1 ∞ ζ λ 1 ∞ ζ λ 1 ∞ ζ λ 1 ∞ ζ λ 1 ∞ ζ λ 1 ∞
Asselah, R, Schapira. Writing up
SLIDE 46 slow d = 1 directed d = 1 biased d = 1 unbiased d = 2 unbiased d ≥ 3 unbiased d ≥ 2 biased fast slow scaling limit ζ λ 1 ∞ ζ λ 1 ∞ ζ λ 1 ∞ ζ λ 1 ∞ ζ λ 1 ∞ ζ λ 1 ∞
Cabezas, R. Writing up
SLIDE 47 slow d = 1 directed d = 1 biased d = 1 unbiased d = 2 unbiased d ≥ 3 unbiased d ≥ 2 biased fast slow scaling limit ζ λ 1 ∞ ζ λ 1 ∞ ζ λ 1 ∞ ζ λ 1 ∞ ζ λ 1 ∞ ζ λ 1 ∞
Hoffman, Richey, R. Writing up
SLIDE 48 slow d = 1 directed d = 1 biased d = 1 unbiased d = 2 unbiased d ≥ 3 unbiased d ≥ 2 biased fast slow scaling limit ζ λ 1 ∞ ζ λ 1 ∞ ζ λ 1 ∞ ζ λ 1 ∞ ζ λ 1 ∞ ζ λ 1 ∞
SLIDE 49
Constructions and main tools
SLIDE 50
Constructions
Harris graphical construction
clocks with marks at each site
SLIDE 51
Constructions
Harris graphical construction
clocks with marks at each site
Site-wise construction
stack of instructions at each site
SLIDE 52
Constructions
Harris graphical construction
clocks with marks at each site
Site-wise construction
stack of instructions at each site
Particle-wise constructions
particles start with a life plan and do pause/resume
SLIDE 53 Main tools
Site-wise representation and Abelianess
Diaconis, Fulton. Rend Semin Mat Torino (1991) | Eriksson. SIAM J Discrete Math (1996)
SLIDE 54 Main tools
Site-wise representation and Abelianess
Diaconis, Fulton. Rend Semin Mat Torino (1991) | Eriksson. SIAM J Discrete Math (1996)
Relate it to the dynamics, preserving monotonicity, etc Reduce fixation-activity question to toppling procedures
R, Sidoravicius. Invent Math (2012)
SLIDE 55 Main tools (cont)
Assuming a particle-wise construction is well-defined: a particle stays active ⇒ sites stay active
Amir, Gurel-Gurevich. Electron Commun Probab (2010)
SLIDE 56 Main tools (cont)
Assuming a particle-wise construction is well-defined: a particle stays active ⇒ sites stay active
Amir, Gurel-Gurevich. Electron Commun Probab (2010)
Well-definedness of the particle-wise construction → ergodicity, mass transport, coupling, surgery Averaged criterion for activity
R, Tournier. Ann Inst H Poincar´ e Probab Statist (2018)
SLIDE 57 Criterion for activity
⊢ Stabilizing a large box forces a large number of particles to visit a specific site, wpp
R, Sidoravicius. Invent Math (2012)
SLIDE 58 Criterion for activity
⊢ Stabilizing a large box forces a large number of particles to visit a specific site, wpp
R, Sidoravicius. Invent Math (2012)
⊢ Stabilizing a large box forces a positive fraction of the particles to leave the box, on average
R, Tournier. Ann Inst H Poincar´ e Probab Statist (2018)
SLIDE 59
Sharpness
SLIDE 60
- Theorem. Given d, λ, p, there exists ζc such that
ζ ζc
fixation activity
for all ergodic initial states with density ζ
R, Sidoravicius, Zindy. Ann Henri Poincar´ e (2019)
SLIDE 61
- Theorem. Given d, λ, p, there exists ζc such that
ζ ζc
fixation activity
for all ergodic initial states with density ζ
R, Sidoravicius, Zindy. Ann Henri Poincar´ e (2019)
– Can drop the previous “i.i.d. Poisson” assumption
SLIDE 62
- Theorem. Given d, λ, p, there exists ζc such that
ζ ζc
fixation activity
for all ergodic initial states with density ζ
R, Sidoravicius, Zindy. Ann Henri Poincar´ e (2019)
– Can drop the previous “i.i.d. Poisson” assumption – Restrictive proofs now yield general theorems
SLIDE 63
- Theorem. Given d, λ, p, there exists ζc such that
ζ ζc
fixation activity
for all ergodic initial states with density ζ
R, Sidoravicius, Zindy. Ann Henri Poincar´ e (2019)
– Can drop the previous “i.i.d. Poisson” assumption – Restrictive proofs now yield general theorems – Contributes to ongoing discussion about some dissipative models mixing better than others
SLIDE 64
Scaling limits and avalanches
SLIDE 65
Critical one-dimensional model
d = 1, directed walks (integrable case) i.i.d. initial condition with critical density ζ = ζc =
λ 1+λ
SLIDE 66
Critical one-dimensional model
d = 1, directed walks (integrable case) i.i.d. initial condition with critical density ζ = ζc =
λ 1+λ
Run the dynamics on Vn = [0, n] until it is stable. Cn := how many particles cross the origin.
SLIDE 67
Critical one-dimensional model
d = 1, directed walks (integrable case) i.i.d. initial condition with critical density ζ = ζc =
λ 1+λ
Run the dynamics on Vn = [0, n] until it is stable. Cn := how many particles cross the origin. Released the active particles at x = n + 1, let them interact with the leftovers of previous step. Cn+1 ≥ Cn
SLIDE 68
Simulation
SLIDE 69
SLIDE 70
Scaling limit
Theorem (Cabezas, R ’19). The counting process (Cn)n0 rescales to (Cρ
x)x0
Pure-jump process constructed from a collection of correlated reflected coalescing Brownian motions
SLIDE 71
Scaling limit
Theorem (Cabezas, R ’19). The counting process (Cn)n0 rescales to (Cρ
x)x0
Pure-jump process constructed from a collection of correlated reflected coalescing Brownian motions Correlations depend on ρ = σs
σp ∈ (0, 1], where
σ2
s = ζ − ζ2 and σ2 p is the variance of initial condition
SLIDE 72
ρ = 1.00, 0.50, 0.30, 0.10, 0.05, 0.00
SLIDE 73
Notes on some selected proofs
SLIDE 74
Toppling procedures
Abelian particle-wise construction mV (x) := odometer at x when V is stabilized lim
k lim V P[mV (x) k] =
0 ⇔ Fixation (B) 1 ⇔ Activity (U) Mn := #Particles which quit when Bn is stabilized lim sup
n
E[Mn] |Bn| > 0 ⇒ Activity (E)
SLIDE 75
Examples
Thm (Stauffer, Taggi). ζc
λ 1+λ
Thm (Taggi; R, Tournier). ζc Fp(λ)
SLIDE 76 Particle-wise construction
Labeled particles → mass transport, ergodicity, surgery Construction: assign to each particle a CTRW+beep MTP Example. Assume particles fixate a.s. ζ =E[start at o] = E[settle at o] 1
- Thm. Site fixation ⇒ ζ < 1
Cabezas, R, Sidoravicius. Probab Theory Relat Fields (2018)
SLIDE 77 Particle-wise construction (cont)
- Thm. The PWC is well-defined
Add particles one by one, updating the whole evolution ⊢ Life of each particle is well-defined through some limit Main step: ∀x, T, the number of particle additions that affect site x by time T has finite expectation
R, Tournier. Ann Inst H Poincar´ e Probab Statist (2018)
SLIDE 78
Open problems and questions
SLIDE 79
Open problems and questions
Proofs of fixation/activity that give different insights
SLIDE 80
Open problems and questions
Proofs of fixation/activity that give different insights Unbiased walks on Z2: ζc < 1 for some λ < ∞
SLIDE 81
Open problems and questions
Proofs of fixation/activity that give different insights Unbiased walks on Z2: ζc < 1 for some λ < ∞ Dichotomy for the slow-fast transition on finite ring
SLIDE 82
Open problems and questions
Proofs of fixation/activity that give different insights Unbiased walks on Z2: ζc < 1 for some λ < ∞ Dichotomy for the slow-fast transition on finite ring Make sense of scaling limit at criticality for d 2
SLIDE 83
Open problems and questions
Proofs of fixation/activity that give different insights Unbiased walks on Z2: ζc < 1 for some λ < ∞ Dichotomy for the slow-fast transition on finite ring Make sense of scaling limit at criticality for d 2 Sharpness when a fraction of particles start active
SLIDE 84
Open problems and questions
Proofs of fixation/activity that give different insights Unbiased walks on Z2: ζc < 1 for some λ < ∞ Dichotomy for the slow-fast transition on finite ring Make sense of scaling limit at criticality for d 2 Sharpness when a fraction of particles start active Many more...