\subsection{Estimates for entire functions.} An entire
function $f(z)$ is said to be of exponential type if there is a
constant $\b$ such that $|f(z)|\leq\const e^{\b |z|}$ everywhere.
The infimum of the set of $\b$ for which such inequality holds is
called the type of $f$.

\no {\bf Definition.}
{\it Let $\cE_\g, \g >0$ denote the set
 of exponential type  functions $f$, which satisfy
\[
\lb{D1} |f(k)|\ge 2|k|\qqq   \qqq \forall  \ k\in \R,
\]
\[
\lb{D2}
\begin{aligned}
&|f(k)-2ik+f_0|\le {Q^2\/|k|_1}e^{ \g (|\Im k|-\Im
k)+{Q\/|k|_1}},\qqq
\qq  \forall \ k\in \C,\\
& Q> 0,\qq  f_0\in \R,\qqq \qqq |f_0|\le Q,
\end{aligned}
\]
where the constants $Q=Q(f), f_0$ depend on $f$ and
$|k|_1=\max\{1,|k|\}$.}

In the proof of Theorem \ref{T1} we will need some properties of zeros of
$f\in \cE_\g $ in terms of the Carleson  measure. Recall that {\it a
positive Borel measure $M$ defined in $\C_-$ is called a
Carleson measure if there is
 a constant $C_M$ such that for all $(r,t)\in \R_+\ts\R$
\[
\lb{1.31}
   M(D_-(t,r))\leq C_Mr,\ \ \ {\rm where}\ \ \
   \ D_-(t,r)\ev\{z\in \C_-: |z-t|<r \},
\]
here $C_M$ is the Carleson  constant independent of $(t,r)$.}

For an entire function $f$ with zeroes $k_n, n\ge 1$
  we define an associated measure by
\[
\lb{MO}
 d\O (k,f)=\sum_{\Im k_n\le 0} \d (k-k_n+i)du dv,\qqq k=u+iv\in \C_-.
\]
We denote the number of zeros of function $f$ having modulus  $<r$ by $\cN (r,f)$,
each zero being counted according to its multiplicity.
In order to prove Theorem \ref{T1} we need

  \begin{theorem}
\lb{TEf}
 Let $f\in \cE_\g $ for some $\g >0$. Then

 i) The number of zeros of $f$ in the disk $\{|k|<r\},r>0, r_1=r+{1\/2}$
satisfies
\[
\lb{L1} \cN (r,f)\le {1\/\log 2}\rt({4r_1\g\/ \pi }+\log (1+{4r_1})+
{|f_0|+8Q\/1+4r_1}\rt).
\]

ii) $\O (\cdot,f)$ is a Carleson measure and for all
$(r,t)\in \R_+\ts \R $ satisfies
\[
\begin{aligned}
\lb{L2}
\O (D_-(t,r),f)\le\cN (r,f(t+\cdot))\le C(f)r ,\\
C(f)={12\/\log 4}\rt({\g \/\pi }+1+{|f_0| +3Q\/4}\rt).
\end{aligned}
\]
iii) Let $k_n, n\ge 1$ be all zeros of $f$. Then for each $p>1$ the following
estimate holds true:
\[
\lb{L3}
\sum_{\pm \Im k_n\le 0} {1\/|k_n-2i|^{p}}
\le   C Y_pC(f),
\]
where $C\le 2^5$ is the absolute  constant and  $Y_p$ is given by
\[
\lb{p1}
  Y_p=\int_\R
(1+x^2)^{-{p\/2}}dx=\sqrt \pi{\G({p-1\/2})\/\G({p\/2})} \qqq \forall \qq p>1.
\]

\end{theorem}
{\bf Proof.} i) Recall the Jensen formula (see p. 2 in \cite{Koo88})
for an entire function $F$:
\[
\lb{PJ}
\log |F(0)|+\int _0^r{\cN (t,F)\/ t}dt= {1\/ 2\pi }\int
_0^{2\pi}\log |F(re^{i\f})|d\f,
\]
 Rewrite $f$ in the form
 \[
 \lb{ef1}
 f=2ik-f_0+f_*, \qqq where \qq |f_*(k)|\le
{Q^2\/|k|_1}e^{{Q\/|k|_1}}e^{2\g v_-},\qqq v_-(k)=
{|\Im k|-\Im k\/2},
\]
for all $k\in \C$.
Consider the case $t\in \R\sm [-{1\/2}, {1\/2}]$. Let $z=t+k, k=re^{i\f}$.
Then \er{D2} and \er{ef1} imply
\[
\begin{aligned}
|f(t+k)|=|2it+2ik-f_0+f_*(z)|\le 2|t|\rt(1+{r\/|t|}\rt)\rt(1+{|f_0|+|f_*(z)|\/2(|t|+r)}\rt)\\
\le 2|t|\rt(1+{r\/|t|}\rt) e^{2\g v_-+{Q\/|z|_1}}\rt(1+{|f_0|+Q^2|z|_1^{-1}\/2(|t|+r)}\rt)\\
\le 2|t|\rt(1+{r\/|t|}\rt) e^{2\g v_-+{Q\/|z|_1}}
\rt(1+{|f_0|\/2(|t|+r)}\rt)  \rt(1+{Q^2\/2(|t|+r)|z|_1}\rt).
\end{aligned}
\]
Let $\cN (s)=\cN (s,f(t+\cdot))$.
Substituting this into \er{PJ} we obtain for the function $F(k)=f(t+k),
k=re^{i\f}, z=t+re^{i\f}$:
\[
\begin{aligned}
\lb{F1}
\log |f(t)|+\int_0^r{\cN (s)\/s}ds={1\/2\pi}\int_0^{2\pi}\log |f(z)|d\f\\
\le \log 2|t|+\log \rt(1+{r\/|t|}\rt)+\log \rt(1+{|f_0|\/2(|t|+r)}\rt)+{1\/2\pi}\int_0^{2\pi}2\g v_-d\f+X(t,r),\\
where \qq X(t,r)={1\/2\pi}\int_0^{2\pi}\rt[{Q\/|z|_1}+\log  \rt(1+{Q^2\/2(|t|+r)|z|_1}\rt) \rt]d\f,
\end{aligned}
\]
where $2v_-=r(|\sin \f| -\sin\f)$ and the simple integration yields
\[
\lb{F12}
\begin{aligned}
{1\/2\pi}\int_0^{2\pi}2\g v_-d\f=
{r\g\/2\pi }\int_0^{2\pi}\rt(|\sin \f|-\sin \f\rt)d\f={2r\g\/\pi}.
\end{aligned}
\]
Substituting the estimate ${|f(t)|\/2|t|}\ge 1$  into \er{F1}
together  with the simple one
$$
\int _0^r{\cN (s)ds\/s}\geq \cN ({r/2})\int _{r\/2}^r{ds\/s}=
\cN ({r/2})\log 2,
$$
we obtain
\[
\begin{aligned}
\lb{F3}
\cN ({r/2})\log 2\le {2r\g\/ \pi }+
\log \rt(1+{r\/|t|}\rt)+{|f_0|\/2(|t|+r)}+X(t,r).
\end{aligned}
\]

We show \er{L1}. Let $t={1\/2}$ and $r_1=r+{1\/2}$ for any $r>0$.
Then
$$
\{|k|<r\}\ss \{|k-{1\/2}|<r_1\},\qqq |z|_1=|(1/2)+k|_1\ge {1\/2}(1+|k|-(1/2))=
{1\/4}(1+2|k|)
$$
and \er{F3} give
\[
\lb{esr}
\begin{aligned}
\cN (r,f)\le \cN (r_1,f({1\/2}+\cdot))\le
{4r_1\g\/ \pi }+\log (1+4r_1)+{|f_0|\/(1+4r_1)}+X(t,r),\\
where \qq X(t,r)\le {1\/2\pi}\int_0^{2\pi}\rt[{4Q\/(1+4r_1)}+\log  \rt(1+{4Q^2\/(1+4r_1)^2}\rt) \rt]d\f\le {8Q\/(1+4r_1)},
\end{aligned}
\]
which yields  \er{L1}.

ii)
Let $r\le 1, t\in \R$. Then by the construction of $\O (\cdot,f)$, we obtain
$\O (D_-(t,r),f) =0$.

  In order to show \er{L2} for thew case  $r> 1, t\in \R$ we need to consider two cases:

Firstly, let  $r>1, t\in \R\sm [-{1\/2}, {1\/2}]$. Then due to \er{F3} and \er{F1},
the measure $\O (\cdot,f)$ satisfies
\[
\begin{aligned}
\lb{F4}
& \O (D_-(t,r),f)\le\cN (r,f(t+\cdot))\le
{1\/\log 2}\rt(   {4r\g\/ \pi }+4r+|f_0| +Q+\log (1+Q)^2\rt)\le C_1r,\\
& \qqq \qqq  C_1={4\/\log 2}\rt({\g\/\pi }+1+{|f_0| +3Q\/4}\rt).
\end{aligned}
\]

Secondly, let  $r>1, t\in [0, {1\/2}]$. The proof for the case $t\in [-{1\/2},0]$
is similar. For two disks $D(t,r)$ and $D({1\/2},r_1), r_1=r+{1\/2}$ we have
\[
\begin{aligned}
\lb{F5}
D(t,r)\cap \{\Im k\le -1\}\ss D\big({1\/2},r_1\big)\cap \{\Im k\le -1\},
\\
r\le r_1=r+{1\/2}\le {3r\/2}.
\end{aligned}
\]
Then due to \er{L1}, the measure $\O (\cdot,f)$ satisfies
\[
\begin{aligned}
\lb{F4a}
\O (D_-(t,r),f)\le \O (D_-(1/2,r_1),f)\le \cN (r_1,f({\small{1\/2}}+\cdot))\le
C_1r_1\le {3\/2}C_1r.
\end{aligned}
\]
Thus $\O (\cdot,f)$ is a Carleson measure with the Carleson constant $C(f)={3\/2}C_1$, where $C_1$ is given by \er{F4}.

iii)  Consider the case $\Im k_n\le 0$. The proof for the case $\Im k_n> 0$
is similar, even simpler.
 In order to show \er{L3} we recall the Carleson result (see p.
63, Theorem 3.9, \cite{G81}):

Let $F$ be analytic on $\C_-$. For $0<p<\iy$  we say $F\in
\mH_p=\mH_p(\C_-)$ if
$$
\sup_{y<0}\int_\R|F(x+iy)|^pdx=\|F\|_{\mH_p}^p<\iy
$$
Note that the definition of the Hardy space $\mH_p$ involve all
$y<0$, instead of small only value of $y$, like say, $y\in (-1,0)$.
We define the Hardy space $\mH_p$ for the case $\C_-$, since below
we work with functions on $\C_-$.

{\it If $M$ is a Carleson measure and satisfies \er{1.31}, then the
following estimate holds:
\[
\lb{Ce}
  \int_{\C_-} |F|^pdM\leq C C_M \|F\|_{\mH_p}^p\qqq

  \forall \qq F\in \mH_p, \ p\in (0,\iy ),
\]
where  $C_M$ is the so-called Carleson constant from \er{1.31} and  $C\le 2^5$ is an
absolute constant.}

 In order to prove \er{L3} we  take the functions $F(k)={1\/k-i}$.
Then estimates \er{Ce}, \er{L2}  yield
\[
\lb{3.19} \sum_{\Im k_n\le 0} {1\/|k_n-2i|^p}=\int_{\C_-}
|F(\l)|^pd\O \leq C C(f) \|F\|_{\mH_p}^p, \qq p\in (1,\iy),
\]
$C$ is an absolute constant and $C(f)$ is defined in \er{L2}. Here
we have the simple identity
\[
\lb{3.19a} \|F\|_{\mH_p}^p=\int_\R {dt\/|t-i|^p}=\int_\R
{dt\/(t^2+1)^{p\/2}}=Y_p.
\]
The function $Y$ is studied in Lemma \ref{Tp}. Combine \er{3.19} and
\er{3.19a} we obtain \er{L3}.
\hfill  \BBox

\

\subsection{Estimates of resonances.} We consider Case 1 and without loss of
generality let $\supp q\ss [0,\g]$. The Schr\"odinger equation
\[
\lb{e1}
  -f''+q(x)f=k^2f,\ \ \  \ \ \ k\in \C \sm \{0\},
\]
has unique solutions $\p_{\pm}(x,k)$  such that $\p_+(x,k )=e^{ixk
},\ \ x\ge \g$ and $\p_-(x,k)=e^{-ik x},\ \ x\leq 0$. Outside of the
support of $q$ any solutions of \er{e1} have to be linear
combinations of $e^{\pm ik x}$. Functions $\p_\pm(x,\cdot),
\p_\pm'(x,\cdot),  x\in \R$ are entire. We define the functions $a,
w, s$ by
\[
\begin{aligned}
\lb{H1a}
w(k)=2ika(k)=\{\p_-(\cdot,k), \p_+(\cdot,k )\},\qqq
  \ s(k)=\{\p_+(\cdot,k ),\p_-(\cdot,-k )\},
\end{aligned}
\]
and $\{f, g\}=fg'-f'g$ denotes the Wronskian.  The functions $w(k),
s(k)$ are entire and the following asymptotic
estimates hold true:
\[
\lb{a1}
 a(k )=1+O(k ^{-1})\qqq {\rm as} \ \ |k |\to \iy ,\ \  k\in\ol\C_+.
\]
The scattering matrix for the operators $H, H_0=-{d^2\/dx^2}$ is given by
\[
S(k)=\ma a(k )^{-1}& r_-(k )\\ r_+(k )&a(k )^{-1}\am ,
\ \ \ \  r_{\pm}={s(\mp k)\/ w(k )}, \ \ \ \  k\in\R,
\lb{1.7}
\]
where $1/a$ is the transmission coefficient and $r_{\pm}$ are the reflection
coefficients. The matrix $S(k), k\in \R$ is unitary, which yields
\[
\lb{ws}
|w(k)|^2=4k^2+|s(k)|^2,\ \  \forall \ k\in \R.
\]
The solution $ \p_+$ of \er{e1}  satisfies the following equation
\[
\lb{3.1}
 \p_+(x,k)=e^{ixk}-\int _x^\g{\sin [k(x-t)]\/k}q(t)
\p_+(t,k)dt \qqq \forall\  (x,k)\in [0, \g]\ts \C.
\]
It is well known that equation \er{3.1} has a unique solution.
 Due to \er{3.1} the function $y(x,k)=e^{-ikx}\p_+(x,k)$ satisfies
the integral equation
\[
\lb{y1}
y(x,k)=1+\int _x^\g G(t-x,k)q(t)y(t,k)dt \qqq  G(t,k)={\sin kt\/k}e^{ikt},
\]
for all $(x,k)\in [0,\g]\ts \C$. We have the standard iterations
\[
\lb{y2}
y(x,k)=1+\sum_{n\ge 1}y_n(x,k),\qqq y_n(x,k)=
\int _x^\g G(t-x,k)q(t)y_{n-1}(t,k)dt,
\qq y_0=1.
\]
We need some properties of the functions introduced above.

\begin{lemma}
\label{Tw}
 Let $q\in L^1(\R)$ and  $\supp q\ss [0,\g]$. Then the functions
 $w, s$ and $\p _{+}(x,\cdot),\p _{+}'(x,\cdot)$, $x\in\R$ are entire and we have
\[
\lb{yn}
|y_n(x,k)|\le {h^n\/n!}e^{2(\g-x)v_-}, \qqq v_-=\ca 0,&  k\in \ol\C_+\\
|\Im k|,  &  k\in \C_-\ac,
\]
for any $n\ge 1, (x,k)\in [0,\g]\ts \C$, where
\[
\lb{v-}
 h={Q\/|k|_1}, \qq Q=\max\{\|q\|,\|q\|_1\},\qq
                                                |k|_1=\max\{1,|k|\};
\]
where $\|q\|=\int_0^\g |q(t)|dt$ and $\|q\|_1=\int_0^\g t|q(t)|dt$, and
\[
\begin{aligned}
\lb{y}
|y(x,k)|\le e^{2(\g-x)v_-+h},\\
|y(x,k)-1|\le h e^{2(\g-x)v_-+h},\\
|y(x,k)-1-y_1|\le {h^2\/2} e^{2(\g-x)v_-+h}.
\end{aligned}
\]
Moreover, $w\in \cE_\g$ and satisfies for any $k\in \C$:
\[
\lb{w1}
w(k)=i2k-q_0 +w_*(k), \qqq w_*(k)=-\int _0^\g q(t)\big(y(t,k)-1\big)dt,
\]
\[
\lb{w2}
|w_*(k)|\le \|q\| he^{h+2\g v_-},
\]
where $q_0=\int_0^\g q(t)dt$.

\end{lemma}
{\bf Proof.} The function $G(t,k)={\sin kt\/k}e^{ikt}$ satisfy
\[
\lb{y3}
|G(t,k)|\le {e^{2 v_-t}\/|k|_1}\ca 1 & if \ |k|\ge 1\\
                                   t & if \ |k|<1 \ac,
 \qqq t\ge 0, \ k\in \C,
\]
Consider the first case in \er{y3}: $|k|\ge 1$, the proof for the second case
$|k|<1$ is similar.
Substituting the estimate in \er{y3} for $|k|\ge 1$ into the identity
$$
y_n(x,k)=\int\limits_{x=t_0<t_1< t_2<...< t_n}
\lt(\prod\limits_{1\le j\le n}
G(t_{j}-t_{j-1},k)q(t_j)\rt)dt_1dt_2...dt_n,
$$
we obtain
\[
\begin{aligned}
|y_n(x,k)|\le {1\/|k|_1^n}\int\limits_{x=t_0<t_1< t_2<...< t_n}
\lt(\prod\limits_{1\le j\le n}e^{2v_-(t_{j}-t_{j-1})}
|q(t_j)|\rt)dt_1dt_2...dt_n\\
= {e^{-2v_- x}\/|k|_1^n}\int\limits_{x=t_0<t_1< t_2<...< t_n}
\lt(\prod\limits_{1\le j\le n}
|q(t_j)|\rt) e^{2v_-t_{n}}dt_1dt_2...dt_n\\
\le {e^{2(\g-x)v_-}\/|k|_1^n}\int\limits_{x=t_0<t_1< t_2<...< t_n}
|q(t_1)q(t_2)....q(t_n)| dt_1dt_2...dt_n\le e^{2(\g-x)v_-}{\|q\|^n\/n! |k|_1^n}.
\end{aligned}
\]

This shows that for each $x\ge 0$
the series \er{y2} converges uniformly on any bounded subset of $\C$.
Each term of this series is an entire function. Hence the sum is
an entire function. Summing the majorants we obtain estimates  \er{y}.
Thus the functions  $w, s$ and $\p _{+}(x,\cdot),\p _{+}'(x,\cdot),x\in\R$
are entire

The identity \er{3.1} gives
\[
\begin{aligned}
\lb{y5}
\p_+(0,k)'=ik-\int _0^\g\cos [kt] \ q(t)
\p_+(t,k)dt, \\
ik\p_+(0,k)=ik+i\int _0^\g\sin [kt] \ q(t)\p_+(t,k)dt,
\end{aligned}
\]
which yields
\[
\begin{aligned}
\lb{y6}
w(k)=\p_-(0,k )\p_+'(0,k )-\p_-'(0,k )\p_+(0,k )
=ik\p_+(0,k ) +\p_+'(0,k )\\
=2ik-q_0+w_*(k), \qqq w_*(k)=-\int _0^\g q(t)\big(y(t,k)-1\big)dt.
\end{aligned}
\]
Due to \er{y} we have the estimate
\[
\begin{aligned}
\lb{y7}
|w_*(k)|\le \int _0^\g|q(t)|h e^{2(\g-t)v_-+h}dt\le
h e^{2\g v_-+h}\int _0^\g|q(t)|dt=h \|q\|e^{2\g v_-+h},
\end{aligned}
\]
which is \er{w2}. Together with \er{ws} this implies that $w\in \cE_\g$.
\hfill \BBox

\

The identity \er{w1} and estimate \er{w2} give
 \[
\lb{w1x}
\begin{aligned}
w(k)=i2k-q_0 +w_*(k), \qq |w_*(k)|\le {Q^2\/|k|_1}e^{{Q\/|k|_1}+2\g v_-},\qq
 \forall \ k\in \C,\\
 where \ \ q_0=\int_0^\g q(t)dt,\qqq Q=\max\{\|q\|,\|q\|_1\}.
 \end{aligned}
\]

\

We show the relations between the Cases 1, 2 and 3. In order to show
this we use a standard trick. For $q\in L^1(\R_+)$ with $\supp q\ss
[0,\g]$ we define the Schr\"odinger operator $\wt
H$ acting in $L^2(\R )$ with an even compactly
supported potential $\wt q$ given by
\[
\lb{wtq} \wt H=-{d^2\/dx^2}+\wt q,\qqq \qqq \wt q(x)=q(|x|), x\in \R.
\]
 The Schr\"odinger equation
$
  -f''+\wt q(x)f=k^2f,\ \ \  \ \ \ k\in \C \sm \{0\},
$ has unique solutions $\wt\p_{\pm}(x,k)$  such that $\wt\p_+(x,k
)=e^{ixk },\ \ x\ge \g$ and $\wt\p_-(x,k)=e^{-ik x},\ \ x\leq -\g$.
Note that the symmetry of the potential $\wt q$ yields
$$
\wt\p_+(x,k)=\p_+(x,k)=\wt\p_-(-x,k) \qqq \forall x\in [0,\g].
$$
This implies that the Wronskian $\wt w(k)$ for the potential $\wt q$
satisfies
\[
\lb{tw} \wt w(k)=\{\wt\p_+(x,k),
\wt\p_-(x,k)\}|_{x=0}=2\p_+(0,k)\p_+'(0,k).
\]
where $\p_+(0,k)$ and $\p_+'(0,k)$ are the Jost functions for the
Cases 2 and 3 respectively with the potential $q$. This identity is
very useful. For example, if we have estimate \er{3} for the Case 1,
then  \er{tw} gives the estimate \er{3} for the Case 2 and 3.

We obtain the estimate of the Wronskian $\wt w(k)$.
We need to define the potential
$q_1=\wt q(x-\g)$, which has the support $\supp q_1\ss [0,2\g]$.
Using Lemma \ref{Tw} and \er{w1}, \er{w2} we obtain
\[
\lb{wtw}
\begin{aligned}
\wt w(k)=i2k-2q_0 +\wt w_*(k), \\
|\wt w_*(k)|\le \|q_1\| {\wt Q^2\/|k|_1}e^{4\g v_-+{\wt Q\/|k|_1}},
\end{aligned}
\]
where the corresponding constant $\wt Q$ is given by
\[
\begin{aligned}
\lb{wtQ}
\wt Q=\max \{\|q_1\|, \|q_1\|_1\}=2\|q\|\max \{1, \g\},
\end{aligned}
\]
since
$$
\|q_1\|=\int_0^{2\g} |q_1(x)|dx=2\|q\|,\qqq \|q_1\|_1=\int_0^{2\g} x|\wt q(x-\g)|dx=
\int_{-\g}^{\g} (t+\g)|\wt q(t)|dt=2\g \|q\|.
$$

We describe few facts about the resonances in the disk $\{|k|<r\}$.

\begin{proposition}
\lb{Tr}
 Let $q\in L^1(\R)$ and  $\supp q\ss [0,\g]$ for some $\g>0$.

 i)  Let $r>0$ and  let $q$ satisfy
\[
\lb{q0}
\|q\|(1+he^{2\g r+h})<2r, \qqq {\rm where} \ h={\max \{\|q\|,\|q\|_1\}\/\max\{1,r\}}.
\]
Then the function $w$ has only one simple zero in the disk $\{|k|<r\}$.

ii) Let $1\le r\le {1\/2\g}$ and let $2\|q\|\le r$.
Then the function $w$ has only one simple zero in the disk $\{|k|<r\}$.

\end{proposition}
\no {\bf Proof}. i) Let $w_0=i2k$ and $|k|=r$.  Then the estimates \er{w1},
\er{w2} imply
\[
|w(k)-w_0(k)|\le \|q\| (1+he^{2\g+h})=|w_0(k)|C_0,\qqq
C_0={\|q\|\/2r}(1+he^{2\g r+h}).
\]
Hence, if $C_0<1$, then by Rouche's theorem, $w(k)$ has only one simple zero,
as $w_0=2ik$ in the disk $\{|k|<r\}$.

ii) If we assume that $f(h)=h+h^2e^{1+h}<2$, then \er{q0} holds true.

We have that the function $f$ is increasing. We get
the estimates $f({1\/2})<2$, since  $h={\|q\|\/r}\le {1\/2}$, which yield that $w$
has only one simple zero in the disk $\{|k|<r\}$.
\hfill\BBox

{\bf Remark.} 1) This proposition shows that for  any $\g, r>0$ and sufficiently
small $\|q\|$  the function $w$ has only one simple zero in the disk $\{|k|<r\}$.

2)  The number of zeros of $w$ in the disk with large $r$  depends on the diameter
$\g$. But it is very interesting to determine the asymptotics of $\cN(r,w)$ as
$r\to \iy$. This problem remains open for  long period.

3) We need to say that \er{q2} holds also true for the Case 2 and 3.
 Then \er{3} for Case 1 yields \er{3} for Cases 2 and 3.

\

 {\bf Proof of Theorem \ref{TrJS}}.
Lemma \ref{Tw} gives that the function $w\in \cE_\g$.

We apply the estimate \er{L1} to the function $w$. Thus  due to the estimate \er{w1x} and the simple fact $|q_0|\le Q$,
we obtain \er{q2}.

We show \er{q3}.
The identity \er{tw} gives $\cN (r,\p_+(0,\cdot))+\cN
(r,\p_+'(0,\cdot))=\cN (r,\wt w)$, where $\wt w$ is the Wronskian
for the potential $\wt q$ given by \er{wtq} and $\wt w$ satisfies \er{wtw}, \er{wtQ}. Then these arguments and  the estimate \er{q2} implies
 \er{q3}. \hfill $\BBox$

 \medskip

{\bf Proof of Theorem \ref{T1}.} We consider Case 1 and let $\supp
q\ss [0,\g]$. By Lemma \ref{Tw}, the function $w$ belongs to
$\cE_\g$ with $Q=\max\{\|q\|,\|q\|_1\}$ and satisfies \er{w1x}.
 Thus  the estimate \er{L3}, \er{L2}  from  Theorem
\ref{TEf} give \er{3} for  Case 1.

Consider Cases 2 and 3.  Then we have the identity $\wt w(k)=2\p_+(0,k)\p_+'(0,k)$,
where $\wt w$ is the Wronskian for the potential $\wt q$ given by \er{wtq}
(see  \er{tw}). Then  \er{3} for $\wt q$ (Case 1) together with  \er{wtQ} implies
\er{3} for Cases 2 and 3.
\hfill $\BBox$

\begin{lemma}
\lb{Tp} The function $Y_p=\int_\R {(1+x^2)^{-{p\/2}}}dx, p>1$ is
strongly monotonic and convex on $(1,\iy)$, since $Y_p'<0, Y_p''>0$,
and  satisfies \er{Yp}  and the identity \er{p1}, i.e., $Y_p=\sqrt \pi{\G({p-1\/2})/\G({p\/2})}$ for all $p>1$.
\end{lemma}
\no {\bf Proof.} We rewrite $Y_p$ in the form $Y_p=\int_\R e^{-p
f(x)}dx$, where $f(x)={1\/2}\log (1+x^2)$. This yields
$Y_p'=-\int_\R f(x)e^{-p f(x)}dx<0$ and $Y_p''=\int_\R f^2(x)e^{-p
f(x)}dx>0$. We rewrite $Y_p$ in another form
$$
Y_p=2\int_0^\iy {(1+x^2)^{-{p\/2}}}dx=\int_0^\iy
{t^{-{1\/2}}(1+t)^{-{p\/2}}}dt=\sqrt \pi{ \G(z)\/\G(z+{1\/2})}, \qq
z={p-1\/2},
$$
where $\G$ is the Gamma function, see identities (2), (5) in Section
1.5 in \cite{EMOT53}.

Asymptotics  ${ \G(z)\/\G(z+{1\/2})}=z^{-{1\/2}}(1+O({1\/z}))$ as
$z\to \iy$, see  (5) in Sect. 1.18 in \cite{EMOT53} give \er{Yp} as
$p\to \iy$. Recall that $\G$ is meromorphic in $\C$ with simple
poles $0,-1,-2,..$. Then the identities $\G(1+z)=z\G(z)$ and
$\G(1)=1, \G(1/2)=\sqrt \pi$ yield  \er{Yp} as $p\to 1$:
$$
Y_p=\sqrt \pi{ \G(z)\/\G(z+{1\/2})}={\sqrt\pi \G(1+z)\/
z\G(z+{1\/2})}={1+O(z)\/z}\qq \as \  \ z={p-1\/2}\to 0,
$$
 \hfill $\BBox$

 \medskip

\no {\bf Acknowledgments.} \footnotesize I am grateful to  Michael
Loss, Atlanta, Grigori Rozenblum, G\"oteborgs, and Sergey Morozov,
Munich, for useful comments about Lieb-Thirring inequalities and
Elliot Lieb, Princeton, for the reference \cite{DFLP06}. Moreover, I
am also grateful to  Alexei Alexandrov  (St. Petersburg) for
stimulating discussions about the Carleson Theorem and about the
absolute constant $C$ in \er{Ce}.

This work was supported by the Ministry of education and science of
the Russian Federation, state contract 14.740.11.0581 and the RFFI
grant "Spectral and asymptotic methods for studying of the
differential operators" No 11-01-00458.


\section{Introduction}\label{sec:1}

Chaotic motion accounts on the one hand for the well-known phenomenon of sensitive dependence on initial conditions, that is, exponentially
fast divergence of nearby orbits, and on the other hand for
the phenomenon of decay of correlations or mixing.
Both properties are intimately related with the observation
that even low-dimensional chaotic systems share common features
with random processes.

This intuitive picture
has been used as a basis to address some of the fundamental
questions arising in nonequilibrium statistical mechanics
\cite{Kryl:79,Kubo_Sci86,FaIsVu_PLA90}.
In fact,
it is a simple exercise to show that topological mixing implies
sensitive dependence on initial conditions. At the  measure-theoretical
level, however, relating Lyapunov exponents, the quantitative measures for sensitive dependence
on initial conditions,
to decay rates of correlation functions is a more involved task.
For instance, it is easy to construct
simple maps with finite Lyapunov exponent and arbitrarily small correlation decay
(see, for example, a Markov model described in \cite{Just_PLA90}).

Thus, at a quantitative level it is tempting to explore in some detail in which way
the rate of correlation decay is linked with Lyapunov exponents, as both
quantities are supposed to have a
common origin. In more general terms and from a wider perspective
this topic can be viewed as belonging to the realm of fluctuation dissipation relations for
nonequilibrium dynamics, where one cause, the underlying detailed dynamical structure, is
responsible for the approach to the stationary state, that is, for the decay of
correlations, but, at the same time, is responsible for fluctuation properties
at a microscopic level, or in our context for the sensitive dependence of initial conditions and
a positive Lyapunov exponent. Furthermore, any relation
between decay of correlations and Lyapunov exponents is of great practical interest,
as the measurement of Lyapunov exponents, unlike the
rate of correlation decay, is notoriously difficult to determine in real world experiments \cite{Abar:96,KaSc:97}.
In \cite{Spro:96} it was even suggested to take correlation decay rates as
a meaningful approximation for Lyapunov exponents.

The problem we want to address can be illustrated by a very basic textbook example,
probably considered for the first time more than two decades ago
\cite{BaHeMePo_PRA88}.
Consider a linear full branch map (see \Fref{fig:1}) on the unit interval $I=[0,1]$, that is,
a map $f:I\rightarrow I$ having a finite partition of $I$ into closed intervals $I_k$ with pairwise
disjoint interior such that
(i) for each $k$, we have $f(I_k)=I$, and
(ii) $f$ has constant slope $\gamma_k$ on each $I_k$.
The physical invariant measure is given by the Lebesgue measure and the
Lyapunov exponent
with respect to this measure
can be expressed in terms of the slopes
\begin{equation}\label{aa}
\Lambda = \sum_k |I_k| \ln|\gamma_k|,
\end{equation}
with $|I_k|=1/|\gamma_k|$ denoting the size of the interval $I_k$.

The exponential rate of decay for correlation
functions
is determined by the negative logarithm of the second largest eigenvalue in modulus of the associated Perron-Frobenius operator.
In this setting it is well known (see, for example, \cite{MoSoOs_PTP81})
that eigenfunctions of the
Perron-Frobenius operator are given by polynomials and
that the corresponding eigenvalues $\nu_m$ can be expressed as
\begin{equation*}\label{ab}
\nu_m = \sum_k \frac{1}{|\gamma_k|} \frac{1}{\gamma_k^{m}} =
\sum_k |I_k| \frac{1}{\gamma_k^{m}} \,\qquad (m \geq 0),
\end{equation*}
with largest eigenvalue $\lambda_0=\nu_0=1$ and the subleading eigenvalue $\lambda_1$
being the second largest in modulus, $|\lambda_1| = \max\{|\nu_1|,\nu_2\}$.
Thus, correlation functions of sufficiently smooth observables decay typically at an
exponential rate $\alpha=-\ln|\lambda_1|$.
Since $|\lambda_1|\geq \nu_2>0$ we obtain an upper bound for the decay rate
$\alpha \leq -\ln \nu_2$ which can now be related to the Lyapunov exponent
\eqref{aa}. If we apply Jensen's inequality to the convex function
$\varphi(x)=-\ln (x)$ we end up with
\begin{equation}\label{ac}
\alpha \leq -\ln \sum_k |I_k| \frac{1}{\gamma_k^{2}} \leq
\sum_k |I_k| (-\ln \frac{1}{\gamma_k^{2}}) = 2 \Lambda \, .
\end{equation}
The estimate of the decay rate in \Eref{ac} has been based on $\nu_2$ which contains
positive terms only, even if the slopes have different signs. As a result the upper
bound is given by twice the Lyapunov exponent. If all slopes have the same sign,
say $\gamma_k>1$, then $\nu_1$ determines the subleading eigenvalue and the Lyapunov
exponent itself yields an upper bound for the decay rate, that is, $\alpha \leq \Lambda$.

\begin{figure}[h!]
  \centering
    \includegraphics[width=.30\textwidth]{fullbranch.eps}
  \caption{Diagrammatic view of a linear full branch map. \label{fig:1}}
\end{figure}

In this article we will address the question
to which extent this simple reasoning can be generalised to a larger class of systems.
From a mathematical perspective there is a considerable amount of
literature on the existence of invariant measures of one-dimensional maps,
the corresponding Lyapunov exponent, estimates for the spectra of the
associated
Perron-Frobenius  operator, and a possible relation between
Lyapunov exponents and the decay of correlation functions \cite{CoEc_JSP04}. Expansiveness of the
underlying map $f:I\to I$
is to some extent the key ingredient to establish a
physical measure $\mu$, that is, an invariant measure which asymptotically characterises a large set of orbits, or in formal terms an invariant ergodic probability measure which is absolutely continuous with respect to Lebesque measure (see, for example,
\cite{Live_AM95,Kell_IJBC99,Bala:00,BrLuSt_ASENS03}).

The Lyapunov exponent for a Lebesque typical point $x\in I$ is given by
\begin{equation}\label{eq:Lyapunov}
\Lambda = \int_I \ln|f'(x)| d\mu
= \lim_{n\to \infty} \frac{1}{n} \ln{|(f^n)^{\prime}(x)|.}
\end{equation}
Regarding decay rates for correlation functions, a substantial part of the problem
consists in giving the notion of the rate of correlation decay a proper meaning
(see, for example, \cite{CoEc_JSP04}). With the standard definition of a
correlation
function for two observables $\varphi$ and $\psi$,
\begin{equation}\label{CorrelFunc}
C_{\varphi,\psi}(n) = \int_I{\varphi(f^n(x)) \psi(x) d\mu} - \int_I{\varphi(x) d\mu}
\int_I{\psi(x) d\mu},
\end{equation}
the exponential decay rate $\alpha_{\varphi,\psi}$
governing the asymptotic behaviour of the correlation function may be
formally introduced by
\begin{equation*}\label{decaycorr}
\alpha_{\varphi,\psi}= \sup\{s:\limsup_{n\rightarrow \infty} |\exp(sn) C_{\varphi,\psi}(n)|<\infty\} \, .
\end{equation*}
The rate $\alpha_{\varphi,\psi}$ is sensitive to the choice of $\varphi$ and
$\psi$ and can be made arbitrarily small or large by special choices for the observables
\cite{CraCa_PHYD83}. However, it is possible to define a rate of decay with respect  to ``typical'' observables from  some linear space~$\mathcal{H}$
\begin{equation*}\label{decayrate}
\alpha_{\mathcal{H}} = \inf\{\alpha_{\varphi,\psi} : \varphi,\psi \in {\mathcal{H}}\} \, .
\end{equation*}
We will call this quantity the
\oemph{mixing rate}. Note that this rate still depends,
for instance, on the degree of smoothness shared by the observables in ${\mathcal H}$.

One approach for determining bounds on the  mixing rate $\alpha_{\mathcal{H}}$ relies on reformulating
\eref{CorrelFunc} in terms of the Perron-Frobenius operator.
The rate $\alpha_{\mathcal{H}}$
is then determined by the subleading eigenvalue or by the essential spectral radius
of the Perron-Frobenius operator. The desired relation between the Lyapunov exponent and the
mixing rate expressed in \Eref{ac} then becomes a lower bound for the subleading eigenvalue.
There is a considerable body of literature on upper bounds
for spectral values
(see, for example, \cite{Live_AM95,Kell_IJBC99,BrLuSt_ASENS03}) providing useful tools to establish ergodic properties
of dynamical systems.
However, to the best of our knowledge hardly any nontrivial lower bounds exist
(see, however, \cite{Naud_AnnHP09},  where an exponential lower bound for correlation functions of suspension semiflows is given).

The main tool to establish the desired inequality consists in studying the properties
of the generalised Perron-Frobenius operator $\mathcal{L}_{\beta}$ with potential $-\beta \ln|f^{\prime}|,$
\begin{equation}\label{eq:RPF}
 (\mathcal{L}_{\beta} h)(x) = \sum_{y \in f^{-1}(x)} \frac{h(y)}{|f^{\prime}(y)|^{\beta}} \, .
\end{equation}
For $\beta=1$ this expression reduces to the Perron-Frobenius operator which has leading
eigenvalue one.
The subleading eigenvalue $\lambda_1$  (or the essential spectral radius, if no subleading
eigenvalue exists) determines the mixing rate $\alpha_{\mathcal{H}}=-\ln|\lambda_1|$.
In addition, the derivative of the largest eigenvalue $\nu_0(\beta)$ of \Eref{eq:RPF}
with respect to $\beta$, that is, the derivative of the topological pressure,
determines the Lyapunov exponent
(see, for example, \cite{BoRa_PD87} for a basic exposition).
Finally, the required estimate follows from the convexity
of the topological pressure and gives a lower bound for the subleading eigenvalue of
the Perron-Frobenius operator.

The main challenge is, of course, to put these ideas into
practice. Restricting to piecewise linear Markov maps considerably reduces
the need to worry about mathematical subtleties, as the operator $\mathcal{L}_\beta$ admits finite-dimensional
matrix representations when considering observables consisting of piecewise analytic functions.
Thus, at a computational level all technical details reduce to
straightforward matrix manipulations (see \Sref{sec:2}), allowing us to keep the presentation
elementary and, at the same time, making the underlying ideas transparent.
To keep the presentation self-contained  basic properties of piecewise linear
Markov maps  together with properties of $\mathcal{L}_\beta$ are summarised in the Appendix.
In \Sref{sec:3} we are addressing the question whether the bounds presented in \Sref{sec:2}
hold for general (non-linear) expanding Markov maps. We have compelling evidence
that the estimate breaks down if analytic observables are considered
(see also \cite{ChPaRu_PRL90}), but
the validity of our proposition can be restored if observables of bounded variation are
considered.

\section{Piecewise linear Markov maps}\label{sec:2}

Let us consider a topologically mixing piecewise linear expanding Markov map $f:I\rightarrow I$
with respect to the partition $\mathcal{P}=\{I_1,\ldots I_N\}$ of the interval $I$
(see Definition~\ref{defn:MarkovMap} for a formal account).
Denote by $f_k:=f_{I_{k}}$ the $k$-th branch of $f$.
A nice property of
such maps is that the corresponding transfer operators have finite matrix
representations. In fact, the space of piecewise polynomial functions of degree less than $M$, that is, the space of functions of the form
$x\mapsto \sum_{k=1}^N \sum_{m=0}^M a_{k m} x^m \chi_{k}(x)$, where $\chi_k$ is the characteristic
function of the interval $I_k$, is an $\mathcal{L}_{\beta}$-invariant subspace.
While it appears to be difficult to trace the earliest reference
for this result,
to the best of our knowledge the first
experimentally relevant application appeared in the context of
power spectra for intermittent dynamics (see \cite{MoSoOs_PTP81,SoYoOkMo_JSP84}).

Now, in the natural basis
of piecewise monomials, the
operator $\mathcal{L}_{\beta}$ restricted to the invariant subspace mentioned above
is represented
by the $(M+1)^N\times (M+1)^N $  block upper triangular matrix

\begin{equation}\label{blockmatrix}
\left(
\begin{array}{cccc}
T^{(00)}(\beta)
 & \quad T^{(01)}(\beta) & \quad \ldots &
T^{(0M)}(\beta) \\
0 & \quad T^{(11)}(\beta) & \quad  \ddots & \vdots \\
\vdots &  \ddots & \quad \ddots & \quad T^{(M-1 M)}(\beta) \\
0 & \cdots & 0 & \quad T^{(M M)}(\beta)
\end{array}
\right) \, .
\end{equation}

A calculation similar to the one used to obtain the matrix representation $T^{(0 0)}(1)$ of $\mathcal{L}_1$ on the space of piecewise constant functions
(see, for example, \cite[p.~176]{BoyaGora:97}),
shows that the matrix elements of the block matrices $T^{(mn)}(\beta)$ are given in terms of the
slopes $\gamma_k$, the intercepts $d_k$ of the branches $f_k$ of the map $f$,
 and the topological transition matrix $(A_{kl})_{1 \leq k,l\leq N}$ induced by $f$ and $\mathcal{P}$
(see \Eref{appeqa}) as follows:
\begin{equation}\label{melements}
T_{kl}^{(mn)}(\beta)=
\frac{A_{lk}}{|\gamma_l|^\beta \gamma_l^n}\cdot
 (-d_l)^{n-m}
{n \choose n-m} \, .
\end{equation}
The eigenvalues are determined by the diagonal blocks
$T^{(mm)}(\beta)$ with matrix elements given by
the first factor in \Eref{melements}.

Note that since the underlying map is assumed to be topologically mixing, the matrices $T^{(mm)}(\beta)$ are
irreducible and aperiodic for even $m$. The Perron-Frobenius Theorem
(see, for example, \cite[p.53]{Gantmacher:2000} or \cite[p.536]{Lancaster:69}),
now guarantees that $T^{(mm)}(\beta)$ has a simple, positive eigenvalue, larger (in modulus) than all other eigenvalues, which we denote by $\nu_m(\beta)$.
Now, $\nu_0(\beta)$ determines the topological pressure, given by
$P(\beta)=\ln\nu_0(\beta)$, which has the following well-known properties.

\begin{lem}\label{prop:TopPres}\noindent
\begin{itemize}
\item[(i)] $P(1) = 0$;
\item[(ii)]  $P(\beta)$  is a convex function of $\beta$;
\item[(iii)] $P(\beta)$ is analytic in $\beta$;
\item[(iv)] $P'(1)=-\Lambda$.
\end{itemize}
\end{lem}

These properties have been established for a rather large class
of dynamical systems (see, for example, \cite{Kell:98}). In our context, however, they follow easily
from the matrix representation above using elementary methods. As an example, we shall derive property (ii).
To simplify notation, we shall denote the first block matrix $T^{(00)}(\beta)$ by $T(\beta)$
for the remainder of the argument.
Observing that $P(\beta)=\ln \nu_0(\beta)$ is obtained from
$\ln \mbox{Tr} \left((T(\beta))^n\right)/n$ as
$n$ tends to infinity, convexity of the topological pressure, that is,
\[
P(\beta_1 (1-t)+ \beta_2 t) \leq
(1-t) P(\beta_1) + t P(\beta_2) \qquad (0\leq t\leq 1),
\]
follows from a simple estimate using the H\"older inequality
\begin{eqnarray*}
 & &\frac{1}{n} \ln  \mbox{Tr} \left(
(T(\beta_1(1-t)+\beta_2 t))^n \right)\nonumber\\
 &=& \frac{1}{n} \ln \sum_{l_0,l_1,\ldots, l_{n-1}}
T_{l_0 l _1}(\beta_1 (1-t) + \beta_2 t)
T_{l_1 l _2}(\beta_1 (1-t) + \beta_2 t)
\cdots
T_{l_{n-1} l _0} (\beta_1 (1-t) + \beta_2 t)
\nonumber \\
 &=&  \frac{1}{n} \ln \sum_{l_0,l_1,\ldots, l_{n-1}}
\left(T_{l_0 l _1}(\beta_1)\right)^{1-t}
\left(T_{l_0 l _1}(\beta_2)\right)^{t}
\cdots
\left( T_{l_{n-1} l _0}(\beta_1)\right)^{1-t}
\left( T_{l_{n-1} l _0}(\beta_2)\right)^{t} \nonumber \\

&\leq& \frac{1}{n} \ln \left( \sum_{l_0,l_1,\ldots, l_{n-1}}
T_{l_0 l _1}(\beta_1)
\cdots
T_{l_{n-1} l _0}(\beta_1) \right)^{1-t}
\left(\sum_{l_0,l_1,\ldots l_{n-1}}
T_{l_0 l _1}(\beta_2)
\cdots
T_{l_{n-1} l _0}(\beta_2) \right)^t \nonumber \\

&=& (1-t) \frac{1}{n} \ln \mbox{Tr} \left((T(\beta_1))^n\right)
+ t \frac{1}{n}\ln \mbox{Tr} \left( (T(\beta_2))^n \right)\, .
\end{eqnarray*}
Here we have used that \Eref{melements} implies
$T_{kl}(\beta_1 +\beta_2)=
T_{kl}(\beta_1)T_{kl}(\beta_2)$ as well as
$T_{kl}(t \beta)=(T_{kl}(\beta))^t$.
The other statements of the
lemma can be derived in a similar way, or can be found in standard textbooks
on ergodic theory.

To establish a relation between the correlation decay, that is,
between the eigenvalues of the operator \eref{eq:RPF} for $\beta=1$, and the
Lyapunov exponent, note that the largest eigenvalue of the Perron-Frobenius operator
is given by $\nu_0(1)=1$, while
$\nu_2(1)$ is a positive eigenvalue, which provides a lower bound for the
subleading eigenvalue of $\mathcal{L}_{\beta}$.
Thus, on the one hand
\begin{equation}\label{ba}
\alpha_\mathcal{H} \leq -\ln \nu_2(1) \, .
\end{equation}
On the other hand $T^{(22)}(\beta)=T^{(00)}(\beta+2)$ by \Eref{melements}, which implies
\begin{equation}\label{bb}
\nu_2(\beta)=\nu_0(\beta+2) \, .
\end{equation}
Hence,  using the properties of the topological pressure
in Lemma~\ref{prop:TopPres}, the relations \eref{ba} and \eref{bb}
yield
\begin{equation}\label{eq:aH_lambda3}
\alpha_\mathcal{H} \leq -P(3)\leq (3-1) \Lambda \, .
\end{equation}
See \Fref{fig:TopPressure} for a graphical illustration of this result.
\begin{figure}[h!]
  \centering
    \includegraphics[angle=0,width=.70\textwidth]{topp.eps}
  \caption{Schematic representation of the topological pressure $P(\beta)$ and
graphical illustration of the estimates (\ref{eq:aH_lambda3}) and (\ref{eq:aH_lambda2}).
\label{fig:TopPressure}}
\end{figure}

Note that if all slopes $\gamma_k$ of $f$ have the same sign, then we have
$T^{(11)}(\beta)
=\mbox{sign}(\gamma_k)T^{(00)}(\beta+1)$. Thus,
we can apply the previous arguments to $|\nu_1(1)|$
to obtain the following improved estimate
\begin{equation}\label{eq:aH_lambda2}
\alpha_\mathcal{H} \leq  -\ln |\nu_1(1)| =  -\ln \nu_0(2) =
-P(2)\leq (2-1) \Lambda \, .
\end{equation}
To summarise, we have shown the following.

\begin{prop}\label{thm:main}
Let $f: I \rightarrow I$ denote a topologically mixing
piecewise linear expanding Markov map.
If $\mathcal{H} = \mathcal{H}(D_R)$ denotes a space of piecewise analytic observables
(see Appendix, Definition~\ref{def:Hspace}) then
the mixing rate is bounded in terms of the Lyapunov exponent $\Lambda$ with respect to the piecewise constant invariant density, by
\begin{equation}\label{ineq}
\alpha_{\mathcal{H}} \leq 2 \Lambda \, .
\end{equation}
If all slopes have the same sign the sharper estimate
\begin{equation}\label{sharp}
\alpha_{\mathcal{H}} \leq  \Lambda
\end{equation}
holds.
\end{prop}

\begin{rem}
The assumption that $f$ be topologically mixing is sufficient but not
necessary. Indeed, there exist piecewise linear expanding Markov
maps $f$ with the following properties: the map $f$ is not topologically
mixing, yet exhibits exponential decay of correlations and
the conclusions of Proposition \ref{thm:main} hold.
\end{rem}

The simplest examples for which the bounds are achieved are
the tent map ($\alpha_{\mathcal{H}}=2\Lambda$) and the doubling map
($\alpha_{\mathcal{H}}=\Lambda$).
It is finally worth mentioning that the condition $P(1)=0$, that is, the absolute continuity
of the reference measure with respect to Lebesgue measure, is not essential
for Proposition~\ref{thm:main} to be valid. The conclusions also hold, for instance, for certain Gibbs measures with respect to piecewise
constant potentials.

\section{Remarks on general expanding Markov maps}\label{sec:3}

The setup of piecewise linear Markov maps is rather special. One may thus be tempted
to ask whether a result like Proposition~\ref{thm:main} extends, say, to Markov maps
with finite curvature. The previous considerations are based
on \Eref{bb}, a result which has been exploited previously in a slightly more
restricted setup \cite{BaHeMePo_PRA88} to conjecture an exact relation between
correlation decay and generalised Lyapunov exponents. It has already been
demonstrated that such an identity breaks down for maps with finite curvature
\cite{ChPaRu_PRL90}. This fact, taken on its own, however, does not prevent the validity of a generalisation of
Proposition~\ref{thm:main}, making a study of how curvature
affects the previous considerations a worthy task.

Using an Ulam-like construction any expanding Markov map can be approximated
by a piecewise linear map (see, for example, \cite{BaIsSc_AIHP95}).
Let $F:I\rightarrow I$ denote a piecewise expanding, but
not necessarily linear Markov map with Markov
partition $\{I_1,\ldots,I_N\}$. Using the cylinder sets $U_{i_0,\ldots, i_{n-1}}
= \cap_{k=0}^{n-1} F^{-k}(I_{i_k})$ we may introduce a piecewise linear approximation
$f_n: I\rightarrow I$ by the following construction. The map $f_n$ linearly interpolates
$F$ on each cylinder set, that is, $f_n(U_{i_0,\ldots, i_{n-1}}) = U_{i_1,\ldots, i_{n-1}}$
(see \Fref{fig:Counterexample}(a)).  It is a straightforward exercise
to show that the Lyapunov exponent of $f_n$ tends to the Lyapunov exponent of $F$ (with respect to the absolutely continuous invariant measure) as $n$ tends to infinity.
However, the analysis of the mixing rate requires greater care.
While Proposition~\ref{thm:main} is still valid for any order $n$ of the approximation, it is far from obvious whether the proposition
is valid for a general expanding map.

To illustrate this point we consider a simple example,
a family of full branch piecewise M\"obius maps $F_c$ defined  on $[-1,1],$
\begin{equation}\label{oscarmap}
F_c(x)=\frac{1-2(c+1)|x|}{1+2c|x|} \, .
\end{equation}
We restrict the parameter to $c \in (-1/4,1/2)$, in order to guarantee expansivity.
\Fref{fig:Counterexample}(a) depicts the map $F_c$ for $c=-0.22$.

The leading part of the spectrum of the corresponding Perron-Frobenius
operator considered
on the space of analytic observables
can be approximated using a spectral approximation method. The basic idea of this method is to approximate
$\mathcal{L}_{\beta}$ by an $n\times n$ square matrix $\Pi_n\mathcal{L}_{\beta}\Pi_n$,
where $\Pi_n$ denotes the projector that sends a function to its Lagrange-Chebyshev interpolating polynomial of degree $n-1$.
This method is easily implemented and, moreover, it is possible to show (see \cite{Oscar}) that
the eigenvalues of $\Pi_n\mathcal{L}_{\beta}\Pi_n$ converge exponentially fast to the
eigenvalues of $\mathcal{L}_{\beta}$.
Using this method the leading eigenvalues
of the Perron-Frobenius operator and their dependence on $c$ are
easily obtained (see \Fref{fig:Counterexample}(b)).
A minimum for the subleading eigenvalue occurs at about
$c=-0.11$.
The corresponding numerical value
reads $\lambda_1 \approx 0.10415 $ resulting in a mixing rate
$\alpha_{\mathcal{H}}=-\ln |\lambda_1| \approx 2.2619$. The corresponding Lyapunov
exponent (which hardly depends on the parameter $c$) is computed using
the numerical approximation
of the invariant density. The numerical value is
$\Lambda \approx 0.685$ so that the inequality \eref{ineq} is clearly violated.

\begin{figure}[h!]

$
\hspace{-.08\textwidth}
\begin{array}{cc}
\includegraphics[width=.65\textwidth]{fig3a.eps} &
\hspace{-.12\paperwidth}
\includegraphics[width=.65\textwidth]{fig3b.eps}
\end{array}$
\caption{Left: (a) Diagrammatic view of the  M\"obius map \eref{oscarmap} for $c<0$ (solid).
Dotted lines indicate the cylinder sets of the first generations. Broken lines
show the first two piecewise linear approximations $f_{1}$ and $f_{2}$. \newline
Right: (b) The largest eigenvalues (in modulus) of the Perron-Frobenius
operator for the map $F_c$ in dependence on $c$, as obtained using the
Lagrange-Chebyshev method with truncation order $n=25$. Positive/negative
eigenvalues are indicated by filled/open symbols. \label{fig:Counterexample}}
\end{figure}

In order to understand why the reasoning at the beginning of this section fails, that is,
why an approximation of the map $F_c$ by a piecewise linear Markov map fails to produce the correct mixing rate, let us consider increasingly finer
piecewise linear approximations $f_n$ of the map \eref{oscarmap},
see \Fref{fig:Counterexample}(a).
For each map $f_n$ the Perron-Frobenius operator restricted
to piecewise polynomial functions has a finite matrix representation (see \eref{blockmatrix}).
For the remainder of the section we shall only deal with the case
$\beta = 1$ and refer to $T(\beta)$ (and $\nu_m(\beta)$) as $T$ (and $\nu_m$).
\Fref{fig:EVs_Efunctions}(a) shows the numerical results
for the leading eigenvalues of the diagonal block $T^{(11)}$
and  $T^{(22)}$, that is, $\nu_{1}$ and
$\nu_{2}$, respectively, for increasing level of approximation $n=1,\ldots,6$. For comparison we display the subleading eigenvalue
$\nu^{(subl)}_0$ of $T^{(00)}$ as well.
For every level $n$ of approximation, $\nu_{2}$ gives the subleading
eigenvalue $\lambda_1$ of the Perron-Frobenius operator of the piecewise
linear approximation, and these values seem to converge as $n$ tends to infinity.
The values and the limit are  larger than $\exp(-2\Lambda)$, meaning that the
inequality \eref{ineq} in the Proposition \ref{thm:main} is satisfied, as expected.
The eigenfunction $u_n:[-1,1]\to \mathbb{R}$ of the corresponding eigenvalue
$\nu_{2}$ is a quadratic polynomial on each element of the partition, but
it develops an increasing number of discontinuities between different intervals of
the increasingly finer partition.
These eigenfunctions do not seem to converge to a smooth limit
(see \Fref{fig:EVs_Efunctions}(b)). In fact, unlike the invariant density there is no
reason why the limit should be smooth. The numerical experiment suggests that we end
up with a function of bounded variation.

\begin{figure}[h!]

\hspace{-0.6\textwidth}
\begin{center}
\includegraphics[width=1\textwidth]{oscarTent_c001_UlamApprox_EV_edit.eps}
\end{center}
 \caption{Left: (a) Leading eigenvalues
$\nu_1$ and $\nu_2$ of the  corresponding matrix blocks
$T^{(11)}$ and $T^{(22)}$
for a piecewise linear approximation of the map \eref{oscarmap} with $c=-0.11$, as a function
of the level of approximation. The subleading eigenvalue of
$T^{(00)}$, $\nu_0^{(subl)}$, is displayed as well.
For comparison, $\exp(-\Lambda)$ and $\exp(-2\Lambda)$ are depicted as well.
The broken lines are a guide for the eye. \newline
Right:  (b) Eigenfunction $u_n$ corresponding to $\nu_2$ for
$n=1,\ldots,6$ with normalisation $\int_I|u_n(x)|\, dx=1$. For clarity, successive
approximations are shifted by $0.5$. The open symbols indicate the discontinuity set of the
eigenfunction, i.e., the increasingly finer Markov partition of the piecewise linear approximaton.
\label{fig:EVs_Efunctions}}
\end{figure}

It is indeed possible to show and perhaps well known that an
estimate like \Eref{sharp} holds for discontinuous observables.
For that purpose let us consider the Perron-Frobenius operator
$\mathcal{L}_1$ on the space  of functions of bounded
variation.
Recall that a function $f:[-1,1] \to \mathbb{R}$ is of bounded variation if it has finite total variation
$\text{var}(f) = \sup \big\{ \sum_{i=1}^{p} |f(x_i)-f(x_{i-1})|:-1\leq x_0\leq \ldots \leq x_p \leq 1 \big \} < \infty$.
In this setup, the spectrum of the Perron-Frobenius operator associated with expanding maps has been studied in detail (see \cite{Kell_CMP84}).
In particular, there is an explicit formula for the essential spectral radius given by
\begin{equation*}
\sigma_{ess} = \lim_{k\to \infty} (\inf \{|(F_c^k)'(x)| : x \in[-1,1] \})^{-1/k} \, .
\end{equation*}
Thus, we have an upper bound for the mixing rate
\begin{equation*}
\alpha_{BV}\leq -\ln \sigma_{ess}=  \lim_{k\to \infty}
\frac{1}{k} \ln \inf \{|(F_c^k)'(x)| : x \in[-1,1] \} \,,
\end{equation*}
which yields the following estimate for the Lyapunov exponent:
\begin{eqnarray*}
\Lambda &=& \frac{1}{k} \int_I \ln |(F_c^k)'(x)| \,d\mu\nonumber \\
& \geq& \frac{1}{k}
\inf \{\ln|(F_c^k)'(x)| : x \in[-1,1] \} \int_I d \mu\nonumber \\
&=& \frac{1}{k} \ln \inf \{|(F_c^k)'(x)| : x \in[-1,1] \} \,.
\end{eqnarray*}
Thus, for observables of bounded variation we have the following result.

\begin{prop}\label{bvcor}
Let $f: I \rightarrow I$ be a piecewise monotonic smooth expanding interval map which is mixing with respect to its unique absolutely continuous
invariant measure.  Then the rate of decay of correlations for
functions of bounded variation is bounded by the Lyapunov exponent
\begin{equation}\label{bvest}
\alpha_{BV}\leq \Lambda \, .
\end{equation}
\end{prop}

In fact, almost identical statements can be found in
\cite{CoEc_JSP04}, for example, Corollary~9.2.

\section{Conclusion}\label{sec:4}

There is no simple, straightforward answer to the question about the relation between
Lyapunov exponents and mixing rates. On formal grounds one may argue
that both quantities probe entirely different and independent aspects of
a dynamical system, and that no particular relation should be expected.
Lyapunov exponents are determined by properties related to the largest eigenvalue
of the Perron-Frobenius operator. By contrast, correlation decay
depends crucially on properties of the observables, with mixing
rates being related to the subleading part of the spectrum.
Thus, abstract operator theory on its own does not seem to provide further insight into the
relation between both quantities. Witness, for example, the doubling map viewed as an analytic map
on the unit circle.
While its Lyapunov exponent is finite, the Perron-Frobenius operator
has no nontrivial eigenvalue when considered on the space of analytic
functions, that is, correlations of analytic observables decay faster than any exponential
(see, e.g., \cite{Bala:00}).

The argument outlined above, however, is a bit too simplistic. In fact,
our results on
piecewise linear expanding Markov maps observed via piecewise analytic functions or  general piecewise smooth expanding maps observed via functions of bounded
variation suggest that bounds on the mixing rate in terms of
Lyapunov exponents can be derived provided that specific properties of
the underlying dynamical system are taken into account.
Estimates of this type rely on nontrivial lower bounds for spectra. As such, they are complementary to
estimates which are available for proving the existence of spectral gaps,
and will thus require completely different approaches.

It turns out that the spirit of the result contained in Proposition \ref{bvcor}
can be understood by considering
an observable with a single discontinuity. In order to substantiate this claim
we have performed numerical simulations on the map \eref{oscarmap} for $c=-0.11$.
We have computed the autocorrelation function \eref{CorrelFunc} for the
observable $\varphi=\psi=\phi_h$ with $\phi_h(x)=x$ if $|x|<1/2$ and
$\phi_h(x)=x-\mbox{sign}(x) h$ if $|x|>1/2$, having a discontinuity of stepsize $h$
at $|x|=1/2$ (see \Fref{fig:4}). Choosing $h=0$, the corresponding observable $\phi_h$ is analytic and the
correlation decay is seen to follow the subleading eigenvalue $\lambda_1$ of the
Perron-Frobenius operator defined on the space of analytic functions
(see \Fref{fig:Counterexample}(b)). In the discontinuous case corresponding to $h\neq 0$
the short time initial decay of the correlations still follows the pattern of the
analytic observable, but
the correlation function now develops an exponential tail which obeys \Eref{bvest}.
The tail becomes more pronounced if the stepsize increases.
In fact, the mixing rate seems to be very close to the Lyapunov exponent.
Revisiting the considerations leading to Proposition \ref{bvcor},
it is tempting to surmise that this coincidence is a consequence of large deviation properties of finite time Lyapunov exponents,
since the expression for the essential spectral radius involves
an extreme value of a finite time Lyapunov exponent.
Thus, the relation between Lyapunov exponents and mixing rates for
observables of bounded variation could be viewed to arise from the same
mechanism already exploited in \cite{BaHeMePo_PRA88} for analysing the simple case mentioned
in the introduction.

\begin{figure}[h!]
\begin{center}
\includegraphics[width=.75\textwidth]{fig4_new.eps}
\end{center}
\caption{Normalised autocorrelation function of the discontinuous variable
$\phi_h$ with different
stepsizes $h$ for the map $F_c$ \eref{oscarmap} with $c=-0.11$. The straight lines represent an
exponential decay with rate $\ln|\lambda_1|$ and $-\Lambda$. Ergodic averages have been
computed as time averages of a series of length $2\times 10^4$ for $5\times 10^8$ uniformly
distributed initial conditions, skipping a transient of length $100$. The horizontal dotted
line indicates the order of magnitude of statistical errors induced by the finite
ensemble size.\label{fig:4}}
\end{figure}

This simple demonstration gives support to the folklore that correlation
decay is linked with Lyapunov exponents.
Even if a real world phenomenon is sufficiently well-modelled by a smooth dynamical
system for which to date no link between correlation decay and
Lyapunov exponents can be established, one should keep in mind that modern data processing
inevitably involves digital devices, which correspond to discontinuous observations.
Therefore, in formal terms observables of bounded variation could be the
relevant class for applications and in these cases Proposition \ref{bvcor} applies.

At an intuitive level it is easy to understand
why discontinuous observations
result in correlation decay related to Lyapunov exponents. A discontinuous
observable is able to distinguish between different ``microstates'' at a ``macroscopic'' level, that is, a discontinuous observation is able to
distinguish two states at infinitesimal distance. As the
distance between two nearby phase space points separated by a discontinuity
grows according to the Lyapunov exponent,
the sensitivity at the
microscale may be transported to the macroscale, that is, it may filter through to the
correlation function by a discontinuous observation.

In our context the mathematical challenge is
to establish a relation between mixing rates and Lyapunov exponents for
natural classes of observables, for example, full branch analytic interval maps
observed via analytic functions. Besides the need for developing tools to
obtain lower bounds for spectra, establishing the relation alluded to above also requires
a deeper understanding of which dynamical feature causes the point spectrum of
the Perron-Frobenius operator, that is, which signature of the microscopic dynamics
survives if viewed via analytic observables. This is reminiscent of coarse-graining approaches in statistical mechanics,
for example, the introduction of collective coordinates and quasi-particles. Thus,
tackling the mathematical problem above may well shed some light
on some of the most fundamental problems in contemporary nonequilibrium statistical physics
of complex systems.

\section*{Acknowledgements}
W.J. gratefully acknowledges support by EPSRC (grant no.\ EP/H04812X/1) and DFG (through SFB910),
as well as the kind hospitality by Eckehard Sch{\"o}ll and his group during the stay
at TU-Berlin.