, , , 25.05.2004

In this paper we study the spectrum of a family

A_{α}

of differential operators in the space

L^{2} (R^{2})

. Each operator

A_{α}

is defined by the same differential expression

\begin{matrix} A U = - U_{x x}^{''} + \frac{1}{2} (- U_{q q}^{''} + q^{2} U) . \end{matrix}

(1.1)

The parameter

α \in R

appears in the “transmission condition” across the line

x = 0

\begin{matrix} U_{x}^{'} (0 +, q) - U_{x}^{'} (0 -, q) = α q U (0, q), q \in R . \end{matrix}

(1.2)

As we shall show in Theorem 4.1 , for any

α

the operator

A_{α}

has a unique natural self-adjoint realization. The replacement

α \mapsto - α

corresponds to the change of variables

q \mapsto - q

which does not affect the spectrum.

For this reason, below we shall discuss only non-negative

α

The family

A_{α}

was suggested by the physicist Smilansky in [12] as a model of an irreversible quantum system. He carried out a formal computation of the scattering matrix for the pair

(A_{0}, A_{α})

and showed that this matrix is unitary only if

α < \sqrt{2}

. The loss of unitarity of the scattering matrix for large values of

α

was interpreted in [12] as irreversibility of the system.

First rigorous mathematical results on the family

A_{α}

were obtained in the paper [13] , inspired by [12] . The present paper can be considered as the second part of [13] , but it can be read independently.

The family

A_{α}

exhibits many unusual features, partly revealed in [12] and [13] . The most important of them is a phase transition at the point

α = \sqrt{2}

: the spectral properties of

A_{α}

for

α < \sqrt{2}

and for

α > \sqrt{2}

are quite different. In what follows we refer to the values

α < \sqrt{2}

as “small” and to

α > \sqrt{2}

as “large”.

The spectrum of the operator

A_{0}

can be easily described via separation of variables. It is absolutely continuous, fills the half-line

[1 / 2, \infty)

, and its multiplicity function is given by eq. 2.4 in Section 2 . It seems natural to study the spectrum

σ (A_{α})

for

α > 0

with the help of the perturbation theory of quadratic forms. The standard assumption in this type of problems is relative compactness (see [11] ) of the perturbation with respect to the unperturbed quadratic form, i.e. to the one of the operator

A_{0}

. However, in our case this property is violated: the perturbation is only form-bounded but not form-compact. This was shown in [13] .

As a rule, for the form-bounded perturbations the quadratic form approach does not give much information. Nevertheless, a rather complete description of the essential spectrum and of the point spectrum of

A_{α}

for small

α

was obtained in [13] and [14] by means of this approach. In particular, it was shown in [14] that for any

α < \sqrt{2}

the spectrum of the operator

A_{α}

below the threshold

1 / 2

consists of a finite number of eigenvalues. This number grows indefinitely as

α ↗ \sqrt{2}

and has regular asymptotics of the order

O ((\sqrt{2} - α)^{- 1 / 2})

. If

α \geq \sqrt{2}

, then the point spectrum of

A_{α}

is empty (see Theorem 5.2 ).

Our goal in this paper is to study the absolutely continuous spectrum

σ_{a . c .}

of the operators

A_{α}

for all values of the parameter

α

. This problem was not dealt with in [13] . A certain Jacobi operator in the space

ℓ^{2}

is involved in the description of

σ_{a . c .} (A_{α})

, namely

J_{0} (μ) : {C_{n}} \mapsto {d_{n + 1} C_{n + 1} + (2 n + 1) μ C_{n} + d_{n} C_{n - 1}}, d_{n} = n^{1 / 2} (n^{2} - 1 / 4)^{1 / 4} .

Here

μ > 0

is an auxiliary parameter. It is often convenient to use it along with

α

. Our main result, Theorem 5.1 , states that

\begin{matrix} \begin{matrix} (1.3) & σ_{a . c .} (A_{α}) = σ_{a . c .} (A_{0}) \cup σ_{a . c .} (J_{0} (\sqrt{2} / α)), m_{a . c .} (E; A_{α}) = m_{a . c .} (E; A_{0}) + m_{a . c .} (E; J_{0} (\sqrt{2} / α)), a . e . E \in R . \end{matrix} \end{matrix}

(1.4)

Here the symbol

m_{a . c .}

stands for the multiplicity function of the absolutely continuous spectrum.

The spectrum of

J_{0} (μ)

is discrete for

μ > 1

and is purely absolutely continuous for

μ \leq 1

. Moreover, we show in Theorem 3.1 that

σ (J_{0} (μ)) = (- \infty, \infty) for μ < 1, σ (J_{0} (1)) = [0, \infty); m_{a . c .} (E; J_{0} (μ)) = 1 a . e . on σ (J_{0} (μ)) .

Thus, the equalities 1.3 and 1.4 give the complete description of the absolutely continuous spectrum of the operators

A_{α}

. Namely, for small

α

it coincides with

σ_{a . c .} (A_{0})

, including equality of the multiplicities.

For large

α

a new branch of the absolutely continuous spectrum of multiplicity

1

adds to

σ_{a . c .} (A_{0})

, its source is the Jacobi matrix

J_{0} (μ)

This birth of an additional branch of the a.c. spectrum is the exact mathematical expression of the effect which was interpreted in [12] as irreversibility.

The family

A_{α}

is a striking example of a problem which, in spite of its seeming simplicity, exhibits many unexpected effects. This refers to both the point spectrum and the absolutely continuous spectrum.

For this reason we believe that the detailed analysis of this family is of general interest. Note that one important question remains unanswered.

Namely, our method does not check whether

σ (A_{α})

has singular continuous component.

In the course of the proof of the equalities 1.3 and 1.4 we use the tools coming from different parts of the spectral theory. In Section 6 we obtain a convenient representation of the operator

\begin{matrix} (A_{α} - Λ)^{- 1} - (A_{0} - Λ)^{- 1} . \end{matrix}

(1.5)

This representation involves some matrix-valued function

J_{0} (Λ; μ)

which arises in a natural way when looking for formal (i.e. not necessarily lying in

L^{2} (R^{2})

) solutions of the equation

A_{α} V = Λ V

. The function

J_{0} (Λ; μ)^{- 1}

is close, in an appropriate sense, to the resolvent of the Jacobi operator

J_{0} (μ)

involved in 1.3 and 1.4 . This allows us to find a connection between the boundary behaviour of these two matrix-valued functions as

Λ

approaches the real line. Technically, this is the most difficult part of the paper. Here we make use of theory of analytic operator-valued functions.

In Section 7 we present Theorem 7.1 which relates the a.c. spectrum of a self-adjoint operator, and also its multiplicity function, to the jump of its ”bordered resolvent” across the real line. We could not find this result in its full generality in the literature. For the reader's convenience, we present its proof in Appendix B .

The above mentioned representation of the operator 1.5 leads to an equality which expresses the jump of the bordered resolvent of the operator

A_{α}

through the similar characteristics of

J_{0} (μ)

. As soon as this is done, Theorem 7.1 applies and gives the equality 1.4 . This scheme is especially transparent for the values

E < 1 / 2

. In order to include the values

E > 1 / 2

, we need an additional technical trick.

Below we briefly describe the structure of the paper. Sections 2 – 4 contain the necessary technical material. Our main result on the absolutely continuous spectrum

σ_{a . c .} (A_{α})

for different values of

α

, Theorem 5.1 , is formulated in Section 5 . Its proof is given in Sections 6 – 8 .We also present Theorem 5.2 on the point spectrum of

A_{α}

. The latter material is mostly borrowed from the papers [13] and [14] . In Section 9 we discuss the possibility to extend the results to operators on quantum graphs. As a matter of fact, the family

A_{α}

was suggested by Smilansky in [12] for this, more general situation.

Two appendices are devoted to the proofs of Theorem 3.1 and Theorem 7.1 . The notation used in the paper is mostly standard. We denote

N_{0} = {0} \cup N

. The symbol

(s) - lim

stands for the strong limit of operators.

Abbreviation ”a.e.” always means almost everywhere with respect to the Lebesgue measure. The symbols

H^{l}

stand for Sobolev spaces. Other necessary notations are introduced in the course of presentation.

2 Reduction to an infinite system of ODE

Equation 1.1 involves the harmonic oscillator in the variable

q

. For this reason, it is convenient to represent the functions

U \in ℌ

U (x, q) = \sum_{n \in N_{0}} u_{n} (x) χ_{n} (q),

where

χ_{n}

are Hermite functions, normalized in

L^{2} (R)

. We often identify a function

U (x, q)

with the sequence

{u_{n} (x)}

and write

U \sim {u_{n}}

. This identification is a unitary mapping of the space

L^{2} (R^{2})

onto the Hilbert space

ℌ = ℓ^{2} (N_{0}, L^{2} (R))

or, equivalently, onto the tensor product

ℓ^{2} \otimes L^{2} (R)

with the natural Hilbert space structure. For

U \sim {u_{n}}

we have

\begin{matrix} A U \sim {(L + n) u_{n}}, (L u) (x) = - u^{''} (x) + u (x) / 2, x \neq 0 . \end{matrix}

(2.1)

The condition 1.2 gives

\sum_{n \in N_{0}} χ_{n} (q) (u_{n}^{'} (0 +) - u_{n}^{'} (0 -)) = α \sum_{n \in N_{0}} q χ_{n} (q) u_{n} (0) .

Taking into account the recurrence equation for the Hermite functions:

\sqrt{n + 1} χ_{n + 1} (q) - \sqrt{2} q χ_{n} (q) + \sqrt{n} χ_{n - 1} (q) = 0,

we come to the system of matching conditions

\begin{matrix} μ (u_{n}^{'} (0 +) - u_{n}^{'} (0 -)) = \sqrt{n + 1} u_{n + 1} (0) + \sqrt{n} u_{n - 1} (0) \end{matrix}

(2.2)

where

μ = \sqrt{2} / α

The operator

A_{0}

admits separation of variables, which leads to the complete description of its spectrum. Let

H

stand for the operator in

L^{2} (R)

, defined as

\begin{matrix} H = - d^{2} / d x^{2} + 1 / 2, D o m H = H^{2} (R) . \end{matrix}

(2.3)

One should distinguish between the self-adjoint operator

H

and the differential expression (formal differential operator)

L

defined in 2.1 .

The operator

A_{0}

splits into the orthogonal sum of the operators

H + n

n \in N_{0}

. An element

U \sim {u_{n}}

belongs to the domain

D_{0} : = D o m A_{0}

if and only if

u_{n} \in H^{2} (R)

for each

n \in N_{0}

and

\sum_{n \in N_{0}} ∥ (H + n) u_{n} ∥^{2} < \infty .

Here and in what follows, unless otherwise explicitly stated, the symbols

(., .), ∥ . ∥

without indication of the space stand for the scalar product and the norm in

L^{2} (R)

The spectrum

σ (H)

is absolutely continuous of multiplicity

2

and coincides with the half-line

[1 / 2, \infty)

. As a consequence, the spectrum of

A_{0}

is also absolutely continuous, and

\begin{matrix} σ (A_{0}) = [1 / 2, \infty), m_{a . c .} (E; A_{0}) = 2 n for E \in (n - 1 / 2, n + 1 / 2), n \in N . \end{matrix}

(2.4)

3 Jacobi matrices $J (Λ; μ)$ and $J_{0} (μ)$

3.1 Preliminaries.

Our main goal in this section is to define the Jacobi matrix

J_{0} (μ)

involved in 1.3 , 1.4 , and another matrix

J (Λ; μ)

, depending on two parameters

Λ \in C

and

μ > 0

. This latter matrix appears when analyzing the homogeneous equation

\begin{matrix} A V = Λ V \end{matrix}

(3.1)

under the matching condition 1.2 . In the representation

V \sim {v_{n}}

equation 3.1 reduces to the infinite system

\begin{matrix} - v_{n}^{''} (x) + (n + 1 / 2 - Λ) v_{n} (x) = 0, x \neq 0; n \in N_{0} \end{matrix}

(3.2)

under the matching conditions 2.2 . Our immediate task is to describe all the formal solutions of this system. This means the following.

Introduce the linear space

\begin{matrix} W (R) = {u \in C (R) : u ↾ R_{\pm} \in H^{2} (R_{\pm})} . \end{matrix}

(3.3)

Note that for any

u \in W (R)

the left-hand side in 2.2 is well-defined.

We seek the solutions

V \sim {v_{n}}

such that

v_{n} \in W (R)

for each

n

, but not necessarily

V \in ℌ

Denote

\begin{matrix} ζ_{n} : = ζ_{n} (Λ) = \sqrt{n + 1 / 2 - Λ} . \end{matrix}

(3.4)

We take the branch of the square root which is analytic in the domain

Ω_{n} : = C \ [n + 1 / 2, \infty)

and such that

ζ_{n} (Λ) > 0

for

Λ = \bar{Λ} < n + 1 / 2

, then

Re ζ_{n} (Λ) > 0, Im ζ_{n} (Λ) \cdot Im Λ < 0, Λ \in Ω_{n} .

The subspace of continuous

L^{2}

-solutions of each equation 3.2 is one-dimensional, it is generated by the function

\begin{matrix} η_{n} (x; Λ) = (n + 1 / 2)^{1 / 4} e^{- ζ_{n} (Λ) | x |} . \end{matrix}

(3.5)

We choose such normalization of the vector-valued functions

η_{n} (.; Λ)

, that each of them is analytic in

Ω_{n}

and

\begin{matrix} c_{1} (Λ) \leq ∥ η_{n} (.; Λ) ∥^{2} \leq c_{2} (Λ), \forall n \in N_{0} . \end{matrix}

(3.6)

We have

\begin{matrix} {\begin{matrix} η_{n} (0; Λ) = (n + 1 / 2)^{1 / 4} & , \\ η_{n}^{'} (0 +; Λ) - η_{n}^{'} (0 -; Λ) = - 2 (n + 1 / 2)^{1 / 4} ζ_{n} (Λ) & . \end{matrix} \end{matrix}

(3.7)

From 3.2 we obtain

v_{n} (x) = C_{n} η_{n} (x; Λ)

, and by 3.7 the matching conditions 2.2 reduce to the recurrence system

(n + 1)^{1 / 2} (n + 3 / 2)^{1 / 4} C_{n + 1} + 2 μ (n + 1 / 2)^{1 / 4} ζ_{n} (Λ) C_{n} + n^{1 / 2} (n - 1 / 2)^{1 / 4} C_{n - 1} = 0, n \in N_{0} .

Taking into account our further needs, we multiply both sides of this equality by

(n + 1 / 2)^{1 / 4}

. As a result, we find that equation 3.1 is equivalent to the system

\begin{matrix} \begin{matrix} (3.8) & V \sim {C_{n} η_{n} (x; Λ)}; d_{n + 1} C_{n + 1} + 2 μ y_{n} (Λ) C_{n} + d_{n} C_{n - 1} = 0 \end{matrix} \end{matrix}

(3.9)

where

\begin{matrix} d_{n} = n^{1 / 2} (n^{2} - 1 / 4)^{1 / 4}, y_{n} (Λ) = (n + 1 / 2)^{1 / 2} ζ_{n} (Λ) . \end{matrix}

(3.10)

The Jacobi matrix

J (Λ; μ)

which corresponds to equation 3.9 is one of our main objects. In our notations we do not distinguish between a Jacobi matrix and the operator which it defines in

ℓ^{2} (N_{0})

. Now we write the operator

J (Λ; μ)

in a more convenient form. Given a number sequence

{ω_{n}}_{n \in N_{0}}

, let

D {ω_{n}}

stand for the diagonal operator in

ℓ^{2} (N_{0})

acting as

D {ω_{n}} : {r_{0}, r_{1}, \dots} \mapsto {ω_{0} r_{0}, ω_{1} r_{1}, \dots} .

Let, in particular,

D = D {d_{n}}, Y (Λ) = D {y_{n} (Λ)} .

Denoting by

S

the operator of the forward shift in

ℓ^{2} (N_{0})

S : {r_{0}, r_{1}, \dots} \mapsto {0, r_{0}, r_{1}, \dots},

we can re-write the operator

J (Λ; μ)

\begin{matrix} J (Λ; μ) = D S + S^{*} D + 2 μ Y (Λ) . \end{matrix}

(3.11)

We also let

\begin{matrix} J_{0} (μ) = D S + S^{*} D + 2 μ Y_{0}, Y_{0} = D {n + 1 / 2} . \end{matrix}

(3.12)

The operator

J_{0} (μ)

, defined initially on the set of all sequences with a finite number of non-zero elements, is essentially self-adjoint in

ℓ^{2}

, and we denote by the same symbol

J_{0} (μ)

its unique self-adjoint extension.

We do not need the explicit description of its domain

D o m J_{0} (μ)

. The next result describes the spectral properties of

J_{0} (μ)

depending on

μ

Theorem 3.1. For

μ > 1

the operator

J_{0} (μ)

is positive definite and its spectrum is discrete.

For

μ \leq 1

the spectrum of

J_{0} (μ)

is purely absolutely continuous and

m_{a . c .} (E; J_{0} (μ)) = 1

a.e. on

σ (J_{0} (μ))

. Moreover,

σ (J_{0} (1)) = [0, \infty); σ (J_{0} (μ)) = R for μ < 1 .

The proof is given in Appendix A. The difference

J (Λ; μ) - (J_{0} (μ) - μ Λ)

is a compact operator. This follows from the equality

\begin{matrix} Ψ (Λ; μ) : = J (Λ; μ) - (J_{0} (μ) - μ Λ) = μ D {ψ_{n} (Λ)} \end{matrix}

(3.13)

where

\begin{matrix} ψ_{n} (Λ) = y_{n} (Λ) - 2 (n + 1 / 2 - Λ / 2) = - Λ^{2} {(4 y_{n} (Λ) + 4 (n + 1 / 2 - Λ / 2))}^{- 1} = O (n^{- 1}) . \end{matrix}

(3.14)

Hence, the operator

J (Λ; μ)

is closed on the domain

D o m J (Λ; μ) = D o m J_{0} (μ)

. Moreover,

Im J (Λ; μ) = 2 μ Im Y (Λ) .

This implies an important property:

Im J (Λ; μ) < 0 for Λ \in C_{+}; Im J (Λ; μ) > 0 for Λ \in C_{-} .

Besides, the operator-valued function

J (Λ; μ)

forms a holomorphic family of type (A) in the variable

Λ

, see [7] , §VII.2. It follows from 3.10 that there is a constant

c > 0

such that

Im y_{n} (- i τ) \geq c \sqrt{τ}, \forall n \in N_{0}, τ > 1 .

Hence, for

τ > 1

we have

Im (J (- i τ; μ) - 2 i μ c \sqrt{τ}) \geq 0 .

By a well known estimate for the dissipative operators, see e.g. [5] , Theorem IV.4.1, this implies a useful inequality

\begin{matrix} ∥ J (- i τ; μ)^{- 1} ∥ \leq (2 μ c \sqrt{τ})^{- 1}, τ > 1 . \end{matrix}

(3.15)

3.2 Birkhoff – Adams theorem.

We base our analysis of the system 3.9 , and also the proof of Theorem 3.1 on the classical result due to Birkhoff and Adams, see e.g. the book [3] , Theorem 8.36. For the reader's convenience, we reproduce the formulation of the part of this theorem which we need below. It concerns the general recurrence system

\begin{matrix} C (n + 1) + p_{1} (n) C (n) + p_{2} (n) C (n - 1) = 0 \end{matrix}

(3.16)

where the functions

p_{1} (n)

and

p_{2} (n)

have asymptotic expansions of the form

\begin{matrix} p_{1} (n) \sim \sum_{j = 0}^{\infty} a_{j} n^{- j}, p_{2} (n) \sim \sum_{j = 0}^{\infty} b_{j} n^{- j}, b_{0} \neq 0 . \end{matrix}

(3.17)

Let

λ_{\pm}

stand for the roots of the equation

λ^{2} + a_{0} λ + b_{0} = 0 .

Proposition 3.2. (a) Let

λ_{+} \neq λ_{-}

, then the system 3.16 has two linearly independent solutions

{C^{\pm}}

with the asymptotics

C^{\pm} (n) \sim λ_{\pm}^{n} n^{d_{\pm}}, d_{\pm} = \frac{a_{1} λ_{\pm} + b_{1}}{a_{0} λ_{\pm} + 2 b_{0}} .

(b) Let

λ_{+} = λ_{-} = λ

but

2 b_{1} \neq a_{0} a_{1}

. Then the system 3.16 has two linearly independent solutions

{C^{\pm} (n)}

with the asymptotics

\begin{matrix} C^{\pm} (n) \sim λ^{n} e^{\pm δ \sqrt{n}} n^{ϰ}, δ = 2 \sqrt{\frac{a_{0} a_{1} - 2 b_{1}}{2 b_{0}}}, ϰ = \frac{1}{4} + \frac{b_{1}}{2 b_{0}} . \end{matrix}

(3.18)

Actually, Theorem 8.36 in [3] describes the complete asymptotic expansions of the solutions

C^{\pm} (n)

of the system 3.9 , but we need only their leading terms. (Note that there is an evident misprint in eq. (8.6.7) in [3] , whose part is reproduced above in 3.18 . In 3.18 this misprint is corrected.) We also need an identity for solutions of recurrence equations with Jacobi matrices, of the type

\begin{matrix} Q_{n + 1} C_{n + 1} + P_{n} C_{n} + Q_{n} C_{n - 1} = 0, n \in N_{0}, \end{matrix}

(3.19)

with

Q_{n}

real and

Q_{0} = 0

. Namely,

\begin{matrix} \sum_{n = 0}^{N} | C_{n} |^{2} Im P_{n} = - Q_{N + 1} Im (C_{N + 1} \bar{C_{N}}), \forall N \in N . \end{matrix}

(3.20)

The proof is straightforward and we skip it.

3.3 Solutions of the system 3.9 .

Here we apply Proposition 3.2 to the system 3.9 that is actually equivalent to equation 3.1 . The system 3.9 can be re-written in the form 3.16 , with the functions

p_{1} (n), p_{2} (n)

admitting the asymptotic expansions of the type 3.17 , where in particular

\begin{matrix} a_{0} = 2 μ, a_{1} = - μ (1 + Λ); b_{0} = 1, b_{1} = - 1 . \end{matrix}

(3.21)

The following Lemma is a direct consequence of Proposition 3.2 .

Lemma 3.3. Let

μ > 0

and

Λ \in Ω_{0}

. Then the system 3.9 has two linearly independent solutions whose asymptotic behaviour is given by

\begin{matrix} C_{n}^{\pm} \sim {\begin{matrix} (μ + i \sqrt{1 - μ^{2}})^{\pm n} n^{- \frac{1}{2} \mp i \frac{Λ μ}{2 \sqrt{1 - μ^{2}}}}, & μ < 1; \\ (- 1)^{n} e^{\pm 2 \sqrt{- Λ n}} n^{- 1 / 4}, & μ = 1; \\ (- μ + \sqrt{μ^{2} - 1})^{\pm n} n^{- \frac{1}{2} \pm \frac{Λ μ}{2 \sqrt{μ^{2} - 1}}}, & μ > 1 . \end{matrix} \end{matrix}

(3.22)

We repeatedly use this Lemma in our further exposition.

4 The self-adjoint operator $A_{α}, α > 0$

Here we describe the domain on which the operator given by equations 1.1 – 1.2 or equivalently, by 2.1 – 2.2 , is self-adjoint in

ℌ

. Consider a linear subset

D_{α} \subset ℌ

: an element

U \sim {u_{n}}

belongs to

D_{α}

if and only if each component

u_{n}

lies in

W (R)

(see 3.3 ), the conditions 2.2 are satisfied, and

\begin{matrix} \sum_{n \in N_{0}} ∥ (L + n) u_{n} ∥^{2} < \infty . \end{matrix}

(4.1)

Along with

D_{α}

, we need its subset

D_{α}^{0}

consisting of all elements

U \in D_{α}

which have only a finite number of non-zero components.

Taking each component equal zero in a vicinity of the point

x = 0

, we obtain a dense subset in

ℌ

. Hence, both

D_{α}^{0}

and

D_{α}

are dense in

ℌ

Define the operators

A_{α}

and

A_{α}^{0}

A_{α} U = A U \sim {(L + n) u_{n}}, D o m A_{α} = D_{α}; A_{α}^{0} = A_{α} ↾ D_{α}^{0} .

Evidently, the operator

A_{α}^{0}

is symmetric in

ℌ

. Our goal is to prove the following result.

Theorem 4.1. The operator

A_{α}

is self-adjoint and coincides with the closure of

A_{α}^{0}

Note that for

α \neq \sqrt{2}

the result is covered by [13] , Theorem 5.1. Nevertheless, below we give the full proof of Theorem. We do this mostly in order to illustrate the usage of Proposition 3.2 . In [13] another, more sophisticated technical tools were used for the proof.

Proof. First of all, we show that

\begin{matrix} A_{α} = {(A_{α}^{0})}^{*} . \end{matrix}

(4.2)

The inclusion “

\subset

” in 4.2 can be easily checked by the direct inspection.

To prove the reverse inclusion, suppose that

V \sim {v_{n}} \in D o m {(A_{α}^{0})}^{*}

and

{(A_{α}^{0})}^{*} V = W \sim {w_{n}}

. According to the definition of the adjoint operator, this means that for any

U \sim {u_{n}} \in D_{α}^{0}

we have

\begin{matrix} \sum_{n \in N_{0}} \int_{R} (L + n) u_{n} \bar{v_{n}} d x = \sum_{n \in N_{0}} \int_{R} u_{n} \bar{w_{n}} d x . \end{matrix}

(4.3)

Take a function

f \in W (R)

such that

f (x) = 0

in a vicinity of

x = 0

and fix a number

n_{0} \in N_{0}

. The element

U \sim {u_{n}}

, such that

u_{n_{0}} = f

and

u_{n} = 0

for

n \neq n_{0}

, belongs to

D_{α}^{0}

. Applying the identity 4.3 to all such

U

, we conclude that if

V \in D o m {(A_{α}^{0})}^{*}

, then

v_{n} \in W (R)

and

w_{n} = (L + n) v_{n}

for all

n \in N_{0}

. Hence, for

{v_{n}}

the condition 4.1 is satisfied. It remains to check that the matching conditions 2.2 are also fulfilled. To this end, we fix a number

n_{0} \in N_{0}

and choose an element

U \sim {u_{n}}

as follows:

u_{n} (x) \equiv 0

for

| n - n_{0} | > 1

; the functions

u_{n_{0}}, u_{n_{0} \pm 1}

are supported in a vicinity of the

x = 0

and in some smaller vicinity are given by

u_{n_{0}} (x) = 1

;

u_{n_{0} \pm 1} (x) = h_{\pm} | x |

where

h_{\pm}

are some appropriate numbers. Then

U \in D_{α}^{0}

, provided that

2 μ h_{+} = \sqrt{n_{0} + 1}, 2 μ h_{-} = \sqrt{n_{0}} .

Now 2.2 for

n = n_{0}

is implied by 4.3 for the element

U

constructed.

So, the equality 4.2 is justified.

The statement of Theorem is equivalent to the fact that both deficiency indexes of the operator

A_{α}^{0}

are equal to zero. Since all the coefficients in the equation and in the matching conditions are real, it is enough to prove that the only solution

V \in ℌ

of the equation

\begin{matrix} A_{α} V - i V = 0 \end{matrix}

(4.4)

V \equiv 0

. Equation 4.4 is nothing but 3.1 for

Λ = i

. Using the representation 3.8 , we reduce the equation to the form 3.9 and can apply Lemma 3.3 .

μ < 1

, then according to 3.22 the system has a pair of linearly independent solutions

{C_{n}^{\pm}}

such that

\begin{matrix} | C_{n}^{\pm} |^{2} \sim n^{- 1 \pm \frac{μ}{\sqrt{1 - μ^{2}}}} . \end{matrix}

(4.5)

In view of 3.6 , only the sequence

{C_{n}^{-}}

may generate a solution

V \in ℌ

of equation 4.4 .

Equation 3.9 is of the form 3.19 . Now we use for it the identity 3.20 , with

C_{n} = C_{n}^{-}

. This gives

\begin{matrix} \sum_{n = 0}^{N} | C_{n}^{-} |^{2} Im ζ_{n} = - \sqrt{N + 1} Im (C_{N + 1}^{-} \bar{C_{N}^{-}}) . \end{matrix}

(4.6)

By 4.5 , the right-hand side of 4.6 vanishes as

N \to \infty

. Since

Im ζ_{n} > 0

for each

n

, we conclude from 4.6 that

V \equiv 0

. This shows that for

α > \sqrt{2}

the operator

A_{α}

is self-adjoint.

α < \sqrt{2}

, then

μ > 1

and according to 3.22 one of the solutions of the system 3.9 exponentially grows and another exponentially decays. Only the latter may give rise to the solution

V \in ℌ

of equation 4.4 . Again, using the identity 3.20 we conclude that this solution is identically zero. Finally, let

α = \sqrt{2}

, then

μ = 1

and the formula 3.22 gives

C_{n}^{\pm} \sim e^{\pm \sqrt{2 n} (1 + i)} n^{- 1 / 4} .

Only the sequence

{C_{n}^{-}}

lies in

ℓ^{2}

and again, we conclude from 3.20 that the deficiency indexes are equal to

0

. □

5 Spectrum of the operators $A_{α}$

5.1 Results.

The following theorem is the central result of the paper.

Theorem 5.1. Let

α > 0

μ = \sqrt{2} / α

and let

J_{0} (μ)

be the Jacobi matrix (operator), defined in 3.12 . Then the a.c. spectrum of the operator

A_{α}

and its multiplicity are described by the equalities 1.3 and 1.4 . In particular, for

α < \sqrt{2}

\begin{matrix} {\begin{matrix} σ_{a . c .} (A_{α}) = [1 / 2, \infty) & , \\ m_{a . c .} (E; A_{α}) = m_{a . c .} (E; A_{0}), & a . e . E \geq 1 / 2 . \end{matrix} \end{matrix}

(5.1)

Further,

\begin{matrix} σ_{a . c .} (A_{\sqrt{2}}) = [0, \infty); σ_{a . c .} (A_{α}) = R, α > \sqrt{2}, \end{matrix}

(5.2)

and for

α \geq \sqrt{2}

and

E \in σ_{a . c .} (A_{α})

we have

\begin{matrix} m_{a . c .} (E; A_{α}) = m_{a . c .} (E; A_{0}) + 1 . \end{matrix}

(5.3)

The equalities 5.1 , 5.2 and 5.3 immediately follow from the relations 1.3 , 1.4 and Theorem 3.1 . So, our main goal for the rest of the paper is to prove 1.3 and 1.4 .

To make the picture more complete, we present also the result concerning the point spectrum

σ_{p} (A_{α})

. Within minor detail, these results were proved in the papers [13] , [14] .

Theorem 5.2. 1. For any

α \geq 0

the operator

A_{α}

has no eigenvalues

E \geq 1 / 2

, and

σ_{p} (A_{α}) = \emptyset

for

α \geq \sqrt{2}

2. For

0 < α < \sqrt{2}

the operator

A_{α}

is positive definite.

Its point spectrum is always non-empty and finite, and the number

N (1 / 2; A_{α})

of its eigenvalues (counted with multiplicities) satisfies the asymptotic formula

N (1 / 2; A_{α}) \sim \frac{1}{4 \sqrt{2 (μ (α) - 1)}}, α ↗ \sqrt{2} (μ (α) = \frac{\sqrt{2}}{\sqrt{α}}) .

It follows from Theorems 5.1 and 5.2 that the essential and the absolutely continuous spectra of the operator

A_{α}

coincide as sets.

However, our approach does not show that

A_{α}

has no singular continuous spectrum.

5.2 Outline of proof of Theorem 5.2 .

1. If

E \geq 1 / 2, V \sim {v_{n}}

and

A_{α} V = E V

, then each component

v_{n}

lies in

W (R)

and satisfies equation 3.2 with

Λ = E

. For

n \leq E - 1 / 2

this yields

v_{n} \equiv 0

. This can be interpreted as the equality

C_{n} = 0

for the coefficients in 3.9 . Then these equations imply that

C_{n} = 0

also for all

n > E - 1 / 2

The absence of eigenvalues

E < 1 / 2

for the operator

A_{α}

with

α > \sqrt{2}

was proved in [13] , Theorem 7.1. The possibility to use Proposition 3.2 simplifies the proof, and also allows one to include the borderline case

α \geq \sqrt{2}

. We leave this to the reader.

The statement 2 is covered by [13] , Theorem 6.2 and [14] , formula (3.10).

6 Representation of the resolvent

6.1 Auxiliary considerations

In this section we derive a convenient representation of the operator

(A_{α} - Λ)^{- 1} - (A_{0} - Λ)^{- 1}, Λ / \in R .

The equality 6.6 which we establish in Theorem 6.1 can be interpreted in terms of the extension theory of symmetric operators. However, formally we do not use this theory in our construction.

Given an element

F \sim {f_{n}} \in ℌ

, denote

U_{α} \sim {u_{α, n}} = (A_{α} - Λ)^{- 1} F, α \geq 0 .

The functions

u_{0, n}

satisfy the equation

- u_{0, n}^{''} + (n + 1 / 2 - Λ) u_{0, n} = f_{n}

and lie in

H^{2} (R)

. Solving this equation, we find that

2 ζ_{n} (Λ) u_{0, n} (x) = \int_{R} e^{- ζ_{n} (Λ) | x - t |} f_{n} (t) d t .

Denote

\begin{matrix} J_{n} = \int_{R} η_{n} (t; Λ) f_{n} (t) d t = (f_{n}, η_{n} (.; \bar{Λ})), \end{matrix}

(6.1)

then

\begin{matrix} 2 ζ_{n} (Λ) (n + 1 / 2)^{1 / 4} u_{0, n} (0) = J_{n}, n \in N_{0} . \end{matrix}

(6.2)

Recall that the numbers

ζ_{n} (Λ)

and the functions

η_{n} (.; Λ)

were defined in 3.4 and 3.5 .

Let now

U_{α} - U_{0} \sim {v_{n}}

. Each function

v_{n}

satisfies the homogeneous equation 3.2 and belongs to

W (R)

and hence,

v_{n} (x) = C_{n} η_{n} (x; Λ)

. The coefficients

C_{n}

are determined by the matching conditions for

U_{α}

. Since in view of 6.2

\begin{matrix} (n + 1 / 2)^{1 / 4} C_{n} = v_{n} (0) = u_{α, n} (0) - u_{0, n} (0) = u_{α, n} (0) - (2 ζ_{n} (Λ) (n + 1 / 2)^{1 / 4})^{- 1} J_{n} \end{matrix}

(6.3)

and the derivative

u_{0, n}^{'}

is continuous at

x = 0

, we get

u_{α, n}^{'} (0 +) - u_{α, n}^{'} (0 -) = v_{n}^{'} (0 +) - v_{n}^{'} (0 -) = - 2 ζ_{n} (Λ) (n + 1 / 2)^{1 / 4} C_{n} = (n + 1 / 2)^{- 1 / 4} J_{n} - 2 ζ_{n} (Λ) u_{α, n} (0) .

It is convenient for us to denote

X_{n} = (n + 1 / 2)^{- 1 / 4} u_{α, n} (0),

then

u_{α, n}^{'} (0 +) - u_{α, n}^{'} (0 -) = (n + 1 / 2)^{- 1 / 4} J_{n} - 2 ζ_{n} (Λ) (n + 1 / 2)^{1 / 4} X_{n}

and the matching conditions 2.2 reduce to

μ ((n + 1 / 2)^{- 1 / 4} J_{n} - 2 ζ_{n} (Λ) (n + 1 / 2)^{1 / 4} X_{n}) = (n + 1)^{1 / 2} (n + 3 / 2)^{1 / 4} X_{n + 1} + n^{1 / 2} (n - 1 / 2)^{1 / 4} X_{n - 1},

d_{n + 1} X_{n + 1} + 2 μ y_{n} (Λ) X_{n} + d_{n} X_{n - 1} = μ (f_{n}, η_{n} (.; \bar{Λ})) .

This is the non-homogeneous counterpart of the recurrence system 3.9 .

It can be re-written in terms of the matrix

J (Λ; μ)

introduced in 3.11 :

\begin{matrix} J (Λ; μ) X = μ (f_{n}, η_{n} (.; \bar{Λ})), X = {X_{n}} . \end{matrix}

(6.4)

6.2 Basic formula

For

Λ / \in R

consider the operator

\begin{matrix} T (Λ) : ℓ^{2} (N_{0}) \to ℌ, T (Λ) {X_{n}} \sim {X_{n} η_{n} (.; Λ)}, \end{matrix}

(6.5)

by 3.6 it is bounded and has bounded inverse. Its adjoint acts from

ℌ

ℓ^{2}

T (Λ)^{*} F = {\int_{R} f_{n} (x) η_{n} (x; \bar{Λ}) d x}, F \sim {f_{n}} .

Substituting

\bar{Λ}

for

Λ

, we obtain that

T (\bar{Λ})^{*} F = {(f_{n}, η_{n} (.; Λ))} .

Both operator-valued functions

T (Λ)

and

T (\bar{Λ})^{*}

are analytic in

C_{\pm}

Now we are in a position to present the basic formula which relates the operator

A_{α}

to the Jacobi matrix

J (Λ; μ), μ = \sqrt{2} / α

Theorem 6.1. Let

α > 0, μ = \sqrt{2} / α

and

Λ / \in R

. Then

\begin{matrix} (A_{α} - Λ)^{- 1} - (A_{0} - Λ)^{- 1} = T (Λ) (μ J (Λ; μ)^{- 1} - (2 Y (Λ))^{- 1}) T (\bar{Λ})^{*} . \end{matrix}

(6.6)

Proof. By 6.1 and 6.3 , we have $C_{n} = X_{n} - (2 y_{n} (Λ))^{- 1} (f_{n}, η_{n} (.; \bar{Λ}))$ , or ${C_{n}} = X - (2 Y (Λ))^{- 1} T (\bar{Λ})^{*},$ whence $U_{α} - U_{0} \sim {v_{n}} = T (Λ) {C_{n}} .$ We also find from 6.4 that $X = μ J (Λ; μ)^{- 1} T (\bar{Λ})^{*}$ . The desired equality 6.6 is an immediate consequence of the three last equations. □

7 Proof of Theorem 5.1 : $E < 1 / 2$

7.1 Absolutely continuous spectrum and jump of the bordered resolvent.

Let

Ω \subset C

be a domain, symmetric with respect to the real axis and containing an interval

I \subset R

. Let

Z (Λ)

be an operator-valued function which is analytic in

Ω \cap C_{+}

and in

Ω \cap C_{-}

. Its jump at a point

E \in I

is defined as

[Z] (E) : = {[Z (Λ)]}_{Λ \to E + i 0} = (s) - {lim}_{Λ \to E + i 0} (Z (Λ) - Z (\bar{Λ})) (non-tangential limit),

provided that the (strong) limit does exist.

One of the ways to investigate the absolutely continuous spectrum of a self-adjoint operator

A

consists in studying the jump of its ”bordered resolvent”

Z_{G} (Λ) : = G^{*} (A - Λ)^{- 1} G

. Here

G

is an appropriate bounded operator. It is well known that for any

A

and any

G \in S_{2}

the non-tangential limits of the bordered resolvent

Z_{G} (Λ)

Λ \to E \pm i 0

exist a.e. and belong to

S_{1}

, see e.g. [1] . Hence, the jump also exists a.e. (with respect to the Lebesgue measure). For a given operator

A

the jump of its bordered resolvent may exist also for a wider class of borderings.

The following statement of a rather general nature plays the key role in our analysis. Its proof is given in Appendix B.

Theorem 7.1. Let

A

be a self-adjoint operator in a separable Hilbert space

ℋ

and

Δ \subset R

a given Borelian subset.

Let

G

be a bounded operator with the dense range, such that the operator

\begin{matrix} V_{A, G} (E) = (s) - {lim}_{Λ \to E + i 0} G^{*} ((A - Λ)^{- 1} - (A - \bar{Λ})^{- 1}) G \end{matrix}

(7.1)

is well-defined a.e. on

Δ

. Then

\begin{matrix} m_{a . c .} (E; A) = rank V_{A, G} (E) a.e. on Δ . \end{matrix}

(7.2)

Corollary 7.2. Let the assumptions of Theorem 7.1 be fulfilled for an interval

Δ \subseteq R

. If

rank V_{A, G} (E) > 0

a.e. on

Δ

, then

Δ \subset σ_{a . c .} (A)

, and if

rank V_{A, G} (E) = 0

a.e. on

Δ

, then

σ_{a . c .} (A) \cap Δ = \emptyset

7.2 Bordered resolvent of the operator $A_{α}$

Along with the operator

Y_{0} = D {n + 1 / 2}

acting in

ℓ^{2}

, let us define

Y_{0} = I \otimes Y_{0}

. The latter operator acts in the space

ℌ

interpreted as the tensor product

ℓ^{2} \otimes L^{2} (R)

, cf. Section 2 . The powers

Y_{0}^{- γ}, γ > 0

are compact operators in

ℓ^{2}

, while the powers

Y_{0}^{- γ}

are only bounded operators. If

U \sim {u_{n}} \in ℌ

, then

Y_{0}^{- γ} U \sim {(n + 1 / 2)^{- γ} u_{n}}

. Note that

\begin{matrix} Y^{- γ} T (Λ) = T (Λ) Y^{- γ}, T (\bar{Λ})^{*} Y^{- γ} = Y^{- γ} T (\bar{Λ})^{*} . \end{matrix}

(7.3)

We shall show that the boundary limits of the bordered resolvent

G^{*} (A_{α} - Λ)^{- 1} G

do exist if we take

G = Y_{0}^{- γ}

with

γ > 1 / 4

, though

G

is non-compact, and even the operator

Y_{0}^{- γ}

lies in

S_{2}

only if

γ > 1 / 2

Let us consider the operator-valued functions

\begin{matrix} Z_{γ} (Λ; μ) = Y_{0}^{- γ} (A_{α} - Λ)^{- 1} Y_{0}^{- γ}, α = \sqrt{2} / μ, \tilde{Z_{γ}} (Λ; μ) = Y_{0}^{- γ} J (Λ; μ)^{- 1} Y_{0}^{- γ}, \tilde{Z_{0, γ}} (Λ; μ) = Y_{0}^{- γ} (J_{0} (μ) - μ Λ)^{- 1} Y_{0}^{- γ} . \end{matrix}

(7.4)

In 7.4

J_{0} (μ)

is the Jacobi matrix introduced in 3.12 .

By 6.6 and 7.3 ,

Z_{γ} (Λ; μ) = Y_{0}^{- γ} (A_{0} - Λ)^{- 1} Y_{0}^{- γ} + T (Λ) (μ \tilde{Z_{γ}} (Λ; μ) - Y_{0}^{- γ} (2 Y (Λ))^{- 1} Y_{0}^{- γ}) T (\bar{Λ})^{*} .

The jumps at a point

E \in R

of all but one terms in the right-hand side are evidently equal to zero, provided that

E < 1 / 2

, and we get, at least formally,

\begin{matrix} {[Z_{γ} (Λ; μ)]}_{Λ \to E + i 0} = μ T (E) {[\tilde{Z_{γ}} (Λ; μ)]}_{Λ \to E + i 0} T (E)^{*}, E < 1 / 2 . \end{matrix}

(7.5)

Taking into account the equalities 3.13 and 3.14 , it is natural to expect that the jumps

[\tilde{Z_{γ}} (.; μ)] (E)

and

[\tilde{Z_{0, γ}} (.; μ)] (E)

are close to each other. Together with 7.5 , this would allow us to express the quantity

[Z_{γ} (.; μ)] (E)

through

[\tilde{Z_{0, γ}} (.; μ)] (E)

, which makes it possible to use Theorem 7.1 . Now we proceed to the successive realization of this program.

We start with the following lemma.

Lemma 7.3. Let

μ > 0

Λ \neq \bar{Λ} \in C

and

γ > 1 / 4

. Then

\tilde{Z_{γ}} (Λ; μ) \in S_{1}

. Besides, the strong non-tangential limits

\tilde{Z_{γ}} (E \pm i 0; μ)

exist for almost all

E < 1 / 2

The same results are valid for

\tilde{Z_{0, γ}} (Λ; μ)

Proof. For definiteness, we consider

Λ \in C_{+}

. The operator-valued function

\tilde{Z_{γ}} (Λ; μ)

is analytic in

Λ

in the half-plane

C_{+}

. Since

Im J (Λ; μ) \leq 0

, we conclude that

Im \tilde{Z_{γ}} (Λ; μ) \geq 0

, i.e. the values of this function are dissipative operators. By a result due to S.Naboko, see [10] , Remark (1) to Theorem 2.2, any such function admits the representation

\begin{matrix} \tilde{Z_{γ}} (Λ; μ) = A + B Λ + R^{*} (I + Λ ℒ) (ℒ - Λ)^{- 1} R, \end{matrix}

(7.6)

where

ℒ

is a self-adjoint operator in an auxiliary Hilbert space

ℌ_{0}

A = A^{*}

and

B = B^{*} \geq 0

are bounded operators in

ℓ^{2}

, and

R

is a bounded operator from

ℓ^{2}

ℌ_{0}

. It immediately follows from 3.15 that in this representation

B = 0

. The operators

ℒ, A

and

R

may depend on

μ

but we do not reflect this dependence in our notations.

Taking in 7.6

Λ = i

, we find that

\begin{matrix} \tilde{Z_{γ}} (i; μ) = A + i R^{*} R . \end{matrix}

(7.7)

Further, let us show that for any

γ > 0

one has

\tilde{Z_{γ}} (i; μ) \in S_{2}

. To this end, consider the inverse matrix

{G_{n, k}} : = J (i; μ)^{- 1}

. This matrix can be expressed through the solutions of the homogeneous system 3.9 for

Λ = i

. Namely, it is symmetric (i.e.

G_{n, k} = G_{k, n}

) and its entries for

k \geq n

are

G_{n, k} = C g^{(1)} (k) g^{(2)} (n)

where

C

is a constant,

{g^{(1)} (k)}

is the decaying solution of the system, and

{g^{(2)} (n)}

is some other solution, a certain linear combination of two basic solutions. The matrix

\tilde{Z_{γ}} (i; μ)

has the entries

(n + 1 / 2)^{- γ} G_{n, k} (k + 1 / 2)^{- γ},

and the standard calculation shows that

\tilde{Z_{γ}} (i; μ) \in S_{2}

for any

γ > 0

and an arbitrary

μ > 0

. This immediately yields that

\tilde{Z_{γ}} (i; μ) \in S_{1}

, provided that

γ > 1 / 4

Now we conclude from 7.7 that

R \in S_{2}

. Hence,

\tilde{Z_{γ}} (Λ; μ) \in S_{1}

for all

Λ \in C_{+}

. Besides, by a result of [1] , Lemma 2.4, this implies the existence a.e. of the non-tangential strong limits of the bordered resolvent

R^{*} (ℒ - Λ)^{- 1} R

. The equality 7.6 shows that such limits do exist also for the operator-valued function

\tilde{Z_{γ}} (Λ; μ)

For the operator-valued function

\tilde{Z_{0, γ}}

the proof is simpler, since we need not the representation 7.6 . Otherwise, the argument remains the same. □

As a result of Lemma 7.3 , the equality 7.5 is justified.

7.3 Jumps of $\tilde{Z_{γ}}$ and of $\tilde{Z_{0, γ}}$ .

A subset

Δ \subset (- \infty, 1 / 2)

of full measure can be selected, such that both these jumps are well-defined for

E \in Δ

. According to the Hilbert identity and using the equality 3.13 , we find that

\begin{matrix} \tilde{Z_{γ}} (Λ; μ) - \tilde{Z_{0, γ}} (Λ; μ) = Y_{0}^{- γ} (J (Λ; μ)^{- 1} - (J_{0} (μ) - μ Λ)^{- 1}) Y_{0}^{- γ} = - Y_{0}^{- γ} J (Λ; μ)^{- 1} Ψ (Λ; μ) (J_{0} (μ) - μ Λ)^{- 1} Y_{0}^{- γ} = - \tilde{Z_{γ}} (Λ; μ) Φ (Λ; μ) \tilde{Z_{0, γ}} (Λ; μ) \end{matrix}

(7.8)

where

\begin{matrix} Φ (Λ; μ) = Y_{0}^{γ} Ψ (Λ; μ) Y_{0}^{γ} . \end{matrix}

(7.9)

More exactly,

Φ (Λ; μ)

is the extension by continuity of this operator, defined originally on the set

R a n Y_{0}^{γ}

The calculations below are carried through for a fixed value of

μ

, and we drop this parameter from our notations. It follows from 3.14 that for

γ < 1 / 2

the operator-valued function

Φ (\cdot)

is analytic in

C_{\pm}

and its values are compact operators in

ℓ^{2}

. Below we always assume

1 / 4 < γ < 1 / 2

, then both operator-valued functions

\tilde{Z_{γ}} (Λ)

and

\tilde{Z_{0, γ}} (Λ)

have boundary limits as

Λ \to E \pm i 0

a.e. on

(- \infty, 1 / 2)

The equality 7.8 and the similar equality for

\bar{Λ}

imply that

\tilde{Z_{0, γ}} (Λ) - \tilde{Z_{0, γ}} (\bar{Λ}) = \tilde{Z_{γ}} (Λ) - \tilde{Z_{γ}} (\bar{Λ}) + (\tilde{Z_{γ}} (Λ) - \tilde{Z_{γ}} (\bar{Λ})) Φ (Λ) \tilde{Z_{0, γ}} (Λ) + \tilde{Z_{γ}} (\bar{Λ}) (Φ (Λ) - Φ (\bar{Λ})) \tilde{Z_{0, γ}} (Λ) + \tilde{Z_{γ}} (\bar{Λ}) Φ (\bar{Λ}) (\tilde{Z_{0, γ}} (Λ) - \tilde{Z_{0, γ}} (\bar{Λ})) .

Letting

Λ \to E + i 0

non-tangentially, we obtain for a.e.

E < 1 / 2

\begin{matrix} [\tilde{Z_{0, γ}}] (E) = [\tilde{Z_{γ}}] (E) + [\tilde{Z_{γ}}] (E) Φ (E) \tilde{Z_{0, γ}} (E + i 0) + \tilde{Z_{γ}} (E - i 0) Φ (E) [\tilde{Z_{0, γ}}] (E) . \end{matrix}

(7.10)

Let us consider the operator-valued functions

\begin{matrix} \begin{matrix} (7.11) & G_{+} (Λ) = I + Φ (Λ) \tilde{Z_{0, γ}} (Λ), Λ \in C_{+}, G_{-} (Λ) = I - \tilde{Z_{γ}} (Λ) Φ (Λ), Λ \in C_{-} . \end{matrix} \end{matrix}

(7.12)

The function

G_{+} (Λ)

can be represented in a different way: using the definitions 7.9 and 7.4 , we find that

\begin{matrix} G_{+} (Λ) = I + Y_{0}^{γ} (J (Λ) - (J_{0} - μ Λ)) Y_{0}^{γ} Y_{0}^{- γ} (J_{0} - μ Λ)^{- 1} Y_{0}^{- γ} = Y_{0}^{γ} J (Λ) (J_{0} - μ Λ)^{- 1} Y_{0}^{- γ} . \end{matrix}

(7.13)

This calculation shows also that the right-hand side in 7.13 is well-defined as a bounded operator in

ℓ^{2}

. For the ”minus” sign, it is more convenient to deal with the adjoint operator, and we get

\begin{matrix} G_{-} (Λ)^{*} = Y_{0}^{γ} (J_{0} - μ \bar{Λ}) J (\bar{Λ})^{- 1} Y_{0}^{- γ} . \end{matrix}

(7.14)

It follows from the definitions 7.11 and 7.12 and from compactness of

Φ (Λ)

that the functions

G_{\pm} (Λ)

are analytic in

C_{\pm}

respectively, and for each

Λ \in C_{\pm}

the operator

G_{\pm} (Λ) - I

is compact. Hence, the image of

G_{\pm} (Λ)

is a closed subset in

ℓ^{2}

Now, we need the following Lemma.

Lemma 7.4. For any

Λ \in C_{\pm}

the operator

G_{\pm} (Λ)

has bounded inverse.

Proof. Throughout the proof,

(., .)

stands for the scalar product in

ℓ^{2}

By 7.13 , we have for any

f \in ℓ^{2}

and any

g \in D o m Y_{0}^{γ}

\begin{matrix} (G_{+} (Λ) f, g) = (J (Λ) (J_{0} - μ Λ)^{- 1} Y_{0}^{- γ} f, Y_{0}^{γ} g) . \end{matrix}

(7.15)

Suppose that

f \in K e r G_{+} (Λ)

for some

f \neq 0

and

Λ \in C_{+}

. Since

R a n Y_{0}^{γ} = ℓ^{2}

, we conclude from 7.15 that then

K e r J (Λ)

is non-trivial which contradicts the dissipativity of

J (Λ)

In a similar way, we derive from 7.14 that for

f \in D o m Y_{0}^{γ}

and any

g \in ℓ^{2}

(G_{-} (Λ) f, g) = (Y_{0}^{γ} f, (J_{0} - μ \bar{Λ}) J (\bar{Λ})^{- 1} Y_{0}^{- γ} g) .

Suppose that

R a n G_{-} (Λ) \neq ℓ^{2}

for some

Λ \in C_{-}

. Then there exists an element

g \in ℓ^{2}

, such that

(G_{-} (Λ) f, g) = 0

for all

f \in ℓ^{2}

. This would imply

(J_{0} - μ \bar{Λ}) J (\bar{Λ})^{- 1} Y_{0}^{- γ} g = 0 .

However, for

g \neq 0

this is impossible, since

J (\bar{Λ})^{- 1}

and

Y_{0}^{- γ}

are inverse operators and hence, have trivial kernels. Further,

K e r (J_{0} - μ \bar{Λ}) = {0}

for

Λ / \in R

, since the operator

J_{0}

is self-adjoint. □

Now we use the following statement, see [9] .

Proposition 7.5. Let

X \in C_{\pm}

be a domain, such that

\bar{X} \cap R

contains an interval

Δ

. Let

G (Λ)

be an analytic operator-valued function in

X

, such that

G (Λ) - I

are compact operators in a Hilbert space

ℋ

. Suppose that for almost all

E \in Δ

the function

G

is non-tangentially bounded at

E

and has a strong non-tangential limit

G (E + i 0)

, such that

G (E + i 0) - I \in S_{\infty}

Suppose also that for at least one point

Λ_{0} \in C_{+}

the operator

G (Λ_{0})

has bounded inverse.

Then for almost all

E \in Δ

the operator

G (E + i 0)

has bounded inverse.

It follows from Lemmas 7.3 and 7.4 that the assumptions of Proposition are fulfilled for the operator-valued functions

G_{\mp} (Λ)

defined in 7.11 and 7.12 . This allows us to conclude from 7.5 and 7.10 that for a.e.

E < 1 / 2

\begin{matrix} [Z_{γ}] (E) = T (E) G_{-} (E - i 0) [\tilde{Z_{0, γ}}] (E) {(G_{+} (E + i 0))}^{- 1} T (E)^{*} . \end{matrix}

(7.16)

The equalities 7.5 and 7.16 yield that

rank [Z_{γ}] (E) = rank [\tilde{Z_{0, γ}}] (E), a . e . E < 1 / 2 .

According to Theorem 7.1 , it follows that

m_{a . c .} (E; A_{α}) = m_{a . c .} (E; J_{0} (μ)) a . e . on (- \infty, 1 / 2) .

Since

m_{a . c .} (E; A_{0}) = 0

for

E < 1 / 2

, the equality 1.4 for such

E

is justified.

Now, we conclude from Corollary 7.2 that for

α < \sqrt{2}

the operator

A_{α}

has no a.c. spectrum below the point

1 / 2

. For

α > \sqrt{2}

we have

(- \infty, 1 / 2) \subset σ_{a . c .} (A_{α})

, and

[0, 1 / 2) \subset σ_{a . c .} (A_{\sqrt{2}})

8 Proof of Theorem 5.1 : $E \geq 1 / 2$

The proof extends to

λ > 1 / 2

with the help of a simple technical trick. It is based on the passage to the subspace

ℌ_{m} = ℓ^{2} (N_{m}, L^{2} (R))

, with

m

large enough. Here

N_{m} = {m, m + 1, \dots}

, so that, in particular,

N_{1} = N

. The subspace

ℌ_{m}

is not invariant for the operator

A_{α}

, however it is invariant for an appropriate operator

{\hat{A}}_{α}^{(m)}

, that can be obtained from

A_{α}

by the perturbation of its resolvent by a finite rank operator.

For the operator

{\hat{A}}_{α}^{(m)}

the scheme developed in Section 7 works for

E < m + 1 / 2

, and return to the original operator

A_{α}

does not change the absolutely continuous spectrum and its multiplicity function.

8.1 Operators $A_{α}^{(m)}$ and ${\hat{A}}_{α}^{(m)}$ .

Denote

N_{m} = {m, m + 1, \dots}

, so that, in particular,

N_{1} = N

. Let us define an operator

A_{α}^{(m)}

, acting in the space

ℌ_{m} = ℓ^{2} (N_{m}, L^{2} (R))

Namely, its domain

D_{α}^{(m)}

consists of the elements

U \sim {u_{n}}_{n \geq m}

, such that each component

u_{n}

lies in

W (R)

(see 3.3 ), the matching conditions

u_{n}^{'} (0 +) - u_{n}^{'} (0 -) = {\begin{matrix} \frac{α}{\sqrt{2}} \sqrt{m + 1} u_{m + 1} (0), & n = m, \\ \frac{α}{\sqrt{2}} (\sqrt{n + 1} u_{n + 1} (0) + \sqrt{n} u_{n - 1} (0)), & n > m \end{matrix}

are satisfied, and

\sum_{n \geq m} ∥ (L + n) u_{n} ∥^{2} < \infty

. For

U \in D_{α}^{(m)}

we let

A_{α}^{(m)} U \sim {(L + n) u_{n}}_{n \geq m},

cf. 2.1 . Evidently,

A_{α}^{(0)} = A_{α}

. Each operator

A_{α}^{(m)}

is self-adjoint, the proof is the same as for

m = 0

Along with

A_{α}^{(m)}

, we need also the operators

{\hat{A}}_{α}^{(m)}

, acting in the original Hilbert space

ℌ

. Namely, we let

\begin{matrix} {\hat{A}}_{α}^{(m)} = ({\sum_{n < m}}^{\oplus} (H + n)) \oplus A_{α}^{(m)} \end{matrix}

(8.1)

where

H

is the operator in

L^{2} (R)

, defined in 2.3 . Note that

{\hat{A}}_{0}^{(m)} = A_{0}

for any

m

Consider also two sequences of Jacobi matrices,

J_{0}^{(m)} (μ)

and

J^{(m)} (Λ; μ)

Each

J_{0}^{(m)} (μ)

is a sub-matrix of the matrix

J_{0} (μ)

given by 3.12 , it is obtained from

J_{0} (μ)

by removing its first

m

rows and

m

columns. The matrix

J^{(m)} (Λ; μ)

is obtained in the same way from the matrix 3.11 .

We also define a sequence

Y^{(m)} (Λ)

of diagonal operators in

ℓ^{2} (N_{m})

, namely

Y^{(m)} (Λ) = D {y_{n} (Λ)}_{n \geq m}

. It is clear that

J_{0}^{(0)} (μ) = J_{0} (μ), J^{(0)} (Λ; μ) = J (Λ; μ), Y^{(0)} (Λ) = Y (Λ) .

Let

T^{(m)} (Λ) : ℓ^{2} (N_{m}) \to ℌ_{m}

be the natural restriction of the operator

T (Λ)

defined in 6.5 . The equality

\begin{matrix} (A_{α}^{(m)} - Λ)^{- 1} - (A_{0}^{(m)} - Λ)^{- 1} = T^{(m)} (Λ) (μ J^{(m)} (Λ; μ)^{- 1} - (2 Y^{(m)} (Λ))^{- 1}) T^{(m)} (\bar{Λ})^{*} \end{matrix}

(8.2)

can be justified in the same way as 6.6 . It is important that all the operator-valued functions appearing in 8.2 are analytic in the domain

Ω_{m} = C \ [m + 1 / 2, \infty)

. It follows from 8.2 that the operator

({\hat{A}}_{α}^{(m)} - Λ)^{- 1} - (A_{0} - Λ)^{- 1}

is equal to the orthogonal sum of the null operator acting in the space

ℓ^{2} ({0, \dots, m - 1}, L^{2} (R))

and the operator

T^{(m)} (Λ) (μ J^{(m)} (Λ; μ)^{- 1} - (2 Y^{(m)} (Λ))^{- 1}) T^{(m)} (\bar{Λ})^{*}

The next statement is an immediate consequence of the above reasonings.

Lemma 8.1. For any

α > 0

and

m \in N

and for any

Λ / \in [1 / 2, \infty)

we have

rank ((A_{α} - Λ)^{- 1} - ({\hat{A}}_{α}^{(m)} - Λ)^{- 1}) = m

and therefore,

\begin{matrix} σ_{a . c .} ({\hat{A}}_{α}^{(m)}) = σ_{a . c .} (A_{α}), m_{a . c .} (E; {\hat{A}}_{α}^{(m)}) = m_{a . c .} (E; A_{α}) a . e . \end{matrix}

(8.3)

Now it is easy to conclude the proof of Theorem 5.1 for

E \geq 1 / 2

Take

m \in N

such that

E < m + 1 / 2

. Theorem 3.1 evidently applies to the matrices

J_{0}^{(m)}

, and the scheme developed in Section 7 works for the operator

{\hat{A}}_{α}^{(m)}

without any change. One only has to keep in mind that now in the corresponding version of Lemma 7.3 one can take

E < m + 1 / 2

. As a result, we obtain that

σ_{a . c .} (A_{α}^{(m)}) = σ_{a . c .} (A_{0}^{(m)}) \cup σ_{a . c .} (J_{0}^{(m)} (μ)), m_{a . c .} (E; A_{α}^{(m)}) = m_{a . c .} (E; A_{0}^{(m)}) + m_{a . c .} (E; J_{0}^{(m)} (μ))

Taking into account the equalities 8.1 and 8.3 , we arrive at the desired result.

Note that for

α < \sqrt{2}

it is easy to prove Theorem 5.1 by means of the quadratic form approach, cf. proof of Theorem 6.2 in [13] . However, in the present paper we decided to give a unified exposition for all values of the parameter.

9 Concluding remarks

In the model suggested by Smilansky in [12] the operators

A_{α}

act in the space

L^{2} (Γ \times R)

where

Γ

is a metric graph (in another terminology, a quantum graph). The model is interpreted as “harmonic oscillator, attached to a graph”. In order to describe the setting of this, more general problem, consider first the case when

Γ = Γ_{m}

is a star graph with

m

bonds, each of infinite length. More precisely,

Γ_{m}

is the union of

m

half-lines

B_{1}, \dots, B_{m}

, emanating from the common vertex

o

, the root of the tree. Let

t \in [0, \infty)

stand for the coordinate along each bond. The value

t = 0

corresponds to the root

o

. Each function

f

Γ_{m}

can be viewed as a family of

m

functions

f_{j} = f |_{B_{j}}

defined on

[0, \infty)

. If each

f_{j}

has the derivative at

t = 0

, we set

[f^{'}] (o) = \sum_{j = 1}^{m} f_{j}^{'} (0) .

In the case considered, the operator

A_{α}

L^{2} (Γ_{m} \times R)

is defined by the differential expression 1.1 for

(x, q) \in B_{j} \otimes R, j = 1, \dots, m

and the matching condition

\begin{matrix} [U_{x}^{'} (., q)] (o) = α q U (o, q), q \in R . \end{matrix}

(9.1)

The real axis

R

with the marked point

o = 0

can be identified with the graph

Γ_{2}

, and in this case the condition 9.1 turns into 1.2 .

All the results of the present paper extend to the star graphs

Γ_{m}

with an arbitrary

m > 0

, with only minor changes: 1) the equality in 2.4 has to be replaced by

m_{a . c .} (E; A_{0}) = m n

; 2) the borderline point between the small and the large values of

α

is now

m / \sqrt{2}

, and the expression for

μ

becomes

μ = m (α \sqrt{2})^{- 1}

. The equalities 1.3 and 1.4 survive. The proofs basically remain the same, but the technical calculations sometimes become rather lengthy. This was the only reason, why we restricted ourselves to the case

m = 2

in the main part of this paper.

In a similar way, the case when

Γ

is a general star graph with

m

bonds can be considered. Some of the bonds (say,

m_{0}

where

0 \leq m_{0} \leq m

) are supposed to be of infinite, and other of finite length. The Dirichlet boundary condition is imposed at the ends of the finite bonds. The point spectrum and the essential spectrum of the operators

A_{α}

were considered in [13] for this case. The a.c. spectrum can be analyzed by means of the same approach as in the present paper, but somewhat more serious changes in the formulations are necessary. They stem from the fact that now in the analog of 2.4 we have

m_{a . c .} (E; A_{0}) = m_{0} n

. In particular, if

m_{0} = 0

, that is if the graph is compact, then the spectrum of

A_{0}

is discrete. It remains discrete for

A_{α}

with

α < m / \sqrt{2}

, but its absolutely continuous component fills the whole of

R

for

α \geq m / \sqrt{2}

The case of an arbitrary metric graph with a finite number of bonds can be also analyzed, but this requires a bit more advanced technical tools. Still, the main ideas remain the same. This material will be presented elsewhere.

A Proof of Theorem 3.1

The proof is based upon the Gilbert – Pearson theory [4] of subordinate solutions, or more exactly upon the version of this theory for Jacobi matrices, see [8] .

First of all, we have to consider the homogeneous equation

J_{0} (μ) h = z h,

\begin{matrix} d_{n + 1} C (n + 1) + ((2 n + 1) μ - z) C (n) + d_{n} C (n - 1) = 0 . \end{matrix}

(A.1)

This equation is similar to 3.9 , it can be written in the form 3.16 , and moreover, for the coefficients in the decompositions 3.17 we have

a_{0} = 2 μ, a_{1} = - (μ + z); b_{0} = 1, b_{1} = - 1 .

Comparing this with 3.21 , we conclude that the asymptotic formulas 3.22 apply to the system A.1 if we take

Λ = z / μ

1. If

μ > 1

, the formula 3.22 shows that for any

E \in R

equation A.1 has a subordinate solution. Namely, this is the solution with the sign “

+

” in the exponent. It follows that the spectrum of

J_{0} (μ)

is discrete.

It is easy to show that for

μ > 1

the matrix

J_{0} (μ)

is positive definite, hence its eigenvalues tend to

+ \infty

2. If

μ < 1

, the formula 3.22 shows that for any

E \in R

equation A.1 has no subordinate solution, since

| C^{\pm} (n) | \sim n^{- 1 / 2}

As for any self-adjoint Jacobi matrix, the spectrum of

J_{0} (μ)

is simple, and the element

e_{0} = {1, 0, 0 \dots}^{⊤}

can be taken as the generating vector. Let

ℰ

stand for the spectral measure of the operator

J_{0} (μ)

, then there exists a non-negative function

τ \in L^{1} (R)

such that

{(ℰ (Δ) e_{0}, e_{0})}_{ℓ^{2}} = \int_{Δ} τ (λ) d λ

for an arbitrary Borelian set

Δ \subset R

. Our aim is to show that

τ (λ) \neq 0

a.e. on

R

Let

m_{\infty} (z)

be the Weyl function for the equation

J (μ) h - z h = 0

. It is a Herglotz function, and therefore it has boundary limits

m_{\infty} (E + i 0)

for a.e.

E \in R

. It follows from [8] , Theorem 1 that

Im m (E + i 0) \neq 0

a.e. According to the formula (5) from [8]

(J_{0} (μ) - z)^{- 1} e_{0}, e_{0}) = {(μ - z + d_{1} m_{\infty} (z))}^{- 1} .

Recall that by 3.10

d_{1} = (3 / 4)^{1 / 4}

. By the Spectral Theorem,

(J_{0} (μ) - z)^{- 1} e_{0}, e_{0}) = \int_{R} \frac{τ (λ) d λ}{λ - z} .

The last two equalities imply (when

z = E + i ɛ

and

ɛ \to + 0

) that

Im {(μ - z + d_{1} m_{\infty} (z))}^{- 1} \to \frac{- d_{1} Im m_{\infty} (E + i 0)}{| μ - E + d_{1} m_{\infty} (E + i 0) |^{2}} \neq 0 .

This limit is equal to

π τ (E)

. Therefore,

τ (E) > 0

a.e. which shows that

m (E; J_{0} (μ)) = 1

a.e.

3. The matrix

J_{0} (1)

is non-negative, therefore its spectrum lies on

[0, \infty)

. The formula 3.22 shows that for

E > 0

the equation

J_{0} (1) h = E h

has no subordinate solution. The equality

m_{a . c .} (E; J_{0} (1)) = 1

a.e.

R_{+}

can be proved in the same way as in the previous case.

This concludes the proof of Theorem 3.1 .

B Proof of Theorem 7.1

B.1 Remarks on the a.c. spectrum and on its multiplicity.

Before giving the proof, we present some remarks of a rather general nature, concerning the notion of a.c. spectrum. We consider this useful because of the dual nature of the a.c. spectrum. Indeed, it combines some features coming from the measure theory with another ones, coming from topology of the real line.

Let

K

be a self-adjoint operator whose spectrum is purely a.c., of multiplicity one. This means that

K

is unitary equivalent to the operator of multiplication,

u (λ) \mapsto λ u (λ)

, in the space

L^{2} (X; d λ)

where

X \subset R

is some Borelian set and

d λ

is the Lebesgue measure. The spectrum of

K

is the closure

\bar{X}

of the set

X

, and it may happen that

meas (\bar{X} \ X) > 0

. For this reason, characterization of the a.c.

spectrum of

K

has to include description of both the set

σ_{a . c .} (K)

and the multiplicity function

m_{a . c .} (λ; K)

. Say, in the above example

σ_{a . c .} (K) = \bar{X}

and the multiplicity function is equal to one a.e. on

X

and to zero a.e. outside

X

. It is clear, how this extends to the case of spectrum of higher, or of varying multiplicity.

Note that any Borelian set

X \subset R

, such that

m_{a . c .} (λ; K) > 0

a.e. on

X

and

m_{a . c .} (λ; K) = 0

a.e. on

R \ X

, in the book [15] is called the core of

σ_{a . c .} (K)

B.2 Preparatory material.

We need some facts from the general theory of spectral measure.

Let

A

be a self-adjoint operator in a separable Hilbert space

ℋ

and

ℰ

be its spectral measure. For elements

φ, ψ \in ℋ

, we denote by

ρ_{φ, ψ}

the scalar complex-valued measure

ρ_{φ, ψ} (\cdot) = (ℰ (\cdot) φ, ψ)

and by

ρ_{φ, ψ}^{'}

its Radon Nikodym derivative with respect to the Lebesgue measure.

This derivative is defined a.e. on

R

. It is equal to zero if either of the elements

φ, ψ

is orthogonal to the absolutely continuous subspace of the operator

A

For any

φ, ψ \in ℋ

there exists a subset

X_{φ, ψ} \subset R

of zero Lebesgue measure, such that

\begin{matrix} {lim}_{Λ \to E + i 0} ((A - Λ)^{- 1} φ - (A - \bar{Λ})^{- 1} φ, ψ) = 2 π i ρ_{φ, ψ}^{'} (E), E / \in X_{φ, ψ} . \end{matrix}

(B.1)

In such cases we do not use the notation

[\cdot] (E)

, since the exceptional set

X_{φ, ψ}

depends on the chosen elements. The set

X_{φ, ψ}

can be chosen in such a way that the relation B.1 is satisfied for the pairs

{φ, φ}

{ψ, ψ}

and

{φ, ψ}

simultaneously, then it is satisfied also for any pair from the linear hull of the elements

φ

and

ψ

Suppose now that

G

is a bounded linear operator, such that the jump 7.1 exists a.e. on

Δ

where

Δ \subset R

is a given Borelian set. This implies that for any pair

φ, ψ \in ℋ

the limit

{lim}_{Λ \to E + i 0} ((A - Λ)^{- 1} G φ - (A - \bar{Λ})^{- 1} G φ, G ψ) = (V_{A, G} (E) φ, ψ)

does exist for

E / \in X_{G}

where

X_{G} \subset Δ

is a set of Lebesgue measure

0

. Unlike B.1 , the set

X_{G}

does not depend on the choice of

φ

and

ψ

It follows from B.1 that necessarily

\begin{matrix} (V_{A, G} (E) φ, ψ) = 2 π i ρ_{G φ, G ψ}^{'} (E) a . e . on Δ . \end{matrix}

(B.2)

Again, here the exceptional set may depend on

φ

and

ψ

. However, it is important that according to our assumption, for

E / \in X_{G}

the expression in the left-hand side of B.2 , and thus the one in the right-hand side, is the sesqui-linear form of a bounded operator.

Take a dense countable set

{φ_{n}}

of elements in

ℋ

, then there exists a subset

X^{'} \subset Δ

of the Lebesgue measure

0

, such that for

E \in Δ \ X^{'}

the equality B.2 is satisfied for all pairs

{φ_{n}, φ_{m}}

, and therefore for any

φ, ψ

from the linear hull

ℳ

of the system

{φ_{n}}

. In other words, there exists a dense linear subspace

ℳ \subset ℋ

, such that B.2 is satisfied for all

E \in Δ \ X^{'}

and for all

φ, ψ \in ℳ

simultaneously.

B.3 Proof of the theorem.

Taking an appropriate partition of the original set, we may assume that

m_{a . c .} (E; A) = ν = c o n s t

a.e. on

Δ

. Suppose first that

ν < \infty

. According to the general theory of spectral measure, see e.g. [2] , there is a subspace

ℋ_{Δ} \subset ℋ

invariant with respect to

A

, isometric to

L^{2} (Δ, ℓ_{ν}^{2}; d x)

and such that on

ℋ_{Δ}

the operator

A

acts as multiplication by

x

. More precisely, let

Π : ℋ_{Δ} \to L^{2} (Δ, ℓ_{ν}^{2}; d x)

stand for the above isometry. Extending it by zero to the orthogonal complement of

ℋ_{Δ}

, we obtain a partially isometric operator

Π^{\circ} : ℋ \to L^{2} (Δ, ℓ_{ν}^{2}; d x)

. If

φ \in D o m A \cap ℋ_{Δ}

, then

A φ \in ℋ_{Δ}

and

(Π A φ) (x) = x (Π φ) (x) .

For any

φ, ψ \in ℋ_{Δ}

and any Borelian subset

\partial \subset Δ

we have

ρ_{φ, ψ} (\partial) = \int_{\partial} {((Π φ) (x), (Π ψ) (x))}_{ℓ_{ν}^{2}} d x

and

ρ_{φ, ψ}^{'} (E) = {((Π φ) (E), (Π ψ) (E))}_{ℓ_{ν}^{2}}, a . e . on Δ .

Moreover,

ρ_{φ, ψ}^{'} (E) = 0

a.e. on

Δ

, provided that either of the elements

φ, ψ

is orthogonal to

ℋ_{Δ}

. Now, the formula B.2 and the above remarks imply that there exists a dense linear subset

ℳ \subset ℋ

, such that for all

φ, ψ \in ℳ

the equality

\begin{matrix} (V_{A, G} (E) φ, ψ) = 2 π i {((Π^{\circ} G φ) (E), (Π^{\circ} G ψ) (E))}_{ℓ_{ν}^{2}} \end{matrix}

(B.3)

is satisfied a.e. on

Δ

, and the exceptional subset does not depend on the choice of

φ, ψ

. The expression in the right-hand side is the sesqui-linear form of an operator in

ℓ_{ν}^{2}

. Necessarily, its rank does not exceed

ν

. Therefore, the restriction of the operator

V_{A, G} (E)

ℳ

also has rank no greater than

ν

. Since the linear set

ℳ

is dense in

ℋ

and the operator

V_{A, G} (E)

is bounded, we find that

rank V_{A, G} (E) \leq ν

In order to obtain the reverse inequality, take elements

f_{1}, \dots, f_{ν} \in ℋ_{Δ}

, such that

(Π f_{j}) (x) = χ_{Δ} (x) e_{j}, j = 1, \dots, ν

where

χ_{Δ}

is the characteristic function of the set

Δ

and the vectors

{e_{j}}

form the natural basis in

ℓ_{ν}^{2}

. Since the range of

G

is assumed dense, we can choose elements

φ_{j} \in ℋ

in such a way that

∥ G φ_{j} - f_{j} ∥_{ℋ} < ɛ, j = 1, \dots, ν,

where

ɛ > 0

is arbitrarily small. Then

\int_{Δ} ∥ (Π^{\circ} G φ_{j}) (x) - χ_{Δ} (x) e_{j} ∥_{ℓ_{ν}^{2}}^{2} d x = ∥ Π^{\circ} (G φ_{j} - f_{j}) ∥_{ℋ}^{2} < ɛ^{2} .

This yields that the vectors

(Π^{\circ} G φ_{j}) (E) \in ℓ_{ν}^{2}

are linearly independent for

E

lying outside a subset

X_{ɛ}

of a small measure

η (ɛ)

, and

η (ɛ) \to 0

ɛ \to 0

. For any

E \in Δ \ X_{ɛ}

we have

rank V_{A, G} (E) = ν

. Letting

ɛ \to 0

, we conclude from B.3 that

rank V_{A, G} = ν

a.e. on

Δ

The equality 7.2 is justified for any

ν < \infty

. Therefore, it remains valid also if

ν = \infty

Acknowledgments. The work on the paper started when S.N. visited the Weizmann Institute in November – December of 2003. The visit was supported in part by the Department of Mathematics, and in part by the Einstein center for theoretical physics. S.N. takes this opportunity to express his gratitude to the Institute for its hospitality and financial support.

M.S. acknowledges partial financial support of the network SPECT of the ESF. References

M.Sh. Birman and S.B. Entina, Stationary approach in abstract scattering theory (Russian) Izv. Akad. Nauk SSSR Ser. Mat. 31, no. 2 (1967), 401–430. English translation in: Math. USSR Izvestija 1, no. 1 (1967), 391–420.
M.Sh. Birman and M. Solomyak, Spectral theory of selfadjoint operators in Hilbert space. D. Reidel Publishing Co., Dordrecht, 1987.
S.N. Elaydi, An introduction to difference equations, Springer, New York, 1999.
D.J. Gilbert and D.B. Pearson, On subordinacy and analysis of the spectrum of one-dimensional Schrdinger operators, J. Math. Anal. Appl. 128 (1987), no. 1, 30–56.
I.C. Gohberg and M.G. Krein, Introduction to the theory of linear non-selfadjoint operators in Hilbert space. Izdat. “Nauka”, Moscow 1965 448 pp. English translation: Amer. Math. Soc., Providence (1969).
J. Janas and S.N. Naboko, Spectral Analysis of Selfadjoint Jacobi matrices with Periodically Modulated Entries, J. Funct. Anal. 191 (2002), no. 2, 318–342.
T. Kato, Perturbation theory for linear operators, Springer, Berlin, 1966.
S. Khan and D.B. Pearson, Subordinacy and spectral theory for infinite matrices, Helv. Phys. Acta 65 (1992), no. 4, 505–527.
S.N. Naboko, Conditions for the existence of wave operators in the non-selfadjoint case, Wave propagation. Scattering theory, Topics in Mathematical Physics 12 (1987) (Russian), 127–149. English translation in: AMS Traslations, Series 2, vol. 157 (1993), 132–155.
S.N. Naboko,Uniqueness theorems for operator-valued functions with positive imaginary part, and the singular spectrum in the self-adjoint Friedrichs model, Arkiv. Math. 25 (1987), 116–140.
M. Reed and B. Simon, Methods of modern mathematical physics. I. Functional analysis, Academic Press, New York-London, 1972.
U. Smilansky, Irreversible quantum graphs, Waves in Random Media, 14 (2004), 143 – 153.
M. Solomyak, On a differential operator appearing in the theory of irreversible quantum graphs, Waves in Random Media 14 (2004), 173 – 185.
M. Solomyak, On the discrete spectrum of a family of differential operators, Functional alalysis and its applications, 38 (2004), 217-223.
D.R. Yafaev, Mathematical scattering theory. General theory, Translations of Mathematical Monographs. 105. Providence, RI: American Mathematical Society, 1992.

Department of Mathematical Physics, St.Petersburg State University, St. Petergoff 198904 St. Peterburg, Russia E-mail address : naboko@snoopy.phys.spbu.ru Department of Mathematics, The Weizmann Institute of Science, Rehovot 76100, Israel E-mail address : michail.solomyak@weizmann.ac.il

25.05.2004

Sergey N. Naboko

Michael Solomyak