pith. sign in

arxiv: 2511.18256 · v2 · submitted 2025-11-23 · ❄️ cond-mat.stat-mech

First-Passage Times for the Space-Fractional Spectral Fokker-Planck Equation

Pith reviewed 2026-05-17 06:45 UTC · model grok-4.3

classification ❄️ cond-mat.stat-mech
keywords first-passage timesspace-fractional Fokker-Planck equationcompounded stepssuperdiffusive processesabsorbing boundaryasymptotic scalingLévy flights
0
0 comments X

The pith

Extending random walks to compounded steps yields first-passage time scalings distinct from Lévy flights.

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

The paper extends the random walk framework to include compounded steps. This provides first-passage time properties for superdiffusive processes governed by the space-fractional spectral Fokker-Planck equation. The approach accounts for space-dependent forces and boundaries hit during jumps. For the one-sided absorbing boundary on the semi-infinite line with no potential, the first-passage time density scales as t to the power -1/(2α)-1 at large times. This agrees with the method of images but differs from the t^{-3/2} scaling for Lévy flights of order 2α, and an optimal α minimizes the mean first-passage time.

Core claim

We extend the random walk framework to include compounded steps, providing first-passage time properties for a new class of superdiffusive processes governed by the space-fractional spectral Fokker-Planck equation. This first-passage process leads to novel FPT properties, different from Lévy flights, that account for space dependent forces and hitting boundaries throughout the path of a jump. For the one-sided absorbing boundary with no potential on the semi-infinite line, we find that the FPT density scales asymptotically as t^{-1/(2α)-1} for large times, where the parameter α ∈ (0,1] relates to the power-law behavior for the distribution of the number of compounded steps. This is in agreem

What carries the argument

Compounded steps in random walks leading to the space-fractional spectral Fokker-Planck equation, which enables first-passage time calculations that incorporate forces and mid-jump boundary hits.

Load-bearing premise

Compounded steps generate processes that obey the space-fractional spectral Fokker-Planck equation while properly handling space-dependent forces and boundaries encountered along jump paths.

What would settle it

Numerical simulation of the random walk with compounded steps on the semi-infinite domain with an absorbing boundary at zero, measuring the tail of the first-passage time distribution to check for the power-law exponent -1/(2α)-1.

Figures

Figures reproduced from arXiv: 2511.18256 by Boris Z. Huang, Bruce I. Henry, Christopher N. Angstmann, Daniel S. Han.

Figure 1
Figure 1. Figure 1: FIG. 1. Example realization of a discrete time random walk [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. The PDF (21) for finding the particle at [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Plot of the survival function (25) (top) and PDF [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Log-log plot for the PDF (26) (solid lines) and its [PITH_FULL_IMAGE:figures/full_fig_p005_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. The expected FPT for the one-sided absorbing diffu [PITH_FULL_IMAGE:figures/full_fig_p006_5.png] view at source ↗
read the original abstract

We extend the random walk framework to include compounded steps, providing first-passage time (FPT) properties for a new class of superdiffusive processes, which are governed by the space-fractional spectral Fokker-Planck equation. This first-passage process leads to novel FPT properties, different from L\'evy flights, that account for space dependent forces and hitting boundaries throughout the path of a jump. The FPT distribution can be derived for different types of barriers and potentials, for which we also provide specific examples. For the one-sided absorbing boundary with no potential on the semi-infinite line, we find that the FPT density scales asymptotically as $t^{-1/(2\alpha)-1}$ for large times, where the parameter $\alpha \in (0,1]$ relates to the power-law behavior for the distribution of the number of compounded steps. This is in agreement with the method of images but different to the Sparre-Andersen scaling $t^{-3/2}$ for corresponding L\'evy flights of order $2\alpha$. In this case, there exists an optimal space-fractional exponent $\alpha$ to minimize the mean FPT.

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 extends the random walk framework to compounded steps whose number follows a power-law distribution parameterized by α ∈ (0,1]. The continuum limit is asserted to be the space-fractional spectral Fokker-Planck equation. First-passage time (FPT) properties are derived for various barriers and potentials; the central explicit result is that, for a one-sided absorbing boundary on the semi-infinite line with no potential, the asymptotic FPT density scales as t^{-1/(2α)-1} for large t. This scaling is stated to agree with the method of images yet to differ from the Sparre-Andersen t^{-3/2} law for Lévy flights of index 2α. An optimal α minimizing the mean FPT is also identified.

Significance. If the boundary condition is shown to arise directly from the discrete compounded-step dynamics and to remain non-standard in the continuum limit, the work supplies a new family of superdiffusive processes whose FPT statistics incorporate continuous boundary monitoring during jumps and space-dependent forces. The α-dependent exponent and the existence of an optimal α constitute concrete, falsifiable predictions that distinguish the model from both standard Lévy flights and ordinary fractional diffusion. The agreement with the method of images is a positive consistency check.

major comments (2)
  1. [Section deriving the one-sided absorbing boundary case and the associated boundary condition] The load-bearing step for the claimed scaling t^{-1/(2α)-1} is the assertion that the space-fractional spectral Fokker-Planck equation, together with the boundary condition inherited from the compounded-step model, encodes hits throughout the path of a jump. The manuscript must supply an explicit derivation (or at least a clear statement) of this boundary condition and demonstrate that it does not reduce to a standard Dirichlet condition on the fractional operator; otherwise the long-time survival probability reverts to the universal t^{-1/2} decay independent of α, eliminating the reported α dependence.
  2. [Results section on the semi-infinite line with one-sided absorber] No derivation steps, error estimates, or direct comparison with discrete simulations are provided for the asymptotic FPT scaling. A short analytic or numerical verification that the continuum FPT density indeed follows t^{-1/(2α)-1} (rather than being fitted post hoc) is required to substantiate the central claim.
minor comments (2)
  1. The abstract states that 'specific examples' are given for different barriers and potentials; the manuscript should ensure these examples are accompanied by explicit formulas, figures, or tables so that readers can reproduce the claimed FPT distributions.
  2. The relation between the discrete parameter α and the fractional order appearing in the spectral Fokker-Planck equation should be stated more explicitly in the model section to avoid ambiguity.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive and detailed report. The comments highlight important points where additional clarity and verification will strengthen the manuscript. We address each major comment below and will incorporate the necessary revisions.

read point-by-point responses
  1. Referee: [Section deriving the one-sided absorbing boundary case and the associated boundary condition] The load-bearing step for the claimed scaling t^{-1/(2α)-1} is the assertion that the space-fractional spectral Fokker-Planck equation, together with the boundary condition inherited from the compounded-step model, encodes hits throughout the path of a jump. The manuscript must supply an explicit derivation (or at least a clear statement) of this boundary condition and demonstrate that it does not reduce to a standard Dirichlet condition on the fractional operator; otherwise the long-time survival probability reverts to the universal t^{-1/2} decay independent of α, eliminating the reported α dependence.

    Authors: We agree that an explicit derivation of the boundary condition is required to support the central claim. In the revised manuscript we will add a dedicated subsection deriving the boundary condition directly from the discrete compounded-step dynamics. The derivation will show that absorption occurs whenever any intermediate point along a random jump trajectory intersects the boundary, yielding a non-local condition on the fractional operator that does not reduce to a standard Dirichlet condition. We will further demonstrate analytically that this condition produces the α-dependent survival decay t^{-1/(2α)} rather than the universal t^{-1/2} behavior. revision: yes

  2. Referee: [Results section on the semi-infinite line with one-sided absorber] No derivation steps, error estimates, or direct comparison with discrete simulations are provided for the asymptotic FPT scaling. A short analytic or numerical verification that the continuum FPT density indeed follows t^{-1/(2α)-1} (rather than being fitted post hoc) is required to substantiate the central claim.

    Authors: We accept that the current presentation lacks sufficient verification steps. The revised version will include a concise analytic derivation of the long-time asymptotic using the method of images for the space-fractional spectral operator subject to the inherited boundary condition, together with the resulting exponent -1/(2α)-1. In addition, we will provide Monte Carlo simulations of the underlying discrete compounded random walk on the semi-infinite line, direct comparison of the empirical FPT density with the continuum prediction, and quantitative error estimates confirming the scaling. revision: yes

Circularity Check

0 steps flagged

Derivation self-contained from compounded-step random walk extension

full rationale

The paper extends the random walk framework to compounded steps whose continuum limit yields the space-fractional spectral Fokker-Planck equation and associated boundary conditions. The asymptotic FPT scaling t^{-1/(2α)-1} is obtained by solving this equation for the one-sided absorbing boundary on the semi-infinite line (no potential), with the result stated to agree with the method of images while differing from Sparre-Andersen. No quoted step reduces the central claim to a fitted parameter, self-definition, or load-bearing self-citation chain; the governing equation and boundary handling are presented as direct consequences of the random-walk construction rather than inputs that presuppose the target scaling. The derivation therefore remains independent of the result it produces.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The central claim rests on the validity of mapping compounded random walks onto the space-fractional spectral Fokker-Planck equation and on the applicability of the method of images for the one-sided boundary problem; these are introduced without independent derivation in the provided abstract.

free parameters (1)
  • α
    Parameter in (0,1] that controls the power-law distribution of the number of compounded steps; its value is chosen to produce the desired superdiffusive behavior and appears in the final scaling exponent.
axioms (1)
  • domain assumption Compounded-step random walks generate processes governed by the space-fractional spectral Fokker-Planck equation that correctly handle space-dependent forces and boundary hits during jumps.
    Invoked as the governing equation for the new class of superdiffusive processes in the abstract.

pith-pipeline@v0.9.0 · 5522 in / 1579 out tokens · 47438 ms · 2026-05-17T06:45:48.038981+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

48 extracted references · 48 canonical work pages · 1 internal anchor

  1. [1]

    The FPT remains a topic of increas- ing significance within numerous other recent studies on machine learning algorithms [13, 14], trading strategies

    and in animal foraging strategies [8, 9], where the for- aging pattern of animals dictates whether they find food and shelter [10–12]. The FPT remains a topic of increas- ing significance within numerous other recent studies on machine learning algorithms [13, 14], trading strategies

  2. [2]

    and quantum tunneling [16–18]. The FPT of random walks exhibiting superdiffusive characteristics, where the mean square displacement (MSD) scales as⟨x 2(t)⟩ ∼t ν forν >1, has been of recent interest in theoretical studies. Examples include: the op- timization of searching for sparse targets [19], searches with external biases [20, 21], searches with short...

  3. [3]

    On the other hand, the expected FPT of the L´ evy flight on the half-line is divergent for all 0< α≤1

    However, forα≥ 1 2, the expected FPT is divergent due to the behavior of the PDF power-law tail. On the other hand, the expected FPT of the L´ evy flight on the half-line is divergent for all 0< α≤1. This is because the PDF of the FPT for the L´ evy flight followsψ(t)∼t −3/2 [41–43] as seen in Figure 4. For the compounded random walk, Figure 5 shows how t...

  4. [4]

    Redner,A guide to first-passage processes(Cambridge University Press, 2001)

    S. Redner,A guide to first-passage processes(Cambridge University Press, 2001)

  5. [5]

    Szabo, K

    A. Szabo, K. Schulten, and Z. Schulten, The Journal of Chemical Physics72, 4350 (1980)

  6. [6]

    R. J. Preston, M. F. Gelin, and D. S. Kosov, The Journal of Chemical Physics154, 10.1063/5.0045652 (2021)

  7. [7]

    S. Park, M. K. Sener, D. Lu, and K. Schulten, The Jour- nal of Chemical Physics119, 1313 (2003)

  8. [8]

    C. Rao, D. Waxman, W. Lin, and Z. Song, The Journal of Chemical Physics162(2025)

  9. [9]

    Iyer-Biswas and A

    S. Iyer-Biswas and A. Zilman, Advances in Chemical Physics160, 261 (2016)

  10. [10]

    M. Reva, D. A. DiGregorio, and D. S. Grebenkov, Scien- tific Reports11, 5377 (2021)

  11. [11]

    James, J

    A. James, J. W. Pitchford, and M. J. Plank, Bulletin of Mathematical Biology72, 896 (2010). 7

  12. [12]

    Z. P. Kilpatrick, J. D. Davidson, and A. El Hady, Journal of The Royal Society Interface18, 20210337 (2021)

  13. [13]

    B´ enichou, M

    O. B´ enichou, M. Coppey, M. Moreau, P.-H. Suet, and R. Voituriez, Phys. Rev. Lett.94, 198101 (2005)

  14. [14]

    B´ enichou, C

    O. B´ enichou, C. Loverdo, M. Moreau, and R. Voituriez, Rev. Mod. Phys.83, 81 (2011)

  15. [15]

    Raposo and M

    E. Raposo and M. Da Luz, Physics of Life Reviews5, 133 (2008)

  16. [16]

    S. Sun, Z. Cao, H. Zhu, and J. Zhao, IEEE Transactions on Cybernetics50, 3668 (2020)

  17. [17]

    Y. Yang, H. Chen, A. A. Heidari, and A. H. Gandomi, Expert Systems with Applications177, 114864 (2021)

  18. [18]

    T. Jin, Y. Zhu, Y. Shu, J. Cao, H. Yan, and D. Jiang, Journal of Intelligent & Fuzzy Systems44, 1585 (2023)

  19. [19]

    A. M. Zheltikov, Opt. Lett.48, 5755 (2023)

  20. [20]

    Ni and Y

    Z. Ni and Y. Zheng, Entropy25, 10.3390/e25081231 (2023)

  21. [21]

    Moosavi, H

    M. Moosavi, H. Zaki Dizaji, and A. Seyyedi, Pramana 99, 1 (2024)

  22. [22]

    M. R. Palyulin VV, Chechkin AV, Proc Natl Acad Sci U S A111, 2931 (2014)

  23. [23]

    V. V. Palyulin, A. V. Chechkin, and R. Metzler, Journal of Statistical Mechanics: Theory and Experiment2014, P11031 (2014)

  24. [24]

    Padash, K

    A. Padash, K. Capa la, H. Kantz, B. Dybiec, B. Shokri, R. Metzler, and A. V. Chechkin, Journal of Physics A: Mathematical and Theoretical58, 185002 (2025)

  25. [25]

    Boyer, G

    D. Boyer, G. Mercado-V´ asquez, S. N. Majumdar, and G. Schehr, Physical Review E109, L022103 (2024)

  26. [26]

    Padash, A

    A. Padash, A. V. Chechkin, B. Dybiec, I. Pavlyukevich, B. Shokri, and R. Metzler, Journal of Physics A: Math- ematical and Theoretical52, 454004 (2019)

  27. [27]

    Padash, T

    A. Padash, T. Sandev, H. Kantz, R. Metzler, and A. V. Chechkin, Fractal and Fractional6, 260 (2022)

  28. [28]

    S. D. Lawley, Journal of Nonlinear Science33, 53 (2023)

  29. [29]

    Pal and S

    A. Pal and S. Reuveni, Physical Review Letters118, 030603 (2017)

  30. [30]

    Mudra and A

    S. Mudra and A. Chakraborty, Chemical Physics Letters 751, 137531 (2020)

  31. [31]

    Wang and D

    D. Wang and D. K. Schwartz, The Journal of Physical Chemistry C124, 19880 (2020)

  32. [32]

    Viswanathan, S

    G. Viswanathan, S. Buldyrev, S. Havlin, M. da Luz, E. Raposo, and H. Stanley, Nature401, 911 (1999)

  33. [33]

    G. M. Viswanathan, M. G. Da Luz, E. P. Raposo, and H. E. Stanley,The physics of foraging: an introduction to random searches and biological encounters(Cambridge University Press, 2011)

  34. [34]

    Levernier, J

    N. Levernier, J. Textor, O. B´ enichou, and R. Voituriez, Physical Review Letters124, 080601 (2020)

  35. [35]

    C. N. Angstmann, D. S. Han, B. I. Henry, B. Z. Huang, and Z. Xu, Phys. Rev. E111, 024136 (2025)

  36. [36]

    Lischke, G

    A. Lischke, G. Pang, M. Gulian, F. Song, C. Glusa, X. Zheng, Z. Mao, W. Cai, M. M. Meerschaert, M. Ainsworth, and G. E. Karniadakis, Journal of Com- putational Physics404, 109009 (2020)

  37. [37]

    I. M. Sokolov, J. Klafter, and A. Blumen, Phys. Rev. E 64, 021107 (2001)

  38. [38]

    Brockmann and I

    D. Brockmann and I. Sokolov, Chemical Physics284, 409 (2002)

  39. [39]

    Sibuya, Annals of the Institute of Statistical Mathe- matics31, 373 (1979)

    M. Sibuya, Annals of the Institute of Statistical Mathe- matics31, 373 (1979)

  40. [40]

    Sibuya and R

    M. Sibuya and R. Shimizu, Annals of the Institute of Statistical Mathematics33, 177 (1981)

  41. [41]

    B. I. Henry, T. Langlands, and P. Straka, Phys. Rev. Lett.105, 170602 (2010)

  42. [42]

    Metzler and J

    R. Metzler and J. Klafter, Physics Reports339, 1 (2000)

  43. [43]

    A. M. Mathai, R. K. Saxena, and H. J. Haubold,The H-function: theory and applications(Springer Science & Business Media, 2009)

  44. [44]

    Zumofen and J

    G. Zumofen and J. Klafter, Physical Review E51, 2805 (1995)

  45. [45]

    A. V. Chechkin, R. Metzler, V. Y. Gonchar, J. Klafter, and L. V. Tanatarov, Journal of Physics A: Mathematical and General36, L537 (2003)

  46. [46]

    V. V. Palyulin, G. Blackburn, M. A. Lomholt, N. W. Watkins, R. Metzler, R. Klages, and A. V. Chechkin, New Journal of Physics21, 103028 (2019)

  47. [47]

    Zaburdaev, S

    V. Zaburdaev, S. Denisov, and J. Klafter, Reviews of Modern Physics87, 483 (2015)

  48. [48]

    Hofert, Computational Statistics & Data Analysis55, 57 (2011)

    M. Hofert, Computational Statistics & Data Analysis55, 57 (2011)