Numerical simulations of the spread from the mean of the SLE and Multiple SLE dynamics
Pith reviewed 2026-06-27 14:35 UTC · model grok-4.3
The pith
Numerical experiments predict bimodal distributions for the spread from the mean in SLE when started near the origin, shifting to bell-shaped farther out, with multiple SLE always bell-shaped.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Numerical simulations of the SLE driven by Brownian motion show that the distributions of |g_t(z) - mean| and Re(g_t(z)) - Re(mean) are bimodal when the starting point is close to the origin and can transition to bell-shaped when started further away. For multiple SLE with Dyson Brownian motion drivers, the distributions are bell-shaped in all cases examined. These findings are obtained by applying Euler's method and provide numerical predictions for future theoretical studies on the spread from the mean behavior.
What carries the argument
The quantities |g_t(z) - \g_t(z)\| and Re(g_t(z)) - Re(\g_t(z)\), which measure the deviation of the conformal map from its sample average at fixed time in the Loewner dynamics.
Load-bearing premise
The Euler method discretization with the chosen time step faithfully approximates the underlying continuous Loewner dynamics without introducing artifacts that alter the distribution shapes, particularly near singularities.
What would settle it
Performing simulations with a significantly smaller time step or an alternative integration method such as Runge-Kutta and checking if the bimodal or bell-shaped nature of the distributions remains unchanged.
Figures
read the original abstract
The Schramm-Loewner Evolution (SLE) describes a family of fractal curves that arise in the study of the scaling limits of many planar Statistical Physics models. These curves are modeled using the Loewner Differential Equation for the conformal maps $g_t(z)$ with a Brownian motion driver. Using Euler's Method, in the current work we performed numerical experiments to study at a fixed time the quantities $|g_t(z) - \overline{g_t(z)}|$ and $Re(g_t(z)) - Re(\overline{g_t(z)})$, where $Re$ denotes the real part and $\overline{g_t(z)}$ refers to the sample average. These random variables measure the 'spread' of the dynamics from the average behavior at fixed time. One of the scopes of this work is to give numerical predictions for future theoretical investigations on these quantities. When investigating these quantities in the SLE case our experiments predict that the distribution is bimodal when the dynamics started close to the origin, and it can become bell-shaped if the dynamics is started further from the origin. In the second part, we performed experiments for a Multiple SLE model whose driver is Dyson Brownian Motion. Due to singularity in the dynamics of the drivers and the many data points needed, this part is challenging from a computational perspective. In the multiple SLE case, our experiments predict that the distribution is bell-shaped in all cases. In addition, we check the changes in the distributions as we vary the parameter $\kappa$ in the SLE case and $\beta$ in the Multiple SLE case.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper performs forward Euler numerical integration of the Loewner ODE for standard SLE (Brownian driver) and multiple SLE (Dyson Brownian driver) at fixed time t. It reports the empirical distributions of |g_t(z) − ar g_t(z)| and Re(g_t(z)) − Re(ar g_t(z)) for varying starting points z, claiming these are bimodal for SLE started near the origin and bell-shaped for SLE started farther away or for all multiple-SLE cases; it also examines dependence on κ and β and presents the shapes as predictions for future analytic work.
Significance. If the reported distribution shapes survive proper numerical validation, the work supplies concrete empirical targets that could guide analytic study of deviation statistics in Loewner evolutions. The computational challenges near singularities are acknowledged, but the lack of documented discretization control and statistical error assessment limits the strength of the predictions.
major comments (2)
- [Numerical experiments / Methods] Numerical experiments (Euler integration of the Loewner ODE): the time step Δt is never stated and no mesh-refinement study, local truncation-error bound, or comparison with a higher-order integrator is supplied. Because the modality claims rest entirely on the observed histograms, this omission is load-bearing; the shapes could be discretization artifacts, especially when drivers approach the real axis or multiple poles interact.
- [Results and figures] Results on distribution shapes: no sample sizes, error bars, or convergence diagnostics with respect to number of Monte Carlo realizations are reported, nor is any benchmark against known SLE properties (e.g., exact moments for small κ or chordal vs. radial cases) provided. Without these, the bimodal-versus-bell-shaped distinction cannot be assessed for statistical reliability.
minor comments (2)
- [Abstract] The abstract states that the multiple-SLE part is “challenging from a computational perspective” due to singularities and data volume, yet no quantitative information on realized sample counts, run times, or singularity-handling strategy is given.
- [Introduction / Notation] Notation: ar g_t(z) is introduced as the sample average; it should be stated explicitly whether this is an ensemble average over independent realizations or a time average, and how many realizations enter the average.
Simulated Author's Rebuttal
We thank the referee for the careful reading and constructive comments on our numerical study of spread distributions in SLE and multiple SLE. We address each major comment below and will revise the manuscript accordingly to improve documentation and statistical validation.
read point-by-point responses
-
Referee: [Numerical experiments / Methods] Numerical experiments (Euler integration of the Loewner ODE): the time step Δt is never stated and no mesh-refinement study, local truncation-error bound, or comparison with a higher-order integrator is supplied. Because the modality claims rest entirely on the observed histograms, this omission is load-bearing; the shapes could be discretization artifacts, especially when drivers approach the real axis or multiple poles interact.
Authors: We agree that the discretization parameters must be explicitly documented. In the revised manuscript we will state the fixed time step Δt used for the forward Euler integration, include a mesh-refinement study for representative starting points and κ values showing that the reported bimodal and bell-shaped features persist under refinement, and briefly justify the choice of Euler scheme given the computational constraints near singularities already noted in the text. revision: yes
-
Referee: [Results and figures] Results on distribution shapes: no sample sizes, error bars, or convergence diagnostics with respect to number of Monte Carlo realizations are reported, nor is any benchmark against known SLE properties (e.g., exact moments for small κ or chordal vs. radial cases) provided. Without these, the bimodal-versus-bell-shaped distinction cannot be assessed for statistical reliability.
Authors: We accept that sample sizes, error estimates, and convergence information are required for assessing reliability. The revised version will report the number of Monte Carlo realizations for each histogram, add error bars, and include diagnostics showing stability with increasing sample count. While the spread quantities |g_t(z)−ḡ_t(z)| and Re(g_t(z))−Re(ḡ_t(z)) lack known closed-form distributions, we will add supporting benchmarks such as the deterministic κ=0 limit and small-κ moment comparisons where analytic results are available. revision: yes
Circularity Check
No circularity: results are direct outputs of numerical integration with no fitted parameters or self-referential derivations
full rationale
The manuscript contains no analytic derivation chain, no parameter fitting, and no self-citations. All reported distribution shapes (bimodal or bell-shaped) are obtained by forward Euler integration of the Loewner ODE driven by Brownian motion or Dyson Brownian motion, followed by direct histogram construction on the simulated trajectories of |g_t(z) - ḡ_t(z)| and Re(g_t(z)) - Re(ḡ_t(z)). These quantities are defined independently of the simulation outputs; the paper simply records their empirical statistics at fixed t. Because no quantity is defined in terms of another or obtained by fitting a subset of the same data, none of the enumerated circularity patterns apply. The work is therefore self-contained against external benchmarks (the continuous Loewner dynamics) and receives the default non-finding.
Axiom & Free-Parameter Ledger
free parameters (4)
- Euler time step
- Number of Monte Carlo samples
- Fixed time t and starting point z
- kappa and beta
axioms (2)
- domain assumption The Loewner differential equation driven by Brownian motion (or Dyson Brownian motion) generates the SLE (multiple SLE) curves.
- standard math Standard properties of Brownian motion and stochastic differential equations hold.
Reference graph
Works this paper leans on
-
[1]
Distribution Fitting with Sum of Square Error (SSE).url:https : / / stackoverflow.com/a/37616966
Timothy Davenport (https://stackoverflow.com/users/2087463/tmthydvnprt). Distribution Fitting with Sum of Square Error (SSE).url:https : / / stackoverflow.com/a/37616966
-
[2]
Pole dynamics and an integral of motion for multiple SLE(0)
Tom Alberts et al. “Pole dynamics and an integral of motion for multiple SLE(0)”. In:Selecta Math. (N.S.)30.5 (2024), Paper No. 92, 77.issn: 1022-1824,1420-9020.doi:10.1007/s00029-024-00980-0.url:https: //doi.org/10.1007/s00029-024-00980-0
-
[3]
On the uniqueness of global multiple SLEs
Vincent Beffara, Eveliina Peltola, and Hao Wu. “On the uniqueness of global multiple SLEs”. In:Ann. Probab.49.1 (2021), pp. 400–434.issn: 0091-1798,2168-894X.doi:10.1214/20- AOP1477.url:https://doi. org/10.1214/20-AOP1477. 10
work page doi:10.1214/20- 2021
-
[4]
Continuity inκ in SLE κ theory using a constructive method and rough path theory
Dmitry Beliaev, Terry J. Lyons, and Vlad Margarint. “Continuity inκ in SLE κ theory using a constructive method and rough path theory”. In:Ann. Inst. Henri Poincar´ e Probab. Stat.57.1 (2021), pp. 455–468. issn: 0246-0203,1778-7017.doi:10 . 1214 / 20 - aihp1084.url:https : //doi.org/10.1214/20-aihp1084
-
[5]
Continuity of zero- hitting times of Bessel processes and welding homeomorphisms of SLE κ
Dmitry Beliaev, Vlad Margarint, and Atul Shekhar. “Continuity of zero- hitting times of Bessel processes and welding homeomorphisms of SLE κ”. In:ALEA Lat. Am. J. Probab. Math. Stat.18.1 (2021), pp. 69–79.issn: 1980-0436.doi:10.30757/alea.v18- 04.url:https://doi.org/10. 30757/alea.v18-04
-
[6]
Global fluctuations in generalβDyson’s Brownian mo- tion
Martin Bender. “Global fluctuations in generalβDyson’s Brownian mo- tion”. In:Stochastic Processes and their Applications118.6 (2008), pp. 1022– 1042.issn: 0304-4149.doi:https://doi.org/10.1016/j.spa.2007.07. 010.url:https://www.sciencedirect.com/science/article/pii/ S0304414907001329
-
[7]
Rate of convergence in multiple SLE using random matrix theory
Andrew Campbell, Kyle Luh, and Vlad Margarint. “Rate of convergence in multiple SLE using random matrix theory”. In:Random Matrices The- ory Appl.14.1 (2025), Paper No. 2450028, 23.issn: 2010-3263,2010-3271. doi:10.1142/S201032632450028X.url:https://doi.org/10.1142/ S201032632450028X
work page doi:10.1142/s201032632450028x.url:https://doi.org/10.1142/ 2025
-
[8]
Stochastic Loewner evolution and Dyson’s circular ensem- bles
John Cardy. “Stochastic Loewner evolution and Dyson’s circular ensem- bles”. In:J. Phys. A36.24 (2003), pp. L379–L386.issn: 0305-4470,1751- 8121.doi:10.1088/0305-4470/36/24/101.url:https://doi.org/10. 1088/0305-4470/36/24/101
work page doi:10.1088/0305-4470/36/24/101.url:https://doi.org/10 2003
-
[9]
Perturbations of multiple Schramm- Loewner evolution with two non-colliding Dyson Brownian motions
Jiaming Chen and Vlad Margarint. “Perturbations of multiple Schramm- Loewner evolution with two non-colliding Dyson Brownian motions”. In: Stochastic Process. Appl.151 (2022), pp. 553–569.issn: 0304-4149,1879- 209X.doi:10.1016/j.spa.2022.06.014.url:https://doi.org/10. 1016/j.spa.2022.06.014
work page doi:10.1016/j.spa.2022.06.014.url:https://doi.org/10 2022
-
[10]
Splitting algorithm and normed con- vergence for drawing the random Loewner curves
Jiaming Chen and Vlad Margarint. “Splitting algorithm and normed con- vergence for drawing the random Loewner curves”. In:Proc. A.481.2320 (2025), Paper No. 20240929, 22.issn: 1364-5021,1471-2946
2025
-
[11]
Commutation relations for Schramm-Loewner evolu- tions
Julien Dub´ edat. “Commutation relations for Schramm-Loewner evolu- tions”. In:Comm. Pure Appl. Math.60.12 (2007), pp. 1792–1847.issn: 0010-3640,1097-0312.doi:10.1002/cpa.20191.url:https://doi.org/ 10.1002/cpa.20191
-
[12]
Courant lecture notes in math- ematics ; 28
L´ aszl´ o Erd˝ os and Horng-Tzer Yau.A dynamical approach to random ma- trix theory / L´ aszl´ o Erd˝ os, Institute of Science and Technology Austria ; Horng-Tzer Yau, Harvard University.eng. Courant lecture notes in math- ematics ; 28. New York, New York: Courant Institute of Mathematical Sciences, New York University, 2017 - 2017.isbn: 9781470436483. 11
2017
-
[13]
Random matrix theory: Dyson Brownian motion
Gianluca Finocchio. “Random matrix theory: Dyson Brownian motion”. en. In:Snapshots of modern mathematics from Oberwolfach;2020(2020), p. 02.doi:10.14760/SNAP-2020-002-EN.url:http://publications. mfo.de/handle/mfo/3736
work page doi:10.14760/snap-2020-002-en.url:http://publications 2020
-
[14]
James Foster, Terry Lyons, and Vlad Margarint. “An asymptotic radius of convergence for the Loewner equation and simulation ofSLE κ traces via splitting”. In:J. Stat. Phys.189.2 (2022), Paper No. 18, 14.issn: 0022-4715,1572-9613.doi:10.1007/s10955-022-02979-3.url:https: //doi.org/10.1007/s10955-022-02979-3
-
[15]
Regularity of SLE in (t, κ) and refined GRR estimates
Peter K. Friz, Huy Tran, and Yizheng Yuan. “Regularity of SLE in (t, κ) and refined GRR estimates”. In:Probab. Theory Related Fields180.1-2 (2021), pp. 71–112.issn: 0178-8051,1432-2064.doi:10 . 1007 / s00440 - 021-01058-0.url:https://doi.org/10.1007/s00440-021-01058-0
-
[16]
N-sided radial Schramm- Loewner evolution
Vivian Olsiewski Healey and Gregory F. Lawler. “N-sided radial Schramm- Loewner evolution”. In:Probab. Theory Related Fields181.1-3 (2021), pp. 451–488.issn: 0178-8051,1432-2064.doi:10 . 1007 / s00440 - 021 - 01033-9.url:https://doi.org/10.1007/s00440-021-01033-9
-
[17]
Absence of Local Conserved Quantity in the Heisenberg Model with Next-Nearest-Neighbor Interaction.J
Ikkei Hotta and Makoto Katori. “Hydrodynamic Limit of Multiple SLE”. In:Journal of Statistical Physics171.1 (Mar. 2018), pp. 166–188.issn: 1572-9613.doi:10.1007/s10955- 018- 1996- y.url:http://dx.doi. org/10.1007/s10955-018-1996-y
-
[18]
Fredrik Johansson Viklund, Steffen Rohde, and Carto Wong. “On the continuity of SLEκ inκ”. In:Probab. Theory Related Fields159.3-4 (2014), pp. 413–433.issn: 0178-8051,1432-2064.doi:10 . 1007 / s00440 - 013 - 0506-z.url:https://doi.org/10.1007/s00440-013-0506-z
-
[19]
Antti Kemppainen.Schramm-Loewner evolution. Vol. 24. SpringerBriefs in Mathematical Physics. Springer, Cham, 2017, pp. ix+145.isbn: 978- 3-319-65327-3; 978-3-319-65329-7.doi:10.1007/978- 3- 319- 65329- 7. url:https://doi.org/10.1007/978-3-319-65329-7
-
[20]
Numerical computations for the Schramm-Loewner evolu- tion
Tom Kennedy. “Numerical computations for the Schramm-Loewner evolu- tion”. In:J. Stat. Phys.137.5-6 (2009), pp. 839–856.issn: 0022-4715,1572- 9613.doi:10.1007/s10955-009-9866-2.url:https://doi.org/10. 1007/s10955-009-9866-2
work page doi:10.1007/s10955-009-9866-2.url:https://doi.org/10 2009
-
[21]
Lawler.Conformally invariant processes in the plane
Gregory F. Lawler.Conformally invariant processes in the plane. Vol. 114. Mathematical Surveys and Monographs. American Mathematical Society, Providence, RI, 2005, pp. xii+242.isbn: 0-8218-3677-3.doi:10.1090/ surv/114.url:https://doi.org/10.1090/surv/114
-
[22]
Conformal in- variance of planar loop-erased random walks and uniform spanning trees
Gregory F. Lawler, Oded Schramm, and Wendelin Werner. “Conformal in- variance of planar loop-erased random walks and uniform spanning trees”. In:The Annals of Probability32.1B (2004), pp. 939–995.doi:10.1214/ aop/1079021469.url:https://doi.org/10.1214/aop/1079021469. 12
-
[23]
Numerical solution of Dyson Brow- nian motion and a sampling scheme for invariant matrix ensembles
Xingjie Helen Li and Govind Menon. “Numerical solution of Dyson Brow- nian motion and a sampling scheme for invariant matrix ensembles”. In: J. Stat. Phys.153.5 (2013), pp. 801–812.issn: 0022-4715,1572-9613.doi: 10 . 1007 / s10955 - 013 - 0858 - x.url:https : / / doi . org / 10 . 1007 / s10955-013-0858-x
2013
-
[24]
Berestycki & J.R
N. Berestycki & J.R. Norris.Lectures on Schramm–Loewner Evolution. Jan. 2016
2016
-
[25]
eprint:https://arxiv.org/pdf/2407.13762
Eveliina Peltola Osama Abuzaid Vivian Olsiewski Healey.Large devia- tions of Dyson Brownian motion on the circle and multiradial SLE0+. eprint:https://arxiv.org/pdf/2407.13762
-
[26]
Scaling limits of loop-erased random walks and uniform spanning trees
Oded Schramm. “Scaling limits of loop-erased random walks and uniform spanning trees”. In:Israel Journal of Mathematics118.1 (Dec. 2000), pp. 221–288.issn: 1565-8511.doi:10 . 1007 / bf02803524.url:http : //dx.doi.org/10.1007/BF02803524
-
[27]
Conformal invariance in random cluster models. I. Holmorphic fermions in the Ising model
Stanislav Smirnov. “Conformal invariance in random cluster models. I. Holmorphic fermions in the Ising model”. In:Annals of Mathematics172.2 (2010), p. 1435.issn: 0003-486X.doi:10.4007/annals.2010.172.1435. url:http://dx.doi.org/10.4007/annals.2010.172.1435
-
[28]
Stanislav Smirnov.Critical percolation in the plane. 2009. arXiv:0909. 4499 [math.PR].url:https://arxiv.org/abs/0909.4499
Pith/arXiv arXiv 2009
-
[29]
Existence and Uniqueness of Nonsimple Multiple SLE
Dapeng Zhan. “Existence and Uniqueness of Nonsimple Multiple SLE”. In:J. Stat. Phys.191.8 (2024), Paper No. 101.issn: 0022-4715,1572-9613. doi:10.1007/s10955-024-03306-8.url:https://doi.org/10.1007/ s10955-024-03306-8. A Numerical experiments studying the spread for the SLE maps dynamics started fromz 0 = 1.02i (a)κ= 1 (b)κ= 2. Figure 8: Histograms forRe(...
work page doi:10.1007/s10955-024-03306-8.url:https://doi.org/10.1007/ 2024
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.