Digital Library of Mathematical Functions
About the Project
NIST
10 Bessel FunctionsBessel and Hankel Functions

§10.21 Zeros

Contents

§10.21(i) Distribution

The zeros of any cylinder function or its derivative are simple, with the possible exceptions of z=0 in the case of the functions, and z=0,\pm\nu in the case of the derivatives.

If \nu is real, then \mathop{J_{{\nu}}\/}\nolimits\!\left(z\right), {\mathop{J_{{\nu}}\/}\nolimits^{{\prime}}}\!\left(z\right), \mathop{Y_{{\nu}}\/}\nolimits\!\left(z\right), and {\mathop{Y_{{\nu}}\/}\nolimits^{{\prime}}}\!\left(z\right), each have an infinite number of positive real zeros. All of these zeros are simple, provided that \nu\geq-1 in the case of {\mathop{J_{{\nu}}\/}\nolimits^{{\prime}}}\!\left(z\right), and \nu\geq-\tfrac{1}{2} in the case of {\mathop{Y_{{\nu}}\/}\nolimits^{{\prime}}}\!\left(z\right). When all of their zeros are simple, the mth positive zeros of these functions are denoted by \mathop{j_{{\nu,m}}\/}\nolimits, \mathop{{j^{{\prime}}_{{\nu,m}}}\/}\nolimits, \mathop{y_{{\nu,m}}\/}\nolimits, and \mathop{{y^{{\prime}}_{{\nu,m}}}\/}\nolimits respectively, except that z=0 is counted as the first zero of {\mathop{J_{{0}}\/}\nolimits^{{\prime}}}\!\left(z\right). Since {\mathop{J_{{0}}\/}\nolimits^{{\prime}}}\!\left(z\right)=-\mathop{J_{{1}}\/}%
\nolimits\!\left(z\right) we have

For an extension see Pálmai and Apagyi (2011).

The positive zeros of any two real distinct cylinder functions of the same order are interlaced, as are the positive zeros of any real cylinder function \mathop{\mathscr{C}_{{\nu}}\/}\nolimits\!\left(z\right) and the contiguous function \mathop{\mathscr{C}_{{\nu+1}}\/}\nolimits\!\left(z\right). See also Elbert and Laforgia (1994).

When \nu\geq-1 the zeros of \mathop{J_{{\nu}}\/}\nolimits\!\left(z\right) are all real. If \nu<-1 and \nu is not an integer, then the number of complex zeros of \mathop{J_{{\nu}}\/}\nolimits\!\left(z\right) is 2\left\lfloor-\nu\right\rfloor. If \left\lfloor-\nu\right\rfloor is odd, then two of these zeros lie on the imaginary axis.

If \nu\geq 0, then the zeros of {\mathop{J_{{\nu}}\/}\nolimits^{{\prime}}}\!\left(z\right) are all real.

For information on the real double zeros of {\mathop{J_{{\nu}}\/}\nolimits^{{\prime}}}\!\left(z\right) and {\mathop{Y_{{\nu}}\/}\nolimits^{{\prime}}}\!\left(z\right) when \nu<-1 and \nu<-\tfrac{1}{2}, respectively, see Döring (1971) and Kerimov and Skorokhodov (1986). The latter reference also has information on double zeros of the second and third derivatives of \mathop{J_{{\nu}}\/}\nolimits\!\left(z\right) and \mathop{Y_{{\nu}}\/}\nolimits\!\left(z\right).

No two of the functions \mathop{J_{{0}}\/}\nolimits\!\left(z\right), \mathop{J_{{1}}\/}\nolimits\!\left(z\right), \mathop{J_{{2}}\/}\nolimits\!\left(z\right),\ldots, have any common zeros other than z=0; see Watson (1944, §15.28).

§10.21(ii) Analytic Properties

If \rho_{\nu} is a zero of the cylinder function

where t is a parameter, then

10.21.5{\mathop{\mathscr{C}_{{\nu}}\/}\nolimits^{{\prime}}}\!\left(\rho_{\nu}\right)=%
\mathop{\mathscr{C}_{{\nu-1}}\/}\nolimits\!\left(\rho_{\nu}\right)=-\mathop{%
\mathscr{C}_{{\nu+1}}\/}\nolimits\!\left(\rho_{\nu}\right).

If \sigma_{\nu} is a zero of {\mathop{\mathscr{C}_{{\nu}}\/}\nolimits^{{\prime}}}\!\left(z\right), then

10.21.6\mathop{\mathscr{C}_{{\nu}}\/}\nolimits\!\left(\sigma_{\nu}\right)=\frac{%
\sigma_{\nu}}{\nu}\mathop{\mathscr{C}_{{\nu-1}}\/}\nolimits\!\left(\sigma_{\nu%
}\right)=\frac{\sigma_{\nu}}{\nu}\mathop{\mathscr{C}_{{\nu+1}}\/}\nolimits\!%
\left(\sigma_{\nu}\right).

For sign properties of the forward differences that are defined by

when t=1,2,3,\ldots, and similarly for \sigma_{\nu}(t), see Lorch and Szego (1963, 1964), Lorch et al. (1970, 1972), and Muldoon (1977).

Some information on the distribution of \rho_{\nu}(t) and \sigma_{\nu}(t) for real values of \nu and t is given in Muldoon and Spigler (1984).

§10.21(iii) Infinite Products

§10.21(iv) Monotonicity Properties

In particular, \mathop{j_{{\nu,m}}\/}\nolimits, \mathop{y_{{\nu,m}}\/}\nolimits, {\mathop{j_{{\nu,m}}\/}\nolimits^{{\prime}}}, and {\mathop{y_{{\nu,m}}\/}\nolimits^{{\prime}}} are increasing functions of \nu when \nu\geq 0. It is also true that the positive zeros j^{{\prime\prime}}_{\nu} and j^{{\prime\prime\prime}}_{\nu} of {\mathop{J_{{\nu}}\/}\nolimits^{{\prime\prime}}}\!\left(x\right) and {\mathop{J_{{\nu}}\/}\nolimits^{{\prime\prime\prime}}}\!\left(x\right), respectively, are increasing functions of \nu when \nu>0, provided that in the latter case j^{{\prime\prime\prime}}_{\nu}>\sqrt{3} when 0<\nu<1.

\mathop{j_{{\nu,m}}\/}\nolimits/\nu and {\mathop{j_{{\nu,m}}\/}\nolimits^{{\prime}}}/\nu are decreasing functions of \nu when \nu>0 for m=1,2,3,\ldots.

For further monotonicity properties see Elbert (2001), Lorch (1990, 1993, 1995), Lorch and Muldoon (2008), Lorch and Szego (1990, 1995), and Muldoon (1981). For inequalities for zeros arising from monotonicity properties see Laforgia and Muldoon (1983).

§10.21(v) Inequalities

For bounds for the smallest real or purely imaginary zeros of \mathop{J_{{\nu}}\/}\nolimits\!\left(x\right) when \nu is real see Ismail and Muldoon (1995).

§10.21(vi) McMahon’s Asymptotic Expansions for Large Zeros

If \nu (\geq 0) is fixed, \mu=4\nu^{2}, and m\to\infty, then

where a=(m+\tfrac{1}{2}\nu-\tfrac{1}{4})\pi for \mathop{j_{{\nu,m}}\/}\nolimits, a=(m+\tfrac{1}{2}\nu-\tfrac{3}{4})\pi for \mathop{y_{{\nu,m}}\/}\nolimits. With a=(t+\tfrac{1}{2}\nu-\tfrac{1}{4})\pi, the right-hand side is the asymptotic expansion of \rho_{\nu}(t) for large t.

where b=(m+\tfrac{1}{2}\nu-\tfrac{3}{4})\pi for {\mathop{j_{{\nu,m}}\/}\nolimits^{{\prime}}}, b=(m+\tfrac{1}{2}\nu-\tfrac{1}{4})\pi for {\mathop{y_{{\nu,m}}\/}\nolimits^{{\prime}}}, and b=(t+\tfrac{1}{2}\nu+\tfrac{1}{4})\pi for \sigma_{\nu}(t).

For the next three terms in (10.21.19) and the next two terms in (10.21.20) see Bickley et al. (1952, p. xxxvii) or Olver (1960, pp. xvii–xviii).

For error bounds see Wong and Lang (1990), Wong (1995), and Elbert and Laforgia (2000). See also Laforgia (1979).

For the mth positive zero j^{{\prime\prime}}_{{\nu,m}} of {\mathop{J_{{\nu}}\/}\nolimits^{{\prime\prime}}}\!\left(x\right) Wong and Lang (1990) gives the corresponding expansion

where c=(m+\tfrac{1}{2}\nu-\tfrac{1}{4})\pi if 0<\nu<1, and c=(m+\tfrac{1}{2}\nu-\tfrac{5}{4})\pi if \nu>1. An error bound is included for the case \nu\geq\tfrac{3}{2}.

§10.21(vii) Asymptotic Expansions for Large Order

Let \mathop{\mathscr{C}_{{\nu}}\/}\nolimits\!\left(x\right), \rho_{\nu}(t), and \sigma_{\nu}(t) be defined as in §10.21(ii) and \mathop{M\/}\nolimits\!\left(x\right), \mathop{\theta\/}\nolimits\!\left(x\right), \mathop{N\/}\nolimits\!\left(x\right), and \mathop{\phi\/}\nolimits\!\left(x\right) denote the modulus and phase functions for the Airy functions and their derivatives as in §9.8.

As \nu\to\infty with t (>0) fixed,

where \alpha is given by

10.21.24\mathop{\theta\/}\nolimits\!\left(-2^{{\frac{1}{3}}}\alpha\right)=\pi t,

and

10.21.25
\alpha_{0}=1,
\alpha_{1}=\alpha,
\alpha_{2}=\tfrac{3}{10}\alpha^{2},
\alpha_{3}=-\tfrac{1}{350}\alpha^{3}+\tfrac{1}{70},
\alpha_{4}=-\tfrac{479}{63000}\alpha^{4}-\tfrac{1}{3150}\alpha,
\alpha_{5}=\tfrac{20231}{80\;85000}\alpha^{5}-\tfrac{551}{1\;61700}\alpha^{2},
10.21.26
\beta_{0}=1,
\beta_{1}=-\tfrac{4}{5}\alpha,
\beta_{2}=\tfrac{18}{35}\alpha^{2},
\beta_{3}=-\tfrac{88}{315}\alpha^{3}-\tfrac{11}{1575},
\beta_{4}=\tfrac{79586}{6\;06375}\alpha^{4}+\tfrac{9824}{6\;06375}\alpha.

In particular, with the notation as below,

and

Here \mathop{a_{{m}}\/}\nolimits, \mathop{b_{{m}}\/}\nolimits, \mathop{a^{{\prime}}_{{m}}\/}\nolimits, \mathop{b^{{\prime}}_{{m}}\/}\nolimits are the mth negative zeros of \mathop{\mathrm{Ai}\/}\nolimits\!\left(x\right), \mathop{\mathrm{Bi}\/}\nolimits\!\left(x\right), {\mathop{\mathrm{Ai}\/}\nolimits^{{\prime}}}\!\left(x\right), {\mathop{\mathrm{Bi}\/}\nolimits^{{\prime}}}\!\left(x\right), respectively (§9.9), \alpha_{k}, \beta_{k}, \alpha^{{\prime}}_{k}, \beta^{{\prime}}_{k} are given by (10.21.25), (10.21.26), (10.21.30), and (10.21.31), with \alpha=-2^{{-\frac{1}{3}}}\mathop{a_{{m}}\/}\nolimits in the case of \mathop{j_{{\nu,m}}\/}\nolimits and {\mathop{J_{{\nu}}\/}\nolimits^{{\prime}}}\!\left(\mathop{j_{{\nu,m}}\/}%
\nolimits\right), \alpha=-2^{{-\frac{1}{3}}}\mathop{b_{{m}}\/}\nolimits in the case of \mathop{y_{{\nu,m}}\/}\nolimits and {\mathop{Y_{{\nu}}\/}\nolimits^{{\prime}}}\!\left(\mathop{y_{{\nu,m}}\/}%
\nolimits\right), \alpha^{{\prime}}=-2^{{-\frac{1}{3}}}\mathop{a^{{\prime}}_{{m}}\/}\nolimits in the case of {\mathop{j_{{\nu,m}}\/}\nolimits^{{\prime}}} and \mathop{J_{{\nu}}\/}\nolimits\!\left({\mathop{j_{{\nu,m}}\/}\nolimits^{{\prime%
}}}\right), \alpha^{{\prime}}=-2^{{-\frac{1}{3}}}\mathop{b^{{\prime}}_{{m}}\/}\nolimits in the case of {\mathop{y_{{\nu,m}}\/}\nolimits^{{\prime}}} and \mathop{Y_{{\nu}}\/}\nolimits\!\left({\mathop{y_{{\nu,m}}\/}\nolimits^{{\prime%
}}}\right).

For error bounds for (10.21.32) see Qu and Wong (1999); for (10.21.36) and (10.21.37) see Elbert and Laforgia (1997). See also Spigler (1980).

For the first zeros rounded numerical values of the coefficients are given by

For numerical coefficients for m=2,3,4,5 see Olver (1951, Tables 3–6).

The expansions (10.21.32)–(10.21.39) become progressively weaker as m increases. The approximations that follow in §10.21(viii) do not suffer from this drawback.

§10.21(viii) Uniform Asymptotic Approximations for Large Order

As \nu\to\infty the following four approximations hold uniformly for m=1,2,\ldots:

Here \mathop{a_{{m}}\/}\nolimits and \mathop{a^{{\prime}}_{{m}}\/}\nolimits denote respectively the zeros of the Airy function \mathop{\mathrm{Ai}\/}\nolimits\!\left(z\right) and its derivative {\mathop{\mathrm{Ai}\/}\nolimits^{{\prime}}}\!\left(z\right); see §9.9. Next, z(\zeta) is the inverse of the function \zeta=\zeta(z) defined by (10.20.3). B_{0}(\zeta) and C_{0}(\zeta) are defined by (10.20.11) and (10.20.12) with k=0. Lastly,

10.21.45h(\zeta)=\left(4\zeta/(1-z^{2})\right)^{{\frac{1}{4}}}.

(Note: If the term z(\zeta)(h(\zeta))^{2}C_{0}(\zeta)/(2\zeta\nu) in (10.21.43) is omitted, then the uniform character of the error term \mathop{O\/}\nolimits(\ifrac{1}{\nu}) is destroyed.)

Corresponding uniform approximations for \mathop{y_{{\nu,m}}\/}\nolimits, {\mathop{Y_{{\nu}}\/}\nolimits^{{\prime}}}\!\left(\mathop{y_{{\nu,m}}\/}%
\nolimits\right), {\mathop{y_{{\nu,m}}\/}\nolimits^{{\prime}}}, and \mathop{Y_{{\nu}}\/}\nolimits\!\left(\mathop{{y^{{\prime}}_{{\nu,m}}}\/}%
\nolimits\right), are obtained from (10.21.41)–(10.21.44) by changing the symbols j, \mathop{J\/}\nolimits, \mathop{\mathrm{Ai}\/}\nolimits, {\mathop{\mathrm{Ai}\/}\nolimits^{{\prime}}}, \mathop{a_{{m}}\/}\nolimits, and \mathop{a^{{\prime}}_{{m}}\/}\nolimits to y, \mathop{Y\/}\nolimits, -\mathop{\mathrm{Bi}\/}\nolimits, -{\mathop{\mathrm{Bi}\/}\nolimits^{{\prime}}}, \mathop{b_{{m}}\/}\nolimits, and \mathop{b^{{\prime}}_{{m}}\/}\nolimits, respectively.

For derivations and further information, including extensions to uniform asymptotic expansions, see Olver (1954, 1960). The latter reference includes numerical tables of the first few coefficients in the uniform asymptotic expansions.

§10.21(ix) Complex Zeros

This subsection describes the distribution in \Complex of the zeros of the principal branches of the Bessel functions of the second and third kinds, and their derivatives, in the case when the order is a positive integer n. For further information, including uniform asymptotic expansions, extensions to other branches of the functions and their derivatives, and extensions to half-integer values of the order, see Olver (1954). (There is an inaccuracy in Figures 11 and 14 in this reference. Each curve that represents an infinite string of nonreal zeros should be located on the opposite side of its straight line asymptote. This inaccuracy was repeated in Abramowitz and Stegun (1964, Figures 9.5 and 9.6). See Kerimov and Skorokhodov (1985a, b) and Figures 10.21.310.21.6.)

See also Cruz and Sesma (1982), Cruz et al. (1991), Kerimov and Skorokhodov (1984c, 1987, 1988), Kokologiannaki et al. (1992), and references supplied in §10.75(iii).

Zeros of \mathop{Y_{{n}}\/}\nolimits\!\left(nz\right) and {\mathop{Y_{{n}}\/}\nolimits^{{\prime}}}\!\left(nz\right)

In Figures 10.21.1, 10.21.3, and 10.21.5 the two continuous curves that join the points \pm 1 are the boundaries of \mathbf{K}, that is, the eye-shaped domain depicted in Figure 10.20.3. These curves therefore intersect the imaginary axis at the points z=\pm ic, where c=0.66274\ldots.

The first set of zeros of the principal value of \mathop{Y_{{n}}\/}\nolimits\!\left(nz\right) are the points z=\mathop{y_{{n,m}}\/}\nolimits/n, m=1,2,\ldots, on the positive real axis (§10.21(i)). Secondly, there is a conjugate pair of infinite strings of zeros with asymptotes \imagpart{z}=\pm ia/n, where

Lastly, there are two conjugate sets, with n zeros in each set, that are asymptotically close to the boundary of \mathbf{K} as n\to\infty. Figures 10.21.1, 10.21.3, and 10.21.5 plot the actual zeros for n=1,5, and 10, respectively.

The zeros of {\mathop{Y_{{n}}\/}\nolimits^{{\prime}}}\!\left(nz\right) have a similar pattern to those of \mathop{Y_{{n}}\/}\nolimits\!\left(nz\right).

Zeros of \mathop{{H^{{(1)}}_{{n}}}\/}\nolimits\!\left(nz\right), \mathop{{H^{{(2)}}_{{n}}}\/}\nolimits\!\left(nz\right), {\mathop{{H^{{(1)}}_{{n}}}\/}\nolimits^{{\prime}}}\!\left(nz\right), {\mathop{{H^{{(2)}}_{{n}}}\/}\nolimits^{{\prime}}}\!\left(nz\right)

In Figures 10.21.2, 10.21.4, and 10.21.6 the continuous curve that joins the points \pm 1 is the lower boundary of \mathbf{K}.

The first set of zeros of the principal value of \mathop{{H^{{(1)}}_{{n}}}\/}\nolimits\!\left(nz\right) is an infinite string with asymptote \imagpart{z}=-id/n, where

The only other set comprises n zeros that are asymptotically close to the lower boundary of \mathbf{K} as n\to\infty. Figures 10.21.2, 10.21.4, and 10.21.6 plot the actual zeros for n=1,5, and 10, respectively.

The zeros of {\mathop{{H^{{(1)}}_{{n}}}\/}\nolimits^{{\prime}}}\!\left(nz\right) have a similar pattern to those of \mathop{{H^{{(1)}}_{{n}}}\/}\nolimits\!\left(nz\right). The zeros of \mathop{{H^{{(2)}}_{{n}}}\/}\nolimits\!\left(nz\right) and {\mathop{{H^{{(2)}}_{{n}}}\/}\nolimits^{{\prime}}}\!\left(nz\right) are the complex conjugates of the zeros of \mathop{{H^{{(1)}}_{{n}}}\/}\nolimits\!\left(nz\right) and {\mathop{{H^{{(1)}}_{{n}}}\/}\nolimits^{{\prime}}}\!\left(nz\right), respectively.

Zeros of \mathop{J_{{0}}\/}\nolimits\!\left(z\right)-i\mathop{J_{{1}}\/}\nolimits\!%
\left(z\right) and \mathop{J_{{n}}\/}\nolimits\!\left(z\right)-i\mathop{J_{{n+1}}\/}\nolimits\!%
\left(z\right)

For information see Synolakis (1988), MacDonald (1989, 1997), and Ikebe et al. (1993).

§10.21(x) Cross-Products

Throughout this subsection we assume \nu\geq 0, x>0, \lambda>1, and we denote 4\nu^{2} by \mu.

The zeros of the functions

and

are simple and the asymptotic expansion of the mth positive zero as m\to\infty is given by

10.21.50\alpha+\frac{p}{\alpha}+\frac{q-p^{2}}{\alpha^{3}}+\frac{r-4pq+2p^{3}}{\alpha^%
{5}}+\cdots,

where, in the case of (10.21.48),

10.21.51
\alpha=\frac{m\pi}{\lambda-1},
p=\frac{\mu-1}{8\lambda},
q=\frac{(\mu-1)(\mu-25)(\lambda^{3}-1)}{6(4\lambda)^{3}(\lambda-1)},
r=\frac{(\mu-1)(\mu^{2}-114\mu+1073)(\lambda^{5}-1)}{5(4\lambda)^{5}(\lambda-1%
)},

and, in the case of (10.21.49),

10.21.52
\alpha=\frac{(m-1)\pi}{\lambda-1},
p=\frac{\mu+3}{8\lambda},
q=\frac{(\mu^{2}+46\mu-63)(\lambda^{3}-1)}{6(4\lambda)^{3}(\lambda-1)},
r=\frac{(\mu^{3}+185\mu^{2}-2053\mu+1899)(\lambda^{5}-1)}{5(4\lambda)^{5}(%
\lambda-1)}.

The asymptotic expansion of the large positive zeros (not necessarily the mth) of the function

is given by (10.21.50), where

10.21.54
\alpha=\frac{(m-\tfrac{1}{2})\pi}{\lambda-1},
p=\frac{(\mu+3)\lambda-(\mu-1)}{8\lambda(\lambda-1)},
q=\frac{(\mu^{2}+46\mu-63)\lambda^{3}-(\mu-1)(\mu-25)}{6(4\lambda)^{3}(\lambda%
-1)},
r=\frac{(\mu^{3}+185\mu^{2}-2053\mu+1899)\lambda^{5}-(\mu-1)(\mu^{2}-114\mu+10%
73)}{5(4\lambda)^{5}(\lambda-1)}.

Higher coefficients in the asymptotic expansions in this subsection can be obtained by expressing the cross-products in terms of the modulus and phase functions (§10.18), and then reverting the asymptotic expansion for the difference of the phase functions.

For further information see Cochran (1963, 1964, 1966a, 1966b), Kalähne (1907), Martinek et al. (1966), Muldoon (1979), and Salchev and Popov (1976).

§10.21(xi) Riccati–Bessel Functions

The Riccati–Bessel functions are (\tfrac{1}{2}\pi x)^{{\frac{1}{2}}}\mathop{J_{{\nu}}\/}\nolimits\!\left(x\right) and (\tfrac{1}{2}\pi x)^{{\frac{1}{2}}}\mathop{Y_{{\nu}}\/}\nolimits\!\left(x\right). Except possibly for x=0 their zeros are the same as those of \mathop{J_{{\nu}}\/}\nolimits\!\left(x\right) and \mathop{Y_{{\nu}}\/}\nolimits\!\left(x\right), respectively. For information on the zeros of the derivatives of Riccati–Bessel functions, and also on zeros of their cross-products, see Boyer (1969). This information includes asymptotic approximations analogous to those given in §§10.21(vi), 10.21(vii), and 10.21(x).

§10.21(xii) Zeros of \alpha\mathop{J_{{\nu}}\/}\nolimits\!\left(x\right)+x{\mathop{J_{{\nu}}\/}%
\nolimits^{{\prime}}}\!\left(x\right)

For properties of the positive zeros of the function \alpha\mathop{J_{{\nu}}\/}\nolimits\!\left(x\right)+x{\mathop{J_{{\nu}}\/}%
\nolimits^{{\prime}}}\!\left(x\right), with \alpha and \nu real, see Landau (1999).

§10.21(xiii) Rayleigh Function

The Rayleigh function \mathop{\sigma_{{n}}\/}\nolimits\!\left(\nu\right) is defined by

10.21.55\mathop{\sigma_{{n}}\/}\nolimits\!\left(\nu\right)=\sum_{{m=1}}^{\infty}(j_{{%
\nu,m}})^{{-2n}},n=1,2,3,\dots.

For properties, computation, and generalizations see Kapitsa (1951b), Kerimov (1999, 2008), and Gupta and Muldoon (2000). See also Watson (1944, §§15.5, 15.51).

§10.21(xiv) \nu-Zeros

For information on zeros of Bessel and Hankel functions as functions of the order, see Cochran (1965), Cochran and Hoffspiegel (1970), Hethcote (1970), Conde and Kalla (1979), and Sandström and Ackrén (2007).