Analytic umbral transmutations and Bessel moments

Technical Summary: Analytic Umbral Transmutations and Bessel Moments

Problem Statement
The paper addresses the evaluation of integrals involving products of Bessel functions of the first kind, specifically the "Bessel moments" defined as $I_n = \int_{-\infty}^{\infty} J_0(x)^n dx$. The authors identify a critical limitation in the classical "formal indicial umbral calculus" approach. While this formalism successfully converts products of Bessel functions into Gaussian integrals involving fractional powers of sums of umbral operators, it relies on hypergeometric series expansions that often fail to converge at physically relevant parameter values.

A primary example is the cubic moment $I_3(1,1,1) = \int_{-\infty}^{\infty} J_0(x)^3 dx$. The naive formal expansion yields an Appell $F_4$ hypergeometric series. However, the convergence domain of this series ($\sqrt{|x|} + \sqrt{|y|} < 1$) explicitly excludes the equilateral case where all parameters are equal. Consequently, the formal method incorrectly suggests the integral diverges or is undefined, despite the integral being finite and well-known to converge to a specific constant. The paper seeks to resolve this obstruction by moving from formal series to a rigorous analytic framework.

Methodology
The authors employ an "analytic umbral framework" (building on previous work [5]) to reinterpret the umbral operator $u$ not merely as a formal symbol, but as a functional in the Borel additive complex variable. The core methodology involves:

1. Analytic Umbral Pairing: The formal operator function $f(\zeta u)$ is mapped to a composed functional $f(\zeta e^u)$. The pairing with a ground state $\phi(t)$ is defined via a Mellin–Barnes integral:
   $$F(\zeta) = \langle f(\zeta e^u), \phi \rangle = \frac{1}{2\pi i} \int_{C_t} J_f(t; \zeta) \phi(t) dt$$
   where $J_f$ is a spectral jump function derived from the Mellin transform.

2. Ramanujan's Master Theorem as a Selection Principle: The authors utilize Ramanujan's Master Theorem to invert the analytic umbral transmutation. This theorem acts as a selection principle for the "spectral ground state" (or "clock") $\Phi(t)$ associated with a given Bessel product. Specifically, if $F(y) = J_0(\sqrt{y})^\alpha$, the ground state is identified as $\Phi_\alpha(t) = \mathcal{M}[F](-t) / \Gamma(-t)$.

3. Clock Factorization: The product of Bessel functions is decomposed into "clocks" (spectral states).

   • The elementary clock for $J_0$ is $\phi_1(t) = 1/\Gamma(1+t)$.
   • The square clock for $J_0^2$ is $\phi_2(t) = \Gamma(2t+1)/\Gamma(t+1)^3$.
   • Higher powers are treated as contractions of these distinct clocks.

4. Mellin–Barnes Decomposition: The spatial integration (Gaussian integral) results in algebraic kernels (e.g., $(e^u + e^v)^{-1/2}$). These kernels are expanded using Mellin–Barnes integrals, fusing the separate clocks into a single effective spectral state or a multivariate Barnes function.

Key Contributions and Results

• Resolution of the Cubic Moment ($n=3$): By factorizing $J_0^3 = J_0 \cdot J_0^2$, the authors combine the elementary clock $\phi_1$ and the square clock $\phi_2$. The resulting integral is reduced to a one-dimensional Barnes integral, equivalent to a Meijer $G$-function:
  $$I_3 = \frac{1}{\sqrt{\pi}} G_{3,3}^{2,1}\left( \frac{1}{4} \middle| \begin{matrix} 1/2, 1/2, 1/2 \\ 0, 0, 0 \end{matrix} \right)$$
  This analytic continuation yields the classical finite value $\frac{\Gamma(1/6)^2}{2^{2/3}\sqrt{3}\pi^{3/2}\Gamma(5/6)}$, proving that the divergent Appell series was merely a local representation of a globally meaningful identity.

• Scaled Cubic Products: The method is extended to $I(a,b) = \int J_0(ax)J_0(bx)^2 dx$. The dependence on the scale ratio is encoded in a single effective one-clock transmutation function $H(\zeta)$, which is also a Meijer $G$-function.

• The Fourth Moment ($n=4$): The symmetric grouping $J_0^4 = J_0^2 \cdot J_0^2$ leads to a contraction of two square clocks. This results in a one-variable Meijer $G$-function. The residue expansion of this function (closing the contour to the right) produces a rapidly convergent series involving central binomial coefficients and harmonic numbers, providing a new numerical representation for $I_4$.

• The Fifth Moment ($n=5$) and Higher Ranks: The paper demonstrates that $n=5$ represents a structural transition. Grouping $J_0^5 = J_0^2 \cdot J_0^2 \cdot J_0$ requires a ternary fusion of clocks. The resulting kernel $(e^u + e^v + e^w)^{-1/2}$ necessitates a bivariate Barnes integral. The effective object is no longer a standard Meijer $G$-function (rank 1) but a two-variable Barnes function. The authors note that while the analytic representation exists, the naive residue expansion at the physical point $(1,1)$ is not numerically efficient, highlighting a separation between the structural analytic representation and the most efficient numerical evaluation methods.

• Fractional Powers ($n=\alpha$): The framework is generalized to real $\alpha > 2$. The authors show that the moment remains the Mellin value of a Ramanujan-selected spectral state. Unlike integer cases where the state may collapse to a finite Gamma-product, the fractional case involves a genuinely selected spectral state that is not a finite product, yet the analytic umbral interpretation persists.

Significance and Claims
The paper claims to provide a rigorous justification for the passage from formal indicial umbral calculus to "Mellin–Barnes umbral transmutation theory." Its primary significance lies in:

1. Global vs. Local Validity: It establishes that umbral identities exist globally, even when their specific hypergeometric realizations (like the Appell $F_4$ series) are restricted to local convergence chambers. The analytic framework recovers the correct finite values where formal series diverge.
2. Structural Hierarchy: It identifies a hierarchy in Bessel moments:
   • $n=3, 4$: Binary fusions lead to one-variable Meijer $G$-functions.
   • $n=5$: Ternary fusions lead to higher-rank (bivariate) Barnes structures.
3. Unification: It unifies the treatment of integer and fractional moments under a single analytic paradigm where Ramanujan's Master Theorem acts as the inverse selection principle for the spectral ground state.

The authors conclude that while the analytic umbral method provides a consistent structural representation for all moments, the transition from $n=4$ to $n=5$ marks the limit of the "Meijer $G$-phenomenon," requiring more complex multivariate structures for higher-order moments.