pith. sign in

arxiv: 2604.24388 · v1 · submitted 2026-04-27 · 🧮 math.AP

From PDEs on standard domains to self-similar particle systems on fractals

Georgi Medvedev , Emmanuel Tr\'elat (LJLL (UMR\_7598) , CaGE) This is my paper

Pith reviewed 2026-05-08 02:08 UTC · model grok-4.3

classification 🧮 math.AP
keywords fractalsPDE transportself-similar particlesisometryGalerkin methodnonlocal operatorsnumerical methods on fractals
0
0 comments X

The pith

Equations on the unit interval can be transported to self-similar fractals to produce approximating particle systems.

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

This paper shows how to transport partial differential equations from the unit interval onto self-similar fractal domains using a measure-preserving isometry between their L2 spaces. The transported equations are then discretized with a Galerkin method on self-similar partitions to obtain explicit self-similar systems of interacting particles. The resulting two-parameter scheme separates approximation error into a nonlocal consistency part and a Galerkin part. A reader cares because the reverse direction turns nonlocal fractal models into ordinary equations on the interval, where standard numerical methods apply directly.

Core claim

We construct transported PDEs on self-similar fractal domains from reference equations posed on the unit interval, and derive explicit self-similar interacting particle systems that approximate the resulting dynamics. The construction combines a measure-preserving isometry between L2-spaces on [0,1] and on the fractal, a nonlocal-to-local approximation of differential operators, and a Galerkin discretization on the canonical self-similar partitions. This yields a two-parameter approximation scheme whose error separates a nonlocal consistency term from a Galerkin network term. We work out examples for the transport, Burgers, and heat equations and outline extensions to local charts and pull-b

What carries the argument

Measure-preserving isometry between L2-spaces on the unit interval and fractal domain, combined with Galerkin discretization on self-similar partitions to produce particle systems.

Load-bearing premise

The measure-preserving isometry between the L2 spaces on the unit interval and the fractal domain allows the PDEs to be transported while retaining their key analytic properties.

What would settle it

For the transported heat equation on a specific fractal, compare the L2 error of the particle approximation against the exact or high-resolution solution as the discretization parameters increase; failure of the error to decrease according to the separated consistency and Galerkin terms would falsify the scheme.

Figures

Figures reproduced from arXiv: 2604.24388 by CaGE), Emmanuel Tr\'elat (LJLL (UMR\_7598), Georgi Medvedev.

Figure 1
Figure 1. Figure 1: The Sierpinski gasket, which serves as a prototype throughout the paper. view at source ↗
Figure 2
Figure 2. Figure 2: The isomorphism between measure spaces ( view at source ↗
Figure 3
Figure 3. Figure 3: a. The heat map of the function f(x) = e −|x1−x2| on SG G ⊂ R 2 (above) and it is representation as a function on [0, 1] (below). b. Graphon representation of W(x, y) = e −∥x−y∥ , (x, y) ∈ G × G ⊂ R 2 × R 2 . Next, define ˜fn : Q → R by ˜fn(x) = X w∈Σn νw(f)1Qw (x). It is shown in [20], that the sequence ( ˜fn) converges to ˜f ∈ L 1 (Q, λ) λ-a.e. and in L 1 (Q, λ). Moreover, ˜f = T f. Under the isomorphism… view at source ↗
read the original abstract

We construct transported PDEs on self-similar fractal domains from reference equations posed on the unit interval, and derive explicit self-similar interacting particle systems that approximate the resulting dynamics. The construction combines a measure-preserving isometry between $L^2$-spaces on $[0,1]$ and on the fractal \cite{Med2026}, a nonlocal-to-local approximation of differential operators \cite{PauTre2025}, and a Galerkin discretization on the canonical self-similar partitions. This yields a two-parameter approximation scheme whose error separates a nonlocal consistency term from a Galerkin network term. We work out the transport, Burgers, and heat equations, discuss the relation with intrinsic operators on fractals, and outline extensions to local charts and to pullbacks of nonlocal equations on fractal domains. Moreover, the reverse mapping transforms a nonlocal evolution equation on a fractal domains into the evolution equation on the unit interval, where the methods of classical numerical analysis can be applied. This suggests a promising direction for the development of numerical methods for nonlocal models on fractals, including fractional heat equation, fractal scattering, and related models.

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 / 1 minor

Summary. The manuscript constructs transported PDEs on self-similar fractal domains by pulling back reference equations from the unit interval [0,1] via a measure-preserving isometry between L² spaces, followed by a nonlocal-to-local approximation of differential operators and Galerkin discretization on canonical self-similar partitions. This produces a two-parameter approximation scheme whose error is claimed to separate into a nonlocal consistency term and a Galerkin network term. Explicit constructions are given for the transport, Burgers, and heat equations, with discussion of relations to intrinsic fractal operators and extensions to local charts and pullbacks of nonlocal equations; the reverse mapping is proposed to enable classical numerical methods for nonlocal fractal models.

Significance. If the central constructions hold, particularly the preservation of explicit self-similar particle systems under transport of nonlinear terms, the work would provide a systematic bridge from standard-domain PDE analysis to fractal domains. The two-parameter error separation (nonlocal consistency plus Galerkin) and the reverse-mapping strategy for applying interval-based numerics to fractional heat or scattering models on fractals represent potentially useful methodological contributions.

major comments (2)
  1. [Burgers equation section] The abstract asserts that explicit self-similar interacting particle systems are derived for the Burgers equation, yet the transported nonlinearity U((U^{-1}v)·D(U^{-1}v)) is not shown to remain inside the span of the Galerkin basis on the self-similar partitions or to inherit the required scaling for an explicit interaction kernel. The manuscript must supply the explicit form of the resulting particle system and verify that the quadratic term does not introduce additional error components that invalidate the claimed two-parameter separation.
  2. [error analysis for the two-parameter scheme] For the nonlinear case, the conjugation of the quadratic term through the isometry may disrupt the clean splitting of the approximation error into nonlocal consistency and Galerkin terms that is stated to hold for the two-parameter scheme. The manuscript should derive the error bound explicitly for Burgers (analogous to the linear transport and heat cases) and identify any cross terms arising from the composition.
minor comments (1)
  1. [Introduction] The citations to Med2026 and PauTre2025 are central; a brief self-contained recap of the isometry properties and the nonlocal approximation operator would improve readability without requiring the reader to consult the prior works.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments. We address each major comment below. We agree that the Burgers equation requires more explicit details on the particle system and error analysis, and these will be incorporated in the revised manuscript.

read point-by-point responses

  1. Referee: [Burgers equation section] The abstract asserts that explicit self-similar interacting particle systems are derived for the Burgers equation, yet the transported nonlinearity U((U^{-1}v)·D(U^{-1}v)) is not shown to remain inside the span of the Galerkin basis on the self-similar partitions or to inherit the required scaling for an explicit interaction kernel. The manuscript must supply the explicit form of the resulting particle system and verify that the quadratic term does not introduce additional error components that invalidate the claimed two-parameter separation.

    Authors: We acknowledge that while the general transported construction for Burgers is outlined via the isometry and Galerkin basis, the explicit form of the resulting self-similar particle system was not written out in full detail. The measure-preserving isometry and the self-similar choice of basis ensure the pulled-back nonlinearity remains in the finite-dimensional span and inherits the scaling for an explicit kernel; no additional error components arise beyond the two-parameter split. In the revision we will insert the explicit interaction kernel and the direct verification that the quadratic term preserves the claimed separation. revision: yes

  2. Referee: [error analysis for the two-parameter scheme] For the nonlinear case, the conjugation of the quadratic term through the isometry may disrupt the clean splitting of the approximation error into nonlocal consistency and Galerkin terms that is stated to hold for the two-parameter scheme. The manuscript should derive the error bound explicitly for Burgers (analogous to the linear transport and heat cases) and identify any cross terms arising from the composition.

    Authors: We agree that the nonlinear error bound must be derived explicitly. The linear cases follow immediately from the isometry and projection properties. For Burgers the quadratic term produces cross terms, but these are controlled by the Lipschitz continuity of the nonlinearity together with L2-orthogonality of the Galerkin basis. In the revision we will add the full error estimate for Burgers, explicitly identifying the cross terms and showing they remain absorbed into the existing nonlocal-consistency and Galerkin-network bounds. revision: yes

Circularity Check

0 steps flagged

No significant circularity; construction uses prior isometry as input but derives new approximations independently

full rationale

The paper's derivation chain starts from the cited measure-preserving isometry (Med2026) and nonlocal approximation (PauTre2025) as external inputs, then applies transport to obtain PDEs on the fractal and performs Galerkin discretization on self-similar partitions to produce explicit interacting particle systems. This does not reduce any claimed result to its inputs by construction: the particle systems and error separation (nonlocal consistency plus Galerkin terms) are obtained through additional steps that are not tautological with the isometry or prior citations. No self-definitional loops, fitted inputs renamed as predictions, or uniqueness theorems imported from the same authors appear. The self-citation supplies a mathematical tool rather than bearing the entire load of the central claim, leaving the derivation self-contained as a new synthesis.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 0 invented entities

The central claim depends on the isometry and approximation from cited literature, plus the application of Galerkin discretization to self-similar partitions; no new entities are postulated.

free parameters (1)
  • two parameters in the approximation scheme
    The scheme has two parameters controlling the nonlocal consistency term and the Galerkin network term in the error separation.
axioms (2)
  • domain assumption There exists a measure-preserving isometry between L2 spaces on [0,1] and the self-similar fractal
    This is used to transport the PDEs and is cited from Med2026.
  • domain assumption Nonlocal-to-local approximation of differential operators is valid
    Cited from PauTre2025 and used to approximate the operators.

pith-pipeline@v0.9.0 · 5505 in / 1557 out tokens · 72914 ms · 2026-05-08T02:08:54.492392+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

28 extracted references · 28 canonical work pages

  1. [1]

    B´ ar´ any, K

    B. B´ ar´ any, K. Simon, and B. Solomyak. Self-similar and self-affine sets and measures , volume 276 of Math. Surv. Monogr. Providence, RI: American Mathematical Society (AMS), 2023

    work page 2023

  2. [2]

    Baudoin and C

    F. Baudoin and C. Lacaux. Fractional gaussian fields on the Sierpi´ nski gasket and related fractals. J. Anal. Math., 146:719–739, 2022

    work page 2022

  3. [3]

    C. J. Bishop and Y. Peres. Fractals in probability and analysis, volume 162 ofCamb. Stud. Adv. Math. Cambridge: Cambridge University Press, 2017

    work page 2017

  4. [4]

    A. M. Caetano, S. N. Chandler-Wilde, X. Claeys, A. Gibbs, D. P. Hewett, and A. Moiola. Integral equation methods for acoustic scattering by fractals. Proceedings of the Royal Society A , 481(2306):20230650, 2025

    work page 2025

  5. [5]

    J. P. Chen, M. Hinz, and A. Teplyaev. From non-symmetric particle systems to non-linear PDEs on fractals. In Stochastic partial differential equations and related fields , volume 229 of Springer Proc. Math. Stat. , pages 503–513. Springer, Cham, 2018

    work page 2018

  6. [6]

    Cipriani and J.-L

    F. Cipriani and J.-L. Sauvageot. Derivations as square roots of Dirichlet forms. Journal of Functional Analysis , 201(1):78–120, 2003

    work page 2003

  7. [7]

    M. G. Crandall and A. Majda. Monotone difference approximations for scalar conservation laws. Math. Comp., 34(149):1–21, 1980

    work page 1980

  8. [8]

    C. M. Dafermos. Hyperbolic conservation laws in continuum physics , volume 325 of Grundlehren der mathema- tischen Wissenschaften. Springer, Berlin, fourth edition, 2016

    work page 2016

  9. [9]

    Fukushima and T

    M. Fukushima and T. Shima. On a spectral analysis for the Sierpi´ nski gasket. Potential Analysis, 1(1):1–35, 1992

    work page 1992

  10. [10]

    Grigor’yan and M

    A. Grigor’yan and M. Yang. Local and non-local Dirichlet forms on the Sierpi´ nski carpet. Transactions of the American Mathematical Society, 372(6):3985–4030, 2019

    work page 2019

  11. [11]

    M. Hinz, D. Koch, and M. Meinert. Sobolev spaces and calculus of variations on fractals. In Analysis, probability and mathematical physics on fractals. Based on the presentations at the 6th conference, Cornell University, Ithaca, NY, USA, June 2017 , pages 419–450. Hackensack, NJ: World Scientific, 2020

    work page 2017

  12. [12]

    Hinz and M

    M. Hinz and M. Meinert. Approximation of partial differential equations on compact resistance spaces. Calculus of Variations and Partial Differential Equations , 61(1):19, 2022

    work page 2022

  13. [13]

    M. Hinz, M. R¨ ockner, and A. Teplyaev. Vector analysis for Dirichlet forms and quasilinear PDE and SPDE on metric measure spaces. Stochastic Processes and their Applications , 123(12):4373–4406, 2013. 23

    work page 2013

  14. [14]

    Hinz and W

    M. Hinz and W. Schefer. On equations of continuity and transport type on metric graphs and fractals. arXiv preprint arXiv:2412.07988, 2024

    work page arXiv 2024

  15. [15]

    Hinz and A

    M. Hinz and A. Teplyaev. Vector analysis on fractals and applications. In Fractal geometry and dynamical systems in pure and applied mathematics II: Fractals in applied mathematics. , pages 147–163. Providence, RI: American Mathematical Society (AMS), 2013

    work page 2013

  16. [16]

    J. E. Hutchinson. Fractals and self-similarity. Indiana Univ. Math. J. , 30(5):713–747, 1981

    work page 1981

  17. [17]

    Kaliuzhnyi-Verbovetskyi and G

    D. Kaliuzhnyi-Verbovetskyi and G. S. Medvedev. Sparse Monte Carlo method for nonlocal diffusion problems. SIAM J. Numer. Anal. , 60(6):3001–3028, 2022

    work page 2022

  18. [18]

    J. Kigami. Analysis on fractals, volume 143 of Cambridge Tracts in Mathematics . Cambridge University Press, Cambridge, 2001

    work page 2001

  19. [19]

    J. Kigami. Measurable Riemannian geometry on the Sierpinski gasket: the Kusuoka measure and the Gaussian heat kernel estimate. Mathematische Annalen, 340(4):781–804, 2008

    work page 2008

  20. [20]

    G. S. Medvedev. Interacting dynamical systems on networks and fractals: discrete and continuous models, mean-field limit, and convergence rates. Preprint, arXiv:2601.23175 [math.DS] (2026), 2026

    work page arXiv 2026

  21. [21]

    G. S. Medvedev and M. S. Mizuhara. The Kuramoto model on the Sierpinski gasket II: twisted states. arXiv preprint arXiv:2506.12940, 2025

    work page arXiv 2025

  22. [22]

    G. S. Medvedev and M. S. Mizuhara. Kuramoto model on sierpinski gasket i: Harmonic maps, 2026

    work page 2026

  23. [23]

    U. Mosco. Analysis and numerics of some fractal boundary value problems. Boll. Unione Mat. Ital. (9) , 6(1):53– 73, 2013

    work page 2013

  24. [24]

    Paul and E

    T. Paul and E. Tr´ elat. Universal approximations of quasilinear PDEs by finite distinguishable particle systems. Preprint, arXiv:2501.11387 [math.AP] (2025), 2025

    work page arXiv 2025

  25. [25]

    R. S. Strichartz. Differential equations on fractals. Princeton University Press, Princeton, NJ, 2006. A tutorial

    work page 2006

  26. [26]

    D. W. Stroock. Mathematics of Probability, volume 149 of Grad. Stud. Math. Providence, RI: American Math- ematical Society (AMS), 2013

    work page 2013

  27. [27]

    Teplyaev

    A. Teplyaev. Harmonic coordinates on fractals with finitely ramified cell structure. Canadian Journal of Math- ematics, 60(2):457–480, 2008

    work page 2008

  28. [28]

    M. Yang. Equivalent semi-norms of non-local Dirichlet forms on the Sierpi´ nski gasket and applications.Potential Analysis, 49(2):287–308, 2018. 24

    work page 2018