Parabolic cylinder function explained

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

This equation is found when the technique of separation of variables is used on Laplace's equation when expressed in parabolic cylindrical coordinates.

The above equation may be brought into two distinct forms (A) and (B) by completing the square and rescaling, called H. F. Weber's equations:and

If f(a,z) is a solution, then so aref(a,-z), f(-a,iz)\textf(-a,-iz).

If f(a,z)\, is a solution of equation, then f(-ia,ze^) is a solution of, and, by symmetry,f(-ia,-ze^), f(ia,-ze^)\textf(ia,ze^)are also solutions of .

Solutions

There are independent even and odd solutions of the form . These are given by (following the notation of Abramowitz and Stegun (1965)):y_1(a;z) = \exp(-z^2/4) \;_1F_1 \left(\tfrac12a+\tfrac14; \;\tfrac12\; ; \; \frac\right)\,\,\,\,\,\, (\mathrm)and y_2(a;z) = z\exp(-z^2/4) \;_1F_1 \left(\tfrac12a+\tfrac34; \;\tfrac32\; ; \; \frac\right)\,\,\,\,\,\, (\mathrm)where

1F1(a;b;z)=M(a;b;z)

is the confluent hypergeometric function.

Other pairs of independent solutions may be formed from linear combinations of the above solutions. One such pair is based upon their behavior at infinity:U(a,z)=\frac\left[\cos(\xi\pi)\Gamma(1/2-\xi)\,y_1(a,z) -\sqrt{2}\sin(\xi\pi)\Gamma(1-\xi)\,y_2(a,z) \right]V(a,z)=\frac\left[\sin(\xi\pi)\Gamma(1/2-\xi)\,y_1(a,z) +\sqrt{2}\cos(\xi\pi)\Gamma(1-\xi)\,y_2(a,z) \right]where \xi = \fraca+\frac .

The function approaches zero for large values of and, while diverges for large values of positive real . \lim_U(a,z)/\left(e^z^\right)=1\,\,\,\,(\text\,\left|\arg(z)\right|<\pi/2) and \lim_V(a,z)/\left(\sqrte^z^\right)=1\,\,\,\,(\text\,\arg(z)=0) .

For half-integer values of a, these (that is, U and V) can be re-expressed in terms of Hermite polynomials; alternatively, they can also be expressed in terms of Bessel functions.

The functions U and V can also be related to the functions (a notation dating back to Whittaker (1902)) that are themselves sometimes called parabolic cylinder functions:\beginU(a,x) &= D_(x), \\V(a,x) &= \frac[\sin(\pi a) D_{-a-\tfrac12}(x)+D_{-a-\tfrac12}(-x)] .\end

Function was introduced by Whittaker and Watson as a solution of eq.~ with \tilde a=-\frac14, \tilde b=0, \tilde c=a+\frac12 bounded at

+infty

. It can be expressed in terms of confluent hypergeometric functions as
D
a(z)=1
\sqrt{\pi
}.Power series for this function have been obtained by Abadir (1993).

Parabolic Cylinder U(a,z) function

Integral representation

Integrals along the real line,U(a,z)=\frac\int_0^\infty e^t^e^dt\,,\; \Re a>-\frac12 \;,U(a,z)=\sqrte^\int_0^\infty \cos\left(zt+\fraca+\frac\right)t^e^dt\,,\; \Re a<\frac12 \;.The fact that these integrals are solutions to equation  can be easily checked by direct substitution.

Derivative

Differentiating the integrals with respect to

z

gives two expressions for

U'(a,z)

,U'(a,z)=-\fracU(a,z)-\frac\int_0^\infty e^t^e^dt=-\fracU(a,z)-\left(a+\frac12\right)U(a+1,z) \;,U'(a,z)=\fracU(a,z)-\sqrte^\int_0^\infty \sin\left(zt+\fraca+\frac\right)t^e^dt= \fracU(a,z)-U(a-1,z) \;.Adding the two gives another expression for the derivative,2U'(a,z) = -\left(a+\frac12\right)U(a+1,z)-U(a-1,z)\;.

Recurrence relation

Subtracting the first two expressions for the derivative gives the recurrence relation,zU(a,z) = U(a-1,z) - \left(a+\frac12\right)U(a+1,z) \;.

Asymptotic expansion

Expandinge^=1-\frac12 t^2+\frac18 t^4 - \dots \;in the integrand of the integral representationgives the asymptotic expansion of

U(a,z)

,U(a,z) = e^z^\left(1- \frac\frac+ \frac\frac- \dots\right) .

Power series

Expanding the integral representation in powers of

z

givesU(a,z)=\frac-\fracz+\fracz^2 - \dots \;.

Values at z=0

From the power series one immediately getsU(a,0)=\frac \;,U'(a,0)=-\frac \;.

Parabolic cylinder Dν(z) function

Parabolic cylinder function

D\nu(z)

is the solution to the Weber differential equation,u+\left(\nu+\frac12-\frac z^2 \right)u=0 \,,that is regular at

\Rez\to+infty

with the asymptoticsD_\nu(z) \to e^z^\nu \,.It is thus given as

D\nu(z)=U(-\nu-1/2,z)

and its properties then directly follow from those of the

U

-function.

Integral representation

D_\nu(z)=\frac\int_0^\infty e^ t^ e^dt\,,\; \Re \nu < 0 \,,\; \Re z > 0\;,D_\nu(z)=\sqrte^\int_0^\infty \cos\left(zt-\nu \frac\right)t^e^dt\,,\; \Re \nu > -1 \;.

Asymptotic expansion

D_\nu(z) = e^z^\left(1- \frac\frac+ \frac\frac- \dots\right)\,,\; \Re z \to +\infty .If

\nu

is a non-negative integer this series terminates and turns into a polynomial, namely the Hermite polynomial,D_n(z) = e^\;2^H_n\left(\frac\right)\,, n=0,1,2,\dots \;.

Connection with quantum harmonic oscillator

Parabolic cylinder

D\nu(z)

function appears naturally in the Schrödinger equation for the one-dimensional quantum harmonic oscillator (a quantum particle in the oscillator potential),\left[-\frac{\hbar^2}{2m}\frac{\partial^2}{\partial x^2}+\frac12 m \omega^2 x^2 \right]\psi(x)=E\psi(x) \;,where

\hbar

is the reduced Planck constant,

m

is the mass of the particle,

x

is the coordinate of the particle,

\omega

is the frequency of the oscillator,

E

is the energy,and

\psi(x)

is the particle's wave-function. Indeed introducing the new quantitiesz=\frac \,,\; \nu=\frac-\frac12 \,,\; b_o=\sqrt \,,turns the above equation into the Weber's equation for the function

u(z)=\psi(zbo)

,u+\left(\nu+\frac12-\frac z^2 \right)u=0 \,.

This article is licensed under the GNU Free Documentation License. It uses material from the Wikipedia article "Parabolic cylinder function".

Except where otherwise indicated, Everything.Explained.Today is © Copyright 2009-2025, A B Cryer, All Rights Reserved. Cookie policy.