Skip to content

Hypergeometric Functions

Hypergeometric functions are a unifying language for many special functions. Legendre, Laguerre, Hermite, Bessel, and Airy functions can all be related to hypergeometric or limiting hypergeometric functions, but quantum-mechanics calculations usually become clearer after the physically adapted descendant has been named.

Use this page as an advanced reference: it explains the main equations, conventions, and termination mechanisms. Use the dedicated pages for the named functions when doing ordinary oscillator, hydrogen, angular-momentum, cylindrical, or turning-point calculations.

The rising factorial, or Pochhammer symbol, is

(a)n=a(a+1)⋯(a+n−1),(a)0=1.(a)_n = a(a+1)\cdots(a+n-1), \qquad (a)_0=1.

Equivalently,

(a)n=Γ(a+n)Γ(a)(a)_n = \frac{\Gamma(a+n)}{\Gamma(a)}

when the gamma functions are finite. The gamma-function identities behind this quotient are collected in Gamma and Beta Functions. This notation is the compact grammar of hypergeometric series.

The generalized hypergeometric function is

pFq(a1,…,apb1,…,bq;z)=∑n=0∞(a1)n⋯(ap)n(b1)n⋯(bq)nznn!.{}_{p}F_q \left( \begin{matrix} a_1,\ldots,a_p\\ b_1,\ldots,b_q \end{matrix} ;z \right) = \sum_{n=0}^{\infty} \frac{ (a_1)_n\cdots(a_p)_n }{ (b_1)_n\cdots(b_q)_n } \frac{z^n}{n!}.

The denominator parameters bjb_j must not make (bj)n(b_j)_n vanish in the denominator unless the expression is being interpreted by a limiting process.

The convergence depends on pp and qq:

  • if p≤qp\le q, the series is entire in zz;
  • if p=q+1p=q+1, the radius of convergence is 11;
  • if p>q+1p>q+1, the series generally diverges for every nonzero zz and is used, if at all, as an asymptotic series.

This convergence statement is one reason branch cuts and analytic continuation matter in hypergeometric notation.

The most important case is the Gauss hypergeometric function

2F1(a,b;c;z)=∑n=0∞(a)n(b)n(c)nznn!.{}_2F_1(a,b;c;z) = \sum_{n=0}^{\infty} \frac{(a)_n(b)_n}{(c)_n} \frac{z^n}{n!}.

The defining series converges for ∣z∣<1\lvert z\rvert\lt1 and, with additional conditions, at some points of ∣z∣=1\lvert z\rvert=1. Outside that disk, one needs analytic continuation.

It satisfies the hypergeometric differential equation

z(1−z)y′′+[c−(a+b+1)z]y′−ab y=0.z(1-z)y'' + \left[ c-(a+b+1)z \right]y' -ab\,y =0.

This equation has regular singular points at

z=0,z=1,z=∞.z=0,\qquad z=1,\qquad z=\infty.

Many second-order equations in mathematical physics can be transformed into this form when their singular-point structure matches these three points.

When

Re⁡c>Re⁡b>0,\operatorname{Re}c > \operatorname{Re}b > 0,

one has Euler’s integral representation

2F1(a,b;c;z)=Γ(c)Γ(b)Γ(c−b)∫01tb−1(1−t)c−b−1(1−zt)−a dt,{}_2F_1(a,b;c;z) = \frac{\Gamma(c)}{\Gamma(b)\Gamma(c-b)} \int_0^1 t^{b-1} (1-t)^{c-b-1} (1-zt)^{-a} \,dt,

with branch choices understood. The formula is useful for analytic continuation, estimates, and links to beta and gamma functions. Outside the stated parameter range it may still hold by analytic continuation, but the conditions should not be silently discarded.

A basic derivative identity is

ddz2F1(a,b;c;z)=abc2F1(a+1,b+1;c+1;z),\frac{d}{dz} {}_2F_1(a,b;c;z) = \frac{ab}{c} {}_2F_1(a+1,b+1;c+1;z),

when cc is not zero or a negative integer. This follows by differentiating the defining series term by term inside the disk of convergence.

If one numerator parameter is a nonpositive integer, the series terminates. For example, if

a=−N,N=0,1,2,…,a=-N, \qquad N=0,1,2,\ldots,

then

(−N)n=0for n>N.(-N)_n=0 \qquad \text{for }n>N.

Thus 2F1(−N,b;c;z){}_2F_1(-N,b;c;z) is a polynomial of degree at most NN. This termination mechanism is central in quantum mechanics. Normalizable bound states often appear when a hypergeometric series stops and becomes a polynomial times a weight factor.

When singular points of the Gauss equation coalesce, one obtains confluent hypergeometric functions. Kummer’s equation is

zy′′+(c−z)y′−ay=0.zy'' +(c-z)y' -ay =0.

One standard solution is

M(a,c,z)=1F1(a;c;z)=∑n=0∞(a)n(c)nznn!.M(a,c,z) = {}_1F_1(a;c;z) = \sum_{n=0}^{\infty} \frac{(a)_n}{(c)_n} \frac{z^n}{n!}.

Another standard solution is Tricomi’s function

U(a,c,z).U(a,c,z).

The MM solution is regular at z=0z=0 when the parameters are ordinary, while UU is often chosen for its large-zz behavior. Which solution is physical depends on the endpoint, boundary condition, and normalization.

If a=−Na=-N, then 1F1(−N;c;z){}_1F_1(-N;c;z) terminates. This is the mechanism behind the generalized Laguerre polynomials in hydrogenic radial bound states.

Hypergeometric notation exposes a family relationship, but the named functions usually carry the more useful physical normalization and orthogonality convention.

FunctionHypergeometric relationPhysics use
Legendre polynomialsPℓ(x)=2F1(−ℓ,ℓ+1;1;(1−x)/2)P_\ell(x)={}_2F_1(-\ell,\ell+1;1;(1-x)/2)central-potential angular equations
Generalized Laguerre polynomialsLn(α)(x)=(α+1)nn!1F1(−n;α+1;x)L_n^{(\alpha)}(x)=\frac{(\alpha+1)_n}{n!}{}_1F_1(-n;\alpha+1;x)hydrogen radial bound states
Bessel functionsJν(x)=(x/2)νΓ(ν+1)0F1(;ν+1;−x2/4)J_\nu(x)=\frac{(x/2)^\nu}{\Gamma(\nu+1)}{}_0F_1(; \nu+1; -x^2/4)cylindrical and spherical radial equations
Hermite polynomialsrelated to Laguerre polynomials of argument x2x^2harmonic-oscillator wavefunctions
Airy functionsexpressible through limiting hypergeometric functionssimple turning points and WKB matching

The right column is usually the best guide to which page to use first.

A common exact-solution pattern is:

  1. separate variables;
  2. factor out endpoint behavior, such as an exponential tail or near-origin power;
  3. reduce the remaining equation to a hypergeometric or confluent hypergeometric form;
  4. impose regularity and normalizability;
  5. obtain polynomial termination, which quantizes an energy or separation constant.

In the hydrogen atom, the radial equation follows this pattern. The final physical answer is written using generalized Laguerre polynomials because their orthogonality and degree labels are adapted to the radial problem.

In angular problems, the hypergeometric form is often hidden behind Legendre or associated Legendre functions. In cylindrical problems, Bessel notation is usually clearer. Near simple WKB turning points, Airy notation is the natural local language.

The defining power series is not the whole function. The branch-choice vocabulary is developed in Branch Cuts. For 2F1{}_2F_1, the standard principal branch has a cut in the zz-plane running from 11 to ∞\infty. Different analytic-continuation paths can produce different branch values.

This matters when hypergeometric functions are used for scattering, tunneling, or Green functions. The physical boundary condition may be encoded not in the local series at one point, but in which analytic continuation or branch is chosen.

  • Treating the power series for 2F1{}_2F_1 as valid beyond its disk of convergence.
  • Forgetting that denominator parameters can make the series ill-defined.
  • Using hypergeometric notation when a named orthogonal polynomial would make normalization and nodes clearer.
  • Missing polynomial termination when a numerator parameter is a nonpositive integer.
  • Assuming every hypergeometric-looking solution is normalizable.
  • Ignoring branch cuts when continuing solutions from one region to another.
  • Confusing M(a,c,z)M(a,c,z) and U(a,c,z)U(a,c,z) in confluent problems; they are adapted to different endpoint behaviors.
  • NIST Digital Library of Mathematical Functions, Chapter 15, Hypergeometric Function.
  • NIST Digital Library of Mathematical Functions, Chapter 13, Confluent Hypergeometric Functions.
  • F. W. J. Olver, D. W. Lozier, R. F. Boisvert, and C. W. Clark, eds., NIST Handbook of Mathematical Functions, Cambridge University Press, 2010.
  • G. E. Andrews, R. Askey, and R. Roy, Special Functions, Cambridge University Press, 1999.
  • L. J. Slater, Confluent Hypergeometric Functions, Cambridge University Press, 1960.
  • M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, Dover, 1965.
  • G. B. Arfken, H. J. Weber, and F. E. Harris, Mathematical Methods for Physicists, 7th ed., Academic Press, 2013.
  1. Show that 2F1(−N,b;c;z){}_2F_1(-N,b;c;z) terminates after finitely many terms when NN is a nonnegative integer.
Solution

The coefficient of znz^n contains

(−N)n=(−N)(−N+1)⋯(−N+n−1).(-N)_n = (-N)(-N+1)\cdots(-N+n-1).

When n=N+1n=N+1, one factor is zero:

(−N+N)=0.(-N+N)=0.

Therefore (−N)n=0(-N)_n=0 for all n>Nn>N, and the series terminates.

  1. Use the hypergeometric expression
Pℓ(x)=2F1(−ℓ,ℓ+1;1;(1−x)/2)P_\ell(x) = {}_2F_1(-\ell,\ell+1;1;(1-x)/2)

to check Pℓ(1)=1P_\ell(1)=1.

Solution

At x=1x=1, the argument is

z=1−x2=0.z=\frac{1-x}{2}=0.

Every hypergeometric series starts with its n=0n=0 term, which is 11. All higher powers of zz vanish at z=0z=0. Hence Pℓ(1)=1P_\ell(1)=1 in this normalization.

  1. Explain why 1F1(−n;α+1;x){}_1F_1(-n;\alpha+1;x) is a polynomial and how this relates to generalized Laguerre polynomials.
Solution

The numerator parameter is −n-n, so the same termination mechanism applies:

(−n)k=0for k>n.(-n)_k=0 \qquad \text{for }k>n.

Thus 1F1(−n;α+1;x){}_1F_1(-n;\alpha+1;x) is a polynomial of degree at most nn. The generalized Laguerre polynomial differs by a normalization factor:

Ln(α)(x)=(α+1)nn!1F1(−n;α+1;x).L_n^{(\alpha)}(x) = \frac{(\alpha+1)_n}{n!} {}_1F_1(-n;\alpha+1;x).
  1. Why is hypergeometric notation often less useful than Laguerre notation in the final hydrogen radial wavefunction?
Solution

The hypergeometric notation displays the differential-equation family and the termination condition. Laguerre notation displays the polynomial degree, parameter, orthogonality weight, and conventional normalization used in radial integrals. Once the series has terminated, the Laguerre form is better adapted to normalization, nodes, and matrix elements.