Confluent hypergeometric function

Last updated
Plot of the Kummer confluent hypergeometric function 1F1(a;b;z) with a=1 and b=2 and input z2 with 1F1(1,2,z2) in the complex plane from -2-2i to 2+2i with colors created with Mathematica 13.1 Plot of the Kummer confluent hypergeometric function 1F1(a;b;z) with a=1 and b=2 and input z2 with 1F1(1,2,z2) in the complex plane from -2-2i to 2+2i with colors created with Mathematica 13.1.svg
Plot of the Kummer confluent hypergeometric function 1F1(a;b;z) with a=1 and b=2 and input z² with 1F1(1,2,z²) in the complex plane from -2-2i to 2+2i with colors created with Mathematica 13.1

In mathematics, a confluent hypergeometric function is a solution of a confluent hypergeometric equation, which is a degenerate form of a hypergeometric differential equation where two of the three regular singularities merge into an irregular singularity. The term confluent refers to the merging of singular points of families of differential equations; confluere is Latin for "to flow together". There are several common standard forms of confluent hypergeometric functions:

Contents

The Kummer functions, Whittaker functions, and Coulomb wave functions are essentially the same, and differ from each other only by elementary functions and change of variables.

Kummer's equation

Kummer's equation may be written as:

with a regular singular point at z = 0 and an irregular singular point at z = ∞. It has two (usually) linearly independent solutions M(a, b, z) and U(a, b, z).

Kummer's function of the first kind M is a generalized hypergeometric series introduced in ( Kummer 1837 ), given by:

where:

is the rising factorial. Another common notation for this solution is Φ(a, b, z). Considered as a function of a, b, or z with the other two held constant, this defines an entire function of a or z, except when b = 0, −1, −2, ... As a function of b it is analytic except for poles at the non-positive integers.

Some values of a and b yield solutions that can be expressed in terms of other known functions. See #Special cases. When a is a non-positive integer, then Kummer's function (if it is defined) is a generalized Laguerre polynomial.

Just as the confluent differential equation is a limit of the hypergeometric differential equation as the singular point at 1 is moved towards the singular point at ∞, the confluent hypergeometric function can be given as a limit of the hypergeometric function

and many of the properties of the confluent hypergeometric function are limiting cases of properties of the hypergeometric function.

Since Kummer's equation is second order there must be another, independent, solution. The indicial equation of the method of Frobenius tells us that the lowest power of a power series solution to the Kummer equation is either 0 or 1 − b. If we let w(z) be

then the differential equation gives

which, upon dividing out z1−b and simplifying, becomes

This means that z1−bM(a + 1 − b, 2 − b, z) is a solution so long as b is not an integer greater than 1, just as M(a, b, z) is a solution so long as b is not an integer less than 1. We can also use the Tricomi confluent hypergeometric function U(a, b, z) introduced by FrancescoTricomi  ( 1947 ), and sometimes denoted by Ψ(a; b; z). It is a combination of the above two solutions, defined by

Although this expression is undefined for integer b, it has the advantage that it can be extended to any integer b by continuity. Unlike Kummer's function which is an entire function of z, U(z) usually has a singularity at zero. For example, if b = 0 and a ≠ 0 then Γ(a+1)U(a, b, z) − 1 is asymptotic to az ln z as z goes to zero. But see #Special cases for some examples where it is an entire function (polynomial).

Note that the solution z1−bU(a + 1 − b, 2 − b, z) to Kummer's equation is the same as the solution U(a, b, z), see #Kummer's transformation.

For most combinations of real or complex a and b, the functions M(a, b, z) and U(a, b, z) are independent, and if b is a non-positive integer, so M(a, b, z) doesn't exist, then we may be able to use z1−bM(a+1−b, 2−b, z) as a second solution. But if a is a non-positive integer and b is not a non-positive integer, then U(z) is a multiple of M(z). In that case as well, z1−bM(a+1−b, 2−b, z) can be used as a second solution if it exists and is different. But when b is an integer greater than 1, this solution doesn't exist, and if b = 1 then it exists but is a multiple of U(a, b, z) and of M(a, b, z) In those cases a second solution exists of the following form and is valid for any real or complex a and any positive integer b except when a is a positive integer less than b:

When a = 0 we can alternatively use:

When b = 1 this is the exponential integral E1(−z).

A similar problem occurs when ab is a negative integer and b is an integer less than 1. In this case M(a, b, z) doesn't exist, and U(a, b, z) is a multiple of z1−bM(a+1−b, 2−b, z). A second solution is then of the form:

Other equations

Confluent Hypergeometric Functions can be used to solve the Extended Confluent Hypergeometric Equation whose general form is given as:

[1]

Note that for M = 0 or when the summation involves just one term, it reduces to the conventional Confluent Hypergeometric Equation.

Thus Confluent Hypergeometric Functions can be used to solve "most" second-order ordinary differential equations whose variable coefficients are all linear functions of z, because they can be transformed to the Extended Confluent Hypergeometric Equation. Consider the equation:

First we move the regular singular point to 0 by using the substitution of A + Bzz, which converts the equation to:

with new values of C, D, E, and F. Next we use the substitution:

and multiply the equation by the same factor, obtaining:

whose solution is

where w(z) is a solution to Kummer's equation with

Note that the square root may give an imaginary or complex number. If it is zero, another solution must be used, namely

where w(z) is a confluent hypergeometric limit function satisfying

As noted below, even the Bessel equation can be solved using confluent hypergeometric functions.

Integral representations

If Re b > Re a > 0, M(a, b, z) can be represented as an integral

thus M(a, a+b, it) is the characteristic function of the beta distribution. For a with positive real part U can be obtained by the Laplace integral

The integral defines a solution in the right half-plane Re z > 0.

They can also be represented as Barnes integrals

where the contour passes to one side of the poles of Γ(−s) and to the other side of the poles of Γ(a + s).

Asymptotic behavior

If a solution to Kummer's equation is asymptotic to a power of z as z → ∞, then the power must be a. This is in fact the case for Tricomi's solution U(a, b, z). Its asymptotic behavior as z → ∞ can be deduced from the integral representations. If z = xR, then making a change of variables in the integral followed by expanding the binomial series and integrating it formally term by term gives rise to an asymptotic series expansion, valid as x → ∞: [2]

where is a generalized hypergeometric series with 1 as leading term, which generally converges nowhere, but exists as a formal power series in 1/x. This asymptotic expansion is also valid for complex z instead of real x, with |arg z| < 3π/2.

The asymptotic behavior of Kummer's solution for large |z| is:

The powers of z are taken using −3π/2 < arg zπ/2. [3] The first term is not needed when Γ(ba) is finite, that is when ba is not a non-positive integer and the real part of z goes to negative infinity, whereas the second term is not needed when Γ(a) is finite, that is, when a is a not a non-positive integer and the real part of z goes to positive infinity.

There is always some solution to Kummer's equation asymptotic to ezzab as z → −∞. Usually this will be a combination of both M(a, b, z) and U(a, b, z) but can also be expressed as ez (−1)a-bU(ba, b, −z).

Relations

There are many relations between Kummer functions for various arguments and their derivatives. This section gives a few typical examples.

Contiguous relations

Given M(a, b, z), the four functions M(a ± 1, b, z), M(a, b ± 1, z) are called contiguous to M(a, b, z). The function M(a, b, z) can be written as a linear combination of any two of its contiguous functions, with rational coefficients in terms of a, b, and z. This gives (4
2
) = 6
relations, given by identifying any two lines on the right hand side of

In the notation above, M = M(a, b, z), M(a+) = M(a + 1, b, z), and so on.

Repeatedly applying these relations gives a linear relation between any three functions of the form M(a + m, b + n, z) (and their higher derivatives), where m, n are integers.

There are similar relations for U.

Kummer's transformation

Kummer's functions are also related by Kummer's transformations:

.

Multiplication theorem

The following multiplication theorems hold true:

Connection with Laguerre polynomials and similar representations

In terms of Laguerre polynomials, Kummer's functions have several expansions, for example

( Erdélyi et al. 1953 , 6.12)

or

Special cases

Functions that can be expressed as special cases of the confluent hypergeometric function include:

(a polynomial if a is a non-positive integer)
for non-positive integer n is a generalized Laguerre polynomial.
for non-positive integer n is a multiple of a generalized Laguerre polynomial, equal to when the latter exists.
when n is a positive integer is a closed form with powers of z, equal to when the latter exists.
for non-negative integer n is a Bessel polynomial (see lower down).
etc.
Using the contiguous relation we get, for example,
This identity is sometimes also referred to as Kummer's second transformation. Similarly
When a is a non-positive integer, this equals 2aθa(x/2) where θ is a Bessel polynomial.
In the second formula the function's second branch cut can be chosen by multiplying with (−1)p.

Application to continued fractions

By applying a limiting argument to Gauss's continued fraction it can be shown that [5]

and that this continued fraction converges uniformly to a meromorphic function of z in every bounded domain that does not include a pole.

Notes

  1. Campos, L.M.B.C. (2001). "On Some Solutions of the Extended Confluent Hypergeometric Differential Equation". Journal of Computational and Applied Mathematics. 137 (1): 177–200. Bibcode:2001JCoAM.137..177C. doi:10.1016/s0377-0427(00)00706-8. MR   1865885.
  2. Andrews, G.E.; Askey, R.; Roy, R. (2001). Special functions. Cambridge University Press. ISBN   978-0521789882..
  3. This is derived from Abramowitz and Stegun (see reference below), page 508, where a full asymptotic series is given. They switch the sign of the exponent in exp(iπa) in the right half-plane but this is immaterial, as the term is negligible there or else a is an integer and the sign doesn't matter.
  4. "Aspects of Multivariate Statistical Theory | Wiley". Wiley.com. Retrieved 2021-01-23.
  5. Frank, Evelyn (1956). "A new class of continued fraction expansions for the ratios of hypergeometric functions". Trans. Am. Math. Soc. 81 (2): 453–476. doi:10.1090/S0002-9947-1956-0076937-0. JSTOR   1992927. MR   0076937.

Related Research Articles

<span class="mw-page-title-main">Bessel function</span> Families of solutions to related differential equations

Bessel functions, first defined by the mathematician Daniel Bernoulli and then generalized by Friedrich Bessel, are canonical solutions y(x) of Bessel's differential equation for an arbitrary complex number , which represents the order of the Bessel function. Although and produce the same differential equation, it is conventional to define different Bessel functions for these two values in such a way that the Bessel functions are mostly smooth functions of .

<span class="mw-page-title-main">Gamma function</span> Extension of the factorial function

In mathematics, the gamma function is one commonly used extension of the factorial function to complex numbers. The gamma function is defined for all complex numbers except the non-positive integers. For every positive integer n,

Lambert <i>W</i> function Multivalued function in mathematics

In mathematics, the Lambert W function, also called the omega function or product logarithm, is a multivalued function, namely the branches of the converse relation of the function f(w) = wew, where w is any complex number and ew is the exponential function. The function is named after Johann Lambert, who considered a related problem in 1758. Building on Lambert's work, Leonhard Euler described the W function per se in 1783.

<span class="mw-page-title-main">Euler's constant</span> Constant value used in mathematics

Euler's constant is a mathematical constant, usually denoted by the lowercase Greek letter gamma, defined as the limiting difference between the harmonic series and the natural logarithm, denoted here by log:

<span class="mw-page-title-main">Legendre function</span> Solutions of Legendres differential equation

In physical science and mathematics, the Legendre functionsPλ, Qλ and associated Legendre functionsPμ
λ
, Qμ
λ
, and Legendre functions of the second kind, Qn, are all solutions of Legendre's differential equation. The Legendre polynomials and the associated Legendre polynomials are also solutions of the differential equation in special cases, which, by virtue of being polynomials, have a large number of additional properties, mathematical structure, and applications. For these polynomial solutions, see the separate Wikipedia articles.

<span class="mw-page-title-main">Generalized hypergeometric function</span> Family of power series in mathematics

In mathematics, a generalized hypergeometric series is a power series in which the ratio of successive coefficients indexed by n is a rational function of n. The series, if convergent, defines a generalized hypergeometric function, which may then be defined over a wider domain of the argument by analytic continuation. The generalized hypergeometric series is sometimes just called the hypergeometric series, though this term also sometimes just refers to the Gaussian hypergeometric series. Generalized hypergeometric functions include the (Gaussian) hypergeometric function and the confluent hypergeometric function as special cases, which in turn have many particular special functions as special cases, such as elementary functions, Bessel functions, and the classical orthogonal polynomials.

<span class="mw-page-title-main">Incomplete gamma function</span> Types of special mathematical functions

In mathematics, the upper and lower incomplete gamma functions are types of special functions which arise as solutions to various mathematical problems such as certain integrals.

<span class="mw-page-title-main">Exponential integral</span> Special function defined by an integral

In mathematics, the exponential integral Ei is a special function on the complex plane.

In mathematics, the associated Legendre polynomials are the canonical solutions of the general Legendre equation

<span class="mw-page-title-main">Parabolic cylinder function</span>

In mathematics, the parabolic cylinder functions are special functions defined as solutions to the differential equation

<span class="mw-page-title-main">Chi distribution</span> Probability distribution

In probability theory and statistics, the chi distribution is a continuous probability distribution over the non-negative real line. It is the distribution of the positive square root of a sum of squared independent Gaussian random variables. Equivalently, it is the distribution of the Euclidean distance between a multivariate Gaussian random variable and the origin. It is thus related to the chi-squared distribution by describing the distribution of the positive square roots of a variable obeying a chi-squared distribution.

<span class="mw-page-title-main">Hypergeometric function</span> Function defined by a hypergeometric series

In mathematics, the Gaussian or ordinary hypergeometric function2F1(a,b;c;z) is a special function represented by the hypergeometric series, that includes many other special functions as specific or limiting cases. It is a solution of a second-order linear ordinary differential equation (ODE). Every second-order linear ODE with three regular singular points can be transformed into this equation.

In mathematics, the local Heun function is the solution of Heun's differential equation that is holomorphic and 1 at the singular point z = 0. The local Heun function is called a Heun function, denoted Hf, if it is also regular at z = 1, and is called a Heun polynomial, denoted Hp, if it is regular at all three finite singular points z = 0, 1, a.

In mathematics, Riemann's differential equation, named after Bernhard Riemann, is a generalization of the hypergeometric differential equation, allowing the regular singular points to occur anywhere on the Riemann sphere, rather than merely at 0, 1, and . The equation is also known as the Papperitz equation.

<span class="mw-page-title-main">Whittaker function</span> In mathematics, a solution to a modified form of the confluent hypergeometric equation

In mathematics, a Whittaker function is a special solution of Whittaker's equation, a modified form of the confluent hypergeometric equation introduced by Whittaker (1903) to make the formulas involving the solutions more symmetric. More generally, Jacquet (1966, 1967) introduced Whittaker functions of reductive groups over local fields, where the functions studied by Whittaker are essentially the case where the local field is the real numbers and the group is SL2(R).

<span class="mw-page-title-main">ARGUS distribution</span>

In physics, the ARGUS distribution, named after the particle physics experiment ARGUS, is the probability distribution of the reconstructed invariant mass of a decayed particle candidate in continuum background.

<span class="mw-page-title-main">Meijer G-function</span> Generalization of the hypergeometric function

In mathematics, the G-function was introduced by Cornelis Simon Meijer as a very general function intended to include most of the known special functions as particular cases. This was not the only attempt of its kind: the generalized hypergeometric function and the MacRobert E-function had the same aim, but Meijer's G-function was able to include those as particular cases as well. The first definition was made by Meijer using a series; nowadays the accepted and more general definition is via a line integral in the complex plane, introduced in its full generality by Arthur Erdélyi in 1953.

<span class="mw-page-title-main">Struve function</span>

In mathematics, the Struve functionsHα(x), are solutions y(x) of the non-homogeneous Bessel's differential equation:

<span class="mw-page-title-main">Jacobi polynomials</span> Polynomial sequence

In mathematics, Jacobi polynomials are a class of classical orthogonal polynomials. They are orthogonal with respect to the weight on the interval . The Gegenbauer polynomials, and thus also the Legendre, Zernike and Chebyshev polynomials, are special cases of the Jacobi polynomials.

In mathematics, the Fox–Wright function (also known as Fox–Wright Psi function, not to be confused with Wright Omega function) is a generalisation of the generalised hypergeometric function pFq(z) based on ideas of Charles Fox (1928) and E. Maitland Wright (1935):

References