REVIEW 2 major objections 4 minor 64 references
Geometry dependence in linear interface growth
T0 review · 2 major / 4 minor · reviewed 2026-08-14 · deepseek-v4-flash
Pith's one-line read This paper establishes that the linear Edwards-Wilkinson and Mullins-Herring universality classes split into flat and radial subclasses, giving universal height-fluctuation variances and covariances that depend on geometry.
desk verdict A solid, useful paper: the new exact 1D radial MH variance is real and the geometry splitting for EW/MH holds up well, but the 2D radial values are torus values until a spherical-mode computation confirms them. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing object is the rescaled height fluctuation $\chi$ defined by the KPZ ansatz $h\simeq v_\infty t+(\Theta t)^\beta \chi$, whose variance $\langle\chi^2\rangle_c$ is the single non-null cumulant of the Gaussian height distribution. For the radial case in 1D the argument runs through the linearized growth equation in polar coordinates, Eq. (A1), with noise amplitude $2D/r$; Fourier decomposition of the radial fluctuations yields the closed formula $\langle\chi^2\rangle_c = c_z(z-1)^{1/z}I_z/(2^{1/z}\pi)$, with $z=2$ for EW and $z=4$ for MH, giving the new radial MH value. The flat-geometry variances come from the standard exact solutions of the EW and MH equations on fixed substrates. Numerically, the radial geometry is realized by stochastically duplicating columns so the substrate size grows as $L(t)=L_0+\omega t$, and universality is checked by comparing discrete models with direct integrations of the growth equations.
What would settle it
Directly integrate the 2D radial EW and MH equations in coordinates adapted to a growing circular or spherical domain and compare the extrapolated $\langle\chi^2\rangle_c$ values with the expanding-substrate numbers in Table III; a mismatch would show the column-duplication mapping fails in 2D.
Extended reading notes
Core claim
The central claim is that the EW and MH universality classes split into flat and radial subclasses, with universal, geometry-dependent one-point fluctuations and two-point correlations. Using the KPZ ansatz $h\simeq v_\infty t+(\Theta t)^\beta \chi$, the authors show that $P(\chi)$ is Gaussian for every case studied, with $\langle\chi\rangle=0$ and all higher cumulants zero; the only distinguishing cumulant is $\langle\chi^2\rangle_c$. For 1D, the exact values are $\langle\chi^2\rangle_c\approx 0.79788$ flat and $\approx 1.25331$ radial for EW, and $\approx 0.64697$ flat and $\approx 0.94573$ radial for MH, the last derived here from the radial linear equation. In 2D, where analytic mode sums are not available, extrapolations of integrations and of discrete models give EW variances $\approx 0.076$ flat versus $\approx 0.160$ radial, and MH variances $\approx 0.199$ flat versus $\approx 0.415$ radial. The rescaled spatial covariance $F(r/\xi)$ is universal for each geometry but different between them; in particular the oscillatory decay of flat MH becomes a monotonic decay in the radial case, and the 1D flat EW curve agrees with the Airy1 curve of flat 1D KPZ interfaces. The temporal covariances match the known analytical expressions in 1D and follow the predicted asymptotic exponents in 2D.
Load-bearing premise
The load-bearing premise is that growth on a substrate whose size increases stochastically and linearly in time, implemented by column duplication, faithfully reproduces true radial growth in two dimensions, since the equivalence is verified against exact results only in 1D.
Editorial extensions
If this is right
- The variance $\langle\chi^2\rangle_c$ can serve as a class indicator for Gaussian interfaces, since EW and MH values differ substantially in the same geometry and dimension.
- Geometry must be reported and controlled when height distributions or covariances are used to identify universality classes in experiments such as thin-film growth or electrodeposition.
- The 1D flat EW spatial covariance coinciding with the Airy1 curve suggests that flat 1D EW interfaces share a deeper statistical description with flat 1D KPZ interfaces.
- In 2D radial EW the temporal covariance does not follow a power-law decay because $\beta=0$, so its logarithmically corrected decay must be accounted for in scaling analyses.
- The column-duplication method for expanding substrates, previously used for nonlinear classes, also produces the radial statistics of linear classes.
Reading between the lines
- I would expect the same flat-versus-radial variance split to appear for the MH class in three dimensions, where $\beta=1/8$ and the growth regime is non-logarithmic; the paper does not test this.
- The closed formula for $\langle\chi^2\rangle_c$ in Appendix A is stated for general even $z$, so a direct check for $z=6$ would test whether the identity extends beyond EW and MH.
- The observed 1D flat EW covariance matching Airy1 hints that exact covariance identities for linear interfaces may be derivable from the same machinery as the KPZ class; the paper only reports the numerical coincidence.
- Since the paper's radial variances come from an expanding-substrate mapping rather than from genuinely curved coordinates in 2D, a direct integration of the linear equations on an expanding sphere would be a decisive independent check.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript studies whether the linear Edwards-Wilkinson (EW) and Mullins-Herring (MH) interface growth universality classes split into distinct flat and radial subclasses, analogous to the known splitting in nonlinear KPZ and VLDS classes. The authors derive exact one-dimensional variances for radial EW and, newly for radial MH (Eq. A12), and confirm these against numerical integrations of the growth equations and simulations of several discrete models. For 2D systems, they report variances extracted from integrations and discrete models on substrates whose size grows linearly in time (Tables II and III), and they analyze spatial and temporal covariances in both geometries and dimensions, finding universal but geometry-dependent scaling functions. The central claim is that the splitting is general to linear classes, including at the 2D EW upper critical dimension.
Significance. If the 2D results are accepted, the paper establishes a new general phenomenon: geometry-induced subclass splitting exists in linear as well as nonlinear universality classes. The 1D exact results, particularly the new radial MH variance, are clean and provide a useful anchor for the numerical method. The identification of universal, geometry-dependent amplitudes is also of practical interest for distinguishing universality classes in thin-film growth experiments. The paper is generally well organized, and the 1D analysis is convincing; the main risk is the unverified identification of the expanding square substrate with the true 2D radial (spherical) geometry, on which the central 2D claim depends.
major comments (2)
- [II, III.B, Table III] The identification of the expanding square lattice (a flat torus) with the 2D radial (spherical) geometry is not validated. In 1D, a circle and a periodic line have the same mode spectrum, so the exact circular EW and MH variances from Appendix A anchor the expanding-substrate data. In 2D, the sphere and the expanding flat torus have different Laplace-Beltrami spectra (spherical harmonics with degeneracy 2l+1 versus plane waves on a torus), and no exact or independent numerical result for the 2D spherical EW/MH variance is provided. The consistency across two parameter sets and two discrete models in Table III shows only that the duplication algorithm is robust, not that it reproduces spherical radial geometry. Because the central claim that 2D EW and MH classes split into flat and radial subclasses rests on the values <chi^2> = 0.160(1) (EW radial) and 0.415(3) (MH radial) in Table III and on the radial covariance curves in Figs. 3(c-d) and 5(c-d), the authors should either compute the variance from the spherical equations using the approach of Ref. [52], or provide an independent numerical test of the expanding-substrate equivalence in 2D, or explicitly reframe the claim as applying to expanding flat substrates and justify why this is the relevant radial subclass.
- [V, Fig. 5(c)] The 2D radial EW temporal covariance does not collapse onto a single scaling curve, and the authors introduce a three-parameter fit A(y) = a + b/(c + ln y) to describe the average data. This is a stated failure for one of the paper's main two-point quantities, and it weakens the universality claim for the temporal covariances in 2D radial geometry. The authors should either provide a theoretical argument for the logarithmic decay form, or explicitly present this case as an exception where universality is not yet established, rather than relying on an ad hoc fit.
minor comments (4)
- [Abstract] The phrase 'university classes' in the abstract should be 'universality classes'.
- [Appendix A, Eq. A9] The Poisson summation is truncated to the k=0 term with the statement that the main contribution comes from k=0. Since the z=4 result is new and used for calibration, it would be helpful to justify this truncation by estimating the k != 0 terms in the large-time limit, or to state explicitly that the result is asymptotic rather than exact.
- [III.B, Eq. 7] For the 2D EW class, the logarithmic ansatz in Eq. (7) is introduced separately from the general KPZ ansatz in Eq. (1). The authors should clarify the relationship between the two, since a reader may wonder whether the logarithmic form follows by a limiting procedure from Eq. (1) when beta = 0.
- [IV, Fig. 3] The fitted values of the exponent gamma for the 2D radial cases (gamma = 0.39 for EW and gamma = 0.61 for MH) are reported without error bars. Including uncertainties would make the claimed universality of the asymptotic decay more persuasive.
Circularity Check
No significant circularity: the central variances are exact solutions or cross-validated numerical integrations, and the expanding-substrate method is independently checked in 1D.
full rationale
The paper's central 1D results are exact solutions of the linear EW and MH equations: the flat variances come from standard results, the radial EW variance from Singha's solution, and the radial MH variance is derived in Appendix A from the polar-coordinate equation. These are parameter-free mathematical consequences, not fitted inputs. In 2D, the variances are obtained by direct numerical integration of the same linear equations on fixed and expanding substrates, with two independent parameter sets agreeing with each other; the discrete-model g2 values provide an independent cross-check because they compare separately simulated roughness amplitudes against the equation-derived variances. The expanding-substrate protocol is motivated by prior works including the authors' own, but it is validated in 1D against exact circular solutions, and in 2D it is an unverified geometric extrapolation (expanding square versus true spherical geometry) that constitutes a correctness risk, not a circular reduction. No fitted parameter is renamed as a prediction, and no load-bearing claim reduces to a self-citation or to a definitional identity. Therefore no circular step is identified.
Assumptions & free parameters
free parameters (10)
- nu_z (Family model, EW) =
0.78(2)
- nu_z (LC1 model, MH) =
0.61(2)
- nu_z (LC2 model, MH) =
0.142(6)
- 2D EW <chi^2>_c (flat) =
0.077(2) (EWI), 0.075(3) (EWII)
- 2D EW <chi^2>_c (radial) =
0.160(1) (EWI), 0.161(2) (EWII)
- 2D MH <chi^2>_c (flat) =
0.200(2) (MHI), 0.198(2) (MHII)
- 2D MH <chi^2>_c (radial) =
0.415(3) (MHI), 0.416(3) (MHII)
- gamma (spatial covariance exponent, 2D EW) =
approx 0.67 (flat), approx 0.39 (radial)
- gamma (spatial covariance exponent, 2D MH) =
approx 0 (flat), approx 0.61 (radial)
- a, b, c (temporal covariance fit, 2D radial EW) =
not reported
assumptions (6)
- domain assumption Growth of the discrete models is described in the continuum by the linear EW or MH equation (Eqs. 3-4) with Gaussian white noise; the models belong to these universality classes.
- domain assumption An interface growing on a substrate whose size increases linearly in time, L(t) = L0 + omega t, realizes the radial (curved) geometry for linear growth equations.
- domain assumption For radial geometry the linear equation can be linearized by replacing r^z with <r>^z = (r0 + F t)^z, and the noise correlator scaled by 1/(r0 + F t).
- standard math The k=0 term dominates the Poisson summation in Eq. A9 (S_z approx c_z sigma).
- domain assumption The KPZ-type ansatz (Eq. 1), h approx v_inf t + (Theta t)^beta chi, applies to linear classes, with Gaussian chi and zero mean.
- domain assumption For 2D EW (beta=0), the width grows as (D/nu_2)<chi^2> ln(nu_2 t/a), so the rescaled ansatz becomes h approx v_inf t + sqrt(Theta ln t) chi (Eq. 7).
Cite this review
Pith. "Pith review of Geometry dependence in linear interface growth." pith.science (2026). https://pith.science/paper/QT6APIT2
@misc{pith2026190810892,
author = {Pith},
title = {Pith review of: Geometry dependence in linear interface growth},
year = {2026},
howpublished = {\url{https://pith.science/paper/QT6APIT2}},
note = {Machine review of arXiv:1908.10892}
}
abstract
The effect of geometry in the statistics of \textit{nonlinear} universality classes for interface growth has been widely investigated in recent years and it is well known to yield a split of them into subclasses. In this work, we investigate this for the \textit{linear} classes of Edwards-Wilkinson (EW) and of Mullins-Herring (MH) in one- and two-dimensions. From comparison of analytical results with extensive numerical simulations of several discrete models belonging to these classes, as well as numerical integrations of the growth equations on substrates of fixed size (flat geometry) or expanding linearly in time (radial geometry), we verify that the height distributions (HDs), the spatial and the temporal covariances are universal, but geometry-dependent. In fact, the HDs are always Gaussian and, when defined in terms of the so-called "KPZ ansatz" $[h \simeq v_{\infty} t + (\Gamma t)^{\beta} \chi]$, their probability density functions $P(\chi)$ have mean null, so that all their cumulants are null, except by their variances, which assume different values in the flat and radial cases. The shape of the (rescaled) covariance curves is analyzed in detail and compared with some existing analytical results for them. Overall, these results demonstrate that the splitting of such university classes is quite general, being not restricted to the nonlinear ones.
Figures
Reference graph
Works this paper leans on
-
[52]
A. A. Masoudi, M. Khorrami, M. Stastna, and M. Ko- handel, Eur. Lett. 100, 16004 (2012)
work page 2012
-
[1]
Ohring, Materials Science of Thin Films-Deposition and Structure, 2nd ed
M. Ohring, Materials Science of Thin Films-Deposition and Structure, 2nd ed. (Academic Press, Waltham, Eng- land, 2001)
work page 2001
-
[2]
J. Krug and T. Michely, Islands, mounds and atoms: (Cambridge University Press, Cambridge, England, 1995)
work page 1995
-
[3]
4) equations are linear, the HDs in these classes are Gaussian, in both 1D and 2D
and MH (Eq. 4) equations are linear, the HDs in these classes are Gaussian, in both 1D and 2D. Thereby, the single non-null cumulant of P (χ) is its variance ⟨ χ2⟩ c, once as noted above ⟨χ⟩ = 0, and ⟨χn⟩c = 0 for n >2 for Gaussian distributions. So, any difference between HDs for flat and radial geometries might appear in ⟨ χ2⟩ c. This contrasts with the n...
-
[4]
K. J. M˚ aløy and J. Schmittbuhl, Phys. Rev. Lett. 87, 105502 (2001)
work page 2001
-
[5]
M. A. C. Huergo, M. A. Pasquale, A. E. Bolz´ an, A. J. Arvia, and P. H. Gonz´ alez, Phys. Rev. E - Stat. Nonlin- ear, Soft Matter Phys. 82, 1 (2010)
work page 2010
- [6]
- [7]
Show all 64 references
-
[8]
Villain, J
J. Villain, J. Phys. I 1, 19 (1991)
1991
-
[9]
Kardar, G
M. Kardar, G. Parisi, and Y.-C. Zhang, Phys. Rev. Lett. 56, 889 (1986)
1986
-
[10]
S. F. Edwards and D. R. Wilkinson, Proc. R. Soc. Lon- don, Ser. A 381, 17 (1982)
1982
-
[11]
Lai and S
Z.-W. Lai and S. Das Sarma, Phys. Rev. Lett. 66, 2348 (1991)
1991
-
[12]
Herring, Journal of Applied Physics 21, 301 (1950)
C. Herring, Journal of Applied Physics 21, 301 (1950)
1950
-
[13]
W. W. Mullins, Journal of Applied Physics 28, 333 (1957)
1957
-
[14]
Krug, Adv
J. Krug, Adv. Phys. 46, 139 (1997)
1997
-
[15]
Barabasi and H
A.-L. Barabasi and H. E. Stanley, Fractal Concepts in Surface Growth (Cambridge University Press, Cam- bridge, England, 1995)
1995
-
[16]
Pr¨ ahofer and H
M. Pr¨ ahofer and H. Spohn, Phys. Rev. Lett. 84, 4882 (2000)
2000
-
[17]
J. Krug, P. Meakin, and T. Halpin-Healy, Phys. Rev. A 45, 638 (1992)
1992
-
[18]
G. Amir, I. Corwin, and J. Quastel, Commun. Pure Appl. Math. 64, 466 (2011)
2011
-
[19]
Sasamoto and H
T. Sasamoto and H. Spohn, Phys. Rev. Lett. 104, 230602 (2010)
2010
-
[20]
S. G. Alves, T. J. Oliveira, and S. C. Ferreira, Eur. Lett. 96, 48003 (2011)
2011
-
[21]
Calabrese and P
P. Calabrese and P. Le Doussal, Phys. Rev. Lett. 106, 250603 (2011)
2011
-
[22]
S. G. Alves, T. J. Oliveira, and S. C. Ferreira, J. Stat. Mech. 2013, P05007 (2013)
2013
-
[23]
T. J. Oliveira, S. C. Ferreira, and S. G. Alves, Phys. Rev. E 85, 010601(R) (2012)
2012
-
[24]
S. N. Santalla, J. Rodr´ ıguez-Laguna, A. Celi, and R. Cuerno, J. Stat. Mech. 2017, P023201 (2017)
2017
-
[25]
Halpin-Healy and Y
T. Halpin-Healy and Y. Lin, Phys. Rev. E 89, 010103(R) (2014)
2014
-
[26]
K. A. Takeuchi, M. Sano, T. Sasamoto, and H. Spohn, Sci. Rep. 1, 34 (2011)
2011
-
[27]
K. A. Takeuchi and M. Sano, Phys. Rev. Lett. 104, 230601 (2010)
2010
-
[28]
I. S. S. Carrasco and T. J. Oliveira, Phys. Rev. E 98, 010102(R) (2018)
2018
-
[29]
Y. T. Fukai and K. A. Takeuchi, Phys. Rev. Lett. 119, 030602 (2017)
2017
-
[30]
Halpin-Healy, Phys
T. Halpin-Healy, Phys. Rev. Lett. 109, 170602 (2012)
2012
-
[31]
I. S. S. Carrasco and T. J. Oliveira, Phys. Rev. E 99, 032140 (2019)
2019
-
[32]
Halpin-Healy, Phys
T. Halpin-Healy, Phys. Rev. E 88, 042118 (2013)
2013
-
[33]
T. J. Oliveira, S. G. Alves, and S. C. Ferreira, Phys. Rev. E 87, 040102(R) (2013)
2013
-
[34]
I. S. S. Carrasco and T. J. Oliveira, Phys. Rev. E 94, 050801(R) (2016)
2016
-
[35]
I. S. S. Carrasco, K. A. Takeuchi, S. C. Ferreira, and T. J. Oliveira, New J. Phys. 14, 123057 (2014)
2014
-
[36]
A. A. Masoudi, S. Hosseinabadi, J. Davoudi, M. Khor- rami, and M. Kohandel, J. Stat. Mech. 2012, L02001 (2012)
2012
-
[37]
S. B. Singha, J. Stat. Mech. Theory Exp. 2005, P08006 (2005)
2005
-
[38]
J. Krug, H. Kallabis, S. N. Majumdar, S. J. Cornell, A. J. Bray, and C. Sire, Phys, Rev. E 56, 2702 (1997)
1997
-
[39]
Majaniemi, T
S. Majaniemi, T. Ala-Nissila, and J. Krug, Phys. Rev. B 53, 8071 (1996)
1996
-
[40]
Meakin, P
P. Meakin, P. Ramanlal, L. M. Sander, and R. C. Ball, Phys. Rev. A 34, 5091 (1986)
1986
-
[41]
Moser, J
K. Moser, J. Kert´ esz, and D. E. Wolf, Physica A 178, 215 (1991)
1991
-
[42]
J. M. Kim and S. Das Sarma, Phys. Rev. Lett. 72, 2903 (1994)
1994
-
[43]
Familyt and F
F. Familyt and F. Family, J. Phys. A Math. Gen 19, 441 (1986)
1986
-
[44]
Halpin-Healy and K
T. Halpin-Healy and K. A. Takeuchi, J. Stat. Phys. 160, 794 (2015)
2015
-
[45]
Krug, Phys
J. Krug, Phys. Rev. Lett. 72, 2907 (1994)
1994
-
[46]
Escudero, Phys
C. Escudero, Phys. Rev. E 74, 021901 (2006)
2006
-
[47]
KPZ ansatz
and, so, Θ = D2/ν2 = 1. Moreover, by construc- tion, for the Family and LC models one has D = 1/2 [14], whereas the parameters νz, with z = 2, 4, are not exactly known for these models and the scarce numeri- cal estimates of them are not so accurate. For instance, the value ν4...
-
[48]
Escudero, Phys
C. Escudero, Phys. Rev. E 73, 020902(R) (2006)
2006
-
[49]
Plischke, Z
M. Plischke, Z. R´ acz, and D. Liu, Phys. Rev. B35, 3485 (1987)
1987
-
[50]
Giacometti and M
A. Giacometti and M. Rossi, Phys. Rev. E 63, 046102 (2001)
2001
-
[51]
Lam and L
C.-H. Lam and L. M. Sander, Phys. Rev. Lett. 71, 561 (1993)
1993
-
[53]
Escudero, J
C. Escudero, J. Stat. Mech. 2009, P07020 (2009)
2009
-
[54]
Escudero, Ann
C. Escudero, Ann. Phys. (N. Y). 324, 1796 (2009)
2009
-
[55]
Pr¨ ahofer and H
M. Pr¨ ahofer and H. Spohn, J. Stat. Phys. 108, 1071 (2002)
2002
-
[56]
Sasamoto, J
T. Sasamoto, J. Phys. A 38, L549 (2005)
2005
-
[57]
Borodin, P
A. Borodin, P. Ferrari, and T. Sasamoto, Commun. Math. Phys. 283, 417 (2008). 11
2008
-
[58]
Kallabis and J
H. Kallabis and J. Krug, Eur. Lett. 45, 20 (1999)
1999
-
[59]
R. A. L. Almeida, S. O. Ferreira, T. J. Oliveira, and F. D. A. Aar˜ ao Reis, Phys. Rev. B89, 045309 (2014)
2014
-
[60]
Halpin-Healy and G
T. Halpin-Healy and G. Palasantzas, Europhys. Lett. 105, 50001 (2014)
2014
-
[61]
R. A. L. Almeida, S. O. Ferreira, I. R. B. Ribeiro, and T. J. Oliveira, Eur. Lett. 109, 46003 (2015)
2015
-
[62]
R. A. L. Almeida, S. O. Ferreira, I. Ferraz, and T. J. Oliveira, Sci. Rep. 7, 3773 (2017)
2017
-
[63]
I. S. Brandt, V. C. Zoldan, V. Stenger, C. C. Pl´ a Cid, A. A. Pasa, T. J. Oliveira, and F. D. A. Aar˜ ao Reis, J. Appl. Phys. 118, 145303 (2015)
2015
-
[64]
P. A. Orrillo, S. N. Santalla, R. Cuerno, L. V´ azquez, S. B. Ribotta, L. M. Gassa, F. J. Mompean, R. C. Salvarezza, and M. E. Vela, Sci. Rep. 7, 17997 (2017)
2017
Reviewed August 14, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.