The change of the coordinates $\sigma \to \alpha$  allows us to rewrite $h_1$  in the form

$$\begin{array}{c}{h}_1\left(z\right)={\int }_{\left|\alpha \right|<\varepsilon }\frac{v\left(\alpha \right)}{|\alpha {|}^2-z}d\alpha ,\text{Im}z>0,\end{array}$$

(21)

where the function $v$  is analytic and

$$\begin{array}{c}v\left(0\right)=r\left({\sigma }^0\right)J\left({\sigma }^0\right)=r\left({\sigma }^0\right)|detA{|}^{-1/2}.\end{array}$$

(22)

Let $v_1$  and $v_2$  be even and odd parts of $v$  :

$$\begin{array}{c}{v}_1\left(\sigma \right)=\frac{v\left(\sigma \right)+v(-\sigma )}{2},v_2\left(\sigma \right)=\frac{v\left(\sigma \right)-v(-\sigma )}{2}.\end{array}$$

(23)

We substitute $v_1+v_2$  for $v$  in ( 23 ). Then $h_1=h_{1,1}+h_{1,2}$  where the terms in the right hand side are given by ( 21 ) with $v_1$  and $v_2,$  respectively, instead of $v.$  Obviously, $h_{1,2}=0,$  i.e. $h_1=h_{1,1}.$  Let $(\rho ,\theta )$  be the polar coordinates in the $\alpha$  plane. We write the integral for $h_{1,1}$  in the polar coordinates and integrate with respect to $\theta .$  Since $v_1$  depends analytically on ${\rho }^2,$  we arrive at

$$\begin{array}{c}{h}_1\left(z\right)={\int }_0^{\varepsilon }\frac{\rho w\left({\rho }^2\right)}{{\rho }^2-z}d\rho ,\text{Im}z>0,\end{array}$$

(24)

where the function $w$  is analytic and

$$\begin{array}{c}w\left(0\right)=2\pi v\left(0\right)=2\pi r\left({\sigma }^0\right)|detA{|}^{-1/2}.\end{array}$$

(25)

Formula ( 24 ) implies that $$h_1\left(z\right)={\int }_0^{\varepsilon }\frac{\rho \left[w\right({\rho }^2)-w(z\left)\right]}{{\rho }^2-z}d\rho +w\left(z\right){\int }_0^{\varepsilon }\frac{\rho }{{\rho }^2-z}d\rho .$$  The first integrand above is analytic in $\rho$  and $z$  when $\rho$  and $\left|z\right|$  are small $,$  and therefore the first integral is an analytic function in a neighborhood of $z=0.$  The second integral is equal to $\frac{1}{2}w\left(z\right)$  log $({\varepsilon }^2-z)-\frac{1}{2}w\left(z\right)$  log $(-z).$  The first term in the last sum is analytic when $\left|z\right|$  is small enough, and the second term is a sum of an analytic function $-\frac{\pi i}{2}w\left(z\right)$  and the function $-\frac{1}{2}w\left(z\right)$  log $z.$  This justifies ( 26 ) with $u_1\left(z\right)=-\frac{1}{2}w\left(z\right).$  The proof of the Lemma  4 is completed Proof. There exist analytic local coordinates $\alpha =({\alpha }_1,{\alpha }_2)$  in a neighborhood of ${\sigma }^0$  in which the function $\psi$  has the form $\psi =$  ${\alpha }_1^2-{\alpha }_2^2.$  Let $$D_1=\{\sigma =\sigma (\alpha ):|{\alpha }_1|<\varepsilon ,|{\alpha }_2|<2\varepsilon \}$$  We represent $h\left(z\right)$  as a sum of two terms $h_1\left(z\right)+h_2\left(z\right)$  where $h_2$  is the integral over $D\backslash D_1$  which is analytic in a neighborhood of $z=0,$  and the function $h_1$  is the integral over $D_1$  which has the following form

$$\begin{array}{c}{h}_1\left(z\right)={\int }_{-\varepsilon }^{\varepsilon }{\int }_{-2\varepsilon }^{2\varepsilon }\frac{v\left(\alpha \right)}{{\alpha }_1^2-{\alpha }_2^2-z}d{\alpha }_{2}d{\alpha }_1,Imz>0,\end{array}$$

(27)

where the function $v$  is analytic and ( 22 ) holds. Lemma  2 implies that the function $h_2$  is analytic in a neighborhood of $z=0.$  Hence, it remains to prove the validity of Lemma  4 forthe function ( 27 ). We represent the function $v$  as the sum $v_1+v_2$  where the function $v_1$  is even with respect to ${\alpha }_1$  and $v_2$  is odd. Obviously, the integral with the function $v_2$  instead of $v$  is equal to zero. Then we represent $v_1$  as the sum of an even with respect to ${\alpha }_2$  function and an odd with respect to ${\alpha }_2$  function. Similarly, only the even part remains. Thus, ( 27 ) can be rewritten in the form

$$\begin{array}{c}{h}_1\left(z\right)={\int }_{-\varepsilon }^{\varepsilon }{\int }_{-2\varepsilon }^{2\varepsilon }\frac{w({\alpha }_1^2,{\alpha }_2^2)}{{\alpha }_1^2-{\alpha }_2^2-z}d{\alpha }_{1}d{\alpha }_2=4{\int }_0^{\varepsilon }{\int }_0^{2\varepsilon }\frac{w({\alpha }_1^2,{\alpha }_2^2)}{{\alpha }_1^2-{\alpha }_2^2-z}d{\alpha }_{1}d{\alpha }_2,\text{Im}z>0,\end{array}$$

(28)

where $w$  is analytic and $$w\left(0\right)=r\left({\sigma }^0\right)|detA{|}^{-1/2}.$$  From ( 28 ) it follows that $$h_1\left(z\right)=4{\int }_0^{\varepsilon }{\int }_0^{2\varepsilon }\frac{w({\alpha }_1^2,{\alpha }_2^2)-w({\alpha }_2^2+z,{\alpha }_2^2)}{{\alpha }_1^2-{\alpha }_2^2-z}d{\alpha }_{1}d{\alpha }_2+4{\int }_0^{\varepsilon }{\int }_0^{2\varepsilon }\frac{w({\alpha }_2^2+z,{\alpha }_2^2)}{{\alpha }_1^2-{\alpha }_2^2-z}d{\alpha }_{1}d{\alpha }_2.$$  The first integrand above is analytic in $\alpha$  and $z,$  and therefore the first integral is an analytic function in a neighborhood of $z=0.$  Thus the statement of the Lemma has to be proved for the second term $g$  which, after integrating over ${\alpha }_1$  and replacing ${\alpha }_2$  by $\tau$  , takes the form

$$\begin{array}{c}g=2{\int }_0^{\varepsilon }\frac{w({\tau }^2+z,{\tau }^2)}{\sqrt{{\tau }^2+z}}\text{log}\frac{{\alpha }_1-\sqrt{{\tau }^2+z}}{{\alpha }_1+\sqrt{{\tau }^2+z}}{|}_{{\alpha }_1=0}^{{\alpha }_1=2\varepsilon }d\tau ,Imz>0.\end{array}$$

(29)

If $\mathbb{\varkappa }=\frac{{\alpha }_1-\sqrt{{\tau }^2+}z}{{\alpha }_1+\sqrt{{\tau }^2+z}}$  and Im $z>0,$  ${\alpha }_1>0,$  then Im $\mathbb{\varkappa }<0.$  The unique branch of the log $\mathbb{\varkappa }$  in ( 29 ) is chosen by the condition that $\pi <arg($  log $\mathbb{\varkappa })<2\pi$  when Im $z>0$  . Thus when Im $z>0,$  we have

$$\begin{array}{c}g=2{\int }_0^{\varepsilon }\frac{w({\tau }^2+z,{\tau }^2)}{\sqrt{{\tau }^2+z}}\text{log}\frac{2\varepsilon -\sqrt{{\tau }^2+z}}{2\varepsilon +\sqrt{{\tau }^2+z}}d\tau -2\pi i{\int }_0^{\varepsilon }\frac{w({\tau }^2+z,{\tau }^2)}{\sqrt{{\tau }^2+z}}d\tau .\end{array}$$

(30)

One can easily check that the function $\frac{1}{u}$  log $\frac{2\varepsilon -u}{2\varepsilon +u}$  is analytic in $u$  in the circle $\left|u\right|<2\varepsilon$  (if the function is defined as zero at $u=0$  ). Since the integrand in the first term in the right-hand side of ( 30 ) is analytic in $\tau$  and $z$  when $\left|z\right|<{\varepsilon }^2$  , we obtain that the first term in the right-hand side of ( 30 ) is analytic in $z,$  $\left|z\right|<{\varepsilon }^2$  . Hence, it remains to prove the statement of the Lemma for the function

$$\begin{array}{c}{g}_1\left(z\right)=-2\pi i{\int }_0^{\varepsilon }\frac{w({\tau }^2+z,{\tau }^2)}{\sqrt{{\tau }^2+z}}d\tau ,\text{Im}z>0.\end{array}$$

(31)

Obviously $g_1$  is an analytic function of $z$  when $z\in U_{\varepsilon }=\{z:|z|<{\varepsilon }^2,$  $z/\in {\overline{R}}_{-}\}$  , where $R_{-}$  is the negative part of the real axis. Let us compare the limit values of $g_1$  when $z\to -t\pm i0,$  $0<\sqrt{t}<\varepsilon .$  We represent $g_1$  in the form $g_1\left(z\right)=f_1\left(z\right)+f_2\left(z\right)$  , where $$f_1\left(z\right)=-2\pi i{\int }_0^{\sqrt{t}}\frac{w({\tau }^2+z,{\tau }^2)}{\sqrt{{\tau }^2+z}}d\tau ,f_2\left(\mathbb{\varkappa }\right)=-2\pi i{\int }_{\sqrt{t}}^{\varepsilon }\frac{w({\tau }^2+z,{\tau }^2)}{\sqrt{{\tau }^2+z}}d\tau ,0<\sqrt{t}<\varepsilon .$$  Then $$f_1\left(-t+i0\right)-f_1\left(-t-i0\right)=-4\pi {\int }_0^{\sqrt{t}}\frac{w({\tau }^2-t,{\tau }^2)}{\sqrt{t-{\tau }^2}}d\tau ,$$  and $$f_2\left(-t+i0\right)-f_2\left(-t-i0\right)=0.$$  Thus, after the substitution $\tau =\sqrt{t}s$  , we get

$$\begin{array}{c}{g}_1\left(-t+i0\right)-g_1\left(-t+i0\right)=-4\pi {\int }_0^{\sqrt{t}}\frac{w({\tau }^2-t,{\tau }^2)}{\sqrt{t-{\tau }^2}}d\tau =-4\pi {\int }_0^1\frac{w\left(t{s}^2-t,t{s}^2\right)}{\sqrt{1-s^2}}ds.\end{array}$$

(32)

We denote the latter function by ${\phi }_1={\phi }_1\left(t\right).$  Obviously, ${\phi }_1\left(z\right)$  is analytic in $z.$  Define

$$\begin{array}{c}{\phi }_2\left(z\right)=g_1\left(z\right)-\frac{1}{2\pi i}{\phi }_1\left(z\right)\text{log}z.\end{array}$$

(33)

Then $${\phi }_2\left(-t+i0\right)-{\phi }_2\left(-t-i0\right)=0,$$  i.e. ${\phi }_2\left(z\right)$  does not have a branch point at $z=0$  . Now in order to prove that ${\phi }_2$  is analytic at $z=0$  it is enough to show that $\left|{\phi }_2\left(z\right)\right|\le c$  log $\frac{1}{\left|z\right|}$  when $\left|z\right|$  is small enough. From ( 33 ) it follows that it is enough to get the corresponding logarithmic estimate for $g_1.$  We represent $g_1$  as the sum of two terms $g_1=g_{1,1}+g_{1,2}$  by writing the function $w$  in the integral ( 31 ) as $w_1+w_2$  where $w_1=w(0,-z),$  $w_2=w({\tau }^2+z,{\tau }^2)-w(0,-z).$  Since $$\left|w_2\right|=\left|w\right({\tau }^2+z,{\tau }^2)-w(0,-z\left)\right|\le C|{\tau }^2+z|\text{for small}|{\tau }^2+z|,$$  the function $\left|g_{1,2}\right|$  is bounded when $\left|z\right|$  is small enough. The integral defining $g_{1,1}$  can be easily evaluated: $$g_{1,1}=-2\pi iw(0,-z)\text{log}\frac{\varepsilon +\sqrt{{\varepsilon }^2+z}}{\sqrt{z}}.$$  Hence, $g_1$  has the logarithmic estimate as $\left|z\right|\to 0,$  and therefore $g_1\left(z\right)$  is analytic at $z=0.$  The latter, together with ( 33 ), proves ( 26 ) for $g_1$  with $u_1\left(z\right)=\frac{1}{2\pi i}{\phi }_1\left(z\right).$  Thus ( 14 ) is valid for $h$  with $$u_1\left(0\right)=2iw\left(0\right){\int }_0^1\frac{ds}{\sqrt{1-s^2}}=\pi iw\left(0\right)=\pi ir\left({\sigma }^0\right)|detA{|}^{-1/2}.$$  The proof of the Lemma is completed.

Proof of the Theorem  1 . Let us fix an arbitrary point $k^0\in \left[-\sqrt{8},\sqrt{8}\right]\backslash \cup k_s$  and put $k=k^0+\mathbb{\varkappa }$  in ( 12 ). Then the function $2\pi G$  defined by ( 12 ) takes the form ( 14 ) with $\psi \left(\sigma \right)=\phi \left(\sigma \right)-(k^0{)}^2$  and $z=2k^0\mathbb{\varkappa }+{\mathbb{\varkappa }}^2.$  One can easily check that the integral ( 14 ), which we got, satisfies all the assumptions of Lemma  2 (with $\Gamma$  being a closed curve). Thus $G(k^0+\mathbb{\varkappa },\xi )$  is analytic in $z$  when $\left|z\right|$  is small enough, and therefore it is analytic in $\mathbb{\varkappa }$  in a neighborhood of the point $\mathbb{\varkappa }=0.$  This proves the first statement of Theorem  1 .

The statement of Theorem  1 concerning the point $k=k_0=0$  is the direct consequence of Lemma  3 , since $2\pi G$  defined by ( 12 ) has the form ( 14 ) with $\psi \left(\sigma \right)=\phi \left(\sigma \right)$  and $z=k^2,$  and all the assumptions of Lemma  3 hold in this case (with ${\sigma }^0=0)$  . In order to prove the statement concerning the points $k=k_{\pm 2},$  we use the periodicity of the integrand in ( 12 ) and replace the square $T$  in ( 12 ) by $T^{\prime }=(0,2\pi )\times (0,2\pi )\subset R_{\sigma }^2.$  After that, we again can apply Lemma  3 with $\psi \left(\sigma \right)=\phi \left(\sigma \right)-8,$  $z=$  $k^2-8$  and ${\sigma }^0=(\pi ,\pi ).$  In order to prove the statement of the Theorem concerning the points $k=k_{\pm 1},$  we replace the square $T$  in ( 12 )by $T^{\prime \prime }=(\frac{-\pi }{2},\frac{3\pi }{2})\times (\frac{-\pi }{2},\frac{3\pi }{2})$  and split $T^{\prime \prime }$  in two parts by the line ${\sigma }_1={\sigma }_2.$  Lemma  4 can be applied to each of these two integrals with $\psi \left(\sigma \right)=\phi \left(\sigma \right)-4,$  $z=$  $k^2-4$  and the point ${\sigma }^0$  equal to $(0,\pi )$  in the first integral, and equal to $(\pi ,0)$  in the second integral.

The proof of the Theorem is completed.

Later, when the large time asymptotic behavior of the solutions to the Cauchy problem is analyzed, we shall need to know some properties of the coefficients $u_{1,s}$  and $u_{2,s}$  in ( 13 ).

To be more exact, we will need the following

Lemma 5 . Functions $u_{1,0}\left(0,\xi \right)$  and $u_{2,0}\left(0,\xi \right)$  have the following properties:

a) $u_{1,0}\left(0,\xi \right)=c_1$  is a constant, b) $u_{2,0}\left(0,\xi \right)=c_{1}ln\left|\xi \right|+c_2+o\left({\left|\xi \right|}^{-1}\right)$  as $\left|\xi \right|\to \infty .$ 

Proof. Let $\varsigma \in C_0^{\infty }\left(R_{\sigma }^2\right),$  and $\varsigma \left(\sigma \right)=1$  if $\left|\sigma \right|<1,$  $\varsigma \left(\sigma \right)=0$  if $\left|\sigma \right|>2.$  From ( 12 ) it follows that

$$\begin{array}{c}2\pi G\left(k,\xi \right)={\int }_{R^2}\varsigma \left(\sigma \right)\frac{e^{i\sigma \xi }}{\phi \left(\sigma \right)-k^2}d\sigma +{\int }_T\left(1-\varsigma \left(\sigma \right)\right)\frac{e^{i\sigma \xi }}{\phi \left(\sigma \right)-k^2}d\sigma .\end{array}$$

(34)

We denote the second term in the right hand side of ( 34 ) by $v=v(k,\xi ).$  Since $\phi \left(\sigma \right)=0$  only at $\sigma =0,$  and $1-\varsigma =0$  in a neighborhood of the origin, then $v$  is analytic in $k$  for small enough $\left|k\right|,$  and $$v=v(0,\xi )={\int }_T\left(1-\varsigma \left(\sigma \right)\right)\frac{e^{i\sigma \xi }}{\phi \left(\sigma \right)}d\sigma .$$  The integrand above is infinitely smooth. We apply the integration by parts to that integral, integrating $e^{i\sigma \xi }$  and differentiating the factor $\frac{1-\varsigma \left(\sigma \right)}{\phi \left(\sigma \right)}$  . Since $1-\varsigma$  is equal to one in a neighborhood of the boundary of $T$  and the other factors in the integrand are periodic, the contributions from the boundary will be cancelled, and the integration by parts leads to the estimate:

$$\begin{array}{c}v\left(0,\xi \right)=O\left({\left|\xi \right|}^{-\infty }\right)\text{as}\left|\xi \right|\to \infty .\end{array}$$

(35)

Consider $$w\left(k,\xi \right)={\int }_{R^2}\varsigma \left(\sigma \right)\frac{e^{i\sigma \xi }}{{\left|\sigma \right|}^2-k^2}d\sigma .$$  Let us study

$$\begin{array}{c}2\pi G-v-w={\int }_{R^2}\varsigma e^{i\sigma \xi }\left[\frac{1}{\phi \left(\sigma \right)-k^2}-\frac{1}{{\left|\sigma \right|}^2-k^2}\right]d\sigma ={\int }_{R^2}\varsigma e^{i\sigma \xi }\frac{{\left|\sigma \right|}^2-\phi \left(\sigma \right)}{\left(\phi \left(\sigma \right)-k^2\right)({\left|\sigma \right|}^2-k^2)}d\sigma .\end{array}$$

(36)

The numerator of the integrand in ( 36 ) is of order $O\left({\left|\sigma \right|}^4\right)$  as $\left|\sigma \right|\to 0$  . Thus, the integral has a limit as $k\to 0:$ 

$$\begin{array}{c}\left(2\pi G-v-w\right)\left(0,\xi \right)={\int }_{R^2}{e}^{i\sigma \xi }\psi \left(\sigma \right)d\sigma ,\text{where}\psi \left(\sigma \right)=\varsigma \frac{{\left|\sigma \right|}^2-\phi \left(\sigma \right)}{\phi \left(\sigma \right){\left|\sigma \right|}^2}.\end{array}$$

(37)

Since $\nabla \psi \left(\sigma \right)\in L^1,$  its Fourier transform is $o|\xi {|}^{-1}$  as $\left|\xi \right|\to \infty .$  From here, ( 37 ), and ( 35 ) it follows that

$$\begin{array}{c}\left|\left(2\pi G-w\right)\left(0,\xi \right)\right|=o\left({\left|\xi \right|}^{-1}\right)\text{as}\left|\xi \right|\to \infty .\end{array}$$

(38)

Consider $$w_1\left(k,\xi \right)={\int }_{R^2}\left(1-\varsigma \right)\frac{e^{i\sigma \xi }}{{\left|\sigma \right|}^2-k^2}d\sigma ,Imk>0.$$  Then $w_1$  is analytic in $k$  for $\left|k\right|<\varepsilon$  and

$$\begin{array}{c}{w}_1\left(0,\xi \right)=O\left({\left|\xi \right|}^{-\infty }\right)\text{as}\left|\xi \right|\to \infty .\end{array}$$

(39)

Thus,

$$\begin{array}{c}\left|\left(2\pi G-\left(w+w_1\right)\right)\right|\left(0,\xi \right)=O\left({\left|\xi \right|}^{-1}\right)\text{as}\left|\xi \right|\to \infty .\end{array}$$

(40)

Note that $(w+w_1)/2\pi$  is the inverse Fourier transform of the function $\frac{1}{{\left|\sigma \right|}^2-k^2}$  in $R_{\sigma }^2$  , i.e. $$w+w_1=c_0{H}_0\left(k\left|\xi \right|\right),$$  where $c_0$  is a constant, and $H_0$  is the Hankel function. Hence,

$$\begin{array}{c}w+w_1=2\pi [c_1\text{log}\left(k\left|\xi \right|\right)+c_2]+\alpha \left(k\left|\xi \right|\right),\end{array}$$

(41)

where the function $\alpha$  is continuous and $\alpha \left(0\right)=0.$  From here and ( 40 ) it follows that $G-c_1$  log $k$  has a limit as $k\to 0$  and

$$\begin{array}{c}limk\to 0\left(G-c_1\text{log}k\right)=c_1\text{log}\left|\xi \right|+c_2+o\left({\left|\xi \right|}^{-1}\right)\text{as}\left|\xi \right|\to \infty .\end{array}$$

(42)

On the other hand, from ( 13 ) it follows that the limit ${lim}_{k\to 0}[G-u_1\left(0,\xi \right)$  log $k]$  exists and is equal to $u_2\left(0,\xi \right).$  It can only happen if $u_1\left(0,\xi \right)=c_1$  and $u_2\left(0,\xi \right)$  is given by ( 42 ).

The proof is completed.

III. Analytic properties of the resolvent of the Schrödinger operator. Let $R_{\lambda }={\left(H-\lambda \right)}^{-1}$  be the resolvent of the perturbed operator $H=\left(-\Delta +q\right)$  , where the support of $q$  belongs to a square

$$\begin{array}{c}S=\{\xi \in Z^2:|{\xi }_i|\le m\}.\end{array}$$

(43)

We are going to study the kernel $R_{k^2}(\xi ,\eta )$  of the operator $R_{k^2}$  (this kernel is Green's function for $H-k^2$  ) when both $\xi ,\eta \in S,$  and therefore we shall consider the truncated resolvent ${\widehat{R}}_{k^2}=\chi R_{k^2}\chi ,$  where $\chi$  is the characteristic function of $S.$  Since the operator $H$  with a potential $q,$  whose support is bounded, may have at most a finite number of eigenvalues, the resolvent $R_{\lambda }$  (and therefore, the operator ${\widehat{R}}_{\lambda })$  is meromorphic in $\lambda$  in the complex $\lambda$  -plane with the cut along the segment $[0,8],$  and it has at most a finite number of poles which are eigenvalues of $H$  . The following Proposition is an immediate consequence of the statement above.

Proposition 6 . The operators $R_{k^2}$  and ${\widehat{R}}_{k^2}$  are meromorphic in $k\in \mathbf{C}\backslash \left[-\sqrt{8},\sqrt{8}\right].$  Outside of the segment $\left[-\sqrt{8},\sqrt{8}\right]$  they have poles of the first order at the points $k=\pm i{\sigma }_j$  for which ${\lambda }_j=-{\sigma }_j^2$  are negative eigenvalues of $H$  and at points $k=\pm {\rho }_j$  for which ${\lambda }_j={\rho }_j^2>8$  are positive eigenvalues of $H.$ 

The following statement was proved in [13] :

Proposition 7 . The operator $H$  does not have eigenvalues on the set $(0,4)\cup (4,8).$ 

It means that $H$  may have at most three eigenvalues imbedded into the continuous spectrum, and those possible eigenvalues are $\lambda =0,4,8.$  We consider the upper half plane ${\mathbf{C}}_{+}=\{k:$  Im $k>0\}$  as the image of the set $\{\lambda :0<$  arg $\lambda <2\pi \}$  under the map $k=\sqrt{\lambda }$  . Thus, the resolvent set of the operator $H$  in the $k$  -plane is ${\overline{\mathbf{C}}}_{+}\backslash \left[-\sqrt{8},\sqrt{8}\right]$  , where the points $k=\pm \sigma$  with the same $\sigma >\sqrt{8}$  correspond to the same value of $\lambda$  . The next statement concerns an analytic extension of the operator ${\widehat{R}}_{k^2}$  from the upper half plane ${\mathbf{C}}_{+}$  through the segment $\left[-\sqrt{8},\sqrt{8}\right].$  We denote by $P_j,$  $j=0,1,2,$  the orthogonal projections (in $L^2\left(Z^2\right)$  ) onto the eigenspaces of the operator $H$  with the eigenvalues $\lambda =0,4,8,$  respectively. If some of those values of $\lambda$  are not eigenvalues then the corresponding operators set to be zero. Let ${\widehat{P}}_j=\chi P_j\chi .$  Let ${\widehat{R}}_{k^2}($  mod $P)$  be the operator function ${\widehat{R}}_{k^2}$  by modulus of polynomial in $k$  operator functions.

Theorem 8 1). The operator ${\widehat{R}}_{k^2},$  defined for $k\in {\overline{\mathbf{C}}}_{+}\backslash \left[-\sqrt{8},\sqrt{8}\right],$  admits an analytic extension on $\left[-\sqrt{8},\sqrt{8}\right]\backslash \cup k_s$  , where $k_0=0,$  $k_{\pm 1}=\pm 2$  , and $k_{\pm 2}=\pm \sqrt{8}.$  2). It has branch points of the logarithmic type at $k=k_s,$  and there are integers ${\alpha }_s$  and ${\beta }_s$  and operators $A_s\ne 0$  and $B_s$  ( in the finite dimensional space $L$  of functions supported on $S$  ) such that the truncated resolvent ${\widehat{R}}_{k^2}$  has the following behavior when $k\to k_s$  , $k\ne k_s-i\sigma ,$  $\sigma \ge 0:$ 

$$\begin{array}{ccc}{\widehat{R}}_{k^2}\left(\text{mod}P\right)& =& {\widehat{P}}_{\left|s\right|}\frac{1}{k_s^2-k^2}+A_s{\left(k-k_s\right)}^{{\alpha }_s}{\text{log}}^{{\beta }_s}\left(k-k_s\right)+B_s{\left(k-k_s\right)}^{{\alpha }_s}{\text{log}}^{{\beta }_s-1}\left(k-k_s\right)\end{array}$$

$$\begin{array}{ccc}& & +O\left({\left(k-k_s\right)}^{{\alpha }_s}{\text{log}}^{{\beta }_s-2}\left(k-k_s\right)\right)\text{,}\end{array}$$

(44)

3). The kernels ${\widehat{R}}_{k^2}(\xi ,\eta )$  of the operators ${\widehat{R}}_{k^2}$  with $k=\pm \sigma +i\varsigma ,$  where $\sigma$  is real, $\varsigma \ge 0,$  are complex adjoint (and therefore, ${\alpha }_s={\alpha }_{\left|s\right|},$  ${\beta }_s={\beta }_{\left|s\right|})$  .

The following inequalities hold: a) ${\alpha }_0\ge -2,$  and ${\beta }_0\le -1$  when ${\alpha }_0=-2;$  b) Let $s=\pm 1,\pm 2.$  Then ${\alpha }_s\ge -1,$  and ${\beta }_s\le -1$  when ${\alpha }_s=-1.$  4). If $q\left(\xi \right)\ge 0$  and $q$  is not equal to zero identically, then ${\widehat{R}}_{k^2}$  is bounded in a neighborhood of $k=0.$  In particular, in this case $\lambda =0$  is not an eigenvalue of $H$  and ${\alpha }_0\ge 0.$ 

Remark. The reason to consider ${\widehat{R}}_{k^2}($  mod $P)$  in ( 44 ), but not ${\widehat{R}}_{k^2}$  , is the following:

the large time asymptotic behavior of the solutions $v=v(t,x)$  to the Cauchy problem ( 1 ) depends on the first singular in $k-k_s$  term in the asymptotic expansion for ${\widehat{R}}_{k^2}-{\widehat{P}}_{\left|s\right|}\frac{1}{k_s^2-k^2}$  as $k\to k_s$  . In fact, ( 44 ) will be proven for ${\widehat{R}}_{k^2},$  but it will be clear from the proof that one can omit integer nonnegative powers of $k-k_s$  and get ( 44 ) for ${\widehat{R}}_{k^2}($  mod $P).$  Proof. Let

$$\begin{array}{c}\left(-\Delta -k^2+q\right)u=f,\text{Im}k>0,\end{array}$$

(45)

where supports of $q$  and $f$  belong to $S$  . The formula

$$\begin{array}{c}u=R_{k^2}^{0}h,\text{Im}k>0,\end{array}$$

(46)

establishes a one-to-one correspondence between solutions $u=R_{k^2}f\in L^2\left(Z^2\right)$  of ( 45 ) and functions $h$  which are supported on $S$  and satisfy the equation

$$\begin{array}{c}h+q{R}_{k^2}^{0}h=f,\text{Im}k>0.\end{array}$$

(47)

In fact, if $u\in L^2\left(Z^2\right)$  is a solution of ( 45 ), then $\left(-\Delta -k^2\right)u=h$  with $h=f-qu.$  Hence, $u=R_{k^2}^{0}h$  and the support of $h$  belongs to $S.$  The substitution of ( 46 ) into ( 45 ) leads to the equation ( 47 ) for $h.$  On the other hand, if $h$  is supported on $S$  and ( 47 ) holds, then $u,$  given by ( 46 ), belongs to $L^2\left(Z^2\right)$  and satisfies ( 45 ). The latter can be checked by substituting ( 46 ) into ( 45 ).

The equation ( 47 ) has the form $T_{k}h=f$  , where the operator $T_k$  acts in the finite dimensional space $L$  of functions supported on $S$  . It also can be written in the form of a linear system:

$$\begin{array}{c}h\left(\xi \right)+\sum \eta \in Sq\left(\xi \right)G(k,\xi -\eta )h\left(\eta \right)=f\left(\xi \right),\xi \in S,\text{Im}k>0,\end{array}$$

(48)

where $G$  is Green's function of the operator $-\Delta -k^2$  (the kernel of $R_{k^2}^0).$  If one chooses the basis in $L,$  which consists of functions equal to one at one point of $S$  and equal to zero everywhere else, then the matrix $\left[T_k\right]$  of the operator $T_k$  in this basis is: $$\left[T_k\right]=\left[T_k\right(\xi ,\eta ){]}_{\xi ,\eta \in S}=E+[q\left(\xi \right)G(k,\xi -\eta )],$$  where $E$  is the identity matrix. Note also that ( 46 ) implies that

$$\begin{array}{c}{\widehat{R}}_{k^2}=\chi R_{k^2}^0{T}_k^{-1},\text{Im}k>0.\end{array}$$

(49)

From the first statement of Theorem  1 it follows that the matrix $\left[T_k\right]$  can be analytically extended on $\left[-\sqrt{8},\sqrt{8}\right]\backslash \cup k_s.$  The equation ( 45 ) is uniquely solvable if $k\ne i{\sigma }_j,$  and therefore the same is true for the equation ( 47 ). Thus, $det\left[T_k\right]\ne 0$  when Im $k>0,$  $k\ne i{\sigma }_j.$  Hence, the operator $T_k^{-1}$  is meromorphic on $\left[-\sqrt{8},\sqrt{8}\right]\backslash \cup k_s.$  From here, ( 49 ), and the first statement of Theorem  1 it follows that ${\widehat{R}}_{k^2}$  has a meromorphic extension on $\left[-\sqrt{8},\sqrt{8}\right]\backslash \cup k_s.$  On the other hand, it was proved in ([13] ) that ${lim}_{k\to k^0+i0}{\widehat{R}}_{k^2}$  exists for any $k^0\in \left[-\sqrt{8},\sqrt{8}\right]\backslash \cup k_s.$  This justifies the first statement of Theorem  8 .

The second statement of Theorem  8 can be proven absolutely similarly. The only difference is that the analytic extension of $\left[T_k\right]$  has branch points of the logarithmic type at $k=k_s,$  $s=0,\pm 1,\pm 2$  . To be more exact, the elements of $\left[T_k\right]$  in neighborhoods of the points $k=k_s$  are linear functions of log $(k-k_s)$  with coefficients analytic in $k$  (see ( 13 )). Let $N$  be the number of points in $S$  . By solving the system $\left[T_k\right]h=f$  using the Kramer rule we get that $T_k^{-1}$  in a neighborhood of the point $k=k_s$  has the form

$$\begin{array}{c}{T}_k^{-1}=\frac{{\widehat{A}}_s(k,\text{log}(k-k_s\left)\right)}{B_s(k,\text{log}(k-k_s\left)\right)},s=0,\pm 1,\pm 2,\end{array}$$

(50)

where the functions $B_s=det\left[T_k\right]$  are polynomials of order $N$  with respect to the second argument (i.e., log $(k-k_s))$  whose coefficients are analytic in $k$  , and ${\widehat{A}}_s$  are polynomials in log $(k-k_s)$  of order $N-1$  , whose coefficients are linear operators on $L$  (matrices) which depend analytically in $k$  . From here, ( 49 ), and Theorem  1 it follows that the formula ( 50 ) is valid for ${\widehat{R}}_{k^2}$  with the only difference that the order of the polynomials in the numerator is $N:$ 

$$\begin{array}{c}{\widehat{R}}_{k^2}=\frac{{\widehat{C}}_s(k,\text{log}(k-k_s\left)\right)}{D_s(k,\text{log}(k-k_s\left)\right)},\end{array}$$

(51)

where the functions $D_s$  are polynomials in log $(k-k_s)$  with analytic in $k$  coefficients, and ${\widehat{C}}_s$  are polynomials in log $(k-k_s)$  whose coefficients are analytic in $k$  operators. Then a similar representation is valid for the difference of ${\widehat{R}}_{k^2}$  and ${\widehat{P}}_{\left|s\right|}\frac{1}{k_s^2-k^2}:$ 

$$\begin{array}{c}{\widehat{R}}_{k^2}-{\widehat{P}}_{\left|s\right|}\frac{1}{k_s^2-k^2}=\frac{{\widehat{C}}_s}{D_s}-\frac{{\widehat{P}}_{\left|s\right|}}{k_s^2-k^2}=\frac{{\widehat{F}}_s(k,\text{log}(k-k_s\left)\right)}{G_s(k,\text{log}(k-k_s\left)\right)},\end{array}$$

(52)

where ${\widehat{F}}_s=\left(k_s^2-k^2\right){\widehat{C}}_s-D_s{\widehat{P}}_{\left|s\right|},$  $G_s=\left(k_s^2-k^2\right)D_s.$  The formula ( 44 ) immediately follows from here. The second statement of Theorem  8 is proved.

The first part of the third statement of the Theorem is an obvious consequence of the fact that

$$\begin{array}{c}{R}_{\overline{\lambda }}=\overline{R_{\lambda }}\text{if}\lambda /\in [0,\infty ).\end{array}$$

(53)

Now we are going to prove the second part. Let $d{E}_{\sigma },$  $\sigma$  is real, be the operator valued spectral measure for the operator $H$  . Let ${\mu }_{ac}\left(d\sigma \right)$  be the absolutely continuous part of the measure $d{E}_{\sigma },$  and let the operator valued function $\nu =\nu \left(\sigma \right)$  be the density of the measure ${\mu }_{ac}\left(d\sigma \right).$  The matrix elements ${\nu }_{\xi ,\eta },$  $\xi ,\eta \in Z^2,$  of the operator $\nu$  are the densities of the absolutely continuous parts of the scalar measures $(d{E}_{\sigma }{\delta }_{\xi },{\delta }_{\eta }),$  where ${\delta }_{\xi }$  is the function on $Z^2$  equal to one at a fixed point $\xi$  and equal zero elsewhere. We also shall consider the restriction $\widehat{\nu }=\chi \nu \left(\sigma \right)\chi$  of the operators $\nu =\nu \left(\sigma \right)$  onto the space $L$  of functions $\psi$  supported on $S.$  From Stone's formula and the first part of Theorem  8 it follows that the density $\widehat{\nu }\left(\sigma \right)$  in a neighborhood of $\sigma ={\lambda }_s$  can be obtained as

$$\begin{array}{c}\widehat{\nu }\left(\sigma \right)=\frac{1}{\pi }{lim}_{\lambda \to \sigma +i0}\text{Im}[{\widehat{R}}_{\lambda }-{\widehat{P}}_{\left|s\right|}\frac{1}{{\lambda }_s-\lambda }].\end{array}$$

(54)

Obviously, the operator function $\nu$  is summable in $\sigma$  :

$$\begin{array}{c}\widehat{\nu }\in L^1.\end{array}$$

(55)

The formula ( 52 ) implies a more exact asymptotic expansion of ${\widehat{R}}_{k^2}$  than the expansion ( 44 ). Namely,

$$\begin{array}{c}{\widehat{R}}_{k^2}-{\widehat{P}}_{\left|s\right|}\frac{1}{k_s^2-k^2}\sim {\sum }_{j\ge 0}(k-k_s{)}^{{\alpha }_s+j}\frac{{\widehat{T}}_{j,s}\left(\text{log}\right(k-k_s\left)\right)}{N_{j,s}\left(\text{log}\right(k-k_s\left)\right)},k\to k_s,\end{array}$$

(56)

where $N_{j,s}$  are polynomials, ${\widehat{T}}_{j,s}$  are operator valued polynomials, and ${\widehat{T}}_{0,s}$  is a non-zero operator function. Consider first the case $s=1$  or $2$  . Then ( 56 ) implies that the following expansion is valid when $\lambda$  is close to ${\lambda }_s$  and Im $\lambda >0:$ 

$$\begin{array}{c}{\widehat{R}}_{\lambda }-{\widehat{P}}_{\left|s\right|}\frac{1}{{\lambda }_s-\lambda }\sim {\sum }_{j\ge 0}(\lambda -{\lambda }_s{)}^{{\alpha }_s+j}\frac{{\widehat{T}}_{j,s}\left(\text{log}\right(\lambda -{\lambda }_s\left)\right)}{N_{j,s}\left(\text{log}\right(\lambda -{\lambda }_s\left)\right)},\end{array}$$

(57)

where the polynomials ${\widehat{T}}_{j,s},$  $N_{j,s}$  are different from those in ( 56 ), and still ${\widehat{T}}_{0,s}\ne 0$  . Together with ( 54 ) this leads to the following representation

$$\begin{array}{c}\widehat{\nu }\left(\sigma \right)={lim}_{\lambda \to \sigma +i0}\text{Im{}(\lambda -{\lambda }_s{)}^{{\alpha }_s}{\text{log}}^{{\beta }_s}(\lambda -{\lambda }_s\left)f\right({\text{log}}^{-1}(\lambda -{\lambda }_s))+O((\lambda -{\lambda }_s{)}^{{\alpha }_s+1}{\text{log}}^{{\gamma }_s}(\lambda -{\lambda }_s\left)\right)\text{},}\end{array}$$

(58)

where the operator function $f$  is analytic, $f\left(0\right)\ne 0$  , and ${\gamma }_s$  is an integer.

Assume that ${\alpha }_s\le -2,$  ${\tilde{\beta }}_s\ne 0.$  From ( 58 ) with $\sigma >{\lambda }_s$  and ( 55 ) we obtain that the operators $f\left(\sigma \right)$  with real $\sigma >{\lambda }_s$  are real valued (i.e. they map any real valued function into a real valued function). In particular, the operators $f^{\left(j\right)}\left(0\right)$  are real valued for all $j\ge 0$  .

Then ( 58 ) with $\sigma <{\lambda }_s$  implies that

$$\begin{array}{ccc}\widehat{\nu }\left(\sigma \right)& =& (\sigma -{\lambda }_s{)}^{{\alpha }_s}\text{Im}[{\text{log}}^{{\beta }_s}\left(\right|\sigma -{\lambda }_s|+i\pi )f\left(0\right)+{\text{log}}^{{\beta }_s-1}\left(\right|\sigma -{\lambda }_s|+i\pi )f^{\prime }\left(0\right)\end{array}$$

$$\begin{array}{ccc}& & +O\left({\text{log}}^{{\beta }_s-2}\right(|\sigma -{\lambda }_s|\left)\right]\end{array}$$

$$\begin{array}{ccc}& =& \pi (\sigma -{\lambda }_s{)}^{{\alpha }_s}[{\beta }_{s}f\left(0\right){\text{log}}^{{\beta }_s-1}|\sigma -{\lambda }_s|+O\left({\text{log}}^{{\beta }_s-2}\right(|\sigma -{\lambda }_s|\left)\right],\sigma <{\lambda }_s.\end{array}$$

(59)

This contradicts ( 55 ), and therefore the assumption is wrong. Let us assume that ${\alpha }_s\le -2,$  ${\beta }_s=0.$  Let $n$  be the smallest $j>0$  for which $f^{\left(j\right)}\left(0\right)\ne 0.$  The formula ( 58 ) leads to the following analog of ( 59 ): $$\widehat{\nu }\left(\sigma \right)=\pi (\sigma -{\lambda }_s{)}^{{\alpha }_s}[\frac{-n}{n!}{f}^{\left(n\right)}\left(0\right){\text{log}}^{-n-1}|\sigma -{\lambda }_s|+O\left({\text{log}}^{-n-2}\right(|\sigma -{\lambda }_s|\left)\right],\sigma <{\lambda }_s,$$  which also contradicts ( 55 ). Thus ${\alpha }_s\le -2$  requires ${\beta }_s=0,$  $f$  is a real valued operator which does not depend on $\lambda$  . Now the same arguments can be applied successively to the next terms in ( 57 ) for which ${\alpha }_s+j<-1$  . When ${\alpha }_s+j=-1,$  we can only get that

$$\begin{array}{c}\frac{{\widehat{T}}_{j,s}}{N_{j,s}}=M+O\left({\text{log}}^{-1}\right(\lambda -{\lambda }_s\left)\right),{\alpha }_s+j=-1,\end{array}$$

(60)

where $M$  is a real valued operator. So, we obtain that ( 57 ) has the form

$$\begin{array}{c}{\widehat{R}}_{\lambda }-{\widehat{P}}_{\left|s\right|}\frac{1}{{\lambda }_s-\lambda }={\sum }_{-{\alpha }_s>j\ge 0}(\lambda -{\lambda }_s{)}^{{\alpha }_s+j}{M}_{j,s}+O((\lambda -{\lambda }_s{)}^{-1}{\text{log}}^{-1}(\lambda -{\lambda }_s\left)\right)\end{array}$$

(61)

where $M_{j,s}$  are real valued and independent of $\lambda$  operators on $L$  . If ${\alpha }_s=-1,$  then the arguments above still lead to ( 60 ) when $j=0.$  Thus ( 61 ) is valid when ${\alpha }_s\le -1.$  Our next goal is to show that ( 61 ) with ${\alpha }_s<-1$  implies that $M_{0,s}$  are zero operators.

Let $\delta >0$  be so small that the interval $\omega =[{\lambda }_s,{\lambda }_s+\delta ]$  does not contain the eigenvalues of $H$  different from ${\lambda }_s.$  Let ${\omega }_{\varepsilon }=$  $[{\lambda }_s+i\varepsilon ,{\lambda }_s+\delta +i\varepsilon ]$  be the shift of $\omega .$  Stone's formula implies that

$$\begin{array}{c}{\int }_{\omega }\widehat{\nu }\left(\sigma \right)d\sigma =\frac{1}{\pi }{lim}_{\varepsilon \to +0}\text{Im}{\int }_{{\omega }_{\varepsilon }}[{\widehat{R}}_{\lambda }-{\widehat{P}}_{\left|s\right|}\frac{1}{{\lambda }_s-\lambda }]d\lambda .\end{array}$$

(62)

Assume that $M_{0,s}\ne 0.$  Then the existence of the limit in the right hand side above contradicts ( 61 ) if ${\alpha }_s<-1$  is odd. It also contradicts ( 61 ) if ${\alpha }_s<-1$  is even and $M_{1,s}\ne 0.$  One can easily show that ${\alpha }_s$  does not depend on the choice of the square $S$  if $S$  contains the support of $q$  (this follows from ( 49 ) and ( 48 )). Hence, one can take a bigger $S$  , so that there is a function $\psi$  for which $M_{0,s}\psi$  is not identically equal to zero on a smaller cube $S^{\prime }=\{\xi :|{\xi }_i|\le m-1\}\subset S$  (see ( 43 )). We apply the operator $$H-\lambda =-\Delta +q\left(\xi \right)-{\lambda }_s-(\lambda -{\lambda }_s)$$  to both sides of ( 61 ). The following relations hold in $S^{\prime }$  for any function $\psi$  defined on $Z^2:$  $$(H-\lambda ){\widehat{R}}_{\lambda }\psi =\psi ,(H-\lambda ){\widehat{P}}_{\left|s\right|}\frac{\psi }{{\lambda }_s-\lambda }=-{\widehat{P}}_{\left|s\right|}\psi ,Im\lambda >0.$$  The right hand sides here are independent of $\lambda .$  Thus, the negative powers of $\lambda -{\lambda }_s$  will be cancelled when $H-\lambda$  is applied to the right hand side of ( 61 ). This leads to the following relation in $S^{\prime }$  : $$(H-{\lambda }_s)M_{1,s}\psi =M_{0,s}\psi ,\text{if}{\alpha }_s\le -2.$$  Hence, the assumption $M_{0,s}\ne 0$  implies that $M_{1,s}\ne 0,$  and this leads to the contradiction discussed above. This proves that ${\alpha }_s\ge -1.$  If ${\alpha }_s=-1,$  then ( 61 ) implies that the limit in the right hand side of ( 62 ) exists and converges to $\frac{1}{2}{M}_{0,s}$  when $\delta =\left|\omega \right|\to 0.$  The left hand side in ( 62 ) converges to zero as $\delta \to 0.$  Thus, $M_{0,s}=0$  and ( 61 ) takes the form $${\widehat{R}}_{\lambda }-{\widehat{P}}_{\left|s\right|}\frac{1}{{\lambda }_s-\lambda }=O\left(\right(\lambda -{\lambda }_s{)}^{-1}{\text{log}}^{-1}(\lambda -{\lambda }_s)),s=1,2.$$  This proves the second part of statement 3 of Theorem  8 in the case $s=1,2.$  If $s=-1$  or $-2$  , the result follows from ( 53 ). The statement for $s=0$  can be justified similarly to the case $s=1,2.$  We leave it for the reader. This completes the proof of statement 3.

Next we shall show that ${\widehat{R}}_{k^2}$  is bounded in a neighborhood of $k=0$  when $q\ge 0$  . From ( 51 ) it follows that

$$\begin{array}{c}{\widehat{R}}_{k^2}=A{k}^{\alpha }{\text{log}}^{\beta }k+B{k}^{\alpha }{\text{log}}^{\beta -1}k+O\left(k^{\alpha }{\text{log}}^{\beta -2}k\right)\text{as}k\to 0.\end{array}$$

(63)

where $A$  is a non zero operator defined on the space $L$  of functions with supports in $S$  . One could get ( 63 ) from ( 44 ) with $\alpha ={\alpha }_0,$  $\beta ={\beta }_0$  (when $P_0=0$  ) $,$  or $\alpha =-2,$  $\beta =0$  (if $P_0\ne 0$  ) $.$  It is also possible to get ( 63 ) from ( 51 ). The last statement of the Theorem will be proved if we show that the first term in the right hand side of ( 63 ) is bounded as $k\to 0.$  Suppose, to the contrary, that $\alpha <0$  or $\alpha =0,$  $\beta >0$  . Let $f$  be a function supported on $S$  , such $Af$  is a non-zero function. We have

$$\begin{array}{c}\left(-\Delta -k^2\right)R_{k^2}f=f-q{R}_{k^2}f,\text{Im}k>0.\end{array}$$

(64)

Since $q$  and $f$  are supported on $S$  , one can replace $R_{k^2}$  in the right hand side above by ${\widehat{R}}_{k^2}.$  Thus, $$R_{k^2}f=R_{k^2}^0(f-q{\widehat{R}}_{k^2}f),\text{Im}k>0.$$  From here, ( 63 ), and Lemma  5 it follows that, for any $\xi \in Z^2,$  Im $k>0$  and $k\to 0,$  the function $\left(R_{k^2}f\right)\left(\xi \right)$  has the form

$$\begin{array}{c}{R}_{k^2}f=[c_1\text{log}k+u_2(0,\xi )+O(k\text{log}k\left)\right]*\{f-q[Af{k}^{\alpha }{\text{log}}^{\beta }k+Bf{k}^{\alpha }{\text{log}}^{\beta -1}k+O\left(k^{\alpha }lo{g}^{\beta -2}k\right)\left]\right\},\end{array}$$

(65)

where the convolution of two functions $g,$  $h$  on $Z^2$  is defined as $$(g*h)\left(\xi \right)={\sum }_{\eta }g(\xi -\eta )h\left(\eta \right).$$  Hence, the main term in the asymptotic expansion of $R_{k^2}f$  as $k\to 0$  is

$$\begin{array}{c}-c_1{k}^{\alpha }{\text{log}}^{\beta +1}k{\sum }_{\eta }\left(qAf\right)\left(\eta \right).\end{array}$$

(66)

On the other hand, $\left(R_{k^2}f\right)\left(\xi \right)=\left({\widehat{R}}_{k^2}f\right)\left(\xi \right)$  when $\xi \in S,$  and this value has a smaller order due to ( 63 ). Thus ( 66 ) is zero. Since $c_1\ne 0,$  we arrive at

$$\begin{array}{c}{\sum }_{\eta }\left(qAf\right)\left(\eta \right)=0.\end{array}$$

(67)

Now from ( 65 ) we get

$$\begin{array}{c}{R}_{k^2}f=v\left(\xi \right)k^{\alpha }{\text{log}}^{\beta }k+O\left(k^{\alpha }{\text{log}}^{\beta -1}k\right),\text{Im}k>0,k\to 0,\end{array}$$

(68)

where $$v\left(\xi \right)=-u_2(0,\xi )*qAf-c_1*qBf$$  Together with Lemma  5 this leads to the following asymptotics for $v\left(\xi \right):$ 

$$\begin{array}{c}v\left(\xi \right)=-c_{1}ln\left|\xi \right|*qAf+c_2+o\left(\right|\xi {|}^{-1}),|\xi |\to \infty .\end{array}$$

(69)

We will need more specific behavior of $v$  at infinity. First of all note that ( 67 ) implies that $$ln\left|\xi \right|*qAf={\sum }_{\eta }(ln|\xi -\eta |-ln|\xi \left|\right)qAf\left(\eta \right)=O\left(\right|\xi {|}^{-1}),|\xi |\to \infty ,$$  and therefore,

$$\begin{array}{c}\left|v\right(\xi \left)\right|<C<\infty .\end{array}$$

(70)

Let us denote by $\frac{\partial }{de},$  $e\in Z^2,$  the difference derivative in the direction of $e:$ 

$$\begin{array}{c}\frac{\partial h\left(\xi \right)}{de}=h(\xi +e)-h\left(\xi \right).\end{array}$$

(71)

Let us show that

$$\begin{array}{c}\frac{\partial v\left(\xi \right)}{de}=o\left(\right|\xi {|}^{-1}),|\xi |\to \infty .\end{array}$$

(72)

Formula ( 67 ) implies that

$$\begin{array}{ccc}\frac{\partial }{de}ln\left|\xi \right|*qAf& =& {\sum }_{\eta }(ln|\xi +e-\eta |-ln|\xi -\eta \left|\right)qAf\left(\eta \right)\end{array}$$

$$\begin{array}{ccc}& =& {\sum }_{\eta }(ln|\xi +e-\eta |-ln|\xi -\eta |-ln|\xi +e|+ln|\xi \left|\right)qAf\left(\eta \right)=O\left(\right|\xi {|}^{-2}).\end{array}$$

The last equality can be easily obtained by evaluating the asymptotics for the combination of the logarithms above. The formula above and ( 69 ) justify ( 72 ).

Now we multiply both sides of ( 64 ) by $k^{-\alpha }lo{g}^{-\beta }k$  and pass to the limit as $k\to 0.$  This leads to

$$\begin{array}{c}\left(-\Delta +q\right)v=0.\end{array}$$

(73)

We apply the Green's Theorem for the difference Schrödinger operator to $v$  over the square $-N\le {\xi }_1,{\xi }_2\le N$  . Estimates ( 70 ), ( 69 ) allow us to pass to the limit as $N\to \infty$  .

This leads to the following result:

$$\begin{array}{c}0={\sum }_{\xi \in Z^2}\left(-\Delta +q\right)v\cdot \overline{v}={\sum }_{\xi \in Z^2}\left(\right|{\nabla }_{\xi }v{|}^2+q\left|v{|}^2\right),\end{array}$$

(74)

where ${\nabla }_{\xi }=(\frac{\partial }{\partial e_1},\frac{\partial }{\partial e_2}),$  $e_1=(1,0),$  $e_2=(0,1)$  [(+,)] and $\frac{\partial }{\partial e}$  is defined in ( 71 ). Since $q\ge 0$  , ( 74 ) implies that $v$  is a constant. Now from ( 73 ) (or ( 74 )) it follows that $v=0,$  since $q$  is a non-zero function. This conclusion contradicts ( 68 ) and ( 63 ). In fact, the later two relations imply that $v\left(\xi \right)=Af$  when $\xi \in S.$  Yet, $f$  was chosen so that $Af$  is a non-zero function on $S$  . The contradiction proves that our assumption of unboundedness of the first term in the right hand side of ( 63 ) is wrong. This completes the proof of the last statement of Theorem  8 .

IV. The large time behavior of solutions to the Cauchy problem. The solution $v=v\left(t,\xi \right)$  to the Cauchy problem ( 1 ) is given by ( 7 ). We assume that $f$  and $q$  are real valued and have bounded supports, and we will study the asymptotic behavior of $v$  when $t\to \infty$  and $\xi \in S$  , where $S$  is a square in $Z^2$  which contains the supports of $q$  and $f.$  In order to keep in mind this restriction on $\xi$  we multiply both sides of ( 7 ) by the characteristic function $\chi$  of $S$  . Since the support of $f$  belongs to $S$  , we can replace $f$  in ( 7 ) by $\chi f$  , and we arrive to

$$\begin{array}{c}\chi v\left(t,\xi \right)=\frac{1}{2\pi }{\int }_{B-i\infty }^{B+i\infty }{\widehat{R}}_{k^2}f{e}^{-ikt}dk,B>A.\end{array}$$

(75)

Let us denote by $P_{\left(j\right)},$  $1\le j\le n_1,$  the projection operators in $L^2\left(Z^2\right)$  onto the eigenspaces of the operator $H$  with the negative eigenvalues $\lambda =-{\sigma }_j^2.$  Let $P_j,$  $1\le j\le n_2,$  be the projection operators onto the eigenspaces of the operator $H$  with the positive eigenvalues $\lambda ={\rho }_j^2.$  We consider here all positive eigenvalues, including those which belong to the continuous spectrum (see the Proposition  7  $)$  . Let $P_0$  be the projection operator onto the eigenspace with the eigenvalue $\lambda =0$  , and let $P_0=0$  if $\lambda =0$  is not an eigenvalue. Let ${\widehat{P}}_{\left(j\right)}=\chi P_{\left(j\right)}\chi ,$  ${\widehat{P}}_j=\chi P_j\chi ,$  ${\widehat{P}}_0=\chi P_0\chi .$ 

Theorem 9 . The solution $v$  of the problem ( 1 ) has the following asymptotic behavior as $t\to \infty :$ 

$$\begin{array}{ccc}\chi v\left(t,\xi \right)& =& {\sum }_j\frac{e^{{\sigma }_{j}t}}{2{\sigma }_j}{\widehat{P}}_{\left(j\right)}f+{\sum }_j\frac{sin{\rho }_{j}t}{{\rho }_j}{\widehat{P}}_{j}f+t{\widehat{P}}_{0}f\end{array}$$

$$\begin{array}{ccc}& & +{\sum }_{0\le s\le 2}{a}_s{t}^{-{\alpha }_s-1}{\left(\text{log}t\right)}^{{\beta }_s-{\gamma }_s}sin(k_{s}t+{\omega }_s)A^{s}f+w\left(t,\xi \right),\end{array}$$

where $a_s\ne 0,$  ${\omega }_s$  are some constants, ${\alpha }_s,$  ${\beta }_s,$  $A^s$  are defined in ( 44 ), ${\gamma }_s=0$  if ${\alpha }_s=-1,$  ${\gamma }_s=1$  if ${\alpha }_s\ge 0,$  $k_0=0,$  $k_1=2,$  $k_2=\sqrt{8},$  and $\frac{d^n}{d{t}^n}w(t,\xi )$  decay as $t\to \infty$  by the factor log ${}^{-1}t$  faster than the slowest term in the last sum in the right hand side above.

Proof. We plan to move the contour of integration in ( 75 ) down and replace it by the following contour $\Gamma$  . Let $c_s$  be the following rays in the complex $k$  -plane $\mathbf{C}:$  $$c_s=\{k=k_s-i\sigma ,\sigma \ge 0\},s=0,\pm 1,\pm 2,$$  and let ${\mathbf{C}}^{\prime }=$  $\mathbf{C}\backslash \cup c_s$  be the complex $k$  -plane $\mathbf{C}$  with the cuts along $c_s.$  Let ${\Gamma }_0$  be the line Im $k=-\delta ,$  $\delta >0,$  without points $k_s-i\delta$  , where the line intersects $c_s.$  Let ${\gamma }_s$  be smooth loops joining the points $k_s-i\delta -0$  (on the left side of $c_s$  ) and $k_s-i\delta +0$  (on the right side of $c_s$  ) going round the cuts $c_s$  and lying in the $\delta$  -neighborhood of the points $k_s.$  Then $\Gamma ={\Gamma }_0\cup {\gamma }_s$  is the union of ${\Gamma }_0$  and all ${\gamma }_s,$  and we need only to chose $\delta$  in order to specify $\Gamma .$  From Proposition  6 it follows that there is a ${\delta }_1>0$  such that the operator ${\widehat{R}}_{k^2}$  defined on $\mathbf{C}\backslash [\sqrt{8},\sqrt{8}]$  does not have poles in the strip $-{\delta }_1\le$  Im $k<0.$  Below we shall always consider ${\widehat{R}}_{k^2}$  in ${\mathbf{C}}^{\prime },$  i.e. we start with ${\widehat{R}}_{k^2}$  in the upper half plane ${\mathbf{C}}_{+}$  and then extend it analytically down. Due to Proposition  6 , this extension exists in the region $|$  Re $k|>\sqrt{8}$  (which does not contain cuts $c_s$  ). According to Theorem  8 , there exists ${\delta }_2>0$  such that ${\widehat{R}}_{k^2}$  has an analytic (without poles) extension onto the ${\delta }_2$  -neighborhood in ${\mathbf{C}}^{\prime }$  of the segment $[\sqrt{8},\sqrt{8}]$  with branch points at $k=k_s$  . In fact, one can show that ${\widehat{R}}_{k^2}$  can be meromorphically extended onto the whole ${\mathbf{C}}^{\prime },$  but we don't need this statement. We choose an arbitrary $\delta <min({\delta }_1,{\delta }_2).$  Then ${\widehat{R}}_{k^2}$  is meromorphic strictly inside the region $\Omega \subset {\mathbf{C}}^{\prime }$  between the contour of integration in ( 75 ) and $\Gamma ,$  and ${\widehat{R}}_{k^2}$  has there poles only at points $k=i{\sigma }_j,$  ${\sigma }_j>0$  , and $k=\pm {\rho }_j$  for which $\lambda =k^2$  are eigenvalues of $H$  , which do not belong to the absolutely continuous spectrum of $H$  .

Note that $$\left|\right|{\widehat{R}}_{k^2}\left|\right|\le \left|\right|R_{k^2}\left|\right|\le 1/d,k\in \mathbf{C}\backslash [\sqrt{8},\sqrt{8}],$$  where $d$  is the distance in the complex $\lambda$  -plane between the point $\lambda =k^2$  and the spectrum of the operator $H$  . Thus,

$$\begin{array}{c}\left|\right|R_{k^2}\left|\right|\le C/|k{|}^2,k\in \Omega ,|k|\to \infty .\end{array}$$

(76)

Since $\left|\right|{\widehat{R}}_{k^2}\left|\right|\le \left|\right|R_{k^2}\left|\right|,$  the estimate ( 76 ) is valid also for $\left|\right|{\widehat{R}}_{k^2}\left|\right|.$  This allows us to use the Cauchy theorem and reduce the integral ( 75 ) to the following form:

$$\begin{array}{c}\chi v\left(t,\xi \right)=-i{\sum }_j\text{Res}k=i{\sigma }_j\text{(}{\widehat{R}}_{k^2}f{e}^{-ikt})-i\widetilde{{\sum }_j}\text{Res}k=\pm {\rho }_j\text{(}{\widehat{R}}_{k^2}f{e}^{-ikt})+\frac{1}{2\pi }{\int }_{\Gamma }{\widehat{R}}_{k^2}f{e}^{-ikt}dk,\end{array}$$

(77)

where the sum does not include points $\pm {\rho }_j\in [-\sqrt{8},\sqrt{8}].$  Let us estimate the last term in ( 77 ). We split ${\Gamma }_0={\Gamma }_1\cup {\Gamma }_2,$  where ${\Gamma }_1$  is the bounded part of ${\Gamma }_0$  for which $|$  Re $k|<\sqrt{8}+1$  and ${\Gamma }_2$  is the part of ${\Gamma }_0$  where $|$  Re $k|>\sqrt{8}+1.$  We denote

$$\begin{array}{c}{u}_s\left(t,\xi \right)=\frac{1}{2\pi }{\int }_{\Gamma s}{\widehat{R}}_{k^2}f{e}^{-ikt}dk,s=1,2,\end{array}$$

(78)

Since $\left|\right|{\widehat{R}}_{k^2}\left|\right|$  is bounded and $\left|e^{-ikt}\right|=e^{-\delta t}$  on ${\Gamma }_1,$  $$\left|\frac{{\partial }^j}{{\partial }^{j}t}{u}_1\left(t,\xi \right)\right|\le \frac{1}{2\pi }{\int }_{\Gamma s}\left|k^j{\widehat{R}}_{k^2}f{e}^{-ikt}\right|dk\le C_j{e}^{-\delta t},j=0,1,...\text{.}$$  A similar estimate for $u_2$  with $j=0$  immediately follows from ( 76 ). To get the estimate for the derivatives we differentiate the equation $$\left(-\Delta +q+k^2\right)u=f,k\in {\Gamma }_2,u=R_{k^2}f,$$  with respect to $k$  : $$\left(-\Delta +q+k^2\right)u_k=-2ku.$$  From here and ( 76 ) it follows that $\left|\right|\frac{\partial }{\partial k}{R}_{k^2}\left|\right|\le C/|k{|}^3,$  $k\in {\Gamma }_2.$  If we differentiate the equation again we can estimate the second derivative of the resolvent, and so on. Thus,

$$\begin{array}{c}\left|\right|\frac{{\partial }^j}{\partial k^j}{R}_{k^2}\left|\right|\le C/|k{|}^{2+j},k\in {\Gamma }_2.\end{array}$$

(79)

We apply the integration by parts $j$  times to the integral ( 78 ) with $s=2:$  $$u_2\left(t,\xi \right)={\sum }_{1\le l\le j}\left[\frac{{\partial }^{l-1}}{\partial k^{l-1}}\right({\widehat{R}}_{k^2}f\left)\frac{e^{-ikt}}{(it{)}^l}\right]|{}_{k=-\sqrt{8}-1-i\delta }^{k=\sqrt{8}+1-i\delta }+\frac{1}{2\pi }{\int }_{{\Gamma }_2}\frac{{\partial }^j}{\partial k^j}\left({\widehat{R}}_{k^2}f\right)\frac{e^{-ikt}}{(it{)}^j}dk.$$  This formula and ( 79 ) allow one to estimate the derivatives of $u_2:$  $$\left|\frac{{\partial }^j}{{\partial }^{j}t}{w}_2\left(t,\xi \right)\right|\le C_j{e}^{-\delta t},j=0,1,...\text{.}$$  Now ( 77 ) takes the form

$$\begin{array}{c}\chi v\left(t,\xi \right)=-i{\sum }_j\text{Res}k=i{\sigma }_j\text{(}{\widehat{R}}_{k^2}f{e}^{-ikt})-i\widetilde{{\sum }_j}\text{Res}k=\pm {\rho }_j\text{(}{\widehat{R}}_{k^2}f{e}^{-ikt})+{\sum }_s\frac{1}{2\pi }{\int }_{{\gamma }_s}{\widehat{R}}_{k^2}f{e}^{-ikt}dk+u,\end{array}$$

(80)

where $u=u_1+u_2,$  and therefore $$\left|\frac{{\partial }^j}{{\partial }^{j}t}u\left(t,\xi \right)\right|\le C_j{e}^{-\delta t},j=0,1,...\text{.}$$  Since the operator $R_{\lambda }$  in a neighborhood of the eigenvalue $\lambda ={\rho }_j^2>8$  has the form $$R_{\lambda }=P_j\frac{1}{{\rho }_j^2-\lambda }+Q_{\lambda },$$  where $P_j$  is the projection operator onto the corresponding eigenspace and $Q_{\lambda }$  is analytic in a neighborhood of ${\rho }_j^2,$  we obtain $$\text{Res}k={\rho }_j\text{(}{\widehat{R}}_{k^2}f{e}^{-ikt})+\text{Res}k=-{\rho }_j\text{(}{\widehat{R}}_{k^2}f{e}^{-ikt})=\frac{isin{\rho }_{j}t}{{\rho }_j}{\widehat{P}}_{j}f\text{, where}{\widehat{P}}_j=\chi P_j\chi .$$  Similarly, $$\text{Res}k=i{\sigma }_j\text{(}{\widehat{R}}_{k^2}f{e}^{-ikt})=\frac{i{e}^{{\sigma }_{j}t}}{2{\sigma }_j}{\widehat{P}}_{\left(j\right)}f.$$  Since the integral of any analytic (in particular, polynomial) function over ${\gamma }_s$  is zero, we can replace ${\widehat{R}}_{k^2}$  in the last term in ( 80 ) by the right hand side in ( 44 ). After that the corresponding integrals involving the first term from the right hand side of ( 44 ) can be evaluated with the help of the residue theorem. Together with the two formulas above this allows us to rewrite ( 80 ) in the form

$$\begin{array}{c}\chi v\left(t,\xi \right)={\sum }_{j=1}^{n_1}\frac{e^{{\sigma }_{j}t}}{2{\sigma }_j}{\widehat{P}}_{\left(j\right)}f+{\sum }_{j=1}^{n_2}\frac{sin{\rho }_{j}t}{{\rho }_j}{\widehat{P}}_{j}f+t{\widehat{P}}_{0}f+{\sum }_s\frac{1}{2\pi }{\int }_{{\gamma }_s}{Q}_s\left(k\right)f{e}^{-ikt}dk+u,\end{array}$$

(81)

where $$Q_s\left(k\right)={\widehat{R}}_{k^2}\left(\text{mod}P\right)-{\widehat{P}}_{\left|s\right|}\frac{1}{k_s^2-k^2}$$  is defined in ( 44 ), and the first three terms in the right hand side of ( 81 ) contain the projections on all eigenspaces of $H$  , including those for which ${\lambda }_j\in [0,8]$  . In order to complete the proof of the theorem, it remains only to evaluate the last term in ( 81 ) using Theorem  8 and the following simple statement which can be found in [14] , Ch X, Lemmas 7, 8:

Lemma 10 Let $$v\left(t\right)={\int }_l\omega \left(k\right)e^{-ikt}dk,$$  where the function $\omega$  is analytic when $\left|k\right|<2\delta$  with a branch point of the logarithmic type at $k=0,$  and $l$  is a smooth contour in the $\delta -$  neighborhood of $k=0$  starting at $k=-i\delta -0,$  ending at $k=-i\delta +0$  and not having common points with the negative imaginary semiaxis except the end points $k=-i\delta \pm 0.$  Let $\omega \left(k\right)$  have the following asymptotic behavior as $k\to 0,$  $\frac{3\pi }{2}<$  arg $k<-\frac{\pi }{2}:$  $$\omega \left(k\right)=k^p{\text{log}}^{q}k+c{k}^p{\text{log}}^{q-1}k+O\left(k^{p}lo{g}^{q-2}k\right)$$  where $p$  and $q$  are integers and $q\ne 0$  if $p\ge 0.$  Then $$v\left(t\right)=a{t}^{-p-1}{\text{log}}^{q-\gamma }t+w\left(t\right),$$  where $\gamma =0$  if $p<0,$  $\gamma =1$  if $p\ge 0,$  $a=a(p,q)$  is a constant and $$\left|\frac{d^j}{d{t}^j}w\left(t\right)\right|<C_j|\frac{d^j}{d{t}^j}\left[t^{-p-1}{\text{log}}^{q-\gamma -1}t\right],t\to \infty ,j=0,1,2....$$ 

References

1. Bourgain, J.: On random Schrödinger operators on $Z^2,$  Discrete and Continuous Dynamical Systems, 8, 1 (2002).
2. Bentosela, F., Briet, Ph., Pastur, L.: On the spectral and wave propagation properties of the surface Maryland model, J. Math. Physics, 44, no 1 (2003), 1-35.
3. Chahrour, A., Sahbani, J.: On the spectral and scattering theory of the Schrödinger operator with surface potential, Rev. Math. Phys., 12, 561 (2000).
4. Hundertmark, D., Kirsh, W.: Spectral theory of the sparse potentials, In Stochastic processes, physics and geometry: new interplays, I (Leipzig, 1999), CMS Conf. Proc., 28, 213 (2000).
5. Jaksic, V., Last, Y.: Surface states and spectra, Comm. Math. Phys., 218, 459 (2001).
6. Jaksic, V., Molchanov, S.: Wave operators for the surface Maryland model, J. Math. Physics, 41, 4452 (2000).
7. Jaksic, V., Molchanov, S., Pastur, L.: On the propagation properties of surface waves, In Wave Propagation in Complex Media, IMA Vol. Math. Appl., 96, 143 (1998).
8. Lax, P., Phillips, R.: Scattering Theory, Academic Press, New York, 1967.
9. Molchanov, S., Vainberg, B.: Scattering on the system of the sparse bumps: multidimensional case. Appl. Analysis, 71, 167 (1999).
10. Molchanov, S., Vainberg, B.: Spectrum of multidimensional Schrödinger operators with sparse potentials, In Analytical and Computational Methods in Scattering and Applied Mathematics, Chapman and Hall/CRC, Res. Notes Math., 417, 231 (2000).
11. Poulin, P.: Molchanov-Vainberg Laplacian, preprint (2004).
12. Rodnianski, I., Schlag, W.: Classical and quantum scattering for a class of long range random potentials. Int. Math. Res. Notes, no 5 (2003), 243-300.
13. Shaban, W., Vainberg, B.: Radiation conditions for the difference Schrödinger operators, Applicable Analysis, 80, 525-556 (2001).
14. Vainberg, B.: Asymptotic Methods in Equations of Mathematics Physics, Gordon and Breach Publishers, 1989.