Biconfluent Heun Form

• Biconfluent Heun form is a confluent reduction of the general Heun equation that merges two regular singularities into one regular and one irregular singularity, enabling both Frobenius expansions and complex asymptotic behaviors.
• It supports diverse solution methods including power-series, polynomial truncation, special function, and integral representations governed by multi-term recurrence relations.
• Its broad applications span mathematical physics and integrable systems, influencing quasi-exact solvability, q-difference equations, and isomonodromic deformation studies.

The biconfluent Heun form is a confluent reduction of the general Heun equation in which two of the four regular singularities are merged to produce a second-order linear ordinary differential equation with one regular singularity and one irregular singularity of rank 2. This form arises in a variety of areas in mathematical physics, special function theory, and the paper of integrable systems, linking classical differential equations, $q$-difference equations, and problems involving quasi-exact solvability. The biconfluent Heun equation (BHE) supports a rich structure of polynomial, power-series, and special-function solutions governed by characteristic multi-term recurrence relations, and its discrete analog, the biconfluent $q$-Heun equation, connects directly to the geometry of confluence in both the $q$-difference and differential settings.

1. Canonical Form and Singular Structure

The classical biconfluent Heun equation is typically written in the normalized canonical form: $$x y''(x) + [1 + \alpha - \beta x - 2x^2] y'(x) + \left[ (\gamma - \alpha - 2) x - \frac{1}{2}(\delta + (1+\alpha)\beta) \right] y(x) = 0.$$ This ODE has:

• a regular singularity at $x=0$,
• an irregular singularity of rank 2 at infinity.

The singularity structure ensures that, near the origin, solutions admit Frobenius (power-series) expansions, while at infinity the dominant behavior is dictated by the irregular singularity, which may entail essential growth or decay depending on parameter values. This structure is a direct consequence of the coalescence (confluence) of singularities in the general Heun or $q$-Heun family (El-Jaick et al., 2012, Choun, 2013, Vieira et al., 2015, Chiang et al., 2013).

2. Derivation and Confluent Limit Procedures

The biconfluent form can be realized as a degeneration of more general equations, both in the differential and $q$-difference settings.

a. Differential Confluence:

One can obtain the biconfluent Heun equation from the general Heun equation by merging two regular singularities and suitable rescaling of parameters—a process that induces the transformation of singularity types (El-Jaick et al., 2012, Choun, 2013).

b. $q$-Difference Degeneration:

On the $q$-difference side, the biconfluent $q$-Heun equation arises through the degeneration: $$f(x) - (B_2 x^2 + B_1 x + B_0) f(x) + (C_2 x^2 + C_1 x + C_0) f(x/q) = 0,$$ where the collapse of certain highest-order coefficients in $x$ signals the confluencing of singularities. Under the limiting process $q = 1 + \varepsilon$ with $\varepsilon \rightarrow 0$, Taylor expansions and suitable scalings yield a biconfluent differential equation in the continuum limit: $$x^2 g''(x) + \{ b_{2,1} x^2 + b_{1,1} x + (b_{0,1} + b_0) \} g'(x) + (b_{2,0} x + b_{1,0}) g(x) = 0,$$ which is equivalent, under gauge transformations, to the standard biconfluent Heun equation (Sato et al., 11 May 2025).

3. Solution Structures and Recurrence Relations

a. Power Series and Polynomial Solutions

A generic solution near the origin can be expressed as a Frobenius series: $y(x) = \sum_{n=0}^{\infty} c_n x^n,$ with coefficients $c_n$ obeying a three-term recurrence relation: $R_n c_n + Q_{n-1} c_{n-1} + P_{n-2} c_{n-2} = 0,$ where $R_n$, $Q_n$, $P_n$ are functions of the equation parameters. For specific parameter choices (e.g., $q_1 = 2N$ in the scaled Hermite expansion), the series truncates, yielding polynomial (Heun polynomial) solutions associated with quantization conditions and spectral problems in quantum mechanics (Vieira et al., 2015, Melikdzhanian et al., 2019).

b. Multi-Term Recurrence in Special Function Expansions

Expansions in terms of incomplete Beta, incomplete Gamma, Hermite, or parabolic cylinder functions correspond to solutions where the expansion coefficients satisfy four- or five-term recurrences:

• Incomplete Beta/Gamma Expansions:

For example, by extracting a derivative and integrating termwise, one gets expansions such as

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

with $c_n$ satisfying a four-term recurrence; or, in Gamma expansions,

$$u(z) = C_0 + e^{-sz} \sum_{n=0}^{\infty} c_n\, \Gamma(1+n+\lambda-\alpha; sz),$$

with $c_n$ generally governed by a five-term relation (Ishkhanyan et al., 2014, Ishkhanyan, 2014).

• Hermite Function Expansions:

Series of the form

$$y(x) = \sum_{n=0}^\infty d_n H_{P_0 + q_1/2 + n}(s x - P_1/2)$$

where $H_\nu$ is the Hermite function, and $d_n$ satisfies a modified three-term relation (Ishkhanyan et al., 2016, Melikdzhanian et al., 2019).

• Parabolic Cylinder Function Expansions:

In degenerate cases, especially for isomonodromic deformations/Painlevé IV, solutions are finite sums:

$y(x) = e^{x^2/4} \sum_{k=0}^N A_k D_{e/2-k}(x),$

with $D_\nu$ the parabolic cylinder function and $A_k$ governed by a three-term recurrence. Such solutions terminate when parameters are tuned according to monodromy and apparent singularity constraints (Chiang et al., 2019).

4. Integral Representations and Generating Functions

a. Integral Representations:

A systematic method to derive integral forms starts from the power series and exploits the representation of beta/gamma factors in the coefficients. For the biconfluent Heun function (as a special case of the grand confluent hypergeometric function), the solution admits integral representations of the type: $$y(x) = y_0(x) + \sum_{n=1}^\infty \left\{ \prod_{j=1}^{n} \int_0^1 dt_j \int_0^1 du_j\, W_j(t_j, u_j) \right\} \phi_n(x),$$ where $W_j$ are kernel functions arising from beta-function integral identities (Choun, 2013).

b. Path-Sum Approach:

An alternative robust approach recasts the second-order ODE as a matrix-valued first-order system, with the solution given by the path-ordered exponential of the coefficient matrix. This yields an explicit Neumann/Vottera series involving only rational and exponential elementary integrands (Giscard et al., 2020): $$H_B(z) = H_0 + H_0 \int_{z_0}^{z} G_1(\zeta, z_0) d\zeta + (H'_0 - H_0) \left( e^{z-z_0} - 1 + \int_{z_0}^{z} (e^{z-\zeta}-1) G_2(\zeta, z_0)d\zeta \right),$$ where $G_{1,2}$ are defined via absolutely convergent Volterra iterated convolution series.

c. Generating Functions:

Based on the three-term recurrence, generating functions can be built for the sequence of (terminated) polynomial solutions, often facilitating orthogonality and further analytic properties.

5. Algebraic Structure and Quasi-Exact Solvability

The biconfluent Heun operator admits an embedding into the enveloping algebra of $\mathfrak{sl}_2(\mathbb{C})$ (or $su(1,1)$ in the real case). This algebraic realization explains both the existence of quasi-exactly solvable models—where a finite portion of the spectrum and eigenfunctions can be constructed exactly—and the analytic termination of power or special function series (Chiang et al., 2013, Kar, 2020).

• su(1,1) and Quasi-Polynomial Solutions:

The biconfluent Heun class supports solutions of the form $z^\alpha P_N(z)$, where $P_N$ is a polynomial, whenever underlying parameters satisfy algebraicity conditions arising from possible finite-dimensional invariant subspaces under the algebra action.

• Connection to QES Potentials:

In quantum systems, such as two-electron quantum dot models or certain Schrödinger and Klein-Gordon equations with polynomial and inverse-polynomial potentials, the exact solvability of a subset of the spectrum is governed by the possibility of terminating the corresponding Heun series (Caruso et al., 2013, Vieira et al., 2015, Arda et al., 2016).

6. Special Function Expansions and Orthogonality

a. Hypergeometric and Generalized-Hypergeometric Expansions:

In many cases, the biconfluent Heun functions can be written as irreducible linear combinations of four generalized hypergeometric functions (not reducible to a single function except in trivial cases). These combinations capture the structure of the solution space, which is generally richer than, for example, the classical confluent hypergeometric case (Melikdzhanian et al., 2019).

b. Orthogonality and Polynomial Analogs:

• For periodic biconfluent Heun equations (PBHE), eigen-solutions can display both single and double orthogonality, paralleling classical results for Lamé equations.
• The Bessel polynomial analogy: The biconfluent Heun equation and Bessel equations both have a regular singularity at the origin and an irregular singularity at infinity, with Bessel polynomials providing a structural model for some cases (Chiang et al., 2013).

7. Connections to Painlevé Equations and Isomonodromic Deformations

There is a deep correspondence between the degenerate monodromy/Stokes data of biconfluent Heun operators and special-function solutions of the Painlevé IV equation. Finite-sum eigen-solutions (in parabolic cylinder functions) of the BHE arise precisely for parameter choices corresponding to solvable or reducible monodromy. Infinite series in these bases converge for generic parameters and terminate only for the special cases matching the Painlevé IV special function regime (Chiang et al., 2019).

8. Discrete Analogs: The Biconfluent $q$-Heun Equation

The biconfluent $q$-Heun equation

$$f(x) - (B_2 x^2 + B_1 x + B_0) f(x) + (C_2 x^2 + C_1 x + C_0) f(x/q) = 0$$

is a discrete counterpart to the classical BHE, arising as a degeneration of the general $q$-Heun equation via coalescence of singularities. Its limit as $q \to 1$ recovers the continuous biconfluent Heun differential equation, mirroring the confluence in the continuous case. Study of these equations reveals deep links to $q$-Painlevé systems, Newton polygon stratification, and the transition between discrete and continuous special function theory (Sato et al., 11 May 2025).

---

This comprehensive structure frames the biconfluent Heun form as a central confluent reduction of the Heun hierarchy, unifying power series, polynomial, special function, and integral solution techniques with direct implications in spectral problems, quasi-exact solvability, isomonodromy theory, and discrete–continuous analytic correspondences.