Oblate spheroidal wave function

In applied mathematics, oblate spheroidal wave functions (like also prolate spheroidal wave functions and other related functions[1]) are involved in the solution of the Helmholtz equation in oblate spheroidal coordinates. When solving this equation, ${\displaystyle \Delta \Phi +k^{2}\Phi =0}$, by the method of separation of variables, ${\displaystyle (\xi ,\eta ,\varphi )}$, with:

${\displaystyle \ z=(d/2)\xi \eta ,}$

${\displaystyle \ x=(d/2){\sqrt {(\xi ^{2}+1)(1-\eta ^{2})}}\cos \varphi ,}$

${\displaystyle \ y=(d/2){\sqrt {(\xi ^{2}+1)(1-\eta ^{2})}}\sin \varphi ,}$

${\displaystyle \ \xi \geq 0{\text{ and }}|\eta |\leq 1.}$

the solution ${\displaystyle \Phi (\xi ,\eta ,\varphi )}$ can be written as the product of a radial spheroidal wave function ${\displaystyle R_{mn}(-ic,i\xi )}$ and an angular spheroidal wave function ${\displaystyle S_{mn}(-ic,\eta )}$ by ${\displaystyle e^{im\varphi }}$. Here ${\displaystyle c=kd/2}$, with ${\displaystyle d}$ being the interfocal length of the elliptical cross section of the oblate spheroid.

The radial wave function ${\displaystyle R_{mn}(-ic,i\xi )}$ satisfies the linear ordinary differential equation:

${\displaystyle \ (\xi ^{2}+1){\frac {d^{2}R_{mn}(-ic,i\xi )}{d\xi ^{2}}}+2\xi {\frac {dR_{mn}(-ic,i\xi )}{d\xi }}-\left(\lambda _{mn}(c)-c^{2}\xi ^{2}-{\frac {m^{2}}{\xi ^{2}+1}}\right){R_{mn}(-ic,i\xi )}=0}$.

The angular wave function satisfies the differential equation:

${\displaystyle \ (1-\eta ^{2}){\frac {d^{2}S_{mn}(-ic,\eta )}{d\eta ^{2}}}-2\eta {\frac {dS_{mn}(-ic,\eta )}{d\eta }}+\left(\lambda _{mn}(c)+c^{2}\eta ^{2}-{\frac {m^{2}}{1-\eta ^{2}}}\right){S_{mn}(-ic,\eta )}=0}$.

It is the same differential equation as in the case of the radial wave function. However, the range of the radial coordinate ${\displaystyle \xi }$ is different from that of the angular coordinate ${\displaystyle \eta }$.

The eigenvalue ${\displaystyle \lambda _{mn}(-ic)}$ of this Sturm–Liouville problem is fixed by the requirement that ${\displaystyle {S_{mn}(-ic,\eta )}}$ be finite for ${\displaystyle |\eta |=1}$.

For ${\displaystyle c=0}$ these two differential equations reduce to the equations satisfied by the associated Legendre polynomials. For ${\displaystyle c\neq 0}$, the angular spheroidal wave functions can be expanded as a series of Legendre functions. Such expansions have been considered by Müller.[2]

The differential equations given above for the oblate radial and angular wave functions can be obtained from the corresponding equations for the prolate spheroidal wave functions by the substitution of ${\displaystyle -ic}$ for ${\displaystyle c}$ and ${\displaystyle i\xi }$ for ${\displaystyle \xi }$. The notation for the oblate spheroidal functions reflects this relationship.

There are different normalization schemes for spheroidal functions. A table of the different schemes can be found in Abramowitz and Stegun.[3] Abramowitz and Stegun (and the present article) follow the notation of Flammer.[4]

Originally, the spheroidal wave functions were introduced by C. Niven,[5] which lead to a Helmholtz equation in spheroidal coordinates. Monographs tying together many aspects of the theory of spheroidal wave functions were written by Strutt,[6] Stratton et al.,[7] Meixner and Schafke,[8] and Flammer.[4]

Flammer[4] provided a thorough discussion of the calculation of the eigenvalues, angular wavefunctions, and radial wavefunctions for both the oblate and the prolate case. Computer programs for this purpose have been developed by many, including Van Buren et al.,[9] King and Van Buren,[10] Baier et al.,[11] Zhang and Jin,[12] and Thompson.[13] Van Buren has recently developed new methods for calculating oblate spheroidal wave functions that extend the ability to obtain numerical values to extremely wide parameter ranges. These results are based on earlier work on prolate spheroidal wave functions.[14][15] Fortran source code that combines the new results with traditional methods is available at http://www.mathieuandspheroidalwavefunctions.com.

Tables of numerical values of oblate spheroidal wave functions are given in Flammer,[4] Hanish et al.,[16][17][18] and Van Buren et al.[19]

Asymptotic expansions of angular oblate spheroidal wave functions for large values of ${\displaystyle c}$ have been derived by Müller.,[20] also similarly for prolate spheroidal wave functions.[21]

The Digital Library of Mathematical Functions http://dlmf.nist.gov provided by NIST is an excellent resource for spheroidal wave functions.