%\pdfoutput=1\relax\pdfpagewidth=8.26in\pdfpageheight=11.69in\pdfcompresslevel=9
\documentclass[twoside]{article}
\usepackage{amsfonts} % used for R in Real numbers
\pagestyle{myheadings}
\markboth{ Center manifold and exponentially bounded solutions }
{ Luis Garcia \&  Hugo Leiva }
\begin{document}
\setcounter{page}{69}
\title{\vspace{-1in}\parbox{\linewidth}{\footnotesize\noindent
Nonlinear Differential Equations, \newline
Electron. J. Diff. Eqns., Conf. 05, 2000, pp. 69--80\newline
http://ejde.math.swt.edu or http://ejde.math.unt.edu
\newline ftp  ejde.math.swt.edu or ejde.math.unt.edu (login: ftp)}
 \vspace{\bigskipamount} \\
%
 Center manifold and exponentially-bounded solutions 
 of a forced Newtonian system with dissipation
% 
\thanks{ {\em Mathematics Subject Classifications:} 34A34, 34C27, 34C30.
 \hfil\break\indent 
{\em Key words:} center manifold, exponentially-bounded solutions.
 \hfil\break\indent
\copyright 2000 Southwest Texas State University. 
\hfil\break\indent Published October 24, 2000. \hfil\break\indent
Partially supported by CDCHT-ULA } } 

\date{}
\author{ Luis Garcia \&  Hugo Leiva }

\maketitle


\begin{abstract}
We study the existence of exponentially-bounded solutions to
the following system of second-order ordinary differential equations with 
dissipation:  
$$
u''+cu'+Au+kH(u) = P(t), \quad u \in {\mathbb R}^n, \quad t \in {\mathbb R},
$$
where $c$ and $k$ are positive constants, $H$ is a globally Lipschitz 
function, and $P$ is a bounded and continuous function. 
$A$ is a $n \times n$ symmetric matrix whose first eigenvalue 
is equal to zero and the others are positive.
Under these conditions, we prove that for some values of $c$, 
and $k$ there exist a continuous manifold such that solutions starting 
in this manifold are exponentially bounded.  
Our results are applied to the spatial discretization of well-known 
second-order partial differential equations with Neumann boundary 
conditions. 
\end{abstract}


\newtheorem{theorem}{Theorem} [section]
\newtheorem{lemma}[theorem]{Lemma} 
\newtheorem{corollary}[theorem]{Corollary}
\newtheorem{example}[theorem]{Example}
\renewcommand{\theequation}{\thesection.\arabic{equation}}

\section{Introduction}

In this note, we study the existence of exponentially bounded 
solutions of the following system of second-order ordinary 
differential equations with a damping force and dissipation
in ${\mathbb R}^n$:
\begin{equation} \label{S1}
u'' + cu' + Au + kH(u) = P(t), \quad u \in {\mathbb R}^n, \quad t \in {\mathbb R},
\end{equation}
where $c$ and $k$ are positive constants, $P \in C_{b}({\mathbb R};{\mathbb R}^n)$
(the space of continuous and bounded functions) and $H: {\mathbb R}^n \to 
{\mathbb R}^n$ is a globally Lipschitz function. i.e., there exists a 
constant $L>0$ such that
\begin{equation}\label{Lip}
\|H(U_1)-H(U_2) \| \leq L \|U_1 - U_2 \|, \quad U_1, U_2 \in {\mathbb R}^n.
\end{equation}
$A$ is a $n \times n$ symmetric matrix whose first eigenvalue 
$\lambda_1$ is equal to zero and the other eigenvalues are positive
$$
0= \lambda_1 < \lambda_{2} < \cdots \lambda_{l}
$$
each one with multiplicity $\gamma_j$ equal to the dimension 
of the corresponding eigenspace.

The equation (\ref{S1}) has been studied in \cite{HL} for the 
case that the first eigenvalue $\lambda_1$ of the matrix $A$ is 
positive ($\lambda_1 >0$); under these conditions they prove 
that for some values of $c$ the equation (\ref{S1}) has 
a bounded solution which is exponentially  stable and, if 
$P(t)$ is almost periodic, this bounded solution is also almost 
periodic.

The fact that the first eigenvalue $\lambda_1$ of the matrix 
$A$ is equal to zero, in general, does not allow
us to prove the existence of bounded solutions of (\ref{S1}).  
However, we prove that for some values of $c$ and $k$ 
there exist a positive number $\eta$ depending on $c$ and a 
continuous manifold ${\cal M}= {\cal M}(c,k,P(\cdot))$ such that any 
solution of the system (\ref{S1}) starting in ${\cal M}$ is exponentially 
bounded. i.e., 
$$
\sup_{t \in {\mathbb R}} e^{-\eta|t|}\left\{ \|u'(t) 
\|^{2} + \|u(t) \|^{2} \right\}^{1/2} < \infty.
$$
Our method is similar to the one used in \cite{HL}, 
we just rewrite the equation (\ref{S1}) as a first order system 
of ordinary differential equations and prove that the linear 
part of this system has an exponential trichotomy with trivial 
unstable space. Next, we use the variation constants formula and some ideas 
from \cite{V} \cite{V1} to find a formula for the exponentially bounded 
solutions of (\ref{S1}).  From this formula we can prove the existence 
of such manifold ${\cal M}= {\cal M}(c,k,P(\cdot))$.  These results are applied to 
the spatial discretization of very well known second order partial 
differential equations with Neumann boundary conditions:  

\begin{example} \rm 
The Sine-Gordon equation with Neumann boundary conditions is 
given by
\begin{eqnarray} \label{S2}
&U_{tt} + cU_{t} -dU_{xx} + k\sin{U} =  p(t,x), \quad 0 < x < L, \quad t 
\in {\mathbb R},& \\
&U_{x}(t,0) = U_{x}(t,L)  = 0,  \quad t \in {\mathbb R}, & \nonumber
\end{eqnarray}
where $c$, $d$ and $k$ are positive constants, $p: {\mathbb R} \times [0,L] 
\to {\mathbb R}$ is continuous and bounded. 


The following paragraph was taken from Temam's book(\cite{T}, 
pg. 184). With Neumann boundary the Sine-Gordon Equation
is physically interesting, the average value of the function $U$ 
is not expected to remain bounded and actually leads to 
nontrivial dynamics. From mathematical point of view this case 
is also interesting.

For each $N \in {\mathbb N}$ the spatial discretization of this equation 
is given by
\begin{eqnarray}
&u_{i}'' + cu_{i}' + d \delta^{-2}(2u_{i} - u_{i+1} - u_{i-1}) + 
k\sin{u_{i}} =  p_{i}(t), \quad 1 \leq i \leq N, \quad t \in {\mathbb R}, &
\nonumber\\
&u_{0} = u_1, \quad u_{N} = u_{N+1}  =  0\,. & \label{s3}
\end{eqnarray}
This equation can be written in the form of 
(\ref{S1}) with 
\begin{equation}\label{A1}
A=d \delta^{-2} \left(
\begin{array}{rrrrr}
1 & -1 & 0 & \cdots & 0 \\
-1 & 2 & -1 & \cdots & 0  \\
\vdots & \ddots & \ddots & \ddots & \vdots \\
0 & \cdots & -1 & 2 & -1 \\
0 & \cdots & 0 & -1 & 1
\end{array}
\right)
\end{equation}
where $\delta = L/(N+1)$.  The eigenvalues of this matrix are 
simple, the first one being zero and the others positive.
\end{example}

\begin{example} \rm
A telegraph equation with Neumman boundary conditions
\begin{eqnarray}\label{s2}
&U_{tt} + cU_{t} -dU_{xx} + \arctan {U} =  p(t,x), \quad 0 < x < L,
 \quad t \in {\mathbb R},& \\
&U_{x}(t,0) = U_{x}(t,L)  = 0,  \quad t \in {\mathbb R}, & \nonumber
\end{eqnarray}
where $c$, $d$ are positive constants, $p: {\mathbb R} \times [0,L] 
\to {\mathbb R}$ is continuous and bounded.

For each $N \in {\mathbb N}$ the spatial discretization of this equation 
is given by
\begin{eqnarray} \label{s4}
&u_{i}'' + cu_{i}' + d \delta^{-2}(2u_{i} - u_{i+1} - u_{i-1}) + 
\arctan {u_{i}} =  p_{i}(t) &\\
% \quad 1 \leq i \leq N, \quad t \in {\mathbb R}, 
&u_{0} = u_1, \quad u_{N} = u_{N+1}  =  0\,.  &\nonumber
\end{eqnarray}
This equation can be written in the form of 
(\ref{S1}) with the the same matrix $A$ as in (\ref{A1}).
\end{example}

\section{Preliminaries} \setcounter{equation}{0}

Most of the ideas presented in this section can be found in 
\cite{HL}. So, we shall prove only the new results.  
Equation (\ref{S1}) can be written as a first order system of 
ordinary differential equations in the space $W = {\mathbb R}^n \times 
{\mathbb R}^n$ as follows:
\begin{equation}\label{S21}
w' + {\cal A}w + k{\cal H}(w) = {\cal P}(t), \quad w \in W, 
\quad t \in {\mathbb R},
\end{equation}
where $v = u'$ and
\begin{equation}\label{S22}
w = \left( \begin{array}{r} u \\  v \end{array} \right), \quad
{\cal H}(w) = \left( \begin{array}{c}0 \\ H(u) \end{array}\right), \quad
{\cal P}(t) = \left( \begin{array}{c}0 \\ P(t)\end{array}\right), \quad
{\cal A} = \left(\begin{array}{rr}0 & -I \\ A & cI\end{array}\right).
\end{equation}
Now, we are ready to study the linear part of (\ref{S21}):
\begin{equation} \label{S23}
w' + {\cal A}w  = 0, \quad w \in W, \quad t \in {\mathbb R}.
\end{equation}
 For the rest of this article, we shall assume that each eigenvalue 
of the  matrix $A$ has multiplicity $\gamma_j$ equal to the 
dimension of the corresponding eigenspace and the first one is 
equal to zero and the others are positive.  Therefore, if  
$0 = \lambda_1 < \lambda_2 < \cdots \lambda_l$ are the 
eigenvalues of $A$, we have the following: 
\begin{description}

\item{\bf a)} There exists a complete orthonormal set 
$\{ \phi_{j,k} \}$ of eigenvector of $A$ in ${\mathbb R}^n$

\item{\bf b)} For all $x \in {\mathbb R}^n$ we have
\begin{equation} \label{S24}
Ax = \sum_{j = 1}^{l} \lambda_{j} \sum_{k = 1}^{\gamma_j}
\langle x, \phi_{j,k}\rangle  \phi_{j,k}
=\sum_{j = 1}^{l} \lambda_{j}E_{j}x,
\end{equation}
where $\langle \cdot, \cdot\rangle$ is the inner product in ${\mathbb R}^n$ and 
\begin{equation} \label{S25}
E_{j}x = \sum_{k = 1}^{\gamma_j} <x, \phi_{j,k}> \phi_{j,k}.
\end{equation}
So, $\{ E_j \}$ is a family of complete orthogonal projections 
in $R^n$ and $x = \sum_{j = 1}^{l}E_{j}x$, $x \in {\mathbb R}^n$

\item{\bf c)} The exponential matrix $e^{-At}$ is given by
\begin{equation} \label{S26}
 e^{-At} = \sum_{j = 1}^{l} e^{-\lambda_{j}t}E_{j}
\end{equation}
\end{description}

\begin{theorem}\label{T1}
Suppose that $c \not= 2\sqrt{\lambda_{j}}$, $j=1,2, \dots,l$. 
Then the exponential matrix $e^{-{\cal A}t}$ of the matrix 
$-{\cal A}$ given by (\ref{S22}) can be written as follow
\begin{equation}\label{S27}
e^{-{\cal A}t}w = \sum_{j = 1}^{l} \left\{ 
e^{\rho_1(j)t}Q_1(j)w + e^{\rho_{2}(j)t}Q_{2}(j)w \right\}, 
\quad w \in W, \quad t \in {\mathbb R},
\end{equation}
where
\begin{equation}\label{S28}
\rho(j) = {-c \underline{+} \sqrt{c^{2}-4 \lambda_{j}} \over 2}, \quad j=1,2, 
\dots,l
\end{equation}
and $\{ Q_{i}(j) : i=1,2 \}_{j=1}^{l}$ is a complete orthogonal 
system of projections in $W$.
\end{theorem}

\begin{corollary}\label{C1}
The spectrum $\sigma(-{\cal A})$ of the matrix $-{\cal A}$ is 
given by
$$
\sigma({-\cal A})= \{{-c \underline{+} \sqrt{c^{2}-4 \lambda_{j}} \over 2},
 \quad j=1,2, \dots,l \}.
$$
\end{corollary}

\begin{corollary}\label{C2}
Under the hypothesis of Theorem \ref{T1}, there exist two 
orthogonal projectors $\pi_{0}, \pi_{s}: W \to W$ and a 
constant $M > 0$ such that
$$\displaylines{
\| e^{-{\cal A}t} \pi_{0} \|  \leq  M, \quad t \in {\mathbb R},\cr
\| e^{-{\cal A}t} \pi_{s} \|  \leq   e^{-\beta t}, \quad t \geq 0,\cr
I_{W} = \pi_{0} + \pi_{s}, \quad  W = W_{0} \oplus W_{s}\,,
}$$
where $W_{0} = \mathop{\rm Ran}(\pi_{0}), \quad W_{s} = 
\mathop{\rm Ran}(\pi_{s})$, and $\beta =\beta(c)>0$ is given by
$$
-\beta = \max \big\{ -c, \quad 
\mathop{\rm Re}(\rho_{j})= \mathop{\rm Re} 
\big( {-c \pm \sqrt{c^{2}-4  \lambda_{j}} \over 2}  
\big): j=2, \dots , l. \; i=1,2. \big\}.
$$
\end{corollary}

\paragraph{Proof}
Since Theorem \ref{T1} and Corollary \ref{C1} are proved in 
\cite{HL}, we shall only prove Corollary \ref{C2}.  From 
Corollary \ref{C1}, we have that $\mathop{\rm Re}(\rho_{i}(j) < 0, \quad j=2, 
\dots,l$.  Therefore $\beta >0$. Since $\lambda_1 = 0$, the formula
 (\ref{S27}) can be written as
$$
e^{-{\cal A}t}w = Q_1(1)w + e^{-ct}Q_{2}(1)w + \sum_{j = 2}^{l} \left\{ 
e^{\rho_1(j)t}Q_1(j)w + e^{\rho_{2}(j)t}Q_{2}(j)w \right\}, 
$$
with $w \in W$ and $t \in {\mathbb R}$. If we put $\pi_{0} = Q_1(1)$ and 
$$
\pi_{s} = I-\pi_{0} = Q_{2}(1) + \sum_{j = 2}^{l} \left\{ Q_1(j) +
Q_{2}(j) \right\},
$$
we obtain the required projections, and
$$
\| e^{-{\cal A}t} \pi_{0}w \| = \| \pi_{0} e^{-{\cal A}t} w \| 
= \| Q_1(1) w \| \leq   \|Q_1(1) \| \|w \|, \quad t \in {\mathbb R}.
$$
Therefore,  $\| e^{-{\cal A}t} \pi_{0} \| \leq \|Q_1(1) \| = 
M$.
In the same way we get
\begin{eqnarray*}
\lefteqn{ \| e^{-{\cal A}t} \pi_{s}w \|^2 }\\
& = & \| \pi_{s} e^{-{\cal A}t} w 
\|^2 = \| e^{-ct}Q_{2}(1)w + \sum_{j = 2}^{l} \left\{ 
e^{\rho_1(j)t}Q_1(j)w + e^{\rho_{2}(j)t}Q_{2}(j)w \right\} 
\|^2\\
& = & e^{-2ct} \|Q_{2}(1)w \|^2 + \sum_{j = 2}^{l} \left\{ 
e^{2\mathop{\rm Re}{\rho_1(j)t}} \|Q_1(j)w \|^2 + 
e^{2\mathop{\rm Re}{\rho_{2}(j)t}} \|Q_{2}(j)w \|^2 \right\} \\
& \leq & e^{-2\beta t} \Big\{ \|Q_{2}(1)w \|^2 + \sum_{j = 2}^{l}  
 \|Q_1(j)w \|^2 + \|Q_{2}(j)w \|^2 \Big\}\\
& \leq & e^{-2\beta t} \| w \|^2, \quad t \geq 0.
\end{eqnarray*}
Therefore,  $\| e^{-{\cal A}t} \pi_{s} \|  \leq   e^{-\beta t}$ for 
$t \geq 0$.\hfill $\diamondsuit$

\begin{corollary}\label{C3}
For each $\epsilon \in [0, \beta)$ there exists some 
$M(\epsilon) >0$ such that
$$\displaylines{
\| e^{-{\cal A}t} \pi_{0} \|  \leq  M(\epsilon)e^{\epsilon 
|t|}, \quad t \in {\mathbb R}, \cr
\| e^{-{\cal A}t} \pi_{s} \|  \leq   M(\epsilon)
e^{-(\beta -\epsilon) t}, \quad t \geq  0\,.
}$$
\end{corollary}

\section{ Main Result} \setcounter{equation}{0}

In this section we shall prove the main Theorem of this paper 
under the hypothesis of Theorem \ref{T1} 
($c \not= 2\sqrt{ \lambda_{j}}, \quad j=1,2, \dots,l$).

The solution of (\ref{S21})  passing through the point $w_0$ at 
time $t=t_0$ is given by the variation constants formula
\begin{equation}\label{S31}
w(t) = e^{-{\cal A}(t-t_0)}w_{0} + \int_{t_0}^{t} e^{-{\cal 
A}(t-s)}\left\{-k{\cal H}(w(s)) + {\cal P}(s) \right\}ds, \quad t 
\in {\mathbb R}.
\end{equation}
We shall use the following notation:  For each $\eta \geq 0$
we denote by $Z_{\eta}$ the Banach space
\begin{equation}
Z_{\eta} = \Big\{ z \in C({\mathbb R};W) :  \|z \|_{\eta} = \sup_{t 
\in {\mathbb R}} e^{-\eta |t|} \|z(t) \| < \infty \Big\}.
\end{equation}
In particular, $Z_{0} = C_{b}({\mathbb R}, W)$ the space of 
bounded and continuous functions defined in ${\mathbb R}$ taking values 
in $W={\mathbb R}^n \times {\mathbb R}^n$.

\begin{theorem}\label{T2}
Suppose that $H$ is a bounded function or 
$H(0) =0$.  Then for some $c$ and $k$ positive there exist 
$\eta = \eta(c) \in (0, \beta)$ and a continuous manifold ${\cal M}= {\cal 
M}(c,k,P)$ such that any solution $u(t)$ of (\ref{S1}) with
$(u(0), u'(0)) \in {\cal M}$ satisfies:
\begin{equation}\label{M2}
\sup_{t \in {\mathbb R}} e^{-\eta |t|} \left\{ \|u(t) \|^2 + \|u'(t) 
\|^2 \right\}^{1/2} < \infty,
\end{equation}
where $\beta$, $W_0$ and $W_s$ are as in corollary \ref{C2}.  Moreover,
\begin{description}
\item{\bf (a)} There exists a globally Lipschitz function $\psi 
:W_{0} \to W_{s}$ depending on $P$ such that
\begin{equation}\label{M3}
{\cal M} = \{ w_{0} + \psi(w_{0}) \quad : w_{0} \in W_{0} \},
\end{equation}
and if we put $\psi(w) = \psi(w,P)$, there exist $M \geq 1$ 
and $0 < \Gamma <1$ such that 
\begin{equation}
\|\psi(w_1, P_1)-\psi(w_2, P_2) \| \leq {kLM(1-\Gamma)^{-1} 
\over \beta - \eta} \|w_1 -w_2 \| +  {1 \over \beta} \|P_1 - 
P_2  \|, 
\end{equation}
for $w_1, w_2  \in W_{0}, \quad P_1, P_2 \in C_{b}({\mathbb R} ,{\mathbb R}^n)$.

\item{\bf (b)} If $H$ is bounded, then $\psi$ is also bounded.

\item{\bf (c)} If $P = 0$ and $H(0) = 0$, then ${\cal M}$ is unique 
and invariant under the equation $w' + {\cal A}w + k{\cal H}(w) 
= 0$. In this case ${\cal M}$ is called center manifold and it 
is tangent to the space $W_{0}$ at $w_{0} = 0$.
\end{description}
\end{theorem}
Before proving  the main theorem, we need some previous results. 

\begin{lemma}\label{L1}
Let $z \in Z_{0}=C_{b}({\mathbb R},W)$.  Then, $z$ is a solution of 
(\ref{S21}) if and only if there exists some $w_{0} \in W_{0}$ 
such that
\begin{eqnarray}\nonumber
z(t) & = & e^{-{\cal A}t}w_{0} +  \int_{0}^{t} e^{-{\cal 
A}(t-\tau)} \pi_{0} \left\{-k{\cal H}(z(\tau)) + {\cal P}(\tau) 
\right\}d\tau\\
&&+ \int_{-\infty}^{t} e^{-{\cal 
A}(t-\tau)} \pi_{s} \left\{-k{\cal H}(z(\tau)) + {\cal P}(\tau) \right\}d\tau, \quad t 
\in {\mathbb R}. \label{S32}
\end{eqnarray}
\end{lemma}

\paragraph{Proof}
Suppose that $z$ is a solution of (\ref{S21}).  Then, from 
Corollary \ref{C2} we get  $z(t) = \pi_{0}z(t) + 
\pi_{s}z(t)$ and from the variation of constants formula (\ref{S31})
we obtain
\begin{equation}\label{S33}
\pi_{0}z(t) = e^{-{\cal A}t}\pi_{0}z(0) + \int_{0}^{t} e^{-{\cal 
A}(t-\tau)} \pi_{0}\left\{-k{\cal H}(z(\tau)) + {\cal P}(\tau) 
\right\}d\tau, \quad t \in {\mathbb R},
\end{equation}
and
\begin{equation}\label{S34}
\pi_{s}z(t) = e^{-{\cal A}(t-t_0)}\pi_{s}z(t_0) + \int_{t_0}^{t} e^{-{\cal 
A}(t-\tau)} \pi_{s}\left\{-k{\cal H}(z(\tau)) + {\cal P}(\tau) 
\right\}d\tau, \quad t \in {\mathbb R}.
\end{equation}
Since $z(t)$ is bounded, there exists $R >0$ such that $\|z(t) 
\| \leq R$, for all $t \in {\mathbb R}$. Then , from corollary \ref{C2} 
we obtain that
$$
\| e^{-{\cal A}(t-t_0)}\pi_{s}z(t_0) \| \leq Re^{-\beta(t-t_0)} 
\to 0 \quad \mbox{as} \quad t_0 \to - \infty.
$$
Now, if we put
$$
l= \sup_{\tau \in {\mathbb R}} \|{\cal H}(z(\tau)) \| \quad \mbox{and} \quad 
L_{p} = \sup_{\tau \in {\mathbb R}} \|{\cal P}(\tau) \|,
$$
then 
\begin{eqnarray*}
\| \int_{-\infty}^{t} e^{-{\cal A}(t-\tau)} 
\pi_{s}\left\{-k{\cal H}(z(\tau)) + {\cal P}(\tau) 
\right\}d\tau \| & \leq & \| \int_{-\infty}^{t} e^{-\beta(t-\tau)} 
\left\{kl + L_{p} \right\}d\tau \| \\
 & = & {kl + L_{p} \over \beta }.
\end{eqnarray*}
Hence, passing to the limit in (\ref{S34}), as $t_0$ approaches  
$-\infty$,  we obtain
\begin{equation}\label{S35}
\pi_{s}z(t) = \int_{-\infty}^{t} e^{-{\cal 
A}(t-\tau)} \pi_{s}\left\{-k{\cal H}(z(\tau)) + {\cal P}(\tau) 
\right\}d\tau, \quad t \in {\mathbb R}.
\end{equation}
Therefore, putting $w_{0} = \pi_{0}z(0)$ we get (\ref{S32}).

Conversely, suppose that $z$ is a solution of (\ref{S32}).  
Then 
\begin{eqnarray*}
z(t) & = & e^{-{\cal A}t}w_{0} +  \int_{0}^{t} e^{-{\cal 
A}(t-\tau)} \pi_{0} \left\{-k{\cal H}(z(\tau)) + {\cal P}(\tau) 
\right\}d\tau\\
&&+ \int_{0}^{t} e^{-{\cal 
A}(t-\tau)} \pi_{s} \left\{-k{\cal H}(z(\tau)) + {\cal P}(\tau) 
\right\}d\tau\\
&&+ \int_{-\infty}^{0} e^{-{\cal 
A}(t-\tau)} \pi_{s} \left\{-k{\cal H}(z(\tau)) + {\cal P}(\tau) 
\right\}d\tau\\
& = & e^{-{\cal A}t} \left\{w_0 + \int_{-\infty}^{0} e^{{\cal 
A} \tau} \pi_{s} \left\{-k{\cal H}(z(\tau)) + {\cal P}(\tau) 
\right\}d\tau \right\}\\
&&+ \int_{0}^{t} e^{-{\cal 
A}(t-\tau)} \left\{-k{\cal H}(z(\tau)) + {\cal P}(\tau) 
\right\}d\tau\\
& = & e^{-{\cal A}t}z(0) + \int_{0}^{t} e^{-{\cal 
A}(t-\tau)} \pi_{0} \left\{-k{\cal H}(z(\tau)) + {\cal P}(\tau) 
\right\}d\tau,
\end{eqnarray*}
where
\begin{equation}\label{S36}
z(0) = w_0 + \int_{-\infty}^{0} e^{{\cal 
A} \tau} \pi_{s} \left\{-k{\cal H}(z(\tau)) + {\cal P}(\tau) 
\right\}d\tau.
\end{equation}
This concludes the proof of the lemma. \hfill $\diamondsuit$

\begin{lemma}\label{L2}
Suppose that $H(0)=0$ and $z \in Z_{\eta}$ for $\eta \in 
[0, \beta)$.  Then, $z$ is a solution of (\ref{S21}) if and only 
if there exists some $w_{0} \in W_{0}$ such that $z$ satisfies 
(\ref{S32}).
\end{lemma}

\paragraph{Proof}
Suppose that $z$ is a solution of (\ref{S21}). Then, in the 
same way as the proof of lemma \ref{L1}, we consider
$$
\pi_{0}z(t) = e^{-{\cal A}t} \pi_{0}z(0) + \int_{0}^{t} e^{-{\cal 
A}(t-\tau)} \pi_{0}\left\{-k{\cal H}(z(\tau)) + {\cal P}(\tau) 
\right\}d\tau, \quad t \in {\mathbb R}.
$$
and
$$
\pi_{s}z(t) = e^{-{\cal A}(t-t_0)} \pi_{s}z(t_0) + \int_{t_0}^{t} e^{-{\cal 
A}(t-\tau)} \pi_{s}\left\{-k{\cal H}(z(\tau)) + {\cal P}(\tau) 
\right\}d\tau, \quad t \in {\mathbb R}.
$$
Since $z$ belongs to $Z_{\eta}$, there exists $R >0$ such that $\|z(t) 
\| \leq Re^{\eta |t|}$, for all $t \in {\mathbb R}$.  Fix some $t \in 
{\mathbb R}$ and let $t_0 \leq \min \{t, 0 \}$; then we have
$$
\| e^{-{\cal A}(t-t_0)}\pi_{s}z(t_0) \| \leq Re^{-\beta(t-t_0)} 
e^{-\eta t_{0}} = Re^{-\beta t} e^{(\beta-\eta) t_{0}} 
\to 0 \quad \mbox{as} \quad t_0 \to - \infty.
$$
On the other hand, we obtain the following estimate
\begin{equation}\label{S37}
\| \int_{-\infty}^{t} e^{-{\cal A}(t-\tau)} 
\pi_{s}\left\{-k{\cal H}(z(\tau)) + {\cal P}(\tau) 
\right\}d\tau \|  \leq e^{\eta |t|} \left\{ {kLR \over 
\beta + \eta } + {kLR \over 
\beta - \eta } + {L_{p} \over 
\beta } \right\},
\end{equation}
where $L$ is the Lipschitz constant of $H$.  Hence, putting $w_{0} = 
\pi_{0}z(0)$  and passing to the limit when 
$t_0$ goes to $-\infty$ we get (\ref{S32}).

The converse follows in the same way as the foregoing lemma.
\hfill $\diamondsuit$


\begin{lemma}\label{L5}
Suppose that $H$ is bounded and $z \in Z_{\eta}$ for $\eta \in 
[0, \beta)$.  Then, $z$ is a solution of (\ref{S21}) if and only 
if there exists some $w_{0} \in W_{0}$ such that $z$ satisfies 
(\ref{S32}).
\end{lemma}

Now, from (\ref{S32}) we only have to prove that the following 
set 
\begin{equation}\label{M1}
{\cal M} = {\cal M}(c,k, P) = \{ z(0) : z \in Z_{\eta}, 
\mbox{ and satisfies  (\ref{S32}) }\}
\end{equation}
is a continuous manifold for some values of $c$ and $\eta 
\in (0, \beta(c))$.
 From (\ref{S36}) we get that
\begin{equation}\label{M4}
{\cal M}= \{w_0 +\pi_{s}z(0) \quad : (w_0, z) \in W_{0} \times Z_{\eta}, \quad 
(w_0, z) \quad \mbox{satisfying} \quad (\ref{S32}) \}
\end{equation}
We shall need the following definition and notations:

\paragraph{Definition} \begin{description}

\item{(a)} For each $w_0 \in W_0$ we define the function 
$Sw_0 : {\mathbb R} \to W$ by:
$$
(Sw_0)(t) = e^{-{\cal A}t}w_0, \quad t \in {\mathbb R};
$$

\item{(b)} for each function $z: {\mathbb R} \to W$ we define 
the non-autonomous Nemytski operator $G(z): {\mathbb R} \to W$ by
$$
G(z)(t) = -k{\cal H}(z(t)) + {\cal P}(t), \quad t \in {\mathbb R};
$$

\item{(c)}  for those functions $z: {\mathbb R} \to W$ for which the 
integrals make sense we define $Kz: {\mathbb R} \to W$ by
$$
Kz(t) = \int_{0}^{t} e^{-{\cal A}(t-\tau)} \pi_{0}z(\tau) d\tau\\
  +  \int_{-\infty}^{t} e^{-{\cal 
A}(t-\tau)} \pi_{s}z(\tau)d\tau, \quad t \in {\mathbb R}.
$$
\end{description}

\noindent
With this notation,  (\ref{S32}) can be written in 
the following equivalent form in $Z_{\eta}$
\begin{equation}\label{E1}
z = Sw_0 + K \circ G(z).
\end{equation}
The proof of the following lemma is not hard:

\begin{lemma}\label{L4} \begin{description}

\item{(a)} $S$ is a bounded operator from $W_0$ into $Z_{\eta}$ for 
each $\eta \geq 0$.

\item{(b)} if ${\cal H}(0) = 0$ or ${\cal H}$ is bounded, 
then $G$ maps $Z_{\eta}$ into itself for $\eta \geq 0$ and
$$
\|G(z_1) - G(z_2) \|_{\eta} \leq kL\|z_1 - z_2 \|_{\eta}, \quad 
z_1, z_2 \in Z_{\eta};
$$
\item{(c)} for $\eta \in (0,\beta)$ the linear operator $K$ is bounded 
from $Z_{\eta}$ into itself and
\begin{equation}\label{K1}
\|K \|_{\eta} \leq R(c) = M(\epsilon) \left\{ {1 \over \eta - 
\epsilon} + {1 \over \beta(c) -\eta  } \right\},
\end{equation}
where $L$ is the Lipschitz constant of $H$, $\beta$ is given by 
corollary \ref{C2}, $0 < \epsilon < \eta < \beta$ and $M(\epsilon)$ is given by 
corollary \ref{C3}.
\end{description}
\end{lemma}

\begin{lemma}\label{LL4}
Let $c >0$, $k >0$ and $\eta \in (0,\beta)$ such that
\begin{equation}\label{H1}
\Gamma = \|K \|_{\eta} kL < 1.
\end{equation}
Then $(I-K \circ G): Z_{\eta} \to Z_{\eta}$ is a homeomorphism 
with inverse $\Psi: Z_{\eta} \to Z_{\eta}$ and the 
manifold ${\cal M}= {\cal M}(c,k,P)$ is given by
\begin{equation}\label{M5}
{\cal M} = \{ w_{0} + \pi_{s}\Psi(Sw_{0})(0) \quad : w_{0} \in W_{0} \}.
\end{equation}
\end{lemma}

\paragraph{Proof} 
It follows from Lemma \ref{L4} that $K \circ G$ maps $Z_{\eta}$ 
into itself for $\eta \in (0,\beta)$ and is globally 
Lipschitzian with Lipschitz constant $\Gamma$.  Then, under the 
condition (\ref{H1}) the map $(I-K \circ G): Z_{\eta} 
\to Z_{\eta}$ is invertible, with inverse $\Psi : Z_{\eta} 
\to Z_{\eta}$ which is also globally Lipschitzian with 
Lipschitz constant $(1- \Gamma)^{-1}$.  In particular $\Psi$ is 
a continuous function.  Therefore, the equation (\ref{E1}) has 
a unique solution given by
\begin{equation}\label{EE1}
z(t) = (I - K \circ G)^{-1}(Sw_0)(t) = \Psi (Sw_0)(t), \quad t 
\in {\mathbb R}.
\end{equation}
Hence, from (\ref{M4}) we get (\ref{M5}).
\hfill $\diamondsuit$

\paragraph{Proof of Theorem \ref{T2}.}  If we take for example $\eta = 
\beta /2$, then $R(c)$ given by (\ref{K1}) can be written as 
follow
$$
R(c) =  M(\epsilon) \left\{ {2 \over \beta - 
2 \epsilon} + {2 \over \beta } \right\},
$$
with $0 < \epsilon < {\beta \over 2}$.  Hence, from $\lim_{c \to 
\infty} R(c) = 0$.  So, we can choose 
$c$ such that 
$$
\Gamma = \|K \|_{\eta}kL \leq R(c)kL <1.
$$
Then, using Lemma \ref{LL4} we obtain the first conclusion of 
the Theorem and defining  $\psi: W_0 \to W_{s}$ by
$$
\psi(w_0) = \pi_{s} \Psi(Sw_0)(0), \quad w_0 \in W_0,
$$
we get (\ref{M3}).  Clearly, the function $\psi$ is globally 
Lipschitzian.

On the other hand, from (\ref{S36}) we have
\begin{eqnarray*}
z(0) & = & w_0 + \pi_{s}z(0) = w_0 + \psi(w_0)\\
& = & w_0 + \int_{-\infty}^{0} e^{{\cal 
A} \tau} \pi_{s} \left\{-k{\cal H}(\Psi(Sw_0)(\tau)) + {\cal P}(\tau) 
\right\}d\tau.
\end{eqnarray*}
Therefore,
\begin{equation}\label{P1}
\psi(w_0,P) = \int_{-\infty}^{0} e^{{\cal 
A} \tau} \pi_{s} \left\{-k{\cal H}(\Psi(Sw_0)(\tau)) + {\cal P}(\tau) 
\right\}d\tau.
\end{equation}
To complete the proof of part (a), let us consider 
$w_1, w_2  \in W, \quad P_1, P_2 \in C_{b}({\mathbb R} ,{\mathbb R}^n)$ and
\begin{eqnarray*}
\psi(w_1 , P_1) - \psi(w_2,P_2) & = & \int_{-\infty}^{0} -ke^{{\cal 
A} \tau} \pi_{s}  \left\{{\cal H}(\Psi(Sw_1)(\tau)) -
{\cal H}(\Psi(Sw_2)(\tau)) \right\}d\tau\\
&&+ \int_{-\infty}^{0} e^{{\cal 
A} \tau} \pi_{s} \left\{{\cal P}_1(\tau) -{\cal P}_{2}(\tau) \right\}d\tau.
\end{eqnarray*}
 From Corollary \ref{C2} we get that
\begin{eqnarray*}
\|Sw_1 - Sw_2 \|_{\eta} & = & \sup_{\tau \in {\mathbb R}} e^{-\eta |\tau |} 
\|e^{{\cal A} \tau}(Sw_1 - Sw_2) \|\\
& \leq & M \|w_1 -w_2 \|,
\end{eqnarray*}
and from Lemma \ref{LL4} we get that
\begin{eqnarray*}
\|\Psi(Sw_1) - \Psi(Sw_2) \|_{\eta} & \leq & (1-\Gamma)^{-1} \|Sw_1 - Sw_2 
\|_{\eta}\\ & \leq & M (1-\Gamma)^{-1} \|w_1 - w_2 \|.
\end{eqnarray*}
Therefore,
\begin{eqnarray*}
\|\psi(w_1 , P_1) - \psi(w_2,P_2)\| & \leq & \int_{-\infty}^{0} 
kLM(1-\Gamma)^{-1} \|w_1 -w_2 \|e^{(\beta - \eta)\tau}d\tau\\
& &+ \int_{-\infty}^{0} \|P_1 -P_2 \|e^{\beta \tau}d\tau.
\end{eqnarray*}
Hence,
$$
\|\psi(w_1, P_1)-\psi(w_2, P_2) \| \leq {kLM(1-\Gamma)^{-1} 
\over \beta - \eta} \|w_1 -w_2 \| +  {1 \over \beta} \|P_1 - 
P_2  \|.
$$
To prove part (b), let us suppose that:
$ \| H(u) \| \leq l, \quad u \in {\mathbb R}^n \quad \mbox{and} \quad L_{p}= 
\sup_{\tau \in {\mathbb R}} \|P(\tau) \|$.  Then, from (\ref{P1}) we get 
that
$$
\| \psi(w_0) \| \leq \int_{-\infty}^{0} e^{\beta \tau} \{ kl + 
L_{P} \}d\tau \leq {kl + L_P \over \beta}, \quad w_0 \in W_0.
$$
Part (c) follows from Theorem 2.1 of \cite{V}.
\hfill $\diamondsuit$

\paragraph{Remark}
The equation (\ref{S32}) may not have bounded solutions in 
${\mathbb R}$.  However, if $H=0$ and ${\cal P}$ satisfies the condition
$$
\sup \left\{ \mid \int_{0}^{t} \| \pi_{0} {\cal P}(\tau) 
\|d \tau \mid : \quad t \in {\mathbb R} \right\} < \infty,
$$
then for each $w_0 \in W_0$ the equation (\ref{S32}) has a 
bounded solution which is given by
$$
z(t)  =  e^{-{\cal A}t}w_{0} +  \int_{0}^{t} e^{-{\cal 
A}(t-\tau)} \pi_{0}{\cal P}(\tau)d\tau
 +  \int_{-\infty}^{t} e^{-{\cal 
A}(t-\tau)} \pi_{s}{\cal P}(\tau)d\tau, \quad t 
\in {\mathbb R}.
$$
An open question, is the following:
What conditions do we have to impose to the functions $H$ and 
$P$ to insure the existence of bounded solutions of the 
equation (\ref{S32})?

 

\begin{thebibliography}{99}

\bibitem {AO2} J.M. Alonso and R. Ortega, ``Global asymptotic stability 
of a forced newtonian system with dissipation'', J. Math. Anal. 
and Applications {\bf 196} (1995), 965-986.

\bibitem {AMO} J.M. Alonso, J. Mawhin and R. Ortega, ``Bounded 
solutions of second order semilinear evolution equations and 
applications to the telegraph equation'', J. Math. Pures Appl., 
{\bf 78}, (1999), p. 49-663.

\bibitem {CJM} Y.S. Choi, K.C. Jen, and P.J. McKenna, `` The 
structure of the solution set for periodic oscillations in a 
suspension bridge model'', IMA J. Appl. Math. {\bf 47} (1991), 
283-306. 

\bibitem {GLM} J. Glover, A.C. Lazer and P.J. McKenna, `` 
Existence and stability of large scale nonlinear oscillations 
in suspension bridges'' J. Appl. Math. Phys. {\bf 40} (1989), 
172-200.

\bibitem {HL} H. Leiva, ``Existence of Bounded Solutions 
of a Second Order System with Dissipation" J. Math. Anal. and 
Applications {\bf 237}, 288-302 (1999).

\bibitem {T} R. Temam, ``Infinite-dimensional dynamical systems 
in mechanics and physics'', Appl. Math. Scien. Vol., 68, 
Springer-Verlag(1988).

\bibitem {V} A. Vanderbauwhede, ``Center Manifolds, 
Normal Forms and Elementary Bifurcations'' Dyns.  Reported {\bf 
2}(1987), 89-170.  

\bibitem {V1} A. Vanderbauwhede, ``Center Manifolds, Normal and 
Contractions on a Scale of B-Spaces'', J. Funct. Analy. {\bf 72}(1987), 
209-224.


\end{thebibliography}\medskip

\noindent{\sc Luis Garcia \&  Hugo Leiva}\\
Universidad de los Andes, Facultad de Ciencias \\
Departamento de Matem\'atica \\
M\'erida 5101-Venezuela \\
e-mail: hleiva@ula.ve

                                                                          

\end{document}
