Pith. sign in

REVIEW 3 major objections 4 minor 22 references

Pointwise upper bound for the fundamental solution of fractional Fokker-Planck equation

T0 review · 3 major / 4 minor · reviewed 2026-08-10 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read The fractional Fokker-Planck kernel satisfies a pointwise polynomial decay bound for all derivatives.

desk verdict The main lemma is false as stated, a key case is not proved, and the theorem is already subsumed by Hou-Zhang's two-sided bounds; this is not ready for review. read the letter →

arxiv 2501.13771 v2 pith:P7ILDAFH submitted 2025-01-23 math.AP

Haina Li , Yiran Xu This is my paper · ORCID
classification math.AP MSC 35A0835Q84
keywords fractionalFokker-PlanckequationKolmogorovfundamentalsolutionpointwiseupperboundLittlewood-PaleydecompositionLaplaciankineticderivativeestimates
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

The paper sets out to prove a pointwise upper bound for the fundamental solution --- the kernel that propagates a point source --- of the one-dimensional fractional kinetic Fokker-Planck equation $$\partial_t f + v\partial_x f + |D_v|^{2s} f = 0, \qquad 0

What carries the argument

The load-bearing object is the phase $$M(\xi,\eta)=\frac{1}{2s+1}\frac{|\xi|^{2s}\xi-|\eta|^{2s}\eta}{\xi-\eta}.$$ Lemma 3.1 bounds its derivatives; Lemma 3.2 transfers those bounds to $e^{-M}$ through the higher-chain-rule formula for derivatives of composed functions; Lemma 3.3 localizes the Fourier integral with dyadic cutoffs $\chi(|\xi|/2^{m_1})\chi(|\eta|/2^{m_2})$ and integrates by parts in $\eta$ --- twice for $0<s\le 1/2$ and three times for $1/2<s<1$ --- to convert powers of $1/v$ into decay. The three dyadic regimes $|m_1-m_2|\le 2$, $m_2\ge m_1+3$, and $m_2\le m_1-3$ are then summed separately. The Fourier representation evaluates the inverse transform at $(x,-x-v)$, which is why the final weight is expressed through $x$ and $x+v$.

What would settle it

Evaluate Lemma 3.2 at $\eta=0$ for $(m_1,m_2)=(1,0)$: compute $\partial_\xi e^{-|\xi|^{2s}/(2s+1)}$ and compare both sides over $\xi>0$; if no uniform constant works, the dyadic sums in Lemma 3.3 cannot be uniform. Separately, computing $|K(1,x,0)|$ for large $x$ at several values of $s$ would test the claimed $\langle x\rangle^{-(2+2s)}$ decay without relying on the omitted $|x|\ge |v|$ calculation.

Watch

Extended reading notes

Core claim

On its own terms, the central claim is Theorem 1.1: for $K$ the fundamental solution of (1.1), for every $b_1,b_2\in\mathbb{N}$ and every $0<s<1$ there is an arbitrarily small $\varepsilon>0$ with $$|\$partial_x^{{b_1}}$\$partial_v^{{b_2}}$K(1,x,v)| \lesssim \frac{1}{\langle x,x+v\$rangle^{{2+2s-2\varepsilon}}$\langle x\$rangle^{{\varepsilon+b_1}}$\langle x+v\$rangle^{{\varepsilon+b_2}}$},$$ where $\langle u,w\rangle=(1+|u|^2+|w|^2)^{1/2}$. The proof derives $K$ by space-velocity Fourier transform and characteristics, obtaining the phase $M(\xi,\eta)=\int_0^1 |(\xi-\eta)\tau+\eta|^{2s}\,d\tau$, then proves derivative bounds for $M$ and for $e^{-M}$ and sums Littlewood-Paley blocks after integration by parts in $v$. The fully written computation covers $0<s\le 1/2$ in the regime $|v|\ge |x|$; the regime $|x|\ge |v|$ is stated to be similar, and the case $1/2<s<1$ is summarized as analogous.

Load-bearing premise

The whole dyadic summation rests on Lemma 3.2's pointwise control of derivatives of $e^{-M}$, and that displayed bound is not homogeneous for pure $\xi$-derivatives: at $\eta=0$ and order $(1,0)$ the left side scales like $|\xi|^{2s-1}$ while the stated right-hand side scales like $|\xi|^{2s-2}$, so the block estimates in Lemma 3.3 inherit whatever defect that bound has.

Editorial extensions

If this is right

  • Because the combined weight $\langle x,x+v\rangle^{-(2+2s)}$ is integrable over $\mathbb{R}^2$, the bound puts $K$ and all its derivatives into natural weighted $L^1$ spaces.
  • The scaling identity $K(t,x,v)=t^{-1-1/s}K(1,x/t^{1+1/(2s)},v/t^{1/(2s)})$ upgrades the time-one estimate to an explicit $t$-dependent decay bound at every time.
  • Since $b_1$ and $b_2$ are arbitrary, the kernel is smooth with controlled tails in both $x$ and $v$, not merely integrable.
  • The bound has the expected heat-kernel shape for fractional Kolmogorov equations, which is the upper estimate the authors point to for future well-posedness and regularity arguments.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • If Lemma 3.2 survives a direct homogeneity check, the dyadic summation is the only place it enters; a short calculation at $\eta=0$ would settle whether the displayed exponents are harmless or need repair.
  • The omitted $|x|\ge |v|$ case is likely reachable by the same scheme with integration by parts in $\xi$ instead of $\eta$, since the final representation shows the decay must be split between $x$ and $x+v$.
  • A two-sided version would need lower bounds for $e^{-M}$ on the same dyadic blocks, so the phase formula turns sharpness into a convexity question; the paper proves only the upper side.
  • The same frequency-block framework should adapt to higher-dimensional kinetic equations by replacing $M$ with the vector analogue of the same time integral and using anisotropic dyadic blocks.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 4 minor

Summary. The paper studies the fundamental solution of the one-dimensional fractional Fokker-Planck equation ∂tf + v∂xf + |Dv|^{2s}f = 0. The fundamental solution is written as an explicit oscillatory integral with phase M(ξ,η)=∫_0^1 |(ξ−η)τ+η|^{2s}dτ. The authors prove derivative estimates for M and e^{-M} (Lemmas 3.1 and 3.2), then via Littlewood-Paley dyadic decomposition and several integrations by parts (Lemma 3.3) they bound each dyadic piece and sum the pieces (Section 4) to obtain the pointwise upper bound in Theorem 1.1. The result is claimed for all s∈(0,1) and all derivative orders.

Significance. If the pointwise bound is correct, it gives a sharp (up to ε) polynomial decay in the kinetic variables for the fundamental solution and its derivatives, with a transparent Fourier/Littlewood-Paley proof. The paper uses no fitting parameters, and the statement is consistent with the known two-sided heat-kernel estimates of [9]. However, the proof as written contains a false key lemma and omits essential cases, so the result is not established by this manuscript.

major comments (3)
  1. [Lemma 3.2, Eq. (3.2)] The estimate (3.2) is false. Take s=1/2, m1=2, m2=0, η=0, ξ=8. Then M(ξ,0)=|ξ|/2 and ∂_ξ^2 e^{-M} = (1/4)e^{-|ξ|/2}, so the left side is (1/4)e^{-4} ≈ 0.00458. The right side of (3.2) is |ξ|^{1+2s-m1}/(|ξ|^{m1+1}+1) = 1/8^3 = 1/512 ≈ 0.00195, so the inequality fails. Since Lemma 3.3 invokes (3.2) in every dyadic block, the dyadic summation in Section 4 is unjustified. The source of the failure is visible in Lemma 3.1, whose final '∼' display drops the |ξ|^{2s-m1-m2} term present in the subcase |ξ|≥4|η|.
  2. [Section 4, case 3 (|x|≥|v|)] The entire case |x|≥|v| is dismissed with the sentence 'Similar to case 2' and the asserted bound ∑_{m1,m2} ~K_{m1,m2} ≲ |x|^{-(2+2s+b1-ε)} |v|^{-(b2+ε)} is not derived. Because the theorem's final bound requires this regime, the proof is incomplete.
  3. [Section 4, case 2.2 (1/2<s<1)] The case 1/2<s<1 is only summarized: the text says 'Given its similarity to case 2.1, we just summarize key findings here' and 'we omit the details here'. The three displayed estimates for H_{m1,m2} are stated without derivation, and the summation to the claimed bounds is not shown. This is a separate case in the main theorem, so the omission is load-bearing.
minor comments (4)
  1. [Lemma 3.3 statement] The integers n1,...,n10 appear in the estimates before being defined; they should be quantified in the statement of the lemma, since the proof later chooses them.
  2. [Theorem 1.1 statement] The phrase 'there exists ε > 0 ... where ε arbitrarily small' is ambiguous; the theorem should state that for every sufficiently small ε>0 the bound holds with a constant depending on ε.
  3. [End of proof, Section 4] The notation ⟨x, x+v⟩ is not defined, and the equality |∂_x^{b1} ∂_v^{b2} K(1,x,v)| = |∂_x^{b1} ∂_v^{b2} K(1,x,−x−v)| in the final display is not justified from the Fourier integral representation.
  4. [Throughout] There are several typos, e.g., 'stvarepsilons' in the proof of Lemma 3.3 and 'FUNDAMENT AL' in the title header of the arXiv source; a careful copyedit is needed.

Circularity Check

1 steps flagged · score 4.0 of 10

Lemma 3.2's proof invokes Lemma 3.2 itself in the |ξ| ≤ |η|/4 subcase; the Fourier/dyadic derivation is otherwise self-contained.

  1. other [Section 3, proof of Lemma 3.2, subcase |ξ| ≤ 1/4 |η|]
    "For the case |ξ| ≤ 1/4 |η|, we apply Lemma 3.2 and (3.3) once more to achieve"

    Lemma 3.2 is the estimate being proved. In this subcase the proof feeds Lemma 3.2 itself into the Faà di Bruno expansion (3.3), so the displayed bound for ∂ξ^m1 ∂η^m2 e^{-M} is assumed rather than derived on that branch. Lemma 3.3 applies Lemma 3.2 in every dyadic block and Theorem 1.1 relies on Lemma 3.3, making the self-reference load-bearing in the written derivation. If the intended reference was Lemma 3.1, the loop would be a typo, but as written the proof is circular here.

full rationale

The main result is not fitted: K(1,x,v) is written as an explicit Fourier integral of e^{-M}, and the dyadic estimates and summations in Lemmas 3.3 and Section 4 are independent computations. No parameter is calibrated to the desired bound, and there is no load-bearing citation to prior work by the same authors; Remark 1.3 explicitly disclaims dependence on the concurrent paper [9]. The only circular step is internal: the proof of Lemma 3.2 says 'we apply Lemma 3.2 and (3.3) once more' in the subcase |ξ| ≤ |η|/4, i.e. it assumes the very estimate being proved. Since Lemma 3.3 invokes Lemma 3.2 block-by-block and Theorem 1.1 depends on Lemma 3.3, this self-reference is load-bearing as written. The skeptic's scaling counterexample is a correctness defect, not a circularity defect: a false estimate is not the same as an estimate used as its own input. The self-reported omissions (case 2.2 'we omit the details here', case 3 'Further details are omitted here', Remark 4.1) are completeness gaps rather than circularity. Score 4 reflects one load-bearing self-referential line with otherwise independent content.

Assumptions & free parameters 0 free parameters · 4 assumptions · 0 invented entities

No free parameters appear, since the proof derives the bound from the explicit symbol without fitting exponents. No new entities are introduced. The axioms are standard analytic tools, plus an implicit regularity assumption on oscillatory integrals; the main risk is not the ledger but the correctness of the derivative estimates.

assumptions (4)
  • domain assumption The kernel K(1,x,v) can be written as an absolutely convergent Fourier integral and repeated integration by parts produces no boundary terms.
    Section 2 defines K as the inverse Fourier transform of exp(-M); the proof in Section 4 integrates by parts in η or ξ without justifying decay at infinity.
  • standard math Littlewood-Paley dyadic blocks can be summed after integration by parts.
    Lemma 3.3 estimates individual blocks χ(|ξ|/2^{m1})χ(|η|/2^{m2}) and then sums over all m1,m2∈Z; the paper relies on standard convergence of the decomposition.
  • standard math The fundamental solution satisfies the stated scaling relation K(t,x,v)=t^{-1-1/s}K(1,x/t^{1+1/(2s)},v/t^{1/(2s)}).
    Used to obtain Corollary 1.2; it follows from homogeneity of the symbol but is not proved in detail.
  • domain assumption The non-smooth power |η|^{2s} and the function e^{-M} are regular enough for the derivative bounds in Lemmas 3.1 and 3.2 for all s∈(0,1), including s<1/2 where |η|^{2s} is not C^1 at the origin.
    The proof applies Faà di Bruno and differentiates |η|^{2s} as if it were smooth; singularities at η=0 are not discussed.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Pointwise upper bound for the fundamental solution of fractional Fokker-Planck equation." pith.science (2026). https://pith.science/paper/P7ILDAFH

@misc{pith2026250113771,
  author       = {Pith},
  title        = {Pith review of: Pointwise upper bound for the fundamental solution of fractional Fokker-Planck equation},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/P7ILDAFH}},
  note         = {Machine review of arXiv:2501.13771}
}
read the original abstract

In this paper, we investigate the fundamental solution of the fractional Fokker-Planck equation. Utilizing the Littlewood-Paley decomposition technology, we present a concise proof of the pointwise estimate for the fundamental solution.

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

22 extracted references · 16 canonical work pages

  1. [9]

    https://doi.org/10.48550/arXiv.2410.18614

    Haojie Hou and Xicheng Zhang, Heat kernel estimates for n onlocal kinetic operators, arxiv preprint arxiv:2410.186 14, (2024). https://doi.org/10.48550/arXiv.2410.18614

    arXiv doi:10.48550/arxiv.2410.18614 2024

  2. [1]

    D. G. Aronson, Bounds for the fundamental solution of a pa rabolic equation, Bull. Amer. Math. Soc. , 73 (1967), 890–896

    work page 1967

  3. [2]

    Auscher, C

    P. Auscher, C. Imbert and L. Niebel, W eak solutions to Kol mogorov-Fokker-Planck equations: regularity, existence and uniqueness, arxiv preprint arXiv:2403.17464, (2024). htt ps://doi.org/10.48550/arXiv.2403.17464

    work page doi:10.48550/arxiv.2403.17464 2024

  4. [3]

    R. M. Blumenthal and R. K. Getoor, Some Theorems on Stable Processes, Trans. Amer. Math. Soc. , 95 (1960), 263–273. https://doi.org/10.2307/1993291

    work page doi:10.2307/1993291 1960

  5. [4]

    E. B. Davies, Explicit Constants for Gaussian Upper Boun ds on Heat Kernels, American Journal of Mathematics , 109 (1987), 319–333. https://doi.org/10.2307/2374577

    work page doi:10.2307/2374577 1987

  6. [5]

    E. B. Fabes, Gaussian upper bounds on fundamental soluti ons of parabolic equations; the method of nash. In: Dell’Ant onio, G., Mosco, U. (eds) Dirichlet Forms, Lecture Notes in Mathem atics, 1563. Springer, Berlin, Heidelberg, (1993). https://doi.org/10.1007/BFb0074089

    work page doi:10.1007/bfb0074089 1993

  7. [6]

    Farkas and L

    B. Farkas and L. Lorenzi, On a class of hypoelliptic opera tors with unbounded coefficients in RN , Communications on Pure and Applied Analysis , 8 (2009), 1159–1201. doi: 10.3934/cpaa.2009.8.1159

    work page doi:10.3934/cpaa.2009.8.1159 2009

  8. [7]

    https://doi.org/10.48550/arXiv.2302.04392

    Zimo Hao and Michael R¨ ockner and Xicheng Zhang, Second o rder fractional mean-field SDEs with singular kernels and measure initial data, arxiv preprint arxiv:2302.04392, (2 023). https://doi.org/10.48550/arXiv.2302.04392

    arXiv doi:10.48550/arxiv.2302.04392

Show all 22 references
  1. [8]

    H¨ ormander, Hypoelliptic second order differential e quations, Acta Math

    L. H¨ ormander, Hypoelliptic second order differential e quations, Acta Math. , 119 (1967), 147–171. https://doi.org/10.1007/BF02392081

    1967 doi

  2. [10]

    Imbert and L

    C. Imbert and L. Silvestre, The weak Harnack inequality for the Boltzmann equation without cut-off, Journal of the European Mathematical Society , 22 (2019), 507–592

    2019

  3. [11]

    Kolmogorov, Zuf¨ allige Bewegungen (zur Theorie der Brownschen Bewegung), Ann

    A. Kolmogorov, Zuf¨ allige Bewegungen (zur Theorie der Brownschen Bewegung), Ann. of Math.(Second Series) , 35 (1934), 116–117. https://doi.org/10.2307/1968123

    1934 doi

  4. [12]

    L. P. Kuptsov, Fundamental solutions of certain degene rate second-order parabolic equations, Mathematical Notes of the Academy of Sciences of the USSR , 31 (1982), 283–289 . https://doi.org/10.1007/BF01138938

    1982 doi

  5. [13]

    Lanconell and S

    E. Lanconell and S. Polidoro, On a class of hypoelliptic evolution operators, Rend. Sem. Mat. , 52 (1994), 29–63. https://www.researchgate.net/publication/243766432

    1994

  6. [14]

    Lanconelli and A

    A. Lanconelli and A. Pascucci, Nash Estimates and Upper Bounds for Non-homogeneous Kolmogorov Equations, Potential Anal, 47 (2017), 461–483 . https://doi.org/10.1007/s11118-017-9 622-1

    2017 doi

  7. [15]

    Loher, Semi-local behaviour of non-local hypoellip tic equations: divergence form, arxiv preprint arxiv:2404 .05612, (2024)

    A. Loher, Semi-local behaviour of non-local hypoellip tic equations: divergence form, arxiv preprint arxiv:2404 .05612, (2024). https://doi.org/10.48550/arXiv.2404.05612

    2024 doi arXiv

  8. [16]

    Nash, Continuity of Solutions of Parabolic and Ellip tic Equations, American Journal of Mathematics , 80 (1958), 931–954

    J. Nash, Continuity of Solutions of Parabolic and Ellip tic Equations, American Journal of Mathematics , 80 (1958), 931–954. https://doi.org/10.2307/2372841

    1958 doi

  9. [17]

    Niebel, Kinetic maximal Lp µ (Lp)-regularity for the fractional Kolmogorov equation with v ariable density, Nonlinear Analysis, 214 (2022), 112517

    L. Niebel, Kinetic maximal Lp µ (Lp)-regularity for the fractional Kolmogorov equation with v ariable density, Nonlinear Analysis, 214 (2022), 112517. https://doi.org/10.1016/j.na.2021.112 517

    2022 doi

  10. [18]

    Niebel and R

    L. Niebel and R. Zacher, Kinetic maximal L2-regularity for the (fractional) Kolmogorov equation, J. Evol. Equ. , 21 (2021), 3585–3612

    2021

  11. [19]

    Pascucci and S

    A. Pascucci and S. Polidoro, A Gaussian upper bound for t he fundamental solutions of a class of ultraparabolic equat ions, J. Math. Anal. Appl. , 282 (2003), 396–409. https://doi.org/10.1016/S0022-247X(0 3)00159-8

    2003 doi

  12. [20]

    Polidoro, On a class of ultraparabolic operators of K olmogorov-Fokker-Planck type, Le Matematiche , 49 (1994), 53–105

    S. Polidoro, On a class of ultraparabolic operators of K olmogorov-Fokker-Planck type, Le Matematiche , 49 (1994), 53–105

    1994

  13. [21]

    P´ olya, On the zeros of an integral function represen ted by Fourier’s integral, Messenger of Math , 52 (1923), 185–188

    G. P´ olya, On the zeros of an integral function represen ted by Fourier’s integral, Messenger of Math , 52 (1923), 185–188

    1923

  14. [22]

    W eber, The fundamental solution of a degenerate part ial differential equation of parabolic type, Transactions of the American Mathematical Society , 71 (1951), 24–37

    M. W eber, The fundamental solution of a degenerate part ial differential equation of parabolic type, Transactions of the American Mathematical Society , 71 (1951), 24–37. School of Mathematics and Statistics, Beijing Institute of Technology, Beijing 100081, China. Email addres...

    1951

Pith tools

Reviewed August 10, 2026 · model on record in the stance chip above.