Exponents and front fluctuations in the quenched Kardar-Parisi-Zhang universality class of one- and two- dimensional interfaces
Pith reviewed 2026-05-16 21:49 UTC · model grok-4.3
The pith
Quenched KPZ interfaces in one and two dimensions show scaling exponents matching the directed percolation depinning class.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Simulations of an automaton discretization of the quenched KPZ equation in one and two dimensions at the depinning transition directly yield the exponents α, β, θ, δ, and z. These values are compatible with the directed percolation depinning class. The probability density function of the front fluctuations is computed in the growth regime and found to be asymptotically non-Gaussian with skewness and kurtosis values that also differ from the KPZ equation with time-dependent noise in both the central region and the tails.
What carries the argument
The automaton discretization of the quenched KPZ equation, together with direct measurement of the height-difference correlation function to obtain the dynamic exponent z.
If this is right
- Interface roughening and depinning in quenched media obey the same scaling relations as directed percolation in one and two dimensions.
- The probability density function of front fluctuations supplies an independent signature that distinguishes quenched noise from time-dependent noise.
- Universality allows direct transfer of known directed percolation results to predict correlation lengths and fluctuation statistics in physical interfaces.
- The measured dynamic exponent z obtained from real-space height correlations closes the set of scaling relations for the universality class.
Where Pith is reading between the lines
- Real experiments on fluid invasion or magnetic domain walls in two dimensions could test the predicted non-Gaussian front statistics without needing to resolve microscopic disorder details.
- If the compatibility persists when the same automaton is run in three dimensions, the directed percolation mapping would extend beyond the physical substrate dimensions studied here.
- The difference between the quenched and time-dependent PDFs suggests that annealed versus quenched noise may belong to distinct universality classes even when the deterministic part of the equation is identical.
Load-bearing premise
The automaton discretization faithfully reproduces the continuum quenched KPZ dynamics and the simulated system sizes and times are large enough to reach the asymptotic scaling regime without crossover effects.
What would settle it
A larger-scale simulation or laboratory experiment on a quenched KPZ interface in one or two dimensions that returns an exponent set, for example the roughness exponent α, clearly outside the range reported for the directed percolation depinning class.
Figures
read the original abstract
We have simulated an automaton version of the quenched Kardar-Parisi-Zhang (qKPZ) equation in one and two dimensions in order to study the scaling properties of the interface at the depinning transition. Specifically, the $\alpha$, $\beta$, $\theta$, and $\delta$ critical exponents characterizing the surface kinetic roughening and depinning behaviors have been directly computed from the simulations. In addition, by studying the height-difference correlation function in real space, we have also been able to directly compute the dynamic correlation length and its associated dynamic critical exponent $z$. The full sets of scaling exponents are largely compatible with those of the Directed Percolation Depinning universality class for one and two dimensional interfaces. Furthermore, we have computed numerically the probability density function (PDF) of the front fluctuations in the growth regime, finding its asymptotic form in one and two dimensions. While the PDF features strongly non-Gaussian skewness and kurtosis values, it also differs from the PDF of the KPZ equation with time-dependent noise for physical substrate dimensions, both in the central part and at the tails of the distribution.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript reports simulations of an automaton discretization of the quenched Kardar-Parisi-Zhang (qKPZ) equation in one and two dimensions at the depinning transition. The authors directly measure the roughness exponent α, growth exponent β, and the additional exponents θ, δ, and z from the height-difference correlation function, concluding that the full set is largely compatible with the Directed Percolation Depinning (DPD) universality class. They further compute the probability density function of front fluctuations in the growth regime and characterize its non-Gaussian features, which differ from those of the KPZ equation with time-dependent noise.
Significance. If the simulations are shown to reach the true asymptotic regime, the work would strengthen the numerical case for qKPZ depinning belonging to the DPD class in low dimensions and supply a concrete characterization of fluctuation statistics. The direct extraction of z from real-space correlations and the PDF analysis are useful additions to the literature on interface depinning.
major comments (3)
- [Numerical methods and scaling analysis] The central claim of compatibility with the DPD class rests on the assumption that the automaton faithfully reproduces the continuum qKPZ dynamics and that the measured exponents are asymptotic. No explicit tests of exponent stability versus system size L or time are presented, nor is there a comparison against alternative discretizations or continuum integrations; this is especially critical in 2D where depinning crossovers are known to be slow.
- [Results on critical exponents] The statement that the exponents are 'largely compatible' with DPD values is not accompanied by quantitative error bars on the differences or by data-collapse quality metrics; without these, it is difficult to assess whether the agreement is within numerical uncertainty or merely pre-asymptotic.
- [Dynamic correlation length] For the dynamic exponent z extracted from the height-difference correlation function, the manuscript does not discuss possible corrections-to-scaling or the fitting window used; this leaves open the possibility that the reported z reflects a transient rather than the true asymptotic value.
minor comments (2)
- [Front fluctuation statistics] The PDF figures would benefit from explicit reporting of skewness and kurtosis values together with their uncertainties.
- [Introduction] Notation for the exponents (α, β, θ, δ, z) should be defined once in a dedicated table or section to avoid ambiguity when comparing to DPD literature.
Simulated Author's Rebuttal
We thank the referee for the thoughtful and detailed report. We address each major comment below and have revised the manuscript to strengthen the presentation of our numerical results and scaling analysis.
read point-by-point responses
-
Referee: [Numerical methods and scaling analysis] The central claim of compatibility with the DPD class rests on the assumption that the automaton faithfully reproduces the continuum qKPZ dynamics and that the measured exponents are asymptotic. No explicit tests of exponent stability versus system size L or time are presented, nor is there a comparison against alternative discretizations or continuum integrations; this is especially critical in 2D where depinning crossovers are known to be slow.
Authors: We agree that explicit checks of asymptotic convergence are important, particularly in 2D. In the revised manuscript we add supplementary figures that plot the effective exponents versus L and versus time for the largest accessible scales; these show clear stabilization of α, β, θ, δ and z within the statistical uncertainties of our runs. We also expand the methods section to recall that the automaton is a standard, previously validated discretization of qKPZ whose exponents have been shown to match continuum results in the literature. While we have not performed new continuum integrations (owing to their substantially higher computational cost), the added convergence data directly address the concern about pre-asymptotic behavior. revision: partial
-
Referee: [Results on critical exponents] The statement that the exponents are 'largely compatible' with DPD values is not accompanied by quantitative error bars on the differences or by data-collapse quality metrics; without these, it is difficult to assess whether the agreement is within numerical uncertainty or merely pre-asymptotic.
Authors: We have revised the results section to report statistical error bars on every exponent, obtained from both ensemble fluctuations across independent realizations and from the covariance of the fitting procedures. We further include data-collapse plots for the height-difference correlation function together with a quantitative goodness-of-collapse metric (mean-squared deviation from the master curve). These additions allow the reader to judge that the observed agreement with DPD values lies within the estimated numerical uncertainties. revision: yes
-
Referee: [Dynamic correlation length] For the dynamic exponent z extracted from the height-difference correlation function, the manuscript does not discuss possible corrections-to-scaling or the fitting window used; this leaves open the possibility that the reported z reflects a transient rather than the true asymptotic value.
Authors: The revised text now explicitly states the fitting windows (in both real-space and time) used to extract z and shows the effective exponent versus scale to demonstrate the absence of strong drift. We also include a brief analysis of leading corrections-to-scaling by fitting an extended form that incorporates a sub-leading term; the leading exponent remains stable and consistent with the DPD value once the fitting window is restricted to the largest scales. revision: yes
Circularity Check
No circularity: all results are direct numerical measurements
full rationale
The paper reports direct numerical simulations of an automaton discretization of the qKPZ equation in 1D and 2D. Scaling exponents (α, β, θ, δ, z) and the PDF of front fluctuations are computed from the raw simulation trajectories without any analytic derivation, parameter fitting that is then re-labeled as a prediction, or load-bearing self-citation of a uniqueness theorem. The compatibility statement with the Directed Percolation Depinning class is a post-hoc numerical comparison, not a reduction of the reported values to the paper's own inputs. The derivation chain is therefore self-contained and non-circular.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption The discrete automaton version belongs to the quenched KPZ universality class
Reference graph
Works this paper leans on
-
[1]
and others) for this and other discrete models in the qKPZ universality class, and for the equation itself, see Ref. [31] for a review. Very recently, further numerical and analytical works [45, 46] have revised previous arguments and exponent values in the literature [42], concluding in particular that the dynamic exponentzdoes not derive directly from t...
- [2]
- [3]
-
[4]
G. Pruessner,Self-Organised Criticality: Theory, Models and Characterisation(Cambridge University Press, Cambridge, 13 UK, 2012)
work page 2012
-
[5]
U. C. T ¨auber,Critical dynamics: a field theory approach to equilibrium and non-equilibrium scaling behavior(Cambridge University Press, Cambridge, UK, 2014)
work page 2014
-
[6]
A.-L. Barab ´asi and H. Stanley,Fractal Concepts in surface growth(Cambrige University Press, Cambridge, UK, 1995)
work page 1995
- [7]
- [8]
-
[9]
T. Kriecherbauer and J. Krug, J. Phys. A: Math. Theor.43, 403001 (2010)
work page 2010
- [10]
-
[11]
K. A. Takeuchi, Physica A504, 77 (2018)
work page 2018
-
[12]
Q. Fontaine, D. Squizzato, F. Baboux, I. Amelio, A. Lema ˆıtre, M. Morassi, I. Sagnes, L. L. Gratiet, A. Harouri, M. Wouters, I. Carusotto, A. Amo, M. Richard, A. Minguzzi, L. Canet, S. Ravets, and J. Bloch, Nature608, 687 (2022)
work page 2022
-
[13]
D. Wei, A. Rubio-Abadal, B. Ye, F. Machado, J. Kemp, K. Srakaew, S. Hollerith, J. Rui, S. Gopalakrishnan, N. Y . Yao, I. Bloch, and J. Zeiher, Science376, 716 (2022)
work page 2022
-
[14]
L. M. Sieberer, M. Buchhold, J. Marino, and S. Diehl, Rev. Mod. Phys.97, 25004 (2025)
work page 2025
-
[15]
F. Caballero, C. Nardini, and M. E. Cates, J. Stat. Mech.2018, 123208 (2018)
work page 2018
- [16]
- [17]
- [18]
- [19]
- [20]
- [21]
- [22]
-
[23]
K. A. Takeuchi and M. Sano, Phys. Rev. Lett.104, 230601 (2010)
work page 2010
-
[24]
R. A. L. Almeida, S. O. Ferreira, T. J. Oliveira, and F. D. A. Aar˜ao-Reis, Phys. Rev. B89, 045309 (2014)
work page 2014
- [25]
-
[26]
P. A. Orrillo, S. N. Santalla, R. Cuerno, L. V ´azquez, S. B. Ri- botta, L. M. Gassa, F. J. Mompean, R. C. Salvarezza, and M. E. Vela, Sci. Rep.7, 17997 (2017)
work page 2017
-
[27]
I. S. S. Carrasco and T. J. Oliveira, Phys. Rev. E94, 050801(R) (2016)
work page 2016
-
[28]
E. Rodr ´ıguez-Fern´andez and R. Cuerno, Phys. Rev. Res.3, L012020 (2021)
work page 2021
-
[29]
J. M. Marcos, P. Rodr ´ıguez-L´opez, J. J. Mel ´endez, R. Cuerno, and J. J. Ruiz-Lorenzo, Phys. Rev. E105, 054801 (2022)
work page 2022
-
[30]
H. Vaquero del Pino and R. Cuerno, Phys. Rev. Res.7, 043192 (2025)
work page 2025
-
[31]
M. Allemand, G. Dupuy, P. Paquiez, N. Dupuis, A. Ranc ¸on, T. Roscilde, T. Chalopin, and D. Cl ´ement, (2025), arXiv:2508.21623 [cond-mat.quant-gas]
-
[32]
K. J. Wiese, Rep. Prog. Phys.85, 086502 (2022)
work page 2022
- [33]
- [34]
-
[35]
K.-W. Moon, D.-H. Kim, S.-C. Yoo, C.-G. Cho, S. Hwang, B. Kahng, B.-C. Min, K.-H. Shin, and S.-B. Choe, Phys. Rev. Lett.110, 107203 (2013)
work page 2013
-
[36]
S. Atis, A. K. Dubey, D. Salin, L. Talon, P. Le Doussal, and K. J. Wiese, Phys. Rev. Lett.114, 234502 (2015)
work page 2015
- [37]
- [38]
- [39]
-
[40]
J. J. Ramasco, J. M. L ´opez, and M. A. Rodr´ıguez, Phys. Rev. E 64, 066109 (2001)
work page 2001
-
[41]
L.-H. Tang, M. Kardar, and D. Dharg, Phys. Rev. Lett.74, 920 (1995)
work page 1995
- [42]
-
[43]
L. A. N. Amaral, A.-L. Barab ´asi, S. V . Buldyrev, S. T. Harring- ton, S. Havlin, R. Sadr-Lahijany, and H. E. Stanley, Phys. Rev. E64, 4655 (1995)
work page 1995
-
[44]
Y . J. Chen, S. Papanikolaou, J. P. Sethna, S. Zapperi, and G. Durin, Phys. Rev. E84, 061103 (2011)
work page 2011
-
[45]
Y . J. Chen, S. Zapperi, and J. P. Sethna, Phys. Rev. E92, 022146 (2015)
work page 2015
-
[46]
G. Mukerjee, J. A. Bonachela, M. A. Mu ˜noz, and K. J. Wiese, Phys. Rev. E107, 054136 (2023)
work page 2023
- [47]
-
[48]
H. S. Song and J. M. Kim, J. Korean Phys. Soc.51, 1630 (2007)
work page 2007
- [49]
- [50]
-
[51]
B. G. Barreales, J. J. Mel ´endez, R. Cuerno, and J. J. Ruiz- Lorenzo, Phys. Rev. E108, 044801 (2023)
work page 2023
-
[52]
E. E. Mozo Luis, T. A. de Assis, and F. A. Oliveira, J. Stat. Mech.2022, 083202 (2022)
work page 2022
-
[53]
B. G. Barreales, J. J. Mel ´endez, R. Cuerno, and J. J. Ruiz- Lorenzo, J. Stat. Mech. , 023203 (2020)
work page 2020
-
[54]
B. Efron,The Jackknife, the Bootstrap, and Other Resam- pling Plans(Society for Industrial and Applied Mathematics, Philadelphia, 1982)
work page 1982
-
[55]
P. Young,Everything you wanted to know about data analysis and fitting but were afraid to ask(Springer, London, 2015)
work page 2015
-
[56]
J. M. L ´opez, M. A. Rodr´ıguez, and R. Cuerno, Phys. Rev. E56, 3993 (1997)
work page 1997
-
[57]
J. J. Ramasco, J. M. L ´opez, and M. A. Rodr ´ıguez, Phys. Rev. Lett.84, 2199 (2000)
work page 2000
-
[58]
R. Cuerno and L. V ´azquez, inAdvances in Condensed Matter and Statistical Physics, edited by E. Korutcheva and R. Cuerno (Nova Science Publishers, New York, 2004)
work page 2004
-
[59]
E. Toivonen, M. Molkkari, E. R ¨as¨anen, and L. Laurson, Phys. Rev. Lett.129, 175701 (2022)
work page 2022
-
[60]
T. Sillanp ¨a¨a, S. Nousiainen, and L. Laurson, (2025), arXiv:2511.16443 [cond-mat.stat-mech]
-
[61]
I. S. S. Carrasco and T. J. Oliveira, Phys. Rev. E100, 042107 (2019)
work page 2019
-
[62]
A. Tajuelo, J. Trujillo-Mulero, J. J. Mel ´endez, R. Cuerno, and J. J. Ruiz-Lorenzo, 10.5281/zenodo.17879472 (2025)
- [63]
- [64]
-
[65]
T. J. Oliveira, S. G. Alves, and S. C. Ferreira, Phys. Rev. E87, 040102 (2013)
work page 2013
- [66]
-
[67]
S. N. Majumdar and G. Schehr, J. Stat. Mech.2014, P01012 (2014)
work page 2014
-
[68]
J. M. Kim, M. A. Moore, and A. J. Bray, Phys. Rev. A44, 2345 (1991)
work page 1991
- [69]
-
[70]
S. V . Buldyrev, A.-L. Barab ´asi, F. Caserta, S. Havlin, H. E. Stanley, and T. Vicsek, Phys. Rev. A45, R8313 (1992)
work page 1992
- [71]
-
[72]
S. V . Buldyrev, S. Havlin, J. Kert ´esz, A. Shehter, and H. E. Stanley, Fractals01, 827 (1993)
work page 1993
- [73]
- [74]
-
[75]
L. H. Tang and H. Leschhorn, Phys. Rev. A45, R8309 (1992)
work page 1992
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.