1991 Mathematics Subject Classification. Primary 47B36, 39A70; Secondary 34L25. The work was partially supported by INTAS grant no. 03-51-6637.
<ph f="cmti"> </ph><ph f="cmbx">The discrete spectrum for complex perturbations of periodic Jacobi matrices</ph>

Iryna Egorova

Leonid Golinskii

Institute for Low Temperature Physics and Engineering, 47, Lenin Ave, Kharkiv, Ukraine E-mail address : egorova@ilt.kharkov.ua Institute for Low Temperature Physics and Engineering, 47, Lenin Ave, Kharkiv, Ukraine E-mail address : golinskii@ilt.kharkov.ua

1 Introduction

The main object under investigation is an infinite complex Jacobi matrix J   , which is a compact perturbation of a real symmetric periodic matrix J 0   of period N   (background). The spectrum of the background as an operator in 2 ( N )   is known to be a union of at most N   bands divided by gaps, each of which may carry at most one eigenvalue. Since compact perturbations preserve the continuous spectrum we are primarily interested in the location of the discrete spectrum of J   . According to the general results of operator theory for small in norm perturbations the spectrum of perturbed operator is settled in a small vicinity of the spectrum of the background.
A specific feature of the given situation is that the spectrum is accumulated near the edges of bands and near the discrete spectrum of the background. Whereas a qualitative picture holds for any value of N   , we find certain quantitative bounds for the size of domain where the discrete spectrum lives only for N = 2   . In this case we also show that the discrete spectrum of J   is empty provided the spectrum of the background is pure continuous and the first moment of perturbation is small enough.
We proceed as follows. In Section 2 we present some basic facts from the spectral theory of periodic Jacobi matrices. The key tool of our investigation the discrete integral equation for the Jost solution of the Jacobi difference equation (three term recurrence relation) for the perturbed matrix J   is introduced in Section 3. In Section 4 we describe the algorithm of reconstruction of a 2-periodic matrix from its Borg sequence (lemma 4.1 ) and obtain a new representation for the Floquet solution (lemma 4.3 ). Finally, in the last Section 5 we prove the main results of the paper.
Let J = ( b 1 c 1 a 1 b 2 c 2 a 2 b 3 c 3 . . . . . . . . . ) , J 0 = ( b 1 0 a 1 0 a 1 0 b 2 0 a 2 0 a 2 0 b 3 0 a 3 0 . . . . . . . . . ) ,   be Jacobi matrices with the entries a n , b n , c n C , a n c n 0 ; a n + 2 0 = a n 0 > 0 , b n + 2 0 = b n 0 R ,   and
( J y ) n = a n 1 y n 1 + b n y n + c n y n + 1 = λ y n , λ C
( J 0 y ) n = a n 1 0 y n 1 + b n 0 y n + a n 0 y n + 1 = λ y n ,
be the Jacobi difference equations, which correspond to the matrices J   and J 0   , respectively. We assume that either the zero moment m d m   or the first moment m m d m   of the perturbation is finite, where
d m : = | b m 0 b m | + | a m 1 0 a m 1 c m 1 a m 1 0 | . (1.1)
With no loss of generality we may suppose that the continuous spectrum of the background J 0   is a set of the form σ c ( J 0 ) = [ d , s ] [ s , d ] , 0 < s < d .   Put ɛ = 1   when J 0   has an eigenvalue ν   , s < ν < s   , and ɛ = 1   otherwise (this is actually the case when ν   is an eigenvalue of the Jacobi matrix with shifted entries a ^ n 0 = a n + 1 0   , b ^ n 0 = b n + 1 0   or ν = ± s   )1 . Let w 2 ( λ )   , | w ( λ ) | 1   , be the Floquet multiplier (see (2.9 )) of J 0   . The following result holds.
Theorem 1.1. Let m = 1 d m <   , where the value d m   is given in ( 1.1 )   .
Then in the domain of variable λ   , defined by the relation
K 2 ( s , d ) | w ( λ ) | | w 2 ( λ ) w 2 ɛ ( ν ) | | 1 w 4 ( λ ) | m = 1 d m < t , (1.2)
with t 0.567 . .   being a unique positive root of the equation x exp ( x ) = 1   and K 2 ( s , d ) = 1024 2 d 5 ( d s ) 5 ( d 2 s 2 ) 1 / 2 ,   the operator J   has no discrete spectrum at all.
It can be seen from (1.2 ), that if perturbation d m   is small enough, the discrete spectrum of J   is concentrated near the edges of bands (which are exactly the points satisfying equation w 4 = 1   ) and in a vicinity of the point ν   , as long as it is an eigenvalue of J 0   (i.e., ɛ = 1   ), or ν = ± s   . We can also observe that even if ν   is not an eigenvalue of J 0   , its position relative to the edge of the gap exerts certain influence on the size of the domain where the discrete spectrum of J   is contained.
In the case when the first moment is small enough we have
Theorem 1.2. Let the background J 0   has pure continuous spectrum and a n 0 a n + 1 0   . Let K 2 ( s , d )   and t   be taken from theorem 1.1. Then under assumption K 2 ( s , d ) 1 | w ( ν ) | 2 m = 1 m d m < t ,   the complex Jacobi matrix J   has no discrete spectrum.
Our work was obviously very strongly influenced by Geronimo and Van Assche [7, wherein the behavior of the Jost solutions for real perturbations of a periodic Jacobi matrix was investigated (see also [1,Section7). More detailed analysis of such solutions in the case N = 2   is based upon the reconstruction algorithm for a periodic Jacobi matrix from its Borg sequence, suggested in [5.
In conclusion let us point out that the similar problems for the case of constant background ( a n 0 = 1   , b n 0 = 0   ) are discussed in [4(see also [6).

1 See Section 2 for more information about periodic problem

2 Preliminaries on periodic spectral problem

Let J 0 = ( b 1 0 a 1 0 a 1 0 b 2 0 a 2 0 a 2 0 b 3 0 a 3 0 . . . . . . . . . ) , a j 0 > 0 , b j 0 = b j 0 ¯   be an infinite Jacobi matrix with N   periodic entries, N 2   , a j + N 0 = a j 0 ; b j + N 0 = b j 0 ; j N = { 1 , 2 , } .   Put a 0 0 = a N 0   and consider the three term equation
a n 1 0 y n 1 + b n 0 y n + a n 0 y n + 1 = λ y n ; n N , λ C . (2.1)
The solutions s n   and c n   of ` sin   ' and ` cos   ' type, defined by the initial data s 0 = c 1 = 0   , c 1 = s 0 = 1   , are polynomials with real coefficients of degree n 1   and n 2   respectively, and also s n ( λ ) = λ n 1 a 1 0 a n 1 0 + .   The Wronskian s , c ( n ) = a n 1 0 ( s n 1 c n s n c n 1 )   of these two solutions does not depend on n   and by periodicity of a n 0   s , c ( N + 1 ) = a N 0 ( s N c N + 1 s N + 1 c N ) = s , c ( 1 ) = a N 0   or
s N + 1 ( λ ) c N ( λ ) s N ( λ ) c N + 1 ( λ ) = 1 , λ C . (2.2)
The following two polynomials with real coefficients
u ( λ ) = 1 2 ( s N + 1 ( λ ) + c N ( λ ) ) = λ N 2 a 1 0 a N 0 + (2.3)
v ( λ ) = 1 2 ( s N + 1 ( λ ) c N ( λ ) ) (2.4)
are of particular interest. The first one, u   , known as the Hill discriminant, characterizes the continuous spectrum of the operator J 0   in the space 2 ( N )   . Indeed, as is known by [2,Chapter4, the spectrum σ ( J 0 ) = σ c ( J 0 ) σ d ( J 0 )   is a union of at most N   disjoint closed intervals (bands) which constitute the continuous spectrum σ c ( J 0 )   σ c ( J 0 ) = [ μ 0 , μ N ] \ M 1 k = 1 Γ k , Γ k = ( μ k , μ k + ) , M N ,   with at most one eigenvalue inside each gap Γ k   (they form the discrete spectrum σ d ( J 0 )   ). The continuous spectrum is related to the Hill discriminant by σ c ( J 0 ) = { λ : 1 u ( λ ) 1 }   (see [2,Section7.1). Throughout the paper we assume that M = N   , i.e., an equation u 2 = 1   has no multiple zeros. Hence μ k < μ k +   and { λ : u ( λ ) = ± 1 } = { μ 0 , μ 1 , μ 1 + , . . . , μ N 1 + , μ N } .   By (2.3 ) u ( μ N ) = 1   and u +   as λ +   . The behavior of u   for λ   is determined by the parity of N   .
Consider the function ( u 2 ( λ ) 1 ) 1 / 2 = 1 2 a 1 0 a n 0 ( ( λ μ 0 ) ( λ μ N ) N 1 k = 1 ( λ μ k ) ( λ μ k + ) ) 1 / 2 ,   where the branch of the square root is taken positive for λ > μ N   . This function is real-valued inside the gaps Γ k   in such a way that sign ( u 2 1 ) 1 / 2 = ( 1 ) N k + 1 ,   λ Γ k   . On the set σ c ( J 0 )   it takes pure imaginary values.
The Weyl function m   for periodic Jacobi matrices J 0   is known to be of the form (see [3,formula(21.3.20)) m ( λ ) = ( J 0 λ I ) 1 e 1 , e 1 = v ( λ ) ( u 2 ( λ ) 1 ) 1 / 2 s N ( λ ) .   Denote by Λ r   the set of all poles of m   , which is clearly a subset of the set Λ = { λ k } k = 1 N 1   of all roots of the polynomials s N   (they all are real and simple and λ k [ μ k , μ k + ] ) .   We assume that all λ k   stay strictly inside the gaps, that is
Λ N 1 k = 1 Γ k . (2.5)
The Weyl solution ψ n   of (2.1 ) is defined in a standard way ψ n ( λ ) = c n ( λ ) + m ( λ ) s n ( λ ) ; ψ 0 = 1 , ψ 1 = m ,   { ψ n } 2 ( N )   for λ / σ c ( J 0 ) Λ r .   Along with the function m   consider the function 2 m ^ ( λ ) : = v ( λ ) + ( u 2 ( λ ) 1 ) 1 / 2 s N ( λ ) .   It is easy to see from (2.2 )-(2.4 ) that
v 2 ( λ k ) = u 2 ( λ k ) 1 , (2.6)
and so one of the functions m , m ^   has a pole at each point λ k   and the other a removable singularity .
Let us introduce the second solution of (2.1 ) ψ ^ n ( λ ) = c n ( λ ) + m ^ ( λ ) s n ( λ ) , λ / Λ \ Λ r .   As the Wronskian ψ , ψ ^ = a N 0 ( ψ 0 ψ ^ 1 ψ 1 ψ ^ 0 ) = 2 a N 0 ( u 2 ( λ ) 1 ) 1 / 2 s N 0 ,   the solutions ψ   and ψ ^   are linearly independent off the endpoints of Γ k   . Note also that ψ ^ n ( λ + i 0 ) = ψ n ( λ + i 0 ) ¯ = ψ n ( λ i 0 )   for λ σ c ( J 0 )   , and both solutions take real values for λ R \ σ ( J 0 )   . It follows from the properties of the Weyl functions m , m ^   that for n k N ,   k Z + = { 0 , 1 , }   the only one of ψ n , ψ ^ n   may have a pole at λ k   . For n = k N   they both are continuous at λ k   .
The typical feature of the periodic Jacobi difference equation is the existence of a special solution of (2.1 ) known as the Floquet solution φ n ( λ )   , which satisfies
φ n + N ( λ ) = ρ ( λ ) φ n ( λ ) , n Z + . (2.7)
Here ρ ( λ ) = u ( λ ) ( u 2 ( λ ) 1 ) 1 / 2   is the Floquet multiplier (see [1), with | ρ ( λ ) | = 1   for λ σ c ( J 0 )   and | ρ ( λ ) | < 1   for λ σ c ( J 0 )   . It is clear that { φ n } 2 ( Z + )   for λ C \ σ ( J 0 )   , and so it agrees up to a constant multiplier with the Weyl solution. We can put φ n ( λ ) ψ n ( λ )   . In addition by ([1,formula(7.19)) for the second solution
ψ ^ n + N ( λ ) = ρ 1 ( λ ) ψ ^ n ( λ ) (2.8)
holds. Following [7we introduce the function w   :
w N ( λ ) = ρ ( λ ) . (2.9)
It defines a conformal mapping of the upper C +   and lower C   half-planes onto the lower W   and upper W +   half-disks with cuts
W ± = { w : | w | < 1 ; ± Im w > 0 } \ N 1 k = 1 β k ± , (2.10)
where β k ± = [ exp ( ± i π ( N k ) N 1 ) , exp ( ± i π ( N k ) N 1 ) | w ( γ k ) | ] ,   and { γ k } k = 1 N 1   are the points of local extremum of the Hill discriminant u   , γ k ( μ k , μ k + )   (see [5,p.110). Moreover, the upper side of the interval [ μ k + , μ k + 1 ]   goes to the circular arc with endpoints exp i π ( k N ) N   and exp i π ( k + 1 N ) N ,   the lower side goes to the symmetric arc on the upper semicircle. The gap Γ k ¯   viewed as a part of the boundary of C +   is mapped on both sides of the cut β k   , in such a way that γ k   corresponds to the vertex of the cut and μ k   , μ k +   to exp ( i π ( k N ) N 1 i 0 )   , respectively. As a part of the boundary of C +   the gap Γ k ¯   goes to the symmetric cut β k +   . The infinite gaps ( , μ 0 )   and ( μ N , + )   are mapped onto ( 1 , 0 )   and ( 0 , 1 )   , respectively; w ( λ ) 0   for λ   . Hence the points μ 0 , μ 1 , . . . , μ N   correspond to w 2 N = 1   and each point from Λ   has two symmetric images on symmetric sides of the cuts.
Consider the functions ψ n ( λ )   and ψ ^ n ( λ )   as functions of variable w   for w W + W ( 1 , 0 ) ( 0 , 1 )   . Then they have continuous boundary values with
ψ n ( w ¯ ) = ψ n ( w ) ¯ , ψ ^ n ( w ¯ ) = ψ ^ n ( w ) ¯ , | w | = 1 , (2.11)
ψ n ( w ¯ ) = ψ n ( w ) R ψ ^ ( w ¯ ) = ψ ^ n ( w ) ¯ R , w β k ± .   Note that by (2.7 ),(2.8 ) and (2.9 ) the equalities
ψ n ( w ) = w n χ n ( w ) , ψ ^ n ( w ) = w n χ ^ ( w ) (2.12)
hold, where χ n ( w ) = χ N + n ( w )   , χ ^ n ( w ) = χ ^ N + n ( w )   are N   periodic functions.

2 The function m ^   appears to be the Weyl function of the operator J ^ 0 : 2 ( Z ) 2 ( Z )   with periodically extended entries a i   , b i   .

3 Integral equation

Take an infinite Jacobi matrix J = ( b 1 c 1 a 1 b 2 c 2 a 2 b 3 c 3 . . . . . . . . . ) , a j , b j , c j C , a j c j 0 ,   and the related three term recurrence relation
a n 1 y n 1 + b n y n + c n y n + 1 = λ y n ; λ C , n N . (3.1)
We assume that the matrix entries { a n }   , { b n }   and { c n }   are close to N   periodic sequences { a n 0 }   { b n 0 }   in the following sense
n = 1 { | a n a n 0 | + | b n b n 0 | + | c n a n 0 | } < . (3.2)
It is convenient to modify (3.1 ) by changing the variables v m = y m j = m a j a j 0 , m = 1 , 2 ,   (the infinite product converges due to (3.2 )). For v n   we have
a n 1 0 v n 1 + b n v n + c n a n a n 0 v n + 1 = λ v n . (3.3)
Under the Green function (for the unperturbed problem) we mean
G ( λ ; n , m ) : = { ψ n ( λ ) ψ ^ m ( λ ) ψ m ( λ ) ψ ^ n ( λ ) ψ , ψ ^ ; m > n 0 ; m n , m , n 0 , (3.4)
where the Weyl solutions ψ n   , ψ ^ n   are defined in Section 2. It is easy to check that a n 1 0 G ( λ ; n 1 , m ) + ( b n 0 λ ) G ( λ ; n , m ) + a n 0 G ( λ ; n + 1 , m ) = δ n , m ,   where δ n , m   is the Kronecker symbol.
The key role throughout the rest of the paper is played by the discrete integral equation (cf. [1,formulae(7.87)-(7.88)
v n ( λ ) = ψ n ( λ ) + m = n + 1 A ( λ ; n , m ) v n ( λ ) , λ / Λ r . (3.5)
with the kernel
A ( λ ; n , m ) = ( b m 0 b m ) G ( λ ; n , m ) + ( a m 1 0 a m 1 c m 1 a m 1 0 ) G ( λ ; n , m 1 ) . (3.6)
Let us show that each solution of (3.5 ) satisfies (3.3 ) and has the following asymptotic behavior for n   :
v n ψ n 0 .   Indeed, by the definition
A ( λ ; n 1 , m ) = b n 0 b n a n 1 0 ,
A ( λ ; n 1 , n + 1 ) = ( b n + 1 0 b n + 1 ) ( λ b n 0 ) a n 1 0 a n 0 + ( a n 0 a n c n a n 0 ) .
We apply the `unperturbed' operator to both parts of (3.5 )
( J 0 λ ) v n = a n 1 0 v n 1 + ( b n 0 λ ) v n + a n 0 v n + 1
= ( J 0 λ ) ψ n ( λ ) + a n 1 0 { A ( λ ; n 1 , n ) v n + A ( λ ; n 1 , n + 1 ) v n + 1 }
+ ( b n 0 λ ) A ( λ ; n , n + 1 ) v n + 1 + m = n + 2 ( J 0 λ ) n A ( λ ; n , m ) v m
or a n 1 0 v n 1 + ( b n 0 λ ) v n + a n 0 v n + 1 = ( b n 0 b n ) v n + ( a n 0 a n c n a n 0 ) v n + 1 ,   which gives exactly (3.3 ).
Multiply (3.5 ) through by w n   , where w   is defined in (2.9 ):
v ~ n : = v n w n = ψ n ( λ ) w n + m = n + 1 A ~ ( λ ; n , m ) v ~ m , (3.7)
A ~ ( λ ; n , m ) : = A ( λ ; n , m ) w m n ,   or
V n ( λ ) = m = n + 1 A ~ ( λ ; n , m ) ψ m ( λ ) w m + m = n + 1 A ~ ( λ ; n , m ) V m ( λ ) , (3.8)
with V n ( λ ) = ( v n ( λ ) ψ n ( λ ) ) w n .   Equation (3.8 ) is analyzed thoroughly in [7,Theorem3, for real symmetric Jacobi matrices ( c n = a n > 0 ,   b n = b n ¯ )   . The existence of the solution V   such that | V n ( λ ) s N ( λ ) | = O ( 1 ) , n   uniformly on compact subsets of the domain C \ { μ k ± }   has been established (the reason for the factor s N   to enter the asymptotic relation is to `suppress' poles of ψ n   ).
We will be interested in the bounds of the form | v n ( λ ) ψ n ( λ ) | K ( λ ) w n κ 0 ( n ) | 1 w 2 N | , λ / Λ r ,   wherein a positive function K ( λ )   is given explicitly in terms of the spectral data of the background and uniformly bounded off certain small neighborhoods of the set Λ r   ,
κ 0 ( n ) : = m = n + 1 d m , d m : = | b m 0 b m | + | a m 1 0 a m 1 c m 1 a m 1 0 | . (3.9)
We solve this problem completely for the case N = 2   , that is performed later on.
In the rest of this section we show how the method of successive approximations works in connection with equation (3.8 ). Appropriate bounds (see lemma 5.2 ) for the kernel A ~   lead to the uniform convergence of the series
V n , 1 ( λ ) : = m = n + 1 A ~ ( λ ; n , m ) ψ m ( λ ) w m (3.10)
on compact subsets of the domain C \ ( { μ k ± } Λ r )   , and moreover | V n , 1 ( λ ) | K ( λ ) | 1 w 2 N | κ 0 ( n ) .   Put V n , j + 1 : = m = n + 1 A ~ ( λ ; n , m ) V m , j ( λ )   and prove by induction that | V n , j ( λ ) | 1 ( j 1 ) ! ( K ( λ ) κ 0 ( n ) | 1 w 2 N | ) j .   Then the required solution of (3.8 ) is given by
V n ( λ ) = j = 1 V n , j ( λ ) , | V n ( λ ) | K ( λ ) κ 0 ( n ) | 1 w 2 N | exp ( K ( λ ) κ 0 ( n ) | 1 w 2 N | ) . (3.11)
In particular for n = 0   thanks to ψ 0 = 1   we have | v 0 ( λ ) 1 | K 1 ( w ) m = 1 d m | 1 w 2 N | exp ( K 1 ( w ) m = 1 d m | 1 w 2 N | ) ,   which makes it possible to describe the region of the variable w   wherein v 0 0   .
The similar argument (with slightly different bounds) holds under the assumption
m = 1 m d m < . (3.12)
It remains only to point out that the set of zeros of the function v 0   agrees with the set of eigenvalues (discrete spectrum) of the operator J   off σ c ( J 0 )   .

4 Reconstruction of matrix of period 2 from spectral data

Let J 0   be a real symmetric Jacobi matrix of period 2 with the matrix entries a j + 2 0 = a j 0   , b j + 2 0 = b j 0   . The first few polynomials s n   , c n   can be easily computed from recurrence relation (2.1 ):
{ s 0 = 0 , s 1 = 1 , s 2 = λ b 1 0 a 1 0 , s 3 = ( λ b 1 0 ) ( λ b 2 0 ) ( a 1 0 ) 2 a 1 0 a 2 0 , c 0 = 1 , c 1 = 0 , c 2 = a 2 0 a 1 0 , c 3 = b 2 0 λ a 1 0 . (4.1)
The Hill discriminant
u ( λ ) = 1 2 ( s 3 ( λ ) + c 2 ( λ ) ) = λ 2 ( b 1 0 + b 2 0 ) λ + b 1 0 b 2 0 ( a 2 0 ) 2 ( a 1 0 ) 2 2 a 1 0 a 2 0 (4.2)
is a quadratic parabola with the positive leading coefficient. Hence the bands of continuous spectrum { λ : 1 u 1 }   are symmetric with respect to the vertex of this parabola and have the same length. By shifting the spectral parameter we can put the vertex to the origin. Denote the endpoints of the spectrum by { d , s , s , d } ,   0 < s < d   . Then
u ( λ ) 1 = C ( λ 2 d 2 ) , u ( λ ) + 1 = C ( λ 2 s 2 ) , u 2 ( λ ) 1 = C 2 ( λ 2 d 2 ) ( λ 2 s 2 ) , (4.3)
whence it follows that C = 2 ( d 2 s 2 ) 1   and
u ( λ ) = 2 d 2 s 2 ( λ 2 d 2 + s 2 2 ) . (4.4)
Let ν   be a root of the polynomial s 2 ,   ν = b 1 0   . By the assumption (2.5 ) we have | ν | < s   . Put ɛ = 1   in the case when the Weyl function m ( λ )   has a pole at ν   and ɛ = 1   otherwise. Our goal here is to restore the matrix J 0   from its Borg sequence { ± s , ± d , ν , ɛ }   and to derive a formula for the Weyl function as a function of variable w   (2.9 ), which is related to u   by the equality
u ( λ ) = w 2 ( λ ) + w 2 ( λ ) 2 . (4.5)
By comparing (4.2 ) and (4.4 ) we come to the relations for the matrix entries:
b 2 0 = b 1 0 = ν (4.6)
and
a 1 0 a 2 0 = d 2 s 2 4 , ν 2 + ( a 1 0 ) 2 + ( a 2 0 ) 2 = d 2 + s 2 2 , (4.7)
so that ( a 1 0 ) 2   , ( a 2 0 ) 2   are roots of the quadratic equation x 2 ( d 2 + s 2 2 ν 2 ) x + ( d 2 s 2 4 ) 2 = 0 .   Hence
4 ( a min 0 ) 2 = d 2 + s 2 2 ν 2 2 ( d 2 ν 2 ) ( s 2 ν 2 ) ) + 1 / 2 ,
4 ( a max 0 ) 2 = d 2 + s 2 2 ν 2 + 2 ( ( d 2 ν 2 ) ( s 2 ν 2 ) ) + 1 / 2 , (4.8)
where a min 0 = min { a 1 0 , a 2 0 }   , a max 0 = max { a 1 0 , a 2 0 }   and ( ) + 1 / 2   designates the arithmetic value of the square root.
Let us show that ɛ = 1 ( ɛ = 1 )   is equivalent to a 1 0 < a 2 0   ( a 2 0 < a 1 0 )   . To this end write
v ( λ ) = 1 2 ( s 3 ( λ ) c 2 ( λ ) ) = λ 2 ν 2 ( a 1 0 ) 2 + ( a 2 0 ) 2 2 a 1 0 a 2 0
= 2 d 2 s 2 ( λ 2 ν 2 ( a 1 0 ) 2 + ( a 2 0 ) 2 ) ,
and so
v ( ν ) = 2 d 2 s 2 ( ( a 2 0 ) 2 ( a 1 0 ) 2 ) . (4.9)
With regard to (4.8 ) and (4.3 ) (see also (2.6 )) the latter gives v 2 ( ν ) = 4 ( d 2 s 2 ) 2 ( d 2 ν 2 ) ( s 2 ν 2 ) = u 2 ( ν ) 1 ,   and therefore | v ( ν ) | = ( u 2 ( ν ) 1 ) + 1 / 2   . In two-band case the Weyl function admits representation m ( λ ) = v ( λ ) + ( u 2 ( λ ) 1 ) + 1 / 2 s 2 ( λ ) , | λ | < s ,   and has a pole if and only if | v ( ν ) | = v ( ν )   , that is equivalent to a 1 0 < a 2 0   by (4.9 ).
Thus we come to the following
Lemma 4.1. The entries of a two-periodic Jacobi matrix can be restored from the Borg sequence { d , s , s , d , ν , ɛ }   by the formulae:
b 2 0 = b 1 0 = ν ,   4 ( a 1 0 ) 2 = d 2 + s 2 2 ν 2 2 ɛ ( ( d 2 ν 2 ) ( s 2 ν 2 ) ) + 1 / 2 ,   4 ( a 2 0 ) 2 = d 2 + s 2 2 ν 2 + 2 ɛ ( ( d 2 ν 2 ) ( s 2 ν 2 ) ) + 1 / 2 .  
Corollary 4.2. The following inequalities hold
d s 2 a j 0 d + s 2 . (4.10)
Our next problem is to find effective bounds for the Weyl solution, more precisely, for the function χ n ( λ ) = ψ n ( λ ) w n   in terms of the variable w   . As χ n + 2 = χ n   , χ 0 = ψ 0 = 1   , we only need estimate the function χ 1 ( λ ) = w 1 ( λ ) ψ 1 ( λ )   . By the definition of w   and in view of (4.5 ) we have
w 2 ( ν ) = u ( ν ) ( u 2 ( ν ) 1 ) 1 / 2 . (4.11)
Remind that in the two-band case ( u 2 ( ν ) 1 ) 1 / 2 < 0   .
Lemma 4.3. The following representation holds
χ 1 ( λ ) = 4 a 1 0 ( λ + ν ) w ( λ ) d 2 s 2 1 w 2 ( λ ) w 2 ɛ ( ν ) . (4.12)
  • Proof. The wronskian of two solutions s n ( λ )   and ψ n ( λ )   does not depend on n   so ψ , s = a 2 0 ( ψ 2 s 3 s 2 ψ 3 ) = a 0 0 ( ψ 0 s 1 ψ 1 s 0 ) = a 0 0 : = a 2 0 ,   that is
    s 3 ( λ ) ψ 2 ( λ ) s 2 ( λ ) ψ 3 ( λ ) = 1 . (4.13)
    Since (see Section 1) ψ 2 = w 2 ψ 0 = w 2   , ψ 3 = w 2 ψ 1   , it follows from (4.13 ), (2.1 ), (4.1 ) and theorem 4.1 that
    ψ 1 = s 3 w 2 1 w 2 s 2 = ( ( λ b 2 0 ) s 2 a 1 0 ) w 2 ( a 2 0 ) 1 1 w 2 s 2
    = λ + ν a 2 0 a 1 0 w 2 + a 2 0 a 2 0 w 2 s 2 = λ + ν a 2 0 a 1 0 ψ 2 + a 0 0 ψ 0 a 2 0 w 2 s 2
    = λ + ν a 2 0 ( λ b 1 0 ) ψ 1 a 2 0 w 2 s 2 = λ + ν a 2 0 a 1 0 a 2 0 w 2 ψ 1 ,
    i.e., ψ 1 ( 1 + a 1 0 w 2 a 2 0 ) = λ + ν a 2 0 ,   and hence
    ψ 1 ( λ ) = ( λ + ν ) w 2 a 2 0 w 2 + a 1 0 = 4 a 1 0 ( λ + ν ) w 2 4 a 1 0 a 2 0 w 2 + 4 ( a 1 0 ) 2 . (4.14)
    But by (4.3 ), (4.4 ), (4.7 ) and theorem 4.1 4 ( a 1 0 ) 2 d 2 s 2 = 2 d 2 s 2 ( d 2 + s 2 2 ν 2 ) + ɛ 2 d 2 s 2 ( ( d 2 ν 2 ) ( s 2 ν 2 ) ) + 1 / 2 ,   that is 4 ( a 1 0 ) 2 d 2 s 2 = u ( ν ) ɛ ( u 2 ( ν ) 1 ) 1 / 2 .   Relations (4.11 ), (4.5 ) imply
    4 ( a 1 0 ) 2 d 2 s 2 = ( w ( ν ) ) 2 ɛ . (4.15)
    Next we plague (4.15 ) into (4.14 ) and take into account (4.7 ):
    ψ 1 ( λ ) = 4 a 1 0 ( λ + ν ) w 2 ( λ ) ( d 2 s 2 ) ( w 2 ( λ ) w 2 ɛ ( ν ) ) , (4.16)
    which proves (4.12 ) since χ 1 ( λ ) = w 1 ( λ ) ψ 1 ( λ )   .

5 Bounds for solutions of integral equation

The bounds in the method of successive approximations displayed in Section 3 for the case N = 2   is based upon explicit expressions for w   and the Weyl functions m , m ^   . We begin with the following simple result.
Lemma 5.1. Let 0 | ν | < s < d   and ɛ = ± 1   . Put
τ ( ν , ɛ ) : = ɛ ( ( d 2 ν 2 ) ( s 2 ν 2 ) ) 1 / 2 ν 2 d 2 , τ = τ ( 0 , ɛ ) = ɛ s d , ξ = d 2 + s 2 2 d 2 . (5.1)
For | z | < 1   consider the functions f ( z , ν , ɛ ) : = 1 + τ ( ν , ɛ ) z + ( ( 1 z ) ( 1 τ z ) ) 1 / 2 ,   g ( z ) : = 1 ξ z + ( 1 z ) ( 1 τ z ) ) 1 / 2 ,   where the branch of the square root is taken positive for 1 < z < 1   .
Then
| f ( z , ν , ɛ ) | ( 1 + τ ( ν , ɛ ) ) 2 4 ( 1 | τ | ) 2 4 = ( d s ) 2 4 d 2 , (5.2)
| g ( z ) | | ξ 2 τ 2 | 4 = ( 1 τ 2 ) 2 4 = ( d 2 s 2 ) 2 4 d 4 . (5.3)
  • Proof. We have f ( z , ν , ɛ ) = f 1 + f 2   , where f 1 2 f 2 2 = ( 1 + τ ( ν , ɛ ) z ) 2 ( 1 z ) ( 1 τ 2 z ) = ( 2 τ ( ν , ɛ ) + τ 2 + 1 ) z + ( τ 2 ( ν , ɛ ) τ 2 ) z 2 ,   and so f ( z , ν , ɛ ) = 2 τ ( ν , ɛ ) + τ 2 + 1 + ( τ 2 ( ν , ɛ ) τ 2 ) z ( f 1 ( z ) f 2 ( z ) ) z 1 .   The function h = f 1 f 2   is analytic in D   , h ( 0 ) = 0   and | h ( z ) | 1 + | τ ( ν , ɛ ) | + ( 2 ( 1 + τ 2 ) ) 1 / 2 .   An elementary analysis of the function τ ( ν , ɛ )   shows that
    | τ ( ν , ɛ ) | | τ | < 1 , (5.4)
    i.e., | h | 4   . By Schwarz's Lemma | h | 4 | z |   for | z | 1   and hence the second inequality in (5.2 ) follows from (5.4 ).
    The below bound for the function g   comes out in a similar manner: g = g 1 + g 2 , g 1 2 g 2 2 = ( 1 ξ z ) 2 ( 1 z ) ( 1 τ 2 z ) = ( ξ 2 τ 2 ) z 2 ,   by Schwarz's Lemma | g 1 g 2 | 4 | z | 2   , which leads to (5.3 ).
Let us next turn to the kernel A ~   (3.6 ), (3.7 ).
Lemma 5.2. The following inequalities take place
| A ~ ( λ ; n , m ) | K 1 ( s , d ) | w ( λ ) | | 1 w 4 ( λ ) | d m , λ C , n , m Z + , (5.5)
where
K 1 ( s , d ) : = 256 d 3 ( d s ) 4 ; d m = | b m b m 0 | + | a m 1 0 a m 1 c m 1 a m 1 0 | . (5.6)
  • Proof. Consider the expression
    G ( λ ; n , m ) w m n = χ n ( λ ) χ ^ m ( λ ) χ ^ n ( λ ) χ m ( λ ) w 2 ( m n ) ψ , ψ ^ . (5.7)
    For the denominator in (5.7 ) we have
    ψ , ψ ^ = 2 a 2 0 s 2 ( u 2 ( λ ) 1 ) 1 / 2 = 2 a 1 0 a 2 0 λ ν ( u 2 ( λ ) 1 ) 1 / 2
    = d 2 s 2 4 ( λ ν ) ( w 2 ( λ ) w 2 ( λ ) ) .
    The functions χ n , χ ^ n   in the numerator are 2-periodic, and also χ 2 n ( λ ) = χ ^ 2 n = 1 , χ 2 n + 1 ( λ ) = m ( λ ) w 1 ( λ ) , χ ^ 2 n + 1 ( λ ) = m ^ ( λ ) w ( λ ) , n Z + .   Thereby we have 4 different expressions for the right hand side (5.7 ):
    G ( λ ; 2 p , 2 q ) w 2 ( q p ) = D ( λ ν ) w 2 ( 1 w 4 ( q p ) ) 1 w 4 , (5.8)
    G ( λ ; 2 p , 2 q + 1 ) w 2 ( q p ) + 1 = D m ^ ( λ ) ( λ ν ) w 2 m ( λ ) ( λ ν ) w 4 ( q p ) + 3 1 w 4 , (5.9)
    G ( λ ; 2 p + 1 , 2 q ) w 2 ( q p ) 1 = D m ( λ ) ( λ ν ) w m ^ ( λ ) ( λ ν ) w 4 ( q p ) + 1 1 w 4 , (5.10)
    G ( λ ; 2 p + 1 , 2 q + 1 ) w 2 ( q p ) = D m ( λ ) m ^ ( λ ) ( λ ν ) w 2 ( 1 w 4 ( q p ) ) 1 w 4 (5.11)
    with D = 4 d 2 s 2 .   We can estimate right hand sides in (5.8 ) -(5.11 ) by using the explicit formulae for w   , m   , m ^   . Assume first that | λ | d   . Then
    w 2 ( λ ) = u ( λ ) ( u 2 ( λ ) 1 ) 1 / 2 = 1 u ( λ ) + ( u 2 ( λ ) 1 ) 1 / 2
    = d 2 s 2 2 1 λ 2 ( s 2 + d 2 ) / 2 + ( ( λ 2 s 2 ) ( λ 2 d 2 ) ) 1 / 2 ,
    so that ( λ ν ) 2 w 2 ( λ ) = d 2 s 2 2 g ( z ) ( 1 ν λ ) 2 , z = d 2 λ 2 .   By lemma 5.1
    | ( λ ν ) w ( λ ) | 2 d 2 ( d 2 s 2 ) 1 / 2 ( 1 + | ν | d ) , | λ | d . (5.12)
    The Weyl function is easily seen to satisfy as z = d 2 λ 2  
    ( λ ν ) m ( λ ) = 2 a 1 0 d 2 s 2 ( λ 2 + d 2 τ ( ν , ɛ ) ( ( λ 2 s 2 ) ( λ 2 d 2 ) ) 1 / 2 )
    = 2 a 1 0 d 2 s 2 2 d 2 τ ( ν , ɛ ) + s 2 + d 2 + ( d 4 τ 2 ( ν , ɛ ) s 2 d 2 ) λ 2 f ( z , ν , ɛ ) , .
    By the definition of τ ( ν , ɛ )   (5.1 ) and the property (5.4 ) it is not hard to derive | 2 d 2 τ ( ν , ɛ ) + s 2 + d 2 | ( s + d ) 2 , | d 4 τ 2 ( ν , ɛ ) s 2 d 2 | 3 ν 2 ( s + d ) 2 ,   and so the application of lemma 5.1 gives
    | m ( λ ) ( λ ν ) | 8 a 1 0 ( s + d ) 2 d 2 s 2 4 d 2 ( d s ) 2 16 ( s + d ) 2 d 2 ( d s ) 3 , | λ | d (5.13)
    (in the later inequality we use 2 a 1 0 s + d   (see (4.10 )). In the same way
    m ^ ( λ ) ( λ ν ) w 2 ( λ ) = 2 a 1 0 d 2 s 2 λ 2 + d 2 τ ( ν , ɛ ) + ( λ 2 s 2 ) ( λ 2 d 2 ) ) 1 / 2 λ 2 d 2 ξ + ( ( λ 2 s 2 ) ( λ 2 d 2 ) ) 1 / 2 d 2 s 2 2
    = a 1 0 f ( z , ν , ɛ ) g ( z ) , z = d 2 λ 2 ,
    and by lemma 5.1
    | m ^ ( λ ) ( λ ν ) w 2 ( λ ) | 16 a 1 0 d 4 ( d 2 s 2 ) 2 8 d 4 ( d s ) ( d 2 s 2 ) , | λ | d . (5.14)
    Finally m ( λ ) m ^ ( λ ) = c 3 ( λ ) s 2 ( λ ) = λ + ν λ ν , m ( λ ) m ^ ( λ ) ( λ ν ) w ( λ ) = ( λ + ν ) w ( λ )   and as above in (5.12 )
    | m ( λ ) m ^ ( λ ) ( λ ν ) w ( λ ) | 2 d 2 ( d 2 s 2 ) 1 / 2 ( 1 + | ν | d ) , | λ | d . (5.15)
    For | λ | d   the argument is much more elementary and is based on | w | 1   :
    | ( λ ± ν ) w | d + | ν | , (5.16)
    | m ( λ ) ( λ ν ) | 2 a 1 0 d 2 s 2 ( d 2 + d 2 | τ | + 2 d 2 ) 4 d 2 d s , (5.17)
    and the same estimate holds for m ^ ( λ ) ( λ ν )   .
    It is clear that the worst bound comes in (5.13 ) and thus we arrive at the following inequalities
    | G ( λ ; n , m ) w m n | K 1 ( s , d ) | w ( λ ) | | 1 w 4 ( λ ) | , λ C , n , m Z + . (5.18)
    The desired result now stems from (5.18 ) and the definition of the kernel A ~   .
An application of (5.12 ) and (4.10 ) to (4.12 ) gives
| χ 1 ( λ ) | 2 ( d + s ) d 2 s 2 2 2 d 2 ( d 2 s 2 ) 1 / 2 1 | w 2 ( λ ) w 2 ɛ ( ν ) |
= 4 2 d 2 ( d s ) ( d 2 s 2 ) 1 / 2 1 | w 2 ( λ ) w 2 ɛ ( ν ) | , (5.19)
and in fact this bound is true for all χ n   .
We are now in a position to estimate the solution of equation (3.8 )
V n ( λ ) = m = n + 1 A ~ ( λ ; n , m ) χ m ( λ ) + m = n + 1 A ~ ( λ ; n , m ) V m ( λ ) (5.20)
with constants expressed explicitly in terms of the spectral data for the 2-periodic background. Assume first that
m = 1 d m < . (5.21)
By (5.19 ) and lemma 5.2 the first series in (5.20 ) converges uniformly on compact subsets of the domain
Ω ( ɛ ) = { C \ { ± s , ± d } , ɛ = 1 , C \ { ± s , ± d , ν } , ɛ = 1 , (5.22)
and its sum is subject to | V n , 1 ( λ ) | K 2 ( s , d ) κ 0 ( n ) | w ( λ ) | | w 2 ( λ ) w 2 ɛ ( ν ) | | 1 w 4 ( λ ) | ,  
K 2 ( s , d ) = K 1 ( s , d ) 4 2 d 2 ( d s ) ( d 2 s 2 ) 1 / 2 , (5.23)
κ 0 ( n )   is defined in (3.9 ), κ ( n ) = o ( 1 )   as n   . Hence for V n   we have (cf.(3.11 ))
| V n ( λ ) | K ( λ ) κ 0 ( n ) | 1 w 4 ( λ ) | exp ( K ( λ ) κ 0 ( n ) | 1 w 4 ( λ ) | ) , K ( λ ) = K 2 ( s , d ) | w ( λ ) | | w 2 ( λ ) w 2 ɛ ( ν ) | , (5.24)
which leads to the following result.
Theorem 5.3. In assumption ( 5.21 )   there exists a solution v n   of equation ( 3.3 )   which satisfies
| v n ( λ ) ψ n ( λ ) | K ( λ ) κ 0 ( n ) | w n ( λ ) | | 1 w 4 ( λ ) | exp ( K ( λ ) κ 0 ( n ) | 1 w 4 ( λ ) | ) , λ Ω ( ɛ ) , (5.25)
where K ( λ )   is given in ( 5.24 )   , ( 5.23 )   and ( 5.6 )   , κ 0 ( n )   is defined in ( 3.9 )   .
 
Inequality (5.25 ) with n = 0   implies | v 0 ( λ ) 1 | K ( λ ) | 1 w 4 ( λ ) | m = 1 d m exp ( K ( λ ) | 1 w 4 ( λ ) | m = 1 d m ) .   Taking into account the spectral interpretation of zeros of the Jost function v 0   we can draw the following conclusion.
Corollary 5.4. Define the domain G   in λ   -plain by inequality
G : = { λ C : K 2 ( s , d ) | w ( λ ) | | w 2 ( λ ) w 2 ɛ ( ν ) | | 1 w 4 ( λ ) | m = 1 d m < t } , (5.26)
where t 0.567 . .   is a unique positive root of the equation x exp ( x ) = 1   .
Then G   is free from the discrete spectrum of J   .  
In the case when the first moment of the perturbation is finite
m = 1 m d m < , (5.27)
the bounds for the Green kernel are performed in a slightly different way by using a simple inequality | 1 w n 1 w | n , | w | 1 .   For the right hand sides in (5.8 ) and (5.11 ) we see that
| G ( λ ; 2 p , 2 q ) w 2 ( q p ) | 4 ( q p ) d 2 s 2 | ( λ ν ) w ( λ ) | | w | ,
| G ( λ ; 2 p + 1 , 2 q + 1 ) w 2 ( q p ) | 4 ( q p ) d 2 s 2 | ( λ + ν ) w ( λ ) | | w | .
Next, for (5.9 ) G ( λ ; 2 p , 2 q + 1 ) w 2 ( q p ) + 1 = G + G ,   where
G : = 4 ( λ ν ) ( d 2 s 2 ) ( w 2 ( λ ) w 2 ( λ ) ) ( m ^ ( λ ) m ( λ ) ) w ( λ )
= s 2 2 a 2 0 ( u 2 1 ) 1 / 2 ± 2 ( u 2 1 ) 1 / 2 s 2 w ( λ ) ,
and so | G | | w ( λ ) | a 2 0 < 2 | w ( λ ) | d s ;   G : = 4 ( λ ν ) d 2 s 2 m ( λ ) w ( λ ) ( 1 w 4 ( p q ) ) w 2 ( λ ) w 2 ( λ ) ,   and | G | 4 ( q p ) d 2 s 2 | ( λ ν ) m w 3 | .   Hence | G ( λ ; 2 p , 2 q + 1 ) w 2 ( q p ) + 1 | 2 | w | d s + 4 | ( λ ν ) m w 3 | d 2 s 2 ( q p ) .   The similar inequality holds for G ( λ ; 2 p + 1 , 2 q ) w 2 q 2 p 1   . Finally, | A ~ ( λ ; n , m ) | K 1 ( s , d ) ( m n ) d m | w ( λ ) | .   The method of successive approximations of Section 3 is applicable and provides the bound for the solution of the basic integral equation in assumption (5.27 ):
| v n ( λ ) ψ n ( λ ) | K ( λ ) κ 1 ( n ) exp ( K ( λ ) κ 1 ( n ) ) ,   where K ( λ )   is defined in (5.24 ), κ 1 ( n ) : = m = n + 1 ( m n ) d m .   We are aimed here at the condition which guarantees the lack of discrete spectrum for the operator J   . That is why it seems reasonable to suppose that the background operator J 0   has pure continuous spectrum, i.e., in our notation ɛ = 1   . Then 1 | w 2 ( λ ) w 2 ( ν ) | = | w 2 ( ν ) | | 1 w 2 ( ν ) w 2 ( λ ) | 1 1 | w 2 ( ν ) | ,   and hence | v 0 ( λ ) 1 | K 2 ( s , d ) 1 | w ( ν ) | 2 m = 1 m d m exp ( K 2 ( s , d ) 1 | w ( ν ) | 2 m = 1 m d m ) ,   that gives the following result.
Theorem 5.5. The operator J   has no discrete spectrum as long as K 2 ( s , d ) 1 | w ( ν ) | 2 m = 1 m d m < t .  
References

  1. G. Teschl, Jacobi Operators and Completely Integrable Nonlinear Lattices Math. Surv. and Mon.,Vol. 72, AMS, Rhode Island, 2000.
  2. M. Toda, Theory of Nonlinear Lattices, Springer, Berlin, 1981.
  3. E. Titchmarsh, Eigenfunction Expansions Associated with Second-Order Differential Equations, Part II, Oxford University Press, Oxford, 1958.
  4. I. Egorova, L. Golinskii, On location of discrete spectrum for complex Jacobi matrices, to appear in Proc. AMS
  5. L. Percolab, Inverse problem for periodic Jacobi matrix,Teor. funk., funk. anal., pril., 42 (1984), 107-121.
  6. V. Lyantse, Spectrum and resolvent of a non-selfconjugate difference operator,Ukr. Math. J., 20, no.4 (1968), 489-503.
  7. J. Geronimo, W. Van Assche, Orthogonal polynomials with asymptotically periodic recurrence coefficients, J.of Approx. Theory , 46 (1986), 251-283.

Institute for Low Temperature Physics and Engineering, 47, Lenin Ave, Kharkiv, Ukraine E-mail address : egorova@ilt.kharkov.ua Institute for Low Temperature Physics and Engineering, 47, Lenin Ave, Kharkiv, Ukraine E-mail address : golinskii@ilt.kharkov.ua