Recognition: unknown
Hydrodynamic limits for TASEP with space-time discontinuities
Pith reviewed 2026-05-14 18:20 UTC · model grok-4.3
The pith
A height-dependent TASEP with space-time discontinuous jump rates has a hydrodynamic limit given by a Lax-Oleinik variational formula for the current.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Combining the law of large numbers for the associated inhomogeneous directed last-passage percolation model with an extension of the variational coupling method, the authors prove a hydrodynamic limit for the height function and associated particle density. The limiting current is characterised by a Lax-Oleinik type variational formula built from the discontinuous last-passage shape function. At points of differentiability and continuity of the coefficient, this current solves a Hamilton-Jacobi equation with discontinuous dependence on space and height. At discontinuities, it satisfies a natural envelope-based discontinuous viscosity formulation. For coefficients with only spatial discontinu
What carries the argument
The Lax-Oleinik variational formula constructed from the discontinuous last-passage shape function, which selects the limiting current and determines the PDE structure.
If this is right
- The particle density follows the maximal-current weak solution of the conservation law with discontinuous flux when discontinuities are only spatial.
- The limiting current is the unique nondecreasing Lipschitz solution to the envelope-based discontinuous viscosity Hamilton-Jacobi equation.
- The microscopic dynamics naturally select the envelope solution at discontinuity curves of the speed function.
Where Pith is reading between the lines
- This method could be adapted to other particle systems or growth models with space-time varying rates.
- Testing the rate of convergence near discontinuity curves in simulations would provide further evidence.
- Applications may include modeling traffic or biological transport with abrupt environmental changes.
Load-bearing premise
The speed function has discontinuities along only locally finitely many curves so that the variational coupling and envelope formulation can be controlled at the jumps.
What would settle it
A simulation of the process with a speed function discontinuous across one space-time curve in which the observed height function deviates from the predicted Lax-Oleinik variational value would falsify the limit.
Figures
read the original abstract
We develop a hydrodynamic theory for a height-dependent version of the totally asymmetric simple exclusion process in which the jump rate at a growth site is sampled from a macroscopic two-dimensional speed function evaluated at the spatial coordinate and the current height level. The speed function is allowed to have discontinuities along locally finitely many curves. Through the TASEP height-function representation, the process is coupled to an inhomogeneous directed last-passage percolation model whose exponential rates vary discontinuously in the two macroscopic LPP coordinates. Combining the law of large numbers for this last-passage model with an extension of the variational coupling method, we prove a hydrodynamic limit for the height function and for the associated particle density. The limiting current is characterised by a Lax-Oleinik type variational formula built from the discontinuous last-passage shape function. We then identify the first-order PDE structure selected by the microscopic dynamics. At points of differentiability of the limiting current and continuity of the sampled coefficient, the current solves a Hamilton-Jacobi equation whose Hamiltonian depends discontinuously on the spatial variable and on the value of the solution itself. At discontinuities, the variational formula leads to a natural envelope-based discontinuous viscosity formulation, and we prove that the limiting current satisfies this formulation. Finally, when the coefficient has only spatial discontinuities, we prove uniqueness of the Hamilton-Jacobi solution in the natural class of nondecreasing Lipschitz currents, and identify its spatial derivative as the maximal-current weak solution of the associated scalar conservation law with discontinuous flux.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper proves a hydrodynamic limit for the height function (and associated density) of a height-dependent TASEP whose jump rates are sampled from a macroscopic speed function that may be discontinuous along locally finitely many curves in space-time. The argument couples the microscopic process to an inhomogeneous directed last-passage percolation model, invokes an external LLN for the LPP shape function, and extends the variational coupling method to pass to the limit; the limiting current is characterized by a Lax-Oleinik variational formula that yields a discontinuous-viscosity Hamilton-Jacobi equation whose Hamiltonian depends on both the spatial variable and the solution value itself. Uniqueness of the HJ solution (and identification of its derivative as the maximal-current weak solution of the associated conservation law) is established only when the discontinuities are purely spatial.
Significance. If the technical controls hold, the result supplies the first hydrodynamic theory for TASEP-type models whose coefficients are allowed to jump in both space and time. The combination of an inhomogeneous LPP law of large numbers with a variational-coupling argument, together with the envelope formulation that selects the correct weak solution at discontinuities, provides a reusable template for other particle systems with discontinuous rates. The explicit identification of the selected PDE structure is a clear strength.
major comments (2)
- [§5] §5 (variational coupling extension): local finiteness of the discontinuity curves is invoked to control the coupling and envelope formulation at jumps, yet the argument supplies no uniform bound on the number of crossings per path under hydrodynamic scaling; repeated crossings could permit o(1) errors to accumulate before the limit is taken, and this bound is load-bearing for the space-time case.
- [§6] §6 (uniqueness): uniqueness of the Hamilton-Jacobi solution in the class of nondecreasing Lipschitz currents is proved only for purely spatial discontinuities; the space-time case therefore rests entirely on the variational characterization without an independent uniqueness statement, leaving open the possibility that other solutions to the envelope formulation exist.
minor comments (2)
- [§2] The notation distinguishing the macroscopic speed function from the microscopic rate function is introduced only in the abstract and could be restated explicitly at the beginning of §2.
- Figure 1 (schematic of discontinuity curves) would benefit from a caption that explicitly labels the local-finiteness assumption.
Simulated Author's Rebuttal
We thank the referee for the careful reading and for identifying these two technical points that require clarification. We address each major comment below.
read point-by-point responses
-
Referee: [§5] §5 (variational coupling extension): local finiteness of the discontinuity curves is invoked to control the coupling and envelope formulation at jumps, yet the argument supplies no uniform bound on the number of crossings per path under hydrodynamic scaling; repeated crossings could permit o(1) errors to accumulate before the limit is taken, and this bound is load-bearing for the space-time case.
Authors: We agree that an explicit uniform bound on the number of crossings is needed for rigor. Because the discontinuity curves are locally finite and the variational paths are 1-Lipschitz, the number of crossings per path in any compact space-time region is bounded by a constant depending only on the local geometry and the Lipschitz constant. We will insert a short lemma establishing this bound under hydrodynamic scaling and verify that the resulting error remains o(1) uniformly; the revised manuscript will contain this estimate. revision: yes
-
Referee: [§6] §6 (uniqueness): uniqueness of the Hamilton-Jacobi solution in the class of nondecreasing Lipschitz currents is proved only for purely spatial discontinuities; the space-time case therefore rests entirely on the variational characterization without an independent uniqueness statement, leaving open the possibility that other solutions to the envelope formulation exist.
Authors: The referee correctly notes that uniqueness of the discontinuous-viscosity solution is proved only when discontinuities are purely spatial. In the space-time setting the limiting current is identified directly as the hydrodynamic limit of the microscopic height function, which is given by the Lax-Oleinik variational formula constructed from the LPP shape function. This variational object is the unique candidate selected by the particle system; however, we do not supply an independent uniqueness theorem for the envelope formulation when the Hamiltonian depends on both space and time. Such a result would require additional analytic tools and lies outside the scope of the present work. revision: no
- Uniqueness of the envelope formulation for the Hamilton-Jacobi equation with space-time discontinuities
Circularity Check
No significant circularity; derivation combines external LLN with variational coupling extension
full rationale
The paper's central result is obtained by combining a law of large numbers for the inhomogeneous last-passage percolation model with an extension of the variational coupling method, then verifying that the limit satisfies a Lax-Oleinik variational formula and a discontinuous viscosity formulation of the Hamilton-Jacobi equation. No step reduces a prediction to a fitted input by construction, nor does any load-bearing claim rest solely on a self-citation whose content is unverified within the paper. The local-finiteness assumption on discontinuity curves is an explicit hypothesis used to control crossings in the coupling, not a self-definitional loop. Uniqueness of the HJ solution is proven only in the spatial-discontinuity case; the space-time case is characterized variationally without claiming uniqueness, which is consistent with the stated scope. The derivation therefore remains self-contained against the external LLN benchmark and does not exhibit any of the enumerated circular patterns.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption Law of large numbers for inhomogeneous directed last-passage percolation with discontinuous rates
Reference graph
Works this paper leans on
-
[1]
Adimurthi and G. D. Veerappa Gowda,Conservation law with discontinuous flux, J. Math. Kyoto Univ.43(2003), no. 1, 27–70. doi:10.1215/kjm/1250283740
-
[2]
D. Aldous and P. Diaconis. Hammersley’s interacting particle process and longest increasing subsequences.Probability Theory and Related Fields, 103:199–213, 1995. doi:10.1007/BF01204214
-
[3]
G. Amir, C. Bahadoran, O. Busani, and E. Saada. Invariant measures for multilane exclusion process.Annales de l’Institut Henri Poincar´ e, Probabilit´ es et Statistiques, 61(3):2184–2234, 2025. doi:10.1214/24-AIHP1482
- [4]
-
[5]
B. Andreianov, K. H. Karlsen, and N. H. Risebro,A theory ofL 1-dissipative solvers for scalar conservation laws with discontinuous flux, Arch. Ration. Mech. Anal.201(2011), 27–86. doi:10.1007/s00205-010-0389-4
-
[6]
B. Andreianov and D. Mitrovi´ c,Entropy conditions for scalar conservation laws with discontinuous flux revisited, Ann. Inst. H. Poincar´ e C Anal. Non Lin´ eaire32(2015), no. 6, 1307–1335. doi:10.1016/j.anihpc.2014.08.002
-
[7]
E. Audusse and B. Perthame,Uniqueness for scalar conservation laws with discontinuous flux via adapted entropies, Proc. Roy. Soc. Edinburgh Sect. A135(2005), no. 2, 253–265. doi:10.1017/S0308210500003863
-
[8]
C. Bahadoran. Hydrodynamical limit for spatially heterogeneous simple exclusion processes.Probability Theory and Related Fields, 110:287–331, 1998. doi:10.1007/s004400050150
-
[9]
C. Bahadoran, H. Guiol, K. Ravishankar, and E. Saada. Euler hydrodynamics for attractive particle systems in random environment.Annales de l’Institut Henri Poincar´ e, Probabilit´ es et Statistiques, 50(2):403–424, 2014. doi:10.1214/12- AIHP510
work page doi:10.1214/12- 2014
-
[10]
J. Baik, P. Deift, and K. Johansson. On the distribution of the length of the longest increasing subsequence of random permutations.Journal of the American Mathematical Society, 12(4):1119–1178, 1999. doi:10.1090/S0894-0347-99-00307-0
-
[11]
Y. Bakhtin. The Burgers equation with Poisson random forcing.Annals of Probability, 41(4):2961–2989, 2013
2013
-
[12]
Bakhtin, E
Y. Bakhtin, E. Cator, and K. Khanin. Space-time stationary solutions for the Burgers equation.Journal of the American Mathematical Society, 27(1):193–238, 2014
2014
-
[13]
G. Barles and E. Chasseigne,(Almost) Everything you always wanted to know about deterministic control problems in stratified domains, Netw. Heterog. Media10(2015), no. 4, 809–836. doi:10.3934/nhm.2015.10.809
-
[14]
R. Basu, V. Sidoravicius, and A. Sly. Last passage percolation with a defect line and the solution of the slow bond problem.arXiv preprint, arXiv:1408.3464, 2016
work page internal anchor Pith review Pith/arXiv arXiv 2016
-
[15]
I. Benjamini, P. A. Ferrari, and C. Landim. Asymmetric conservative processes with random rates.Stochastic Processes and their Applications, 61(2):181–204, 1996. doi:10.1016/0304-4149(95)00077-1
-
[16]
L. Bertini and G. Giacomin. Stochastic Burgers and KPZ equations from particle systems.Communications in Mathe- matical Physics, 183:571–607, 1997. doi:10.1007/s002200050044
-
[17]
E. Bisi, Y. Liao, A. Saenz, and N. Zygouras. Non-intersecting path constructions for TASEP with inhomogeneous rates and the KPZ fixed point.Communications in Mathematical Physics, 402:285–333, 2023. doi:10.1007/s00220-023-04723-8
-
[18]
A. Borodin and I. Corwin. Macdonald processes.Probability Theory and Related Fields, 158:225–400, 2014. doi:10.1007/s00440-013-0482-3
-
[19]
A. Borodin, I. Corwin, and V. Gorin. Stochastic six-vertex model.Duke Mathematical Journal, 165(3):563–624, 2016. doi:10.1215/00127094-3166843
-
[20]
A. Borodin, I. Corwin, and T. Sasamoto. From duality to determinants forq-TASEP and ASEP.Annals of Probability, 42(6):2314–2382, 2014. doi:10.1214/13-AOP868
-
[21]
J. Calder. Directed last passage percolation with discontinuous weights.Journal of Statistical Physics, 158:903–949, 2015. doi:10.1007/s10955-014-1146-0
-
[22]
G.-Q. Chen, N. Even, and C. Klingenberg. Hyperbolic conservation laws with discontinuous fluxes and hydrodynamic limit for particle systems.Journal of Differential Equations, 245(11):3095–3126, 2008. doi:10.1016/j.jde.2008.07.036
-
[23]
X. Chen and B. Hu, Viscosity solutions of discontinuous Hamilton–Jacobi equations,Interfaces and Free Boundaries10 (2008), 339–359. doi:10.4171/IFB/192
-
[24]
D. Chowdhury, L. Santen, and A. Schadschneider. Statistical physics of vehicular traffic and some related systems. Physics Reports (2000), 329(4-6):199–329. doi:10.1016/S0370-1573(99)00117-9
-
[25]
F. Ciech and N. Georgiou. Last passage percolation in an exponential environment with discontinuous ratesAnnales de l’Institut Henri Poincar´ e, Probabilit´ es et Statistiques(2021), 57(4):2165–2188. doi:10.1214/21-AIHP1172
-
[26]
G. M. Coclite and N. H. Risebro,Viscosity solutions of Hamilton–Jacobi equations with discontinuous coefficients, J. Hyperbolic Differ. Equ.4(2007), no. 4, 771–795. doi:10.1142/S0219891607001355
-
[27]
I. Corwin. The Kardar–Parisi–Zhang equation and universality class.Random Matrices: Theory and Applications, 1(1):1130001, 2012. doi:10.1142/S2010326311300014
-
[28]
I. Corwin and L.-C. Tsai. SPDE limit of weakly inhomogeneous ASEP.Electronic Journal of Probability, 25:paper no. 156, 1–55, 2020. doi:10.1214/20-EJP565. SPACE-TIME DISCONTINUOUS TASEP 61
-
[29]
P. Covert and F. Rezakhanlou. Hydrodynamic limit for particle systems with nonconstant speed parameter.Journal of Statistical Physics, 88:383–426, 1997. doi:10.1007/BF02508477
-
[30]
G. Crasta, V. De Cicco, and G. De Philippis,Kinetic formulation and uniqueness for scalar conservation laws with discontinuous flux, Comm. Partial Differential Equations40(2015), no. 4, 694–726. doi:10.1080/03605302.2014.979998
-
[31]
D. Dauvergne, J. Ortmann, and B. Vir´ ag. The directed landscape.Acta Mathematica, 229(2):201–285, 2022. doi:10.4310/ACTA.2022.v229.n2.a1
-
[32]
K. Deckelnick and C. M. Elliott,Uniqueness and error analysis for Hamilton–Jacobi equations with discontinuities, Interfaces Free Bound.6(2004), no. 3, 329–349. doi:10.4171/IFB/103
-
[33]
C. De Zan and P. Soravia,Cauchy problems for noncoercive Hamilton–Jacobi–Isaacs equations with discontinuous coef- ficients, Interfaces Free Bound.12(2010), no. 3, 347–368. doi:10.4171/IFB/238
-
[34]
W. E, K. Khanin, A. Mazel, and Y. G. Sinai. Invariant measures for Burgers equation with stochastic forcing.Annals of Mathematics, 151(3):877–960, 2000
2000
-
[35]
E. Emrah. Limit shapes for inhomogeneous corner growth models with exponential and geometric weights.Electronic Communications in Probability, 21:paper no. 42, 16 pp., 2016. doi:10.1214/16-ECP4
-
[36]
E. Emrah, C. Janjigian, and T. Sepp¨ al¨ ainen. Flats, spikes and crevices: the evolving shape of the inhomogeneous corner growth model.Electronic Journal of Probability, 26:paper no. 33, 45 pp., 2021. doi:10.1214/21-EJP595
-
[37]
E. Emrah, C. Janjigian, and T. Sepp¨ al¨ ainen. Anomalous geodesics in the inhomogeneous corner growth model.Commu- nications in Mathematical Physics, 406:paper no. 316, 2025. doi:10.1007/s00220-025-05470-8
-
[38]
L. C. Evans,Partial Differential Equations, 2nd ed., Graduate Studies in Mathematics, vol. 19, American Mathematical Society, Providence, RI, 2010
2010
-
[39]
P. A. Ferrari. TASEP hydrodynamics using microscopic characteristics.Probability Surveys, 15:1–27, 2018. doi:10.1214/17- PS284
work page doi:10.1214/17- 2018
-
[40]
N. Forcadel, C. Imbert, and R. Monneau. Germs for scalar conservation laws: the Hamilton–Jacobi equation point of view.arXiv preprint,arXiv:2407.04318, 2024
-
[41]
T. Franco and C. Landim, Hydrodynamic limit of gradient exclusion processes with conductances,Archive for Rational Mechanics and Analysis195(2010), 409–439. doi:10.1007/s00205-008-0206-5
-
[42]
N. Gantert, N. Georgiou, and D. Schmid. The TASEP on Galton–Watson trees.Electronic Journal of Probability, 26:paper no. 159, 38 pp., 2021. doi:10.1214/21-EJP725
-
[43]
M. Garavello, R. Natalini, B. Piccoli, and A. Terracina,Conservation laws with discontinuous flux, Netw. Heterog. Media 2(2007), no. 1, 159–179. doi:10.3934/nhm.2007.2.159
-
[44]
Georgiou, R
N. Georgiou, R. Kumar, and T. Sepp¨ al¨ ainen. TASEP with discontinuous jump rates.ALEA. Latin American Journal of Probability and Mathematical Statistics, 7:293–318, 2010
2010
-
[45]
Y. Giga, P. G´ orka, and P. Rybka,A comparison principle for Hamilton–Jacobi equations with discontinuous Hamiltonians, Proc. Amer. Math. Soc.139(2011), no. 5, 1777–1785. doi:10.1090/S0002-9939-2010-10630-5
-
[46]
Hairer, Solving the KPZ equation,Annals of Mathematics178(2013), no
M. Hairer, Solving the KPZ equation,Annals of Mathematics178(2013), no. 2, 559–664. doi:10.4007/annals.2013.178.2.4
-
[47]
T. E. Harris, Nearest-neighbor Markov interaction processes on multidimensional lattices,Advances in Mathematics9 (1972), no. 1, 66–89. doi: 10.1016/0001-8708(72)90030-8
-
[48]
D. Helbing, Traffic and related self-driven many-particle systems,Reviews of Modern Physics, 73(4):106, 2001. doi:10.1103/RevModPhys.73.1067
-
[49]
C. Imbert and R. Monneau. Flux-limited solutions for quasi-convex Hamilton–Jacobi equations on networks.Annales scientifiques de l’ ´Ecole Normale Sup´ erieure, 50(2):357–448, 2017. doi:10.24033/asens.2323
-
[50]
C. Imbert, R. Monneau, and H. Zidani,A Hamilton–Jacobi approach to junction problems and application to traffic flows, ESAIM Control Optim. Calc. Var.19(2013), no. 1, 129–166. doi:10.1051/cocv/2012002
-
[51]
Ishii,Hamilton–Jacobi equations with discontinuous Hamiltonians on arbitrary open sets, Bull
H. Ishii,Hamilton–Jacobi equations with discontinuous Hamiltonians on arbitrary open sets, Bull. Facul. Sci. Eng. Chuo University.28(1985), 33–77
1985
-
[52]
K. Johansson. Shape fluctuations and random matrices.Communications in Mathematical Physics, 209(2):437–476, 2000. doi:10.1007/s002200050027
-
[53]
Johansson, Non-intersecting paths, random tilings and random matrices
K. Johansson, Non-intersecting paths, random tilings and random matrices. Probability theory and related fields (2002), 123(2):225–280. doi:10.1007/s004400100187
-
[54]
M. Kardar, G. Parisi, and Y.-C. Zhang, Dynamic scaling of growing interfaces, Physical Review Letters, 56(9):889, 1986. 10.1103/PhysRevLett.56.889
-
[55]
T. Kriecherbauer, and J. Krug. A pedestrian’s view on interacting particle systems, KPZ universality and random matri- ces. Journal of Physics A: Mathematical and Theoretical,(2010): 43(40):403001. doi:10.1088/1751-8113/43/40/403001
-
[56]
S. N. Kruˇ zkov, (1970). First order quasilinear equations in several independent variables. Mathematics of the USSR- Sbornik, 10(2):217
1970
-
[57]
T. M. Liggett, An improved subadditive ergodic theorem,The Annals of Probability13(1985), no. 4, 1279–1285. doi: 10.1214/aop/1176992811. 62 JACOB BUTT, NICOS GEORGIOU, AND ENRICO SCALAS
-
[58]
K. Matetski, J. Quastel, and D. Remenik. The KPZ fixed point.Acta Mathematica, 227(1):115–203, 2021. doi:10.4310/ACTA.2021.v227.n1.a3
-
[59]
K. Matetski and D. Remenik. Exact solution of TASEP and variants with inhomogeneous speeds and memory lengths. Electronic Journal of Probability, 30:paper no. 167, 54 pp., 2025. doi:10.1214/25-EJP1430
-
[60]
F. Rezakhanlou. Continuum limit for some growth models.Stochastic Processes and their Applications, 101(1):1–41, 2002. doi:10.1016/S0304-4149(02)00100-X
-
[61]
L. T. Rolla and A. Q. Teixeira. Last passage percolation in macroscopically inhomogeneous media.Electronic Commu- nications in Probability, 13:131–139, 2008. doi:10.1214/ECP.v13-1287
-
[62]
H. Rost. Non-equilibrium behaviour of a many particle process: density profile and local equilibria.Zeitschrift f¨ ur Wahrscheinlichkeitstheorie und Verwandte Gebiete, 58:41–53, 1981. doi:10.1007/BF00536194
-
[63]
T. Sepp¨ al¨ ainen. A microscopic model for the Burgers equation and longest increasing subsequences.Electronic Journal of Probability, 1(5):1–51, 1996. doi:10.1214/EJP.v1-5
-
[64]
T. Sepp¨ al¨ ainen, Existence of hydrodynamics for the totally asymmetric simpleK-exclusion process,The Annals of Probability27(1999), no. 1, 361–415. doi: 10.1214/aop/1022677266
-
[65]
T. Sepp¨ al¨ ainen, Hydrodynamic profiles for the totally asymmetric exclusion process with a slow bond,Journal of Statis- tical Physics102(2001), no. 1–2, 69–96. doi: 10.1023/A:1026508625058
-
[66]
T. Sepp¨ al¨ ainen. Second class particles as microscopic characteristics in totally asymmetric nearest-neighbor-exclusion processes. Transactions of the American Mathematical Society, 353(12):4801–4829, 2001. 10.1090/S0002-9947-01-02872- 0
-
[67]
Sepp¨ al¨ ainen, The corner growth model with exponential weights, inRandom Growth Models, Proc
T. Sepp¨ al¨ ainen, The corner growth model with exponential weights, inRandom Growth Models, Proc. Sympos. Appl. Math., vol. 75, Amer. Math. Soc., Providence, RI, 2018, pp. 133–201
2018
-
[68]
T. Sepp¨ al¨ ainen and J. Krug. Hydrodynamics and platoon formation for a totally asymmetric exclusion model with particlewise disorder.Journal of Statistical Physics, 95:525–567, 1999. doi:10.1023/A:1007535124155
-
[69]
C. A. Tracy and H. Widom. Integral formulas for the asymmetric simple exclusion process.Communications in Mathe- matical Physics, 279:815–844, 2008. doi:10.1007/s00220-008-0443-3. Jacob Butt, University of Rome ‘La Sapienza’, Rome, Italy Email address:jacob.butt@uniroma1.it Nicos Georgiou, University of Sussex, Brighton, UK. Email address:N.Georgiou@susse...
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.