Uniform Asymptotic Approximation and Numerical Evaluation of the Reverse Generalized Bessel Polynomial Zeros

Basic Information

• Paper ID: 2510.14295
• Title: Uniform Asymptotic approximation and numerical evaluation of the Reverse Generalized Bessel Polynomial zeros
• Authors: T. M. Dunster, A. Gil, D. Ruiz-Antolin, J. Segura
• Classification: math.CA (Classical Analysis and ODEs)
• Publication Date: October 16, 2025
• Paper Link: https://arxiv.org/abs/2510.14295

Abstract

This paper derives uniform asymptotic expansions for the zeros of reverse generalized Bessel polynomials for large degree $n$ and real parameter $a$. Assuming the parameter range $-\Delta_1 n + \frac{3}{2} \leq a \leq \Delta_2 n$, where $\Delta_1 \in (0,1)$ is an arbitrary fixed constant and $\Delta_2$ is a bounded positive constant, the new expansions are uniformly valid for all zeros. Within this parameter range, at most one zero is real, with the remainder being complex conjugate pairs. The expansions demonstrate high accuracy for moderate to large values of $n$ and serve as initial values for an efficient numerical algorithm that utilizes Taylor series to compute the remaining complex zeros.

Research Background and Motivation

Problem Description

Generalized Bessel polynomials are defined as: $y_n(z; a) = \sum_{k=0}^{n} \binom{n}{k} (n+a-1)_k \left(\frac{1}{2z}\right)^k$

where $(α)_k = Γ(α + k)/Γ(α)$ is the Pochhammer symbol. The zeros of these polynomials are typically complex-valued and have numerous applications in applied mathematics and engineering.

Research Motivation

1. Practical Application Needs: The zeros of generalized Bessel polynomials have important applications in applied mathematics and engineering
2. Computational Challenges: Existing methods for computing these zeros have limitations in efficiency and accuracy
3. Theoretical Refinement: Establishing uniformly valid asymptotic theory to describe the behavior of all zeros is needed
4. Numerical Algorithm Optimization: Efficient numerical methods for computing complex zeros are required

Limitations of Existing Methods

• Previous work 2, 13 has restrictions on computational precision and efficiency
• Lack of uniformly valid asymptotic expansions for all zeros
• Suboptimal initial value selection for numerical algorithms

Core Contributions

1. Derivation of Uniform Asymptotic Expansions: Established asymptotic expansions uniformly valid for all zeros of reverse generalized Bessel polynomials
2. Establishment of Turning Point Theory: Based on WKB theory and Liouville-Green methods, handled differential equations with turning points
3. Construction of Efficient Numerical Algorithm: Combined asymptotic expansions and Taylor series to develop an efficient algorithm for computing complex zeros
4. Provision of High-Precision Approximations: Expansions achieve relative errors on the order of $10^{-15}$ for moderate to large values of $n$

Detailed Methodology

Problem Formulation

Finding the zeros $t_m(u, a)$ of the reverse generalized Bessel polynomial $θ_n(z; a) = z^n y_n(z^{-1}; a)$, where $u = n + \frac{1}{2}$, within the parameter range $-\Delta_1 n + \frac{3}{2} \leq a \leq \Delta_2 n$.

Core Theoretical Framework

1. Differential Equation Transformation

Through the scaling function: $w_n^{(0)}(z; a) = 2^{-n-a+1} z^{1-n-a/2} e^{-z} θ_n(z; a)$

the problem is transformed into solving the differential equation: $\frac{d^2w}{dz^2} = \left\{1 + \frac{a-2}{z} + \frac{(2n+a)(2n+a-2)}{4z^2}\right\} w$

2. Liouville-Green Theory

Define parameters: $u = n + \frac{1}{2}, \quad α = \frac{a-2}{u}$

The differential equation is rewritten as: $\frac{d^2w}{dz^2} = \{u^2 f(α, z) + g(z)\} w$

where: $f(α, z) = \frac{(z + \frac{1}{2}α)^2 + 1 + α}{z^2}, \quad g(z) = -\frac{1}{4z^2}$

3. Turning Point Analysis

Turning points are located at: $z_{1,2}(α) = ±iσ - \frac{1}{2}α, \quad σ = \sqrt{1 + α}$

Through the Liouville variable $ζ$ and LG variable $ξ$: $\frac{2}{3}ζ^{3/2} = ξ = \int_{z_1(α)}^z f^{1/2}(α, t) dt$

Asymptotic Expansion of Zeros

The zeros satisfy the implicit equation: $Z(u, a, u^{-1}t_m(u, a)) = u^{-2/3}a_m$

where $a_m$ is the $m$-th negative zero of the Airy function $Ai(x)$.

The expansion has the form: $t_m(u, a) \sim u \sum_{s=0}^{\infty} \frac{τ_{m,s}(α)}{u^{2s}} \quad (u → ∞)$

Main Coefficient Calculations

• Leading term $τ_{m,0}$: Determined by the complex implicit equation (3.16)
• Subsequent terms: Computed through recurrence relations (3.20)-(3.23):

$τ_{m,1} = -\frac{Υ_{m,1}}{ζ'_{m,0}}$

$τ_{m,2} = -\frac{1}{2ζ'_{m,0}}\{τ_{m,1}^2 ζ''_{m,0} + 2τ_{m,1} Υ'_{m,1} + 2Υ_{m,2}\}$

Numerical Algorithm Design

1. Iteration Function

$T_n(a, z) = z - \frac{1}{\sqrt{Ω_n(a, z)}} \arctan\left(\frac{\sqrt{Ω_n(a, z)}w_n^{(0)}(z; a)}{\partial w_n^{(0)}(z; a)/\partial z}\right)$

2. Step Function

$H_n(a, z) = z + \frac{π}{\sqrt{Ω_n(a, z)}}$

3. Taylor Series Expansion

High-order derivatives are computed using recurrence relations (4.7), avoiding direct function evaluation.

Experimental Setup

Test Parameters

• Degree Range: $n = 15, 30, 50, 200, 500, 1000, 2000$
• Parameter Values: $a = 1.01, 20.2$, etc.
• Zero Indices: $m = 1, 3, 10, 15, 25$, etc.

Evaluation Metrics

• Relative Error: Comparison with Maple high-precision implementation
• Computation Time: CPU execution time
• Convergence: Iteration algorithm convergence precision set to $ε = 10^{-12}$

Implementation Details

• Asymptotic expansions use the first 5 coefficient terms
• Maple implementation uses 60-digit precision
• MATLAB algorithm tested on Dell Latitude 7410

Experimental Results

Main Results

1. High-Precision Approximation

Relative errors shown in Tables 1 and 2:

• For $a = 1.01$: Relative errors on the order of $10^{-18}$ to $10^{-15}$
• For $a = 20.2$: Relative errors on the order of $10^{-18}$ to $10^{-15}$
• All test cases achieve relative errors near or better than $10^{-15}$

2. Parameter Dependency Analysis

Figure 1 shows that relative error decreases with increasing parameter $a$, with maximum relative error less than $4×10^{-15}$.

3. Algorithm Efficiency

Computation time results in Table 3:

• $n = 30$ (15 zeros): $2.8×10^{-3}$ seconds
• $n = 2000$ (1000 zeros): $3.7×10^{-2}$ seconds
• With 67-fold increase in zero count, computation time increases only 13-fold, demonstrating algorithm efficiency

Convergence Analysis

Figure 4 demonstrates that relative error remains well-controlled even when computing large numbers of zeros, indicating excellent numerical stability of the algorithm.

Case Verification

Figure 2 presents visualization of solving the equation $F(w) = 0$, showing that the Newton method accurately locates zeros.

Related Work

Historical Development

1. Carpenter (1992) 2: Early asymptotic analysis of generalized Bessel polynomial zeros
2. Pasquini (2000) 13: Numerical methods for precise computation of generalized Bessel polynomial zeros
3. Segura (2013) 14: General computational methods for complex zeros of special functions
4. de Bruin et al. (1981) 3,4: Study of distribution domains of generalized Bessel polynomial zeros

Advantages of This Work

• Uniformity: Expansions are uniformly valid for all zeros
• Precision: Achieves machine precision level
• Efficiency: Algorithm complexity significantly superior to existing methods
• Applicability: Covers a broader parameter range

Conclusions and Discussion

Main Conclusions

1. Theoretical Contribution: Established complete asymptotic theory for zeros of reverse generalized Bessel polynomials
2. Numerical Achievement: Developed high-precision, efficient zero computation algorithm
3. Practical Value: Provides reliable computational tools for related applications

Limitations

1. Parameter Restrictions: Method applies to parameter range $-\Delta_1 n + \frac{3}{2} \leq a \leq \Delta_2 n$
2. Complexity: Asymptotic expansion derivation is mathematically involved
3. Implementation Requirements: Requires high-precision arithmetic support

Future Directions

1. Extension to more general parameter ranges
2. Investigation of similar methods for other special function zeros
3. Algorithm efficiency optimization
4. Exploration of applications in practical engineering problems

In-Depth Evaluation

Strengths

1. Theoretical Rigor: Rigorous mathematical derivation based on turning point theory and WKB methods
2. Methodological Innovation: Clever combination of asymptotic expansions with numerical algorithms
3. Excellent Precision: Achieves relative errors on the order of $10^{-15}$
4. Outstanding Efficiency: Algorithm complexity superior to traditional methods
5. Strong Practicality: Provides complete implementation and open-source code

Weaknesses

1. Theoretical Complexity: Mathematical foundations are quite sophisticated with high comprehension threshold
2. Parameter Dependence: Method validity depends on specific parameter range assumptions
3. Implementation Complexity: Requires computation of multiple complex coefficients and derivatives

Impact

1. Academic Value: Provides important theoretical tools for special function theory
2. Application Prospects: Potential applications in signal processing, control theory, etc.
3. Methodological Significance: Provides a paradigm for research on similar problems

Applicable Scenarios

1. Scientific Computing: Situations requiring high-precision computation of Bessel polynomial zeros
2. Engineering Applications: Filter design, control system analysis, etc.
3. Theoretical Research: Further investigation of special functions and asymptotic analysis

References

This paper cites 14 important references covering classical and cutting-edge work in asymptotic analysis, special functions, and numerical algorithms, providing a solid theoretical foundation for the research.

Overall Assessment: This is a high-quality mathematical analysis paper achieving high standards in both theoretical derivation and numerical implementation. The innovation of the method, precision of results, and efficiency of the algorithm are all commendable, providing valuable tools for research and applications in related fields.