Pith. sign in

REVIEW 4 major objections 4 minor 24 references

Product formulas for multivariate and matrix-variate confluent hypergeometric functions

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

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

Pith's one-line read An old product formula for ${}_1F_1$ now covers multivariate $\Phi_2^{(k)}$, $\Psi_2^{(k)}$, and matrix-variate confluent hypergeometric functions.

desk verdict A clean but modest extension of Erdélyi's product formula; the formulas look right, with one genuinely missing convergence justification in the contour cases. read the letter →

arxiv 2608.19270 v1 pith:22FEZXNL submitted 2026-08-18 math.CA

Min-Jie Luo , Shu-Ying Shen This is my paper · ORCID
classification math.CA MSC 33C1533C65
keywords confluenthypergeometricfunctionmultivariateWhittakerproductformulazonalpolynomialmatrixargumentHumbertfunctionsΦ2andΨ2Erdélyifractionalcalculus
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 note extends Erdélyi's 1939 product formula for the univariate confluent hypergeometric function ${}_1F_1$ to two higher-dimensional families: the multivariate Humbert functions $\Phi_2^{(k)}$ and $\Psi_2^{(k)}$, and the confluent hypergeometric function of a matrix argument. In the multivariate setting the product of two $\Phi_2^{(k)}$'s is written as an integral of a ${}_1F_1$ over a simplex product and a one-dimensional $\beta$ law, or as an integral of a $\Psi_2^{(k)}$ over three simplices, with contour-integral analogues for both $\Phi_2^{(k)}$ and $\Psi_2^{(k)}$. For matrices, the product of two matrix-variate ${}_1F_1$'s is shown to equal a threefold matrix-$\beta$ integral of a third ${}_1F_1$. These identities matter because product formulas are the standard bridge from special functions to convolution structures, addition formulas, and bounds for orthogonal polynomials, and the matrix-argument case feeds directly into multivariate statistical distribution theory. A product formula for Humbert's multivariate Whittaker function and a connection to fractional calculus are recorded as corollaries.

What carries the argument

The machinery is a family of coupling identities that rewrite the exponential as an integral of a confluent hypergeometric function: $e^z = \int_0^1 {}_1F_1(g;h;zt)\,d\mu_{h,g-h}(t)$, its simplex analogue $e^{\langle z\rangle} = \int_{E^k} \Psi_2^{(k)}[g;h;z\circ t]\,d\mu(h,g-\langle h\rangle)(t)$, the contour analogue $e^{\langle z\rangle} = \int_{L^k} \Phi_2^{(k)}[h;g;z/t]\,d\widetilde{\mu}(h-1,g+k-\langle h\rangle)(t)$, and the matrix-argument version $\mathrm{etr}(X) = \int_{O<T<I} {}_1F_1(h;g;XT)\,d\mu^{(m)}_{g,h-g}(T)$. Substituting one of these into the double integral obtained from the Euler representations of the two factors turns the product into a single integral of a hypergeometric function. The parameter bookkeeping is carried by $\beta$ and Dirichlet measures in the simplex cases and by the complex measure $\widetilde{\mu}$ on the contour $L^k$ in the other cases, with the matrix case resting on the zonal polynomial integral $\int_{O<X<I} C_\kappa(XT)\,d\mu^{(m)}_{a,b}(T) = \frac{[a]_\kappa}{[a+b]_\kappa} C_\kappa(X)$ and the identity $\sum_{\langle\kappa\rangle=\ell} C_\kappa(X) = (\mathrm{tr}\,X)^\ell$.

What would settle it

Evaluate identity (15) for $k=1$ numerically: with $z=1$, $h=3/2$, $g=5/2$, compute the right-hand side by quadrature over the vertical line $L=\{1/2+is: s\in\mathbb{R}\}$, truncating the ${}_1F_1$ series with a rigorous tail bound; if the truncated sum and integral do not commute to within the error tolerance, the term-by-term interchange fails and formula (14) cannot hold as stated.

Watch

Extended reading notes

Core claim

The central claim is that formulas (9), (10), (14), (16), and (19) hold under the stated parameter conditions, each expressing a product of two hypergeometric functions as an integral of one hypergeometric function against a product of $\beta$ or Dirichlet measures, with the complex measure $\widetilde{\mu}$ supported on $L^k$ in the contour cases. Formula (19) is representative: for positive definite $X,Y \in \mathrm{Sym}_m(\mathbb{R})$ and parameters with $\min\{\Re(a),\Re(b),\Re(g),\Re(c-a),\Re(b-d),\Re(h-g)\} > \tfrac12(m-1)$, the product ${}_1F_1(a;c;X)\,{}_1F_1(b;d;Y)$ equals a triple integral over matrix $\beta$ measures of ${}_1F_1\bigl(h;g;(U^{1/2}XU^{1/2}+V^{1/2}YV^{1/2})T\bigr)$. The multivariate formulas similarly couple two $\Phi_2^{(k)}$ or $\Psi_2^{(k)}$ factors through a one-dimensional or simplex integral of a ${}_1F_1$ or $\Psi_2^{(k)}$. The paper presents two proofs of the original Erdélyi formula and adopts the more instructive one, inserting an exponential representation of $e^z$ into the Euler integral representation of the two factors, and the Whittaker function product formula follows by specializing (16).

Load-bearing premise

The load-bearing premise is that the series defining $\Phi_2^{(k)}$ and $\Psi_2^{(k)}$ may be integrated term by term over the noncompact contour $L^k$ in formulas (14) and (16); if that interchange of sum and integral fails for the stated parameter ranges, those two product formulas are not established.

Editorial extensions

If this is right

  • For $\Phi_2^{(k)}$ and $\Psi_2^{(k)}$, the product of two functions with different parameter sets becomes a single hypergeometric function integrated against an explicit measure, giving ready-made convolution kernels for multivariate Laguerre and Whittaker expansions.
  • From (16) the paper derives a product formula for Humbert's multivariate Whittaker function $M_{\kappa,\mu}$, expressing the product $M_{\kappa,\nu}(x)M_{\kappa,\beta}(y)$ as an integral of a single $M_{\kappa,\lambda}$ times an explicit prefactor.
  • The matrix-variate formula (19) holds for positive definite matrix arguments under the analogue $\min\{\Re(a),\Re(b),\Re(g),\Re(c-a),\Re(b-d),\Re(h-g)\} > \tfrac12(m-1)$ of the univariate parameter conditions, providing the matrix-argument product identity used in random matrix distribution theory.
  • The contour formulas (14) and (16) cover parameter ranges where the simplex-based formulas (9) and (10) do not apply, provided the stated dominance conditions such as $\Re(g+k)>\Re(\langle h\rangle)$ hold.

Reading between the lines

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

  • The same coupling scheme should yield product formulas for other Humbert or Appell-type series, for example $\Phi_1$, $\Phi_3$, or the two-variable Appell functions, by writing their Euler kernels as exponentials and applying identity (8); this is a direct test of the method's scope.
  • Applying the matrix-argument argument to the matrix-variate $\Phi_2^{(k)}$ should give a threefold or fourfold matrix-Dirichlet integral; the authors explicitly leave that extension to the reader, so it is a natural check on the method.
  • If the term-by-term interchange over $L^k$ can be justified, formulas (14) and (16) may extend by analytic continuation to a larger parameter region, and the boundary $\Re(g+k)=\Re(\langle h\rangle)$ would mark the true domain.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Request a human review

A listed scientist reviews the paper for a fee and the review publishes here regardless of verdict. See the reviewers or get listed.

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

4 major / 4 minor

Summary. The paper revisits Erdélyi's 1939 product formula for the univariate confluent hypergeometric function and proposes two higher-dimensional generalizations. Section 2 gives two proofs of the univariate formula (4). Section 3 derives product formulas for the multivariate functions Φ_2^{(k)} and Ψ_2^{(k)}: formulas (9) and (10) use integrals over simplexes, while formulas (14) and (16) use contour integrals over L^k. A corollary product formula for multivariate Whittaker functions is also given. Section 4 states a matrix-variate analogue, Proposition 4.1, expressing the product of two matrix-variate 1F1's as a threefold matrix-beta integral. The derivations are reductions to standard integral representations of the exponential function, (2), (8), (12), (13), and (21), together with zonal-polynomial identities.

Significance. If the identified gaps are repaired, the paper provides clean, explicit product formulas for two classes of multivariate confluent hypergeometric functions and for the matrix-variate 1F1. The use of Oshima's contour transformations to obtain Φ_2-type and Ψ_2-type product formulas is an interesting and potentially useful technique, and the paper is self-contained up to standard references. The formulas are concrete and falsifiable, and the univariate and matrix-variate cases are checked against known integral representations. The main weakness is that two load-bearing steps are asserted rather than proved: the noncompact contour interchange behind Propositions 3.3 and 3.4, and the well-definedness of the matrix-variate integrand in (19). The Whittaker corollary also has an incorrect parameter condition. These issues are local and repairable, so the central ideas appear sound.

major comments (4)
  1. [Sec. 3, Eq. (15), Props. 3.3 and 3.4] The contour formula e^{⟨z⟩} = ∫_{L^k} Φ_2^{(k)}[h;g;z/t] dẽµ(h-1,g+k-⟨h⟩)(t) is asserted by termwise integration, but L^k is noncompact and no domination argument is supplied. This interchange is load-bearing because formulas (14) and (16) depend on it. The stated hypotheses ℜ(h-1)>0 and ℜ(g+k)>ℜ(⟨h⟩) likely suffice: on L one has |t| ≥ (k+1)^{-1}, and the ratio |(h)_m/(g)_{⟨m⟩}| grows at most polynomially, so an absolutely convergent double-series estimate can be written down. However, the paper does not provide it, and the proof is therefore incomplete as written.
  2. [Sec. 4, Prop. 4.1] The hypothesis states min{ℜ(a),ℜ(b),ℜ(g),ℜ(c-a),ℜ(b-d),ℜ(h-g)} > (m-1)/2. But the measure dµ_{b,d-b} in the integrand requires ℜ(d-b) > (m-1)/2, not ℜ(b-d) > (m-1)/2. With the printed sign, the measure dµ_{b,d-b} is not covered by (18), and Proposition 4.1 is not stated under valid conditions.
  3. [Sec. 4, Eq. (19)] The integrand 1F1(h/g; (U^{1/2}XU^{1/2}+V^{1/2}YV^{1/2})T) has a matrix argument W T that is not symmetric in general, because the positive definite matrices W and T need not commute. Since 1F1 of matrix argument is defined via zonal polynomials on symmetric matrices, the expression is not well-defined without an additional convention. The proof implicitly uses the identity Cκ(W T)=Cκ(W^{1/2} T W^{1/2}), which follows from the fact that W T and W^{1/2} T W^{1/2} have the same eigenvalues. This convention should be stated explicitly in the manuscript.
  4. [Sec. 3, Whittaker corollary] The stated parameter conditions for the Whittaker product formula are incomplete. Applying (16) to the Whittaker representation requires the t-measure dµ(2λ+1, -⟨λ⟩-κ-k/2) to satisfy the condition ℜ(g)>ℜ(⟨h⟩), i.e., ℜ(κ)<-ℜ(⟨λ⟩)-k/2. The printed condition ℜ(κ)<min{0, k/2-ℜ(⟨ν⟩), k/2-ℜ(⟨β⟩)} does not imply this estimate. As written, the corollary is stated for a parameter range that may violate the hypotheses of Proposition 3.4.
minor comments (4)
  1. [Abstract and Section 1] There are small typos, such as "Erdélyi' s" in the abstract and inconsistent spacing in the notation Φ(k)2, Ψ(k)2; these should be cleaned up.
  2. [Sec. 3, after Eq. (12)] The statement that (12) 'can be easily verified by expanding' is acceptable because the integration is over the compact simplex and the series is locally uniformly convergent, but a short justification would improve the exposition.
  3. [Sec. 4, Eq. (18)] In the definition of the matrix beta measure, the positivity condition is written as min{ℜ(a),ℜ(c-a)} > (m-1)/2, which is correct only up to the usual convention that the real parts are positive; this is fine, but the subsequent use in Proposition 4.1 should match the signs exactly.
  4. [References] References [10] and [13] are self-citations to unpublished or very recent work; this is acceptable since they are used only for notation and motivation, but the authors might note that the present results do not depend on those papers.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity found: the claimed product formulas are derived from independent integral identities and are not equivalent to their own inputs.

full rationale

The paper's central derivations are self-contained reductions to external, verifiable integral representations rather than circular chains. Proposition 3.1 starts from the known Euler-type representation of Phi_2 (cited to Erdelyi [7]) and from the elementary identities (8) and (12), both of which are beta/Dirichlet moment identities that are checked by expanding the hypergeometric series and integrating term by term. Propositions 3.3 and 3.4 similarly use the contour moment identity (13) to verify the representations (15) and the analogous Psi_2 representation; these representations are independent of the target product formulas, and the product formulas follow by substituting them into the product of two integral representations. Proposition 4.1 derives the matrix-variate identity (21) directly from standard zonal polynomial facts in the NIST Handbook [16] and then substitutes it into (20); no fitted parameter is relabeled as a prediction, and no uniqueness claim is imported from the authors' own prior work. The self-citations ([10], [13]) are used only for notation and motivation and are not load-bearing. The reader's noted weakness, the unproved interchange of a series with an integral over the noncompact contour L^k in (15), is a genuine convergence/rigor concern about the proof as written, but it is not circularity: a failure of that interchange would affect the validity of (14) and (16), not reduce those formulas to the assumptions used to derive them. Accordingly, the appropriate circularity score is 0.

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

No free parameters, no fitted values, and no new postulated entities. The paper's claims rest on prior integral representations and standard interchange theorems, with the contour interchange the least documented.

assumptions (6)
  • standard math Euler integral representation of 1F1, Eq (2), giving 1F1 as an integral against a beta measure.
    Used throughout, including the univariate proposition and as template for multivariate and matrix analogues; cited to standard references.
  • domain assumption Integral representation of Φ2^(k) over the simplex, cited to Erdélyi [7, p. 449, Eq. (8.5)].
    Basis for Propositions 3.1 and 3.3; not proved in this paper.
  • domain assumption Exponential representation e^{<z>} = ∫_{E_k} Ψ2^(k)[g;h;z∘t] dµ(h,g-<h>)(t), Eq (12), cited to Exton [9, Eq. (2.7.10)].
    Used to create coupling terms in the Φ2 and Ψ2 product formulas.
  • domain assumption Oshima's contour integral identity, Eq (13), giving moments of t^{-m} over L^k.
    Used to derive Eq (15) and Prop 3.3; relies on complex line integrals from Oshima [17].
  • standard math Interchange of infinite series with integrals (over compact simplex and unit intervals, and over complex contours) is permissible.
    The paper asserts permissibility but does not provide detailed convergence arguments; this is the main technical premise for the multivariate formulas.
  • standard math Matrix-variate beta measure integral and zonal polynomial integration formula, cited to Muirhead [15] and NIST [16].
    Used in the proof of Eq (21) and Proposition 4.1.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Product formulas for multivariate and matrix-variate confluent hypergeometric functions." pith.science (2026). https://pith.science/paper/22FEZXNL

@misc{pith2026260819270,
  author       = {Pith},
  title        = {Pith review of: Product formulas for multivariate and matrix-variate confluent hypergeometric functions},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/22FEZXNL}},
  note         = {Machine review of arXiv:2608.19270}
}
abstract

In this note, we revisit one of Erd\'{e}lyi's product formulas for the univariate confluent hypergeometric function ${}_{1}F_{1}$ and extend it to the multivariate confluent hypergeometric functions $\Phi_2^{(k)}$, $\Psi_2^{(k)}$, as well as to a certain matrix-variate confluent hypergeometric function. A useful connection related to fractional calculus is also given.

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

24 extracted references · 19 canonical work pages

  1. [1]

    Amri, About integral product formula for Jack polynomials of two variables, Afr

    B. Amri, About integral product formula for Jack polynomials of two variables, Afr. Mat. 36 (4) (2025), Paper No. 160, 14p

    work page 2025

  2. [2]

    Beals, Y

    R. Beals, Y. Kannai, Inverse Laplace transforms of products of Whittaker functions, Proc. R. Soc. Lond., Ser. A, Math. Phys. Eng. Sci. 464 (2092) (2008), 795–806

    work page 2008

  3. [3]

    Carlitz, An integral for the product of two Laguerre polynomials, Boll

    L. Carlitz, An integral for the product of two Laguerre polynomials, Boll. Unione Mat. Ital., III. Ser. 17 (1962), 25–28. Product formulas for multivariate and matrix-variate confluent hypergeometric functions 9

    work page 1962

  4. [4]

    Carlson, Special Functions of Applied Mathematics, Academic Press, New York, 1977

    B.C. Carlson, Special Functions of Applied Mathematics, Academic Press, New York, 1977

    1977

  5. [5]

    Chatterjea, An integral representation for the product of two generalized Bessel polynomials, Boll

    S. Chatterjea, An integral representation for the product of two generalized Bessel polynomials, Boll. Unione Mat. Ital., III. Ser. 18 (1963), 377–381

    work page 1963

  6. [6]

    Chatterjea, Integral representation for the product of two Jacobi polynomials, J

    S.K. Chatterjea, Integral representation for the product of two Jacobi polynomials, J. Lond. Math. Soc. 39 (1964), 753–756

    work page 1964

  7. [7]

    Erdélyi, Beitrag zur Theorie der konfluenten hypergeometrischen Funktionen von mehreren Veränderlichen, Sitzungsber

    A. Erdélyi, Beitrag zur Theorie der konfluenten hypergeometrischen Funktionen von mehreren Veränderlichen, Sitzungsber. Akad. Wiss. Wien, Math.-Naturw. Kl., Abt. IIa 146 (1937), 431–467

    1937

  8. [8]

    Erdélyi, An integral representation for the product of two Whittaker functions, J

    A. Erdélyi, An integral representation for the product of two Whittaker functions, J. Lond. Math. Soc. 14 (1939), 23–30

    work page 1939

Show all 24 references
  1. [9]

    Exton, Multiple Hypergeometric Functions and Applications, Halsted Press, London (1976)

    H. Exton, Multiple Hypergeometric Functions and Applications, Halsted Press, London (1976)

    1976

  2. [10]

    Guo, M.-J

    L.-J. Guo, M.-J. Luo, R.K. Raina, J.-J. Wang, Multivariate Laguerre polynomials: new results and insights, 2026. Available online athttps://arxiv.org/abs/2604.18629

    2026 arXiv

  3. [11]

    Humbert, La fonctionW k,µ1,µ2,···,µ n (x1,x 2,· · ·,x n), C

    P . Humbert, La fonctionW k,µ1,µ2,···,µ n (x1,x 2,· · ·,x n), C. R. 171 (1920), 428–430

    1920

  4. [12]

    Koornwinder, A.L

    T .H. Koornwinder, A.L. Schwartz, Product formulas and associated hypergroups for orthogonal polynomials on the simplex and on a parabolic biangle, Constr. Approx. 13 (4) (1997), 537–567

    1997

  5. [13]

    Luo, R.K

    M.-J. Luo, R.K. Raina, Upper bounds for Erdélyi’ s multivariate Laguerre polynomials, J. Math. Inequal. 20 (1) (2026), 147–154

    2026

  6. [14]

    Mathai, Jacobians of Matrix Transformations and Functions of Matrix Argument, World Scientific Publishing, New York, 1997

    A.M. Mathai, Jacobians of Matrix Transformations and Functions of Matrix Argument, World Scientific Publishing, New York, 1997

    1997

  7. [15]

    Muirhead, Aspects of Multivariate Statistical Theory, John Wiley & Sons, New York, 1982

    R.J. Muirhead, Aspects of Multivariate Statistical Theory, John Wiley & Sons, New York, 1982

    1982

  8. [16]

    Olver, D.W

    F .W .J. Olver, D.W . Lozier, R.F . Boisvert, C.W . Clark (Eds.), NIST Handbook of Mathematical Functions, Cambridge University Press, New York, 2010

    2010

  9. [17]

    T . Oshima, Integral transformations of hypergeometric functions with several variables, InSymmetry in Geometry and Analysis, Vol 2, Festschrift in Honor of Toshiyuki Kobayashi, Progress in Mathematics, Vol. 358, pages 551–586. Birkhäuser, 2025

    2025

  10. [18]

    Srivastava, P .W

    H.M. Srivastava, P .W . Karlsson, Multiple Gaussian Hypergeometric Series, Ellis Horwood Ltd., Chichester, 1985

    1985

  11. [19]

    Srivastava, H.L

    H.M. Srivastava, H.L. Manocha, A Treatise on Generating Functions, Halsted Press (Ellis Horwood Limited, Chich- ester), John Wiley and Sons, New York, Chichester, Brisbane, and Toronto, 1984

    1984

  12. [20]

    Srivastava, R

    H.M. Srivastava, R. Panda, An integral representation for the product of two Jacobi polynomials, J. Lond. Math. Soc., II. Ser. 12 (1976), 419–425

    1976

  13. [21]

    Veestraeten, An integral representation for the product of parabolic cylinder functions, Integral Transforms Spec

    D. Veestraeten, An integral representation for the product of parabolic cylinder functions, Integral Transforms Spec. Funct. 28 (1) (2017), 15–21

    2017

  14. [22]

    Veestraeten, An alternative integral representation for the product of two parabolic cylinder functions, Integral Transforms Spec

    D. Veestraeten, An alternative integral representation for the product of two parabolic cylinder functions, Integral Transforms Spec. Funct. 28 (12) (2017), 915–922

    2017

  15. [23]

    Watson, A note on the polynomials of Hermite and Laguerre, J

    G.N. Watson, A note on the polynomials of Hermite and Laguerre, J. Lond. Math. Soc. 13 (1938), 204–209

    1938

  16. [24]

    Watson, Another note on Laguerre polynomials, J

    G.N. Watson, Another note on Laguerre polynomials, J. London Math. Soc., 14 (1939), 19–22

    1939

Pith tools

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