, , March 7, 2005

In this paper, we are going to consider semi-infinite Jacobi matrices

\begin{matrix} J = (\begin{matrix} b_{1} & a_{1} & 0 & . . . \\ a_{1} & b_{2} & a_{2} & . . . \\ 0 & a_{2} & b_{3} & . . . \\ . . . & . . . & . . . & . . . \end{matrix}) \end{matrix}

(1.1)

whose Jacobi parameters have exponential decay (i.e.,

{limsup}_{n \to \infty} (| b_{n} | + | a_{n} - 1 |)^{1 / n} < 1

. As explained in the first two papers of this series [2, 3] (and well-known earlier), such a

J

has an associated Jost function,

u

, defined and analytic in a neighborhood of

\bar{D}

where

D = {z \in C | | z | < 1}

As is standard,

J

describes the recursion relations for orthogonal polynomials on the real line (OPRL). There is a probability measure,

γ

, so that the orthonormal polynomials,

p_{n} (x)

[14, 10, 11] , defined by

γ

obey

\begin{matrix} x p_{n} (x) = a_{n + 1} p_{n + 1} (x) + b_{n + 1} p_{n} (x) + a_{n} p_{n - 1} (x) \end{matrix}

(1.2)

γ

is the spectral measure for

J

and vector

(100 . . .)^{t}

and the

a

's and

b

's can be obtained from

γ

by Gram-Schmidt on the moments.

u

is defined by

γ

via the following three facts:

(i) $u (z) = 0$ for $z \in D$ if and only if $z + z^{- 1}$ is an eigenvalue of $J$ .

(ii) The support of

d γ_{s}

, the singular part of

d γ

, is a finite set of eigenvalues in

R \ [- 2, 2]

and

\begin{matrix} d γ ↾ [- 2, 2] = f (x) d x \end{matrix}

(1.3)

where for any

θ \in [0, 2 π)

\begin{matrix} f (2 cos θ) = \frac{1}{π} [\frac{sin θ}{| u (e^{i θ}) |^{2}}] \end{matrix}

(1.4)

(iii) $u (0) > 0$ .

These determine

u

by standard theory of nice analytic functions on

D

u

does not determine

γ

in many cases. For by 1.3 / 1.4 ,

u

determines the a.c. part of

γ

and the positions of the pure points but not their weights. We prefer to normalize the weights by looking at

\begin{matrix} M (z) = - \int \frac{d ρ (x)}{x - (z + z^{- 1})} \end{matrix}

(1.5)

and looking at the residues of

M

at the points where

u (z) = 0

In any event,

u

plus the weights are spectral data, and our goal here is to produce equivalences between this spectral data side and the recursion coefficient side.

To state our main theorems, we define

Definition. A sequence,

(x_{0}, . . ., x_{n}, . . .)

, of complex numbers is said to have an asymptotic series up to

R > 1

if and only if there exists

μ_{1}, . . ., μ_{K (R)}

{z | 1 < | z | < R}

and polynomials

p_{1}, . . ., p_{K (R)}

so that

\begin{matrix} {limsup}_{n \to \infty} {| x_{n} - \sum_{j = 1}^{K (R)} p_{j} (n) μ_{j}^{- n} |}^{1 / n} \leq R^{- 1} \end{matrix}

(1.6)

We say

(x_{0}, . . .)

has a complete asymptotic series if it has one for each

R > 1

It is easy to see that the

x

's uniquely determine the

p

's and

μ

's and that

Theorem 1.1.

{x_{n}}_{n = 0}^{\infty}

has an asymptotic series up to

R

if and only if

\begin{matrix} f (z) \equiv \sum_{n = 0}^{\infty} x_{n} z^{n} \end{matrix}

(1.7)

is meromorphic in

{z | | z | < R}

with no singularities in a neighborhood of

\bar{D}

and finitely many poles in the region.

{x_{n}}_{n = 1}^{\infty}

has a complete asymptotic series if and only if

f

is entire meromorphic.

Indeed, the poles are at the

μ_{j}

and their orders are one plus the degrees of the

p_{j}

We say a set of Jacobi parameters has an asymptotic series up to

R

if and only if the sequence

\begin{matrix} (1, - b_{1}, 1 - a_{1}^{2}, - b_{2}, 1 - a_{2}^{2}, . . .) \end{matrix}

(1.8)

has an asymptotic series up to

R

. Thus, the function

f

\begin{matrix} B (z) = 1 - \sum_{n = 0}^{\infty} [b_{n + 1} z^{2 n + 1} + (a_{n + 1}^{2} - 1) z^{2 n + 2}] \end{matrix}

(1.9)

B

will enter naturally below, but we note the following interpretation: If

J_{0}

is the Jacobi matrix with

a_{n} = 1

b_{n} \equiv 0

, and

δ J = J - J_{0}

, then (see Lemma 6.2 of [2] ):

\begin{matrix} Tr (δ J (J_{0} - (z + z^{- 1}))^{- 1}) = - (z^{- 1} - z)^{- 1} {\sum_{n = 1}^{\infty} b_{n} (1 - z^{2 n}) + 2 \sum_{n = 1}^{\infty} (a_{n} - 1) (z - z^{2 n + 1})} \end{matrix}

(1.10)

Moreover (see Theorem 2.16 of [9] ),

\begin{matrix} u (z) = {(^{\infty} \prod_{j = 1} a_{j})}^{- 1} det (1 + δ J (J_{0} - (z + z^{- 1}))^{- 1}) \end{matrix}

(1.11)

Taking into account that

a_{n}^{2} - 1 = 2 (a_{n} - 1) + O ((a_{n} - 1)^{2})

and

det (1 + A) = 1 + Tr (A) + O (∥ A ∥_{1}^{2})

, we see that if

δ J

is trace class, then

\begin{matrix} - (z^{- 1} - z) (^{\infty} \prod_{j = 1} a_{j}) u (z) = c + z B (z) + O (∥ δ J ∥_{1}^{2}) \end{matrix}

(1.12)

for a constant,

c

. Thus,

B (z)

is a kind of first-order (Born) approximation to

u

In some ways, our main result in this paper is

Theorem 1.2. The Jacobi parameters have a complete asymptotic series if and only if

u

is an entire meromorphic function. Equivalently,

B (z)

is entire meromorphic if and only if

u (z)

is.

Of course, one wants to understand the relation between the poles of

u

and those of

B

. Both for that understanding and because we will actually use them in our proofs in Section 3 , it pays to review our recent results [12] on the analogous problem for orthogonal polynomials on the unit circle (OPUC). The basics (see [10, 11] for background) associate to a nontrivial probability measure,

μ

, on

\partial D

a sequence of Verblunsky coefficients defined by

\begin{matrix} Φ_{n + 1} (z) = z Φ_{n} (z) - {\bar{α}}_{n} Φ_{n}^{*} (z) \end{matrix}

(1.13)

where

Φ_{n}

are the monic orthogonal polynomials for

μ

and

\begin{matrix} Φ_{n}^{*} (z) = z^{n} \bar{Φ_{n} (1 / \bar{z})} \end{matrix}

(1.14)

In place of

B

, [12] uses

\begin{matrix} S (z) = 1 - \sum_{j = 1}^{\infty} α_{j - 1} z^{j} \end{matrix}

(1.15)

and, in place of

u

, the Szegő function

\begin{matrix} D (z) = exp (\int \frac{e^{i θ} + z}{e^{i θ} - z} log (w (θ)) \frac{d θ}{4 π}) \end{matrix}

(1.16)

where

d μ = w (θ) \frac{d θ}{2 π} + d μ_{s}

. One also defines

\begin{matrix} r (z) = \frac{D^{- 1} (z)}{\bar{D^{- 1} (1 / \bar{z})}} \end{matrix}

(1.17)

The main theorems of [12] are:

Theorem 1.3 ([4] ). If

limsup | α_{n} |^{1 / n} = R^{- 1} < 1

, then

r (z) - S (z)

is analytic in

{z | 1 - δ < | z | < R^{3}}

for some

δ > 0

Remarks. 1. This result is due to Deift-Ostensson [4] , but [12] has a new proof.

Earlier, [10] proved the weaker result when

R^{3}

is replaced by

R^{2}

2. The point is that

r

and

S

both have singularities on

| z | = R

. This theorem says they cancel, as do other singularities in

{z | R < | z | < R^{3}}

3. [12] has explicit examples where

r (z) - S (z)

has singularities on

{z | | z | = R^{3}}

and shows that this is the case generically. So

R^{3}

is best possible.

Given a discrete set,

Ω \subset {z | | z | > 1}

, with limit points only at

\infty

, we define

\begin{matrix} G^{2 j - 1} (Ω) & = {μ_{1} . . . μ_{j} {\bar{μ}}_{j + 1} . . . {\bar{μ}}_{2 j - 1} | μ_{k} \in Ω} \end{matrix}

(1.18)

\begin{matrix} G (Ω) & =^{\infty} ⋃_{j = 1} G^{2 j - 1} (Ω) \end{matrix}

(1.19)

\begin{matrix}  \end{matrix}

Theorem 1.4 ([12] ).

S

is entire analytic if and only if

D^{- 1}

is. If

T

is the set of poles of

S (z)

and

P

the poles of

D^{- 1} (z)

, then

\begin{matrix} T \subset G (P) P \subset G (T) \end{matrix}

(1.20)

Analogously to Theorem 1.3 , we will prove in Section 2 that

Theorem 1.5. Suppose

\begin{matrix} {limsup}_{n \to \infty} (| a_{n}^{2} - 1 | + | b_{n} |)^{1 / 2 n} = R^{- 1} < 1 \end{matrix}

(1.21)

Then

\begin{matrix} (1 - z^{2}) u (z) + z^{2} \bar{u (1 / \bar{z})} B (z) \end{matrix}

(1.22)

is analytic in

{z | R^{- 1} < | z | < R^{2}}

Remarks. 1. [3] has necessary and sufficient conditions on

{u, weights}

for 1.21 to hold. If there are no eigenvalues of

J

outside

[- 2, 2]

, the condition is that

u

is analytic in

{z | | z | < R}

u

is real on

R

\bar{u (\bar{z})} = u (z)

and thus, 1.22 could be written

(1 - z^{2}) u (z) + u (1 / z) B (z)

; we write it as we do for analogy with the OPUC case.

3. The point, of course, is that

B

has singularities on

{z | | z | = R}

, so this theorem implies a cancellation either via zeros of

\bar{u (1 / \bar{z})}

or singularities of

u

. Since

\bar{u (1 / \bar{z})}

can have zeros in

| z | > 1

(while

\bar{D (\bar{z})^{- 1}}

cannot), the situation is somewhat different from OPUC. We will discuss this further in Section 2 .

4. As we will show in Section 3 , the function in 1.22 often has a singularity at

z = R^{2}

, so one cannot increase the

R^{2}

R^{3}

as one can in the OPUC case. The reason for this difference will become clear in Section 3 .

For the analog of Theorem 1.4 , we need to define a larger set than

G

. In our situation,

u

and

B

are real on

R

so their poles are symmetric about

R

. So for this, we will suppose

Ω \subset {z | | z | > 1}

with limit point only at infinity, and

\begin{matrix} \bar{Ω} = Ω \end{matrix}

(1.23)

In that case, for any

m

, we define

\begin{matrix} G^{(m)} (Ω) = {μ_{1} . . . μ_{m} | μ_{k} \in Ω} \end{matrix}

(1.24)

When 1.23 holds, this agrees with the previous definition if

m = 2 k - 1

\begin{matrix} \tilde{G} (Ω) = [^{\infty} ⋃_{m = 1} G^{(m)} (Ω)] \cup [-^{\infty} ⋃_{m = 1} G^{(m)} (Ω)] \end{matrix}

(1.25)

Our main results refine Theorem 1.2 :

Theorem 1.6. Let

J

have no spectrum outside

[- 2, 2]

and let

u

be entire meromorphic and nonvanishing at

z = \pm 1

. Let

P

be the poles of

u

and

T

the poles of

B

. Then

\begin{matrix} P \subset \tilde{G} (T) T \subset \tilde{G} (P) \end{matrix}

(1.26)

To state the result when there are bound states, we recall and extend a notion from [3] .

Definition. Let

u

be a meromorphic function and

z_{0} \in D

a point with

u (z_{0}) = 0

(so

z_{0}

is real and

z_{0} + z_{0}^{- 1} \in σ (J)

z_{0}

is called a noncanonical zero for

J

if and only if

1 / z_{0}

is not a pole of

u

and

\begin{matrix} {lim}_{z \to z_{0}} (z - z_{0}) M (z) \neq - (z_{0} - z_{0}^{- 1}) {[u^{'} (z_{0}) u (\frac{1}{z_{0}})]}^{- 1} \end{matrix}

(1.27)

Thus,

z_{0}

is not noncanonical (which we will call canonical) if

u

is regular at

1 / z_{0}

and equality holds in 1.27 . Here is what we will prove in case there are bound states or

u (\pm 1) = 0

Theorem 1.7. Suppose

u

is entire meromorphic. Let

T

be the poles of

B

. Let

P_{1}

be the poles of

u

and

P_{2}

the

{z^{- 1} | z

is a noncanonical zero for

J}

. Let

P = P_{1} \cup P_{2}

. Then 1.26 holds.

As in [12] , one can easily prove results relating meromorphicity of

u

{z | | z | < R^{2 ℓ - 1}}

to meromorphicity of

B

there.

In Section 2 , we use the Geronimo-Case equations to prove Theorem 1.5 . In Section 3 , we use the second Szegő map from OPRL to OPUC to prove Theorem 1.6 .

In Section 4 , we extend the analysis of [3] to obtain Theorem 1.7 from Theorem 1.6 .

This research was completed during my stay as a Lady Davis Visiting Professor at Hebrew University, Jerusalem. I'd like to thank H. Farkas and Y. Last for the hospitality of the Mathematics Institute at Hebrew University.

2 The Geronimo-Case Equations and the $R^{- 2}$ Result

In this section, we will prove Theorem 1.5 using a strategy similar to that used in [12] to prove Theorem 1.3 . There the critical element was the use of Szegő recursion (1.13) and its

^{*}

, that is,

\begin{matrix} Φ_{n + 1}^{*} (z) = Φ_{n}^{*} (z) - α_{n} z Φ_{n} (z) \end{matrix}

(2.1)

z

and

1 / \bar{z}

Here we will instead use the Geronimo-Case equations [5] in the form introduced in [3] . Define

\begin{matrix} C_{n} (z) = z^{n} P_{n} (z + \frac{1}{z}) \end{matrix}

(2.2)

The equations

\begin{matrix} C_{n} (z) & = (z^{2} - b_{n} z) C_{n - 1} (z) + G_{n - 1} (z) \end{matrix}

(2.3)

\begin{matrix} G_{n} (z) & = G_{n - 1} (z) + [(1 - a_{n}^{2}) z^{2} - b_{n} z] C_{n - 1} (z) \end{matrix}

(2.4)

\begin{matrix}  \end{matrix}

are the unnormalized GC equations. With initial condition

G_{0} (z) = C_{0} (z) = 1

, they define monic polynomials of degree at most

n

C_{n}

has the form 2.2 , and if

\begin{matrix} \sum_{n = 1}^{\infty} | a_{n}^{2} - 1 | + | b_{n} | < \infty \end{matrix}

(2.5)

then for

| z | < 1

\begin{matrix} {lim}_{n \to \infty} G_{n} (z) = (^{\infty} \prod_{j = 1} a_{j}) u (z) \end{matrix}

(2.6)

(see Theorem A.3 of [3] ).

2.3 / 2.4 have a structure somewhat like 1.13 / 2.1 . The difference is that 1.14 is replaced by

\begin{matrix} C_{n} (z) = z^{2 n} C_{n} (\frac{1}{z}) \end{matrix}

(2.7)

as is obvious from 2.2 . We introduce

f = \tilde{O} (g)

if and only if for all

ɛ > 0

| f | / | g |^{1 - ɛ} \to 0

Lemma 2.1. If 1.21 holds, then for

z \in D

\begin{matrix} (i) & | G_{n} (z) - u (z) | & \leq \tilde{O} (R^{- 2 n}) \end{matrix}

(2.8)

\begin{matrix} (ii) & | C_{n} (z) - \frac{u (z)}{1 - z^{2}} | & \leq \tilde{O} ([max (| z |, R^{- 1})]^{2 n}) \end{matrix}

(2.9)

\begin{matrix}  \end{matrix}

Proof. (i) By Theorem A.3 of [3] ,

\begin{matrix} {lim}_{n \to \infty} C_{n} (z) = \frac{u (z)}{1 - z^{2}} \end{matrix}

(2.10)

By 2.4 and

{sup}_{n} | C_{n} (z) | < \infty

, we see

\begin{matrix} | G_{n} (z) - u (z) | & \leq \sum_{m = n}^{\infty} | G_{m + 1} (z) - G_{m} (z) | \end{matrix}

\begin{matrix} \leq ({sup}_{n} | C_{n} (z) |) \sum_{m = 1}^{\infty} (| 1 - a_{n + m}^{2} | + | b_{n + m} |) \end{matrix}

\begin{matrix} = \tilde{O} (R^{- 2 n}) \end{matrix}

\begin{matrix}  \end{matrix}

since the series of bounds converges exponentially.

(ii) By 2.3 ,

| C_{n} - G_{n - 1} - z^{2} C_{n - 1} | \leq {sup}_{n} | C_{n} (z) | | b_{n} |

so iterating,

\begin{matrix} | C_{n} - \sum_{j = 0}^{n - 1} G_{n - j - 1} z^{2 j} | & \leq | z |^{2 n} + {sup}_{n} | C_{n} (z) | \sum_{j = 0}^{n - 1} | b_{n - j} | | z^{2 j} | \end{matrix}

\begin{matrix} \leq \tilde{O} (max (| z |, R^{- 1})^{2 n}) \end{matrix}

\begin{matrix}  \end{matrix}

By 2.8 ,

| \sum_{j = 0}^{n - 1} (G_{n - j - 1} - u) z^{2 j} | \leq \tilde{O} (max (| z |, R^{- 1})^{2 n})

Since

\sum_{j} z^{2 j} u = (1 - z^{2})^{- 1} u

, we have 2.9 . □

Proof of Theorem 1.5 . By 2.4 and 2.7 for

| z | > 1

\begin{matrix} | G_{n + 1} - G_{n} | \leq [{sup}_{n} | C_{n} (\frac{1}{z}) |] | z |^{2 n + 2} [| 1 - a_{n}^{2} | + | b_{n} |] \end{matrix}

(2.11)

which proves that for

1 < | z | < R

G_{n}

converges uniformly, so by the maximum principle, we have convergence for

| z | < R

, so

u

has an analytic continuation to that region. In that region,

\begin{matrix} u (z) & = 1 + \sum_{n = 0}^{\infty} (G_{n + 1} (z) - G_{n} (z)) \end{matrix}

\begin{matrix} = 1 + \sum_{n = 0}^{\infty} ((1 - a_{n + 1}^{2}) z^{2} - b_{n + 1} z) C_{n} (z) \end{matrix}

(2.12)

\begin{matrix} = (\frac{u (\frac{1}{z})}{1 - \frac{1}{z^{2}}}) (B (z) - 1) + 1 + \sum_{n = 0}^{\infty} f_{n} (z) \end{matrix}

(2.13)

\begin{matrix}  \end{matrix}

where

\begin{matrix} f_{n} (z) = ((1 - a_{n + 1}^{2}) z^{2} - b_{n + 1} z) z^{2 n} (C_{n} (\frac{1}{z}) - \frac{u (\frac{1}{z})}{1 - \frac{1}{z^{2}}}) \end{matrix}

(2.14)

Thus

\begin{matrix} (1 - z^{2}) u (z) + u (\frac{1}{z}) z^{2} B (z) = u (\frac{1}{z}) z^{2} + (1 - z^{2}) + \sum_{n = 0}^{\infty} (1 - z^{2}) f_{n} (z) \end{matrix}

(2.15)

Each function

f_{n}

is analytic in

{z | | z | > 1}

, so if we can prove that the sum converges uniformly in

{z | 1 < | z | < R^{2}}

, we know the left side of 2.15 has an analytic continuation in that region.

By 2.9 , for

| z | > 1

| C_{n} (\frac{1}{z}) - \frac{u (\frac{1}{z})}{1 - \frac{1}{z^{2}}} | \leq \tilde{O} (max {(\frac{1}{| z |}, R^{- 1})}^{2 n})

| z^{2 n} [C_{n} (\frac{1}{z}) - \frac{u (\frac{1}{z})}{1 - \frac{1}{z^{2}}}] | \leq \tilde{O} (max (1, | z | R^{- 1})^{2 n})

and thus,

| f_{n} (z) | \leq \tilde{O} (R^{- 2 n}) \tilde{O} (max (1, | z | R^{- 1})^{2 n})

For

1 < | z | < R

, this is

\tilde{O} (R^{- 2 n})

and so summable. For

R \leq | z | < R^{- 2}

, it is

\tilde{O} ((| z | R^{- 2})^{2 n})

and so also summable. □

u (\pm R^{- 1}) \neq 0

, 1.22 tells us that since

B

has a singularity on the circle of radius

R

, so must

u

. However, if

u (R^{- 1}) = 0

and/or

u (- R^{- 1}) = 0

, that zero can compensate for a pole in

B

and

u

can have a larger region of analyticity than

B

This is exactly what happens in the case of noncanonical weights, as explained in [3] .

3 The Second Szego Map and Jost Functions With No Bound States

In [13, 14] , Szegő defined two maps from the probability measures on

\partial D

invariant under

z \to \bar{z}

to the probability measures on

[- 2, 2]

; let us call them

{Sz}_{1}

and

{Sz}_{2}

Both are injective, but only

{Sz}_{1}

is surjective — and for this reason,

{Sz}_{1}

is the one most often used and studied (see [11,Section 13.1] ). Here we will see that

{Sz}_{2}

is also exceedingly useful, especially for studying Jost functions analytic in a neighborhood of

\bar{D}

and nonvanishing on

\bar{D}

(i.e.,

J

has no bound states and no resonance at

\pm 2

For a.c. measures, the relations are

\begin{matrix} d μ = w (θ) \frac{d θ}{2 π} {Sz}_{1} (d μ) = f_{1} (x) d x {Sz}_{2} (d μ) = f_{2} (x) d x \end{matrix}

(3.1)

where

w (θ) = w (- θ)

and (formulae (13.1.6) and (13.2.22) of [11] )

\begin{matrix} f_{1} (x) & = π^{- 1} (4 - x^{2})^{- 1 / 2} w (arccos (\frac{x}{2})) \end{matrix}

(3.2)

\begin{matrix} f_{2} (x) & = π^{- 1} c^{2} (4 - x^{2})^{1 / 2} w (arccos (\frac{x}{2})) \end{matrix}

(3.3)

\begin{matrix}  \end{matrix}

where

\begin{matrix} c = [2 (1 - | α_{0} |^{2}) (1 - α_{1})]^{- 1 / 2} \end{matrix}

(3.4)

Taking into account that

{Sz}_{1}

is a bijection of even measures on

\partial D

and all measures on

[- 2, 2]

, we see that

\begin{matrix} d γ \in ran ({Sz}_{2}) \Leftrightarrow \int_{- 2}^{2} (4 - x^{2})^{- 1} d γ (x) < \infty \end{matrix}

(3.5)

Proposition 3.1. If

d γ

has a Jost function

u

analytic in a neighborhood of

\bar{D}

and nonvanishing on

\bar{D}

, then

d γ \in ran ({Sz}_{2})

Proof. Since

d (2 cos θ) = - 2 sin θ d θ

and

(4 - 2 {cos}^{2} θ) = 4 {sin}^{2} θ

, by 1.4 the right side of 3.5 is equivalent to

\begin{matrix} \frac{1}{π} \int_{0}^{π} \frac{d θ}{| u (e^{i θ}) |^{2}} & = \int_{0}^{π} ({sin}^{2} θ)^{- 1} f (2 cos θ) sin θ d θ \end{matrix}

\begin{matrix} = 2 \int_{- 2}^{2} (4 - x^{2})^{- 1} f (x) d x < \infty \end{matrix}

\begin{matrix}  \end{matrix}

which is true if

| u |

is bounded away from zero. □

For our purposes, what is critical is:

Theorem 3.2. Let

d μ

be a nontrivial probability measure on

\partial D

obeying the Szegő condition with Verblunsky coefficients

{α_{n}}_{n = 0}^{\infty}

and Szegő function

D (z_{0})

. Let

d γ = {Sz}_{2} (d μ)

have Jost function

u

and Jacobi parameters

{a_{n}, b_{n}}_{n = 1}^{\infty}

. Then

\begin{matrix} b_{n + 1} & = α_{2 n} - α_{2 n + 2} - α_{2 n + 1} (α_{2 n} + α_{2 n + 2}) \end{matrix}

(3.6)

\begin{matrix} a_{n + 1}^{2} - 1 & = α_{2 n + 1} - α_{2 n + 3} - α_{2 n + 2}^{2} (1 - α_{2 n + 3}) (1 + α_{2 n + 1}) - α_{2 n + 3} α_{2 n + 1} \end{matrix}

(3.7)

\begin{matrix} u (z) & = (1 - | α_{0} |^{2}) (1 - α_{1}) D (z)^{- 1} \end{matrix}

(3.8)

\begin{matrix}  \end{matrix}

Remark. Formulae of the form 3.6 / 3.7 for

{Sz}_{1}

go back to Geronimus [6, 7] . For

{Sz}_{2}

, the earliest reference I am aware of is Berriochoa, Cachafeiro, and García-Amor [1] ; see also [8] .

Proof. 3.6 / 3.7 are (13.2.20)/(13.2.21) of [11] . To see 3.8 , note that, by 1.4 and 3.3 ,

\begin{matrix} | u (e^{i θ}) |^{- 2} & = π f_{2} (2 cos θ) (sin θ)^{- 1} \end{matrix}

\begin{matrix} = 2 c^{2} w \end{matrix}

\begin{matrix} = 2 c^{2} | D |^{2} \end{matrix}

\begin{matrix}  \end{matrix}

Thus, the absolute value of 3.8 holds if

| z | = e^{i θ}

. Since both sides are analytic, nonvanishing on

D

, and positive at

z = 0

, 3.8 holds for all

z

. □

3.6 / 3.7 first of all provide a second proof of Theorem 1.5 in case

u

is nonvanishing on

\bar{D}

and, more importantly, show generically that

R^{2}

is optimal. We note first:

Proposition 3.3. We have that

\begin{matrix} B (z) = α_{0} z^{- 1} + α_{1} + 1 + (S (z) - 1) (1 - z^{- 2}) + Q (z) \end{matrix}

(3.9)

where, if

\begin{matrix} limsup | α_{n} |^{1 / n} = R^{- 1} \end{matrix}

(3.10)

then

Q

is analytic in

{z | | z | < R^{2}}

Proof. By 1.9 , 1.15 , and 3.6 / 3.7 , we have 3.9 where

\begin{matrix} \begin{matrix} Q (z) = \sum_{n = 0}^{\infty} & α_{2 n + 1} (α_{2 n} + α_{2 n + 2}) z^{2 n + 1} & + {α_{2 n + 2}^{2} (1 - α_{2 n + 3}) (1 + α_{2 n + 1}) + α_{2 n + 3} α_{2 n + 1}} z^{2 n + 2} \end{matrix} \end{matrix}

(3.11)

By 3.10 ,

Q (z)

is analytic in

{z | | z | < R^{2}}

. □

Second proof of Theorem 1.5 when

u

is nonvanishing on

\bar{D}

. As we will show below (see Lemma 3.5 ), 1.26 implies 3.10 . By Theorem 1.3 and 3.8 , we conclude that

\begin{matrix} (z^{2} - 1) [u (z) - \bar{u (1 / \bar{z})} S (z)] \end{matrix}

(3.12)

is analytic in

{z | R^{- 1} < | z | < R^{3}}

. By Proposition 3.3 ,

\begin{matrix} z^{2} B (z) - (z^{2} - 1) S (z) \end{matrix}

(3.13)

is analytic in

{z | | z | < R^{2}}

, so by 3.12 , the function in 1.22 is analytic in

{z | R^{- 1} < | z | < R^{2}}

. □

Example 3.4. Suppose

α_{2 n} \equiv 0

(true if and only if

b_{n} \equiv 0

) and

α_{2 n + 1} = R^{- (2 n + 1)}

. Then, by 3.10 ,

\begin{matrix} Q (z) & = \sum_{n = 0}^{\infty} z^{2 n + 2} R^{- 4 n - 4} \end{matrix}

\begin{matrix} = z^{2} R^{- 4} (1 - z^{2} R^{- 4})^{- 1} \end{matrix}

\begin{matrix}  \end{matrix}

has poles at

z = \pm R^{2}

. This shows that 1.22 may not be analytic in any larger annulus than

{z | R^{- 1} < | z | < R^{2}}

. It is also clear that by a similar analysis, if

B

is meromorphic in

{z | | z | < R^{1 + ɛ}}

, then generically 1.22 will have singularities on the circle of radius

R^{2}

. The change from

R^{3}

R^{2}

in going from Theorem 1.3 to Theorem 1.5 is due to the quadratic terms in 3.6 / 3.7 . □

Above we used and below we will need:

Lemma 3.5. Suppose

{α_{n}}_{n = 0}^{\infty}

and

{a_{n}, b_{n}}_{n = 1}^{\infty}

are related by 3.6 / 3.7 , and

\begin{matrix} limsup | α_{n} |^{1 / n} \equiv R_{1}^{- 1} < 1 {limsup}_{n \to \infty} (| a_{n}^{2} - 1 | + | b_{n} |)^{1 / 2 n} \equiv R_{2}^{- 1} < 1 \end{matrix}

(3.14)

Then

R_{1} = R_{2}

. Moreover,

α_{n}

has a complete asymptotic series if and only if

{a_{n}, b_{n}}

do, and if

T

is the set of powers that enter for

{a_{n}, b_{n}}_{n = 1}^{\infty}

(i.e.,

T

are the poles of

B

) and

\tilde{T}

for

{α_{n}}_{n = 0}^{\infty}

(i.e.,

\tilde{T}

are the poles of

S

), then

\begin{matrix} T \subset \tilde{G} (\tilde{T}) \tilde{T} \subset \tilde{G} (T) \end{matrix}

(3.15)

Remark. For

{Sz}_{1}

, there are equations similar to 3.6 / 3.7 which have solutions where

{a_{n}, b_{n}}_{n = 0}^{\infty}

has rapid decay while

α_{2 n + 1} \sim n^{- 1}

at infinity. (Indeed, for

{Sz}_{1}

but not

{Sz}_{2}

, this happens for

J_{0}

; see Example 13.1.3 revisited in [11] .) In fact, the results in this paper plus [12] imply that

R_{1}^{- 1} < 1

if and only if

R_{2}^{- 1} < 1

Proof. It follows from 3.6 that if

R_{1}^{- 1}, R_{2}^{- 1} < 1

, then

R_{2} = R_{1}

since the nonleading terms are exponentially small. In addition, if

α_{n}

has a complete asymptotic series, one gets that

b_{n + 1}

and

a_{n + 1}^{2} - 1

individually have asymptotic series in

μ_{k}^{- 2 n}

with

μ_{k} \in \cup_{j = 1}^{\infty} G^{(j)} (\tilde{T})

. Since

\begin{matrix} c_{1} μ_{k}^{- 2 n} & = \frac{1}{2} (c_{1} + c_{2}) μ_{k}^{- 2 n} + \frac{1}{2} (c_{1} - c_{2}) (- μ_{k})^{- 2 n} \end{matrix}

(3.16)

\begin{matrix} c_{2} μ_{k}^{- 2 n - 1} & = \frac{1}{2} (c_{1} + c_{2}) μ_{k}^{- 2 n - 1} + \frac{1}{2} (c_{1} - c_{2}) (- μ_{k})^{- 2 n - 1} \end{matrix}

(3.17)

\begin{matrix}  \end{matrix}

we can combine into a single series by taking

- μ

's as well as

μ

's.

For the converse, note that since the

α

's decay exponentially,

α_{2 n} = \sum_{m = 0}^{\infty} b_{n + m + 1} + O (R^{- 2 n})

and similarly for

α_{2 n + 1}

and

\sum_{m = 0}^{\infty} (a_{n + m + 1}^{2} - 1)

. Plugging this into 3.6 and summing yields

α_{2 n}

and

α_{2 n + 1}

as explicit sums of products of four or fewer

b

's and

(1 - a^{2})

's plus an error of

O (R^{- 3 n})

. Iterating gives explicit formulae for

α

's as “polynomials” in

b

and

1 - a^{2}

of degree

k

plus an error of order

O (R^{- (k + 2) n})

This shows that if

a

and

b

have asymptotic series to order

R^{- (k + 1) n}

, so do

α_{2 n}

and

α_{2 n - 1}

with rates in

\cup_{j = 1}^{\infty} G^{(j)} (T)

. Using formulae like 3.16 / 3.17 , we can combine to a single series by using

- μ

's, so

\tilde{T} \subset \tilde{G} (T)

. □

Proof of Theorem 1.6 and Theorem 1.2 when $u$ is nonvanishing on $\bar{D}$ . Since $u$ is nonvanishing on $\bar{D}$ , $γ \in ran ({Sz}_{2})$ , so we can define $S, α_{n}$ , etc. If $u$ is entire meromorphic, by 3.8 , so is $D^{- 1}$ . Thus, by Theorem 1.4 , $S$ is entire meromorphic, and if $\tilde{T}$ are the poles of $S$ , then $\tilde{T} \subset G (P)$ By Lemma 3.5 , $B$ is meromorphic and $T \subset \tilde{G} (\tilde{T}) \subset \tilde{G} (G (P)) = \tilde{G} (P)$ Conversely, if $B$ is entire meromorphic, by Lemma 3.5 , so is $S$ , and if $\tilde{T}$ are the poles of $S$ , then $\tilde{T} \subset \tilde{G} (T)$ By Theorem 1.4 , $D^{- 1}$ , and so $u$ , is entire meromorphic and $P \subset G (\tilde{T}) \subset G (\tilde{G} (T)) = \tilde{G} (T) Λ$ □

4 Coefficient Stripping and Jost Functions With Bound States

As in [3] , we will go from the no bound state theorem to the general case (i.e., in our situation, Theorem 1.6 to Theorem 1.7 ) by coefficient stripping, that is, pass from

J

to the Jacobi matrix

J^{(m)}

with Jacobi parameters

{a_{n + m}, b_{n + m}}_{n = 1}^{\infty}

. By Theorem 3.1 of [2] , if

J

has a Jost function analytic in a neighborhood of

\bar{D}

, there exists a

k

with

σ (J^{(k)}) = [- 2, 2]

, and by a slight extension of the argument, we can also suppose its Jacobi function obeys

u^{(k)} (\pm 1) \neq 0

(for if

σ (J^{(k - 1)}) = [- 2, 2]

, if

u^{(k + 1)}

vanishes at

\pm 1

M^{(k - 1)} (z)

has a pole there and

u^{(k)} = u^{(k - 1)} M^{(k - 1)}

is nonvanishing). Thus, we claim that we need only prove (as we shall do below) that

Theorem 4.1. If

P = P_{1} \cup P_{2}

as in Theorem 1.7 and we make the

J

-dependence explicit, then

\begin{matrix} P (J) = P (J^{(1)}) \end{matrix}

(4.1)

Proof of Theorems 1.1 and 1.7 given Theorems 4.1 and 1.6 . Theorem 1.7 implies Theorem 1.1 . By 4.1 and induction, $P (J) = P (J^{(k)})$ where $k$ is chosen as above so Theorem 1.6 is applicable. 1.26 for $J^{(k)}$ implies it for $J$ . □

As in [3] , we will make use of the

M

-function, its connection to

u

, and the update relations. We define (consistently with 1.5 )

\begin{matrix} M^{(k)} (z) = 〈 δ_{1}, (z + z^{- 1} - J^{(k)})^{- 1} δ_{1} 〉 \end{matrix}

(4.2)

for

z \in D \ {w | w + w^{- 1} \in σ (J^{(k)})}

M^{(k)}

has poles at the set in

D

with

w + w^{- 1} \in σ (J^{(k)})

, and

u^{(k)}

has zeros there. The update equations ((2.4)/(2.5) of [3] ) are (initially for

z \in D

)

\begin{matrix} u^{(k + 1)} (z) & = a_{k + 1} z^{- 1} u^{(k)} (z) M^{(k)} (z) \end{matrix}

(4.3)

\begin{matrix} M^{(k)} (z)^{- 1} & = z + z^{- 1} - b_{k + 1} - a_{k + 1}^{2} M^{(k + 1)} (z) \end{matrix}

(4.4)

\begin{matrix}  \end{matrix}

Moreover, we have the analytic continuation of 1.4 plus

π f (2 cos θ) = Im M (e^{i θ})

for

θ \in [0, π]

\begin{matrix} [M (z) - \bar{M} (1 / \bar{z})] \bar{u (1 / \bar{z})} u (z) = z - z^{- 1} \end{matrix}

(4.5)

4.5 can be used to meromorphically continue

M

from

D

C

u

is entire meromorphic. Once one makes these continuations, 4.3 / 4.4 extend to all

z \in C

(as meromorphic relations including possible cancellations of poles and zeros). 4.3 / 4.4 also show that if

u

is entire meromorphic, so is

u^{(1)}

We begin by rephrasing the set

P_{2}

Proposition 4.2.

z_{0} \in C \ \bar{D}

is in

P_{2}

if and only if

(i) $z_{0}$ is not a pole of $u$ .
(ii) Both $z_{0}$ and $z_{0}^{- 1}$ are poles of $M$ .

Remark. In

D

, all poles of

M

are real, so (ii) implies that

z_{0}

is real.

Proof. By definition, $z_{0} \in P_{2}$ if and only if $u (z_{0}^{- 1}) = 0$ , $z_{0}$ is not a pole of $u$ , and 1.27 holds. Since $z_{0}^{- 1} \in D$ , $u (z_{0}^{- 1}) = 0 \Leftrightarrow z_{0} + z_{0}^{- 1} \in σ (J) \Leftrightarrow z_{0}^{- 1}$ is a pole of $M (z)$ . As shown in [3] , by 4.5 , if $u (z_{0}) = 0$ , $z_{0}$ has a pole of $M$ of order two or more and, of course, 1.27 holds at $z_{0}^{- 1}$ since the left side is infinite. If $u (z_{0}) \neq 0$ , 1.27 is precisely the condition, via 4.5 , that $M (z)$ has a pole at $z_{0}$ . □

We have been careful in considering situations where

z_{0}

is a pole of

u

and

z_{0}^{- 1}

is a zero of

u

. We need to consider that case separately:

Proposition 4.3. If

z_{0} \in C \ D

is a pole of

u

and

z_{0}^{- 1}

is a zero of

u

, then

z_{0}

is a pole of

u^{(1)}

Proof. Consider 4.5 near $z = z_{0}$ . Zeros of $u$ in $D$ are simple, so $u (z) \bar{u (1 / \bar{z})}$ either has a pole at $z_{0}$ or a finite nonzero limit. Thus, 4.5 shows $M (z) - \bar{M (1 / \bar{z})}$ must be regular (perhaps even zero) at $z_{0}$ . Since $\bar{M (1 / \bar{z})}$ has a pole at $z_{0}$ , $M (z)$ must have a pole there also. It follows that $u^{(1)} = u M$ has a pole (indeed, at least a second-order pole) at $z_{0}$ . □

Proposition 4.4. If

z_{0} \in P_{2} (J)

, then

z_{0} \in P_{1} (J^{(1)})

Remarks. 1. For

z_{0}

within the disk of analyticity of

u

, this result is in [3] . The proof here is essentially identical.

z_{0} \in P_{2} (J)

is essentially a statement of the vanishing of a “resonance eigenfunction,” so this says that such eigenfunctions cannot have successive zeros because of a second-order difference equation.

Proof. Suppose first that $u (z_{0}) \neq 0$ . By Proposition 4.2 , $M$ has a pole at $z_{0}$ , so $u^{(1)} = u M$ has a pole at $z_{0}$ .
If $u$ has a $k$ th-order zero, $k \geq 1$ , at $z_{0}$ , $\bar{u (1 / \bar{z})} u (z)$ has a $(k + 1)$ st-order zero, so $M (z) - \bar{M (1 / \bar{z})}$ has a $(k + 1)$ st-order pole $z_{0}$ by 4.5 . Since $M$ has simple poles at points in $D$ like $1 / z_{0}$ , $M$ has to have a $(k + 1)$ st-order pole at $z_{0}$ . Thus, $u^{(1)} = u M$ has a pole at $z_{0}$ . □

Proposition 4.5. If

z_{0} \in P_{1} (J)

and

z_{0} / \in P_{1} (J^{(1)})

, then

z_{0} \in P_{2} (J^{(1)})

Proof. By 4.4 , poles of

M^{(1)} (z)

are precisely at zeros of

M (z)

. Thus, by Proposition 4.2 , we need to prove that

\begin{matrix} z_{0} \in P_{1} (J) + z_{0} / \in P_{1} (J^{(1)}) \Rightarrow M (z_{0}) = M (\frac{1}{z_{0}}) = 0 \end{matrix}

(4.6)

Since

u

has a pole at

z_{0}

and

u^{(1)} = u M

does not,

M (z_{0}) = 0

. By Proposition 4.3 ,

z_{0} / \in P_{1} (J^{(1)})

implies

u (1 / \bar{z}) \neq 0

. Thus,

u (z) \bar{u (1 / \bar{z})}

has a pole at

z_{0}

. 4.5 then implies that

M (z) - \bar{M (1 / \bar{z})} |_{z = z_{0}} = 0

. Since

M (z_{0}) = 0

, we conclude

M (1 / {\bar{z}}_{0}) = 0

This proves 4.6 . □

We also need some results that go back from

z_{0} \in P (J^{(1)})

Proposition 4.6. If

z_{0} \in P_{2} (J^{(1)})

, then

z_{0} \in P_{1} (J)

Proof. By Proposition 4.2 , $z_{0}$ and $z_{0}^{- 1}$ are poles of $M^{(1)}$ , so by 4.4 , they are zeros of $M$ . As in the proof of Proposition 4.4 , if $u^{(1)}$ has a $k$ th-order zero (including $k = 0$ , i.e., $u^{(1)} (z_{0}) \neq 0$ ), then $M^{(1)} (z)$ has a $(k + 1)$ st-order pole there, and so $M (z)$ has a $(k + 1)$ st-order zero. This is only consistent with $u^{(1)} = u M$ if $u$ has a pole at $z_{0}$ . □

Proposition 4.7. If

z_{0} \in P_{1} (J^{(1)})

and

z_{0} / \in P_{1} (J)

, then

z_{0} \in P_{2} (J)

Proof. By hypothesis, $z_{0}$ is not a pole of $u$ , and so (i) of Proposition 4.2 holds. So we need only show that $M (z)$ has poles at $z_{0}$ and $z_{0}^{- 1}$ . Suppose $u$ has a $k$ th-order zero at $z_{0}$ (including $k = 0$ , i.e., $u (z_{0}) \neq 0$ ). By 4.3 and the fact that $z_{0}$ is a pole of $u^{(1)}$ , we see that $z_{0}$ is a $(k + 1)$ st-order pole of $M (z)$ and, in particular, a pole of $M (z)$ (since $k + 1 \geq 1$ ).
If $k \geq 1$ , this is only consistent with 4.5 if $\bar{u (1 / \bar{z})}$ has a zero at $z_{0}$ since the possible pole of $M$ at $z_{0}^{- 1}$ is of order $1$ and cannot cancel the $(k + 1)$ st-order pole at $z_{0}$ .
Thus, $\bar{u (1 / \overset{0}{\bar{z}})} = 0$ , so $z_{0}^{- 1}$ is real and a pole of $M (z)$ , that is, (ii) of Proposition 4.2 holds and $z_{0} \in P_{2} (J)$ .
If $k = 1$ and $M (z)$ does not have a pole at ${\bar{z}}_{0}^{- 1}$ , then $M (z) - \bar{M} (1 / \bar{z})$ has a pole at $z_{0}$ , while $\bar{u (1 / \overset{0}{\bar{z}})} \neq 0 \neq u (z_{0})$ (since $k = 1$ and $1 / {\bar{z}}_{0}$ is not a pole), violating 4.5 .
Thus, $M$ must have a pole at $z_{0}$ and $z_{0} \in P_{2} (J)$ by Proposition 4.2 . □

Proof of Theorem 4.1 . If $z_{0} \in P (J^{(1)})$ , either $z_{0} \in P_{2} (J^{(1)}) \Rightarrow z_{0} \in P_{1} (J)$ (by Proposition 4.6 ) or $z_{0} \in P_{1} (J^{(1)}) \Rightarrow z_{0} \in P_{1} (J) \cup P_{2} (J)$ (by Proposition 4.7 ). Thus, $P (J^{(1)}) \subset P (J)$ .
If $z_{0} \in P (J)$ , either $z_{0} \in P_{2} (J) \Rightarrow z_{0} \in P_{1} (J^{(1)})$ (by Proposition 4.4 ) or $z_{0} \in P_{1} (J) \Rightarrow z_{0} \in P_{1} (J^{(1)}) \cup P_{2} (J^{(1)})$ (by Proposition 4.5 ). Thus, $P (J) \subset P (J^{(1)})$ . □

References

E. Berriochoa, A. Cachafeiro, and J. García-Amor, Generalizations of the Szegő transformation interval-unit circle, preprint.
D. Damanik and B. Simon, Jost functions and Jost solutions for Jacobi matrices, I. A necessary and sufficient condition for Szegő asymptotics, preprint.
D. Damanik and B. Simon, Jost functions and Jost solutions for Jacobi matrices, II. Decay and analyticity, preprint.
P. Deift and J. Ostensson, A Riemann-Hilbert approach to some theorems on Toeplitz operators and orthogonal polynomials, in preparation.
J.S. Geronimo and K.M. Case, Scattering theory and polynomials orthogonal on the unit circle, J. Math. Phys. 20 (1979), 299–310.
Ya. L. Geronimus, On the trigonometric moment problem, Ann. of Math. (2) 47 (1946), 742–761.
Ya. L. Geronimus, Polynomials Orthogonal on a Circle and Their Applications, Amer. Math. Soc. Translation 1954 (1954), no. 104, 79 pp.
R. Killip and I. Nenciu, Matrix models for circular ensembles, Inter. Math. Res. Not. (2004), no. 50, 2665–2701.
R. Killip and B. Simon, Sum rules for Jacobi matrices and their applications to spectral theory, Ann. of Math. (2) 158 (2003), 253–321.
B. Simon, Orthogonal Polynomials on the Unit Circle, Part 1: Classical Theory, AMS Colloquium Series, American Mathematical Society, Providence, R.I., 2005.
B. Simon, Orthogonal Polynomials on the Unit Circle, Part 2: Spectral Theory, AMS Colloquium Series, American Mathematical Society, Providence, R.I., 2005.
B. Simon, Meromorphic Szegő functions and asymptotic series for Verblunsky coefficients, preprint.
G. Szegő, Über den asymptotischen Ausdruck von Polynomen, die durch eine Orthogonalitätseigenschaft definiert sind, Math. Ann. 86 (1922), 114–139.
G. Szegő, Orthogonal Polynomials, Amer. Math. Soc. Colloq. Publ., Vol. 23, American Mathematical Society, Providence, R.I., 1939; 3rd edition, 1967.

March 7, 2005

Barry Simon *