What's New
About the Project
NIST
36 Integrals with Coalescing SaddlesProperties

§36.7 Zeros

Contents

§36.7(i) Fold Canonical Integral

This is the Airy function Ai9.2).

§36.7(ii) Cusp Canonical Integral

This is (36.2.4) and (36.2.1) with K=2.

The zeros in Table 36.7.1 are points in the x=(x,y) plane, where phΨ2(x) is undetermined. All zeros have y<0, and fall into two classes. Inside the cusp, that is, for x2<8|y|3/27, the zeros form pairs lying in curved rows. Close to the y-axis the approximate location of these zeros is given by

36.7.1 ym =-2π(2m+1),
m=1,2,3,,
xm,n± =2-ym(2n+12+(-1)m12±14)π,
m=1,2,3,, n=0,±1,±2,.
Table 36.7.1: Zeros of cusp diffraction catastrophe to 5D.
Zeros {xy} inside, and zeros [xy] outside, the cusp x2=827|y|3.
{±0.52768-4.37804} [±2.35218-1.74360]
{±1.41101-5.55470} {±2.36094-5.52321} [±4.42707-3.05791]
{±0.43039-6.64285} {±3.06389-6.44624} {±3.95806-6.40312} [±6.16185-4.03551]
{±1.21605-7.49906} {±2.02922-7.48629} {±4.56537-7.19629} {±5.42206-7.14718} [±7.72352-4.84817]
{±0.38488-8.31916} {±2.71193-8.22315} {±3.49286-8.20326} {±5.96669-7.85723} {±6.79538-7.80456} [±9.17308-5.55831]

More general asymptotic formulas are given in Kaminski and Paris (1999). Just outside the cusp, that is, for x2>8|y|3/27, there is a single row of zeros on each side. With n=0,1,2,, they are located approximately at

36.7.2 xn =±(827)1/2|yn|3/2(1+ξn),
yn =-(3π(8n+5)9+8ξn)1/2,

where ξn is the real solution of

36.7.3 3π(8n+5)9+8ξnξn3/2=2716(32)1/2(ln(1ξn)+3ln(32)).

For a more extensive asymptotic analysis and further tabulations, see Kaminski and Paris (1999).

§36.7(iii) Elliptic Umbilic Canonical Integral

This is (36.2.5) with (36.2.2). The zeros are lines in x=(x,y,z) space where phΨ(E)(x) is undetermined. Deep inside the bifurcation set, that is, inside the three-cusped astroid (36.4.10) and close to the part of the z-axis that is far from the origin, the zero contours form an array of rings close to the planes

36.7.4 zn=±3(14π(2n-12))1/3=3.48734(n-14)1/3,
n=1,2,3,.

Near z=zn, and for small x and y, the modulus |Ψ(E)(x)| has the symmetry of a lattice with a rhombohedral unit cell that has a mirror plane and an inverse threefold axis whose z and x repeat distances are given by

36.7.5 Δz =9π2zn2,
Δx =6πzn.

The zeros are approximated by solutions of the equation

36.7.6 exp(-2πi(z-znΔz+2xΔx))(2exp(-6πixΔx)cos(23πyΔx)+1)=3.

The rings are almost circular (radii close to (Δx)/9 and varying by less than 1%), and almost flat (deviating from the planes zn by at most (Δz)/36). Away from the z-axis and approaching the cusp lines (ribs) (36.4.11), the lattice becomes distorted and the rings are deformed, eventually joining to form “hairpins” whose arms become the pairs of zeros (36.7.1) of the cusp canonical integral. In the symmetry planes (e.g., y=0), the number of rings in the mth row, measured from the origin and before the transition to hairpins, is given by

36.7.7 nmax(m)=25613m-26952.

Outside the bifurcation set (36.4.10), each rib is flanked by a series of zero lines in the form of curly “antelope horns” related to the “outside” zeros (36.7.2) of the cusp canonical integral. There are also three sets of zero lines in the plane z=0 related by 2π/3 rotation; these are zeros of (36.2.20), whose asymptotic form in polar coordinates (x=rcosθ, y=rsinθ) is given by

36.7.8 r=3((2n-1)π4|sin(32θ)|)2/3(1+O(n-1)),
n.

§36.7(iv) Swallowtail and Hyperbolic Umbilic Canonical Integrals

The zeros of these functions are curves in x=(x,y,z) space; see Nye (2007) for Φ3 and Nye (2006) for Φ(H).