Pith. sign in

REVIEW 3 major objections 3 minor 16 references

A novel advancement in the study of Appell polynomials via Pad\`e rational approximants

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

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

Pith's one-line read This paper claims that replacing the Appell amplitude A(t) by its Padé rational approximant, and then substituting the differential operator ∂_x for t in the operational identity a_n(x)=A(∂_x){x^n}, yields explicit polynomial approximations

desk verdict The Hermite identities in Theorems 1–3 are new and correct, but the Euler [2|1] example is invalid—that Padé approximant does not exist—and the paper needs repair before acceptance. read the letter →

arxiv 2509.03178 v1 pith:TGPNPVW2 submitted 2025-09-03 math.CA

Giuseppe Dattoli , Subuhi Khan , Ujair Ahmad This is my paper · ORCID
classification math.CA MSC 32E3041A2111B8333C4505A4041-04
keywords AppellpolynomialsPadéapproximantsHermitetruncatedexponentialChebyshevumbralcalculusmonomialityprincipleoperationalmethods
topics The Collatz Conjecture
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

This paper sets out to show that the whole family of Appell polynomials a_n(x), defined by a single generating amplitude A(t) through the generating function Σ t^n/n! a_n(x)=A(t)e^{xt}, can be approximated order by order by replacing A(t) with a Padé rational function. Because a_n(x) can be recovered as A(∂_x) acting on x^n, a rational approximant of A becomes a rational differential operator, and expanding its denominator yields explicit polynomial approximations of a_n. The paper demonstrates the recipe on Hermite polynomials, obtaining representations through truncated exponential polynomials and, at higher orders, through two-variable Chebyshev polynomials; it then applies the same mechanism to Euler and Bernoulli polynomials using umbral notation and sketches Bessel functions as a further example. If the approach is valid, it gives a general way to connect Appell sequences to other named special-polynomial families and to compute recurrences and differential equations for the approximants.

What carries the argument

The central object is the Padé approximant of the Appell amplitude, [r|s]A(t)=P_r(t)/Q_s(t), promoted to an operator by t→∂_x in the operational identity a_n(x)=A(∂_x){x^n}. Expanding Q_s(∂_x)^{-1} as a formal series converts the action on x^n into finite sums of truncated exponential polynomials e_n(x,y) or e_n^{(2)}(x,y); when Q_s is a quadratic polynomial, the expansion is instead carried by two-variable Chebyshev polynomials U_n(a,b). This substitution is what makes each identity in the paper follow from the Padé table of one function A(t).

What would settle it

Run the Padé matching equations (1.4) for A(t)=2/(e^t+1) with [2|1]. At the t^3 step the condition is c_3 + b_1 c_2 = a_3, i.e. 1/24 = 0, so no such rational approximant exists; Theorem 5's formula (3.14) therefore cannot be derived as claimed from a [2|1] approximant. This is a concrete check the reader can repeat by hand.

Watch

Extended reading notes

Core claim

The discovery is an operational substitution rule: if Σ t^n/n! a_n(x)=A(t)e^{xt}, then a_n(x)=A(∂_x)x^n. Writing the Padé approximant [r|s]A(t)=P_r(t)/Q_s(t) and replacing t by ∂_x gives [r|s]a_n(x)=P_r(∂_x)Q_s(∂_x)^{-1}x^n, an Appell-type approximation of a_n. Expanding the inverse denominator as a formal series turns this into combinations of truncated exponential polynomials; for Hermite amplitudes it yields, for instance, [1|1]He_n(x)=e_n^{(2)}(x,-1/4)-n(n-1)/4 e_{n-2}^{(2)}(x,-1/4), and for [0|2] and [3|2] the expansions are expressed through second-kind Chebyshev polynomials U_r(a,b). The paper further shows the approximated polynomials inherit quasi-monomial operators, so recurrences

Load-bearing premise

The load-bearing premise is that every Padé order [r|s] used for the amplitude A(t) actually has a solution, and the paper never proves this; in the Euler case A(t)=2/(e^t+1) the [2|1] approximant does not exist, so Theorem 5 rests on a nonexistent object.

Editorial extensions

If this is right

  • Every Appell polynomial whose amplitude admits a Padé approximant gets an explicit polynomial [r|s]a_n(x) that is still Appell and is expressible through truncated exponentials or Chebyshev polynomials.
  • Increasing the Padé order [r|s] improves the match with a_n(x), as illustrated for Hermite orders [1|1], [2|1], and [3|2].
  • Approximants inherit quasi-monomiality, so recurrences and differential equations follow automatically; Theorem 4 and Corollary 1 give the [1|1] Hermite example.
  • Umbral notation lets the same rational approximations carry over to Bernoulli, Genocchi, and Bessel-type functions and to second-order Appell families without solving new linear systems each time.

Reading between the lines

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

  • The paper's identities are stated order by order; a natural extension is to prove an existence criterion for the Padé system (1.4) per Appell amplitude, so the method can be applied safely to any A(t).
  • If the Hermite identities are taken at face value, they suggest a dictionary: Hermite approximants live in the truncated-exponential/Chebyshev hierarchy, so recurrences or generating functions of one family can be transported to rational approximants of another.
  • A testable next step, mentioned only as a prospect in the paper, is to apply the same rational-approximant substitution to Sheffer sequences, for instance Laguerre-type generators, and check whether their approximants again collapse to truncated-exponential sums.
  • For the Euler amplitude, a defensible repair would be to define a modified approximant by solving only the consistent matching equations and dropping the impossible t^3 condition; the resulting formula would differ from (3.14) and would need separate numerical testing.
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 / 3 minor

Summary. The paper proposes to approximate Appell polynomial sequences by replacing the amplitude A(t) in the generating function with Padé approximants and using the operational representation a_n(x)=A(∂_x)x^n. For Hermite polynomials, explicit identities are derived expressing the resulting polynomials in terms of truncated exponential polynomials and two-variable Chebyshev polynomials; monomiality is used to obtain multiplicative operators and differential equations. The final sections extend the formalism to Euler, Bernoulli, and Bessel families via umbral notation.

Significance. If correct, the construction gives a systematic way to obtain closed-form rational approximations of Appell sequences and reveals concrete connections to truncated exponential and Chebyshev polynomials. The algebraic identities in Theorems 1–3 are checkable by direct expansion and appear sound, and the umbral treatment in Theorems 6–7 is a useful computational device. However, the Euler example contains a nonexistent Padé approximant, so the central claim is not uniformly established and requires repair.

major comments (3)
  1. [§3, Theorem 5, Eq. (3.14)] The [2|1] Padé approximant of A(t)=2/(e^t+1) does not exist under the paper's own matching conditions (1.4). The Maclaurin coefficients are c0=1, c1=-1/2, c2=0, c3=1/24. For [2|1], the k=3 equation is c3 + b1 c2 = 0, i.e. 1/24 = 0, which is impossible. The rational function displayed in the proof matches only up to order t^2; its t^3 coefficient is 0, not 1/24. Thus Eq. (3.14) is not derived from a Padé approximant, and the claim 'third-order PA' is false.
  2. [§3, Eq. (3.14)] Even if the displayed rational operator is used as a formal approximation, the last term in (3.14) is misindexed. Applying ∂_x^2 to e_n(x,-1/12) gives n(n-1)e_{n-2}(x,-1/12), with e_{n-2} the same first-order truncated exponential polynomial of Eq. (2.5); it is not e^{(2)}_{n-2}(x,-1/12), the second-order polynomial of Eq. (2.15). The proof skips this step, and the two polynomial families are different, so the identity must be corrected.
  3. [Abstract and §1] The method's validity is stated as unconditional, but it depends on the solvability of the linear system (1.4). The Euler [2|1] example is a concrete counterexample. The authors should either restrict the claims to approximants whose existence is verified (as in Theorems 1–3 and 6) or include a criterion for existence. As written, the abstract's promise that the approach 'yields the possibility of determining the approximation' overstates what is proved.
minor comments (3)
  1. [§2, Theorems 1 and 3] The notation [1|1] and [3|2] is nonstandard: these rational functions are Padé approximants of e^{-z} with z=t^2, corresponding to [2|2] and [6|4] approximants in the original variable t. Please clarify that the order is taken in the squared variable, or use the standard [m|n] order in t, to avoid ambiguity.
  2. [§4, Eq. (4.12)] The summation index is inconsistent: the left side uses n, while the summand contains r! and U_r. It should be sum over r from 0 to infinity.
  3. [Throughout] There are several presentation issues: 'exponetial' in §4; the caption 'Figure 3' is repeated before the actual Figure 4; and the term 'second order' is used for different orders in different variables. A careful editing pass is needed.

Circularity Check

0 steps flagged · score 2.0 of 10

No significant circularity: central Hermite approximants are algebraic rewrites checked against exact polynomials; minor self-citations and a definitional umbral step do not make the derivation circular.

full rationale

The claimed derivation chain is not circular in the sense defined here. The central Hermite results (Theorems 1–3) are obtained by taking a rational function (e.g., [1|1]e^{-x} = (1-x/2)/(1+x/2)), substituting x -> ∂²/2, and expanding the resulting rational differential operator on x^n using the generating functions of truncated exponential or Chebyshev polynomials. No free parameter is fitted to the target Hermite polynomials; the identities are exact algebraic rewrites of the rational operator, and the figures compare them against the exact He_n(x) as external benchmarks. The operational representation (2.8) is an elementary consequence of the Appell generating function, so the citation to [10] is not load-bearing. The umbral extension in Section 3 is the paper's own formalism (refs [13], [16] have overlapping authorship), but the Bernoulli/Euler examples are secondary and are not used to fit anything; however, they appear definitional rather than derived from the Padé matching system (1.4). I also flag, as a correctness risk rather than circularity, that the claimed [2|1] Padé approximant in Theorem 5 does not satisfy the paper's own matching equations: for A(t)=2/(e^t+1), c0=1, c1=-1/2, c2=0, c3=1/24, so the k=3 condition c3+b1 c2=0 would require 1/24=0. Thus Theorem 5 and equation (3.14) lack a valid Padé derivation. Similarly, substituting x=t²/2 into [1|1]e^{-x} does not yield the [1|1] Padé of e^{-t²/2} under (1.4). These are mathematical defects, not circularity. Overall score 2 reflects minor reliance on the authors' own umbral references and a definitional umbral step; the central Hermite content is self-contained and externally checkable.

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

No free parameters are fitted; the Padé coefficients are determined by the matching conditions. The main axioms are standard special-function tools or prior results cited by the authors. The most fragile assumption is that formal umbral Padé substitution preserves the validity of standard Padé approximation, which fails in the Euler [2|1] example.

assumptions (6)
  • standard math Padé approximant matching conditions (1.3)-(1.4) require solving m+n+1 linear equations with b0=1
    Used throughout to derive [r|s] approximants of the amplitude A(t).
  • domain assumption Operational representation a_n(x)=A(∂_x){x^n}
    Cited from [10]; converts an amplitude approximation into a polynomial approximation by substituting ∂_x for t.
  • domain assumption Monomiality operators defined by P=∂_x and M=x+A'(∂_x)/A(∂_x)
    Used in Section 3 to derive differential equations for the approximated polynomials.
  • domain assumption Bernoulli amplitude is represented umbrally as A(t)=e^{hat B t}φ0 with hat B^r φ0 = B_r
    Cited from [16]; underlies the umbral derivation in Theorem 7.
  • ad hoc to paper Substituting a Padé approximant of the umbral image yields a valid Padé approximation of the original amplitude
    Assumed in Theorems 5 and 7; not demonstrated and inconsistent with standard PA conditions in the Euler [2|1] case.
  • standard math Truncated exponential polynomial generating function and derivative identity ∂_x e_n = n e_{n-1}
    Used for the explicit representations in terms of e_n and e_n^{(2)}.

how reviews work

0 comments
Cite this review

Pith. "Pith review of A novel advancement in the study of Appell polynomials via Pad\`e rational approximants." pith.science (2026). https://pith.science/paper/TGPNPVW2

@misc{pith2026250903178,
  author       = {Pith},
  title        = {Pith review of: A novel advancement in the study of Appell polynomials via Pad\`e rational approximants},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/TGPNPVW2}},
  note         = {Machine review of arXiv:2509.03178}
}
abstract

The use of approximants of Pad\`e type are employed to develop a method aimed at opening new perspectives in the theory of Appell polynomials $a_n(x)$, specified by the generating function \sum_{n=0}^{\infty} \frac{t^n}{n!} a_n(x) = A(t) e^{xt}. In this article, the expansion of amplitude $A(t)$ of the Appell polynomials family in terms of rational approximants yields the possibility of determining the approximation of the $a_n(x)$ in terms of other special polynomials. Application of this approach to Hermite polynomials yields highly accurate approximations in terms of truncated exponential polynomials. Further, monomiality conditions are explored and formalism is extended to consider the Pad\'e approximants within the context of umbral notation.

Figures

Figures reproduced from arXiv: 2509.03178 by the authors.

Figure 1
Figure 1. Comparison of decaying exponential with Maclaurin [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 3
Figure 3. (a): Red H2P(x, 2) = [1|1]He2(x), Dot H2  x, −1 2 , 2  = He2(x) figure 3(b): Red H2P(x, 3) = [1|1]He3(x), Dot H2  x, −1 2 , 3  = He3(x) [PITH_FULL_IMAGE:figures/full_fig_p008_3.png] view at source ↗
Figure 3
Figure 3. Hermite polynomials Hen(x) (dot) and their Pad´e Approximation (red) for n = 11 and 12. The comparison between Hen(x) and its approximated forms is given in [PITH_FULL_IMAGE:figures/full_fig_p011_3.png] view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

16 extracted references · 16 canonical work pages

  1. [1]

    L. C. Andrews, Special Functions for Engineers and Applied Mathematicians, Macmillan, New York (1985)

    work page 1985

  2. [2]

    Appell, Sur une classe de polynˆomes, Ann

    P. Appell, Sur une classe de polynˆomes, Ann. Sci. ´Ecole Norm. Sup., 9(2) 119-144 (1880)

    work page

  3. [3]

    Appell, J

    P. Appell, J. Kamp´ e de F´ riet, Fonctions Hyperg´ eom´ etriqueset Hypersph´ eriques: Polynˆomes d’Hermite, Gauthier-Villars, Paris, France (1926)

    work page 1926

  4. [4]

    G. A. Baker, Jr., Essentials of Pad´ e Approximants, Academic Press, New York (1975)

    work page 1975

  5. [5]

    Babusci, G

    D. Babusci, G. Dattoli, S. Licciardi, E. Sabia, Mathematical Methods for Physicists, World Scientific Publishing, Singapore (2019)

    work page 2019

  6. [6]

    R. L. Burden and J. D. Faires, Numerical Analysis, Brooks/Cole Publishing Co., Pacific Grove (1992)

    work page 1992

  7. [7]

    Dattoli, Modified Pad´ e approximants and perturbative series in quantum field theory, Nuov Cim A 109 1655–1667 (1996)

    G. Dattoli, Modified Pad´ e approximants and perturbative series in quantum field theory, Nuov Cim A 109 1655–1667 (1996)

    work page 1996

  8. [8]

    Dattoli, Integral transforms and Chebyshev-like polynomials, Appl

    G. Dattoli, Integral transforms and Chebyshev-like polynomials, Appl. Math. Comput. 148(1) 225-234 (2004)

    work page 2004

Show all 16 references
  1. [9]

    Dattoli, C

    G. Dattoli, C. Cesarano, D. Sacchetti, A note on truncated polynomials, Appl. Math. Comput. 134 595-605 (2003)

    2003

  2. [10]

    Dattoli, B

    G. Dattoli, B. Germano, M. R. Martinelli, P. E. Ricci, Monomiality and partial differential equations, Math. Comput. Modelling 50(9-10) 1332-1337 (2009)

    2009

  3. [11]

    Dattoli, M

    G. Dattoli, M. Migliorati, H. M. Srivastava, Sheffer polynomials, monomiality principle, algebraic methods and the theory of classical polynomials, Math. Comput. Modelling 45 1033-1041 (2007)

    2007

  4. [12]

    Gilmore, Lie Groups, Physics and Geometry, Cambridge University Press, Cambridge, UK (2008)

    R. Gilmore, Lie Groups, Physics and Geometry, Cambridge University Press, Cambridge, UK (2008)

    2008

  5. [13]

    Licciardi, G

    S. Licciardi, G. Dattoli, Guide to the Umbral Calculus, A Different Mathematical Language, World Scientific (2022)

    2022

  6. [14]

    L. M. Milne Thomson, Calculus of finite differences, Mac Millan & C. Ltd, London (1951)

    1951

  7. [15]

    Roman, The Umbral Calculus, in: Series on Pure and Applied Mathematics, vol

    S. Roman, The Umbral Calculus, in: Series on Pure and Applied Mathematics, vol. 111, Academic Press, New York, Orlando (Florida) (1984)

    1984

  8. [16]

    Subuhi Khan, Ujair Ahmad, Mehnaz Haneef, Gauss-Appell polynomials: An umbral calcu- lus approach, arXiv:2504.05737v3 (2025). 19

    2025 arXiv

Pith tools

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