pith. sign in

arxiv: 2604.13821 · v1 · submitted 2026-04-15 · ❄️ cond-mat.mtrl-sci · nlin.CG

Step Bunching and Meandering as Common Growth Modes: A Discrete Model and a Continuum Description

Vassil Ivanov , Vesselin Tonchev , Marta A. Chabowska , Hristina Popova , Magdalena A. Za{\l}uska-Kotur This is my paper

Pith reviewed 2026-05-10 13:15 UTC · model grok-4.3

classification ❄️ cond-mat.mtrl-sci nlin.CG
keywords step bunchingstep meanderingVicCAcontinuum modelcrystal growthsurface patternsEhrlich-Schwoebel effectvicinal surfaces
0
0 comments X

The pith

A properly shaped potential energy landscape in the discrete model allows step bunching and meandering to coexist, matching patterns from the continuum description.

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

The paper confronts the apparent contradiction that step bunching and step meandering cannot occur together because they require opposite Ehrlich-Schwoebel effects. It does this by comparing results from a two-dimensional-plus-time Vicinal Cellular Automaton model and a differential-difference partial differential equation model that adds a relaxation term to a bunching description. When the potential energy landscape in the automaton model is given an appropriate shape, both approaches generate similar surface patterns. This establishes a parameter-level link between the discrete and continuum frameworks. Sympathetic readers would care as it offers a route to modeling mixed instability regimes observed in real crystal growth processes.

Core claim

By introducing a proper shape of the potential energy landscape in the VicCA model, similar patterns to those from the PDE-based continuum description are produced, thereby linking the two models on the level of parameters and showing that bunching and meandering can coexist as common growth modes.

What carries the argument

The potential energy landscape in the VicCA discrete model, tuned to match the relaxation term added to the bunching model in the continuum description.

If this is right

  • The continuous framework can explore long-time scales to reveal a large variety of surface patterns.
  • Both instabilities can be captured together without being treated as a simple superposition.
  • The models become connected through shared parameters once the potential shape is chosen appropriately.
  • Coexistence becomes possible in the same growth process despite the usual association with opposing effects.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • This suggests that similar tuning of energy landscapes could unify other pairs of apparently conflicting surface instabilities.
  • Experimentalists might look for growth conditions where both bunching and meandering appear together to test the linked parameters.
  • Extending the relaxation term to two dimensions or other geometries could reveal even richer pattern formation.

Load-bearing premise

That the opposing Ehrlich-Schwoebel effects do not prevent simultaneous bunching and meandering once the potential landscape shape is adjusted in the discrete model and a relaxation term is added in the continuum one.

What would settle it

If simulations with the tuned potential in VicCA show that one instability always dominates and prevents the other from appearing, or if no matching patterns are found between the two models, the unification claim would be falsified.

Figures

Figures reproduced from arXiv: 2604.13821 by Hristina Popova, Magdalena A. Za{\l}uska-Kotur, Marta A. Chabowska, Vassil Ivanov, Vesselin Tonchev.

Figure 2
Figure 2. Figure 2: FIG. 2: 3D surface reconstruction based on Eq. 26 [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 1
Figure 1. Figure 1: FIG. 1: First moments of an already destabilized vicinal [PITH_FULL_IMAGE:figures/full_fig_p005_1.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3: Structures obtained for [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4: a) Crystal structure with the diffusing adatoms [PITH_FULL_IMAGE:figures/full_fig_p007_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5: Structures obtained for [PITH_FULL_IMAGE:figures/full_fig_p008_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6: Comparison of surface morphologies obtained [PITH_FULL_IMAGE:figures/full_fig_p009_6.png] view at source ↗
read the original abstract

The coexistence of step bunching and step meandering remains contradictory in the understanding of the unstable step-flow growth. Considered separately, the two instabilities have generated rich but largely independent modeling traditions. Especially, the one-dimensional framework faces a fundamental difficulty once bunching and meandering occur simultaneously -- step bunching is usually associated with an inverted Ehrlich--Schwoebel effect, whereas step meandering is associated with a direct one. The key experiments also focus mainly on the two basic limiting cases. How, then, can both instabilities coexist within the same growth process once the simultaneous occurrence of bunching and meandering cannot be adequately captured as a simple superposition of the two? In this work, we confront results from two substantially different approaches: a (2+1)D Vicinal Cellular Automaton based model (VicCA) and a differential-difference PDE-based description combining a model of step bunching with a relaxation term in the perpendicular direction. The continuous framework enables to explore long-time scales evolution to find large variety of surface patterns. Introducing a proper shape of the potential energy landscape in the VicCA model produces similar patterns and links both models on the level of parameters.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript claims that the apparent contradiction in modeling simultaneous step bunching and meandering—arising from opposing Ehrlich-Schwoebel effects—can be resolved by augmenting a continuum bunching model with a perpendicular relaxation term in a differential-difference PDE framework, and that an appropriately shaped potential energy landscape in a (2+1)D Vicinal Cellular Automaton (VicCA) discrete model generates matching mixed patterns, thereby establishing a direct parameter-level correspondence between the two approaches.

Significance. If the pattern reproduction and parameter linkage are rigorously demonstrated, the work bridges independent discrete and continuum traditions for step-flow instabilities, enabling long-time-scale exploration of complex surface morphologies that neither limiting case captures alone. This could inform experimental interpretation of mixed bunching-meandering regimes and reduce reliance on ad-hoc superpositions of instabilities.

major comments (2)
  1. [Abstract] Abstract and introduction: the claim that the potential energy landscape shape 'links both models on the level of parameters' requires explicit mapping (e.g., how the landscape parameters correspond to the relaxation coefficient or effective ES barrier signs). Without this, the correspondence risks being qualitative pattern-matching rather than a quantitative parameter equivalence, which is load-bearing for the central resolution of the ES contradiction.
  2. [Introduction] The weakest assumption noted in the review—that opposing ES effects do not prevent simultaneous occurrence—is addressed by the non-standard potential, but the manuscript must show (via stability analysis or simulation sweeps) that the chosen landscape permits both instabilities without one suppressing the other, rather than assuming coexistence from pattern similarity alone.
minor comments (2)
  1. The abstract mentions 'long-time scales evolution' and 'large variety of surface patterns' in the continuum model; include representative time-series or phase diagrams in the results to illustrate the claimed variety.
  2. Notation for the potential energy landscape and the relaxation term should be unified or cross-referenced between the discrete and continuum sections to aid comparison.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the detailed and constructive report. We address the major comments point by point below, indicating the revisions we will incorporate to strengthen the manuscript.

read point-by-point responses

  1. Referee: [Abstract] Abstract and introduction: the claim that the potential energy landscape shape 'links both models on the level of parameters' requires explicit mapping (e.g., how the landscape parameters correspond to the relaxation coefficient or effective ES barrier signs). Without this, the correspondence risks being qualitative pattern-matching rather than a quantitative parameter equivalence, which is load-bearing for the central resolution of the ES contradiction.

    Authors: We agree that an explicit parameter mapping is necessary to elevate the correspondence beyond pattern similarity. In the revised manuscript we will add a dedicated subsection (and supporting table) that derives the direct mapping: the curvature and asymmetry parameters of the VicCA potential are shown to control the sign and magnitude of the effective Ehrlich-Schwoebel barrier, while the overall depth sets the strength of the perpendicular relaxation term in the differential-difference PDE. The mapping is obtained by matching the linear instability thresholds and the adatom current expressions between the two models, thereby providing the quantitative link requested. revision: yes

  2. Referee: [Introduction] The weakest assumption noted in the review—that opposing ES effects do not prevent simultaneous occurrence—is addressed by the non-standard potential, but the manuscript must show (via stability analysis or simulation sweeps) that the chosen landscape permits both instabilities without one suppressing the other, rather than assuming coexistence from pattern similarity alone.

    Authors: We accept the need for explicit verification. The revised version will include (i) a linear stability analysis of the continuum model with the augmented relaxation term that identifies the parameter window in which both bunching and meandering modes are simultaneously unstable, and (ii) a systematic sweep of the VicCA potential parameters (with corresponding continuum equivalents) that demonstrates the coexistence regime without suppression. These results will be presented as new figures and a short subsection, confirming that the chosen landscape indeed allows both instabilities to develop concurrently. revision: yes

Circularity Check

0 steps flagged

No significant circularity detected

full rationale

The paper compares two independent modeling frameworks—a (2+1)D VicCA discrete automaton and a differential-difference continuum PDE that augments bunching with perpendicular relaxation—and demonstrates that a suitably chosen potential-energy landscape in the discrete model can reproduce the mixed bunching+meandering patterns obtained from the continuum description. This establishes a parameter-level correspondence but does not constitute a derivation in which any claimed prediction or result is forced by construction to equal its own inputs. No load-bearing self-citation, self-definitional step, or renaming of a known result appears in the provided abstract or construction; the potential shape is introduced as an adjustable feature that permits coexistence of the two instabilities, not as a tautological fit. The central claim therefore remains self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Abstract-only review yields no identifiable free parameters, axioms, or invented entities; the potential energy landscape shape is referenced but not specified.

pith-pipeline@v0.9.0 · 5534 in / 1024 out tokens · 32923 ms · 2026-05-10T13:15:58.103725+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

40 extracted references · 40 canonical work pages

  1. [1]

    orthogonal

    to obtain the general form of the model of interest. Let {u(t, x)}N−1 n=0 be the positions of N steps at time t, where x is the coordinate along the step edge (perpendicular to the step flow direction). The steps are ordered such that un < u n+1 for all x at t = 0. The system Hamiltonian is then given by: H({u n}) = Z L 0 N−1X n=0 " ˜β 2 ∂un ∂x 2 +V(∆u n)...

    work page 2023

  2. [2]

    Research, Innovation and Digitali- sation for Smart Transformation

    and HPC resources of the GATE Institute, Bulgaria, via the programme “Research, Innovation and Digitali- sation for Smart Transformation” 2021-2027 (PRIDST) (grant agreement no. BG16RFPR002-1.014-0010-C01). V.I. would like to further thank to Boris Kraychev from GATE Institute, Bulgaria for the support with the compu- tational resources. V.T. and V.I. tha...

    work page 2021

  3. [3]

    Kangawa, A

    Y. Kangawa, A. Kusaba, T. Akiyama, S. Nitta, M. Iwaya, H. Miyake, P. Kempisty, M. A. Chabowska, H. Popova, and M. A. Za luska-Kotur, Chemical vapor deposition of compound semiconductors: process simulation and exper- imental validation, Applied Physics Express19, 010101 (2026)

    work page 2026

  4. [4]

    Misbah, O

    C. Misbah, O. Pierre-Louis, and Y. Saito, Crystal surfaces in and out of equilibrium: A modern view, Reviews of Modern Physics82, 981–1040 (2010)

    work page 2010

  5. [5]

    Ohtani, M

    N. Ohtani, M. Katsuno, T. Aigo, T. Fujimoto, H. Tsuge, H. Yashiro, and M. Kanaya, Step bunching behaviour on the {0001} surface of hexagonal SiC, Journal of crystal growth210, 613 (2000)

    work page 2000

  6. [6]

    N´ eel, T

    N. N´ eel, T. Maroutian, L. Douillard, and H.-J. Ernst, From meandering to faceting, is step flow growth ever stable?, Physical Review Letters91, 226103 (2003)

    work page 2003

  7. [7]

    H. Omi, Y. Homma, V. Tonchev, and A. Pimpinelli, New types of unstable step-flow growth on Si (111)-(7 × 7) during molecular beam epitaxy: Scaling and universality, Physical Review Letters95, 216101 (2005)

    work page 2005

  8. [8]

    Galiana, M

    B. Galiana, M. Benedicto, and P. Tejedor, Step-step in- teractions on gaas (110) nanopatterns, Journal of Applied Physics113, 10.1063/1.4774215 (2013)

    work page doi:10.1063/1.4774215 2013

  9. [9]

    L. Huo, R. Lingaparthi, N. Dharmarasu, and K. Radhakr- ishnan, Growth instability of N-polar GaN on vicinal SiC substrate using plasma-assisted molecular beam epitaxy, Thin Solid Films808, 140572 (2024)

    work page 2024

  10. [10]

    de Theije, J

    F. de Theije, J. Schermer, and W. van Enckevort, Effects of nitrogen impurities on the cvd growth of diamond: step bunching in theory and experiment, Diamond and Related Materials9, 1439–1449 (2000)

    work page 2000

  11. [11]

    Liu and J

    D.-J. Liu and J. D. Weeks, Quantitative theory of current- induced step bunching on Si(111), Physical Review B57, 14891–14900 (1998)

    work page 1998

  12. [12]

    Tersoff, Y

    J. Tersoff, Y. H. Phang, Z. Zhang, and M. G. Lagally, Step-bunching instability of vicinal surfaces under stress, Physical Review Letters75, 2730–2733 (1995)

    work page 1995

  13. [13]

    Krasteva, H

    A. Krasteva, H. Popova, N. Akutsu, and V. Tonchev, Time 11 scaling relations for step bunches from models with step- step attractions (B1-type models), inAIP Conference Proceedings, Vol. 1722 (AIP Publishing LLC, 2016) p. 220015

    work page 2016

  14. [14]

    Stoyanov and V

    S. Stoyanov and V. Tonchev, Properties and dynamic interaction of step density waves at a crystal surface during electromigration affected sublimation, Physical Review B58, 1590–1600 (1998)

    work page 1998

  15. [15]

    Stoyanov, New type of step bunching instability at vicinal surfaces in crystal evaporation affected by electro- migration, Surface Science416, 200–213 (1998)

    S. Stoyanov, New type of step bunching instability at vicinal surfaces in crystal evaporation affected by electro- migration, Surface Science416, 200–213 (1998)

    work page 1998

  16. [16]

    Pimpinelli, V

    A. Pimpinelli, V. Tonchev, A. Videcoq, and M. Vladimirova, Scaling and universality of self- organized patterns on unstable vicinal surfaces, Physical Review Letters88, 206103 (2002)

    work page 2002

  17. [17]

    J. Krug, V. Tonchev, S. Stoyanov, and A. Pimpinelli, Scaling properties of step bunches induced by sublimation and related mechanisms, Physical Review B71, 045412 (2005)

    work page 2005

  18. [18]

    crea- tures

    M. Za luska-Kotur, H. Popova, and V. Tonchev, Step bunches, nanowires and other vicinal “crea- tures”—Ehrlich–Schwoebel effect by cellular automata, Crystals11, 1135 (2021)

    work page 2021

  19. [19]

    M. A. Chabowska and M. A. Za luska-Kotur, Engineering 2D surface patterns with the VicCa model, Journal of Crystal Growth667, 128245 (2025)

    work page 2025

  20. [20]

    M. A. Chabowska, H. Popova, and M. A. Za luska-Kotur, Step meandering: Balance between the potential well and the Ehrlich-Schwoebel barrier, Phys. Rev. B111, 155401 (2025)

    work page 2025

  21. [21]

    Chabowska and M

    M. Chabowska and M. Za luska-Kotur, Surface patterns shaped by additives in crystals, Vacuum240, 114431 (2025)

    work page 2025

  22. [22]

    Kandel and J

    D. Kandel and J. D. Weeks, Theory of impurity-induced step bunching, Physical Review B49, 5554–5564 (1994)

    work page 1994

  23. [23]

    J. D. Weeks, D.-J. Liu, and H.-C. Jeong, Two-dimensional models for step dynamics, inDynamics of Crystal Surfaces and Interfaces(Kluwer Academic Publishers, 1997) p. 199–216

    work page 1997

  24. [24]

    T. J. Stasevich and T. L. Einstein, Analytic formulas for the orientation dependence of step stiffness and line ten- sion: Key ingredients for numerical modeling, Multiscale Modeling and Simulation6, 90–104 (2007)

    work page 2007

  25. [25]

    C. R. Connell, Coarsening of step bunches in step flow growth: a reaction–diffusion model and its travelling wave solutions, Physica D: Nonlinear Phenomena189, 287–316 (2004)

    work page 2004

  26. [26]

    Margetis and R

    D. Margetis and R. E. Caflisch, Anisotropic step stiff- ness from a kinetic model of epitaxial growth, Multiscale Modeling & Simulation7, 242–273 (2008)

    work page 2008

  27. [27]

    Y.-M. Yu, A. Voigt, X. Guo, and Y. Liu, Simultaneous step meandering and bunching instabilities controlled by ehrlich-schwoebel barrier and elastic interaction, Applied Physics Letters99, 10.1063/1.3666781 (2011)

    work page doi:10.1063/1.3666781 2011

  28. [28]

    Krukowski, K

    S. Krukowski, K. Sakowski, P. Strak, P. Kempisty, J. Piechota, and I. Grzegory, Macrosteps dynamics and the growth of crystals and epitaxial layers, Progress in Crystal Growth and Characterization of Materials68, 100581 (2022)

    work page 2022

  29. [29]

    Higham, What is the second difference matrix? (2022)

    N. Higham, What is the second difference matrix? (2022)

    work page 2022

  30. [30]

    Saad and M

    Y. Saad and M. H. Schultz, GMRES: A generalized mini- mal residual algorithm for solving nonsymmetric linear systems, SIAM Journal on Scientific and Statistical Com- puting7, 856–869 (1986)

    work page 1986

  31. [31]

    Bradbury, R

    J. Bradbury, R. Frostig, P. Hawkins, M. J. Johnson, C. Leary, D. Maclaurin, G. Necula, A. Paszke, J. Van- derPlas, S. Wanderman-Milne, and Q. Zhang, JAX: com- posable transformations of Python+NumPy programs (2018)

    work page 2018

  32. [32]

    Rader, T

    J. Rader, T. Lyons, and P. Kidger, Lineax: uni- fied linear solves and linear least-squares in JAX and Equinox, AI for science workshop at Neural In- formation Processing Systems 2023, arXiv:2311.17283 10.48550/ARXIV.2311.17283 (2023)

    work page doi:10.48550/arxiv.2311.17283 2023

  33. [33]

    2024, arXiv:2402.09983

    J. Rader, T. Lyons, and P. Kidger, Optimistix: modu- lar optimisation in JAX and Equinox, arXiv:2402.09983 10.48550/ARXIV.2402.09983 (2024)

    work page doi:10.48550/arxiv.2402.09983 2024

  34. [34]

    Ivanov, Source code for the finite difference scheme for the continuum model (2026)

    V. Ivanov, Source code for the finite difference scheme for the continuum model (2026)

    work page 2026

  35. [35]

    F. C. Frank,Growth and Perfection of Crystals: Proceed- ings of an International Conference on Crystal Growth, Cooperstown, New York, August 1958(International Union of Crystallography (IUCr), Chester, England, 1958) pp. 411–421, proceedings of an International Conference on Crystal Growth, Cooperstown, New York, August 1958

    work page 1958

  36. [36]

    A. Stukowski, Visualization and analysis of atomistic sim- ulation data with OVITO-the Open Visualization Tool, Modelling and Simulation in Materials Science and Engi- neering18, 015012 (2010)

    work page 2010

  37. [37]

    Pimpinelli, J

    A. Pimpinelli, J. Villain, D. Wolf, J. M´ etois, J. Heyraud, I. Elkinani, and G. Uimin, Equilibrium step dynamics on vicinal surfaces, Surface Science295, 143–153 (1993)

    work page 1993

  38. [38]

    N. C. Bartelt, J. L. Goldberg, T. L. Einstein, E. D. Williams, J. C. Heyraud, and J. J. M´ etois, Brown- ian motion of steps on Si(111), Physical Review B48, 15453–15456 (1993)

    work page 1993

  39. [39]

    T. Ihle, C. Misbah, and O. Pierre-Louis, Equilibrium step dynamics on vicinal surfaces revisited, Physical Review B58, 2289–2309 (1998)

    work page 1998

  40. [40]

    Gillet, O

    F. Gillet, O. Pierre-Louis, and C. Misbah, Non-linear evo- lution of step meander during growth of a vicinal surface with no desorption, The European Physical Journal B18, 519–534 (2000)

    work page 2000