Eigenfunction Expansion and the Decomposition of Jacobi Operators on $\mathbb{Z}$

Abstract

Let $J$ be a Jacobi operator on $\ell^{2}\left(\mathbb{Z}\right)$. We prove an eigenfunction expansion theorem for the singular part of $J$ using subordinate solutions to the eigenvalue equation. We exploit this theorem in order to show that $J$ can be decomposed as a direct integral of half-line operators.

1 Introduction

Let $J$ be a Jacobi operator on $\ell^{2}\left(\mathbb{Z}\right)$, by which we mean an operator of the form

$$\left(J\psi\right)\left(n\right)=a_{n-1}\psi\left(n-1\right)+a_{n}\psi\left(n+1\right)+b_{n}\psi\left(n\right),$$ (1.1)

where $\left(a_{n}\right)_{n\in\mathbb{Z}}$ and $\left(b_{n}\right)_{n\in\mathbb{Z}}$ are real-valued and bounded and $\left(a_{n}\right)_{n\in\mathbb{Z}}$ is positive. In this work, we are inerested in studying the connection between the characterization of the singular part of $J$ via the notion of subordinacy, and the eigenfunction expansion of $J$. Throughout, we will study solutions to the eigenvalue equation, i.e. functions $u:\mathbb{Z}\to\mathbb{C}$ which satisfy

$$a_{n-1}u\left(n-1\right)+a_{n}u\left(n+1\right)+b_{n}u\left(n\right)=Eu\left(n\right)$$ (1.2)

for some $E\in\mathbb{R}$ and for every $n\in\mathbb{Z}$. We say that a Borel measure $\mu$ is in the spectral measure class of $J$ if it is equivalent to $P$, the spectral measure of $J$, namely $P\left(A\right)=0\iff\mu\left(A\right)=0$ for any Borel set $A$ (see Definition 2.5).

It is well-known (see, e.g. [1, Chapter V, Theorem 2.1 ]) that for any such $\mu$, there exists a function $N:\mathbb{R}\to\mathbb{N}$ which is defined $\mu$-almost everywhere, and for $\mu$-almost every $E\in\mathbb{R}$ there exists a collection of solutions of (1.2), $\left(u_{E}^{k}\right)_{k=1}^{N\left(E\right)}$, such that for any compactly supported $\psi\in\ell^{2}\left(\mathbb{Z}\right)$ and for any Borel set $A$,

$$P\left(A\right)\psi=\int_{\mathbb{R}}\sum\limits_{k=1}^{N\left(E\right)}\langle u_{E}^{k},\psi\rangle u_{E}^{k}d\mu\left(E\right).$$ (1.3)

A few remarks are in order

1. Both the choice of solutions and the function $N$ are defined $\mu$-almost everywhere.

2. $N$ is called the multiplicity function of $J$.

3. This is a generalization of the finite-dimensional version of the spectral theorem.

4. We are citing here a very partial version of the eigenfunction expansion theorem. For example, (1.3) actually holds for a larger class of sequences. We refer the reader to [1, Chapter V] for more information on eigenfunction expansions.

It is easy to see that for any $E\in\mathbb{R}$, any solution $u$ of (1.2) is determined by $u\left(0\right)$ and $u\left(1\right)$, and vice-versa: for any $\alpha_{0},\alpha_{1}\in\mathbb{C}$, there exists a solution $u$ such that $u\left(0\right)=\alpha_{0},u\left(1\right)=\alpha_{1}$. This implies that the space of solutions of (1.2) is two-dimensional, and so $N\left(E\right)\leq 2$ for any $E\in\mathbb{R}$. Furthermore, if we restrict our attention to the singular part of $J$, then $N\left(E\right)=1$. More formally, in [14] it was shown that $N\left(E\right)=1$ for $\mu_{s}$-almost every $E\in\mathbb{R}$. This was also proved later in [23] using the theory of finite rank perturbations, and in [9, 18] using subordinacy theory, which relates singularity of $\mu$ to asymptotic properties of solutions of (1.2).

Subordinacy theory was first introduced in [10], and was later extended and developed in many ways (for example [7, 8, 9, 12, 13, 16, 18]). Roughly speaking, a solution to (1.2) will be called subordinate if it has slower growth than any other solution at both $\pm\infty$. A formal definition is given in Subection 2.2. Now, a support to the singular part of $\mu$ can be given in terms of subordinacy in the following way:

Theorem 1.1.

([6, Theorem 2.7.11]) Let $\mu_{s}$ be the part of $\mu$ which is singular w.r.t. the Lebesgue measure. Then $\mu_{s}$ is supported on the set

$S=\left\{E\in\mathbb{R}:\text{ there exists a subordinate solution to (\ref{ev_eq_intro})}\right\}$.

For the rest of this section, the measure $\mu$ is taken to be $\mu_{\delta_{0}}+\mu_{\delta_{1}}$, where $\mu_{\psi}$ is the spectral measure of $\psi$ w.r.t. $J$. It is not hard to see that $\mu$ is in the spectral measure class of $J$. It is also not hard to show that if a subordinate solution exists, then it is unique (up to scalar multiplication). However, as mentioned before, the space of solutions is two-dimensional, and so since $N\left(E\right)=1$ for $\mu_{s}$-almost $E\in\mathbb{R}$, it is natural to expect that the solution which is chosen in the RHS of (1.3) is the subordinate one. Given a function $f:\mathbb{Z}\to\mathbb{R}_{>0}$, we denote

$\ell^{2}\left(\mathbb{Z};f\right)\coloneqq\left\{u:\mathbb{Z}\to\mathbb{C}:\sum\limits_{n\in\mathbb{Z}}\frac{\left|u\left(n\right)\right|^{2}}{f\left(n\right)}<\infty\right\}$.

Let us also denote by $P_{s}$ the part of $P$ which is singular w.r.t the Lebesgue measure (we refer the reader to Subsection 2.1 for precise definitions). In this work, we show

Theorem 1.2.

Let $\mu\coloneqq\mu_{\delta_{0}}+\mu_{\delta_{1}}$. Then there exists a choice, for $\mu_{s}$-almost every $E\in\mathbb{R}$, of a subordinate solution $u_{E}$ which satisfies $\left|u_{E}\left(0\right)\right|^{2}+\left|u_{E}\left(1\right)\right|^{2}=1$, such that for any Borel set $B$ and for any $\psi\in\ell^{2}\left(\mathbb{Z};n^{2}\right)$,

$P_{s}\left(B\right)\psi=\int_{B\cap S}\langle u_{E},\psi\rangle u_{E}d\mu_{s}\left(E\right)$.

Remark 1.3.

A certain version of Theorem 1.2 remains valid for a general Borel measure $\rho$ which is in the spectral measure class of $J$. The only modification is in the normalization of the subordinate solutions.

We will exploit Theorem 1.2 in order to obtain a decomposition of $J$ into half-line operators. In general, the decomposition of Jacobi operators on graphs has drawn a lot of attention and proved useful in many cases. A very non-comprehensive list of works in that area includes [2, 3, 4, 5, 17, 19, 24]. Generally, such decompositions are obtained when the underlying operator is the adjacency matrix (namely $a_{n}=1$ and $b_{n}=0$ for all $n\in\mathbb{Z}$), and the underlying graph possesses some strong symmetry properties. In these cases, the operator is decomposed into a direct sum of Jacobi operators on $\mathbb{N}$. In this work, we describe a decomposition for general bounded self-adjoint Jacobi operators on $\mathbb{Z}$ with no symmetry assumptions on the potential. This is the first result of this kind as far as we know. The price that is paid for this generality is that the decomposition is no longer given by a direct sum, but as a certain integral of operators.

Let us denote by $J_{0}^{\pm}$ the restriction of $J$ to $\mathbb{Z}_{\pm}$, where $\mathbb{Z}_{+}=\left\{1,2,\ldots\right\}$ and $\mathbb{Z}_{-}=\mathbb{Z}\setminus\mathbb{Z}_{+}$. For any $\theta\in\left[0,\pi\right)$, we denote

$J_{\theta}\coloneqq\left(J_{0}^{-}-\cot\theta\langle\delta_{0},\cdot\rangle\delta_{0}\right)\oplus\left(J_{0}^{+}-\tan\theta\langle\delta_{1},\cdot\rangle\delta_{1}\right)$.

Note that for any $\theta\in\left[0,\pi\right)$, $J_{\theta}$ is the direct sum of two half-line Jacobi operators. For any self-adjoint operator $T$, we denote by $T_{s}$ the singular part of $T$ (again, we refer the reader to Subsection 2.1 for precise definitions). Our first result is the following decomposition theorem:

Theorem 1.4.

There exists a projection-valued measure $Q$ which depends on $J$ and which is defined on Borel subsets of $\left[0,\pi\right)$, such that

$$\int_{0}^{\pi}J_{\theta}dQ\left(\theta\right)=J_{s}.$$ (1.4)

Theorem 1.4 follows from the following

Theorem 1.5.

The measure $Q$ commutes with $J$, and for any $a,b\in\left(0,\frac{\pi}{2}\right)$ there exists $\gamma>0$ such that for any $\alpha,\beta\in\left(a,b\right)$ and for any $\psi\in\ell^{2}\left(\mathbb{Z};n^{2}\right)$,

$$\|J_{s}Q\left(\left[\alpha,\beta\right)\right)\psi-J_{\theta}Q\left(\left[\alpha,\beta\right)\right)\psi\|<\gamma\left|\beta-\alpha\right|\|Q\left(\left[\alpha,\beta\right)\right)\psi\|$$ (1.5)

Remark 1.6.

1. The decomposition (1.4) is given by the strong limit of the Riemann-Stieltjes sums.

2. Theorems 1.4 and 1.5 actually hold when replacing $J$ and $J_{\theta}$ with $F\left(J\right)$ and $F\left(J_{\theta}\right)$ for certain functions $F$. The stronger versions of these theorems are given in Theorems 4.3 and 4.6.

The rest of the paper is structured as follows. In Section $2$, we present some preliminary background in the one-dimensional theory of Jacobi operators and some basic measure-theoretic background. In Section $3$, we present the eigenfunction expansion of Jacobi operators in $\ell^{2}\left(\mathbb{Z}\right)$ and in particular, prove Theorem 1.2. In Section $4$, we define the measure $Q$ and the collection $\left(J_{\theta}\right)_{\theta\in\left[0,\pi\right)}$ mentioned in Theorems 1.5 and 1.4, and prove these theorems. Section $5$ contains some examples and further discussion.

Acknowledgments I would like to thank Jonathan Breuer for useful discussions.

2 Preliminaries

2.1 The spectral theorem and operator-valued measures

Let $\mathcal{H}$ be a Hilbert space and let $\Theta:\text{Borel}\left(\mathbb{R}\right)\to\mathcal{B}\left(\mathcal{H}\right)$. $\Theta$ will be called a nonnegative operator-valued measure (NOVM), if $\Theta\left(\emptyset\right)=0$, $\Theta\left(B\right)>0$ for any Borel set $B$, and for any collection of disjoint Borel set $\left(B_{j}\right)_{j\in\mathbb{N}}$, one has

$\Theta\left(\underset{n\in\mathbb{N}}{\bigcup}B_{j}\right)=\sum\limits_{j=1}^{\infty}\Theta\left(B_{j}\right)$

where the RHS converges weakly to the LHS. If for any two Borel sets $B_{1},B_{2}$, $\Theta\left(B_{1}\right)$ is an orthogonal projection and in addition, $\Theta\left(B_{1}\cap B_{2}\right)=\Theta\left(B_{1}\right)\Theta\left(B_{2}\right)$, then we will say that $\Theta$ is a projection-valued measure (PVM). We will also be interested in NOVMs which have bounded trace.

Definition 2.1.

We say that a NOVM $\Theta$ has bounded trace if and only if for any bounded Borel set $B\subseteq\mathbb{R}$, for any orthonormal basis of $\mathcal{H}$, $\left(e_{n}\right)_{n=1}^{\infty}$

$\operatorname{tr}\Theta\left(B\right)\coloneqq\sum\limits_{n=1}^{\infty}\langle\Theta\left(B\right)e_{j},e_{j}\rangle<\infty$.

Remark 2.2.

The trace of a nonnegative operator does not depend on the choice of the orthonormal basis.

Throughout this work we will discuss certain limits which are called Riemann-Stieltjes integrals. These are integrals of operator-valued functions w.r.t. either PVMs or standard Borel measures. We present the case of a PVM here. Let us fix a $PVM$ $P$ and a $\mathcal{B}\left(\mathcal{H}\right)$-valued function $F:\left[a,b\right]\to\mathcal{B}\left(\mathcal{H}\right)$. Given a partition of $\left[a,b\right]$ $\Delta=\left(x_{0}=a,x_{1},\ldots,x_{n}=b\right)$ and a selection of points $\Delta^{*}=\left(t_{1},\ldots,t_{n}\right)$ such that $t_{i}\in\left[x_{i-1},x_{i}\right)$, The Riemann Stieltjes sum of $F$ w.r.t. $\Delta$ and $\Delta^{*}$ is defined by

$RS\left(F,\Delta,\Delta^{*}\right)\coloneqq\sum\limits_{k=1}^{n}F\left(t_{j}\right)P\left(\left[x_{j-1},x_{j}\right)\right)$.

We will say that $S\in\mathcal{B}\left(\mathcal{H}\right)$ is the Riemann-Stieltjes integral of $F$ w.r.t. $P$ if the Riemann-Stieltjes sums converge (in operator norm) to $S$ as $\left|\Delta\right|\to 0$, where $\left|\Delta\right|=\underset{i=1,\ldots,n}{\max}\left|x_{i}-x_{i-1}\right|$. If this is the case, we also say that $F$ is integrable on $\left[a,b\right]$ and we denote $S=\int_{a}^{b}F\left(t\right)dP\left(t\right)$.

Remark 2.3.

Given $a,b\in\mathbb{R}$, if $F$ is integrable on any interval $\left[a^{\prime},b^{\prime}\right]\subset\left(a,b\right)$, then we define

$\int_{a}^{b}F\left(t\right)dP\left(t\right)\coloneqq\underset{a^{\prime}\to a,b^{\prime}\to b}{\operatorname{s-\lim}}\int_{a^{\prime}}^{b^{\prime}}F\left(t\right)dP\left(t\right)$

if the limit on the RHS exists.

The following is the well known spectral theorem (see, e.g. [21, Theorem VII.8]).

Theorem 2.4.

Let $T:\mathcal{H}\to\mathcal{H}$ be a self-adjoint operator. Then there exists a projection-valued measure, $P_{T}$, such that for any continuous function $f:\sigma\left(T\right)\to\mathbb{C}$,

$$\int_{\sigma\left(T\right)}f\left(\lambda\right)dP_{T}\left(\lambda\right)=f\left(T\right),$$ (2.1)

where the RHS is taken to be the limit, in operator norm, of the Riemann-Stieltjes sums. $P_{T}$ is called the spectral measure of $T$.

Given $\varphi,\psi\in\mathcal{H}$, their joint spectral measure w.r.t. $T$ is a Borel measure which is given by

$\mu_{\varphi,\psi}\left(B\right)=\langle P_{T}\left(B\right)\varphi,\psi\rangle$

for any Borel set $B$. $\mu_{\varphi,\psi}^{T}$. We will write $\mu_{\varphi}^{T}$ when $\varphi=\psi$. and we will omit the $T$ when it is clear from context. Let us now define the singular and absolutely continuous parts of a self-adjoint operator and of its spectral measure. To that end, let

$\mathcal{H}_{\text{s}}=\left\{u\in\mathcal{H}:\mu_{u}\text{ is singular w.r.t.\ the Lebesgue measure}\right\}$,

$\mathcal{H}_{\text{ac}}=\left\{u\in\mathcal{H}:\mu_{u}\text{ is absolutely continuous w.r.t.\ the Lebesgue measure}\right\}$.

It is well known (see, e.g. [21, Theorem VII.4]) that $\mathcal{H}=\mathcal{H}_{\text{s}}\oplus\mathcal{H}_{\text{ac}}$, and that each of these subspaces is invariant under $T$. We will denote by $T_{\bullet}$ the operator $T|_{\mathcal{H}_{\bullet}}$ and by $\left(P_{T}\right)_{\bullet}$ its spectral measure, where $\bullet\in\left\{\text{s},\text{ac}\right\}$. We will omit the $T$ when it is clear from context. We will also be interested in scalar spectral measures of $T$, which are measures that are in the spectral measure class of $T$.

Definition 2.5.

Let $\Theta$ be a NOVM and let $\rho$ be a Borel measure. We will say that $\rho$ is in the measure class of $\Theta$ if for any Borel set $B$, $\rho\left(B\right)=0\iff\Theta\left(B\right)=0$. Given a self-adjoint operator $T$, we will say that $\rho$ is in the spectral measure class of $T$ if it is in the measure class of $P_{T}$.

We will also need the notion of a set which supports a measure. We define it here for scalar Borel measures. The definition for NOVMs is similar.

Definition 2.6.

Let $\mu$ be a Borel measure. We will say that a Borel set $B\subseteq\mathbb{R}$ supports $\mu$ (or that $\mu$ is supported by $B$) if $\mu\left(\mathbb{R}\setminus B\right)=0$.

Lemma 2.7.

Let $T:\mathcal{H}\to\mathcal{H}$ be a self-adjoint operator and let $\mu$ be a Borel measure which is in the spectral measure class of $T$. Let $A\subseteq\mathbb{R}$ be a Borel set of Lebesgue measure zero which supports $\mu_{s}$. Then $T_{s}=TP\left(A\right)$.

Proof.

This follows immediately from the fact that for any $\psi\in\mathcal{H}_{\text{ac}}$,

$\|P\left(A\right)\psi\|^{2}=\langle P\left(A\right)\psi,P\left(A\right)\psi\rangle=\langle P\left(A\right)\psi,\psi\rangle=\int_{A}xd\mu_{\psi}=0$

∎

2.2 Subordinacy theory and the Borel transform

2.2.1 The half-line case

Let $J^{+}$ be a bounded self-adjoint Jacobi operator on $\ell^{2}\left(\mathbb{N}\right)$, namely an operator of the form

$$\left(J^{+}u\right)\left(n\right)=\begin{cases}a_{n-1}u\left(n-1\right)+a_{n}u\left(n+1\right)+b_{n}u\left(n\right)&n>1\\ a_{1}u\left(2\right)+b_{1}u\left(1\right)&n=1\end{cases}$$ (2.2)

for bounded real-valued sequences $\left(a_{n}\right)_{n\in\mathbb{N}}$ and $\left(b_{n}\right)_{n\in\mathbb{N}}$ such that $\left(a_{n}\right)_{n}\in\mathbb{N}$ is positive. For any $\theta\in\left[0,\pi\right)$, define the operator $J_{\theta}^{+}$ by

$J_{\theta}^{+}\coloneqq J^{+}-\tan\theta\langle\delta_{1},\cdot\rangle\delta_{1}$.

It is not hard to verify that for any $\theta\in\left[0,\pi\right)$, $\delta_{1}$ is a cyclic vector for $J^{+}_{\theta}$, namely $\ell^{2}\left(\mathbb{N}\right)=\overline{\operatorname{sp}\left\{\delta_{1},J_{\theta}^{+}\delta_{1},\left(J_{\theta}^{+}\right)^{2}\delta_{1},\ldots\right\}}$. In that case, the spectral measure of $\delta_{1}$, which we will denote by $\mu_{\theta}$, is in the spectral measure class of $J_{\theta}^{+}$. The purpose of subordinacy theory is to find a set $S$ which supports the singular part of $\mu_{\theta}$, and is defined through asymptotic properties of solutions of the eigenvalue equation,

$$a_{n-1}u\left(n-1\right)+a_{n}u\left(n+1\right)+b_{n}u\left(n\right)=Eu\left(n\right)$$ (2.3)

subject to the boundary condition

$$u\left(0\right)\cos\theta+u\left(1\right)\sin\theta=0,$$ (2.4)

by which we mean that $u$ is defined on $\mathbb{N}\cup\left\{0\right\}$, satisfies (2.3) for $n\in\mathbb{N}$, and satisfies (2.4).

Given $u:\mathbb{N}\to\mathbb{C}$ and $L>1$, we define

$\|u\|_{L}\coloneqq\left[\sum\limits_{n=1}^{\left[L\right]}\left|u\left(n\right)\right|^{2}+\left(L-\left[L\right]\left|u\left(\left[L\right]+1\right)\right|^{2}\right)\right]^{\frac{1}{2}}$.

In other words, $\|\cdot\|_{L}$ denotes the norm of a vector $u$ over an interval of length $L$. Given $E\in\mathbb{R}$ and $\theta\in\left[0,\pi\right)$, we will denote by $u_{\theta,E}$ the solution of (2.3) which satisfies (2.4), normalized so that $u_{\theta,E}\left(1\right)=\cos\theta$. With that in mind, we can present the definition of subordinacy.

Definition 2.8.

Given $E\in\mathbb{R}$ and $\theta\in\left[0,\pi\right)$, $u_{\theta,E}$ will be called subordinate if for any $\eta\neq\theta$, we have

$$\underset{L\to\infty}{\lim}\frac{\|u_{\theta,E}\|_{L}}{\|u_{\eta,E}\|_{L}}=0.$$ (2.5)

Remark 2.9.

For any $E\in\mathbb{R}$, the collection $\left\{\alpha u_{\theta,E}:\alpha\in\mathbb{C},\theta\in\left[0,\pi\right)\right\}$ is a two-dimensional vector space. This implies that $u_{\theta,E}$ is subordinate if and only if there exists $\eta\neq\theta$ such that (2.5) holds.

We have the following

Theorem 2.10.

([16]) The part of $\mu_{\theta}$ which is singular w.r.t. the Lebesgue measure is supported on the set

$S_{\theta}\coloneqq\left\{E\in\mathbb{R}:u_{\theta,E}\text{ is subordinate}\right\}$.

Remark 2.11.

Originally, subordinacy theory was proved in [10] for continuum Schrödinger operators which act on a subset of $L^{2}\left(\mathbb{R}_{\geq 0}\right)$. We state here the discrete setting which was proved in [16] as it is more relevant for this work.

2.2.2 The line case

Subordinacy theory was extended in several ways (see, e.g. [7, 8, 12, 18]). In this work, we will use its extension to Jacobi operators on $\ell^{2}\left(\mathbb{Z}\right)$. To that end, let $J:\ell^{2}\left(\mathbb{Z}\right)\to\ell^{2}\left(\mathbb{Z}\right)$ be given by (1.1). $\delta_{1}$ is no longer necessarily a cyclic vector, but now the pair $\left\{\delta_{0},\delta_{1}\right\}$ is cyclic, i.e. $\ell^{2}\left(\mathbb{Z}\right)=\overline{\operatorname{sp}\left\{J^{n}\delta_{j}:n\in\mathbb{N}\cup\left\{0\right\},j\in\left\{0,1\right\}\right\}}$. This implies that $\mu\coloneqq\mu_{\delta_{0}}+\mu_{\delta_{1}}$ is in the spectral measure class of $J$. In the case of $\mathbb{Z}$, a solution $u$ of the equation

$$u\left(n-1\right)+u\left(n+1\right)+V\left(n\right)u\left(n\right)=Eu\left(n\right)$$ (2.6)

for some $E\in\mathbb{R}$ will be called subordinate if and only if its restrictions (under an appropriate normalization) to $\mathbb{Z}_{+}\coloneqq\left\{1,2,\ldots\right\}$ and to $\mathbb{Z}_{-}\coloneqq\left\{0,-1,\ldots\right\}$ are subordinate. We have the following

Theorem 2.12.

The singular part of $\mu$ is supported on the set

$$S\coloneqq\left\{E\in\mathbb{R}:\text{(\ref{ev_eq_line}) has a subordinate solution}\right\}.$$ (2.7)

This theorem was proved for the continuum case in [8]. A proof for the discrete case can be found in [6, 18].

2.2.3 The Borel transform of a measure

Given a Borel measure $\sigma$, its Borel transform is defined by

$M_{\sigma}\left(z\right)=\int_{\mathbb{R}}\frac{d\sigma\left(E\right)}{E-z}$

for any $z\in\mathbb{C}$ such that $\operatorname{Im}z>0$. The boundary behavior of the Borel transform of a measure $\mu$ is related to continuity properties of $\mu$ by the following theorem (see, for example, [22, Theorem 1.6]).

Theorem 2.13.

Let $\mu$ be a positive measure satisfying $\int_{\mathbb{R}}\frac{d\mu\left(x\right)}{\left|x\right|+1}<\infty$. Denote by $\mu_{ac},\mu_{s}$ the absolutely continuous and singular parts of $\mu$ (with respect to the Lebesgue measure) respectively. Then

1. $\mu_{ac}$ is supported on the set $\left\{E\in\mathbb{R}:0<\underset{\epsilon\to 0}{\lim}\operatorname{Im}m_{\mu}(E+i\epsilon)<\infty\right\}$.

2. $\mu_{s}$ is supported on the set $\left\{E\in\mathbb{R}:\underset{\epsilon\to 0}{\lim}\left|m_{\mu}(E+i\epsilon)\right|=\infty\right\}$.

Remark 2.14.

The spectral measure of any vector w.r.t. any bounded self-adjoint operator satisfies the condition required in Theorem 2.13. With that in mind, there is a useful reformulation of Theorem 2.10: Denote by $m_{\theta}$ the Borel transform of $\mu_{\theta}$. Then, using the formula

$$m_{\theta}=\frac{m_{0}}{1-\tan\theta m_{0}},$$ (2.8)

and an inequality which relates the growth rate of solutions to the behavior of the Borel transoform of $m_{0}$ (see [12, Theorem 1.1]), one can show that $\underset{\epsilon\to 0}{\lim}\left|m_{\theta}\left(E+i\epsilon\right)\right|=\infty$ if and only if $u_{\theta,E}$ is subordinate. In other words, we have

$$u_{\theta,E}\text{ is subordinate}\iff\underset{\epsilon\to 0}{\lim}m_{0}\left(E+i\epsilon\right)=\cot\theta.$$ (2.9)

We will use the following theorem of Poltoratskii:

Theorem 2.15.

([11, Theorem 1.1]) Let $\nu,\sigma$ be two Borel measures such that $\sigma$ is positive. Suppose that $\nu\ll\sigma$. Then for $\sigma_{s}$-almost every $E\in\mathbb{R}$, we have

$\underset{\epsilon\to 0}{\lim}\frac{M_{\nu}\left(E+i\epsilon\right)}{M_{\sigma}\left(E+i\epsilon\right)}=\frac{d\nu}{d\sigma}\left(E\right)$,

where $\sigma_{s}$ denotes the part of $\sigma$ which is singular w.r.t. the Lebesgue measure.

Remark 2.16.

This theorem was originally proved in [20] for measures on the unit circle. Although the transition between the two versions is elementary, we cite here the real-valued version, which was proved in [11].

We will also use the following

Lemma 2.17.

Let $T$ be a bounded self-adjoint operator on a Hilbert space $\mathcal{H}$ and let $\varphi_{1},\varphi_{2}\in\mathcal{H}$. Let $\mu=\mu_{\varphi_{1}}+\mu_{\varphi_{2}}$. Then for any $k,j\in\left\{1,2\right\}$, $\mu_{\varphi_{k},\varphi_{j}}\ll\mu$, and for $\mu$-almost every $E\in\mathbb{R}$, $\left|\frac{d\mu_{\varphi_{k}\varphi_{j}}}{d\mu}\left(E\right)\right|\leq 1$.

Proof.

For any Borel set $B\subseteq\mathbb{R}$ and $\psi,\varphi\in\mathcal{H}$, we have that $\mu_{\psi,\varphi}\left(B\right)=\langle P_{T}\left(B\right)\psi,\varphi\rangle$. This implies that

$\mu_{\varphi_{k},\varphi_{j}}^{2}\left(B\right)=\left|\langle P_{T}\left(B\right)\varphi_{k},\varphi_{j}\rangle\right|^{2}\leq\langle P_{T}\left(B\right)\varphi_{k},\varphi_{k}\rangle\cdot\langle P_{T}\left(B\right)\varphi_{j},\varphi_{j}\rangle\leq\mu^{2}\left(B\right)$,

which implies both that $\mu_{\varphi_{k},\varphi_{j}}\ll\mu$ and that $\left|\frac{d\mu_{\varphi_{k}\varphi_{j}}}{d\mu}\right|\leq 1$ $\mu$-almost everywhere, as required. ∎

It will be useful to show that the subordinate solution is actually given in terms of the Radon-Nikodym derivative of certain spectral measures. More precisely, we will use the following

Proposition 2.18.

Let $J$ be a Jacobi operator on $\ell^{2}\left(\mathbb{Z}\right)$, let $\mu\coloneqq\mu_{\delta_{0}}+\mu_{\delta_{1}}$ and denote by $M$ the Borel transform of $\mu$. For any $k,j\in\mathbb{Z}$, let $\mu_{kj}$ be the joint spectral measure of $\delta_{k}$ and $\delta_{j}$ w.r.t. $J$ and let $M_{kj}$ be its Borel transform. There exists a set $\widetilde{S}\subseteq S$ (where $S$ is given by (2.7)) which supports $\mu_{s}$ such that for any $E\in\widetilde{S}$ and for any $k,j\in\mathbb{Z}$, we have

$$\underset{\epsilon\to 0}{\lim}\frac{\mu_{kj}\left(\left(E-\epsilon,E+\epsilon\right)\right)}{\mu\left(E-\epsilon,E+\epsilon\right)}=\underset{\epsilon\to 0}{\lim}\frac{M_{kj}\left(E+i\epsilon\right)}{M\left(E+i\epsilon\right)}=u_{E}\left(k\right)u_{E}\left(j\right),$$ (2.10)

where $u_{E}$ is the subordinate solution of (2.6) which satisfies $\left|u_{E}\left(0\right)\right|^{2}+\left|u_{E}\left(1\right)\right|^{2}=1$.

The proof of Proposition 2.18 mostly follows the proof of [18, Theorem 1.4], and so it is given in the Appendix.

3 Eigenfunction expansion for Jacobi operators on $\ell^{2}\left(\mathbb{Z}\right)$

Let $J$ be a Jacobi operator on $\ell^{2}\left(\mathbb{Z}\right)$. Our goal in this section is to obtain an eigenfunction expansion of the singular part of $J$ following the lines of Berezanskiĭ [1] with slight modifications. Throughout this section, we denote the spectral measure of $J$ by $P$. We will also denote $\mu=\mu_{\delta_{0}}+\mu_{\delta_{1}}$, and $S$ the support of $\mu$ which is given by Theorem 2.12. We will use the following theorem:

Proposition 3.1.

([1, Chapter V, Theorem 1.1]) Let $\Theta$ be a NOVM which has bounded trace, and let $\rho$ be a positive Borel measure which is in the measure class of $\Theta$. Then there exists an operator-valued function $\Psi:\mathbb{R}\to\mathcal{B}\left(\mathcal{H}\right)$ which is defined $\rho$-almost everywhere, such that for any Borel set $B\subseteq\mathbb{R}$,

$$\Theta\left(B\right)=\int_{B}\Psi\left(E\right)d\rho\left(E\right),$$ (3.1)

where the Riemann-Stieltjes integral in (3.1) converges in operator norm. For $\rho$-almost every $E\in\mathbb{R}$, The operator $\Psi$ can be taken to be the weak limit of $\frac{\Theta\left(I_{E,n}\right)}{\rho\left(I_{E,n}\right)}$, where $I_{E,n}=\left(E-\frac{1}{n},E+\frac{1}{n}\right)$. In particular, for $\rho$-almost every $E\in\mathbb{R}$, the weak limit exists.

Remark 3.2.

In [1], the theorem is proven for $\rho=\operatorname{tr}\Theta$. The proof remains exactly the same when one takes $\rho$ to be an arbitrary measure which is in the measure class of $\Theta$.

Our goal is to obtain a formula similar to (3.1) with $P_{s}$ instead of $\Theta$, where $P_{s}$ is the singular part of the spectral measure of $J$. To that end, let us introduce the operator $A:D\left(A\right)\subset\ell^{2}\left(\mathbb{Z}\right)\to\ell^{2}\left(\mathbb{Z}\right)$ which is given by

$\left(A\psi\right)\left(n\right)=n\psi\left(n\right)$

and is defined on the maximal domain

$D\left(A\right)\coloneqq\left\{\psi\in\ell^{2}\left(\mathbb{Z}\right):A\psi\in\ell^{2}\left(\mathbb{Z}\right)\right\}$.

Note that $D\left(A\right)$ contains all of the sequences which are compactly supported. Also note that $A$ has a bounded and self-adjoint inverse.

Claim 3.3.

For any Borel set $B\subseteq\mathbb{R}$, let $\Theta\left(B\right)\coloneqq A^{-1}P\left(B\right)A^{-1}$. Then $\Theta$ is a NOVM with a bounded trace.

Proof.

The fact that $P$ is a PVM along with $A^{-1}$ being a positive operator imply that $\Theta$ is a NOVM. We turn to show that $\Theta$ has bounded trace. Let $B\subseteq\mathbb{R}$ be a Borel set.

$$\begin{split}\sum\limits_{n\in\mathbb{Z}}\langle\Theta\left(B\right)\delta_{n},\delta_{n}\rangle=\sum\limits_{n\in\mathbb{Z}}\langle A^{-1}P\left(B\right)A^{-1}\delta_{n},\delta_{n}\rangle=\sum\limits_{n\in\mathbb{Z}}\frac{1}{n^{2}}\langle P\left(B\right)\delta_{n},\delta_{n}\rangle\leq\sum\limits_{n\in\mathbb{Z}}\frac{1}{n^{2}}<\infty,\end{split}$$

where the last equality follows from the fact that $P\leq\text{Id}$. ∎

By Claim 3.3, $\Theta$ satisfies the conditions of Proposition 3.1. In addition, $\mu$ is in the measure class of $\Theta$ and so (3.1) (with $\rho=\mu$) holds for some operator-valued function $\Psi$ which is defined $\mu$-almost everywhere. Furthermore, $\Psi$ is given by

$\Psi=\underset{n\to\infty}{\operatorname{w-\lim}}\frac{\Theta\left(I_{E,n}\right)}{\mu\left(I_{E,n}\right)}$,

and so for $\mu$-almost every $E\in\mathbb{R}$, we have

$$\langle\Psi\left(E\right)\delta_{k},\delta_{j}\rangle=\underset{n\to\infty}{\lim}\frac{1}{\mu\left(I_{E,n}\right)}\langle T^{-1}\Theta\left(I_{E,n}\right)T^{-1}\delta_{k},\delta_{j}\rangle=\underset{n\to\infty}{\lim}\frac{1}{k^{2}}\frac{\mu_{kj}\left(I_{E,n}\right)}{\mu\left(I_{E,n}\right)}=\frac{1}{k^{2}}\frac{d\mu_{kj}}{d\mu}\left(E\right).$$ (3.2)

By Theorem 2.15 and Proposition 2.18, there exists a Borel set $\widetilde{S}\subseteq S$ such that $\mu_{s}\left(S\setminus\widetilde{S}\right)=0$, and for any $E\in\widetilde{S}$, $k,j\in\mathbb{Z}$, we have

$$\underset{n\to\infty}{\lim}\frac{\mu_{kj}\left(I_{E,n}\right)}{\mu\left(I_{E,n}\right)}=\underset{\epsilon\to 0}{\lim}\frac{M_{kj}\left(E+i\epsilon\right)}{M\left(E+i\epsilon\right)}=u_{E}\left(k\right)u_{E}\left(j\right),$$ (3.3)

where $u_{E}$ is the subordinate solution of (2.6) which satisfies $\left|u_{E}\left(0\right)\right|^{2}+\left|u_{E}\left(1\right)\right|^{2}=1$. We can now prove Theorem 1.2.

Proof.

Define $\mathcal{H}_{\pm}\coloneqq\ell^{2}\left(\mathbb{Z};\left|n\right|^{\pm 2}\right)$. For any $\varphi\in\mathcal{H}$, $A\varphi\in\mathcal{H}_{-}$ and in addition, $\mathcal{H}_{+}=D\left(A\right)$. Thus, for any $E\in\widetilde{S}$, the operator $\Phi\left(E\right)\coloneqq A\Psi\left(E\right)A$ is defined from $\mathcal{H}_{+}$ to $\mathcal{H}_{-}$. Furthermore, for any $u,v\in\mathcal{H}_{+}$, we have

$\langle\Psi\left(E\right)Au,Av\rangle=\langle A\Psi\left(E\right)Au,v\rangle=\langle\Phi\left(E\right)u,v\rangle$

and so for any Borel set $B\subseteq\mathbb{R}$,

$$\begin{split}\langle P\left(B\right)u,v\rangle=\langle P\left(B\right)A^{-1}Au,A^{-1}Av\rangle=\langle A^{-1}P\left(B\right)A^{-1}Au,Av\rangle=\\ =\langle\Theta\left(B\right)Au,Av\rangle=\int_{B}\langle\Psi\left(E\right)Au,Av\rangle d\mu\left(E\right)=\int_{B}\langle\Phi\left(E\right)u,v\rangle d\mu\left(E\right).\end{split}$$ (3.4)

Thus, the integral

$$P\left(B\right)=\int_{B}\Phi\left(E\right)d\mu\left(E\right)$$ (3.5)

converges weakly. Note that $\Phi\left(E\right)\in\mathcal{B}\left(\mathcal{H}_{+},\mathcal{H}_{-}\right)$ and in addition, $\|\Phi\left(E\right)\|\leq 1$. Thus, the weak convergence of (3.5) implies convergence in the operator norm of $\mathcal{B}\left(\mathcal{H}_{+},\mathcal{H}_{-}\right)$. Furthermore, by (3.2) and (3.3), for $\mu_{s}$-almost every $E\in\mathbb{R}$, $\Phi\left(E\right)=\langle u_{E},\cdot\rangle u_{E}$, and so for any Borel set $B$ and $u\in\mathcal{H}_{+}$, we get

$P_{s}\left(B\right)u=\int_{B\cap S}\Phi\left(E\right)ud\mu_{s}\left(E\right)=\int_{B\cap S}\langle u_{E},u\rangle u_{E}d\mu_{s}\left(E\right)$,

as required. ∎

4 $J_{s}$ as an integral of half-line operators

The notation used in the prior section applies throughout this section as well. We will also denote $I\coloneqq\left[0,\pi\right)$.

4.1 The correlative spectral resolution

In this subsection, we will define a projection-valued measure $Q$ on $\left[0,\pi\right)$ which, in some sense, measures the correlation between the restrictions of $J$ to the half-lines $\mathbb{Z}_{\pm}$. Recall that $S\subseteq\mathbb{R}$ is the set of all $E\in\mathbb{R}$ for which there exists a subordinate solution of (2.6), and that $P_{s}$ is supported on the set $\widetilde{S}\subseteq S$ which was defined in the previous section. Denote by $J_{\pm}$ the operators

$J_{\pm}=P_{\mathbb{Z}_{\pm}}JP_{\mathbb{Z}_{\pm}}$

where $\mathbb{Z}_{+}=\left\{1,2,\ldots\right\}$ and $\mathbb{Z}_{-}=\mathbb{Z}_{+}^{c}$. For any $\theta\in\left[0,\pi\right)$, we define $J_{+}^{\theta}=J_{+}-\tan\theta\langle\delta_{1},\cdot\rangle\delta_{1}$, and $J_{-}^{\theta}=J_{-}-\cot\theta\langle\delta_{0},\cdot\rangle\delta_{0}$. Recall that for any $E\in\mathbb{R}$ and $\theta\in\left[0,\pi\right)$, $u_{\theta,E}$ is the solution of (2.3) which satisfies (2.4), and let

$S_{\theta}^{\pm}=\left\{E\in\mathbb{R}:u_{\theta,E}^{\pm}\text{ is subordinate}\right\}$.

Let us denote $S_{\theta}\coloneqq S_{-}^{\theta}\cap S_{+}^{\theta}$. It is not hard to see that $S=\underset{\theta\in I}{\bigcup}S_{\theta}$. In addition, by (3.3) and (2.4), for any $\theta\in\left[0,\pi\right)$ and for any $E\in\widetilde{S}$,

$$E\in S_{\theta}\iff\underset{\epsilon\to 0}{\lim}\frac{M_{0}\left(E+i\epsilon\right)}{M\left(E+i\epsilon\right)}=\frac{d\mu_{00}}{d\mu}\left(E\right)=\cos\theta.$$ (4.1)

Proposition 4.1.

For any Borel set $B\subseteq I$, the set $S_{B}\coloneqq\underset{\theta\in B}{\bigcup}S_{\theta}\cap\widetilde{S}$ is Borel measurable.

Proof.

Denote $U=\cos\left(B\right)$ and let $f=\frac{d\mu_{00}}{d\mu}\left(E\right)$ given by (4.1). Then $S_{B}=f^{-1}\left(U\right)$, which is a measurable set. ∎

Now, for any Borel set $B\subseteq I$, we can define

Q$\left(B\right)=P\left(S_{B}\right)$.

Claim 4.2.

$Q$ is a projection-valued measure which commutes with $J$.

Proof.

This is immediate by the fact that $P$ is a projection-valued measure which commutes with $J$. ∎

One of the main reasons to consider the projection-valued measure $Q$ is that it enables one to split $\ell^{2}\left(\mathbb{Z}\right)$ into a direct sum of subspaces such that on each subspace, $J$ looks like the direct sum of half-line operators. For any $\theta\in\left[0,\pi\right)$, we define $J_{\theta}\coloneqq J^{\theta}_{-}\oplus J^{\theta}_{+}$.

Theorem 4.3.

Let $a,b\in\left(0,\frac{\pi}{2}\right)$ and let $\mathcal{F}=\operatorname{sp}\left(\mathcal{F}_{1}\cup\mathcal{F}_{2}\cup\mathcal{F}_{2}\right)\subseteq C\left(\mathbb{R}\right)$, where

$\mathcal{F}_{1}=\left\{x^{n}:n\in\mathbb{N}\right\}$,

$\mathcal{F}_{2}=\left\{e^{itx}:t\in\mathbb{R}\right\}$

$\mathcal{F}_{3}=\left\{\left(x-z\right)^{-1}:z\in\mathbb{C},\operatorname{Im}z>0\right\}$.

For any $F\in\mathcal{F}$ there exists $\gamma>0$ such that for any $a\leq\alpha<\beta\leq b$, $\psi\in\mathcal{H}_{+}$ and $\theta\in\left[\alpha,\beta\right)$,

$$\|F\left(J_{\theta}\right)Q\left(\left[\alpha,\beta\right)\right)\psi-F\left(J\right)Q\left(\left[\alpha,\beta\right)\right)\psi\|<\gamma\left|\beta-\alpha\right|\|Q\left(\left[\alpha,\beta\right)\right)\psi\|.$$ (4.2)

Proof.

We will prove the theorem for basis elements. The general case then easily follows using the triangle inequality. As before, for any $E\in S$, we denote by $u_{E}$ the subordinate solution of (2.6) which satisfies (2.4). Let $0<M\in\mathbb{R}$ be some constant such that $\|J\|\leq M$ and for any $\theta\in\left[a,b\right]$, $\|J_{\theta}\|\leq M$. We treat each $\mathcal{F}_{i}$ separately.
First case: For $F\in\mathcal{F}_{1}$, we will prove the theorem by induction on $n$. Suppose $F\in\mathcal{F}$ is given by $F\left(x\right)=x$. A straightforward computation shows that for any $\psi\in\ell^{2}\left(\mathbb{Z}\right)$,

$$\left(J_{\theta}-J\right)\psi=\left(\psi\left(1\right)+\cot\theta\psi\left(0\right)\right)\delta_{0}+\left(\psi\left(0\right)+\tan\theta\psi\left(1\right)\right)\delta_{1}.$$ (4.3)

In addition, let $A=\underset{\theta\in\left[\alpha,\beta\right)}{\bigcup}S_{\theta}$. By Theorem 1.2, for any $\psi\in\mathcal{H}_{+}$ and for any $n\in\mathbb{Z}$,

$$\left(Q\left(\left[\alpha,\beta\right)\right)\psi\right)\left(n\right)=\int_{A}\langle\psi,\psi_{E}\rangle\psi_{E}\left(n\right)d\mu\left(E\right)$$ (4.4)

and in addition,

$\|J_{\theta}Q\left(\left[\alpha,\beta\right)\right)\psi-JQ\left(\left[\alpha,\beta\right)\right)\psi\|^{2}=\|\left(J_{\theta}-J\right)Q\left(\left[\alpha,\beta\right)\right)\psi\|^{2}=$

$=\left|\left(J-J_{\theta}\right)Q\left(\left[\alpha,\beta\right)\right)\psi\left(0\right)\right|^{2}+\left|\left(J-J_{\theta}\right)Q\left(\left[\alpha,\beta\right)\right)\psi\left(1\right)\right|^{2}.$

where the last equality follows from (4.3). Recall that $S_{\left[\alpha,\beta\right]}=\underset{\theta\in\left[\alpha,\beta\right]}{S_{\theta}}$ and for any $E\in S_{\left[\alpha,\beta\right]}$, let $\theta\left(E\right)$ be the unique $\theta\in\left[\alpha,\beta\right)$ such that $E\in S_{\theta}$. Now, note that

$\left|\left(J-J_{\theta}\right)Q\left(\left[\alpha,\beta\right]\right)\psi\left(1\right)\right|^{2}=\left|\int_{S_{\left[\alpha,\beta\right]}}\langle\psi,\psi_{E}\rangle\psi_{E}\left(0\right)d\mu\left(E\right)+\tan\theta\int_{S_{\left[\alpha,\beta\right]}}\langle\psi,\psi_{E}\rangle\psi_{E}\left(1\right)d\mu\left(E\right)\right|^{2}=$

$=\left|\int_{S_{\left[\alpha,\beta\right]}}\left(\tan\theta-\tan\theta\left(E\right)\right)\langle\psi,\psi_{E}\rangle\psi_{E}\left(1\right)d\mu\left(E\right)\right|^{2}\leq\underset{\eta\in\left[\alpha,\beta\right]}{\max}\left|\tan\theta-\tan\eta\right|^{2}\left|\int_{S_{\left[\alpha,\beta\right]}}\langle\psi,\psi_{E}\rangle\psi_{E}\left(1\right)d\mu\left(E\right)\right|^{2}=$

$=\underset{\eta\in\left[\alpha,\beta\right]}{\max}\left|\tan\theta-\tan\eta\right|^{2}\|Q\left(\left[\alpha,\beta\right]\right)\psi\|^{2}$.

In addition, $\tan$ satisfies the Lipschitz condition on $\left[a,b\right]$ and so there exists $\gamma>0$ such that $\left|\tan\theta-\tan\eta\right|<\gamma\left|\theta-\eta\right|$. This immediately implies the required inequality. The analysis of the second summand is done in a similar way, noting that $\cot\theta$ also satisfies the Lipschitz condition on $\left[a,b\right]$. Now suppose that the theorem holds for $x^{n-1}$ and let $F\in\mathcal{F}$ be given by $F\left(x\right)=x^{n}$. We have

$\|J_{\theta}^{n}Q\left(\left[\alpha,\beta\right)\right)\psi-J^{n}P\left(\left[\alpha,\beta\right)\right)\psi\|=\|J_{\theta}^{n-1}J_{\theta}Q\left(\left[\alpha,\beta\right)\right)\psi-J^{n-1}JQ\left(\left[\alpha,\beta\right)\right)\psi\|\leq$

$\leq\|J_{\theta}^{n-1}J_{\theta}^{n}Q\left(\left[\alpha,\beta\right)\right)\psi-J_{\theta}^{n-1}JQ\left(\left[\alpha,\beta\right)\right)\psi\|+\|J_{\theta}^{n-1}JQ\left(\left[\alpha,\beta\right)\right)\psi-J^{n-1}JQ\left(\left[\alpha,\beta\right)\right)\psi\|=$

$=\|J_{\theta}^{n-1}\big{[}J_{\theta}Q\left(\left[\alpha,\beta\right)\right)\psi-JQ\left(\left[\alpha,\beta\right)\right)\psi\big{]}\|+\|J_{\theta}^{n-1}Q\left[\alpha,\beta\right)J\psi-J^{n-1}Q\left[\alpha,\beta\right)J\psi\|\leq$

$\leq\|J_{\theta}^{n-1}\|\cdot\gamma_{1}\left|\beta-\alpha\right|\|Q\left(\left[\alpha,\beta\right)\psi\right)\|+\gamma_{2}\left|\beta-\alpha\right|\|Q\left(\left[\alpha,\beta\right)\right)J\psi\|$

where the last inequality follows from the induction hypothesis for suitable constants $\gamma_{1},\gamma_{2}$. Note that $J$ and $Q$ commute, and so we may proceed

$\leq M^{n-1}\cdot\gamma_{1}\left|\beta-\alpha\right|\|Q\left(\left[\alpha,\beta\right)\right)\psi\|+\gamma_{2}\left|\beta-\alpha\right|\|JQ\left(\left[\alpha,\beta\right)\right)\psi\|\leq$

$\leq\left(\underset{=\gamma}{\overset{\left(\star\right)}{\underbrace{M^{n-1}\gamma_{1}+M\gamma_{2}}}}\right)\left|\beta-\alpha\right|\|Q\left(\left[\alpha,\beta\right)\right)\psi\|$

as required. Note also that using ($\star$), if we denote by $\gamma_{n}$ the constant which corresponds with $x^{n}$, then it can easily be shown by induction that there exists $C>0$ such that $\gamma_{n}\leq CM^{n}$.
Second case: Let $t\in\mathbb{R}$, $\psi\in\mathcal{H}_{+}$. We have

$\|e^{itJ_{\theta}}Q\left(\left[\alpha,\beta\right)\right)\psi-e^{itJ}Q\left(\left[\alpha,\beta\right)\right)\psi\|=\|\sum\limits_{n=0}^{\infty}\frac{\left(itJ_{\theta}\right)^{n}}{n!}Q\left(\left[\alpha,\beta\right)\right)\psi-\frac{\left(itJ\right)^{n}}{n!}Q\left(\left[\alpha,\beta\right)\right)\psi\|\leq\sum\limits_{n=0}^{\infty}\frac{\left|t\right|^{n}\gamma_{n}}{n!}\left|\beta-\alpha\right|\|Q\left(\left[\alpha,\beta\right)\right)\psi\|<C\left|\beta-\alpha\right|\|Q\left(\left[\alpha,\beta\right)\right)\psi\|\sum\limits_{n=0}^{\infty}\frac{\left|tM\right|^{n}}{n!}=Ce^{tM}\left|\beta-\alpha\right|\|Q\left(\left[\alpha,\beta\right)\right)\psi\|$.

Setting $\gamma=Ce^{tM}$, we get the desired result.
Third case: For functions in $\mathcal{F}_{3}$, we have the resolvent formula

$$\left(T-z\right)^{-1}-\left(S-z\right)^{-1}=\left(T-z\right)^{-1}\left(S-T\right)\left(S-z\right)^{-1}$$ (4.5)

for self-adjoint operators $T,S\in\mathcal{B}\left(\mathcal{H}\right)$ and $z\in\mathbb{C}\setminus\mathbb{R}$. In addition, $\mathcal{H}_{+}$ is invariant under $\left(J-z\right)^{-1}$. Thus, we have

$\|\left(J_{\theta}-z\right)^{-1}Q\left(\left[\alpha,\beta\right)\right)\psi-\left(J-z\right)^{-1}Q\left(\left[\alpha,\beta\right)\right)\psi\|=$

$=\|\left(J_{\theta}-z\right)^{-1}\left(J-J_{\theta}\right)\left(J-z\right)^{-1}Q\left(\left[\alpha,\beta\right)\right)\psi\|=\|\left(J_{\theta}-z\right)^{-1}\left(J-J_{\theta}\right)Q\left(\left[\alpha,\beta\right)\right)\left(J-z\right)^{-1}\psi\|\leq\|\left(J_{\theta}-z\right)^{-1}\|\cdot\gamma\left|\beta-\alpha\right|\|Q\left(\left[\alpha,\beta\right)\right)\left(J-z\right)^{-1}\psi\|=\|\left(J_{\theta}-z\right)^{-1}\|\gamma\left|\beta-\alpha\right|\|\left(J-z\right)^{-1}Q\left(\left[\alpha,\beta\right)\right)\psi\|\leq$

$\leq M^{2}\gamma\left|\beta-\alpha\right|\|Q\left(\left[\alpha,\beta\right)\right)\psi\|$

which implies the desired result. ∎

Remark 4.4.

Note that Theorem 4.3 holds also when taking $a,b\in\left(\frac{\pi}{2},\pi\right)$.

Theorem 1.5 now follows immediately from Theorem 4.3 and Lemma 2.7.

4.2 Integrating w.r.t. a projection-valued measure

Our goal in this section is to prove a (stronger version of) Theorem 1.5. We will use the following lemma:

Lemma 4.5.

Let $P$ be a PVM, $T\in\mathcal{B}\left(\mathcal{H}\right)$ self-adjoint and $F:\left[a,b\right]\to\mathcal{B}\left(\mathcal{H}\right)$. Suppose that there exists $M\in\mathbb{R}$ such that $\|T\|\leq M$ and for any $t\in\left[a,b\right]$, $\left\|F\left(t\right)\right\|\leq M$. In addition, suppose that there exists a dense subspace $\mathcal{H}_{+}\subseteq\mathcal{H}$ and $\gamma>0$ such that for any $a\leq\alpha<\beta\leq b$ and for any $\psi_{0}\in\mathcal{H}_{+}$,

$$\|F\left(t\right)P\left(\left[\alpha,\beta\right)\right)\psi_{0}-TP\left(\left[\alpha,\beta\right)\right)\psi\|<\gamma\left|\beta-\alpha\right|\|P\left(\left[\alpha,\beta\right)\right)\psi_{0}\|.$$ (4.6)

Then $F$ is integrable w.r.t. $P$ and

$\int_{a}^{b}F\left(t\right)dP\left(t\right)=TP\left(\left[a,b\right)\right)$.

Proof.

We first claim that (4.6) holds for any $\psi\in\mathcal{H}$. Indeed, given $\psi\in\mathcal{H}$ and $\psi_{0}\in\mathcal{H}_{+}$, we have

$\|F\left(t\right)P\left(\left[\alpha,\beta\right)\right)\psi-TP\left(\left[\alpha,\beta\right)\right)\psi\|\leq\|F\left(t\right)P\left(\left[\alpha,\beta\right)\right)\psi-F\left(t\right)P\left(\left[\alpha,\beta\right)\right)\psi_{0}\|+\|F\left(t\right)P\left(\left[\alpha,\beta\right)\right)\psi_{0}-TP\left(\left[\alpha,\beta\right)\right)\psi_{0}\|+\|TP\left(\left[\alpha,\beta\right)\right)\psi_{0}-TP\left(\left[\alpha,\beta\right)\right)\psi\|<\|F\left(t\right)\|\|\psi-\psi_{0}\|+\gamma\left|\beta-\alpha\right|\|P\left(\left[\alpha,\beta\right)\right)\psi_{0}\|+\|T\|\|\psi-\psi_{0}\|\leq\gamma\left|\beta-\alpha\right|\|P\left(\left[\alpha,\beta\right)\right)\psi_{0}\|+\left(\|F\left(t\right)\|+\|T\|+\|P\left(\left[\alpha,\beta\right)\right)\|\right)\|\psi-\psi_{0}\|\leq\gamma\left|\beta-\alpha\right|\|P\left(\left[\alpha,\beta\right)\right)\psi_{0}\|+\left(2M+1\right)\|\psi-\psi_{0}\|$.

Letting $\psi_{0}\to\psi$, we get (4.6). Now, let $\psi\in\mathcal{H}$, let $\Delta=\left(a=x_{0},\ldots,x_{n}=b\right)$ be any partition and let $\left(t_{j}\right)_{j=1}^{n}$ be any choice of points.

$\|\sum\limits_{j=1}^{n}F\left(t_{j}\right)P\left(\left[x_{j-1},x_{j}\right)\right)\psi-T\psi\|=\|\sum\limits_{j=1}^{n}F\left(t_{j}\right)P\left(\left[x_{j-1},x_{j}\right)\right)\psi-TP\left(\left[x_{j-1},x_{j}\right)\right)\psi\|<\sum\limits_{j=1}^{n}\gamma\left|x_{j}-x_{j-1}\right|\|P\left(\left[x_{j-1},x_{j}\right)\right)\psi\|\leq\gamma\left(\sum\limits_{j=1}^{n}\left|x_{j}-x_{j-1}\right|^{2}\right)^{\frac{1}{2}}\left(\sum\limits_{j=1}^{n}\|P\left(\left[x_{j-1},x_{j}\right)\right)\psi\|^{2}\right)^{\frac{1}{2}}\leq\gamma\left|\Delta\right|\sqrt{b-a}\|\psi\|\underset{\left|\Delta\right|\to 0}{\longrightarrow}0$,

as required. ∎

Lemma 4.5 and Theorem 4.3 imply

Theorem 4.6.

Let $\mathcal{F},\,Q$ be as in Theorem 4.3. Then for any $F\in\mathcal{F}$,

$$\int_{I}F\left(J_{\theta}\right)dQ\left(\theta\right)=F\left(J\right)Q\left(I\right),$$ (4.7)

where $I=\left[0,\pi\right)$.

Proof.

Let $a_{n}=\frac{1}{n}$, $b_{n}=\frac{\pi}{2}-\frac{1}{n}$. Then, by Theorem 4.3, (4.6) holds for $F\left(J_{\theta}\right)$ and $F\left(J\right)Q\left(\left[a_{n},b_{n}\right]\right)$. Thus, by Lemma 4.5 we get that

$\int_{a_{n}}^{b_{n}}F\left(J_{\theta}\right)dQ\left(\theta\right)=F\left(J\right)Q\left(\left[a_{n},b_{n}\right]\right)$.

Letting $n\to\infty$ and noting that $F\left(J\right)Q\left(\left[a_{n},b_{n}\right]\right)\to F\left(J\right)Q\left(\left(0,\frac{\pi}{2}\right)\right)$ as $n\to\infty$, we get that

$\int_{I_{1}}F\left(J_{\theta}\right)dQ\left(\theta\right)=F\left(J\right)Q\left(I_{1}\right)$,

where $I_{1}=\left(0,\frac{\pi}{2}\right)$. In a similar way, for $I_{2}=\left(\frac{\pi}{2},\pi\right)$,

$\int_{I_{2}}F\left(J_{\theta}\right)dQ\left(\theta\right)=F\left(J\right)Q\left(I_{2}\right)$.

and so

$\int_{I_{1}\cup I_{2}}F\left(J_{\theta}\right)dQ\left(\theta\right)=F\left(J\right)Q\left(I_{1}\cup I_{2}\right)$.

Now, note that for any $\theta\in\left[0,\pi\right)$ such that $Q\left(\left\{\theta\right\}\right)\neq 0$, $J|_{\text{Ran}Q\left(\left\{\theta\right\}\right)}=J_{\theta}|_{\text{Ran}Q\left(\left\{\theta\right\}\right)}$ and so for any $F\in\mathcal{F}$, $F\left(J\right)|_{\text{Ran}Q\left(\left\{\theta\right\}\right)}=F\left(J_{\theta}\right)|_{\text{Ran}Q\left(\left\{\theta\right\}\right)}$. Thus, trivially,

$\int_{\left\{\theta\right\}}F\left(J_{\theta}\right)dQ\left(\theta\right)=F\left(J\right)Q\left(\left\{\theta\right\}\right)$.

Thus, since $I=I_{1}\cup I_{2}\cup\left\{0\right\}\cup\left\{\frac{\pi}{2}\right\}$, we get (4.7), as required. ∎

Theorem 1.4 now follows from Theorem 4.6 and Lemma 2.7, along with the fact that $Q\left(I\right)=P\left(S\right)$, and $S$ supports $\mu_{s}$.

5 Examples and further discussion

5.1 An operator with a symmetric potential

Let us discuss a special case of Jacobi operators in which the potential $b$ satisfies some symmetry condition. For simplicity, in this example we will assume that $a\equiv 1$. We say that $b$ is symmetric around $\frac{1}{2}$ if for any $n\in\mathbb{Z}$, $b_{n}=b_{-n+1}$.

Proposition 5.1.

Suppose $b$ is symmetric around $\frac{1}{2}$. Then $Q$ is given by $Q=Q_{\text{Sym}}+Q_{\text{Sym}^{\perp}}$, where

$Q_{\text{Sym}}\left(A\right)=\begin{cases}P_{\text{Sym}}&-\frac{\pi}{4}\in A\\ 0&\text{otherwise}\end{cases},\,\,\,\,Q_{\text{Sym}^{\perp}}\left(A\right)=\begin{cases}\text{Id}-P_{\text{Sym}}&\frac{\pi}{4}\in A\\ 0&\text{otherwise}\end{cases}$

and $P_{\text{Sym}}$ is the orthogonal projection onto $\mathcal{H}_{\text{Sym}}\coloneqq\left\{\psi\in\mathcal{H}:\psi\left(n\right)=\psi\left(-n+1\right)\right\}$.

Proof.

By the symmetry of the potential, we get that $J_{+}^{0}=J_{-}^{\frac{\pi}{2}}$. For any $\theta\in\left[0,\pi\right)$, Let $\mu_{+}^{\theta}$ ($\mu_{-}^{\theta}$) be the spectral measure of $\delta_{1}$ ($\delta_{0}$) w.r.t. $J_{+}^{\theta}$ ($J_{-}^{\theta}$). Also let $m_{\pm}^{\theta}$ be the Borel transform of $\mu_{\pm}^{\theta}$. We claim that for any $\theta\in\left[0,\pi\right)$, for any $E\in\mathbb{R}$, $\underset{\epsilon\to 0}{\lim}\left|m_{\pm}^{\theta}\left(E+i\epsilon\right)\right|=\infty\iff\theta\in\left\{\frac{\pi}{4},\frac{3\pi}{4}\right\}$. The fact that $J_{+}^{0}=J_{-}^{\frac{\pi}{2}}$ along with (2.8) imply that

$$m_{-}^{\theta}=\frac{m_{-}^{\frac{\pi}{2}}}{1-\cot\theta m_{-}^{\frac{\pi}{2}}}=\frac{m_{+}^{0}}{1-\cot\theta m_{+}^{0}}$$ (5.1)

Now, suppose that $\underset{\epsilon\to 0}{\lim}\left|m_{\pm}^{\theta}\left(E+i\epsilon\right)\right|=\infty$. Then by (2.8) we get that $\underset{\epsilon\to 0}{\lim}m_{+}^{0}\left(E+i\epsilon\right)=\cot\theta$. Plugging this in (5.1), we get

$\underset{\epsilon\to 0}{\lim}\left|m_{-}^{\theta}\left(E+i\epsilon\right)\right|=\begin{cases}\left|\frac{\cot\theta}{1-\cot^{2}\theta}\right|&\theta\in\left[0,\pi\right)\setminus\left\{\frac{\pi}{4},\frac{3\pi}{4}\right\}\\ \infty&\theta\in\left\{\frac{\pi}{4},\frac{3\pi}{4}\right\}\end{cases}$,

as required. For any $\theta\in\left[0,\pi\right)$, let

$S_{\pm}^{\theta}=\left\{E\in\mathbb{R}:u_{\theta,E}^{\pm}\text{ is subordinate as a solution to (\ref{ev_eq1})}\right\}$.

Recall that by Theorem 2.12, $\mu_{s}$ is supported on the set $S$ which satisfies $S=\underset{\theta\in\left[0,\pi\right)}{\bigcup}S_{-}^{\theta}\cap S_{+}^{\theta}$. Now, by Remark 2.14, we have

$S_{\pm}^{\theta}=\left\{E\in\mathbb{R}:\underset{\epsilon\to 0}{\lim}\left|m_{\pm}^{\theta}\left(E+i\epsilon\right)\right|=\infty\right\}$,

and by what we’ve shown so far, if $\theta\in\left[0,\pi\right)\setminus\left\{\frac{\pi}{4},\frac{3\pi}{4}\right\}$ then $S_{+}^{\theta}\cap S_{-}^{\theta}=\emptyset$, and otherwise $S_{-}^{\theta}=S_{+}^{\theta}$. Thus, we get that $Q$ is supported on $\left\{\frac{\pi}{4},\frac{3\pi}{4}\right\}$. By the properties of the potential we get that for any $E\in S_{-}^{\frac{\pi}{4}}\cap S_{+}^{\frac{\pi}{4}}$, any subordinate solution $u$ of (2.6) satisfies $u\left(n\right)=-u\left(-n+1\right)$, and so $\text{Ran}Q\left(\left\{\frac{3\pi}{4}\right\}\right)=\mathcal{H}_{\text{Sym}}^{\perp}$. Similarly, we get that $\text{Ran}Q\left(\left\{\frac{3\pi}{4}\right\}\right)=\mathcal{H}_{\text{Sym}}$, which implies the desired result. ∎

5.2 An operator with pure point spectrum

Suppose now that $J$ has pure point spectrum. Recall that $S_{\theta}\coloneqq S_{-}^{\theta}\cap S_{-}^{\theta}$ is the set of energies for which there exists a subordinate solution to (2.6) which satisfies (2.4). In addition, in the case of pure point spectrum, we may assume that a solution is subordinate if and only if it is $\ell^{2}$ at $\pm\infty$, namely we can take a support $S$ to be the set of eigenvalues of $J$. Then, the set $S_{\theta}$ which is defined in Section 4 consists of the eigenvalues of $J$ for which the corresponding eigenvector $\psi$ satisfies

$\psi\left(0\right)\cos\theta+\psi\left(1\right)\sin\theta=0$.

For each $\theta\in\left[0,\pi\right)$ such that the set $S_{\theta}$ is non-empty, we have that $Q\left(S_{\theta}\right)\neq 0$, and so $Q$ is supported on a countable set $\left\{\theta_{n}\right\}_{n}\in\mathbb{N}$. Thus, for any $F\in\mathcal{F}$, where $\mathcal{F}$ is the family of functions given in Theorem 4.3, the integral formula (4.7) is given by

$\int_{I}F\left(J_{\theta}\right)dQ\left(\theta\right)=\underset{n\in\mathbb{N}}{\oplus}F\left(J_{\theta_{n}}\right)Q\left(\left\{\theta_{n}\right\}\right)=F\left(J\right)$,

where $I=\left[0,\pi\right)$.

6 Appendix - Proof of Proposition 2.18

Our goal in this appendix is to prove Proposition 2.18. For any $n\in\mathbb{N}$, we denote $\left[-n,n\right]\coloneqq\left\{-n,-\left(n-1\right),\ldots,0,\ldots,n-1,n\right\}$. We will use the following theorem, which is a special case of [18, Lemma 3.2]:

Lemma 6.1.

Let $J:\ell^{2}\left(\mathbb{Z}\right)\to\ell^{2}\left(\mathbb{Z}\right)$ be a Jacobi operator. Given $n\in\mathbb{N}$ and $z\in\mathbb{C}$, let $F\left(n,z\right)\in\mathcal{M}_{2n}\left(\mathbb{C}\right)$ be defined by

$\left(F\left(n,z\right)\right)_{kj}=\langle\delta_{k},\left(J-z\right)^{-1}\delta_{j}\rangle,\,\,\,\,\,\,\,\,\,k,j\in\left[-n,n\right]$.

Then $F\left(n,z\right)$ is invertible, and we have

$$\left(F\left(n,z\right)\right)^{-1}=A\left(n\right)+\operatorname{diag}\left(\frac{1}{m_{-n}\left(z\right)},\ldots,\frac{1}{m_{+n}\left(z\right)}\right)$$ (6.1)

where $A\left(n\right)\in\mathcal{M}_{2n}$ is given by

$A\left(n,z\right)=\left(\begin{matrix}\frac{1}{m_{-n}\left(z\right)}&a_{-n}&0&0&\cdots&0&0\\ a_{-n}&b_{-n+1}-z&a_{-n+1}&0&\cdots&0&0\\ 0&a_{-n+1}&b_{-n+2}-z&a_{-n+2}&\cdots&0&0\\ \vdots&\vdots&\vdots&\vdots&\ddots&\vdots&\vdots\\ 0&0&0&0&\cdots&b_{n-1}-z&a_{n-1}\\ 0&0&0&0&\cdots&a_{n-1}&\frac{1}{m_{+n}\left(z\right)}\end{matrix}\right)$

and $m_{\pm n}$ is the Borel transform of the spectral measure of $\delta_{\pm n}$ w.r.t. the restriction of $J$ to $\left\{\pm n,\pm n+1,\ldots\right\}$.

We now turn to prove Proposition 2.18.

Proof.

For any $k,j\in\mathbb{Z}$, $\mu_{kj}\ll\mu$. Let $f_{kj}=\frac{d\mu_{kj}}{d\mu}$. It is well known that for $\mu$-almost every $E\in\mathbb{R}$, we have

$f_{kj}\left(E\right)=\underset{\epsilon\to 0}{\lim}\frac{\mu_{kj}\left(\left(E-\epsilon,E+\epsilon\right)\right)}{\mu\left(\left(E-\epsilon,E+\epsilon\right)\right)}$

and so by Theorem 2.15, there exists a set $\widetilde{S}_{kj}\subseteq\mathbb{R}$ which supports $\mu_{s}$, such that for any $E\in\widetilde{S}_{kj}$, we have

$$\underset{\epsilon\to 0}{\lim}\frac{\mu_{kj}\left(\left(E-\epsilon,E+\epsilon\right)\right)}{\mu\left(\left(E-\epsilon,E+\epsilon\right)\right)}=\underset{\epsilon\to 0}{\lim}\frac{M_{kj}\left(E+i\epsilon\right)}{M\left(E+i\epsilon\right)}.$$ (6.2)

Let $\widetilde{S}=S\cap\underset{k,j\in\mathbb{Z}}{\bigcap}\widetilde{S}_{kj}$. It is easy to see that $\widetilde{S}$ supports $\mu_{s}$, and that (6.2) holds for any $E\in\widetilde{S}$ and for any $k,j\in\mathbb{Z}$. We denote the limit by $\alpha_{kj}$. By Theorem 2.13, we may assume that for any $E\in\widetilde{S}$,

$$\underset{\epsilon\to 0}{\lim}\left|M\left(E+i\epsilon\right)\right|=\infty.$$ (6.3)

Fix $n\in\mathbb{Z}$, $E\in\mathbb{R}$ and $k\in\left[-n.n\right]$. Let $e_{k}$ be the $k$’th element in the standard basis of $\mathbb{C}^{2n}$. By (6.3), we have

$0=\underset{\epsilon\to 0}{\lim}\frac{1}{\left|M\left(E+i\epsilon\right)\right|}=\underset{\epsilon\to 0}{\lim}\left\|\frac{e_{k}}{\left|M\left(E+i\epsilon\right)\right|}\right\|=\underset{\epsilon\to 0}{\lim}\left\|F\left(n,E+i\epsilon\right)^{-1}\left(F\left(n,E+i\epsilon\right)\frac{e_{k}}{M\left(E+i\epsilon\right)}\right)\right\|$.

Now, note that $F\left(n,E+i\epsilon\right)=\left(\begin{matrix}M_{1k}\left(E+i\epsilon\right)\\ \vdots\\ M_{nk}\left(E+i\epsilon\right)\end{matrix}\right)$ and so we get

$$0=\underset{\epsilon\to 0}{\lim}F\left(n,E+i\epsilon\right)^{-1}\left(\begin{matrix}\frac{M_{-nk}\left(E+i\epsilon\right)}{M\left(E+i\epsilon\right)}\\ \vdots\\ \frac{M_{nk}\left(E+i\epsilon\right)}{M\left(E+i\epsilon\right)}\end{matrix}\right)$$ (6.4)

Taking into account (6.1), we get that for any $j\in\left[-\left(n-1\right),n-1\right]$,

$$\underset{\epsilon\to 0}{\lim}a_{j-1}\frac{M_{\left(j-1\right)k}\left(E+i\epsilon\right)}{M\left(E+i\epsilon\right)}+a_{j}\frac{M_{\left(j+1\right)k}\left(E+i\epsilon\right)}{M\left(E+i\epsilon\right)}+\frac{M_{jk}\left(E+i\epsilon\right)}{M\left(E+i\epsilon\right)}\left(b_{j}-\left(E+i\epsilon\right)\right)=0,$$ (6.5)

Define $\widetilde{u}$ on $\left[-n,n\right]$ by $\widetilde{u}\left(j\right)=\alpha_{kj}$. By (6.5), we get that for any $j\in\left[-\left(n-1\right),n-1\right]$,

$a_{j-1}\widetilde{u}\left(j-1\right)+a_{j}\widetilde{u}\left(j+1\right)+b_{j}\widetilde{u}\left(j\right)=0$.

This implies that there exists a solution $u$ of (2.6) such that $u|_{\left[-n,n\right]}=\widetilde{u}|_{\left[-n,n\right]}$. We claim that $u$ is either a subordinate solution or $u$ vanishes identically. We will show this for $u|_{\left[n,\infty\right)}$, the proof for $u|_{\left(-\infty,-n\right]}$ is similar. Suppose $\alpha_{n}=0$. In addition, suppose that there exists a sequence $\left(\epsilon_{n}\right)_{n=1}^{\infty}$ such that $\underset{n\to\infty}{\lim}\epsilon_{n}=0$ and

$\underset{n\to\infty}{\lim}m_{+n}\left(E+i\epsilon_{n}\right)\neq 0$.

By (6.1) and (6.4), this implies that $a_{n-1}\alpha_{n-1}=0$. Thus, $u$ is a solution of (2.6) which satisfies $u\left(n-1\right)=u\left(n\right)=0$, which implies that it vanishes identically. If no such sequence $\left(\epsilon_{n}\right)_{n=1}^{\infty}$ exists, then $\underset{\epsilon\to 0}{\lim}m_{+n}\left(E+i\epsilon\right)=0$. The half-line subordinacy now implies (see Remark 2.14) that any non-zero solution which satisfies $u\left(n\right)=0$ is subordinate. This concludes the proof for the case $\alpha_{n}=0$. Suppose now that $\alpha_{n}\neq 0$. Again, by (6.1) and (6.4), we get that the limit

$\underset{\epsilon\to 0}{\lim}\frac{1}{m_{+n}\left(E+i\epsilon\right)}$

exists and is real, and so there exists $\theta\in\left[0,\pi\right)$ such that $\underset{\epsilon\to 0}{\lim}m_{+n}\left(E+i\epsilon\right)=\cot\theta$. By (2.9), we get $u|_{\left[n,\infty\right)}$ is subordinate if and only if $a_{n-1}u\left(n-1\right)\cos\theta+u\left(n\right)\sin\theta=0$. By (6.1) and (6.4), we get

$$a_{n-1}\alpha_{n-1}+\tan\theta\alpha_{n}=\underset{\epsilon\to 0}{\lim}a_{n-1}\frac{M_{\left(n-1\right)k}\left(E+i\epsilon\right)}{M\left(E+i\epsilon\right)}+\frac{M_{nk}\left(E+i\epsilon\right)}{M\left(E+i\epsilon\right)}\frac{1}{m_{+n}\left(E+i\epsilon\right)}=0$$ (6.6)

which implies that $u$ is subordinate, as required. For any $E\in\widetilde{S}$, let us denote by $\Gamma\left(E\right)$ the $\mathbb{Z}\times\mathbb{Z}$ matrix which is given by $\Gamma\left(E\right)_{kj}=\alpha_{kj}$. We conclude that for any $j\in\mathbb{Z}$ and for any $E\in\widetilde{S}$, the function $u:\mathbb{Z}\to\mathbb{C}$ which is defined by $u\left(k\right)=\Gamma\left(E\right)_{kj}$ is either a subordinate solution of (2.6) or identically zero. Thus, for any $n\in\mathbb{Z}$, the matrix $\Gamma_{n}\left(E\right)\in M_{2n}\left(\mathbb{C}\right)$ which is given by restricting $\Gamma\left(E\right)$ to $\left[-n,n\right]\times\left[-n,n\right]$ has rank one. We claim that it is also self-adjoint. For any $\epsilon>0$, let $I_{\epsilon}^{E}\coloneqq\left(E-\epsilon,E+\epsilon\right)$. For any $k,j\in\left[-n,n\right]$, we have

$\Gamma\left(E\right)_{kj}=\underset{\epsilon\to 0}{\lim}\frac{\mu_{kj}\left(I_{\epsilon}^{E}\right)}{\mu\left(I_{\epsilon}^{E}\right)}=\underset{\epsilon\to 0}{\lim}\frac{1}{\mu\left(I_{\epsilon}^{E}\right)}\langle P\left(I_{\epsilon}^{E}\right)\delta_{k},\delta_{j}\rangle=\underset{\epsilon\to 0}{\lim}\frac{1}{\mu\left(I_{\epsilon}^{E}\right)}\langle\delta_{k},P\left(I_{\epsilon}^{E}\right)\delta_{j}\rangle\overset{\left(\star\right)}{=}\underset{\epsilon\to 0}{\lim}\frac{1}{\mu\left(I_{\epsilon}^{E}\right)}\langle P\left(I_{\epsilon}^{E}\right)\delta_{j},\delta_{k}\rangle=\underset{\epsilon\to 0}{\lim}\frac{\mu_{jk}\left(I_{\epsilon}^{E}\right)}{\mu\left(I_{\epsilon}^{E}\right)}=\Gamma\left(E\right)_{jk}$.

Let us justify $\left(\star\right)$. First, note that $P\left(I_{\epsilon}^{E}\right)=\underset{n\to\infty}{\lim}q_{n}\left(J\right)$, where $\left(q_{n}\right)_{n=1}^{\infty}\subseteq\mathbb{R}\left[X\right]$. This implies that $\langle\delta_{k},P\left(I_{\epsilon}^{E}\right)\delta_{j}\rangle=\underset{n\to\infty}{\lim}\langle\delta_{k},q_{n}\left(J\right)\delta_{j}\rangle$, and this implies that $\langle\delta_{k},P\left(I_{\epsilon}^{E}\right)\delta_{j}\rangle\in\mathbb{R}$ as a limit of a sequence of real numbers, as required.

Since $\Gamma_{n}\left(E\right)$ is a self-adjoint matrix of rank-one, there exists $u_{n}\in\mathbb{C}^{n}$ such that $\Gamma\left(E\right)_{kj}=u_{n}\left(k\right)u_{n}\left(j\right)$. Furthermore, it is clear that for every $n,m\in\mathbb{N}$ such that $n>m$, $u_{n}|_{\left[-m,m\right]}=u_{m}$ and so there exists $u_{E}:\mathbb{Z}\to\mathbb{C}$ such that for any $k,j\in\mathbb{Z}$, $\left(\Gamma\left(E\right)\right)_{kj}=u_{E}\left(k\right)u_{E}\left(j\right)$. Finally, by what we have shown, $u$ must be a subordinate solution of (2.6), and

$u_{E}\left(0\right)^{2}+u_{E}\left(1\right)^{2}=\Gamma\left(E\right)_{00}+\Gamma\left(E\right)_{11}=\underset{\epsilon\to 0}{\lim}\frac{\mu_{00}\left(I_{\epsilon}^{E}\right)}{\mu\left(I_{\epsilon}^{E}\right)}+\frac{\mu_{11}\left(I_{\epsilon}^{E}\right)}{\mu\left(I_{\epsilon}^{E}\right)}=\underset{\epsilon\to 0}{\lim}\frac{\mu\left(I_{\epsilon}^{E}\right)}{\mu\left(I_{\epsilon}^{E}\right)}=1$,

as required. ∎