pith. sign in

arxiv: 2312.14883 · v3 · submitted 2023-12-22 · 🧮 math.PR

Roots of polynomials under repeated differentiation and repeated applications of fractional differential operators

Pith reviewed 2026-05-24 05:00 UTC · model grok-4.3

classification 🧮 math.PR
keywords random polynomialsroot distributionsfractional differential operatorslimiting measurestransport mapslogarithmic potentialfree probability
0
0 comments X

The pith

The limiting root distribution after repeated applications of a fractional differential operator is the push-forward of the initial distribution under a transport map from PDE characteristics.

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

The paper starts from random polynomials of degree N with independent coefficients that converge in root distribution to some measure as N grows. It applies an operator of the form z^a (d/dz)^b a total of roughly Nt times and derives the new limiting root measure μ_t. This measure equals the image of the initial limiting measure under an explicit transport map T_t. The map is constructed by solving a PDE for the logarithmic potential of μ_t and integrating its characteristic curves. The construction recovers radial inward motion at constant speed for pure repeated differentiation and supplies a free-probability multiplication rule for general parameters.

Core claim

Starting from a random polynomial P^N of degree N with independent coefficients that possesses a limiting root distribution ν, the polynomial obtained after roughly Nt applications of z^a (d/dz)^b has limiting root distribution μ_t equal to the push-forward of ν under the map T_t, where T_t is the flow along the characteristic curves of the PDE satisfied by the logarithmic potential of μ_t.

What carries the argument

The transport map T_t obtained by flowing along the characteristic curves of the PDE satisfied by the logarithmic potential of the evolving measure.

If this is right

  • For repeated differentiation the roots move radially inward at constant speed until they reach the origin and disappear.
  • The transport map admits an interpretation in free probability as multiplication of an R-diagonal operator by an R-diagonal transport operator.
  • The construction supplies a push-forward characterization of the free self-convolution semigroup of radial measures on the complex plane.
  • When the operator involves integration the root dynamics become more complicated than simple transport.

Where Pith is reading between the lines

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

  • The same characteristic-flow method could be applied to other families of linear operators on polynomials whose action on the log-potential yields a closed PDE.
  • The free-probability multiplication rule may connect the root evolution to the multiplicative free convolution of circular elements.
  • Numerical checks on finite-N polynomials with Gaussian coefficients could confirm the predicted speed of radial motion before the large-N limit is taken.

Load-bearing premise

The initial random polynomial of degree N with independent coefficients possesses a limiting root distribution as N tends to infinity.

What would settle it

For a concrete initial distribution such as the circular law, compute the empirical roots of large-N polynomials after exactly floor(Nt) applications of the operator and test whether their empirical measure converges to the predicted push-forward.

Figures

Figures reproduced from arXiv: 2312.14883 by Brian C. Hall, Ching-Wei Ho, Jonas Jalowy, Zakhar Kabluchko.

Figure 1
Figure 1. Figure 1: The smaller roots (left) travel radially at constant speed until they hit the origin and die before time t. The larger roots (right) travel radially at constant speed without hitting the origin. The blue dots show the roots of all the polynomials with time s < t, while the red dots show the roots with time t. Shown for t = 0.4 starting from a Weyl polynomial with N = 300. transform of µt, as in (1.2). Sinc… view at source ↗
Figure 2
Figure 2. Figure 2: The degree-increasing case with a = 5/2 and b = 3/2. The roots move radially inward without reaching the origin. Shown for N = 150 and t = 1/5, starting from a Weyl polynomial. ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● … view at source ↗
Figure 3
Figure 3. Figure 3: The repeated integration case (a = 0, b = −1) starting from a Weyl polynomial. The blue dots show the roots of all the polynomials with time s < t, while the red dots show the roots with time t. Shown for t = 0.15 and N = 300 [PITH_FULL_IMAGE:figures/full_fig_p008_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: The roots of QN t (left) and RN t (right) with a = 0 and b = 0.999, starting from a Weyl polynomial. Shown with N = 1000 and t = 1/2. of the original polynomial, together with N t roots at the origin. By contrast, RN t has a positive fraction of its roots concentrated near a circle of some positive radius r. (The roots having magnitude greater than r are almost the same for the two polynomials.) Let us now… view at source ↗
Figure 5
Figure 5. Figure 5: The exponential profile (solid) and its concave majo￾rant (dashed) in the Weyl case, for t = 0.15 Thus, the limiting root distribution of QN t will be concentrated entirely on the ring of radius e −Mt = (1 + t) 1+t t t = (1 + t)  1 + 1 t t . When t is large, e −Mt ≈ (1 + t)e. 8.2. Singular behavior: the Weyl polynomial case. As our next example, we consider the Weyl polynomials, which correspond to g0(α)… view at source ↗
Figure 6
Figure 6. Figure 6: The roots of the original exponential polynomial, di￾lated by a factor of 1 + t (left), and the roots of P N t (right). The inner circle has radius t. Shown for t = 0.3 and N = 1, 000. the exponential polynomials is easily computed using Theorem 3.4 and µ0 assigns mass r to the disk of radius r, for all 0 ≤ r ≤ 1. In that case, Example 3.11 with a = 0 and b = −1 applies. The limiting root distribution µt o… view at source ↗
read the original abstract

We start with a random polynomial $P^{N}(z)$ of degree $N$ with independent coefficients. We then consider a new polynomial $P_{t}^{N}$ obtained by $\lceil Nt\rceil$ applications of a fractional differential operator of the form $z^{a} (d/dz)^{b},$ where $a$ and $b$ are real numbers. When $b>0,$ we compute the limiting root distribution $\mu_{t}$ of $P_{t}^{N}$ as $N\rightarrow\infty.$ We show that $\mu_{t}$ is the push-forward of the limiting root distribution of $P^{N}$ under a transport map $T_{t}$. The map $T_{t}$ is defined by flowing along the characteristic curves of a PDE satisfied by the log potential of $\mu_{t}.$ In the special case of repeated differentiation, our results may be interpreted as saying that the roots evolve radially \textit{with constant speed} until they hit the origin, at which point, they cease to exist. For general $a$ and $b,$ the transport map $T_{t}$ has a free probability interpretation as multiplication of an $R$-diagonal operator by an $R$-diagonal \textquotedblleft transport operator.\textquotedblright As an application, we obtain a push-forward characterization of the free self-convolution semigroup $\oplus$ of radial measures on $\mathbb{C}$. We also consider the case $b<0,$ which includes the case of repeated integration. More complicated behavior of the roots can occur in this case.

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 paper starts with a random polynomial P^N(z) of degree N with independent coefficients and considers the polynomial obtained after applying the fractional differential operator z^a (d/dz)^b a total of ⌈Nt⌉ times. It claims to compute the limiting root distribution μ_t as N→∞, showing that μ_t is the push-forward of the initial limiting root distribution under a transport map T_t obtained by flowing along characteristics of a PDE satisfied by the log potential of μ_t. Special cases include repeated differentiation (roots move radially at constant speed until reaching the origin) and a free-probability interpretation as multiplication by an R-diagonal transport operator; an application yields a push-forward characterization of the free self-convolution semigroup of radial measures on ℂ. The case b<0 (including repeated integration) is also treated.

Significance. If the claims hold under appropriate conditions, the work supplies a dynamical description of root evolution under fractional operators that links random polynomials to free probability and gives an explicit characterization of the radial free self-convolution semigroup. The transport-map construction via PDE characteristics is a potentially useful technical contribution when the initial limiting measure exists.

major comments (2)
  1. [Abstract] Abstract and opening paragraphs: the central claim that μ_t is computed for general independent coefficients presupposes the existence of a limiting root distribution μ_0 for the initial P^N. No moment or distributional hypotheses are stated, yet existence of such limits is a non-trivial prerequisite (typically requiring e.g. finite log-moments or Gaussianity) that is not proved or even explicitly assumed in the manuscript; without it the push-forward construction cannot be applied.
  2. [Abstract] The derivation that the evolved measure is exactly the push-forward under the characteristic flow T_t of the log-potential PDE is presented as the main result, but the manuscript supplies no error estimates, tightness arguments, or verification that the limiting empirical measure converges to this transported measure; the soundness of the transport step therefore cannot be assessed from the given text.
minor comments (1)
  1. The notation for the fractional operator z^a (d/dz)^b should include an explicit definition or reference for non-integer exponents a and b.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments on the assumptions and rigor. We address the two major comments point by point below.

read point-by-point responses
  1. Referee: [Abstract] Abstract and opening paragraphs: the central claim that μ_t is computed for general independent coefficients presupposes the existence of a limiting root distribution μ_0 for the initial P^N. No moment or distributional hypotheses are stated, yet existence of such limits is a non-trivial prerequisite (typically requiring e.g. finite log-moments or Gaussianity) that is not proved or even explicitly assumed in the manuscript; without it the push-forward construction cannot be applied.

    Authors: We agree that existence of the initial limiting measure μ_0 is a prerequisite not proved in the paper. The abstract and main text presuppose this limit exists when they refer to computing μ_t as its push-forward under T_t. The manuscript's focus is the evolution step assuming μ_0 is given, rather than proving existence (which is known in the literature under conditions such as i.i.d. Gaussian coefficients or finite log-moments). In revision we will explicitly state the assumption in the abstract and introduction, with a brief reference to standard sufficient conditions from the random-polynomial literature. revision: yes

  2. Referee: [Abstract] The derivation that the evolved measure is exactly the push-forward under the characteristic flow T_t of the log-potential PDE is presented as the main result, but the manuscript supplies no error estimates, tightness arguments, or verification that the limiting empirical measure converges to this transported measure; the soundness of the transport step therefore cannot be assessed from the given text.

    Authors: The derivation obtains the explicit form of μ_t by evolving the log-potential, deriving the associated PDE, and solving along characteristics to produce the transport map T_t; the limiting distribution is then defined as the push-forward of μ_0. We acknowledge that the manuscript contains no quantitative error bounds, tightness proofs, or direct verification that the empirical measures converge to this transported limit. The result is therefore formal in the sense that it describes the candidate limit under the standing assumption that the initial empirical measures converge to μ_0 and that the transport applies in the limit. In revision we will add a clarifying remark in the introduction and main theorem statement to this effect, and we will indicate the additional analytic work that would be needed for a fully rigorous convergence proof. revision: yes

Circularity Check

0 steps flagged

No circularity; forward evolution from assumed initial limit

full rationale

The paper assumes existence of the initial limiting root measure μ_0 for P^N with independent coefficients (a standard external prerequisite in random polynomial theory) and derives the evolved measure μ_t as its push-forward under the characteristic flow of the log-potential PDE. No equation reduces the transport map T_t or μ_t to a fitted quantity or self-definition; the construction is explicit from the PDE and free-probability interpretation. No self-citations are invoked as load-bearing uniqueness results. The derivation chain is therefore self-contained once the initial convergence (independent of the present work) is granted.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on the existence of a limiting root distribution for the initial polynomial and on the well-posedness of the characteristic flow for the log-potential PDE; no free parameters or invented entities are visible in the abstract.

axioms (2)
  • domain assumption A random polynomial of degree N with independent coefficients possesses a limiting empirical root measure as N tends to infinity.
    The paper begins with this measure and evolves it; the assumption is invoked at the outset of the abstract.
  • domain assumption The fractional differential operator z^a (d/dz)^b generates a well-defined evolution of the root measure via the stated PDE.
    The transport map is defined by flowing along characteristics of that PDE.

pith-pipeline@v0.9.0 · 5836 in / 1331 out tokens · 21887 ms · 2026-05-24T05:00:04.974691+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

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

  1. [1]

    Alazard, O

    T. Alazard, O. Lazar, Q.-H. Nguyen, On the dynamics of the roots of polynomials under differentiation, arXiv:2104.06921 [math.AP]

  2. [2]

    Angst, D

    J. Angst, D. Malicet, G. Poly, Almost sure behavior of the critical points of random polyno- mials, arXiv:2301.06973 [math.PR]

  3. [3]

    Arizmendi, J

    O. Arizmendi, J. Garza-Vargas, and D. Perales, Finite free cumulants: Multiplicative convo- lutions, genus expansion and infinitesimal distributions, arXiv:2108.08489 [math.CO]

  4. [4]

    S. T. Belinschi and H. Bercovici, Partially defined semigroups relative to multiplicative free convolution, Int. Math. Res. Not. (2005), 65–101

  5. [5]

    E. W. Benz, ¨Uber lineare verschiebungstreue Funktionaloperatoren und die Nullstellen ganzer Funktionen, Comment. Math. Helv. 7 (1934), 243–289

  6. [6]

    Bercovici and D

    H. Bercovici and D. Voicolescu, L´ evy–Hinˇ cin type theorems for multiplicative and additive free convolution, Pacific J. Math. 153 (1992), 217–248

  7. [7]

    Bercovici and D

    H. Bercovici and D. Voicolescu, Free convolution of measures with unbounded support, In- diana Univ. Math. J. 42 (1993), 733–773

  8. [8]

    Bercovici and D

    H. Bercovici and D. Voicolescu, Superconvergence to the central limit and failure of the Cramer theorem for free random variables, Probab. Theory Related Fields 103 (1995), 215– 222

  9. [9]

    Bøgvad, C

    R. Bøgvad, C. H¨ agg, B. Shapiro, Rodrigues’ descendants of a polynomial and Boutroux curves, Constr. Approx. 59 (2024), 737–798

  10. [10]

    S.-S. Byun, J. Lee, and T. R. Reddy. Zeros of random polynomials and their higher derivatives, Trans. Amer. Math. Soc., 375 (2022), 6311–6335

  11. [11]

    Campbell, S

    A. Campbell, S. O’Rourke, and D. Renfrew, The fractional free convolution of R-diagonal op- erators and random polynomials under repeated differentiation, arXiv:2307.11935 [math.PR]

  12. [12]

    Castin, Z

    Y. Castin, Z. Hadzibabic, S. Stock, J. Dalibard, and S. Stringari, Quantized vortices in the ideal Bose gas: a physical realization of random polynomials, Phys. Rev. Lett. 96 (2006), Article 040405

  13. [13]

    B. K. Driver, B. C. Hall, and T. Kemp, The Brown measure of the free multiplicative Brow- nian motion, Probab. Theory Related Fields 184 (2022), 209–273

  14. [14]

    Dennis and J

    M. Dennis and J. Hannay, Saddle points in the chaotic analytic function and Ginibre char- acteristic polynomial, J. Phys. A 36 (2003), 3379–3384

  15. [15]

    L. C. Evans, Partial differential equations. Second edition. Graduate Studies in Mathematics,

  16. [16]

    xxii+749 pp

    American Mathematical Society, Providence, RI, 2010. xxii+749 pp

  17. [17]

    Feng and D

    R. Feng and D. Yao, Zeros of repeated derivatives of random polynomials, Anal. PDE 12 (2019), 1489-1512

  18. [18]

    A. Galligo. Modeling complex root motion of real random polynomials under differentiation. In ISSAC 22—Proceedings of the 2022 International Symposium on Symbolic and Algebraic Computation, pp. 245–253. ACM, New York, 2022

  19. [19]

    B. C. Hall and C.-W. Ho, The Brown measure of a family of free multiplicative Brownian motions, Probab. Theory Related Fields 186 (2023), 1081–1166

  20. [20]

    B. C. Hall and C.-W. Ho, The heat flow conjecture for random matrices, preprint arXiv:2202.09660 [math.PR]

  21. [21]

    Haagerup and F

    U. Haagerup and F. Larsen, Brown’s spectral distribution measure for R-diagonal elements in finite von Neumann algebras, J. Funct. Anal. 176 (2000), 331–367

  22. [22]

    B. C. Hall, C.-W. Ho, J. Jalowy, and Z. Kabluchko, The heat flow, GAF, and SL(2;R), Indiana Univ. Math. J. , to appear, arXiv:2304.06665 [math.PR]. ROOTS OF POLYNOMIALS UNDER REPEATED DIFFERENTIATION 47

  23. [23]

    B. C. Hall, C.-W. Ho, J. Jalowy, and Z. Kabluchko, Zeros of random polynomials undergoing the heat flow, arXiv:2308.11685 [math.PR]

  24. [24]

    B. Hanin. Correlations and pairing between zeros and critical points of Gaussian random polynomials, Int. Math. Res. Not. IMRN (2015), 381–421

  25. [25]

    B. Hanin. Pairing of zeros and critical points for random polynomials, Ann. Inst. H. Poincar´ e Probab. Statist., 53 (2017)1498–1511

  26. [26]

    Hoskins and Z

    J. Hoskins and Z. Kabluchko, Dynamics of zeroes under repeated differentiation, Exp. Math. 32 (2023), 573–599

  27. [27]

    Kabluchko, Critical points of random polynomials with independent identically distributed roots, Proc

    Z. Kabluchko, Critical points of random polynomials with independent identically distributed roots, Proc. Amer. Math. Soc. , 143 (2015), 695–702

  28. [28]

    Kabluchko, Repeated differentiation and free unitary Poisson process, arXiv:2112.14729 [math.PR]

    Z. Kabluchko, Repeated differentiation and free unitary Poisson process, arXiv:2112.14729 [math.PR]

  29. [29]

    Kabluchko, Lee-Yang zeroes of the Curie-Weiss ferromagnet, unitary Hermite polynomials, and the backward heat flow, arXiv:2203.05533 [math.PR]

    Z. Kabluchko, Lee-Yang zeroes of the Curie-Weiss ferromagnet, unitary Hermite polynomials, and the backward heat flow, arXiv:2203.05533 [math.PR]

  30. [30]

    Kabluchko and H

    Z. Kabluchko and H. Seidel, Distances between zeroes and critical points for random poly- nomials with i.i.d. zeroes, Electron. J. Probab., 24 (2019), Paper No. 34

  31. [31]

    Kabluchko and D

    Z. Kabluchko and D. Zaporozhets, Asymptotic distribution of complex zeros of random an- alytic functions, Ann. Probab. 42 (2014), 1374–1395

  32. [32]

    A. A. Kilbas, H. M. Srivastava and J. J. Trujillo, Theory and applications of fractional differential equations (Vol. 204). elsevier. (2006)

  33. [33]

    Kiselev and C

    A. Kiselev and C. Tan, The flow of polynomial roots under differentiation, Ann. PDE 8 (2022), Paper No. 16, 69 pp

  34. [34]

    Limiting Spectral Distributions of Sums of Products of Non-Hermitian Random Matrices

    H. K¨ osters and A. Tikhomirov, Limiting spectral distributions of sums of products of non- Hermitian random matrices, arXiv:1506.04436

  35. [35]

    Marcus, Polynomial convolutions and (finite) free probability, arXiv:2108.07054 [math.CO]

    A. Marcus, Polynomial convolutions and (finite) free probability, arXiv:2108.07054 [math.CO]

  36. [36]

    Marcus, D

    A. Marcus, D. A. Spielman, and N. Srivastava, Finite free convolutions of polynomials, Probab. Theory Related Fields 182 (2022), 807–848

  37. [37]

    Michelen and X.-T

    M. Michelen and X.-T. Vu., Zeros of a growing number of derivatives of random polynomials with independent roots, arXiv:2212.11867 [math.PR]

  38. [38]

    Mingo and R

    J. Mingo and R. Speicher, Free probability and random matrices, Springer, 2016

  39. [39]

    Nica and R

    A. Nica and R. Speicher. On the multiplication of free N-tuples of noncommutative random variables, Amer. J. Math. , 118 (1996), 799–837

  40. [40]

    Nica and R

    A. Nica and R. Speicher. Lectures on the combinatorics of free probability, London Math. Soc. Lecture Note Ser., 335, Cambridge University Press (2006)

  41. [41]

    Olver, ed

    F. Olver, ed. NIST handbook of mathematical functions, Cambridge university press. (2010)

  42. [42]

    O’Rourke and S

    S. O’Rourke and S. Steinerberger, A nonlocal transport equation modeling complex roots of polynomials under differentiation, Proc. Amer. Math. Soc. 149 (2021), 1581–1592

  43. [43]

    O’Rourke, Critical points of random polynomials and characteristic polynomials of random matrices, Int

    S. O’Rourke, Critical points of random polynomials and characteristic polynomials of random matrices, Int. Math. Res. Not. 18 (2016), 5616–5651

  44. [44]

    O’Rourke and N

    S. O’Rourke and N. Williams, On the local pairing behavior of critical points and roots of random polynomials, Electron. J. Probab. 25(2020), Paper No. 100

  45. [45]

    Pemantle and I

    R. Pemantle and I. Rivin, The distribution of zeros of the derivative of a random polynomial, in: Advances in combinatorics (2013), 259–273, in part based on the 3rd Waterloo workshop on computer algebra (WWCA, W80), Waterloo, Canada, 2011. Berlin: Springer

  46. [46]

    Petz and J

    D. Petz and J. R´ effy, Large deviation for the empirical eigenvalue density of truncated Haar unitary matrices, Probab. Theory Related Fields, 133 (2005), 175–189

  47. [47]

    Shlyakhtenko and T

    D. Shlyakhtenko and T. Tao, Fractional free convolution powers, Indiana Univ. Math. J. 71(2022), 2551–2594

  48. [48]

    Steinerberger, A nonlocal transport equation describing roots of polynomials under differ- entiation, Proc

    S. Steinerberger, A nonlocal transport equation describing roots of polynomials under differ- entiation, Proc. Amer. Math. Soc. 147 (2019), 4733–4744

  49. [49]

    S. D. Subramanian, On the distribution of critical points of a polynomial, Electron. Commun. Probab. 17 (2012), Paper No. 37

  50. [50]

    Voiculescu, Multiplication of certain noncommuting random variables, J

    D. Voiculescu, Multiplication of certain noncommuting random variables, J. Operator Theory 18 (1987), 223–235

  51. [51]

    Bercovici and D

    H. Bercovici and D. Voiculescu, Free convolution of measures with unbounded support, In- diana Univ. Math. J. 42 (1993), 733–773. 48 BRIAN C. HALL, CHING-WEI HO, JONAS JALOWY, AND ZAKHAR KABLUCHKO

  52. [52]

    ˙Zyczkowski and H-J

    K. ˙Zyczkowski and H-J. Sommers, Truncations of random unitary matrices, J. Phys. A , 33 (2000), 2045–2057. Brian C. Hall: University of Notre Dame, Notre Dame, IN 46556, USA Email address: bhall@nd.edu Ching-Wei Ho: Institute of Mathematics, Academia Sinica, Taipei 10617, Taiwan Email address: chwho@gate.sinica.edu.tw Jonas Jalowy: Institut f¨ur Mathemat...