, , ,

int In the paper [7] U. Smilansky suggested a mathematical model to which he gave the name “Irreversible quantum graph”. In this model one studies the interaction between a quantum graph and a finite system of one-dimensional oscillators attached at several different points in the graph. Recall that the term “quantum graph” usually stands for a metric graph

Γ

equipped with a self-adjoint differential operator acting on

L^{2} (Γ)

; see the survey paper [4] and references therein. In our case this operator will be the Laplacian

- Δ

In Smilansky's model one initially deals with two independent dynamical systems. One of the systems acts in

L^{2} (Γ)

and its Hamiltonian is the Laplacian. Another system acts in the space

L^{2} (R^{K}), K \geq 1

and is generated by the Hamiltonian

H_{o s c} = \sum_{k = 1}^{K} h_{k}

where

h_{k} = \frac{ν_{k}^{2}}{2} (- \frac{\partial^{2}}{\partial {q_{k}}^{2}} + q_{k}^{2}), k = 1, \dots, K;

in [7] the oscillators are written in a slightly different form; one form reduces to another by scaling. In what follows the points in

Γ

are denoted by

x

and the points in

R^{K}

q = (q_{1}, \dots, q_{K})

Consider now the operator

\begin{matrix} A_{0} = - Δ \otimes I + I \otimes H_{o s c} 1.1 \end{matrix}

(1.1)

in the space

L^{2} (Γ \times R^{K})

. It is defined by the differential expression

\begin{matrix} A U = - Δ_{x} U + \frac{1}{2} \sum_{k = 1}^{K} ν_{k}^{2} (- \frac{\partial^{2} U}{\partial {q_{k}}^{2}} + q_{k}^{2} U) 1.2 \end{matrix}

(1.2)

and is self-adjoint on the natural domain. The terms in 1.1 do not interact with each other.

Interaction is introduced with the help of a system of “matching conditions” on the derivative

U_{x}^{'}

at some points

o_{1}, \dots, o_{K} \in Γ

. One says that the

k

-th oscillator is attached to the graph at the point

o_{k}

. The condition at the point

o_{k}

\begin{matrix} [U_{x}^{'}] (o_{k}, q) = α_{k} q_{k} U (o_{k}, q), k = 1, \dots, K, 1.4 \end{matrix}

(1.3)

where

[f_{x}^{'}] (.)

stands for the expression appearing in the Kirchhoff condition, well known in the theory of electric networks. When

Γ = R

(which is the only case we deal with in the main body of the paper),

[f_{x}^{'}] (.)

is the jump of the derivative,

\begin{matrix} [f_{x}^{'}] (o) = f_{x}^{'} (o +) - f_{x}^{'} (o -) . 1.5 \end{matrix}

(1.4)

The real parameter

α_{k}

in 1.3 expresses the strength of interaction between the quantum graph and the oscillator

h_{k}

. The case

α_{1} = \dots = α_{K} = 0

corresponds to the operator

A_{0}

as in 1.1 .

Sometimes we shall denote by

α, ν

the multi-dimensional parameters

α = {α_{1}, \dots, α_{K}}

ν = {ν_{1}, \dots, ν_{K}}

. Let

A_{α; ν} = A_{α_{1}, \dots, α_{K}; ν_{1}, \dots, ν_{K}}

stand for the operator defined by the differential expression 1.2 and the conditions 1.3 .

Usually, the values of

ν_{k}

are fixed and we exclude them from the notation.

On the other hand, we use the notation

A_{Γ; α; ν}

for this operator when it is necessary to reflect its dependence on the graph.

The problem to be considered is the description of the spectrum of the dynamical system generated by the Hamiltonian

A_{α; ν}

. More specifically, it is to construct the self-adjoint realization of

A_{α; ν}

as an operator in the Hilbert space

L^{2} (Γ \times R^{K})

and to describe its spectrum.

Up until now, the problem has only been investigated for the simplest case

K = 1

. The first results were obtained in the paper [7] by Smilansky. Then a detailed study of the problem was carried out in the papers [8] , [9] and [5] . In [10] , along with some new results, a detailed survey of the current state of the problem is given.

On first sight, the problem might seem amenable to the perturbation theory of quadratic forms. Indeed, the spectrum

σ (A_{0})

can be easily described by separation of variables and the perturbation in the quadratic form, which appears when passing from

A_{0}

A_{α}

with

α \neq 0

, seems not to be too strong.

However, this is not so: this perturbation turns out to be only form-bounded but not form-compact, which makes it impossible to apply the standard techniques.

So, the problem requires certain specific tools which were developed in [8] – [10] and [5] . The most important of these tools is the systematic use of Jacobi matrices.

It was found in the above mentioned papers on the one-oscillator problem that the character of the spectrum strongly depends on the size of

α

: there exists some

α^{*} > 0

such that the absolutely continuous spectrum

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

coincides with

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

| α | < α^{*}

(in particular, it is absent if the graph is compact) and fills the whole of

R

| α | > α^{*}

. The dependence of the structure of the point spectrum

σ_{p} (A_{α})

α

is also well understood.

This is the first of two papers on the problem for

K > 1

and in it we study the absolutely continuous spectrum; in our other paper [1] the point spectrum is investigated. This division is natural, since the technical tools used in each part are different. We address the simplest situation, when

Γ = R

and

K = 2

, but our argument is of a rather general character and we firmly believe that it applies to a wide class of graphs and to any

K

. However, in the general case, the calculations become more complicated and this obscures the main features of the argument.

We first describe informally the main ideas lying behind our approach.

The effect of adding one more oscillator to a system with

K

oscillators is twofold. Firstly, the total dimension of the set

Γ \times R^{K}

increases by one which certainly affects the spectrum. Secondly, there is some effect coming from the additional matching condition 1.3 at the point

o_{K + 1}

. This second effect disappears if we take

α_{K + 1} = 0

. Indeed, then the variable

q_{K + 1}

can be separated and the operator decomposes into the orthogonal sum of simpler operators. More exactly, denote by

\tilde{A}

the operator which corresponds to the configuration with the

(K + 1)

-th oscillator removed,

\tilde{A} = A_{α_{1}, \dots, α_{K}; ν_{1}, \dots, ν_{K}} .

Then it is easy to see that

\begin{matrix} A_{α_{1}, \dots, α_{K}, 0; ν_{1}, \dots, ν_{K}, ν_{K + 1}} = {\sum_{n \in N_{0}}}^{\oplus} (\tilde{A} + ν_{K + 1}^{2} (n + 1 / 2)) . 1.6 \end{matrix}

(1.5)

This orthogonal decomposition yields the complete description of the spectrum of the operator on the left-hand side, provided that the spectrum of

\tilde{A}

is known.

The key observation which allows one to solve the general problem is that the interaction between the oscillators attached at different points is weak. For

K = 2

this observation leads to the conclusion that the study of

σ (A_{α_{1}, α_{2}; ν_{1}, ν_{2}})

can be reduced to the same problem for the operators

A_{α_{1}, 0; ν_{1}, ν_{2}}

and

A_{0, α_{2}; ν_{1}, ν_{2}}

Due to the equality 1.5 this reduces the problem to the study of the spectra of two operators,

A_{α_{1}; ν_{1}}

and

A_{α_{2}; ν_{2}}

, each corresponding to the case of only one oscillator. Since the latter case is already well understood, we obtain the desired results for our more complicated case.

An accurate realization of this idea is different for the point spectrum and for the absolutely continuous spectrum. In the present paper we concentrate on the absolutely continuous spectrum. Here an important correction to the above scheme is necessary: the study of

σ_{a . c .} (A_{α_{1}, α_{2}; ν_{1}, ν_{2}})

does not reduce to the study of

σ_{a . c .} (A_{α_{1}; ν_{1}})

and

σ_{a . c .} (A_{α_{1}; ν_{1}})

for the same graph

Γ

. Rather, we have to divide

Γ

into two parts,

Γ = Γ_{1} \cup Γ_{2}

in such a way that

o_{j} \in Γ_{j}

and

o_{j} / \in Γ_{3 - j}

. Then

σ_{a . c .} (A_{Γ; α_{1}, α_{2}; ν_{1}, ν_{2}})

can be expressed in terms of

σ_{a . c .} (A_{Γ_{j}; α_{j}; ν_{j}})

j = 1, 2

. The paper [1] is devoted to the study of the point spectrum. There such a partition of

Γ

is unnecessary.

We use the following notation. We write

N_{0}

for the set

{0, 1, \dots}

. The diagonal operator in an appropriate

ℓ^{2}

-space, with the diagonal elements

a_{0}, a_{1}, \dots

, is denoted by

d i a g {a_{n}}

. We apply similar notation for the block-diagonal operators. The notation

J ({a_{n}}, {b_{n}})

stands for the Jacobi matrix whose non-zero entries are

j_{n, n} = a_{n}

and

j_{n, n + 1} = j_{n + 1, n} = b_{n}

. If

A

is a self-adjoint operator in a Hilbert space, then

σ (A), σ_{a . c} (A), σ_{p} (A)

stand for its spectrum, absolutely continuous (a.c.) spectrum and point spectrum respectively. We use the symbol

m_{a . c} (λ; A)

for the multiplicity function of the a.c. spectrum.

The symbol

S_{1}

stands for the trace class of compact operators.

Other necessary notations are introduced in the course of the presentation.

2 Statement of the problem. Results

2.1 The operator $A_{α}$ .

st1 As was mentioned in the introduction, we present our argument for the graph

Γ = R

and

K = 2

. We choose the points

o_{1} = 1, o_{2} = - 1

and denote the coordinates in

R^{2}

q_{+}, q_{-}

and the parameters by

α = {α_{+}, α_{-}}

ν = {ν_{+}, ν_{-}}

The Laplacian on

Γ

is just the operator

- d^{2} / d x^{2}

with the Sobolev space

H^{2} (R)

as the operator domain. The operator

A_{α, ν}

acts in the Hilbert space

ℌ = L^{2} (R^{3})

and is defined by the differential expression

\begin{matrix} A U = A_{ν} U = - U_{x^{2}}^{''} + \frac{ν_{+}^{2}}{2} (- U_{q_{+}^{2}}^{''} + q_{+}^{2} U) + \frac{ν_{-}^{2}}{2} (- U_{q_{-}^{2}}^{''} + q_{-}^{2} U) st.1 \end{matrix}

(2.1)

and the matching conditions (cf. 1.4 )

\begin{matrix} [U_{x}^{'}] (\pm 1, q_{+}, q_{-}) = α_{\pm} q_{\pm} U (\pm 1, q_{+}, q_{-}) . st2 \end{matrix}

(2.2)

So, in the notation of the introduction, we are dealing with the operator

A_{R; α_{+}, α_{-}; ν_{+}, ν_{-}} .

However, as a rule we use the shortened notation

A_{α}

. Note that the replacement

α_{\pm} \mapsto - α_{\pm}

corresponds to the change of variables

q_{\pm} \mapsto - q_{\pm}

which does not affect the spectrum. For this reason, we discuss only

α_{\pm} \geq 0

The structure of the differential expression

A

makes it natural to decompose the function

U

in a double series in terms of the normalized Hermite functions

χ_{n}

, namely

\begin{matrix} U (x, q_{+}, q_{-}) = \sum_{m, n \in N_{0}} u_{m, n} (x) χ_{m} (q_{+}) χ_{n} (q_{-}), sa.0 \end{matrix}

(2.3)

which is hereafter represented by

U \sim {u_{m, n}}

. The mapping

U \mapsto {u_{m, n}}

is an isometry of the space

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

onto the Hilbert space

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

. We evidently have

A U \sim {L_{m, n} u_{m, n}}

where

\begin{matrix} (L_{m, n} u) (x) = - u^{''} (x) + r_{m, n} u (x), x \neq \pm 1; sa.1 \end{matrix}

(2.4)

\begin{matrix} r_{m, n} = ν_{+}^{2} (m + 1 / 2) + ν_{-}^{2} (n + 1 / 2), m, n \in N_{0} . sa.2 \end{matrix}

(2.5)

The conditions at

x = \pm 1

reduce to

\begin{matrix} \begin{matrix} [u_{m, n}^{'}] (1) & = \frac{α_{+}}{\sqrt{2}} (\sqrt{m + 1} u_{m + 1, n} (1) + \sqrt{m} u_{m - 1, n} (1)); & [u_{m, n}^{'}] (- 1) & = \frac{α_{-}}{\sqrt{2}} (\sqrt{n + 1} u_{m, n + 1} (- 1) + \sqrt{n} u_{m, n - 1} (- 1)) . \end{matrix} sa.3 \end{matrix}

(2.6)

To derive the conditions 2.6 from 2.2 , one uses the recurrency equation for the functions

χ_{n}

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

2.2 Operator $A_{0}$ .

op0 The operator

A_{0} : = A_{0, 0; ν_{+}, ν_{-}}

admits separation of variables and we get

\begin{matrix} A_{0} = \sum_{m, n}^{\oplus} (H_{0} + r_{m, n}), add1 \end{matrix}

(2.7)

where

H_{0}

is the self-adjoint operator

- d^{2} / d x^{2}

L^{2} (R)

. This leads to the complete description of the spectrum

σ (A_{0})

, namely, that it is purely a.c.

and fills the half-line

[r_{0, 0}, \infty) = [(ν_{+}^{2} + ν_{-}^{2}) / 2, \infty)

. The expression for the multiplicity function

m_{a . c .} (λ; A_{0})

immediately follows from 2.7 , but is omitted.

2.3 Domain of $A_{α}$ .

opa It is convenient to describe the domain of the self-adjoint realization of the operator

A_{α}

in terms of the decomposition 2.3 . Define the set

D_{α}

as follows.

Definition 2.1. An element

U \sim {u_{m, n}}

lies in

D_{α}

if and only if d 1.

u_{m, n} \in H^{1} (R)

for all

m, n

2. For all

m, n

the restriction of

u_{m, n}

to each interval

(- \infty, - 1), (- 1, 1)

(1, \infty)

lies in

H^{2}

and moreover,

\sum_{m, n} \int_{R} | L_{m, n} u_{m, n} |^{2} d x < \infty .

3. The conditions 2.6 are satisfied.

Along with the set

D_{α}

, define its subset

D_{α}^{∙} = {U \in D_{α} : U \sim {u_{m, n}} f i n i t e}

where by finite we mean that the sequence has only a finite number of non-zero components. Denote by

A_{α}^{∙}

the operator in

ℌ = L^{2} (R^{3})

, defined by the system 2.4 on the domain

D_{α}^{∙}

Lemma 2.2. The operator

A_{α}^{∙}

is symmetric in

ℌ

. Its adjoint l1 coincides with the operator

A_{α}

considered on the domain

D_{α}

(A_{α}^{∙})^{*} = A_{α} .

The proof is a straightforward modification of that for (5.2) in [5] .

Theorem 2.3. For any

α_{+}, α_{-} \geq 0

the operator

A_{α}

is self-adjoint.

sa.thm

The proof is given in section 4 . Theorem 2.3 and Lemma 2.2 show that

A_{α}

is the unique natural self-adjoint realization of the operator, defined by the differential expression 2.1 and the matching conditions 2.2 .

2.4 Absolutely continuous spectrum of the operator $A_{α}$ .

circ Below we construct an operator

A_{α}^{\circ}

whose a.c. spectrum admits a complete description. Then we show that the a.c. spectra of both operators

A_{α}

and

A_{α}^{\circ}

coincide, including the multiplicities.

As a first step, let us consider two operators,

A_{α_{+}}^{+}

and

A_{α_{-}}^{-}

. The operator

A_{α_{+}}^{+}

, say, acts in the space

L^{2} (R_{+} \times R^{2})

and is defined by the differential expression 2.1 , the matching condition 2.2 at the point

o_{1} = 1

and the Dirichlet condition

U (0, q_{+}, q_{-}) = 0

. The definition of

A_{α_{-}}^{-}

is similar, with

R_{+}

replaced by

R_{-}

and the point

o_{1} = 1

o_{2} = - 1

. By separation of variables, the operators

A_{α_{\pm}}^{\pm}

can be identified with the orthogonal sum of simpler operators:

\begin{matrix} \begin{matrix} A_{α_{+}}^{+} = {\sum_{n \in N_{0}}}^{\oplus} (A_{R_{+}; α_{+}; ν_{+}} + ν_{-}^{2} (n + 1 / 2)) & , & A_{α_{-}}^{-} = {\sum_{m \in N_{0}}}^{\oplus} (A_{R_{-}; α_{-}; ν_{-}} + ν_{+}^{2} (m + 1 / 2)) & . \end{matrix} one.2 \end{matrix}

(2.8)

Hence, both operators are self-adjoint. The direct sum

\begin{matrix} A_{α}^{\circ} = A_{α_{+}, α_{-}; ν_{+}, ν_{-}}^{\circ} : = A_{α_{+}}^{+} \oplus A_{α_{-}}^{-} one.3 \end{matrix}

(2.9)

is a self-adjoint operator in the original Hilbert space

ℌ

The following theorem is the main result of the paper. Its formulation involves the notion of wave operator, which is one of the basic notions in mathematical scattering theory; see e.g. [3] , [6] or [11] .

Theorem 2.4. For each of the pairs

(A_{α}, A_{α}^{\circ}), (A_{α}^{\circ}, A_{α})

, there exist one.t1 complete isometric wave operators. In particular, the absolutely continuous parts of

A_{α}

and

A_{α}^{\circ}

are unitarily equivalent.

Theorem 2.4 and the formulae 2.8 , 2.9 reduce the study of

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

to the similar problem for the case of only one oscillator. The latter problem was solved in [5] and [10] . The next statement collects, for the particular case we need, the results of section 3 in [10] ; see also Theorem 5.1 and remarks in section 9 of [5] . In both papers it was assumed that

ν = 1

, and we arrive at the formulation below via scaling. By default, we take

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

λ / \in σ_{a . c .} (A)

Proposition 2.5. (The case of one oscillator.) Let

Γ = R_{+}

and

o = 1

, prop or

Γ = R_{-}

and

o = - 1

. Then 1)

σ_{a . c .} (A_{0; ν}) = [ν^{2} / 2, \infty)

;

m_{a . c .} (λ; A_{0; ν}) = n f o r - ν^{2} / 2 \leq λ - ν^{2} n < ν^{2} / 2, n \in N;

2) if

0 < α < ν \sqrt{2}

, then

σ_{a . c .} (A_{α; ν}) = σ_{a . c .} (A_{0; ν}) = [ν^{2} / 2, \infty);

m_{a . c .} (λ; A_{α; ν}) = m_{a . c .} (λ; A_{0; ν});

3) if

α = ν \sqrt{2}

, then

σ_{a . c .} (A_{α; ν}) = [0, \infty); m_{a . c .} (λ; A_{α; ν}) = m_{a . c .} (λ; A_{0; ν}) + 1, \forall λ \geq 0;

4) if

ν \sqrt{2} < α < \infty

, then

σ_{a . c .} (A_{α; ν}) = R; m_{a . c .} (λ; A_{α; ν}) = m_{a . c .} (λ; A_{0; ν}) + 1, \forall λ \in R .

Now we are in a position to present the final formula for the function

m_{a . c .} (λ; A_{α}^{\circ})

, and thus for our original operator

A_{α}

\begin{matrix} \begin{matrix} m_{a . c .} (λ; A_{α}) = & \sum_{n \in N_{0}} m_{a . c .} (λ - ν_{-}^{2} (n + 1 / 2); A_{R_{+}; α_{+}; ν_{+}}) & + & \sum_{m \in N_{0}} m_{a . c .} (λ - ν_{+}^{2} (m + 1 / 2); A_{R_{-}; α_{-}; ν_{-}}) . \end{matrix} add3 \end{matrix}

(2.10)

Combining the equality 2.10 with Proposition 2.5 , we obtain the following description of the a.c. spectrum of the operator

A_{α}

for any

α_{+}, α_{-} \geq 0

Theorem 2.6. Let

A_{α} = A_{α; ν}

be the self-adjoint operator defined one.t2a by the differential expression 2.1 on the operator domain

D_{α}

1) If

α_{\pm} / ν \pm < \sqrt{2}

, then

σ_{a . c .} (A_{α}) = [r_{0, 0}, \infty) = [(ν_{+}^{2} + ν_{-}^{2}) / 2, \infty) .

2) Let

α_{+} / ν_{+} = \sqrt{2}

and

α_{-} / ν_{-} < \sqrt{2}

, or

α_{-} / ν_{-} = \sqrt{2}

and

α_{+} / ν_{+} < \sqrt{2}

. Then

σ_{a . c .} (A_{α}) = [ν_{-}^{2} / 2, \infty) o r σ_{a . c .} (A_{α}) = [ν_{+}^{2} / 2, \infty)

respectively.

3) Let

α_{+} / ν_{+} = α_{-} / ν_{-} = \sqrt{2}

, then

σ_{a . c .} (A_{α}) = [0, \infty) .

In all the cases 1 – 3 the multiplicity function

m_{a . c .} (λ; A_{α})

, given by the equality 2.10 , is finite for all

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

4) Let

m a x (α_{+} / ν_{+}, α_{-} / ν_{-}) > \sqrt{2}

. Then

σ_{a . c .} (A_{α}) = R, m_{a . c .} (λ; A_{α}) \equiv \infty .

In connection with this theorem, we would like to emphasize that the existence of the wave operators established in Theorem 2.4 gives much more information about the operator

A_{α}

than just the description of its a.c. spectrum.

For the proof of Theorem 2.4 we use the following classical result due to Kato, see Theorem 6.5.1 and Remark 6.5.2 in [11] .

Proposition 2.7. Let

A, A^{\circ}

be self-adjoint operators in a Hilbert kato space. Suppose that for some natural number

p

the inclusion

(A^{\circ} - Λ)^{- p} - (A^{\circ} - Λ)^{- p} \in S_{1}

is satisfied for all non-real

Λ \in C

. Then the complete isometric wave operators exist for both pairs

A, A^{\circ}

and

A^{\circ}, A

In our case the conditions of Proposition 2.7 turn out to be fulfilled with

p = 3

. This is the result of the following statement whose proof is our main technical goal in this paper.

Theorem 2.8. For any non-real

Λ \in C

one has

\begin{matrix} one.t2 (A_{α}^{\circ} - Λ)^{- 3} - (A_{α} - Λ)^{- 3} \in S_{1} . goal \end{matrix}

(2.11)

Theorem 2.4 is a direct consequence of Theorem 2.8 .

The proof of Theorem 2.8 is rather long and requires some preparatory work.

3 Auxiliary material

au In this section we present some elementary technical material concerning the equations

\begin{matrix} \begin{matrix} (3.1) & - u^{''} + ζ^{2} u = f, res.1g - v^{''} + ζ^{2} v = 0 . res.1go \end{matrix} \end{matrix}

(3.2)

where

ζ = γ + i δ

is a complex parameter. We need this material for the proofs of both our main technical results, Theorem 2.3 and Theorem 2.8 .

We assume that

γ > 0

, and are mainly interested in estimates which are uniform with respect to

ζ

3.1 Homogeneous equation.

he Let

ℱ_{ζ}

be the two-dimensional space of functions on

R

which are continuous, vanish as

| x | \to \infty

, and for

x \neq \pm 1

satisfy the equation 3.2 . We choose the following basis

φ_{ζ}^{+}, φ_{ζ}^{-}

ℱ_{ζ}

φ_{ζ}^{+} (x) = {\begin{matrix} 0, & x < - 1, \\ \frac{sinh ζ (x + 1)}{sinh 2 ζ}, & - 1 \leq x \leq 1; \\ e^{- ζ (x - 1)}, & x > 1; \end{matrix} φ_{ζ}^{-} (x) = φ_{ζ}^{+} (- x) .

Then

\begin{matrix} φ_{ζ}^{+} (1) = φ_{ζ}^{-} (- 1) = 1, φ_{ζ}^{+} (- 1) = φ_{ζ}^{-} (1) = 0.1 \end{matrix}

(3.3)

Just for this reason this basis is more convenient than the “natural” basis consisting of the functions

e^{- ζ | x \pm 1 |}

. Note also that

[(φ_{ζ}^{\pm})^{'}] (\pm 1) = - \frac{2 ζ}{1 - e^{- 4 ζ}}, [(φ_{ζ}^{\pm})^{'}] (\mp 1) = \frac{2 ζ e^{- 2 ζ}}{1 - e^{- 4 ζ}}

and hence, for all

v \in ℱ_{ζ}

\begin{matrix} \begin{matrix} [v^{'}] (1) = & - \frac{2 ζ}{1 - e^{- 4 ζ}} (v (1) - e^{- 2 ζ} v (- 1)), & [v^{'}] (- 1) = & - \frac{2 ζ}{1 - e^{- 4 ζ}} (v (- 1) - e^{- 2 ζ} v (1)) . \end{matrix} 1c \end{matrix}

(3.4)

A standard calculation shows that for the norm and scalar product in

L^{2} (R)

\begin{matrix} ∥ φ_{ζ}^{+} ∥^{2} = ∥ φ_{ζ}^{-} ∥^{2} = \frac{1}{2 γ} + \frac{γ^{- 1} sinh 4 γ - δ^{- 1} sin 4 δ}{4 ({sinh}^{2} 2 γ + {sin}^{2} 2 δ)} = γ^{- 1} + o (e^{- 4 γ}), 1f (φ_{ζ}^{+}, φ_{ζ}^{-}) = O (e^{- 2 γ}), γ \to \infty . \end{matrix}

(3.5)

This shows that for

γ

large the chosen basis is “almost orthogonal”. It follows that the two-sided estimate

\begin{matrix} c_{0}^{- 1} γ ∥ v ∥^{2} \leq | C_{+} |^{2} + | C_{-} |^{2} \leq c_{0} γ ∥ v ∥^{2}, v = C_{+} φ_{ζ}^{+} + C_{-} φ_{ζ}^{-} \in ℱ_{ζ} 1x \end{matrix}

(3.6)

with some

c_{0} > 1

is satisfied uniformly in any half-plane

γ = Re ζ \geq γ_{0} > 0

Now we turn to the subspace

ℱ_{ζ}^{\circ}

formed by the functions

v \in ℱ_{ζ}

, satisfying an additional condition

v (0) = 0

. The functions

φ_{ζ}^{\circ, +} (x) = {\begin{matrix} 0, & x < 0, \\ \frac{sinh ζ x}{sinh ζ}, & 0 \leq x \leq 1, \\ e^{- ζ (x - 1)}, & x > 1; \end{matrix} φ_{ζ}^{\circ, -} (x) = φ_{ζ}^{\circ, +} (- x)

form a natural basis in

ℱ_{ζ}^{\circ}

. For the functions

φ_{ζ}^{\circ, \pm}

the equalities 3.3 are satisfied, and instead of 3.4 we have

\begin{matrix} [v^{'}] (\pm 1) = - \frac{2 ζ}{1 - e^{- 2 ζ}} v (\pm 1), \forall v \in ℱ_{ζ}^{\circ} . 1ac \end{matrix}

(3.7)

Similarly to 3.5 , we find that

∥ φ_{ζ}^{\circ, \pm} ∥^{2} = γ^{- 1} + O (e^{- 2 γ}), (φ_{ζ}^{\circ, +}, φ_{ζ}^{\circ, -}) = O (e^{- γ}), γ \to \infty .

As a consequence, we conclude that an analogue of 3.6 , with the functions

φ_{ζ}^{\pm}

replaced by

φ_{ζ}^{\circ, \pm}

, is valid for

v \in ℱ_{ζ}^{\circ}

A straightforward calculation shows also that

\begin{matrix} ∥ φ_{ζ}^{\circ, \pm} - φ_{ζ}^{\pm} ∥ = O (e^{- γ}), γ \to \infty . 2c \end{matrix}

(3.8)

3.2 Non-homogeneous equation.

nhom Here we discuss the equation 3.1 without the matching conditions at

x = \pm 1

or, equivalently, under the conditions of the type 2.6 with

α = 0

. Then the solution is given by

\begin{matrix} u_{ζ} (x) = (2 ζ)^{- 1} \int_{R} e^{- ζ | x - t |} f (t) d t . res.1 \end{matrix}

(3.9)

The solution of the same equation 3.1 subject to the condition

u (0) = 0

\begin{matrix} u_{ζ}^{\circ} (x) = {\begin{matrix} (2 ζ)^{- 1} \int_{R_{+}} (e^{- ζ | x - t |} - e^{- ζ (x + t)}) f (t) d t, x > 0; \\ (2 ζ)^{- 1} \int_{R_{-}} (e^{- ζ | x - t |} - e^{ζ (x + t)}) f (t) d t, x < 0 . \end{matrix} res.1b \end{matrix}

(3.10)

The difference

u_{ζ}^{\circ} - u_{ζ}

is given by a rank one operator,

\begin{matrix} u_{ζ}^{\circ} (x) - u_{ζ} (x) = - (2 ζ)^{- 1} g_{ζ} (x) \int_{R} f (t) g_{ζ} (t) d t, g_{ζ} (x) = e^{- ζ | x |} . res.1f \end{matrix}

(3.11)

Note that

∥ g_{ζ} ∥^{2} = γ^{- 1}

3.3 Dependence on the additional parameters.

add We are particularly interested in the case when

ζ

depends on two parameters

r \in R

and

Λ \in C

, where

r \geq r_{0} > 0, Λ / \in R_{+}

\begin{matrix} ζ = ζ_{r} (Λ) : = γ_{r} (Λ) + i δ_{r} (Λ) = \sqrt{r - Λ} . 2 \end{matrix}

(3.12)

We select the branch of the square root in 3.12 to have

Re ζ_{r} (Λ) > 0, Im Λ \cdot Im ζ_{r} (Λ) \leq 0 .

For

Λ

fixed all the points

ζ_{r} (Λ)

lie in some half-plane

Re ζ_{r} (Λ) \geq γ_{0} (Λ) > 0

, hence 3.6 is satisfied. It is clear that

γ_{r} (Λ) \sim r^{1 / 2}

r \to \infty

. Therefore, for any

Λ / \in R

there exists a constant

c_{1} = c_{1} (Λ) > 1

such that

\begin{matrix} \begin{matrix} c_{1}^{- 1} r^{1 / 2} ∥ v ∥^{2} \leq | C_{+} |^{2} + | C_{-} |^{2} & \leq c_{1} r^{1 / 2} ∥ v ∥^{2}, & \forall v = C_{+} φ_{ζ}^{+} + C_{-} φ_{ζ}^{-} \in ℱ_{ζ}, & ζ = ζ_{r} (Λ) . \end{matrix} 2r \end{matrix}

(3.13)

4 Self-adjointness: proof of Theorem 2.3

pr According to the general theory of self-adjoint operators, we must show that the equation

A_{α} V = Λ V

has only the trivial solution for some (and then all)

Λ \in C_{\pm}

. To simplify our notation, we shall denote

\begin{matrix} ζ_{m, n} (Λ) = ζ_{r_{m, n}} (Λ) = (r_{m, n} - Λ)^{1 / 2}, φ_{m, n}^{\pm} (x; Λ) = φ_{ζ_{m, n} (Λ)}^{\pm} (x) . sa.4 \end{matrix}

(4.1)

V \sim {v_{m, n}}

, then each function

v_{m, n}

can be written as

\begin{matrix} v_{m, n} (x) = r_{m, n}^{1 / 4} (C_{m, n}^{+} φ_{m, n}^{+} (x; Λ) + C_{m, n}^{-} φ_{m, n}^{-} (x; Λ)), sa.4y \end{matrix}

(4.2)

with some coefficients

C_{m, n}^{\pm}

. We have inserted the factor

r_{m, n}^{1 / 4}

in order that (cf. 3.13 )

\begin{matrix} {C_{m, n}^{+}, C_{m, n}^{-}} \in ℓ^{2} ⟺ {V \in ℌ} . sa.4a \end{matrix}

(4.3)

The matching conditions 2.6 at

x = \pm 1

yield an infinite system of homogeneous linear equations for the unknown coefficients

C_{m, n}^{\pm}

. Below we set

μ_{\pm} = \sqrt{2} / α_{\pm}

Taking 3.4 into account, we get from the condition at

x = 1

\begin{matrix} r_{m + 1, n}^{1 / 4} (m + 1)^{1 / 2} C_{m + 1, n}^{+} + \frac{2 μ_{+} ζ_{m, n} (Λ) r_{m, n}^{1 / 4}}{1 - e^{- 4 ζ_{m, n} (Λ)}} (C_{m, n}^{+} - C_{m, n}^{-} e^{- 2 ζ_{m, n} (Λ)}) + r_{m - 1, n}^{1 / 4} m^{1 / 2} C_{m - 1, n}^{+} = 0 . sa.4b \end{matrix}

(4.4)

It is convenient to multiply each equation by the factor

r_{m, n}^{1 / 4}

. Let us also denote

\begin{matrix} \begin{matrix} q_{m, n}^{+} = m^{1 / 2} r_{m, n}^{1 / 4} r_{m - 1, n}^{1 / 4}, & q_{m, n}^{-} = n^{1 / 2} r_{m, n}^{1 / 4} r_{m, n - 1}^{1 / 4}; & p_{m, n} (Λ) & = ζ_{m, n} (Λ) r_{m, n}^{1 / 2} . \end{matrix} sa.5 \end{matrix}

(4.5)

The equation 4.4 and the similar equation coming from the condition 2.6 at

x = - 1

yield

\begin{matrix} \begin{matrix} q_{m + 1, n}^{+} C_{m + 1, n}^{+} + \frac{2 μ_{+} p_{m, n} (Λ)}{1 - e^{- 4 ζ_{m, n} (Λ)}} (C_{m, n}^{+} - C_{m, n}^{-} e^{- 2 ζ_{m, n} (Λ)}) + q_{m, n}^{+} C_{m - 1, n}^{+} = 0, & q_{m, n + 1}^{-} C_{m, n + 1}^{-} + \frac{2 μ_{-} p_{m, n} (Λ)}{1 - e^{- 4 ζ_{m, n} (Λ)}} (C_{m, n}^{-} - C_{m, n}^{+} e^{- 2 ζ_{m, n} (Λ)}) + q_{m, n}^{-} C_{m, n - 1}^{-} = 0 . \end{matrix} sa.6 \end{matrix}

(4.6)

Denote by

ℛ = ℛ (Λ)

the infinite matrix which corresponds to this system.

In view of 4.3 , we consider

ℛ

as an operator in the space

G = ℓ^{2} (N_{0}^{2}; C^{2}) .

Removing in 4.6 the exponentially small terms, we come to a simpler system

\begin{matrix} q_{m + 1, n}^{+} C_{m + 1, n}^{+} + 2 μ_{+} p_{m, n} (Λ) C_{m, n}^{+} + q_{m, n}^{+} C_{m - 1, n}^{+} = 0; sa.7 \end{matrix}

(4.7)

\begin{matrix} q_{m, n + 1}^{-} C_{m + 1, n}^{-} + 2 μ_{-} p_{m, n} (Λ) C_{m, n}^{-} + q_{m, n}^{-} C_{m, n - 1}^{-} = 0 . sa.8 \end{matrix}

(4.8)

Let

ℛ^{'} = ℛ^{'} (Λ)

stand for the matrix which corresponds to the system 4.7 – 4.8 , and also for the operator in

G

generated by this matrix. The operator

ℛ^{'}

decomposes into an infinite family of simpler operators. First of all, the equations 4.7 (for

C_{m, n}^{+}

) and 4.8 (for

C_{m, n}^{-}

) are mutually independent. Further, fix any

n \in N_{0}

. The equations in 4.7 which correspond to the chosen value of

n

form a linear system in

ℓ^{2} (N_{0})

with the Jacobi matrix

J_{n}^{+} (Λ) = J ({2 μ_{+} p_{m, n} (Λ)}, {q_{m, n}^{+}}) .

In the same way, the equations in 4.8 , which correspond to the chosen value of

m

, form a linear system in

ℓ^{2} (N_{0})

with the Jacobi matrix

J_{m}^{-} (Λ) = J ({2 μ_{-} p_{m, n} (Λ)}, {q_{m, n}^{-}}) .

The above reasoning shows that

\begin{matrix} ℛ^{'} (Λ) = {\sum_{n}}^{\oplus} J_{n}^{+} (Λ) \oplus {\sum_{m}}^{\oplus} J_{m}^{-} (Λ) . sa.11 \end{matrix}

(4.9)

The original operator

ℛ

can be written as

\begin{matrix} ℛ (Λ) = ℛ^{'} (Λ) + N (Λ) sa.11p \end{matrix}

(4.10)

where

N = N (Λ)

is a block-diagonal matrix with

2 \times 2

-blocks:

\begin{matrix} N_{m, n} (Λ) = \frac{2 p_{m, n} (Λ) e^{- 2 ζ_{m, n} (Λ)}}{1 - e^{- 4 ζ_{m, n} (Λ)}} (\begin{matrix} μ_{+} e^{- 2 ζ_{m, n} (Λ)} & - μ_{+} \\ - μ_{-} & e^{- 2 ζ_{m, n} (Λ)} \end{matrix}) . sa.11n \end{matrix}

(4.11)

The last two equations elucidate the structure of the matrix

ℛ (Λ)

In the rest of the section we take

Λ = i τ \in i R

. We will show that each term on the right-hand side of 4.9 is an invertible operator in

ℓ^{2}

and that the norms of

∥ J_{k}^{\pm} (i τ)^{- 1} ∥

are uniformly bounded. For this purpose, we note that

(2 μ_{+})^{- 1} Im J_{n}^{+} (i τ) = (2 μ_{-})^{- 1} Im J_{m}^{-} (i τ) = d i a g {p_{m, n} (i τ)} .

We have

p_{m, n} (i τ) = \sqrt{r_{m, n}^{2} - i r_{m, n} τ} = X + i Y

where

2 Y^{2} = (r_{m, n}^{4} + r_{m, n}^{2} τ^{2})^{1 / 2} - r_{m, n}^{2} = \frac{r_{m, n} τ^{2}}{(r_{m, n}^{2} + τ^{2})^{1 / 2} + r_{m, n}} \geq 2 c^{2} | τ | .

The last inequality, with some constant

c > 0

, is valid for

| τ | \geq τ_{0}

and for any

m, n \geq 0

; we have taken into account that

r_{m, n} \geq r_{0, 0} = (ν_{+}^{2} + ν_{-}^{2}) / 2

. By the well known estimate for the operators with sign-defined imaginary part, see e.g. Theorem IV.4.1 in [2] , this implies that

∥ J_{k}^{\pm} (i τ)^{- 1} ∥ \leq {(c \sqrt{| τ |})}^{- 1}, \forall k \in N_{0},

and therefore

\begin{matrix} ∥ ℛ^{'} (τ)^{- 1} ∥ = sup ∥ J_{k}^{\pm} (τ)^{- 1} ∥ \leq {(c \sqrt{| τ |})}^{- 1} . sa.13 \end{matrix}

(4.12)

The norms of the blocks

N_{m, n} (i τ)

in 4.11 are controlled by

| p_{m, n} (i τ) | e^{- 2 γ_{m, n} (i τ)}

and hence, are bounded uniformly in

m, n \in N_{0}

. Therefore,

∥ ℛ (\pm i τ) - ℛ^{'} (\pm i τ) ∥ \leq C = C (τ_{0}) .

Choosing

| τ |

large enough, we conclude from 4.12 that

∥ ℛ (\pm i τ) - ℛ^{'} (\pm i τ) ∥ < ∥ ℛ^{'} (\pm i τ)^{- 1} ∥^{- 1} .

It follows that the operator

ℛ (\pm i τ)

has bounded inverse, and, in particular, the system 4.6 has only the trivial solution in

G

The proof of Theorem 2.3 is complete.

5 Representation of the resolvent $(A_{α} - Λ)^{- 1}$

res

5.1 Resolvent.

r In order to prove Theorem 2.8 , we need a convenient representation for both resolvents involved in 2.11 . Here we do this for the operator

A_{α}

. We derive an analogue of the formula (6.6) (the basic formula) in [5] or, equivalently, (6.4) in [10] . However, there is an important difference between the techniques we employ here and those in [5] and [10] . The main goal in both cited papers was the direct study of the a.c. spectrum of the operator corresponding to the case of one oscillator. To achieve this objective, the behaviour of the resolvent as the spectral parameter approaches the real axis was carefully studied. What we do here is to apply scattering theory, and use the already known results of [5] and [10] . This makes our analysis much easier. We are able to work with the resolvents for a fixed value of the parameter

Λ

. We exclude

Λ

from the notation, unless to do so would be confusing.

Let a function

F \in ℌ

have the decomposition

F \sim {f_{m, n}}

. For any

α \geq 0

let us denote

\begin{matrix} U_{α} = {(A_{α} - Λ)}^{- 1} F \sim {u_{α; m, n}}, V \sim {v_{m, n}} = {u_{α; m, n} - u_{0; m, n}} . au0 \end{matrix}

(5.1)

In the notation for

v

we do not reflect dependence on the parameter

α

The operator

(A_{0} - Λ)^{- 1}

can be written in an explicit form. The functions

u_{0; m, n}

are given by the formula 3.9 , with

ζ = ζ_{m, n} (Λ)

, see 4.1 . It follows that

\begin{matrix} u_{0; m, n} (\pm 1) = (2 ζ_{m, n})^{- 1} \int_{R} e^{- ζ_{m, n} | t \mp 1 |} f_{m, n} (t) d t; [u_{0; m, n}^{'}] (\pm 1) = 0 . res.1amn \end{matrix}

(5.2)

Each function

v_{m, n} (x)

belongs to the space

ℱ_{ζ_{m, n}}

. Hence, the equalities 3.4 are satisfied for it, again with

ζ = ζ_{m, n} (Λ)

. Using these equalities and taking into account that

[u_{α; m, n}^{'}] (\pm 1) = [v_{m, n}^{'}] (\pm 1)

, we find from the matching condition in 2.6 at the point

x = 1

that

\begin{matrix} \begin{matrix} - \frac{2 ζ_{m, n}}{1 - e^{- 4 ζ_{m, n}}} & (v_{m, n} (1) - e^{- 2 ζ_{m, n}} v_{m, n} (- 1)) & = \frac{α_{+}}{\sqrt{2}} & (\sqrt{m + 1} u_{α; m + 1, n} (1) + \sqrt{m} u_{α; m - 1, n} (1)) . \end{matrix} res.3 \end{matrix}

(5.3)

As in section 4 , we let

μ_{\pm} = \sqrt{2} / α_{\pm}

. Since

v_{m, n} = u_{α; m, n} - u_{0; m, n}

, the equation 5.3 yields

\begin{matrix} \begin{matrix} \sqrt{m + 1} u_{α; m + 1, n} (1) + \frac{2 μ_{+} ζ_{m, n}}{1 - e^{- 4 ζ_{m, n}}} (u_{α; m, n} (1) - e^{- 2 ζ_{m, n}} u_{α; m, n} (- 1)) & + \sqrt{m} u_{α; m - 1, n} (1) = \frac{2 μ_{+} ζ_{m, n}}{1 - e^{- 4 ζ_{m, n}}} (u_{0; m, n} (1) - e^{- 2 ζ_{m, n}} u_{0; m, n} (- 1)) . \end{matrix} res.3q \end{matrix}

(5.4)

The next step is the same normalization as in section 4 . Denote

\begin{matrix} X_{m, n}^{\pm} = r_{m, n}^{- 1 / 4} u_{0; m, n} (\pm 1), Z_{m, n}^{\pm} = r_{m, n}^{- 1 / 4} u_{α; m, n} (\pm 1) . xz \end{matrix}

(5.5)

Each function

v_{m, n}

can be represented as in 4.2 , with

C_{m, n}^{\pm} = Z_{m, n}^{\pm} - X_{m, n}^{\pm}

We use a shortened notation for the corresponding elements in

G

X = {X_{m, n}^{+}, X_{m, n}^{-}}, Z = {Z_{m, n}^{+}, Z_{m, n}^{-}}, C = {C_{m, n}^{+}, C_{m, n}^{-}}, m, n \in N_{0} .

Multiplying each equation in 5.4 by

r_{m, n}^{1 / 4}

and writing out the similar equations coming from the matching conditions at

x = - 1

, we reduce the system to the form

\begin{matrix} \begin{matrix} q_{m + 1, n}^{+} Z_{m + 1, n}^{+} + \frac{2 μ_{+} p_{m, n}}{1 - e^{- 4 ζ_{m, n}}} (Z_{m, n}^{+} - e^{- 2 ζ_{m, n}} Z_{m, n}^{-}) + q_{m, n}^{+} Z_{m - 1, n}^{+} & = \frac{2 μ_{+} p_{m, n}}{1 - e^{- 4 ζ_{m, n}}} (X_{m, n}^{+} - e^{- 2 ζ_{m, n}} X_{m, n}^{-}); & q_{m, n + 1}^{-} Z_{m, n + 1}^{-} + \frac{2 μ_{-} p_{m, n}}{1 - e^{- 4 ζ_{m, n}}} (Z_{m, n}^{-} - e^{- 2 ζ_{m, n}} Z_{m, n}^{+}) + q_{m, n}^{-} Z_{m - 1, n}^{-} & = \frac{2 μ_{-} p_{m, n}}{1 - e^{- 4 ζ_{m, n}}} (X_{m, n}^{-} - e^{- 2 ζ_{m, n}} X_{m, n}^{+}) . \end{matrix} res.3r \end{matrix}

(5.6)

This is the non-homogeneous counterpart of the system 4.6 . In order to write it more conveniently, we need more notation. All the operators introduced below depend on

Λ

and we always assume that

Λ / \in R

Define the operator

S = S (Λ) : ℌ \to G

S : F \mapsto {\frac{r_{m, n}^{1 / 4}}{2} \int_{R} e^{- ζ_{m, n} | t - 1 |} f_{m, n} (t) d t, \frac{r_{m, n}^{1 / 4}}{2} \int_{R} e^{- ζ_{m, n} | t + 1 |} f_{m, n} (t) d t} .

According to 5.2 and 5.5 , this can be written as

\begin{matrix} S : F \mapsto {p_{m, n} X_{m, n}^{+}, p_{m, n} X_{m, n}^{-}} . cs \end{matrix}

(5.7)

It follows from the Cauchy-Schwartz inequality and 3.13 that the operator

S

is bounded.

The diagonal operator

P = P (Λ) : {X_{m, n}^{+}, X_{m, n}^{-}} \mapsto {p_{m, n} X_{m, n}^{+}, p_{m, n} X_{m, n}^{-}}

acts in

G

and is unbounded. Its inverse

P^{- 1}

is a bounded operator.

Further, let

ℳ = ℳ (Λ)

be the operator generated by the block-diagonal matrix,

ℳ = d i a g {M_{m, n}}

, where

\begin{matrix} M_{m, n} = (1 - e^{- 4 ζ_{m, n}})^{- 1} (\begin{matrix} μ_{+} & - μ_{+} e^{- 2 ζ_{m, n}} \\ - μ_{-} e^{- 2 ζ_{m, n}} & μ_{-} \end{matrix}) . 3e \end{matrix}

(5.8)

Evidently,

ℳ

is bounded in

G

Finally, we let

\begin{matrix} T = T (Λ) : {C_{m, n}^{+}, C_{m, n}^{-}} \mapsto {r_{m, n}^{1 / 4} (C_{m, n}^{+} φ_{m, n}^{+} + C_{m, n}^{-} φ_{m, n}^{-})} . t \end{matrix}

(5.9)

This is a bounded operator acting from

G

into

ℌ

. The system 5.6 can be written in the operator form

ℛ Z = 2 ℳ S F,

whence

Z = 2 ℛ^{- 1} ℳ S F

. Here

ℛ = ℛ (Λ)

is the operator in

G

which corresponds to the left-hand side of the system 5.6 , or, equivalently, of the homogeneous system 4.6 . We also have

X = P^{- 1} S F

, so that

C = Z - X = (2 ℛ^{- 1} ℳ - P^{- 1}) S F .

C

is found from this equation, then evidently

U_{α} - U_{0} = T C

Now it follows from the construction that

\begin{matrix} (A_{α} - Λ)^{- 1} - (A_{0} - Λ)^{- 1} = T (2 ℛ^{- 1} ℳ - P^{- 1}) S . 3f \end{matrix}

(5.10)

This is the desired representation of the resolvent of the operator

A_{α}

5.2 On the matrix $ℛ (Λ)$ .

matr It was shown in section 4 that for

τ

large enough the matrix

ℛ (i τ)

has a bounded inverse. This allowed us to conclude that the operator

A_{α} - Λ

has a bounded inverse for all

Λ / \in R

, and hence

ker ℛ (Λ) = {0}

for all such

Λ

. So, the operator

ℛ (Λ)^{- 1}

is well-defined. However, this does not imply automatically that this operator is bounded in

G

. We now show that this property is a direct consequence of the representation 5.10 . Indeed, 5.10 implies that

2 ℛ^{- 1} ℳ - P^{- 1} = (T^{*} T)^{- 1} T^{*} ((A_{α} - Λ)^{- 1} - (A_{0} - Λ)^{- 1}) S^{*} (S S^{*})^{- 1} .

It is easy to show that for

Λ / \in R_{+}

the operators

ℳ

and

T^{*} T

and

S S^{*}

(acting in

G

) have bounded inverses. This yields the desired result.

6 Representation of the resolvent $(A_{α}^{\circ} - Λ)^{- 1}$

reso Our aim here is to derive an analogue of the representation 5.10 for the operator

A_{α}^{\circ}

. One possible way to proceed is to use the decompositions 2.8 , 2.9 . However, we prefer another way, one which parallels our argument in section 5 . The calculations are easier for

A_{α}^{\circ}

than for

A_{α}

For the objects related to the operator

A_{α}^{\circ}

we use the notation

\begin{matrix} U_{α}^{\circ} = {(A_{α}^{\circ} - Λ)}^{- 1} F \sim {u_{α; m, n}^{\circ}}, V^{\circ} \sim {v_{m, n}^{\circ}} = {u_{α; m, n}^{\circ} - u_{0; m, n}^{\circ}} . au00 \end{matrix}

(6.1)

An analogue of 5.2 is given by

\begin{matrix} u_{0; m, n}^{\circ} (\pm 1) = (2 ζ_{m, n})^{- 1} \int_{R_{\pm}} (e^{- ζ_{m, n} | t \mp 1 |} - e^{- ζ_{m, n} (1 \pm t)}) f_{m, n} (t) d t . u1a \end{matrix}

(6.2)

Next, we derive an analogue of 5.3 . Taking 3.7 and

[{u_{0; m, n}^{\circ}}^{'}] (\pm 1) = 0

into account, we get from the matching condition at

x = 1

- \frac{2 ζ_{m, n}}{1 - e^{- 2 ζ_{m, n}}} v_{m, n}^{\circ} (1) = \frac{α_{+}}{\sqrt{2}} (\sqrt{m + 1} u_{α; m + 1, n}^{\circ} (1) + \sqrt{m} u_{α; m - 1, n}^{\circ} (1)) .

Since

v_{m, n}^{\circ} = u_{α; m, n}^{\circ} - u_{0; m, n}^{\circ}

, we find, taking, as before,

μ_{\pm} = \sqrt{2} / α_{\pm}

\begin{matrix} \sqrt{m + 1} u_{α; m + 1, n}^{\circ} (1) & + & \frac{2 μ_{+} ζ_{m, n}}{1 - e^{- 2 ζ_{m, n}}} u_{α; m, n}^{\circ} (1) + \sqrt{m} u_{α; m - 1, n}^{\circ} (1) \end{matrix}

\begin{matrix} = & \frac{2 μ_{+} ζ_{m, n}}{1 - e^{- 2 ζ_{m, n}}} u_{0; m, n}^{\circ} (1) . \end{matrix}

(6.3)

This is much simpler than the system 5.4 , which, of course, merely reflects the special structure of the operator

A_{α}^{\circ}

as given by 2.9 .

The normalization, as in section 5 , reduces 6.3 and the similar equations for

x = - 1

to the form

\begin{matrix} \begin{matrix} q_{m + 1, n}^{+} Z_{m + 1, n}^{\circ, +} + \frac{2 μ_{+} p_{m, n}}{1 - e^{- 2 ζ_{m, n}}} Z_{m, n}^{\circ, +} + q_{m - 1, n}^{+} Z_{m - 1, n}^{\circ, +} = \frac{2 μ_{+} p_{m, n}}{1 - e^{- 2 ζ_{m, n}}} X_{m, n}^{\circ, +}, & q_{m, n + 1}^{-} Z_{m, n + 1}^{\circ, -} + \frac{2 μ_{-} p_{m, n}}{1 - e^{- 2 ζ_{m, n}}} Z_{m, n}^{\circ, -} + q_{m, n + 1}^{-} Z_{m, n - 1}^{\circ, -} = \frac{2 μ_{-} p_{m, n}}{1 - e^{- 2 ζ_{m, n}}} X_{m, n}^{\circ, -} . \end{matrix} 3m \end{matrix}

(6.4)

Here

X_{m, n}^{\circ, \pm} = r_{m, n}^{- 1 / 4} u_{0; m, n}^{\circ} (\pm 1), Z_{m, n}^{\circ, \pm} = r_{m, n}^{- 1 / 4} u_{α; m, n}^{\circ} (\pm 1) .

The coefficients

q_{m, n}^{\pm}, p_{m, n}

are the same as in 5.6 , being defined in 4.5 . By 6.2 , we have

\begin{matrix} 2 p_{m, n} X_{m, n}^{\circ, +} & = \int_{R_{+}} (e^{- ζ_{m, n} | t - 1 |} - e^{- ζ_{m, n} (t + 1)}) f_{m, n} (t) d t, & 2 p_{m, n} X_{m, n}^{\circ, -} & = \int_{R_{-}} (e^{- ζ_{m, n} | t + 1 |} - e^{ζ_{m, n} (t - 1)}) f_{m, n} (t) d t . \end{matrix}

Now we define analogues of the operators involved in the equality 5.6 . First of all,

ℛ^{\circ} = ℛ^{\circ} (Λ)

is the operator in

G

, defined by the infinite matrix which corresponds to the left-hand side of 6.4 . The operator

ℛ^{\circ}

can be written in the form similar to 4.10 :

\begin{matrix} ℛ^{\circ} (Λ) = ℛ^{'} (Λ) + N^{\circ} (Λ) ro \end{matrix}

(6.5)

where

\begin{matrix} N^{\circ} (Λ) = d i a g {N_{m, n}^{\circ}}, N_{m, n}^{\circ} = 2 \frac{p_{m, n} e^{- 2 ζ_{m, n}}}{1 - e^{- 2 ζ_{m, n}}} (\begin{matrix} μ_{+} & 0 \\ 0 & μ_{-} \end{matrix}) . ro0 \end{matrix}

(6.6)

The self-adjointness of the operator

A_{α}^{\circ}

ℌ

implies that the operator

ℛ^{\circ} (Λ)

is invertible for any

Λ / \in R

The operator

S^{\circ} = S^{\circ} (Λ) : ℌ \to G

is a bounded operator defined by

S^{\circ} : F \sim {f_{m, n}} \mapsto {p_{m, n} X_{m, n}^{\circ, +}, p_{m, n} X_{m, n}^{\circ, -}};

cf. 5.7 . The operator

ℳ^{\circ} = ℳ^{\circ} (Λ)

is the bounded operator on

G

of block-diagonal form

\begin{matrix} ℳ^{\circ} = d i a g {M_{m, n}}, M_{m, n} = (1 - e^{- 2 ζ_{m, n}})^{- 1} (\begin{matrix} μ_{+} & 0 \\ 0 & μ_{-} \end{matrix}) . 3ec \end{matrix}

(6.7)

Finally, let

T^{\circ} = T^{\circ} (Λ) : {C_{m, n}^{+}, C_{m, n}^{-}} \mapsto {r_{m, n}^{1 / 4} (C_{m, n}^{+} φ_{m, n}^{\circ, +} + C_{m, n}^{+} φ_{m, n}^{\circ, -})}

This is a bounded operator acting from

G

into

ℌ

As in section 5 , we can re-write the system 6.4 as

\begin{matrix} (A_{α}^{\circ} - Λ)^{- 1} - (A_{0}^{\circ} - Λ)^{- 1} = T^{\circ} (2 (ℛ^{\circ})^{- 1} ℳ^{\circ} - P^{- 1}) S^{\circ} . 3fc \end{matrix}

(6.8)

Note that for any

Λ \neq \bar{Λ}

the operator

ℛ^{\circ} (Λ)^{- 1}

is bounded. The proof is the same as for the operator

ℛ^{- 1} (Λ)

, see section 5.2 .

7 Proof of Theorem 2.8 .

proof

7.1 The case $α = 0$ .

00 Here we show that

\begin{matrix} (A_{0}^{\circ} - Λ)^{- 3} - (A_{0} - Λ)^{- 3} \in S_{1} . goal0 \end{matrix}

(7.1)

Recall that in the notation of 5.1 and 6.1

(A_{0} - Λ)^{- 1} F = {u_{0, m, n}}, (A_{0}^{\circ} - Λ)^{- 1} F = {u_{0, m, n}^{\circ}}

where the functions

u_{0, m, n}, u_{0, m, n}^{\circ}

are given by the equations 3.9 and 3.10 respectively, with

ζ = ζ_{m, n}

. So, both operators are diagonal. Denote by

Φ_{m, n}, Φ_{m, n}^{\circ}

their components, and let

\begin{matrix} \begin{matrix} Q = d i a g {Q_{m, n}} = & (A_{0}^{\circ} - Λ)^{- 1} - (A_{0} - Λ)^{- 1}, & Q_{m, n} = & Φ_{m, n}^{\circ} - Φ_{m, n} . \end{matrix} pr0.qq \end{matrix}

(7.2)

According to 3.11 , each

Q_{m, n}

is a rank one operator:

\begin{matrix} Q_{m, n} : f_{m, n} (x) \mapsto - (2 ζ_{m, n})^{- 1} g_{m, n} (x) \int_{R} g_{m, n} (t) f_{m, n} (t) d t p2 \end{matrix}

(7.3)

where

g_{m, n} (x) = e^{- ζ_{m, n} | x |}

. It follows from 7.3 that

∥ Q_{m, n} ∥ = (2 | ζ_{m, n} | Re ζ_{m, n})^{- 1} \leq C r_{m, n}^{- 1}, C = C (Λ) .

The norms of

Φ_{m, n}

and

Φ_{m, n}^{\circ}

can be easily estimated (actually,

∥ Φ_{m, n} ∥

can be calculated explicitly, since this is a convolution operator). By the “Schur test”, the norm

∥ K ∥

of an integral operator in

L^{2}

with the kernel

K (x, t)

can be estimated as

∥ K ∥^{2} \leq {sup}_{t} \int | K (x, t) | d x {sup}_{x} \int | K (x, t) | d t .

Applying this to the operators 3.9 and 3.10 , we find that

∥ Φ_{m, n} ∥, ∥ Φ_{m, n}^{\circ} ∥ \leq C r_{m, n}^{- 1} .

Furthermore, the components of the operator

(A_{0}^{\circ} - Λ)^{- 3} - (A_{0} - Λ)^{- 3}

are

Φ_{m, n}^{2} Q_{m, n} + Φ_{m, n} Q_{m, n} Φ_{m, n}^{\circ} + Q_{m, n} (Φ_{m, n}^{\circ})^{2} .

The norm of this operator does not exceed

3 C^{3} r_{m, n}^{- 3}

and, since its rank is no greater than

3

, its trace class norm does not exceed

9 C^{3} r_{m, n}^{- 3}

. By 2.5 , these numbers form a convergent double series, and hence, 7.1 is established.

Note that the exponent

3

in 7.1 can not be replaced by

2

7.2 Difference between the right-hand sides in 5.10 , 6.8 .

e To shorten our notation, let us denote

\begin{matrix} H = T (2 ℛ^{- 1} ℳ - P^{- 1}) S, H^{\circ} = T^{\circ} ((2 ℛ^{\circ})^{- 1} ℳ^{\circ} - P^{- 1}) S^{\circ} . end.0 \end{matrix}

(7.4)

Here we show that

\begin{matrix} Ψ : = H^{\circ} - H = \in S_{1} . end.1 \end{matrix}

(7.5)

Since all the inverse operators appearing in 7.4 are bounded, we only need to check that

T^{\circ} - T, S^{\circ} - S, ℳ^{\circ} - ℳ, ℛ^{\circ} - ℛ \in S_{1} .

Here each operator has a block-diagonal structure, with

(2 \times 2)

-blocks, and it is sufficient to estimate the operator norm of each block and to verify that the corresponding series converge.

For the operators

T^{\circ} - T

and

S^{\circ} - S

the result immediately follows from the definitions of the operators involved and the estimate 3.8 . For the operator

ℳ^{\circ} - ℳ

the result is evident from the comparison of 5.8 and 6.7 . Finally, for

ℛ^{\circ} - ℛ

the result is implied by 4.10 and 6.5 , if we take into account evident estimates of the norms of blocks

N_{m, n}

in 4.11 and

N_{m, n}^{\circ}

in 6.6 .

7.3 End of the proof

end Unfortunately, the desired inclusion 2.11 is not implied by 7.1 and 7.5 automatically and we need an extra argument in order to finalize the proof.

Let us denote

G = (A_{0} - Λ)^{- 1}

. Using also the notation

Q

Ψ

as in 7.2 , 7.5 , we can re-write the equalities 5.10 and 6.8 as

(A_{α} - Λ)^{- 1} = G + H, (A_{α}^{\circ} - Λ)^{- 1} = G + H + Q + Ψ .

We already know that

Ψ \in S_{1}

and

(G + Q)^{3} - G^{3} \in S_{1}

. Therefore, the following equality is satisfied modulo a trace class correction:

(A_{α}^{\circ} - Λ)^{- 1} - (A_{α} - Λ)^{- 1} = (G + Q + H + Ψ)^{3} - (G + H)^{3} = (G + Q + H)^{3} - (G + H)^{3} (m o d S_{1}) = ((G + Q)^{2} - G^{2}) H + (G + Q) H (G + Q) - G H G + H ((G + Q)^{2} - G^{2}) + Q H^{2} + H Q H + H^{2} Q (m o d S_{1}) .

Removing the parentheses, we come to the sum where each term involves one of the products

Q H

H Q

Q G H

H G Q

. Taking into account the structure of the operator

H

, we see that it is sufficient for us to prove that the operators

\begin{matrix} Q T, S Q, Q (A_{0} - Λ)^{- 1} T, S (A_{0} - Λ)^{- 1} Q fin \end{matrix}

(7.6)

are trace class. All these operators have block-diagonal form, with the blocks given by explicit formulae implied by the corresponding definitions. For instance, according to 5.9 the operator

Q T

transforms the number sequence

{C_{m, n}^{+}, C_{m, n}^{-}}

into the sequence of functions

{w_{m, n}}

where

w_{m, n} (x) = - \frac{r_{m, n}^{1 / 4}}{2 ζ_{m, n}} g_{m, n} (x) \int_{R} (C_{m, n}^{+} φ_{m, n}^{+} (t) +, C_{m, n}^{-} φ_{m, n}^{-} (t)) g_{m, n} (t) d t .

An elementary calculation shows that the integral here is of order

O (e^{- γ_{m, n}})

This happens because the function

g_{m, n} (t)

, see 7.3 , is concentrated around the point

t = 0

, while

φ_{m, n}^{\pm} (t)

is concentrated around

t = \pm 1

, and all the three functions decay exponentially when

t

moves away from the corresponding center. Clearly, this estimate implies that

Q T \in S_{1}

. The proofs for the other operators in 7.6 are similar.

The proof of is complete.

8 Acknowledgments

The work on the paper started in the Summer of 2004 when one of the authors (M.S.) was a guest of the School of Mathematics, Cardiff University. M.S. takes this opportunity to express his gratitude to the University for its hospitality and to the EPSRC for financial support under grant GR/T01556. References

W. D. Evans, M. Solomyak, Smilansky's model of irreversible quantum graphs, II: the ES point spectrum, in preparation.
I.C. Gohberg and M.G. Krein, Introduction to the theory of linear non-selfadjoint GK operators in Hilbert space. Izdat. “Nauka”, Moscow 1965. English translation: Amer. Math. Soc., Providence (1969).
T. Kato, Perturbation theory for linear operators, Springer-Verlag, 1966. K
P. Kuchment, Graph models for waves in thin structures, Waves Random Media 12 Ku (2002), no. 4, R1–R24.
S. N. Naboko and M. Solomyak, On the absolutely continuous spectrum in a model NS of irreversible quantum graph, Proc. LMS (in press).
M. Reed and B. Simon, Methods of modern mathematical physics, III: Scattering RS theory, Academic Press, 1979.
U. Smilansky, Irreversible quantum graphs, Waves in Random Media, 14 (2004), 143 SM – 153.
M. Solomyak, On a differential operator appearing in the theory of irreversible S1 quantum graphs, Waves in Random Media, 14(2004), 173-185.
M. Solomyak, On the discrete spectrum of a family of differential operators, Funct. S2 Analysis and its appl., 38 (2004), 217-223.
M. Solomyak, On a mathematical model of the irreversible quantum graph, S3 St.-Petersburg Math. J., submitted.
D. R. Yafaev, Mathematical scattering theory. General theory Translations of Y mathematical monographs, 105. American mathematical Society, Providence, RI, 1993.

School of Mathematics, Cardiff University, 23 Senghennydd Road, Cardiff CF24 4AG, UK E-mail address : EvansWD@cardiff.ac.uk Department of Mathematics, The Weizmann Institute of Science, Rehovot 76100, Israel E-mail address : michail.solomyak@weizmann.ac.il

W.D. Evans

M. Solomyak