Biconfluent Heun Equation Overview

• Biconfluent Heun Equation is a second-order linear ODE characterized by a regular singularity at zero and a rank-two irregular singularity at infinity.
• It supports diverse analytic representations including Frobenius, incomplete Beta/Gamma, and Hermite series, enabling robust solution strategies.
• The equation underpins physical models in quantum mechanics and solvable potentials through polynomial truncation and quasi-exact solvability frameworks.

The biconfluent Heun equation (BHE) is a canonical second-order linear ordinary differential equation with one regular singularity and one irregular singularity of rank two. It emerges as a confluent limit of the general Heun equation and serves as a unifying framework for several families of special functions and exactly solvable models in mathematical physics, particularly quantum mechanics, general relativity, and integrable systems. The BHE admits rich analytic structures, including expansions in terms of incomplete Beta and Gamma functions, Hermite functions, generalized hypergeometric series, and supports quasi-exact solvability, orthogonal polynomial solutions, and connections to isomonodromic deformations and Painlevé equations.

1. Canonical Form, Singularity Structure, and Recurrence

The standard canonical form of the biconfluent Heun equation (BHE) for a function $u(z)$ is

$$u''(z) + \left(\frac{1+\alpha}{z} - \beta - 2z\right) u'(z) + \left(\frac{\gamma-\alpha-2}{z} - \frac{1}{2}\left[\delta+(1+\alpha)\beta\right] z \right) u(z) = 0,$$

where $\alpha, \beta, \gamma, \delta \in \mathbb{C}$ are independent parameters. The singularity structure consists of a regular singularity at $z=0$ (exponents $0$ and $-\alpha$) and an irregular singularity of rank two at $z=\infty$. The local solutions around $z=0$ may be constructed via a Frobenius series $u(z) = \sum_{n=0}^\infty a_n z^{n+\rho}$ with $\rho \in \{0,-\alpha\}$, leading to a three-term recurrence for the coefficients: $$a_n = -\frac{ [\beta(n-1)+\delta] a_{n-1} + [\gamma-2n+4] a_{n-2} }{ n(n-1+\alpha) },\quad (n\geq 2),$$ and $a_1 = -\delta/\alpha\; a_0$ when $\alpha \neq 0$ (Melikdzhanian et al., 2019, &&&1&&&, Sato et al., 11 May 2025). The general solution space is two-dimensional and may be described by local analytic expansions and global connection formulae involving Stokes phenomena at infinity.

2. Series Expansions: Power Series, Incomplete Beta, Gamma, and Hermite Expansions

Solutions to the BHE admit several analytic representations:

Power Series (Frobenius)

The standard power-series expansion about $z=0$ leads to entire solutions for generic parameter values, with infinite radius of convergence due to the absence of other finite singularities (Melikdzhanian et al., 2019, Choun, 2013).

Incomplete Beta and Gamma Function Expansions

A systematic approach, based on differential equations satisfied by auxiliary functions involving derivatives of the Heun function, enables the construction of expansions with terms involving incomplete Beta or Gamma functions (Ishkhanyan et al., 2014, Ishkhanyan, 2014):

• First Beta Expansion (single Beta terms):

$$u(z) = C_0 + \sum_{n=0}^\infty c_n (-z_0)^{n+\mu} B(1-\alpha, n+\mu+1; z/z_0),$$

where $B(a, b; x)$ is the incomplete Beta function, and the coefficients satisfy a four-term recurrence.

• Second Beta Expansion (pairs of Beta terms):

$$u(z) = C_0 + C_1 z + \sum_{n=0}^\infty a_n (-z_1)^{n+\nu-2} \left\{ B(1-\alpha, n+\nu+1; z/z_1) - B(1-\alpha, n+\nu-1; z/z_1) \right\},$$

with a five-term recurrence for coefficients.

• Gamma Function Expansions are obtained analogously:

$$u(z) = C_0 - \sum_{n=0}^\infty d_n \; \Gamma(n+\mu+1-\alpha, \beta z),$$

where $\Gamma(s,x)$ is the incomplete Gamma function. Four- and five-term recurrences for expansion coefficients are derived for various auxiliary equations (Ishkhanyan et al., 2014, Ishkhanyan, 2014).

Hermite and Generalized Hypergeometric Expansions

The BHE also admits expansions in terms of Hermite functions $H_\nu(\xi)$: $u(z) = \sum_{n=0}^\infty d_n H_{B+n}(\xi),$ with coefficients obeying three-term recursions (Melikdzhanian et al., 2019, Ishkhanyan et al., 2016, Ishkhanyan et al., 2017). If the Hermite series terminates (i.e., $d_{N+1}=0$), the solution is a finite sum reducible to an irreducible linear combination of four generalized hypergeometric functions of the type ${}_{p}F_{p}$ with $p$ determined by the truncation order.

3. Polynomial and Quasi-Exact Solutions: Termination, Orthogonality, and Invariant Subspaces

Polynomial (finite-sum) solutions, often termed "Heun polynomials," are obtained when series truncation conditions are met. Specifically, truncation occurs if: $\gamma-\alpha-2=2n; \quad \text{and}\ A_{n+1}=0,$ where $A_n$ are the Frobenius or associated expansion coefficients (Vieira et al., 2015, Caruso et al., 2013, Chiang et al., 2019). Such finite-dimensional solution spaces correspond to invariant subspaces of the BHE operator, often spanned by parabolic cylinder or Hermite functions, and underpin connections to quasi-exact solvability (QES) and orthogonal polynomial theory.

Orthogonality: Explicit Heun-polynomial solutions with specified parameter constraints form orthogonal systems with respect to weights such as $w(z) = z^{\alpha} e^{-z^2}$ or Bessel-type weights in limiting regimes. The biconfluent Heun equation admits analogs of both classical Hermite and Bessel polynomial orthogonalities (Chiang et al., 2013).

QES and Lie-Algebraic Structure: The BHE is gauge-equivalent to Schrödinger operators preserving finite-dimensional polynomial spaces, typically constructed via representations of the $\mathfrak{sl}_2(\mathbb{C})$ algebra. These QES instances produce all solutions in terms of generalized Hermite, Laguerre, or Bessel functions for appropriate parameter values (Chiang et al., 2013, Chiang et al., 2019).

4. Physical Applications: Quantum Models, Potentials, and Painlevé Correspondence

The BHE underlies a wide variety of physical models, including:

• Quantum Newtonian Cosmology: The radial Schrödinger equation for a galaxy in a Newtonian universe with cosmological constant reduces to a BHE, leading to quantized energy levels and polynomial wavefunctions (Vieira et al., 2015).
• Two-Electron Quantum Dots: Exact energy levels and wavefunctions are realized for two electrons interacting via Coulomb and oscillator potentials through polynomial truncation of the BHE (Caruso et al., 2013).
• Exactly/Conditionally Solvable Schrödinger Potentials: Multi-parameter Schrödinger problems, including quartic and exponential wells, can be mapped to the BHE, with physical states expressed in terms of Hermite or confluent hypergeometric functions. Conditioned parameter choices guarantee finite-sum solutions and closed-form spectra (Ishkhanyan et al., 2016, Ishkhanyan et al., 2017).
• Isomonodromic Deformations and Painlevé IV: The BHE emerges as a scalar reduction of the isomonodromic deformation equations, with degeneration loci corresponding to special solutions (rational, classical) of the Painlevé IV equation. Explicit finite-sum solutions correspond to special parabolic cylinder or Hermite solutions of PIV (Chiang et al., 2019).

5. Orthogonality, Integral Representations, and Indefinite Integrals

Integral and orthogonality structures play an important role in the analytic theory:

• Integral Representations: Under suitable parameter conditions, the BHE admits Euler-type integral representations, for example,

$$y(z) = \int_C t^{-1-\alpha} \exp(-t^2 -\beta t + zt) dt,$$

with contours $C$ chosen according to Stokes data (Sato et al., 11 May 2025, Batic et al., 2018).

• Fredholm Integral Equations: The BHE eigenproblem can be equivalently posed as a Fredholm integral equation of the second kind, with kernel constructed from a fundamental solution system and the “potential” function in Sturm-Liouville form (Chiang et al., 2013).
• Indefinite Integrals and Derivative Formulae: The BHE admits new indefinite integrals and explicit derivative formulae, often derived via adjoint or Lagrangian identities. For instance, for $\gamma = a+2$,

$$H_B'(a,B,a+2,\delta; x) = \frac{\delta + B(a+1)}{2(1+a)} H_B(a+1,B+2,a-1,3+\delta; x),\quad a\neq -1.$$

Several classes of indefinite integrals involving products of biconfluent Heun functions and elementary, error, or exponential functions can be constructed directly (Batic et al., 2018).

6. Degenerations, $q$-Difference Generalizations, and Special Reductions

The BHE is obtained as a confluent limit of more general $q$-Heun difference equations. In the limit $q\to 1$, the biconfluent $q$-Heun operator degenerates to the differential BHE, with the standard singularity structure and recurrence properties preserved (Sato et al., 11 May 2025). Further parameter tuning recovers the confluent Heun, Hermite, Laguerre, or Bessel equations as subcases, with BHE solutions reducing to those of the classical special functions.

Table: Representative Expansions for the BHE

Expansion Type	Basis Functions	Recursion (terms)
Power series (Frobenius)	$z^n$	3
Incomplete Beta	$B(a,b;z)$	4 or 5
Incomplete Gamma	$\Gamma(s, z)$	4 or 5
Hermite (parabolic cylinder)	$H_\nu(\xi)$	3
Generalized hypergeometric	${}_{p}F_{p}$	N/A, arises from finite Hermite sum

7. Further Directions and Generalizations

Auxiliary function methods applied to the BHE can generate expansions in numerous other special functions, including the Appell, Lauricella, incomplete elliptic integrals, and Fox–Wright functions, through the strategic choice of weight functions in the auxiliary ODEs (Ishkhanyan et al., 2014). The structure of recurrence relations (three-, four-, or five-term) enables both numerical implementations and analytic classification of spectral problems in associated Schrödinger or isomonodromic deformation theory.

The BHE thus serves as a central object in the modern theory of special functions, exactly solvable systems, and integrable models, bridging classical function theory with contemporary applications in mathematical and mathematical-physics research (Ishkhanyan et al., 2014, Chiang et al., 2013, Chiang et al., 2019, Sato et al., 11 May 2025, Batic et al., 2018).