1. Introduction
Plancherel–Rotach asymptotic expansions for the orthogonal polynomials as the polynomial degree tends to infinity have been studied extensively in the literature. There are many standard asymptotic techniques including the steepest descent method for integrals [Reference Wong12], the WKB method for differential equations [Reference Olver8], the Deift–Zhou method for Riemann–Hilbert problems [Reference Deift, Kriecherbauer, McLaughlin, Venakides and Zhou4, Reference Deift and Zhou5], asymptotic theory of difference equations [Reference Wong13] and Darboux's method [Reference Wong and Zhao14]. The error term for the truncated asymptotic expansion is usually of the same order of magnitude as the first neglected term. In other words, the ratio of the error term and the first omitted term is bounded by a constant independent of the polynomial degree. The existence of such a bound can be proved by standard ‘soft’ analysis. However, as far as we know, the quantitative information for this error bound is unknown even for the classical orthogonal polynomials. One may expect a rather complicated ‘hard’ analysis on finding explicit and computable expressions for error bounds. Various reasons why we need error bounds for asymptotic expansions are discussed in Ref. [Reference Olver9]. As mentioned in Ref. [Reference Wong11], the upper bounds of the error term obtained from standard asymptotic techniques are difficult to compute, and thus may not be realistic. In this paper, we will derive explicit and computable error bounds for the asymptotic expansions of the Hermite polynomials.
The main challenge is to find a convenient integral representation of the error term on an appropriate contour so that the error bound is computable. Berry and Howls [Reference Berry and Howls2] proposed the so-called ‘adjacent saddles’ method which was further developed in [Reference Bennett, Howls, Nemes and Olde Daalhuis1, Reference Boyd3]. The key idea of their method is to express the error term as an integral over the ‘adjacent contours’ passing through the ‘adjacent saddles’. However, the technique of ‘adjacent saddles’ cannot be applied directly to the Hermite polynomials and many other orthogonal polynomials because the integrand (or the phase function) of the integral representation for the Hermite polynomials has a branch point in the complex plane. It is thus difficult to express the error term as an integral over the ‘adjacent contours’. Moreover, for the turning point case, when two saddle points coincide with each other, there does not exist any ‘adjacent saddle’. To resolve these two difficulties, we introduce a new ‘branch cut’ technique which deforms the contour of integration for the error term to the branch cut of the phase function. More precisely, the contour is defined as the limit of the curves approaching one side of the branch cut. By virtue of the ‘branch cut’ technique, we are able to find simple error bounds in terms of elementary functions.
The Hermite polynomials can be expressed as contour integrals:
where $\Gamma$ is a counter-clockwisely oriented contour encircling the negative real line (cf. [7, § 18.10(iii)]). To be more specific, we may choose
where $\delta >0$ is any fixed positive number; see figure 1.
After introducing the Plancherel–Rotach scale, we rewrite the integral representation as
where $N=2n+1$. By symmetry, we may assume that $x\ge 0$. We shall denote the phase function in the above integral by
The zeros of $f'(t;x)$ are called saddle points and can be calculated as
Note that the two saddle points coincide ($t_+=t_-$) at the turning point $x=1$. For simplicity, we shall drop the dependence of $f$ on $x$. There are three cases to be considered.
Case I: $x=\cosh \beta >1$ with $\beta >0$. This is called the outer interval.
Case II: $x=\cos \alpha \in [0,\,1)$ with $\alpha \in (0,\,\pi /2]$. This is called the oscillatory interval.
Case III: $x=1$. This is called the turning point.
We will consider these three cases separately. We mention here that an asymptotic expansion with error bounds for the slightly differently scaled $H_n(\sqrt {N+1}\cos \alpha )$ was given earlier by van Veen [Reference van Veen10].
The rest of the paper is organized as follows. In § 2–4, we derive the asymptotic expansions with error bounds for Case I, II and III, respectively. The main results are stated at the end of each section. In § 5, we demonstrate the accuracy of the error bounds through numerical examples. Recursive procedures for the computation of the coefficients appearing in the asymptotic expansions are given in appendix A.
2. Case I: $x=\cosh \beta >1$ with $\beta >0$
The saddle points are $t_\pm =e^{\pm \beta }/2$ with $0< t_-< t_+$. Moreover, $f''(t_-)=2(1-e^{2\beta })<0$ and $f''(t_+)=2(1-e^{-2\beta })>0$. We shall deform the contour of integration $\Gamma$ to the path of steepest descent passing through the saddle point $t_-$. To describe this steepest descent contour, we shall introduce the analytic function
where the branch of the square root function is chosen so that
for $z$ in a small complex neighbourhood of $t_-$. Let $z=re^{i\theta }$ with $r>0$ and $\theta \in (-\pi,\,\pi )$. We investigate the equation ${\rm Im} [f(z)-f(t_-)] = {\rm Im}\ f(z)=0$; namely,
Clearly, $\theta =0$ satisfies the above equation. By symmetry, we can restrict $\theta \in (0,\,\pi )$ and find two solutions
It is readily seen that $r_\pm (\theta )\to t_\pm =e^{\pm \beta }/2$ as $\theta \to 0$. For consistency, we define $r_\pm (0)=t_\pm$. Consequently, we can regard $r_-(\theta )$ as a symmetric function on $(-\pi,\,\pi )$ and $r_+(\theta )$ a symmetric function on $(-\pi /2,\,\pi /2)$. Define two contours
It follows that the solution of ${\rm Im}\ f(z)=0$ is the union $\Gamma _-\cup \Gamma _+\cup \mathbb {R}^{+}$, where $\mathbb {R}^{+}$ is the positive real line. The contours $\Gamma _\pm$ are depicted in figure 2.
Since $f'(z)$ has no zeros other than $t_\pm$, it follows from the asymptotic behaviour $f(z)\sim z^{2}$ as $z\to \infty$ that $\Gamma _-$ is a steepest descent contour passing through $t_-$, while $\Gamma _+$ is a steepest ascent contour passing through $t_+$. Moreover, the solutions of $f(z)-f(t_-)=0$ are real and positive. Consider $f(t)=t^{2}-2t\cosh \beta +(\ln t)/2$ for $t>0$. Since $f(t)\to -\infty$ as $t\to 0^{+}$ and $f(t)\to \infty$ as $t\to \infty$, we obtain from $f'(t_\pm )=0$ that the equation $f(t)-f(t_-)=0$ has exactly two positive roots: one is $t_-$; and the other, denoted by $t_c$, is larger than $t_+$. Hence, we can choose $[t_c,\,\infty )$ as a branch cut and extend the definition of $w(z)$ on $\mathbb {C}\setminus (-\infty,\,0]\cup [t_c,\,\infty )$ by analytic continuation.
Now, we deform the contour of integration from $\Gamma$ to the steepest descent contour $\Gamma _-$ with counter-clockwise orientation and then make a change of variable $u=w(t)$:
Recall that $\Gamma _+$ is the steepest ascent contour passing through $t_+$. We also denote by $\Gamma _0$ a clockwise oriented contour encircling $(-\infty,\,0]$; namely,
where $\delta >0$ is a small positive number such that $\delta < t_c$. In fact, $\Gamma _0$ is identical to $\Gamma$ but has the opposite orientation. Since $|w(z)|/|z|\to 1$ as $z\to \infty$, we have from the Cauchy integral formula
for any $u\in \mathbb {R}$. Substituting this expression into (2.4) gives a double integral representation
For any $p\ge 0$, it follows from
and
that
where
We show in appendix A.1 that the coefficients $A_j(\coth \beta )$ are polynomials in $\coth \beta$ of degree $3j$ with rational coefficients, and can be calculated using a recursive formula.
To find an upper bound for the remainder $\varepsilon _p(N,\,\beta )$, we need to estimate the integral
If there exists $C_p>0$ such that $|c_p(u)|\le C_p$ for all $u\in \mathbb {R}$, then we obtain from (2.5) with $k=2p+1$ and (2.8) that
First, we estimate the integral
Recall that the contour $\Gamma _0$ encircles the negative real line having a distance $\delta >0$ from it. By analyticity of the integrand, the value $I_0$ does not change if we let $\delta \to 0^{+}$. Hence,
where $w(se^{\pm i\pi })$ is defined as the limit of $w(se^{i\theta })$ as $\theta \to \pm \pi$. Note from (1.2) and (2.1) that
Thus, we obtain
where the last equality defines $R(s)$, the modulus of $w(se^{\pm i\pi })$. Let $\theta _\pm (s)$ be the phase of $w(se^{\pm i\pi })$. We then have
which implies that
Assume that $p\geq 1$. Substituting this inequality into (2.11) yields
To estimate the integral on the right-hand side of the above inequality, we note from (2.12) that $R(s)\ge \sqrt {\pi /2}$ for $s\in [0,\,1]$ and $R(s)\ge \sqrt {2s\cosh \beta }$ for $s\ge 1>t_-$. Thus,
where we have assumed $p\ge 2$ to ensure convergence of the second integral.
It remains to estimate the integral
For $z\in \Gamma _+$, we have ${\rm Im}\ f(z)=0$ and ${\rm Re} [f(z)-f(t_-)]< f(t_+)-f(t_-)<0$. It then follows from (2.1) that $w(z)$ is purely imaginary and $|w(z)-u|\ge |w(z)|$ for $z\in \Gamma _+$ and $u\in \mathbb {R}$. Now, we parameterize $\Gamma _+$ as $z=z_+(\theta )=r_+(\theta )e^{i\theta }$ with $\theta \in (-\pi /2,\,\pi /2)$, where $r_+(\theta )$ is given in (2.3). Since $dz/d\theta =[r_+'(\theta )+ir_+(\theta )]e^{i\theta }$ and $|w(z_+(\theta ))|^{2}=|f(z_+(\theta ))-f(t_-)|=f(t_-)-f(z_+(\theta ))$, we deduce
As both $\theta /\tan \theta$ and $\cos \theta$ are decreasing functions for $\theta \in (0,\,\pi /2)$, it follows from (2.3) that $r_+(\theta )$ is an increasing function for $\theta \in (0,\,\pi /2)$; namely, $r_+'(\theta )\ge 0$. To estimate $r_+'(\theta )$, we apply an implicit differentiation on (2.2) and obtain
A simple calculation together with (2.3) gives
and
Consequently,
and
Let $\theta _0=\arccos (1/4)$. If $\theta \in (0,\,\theta _0)$, then $\cos \theta >1/4$, $f(t_-)-f(z_+(\theta ))>f(t_-)-f(t_+)=\sinh (2\beta )/2-\beta$ and $r_+(\theta )<4\cosh \beta$. Hence,
If $\theta \in (\theta _0,\,\pi /2)$, then $r_+(\theta )>2e^{\beta }=4t_+>2\cosh \beta >2$. This together with (2.14) implies that
and
Coupling the above two inequalities gives
Therefore,
where we have assumed $p\ge 1$. Substituting (2.16) and (2.17) into (2.15) yields
A combination of (3.5), (2.13) and (2.18) gives
for $p\ge 2$. Let $C_p$ be the number on the right-hand side of the above inequality. The error bound (2.10) holds for $p\ge 2$. Note from (2.6) that
Hence, we may improve the error bound (2.10) to
for $p\ge 1$, where
Using the explicit value
in (2.6), we can summarize our results as follows.
Theorem 2.1 Let $\beta >0$ and $N=2n+1 \geq 1$. Then for any $p\ge 1,$ we have
where
with
and
In particular, $\widetilde C_p\to |A_p(\coth \beta )|$ as $N\to +\infty$; namely, the error bound is close to the absolute value of the first neglected term when $N$ is large. The coefficients $A_j(\coth \beta )$ are polynomials in $\coth \beta$ of degree $3j$ with rational coefficients, and these polynomials can be calculated using a recursive formula given in appendix A.1.
3. Case II: $x=\cos \alpha \in [0,\,1)$ with $\alpha \in (0,\,\pi /2]$
The saddle points are $t_\pm =e^{\pm i\alpha }/2$, and we have
Now, we define an analytic function $w(z)$ via
where the branch of the square root function is chosen so that
for $z$ in a small complex neighbourhood of $t_+$. The steepest descent contour passing through $t_+$ is described via the equation ${\rm Im}\ [f(z)-f(t_+)]=0$. By using polar coordinates $z=re^{i\theta }$, this equation becomes
with
We are only interested in solutions in the upper half-plane; namely, $\theta \in (0,\,\pi )$. By the quadratic formula, there are exactly two solutions:
and
where $r_{d+}(\theta )$ and $r_{a+}(\theta )$ are defined as piecewise functions so that they are differentiable at $\theta =\alpha$. This is because the function $\sqrt {\cos ^{2}\alpha -(\theta -\theta _0)/\tan \theta }$ is not differentiable at $\theta =\alpha$. For any $\theta$ in a small real neighbourhood of $\alpha$, we have, by Taylor expansion, $\cos ^{2}\alpha -(\theta -\theta _0)/\tan \theta =(\theta -\alpha )^{2}+\mathcal {O}((\theta -\alpha )^{3})$ and $\sqrt {\cos ^{2}\alpha -(\theta -\theta _0)/\tan \theta }=|\theta -\alpha |+\mathcal {O}((\theta -\alpha )^{2})$. Hence,
In particular, we have
Now, we define two contours
Recall that $-Nf(t)$ is the phase function in the integral representation of $H_n(\sqrt N \cos \alpha )$. It is readily seen from the asymptotic behaviours $f(z)\sim z^{2}$ as $z\to \infty$ and $2f(z)\sim \ln z$ as $z\to 0$ that $\Gamma _{d+}$ is a steepest descent contour and $\Gamma _{a+}$ is a steepest ascent contour passing through the saddle point $t_+=e^{i\alpha }/2$. Moreover, the only zero of $f(z)-f(t_+)$ in the upper half-plane is $t_+$. Hence, $w(z)$ defined in (3.1) is analytic for all ${\rm Im}\ z>0$ and positive $z>0$. For $z=se^{i\pi }$ with $s>0$, we define
We can define in a similar manner the contours $\Gamma _{d-}$ and $\Gamma _{a-}$ which, together with $\Gamma _{d+}$ and $\Gamma _{a+}$, are illustrated in figure 3.
Now, we return to the integral representation of $H_n(\sqrt N \cos \alpha )$ in (1.1) and deform the contour of integration to the steepest descent contours $\Gamma _{d+}\cup \Gamma _{d-}$ oriented with increasing parameter $\theta$ (i.e., counter-clockwise direction):
where the last equality defines the two integrals $I_+$ and $I_-$, respectively. Since the contours $\Gamma _{d+}$ and $\Gamma _{d-}$ are symmetric with respect to the real axis, have opposite orientations, and $f(t) = \overline {f(\bar t)}$, the integrals $I_+$ and $I_-$ are complex conjugates of each other. Therefore, it suffices to study the integral $I_+$. First, we introduce the change of variable $u=w(t)$:
Since $w(z)$ is analytic in the upper half-plane and $|w(z)|/|z|\to 1$ as $z\to \infty$, we can infer from the Cauchy integral formula that
where $\mathbb {R}^{+}$ is the positive real line and $\Gamma _0^{+}$ is the negative real line with argument $\pi$:
Both $\mathbb {R}^{+}$ and $\Gamma _0^{+}$ are oriented from left to right. Now, we have the following double integral representation:
For any $p\ge 0$, we use the identities (2.5) and
to derive
where
The coefficients $A_j(i\cot \alpha )$ are polynomials in $i\cot \alpha$ and these polynomials are identical to those appearing in the expansion (2.19) (see appendix A.1 for details).
To find an upper bound for the remainder $\widetilde {\varepsilon }_p(N,\,\alpha )$, we need to estimate the integral
By expressing $w(z)$ in polar coordinates, one can easily show that
Note from (3.1) that
Assume that $p\geq 1$. Substituting the above two formulas into (2.9) yields
For $s\in (0,\,2)$, we have
and for $s>2$, we have
Accordingly,
Denote the right-hand side of the above inequality by $C_p$. It then follows from (2.5), (3.4) and (2.9) that
provided $p \geq 1$. We define
and observe that
for any $p\geq 0$. Using the explicit value
in (3.2), and the fact that $H_n(\sqrt N \cos \alpha )=I_+ + I_- = I_+ +\overline {I_+}=2{\rm Re}\ I_+$, we can summarize our results as follows.
Theorem 3.1 Let $\alpha \in (0,\,\pi /2],$ $N=2n+1 \geq 1$ and $\theta _0 =\alpha -\sin (2\alpha )/2$. Then for any $p\ge 0,$ we have
where
with
and
In particular, $\widetilde C_p\to |{\rm Re} (A_p(i\cot \alpha )e^{-i(\theta _0 N-\pi /2)/2})|$ as $N\to +\infty ;$ namely, the error bound is close to the absolute value of the first omitted term for large $N$. The coefficients $A_j(i\cot \alpha )$ are polynomials in $i\cot \alpha$ of degree $3j$ with rational coefficients, and these polynomials can be calculated using a recursive formula given in appendix A.1.
4. Case III: $x=1$
The saddle points coincide $t_s=t_\pm =1/2$. We have $f(t_s)=-3/4-(\ln 2)/2$, $f'(t_s)=f''(t_s)=0$ and $f'''(t_s)=8$. We introduce the analytic function
with the branch of the cube root being chosen so that
for $z$ in a small complex neighbourhood of $t_s$. We shall extend the definition of $w(z)$ by analytic continuation to $\mathbb {C}\setminus (-\infty,\,0]$. First, we investigate the equation ${\rm Im} [f(z)-f(t_s)]={\rm Im}\ f(z)=0$, which in polar coordinates $z=re^{i\theta }$ is
This equation is clearly satisfied by all $z>0$. There are two further solutions:
and
where we have defined $r_\pm (\theta )$ as piecewise functions so that they are differentiable at $\theta =0$. Note that since $1-\theta \cot \theta =\theta ^{2}/3+\mathcal {O}(\theta ^{3})$ for small $|\theta |>0$, we have
and
for small values of $\theta$. Especially,
Now, we introduce two contours
which we illustrate in figure 4.
We have the following lemma.
Lemma 4.1 Let $z_+(\theta )=r_+(\theta )e^{i\theta }$ for $\theta \in (-\pi /2,\,\pi );$ the function $f(z_+(\theta ))-f(t_s)$ increases from $-\infty$ to $+\infty$ as $\theta$ increases from $-\pi /2$ to $\pi$. Let $z_-(\theta )=r_-(\theta )e^{i\theta }$ for $\theta \in (-\pi,\,\pi /2);$ the function $f(z_-(\theta ))-f(t_s)$ decreases from $+\infty$ to $-\infty$ as $\theta$ increases from $-\pi$ to $\pi /2$.
Proof. For small $\theta$, we have from $f'(t_s)=f''(t_s)=0$, $f'''(t_s)=8$ and (4.4) that
and
Hence, the function $f(z_+(\theta ))-f(t_s)$ is increasing for small $\theta$. We claim that $f(z_+(\theta ))-f(t_s)$ is increasing for all $\theta \in (-\pi /2,\,\pi )$. Assume to the contrary that
for some $\theta \in (-\pi /2,\,0)\cup (0,\,\pi )$. Since $f'(z)\neq 0$ for any $z\neq t_s$, we have $r_+'(\theta )=r_+(\theta )=0$. This contradicts to the equation (4.3). Hence, the function $f(z_+(\theta ))-f(t_s)$ is increasing for all $\theta \in (-\pi /2,\,\pi )$. Note that $r_+(\theta )\to \infty$ as $\theta$ tends to $-\pi /2$ from the right or $\pi$ from the left. It is readily seen that $f(z_+(\theta ))-f(t_s)=r_+(\theta )^{2}\cos (2\theta )-2r_+(\theta )\cos \theta +(1/2)\ln r_+(\theta )-f(t_s)$ increases from $-\infty$ to $+\infty$ as $\theta$ increases from $-\pi /2$ to $\pi$.
The second statement of the lemma follows in a similar manner.
Since the solutions of ${\rm Im}\ f(z)=0$ are the union $\Gamma _1\cup \Gamma _2\cup \mathbb {R}^{+}$, the equation $f(z)-f(t_s)=0$ has a unique solution $z=t_s$ in $\mathbb {C}\setminus (-\infty,\,0]$. This implies that $w(z)$ defined in (4.1) can be analytically continued to $\mathbb {C}\setminus (-\infty,\,0]$. Moreover, it follows from (4.2) that $w(z_+(\theta ))$ increases from $-\infty$ to $\infty$ as $\theta$ increases from $-\pi /2$ to $\pi$.
Next, we let $\Gamma _{1+}$ be the portion of $\Gamma _1$ in the upper half-plane and $\Gamma _{2-}$ be the portion of $\Gamma _2$ in the lower half-plane; namely,
Both $\Gamma _{1+}$ and $\Gamma _{2-}$ are oriented with increasing $\theta$ (counter clockwise direction). We now deform the contour of integration for $H_n(\sqrt N)$ in (1.1) from $\Gamma$ to $\Gamma _{1+}\cup \Gamma _{2-}$:
where the last equality defines the two integrals $I_+$ and $I_-$, respectively. Since the contours $\Gamma _{1+}$ and $\Gamma _{2-}$ are symmetric with respect to the real axis, have opposite orientations, and $f(t) = \overline {f(\bar t)}$, the integrals $I_+$ and $I_-$ are complex conjugates of one another. Thus, it is sufficient to study the integral $I_+$. First, we introduce the change of variable $u = w(t)$:
We denote by $\Gamma _0^{\pm }$ the negative real lines with arguments $\pm \pi$ and orientation from left to right. For convenience, we denote by $\Gamma _*$ the union of $\Gamma _0^{+}$ oriented to the right and $\Gamma _0^{-}$ oriented to the left. Since $w(z)$ is analytic in $\mathbb {C}\setminus (-\infty,\,0]$ and $|w(z)|/|z|^{2/3}\to 1$ as $z\to \infty$, we can assert from the Cauchy integral formula that
Substituting this into (4.5) gives
For any $p\ge 1$, we use the identities
and
to find
where
It is shown in appendix A.2 that the coefficients $D_j$ are rational numbers which can be calculated via a recurrence relation.
To find an upper bound for the remainder $\widetilde {\varepsilon }_p(N)$, we need to estimate the integral
By expressing $w(z)$ in polar coordinates, one can easily show that
We also note from (4.1) that
for $z=se^{\pm i\pi }\in \Gamma _*$. Assume that $p\ge 4$. Substituting the above two formulas into (4.10) yields
For $z=se^{\pm i\pi }\in \Gamma _*$ with $s=|z|\in (0,\,1)$, we have
For $z=se^{\pm i\pi }\in \Gamma _*$ with $s=|z|>1$, we have
Consequently,
Denote the right-hand side of the above inequality by $C_p$. It then follows from (4.6), (4.9) and (4.10) that
provided $p \geq 4$. We define
and observe that
for any $p\ge 3$. Using the explicit value $-f(t_s) = 3/4+(\ln 2)/2$ in (4.7), and the fact that $H_n(\sqrt N)=I_+ + I_- = I_+ +\overline {I_+}=2{\rm Re}\ I_+$, we can summarize our results as follows.
Theorem 4.2 Denote $N=2n+1 \geq 1$. Then for any $p\ge 3,$ we have
where
with
and
In particular, if $p$ is not divisible by $3$, $\widetilde C_p \to \left |(3/4)^{p/3}D_p\sin (2\pi p/3)\right |\Gamma (p/3)$ as $N \to +\infty ;$ namely, the error bound is close to the absolute value of the first neglected term when $N$ is large. The coefficients $D_j$ are rational numbers and can be calculated using a recurrence relation given in appendix A.2.
5. Numerical examples
In tables 1, 2 and 3, we present some numerical results that demonstrate the accuracy of our error bounds given in theorems 2.1, 3.1 and 4.2, respectively. We can infer from the tables that the bounds are rather realistic, that is, they do not seriously overestimate the actual error.
Acknowledgments
We are very grateful to the anonymous referee for his/her careful reading and valuable suggestions which have helped to improve the presentation of this paper. The second author was supported by the JSPS Postdoctoral Research Fellowship No. P21020.
Appendix A. The computation of the coefficients
Appendix A.1 The coefficients $A_j$
Our starting point is the integral representation (2.7). We first reverse the orientation of the path $\Gamma _0 \cup \Gamma _+$ and thereby remove the minus sign in front of the integral. Next, by an appeal to Cauchy's theorem, the contour of integration can be shrunk into a small, positively oriented circle around $t_-$. After performing the change of variable $u=ze^{\beta }$, we arrive at
The contour of integration is a small loop surrounding $1/2$ in the positive sense. We would like to obtain a recursive scheme for the evaluation of the integrals in (A.1). To this end, we employ a method of Lauwerier [Reference Lauwerier6]. The method requires the power series expansions about $1/2$ of the functions appearing in the integrand of (A.1), namely
and
respectively. An application of Lauwerier's method [Reference Lauwerier6, eqs. (8) and (9)] (with the slightly different notation $P_j (t):=(-1)^{j} q_j (-t/2)$) then yields
where $P_0(t)=1$ and
for $j\geq 1$. It is readily seen from (A.3) that $P_j(t)$ is a polynomial in $t$ of degree $j$ with rational coefficients. With the substitution $t=\coth \beta$, (A.2) reduces to the more compact form
The $A_j(t)$ is the product of $(1+t)^{j}$ and a polynomial in $t$ of degree $2j$, and has rational coefficients. In particular, $A_0(t)=1$ and
Consider now the right-hand side of (3.3). We deform the contour of integration into a small, positively oriented circle around $t_+$, employ the identity $\sqrt {2\sin \alpha }=e^{-\pi i/4}e^{-i\alpha /2}\sqrt {e^{2i\alpha }-1}$, and perform the change of variable $u=ze^{-i\alpha }$, to obtain the equivalent form
Since $\cos \alpha =\cosh (-i\alpha )$, we see that (A.1) is identical to (A.4) when $\beta$ is replaced by $-i\alpha$. This confirms that the right-hand side of (3.3) is indeed $A_j(\coth (-i\alpha ))=A_j(i\cot \alpha )$.
Appendix A.2 The coefficients $D_j$
We proceed similarly as in the case of the coefficients $A_j$. By Cauchy's theorem, we can deform the contour of integration in (4.8) into a small, positively oriented loop contour around the saddle point $t_s=1/2$, to obtain
The square and cube roots in the second integral are defined to be positive for positive real $z$ and are defined by continuity elsewhere. To apply Lauwerier's method, we expand the functions in the rightmost integral in (A.5) into Taylor series about $1/2$:
Lauwerier's method [Reference Lauwerier6, eqs. (8) and (9)] (with the notation $Q_j(t) :=(-1)^{j-1} q_{j-1}(3t/4)$) then gives
where $Q_1(t)=1$ and
for $j\geq 1$. It is easy to verify that $Q_j(t)$ is a polynomial in $t$ of degree $j-1$ with rational coefficients. When integrating $Q_j (s)$ term-by-term in (A.6), the factor $\Gamma (j/3)$ in the denominator cancels, showing that $D_j$ is always a rational number. In particular,