2025-11-10T02:38:06.672947

Summing Sneddon-Bessel series explicitly

Durán, Pérez, Varona
We sum in a close form the Sneddon-Bessel series \[ \sum_{m=1}^\infty \frac{J_α(x j_{m,ν})J_β(y j_{m,ν})} {j_{m,ν}^{2n+α+β-2ν+2} J_{ν+1}(j_{m,ν})^2}, \] where $0<x$, $0<y$, $x+y<2$, $n$ is an integer, $α,β,ν\in \mathbb{C}\setminus \{-1,-2,\dots \}$ with $2\operatorname{Re} ν< 2n+1 + \operatorname{Re} α+ \operatorname{Re} β$ and $\{j_{m,ν}\}_{m\geq 0}$ are the zeros of the Bessel function $J_ν$ of order $ν$. As an application we prove some extensions of the Kneser-Sommerfeld expansion.
academic

Summing Sneddon-Bessel series explicitly

Basic Information

  • Paper ID: 2207.08709
  • Title: Summing Sneddon-Bessel series explicitly
  • Authors: Antonio J. Durán, Mario Pérez, Juan L. Varona
  • Classification: math.CA (Classical Analysis and ODEs)
  • Publication Date: July 15, 2022
  • Paper Link: https://arxiv.org/abs/2207.08709

Abstract

This paper provides explicit summation of Sneddon-Bessel series: m=1Jα(xjm,ν)Jβ(yjm,ν)jm,ν2n+α+β2ν+2Jν+1(jm,ν)2\sum_{m=1}^\infty \frac{J_\alpha(x j_{m,\nu})J_\beta(y j_{m,\nu})}{j_{m,\nu}^{2n+\alpha+\beta-2\nu+2} J_{\nu+1}(j_{m,\nu})^2}

where 0<x0<x, 0<y0<y, x+y<2x+y<2, nn is an integer, α,β,νC{1,2,}\alpha,\beta,\nu\in \mathbb{C}\setminus \{-1,-2,\ldots \} satisfy 2Reν<2n+1+Reα+Reβ2\operatorname{Re} \nu< 2n+1 + \operatorname{Re} \alpha+ \operatorname{Re} \beta, and {jm,ν}m0\{j_{m,\nu}\}_{m\geq 0} are the zeros of the Bessel function JνJ_\nu of order ν\nu. As applications, several extensions of the Kneser-Sommerfeld expansion are established.

Research Background and Motivation

  1. Problem Background: Explicit summation of Bessel series is an important problem in classical analysis, particularly for infinite series involving zeros of Bessel functions. Sneddon considered bivariate Bessel series in his monograph, but lacked explicit closed-form expressions.
  2. Problem Significance: Summation of Bessel series has widespread applications in applied mathematics, mathematical physics, and engineering, particularly in solving boundary value problems, heat conduction equations, and wave equations.
  3. Limitations of Existing Methods: Currently, there is a lack of unified closed-form expressions for Sneddon-Bessel series with general parameters, especially for complex parameter cases.
  4. Research Motivation: To systematically resolve the explicit summation of Sneddon-Bessel series in commemoration of the centennial of the first edition of Watson's monograph on the theory of Bessel functions.

Core Contributions

  1. Main Theorem: Provides complete closed-form expressions for Sneddon-Bessel series, covering both cases n0n \geq 0 and n<0n < 0
  2. Technical Innovation: Establishes partial fraction decomposition of Bessel functions using the residue theorem
  3. Theoretical Extension: Proves multiple generalizations of the Kneser-Sommerfeld expansion
  4. Unified Framework: Provides a unified framework for handling Bessel series with complex parameters

Detailed Methodology

Problem Definition

The research objective is to find the closed form of the series: Sα,β,νqn(x,y)=m=1Jα(xjm,ν)Jβ(yjm,ν)jm,νqnJν+1(jm,ν)2S_{\alpha,\beta,\nu}^{q_n}(x,y) = \sum_{m=1}^\infty \frac{J_\alpha(xj_{m,\nu})J_\beta(yj_{m,\nu})}{j_{m,\nu}^{q_n}J_{\nu+1}(j_{m,\nu})^2}

where qn=2n+α+β2ν+2q_n = 2n+\alpha+\beta-2\nu+2.

Core Technical Framework

1. Function Transformation

Introduce the entire function: Φν(z)=2νΓ(ν+1)Jν(z)zν\Phi_\nu(z) = 2^\nu\Gamma(\nu+1)\frac{J_\nu(z)}{z^\nu}

Define the generating function: Φα(xz)Φβ(yz)Φν(z)2=n=0δα,β,ν(n)(x,y)z2n\frac{\Phi_\alpha(xz)\Phi_\beta(yz)}{\Phi_\nu(z)^2} = \sum_{n=0}^\infty \delta_{\alpha,\beta,\nu}^{(n)}(x,y)z^{2n}

2. Application of Residue Theorem

Theorem 1 (Partial Fraction Decomposition): Let ff be an entire function satisfying growth conditions. Then: 1n!dndtn(f(t)Φν(t)2)=mZ{0}4(ν+1)2[(2ν+1)t(2νn)jm]f(jm)jm(jmt)f(jm)(jmt)n+2jm3Φν+1(jm)2\frac{1}{n!}\frac{d^n}{dt^n}\left(\frac{f(t)}{\Phi_\nu(t)^2}\right) = \sum_{m\in\mathbb{Z}\setminus\{0\}} \frac{4(\nu+1)^2[(2\nu+1)t-(2\nu-n)j_m]f(j_m)-j_m(j_m-t)f'(j_m)}{(j_m-t)^{n+2}j_m^3\Phi_{\nu+1}(j_m)^2}

3. Partial Differential Equation Method

Via the residue theorem, establish the partial differential equation: δα,β,ν(n)(x,y)8(ν+1)2=(2ν2n)ξn,α,β,ν(x,y)+xξn,α,β,νx+yξn,α,β,νy\frac{\delta_{\alpha,\beta,\nu}^{(n)}(x,y)}{-8(\nu+1)^2} = (2\nu-2n)\xi_{n,\alpha,\beta,\nu}(x,y) + x\frac{\partial\xi_{n,\alpha,\beta,\nu}}{\partial x} + y\frac{\partial\xi_{n,\alpha,\beta,\nu}}{\partial y}

Main Results

Case 1: n0n \geq 0

Sα,β,νqn(x,y)=Γ(ν+1)2xαyβ2q0Γ(α+1)Γ(β+1)(x2n2νφα,β,ν(n)(y/x)2j+2k2nA2j,2k,nα,β,νj+k+νnx2jy2k)S_{\alpha,\beta,\nu}^{q_n}(x,y) = \frac{\Gamma(\nu+1)^2x^\alpha y^\beta}{2^{q_0}\Gamma(\alpha+1)\Gamma(\beta+1)}\left(x^{2n-2\nu}\varphi_{\alpha,\beta,\nu}^{(n)}(y/x) - \sum_{2j+2k\leq 2n}\frac{A_{2j,2k,n}^{\alpha,\beta,\nu}}{j+k+\nu-n}x^{2j}y^{2k}\right)

where the function φα,β,ν(n)\varphi_{\alpha,\beta,\nu}^{(n)} is determined by the recurrence relation: φα,β,ν(0)(t)=1ν(αν)2F1(να,νβ+1;t2)\varphi_{\alpha,\beta,\nu}^{(0)}(t) = \frac{1}{\nu}\binom{\alpha}{\nu}{}_2F_1\left(\frac{\nu-\alpha,\nu}{\beta+1};t^2\right)

φα,β,ν(n)(t)=12ν2n(12(α+1)φα+1,β,ν(n1)(t)+t22(β+1)φα,β+1,ν(n1)(t))\varphi_{\alpha,\beta,\nu}^{(n)}(t) = \frac{1}{2\nu-2n}\left(\frac{1}{2(\alpha+1)}\varphi_{\alpha+1,\beta,\nu}^{(n-1)}(t) + \frac{t^2}{2(\beta+1)}\varphi_{\alpha,\beta+1,\nu}^{(n-1)}(t)\right)

Case 2: n<0n < 0

Sα,β,νqn(x,y)=xα2ν+2nyβΓ(νn)2qnΓ(β+1)Γ(n+αν+1)2F1(νn,ναnβ+1;y2x2)S_{\alpha,\beta,\nu}^{q_n}(x,y) = \frac{x^{\alpha-2\nu+2n}y^\beta\Gamma(\nu-n)}{2^{q_n}\Gamma(\beta+1)\Gamma(n+\alpha-\nu+1)}{}_2F_1\left(\frac{\nu-n,\nu-\alpha-n}{\beta+1};\frac{y^2}{x^2}\right)

Applications: Extensions of the Kneser-Sommerfeld Expansion

Classical Kneser-Sommerfeld Expansion

m=1Jν(xjm,ν)Jν(yjm,ν)(jm,ν2z2)Jν+1(jm,ν)2=πJν(yz)4Jν(z)[Yν(z)Jν(xz)Jν(z)Yν(xz)]\sum_{m=1}^\infty \frac{J_\nu(xj_{m,\nu})J_\nu(yj_{m,\nu})}{(j_{m,\nu}^2-z^2)J_{\nu+1}(j_{m,\nu})^2} = \frac{\pi J_\nu(yz)}{4J_\nu(z)}[Y_\nu(z)J_\nu(xz)-J_\nu(z)Y_\nu(xz)]

New Extended Forms

m=1jm,ννβJν(xjm,ν)Jβ(yjm,ν)(jm,ν2z2)Jν+1(jm,ν)2=πJβ(yz)4zβνJν(z)[Yν(z)Jν(xz)Jν(z)Yν(xz)]\sum_{m=1}^\infty \frac{j_{m,\nu}^{\nu-\beta}J_\nu(xj_{m,\nu})J_\beta(yj_{m,\nu})}{(j_{m,\nu}^2-z^2)J_{\nu+1}(j_{m,\nu})^2} = \frac{\pi J_\beta(yz)}{4z^{\beta-\nu}J_\nu(z)}[Y_\nu(z)J_\nu(xz)-J_\nu(z)Y_\nu(xz)]

Applicable conditions: Reν<Reβ+1\operatorname{Re}\nu < \operatorname{Re}\beta + 1, 0yx10 \leq y \leq x \leq 1.

Technical Innovations

  1. Systematic Application of Residue Methods: First systematic application of the residue theorem to partial fraction decomposition of Bessel functions
  2. Generating Function Techniques: Cleverly utilizes generating functions to transform infinite series problems into finite polynomial recurrences
  3. Integral Transform Methods: Employs integral transforms such as Sonin's formula to handle complex parameter cases
  4. Analytic Continuation Techniques: Extends results to broader parameter ranges through analytic continuation

Mathematical Analysis

Convergence Conditions

The key condition for series convergence is: 2Reν<2n+1+Reα+Reβ2\operatorname{Re}\nu < 2n+1+\operatorname{Re}\alpha+\operatorname{Re}\beta

Special Case Treatment

  • When ν{0,1,,n}\nu \in \{0,1,\ldots,n\}, handled through limiting processes
  • When Reν<1\operatorname{Re}\nu < -1, results are extended using the integral transform Tμ,η,hT_{\mu,\eta,h}

Boundary Cases

  • The boundary case x+y=2x+y=2 requires stronger conditions: 2Reν<2n+Reα+Reβ2\operatorname{Re}\nu < 2n+\operatorname{Re}\alpha+\operatorname{Re}\beta
  1. Classical Results: Watson's theory of Bessel functions, Sneddon's mixed boundary value problems
  2. Modern Developments: Martin's research on Fourier-Bessel series, Grebenkov's spectral sum formulas
  3. Contribution of This Paper: First complete closed-form solution for general Sneddon-Bessel series

Conclusions and Discussion

Main Conclusions

The paper successfully resolves the explicit summation problem for Sneddon-Bessel series, providing complete closed-form expressions and establishing important generalizations of the Kneser-Sommerfeld expansion.

Theoretical Significance

  1. Perfects the theoretical framework of Bessel series
  2. Provides powerful tools for related physics and engineering problems
  3. Demonstrates the power of complex analysis methods in the theory of special functions

Application Prospects

The results have important application value in solving partial differential equations, integral equations, and mathematical physics problems.

In-Depth Evaluation

Strengths

  1. Theoretical Completeness: Provides a complete theoretical framework with rigorous mathematical proofs
  2. Methodological Innovation: Innovatively combines the residue theorem, generating functions, and integral transforms
  3. Practical Value of Results: The provided formulas have strong practical utility
  4. Clarity of Exposition: The paper has clear structure and rigorous mathematical presentation

Technical Highlights

  1. Skillfully handles analytic continuation for complex parameter cases
  2. Systematically resolves convergence issues across different parameter ranges
  3. Provides concrete computational methods and recurrence relations

Impact

This work makes important contributions to Bessel function theory and is expected to have broad impact in mathematical physics, engineering mathematics, and related fields. The methods and results provide important theoretical foundations and technical tools for subsequent related research.