# Heteroclinic Solutions in Singularly Perturbed Discontinuous Differential
Equations

Flaviano Battelli
Department of Industrial Engineering and Mathematics
Marche Polytecnic University
Ancona - Italy  Michal Fečkan
Department of Mathematical Analysis and Numerical Analysis
Comenius University in Bratislava
Mlynská dolina, 842 48 Bratislava – Slovakia
and Institute of Mathematics, Slovak Academy of Sciences
Štefánikova 49, 814 73 Bratislava – Slovakia Partially supported by the
National Natural Science Foundation of China (12161015), the Slovak Research
and Development Agency under the contract No. APVV-18-0308 and by the Slovak
Grant Agency VEGA No. 1/0084/23 and No. 2/0127/20 JinRong Wang
Department of Mathematics
Guizhou University, Guiyang
Guizhou 550025 – China

###### Abstract

We derive Melnikov type conditions for the persistence of heteroclinic
solutions in perturbed slowly varying discontinuous differential equations
extending to these equations similar results for continuous differential
equations.

Key Words: discontinuous differential equations; heteroclinic solutions;
Melnikov conditions; persistence.

##  1 Introduction

In recent years there has been a great deal in the study of singularly
perturbed equations such as

$\begin{array}[]{l}\dot{x}=f(x,y)\\\ \dot{y}=\varepsilon g(x,y).\end{array}$
(1.1)

where $f(x,y)$ and $g(x,y)$ are sufficiently smooth functions. To the best of
our knowledge this study started with remarkable papers by Tichonov [24],
Vasil’eva et al. [25, 26]. All these results concern the behaviour of the
solutions of (1.1) on a finite time scale. Later, Hoppensteadt [15, 16]
extended the result on an infinite time scale. Then Fenichel [12] developed
his geometric singularly perturbed theory that has been widely used by other
authors as, for example, Szmolyan [23]. These last results essentially prove
that the $x$ component of the solution is asymptotic to invariant manifolds,
sometimes called centre manifolds, described by equations like $x=X_{\pm}(y)$
(the case $X_{+}(y)=X_{-}(y)$ is allowed). In [18] a combination of the
Melnikov method and geometric singular perturbation theory is presented. This
approach has been then extended in [14] to investigate a mechanism of chaos
near resonances both in the dissipative and the Hamiltonian context of
differential equations.

The basic starting point for this kind of results is that, for all $y$ the
frozen equation

$\dot{x}=f(x,y)$ (1.2)

has hyperbolic fixed points, say $u_{\pm}(y)$, that are bounded functions of
$y$ together with their derivatives, and solutions $u_{\pm}(t,y)$, defined for
$t\geq 0$ and $t\leq 0$, resp., such that

$\lim_{t\to\infty}u_{+}(t,y)-u_{+}(y)=0,\;\lim_{t\to-\infty}u_{-}(t,y)-u_{-}(y)=0$
(1.3)

uniformly with respect to $y$ and there exists $y_{0}$ such that
$u_{-}(0,y_{0})=u_{+}(0,y_{0})$. Hence, for $y=y_{0}$ the frozen equation
(1.2) has a heteroclinic orbit connecting $u_{-}(y_{0})$ to $u_{+}(y_{0})$.
Note that in [6, 7, 20] it is also assumed that $u_{\pm}(t,y)$ are bounded
together with their derivatives. However, in Lemma 4.2, we prove that this
condition is automatically satisfied when $u_{\pm}(y)$ are bounded together
with their derivatives and (1.3) holds. The proof of this result is rather
technical and may be skipped at a first reading of this paper.

Other relevant results concerning (1.1) appeared in [6, 7, 8, 17, 20]. In
[17], the second order equation with slowly varying coefficients

$\ddot{x}+x(x-a(\varepsilon t+\theta))(1-x)=0,\quad\varepsilon>0$ (1.4)

is studied, where $a(y)$ is a $C^{1}$, $1$-periodic function, such that
$0<a(y)<1$. When $\varepsilon=0$ and $a(y)=\frac{1}{2}$, (1.4) has the
heteroclinic solution

$x_{het}(t)=\frac{1}{1+e^{-\frac{t}{\sqrt{2}}}}.$

However this heteroclinic solution is broken when $a(y)\neq\frac{1}{2}$. In
[17] it has been proved that if $a(y)$ has a transversal intersection with
$a=\frac{1}{2}$ at $y=y_{0}$, i.e. $2a(y_{0})=1$ and $a^{\prime}(y_{0})\neq
0$, then equation (1.4) has a heteroclinic solution. Writing

$x_{1}=x,\quad x_{2}=\dot{x},\quad y=\varepsilon t+\alpha$

equation (1.4) reads:

$\begin{array}[]{l}\dot{x}=f(x,y)\\\ \dot{y}=\varepsilon\end{array}$ (1.5)

where

$f(x,y)=\begin{pmatrix}x_{2}\\\ x_{1}(x_{1}-1)(x_{1}-a(y))\end{pmatrix}$

and $u(t,y_{0})=(x_{het}(t),\dot{x}_{het}(t)$. This solution is a part of a
family of solutions $u_{\pm}(t,y)$ such that

$\lim_{t\to-\infty}u_{-}(t,y)=(0,0)\quad\lim_{t\to\infty}u_{+}(t,y)=(1,0)$

and it can be proved that the transversality condition on $a(y)$ corresponds
to the fact that the real-valued function

$\psi^{*}[u_{-}(0,y)-u_{+}(0,y)]$

where $\psi\in\mathbb{R}^{n}$ is a suitable vector, has a simple zero at
$y=y_{0}$.

In Theorem 6.2, we extend these results, as far as the asymptotic behaviour of
the $x$-component of the solution is concerned, to discontinuous differential
equations that is a differential equation which is described in different ways
according to the region the solutions belong to. The first difference with the
continuous case is that the way how $u_{\pm}(t,y)$ passes from one domain into
another is important. In this paper we assume that $u_{\pm}(t,y)$ intersects
transversally the boundary of the domain when they pass from the domain into
another (see assumption $A_{1})$ in Section 2).

Problems like the above where $u(t,\alpha,y)$ is, instead, a family of
periodic solutions depending on some further parameters $\alpha\in\mathbb{R}$
and $y\in\mathbb{R}^{m}$, have been studied in [27, 4, 5]. Note that, in [27],
$f(x,y)$ is replaced by $f_{0}(x,y)+\varepsilon f_{1}(x,y,t,\varepsilon)$ and
it is assumed that system

$\dot{x}=f_{0}(x,y)$ (1.6)

has a one-parameter family of periodic solutions $q(t-\theta,y,\alpha)$ with
period $T(y,\alpha)$ being $C^{r}$ in $(y,\alpha)$ and it is because of the
$t$-dependence of the perturbed equation that the extra variable $\theta$ has
been introduced. Then a vector valued function $M^{m/n}(y,\alpha,\theta)$ is
constructed, that they called _subharmonic Melnikov function_ , which is a
measure of the difference between the starting value and the value of the
solution at the time $\frac{m}{n}T$ in a direction transverse to the
unperturbed vector field at the starting point and proved that periodic
solutions of the perturbed vector field arise near the simple zeros of
$M^{m/n}(y,\alpha,\theta)$. In [5] this result concerning the existence of
periodic solutions has been extended to discontinuous systems of differential
equations as the one we consider here.

Now, let us be more precise and give the definition of the discontinuous
differential equation we are studying in this paper.

Let $h(x,y):\mathbb{R}^{n}\times\mathbb{R}^{m}\to\mathbb{R}$ be a $C^{r}$
function, $r\geq 2$, with bounded derivatives,
$f_{\ell}(x,y):\mathbb{R}^{n}\times\mathbb{R}^{m}\to\mathbb{R}$,
$\ell=-M,\ldots,N$, be $C^{r}$-functions, bounded on
$\mathbb{R}^{n}\times\mathbb{R}^{m}$ together with their derivatives, and
$c_{-M}<\ldots<c_{0}\ldots<c_{N}$ be real numbers. By discontinuous
differential equation we mean an equation like

$\begin{array}[]{l}\dot{x}=f(x,y)\\\ \dot{y}=\varepsilon
g(x,y,\varepsilon)\end{array}$ (1.7)

where $x\in\mathbb{R}^{n}$, $y\in\mathbb{R}^{m}$, $\varepsilon\in\mathbb{R}$,
$\varepsilon>0$ and

$f(x,y):=\left\\{\begin{array}[]{ll}f_{-M}(x,y)&\hbox{if $h(x,y)<c_{-M}$}\\\
f_{\ell}(x,y)&\hbox{if $c_{\ell-1}<h(x,y)<c_{\ell}$}\\\
&\ell=-M+1,\ldots,N\end{array}\right.$ (1.8)

It is assumed that the equation

$\dot{x}=f_{N}(x,y)\quad\hbox{[resp. $\dot{x}=f_{-M}(x,y)$]}$ (1.9)

has a hyperbolic fixed point $x=w_{+}(y)$ [resp. $x=w_{-}(y)$] together with
continuous, piecewise $C^{2}$ solutions, $u_{+}(t,y)$, $t\geq 0$,
$u_{-}(t,y)$, $t\leq 0$, such that

$\lim_{t\to\pm\infty}|u_{\pm}(t,y)-w_{\pm}(y)|\to 0$

uniformly with respect to $y$, and

$h(u_{\pm}(t_{\ell},y),y)=c_{\ell}\Rightarrow\inf_{y\in\mathbb{R}^{m}}|h_{x}(u_{\pm}(t_{\ell},y),y)\dot{u}_{\pm}(t_{\ell}^{\pm},y)|>0.$
(1.10)

for any $-M\leq\ell\leq N$, $\ell\neq 0$.

Here we used the shorthand $u_{\pm}(t,y)-w_{\pm}(y)$ for
$u_{+}(t,y)-w_{+}(y)$, when $t\geq 0$, and $u_{-}(t,y)-w_{-}(y)$, when $t\leq
0$. We will use such a shorthand throughout the whole paper.

Note that (1.10) implies that $u_{\pm}(t,y)$ intersects transversally the set

${\cal S}_{\ell}(y)=\\{x\in\mathbb{R}^{n}\mid h(x,y)=c_{\ell}\\},$

at $u_{\pm}(t_{\ell},y)$. Then we look for solutions
$(x(t,\varepsilon),y(t,\varepsilon))$ of (1.7) such that

$\sup_{t\in\mathbb{R}_{\pm}}|x(t,\varepsilon)-u_{\pm}(t,y(t,\varepsilon))|\ll
1.$ (1.11)

Note that, from (1.10) it follows that, for any $y\in\mathbb{R}^{m}$, the set
${\cal S}_{\ell}(y)$, $\ell=-M,\ldots,N$, $\ell\neq 0$, is an hypersurface in
a ball around $u_{\pm}(t_{\ell},y)$ whose radius does not depend on $y$. We
emphasize that (1.10) is all that we need on ${\cal S}_{\ell}(y)$ for our
analysis.

To prove the existence of solutions of (1.7) satisfying (1.11) we use
Lyapunov-Schimtd method, together with a combination of singularly perturbed
analysis and a technique for discontinuous dynamical systems, to construct a
bifurcation function whose zeros are associated to solutions of the perturbed
equation whose $x$-component satisfies (1.11). This is the content of Theorem
6.2 where we prove that if a certain generic condition is satisfied then there
is a manifold of such solutions.

According to $A_{3})$ (see Section 2) $x=w_{\pm}(y)$ are normally hyperbolic
invariant manifolds for the unperturbed system $\dot{x}=f(x,y)$, $\dot{y}=0$.
From [1, 22] these manifolds perturb to invariant manifolds
$x=\tilde{w}_{\pm}(y,\varepsilon)$ such that

$\sup_{y\in\mathbb{R}^{m}}|\tilde{w}_{\pm}(y,\varepsilon)-w_{\pm}(y)|\to
0,\hbox{as $\varepsilon\to 0$.}$

From Remark 6.3 it follows that $x(t,\varepsilon)$ is asymptotic to the
manifold $\tilde{w}_{+}(y,\varepsilon)$ as $t\to\infty$ and to
$\tilde{w}_{-}(y,\varepsilon)$ as $t\to-\infty$. Hence,
$(x(t,\varepsilon),y(t,\varepsilon))$ behaves as a heteroclinic solution of
equation (1.7) connecting the invariant manifold
$\tilde{w}_{-}(y,\varepsilon)$ to $\tilde{w}_{+}(y,\varepsilon)$.

The problem studied in this paper has been motivated by [28], where existence
and bifurcation theorems are derived for homoclinic orbits in three-
dimensional flows that are perturbations of families of planar Hamiltonian
system. In this paper we study the problem of persistence of bounded solutions
in the discontinuous case (1.7) where $x\in\mathbb{R}^{n}$,
$y\in\mathbb{R}^{m}$, all functions considered are sufficiently smooth and
$\varepsilon\in\mathbb{R}$ is a small parameter, assuming the existence of
such an orbit in the unperturbed equation (1.9).

Then in Section 7 we apply Theorem 6.2 to extend [17, Theorem 1] to the
discontinuous equation (1.7). Following [6] we derive a bifurcation function
characterizing the persistence of homoclinic solutions of (1.7) from a generic
homoclinic solution of the unperturbed system (1.6). It can be easily checked
that the results of this paper easily extend if we replace $f(x,y)$ with
$f(x,y,\varepsilon)=f_{0}(x,y)+\varepsilon f_{1}(x,y,\varepsilon)$ (with the
same $h(x,y)$). In this case the unperturbed system will be

$\dot{x}=f_{0}(x,y),\;\dot{y}=0.$

The next step is the study of a degenerate case where, for any
$y\in\mathbb{R}^{m}$ the unperturbed discontinuous equation (1.6) has a
piecewise $C^{1}$ solution $u(t,y)$ heteroclinic to the hyperbolic fixed
points $x=w_{\pm}(y)$. We plan to perform this study in a forthcoming paper as
it is necessary to go into a deeper analysis of the bifurcation function.

We now briefly sketch the content of this paper. In Section 2 we provide basic
assumptions and define the piecewise smooth heteroclinic solution of the
unperturbed system. In Section 3 we recall the definition of exponential
dichotomy and extend this notion to discontinuous, piecewise linear, systems
with a jump at some points. We also extend to these systems some results
concerning existence of bounded solutions on either $t\geq 0$ and $t\leq 0$.
In our opinion, these results are theirselves interesting as they give the
form of the projection of the dichotomy of a linear discontinuous system.
Although this is an important point in the proof of Theorem 6.2, we think that
the proofs in this section can be skipped at a first reading.

Next, in Section 4, we construct families of bounded solutions and describe
them in terms of some parameters. These solutions are continuous and piecewise
smooth and give the bounded solutions we look for, when they assume the same
value at $t=0$. Then, after having defined the variational equation in Section
5, in Section 6 we study the _joining condition at $t=0$_ which is the
bifurcation condition and give a Melnikov-type condition assuring that the
bifurcation equation has a manifold of solutions. Finally, in Section 7 we
first state some general facts concerning two-dimensional discontinuous
equations depending on a slowly varying parameter and give an example of
application of the main result of this paper. In Section 8 we give a hint on a
possible extension of the example.

In the whole paper we will use the following notation. Given a vector $v$ or a
matrix $A$ with $v^{T}$, (resp. $A^{T}$) we denote the transpose of $v$ (resp.
$A$).

##  2 Notation and basic assumptions

Let $\Omega\subset\mathbb{R}^{n}$ be a bounded domain,

$c_{-M}<\ldots<c_{0}<\ldots<c_{N}$

be real numbers and $h:\Omega\times\mathbb{R}^{m}\to\mathbb{R}$ be a
$C^{r}$-functions, $r\geq 2$, with bounded derivatives. For
$\ell=-M,\ldots,N$, we set

$\Omega_{\ell}=\\{(x,y)\in\Omega\times\mathbb{R}^{m}\mid c_{i-1}\leq
h(x,y)<c_{i}\\},$

where we set for simplicity, $c_{-M-1}=-\infty$ and let
$f_{\ell}:\Omega\times\mathbb{R}^{m}\to\mathbb{R}^{n}$ be $C^{r}$-function,
bounded together with their derivatives in $\Omega\times\mathbb{R}^{m}$. We
are looking for solutions of equation

$\dot{x}=f_{\ell}(x,y),\quad(x,y)\in\Omega_{\ell}$ (2.1)

which are contained in a compact subset $H$ of $\Omega$. Hence it is not
restrictive to assume that $\Omega=\mathbb{R}^{n}$. So, from now on, we
suppose $\Omega=\mathbb{R}^{n}$.

First we give the definition of solutions of equation (2.1) we are considering
in this paper

######  Definition 2.1.

A continuous, piecewise smooth function $u(t,y)$ is a solution of equation
(2.1) on $t\geq 0$ intersecting transversally the sets ${\cal
S}_{i}(y)=\\{x\in\Omega\mid h(x,y)=c_{i}\\}$, $i=0,\ldots,N-1$, if there exist
$\eta>0$ and $0<t_{1}(y)<\ldots<t_{N}(y)$ such that the following conditions
hold for $1\leq i\leq N-1$ (note that we set $t_{0}(y)=0$)

  * $a_{1})$

$\dot{u}(t,y)=f_{i-1}(u(t,y),y)$ for $t_{i-1}(y)<t<t_{i}(y)$ and
$\dot{u}(t,y)=f_{N}(u(t,y),y)$ for $t>t_{N}(y)$;

  * $a_{2})$

$h(u(t_{i}(y),y),y)=c_{i-1},\quad\hbox{and}\quad
h_{x}(u(t_{i}(y),y),y)\dot{u}(t_{i}(y)^{\pm},y)>2\eta$;

  * $a_{3})$

$c_{i-1}<h(u(t,y),y)<c_{i}$, for $t_{i}(y)<t<t_{i+1}(y)$ and
$h(u(t,y),y)>c_{N-1}$, for $t>t_{N}(y)$.

Similarly, a continuous, piecewise smooth function $u(t,y)$ is a solution of
equation (2.1) on $t\leq 0$ intersecting transversally the sets ${\cal
S}_{j}(y)=\\{x\in\Omega\mid h(x,y)=c_{j}\\}$, $j=-1,\ldots,-M$, if there exist
$\eta>0$ and $t_{-M}(y)<\ldots<t_{-1}(y)<0$ such that the following conditions
hold for any $-M+1\leq j\leq-1$:

  * $a^{\prime}_{1})$

$\dot{u}(t,y)=f_{j+1}(u(t,y),y)$ for $t_{j}(y)<t<t_{j+1}(y)$ and
$\dot{u}(t,y)=f_{-M}(u(t,y),y)$ for $t<t_{-M}(y)$;

  * $a^{\prime}_{2})$

$h(u(t_{j}(y),y),y)=c_{j},\quad\hbox{and}\quad
h_{x}(u(t_{j}(y),y),y)\dot{u}(t_{j}(y)^{\pm},y)>2\eta$;

  * $a^{\prime}_{3})$

$c_{j}<h(u(t,y),y)<c_{j+1}$, for $t_{j}(y)<t<t_{j+1}(y)$ and
$h(u(t,y),y)<c_{-M}$, for $t<t_{-M}(y)$.

In this paper we assume that continuous, piecewise smooth solution
$u_{+}(t,y)$ and $u_{-}(t,y)$ of equation (2.1) exist, for $t\geq 0$, resp.
$t\leq 0$, such that the following conditions also hold.

  * $A_{1})$

$w_{0}^{\pm}(y):=u_{\pm}(0,y)$ and their derivatives are bounded functions and
belong to an open and bounded subset $B\subset\mathbb{R}^{n}$ such that
$\overline{B}\times\mathbb{R}^{m}\subset\Omega_{0}$.

  * $A_{2})$

There exist smooth and bounded functions $w_{\pm}(y)$ and $\mu_{0}>0$, such
that

$\begin{array}[]{l}f_{N}(w_{+}(y),y)=f_{-M}(w_{-}(y),y)=0\\\
h(w_{+}(y),y)-c_{N-1}>\mu_{0},\quad h(w_{-}(y),y)-c_{-M}<-\mu_{0}\end{array}$

for any $y\in\mathbb{R}^{m}$ and

$\lim_{t\to\pm\infty}u_{\pm}(t,y)-w_{\pm}(y)=0$

uniformly with respect to $y\in\mathbb{R}^{m}$.

  * $A_{3})$

For any $y\in\mathbb{R}^{m}$, $f_{N,x}(w_{+}(y),y)$ and $f_{-M,x}(w_{-}(y),y)$
have $k$ eigenvalues with negative real parts and $n-k$ eigenvalues with
positive real parts, counted with multiplicities and there exists
$\delta_{0}>0$ such that all these eigenvalues satisfy

$|Re\lambda(y)|>\delta_{0}.$

  * $A_{4})$

There exists $y_{0}\in\mathbb{R}^{m}$ such that
$w_{0}^{-}(y_{0})=w_{0}^{+}(y_{0})=x_{0}$.

Figure 1: The piecewise continuous bounded solution of (1.2). For simplicity
we write $w^{\pm}_{j}$ instead of $w^{\pm}_{j}(y)$. The solutions
$u_{\pm}(t,y)$ go from $\Omega_{-M}$ to $\Omega_{N}$ as $t$ goes from
$-\infty$ to $\infty$ and may have a discontinuity at $t=0$. In this figure we
actually draw the intersection of $\Omega_{i}$ with $y=$constant.

######  Remark 2.2.

From $A_{2})-A_{4})$ it follows that $u(t,y_{0})$ is a continuous, piecewice
$C^{1}$ solution of $\dot{x}=f(x,y_{0})$ such that

$\lim_{t\to\pm\infty}u(t,y_{0})=w_{\pm}(y_{0}).$

Moreover, from $a_{2})-a_{3})$ it follows that, for $i=1,\ldots,N$:

$\begin{array}[]{l}\frac{\partial}{\partial
t}h(u_{+}(t,y),y)_{|t=t_{i}(y)^{\pm}}\geq 0\end{array}$

that is

$\begin{array}[]{l}h_{x}(u_{+}(t_{i}(y),y)\dot{u}_{+}(t_{i}(y)^{\pm},y)\geq
0,\end{array}$

Similarly we see that

$\begin{array}[]{l}h_{x}(u_{-}(t_{i}(y),y)\dot{u}_{-}(t_{i}(y)^{\pm},y)\geq
0,\end{array}$

So $a_{2})-a^{\prime}_{2})$ are a kind of transversality assumption on the
solutions $u_{\pm}(t,y)$.

ii) All results in this paper can be generalised to the case where the
solutions exit transversally $\Omega_{i}$ and enter into either $\Omega_{i+1}$
or $\Omega_{i-1}$ transversally in the sense that (1.10) holds at the
intersection points of the solution with the boundary of $\Omega_{i}$. More
precisely suppose, to fix ideas, $t\geq 0$. Then we assume the following (see
Fig. 2). There exists $(i_{0}=0,i_{1}\ldots,i_{M})$ such that given $i_{h}$
then $i_{h+1}$ is either $i_{h}-1$ or $i_{h}+1$ and for
$t_{h}(y)<t<t_{h+1}(y)$ we have

$c_{i_{h}-1}<h(u_{+}(t,y),y)<c_{i_{h}},\quad\hbox{for any $h=0,\ldots,N-1$}.$

Moreover

$\begin{array}[]{l}|h_{x}(u_{+}(t_{i}(y),y),y)\dot{u}_{+}(t_{i}(y)^{\pm},y)|>2\eta\end{array}$

for any $h=1,\ldots,N$.

A similar generalization can be made for $t\leq 0$ and all other assumptions
are changed accordingly.

Another possible generalization is the homoclinic case. This is the case, for
example, where $N=M$ and assumptions $a^{\prime}_{1})-a^{\prime}_{3})$ are
changed to the following where we write $i\in\\{1,\ldots,N\\}$ instead of $-j$
(see Fig. 3):

  * $a^{\prime}_{1})$

$\dot{u}(t,y)=f_{i-1}(u(t,y),y)$ for $t_{-i}(y)<t<t_{-(i-1)}(y)$ and
$\dot{u}(t,y)=f_{N}(u(t,y),y)$ for $t<t_{-N}(y)$;

  * $a^{\prime}_{2})$

$h(u(t_{-i}(y),y),y)=c_{i-1}$ and
$h_{x}(u(t_{-i}(y),y),y)\dot{u}(t_{-i}(y)^{\pm},y)<-2\eta$;

  * $a^{\prime}_{3})$

$c_{i-2}<h(u(t,y),y)<c_{i-1}$, $t_{-i}(y)<t<t_{-(i-1)}(y)$ and
$h(u(t,y),y)>c_{N-1}$, for $t<t_{-N}(y)$.

where, we set for simplicity, $c_{-1}=-\infty$.

Figure 2: An example of the situation described in Remark 2.2 with $N=3$ and
$M=4$. Here we have
$(i_{-1},i_{-2},i_{-3},i_{-4},i_{-5},i_{-6},i_{-7})=(-1,-2,-1,-2,-3,-4,-3)$
and $(i_{0},i_{1},i_{2},i_{3},i_{4},i_{5},i_{6},i_{7})=(0,1,2,1,0,1,2,3)$.
Again in this figure we draw the intersection of $\Omega_{i}$ with
$y=$constant. Figure 3: An example of homoclinic orbit as described in Remark
2.2. The homoclinic orbit arises at $y=y_{0}$ when
$w^{+}_{0}(y_{0})=w^{-}_{0}(y_{0})$. Again in this figure we draw the
intersection of $\Omega_{i}$ with $y=$constant

We have the following

######  Lemma 2.3.

For $i=0\ldots,N$, $j=0,\ldots,-M$, let $w_{i}^{+}(y)=u_{+}(t_{i}(y),y)$,
$w^{-}_{j}(y)=u_{-}(t_{j}(y),y)$. Then $w_{i}^{+}(y)$ and $w^{-}_{j}(y)$ are
$C^{r}$-functions bounded together with their derivatives.

###### Proof.

We know that $w_{0}^{\pm}(y)$ is $C^{r}$ and bounded together with its
derivatives. Now,

$w_{1}^{+}(y)=u_{+}(t_{1}(y),y)=w_{0}^{+}(y)+\int_{0}^{t_{1}(y)}f_{1}(u_{+}(t,y),y)dt.$

Hence $w_{1}^{+}(y)$ is $C^{r}$ and bounded since so are $w_{0}^{+}(y)$,
$t_{1}(y)$ and $f_{1}(x,y)$. Next, from

$h_{x}(w_{1}^{+}(y),y)f_{1}(w_{1}^{+}(y),y)>2\eta$

and the fact that $f_{1}(x,y)$ is bounded, we see that $\bar{h}>0$ exists such
that $\inf_{y\in\mathbb{R}^{m}}|h_{x}(w_{1}^{+}(y),y)|>\bar{h}$. Then,
differentiating $h(w_{1}^{+}(y),y)=c_{0}$ we see that

$\frac{dw_{1}^{+}}{dy}(y)=-\frac{h_{y}(w_{1}^{+}(y),y)}{h_{x}(w_{1}^{+}(y),y)}$

is bounded. More arguments of similar nature show that all derivatives of
$w_{1}^{+}(y)$ are bounded.

Suppose now that $w^{+}_{i}(y)$ is bounded with its derivative. Then, since

$w^{+}_{i+1}(y)=w^{+}_{i}(y)+\int_{t_{i}(y)}^{t_{i+1}(y)}f_{i}(u_{+}(t,y),y)dt$

we see that $w^{+}_{i+1}(y)$ is a bounded function and a similar argument as
before shows that the derivatives of $w^{+}_{i+1}(y)$ are also bounded. The
proof that $w^{-}_{j}(y)$ is a $C^{r}$ function bounded with its derivatives
for any $j=-1,\ldots,-M$ is similar. ∎

Let $i=0,1,\ldots,N$. For $t\geq 0$, let $u_{i}^{+}(t,y)$ be the solution of
$\dot{x}=f_{i}(x,y)$ such that $u_{i}^{+}(t_{i}(y),y)=w^{+}_{i}(y)$.
Similarly, for $j=-M,\ldots,-1,0$ let $u_{j}^{-}(t,y)$ be the solution of
$\dot{x}=f_{j}(x,y)$ such that $u_{j}^{-}(t_{i}(y),y)=w^{-}_{j}(y)$. Then

$\begin{array}[]{l}u(t,y):=\left\\{\begin{array}[]{ll}u_{+}(t,y)&\hbox{$t\geq
0$}\\\ u_{-}(t,y)&\hbox{for $t\leq 0$}\end{array}\right.\\\
\qquad=\left\\{\begin{array}[]{ll}u_{j}^{-}(t,y)&\hbox{for $t_{j-1}(y)\leq
t\leq t_{j}(y)$, $j=0,-1,\ldots-M$}\\\ u_{i}^{+}(t,y)&\hbox{for $t_{i}(y)\leq
t\leq t_{i+1}(y)$, $i=0,\ldots,N$}\end{array}\right.\end{array}$ (2.2)

where, for simplicity, we set $t_{N+1}(y)=\infty$ and $t_{-M-1}(y)=-\infty$

$u_{i}^{+}(t_{i+1}(y),y)=u_{+}(t_{i+1}(y),y)=w^{+}_{i+1}(y)=u_{i+1}^{+}(t_{i}(y),y)$

and similarly,

$u_{j}^{-}(t_{j-1}(y),y)=w^{-}_{j-1}(y)=u^{-}_{j-1}(t_{j-1}(y),y).$

Here we emphasize a convention that we will use throughout the whole paper.
When we use the index $i$, such as in $t_{i}$ or $w^{+}_{i}(y)$, we always
mean $i=1,\ldots,N$, (sometimes also $i=0$) while, when we use the index $j$,
as in $t_{j}$ or $w^{-}_{j}(y)$, we always mean $j=-1,\ldots,-M$ (sometimes
also $j=0$).

##  3 Exponential dichotomy for piecewise discontinuous systems

A basic tool in this paper is the notion of exponential dichotomy, whose
definition we recall here. Let $J$ be either $[a,\infty),(-\infty,a]$, or
$\mathbb{R}$ and $A(t)$, $t\in J$, be a $n\times n$ continuous matrix. We say
that the linear system

$\dot{x}=A(t)x,\quad x\in\mathbb{R}^{n}$ (3.1)

has an exponential dichotomy on $J$ if there exist a projection
$P:\mathbb{R}^{n}\to\mathbb{R}^{n}$ and constants $\delta>0$ and $K\geq 1$
such that the fundamental matrix $X(t)$ of (3.1) satisfying $X(a)=\mathbb{I}$,
when $J=[a,\infty),(-\infty,a]$, or $X(0)=\mathbb{I}$ when $J=\mathbb{R}$,
satisfies

$\begin{array}[]{l}|X(t)PX(s)^{-1}|\leq Ke^{-\delta(t-s)},\quad\hbox{for
$s\leq t,\;s,t\in J$}\\\ |X(s)(\mathbb{I}-P)X(t)^{-1}|\leq
Ke^{-\delta(t-s)},\quad\hbox{for $s\leq t,\;s,t\in J$.}\end{array}$

In this section we extend the definition of exponential dichotomy to systems
with discontinuities. To allow more generality we consent the solutions to
have jump discontinuities at the discontinuity points of the coefficient
matrix.

So, let $t_{0}<t_{1}<\ldots<t_{N}$ be real numbers, $B_{1},\ldots,B_{N}$ be
invertible $n\times n$ matrices and ${\cal A}(t)$, $t\geq t_{0}$ be a
piecewise continuous matrix with possible discontinuity jumps at
$t=t_{1},\ldots,t_{N}$, that is

${\cal A}(t)=\left\\{\begin{array}[]{ll}A_{i-1}(t)&\hbox{if $t_{i-1}\leq
t<t_{i}$,}\\\ A_{N}(t)&\hbox{if $t\geq t_{N}$}\end{array}\right.$ (3.2)

where $A_{0}(t),\ldots,A_{N}(t)$ are continuous matrices. Note that ${\cal
A}(t)$ is continuous for $t\geq t_{0}$, $t\neq t_{i}$, and right-continuous at
$t=t_{i}$, with possible jumps at $t=t_{i}$ given by the matrix
$A_{i}(t_{i})-A_{i-1}(t_{i})$. For $t\geq t_{0}$ we consider the linear,
discontinuous, system

$\begin{array}[]{l}\dot{x}={\cal A}(t)x\\\
x(t_{i}^{+})=B_{i}x(t_{i}^{-}),\end{array}$ (3.3)

Similarly, if $t_{-M}<\ldots t_{-1}<t_{0}$ and

${\cal A}(t)=\left\\{\begin{array}[]{ll}A_{-M}(t)&\hbox{if $t\leq t_{-M}$}\\\
A_{j+1}(t)&\hbox{if $t_{j}<t\leq t_{j+1}$,}\end{array}\right.$ (3.4)

where $A_{0}(t),A_{-1}(t),\ldots A_{-M}(t)$ are continuous matrices, we
consider, for $t\leq t_{0}$, the linear, discontinuous, system

$\begin{array}[]{l}\dot{x}={\cal A}(t)x\\\
x(t_{j}^{+})=B_{j}x(t_{j}^{-})\end{array}$ (3.5)

Note that ${\cal A}(t)$ is continuous for $t\leq t_{0}$, $t\neq t_{j}$, left-
continuous at $t=t_{j}$, with possible jumps at $t=t_{j}$, given by the matrix
$A_{j+1}(t_{j})-A_{j}(t_{j})$.

######  Remark 3.1.

i) As a matter of facts, for $t\geq t_{0}$, we will consider

${\cal A}(t)=\left\\{\begin{array}[]{ll}A_{i-1}(t)&\hbox{if $t_{i-1}\leq t\leq
t_{i}$,}\\\ A_{N}(t)&\hbox{if $t\geq t_{N}$}\end{array}\right.$

and similarly for $t\leq t_{0}$. This may cause a duplicate definition of
${\cal A}(t)$ at $t=t_{i}$, however it will be always clear which among
$A_{i}(t)$ will be taken into account at that point.

ii) The results of this section will be applied to the linear system
$\dot{x}=A(t,y)x$ where $A(t,y)$ is given by

$A(t,y):=\left\\{\begin{array}[]{ll}f_{-M,x}(u(t,y),y)&\hbox{if $t\leq
t_{-M}(y)$.}\\\ f_{j+1,x}(u(t,y),y)&\hbox{if $t_{j}(y)<t\leq t_{j+1}(y)$}\\\
f_{i-1,x}(u(t,y),y)&\hbox{if $t_{i-1}(y)\leq t<t_{i}(y)$}\\\
f_{N,x}(u(t,y),y)&\hbox{if $t\geq t_{N}(y)$.}\end{array}\right.$ (3.6)

Note that, being $u(t,y)$ continuous for $t\neq 0$, $A(t,y)$ is continuous for
$t\neq 0$, $t\neq t_{i}(y),t_{j}(y)$ with jump discontinuities at
$t=0,t_{i}(y),t_{j}(y)$. More precisely

$\begin{array}[]{l}A(0^{-},y)=f_{0,x}(w_{0}^{-}(y),y)\\\
A(0^{+},y)=f_{0,x}(w_{0}^{+}(y),y),y)\\\
A(t_{i}(y)^{-},y)=f_{i-1,x}(w_{i}^{+}(y)),y)\\\
A(t_{i}(y)^{+},y)=f_{i,x}(w_{i}^{+}(y)),y)\\\
A(t_{j}(y)^{+},y)=f_{j+1,x}(w^{-}_{i}(y),y)\\\
A(t_{j}(y)^{-},y)=f_{j,x}(w^{-}_{i}(y),y).\end{array}$

Without loss of generality we may assume that $t_{0}=0$, so in the remaining
part of this section we will take $t_{0}=0$.

Let $U_{i}(t)$ be the fundamental matrix of the linear systems

$\dot{x}=A_{i}(t)x$

on $\mathbb{R}$, that is $\dot{U}_{i}(t)=A_{i}(t)U_{i}(t)$, $t\in\mathbb{R}$,
and $U_{i}(0)=\mathbb{I}$. The fundamental matrix of (3.3), where $t\geq 0$,
is given by the (discontinuous) invertible matrix

$X_{+}(t)=\left\\{\begin{array}[]{ll}U_{0}(t)&\hbox{if $0\leq t<t_{1}$}\\\
U_{i}(t)U_{i}(t_{i})^{-1}B_{i}X_{+}(t_{i}^{-})&\hbox{if $t_{i}\leq
t<t_{i+1}$}\\\ U_{N}(t)U_{N}(t_{N})^{-1}B_{N}X_{+}(t_{N}^{-})&\hbox{if $t\geq
t_{N}$}\end{array}\right.$ (3.7)

that is

$X_{+}(t)=\left\\{\begin{array}[]{ll}U_{0}(t)&\hbox{for $0\leq t<t_{1}$}\\\
U_{1}(t)U_{1}(t_{1})^{-1}B_{1}U_{0}(t_{1})&\hbox{for $t_{1}\leq t<t_{2}$}\\\
U_{i}(t)U_{i}(t_{i})^{-1}B_{i}U_{i-1}(t_{i})U_{i-1}(t_{i-1})^{-1}\ldots
B_{1}U_{0}(t_{1})&\hbox{for $t_{i}\leq t<t_{i+1}$}\\\ \vdots\\\
U_{N}(t)U_{N}(t_{N})^{-1}B_{N}U_{N-1}(t_{N})U_{N-1}(t_{N-1})^{-1}\ldots
B_{1}U_{0}(t_{1})&\hbox{for $t\geq t_{N}$.}\end{array}\right.$

Similarly, the fundamental matrix of (3.3), where $t\leq 0$, is given by the
(discontinuous) invertible matrix

$X_{-}(t)=\left\\{\begin{array}[]{ll}U_{0}(t)&\hbox{if $t_{-1}<t\leq 0$}\\\
U_{j}(t)U_{j}(t_{j})^{-1}B_{j}^{-1}X_{-}(t_{j}^{+})&\hbox{if $t_{j-1}<t\leq
t_{j}$}\\\ U_{-M}(t)U_{-M}(t_{-M})^{-1}B_{-M}^{-1}X_{-}(t_{-M}^{+})&\hbox{if
$t\leq t_{-M}.$}\end{array}\right.$ (3.8)

Note that $X_{+}(t)$ is continuous for $t\neq t_{1},\ldots,t_{N}$ and right-
continuous at $t=t_{1},\ldots,t_{N}$ and $X_{-}(t)$ is continuous for $t\neq
t_{-1},\ldots,t_{-M}$ and left-continuous at $t\neq t_{-1},\ldots,t_{-M}$.

It is clear that $\dot{X}_{+}(t)={\cal A}(t)X_{+}(t)$, for $t\geq 0$, $t\neq
t_{i}$, $\dot{X}_{-}(t)={\cal A}(t)X_{-}(t)$, for $t\leq 0$, $t\neq t_{j}$,
$X_{\pm}(0)=\mathbb{I}$, the identity matrix, and

$\begin{array}[]{l}X_{+}(t_{i}^{+})=B_{i}X_{+}(t_{i}^{-})\\\
X_{-}(t_{j}^{+})=B_{j}X_{-}(t_{j}^{-}).\end{array}$

Actually we can write

$X_{+}(t_{i})=B_{i}X_{+}(t_{i}^{-}),\quad
X_{-}(t_{j})=B_{j}^{-1}X_{-}(t_{j}^{+})$

since $X_{+}(t)$ is right-continuous and $X_{-}(t)$ is left-continuous. Then,
from (3.7) and (3.8) we see that

$X_{+}(t)X_{+}(t_{N}^{+})=U_{N}(t)U_{N}(t_{N})^{-1},\forall t\geq t_{N}$ (3.9)

and, similarly,

$X_{-}(t)X_{-}(t_{N}^{-})=U_{N}(t)U_{N}(t_{N})^{-1},\forall t\leq t_{-M}.$
(3.10)

######  Remark 3.2.

i) Let $\tau\geq 0$ be a fixed number. For $t\geq 0$,
$x(t)=X_{+}(t)X_{+}(\tau)^{-1}\tilde{x}$ is the right-continuous solution of

$\left\\{\begin{array}[]{l}\dot{x}={\cal A}(t)x,\quad\hbox{for $t\geq 0$,
$t\neq t_{1},\ldots,t_{N}$}\\\ x(t_{i}^{+})=B_{i}x(t_{i}^{-})\\\
x(\tau^{+})=\tilde{x}.\end{array}\right.$ (3.11)

Indeed, it is obvious that $\dot{x}(t)={\cal A}(t)x(t)$ for $t\geq 0$, $t\neq
t_{1},\ldots,t_{N}$ and that $x(t_{i}^{+})=B_{i}x(t_{i}^{-})$, since
$X_{+}(t_{i}^{+})=B_{i}X_{+}(t_{i}^{-})$. Moreover
$x(\tau^{+})=X_{+}(\tau^{+})X_{+}(\tau)^{-1}\tilde{x}=X_{+}(\tau)X_{+}(\tau)^{-1}\tilde{x}=\tilde{x}$,
since $X_{+}(t)$ is right-continuous at any $t\geq 0$.

Similarly, for $t\leq 0$ and any fixed $\tau\leq 0$,
$x(t)=X_{-}(t)X_{-}(\tau)^{-1}\tilde{x}$ is the left-continuous solution of

$\left\\{\begin{array}[]{l}\dot{x}={\cal A}(t)x,\quad\hbox{for $t\leq 0$,
$t\neq t_{-1},\ldots,t_{-N}$}\\\ x(t_{-i}^{-})=B_{i}^{-1}x(t_{-i}^{+})\\\
x(\tau^{-})=\tilde{x}.\end{array}\right.$ (3.12)

ii) As $X_{+}(t)X_{+}(\tau)^{-1}\tilde{x}$ is right-continuous at
$t=t_{1},\ldots,t_{N}$ it is also clear that
$X_{+}(t)X_{+}(\tau)^{-1}\tilde{x}$ satisfies

$\left\\{\begin{array}[]{l}\dot{x}={\cal A}(t)x,\quad\hbox{for $t\geq 0$,
$t\neq t_{1},\ldots,t_{N}$}\\\ x(t_{i})=B_{i}x(t_{i}^{-})\\\
x(\tau)=\tilde{x}\end{array}\right.$ (3.13)

and, similarly, $X_{-}(t)X_{-}(\tau)^{-1}\tilde{x}$ satisfies

$\left\\{\begin{array}[]{l}\dot{x}={\cal A}(t)x,\quad\hbox{for $t\leq 0$,
$t\neq t_{1},\ldots,t_{N}$}\\\ x(t_{-i})=B_{i}^{-1}x(t_{-i}^{+})\\\
x(\tau)=\tilde{x}.\end{array}\right.$ (3.14)

We have the following

######  Lemma 3.3.

Suppose that the linear system

$\dot{x}=A_{N}(t)x\quad[resp.\dot{x}=A_{-M}(t)x]$

has an exponential dichotomy on $t\geq t_{N}$ (resp. $t\leq t_{-M}$) with
constant $K$, exponent $\delta$ and projection ${\cal P}_{+}$ (resp. ${\cal
P}_{-}$). Then, the linear system (3.3) (resp. (3.5)) with ${\cal A}(t)$ as in
(3.2) (resp. (3.4)) has an exponential dichotomy on $\mathbb{R}_{+}$, (resp.
$\mathbb{R}_{-}$) with the same exponent $\delta$, constant $\tilde{K}\geq K$
and projection

$\begin{array}[]{c}\tilde{\cal P}_{+}=X_{+}(t_{N})^{-1}{\cal
P}_{+}X_{+}(t_{N})\\\ (\hbox{resp. }\tilde{\cal P}_{-}=X_{-}(t_{-M})^{-1}{\cal
P}_{-}X_{-}(t_{-M})).\end{array}$ (3.15)

###### Proof.

As $X_{+}(t)$ is right continuous at $t=t_{N}$, for $t_{N}\leq s\leq t$ we
have

$\begin{array}[]{l}|X_{+}(t)\tilde{\cal P}_{+}X_{+}(s)^{-1}|=\\\
|U_{N}(t)U_{N}(t_{N})^{-1}B_{N}X_{+}(t_{N}^{-})X_{+}(t_{N}^{+})^{-1}{\cal
P}_{+}X_{+}(t_{N}^{+})X_{+}(t_{N}^{-})^{-1}B_{N}^{-1}U_{N}(t_{N})U_{N}(s)^{-1}|\\\
=|U_{N}(t)U_{N}(t_{N})^{-1}{\cal P}_{+}U_{N}(t_{N})U_{N}(s)^{-1}|\leq
Ke^{-\delta(t-s)}\end{array}$

since $U_{N}(t)U_{N}(t_{N})^{-1}$ is the fundamental matrix of
$\dot{x}=A_{N}(t)x$ on $t\geq t_{N}$. Similarly we see that

$|X_{+}(s)(\mathbb{I}-\tilde{\cal P}_{+})X_{+}(t)^{-1}|\leq Ke^{-\delta(t-s)}$

for $t_{N}\leq s\leq t$. Being $0\leq s\leq t\leq t_{N}$ compact and
$X_{+}(t)\tilde{\cal P}_{+}X_{+}(s)^{-1}$ piecewise continuous with right and
left limits at the discontinuity points, there exists $\hat{K}\geq K$ such
that

$\begin{array}[]{l}|X_{+}(t)\tilde{\cal
P}_{+}X_{+}(s)^{-1}|\leq\hat{K}e^{-\delta(t-s)}\\\
|X_{+}(s)(\mathbb{I}-\tilde{\cal
P}_{+})X_{+}(t)^{-1}|\leq\hat{K}e^{-\delta(t-s)}\end{array}$

for $0\leq s\leq t\leq t_{N}$. Finally, for $0\leq s\leq t_{N}<t$, we have,
using right-continuity of $X_{+}(t)$:

$\begin{array}[]{l}|X_{+}(t)\tilde{\cal
P}_{+}X_{+}(s)^{-1}|\leq|X_{+}(t)\tilde{\cal
P}_{+}X_{+}(t_{N})^{-1}||B_{N}X_{+}(t_{N}^{-})\tilde{\cal
P}_{+}X_{+}(s)^{-1}|\\\ \qquad\leq
Ke^{-\delta(t-t_{N})}|B_{N}|\hat{K}e^{-\delta(t_{N}-s)}\leq
K|B_{N}|\hat{K}e^{-\delta(t-s)}\\\ |X_{+}(s)(\mathbb{I}-\tilde{\cal
P}_{+})X_{+}(t)^{-1}|\\\ \qquad\leq|X_{+}(s)(\mathbb{I}-\tilde{\cal
P}_{+})X_{+}(t_{N}^{-})^{-1}||B_{N}^{-1}X_{+}(t_{N}^{+})(\mathbb{I}-\tilde{\cal
P}_{+})X_{+}(t)^{-1}|\\\ \qquad\leq
Ke^{-\delta(s-t_{N})}|B_{N}^{-1}|\hat{K}e^{-\delta(t_{N}-t)}\leq
K|B_{N}^{-1}|\hat{K}e^{-\delta(t-s)}\end{array}$

proving the result for $t\geq 0$ with
$\tilde{K}=\max\\{K,\hat{K},K|B_{N}^{-1}|\hat{K}\\}$. A similar argument works
when $t\leq 0$. The Lemma is proved. ∎

The following result characterises ${\cal R}\tilde{P}_{+}$, resp. ${\cal
N}\tilde{P}_{-}$, in terms of bounded solutions of system (3.3), resp. (3.5),
extending to the discontinuous case a simular result for continuous equations.

######  Lemma 3.4.

Let ${\cal A}(t)$ be either as in (3.2) or (3.4). Suppose that the condition
of Lemma 3.3 hold and let $\tilde{\cal P}_{\pm}$ be as in (3.15). Then
$\xi_{+}\in{\cal R}\tilde{\cal P}_{+}$ if and only if the solution of the
discontinuous system (3.3) such that $x(0)=\xi_{+}$ is bounded for $t\geq 0$.
Similarly, $\xi_{-}\in{\cal N}\tilde{\cal P}_{-}$ if and only if the solution
of the discontinuous system (3.3) such that $x(0)=\xi_{-}$ is bounded for
$t\leq 0$.

###### Proof.

If $\xi_{+}\in{\cal R}\tilde{\cal P}_{+}$ we have

$|X_{+}(t)\xi_{+}|=|X_{+}(t)\tilde{\cal P}_{+}\xi_{+}|\leq Ke^{-\delta
t}|\xi_{+}|$

that is the solution of (3.3) starting from $\xi_{+}$ is bounded. Vice versa,
suppose that $x(t)$ is a solution of (3.3) bounded on $t\geq 0$. We have

$\begin{array}[]{l}|[\mathbb{I}-\tilde{\cal
P}_{+}]x(0)|=|[\mathbb{I}-\tilde{\cal P}_{+}]X_{+}(t)^{-1}x(t)|\leq
Ke^{-\delta t}\sup_{t\geq 0}|x(t)|\to 0\end{array}$

as $t\to\infty$. Then $[\mathbb{I}-\tilde{\cal P}_{+}]x(0)=0$ and hence
$x(0)\in{\cal R}\tilde{\cal P}_{+}(0)$. By a similar argument we prove the
thesis when $\xi_{-}\in{\cal N}\tilde{\cal P}_{-}$ is concerned. ∎

We conclude this section with the following

######  Lemma 3.5.

Let $B_{i}$, $B_{j}$, be invertile $n\times n$ matrices and $k(t)$ be a
bounded integrable function for $t\geq 0$, (resp. $t\leq 0$). Suppose the
condition of Lemma 3.3 hold and set

$\begin{array}[]{l}\tilde{\cal P}_{+}^{\tau}=X_{+}(\tau)\tilde{\cal
P}_{+}X_{+}(\tau)^{-1}\\\ \tilde{\cal P}_{-}^{\tau}=X_{-}(-\tau)\tilde{\cal
P}_{-}X_{-}(-\tau)^{-1}\end{array}$

where $\tilde{\cal P}_{\pm}$ is as in (3.15) and $0\leq\tau\in\mathbb{R}$ is a
fixed number. Then, for any $\xi_{+}\in{\cal R}\tilde{\cal P}_{+}^{\tau}$
(resp. $\xi_{-}\in{\cal N}\tilde{\cal P}_{-}^{\tau}$) the linear inhomogeneous
system

$\begin{array}[]{l}\dot{x}={\cal A}(t)x+k(t)\\\
x(t_{i}^{+})=B_{i}x(t_{i}^{-}),\\\ \tilde{\cal
P}_{+}^{\tau}x(\tau)=\xi_{+}\end{array}$ (3.16)

with $t\geq 0$, [resp.

$\begin{array}[]{l}\dot{x}={\cal A}(t)x+k(t)\\\
x(t_{j}^{-})=B_{j}^{-1}x(t_{j}^{+})\\\ (\mathbb{I}-\tilde{\cal
P}_{-}^{\tau})x(-\tau)=\xi_{-}\end{array}$

when $t\leq 0$] has the unique right-continuous, [resp. left-continuous when
$t\leq 0$] bounded solution

$\begin{array}[]{l}\displaystyle x(t)=X_{+}(t)\tilde{\cal
P}_{+}X_{+}(\tau)^{-1}\xi_{+}+\int_{\tau}^{t}X_{+}(t)\tilde{\cal
P}_{+}X_{+}(s)^{-1}k(s)ds\\\
\displaystyle\qquad-\int_{t}^{\infty}X_{+}(t)(\mathbb{I}-\tilde{\cal
P}_{+})X_{+}(s)^{-1}k(s)ds\end{array}$ (3.17)

[resp.

$\begin{array}[]{l}\displaystyle x(t)=X_{-}(t)(\mathbb{I}-\tilde{\cal
P}_{-})X_{-}(-\tau)^{-1}\xi_{-}+\int_{-\infty}^{t}X_{-}(t)\tilde{\cal
P}_{-}X_{-}(s)^{-1}k(s)ds\\\
\displaystyle\qquad-\int_{t}^{-\tau}X_{-}(t)(\mathbb{I}-\tilde{\cal
P}_{-})X_{-}(s)^{-1}k(s)ds\end{array}$ (3.18)

if $t\leq 0$]. Moreover such a solution satisfies

$\sup_{t\geq\tau}|x(t)|\leq K[|\xi_{+}|+2\delta^{-1}\sup_{t\geq 0}|k(t)|]$
(3.19)

if $t\geq 0$ [resp.

$\sup_{t\leq-\tau}|x(t)|\leq K[|\xi_{-}|+2\delta^{-1}\sup_{t\leq 0}|k(t)|]$
(3.20)

if $t\leq 0$].

###### Proof.

We only give the proof for $t\geq 0$, the proof for $t\leq 0$ being similar.
We prove uniqueness, first. Suppose that $x_{1}(t),x_{2}(t)$ are two solutions
of (3.16), right-continuous and bounded for $t\geq 0$. Then
$x(t)=x_{1}(t)-x_{2}(t)$ is a right-continuous, bounded solution of

$\begin{array}[]{l}\dot{x}={\cal A}(t)x\\\ x(t_{i}^{+})=B_{i}x(t_{i}^{-})\\\
\tilde{\cal P}_{+}^{\tau}x(\tau)=0.\end{array}$

Then, as we have observed in Remark 3.2-i),
$x(t)=X_{+}(t)X_{+}(\tau)^{-1}x(\tau)$, so:

$\begin{array}[]{l}x(\tau)=(\mathbb{I}-\tilde{\cal
P}_{+}^{\tau})x(\tau)=(\mathbb{I}-\tilde{\cal
P}_{+}^{\tau})X_{+}(\tau)X_{+}(t)^{-1}x(t)\\\
\qquad=X_{+}(\tau)(\mathbb{I}-\tilde{\cal P}_{+})X_{+}(t)^{-1}x(t)\end{array}$

from which we get

$|x(\tau)|\leq Ke^{-\delta(t-\tau)}\sup_{t\geq\tau}|x(t)|\to 0$

as $t\to\infty$, since $x(t)$ is bounded. This proves that $x(\tau)=0$ and
then $x(t)=0$. Hence we have uniqueness. To show the existence we observe that
the function given in (3.17) is right-continuous, bounded for $t\geq\tau$ (and
hence also for $t\geq 0$) and satisfies (3.16). Finally, from (3.17) we get:

$|x(t)|\leq
Ke^{-\delta(t-\tau)}|\xi_{+}|+\int_{0}^{\infty}Ke^{-\delta|t-s|}ds\,\sup_{t\geq
0}|k(t)|$

from which (3.19) easily follows. ∎

##  4 Bounded solutions on the half lines

In this section we prove the existence of continuous solutions $(x(t),y(t))$,
with $t\geq 0$, of the perturbed linear system (1.7) such that

$\sup_{t\geq 0}|x(t)-u(t,y(t))|<\rho$

where $\rho>0$ is a sufficiently small positive real number. By a similar
argument we can also prove the existence of continuous solutions of (1.7),
with $t\leq 0$, such that

$\sup_{t\leq 0}|x(t)-u(t,y(t))|<\rho.$

From $A_{3})$ it follows that the number of the eigenvalues of
$f_{N,x}(w_{+}(y),y)$, $f_{-M,x}(w_{-}(y),y)$ with negative (and then also
positive) real parts, counted with multiplicities, is independent of
$y\in\mathbb{R}^{m}$. Moreover it also follows that all eigenvalues are
bounded functions of $y\in\mathbb{R}^{m}$. Indeed, since $f_{N,x}(w_{+}(y),y)$
is bounded, the matrix $\mathbb{I}-\lambda^{-1}f_{N,x}(w_{+}(y),y)$ is
invertible for $|\lambda|>R$, sufficiently large and independent of $y$. Hence
all eigenvalues have to satisfy $|\lambda|\leq R$. The same arguments work as
far as the eigenvalues of $f_{-M,x}(w_{-}(y),y)$ are concerned.

As in $A_{3})$, let $k<n$ be the number of eigenvalues with negative real
parts, counted with multiplicities, of the matrix $f_{N,x}(w_{+}(y),y)$ and
$\delta_{0}$ be any positive number strictly less than
$\min\\{|{Re}\lambda(y)|\\}$, where $\lambda(y)$ are the eigenvalues of
$f_{N,x}(w_{+}(y),y)$. According to [11] the system
$\dot{x}=f_{N,x}(w_{+}(y),y)x$ has an exponential dichotomy on $\mathbb{R}$
with exponent $\delta_{0}$ and (spectral) projection

$\begin{array}[]{l}\displaystyle P^{0}_{+}(y)=\frac{1}{2\pi
i}\int_{\Gamma}(z\mathbb{I}-f_{N,x}(w_{+}(y),y))^{-1}dz\\\
\qquad\displaystyle=\sum_{Re\lambda(y)<0}{\rm
Res}((z\mathbb{I}-f_{N,x}(w_{+}(y),y))^{-1},z=\lambda(y))\end{array}$

where ${\rm Res}(F(z),z=z_{0})$ is the residual of the meromorphic function
$F(z)$ at $z_{0}$ and $\Gamma$ is a closed curve that contains in its interior
all eigenvalues of $f_{N,x}(0,y)$ with negative real parts, but none of those
with positive real parts. Hence $\sup_{y\in\mathbb{R}^{m}}|P_{+}^{0}(y)|\leq
M$, for some $M\geq 1$. Similarly we see that $\dot{x}=f_{-M,x}(w_{-}(y),y)x$
has an exponential dichotomy on $\mathbb{R}$ with exponent $\delta_{0}$ and
projection $P^{0}_{-}(y)$ such that
$\sup_{y\in\mathbb{R}^{m}}|P_{+}^{0}(y)|\leq M$, for some $M\geq 1$.

Now, from $A_{2})$ we know that

$\lim_{t\to\pm\infty}u^{\pm}_{N}(t,y)=w_{\pm}(y)$

uniformly with respect to $y\in\mathbb{R}^{m}$.

Let $T_{+}>\sup_{y\in\mathbb{R}^{m}}t_{N}(y)$,
$T_{-}<\inf_{y\in\mathbb{R}^{m}}t_{-M}(y)$ and take $0<\delta<\delta_{0}$.
From the roughness of exponential dichotomies (cfr. [11, Proposition 2, p.
34]) the linear systems

$\dot{x}=f_{N,x}(u^{+}_{N}(t+T_{+},y),y)x$ (4.1)

and

$\dot{x}=f_{-M,x}(u^{-}_{-M}(t+T_{-},y),y)x$ (4.2)

have an exponential dichotomy on $\mathbb{R}_{+}$, $\mathbb{R}_{-}$ resp.,
uniform with respect to $y\in\mathbb{R}^{m}$, with projections $P_{+}(y)$,
resp. $P_{-}(y)$, of rank $k$, constant $K$ and exponent $\delta$. Moreover,
according to [20, Proposition 2.3], it can be assumed that, for $|y-y_{0}|$
sufficiently small it results: ${\cal N}P_{+}(y)={\cal N}P_{+}(y_{0})$, ${\cal
R}P_{-}(y)={\cal R}P_{-}(y_{0})$ and in this case the projections are smooth
with respect to $y$. Note that, ${\cal N}P_{+}(y)={\cal N}P_{+}(y_{0})$ and
${\cal R}P_{-}(y)={\cal R}P_{-}(y_{0})$ are equivalent to

$\begin{array}[]{l}P_{+}(y)=P_{+}(y)P_{+}(y_{0}),\quad
P_{+}(y_{0})=P_{+}(y_{0})P_{+}(y)\\\ P_{-}(y)=P_{-}(y_{0})P_{-}(y),\quad
P_{-}(y_{0})=P_{-}(y)P_{-}(y_{0}).\end{array}$ (4.3)

For $i=0,1,\ldots,N$, $j=0,-1,\ldots,-M$, let $U_{i}^{+}(t,y)$, [resp.
$U_{j}^{-}(t,y)$] be the fundamental matrix of

$\dot{x}=f_{i,x}(u_{i}^{+}(t,y),y)x\quad\hbox{[resp.
}\dot{x}=f_{j,x}(u_{j}^{-}(t,y),y)x\hbox{]}$

in $\mathbb{R}_{+}$ resp., $\mathbb{R}_{-}$ that is

$\begin{array}[]{l}\dot{U}_{i}^{+}(t,y)=f_{i,x}(u_{i}^{+}(t,y),y)U_{i}^{+}(t,y),\quad
t\geq 0,\\\ U_{i}^{+}(0,y)=\mathbb{I}.\end{array}$

resp.

$\begin{array}[]{l}\dot{U}_{j}^{-}(t,y)=f_{j,x}(u_{j}^{-}(t,y),y)U_{j}^{-}(t,y),\quad
t\leq 0,\\\ U_{j}^{-}(0,y)=\mathbb{I}.\end{array}$

We have the following

######  Lemma 4.1.

For any $\tau\in\mathbb{R}$ the linear system

$\dot{x}=f_{N,x}(u^{+}_{N}(t,y),y)x,$ (4.4)

resp.

$\dot{x}=f_{-M,x}(u^{-}_{-M}(t,y),y)x,$ (4.5)

has an exponential dichotomy on $t\geq\tau$, resp. $t\leq\tau$, with exponent
$\delta$, constant $\tilde{K}$ independent on $y$ and projections

$\begin{array}[]{l}Q_{+}(y)=U^{+}_{N}(\tau,y)U^{+}_{N}(T_{+},y)^{-1}P_{+}(y)U^{+}_{N}(T_{+},y)U^{+}_{N}(\tau,y)^{-1}\\\
Q_{-}(y)=U^{-}_{-M}(\tau,y)U^{-}_{-M}(T_{-},y)^{-1}P_{-}(y)U^{-}_{-M}(T_{-},y)U^{-}_{-M}(\tau,y)^{-1}.\end{array}$

In particular, if $\tau=T_{+}$, resp. $\tau=T_{-}$, then $Q_{+}(y)=P_{+}(y)$,
resp. $Q_{-}(y)=P_{-}(y)$, and $\tilde{K}=K$.

###### Proof.

As $U_{N}^{+}(t,y)$ is the fundamental matrix of (4.4) in $\mathbb{R}_{\pm}$,
the fundamental matrix of (4.1) on $\mathbb{R}_{+}$ is

$\tilde{U}(t,y):=U^{+}_{N}(t+T_{+},y)U^{+}_{N}(T_{+},y)^{-1}$

while the fundamental matrix of (4.4) on $t\geq\tau$ is

$\bar{U}(t,y)=U^{+}_{N}(t,y)U^{+}_{N}(\tau,y)^{-1}.$

Then

$\bar{U}(t,y)U^{+}_{N}(\tau,y)=U^{+}_{N}(t,y)=\tilde{U}(t-T_{+},y)U^{+}_{N}(T_{+},y).$

So we get, for $t\geq s\geq T_{+}$:

$\begin{array}[]{l}|\bar{U}(t,y)Q_{+}(y)\bar{U}(s,y)^{-1}|=\\\
|\tilde{U}(t-T_{+},y)U^{+}_{N}(T_{+},y)U^{+}_{N}(\tau,y)^{-1}Q_{+}(y)U^{+}_{N}(\tau,y)U^{+}_{N}(T_{+},y)^{-1}\tilde{U}(s-T_{+},y)^{-1}|=\\\
|\tilde{U}(t-T_{+},y)P_{+}(y)\tilde{U}(s-T_{+},y)^{-1}|\leq
Ke^{-\delta(t-s)}\end{array}$

and similarly

$\begin{array}[]{l}|\bar{U}(s,y)(\mathbb{I}-Q_{+}(y))\bar{U}(t,y)^{-1}|=\\\
|\tilde{U}(s-T_{+},y)(\mathbb{I}-P_{+}(y))\tilde{U}(t-T_{+},y)^{-1}|\leq
Ke^{-\delta(t-s)}.\end{array}$

Hence (4.4) has an exponential dichotomy on $t\geq T_{+}$ with projection
$Q_{+}(y)$, constant $\tilde{K}=K$ and exponent $\delta$. If $\tau\geq T_{+}$
the proof is finished. Moreover if $\tau=T_{+}$ we have

$Q_{+}(y)=P_{+}(y)\quad\hbox{and}\quad Q_{-}(y)=P_{-}(y).$

If $\tau<T_{+}$, for $\tau\leq s\leq t\leq T_{+}$ let

$v(t):=\bar{U}(t,y)Q_{+}(y)\bar{U}(s,y)^{-1}=\tilde{U}(t-T_{+},y)P_{+}(y)\tilde{U}(s-T_{+},y)^{-1}.$

Then $v(t)$ satisfies

$\dot{v}(t)=f_{N,x}(u^{+}_{N}(t,y),y)v(t)$

with $v(T_{+})=P_{+}(y)\tilde{U}(s-T_{+},y)^{-1}$. Hence

$v(t)=P_{+}(y))\tilde{U}(s-T_{+},y)^{-1}-\int_{t}^{T_{+}}f_{N,x}(u^{+}_{N}(s,y),y)v(s)ds.$

So, using $|P_{+}(y)\tilde{U}(s-T_{+},y)^{-1}|\leq Ke^{-\delta(T_{+}-s)}$:

$|v(t)|\leq Ke^{-\delta(T_{+}-s)}+\int_{t}^{T_{+}}\bar{F}|v(s)|ds.$

where $\bar{F}$ is an upper bound for $|f_{N,x}(u^{+}_{N}(s,y),y)|$. Setting
$r(t):\int_{t}^{T_{+}}|v(s)|ds$ we get

$\frac{d}{dt}[r(t)e^{\bar{F}t}]=[r^{\prime}(t)+\bar{F}r(t)]e^{\bar{F}t}\geq-
Ke^{\bar{F}t}.$

Integrating on $[t,T_{+}]$ we get, since $r(T_{+})=0$,
$-r(t)e^{\bar{F}t}\geq-K\bar{F}^{-1}[e^{\bar{F}T_{+}}-e^{\bar{F}t}]$. Hence

$\int_{t}^{T_{+}}\bar{F}|v(s)|ds=\bar{F}r(t)\leq K[e^{\bar{F}(T_{+}-t)}-1]$

and then

$\begin{array}[]{l}|v(t)|\leq Ke^{-\delta(T_{+}-s)}+Ke^{\bar{F}(T_{+}-t)}\\\
=K[e^{-\delta(T_{+}-t)}+e^{\bar{F}(T_{+}-t)+\delta(t-s)}]e^{-\delta(t-s)}\leq\hat{K}e^{-\delta(t-s)}\end{array}$

where

$\hat{K}=K[1+e^{(\bar{F}+\delta)(T_{+}-\tau)}].$

Note that $\hat{K}$ is independent of $y$. Finally, for $\tau\leq s\leq
T_{+}\leq t$ we have

$\begin{array}[]{l}|\bar{U}(t,y)Q_{+}(y)\bar{U}(s,y)^{-1}|\leq|\bar{U}(t,y)Q_{+}(y)\bar{U}(T_{+},y)^{-1}||\bar{U}(T_{+},y)Q_{+}(y)\bar{U}(s,y)^{-1}|\\\
\leq
Ke^{-\delta(t-T_{+})}\hat{K}e^{\delta(s-T_{+})}=K\hat{K}e^{\delta(t-s)}.\end{array}$

By a similar argument we prove that

$|\tilde{U}(t-T_{+},y)(\mathbb{I}-P_{+}(y))\tilde{U}(s-T_{+},y)^{-1}|\leq\tilde{K}e^{\delta(t-s)}$

for $\tau\leq t\leq s$.

The proof that (4.5) has an exponential dichotomy on $t\leq\tau$, with
projection $Q_{-}(y)$, exponent $\delta$ and a constant independent of
$y\in\mathbb{R}^{m}$ is similar. ∎

From the proof of the parametric stable (resp. unstable) Theorem (see [19, p.
18]) it follows that, all solutions $w(t,y)$ of $\dot{w}=f_{N}(w,y)$ such that
$\sup_{t\geq 0}|w(t,y)-w^{+}(y)|\leq\Delta\ll 1$ (resp. $\sup_{t\leq
0}|w(t,y)-w^{-}(y)|\leq\Delta\ll 1$) can be obtained as fixed points of a
uniform contraction, $C^{r}$ with respect to $y$, on the space of bounded
functions. It follows, then, that these solutions are bounded together with
their derivatives.

The following Lemma states that the bounds for $u^{+}_{N}(t,y)$,
$u^{-}_{-M}(t,y)$ and their derivatives with respect to $y$ can be taken
independent of $y$. Due to its technical character the proof of Lemma 4.2 is
postponed in the Appendix.

######  Lemma 4.2.

Assume $A_{1})-A_{3})$. Then $u^{+}_{N}(t,y)$ and its derivatives with respect
to $y$ are bounded uniformly with respect to $y$, on $t\geq t_{N}(y)$.
Similarly, $u^{-}_{-M}(t,y)$ and its derivatives with respect to $y$ are
bounded, uniformly with respect to $y$, on $t\leq t_{-M}(y)$.

From [6, Theorem 2] we get the following

######  Theorem 4.3.

Suppose that $A_{1})-A_{4})$ hold and let $\gamma,\beta>0$ be such that
$r\beta<\gamma<\delta$. Then there exist $\rho>0$, $\bar{\varepsilon}>0$ and
$\bar{\sigma}>0$ such that for
$(\alpha,\varepsilon)\in\mathbb{R}^{m}\times\mathbb{R}$ and $\xi_{+}\in{\cal
R}P_{+}(\alpha)$, with $|\alpha-y_{0}|<\bar{\sigma}$,
$|\xi_{+}|<\bar{\sigma}$, and $|\varepsilon|<\bar{\varepsilon}$, system

$\begin{array}[]{l}\dot{x}=f_{N}(x,y)\\\ \dot{y}=\varepsilon
g(x,y,\varepsilon)\\\ P_{+}(\alpha)[x(T_{+})-u(T_{+},\alpha)]=\xi_{+}\\\
y(T_{+})=\alpha,\end{array}$ (4.6)

has a unique solution
$(x^{+}_{N}(t,\xi_{+},\alpha,\varepsilon),y^{+}_{N}(t,\xi_{+},\alpha,\varepsilon))$,
$t\geq T_{+}$ such that

$\sup_{t\geq
T_{+}}|x^{+}_{N}(t,\xi_{+},\alpha,\varepsilon)-u^{+}_{N}(t,y_{+}(t,\xi_{+},\alpha,\varepsilon))|<\rho.$
(4.7)

Moreover,

$\sup_{t\geq
T_{+}}|x^{+}_{N}(t,\xi_{+},\alpha,\varepsilon)-u^{+}_{N}(t,y_{+}(t,\xi_{+},\alpha,\varepsilon))|\to
0$ (4.8)

as $|\xi_{+}|+|\varepsilon|\to 0$, uniformly with respect to $\alpha$.
Similarly for $(\alpha,\varepsilon)\in\mathbb{R}^{m}\times\mathbb{R}$ and
$\xi_{-}\in{\cal N}P_{-}(\alpha)$ such that $|\alpha-y_{0}|<\bar{\sigma}$,
$|\xi_{-}|<\bar{\sigma}$ system

$\begin{array}[]{l}\dot{x}=f_{-M}(x,y)\\\ \dot{y}=\varepsilon
g(x,y,\varepsilon)\\\
(\mathbb{I}-P_{-}(\alpha))[x(T_{-})-u(T_{-},\alpha)]=\xi_{-}\\\
y(T_{-})=\alpha\end{array}$

has a unique solution
$(x^{-}_{-M}(t,\xi_{-},\alpha,\varepsilon),y^{-}_{-M}(t,\xi_{-},\alpha,\varepsilon))$,
$t\leq T_{-}$ such that

$\sup_{t\leq
T_{-}}|x^{-}_{-M}(t,\xi_{-},\alpha,\varepsilon)-u^{-}_{-M}(t,y_{-}(t,\xi_{-},\alpha,\varepsilon))|<\rho.$
(4.9)

Moreover

$\sup_{t\leq
T_{-}}|x^{-}_{-M}(t,\xi_{-},\alpha,\varepsilon)-u^{-}_{-M}(t,y_{-}(t,\xi_{-},\alpha,\varepsilon))|\to
0$ (4.10)

as $|\xi_{-}|+|\varepsilon|\to 0$, uniformly with respect to $\alpha$.
Moreover, for any $p\leq r-1$, the function

$(\xi_{+},\alpha,\varepsilon)\mapsto(x^{+}_{N}(t,\xi_{+},\alpha,\varepsilon),y^{+}_{N}(t,\xi_{+},\alpha,\varepsilon))$

is of class $C^{p}$ and, for $t\geq T_{+}$, the $p$-th order derivatives are
bounded above in absolute value by $K^{\prime}_{p}e^{(p+1)\beta|t-t_{N}(y)|}$,
where $K^{\prime}_{p}$ is a suitable constant. Similarly, for $t\leq T_{-}$,
the map

$(\xi_{-},\alpha,\varepsilon)\mapsto(x^{-}_{-M}(t,\xi_{-},\alpha,\varepsilon),y^{-}_{-M}t,\xi_{-},\alpha,\varepsilon))$

is of class $C^{p}$ and the $p$-th order derivatives are bounded above in
absolute value by $K^{\prime\prime}_{p}e^{(p+1)\beta|t-t_{-N}(y)|)}$, where
$K^{\prime\prime}_{p}$ is a suitable constant.

###### Proof.

In [6] the equations are considered for $t\geq 0$ or $t\leq 0$ and instead of
$P_{\pm}(\alpha)$ it is considered $P_{\pm}(y_{0})$. To obtain the result as
in Theorem 4.3 for $\xi_{+}\in{\cal R}P_{+}(\alpha)$, (resp. $\xi_{-}\in{\cal
N}P_{-}(\alpha)$) we take $P_{+}(y_{0})\xi_{+}\in{\cal R}P_{+}(y_{0})$ (resp.
$(\mathbb{I}-P_{-}(y_{0}))\xi_{-}\in{\cal N}P_{-}(y_{0})$) and apply [6,
Theorem 1] with $u^{+}_{N}(t+T_{+},y)$, resp. $u^{-}_{-M}(t+T_{-},y)$, instead
of $u(t,y)$. Suppose, to fix ideas, that $t\geq 0$. We know that (4.1) has an
exponential dichotomy on $\mathbb{R}_{+}$ with projection $P_{+}(y)$ and, from
Lemma 4.2 that $u^{+}_{N}(t+T_{+},y)$ and its derivatives with respect to
$y\in\mathbb{R}^{m}$ are bounded uniformly with respect to
$y\in\mathbb{R}^{m}$. From [6, Theorem 1] we obtain then a unique solution

$(\tilde{x}(t,\xi_{+},\alpha,\varepsilon),\tilde{y}(t,\xi_{+},\alpha,\varepsilon))$

of (1.7) such that

$\begin{array}[]{l}P_{+}(y_{0})[\tilde{x}(0)-u^{+}_{N}(T_{+},\alpha)]=P_{+}(y_{0})\xi_{+}\\\
\tilde{y}(0)=\alpha\end{array}$

and

$\begin{array}[]{l}\sup_{t\geq
0}|\tilde{x}(t,\xi_{+},\alpha,\varepsilon)-u^{+}_{N}(t+T_{+},\tilde{y}(t,\xi_{+},\alpha,\varepsilon))|<\rho,\\\
\sup_{t\geq
0}|\tilde{x}(t,\xi_{+},\alpha,\varepsilon)-u^{+}_{N}(t+T_{+},\tilde{y}(t,\xi_{+},\alpha,\varepsilon))|\to
0\end{array}$

as $|\xi_{+}|+|\varepsilon|\to 0$, uniformly with respect to $\alpha$. Setting

$\begin{array}[]{l}x^{+}_{N}(t,\xi_{+},\alpha,\varepsilon)=\tilde{x}(t-T_{+},\xi_{+},\alpha,\varepsilon)\\\
y^{+}_{N}(t,\xi_{+},\alpha,\varepsilon)=\tilde{y}(t-T_{+},\xi_{+},\alpha,\varepsilon)\end{array}$

we see that
$(x^{+}_{N}(t,\xi_{+},\alpha,\varepsilon),y^{+}_{N}(t,\xi_{+},\alpha,\varepsilon))$
satisfies

$\begin{array}[]{l}\dot{x}=f_{N}(x,y)\\\ \dot{y}=\varepsilon
g(x,y,\varepsilon)\\\ y(T_{+})=\alpha\end{array}$

and

$P_{+}(\alpha)[x(T_{+})-u^{+}_{N}(T_{+},\alpha)]=P_{+}(\alpha)P_{+}(y_{0})\xi_{+}=\xi_{+}$

since $P_{+}(\alpha)P_{+}(y_{0})=P_{+}(\alpha)$. A similar argument works for
$t\leq T_{-}$. Finally, we observe that, although in [6] it is not explicitly
stated that (4.7), (4.9) hold uniformly with respect to $\alpha$, this fact
easily follow from [6, (20)-(23)] and
$\sup_{t\in\mathbb{R}_{\pm}}|\dot{y}_{\pm}(t,\xi_{\pm},\alpha,\varepsilon)|=O(\varepsilon)$.
∎

######  Remark 4.4.

According to assumption $A_{3})$, $x=w_{\pm}(y)$ are normally hyperbolic
manifolds for the system

$\dot{x}=f_{N}(x,y),\quad\dot{y}=0.$

These manifolds perturb to normally hyperbolic invariant manifolds
$x=\tilde{w}_{\pm}(y,\varepsilon)$ for the system

$\begin{array}[]{l}\dot{x}=f_{N}(x,y),\\\ \dot{y}=\varepsilon
g(x,y,\varepsilon)\end{array}$

(see, for example, [1, 22]). Let $y_{c}(t,\alpha,\varepsilon)$ be the solution
of $\dot{y}=g(\tilde{w}_{\pm}(y,\varepsilon),y,\varepsilon)$ such that
$y_{c}(0,\alpha,\varepsilon)=\alpha$.

It follows from [9, Theorem 1] that, for $|\alpha-y_{0}|$ and $|\xi^{0}_{+}|$
sufficiently small, with $\alpha\in\mathbb{R}^{m}$ and $\xi^{0}_{+}\in{\cal
R}P_{+}^{0}(y_{0})$, there exists a unique solution of

$\begin{array}[]{l}\dot{x}=f_{N}(x,y)\\\ \dot{y}=\varepsilon
g(x,y,\varepsilon)\end{array}$

such that

$\begin{array}[]{l}\sup_{t\geq T_{+}}e^{\beta
t}|x(t,\xi^{0}_{+},\alpha,\varepsilon)-\tilde{w}_{+}(y(t,\xi^{0}_{+},\alpha,\varepsilon),\varepsilon)|\leq\mu_{1}\\\
\sup_{t\geq T_{+}}e^{\beta
t}|y(t,\xi^{0}_{+},\alpha,\varepsilon)-y_{c}(\varepsilon(t-T_{+}),\alpha,\varepsilon)|\leq\mu_{2}\end{array}$

and
$P_{+}^{0}(y_{0})[x(T_{+},\xi^{0}_{+},\alpha,\varepsilon)-\tilde{w}_{+}(y(0,\xi^{0}_{+},\alpha,\varepsilon),\varepsilon)]=\xi^{0}_{+}$.

As $\lim_{\varepsilon\to 0}|\tilde{w}_{+}(y,\varepsilon)-w_{+}(y)|=0$ and
$\lim_{t\to\infty}|u_{N}(t,y)-w_{+}(y)|=0$ uniformly with respect to $y$ we
conclude that

$\sup_{t\geq
T}|x(t,\xi^{0}_{+},\alpha,\varepsilon)-u_{N}(t,y(t,\xi^{0}_{+},\alpha,\varepsilon))|<\rho$

where $T\geq T_{+}$ is sufficiently large. From continuous dependence we see
that the above estimate holds with $\sup_{t\geq T_{+}}$ instead of
$\sup_{t\geq T}$, provided $|\xi^{0}_{+}|+|\varepsilon|\ll 1$ is sufficiently
small, and then

$x^{+}_{N}(t,\xi_{+},\alpha,\varepsilon)=x(t,\xi^{0}_{+},\alpha_{0},\varepsilon),\quad
y^{+}_{N}(t,\xi_{+},\alpha,\varepsilon)=y(t,\xi^{0}_{+},\alpha_{0},\varepsilon),$

for some $\xi^{0}_{+}$ and $\alpha_{0}$. As a consequence the solutions given
in Theorem 4.3 satisfy

$\begin{array}[]{l}\sup_{t\geq T_{+}}e^{\beta
t}|x^{+}_{N}(t,\xi_{+},\alpha,\varepsilon)-\tilde{w}_{+}(y_{N}(t,\xi^{0}_{+},\alpha_{0},\varepsilon),\varepsilon)|\leq\mu_{1}\\\
\sup_{t\geq T_{+}}e^{\beta
t}|y^{+}_{N}(t,\xi_{+},\alpha,\varepsilon)-y_{c}(\varepsilon(t-T_{+}),\xi^{0}_{+},\varepsilon)|\leq\mu_{2}.\end{array}$

Similarly we see that

$\begin{array}[]{l}\sup_{t\leq T_{-}}e^{-\beta
t}|x^{-}_{-M}(t,\xi_{-},\alpha,\varepsilon)-\tilde{w}_{-}(y_{-M}(t,\xi^{0}_{-},\alpha_{0},\varepsilon),\varepsilon)|\leq\mu_{1}\\\
\sup_{t\leq T_{-}}e^{-\beta
t}|y^{-}_{-M}(t,\xi_{-},\alpha,\varepsilon)-y_{c}(\varepsilon(t+T_{-}),\xi^{0}_{-},\varepsilon)|\leq\mu_{2}.\end{array}$

In the remaining part of this section we extend the solutions obtained in
Theorem 4.3 to continuous, piecewise $C^{1}$ solutions of (1.7) in
$\mathbb{R}_{+}$ resp. $\mathbb{R}_{-}$. We have the following

######  Theorem 4.5.

There exist $\rho>0$, bounded $C^{r}$-functions

$t^{*}_{-M}(\xi_{-},\alpha,\varepsilon)<\ldots<t^{*}_{-1}(\xi_{-},\alpha,\varepsilon)<0<t^{*}_{1}(\xi_{+},\alpha,\varepsilon)<\ldots<t^{*}_{N}(\xi_{+},\alpha,\varepsilon),$

and continuous, piecewise $C^{r}$ solutions of (1.7)

$(x_{\pm}(t,\xi_{\pm},\alpha,\varepsilon),y_{\pm}(t,\xi_{\pm},\alpha,\varepsilon))$

defined for $t\geq 0$ and $t\leq 0$ resp., such that

$\begin{array}[]{l}\lim_{(\xi_{+},\varepsilon)\to
0}|t^{*}_{i}(\xi_{+},\alpha,\varepsilon)-t_{i}(\alpha)|=0,\\\
\lim_{(\xi_{-},\varepsilon)\to
0}|t^{*}_{j}(\xi_{-},\alpha,\varepsilon)-t_{j}(\alpha)|=0\end{array}$

($i=1,\ldots,N$, $j=-1,\ldots,-M$) uniformly with respect to
$\alpha\in\mathbb{R}^{m}$ and

$\begin{array}[]{l}c_{i-2}<h(x_{+}(t,\xi_{+},\alpha,\varepsilon),y_{+}(t,\xi_{+},\alpha,\varepsilon))<c_{i-1},\quad\hbox{
for
$t^{*}_{i-1}(\xi_{+},\alpha,\varepsilon)<t<t^{*}_{i}(\xi_{+},\alpha,\varepsilon)$}\\\
h(x_{+}(t,\xi_{+},\alpha,\varepsilon),y_{+}(t,\xi_{+},\alpha,\varepsilon))>c_{N-1},\quad\hbox{
for $t>t^{*}_{N}(\xi_{+},\alpha,\varepsilon)$}\\\
c_{j}<h(x_{-}(t,\xi_{-},\alpha,\varepsilon),y_{-}(t,\xi_{-},\alpha,\varepsilon))<c_{j+1},\quad\hbox{
for
$t^{*}_{j}(\xi_{-},\alpha,\varepsilon)<t<t^{*}_{j+1}(\xi_{-},\alpha,\varepsilon)$}\\\
h(x_{-}(t,\xi_{-},\alpha,\varepsilon),y_{-}(t,\xi_{-},\alpha,\varepsilon))<c_{-M},\quad\hbox{
for $t<t^{*}_{-M}(\xi_{-},\alpha,\varepsilon)$}\\\
h(x_{+}(t^{*}_{i}(\xi_{+},\alpha,\varepsilon),\xi_{+},\alpha,\varepsilon),y_{+}(t^{*}_{i}(\xi_{+},\alpha,\varepsilon),\xi_{+},\alpha,\varepsilon))=c_{i-1}\\\
h(x_{-}(t^{*}_{j}(\xi_{-},\alpha,\varepsilon),\xi_{-},\alpha,\varepsilon),y_{-}(t^{*}_{j}(\xi_{-},\alpha,\varepsilon),\xi_{-},\alpha,\varepsilon))=c_{j}\\\
\frac{\partial}{\partial
t}h(x_{+}(t,\xi_{+},\alpha,\varepsilon),y_{+}(t,\xi_{+},\alpha,\varepsilon))_{|t=t^{*}_{i}(\xi_{+},\alpha,\varepsilon)}>\eta\\\
\frac{\partial}{\partial
t}h(x_{-}(t,\xi_{-},\alpha,\varepsilon),y_{-}(t,\xi_{-},\alpha,\varepsilon))_{|t=t^{*}_{j}(\xi_{-},\alpha,\varepsilon)}>\eta\end{array}$

where $t^{*}_{0}(\xi_{\pm},\alpha,\varepsilon)=0$. Moreover

$\begin{array}[]{ll}(x_{+}(t,\xi_{+},\alpha,\varepsilon),y_{+}(t,\xi_{+},\alpha,\varepsilon))=(x^{+}_{N}(t,\xi_{+},\alpha,\varepsilon),y^{+}_{N}(t,\xi_{+},\alpha,\varepsilon))&\hbox{for
$t\geq T_{+}$}\\\
(x_{-}(t,\xi_{+},\alpha,\varepsilon),y_{-}(t,\xi_{+},\alpha,\varepsilon))=(x^{-}_{-M}(t,\xi_{+},\alpha,\varepsilon),y^{-}_{-M}(t,\xi_{+},\alpha,\varepsilon))&\hbox{for
$t\leq T_{-}$}\end{array}$

and

$\begin{array}[]{ll}\rho>\sup_{t\geq
0}|x_{+}(t,\xi_{+},\alpha,\varepsilon)-u(t,y_{+}(t,\xi_{+},\alpha,\varepsilon))|\to
0&\hbox{as $|\xi_{+}|+|\varepsilon|\to 0$}\\\ \rho>\sup_{t\leq
0}|x_{-}(t,\xi_{-},\alpha,\varepsilon)-u(t,y_{-}(t,\xi_{-},\alpha,\varepsilon))|\to
0&\hbox{as $|\xi_{-}|+|\varepsilon|\to 0$}\end{array}$ (4.11)

uniformly with respect to $\alpha$.

###### Proof.

Suppose $t\geq 0$. Since $y^{+}_{N}(T_{+},\xi_{+},\alpha,\varepsilon)=\alpha$,
from (4.7) we see that

$|x^{+}_{N}(T_{+},\xi_{+},\alpha,\varepsilon)-u^{+}_{N}(T_{+},\alpha)|\to 0$

as $|\xi_{+}|+|\varepsilon|\to 0$ (uniformly with respect to $\alpha$). Hence
$h(x^{+}_{N}(T_{+},\xi_{+},\alpha,\varepsilon),\alpha)>c_{N}$, provided
$|\xi_{+}|+|\varepsilon|$ is sufficiently small and uniformly with respect to
$\alpha$. Then
$(x^{+}_{N}(t,\xi_{+},\alpha,\varepsilon),y^{+}_{N}(t,\xi_{+},\alpha,\varepsilon)$
can be extended to a solution of

$\left\\{\begin{array}[]{l}\dot{x}=f_{N}(x,y)\\\ \dot{y}=\varepsilon
g(x,y)\end{array}\right.$

which is defined for $t\geq 0$ and such that

$\sup_{0\leq t\leq
T_{+}}|x^{+}_{N}(t,\xi_{+},\alpha,\varepsilon)-u^{+}_{N}(t,y^{+}_{N}(t,\xi_{+},\alpha,\varepsilon))|\to
0$

as $|\xi_{+}|+|\varepsilon|\to 0$, uniformly with respect to $\alpha$. Note
that

$\sup_{0\leq t\leq T_{+}}|y^{+}_{N}(t,\xi_{+},\alpha,\varepsilon)-\alpha|\to
0$

as $\varepsilon\to 0$, uniformly with respect to $(\xi_{+},\alpha)$, since
$y^{+}_{N}(T_{+},\xi_{+},\alpha,\varepsilon)=\alpha$ and
$\dot{y}^{+}_{N}(t,\xi_{+},\alpha,\varepsilon)=O(|\varepsilon|)$ uniformly
with respect to $(\xi_{+},\alpha)$. Now, from $A_{1})$ and the implicit
function theorem it follows that there exists a $C^{r}$-function
$t^{*}_{N}(\xi_{+},\alpha,\varepsilon)$, bounded together with its
derivatives, such that

$\begin{array}[]{l}|t^{*}_{N}(\xi_{+},\alpha,\varepsilon)-t_{N}(\alpha)|\to
0,\quad\hbox{as $|\xi_{+}|+|\varepsilon|\to 0$}\\\
h(x^{+}_{N}(t,\xi_{+},\alpha,\varepsilon),y^{+}_{N}(t,\xi_{+},\alpha,\varepsilon))>c_{N},\quad\hbox{for
$t>t^{*}_{N}(\xi_{+},\alpha,\varepsilon)$}\\\
h(x^{+}_{N}(t^{*}_{N}(\xi_{+},\alpha,\varepsilon),\xi_{+},\alpha,\varepsilon),y^{+}_{N}(t^{*}_{N}(\xi_{+},\alpha,\varepsilon),\xi_{+},\alpha,\varepsilon))=c_{N}\\\
\frac{\partial}{\partial
t}h(x^{+}_{N}(t,\xi_{+},\alpha,\varepsilon),y^{+}_{N}(t,\xi_{+},\alpha,\varepsilon))_{|t=t^{*}_{N}(\xi_{+},\alpha,\varepsilon)}>\eta\end{array}$

provided $|\xi_{+}|+|\varepsilon|$ is sufficiently small, uniformly with
respect to $\alpha$. Next, from the continuous dependence on the data, we see
that the system

$\begin{array}[]{l}\dot{x}=f_{N-1}(x,y)\\\ \dot{y}=\varepsilon
g(x,y,\varepsilon)\\\
x(t^{*}_{N}(\xi_{+},\alpha,\varepsilon))=x^{+}_{N}(t^{*}_{N}(\xi_{+},\alpha,\varepsilon),\xi_{+},\alpha,\varepsilon)\\\
y(t^{*}_{N}(\xi_{+},\alpha,\varepsilon))=y^{+}_{N}(t^{*}_{N}(\xi_{+},\alpha,\varepsilon),\xi_{+},\alpha,\varepsilon)\end{array}$

has a unique solution
$(x^{+}_{N-1}(t,\xi_{+},\alpha,\varepsilon),y^{+}_{N-1}(t,\xi_{+},\alpha,\varepsilon))$,
defined for $-1\leq t\leq t^{*}_{N}(\xi_{+},\alpha,\varepsilon)$, such that

$\sup_{-1\leq t\leq
t^{*}_{N}}|x^{+}_{N-1}(t,\xi_{+},\alpha,\varepsilon)-u^{+}_{N-1}(t,y^{+}_{N-1}(t,\xi_{+},\alpha,\varepsilon))|\leq\rho$

(where $t^{*}_{N}=t^{*}_{N}(\xi_{+},\alpha,\varepsilon)$) for
$|\xi_{+}|+|\varepsilon|$ sufficiently small and the following holds

$\begin{array}[]{l}\sup_{-1\leq t\leq
t^{*}_{N}}|x^{+}_{N-1}(t,\xi_{+},\alpha,\varepsilon)-u^{+}_{N-1}(t,y^{+}_{N-1}(t,\xi_{+},\alpha,\varepsilon))|\to
0\\\ \sup_{-1\leq t\leq
t^{*}_{N}}|y^{+}_{N-1}(t,\xi_{+},\alpha,\varepsilon)-\alpha|\to 0\\\
\end{array}$

as $|\xi_{+}|+|\varepsilon|\to 0$, uniformly with respect to $\alpha$. Now, as

$\begin{array}[]{l}y^{+}_{N-1}(t,\xi_{+},\alpha,0)=y^{+}_{N}(t^{*}_{N}(\xi_{+},\alpha,0),\xi_{+},\alpha,0)=\alpha\\\
c_{N-2}<h(u^{+}_{N-1}(t,\alpha),\alpha)<c_{N-1}\hbox{ for
$t_{N-1}(\alpha)<t<t_{N}(\alpha)$,}\\\ \frac{\partial}{\partial
t}h(u^{+}_{N}(t,\alpha),\alpha)_{|t=t_{N-1}(\alpha)}>2\eta\end{array}$

from the implicit function theorem we see that a $C^{r}$-function
$t^{*}_{N-1}(\xi_{+},\alpha,\varepsilon)$, bounded together with its
derivatives, exists such that

$\lim_{(\xi_{+},\varepsilon)\to(0,0)}t^{*}_{N-1}(\xi_{+},\alpha,\varepsilon)-t_{N-1}(\alpha)=0$

uniformly with respect to $\alpha$ and the following holds:

$\begin{array}[]{l}c_{N-2}<h(x^{+}_{N-1}(t,\xi_{+},\alpha,\varepsilon),y^{+}_{N-1}(t,\xi_{+},\alpha,\varepsilon))<c_{N-1},\\\
\qquad\hbox{for $t^{*}_{N-1}(\xi_{+},\alpha,\varepsilon)\leq
t<t^{*}_{N}(\xi_{+},\alpha,\varepsilon)$}\\\
h(x^{+}_{N-1}(t^{*}_{N-1}(\xi_{+},\alpha,\varepsilon),\xi_{+},\alpha,\varepsilon),y^{+}_{N-1}(t^{*}_{N-1}(\xi_{+},\alpha,\varepsilon),\xi_{+},\alpha,\varepsilon))=c_{N-2}\\\
\frac{\partial|}{\partial
t}h(x^{+}_{N-1}(t,\xi_{+},\alpha,\varepsilon),y^{+}_{N-1}(t,\xi_{+},\alpha,\varepsilon))_{t=t^{*}_{N-1}(\xi_{+},\alpha,\varepsilon)}>\eta.\end{array}$

Proceeding this way we construct the solution
$(x_{+}(t,\xi_{+},\alpha,\varepsilon),y_{+}(t,\xi_{+},\alpha,\varepsilon)$
with the properties stated in the Theorem. A similar argument works for $t\leq
0$. The proof is complete ∎

######  Remark 4.6.

According to Theorem 4.5 we have

$\begin{array}[]{l}h(x_{+}(t^{*}_{i}(\xi_{+},\alpha,\varepsilon),\xi_{+},\alpha,\varepsilon),y_{+}(t^{*}_{i}(\xi_{+},\alpha,\varepsilon),\xi_{+},\alpha,\varepsilon))=c_{i-1}\\\
h(x_{-}(t^{*}_{j}(\xi_{-},\alpha,\varepsilon),\xi_{-},\alpha,\varepsilon),y^{*}_{j}(t_{*}(\xi_{-},\alpha,\varepsilon),\xi_{-},\alpha,\varepsilon))=c_{j}.\end{array}$
(4.12)

Differentiating the above equalities with respect to $\xi_{+}$, $\xi_{-}$, at
$\varepsilon=0$ we obtain a formula for the derivatives

$\frac{\partial
t^{*}_{i}}{\partial\xi_{+}}(\xi_{+},\alpha,0),\quad\frac{\partial
t^{*}_{j}}{\partial\xi_{-}}(\xi_{-},\alpha,0).$

However we have to distinguish when $t\to t^{*}_{i}(\xi_{+},\alpha,0)^{+}$ or
$t\to t^{*}_{i}(\xi_{+},\alpha,0)^{-}$ (resp. $t\to
t^{*}_{j}(\xi_{-},\alpha,0)^{+}$ or $t\to t^{*}_{j}(\xi_{-},\alpha,0)^{-}$).
For example if $t\to t^{*}_{i}(\xi_{+},\alpha,0)^{+}$,
$x_{+}(t,\xi_{+},\alpha,0)$ is the solution of $\dot{x}=f_{i}(x,\alpha)$ and
then, differentiating (4.12) with respect to $\xi_{+}$, we get, with
$t^{*}_{i}=t^{*}_{i}(\xi_{+},\alpha,0)$:

$h_{x}(x_{+}(t^{*},\xi_{+},\alpha,0),\alpha)[f_{i}(x_{+}(t^{*}_{i},\xi_{+},\alpha,0),\alpha)\frac{\partial
t^{*}_{i}}{\partial\xi_{+}}(\xi_{+},\alpha,0)+x_{+,\xi_{+}}(t^{*+}_{i},\xi_{+},\alpha,0)]=0.$

Vice versa, when $t\to t^{*}_{i}(\xi_{+},\alpha,0)^{-}$,
$x_{+}(t,\xi_{+},\alpha,0)$ is the solution of $\dot{x}=f_{i-1}(x,\alpha)$ and
then

$h_{x}(x_{+}(t^{*}_{i},\xi_{+},\alpha,0),\alpha)[f_{i-1}(x_{+}(t^{*}_{i},\xi_{+},\alpha,0),\alpha)\frac{\partial
t^{*}_{i}}{\partial\xi_{+}}(\xi_{+},\alpha,0)+x_{+,\xi_{+}}(t^{*-}_{i},\xi_{+},\alpha,0)]=0.$

Similarly we get

$h_{x}(x_{-}(t^{*}_{j},\xi_{-},\alpha,0),\alpha)[f_{j}(x_{-}(t^{*}_{j},\xi_{+},\alpha,0),\alpha)\frac{\partial
t^{*}_{j}}{\partial\xi_{-}}(\xi_{-},\alpha,0)+x_{-,\xi_{-}}(t^{*-}_{j},\xi_{-},\alpha,0)]=0$

and

$h_{x}(x_{-}(t^{*}_{j},\xi_{-},\alpha,0),\alpha)[f_{j+1}(x_{+}(t^{*}_{j},\xi_{-},\alpha,0),\alpha)\frac{\partial
t^{*}_{j}}{\partial\xi_{-}}(\xi_{-},\alpha,0)+x_{-,\xi_{-}}(t^{*+}_{j},\xi_{-},\alpha,0)]=0.$

We will use this remark in the next section.

##  5 The variational equation

Let $\ell\in\\{-M,\ldots,-1,1,\ldots,N\\}$. For fixed
$\alpha\in\mathbb{R}^{m}$ we define linear operators
$B_{\ell}(\alpha):\mathbb{R}^{n}\to\mathbb{R}^{n}$ as follows:

$B_{\ell}(\alpha)x=x-\frac{h_{x}(u(t_{\ell}(\alpha),\alpha),\alpha)x}{h_{x}(u(t_{\ell}(\alpha),\alpha),\alpha)\dot{u}(t_{\ell}(\alpha)^{-},\alpha)}[\dot{u}(t_{\ell}(\alpha)^{-},\alpha)-\dot{u}(t_{\ell}(\alpha)^{+},\alpha)].$
(5.1)

Note that (recall that $i\in\\{1,\ldots,N\\}$, $j\in\\{-1,\ldots,-M\\}$):

$\begin{array}[]{l}h_{x}(u(t_{i}(\alpha),\alpha),\alpha)\dot{u}(t_{i}(\alpha)^{-},\alpha)=h_{x}(u(t_{i}(\alpha),\alpha),\alpha)f_{i-1}(u(t_{i}(\alpha),\alpha),\alpha)\neq
0\\\
h_{x}(u(t_{j}(\alpha),\alpha),\alpha)\dot{u}(t_{j}(\alpha)^{-},\alpha)=h_{x}(u(t_{j}(\alpha),\alpha),\alpha)f_{j}(u(t_{j}(\alpha),\alpha),\alpha)\neq
0\end{array}$

We have the following:

######  Proposition 5.1.

For any $\ell$ and $\alpha\in\mathbb{R}^{m}$, $x\mapsto B_{\ell}(\alpha)x$ are
invertible linear maps. Moreover $x_{+,\xi_{+}}(t,0,\alpha,0)$ is a solution
of

$\begin{array}[]{l}\dot{x}=A(t,\alpha)x:=\left\\{\begin{array}[]{ll}f_{i-1,x}(u(t,\alpha),\alpha)x&\hbox{if
$t_{i-1}(\alpha)<t<t_{i}(\alpha)$}\\\ f_{N,x}(u(t,\alpha),\alpha)x&\hbox{if
$t>t_{N}(\alpha)$}\end{array}\right.\\\
x(t_{i}(\alpha)^{+})=B_{i}(\alpha)x(t_{i}(\alpha)^{-}),\end{array}$ (5.2)

which is $C^{1}$ for $t\neq t_{i}(\alpha)$, bounded for $t\geq 0$ and can be
assumed to be right-continuous at $t=t_{i}(\alpha)$. Similarly
$x_{-,\xi_{-}}(t,0,\alpha,0)$ is a solution of

$\begin{array}[]{l}\dot{x}=A(t,\alpha)x:=\left\\{\begin{array}[]{ll}f_{j+1,x}(u(t,\alpha),\alpha)x&\hbox{if
$t_{j}(\alpha)<t<t_{j+1}(\alpha)$}\\\ f_{-M,x}(u(t,\alpha),\alpha)x&\hbox{if
$t<t_{-M}(\alpha)$}\end{array}\right.\\\
x(t_{j}(\alpha)^{+})=B_{j}(\alpha)x(t_{j}(\alpha)^{-}),\end{array}$ (5.3)

which is $C^{1}$ for $t\neq t_{j}(\alpha)$,bounded for $t\leq 0$ and can be
assumed to be left-continuous at $t=t_{j}(\alpha)$.

###### Proof.

First we prove that $B_{\ell}(\alpha):\mathbb{R}^{n}\to\mathbb{R}^{n}$ is
invertible. If
$\dot{u}(t_{\ell}(\alpha)^{-},\alpha)=\dot{u}(t_{\ell}(\alpha)^{+},\alpha)$,
there is nothing to prove since $B_{\ell}(\alpha)x=x$. So, suppose that
$\dot{u}(t_{\ell}(\alpha)^{-},\alpha)\neq\dot{u}(t_{\ell}(\alpha)^{+},\alpha)$
and $B_{\ell}(\alpha)x=0$. Then

$x=\mu[\dot{u}(t_{\ell}(\alpha)^{-},\alpha)-\dot{u}(t_{\ell}(\alpha)^{+},\alpha)]$

where
$\mu=\frac{h_{x}(u(t_{\ell}(\alpha),\alpha),\alpha)x}{h_{x}(u(t_{\ell}(\alpha),\alpha),\alpha)\dot{u}(t_{\ell}(\alpha)^{-},\alpha)}$.
Hence:

$\begin{array}[]{l}B_{\ell}(\alpha)x=\mu[\dot{u}(t_{\ell}(\alpha)^{-},\alpha)-\dot{u}(t_{\ell}(\alpha)^{+},\alpha)]\\\
\qquad-\mu\frac{h_{x}(u(t_{\ell}(\alpha),\alpha),\alpha)[\dot{u}(t_{\ell}(\alpha)^{-},\alpha)-\dot{u}(t_{\ell}(\alpha)^{+},\alpha)]}{h_{x}(u(t_{i}(\alpha),\alpha),\alpha)\dot{u}(t_{\ell}(\alpha)^{-},\alpha)}[\dot{u}(t_{\ell}(\alpha)^{-},\alpha)-\dot{u}(t_{\ell}(\alpha)^{+},\alpha)]\\\
=\mu\frac{h_{x}(u(t_{\ell}(\alpha),\alpha),\alpha)\dot{u}(t_{\ell}(\alpha)^{+},\alpha)}{h_{x}(u(t_{\ell}(\alpha),\alpha),\alpha)\dot{u}(t_{\ell}(\alpha)^{-},\alpha)}[\dot{u}(t_{\ell}(\alpha)^{-},\alpha)-\dot{u}(t_{\ell}(\alpha)^{+},\alpha)].\end{array}$

But then $\mu=0$ since

$\begin{array}[]{l}h_{x}(u(t_{\ell}(\alpha),\alpha),\alpha)\dot{u}(t_{\ell}(\alpha)^{+},\alpha)\\\
=\left\\{\begin{array}[]{ll}h_{x}(u(t_{i}(\alpha),\alpha),\alpha)f_{i}(u(t_{i}(\alpha),\alpha),\alpha)&\hbox{if
$\ell=i>0$}\\\
h_{x}(u(t_{j}(\alpha),\alpha),\alpha)f_{j+1}(u(t_{j}(\alpha),\alpha),\alpha)&\hbox{if
$\ell=j<0$}\end{array}\right.\neq 0.\end{array}$

Next, $x_{+}(t,\xi_{+},\alpha,0)$ is a continuous, piecewise $C^{r}$, solution
of the differential equation

$\dot{x}=\left\\{\begin{array}[]{ll}f_{i-1}(x,\alpha)&\hbox{if
$t^{*}_{i-1}(\xi_{+},\alpha,0)<t<t^{*}_{i}(\xi_{+},\alpha,0)$}\\\
f_{N}(x,\alpha)&\hbox{if $t>t^{*}_{N}(\xi_{+},\alpha,0)$.}\end{array}\right.$

Hence, for $t^{*}_{i}(\xi_{+},\alpha,0)<t<t^{*}_{i+1}(\xi_{+},\alpha,0)$, when
$i\leq N-1$, or $t\geq t^{*}_{N}(\xi_{+},\alpha,0)$ when $i=N$, we have

$x_{+}(t,\xi_{+},\alpha,0)=x_{+}(t^{*}_{i}(\xi_{+},\alpha,0)^{-},\xi_{+},\alpha,0)+\int_{t^{*}_{i}(\xi_{+},\alpha,0)}^{t}f_{i}(x_{+}(s,\xi_{+},\alpha,0),\alpha)ds.$

Differentiating with respect to $\xi_{+}$ we get, for the same values of $t$:

$\begin{array}[]{l}x_{+,\xi_{+}}(t,\xi_{+},\alpha,0)=x_{+,\xi_{+}}(t^{*}_{i}(\xi_{+},\alpha,0)^{-},\xi_{+},\alpha,0)\\\
\displaystyle[\dot{x}_{+}(t^{*}_{i}(\xi_{+},\alpha,0)^{-},\xi_{+},\alpha,0)-\dot{x}_{+}(t^{*}_{i}(\xi_{+},\alpha,0)^{+},\xi_{+},\alpha,0)]\frac{\partial
t^{*}_{i}}{\partial\xi_{+}}(\xi_{+},\alpha,0)\\\
\displaystyle+\int_{t_{i}^{*}(\xi_{+},\alpha,0)}^{t}f_{i,x}(x_{+}(s,\xi_{+},\alpha,0),\alpha)x_{\xi_{+}}(s,\xi_{+},\alpha,0)ds.\end{array}$

and then,

$\begin{array}[]{l}x_{+,\xi_{+}}(t^{*}_{i}(\xi_{+},\alpha,0)^{+},\xi_{+},\alpha,0)=x_{+,\xi_{+}}(t^{*}_{i}(\xi_{+},\alpha,0)^{-},\xi_{+},\alpha,0)+\\\
\displaystyle+[f_{i-1}(\hat{x}^{*}_{i}(\xi_{+},\alpha),\alpha)-f_{i}(\hat{x}^{*}_{i}(\xi_{+},\alpha),\alpha)]\frac{\partial
t^{*}_{i}}{\partial\xi_{+}}(\xi_{+},\alpha,0)\end{array}$ (5.4)

where we write for simplicity

$\hat{x}^{*}_{i}(\xi_{+},\alpha):=x_{+}(t^{*}_{i}(\xi_{+},\alpha,0),\xi_{+},\alpha,0).$
(5.5)

Now, from Remark 4.6 we see that

$\frac{\partial
t^{*}_{i}}{\partial\xi_{+}}(\xi_{+},\alpha,0)=-\frac{h_{x}(\hat{x}^{*}_{i}(\xi_{+},\alpha),\alpha)x_{+,\xi_{+}}(t^{*}(\xi_{+},\alpha,0)^{-},\alpha,0)}{h_{x}(\hat{x}^{*}_{i}(\xi_{+},\alpha),\alpha)f_{i-1}(\hat{x}^{*}_{i}(\xi_{+},\alpha),\alpha)}.$

Hence

$x_{\xi_{+}}(t^{*}_{i}(\xi_{+},\alpha,0)^{+},\xi_{+},\alpha,0)=B_{i}(\xi_{+},\alpha)x_{\xi_{+}}(t^{*}_{i}(\xi_{+},\alpha,0)^{-},\xi_{+},\alpha,0)$
(5.6)

where

$\begin{array}[]{l}B_{i}(\xi_{+},\alpha)x=\\\
x-[f_{i-1}(\hat{x}^{*}(\xi_{+},\alpha),\alpha)-f_{i}(\hat{x}^{*}(\xi_{+},\alpha),\alpha)]\frac{h_{x}(\hat{x}^{*}_{i}(\xi_{+},\alpha),\alpha)x}{h_{x}(\hat{x}^{*}_{i}(\xi_{+},\alpha),\alpha)f_{i-1}(\hat{x}^{*}_{i}(\xi_{+},\alpha),\alpha)}.\end{array}$

Taking $\xi_{+}=0$, we see that $x_{\xi_{+}}(t,0,\alpha,0)$ is a solution, for
$t\geq 0$ of (5.2) where $B_{i}(\alpha)$ is as in (5.1).

Following a similar argument we see that $x_{\xi_{-}}(t,0,\alpha,0)$ is a
solution, for $t\leq 0$ of (5.3) where $B_{j}(\alpha)$ is as in (5.1).

Finally we prove that $x_{\xi_{+}}(t,0,\alpha,0)$ is bounded for $t\geq 0$. It
is enough to prove this for $t\geq T_{+}$. From Theorem 4.3 we know that
$x_{\xi_{+}}(t,0,\alpha,0)$ is a solution of

$\dot{x}=f_{N,x}(u_{N}^{+}(t,\alpha),\alpha)x$

such that

$\sup_{t\geq T_{+}}|x(t)|e^{-\beta(t-T_{+})}<\infty.$

So, $x_{\xi_{+}}(t,0,\alpha,0)e^{-\beta(t-T_{+})}$ is a bounded solution of
the linear system

$\dot{x}=[f_{N,x}(u_{N}^{+}(t,\alpha),\alpha)-\beta\mathbb{I}]x$ (5.7)

whose fundamental matrix on $t\geq T_{+}$ is
$U^{+}_{N}(t,\alpha)U^{+}_{N}(T_{+},\alpha)^{-1}e^{-\beta(t-T_{+})}$.
According to Lemma 4.1 (5.7) has an exponential dichotomy on $t\geq T_{+}$
with projection $P_{+}(\alpha)$ and exponent $\delta-\beta$. Then we have
$x_{\xi_{+}}(T_{+},0,\alpha,0)\in{\cal R}P_{+}(\alpha)$. But then
$x_{\xi_{+}}(t,0,\alpha,0)$ is bounded for $t\geq T_{+}$ because ${\cal
R}P_{+}(\alpha)$ is the space of initial conditions of solutions of
$\dot{x}=f_{N,x}(u_{N}^{+}(t,\alpha),\alpha)x$ that are bounded for $t\geq
T_{+}$. A similar argument shows that $x_{\xi_{-}}(t,0,\alpha,0)$ is bounded
for $t\leq 0$. The proof is complete. ∎

In the next proposition we show that $\dot{u}(t,\alpha)$ is a nontrivial
bounded solution of (5.2) for $t\geq 0$ (resp. (5.3) for $t\leq 0$).

######  Proposition 5.2.

For $t\geq 0$, resp. $t\leq 0$, the function

$\dot{u}(t,\alpha)=\left\\{\begin{array}[]{ll}\dot{u}_{i-1}^{+}(t,\alpha)&\hbox{for
$t_{i-1}(\alpha)<t<t_{i}(\alpha)$}\\\ \dot{u}^{+}_{N}(t,\alpha)&\hbox{for
$t>t_{N}(\alpha)$,}\end{array}\right.$

resp.

$\dot{u}(t,\alpha)=\left\\{\begin{array}[]{ll}\dot{u}_{j+1}^{-}(t,\alpha)&\hbox{for
$t_{j}(\alpha)<t<t_{j+1}\alpha)$}\\\ \dot{u}^{-}_{-M}(t,\alpha)&\hbox{for
$t<t_{-M}(\alpha)$,}\end{array}\right.$

is a solution of (5.2) (resp. (5.3)) bounded on $t\geq 0$ (resp, $t\leq 0$)
where $B_{\ell}(\alpha)$ is as in (5.1).

###### Proof.

We already know that $\dot{u}(t)$ satisfies (5.2), for $t\geq 0$ and (5.3) for
$t\leq 0$, $t\neq t_{\ell}$. We prove that
$\dot{u}(t_{\ell}(\alpha)^{+},\alpha)=B_{\ell}(\alpha)\dot{u}(t_{\ell}(\alpha)^{-},\alpha)$.
We have

$\begin{array}[]{l}B_{\ell}(\alpha)\dot{u}(t_{\ell}(\alpha)^{-},\alpha)=\dot{u}(t_{\ell}(\alpha)^{-},\alpha)\\\
\qquad-\frac{h_{x}(u(t_{\ell}(\alpha),\alpha),\alpha)\dot{u}(t_{\ell}(\alpha)^{-},\alpha)}{h_{x}(u(t_{\ell}(\alpha),\alpha),\alpha)\dot{u}(t_{\ell}(\alpha)^{-},\alpha)}[\dot{u}(t_{\ell}(\alpha)^{-},\alpha)-\dot{u}(t_{\ell}(\alpha)^{+},\alpha)]\\\
=\dot{u}(t_{\ell}(\alpha)^{+},\alpha)\end{array}$

The proof is complete. ∎

##  6 The Melnikov condition

First we recall that $P_{+}(y)$ is the projections of the exponential
dichotomy on $t\geq 0$, of the linear system (4.1) with constant $K$ and
exponent $\delta$. Then, from Lemma 4.1, we see that (4.4) has an exponential
dichotomy on $t\geq t_{N}(y)$ with exponent $\delta$ and projection

$\begin{array}[]{c}U^{+}_{N}(t_{N}(y),y)U^{+}_{N}(T_{+},y)^{-1}P_{+}(y)U^{+}_{N}(T_{+},y)U^{+}_{N}(t_{N}(y),y)^{-1}\\\
=X_{+}(t_{N}(y),y)X_{+}(T_{+},y)^{-1}P_{+}(y)X_{+}(T_{+},y)X_{+}(t_{N}(y),y)^{-1}\end{array}$
(6.1)

the equality following from (3.9) and $T_{+}>t_{N}(y)$.

Similarly, the linear system (4.5) has an exponential dichotomy on $t\leq
t_{-M}(y)$ with exponent $\delta$ and projection

$\begin{array}[]{c}U^{-}_{-M}(t_{-M}(y),y)U^{-}_{-M}(T_{-},y)^{-1}P_{-}(y)U^{-}_{-M}(T_{-},y)U^{-}_{-M}(t_{-M}(y),y)^{-1}\\\
=X_{-}(t_{-M}(y),y)X_{-}(T_{-},y)^{-1}P_{-}(y)X_{-}(T_{-},y)X_{-}(t_{-M}(y),y)^{-1}\end{array}$
(6.2)

where

$X_{+}(t,y)=\left\\{\begin{array}[]{ll}U_{0}(t,y)&\hbox{if $0\leq
t<t_{1}(y)$}\\\
U_{i}(t,y)U_{i}(t_{i},y)^{-1}B_{i}(y)X_{+}(t_{i}^{-},y)&\hbox{if $t_{i}(y)\leq
t<t_{i+1}(y)$}\\\
U_{N}(t,y)U_{N}(t_{N},y)^{-1}B_{N}(y)X_{+}(t_{N}^{-},y)&\hbox{if $t\geq
t_{N}(y)$}\end{array}\right.$

and

$X_{-}(t,y)=\left\\{\begin{array}[]{ll}U_{0}(t,y)&\hbox{if $t_{-1}(y)<t\leq
0$}\\\ U_{j}(t,y)U_{j}(t_{j},y)^{-1}B_{j}^{-1}(y)X_{-}(t_{j}^{+},y)&\hbox{if
$t_{j-1}(y)<t\leq t_{j}(y)$}\\\
U_{-M}(t,y)U_{-M}(t_{-M},y)^{-1}B_{-M}(y)^{-1}X_{-}(t_{-M}^{+},y)&\hbox{if
$t\leq t_{-M}(y)$}\end{array}\right.$

are the fundamental matrix of $\dot{x}=A(t,y)x$, where $A(t,y)$ is as in
(3.6).

From Lemma 3.3–3.4 and (6.1)-(6.2) we obtain the following

######  Proposition 6.1.

For any $\alpha\in\mathbb{R}^{m}$, the discontinuous linear system (5.2)
(resp. (5.3)) has an exponential dichotomy on $\mathbb{R}_{+}$, (resp.
$\mathbb{R}_{-}$) with projections $Q_{+}(\alpha)$, resp. $Q_{-}(\alpha)$,
given by

$\begin{array}[]{l}Q_{+}(\alpha)=X_{+}(T_{+},\alpha)^{-1}P_{+}(\alpha)X_{+}(T_{+},\alpha)\\\
\\\
Q_{-}(\alpha)=X_{-}(T_{-},\alpha)^{-1}P_{-}(\alpha)X_{-}(T_{-},\alpha).\end{array}$

Moreover ${\cal R}Q_{+}(\alpha)$ (resp. ${\cal N}Q_{-}(\alpha)$) is the space
of initial conditions of solutions of (5.2), resp. (5.3), right-continuous,
when $t\geq 0$ (resp. left-continuous, when $t\leq 0$) and bounded on
$\mathbb{R}_{+}$, (resp, on $\mathbb{R}_{-}$).

For simplicity we write $Q_{\pm}=Q_{\pm}(y_{0})$.

We assume the following condition holds:

  * $A_{5})$

$\dim{\cal R}Q_{+}\cap{\cal N}Q_{-}=d\leq m$.

From Proposition 5.2 we see that $\dot{u}(0,y_{0})\in{\cal R}Q_{+}\cap{\cal
N}Q_{-}$ so

$1\leq\dim[{\cal R}Q_{+}+{\cal N}Q_{-}]^{\perp}=d.$

Next, from $A_{3})$ we know that $\dim{\cal R}Q_{+}=k$ and $\dim{\cal
N}Q_{-}=n-k$, hence $d\leq\min\\{k,n-k\\}$.

Let $\psi_{1},\ldots,\psi_{d}\in\mathbb{R}^{n}$ be such that $[{\cal
R}Q_{+}+{\cal N}Q_{-}]^{\perp}=\hbox{\rm span}\\{\psi_{1},\ldots,\psi_{d}\\}$.
Without loss of generality we assume that $(\psi_{1},\ldots,\psi_{d})$ is an
orthonormal set.

The purpose of this section is the to prove the following

######  Theorem 6.2.

Suppose that $A_{1})-A_{5})$ hold. Suppose further that the matrix
$[\psi_{j}^{T}[w^{-}_{0,y}(y_{0})-w^{+}_{0,y}(y_{0})]_{j=1,\ldots,d}$ has rank
$d$. Then there exists $\rho>0$ and $\varepsilon_{0}>0$ such that for
$0\leq\varepsilon\leq\varepsilon_{0}$ system (1.7) has a $(m-d)$-dimensional
manifold of continuous, piecewise $C^{r}$ solutions
$(x(t,\varepsilon),y(t,\varepsilon))$ such that

$\begin{array}[]{l}\sup_{t\in\mathbb{R}}|x(t,\varepsilon)-u(t,y(t,\varepsilon))<\rho,\\\
\displaystyle\sup_{t\in\mathbb{R}}|x(t,\varepsilon)-u(t,y(t,\varepsilon))|\to
0\end{array}$

as $\varepsilon\to 0$.

###### Proof.

First, we apply Lemma 3.5 to obtain another expression of
$x_{+}(t,\xi_{+},\alpha,\varepsilon)$ (resp.
$x_{-}(t,\xi_{-},\alpha,\varepsilon)$). We know that, for $t\geq 0$,

$z_{+}(t)=x_{+}(t,\xi_{+},\alpha,\varepsilon)-u(t,y_{+}(t,\xi_{+},\alpha,\varepsilon))$

is a bounded and continuous solution of the differential equation

$\begin{array}[]{l}\dot{z}=f(x_{+}(t,\xi_{+},\alpha,\varepsilon),y_{+}(t,\xi_{+},\alpha,\varepsilon))-f(u(t,y_{+}(t,\xi_{+},\alpha,\varepsilon)),y_{+}(t,\xi_{+},\alpha,\varepsilon))\\\
\hskip 20.0pt-\varepsilon
u_{y}(t,y_{+}(t,\xi_{+},\alpha,\varepsilon))g(x_{+}(t,\xi_{+},\alpha,\varepsilon),y_{+}(t,\xi_{+},\alpha,\varepsilon))\end{array}$

that we write:

$\left\\{\begin{array}[]{l}\dot{z}-A(t,\alpha)z=b_{+}(t)\\\
z(t_{i}(\alpha)^{+})=z(t_{i}(\alpha)^{-}),\quad
i=1,\ldots,N\end{array}\right.$

where $A(t,y)$ has been defined in (3.6) and

$\begin{array}[]{l}b_{+}(t)=f(x_{+}(t,\xi_{+},\alpha,\varepsilon),y_{+}(t,\xi_{+},\alpha,\varepsilon))-f(u(t,y_{+}(t,\xi_{+},\alpha,\varepsilon)),y_{+}(t,\xi_{+},\alpha,\varepsilon))\\\
\hskip 20.0pt-A(t,\alpha)x_{+}(t,\xi_{+},\alpha,\varepsilon)-\varepsilon
u_{y}(t,y_{+}(t,\xi_{+},\alpha,\varepsilon))g(x_{+}(t,\xi_{+},\alpha,\varepsilon),y_{+}(t,\xi_{+},\alpha,\varepsilon)).\end{array}$
(6.3)

Note that

$(u(t,y_{+}(t,\xi_{+},\alpha,0)),y_{+}(t,\xi_{+},\alpha,0))=(u(t,\alpha),\alpha)$

and, according to $a_{2})$,

$\begin{array}[]{l}h_{x}(u(t_{i}(\alpha),\alpha)f_{i-1}(u(t_{i}(\alpha),\alpha),\alpha)>2\eta,\\\
h_{x}(u(t_{i}(\alpha),\alpha)f_{i}((u(t_{i}(\alpha),\alpha)>2\eta.\end{array}$

Hence, for $\varepsilon$ sufficiently small there exist $C^{r}$-functions
$\tilde{t}_{i}(\xi_{+},\alpha,\varepsilon)$ such that

$\begin{array}[]{l}|\tilde{t}_{i}(\xi_{+},\alpha,\varepsilon)-t_{i}(\alpha)|\to
0,\hbox{ as $\varepsilon\to 0$}\\\
h(u(\tilde{t}_{i}(\xi_{+},\alpha,\varepsilon),y_{+}(\tilde{t}_{i}(\xi_{+},\alpha,\varepsilon),\xi_{+},\alpha,\varepsilon)),y_{+}(\tilde{t}_{i}(\xi_{+},\alpha,\varepsilon),\xi_{+},\alpha,\varepsilon))=c_{i}\\\
c_{i-1}<h(u(t,y_{+}(t,\xi_{+},\alpha,\varepsilon)),y_{+}(t,\xi_{+},\alpha,\varepsilon))<c_{i},\\\
\qquad\hbox{ for
$\tilde{t}_{i-1}(\xi_{+},\alpha,\varepsilon)<t<\tilde{t}_{i}(\xi_{+},\alpha,\varepsilon)$}\\\
h_{x}(x,y)f_{i-1}(x,y)_{\vrule_{x=u(\tilde{t}_{i}(\xi_{+},\alpha,\varepsilon),y_{+}(\tilde{t}_{i}(\xi_{+},\alpha,\varepsilon),\xi_{+},\alpha,\varepsilon))\atop
y=y_{+}(\tilde{t}_{i}(\xi_{+},\alpha,\varepsilon)\hfill}}>\eta\\\
h_{x}(x,y)f_{i}(x,y)_{\vrule_{x=u(\tilde{t}_{i}(\xi_{+},\alpha,\varepsilon),y_{+}(\tilde{t}_{i}(\xi_{+},\alpha,\varepsilon),\xi_{+},\alpha,\varepsilon))\atop
y=y_{+}(\tilde{t}_{i}(\xi_{+},\alpha,\varepsilon)\hfill}}>\eta\end{array}$

for any $i=1,\ldots,N$ and uniformly with respect to $(\xi_{+},\alpha)$. Then

$\begin{array}[]{l}f(x_{+}(t,\xi_{+},\alpha,\varepsilon),y_{+}(t,\xi_{+},\alpha,\varepsilon))\\\
=\left\\{\begin{array}[]{ll}f_{i-1}(x_{+}(t,\xi_{+},\alpha,\varepsilon),y_{+}(t,\xi_{+},\alpha,\varepsilon))&\hbox{if
$t^{*}_{i-1}(\xi_{+},\alpha,\varepsilon)\leq
t<t^{*}_{i}(\xi_{+},\alpha,\varepsilon)$}\\\
f_{N}(x_{+}(t,\xi_{+},\alpha,\varepsilon),y_{+}(t,\xi_{+},\alpha,\varepsilon))&\hbox{if
$t\geq t^{*}_{N}(\xi_{+},\alpha,\varepsilon)$}\end{array}\right.\end{array}$

and

$\begin{array}[]{l}f(u(t,\xi_{+},\alpha,\varepsilon),y_{+}(t,\xi_{+},\alpha,\varepsilon))\\\
=\left\\{\begin{array}[]{ll}f_{i-1}(u(t,\xi_{+},\alpha,\varepsilon),y_{+}(t,\xi_{+},\alpha,\varepsilon))&\hbox{if
$\tilde{t}_{i-1}(\xi_{+},\alpha,\varepsilon)\leq
t<\tilde{t}_{i}(\xi_{+},\alpha,\varepsilon)$}\\\
f_{N}(u(t,\xi_{+},\alpha,\varepsilon),y_{+}(t,\xi_{+},\alpha,\varepsilon))&\hbox{if
$t\geq\tilde{t}_{N}(\xi_{+},\alpha,\varepsilon)$.}\end{array}\right.\end{array}$

According to Lemma 3.5, with $\tau=0$, we see that

$\begin{array}[]{l}x_{+}(t,\xi_{+},\alpha,\varepsilon)=u(t,y_{+}(t,\xi_{+},\alpha,\varepsilon))+X_{+}(t,\alpha)\tilde{\xi}_{+}\\\
\qquad\displaystyle+\int_{0}^{t}X_{+}(t,\alpha)Q_{+}(\alpha)X_{+}(s,\alpha)^{-1}b_{+}(s)ds\\\
\qquad\displaystyle-\int_{t}^{\infty}X_{+}(t,\alpha)(\mathbb{I}-Q_{+}(\alpha))X_{+}(s,\alpha)^{-1}b_{+}(s)ds\end{array}$
(6.4)

where

$\tilde{\xi}_{+}=Q_{+}(\alpha)[x_{+}(0,\xi_{+},\alpha,\varepsilon)-u(0^{+},y_{+}(0,\xi_{+},\alpha,\varepsilon))]\in{\cal
R}Q_{+}(\alpha).$

Note that

$|y_{+}(0,\xi_{+},\alpha,\varepsilon)-\alpha|\leq|\varepsilon|\int_{0}^{T_{+}}g(x_{+}(t,\xi_{+},\alpha,\varepsilon),y_{+}(t,\xi_{+},\alpha,\varepsilon),\varepsilon)|dt\to
0$

as $\varepsilon\to 0$, uniformly with respect to $(\xi_{+},\alpha)$. We prove
the following

Claim: For $\varepsilon$ sufficiently small, the map
$(\xi_{+},\alpha)\mapsto(\tilde{\xi}_{+},\tilde{\alpha})$, with
$\tilde{\alpha}=y_{+}(0,\xi_{+},\alpha,\varepsilon)$ from ${\cal
R}P_{+}(\alpha)\times\mathbb{R}^{m}$ into ${\cal
R}Q_{+}(\alpha)\times\mathbb{R}^{m}$ is linearly invertible.

Indeed, for $\varepsilon=0$ the above map reduces to

$(\xi_{+},\alpha)\mapsto(\tilde{\xi}_{+},\alpha)$

where

$\tilde{\xi}_{+}=Q_{+}(\alpha)[x_{+}(0,\xi_{+},\alpha,0)-u(0^{+},\alpha)].$

Hence

$\frac{\partial\tilde{\xi}_{+}}{\partial\xi_{+}}(0,\alpha,0)=Q_{+}(\alpha)x_{+,\xi_{+}}(0,0,\alpha,0).$

Now, from Proposition 5.1 we know that, for any $\xi_{+}\in{\cal
R}P_{+}(\alpha)$,

$z(t):=x_{+,\xi_{+}}(t,0,\alpha,0)\xi_{+}$

is a right-continuous solution, bounded on $t\geq 0$, of

$\begin{array}[]{l}\dot{z}(t)=A(t,\alpha)z(t)\\\
z(t_{i}(\alpha)^{+})=B_{i}(\alpha)z(t_{i}^{-}(\alpha))\\\
P_{+}(\alpha)z(T_{+})=\xi_{+}\end{array}$

the last relation following differentiating the equality

$P_{+}(\alpha)[x_{+}(T_{+},\xi_{+},\alpha,0)-u(T_{+},\alpha)]=\xi_{+}$

with respect to $\xi_{+}$ at $\xi_{+}=0$ (see Theorem (4.3)). In particular we
have $z(t)=X_{+}(t,\alpha)X_{+}(\tau,\alpha)^{-1}z(\tau)$ for any $\tau\geq 0$
(see Remark 3.2). From Lemma 4.1 it follows that the linear system (4.4), with
$y=\alpha$, has an exponential dichotomy on $t\geq t_{N}(\alpha)$ with
projection (see also (3.9))

$\begin{array}[]{l}{\cal
P}_{+}(\alpha)=U^{+}_{N}(t_{N}(\alpha),\alpha)U^{+}_{N}(T_{+},\alpha)^{-1}P_{+}(\alpha)U^{+}_{N}(T_{+},\alpha)U^{+}_{N}(t_{N}(\alpha),\alpha)^{-1}\\\
=X_{+}(t_{N}(\alpha),\alpha)X_{+}(T_{+},\alpha)^{-1}P_{+}(\alpha)X_{+}(T_{+},\alpha)X_{+}(t_{N}(\alpha),\alpha)^{-1}\end{array}$

as $T_{+}<t_{N}(y)$. Hence

$\begin{array}[]{l}X_{+}(T_{+},\alpha)X_{+}(t_{N}(\alpha),\alpha)^{-1}{\cal
P}_{+}(\alpha)=P_{+}(\alpha)X_{+}(T_{+},\alpha)X_{+}(t_{N}(\alpha),\alpha)^{-1}.\end{array}$

As a consequence

$\begin{array}[]{l}X_{+}(T_{+},\alpha)X_{+}(t_{N}(\alpha),\alpha)^{-1}{\cal
P}_{+}(\alpha)z(t_{N}(\alpha))\\\
=P_{+}(\alpha)X_{+}(T_{+},\alpha)X_{+}(t_{N}(\alpha),\alpha)^{-1}z(t_{N}(\alpha))=P_{+}(\alpha)z(T_{+})=\xi_{+}.\end{array}$

So, $z(t):=x_{+,\xi_{+}}(t,0,\alpha,0)\xi_{+}$ is a right-continuous solution,
bounded on $t\geq 0$

$\begin{array}[]{l}\dot{z}(t)=A(t,\alpha)z(t)\\\
z(t_{i}(\alpha)^{+})=B_{i}(\alpha)z(t_{i}^{-}(\alpha))\\\ {\cal
P}_{+}(\alpha)z(t_{N}(\alpha))=X_{+}(t_{N}(\alpha),\alpha)X_{+}(T_{+},\alpha)^{-1}\xi_{+}\end{array}$

for any $\xi_{+}\in{\cal R}P_{+}(\alpha)$. Now, with reference to Lemmas 3.3,
3.5 with $\tau=t_{N}(\alpha)$, we have

$\begin{array}[]{l}\tilde{\cal
P}_{+}^{t_{N}(\alpha)}(\alpha)=X_{+}(t_{N}(\alpha),\alpha)Q_{+}(\alpha)X_{+}(t_{N}(\alpha),\alpha)^{-1}\\\
=X_{+}(t_{N}(\alpha),\alpha)X_{+}(T_{+},\alpha)^{-1}P_{+}(\alpha)X_{+}(T_{+},\alpha)X_{+}(t_{N}(\alpha),\alpha)^{-1}\\\
={\cal P}_{+}(\alpha).\end{array}$

Hence $x_{+,\xi_{+}}(t,0,y_{0},0)$ is a bounded solution of

$\begin{array}[]{l}\dot{z}(t)=A(t,y_{0})z(t)\\\
z(t_{i}(\alpha)^{+})=B_{i}(\alpha)z(t_{i}^{-}(\alpha))\\\ \tilde{\cal
P}_{+}^{t_{N}(\alpha)}(\alpha)z(t_{N}(\alpha))=X_{+}(t_{N}(\alpha),\alpha)X_{+}(T_{+},\alpha)^{-1}\xi_{+}.\end{array}$

From Lemma 3.5, equation (3.17), we get then, for any $\xi_{+}\in{\cal
R}P_{+}(\alpha)$, using again the right-continuity of $X_{+}(t,\alpha)$:

$\begin{array}[]{l}\displaystyle
x_{+,\xi_{+}}(t,0,\alpha,0)\xi_{+}=X_{+}(t,\alpha)Q_{+}(\alpha)X_{+}(t_{N}(\alpha),\alpha)^{-1}X_{+}(t_{N}(\alpha),\alpha)X_{+}(T_{+},\alpha)^{-1}\xi_{+}\\\
\displaystyle=X_{+}(t,\alpha)Q_{+}(\alpha)X_{+}(T_{+},\alpha)^{-1}\xi_{+}\end{array}$

So,

$\frac{\partial\tilde{\xi}_{+}}{\partial\xi_{+}}(0,\alpha,0)=Q_{+}(\alpha)x_{+,\xi_{+}}(0,0,\alpha,0)=Q_{+}(\alpha)X_{+}(T_{+},\alpha)^{-1}\xi_{+}=X_{+}(T_{+},\alpha)^{-1}P_{+}(\alpha)\xi_{+}.$

Hence $\frac{\partial\tilde{\xi}_{+}}{\partial\xi_{+}}(0,\alpha,0)$ is an
isomorphism from ${\cal R}P_{+}(\alpha)$ into ${\cal R}Q_{+}(\alpha)$ and the
Claim is proved.

Similarly we see that

$\begin{array}[]{l}x_{-}(t,\xi_{-},\alpha,\varepsilon)=u(t,y_{-}(t,\xi_{-},\alpha_{-},\varepsilon))+X_{-}(t,\alpha)\tilde{\xi}_{-}\\\
\quad\displaystyle+\int_{-\infty}^{t}X_{-}(t,\alpha)Q_{-}(\alpha)X_{-}(s,\alpha)^{-1}b_{-}(s)ds\\\
\quad\displaystyle-\int_{t}^{0}X_{-}(t,\alpha)(\mathbb{I}-Q_{-}(\alpha))X_{-}(s,\alpha)^{-1}b_{-}(s)ds\end{array}$
(6.5)

where

$\begin{array}[]{l}b_{-}(t)=f(x_{-}(t,\xi_{-},\alpha,\varepsilon),y_{-}(t,\xi_{-},\alpha,\varepsilon))-f(u_{0}(t,y_{-}(t,\xi_{-},\alpha,\varepsilon)),y_{-}(t,\xi_{-},\alpha,\varepsilon))\\\
\hskip 20.0pt-{\cal
A}(t,\alpha)x_{-}(t,\xi_{-},\alpha,\varepsilon)-\varepsilon
u_{0,y}(t,y_{-}(t,\xi_{-},\alpha,\varepsilon))g(x_{-}(t,\xi_{-},\alpha,\varepsilon),y_{-}(t,\xi_{-},\alpha,\varepsilon))\end{array}$

and

$\tilde{\xi}_{-}=[\mathbb{I}-Q_{-}(\alpha)][x_{-}(0,\xi_{-},\alpha,\varepsilon)-u_{0}(0^{-},y_{-}(0,\xi_{-},\alpha,\varepsilon))]\in{\cal
N}Q_{-}(\alpha).$

Again we see that

$|y_{-}(0,\xi_{-},\alpha,\varepsilon)-\alpha|\leq|\varepsilon|\int_{0}^{1}g(x_{-}(t,\xi,\alpha_{-},\varepsilon),y_{-}(t,\xi_{-},\alpha,\varepsilon),\varepsilon)|dt\to
0$

as $\varepsilon\to 0$, uniformly with respect to $(\xi_{-},\alpha)$ and the
map
$(\xi_{-},\alpha)\mapsto(\tilde{\xi}_{-},y_{-}(0,\xi_{-},\alpha,\varepsilon))$
from ${\cal N}P_{-}(\alpha)\times\mathbb{R}^{m}$ into ${\cal
N}Q_{-}(\alpha)\times\mathbb{R}^{m}$ is linearly invertible.

From (6.4)-(6.5) we get, for $|\alpha_{+}-y_{0}|+|\alpha_{-}-y_{0}|$
sufficiently small

$\begin{array}[]{l}x_{+}(0,\xi_{+},\alpha_{+},\varepsilon)-x_{-}(0,\xi_{-},\alpha_{-},\varepsilon)\\\
=u_{0}(0^{+},y_{+}(0,\xi_{+},\alpha_{+},\varepsilon))-u_{0}(0^{-},y_{-}(0,\xi_{-},\alpha_{-},\varepsilon))+\tilde{\xi}_{+}-\tilde{\xi}_{-}\\\
\displaystyle-\int_{0}^{\infty}(\mathbb{I}-Q_{+}(\alpha_{+}))X_{+}(s,\alpha_{+})^{-1}b_{+}(s)ds-\int_{-\infty}^{0}Q_{-}(\alpha_{-})X_{-}(s,\alpha_{-})^{-1}b_{-}(s)ds.\end{array}$
(6.6)

Hence the system

$\begin{array}[]{l}x_{+}(0,\xi_{+},\alpha_{+},\varepsilon)=x_{-}(0,\xi_{-},\alpha_{-},\varepsilon)\\\
y_{+}(0,\xi_{+},\alpha_{+},\varepsilon)=y_{-}(0,\xi_{-},\alpha_{-},\varepsilon)\end{array}$

is equivalent to

$\left\\{\begin{array}[]{l}\tilde{\xi}_{+}-\tilde{\xi}_{-}=u(0^{-},y_{-}(0,\xi_{-},\alpha_{-},\varepsilon))-u(0^{+},y_{+}(0,\xi_{+},\alpha_{+},\varepsilon))\\\
\displaystyle\qquad+\int_{0}^{\infty}(\mathbb{I}-Q_{+}(\alpha_{+}))X_{+}(s,\alpha_{+})^{-1}b_{+}(s)ds+\int_{-\infty}^{0}Q_{-}(\alpha_{-})X_{-}(s,\alpha_{-})^{-1}b_{-}(s)ds\\\
y_{+}(0,\xi_{+},\alpha_{+},\varepsilon)-y_{-}(0,\xi_{-},\alpha_{-},\varepsilon)=0.\end{array}\right.$
(6.7)

Let

$\begin{array}[]{l}\displaystyle
k(\xi_{+},\xi_{-},\alpha_{+},\alpha_{-},\varepsilon)=\int_{0}^{\infty}(\mathbb{I}-Q_{+}(\alpha_{+}))X_{+}(s,\alpha_{+})^{-1}b_{+}(s)ds\\\
\qquad\displaystyle+\int_{-\infty}^{0}Q_{-}(\alpha_{-})X_{-}(s,\alpha_{-})^{-1}b_{-}(s)ds.\end{array}$

Differentiating $b_{+}(t)=b_{+}(t,\xi_{+},\alpha_{+},\varepsilon)$ with
respect to $\xi_{+}$ at $\xi_{+}=0$, $\varepsilon=0$ and using
$x_{+}(t,0,\alpha_{+},0)=u(t,\alpha_{+})$,
$y_{+}(t,0,\alpha_{+},0)=\alpha_{+}$, we see that, for
$t_{i-1}(\alpha_{+})<t<t_{i}(\alpha_{+})$, we have

$\begin{array}[]{l}\frac{\partial
b_{+}}{\partial\xi_{+}}(t,0,\alpha_{+},0)=[f_{i-1,x}(u(t,\alpha_{+}),\alpha_{+})-A(t,\alpha_{+})]x_{+,\xi_{+}}(t,0,\alpha_{+},0)=0.\end{array}$

and for $t>t_{N}(\alpha_{+})$:

$\begin{array}[]{l}\frac{\partial
b_{+}}{\partial\xi_{+}}(t,\xi_{+},\alpha_{+},0)=[f_{N,x}(u(t,\alpha_{+}),\alpha_{+})-A(t,\alpha_{+})]x_{+,\xi_{+}}(t,0,\alpha_{+},0)=0.\end{array}$

Then

$\frac{\partial}{\partial\xi_{+}}\int_{0}^{\infty}(\mathbb{I}-P_{+})X_{+}(s)^{-1}b_{+}(s)ds=0$

and similarly

$\displaystyle\frac{\partial}{\partial\xi_{-}}\int_{-\infty}^{0}P_{-}X_{-}(s)^{-1}b_{-}(s)ds=0.$

As a consequence, on account of $u(0^{\pm},y)=w^{\pm}_{0}(y)$, (6.7) reads:

$\begin{array}[]{l}\tilde{\xi}_{+}-\tilde{\xi}_{-}=w^{-}_{0}(\alpha_{-})-w^{+}_{0}(\alpha_{+})+R_{1}(\tilde{\xi}_{+},\tilde{\xi}_{-},\alpha_{+},\alpha_{-},\varepsilon)\\\
\alpha_{+}=\alpha_{-}+R_{2}(\tilde{\xi}_{+},\tilde{\xi}_{-},\alpha_{+},\alpha_{-},\varepsilon)\end{array}$

where
$R_{1}(\tilde{\xi}_{+},\tilde{\xi}_{-},\alpha_{+},\alpha_{-},\varepsilon)=O(|\xi_{+}|^{2}+|\xi_{-}|^{2}+|\varepsilon|)$
and
$R_{2}(\tilde{\xi}_{+},\tilde{\xi}_{-},\alpha_{+},\alpha_{-},\varepsilon)=O(|\varepsilon|)$,
uniformly wth respect to $(\xi_{+},\xi_{-},\alpha_{+},\alpha_{-})$. Since
$(\xi_{+},\alpha_{+})\mapsto(\tilde{\xi}_{+},\tilde{\alpha}_{+})$ and
$(\xi_{-},\alpha_{-})\mapsto(\tilde{\xi}_{-},\tilde{\alpha}_{-})$ are linearly
invertible we see that $|\xi_{\pm}|=O(|\tilde{\xi}_{\pm}|)$ and hence (6.7)
reads:

$\begin{array}[]{l}\tilde{\xi}_{+}-\tilde{\xi}_{-}=w^{-}_{0}(\alpha_{-})-w^{+}_{0}(\alpha_{+})+\tilde{R}_{1}(\tilde{\xi}_{+},\tilde{\xi}_{-},\alpha_{+},\alpha_{-},\varepsilon)\\\
\alpha_{+}=\alpha_{-}+\tilde{R}_{2}(\tilde{\xi}_{+},\tilde{\xi}_{-},\alpha_{+},\alpha_{-},\varepsilon)\end{array}$
(6.8)

where
$\tilde{R}_{1}(\tilde{\xi}_{+},\tilde{\xi}_{-},\alpha_{+},\alpha_{-},\varepsilon)=O(|\tilde{\xi}_{+}|^{2}+|\tilde{\xi}_{-}|^{2}+|\varepsilon|)$
and
$\tilde{R}_{1}(\tilde{\xi}_{+},\tilde{\xi}_{-},\alpha_{+},\alpha_{-},\varepsilon)=O(|\varepsilon|)$,
uniformly wth respect to
$(\tilde{\xi}_{+},\tilde{\xi}_{-},\alpha_{+},\alpha_{-})$. Now we can write
$\tilde{\xi}_{+}=Q_{+}(\alpha)\tilde{\xi}_{+}=Q_{+}\tilde{\xi}_{+}+(Q_{+}(\alpha)-Q_{+})\tilde{\xi}_{+}$
and hence

$\frac{1}{2}|\tilde{\xi}_{+}|\leq|Q_{+}\tilde{\xi}_{+}|\leq
2|\tilde{\xi}_{+}|$

provided $|\alpha_{+}-y_{0}|$ is sufficiently small. Similarly, for
$|\alpha_{-}-y_{0}|$ is sufficiently small,
$\frac{1}{2}|\tilde{\xi}_{-}|\leq|(\mathbb{I}-Q_{-})\tilde{\xi}_{+}|\leq
2|\tilde{\xi}_{-}|$. In particular the map $\tilde{\xi}_{+}\mapsto
Q_{+}\tilde{\xi}_{+}$ from ${\cal R}Q_{+}(\alpha)$ into ${\cal R}Q_{+}$, and
$\tilde{\xi}_{+}\mapsto(\mathbb{I}-Q_{-})\tilde{\xi}_{+}$ from ${\cal
N}Q_{+}(\alpha)$ into ${\cal N}Q_{+}$ are linearly invertible. Then, setting

$\bar{\xi}_{+}=Q_{+}\tilde{\xi}_{+},\quad\bar{\xi}_{-}=(\mathbb{I}-Q_{-})\tilde{\xi}_{-},$

(6.8) can be written as

$\begin{array}[]{l}\bar{\xi}_{+}-\bar{\xi}_{-}=w^{-}_{0}(\alpha_{-})-w^{+}_{0}(\alpha_{+})+\bar{R}_{1}(\bar{\xi}_{+},\bar{\xi}_{-},\alpha_{+},\alpha_{-},\varepsilon)\\\
\alpha_{+}=\alpha_{-}+\bar{R}_{2}(\bar{\xi}_{+},\bar{\xi}_{-},\alpha_{+},\alpha_{-},\varepsilon)\end{array}$
(6.9)

where
$\bar{R}_{1}(\bar{\xi}_{+},\bar{\xi}_{-},\alpha_{+},\alpha_{-},\varepsilon)=O(|\bar{\xi}_{+}|^{2}+|\bar{\xi}_{-}|^{2}+|\varepsilon|)$
and
$\bar{R}_{2}(\bar{\xi}_{+},\bar{\xi}_{-},\alpha_{+},\alpha_{-},\varepsilon)=O(|\varepsilon|)$,
uniformly with respect to $(\xi_{+},\xi_{-},\alpha_{+},\alpha_{-})$. Now the
map $(\bar{\xi}_{+},\bar{\xi}_{-})\mapsto\bar{\xi}_{+}-\bar{\xi}_{-}$ is a
linear map from ${\cal R}Q_{+}\times{\cal N}Q_{-}$ into ${\cal
R}Q_{+}\times{\cal N}Q_{-}$ whose kernel is ${\cal R}Q_{+}\cap{\cal N}Q_{-}$
which, by assumption $A_{5})$, is $d$-dimensional.

Let $W\subset{\cal R}Q_{+}$ be a complement of ${\cal R}Q_{+}\cap{\cal
N}Q_{-}$ in ${\cal R}Q_{+}$, so that

${\cal R}Q_{+}+{\cal N}Q_{-}=W\oplus{\cal N}Q_{-}.$

Note that $\dim W=k-d$ and $\mathbb{R}^{n}=[{\cal R}Q_{+}+{\cal
N}Q_{-}]\oplus\hbox{\rm span}\\{\psi_{1},\ldots,\psi_{d}\\}$. Recall that we
assumed that $(\psi_{1},\ldots,\psi_{d})$ is orthonormal. Then, let
$Q:\mathbb{R}^{n}\to\mathbb{R}^{n}$ be the orthogonal projection such that
${\cal R}Q={\cal R}Q_{+}+{\cal N}Q_{-}$ and ${\cal N}Q=\hbox{\rm
span}\\{\psi_{1},\ldots,\psi_{d}\\}$. Since $(\mathbb{I}-Q)x\in{\cal
N}Q=\hbox{\rm span}\\{\psi_{1},\ldots,\psi_{d}\\}$ and
$(\psi_{1},\ldots,\psi_{d})$ is orthonormal we get

$(\mathbb{I}-Q)x=\sum_{j-1}^{d}\langle\psi_{j},(\mathbb{I}-Q)x\rangle\psi_{j}=\sum_{j-1}^{d}\langle(\mathbb{I}-Q)\psi_{j},x\rangle\psi_{j}=\sum_{j-1}^{d}(\psi_{j}^{T}x)\psi_{j}.$

Hence we replace (6.9) with

$\begin{array}[]{l}\bar{\xi}_{+}-\bar{\xi}_{-}=Q[w^{-}_{0}(\alpha_{-})-w^{+}_{0}(\alpha_{+})+\bar{R}_{1}(\bar{\xi}_{+},\bar{\xi}_{-},\alpha_{+},\alpha_{-},\varepsilon)],\\\
\alpha_{+}-\alpha_{-}=\bar{R}_{2}(\bar{\xi}_{+},\bar{\xi}_{-},\alpha_{+},\alpha_{-},\varepsilon)\\\
\psi_{j}^{T}[w^{-}_{0}(\alpha_{-})-w^{+}_{0}(\alpha_{+})+\bar{R}_{1}(\bar{\xi}_{+},\bar{\xi}_{-},\alpha_{+},\alpha_{-},\varepsilon)]=0.\end{array}$
(6.10)

Since ${\rm dim}[{\cal R}Q_{+}+{\cal N}Q_{-}]=n-d$, for any $\varepsilon$

$\begin{array}[]{l}\bar{\xi}_{+}-\bar{\xi}_{-}-Q[w^{-}_{0}(\alpha_{-})-w^{+}_{0}(\alpha_{+})]=Q\bar{R}_{1}(\bar{\xi}_{+},\bar{\xi}_{-},\alpha_{+},\alpha_{-},\varepsilon)],\\\
\alpha_{+}-\alpha_{-}=\bar{R}_{2}(\bar{\xi}_{+},\bar{\xi}_{-},\alpha_{+},\alpha_{-},\varepsilon)\end{array}$
(6.11)

is essentially a system of $n-d+m$ equations in the $n-d+2m$ variables
$(\bar{\xi}_{+},\bar{\xi}_{-},\alpha_{+},\alpha_{-})$ such that, when
$\varepsilon=0$, has the solution

$(\bar{\xi}_{+},\bar{\xi}_{-})=(0,0),\quad\alpha_{+}=\alpha_{-}=y_{0}.$

The Jacobian matrix at this point is

$J=\begin{pmatrix}L&-Qw^{-}_{0,y}(y_{0})&Qw^{+}{0,y}(y_{0})\\\
0&\mathbb{I}_{\mathbb{R}^{m}}&-\mathbb{I}_{\mathbb{R}^{m}}\end{pmatrix}$

where $L:W\times{\cal N}Q_{-}\to W\oplus{\cal N}Q_{-}$ is the invertible
linear map given by
$L(\bar{\xi}_{+},\bar{\xi}_{-})=\bar{\xi}_{+}-\bar{\xi}_{-}$. We have

$\operatorname{rank}J=\operatorname{rank}\begin{pmatrix}L&Q[w^{+}_{0,y}(y_{0})-w^{-}_{0,y}(y_{0})]&Qw^{+}{0,y}(y_{0})\\\
0&0&-\mathbb{I}_{\mathbb{R}^{m}}\end{pmatrix}=n-d+m$

hence, for $\varepsilon\neq 0$ and sufficiently small (6.11) has a
$m$-dimensional manifold of solutions

$\bar{\xi}_{+}=\bar{\xi}_{+}(\alpha_{+},\varepsilon),\quad\bar{\xi}_{-}=\bar{\xi}_{-}(\alpha_{+},\varepsilon),\quad\alpha_{-}=\alpha_{-}(\alpha_{+},\varepsilon)$

where

$\begin{array}[]{l}|\bar{\xi}_{\pm}(\alpha_{+},\varepsilon)|=O(|\alpha_{+}-y_{0}|+|\varepsilon|)\\\
|\alpha_{-}(\alpha_{+},\varepsilon)-y_{0}|=O(|\alpha_{+}-y_{0}|+|\varepsilon|)\end{array}$

Note also that

$\lim_{\varepsilon\to 0}|\alpha_{-}(\alpha_{+},\varepsilon)-\alpha_{+}|=0$

uniformly with respect to $\alpha_{+}$. Then we plug this solution in the
third equation in (6.10) and obtain the system of equations

$\psi_{j}^{T}[w^{-}_{0}(\alpha_{+})-w^{+}_{0}(\alpha_{+})+O(|\alpha_{+}-y_{0}|^{2}+|\varepsilon|])=0,\quad
j=1,\ldots,d.$

Let

${\cal
M}(\alpha_{+},\varepsilon)=\left(\psi_{j}^{*}[w^{-}_{0}(\alpha_{+})-w^{+}_{0}(\alpha_{+})+O(|\alpha_{+}-y_{0}|^{2}+|\varepsilon|)]\right)_{j=1,\ldots,d}.$

We have ${\cal M}:\mathbb{R}^{m}\times\mathbb{R}\to\mathbb{R}^{d}$, ${\cal
M}(y_{0},0)=0$ and

$\begin{array}[]{l}{\cal
M}_{\alpha_{+}}(y_{0},0)=\left(\psi_{j}^{T}[w^{-}_{0,y}(y_{0})-w^{+}_{0,y}(y_{0})]\right)_{j=1,\ldots,d}.\end{array}$

Hence from the Implicit Functions Theorem the existence follows of
$\varepsilon_{0}>0$ such that for any $|\varepsilon|<\varepsilon_{0}$ there
exists a $(m-d)$-dimensional submanifold ${\cal S}$ of $\mathbb{R}^{m}$ such
that when $\alpha_{+}\in{\cal S}$ we have ${\cal
M}(\alpha_{+},\varepsilon)=0$. For $\alpha_{+}\in{\cal S}$, we take

$\alpha_{-}=\alpha_{-}(\alpha_{+},\varepsilon),\quad\bar{\xi}_{+}=\bar{\xi}_{+}(\alpha_{+},\varepsilon),\quad\bar{\xi}_{-}=\bar{\xi}_{-}(\alpha_{+},\varepsilon),$

and $\tilde{\xi}_{\pm}=\tilde{\xi}_{\pm}(\alpha_{+},\varepsilon)$, so that

$\bar{\xi}_{+}(\alpha_{+},\varepsilon)=Q_{+}\tilde{\xi}_{+}(\alpha_{+},\varepsilon),\quad\bar{\xi}_{-}(\alpha_{+},\varepsilon)=(\mathbb{I}-Q_{-})\tilde{\xi}_{-}(\alpha_{+},\varepsilon).$

Then

$\begin{array}[]{l}x(t,\varepsilon)=x(t,\tilde{\xi}_{+}(\alpha_{+},\varepsilon),\tilde{\xi}_{-}(\alpha_{+},\varepsilon),\alpha_{+},\alpha_{-}(\alpha_{+},\varepsilon),\varepsilon)\\\
y(t,\varepsilon)=y(t,\tilde{\xi}_{+}(\alpha_{+},\varepsilon),\tilde{\xi}_{-}(\alpha_{+},\varepsilon),\alpha_{+},\alpha_{-}(\alpha_{+},\varepsilon),\varepsilon)\end{array}$

with $\alpha_{+}\in{\cal S}$ satisfies the conclusion of the Theorem. The
proof is complete. ∎

######  Remark 6.3.

i) According to Remark 4.4 we see that $(x(t,\varepsilon),y(t,\varepsilon))$
satisfies

$\begin{array}[]{l}\sup_{t\geq T_{+}}e^{\beta
t}|x(t,\varepsilon)-\tilde{w}_{+}(y(t,\varepsilon),\varepsilon)|\leq\tilde{\mu}_{1}\\\
\sup_{t\leq T_{-}}e^{-\beta
t}|x(t,\varepsilon)-\tilde{w}_{-}(y(t,\varepsilon),\varepsilon)|\leq\tilde{\mu}_{1}\\\
\end{array}$

where $\tilde{\mu}_{1},\tilde{\mu}_{2}$ do not depend on $\varepsilon$.

ii) We can replace the orthonormal basis $(\psi_{1},\ldots,\psi_{d})$ of
$[{\cal R}Q_{+}+{\cal N}Q_{-}]^{\perp}$ with any independent set
$(\tilde{\psi}_{1},\ldots,\tilde{\psi}_{d})$ such that

$\mathbb{R}^{n}=[{\cal R}Q_{+}+{\cal N}Q_{-}]\oplus\hbox{\rm
span}\\{\tilde{\psi}_{1},\ldots,\tilde{\psi}_{d}\\}.$

Indeed, let $\langle\cdot,\cdot\rangle$ be a scalar product on
$\mathbb{R}^{n}$ such that

$[{\cal R}Q_{+}+{\cal N}Q_{-}]^{\perp}=\hbox{\rm
span}\\{\tilde{\psi}_{1},\ldots,\tilde{\psi}_{d}\\}.$

and let $(\psi_{1},\ldots,\psi_{d})$ be an orthonormal basis of $\hbox{\rm
span}\\{\tilde{\psi}_{1},\ldots,\tilde{\psi}_{d}\\}$. Then an invertible
$d\times d$ matrix $M$ exists such that

$(\tilde{\psi}_{1}\ldots\tilde{\psi}_{d})=(\psi_{1}\ldots\psi_{d})M.$

Hence

$\begin{array}[]{l}[\tilde{\psi}_{j}^{T}[w^{-}_{0,y}(y_{0})-w^{+}_{0,y}(y_{0})]_{j=1,\ldots,d}=(\tilde{\psi}_{1}\ldots\tilde{\psi}_{d})^{T}[w^{-}_{0,y}(y_{0})-w^{+}_{0,y}(y_{0})]\\\
=M^{T}(\psi_{1}\ldots\psi_{d})^{T}[w^{-}_{0,y}(y_{0})-w^{+}_{0,y}(y_{0})]=M^{T}[\psi_{j}^{T}[w^{-}_{0,y}(y_{0})-w^{+}_{0,y}(y_{0})]_{j=1,\ldots,d}\end{array}$

that is $[\tilde{\psi}_{j}^{T}[w^{-}_{0,y}(y_{0})-w^{+}_{0,y}(y_{0})]$ has
rank $d$ if and only if
$[\psi_{j}^{T}[w^{-}_{0,y}(y_{0})-w^{+}_{0,y}(y_{0})]_{j=1,\ldots,d}$ has rank
$d$.

We conclude this section giving another expression for ${\cal
M}_{\alpha_{+}}(y_{0},0)$ that can be useful in the applications of Theorem
6.2.

######  Proposition 6.4.

Let $u(t,y)$ be the $C^{1}_{b}$-function defined in (2.2) and let
$\psi\in[{\cal R}Q_{+}+{\cal N}Q_{-}]^{\perp}={\rm
span}\\{\psi_{1},\ldots,\psi_{d}\\}$. Then

$\psi^{T}[w_{0,y}^{-}(y_{0})-w_{0,y}^{+}(y_{0})]=\int_{-\infty}^{\infty}\psi(t)^{T}f_{y}(u(t,y_{0}),y_{0})dt$

where

$\begin{array}[]{l}\psi(t)=\left\\{\begin{array}[]{ll}(X_{-}(t,y_{0})^{T})^{-1}Q_{-}^{T}\psi&\hbox{for
$t\leq 0$}\\\ (X_{+}^{T}(t,y_{0}))^{-1}(\mathbb{I}-Q_{+}^{T})\psi&\hbox{for
$t\geq 0$.}\end{array}\right.\end{array}$ (6.12)

Hence, the Melnikov conditon in Theorem 6.2 reads

${\rm
rank}\left[\int_{-\infty}^{\infty}\psi_{j}(t)^{T}f_{y}(u(t,y_{0}),y_{0})dt\right]_{j=1,\ldots,d}=d$
(6.13)

where $\psi_{j}(t)$ is as in (6.12) with $\psi_{j}$ instead of $\psi$.

###### Proof.

As $u_{y}(t,y_{0})$ is a bounded solution of

$\dot{x}=A(t)x+f_{y}(u(t,y_{0}),y_{0})$

where $A(t)$ is as in (3.6) with $y=y_{0}$, and

$f_{y}(u(t,y_{0}),y_{0})=\left\\{\begin{array}[]{ll}f_{-M,y}(u(t,y_{0}),y_{0})&\hbox{if
$t<t_{-M}(y_{0})$}\\\ f_{j+1,y}(u(t,y_{0}),y_{0})&\hbox{if
$t_{j}(y_{0})<t<t_{j+1}(y_{0})$}\\\ f_{i-1,y}(u(t,y_{0}),y_{0})&\hbox{if
$t_{i-1}(y_{0})<t<t_{i}(y_{0})$}\\\ f_{N,y}(u(t,y_{0}),y_{0})&\hbox{if
$t>t_{N}(y_{0})$}\\\ \end{array}\right.$

from to Lemma 3.5, equation (3.18), with $\tau=0$ we get:

$\begin{array}[]{l}\displaystyle
u_{y}(t,y_{0})=X_{-}(t,y_{0})(\mathbb{I}-Q_{-})u_{y}(0^{-},y_{0})\\\
\displaystyle+\int_{-\infty}^{t}X_{-}(t,y_{0})Q_{-}X_{-}(s,y_{0})^{-1}f_{y}(u(s,y_{0}),y_{0})ds\\\
\displaystyle-\int_{t}^{0}X_{-}(t,y_{0})(\mathbb{I}-Q_{-})X_{-}(s,y_{0})^{-1}f_{y}(u(s,y_{0}),y_{0})ds.\end{array}$

Taking the $\lim_{t\to 0^{-}}$ and using
$u_{y}(0^{-},y_{0})=w^{-}_{0,y}(y_{0})$, we get

$w^{-}_{0,y}(y_{0})=(\mathbb{I}-Q_{-})w^{-}_{0,y}(y_{0})+\int_{-\infty}^{0}Q_{-}X_{-}(s)^{-1}f_{y}(u(s,y_{0}),y_{0})ds$

that is

$Q_{-}w^{-}_{0,y}(y_{0})=\int_{-\infty}^{0}Q_{-}X_{-}(s)^{-1}f_{y}(u(s,y_{0}),y_{0})ds.$

Similarly

$(\mathbb{I}-Q_{+})w^{+}_{0,y}(y_{0})=-\int_{0}^{\infty}(\mathbb{I}-Q_{+})X_{+}(s)^{-1}f_{y}(u(s,y_{0}),y_{0})ds.$

Now, since

$[{\cal R}Q_{+}+{\cal N}Q_{-}]^{\perp}=[{\cal R}Q_{+}]^{\perp}\cap[{\cal
N}Q_{-}]^{\perp}={\cal N}Q_{+}^{T}\cap{\cal R}Q_{-}^{T}$

for any $\psi\in[{\cal R}Q_{+}+{\cal N}Q_{-}]^{\perp}$ we get

$\begin{array}[]{l}\psi^{T}(\mathbb{I}-Q_{-})=0\\\
\psi^{T}Q_{+}=0.\end{array}$

Hence, for any $\psi\in[{\cal R}Q_{+}+{\cal N}Q_{-}]^{\perp}$ we have:

$\begin{array}[]{l}\psi^{T}[w^{-}_{0,y}(y_{0})-w^{+}_{0,y}(y_{0})]=\psi_{j}^{T}[Q_{-}w^{-}_{0,y}(y_{0})-(\mathbb{I}-Q_{+})w^{+}_{0,y}(y_{0})]\\\
\qquad\displaystyle=\int_{-\infty}^{\infty}\psi(t)^{T}f_{y}(u(t,y_{0}),y_{0})dt\end{array}$
(6.14)

where $\psi(t)$ is as in (6.12) with $\psi\in[{\cal R}Q_{+}+{\cal
N}Q_{-}]^{\perp}$ instead of $\psi_{j}$. The proof is complete. ∎

The adjoint system to (5.2) and (5.3) is given by [3]

$\begin{array}[]{l}\dot{w}=-A^{T}(t,\alpha)w\quad\hbox{if $t\geq 0$}\\\
B_{\ell}(\alpha)^{T}w(t_{\ell}(\alpha)^{+})=w(t_{\ell}(\alpha)^{-})\end{array}$
(6.15)

where $\ell\in\\{-M,\ldots,N\\}\setminus\\{0\\}$.

It is easy to check that, if $\psi\in[{\cal R}Q_{+}+{\cal N}Q_{-}]^{\perp}$,
the function $\psi(t)$ defined in (6.12) is a bounded solution of (6.15) for
$\alpha=y_{0}$. We prove that if ${\rm
span}\\{\psi_{1},\ldots,\psi_{d}\\}=[{\cal R}Q_{+}+{\cal N}Q_{-}]^{\perp}$
then $\\{\psi_{1}(t),\ldots,\psi_{d}(t)\\}$ is a basis for the space of the
bounded solutions of (6.15). Indeed, the fundamental matrix of (6.15) on
$t\geq 0$ is $[X_{+}(t)^{T}]^{-1}$, and the fundamental matrix of (6.15) on
$t\leq 0$ is $[X_{-}(t)^{T}]^{-1}$. As a consequence (6.15) has an exponential
dichotomy on $\mathbb{R}_{+}$ and $\mathbb{R}_{-}$ with projections
$(\mathbb{I}-Q_{+}^{T})$ and $(\mathbb{I}-Q_{-}^{T})$ respectively. So, the
space of bounded solutions of (6.15), $C^{1}$ for $t\neq t_{\ell}(\alpha)$,
are those whose initial conditions belong to

${\cal R}(\mathbb{I}-Q_{+}^{T})\cap{\cal N}(\mathbb{I}-Q_{-}^{T})=({\cal
R}Q_{+})^{\perp}\cap({\cal N}Q_{-})^{\perp}=({\cal R}Q_{+}+{\cal
N}Q_{-})^{\perp}.$

Then the dimension of the space of solutions of (6.15), bounded on
$\mathbb{R}$, is $d$ and $(\psi_{1}(t),\ldots,\psi_{d}(t))$ span this space.

Now suppose that $x(t)$ and $\psi(t)$ are bounded solution on $\mathbb{R}$ of
(5.2)-(5.3) and (6.15) resp., both continuous for $t\neq t_{\ell}(y_{0})$. For
$t\neq t_{\ell}(y_{0})$ we have

$\frac{d}{dt}[\psi(t)^{T}x(t)]=\dot{\psi}(t)^{T}x(t)+\psi(t)^{T}\dot{x}(t)=-\psi(t)^{T}A(t,\alpha)x(t)+\psi(t)^{T}\dot{x}(t)=0.$

Moreover

$\begin{array}[]{l}\psi(t_{\ell}(\alpha)^{+})^{T}x(t_{\ell}(\alpha)^{+}))=[(B_{\ell}(\alpha)^{-1})^{T}\psi(t_{\ell}(\alpha)^{-})]^{T}x(t_{\ell}(\alpha)^{+})\\\
=\psi(t_{\ell}(\alpha)^{-})^{T}[B_{\ell}(\alpha)^{-1}x(t_{\ell}(\alpha)^{+}]=\psi(t_{\ell}(\alpha)^{-})^{T}x(t_{\ell}(\alpha)^{-})\end{array}$

Thus we conclude that $\psi(t)^{T}x(t)$ is constant on $\mathbb{R}$ (see also
[21]).

##  7 An example

An interesting application of Theorem 6.2 is when $d=1$ that is when

${\cal R}Q_{+}\cap{\cal N}Q_{-}=\hbox{\rm span}\\{\dot{u}(0,y_{0})\\}.$

This condition is trivially satisfied when $n=2$ since in this case $k=n-k=1$.
Moreover, when $n=2$, we also have $\dim{\cal R}Q_{+}=\dim{\cal N}Q_{-}=1$ and
hence

${\cal R}Q_{+}={\cal N}Q_{-}=\hbox{\rm span}\\{\dot{u}(0,y_{0})\\}.$ (7.1)

In this section we consider examples of applications of Theorem 6.2 with
$n=2$, $m=1$ and $d=1$. First we prove some general facts concerning two-
dimensional differential equations depending on a slowly varying variable. So
the system is

$\begin{array}[]{l}\dot{x}_{1}=F_{1}(x_{1},x_{2},y)\\\
\dot{x}_{2}=F_{2}(x_{1},x_{2},y)\\\ \dot{y}=\varepsilon
g(x_{1},x_{2},y).\end{array}$ (7.2)

Suppose $u(t,y)=(u_{1}(t,y),u_{2}(t,y))$ is a piecewise smooth solution of
(7.2) for $\varepsilon=0$ satisfying assumptions $A_{1})-A_{4})$. To write the
Melnikov condition (6.13), that in this case reads

$\int_{-\infty}^{\infty}\psi(t)^{T}\begin{pmatrix}F_{1,y}(u(t,y_{0}),y_{0})\\\
F_{2,y}(u(t,y_{0}),y_{0})\end{pmatrix}dt\neq 0,$

we need to know the (unique) bounded solution $\psi(t)$ of the adjoint system
(6.15).

Let

$J=\begin{pmatrix}0&-1\\\ 1&0\end{pmatrix}$

and, again, $\ell\in\\{-M,\ldots,-1,1,\ldots,N\\}$.

We prove the following

######  Proposition 7.1.

Let $A(t)=[a_{jk}(t)]_{1\leq j,k\leq
2}:=[F_{j,x_{k}}(u_{1}(t),u_{2}(t),y_{0})]_{1\leq j,k\leq 2}$, $B_{\ell}$ as
in (5.1), with $\alpha=y_{0}$ and

$v(t):=e^{-\int_{0}^{t}a_{11}(s)+a_{22}(s)ds}J\dot{u}(t,y_{0})=e^{-\int_{0}^{t}a_{11}(s)+a_{22}(s)ds}\begin{pmatrix}-\dot{u}_{2}(t,y_{0})\\\
\dot{u}_{1}(t,y_{0})\end{pmatrix}.$

Then the space of bounded solution of the adjoint variational system are of
the form

$\psi(t)=\left\\{\begin{array}[]{ll}\mu_{-M}v(t)&\hbox{for $t\leq t_{-M}$}\\\
\mu_{j+1}v(t),&\hbox{for $t_{j}<t\leq t_{j+1}$}\\\ \mu_{i}v(t),&\hbox{for
$t_{i}\leq t<t_{i+1}$}\\\ \mu_{N}v(t)&\hbox{for $t\geq
t_{N}$}\end{array}\right.$

where $\mu_{-M}\neq 0$ is arbitrary and

$\begin{array}[]{l}\mu_{j+1}v(t_{j}^{+})=\mu_{j}[B_{j}^{T}]^{-1}v(t_{j}^{-})\\\
\mu_{i}v(t_{i}^{+})=\mu_{i-1}[B_{i}^{T}]^{-1}v(t_{i}^{-}).\end{array}$ (7.3)

###### Proof.

First we show that constants satisfying (7.3) exist. Indeed recall that for
any $i=1,\ldots,N$

$\begin{array}[]{l}[[B_{i}^{T}]^{-1}v(t_{i}^{-})]^{T}\dot{u}(t_{i}^{+},y_{0})=v(t_{i}^{-})^{T}[B_{i}^{-1}\dot{u}(t_{i}^{+},y_{0})]=v(t_{i}^{-})^{T}\dot{u}(t_{i}^{-},y_{0})=0\\\
v(t_{i}^{+})^{T}\dot{u}(t_{i}^{+},y_{0})=0,\end{array}$

where $t_{i}=t_{i}(y_{0})$, and then

$v(t_{i}^{+})=\tilde{\mu}_{i}[B_{i}^{T}]^{-1}v(t_{i}^{-})$

for some $\tilde{\mu}_{i}\in\mathbb{R}$, since both vectors are orthogonal to
$\dot{u}(t_{i}^{+},y_{0})$. In a similar way we see that constants satisfying
(7.3) for $j=-1,\ldots,-M$ exist.

Next, using $JA(t)={\rm det}A(t)[A(t)^{T}]^{-1}J$:

$\begin{array}[]{l}\dot{v}(t)=e^{-\int_{0}^{t}a_{11}(s)+a_{22}(s)ds}JA(t)\begin{pmatrix}\dot{u}_{1}(t,y_{0})\\\
\dot{u}_{2}(t,y_{0})\end{pmatrix}-(a_{11}(t)+a_{22}(t))v(t)\\\
=e^{-\int_{0}^{t}a_{11}(s)+a_{22}(s)ds}{\rm
det}A(t)[A(t)^{T}]^{-1}J\begin{pmatrix}\dot{u}_{1}(t,y_{0})\\\
\dot{u}_{2}(t,y_{0})\end{pmatrix}-(a_{11}(t)+a_{22}(t))v(t)\\\ ={\rm
det}A(t)[A(t)^{T}]^{-1}v(t)-(a_{11}(t)+a_{22}(t))v(t)=-A(t)^{T}v(t).\end{array}$

Hence

$\dot{\psi}(t)=-A(t)^{T}\psi(t),$

for any $t\neq t_{\ell}$. Moreover

$\begin{array}[]{l}\psi(t_{j}^{-})=\mu_{j}v(t_{j}^{-})=B_{j}^{T}\mu_{j+1}v(t_{j}^{+})=B_{j}^{T}\psi(t_{j}^{+})\\\
\psi(t_{i}^{+})=\mu_{i}v(t_{i}^{+})=\mu_{i-1}[B_{i}^{T}]^{-1}v(t_{i}^{-})=[B_{i}^{T}]^{-1}\psi(t_{i}^{-}).\end{array}$

Finally, we prove that $\psi(t)$ is bounded. As the adjoint system has an
exponential dichotomy on $\mathbb{R}_{+}$, resp. $\mathbb{R}_{-}$, with
projections $\mathbb{I}-Q_{+}^{T}$, resp. $I-Q_{-}^{T}$, it is enough to prove
that

$\psi(0)\in{\cal R}(\mathbb{I}-Q_{+}^{T})\cap{\cal N}(I-Q_{-}^{T})={\cal
R}(Q_{+})^{\perp}\cap{\cal N}(Q_{-})^{\perp}=\\{\dot{u}(0)\\}^{\perp}$

the last equality following from (7.1). But

$\langle\psi(0),\dot{u}(0,y_{0})\rangle=\langle
J\dot{u}(0,y_{0}),\dot{u}(0,y_{0})\rangle=0$

and then $\psi(t)$ is bounded concluding the proof of the Proposition. ∎

######  Remark 7.2.

i) From (7.3) we have

$\begin{array}[]{l}\mu_{i}J\dot{u}(t_{i}^{+},y_{0})=\mu_{i-1}[B_{i}^{T}]^{-1}J\dot{u}(t_{i}^{-},y_{0})\\\
\mu_{j+1}J\dot{u}(t_{j}^{+},y_{0})=\mu_{j}[B_{j}^{T}]^{-1}J\dot{u}(t_{j}^{-},y_{0})\end{array}$

and then

$\mu_{i}\|\dot{u}(t_{i}^{+},y_{0})\|^{2}=\mu_{i}\langle
J\dot{u}(t_{i}^{+},y_{0}),J\dot{u}(t_{i}^{+},y_{0})\rangle=\mu_{i-1}\langle[B_{i}^{T}]^{-1}J\dot{u}(t_{i}^{-},y_{0}),J\dot{u}(t_{i}^{+},y_{0})\rangle$

and similarly

$\mu_{j+1}\|\dot{u}(t_{j}^{+},y_{0})\|^{2}=\mu_{j+1}\langle
J\dot{u}(t_{j}^{+},y_{0}),J\dot{u}(t_{j}^{+},y_{0})\rangle=\mu_{j}\langle[B_{j}^{T}]^{-1}J\dot{u}(t_{j}^{-},y_{0}),J\dot{u}(t_{j}^{+},y_{0})\rangle.$

Hence all $\mu_{\ell}$’s can be computed in terms of
$\dot{u}(t_{\ell}^{\pm},y_{0})$.

ii) Since $\mu_{-M}\neq 0$ and all $B_{j}$, $B_{i}$ are invertible, we see
that $\mu_{\ell}\neq 0$ for all $\ell$.

The case where all $\mu_{\ell}$ are equal is of particular interest, since in
this case we can take $\psi(t)=v(t)$ and the Melnikov condition reads

$\Delta:=\int_{-\infty}^{\infty}e^{-\int_{0}^{t}a_{11}(s)+a_{22}(s)ds}\begin{pmatrix}-\dot{u}_{2}(t,y_{0})\\\
\dot{u}_{1}(t,y_{0})\end{pmatrix}^{T}\begin{pmatrix}F_{1,y}(u(t,y_{0}),y_{0})\\\
F_{2,y}(u(t,y_{0}),y_{0})\end{pmatrix}dt\neq 0.$

If, moreover, $a_{11}(t)+a_{22}(t)=0$ we have

$\begin{array}[]{l}\displaystyle\Delta=\int_{-\infty}^{\infty}F_{2,y}(u(t,y_{0}),y_{0})\dot{u}_{1}(t,y_{0})-F_{1,y}(u(t,y_{0}),y_{0})\dot{u}_{2}(t,y_{0})dt\\\
\displaystyle=\int_{\Gamma}F_{2,y}(u_{1},u_{2},y_{0})du_{1}-F_{1,y}(u_{1},u_{2},y_{0})du_{2}\end{array}$
(7.4)

where $\Gamma=\\{(u_{1}(t),u_{2}(t))|t\in\mathbb{R}\\}$.

We have the following

######  Proposition 7.3.

Equations (7.3) are satisfied with $\mu_{0}=1$ and $\mu_{\ell}=1$, if and only
if there exist $\nu_{\ell}$ such that

$J[\dot{u}(t_{\ell}(y_{0})^{+},y_{0})-\dot{u}(t_{\ell}(y_{0})^{-},y_{0})]=\nu_{\ell}h_{x}(u(t_{\ell}(y_{0}),y_{0}),y_{0})^{T}.$
(7.5)

###### Proof.

We have $\mu_{0}=1$, $\mu_{\ell}=1$ for all $\ell$, if and only if the
following holds.

$B_{\ell}^{T}J\dot{u}(t_{\ell}(y_{0})^{+},y_{0})=J\dot{u}(t_{\ell}(y_{0})^{-},y_{0})$
(7.6)

We check that

$\begin{array}[]{l}\langle x,B_{\ell}y\rangle=\langle
x-\frac{\langle\dot{u}(t_{\ell}(y_{0})^{-},y_{0})-\dot{u}(t_{\ell}(y_{0})^{+},y_{0}),x\rangle}{h_{x}(u(t_{\ell},y_{0}),y_{0})\dot{u}(t_{\ell}(y_{0})^{-},y_{0})}h_{x}(u(t_{\ell}(y_{0}),y_{0}),y_{0})^{T},y\rangle\end{array}$

so that

$B_{\ell}^{T}x=x-\frac{\langle\dot{u}(t_{\ell}(y_{0})^{-},y_{0})-\dot{u}(t_{\ell}(y_{0})^{+},y_{0}),x\rangle}{h_{x}(u(t_{\ell},y_{0}),y_{0})\dot{u}(t_{\ell}(y_{0})^{-},y_{0})}h_{x}(u(t_{\ell}(y_{0}),y_{0}),y_{0})^{T}$

Then (7.6) is equivalent to :

$\begin{array}[]{l}J\dot{u}(t_{\ell}^{-},y_{0})=J\dot{u}(t_{\ell}^{+},y_{0})\\\
-\frac{\langle\dot{u}(t_{\ell}(y_{0})^{-},y_{0})-\dot{u}(t_{\ell}(y_{0})^{+},y_{0}),J\dot{u}(t_{\ell}(y_{0})^{+},y_{0})\rangle}{h_{x}(u(t_{\ell},y_{0}),y_{0})\dot{u}(t_{\ell}(y_{0})^{-},y_{0})}h_{x}(u(t_{\ell}(y_{0}),y_{0}),y_{0})^{T}\end{array}$
(7.7)

or else, as $\langle Jx,x\rangle=0$,

$J\dot{u}(t_{\ell}^{+},y_{0})-J\dot{u}(t_{\ell}^{-},y_{0})=\frac{\langle\dot{u}(t_{\ell}(y_{0})^{-},y_{0}),J\dot{u}(t_{\ell}(y_{0})^{+},y_{0})\rangle}{h_{x}(u(t_{\ell}(y_{0}),y_{0})\dot{u}(t_{\ell}(y_{0})^{-},y_{0})}h_{x}(u(t_{\ell}(y_{0}),y_{0})^{T}$
(7.8)

which is (7.5) with

$\nu_{\ell}=\frac{\langle\dot{u}(t_{\ell}(y_{0})^{-},y_{0}),J\dot{u}(t_{\ell}(y_{0})^{+},y_{0})\rangle}{h_{x}(u(t_{\ell}(y_{0}),y_{0})\dot{u}(t_{\ell}(y_{0})^{-},y_{0})}.$

On the other hand, if (7.5) holds, taking the scalar product with
$\dot{u}(t_{\ell}(y_{0}),y_{0})$ we get

$\begin{array}[]{l}\nu_{\ell}\langle
h_{x}(u(t_{\ell}(y_{0}),y_{0})^{T},\dot{u}(t_{\ell}(y_{0})^{-},y_{0})\rangle=\langle
J[\dot{u}(t_{\ell}(y_{0})^{+},y_{0})-\dot{u}(t_{\ell}(y_{0})^{-},y_{0})],\dot{u}(t_{\ell}(y_{0})^{-},y_{0})\rangle\end{array}$

that is

$\nu_{\ell}=\frac{\langle
J\dot{u}(t_{\ell}(y_{0})^{+},y_{0}),\dot{u}(t_{\ell}(y_{0})^{-},y_{0})\rangle}{h_{x}(u(t_{\ell}(y_{0}),y_{0})\dot{u}(t_{\ell}(y_{0})^{-},y_{0})}$

and then (7.6) follows, given the equivalence between (7.6) and (7.8) ∎

######  Remark 7.4.

As $J[\dot{u}(t_{\ell}(y_{0})+,y_{0})-\dot{u}(t_{\ell}(y_{0})^{-},y_{0})]$ is
orthogonal to
$\dot{u}(t_{\ell}(y_{0})^{+},y_{0})-\dot{u}(t_{\ell}(y_{0})^{-},y_{0})$ and
$h_{x}(u(t_{\ell}(y_{0}),y_{0})^{T}$ is orthogonal to the tangent space to
${\cal S}_{\ell}(y_{0})$ at $u(t_{\ell}(y_{0}),y_{0})$, say
$T_{u(t_{\ell}(y_{0}),y_{0})}{\cal S}_{\ell}(y_{0})$, condition (7.5) is
equivalent to the fact that
$\dot{u}(t_{\ell}(y_{0})^{+},y_{0})-\dot{u}(t_{\ell}(y_{0})^{-},y_{0})$
belongs to $T_{(u(t_{\ell}(y_{0}),y_{0},y_{0})}{\cal S}_{\ell}(y_{0})$.

For example, suppose

$h(x,y)=x_{k}$

where either $k=1$ or $k=2$. Recalling that

$f_{\ell}(x,y)=\begin{pmatrix}F_{1,\ell}(x_{1},x_{2},y)\\\
F_{2,\ell}(x_{1},x_{2},y)\end{pmatrix}$

we get, omitting the argument $u(t_{\ell}(y_{0}),y_{0})$ for simplicity:

$J[\dot{u}(t_{\ell}^{+}(y_{0}),y_{0})-\dot{u}(t_{\ell}^{-}(y_{0}),y_{0})]=\left\\{\begin{array}[]{ll}\begin{pmatrix}F_{2,i-1}-F_{2,i}\\\
F_{1,i}-F_{1,i-1}\end{pmatrix}&\hbox{if $\ell=i$}\\\
\begin{pmatrix}F_{2,j}-F_{2,j+1}\\\ F_{1,j+1}-F_{1,j}\end{pmatrix}&\hbox{if
$\ell=j$}\end{array}\right.$

and then (7.5) holds if and only if

$\begin{array}[]{l}F_{k,i}(u(t_{\ell}),y_{0})=F_{k,i-1}(u(t_{\ell}),y_{0}),\quad
i=1,\ldots,N-1\\\
F_{k,j}(u(t_{\ell}),y_{0})=F_{k,j+1}(u(t_{\ell}),y_{0}),\quad
j=-1,\ldots,-M\end{array}$

To give a specific example, consider the second order, discontinuous equation
with slowly varying coefficients:

$\begin{array}[]{ll}\ddot{u}+u(u-a_{-}(\varepsilon t+\bar{y}))(1-u)=0&\hbox{if
$u<c$}\\\ \ddot{u}+u(u-a_{+}(\varepsilon t+\bar{y}))(1-u)=0&\hbox{if $u>c$}\\\
\end{array}$ (7.9)

where
$0<\inf\\{a_{\pm}(y),y\in\mathbb{R}\\}<\sup\\{a_{\pm}(y),y\in\mathbb{R}\\}<1$.

We prove the following

######  Proposition 7.5.

Let $a_{\pm}(y)$ be $C^{1}$-functions such that

$0<a_{+}(y)<\sup\\{a_{+}(y):y\in\mathbb{R}\\}\leq\frac{1}{2}\leq\inf\\{a_{-}(y):y\in\mathbb{R}\\}<a_{-}(y)<1$

and let $0<c<1$ be a fixed number. Then equation (7.9) with $\varepsilon=0$
has a family of $C^{1}$-solutions $u_{\pm}(t,y)$ defined for $t\leq 0$, $t\geq
0$ resp., bounded together with their derivatives, such that
$u_{-}(0,y)=u_{+}(0,y)=c$ and

$\lim_{t\to-\infty}u_{-}(t,y)=0,\quad\lim_{t\to\infty}u_{+}(t,y)=1$

uniformly with respect to $y$. Next, let
$D(y)=c^{2}(2c-3)[a_{+}(y)-a_{-}(y)]+a_{+}(y)-\frac{1}{2}$ and suppose that
$y_{0}$ and $0<c<1$ exist such that

$\begin{array}[]{l}D(y_{0})=0,\quad D^{\prime}(y_{0})\neq 0.\end{array}$
(7.10)

Then $\dot{u}_{-}(0,y_{0})=\dot{u}_{+}(0,y_{0})$ and there exists
$\varepsilon_{0}>0$ such that for $0<\varepsilon<\varepsilon_{0}$ there exists
a $C^{1}$-function $y_{0}(\varepsilon)$, such that $\lim_{\varepsilon\to
0}y_{0}(\varepsilon)=y_{0}$, and a $C^{1}$-solution $u(t,\varepsilon)$ of
equation (7.9) with $\bar{y}=y_{0}(\varepsilon)$, bounded with its derivatives
and such that

$\begin{array}[]{l}\displaystyle\lim_{\varepsilon\to
0}\sup_{t\in\mathbb{R}_{\pm}}|u(t,\varepsilon)-u_{\pm}(t,\varepsilon
t+y_{0}(\varepsilon))|=0\\\ \displaystyle\lim_{\varepsilon\to
0}\sup_{t\in\mathbb{R}_{\pm}}|\dot{u}(t,\varepsilon)-\dot{u}_{\pm}(t,\varepsilon
t+y_{0}(\varepsilon))|=0\end{array}$

###### Proof.

Writing $x=\begin{pmatrix}u\\\ \dot{u}\end{pmatrix}$, $y=\varepsilon
t+\bar{y}$, equation (7.9) reads

$\dot{x}=f(x,y)=\begin{pmatrix}x_{2}\\\
x_{1}(x_{1}-a(x_{1},y))(x_{1}-1)\end{pmatrix},\quad\dot{y}=\varepsilon$ (7.11)

where

$a(x_{1},y)=\left\\{\begin{array}[]{ll}a_{-}(y)&\hbox{if $x_{1}<c$}\\\
a_{+}(y)&\hbox{if $x_{1}>c$.}\end{array}\right.$

Then $h(x,y)=x_{1}$ and

$f_{\pm}(x,y)=\begin{pmatrix}x_{2}\\\
x_{1}(x_{1}-a_{\pm}(y))(x_{1}-1)\end{pmatrix}.$

Note that

$f_{+}(x,y)-f_{-}(x,y)=\begin{pmatrix}0\\\
(a_{-}(y)-a_{+}(y))x_{1}(x_{1}-1)\end{pmatrix}$

is tangent to the manifold $h(x,y)=c$ and (7.9) has the fixed points
$w_{-}(y)=(0,\;0)^{T}$ and $w_{+}(y)=(1,\;0)^{T}$.

Let (7.11)0 be equation (7.11) with $\varepsilon=0$ and
$(u_{\pm}(t,y),\dot{u}_{\pm}(t,y))$, $t\in\mathbb{R}_{\pm}$, be solutions of
(7.11)0 such that

$\begin{array}[]{l}\displaystyle\lim_{t\to-\infty}(u_{-}(t,y),\dot{u}_{-}(t,y))=(0,\;0),\\\
\displaystyle\lim_{t\to\infty}(u_{+}(t,y),\dot{u}_{+}(t,y))=(1,\;0).\end{array}$

Multiplying (7.9) (with $\varepsilon=0$) by $\dot{u}_{-}(t,y)$ and integrating
from $-\infty$ to $t<0$ we get

$\begin{array}[]{l}\displaystyle\frac{1}{2}\dot{u}_{-}(t,y)^{2}=\int_{0}^{u_{-}(t,y)}u(u-a_{-})(u-1)du\\\
\displaystyle\qquad=\frac{u_{-}(t,y)^{2}(3u_{-}(t,y)^{2}-4(a_{-}(y)+1)u_{-}(t,y)+6a_{-}(y)}{12}\end{array}$

hence, if $u_{-}(t,y)>0$, it satisfies

$\frac{du}{dt}=\frac{u\sqrt{3u^{2}-4(a_{-}(y)+1)u+6a_{-}(y)}}{\sqrt{6}}.$
(7.12)

It is easy to check that, for $1<2a_{-}(y)<2$ we have
$3u^{2}-4(a_{-}(y)+1)u+6a_{-}(y)>0$.

To simplify notation in the following we write $a_{-}$ for $a_{-}(y)$.
Integrating (7.12) we get (see [13, eq. 2.266])

$\begin{array}[]{l}{\rm
sinh}^{-1}\left(\frac{\sqrt{2}[3a_{-}-(a_{-}+1)u_{-}(t,y)]}{u_{-}(t,y)\sqrt{(2-a_{-})(2a_{-}-1)}}\right)={\rm
sinh}^{-1}\left(\frac{\sqrt{2}[3a_{-}-(a_{-}+1)c]}{c\sqrt{(2-a_{-})(2a_{-}-1)}}\right)-t\sqrt{a_{-}}.\end{array}$
(7.13)

where $c=u_{-}(0,y)$. Then, using ${\rm sinh}(\alpha-\beta)={\rm
sinh}\alpha{\rm cosh}\beta-{\rm sinh}\beta{\rm cosh}\alpha$ and ${\rm
cosh}\alpha=\sqrt{1+{\rm sinh}^{2}\alpha}$:

$\begin{array}[]{l}\frac{3a_{-}\sqrt{2}}{u_{-}(t,y)\sqrt{(2-a_{-})(2a_{-}-1)}}={\rm
sinh}\left({\rm
sinh}^{-1}\left(\frac{\sqrt{2}[3a_{-}-(a_{-}+1)c]}{c\sqrt{(2-a_{-})(2a_{-}-1)}}\right)-t\sqrt{a_{-}}\right)+\frac{(a_{-}+1)\sqrt{2}}{\sqrt{(2-a_{-})(2a_{-}-1)}}\\\
={\rm sinh}\left({\rm
sinh}^{-1}\left(\frac{\sqrt{2}[3a_{-}-(a_{-}+1)c]}{c\sqrt{(2-a_{-})(2a_{-}-1)}}\right)\right){\rm
cosh}(t\sqrt{a_{-}})\\\ \qquad-{\rm sinh}(t\sqrt{a_{-}}){\rm cosh}\left({\rm
sinh}^{-1}\left(\frac{\sqrt{2}[3a_{-}-(a_{-}+1)c]}{c\sqrt{(2-a_{-})(2a_{-}-1)}}\right)\right)+\frac{(a_{-}+1)\sqrt{2}}{\sqrt{(2-a_{-})(2a_{-}-1)}}\\\
=\frac{\sqrt{2}[3a_{-}-(a_{-}+1)c]}{c\sqrt{(2-a_{-})(2a_{-}-1)}}{\rm
cosh}(t\sqrt{a_{-}})-\sqrt{1+\left(\frac{\sqrt{2}[3a_{-}-(a_{-}+1)c]}{c\sqrt{(2-a_{-})(2a_{-}-1)}}\right)^{2}}{\rm
sinh}(t\sqrt{a_{-}})+\frac{(a_{-}+1)\sqrt{2}}{\sqrt{(2-a_{-})(2a_{-}-1)}}\\\
=\frac{\sqrt{2}[3a_{-}-(a_{-}+1)c]}{c\sqrt{(2-a_{-})(2a_{-}-1)}}{\rm
cosh}(t\sqrt{a_{-}})-\sqrt{\frac{3a_{-}(3c^{2}-4(a_{-}+1)c+6a_{-})}{c^{2}(2-a_{-})(2a_{-}-1)}}{\rm
sinh}(t\sqrt{a_{-}})+\frac{(a_{-}+1)\sqrt{2}}{\sqrt{(2-a_{-})(2a_{-}-1)}}\end{array}$

Hence

$\frac{1}{u_{-}(t,y)}=\mu_{1}(a_{-}){\rm
cosh}(t\sqrt{a_{-}})-\mu_{2}(a_{-}){\rm
sinh}(t\sqrt{a_{-}})+\frac{a_{-}+1}{3a_{-}}.$ (7.14)

where

$\begin{array}[]{l}\mu_{1}(a)=\frac{3a-(a+1)c}{3ac}=\frac{1}{c}-\frac{1}{3}\left(1+\frac{1}{a}\right).\\\
\mu_{2}(a)=\sqrt{\frac{3c^{2}-4(a+1)c+6a}{6ac^{2}}}=\sqrt{\frac{1}{a}\left(\frac{1}{2}+\frac{2}{3c}\right)+\frac{1}{c^{2}}-\frac{2}{3c}}\end{array}$
(7.15)

We pause for a while to observe that $\mu_{2}(a)$ is well defined if
$3c^{2}-4(a+1)c+6a\geq 0$. This holds for sure if $\frac{1}{2}\leq a<1$.
However, if $0<a<\frac{1}{2}$ we need that

$3c<2(a+1)-\sqrt{4a^{2}-10a+4}.$

We will consider this issue in Section 8.

As $a_{-}(y)>0$ and its derivative are bounded on $\mathbb{R}$, so are
$\mu_{1}(a_{-}(y))$ and $\mu_{2}(a_{-}(y))$, Moreover, from (7.14) it follows
that $u_{-}(t,y)$ is $C^{1}$ in $(t,y)$.

We have

$\begin{array}[]{l}\frac{1}{|u_{-}(t,y)|}\geq|\mu_{1}(a_{-}){\rm
cosh}(t\sqrt{a_{-}})-\mu_{2}(a_{-}){\rm sinh}(t\sqrt{a_{-}})|\\\
=\frac{1}{2}|(\mu_{1}(a_{-})-\mu_{2}(a_{-}))e^{t\sqrt{a_{-}}}+(\mu_{1}(a_{-})+\mu_{2}(a_{-}))e^{-t\sqrt{a_{-}}})|\\\
\geq\frac{1}{2}|(\mu_{1}(a_{-})+\mu_{2}(a_{-}))e^{-t\sqrt{a_{-}}})|-\frac{1}{2}|(\mu_{1}(a_{-})-\mu_{2}(a_{-}))e^{t\sqrt{a_{-}}}|\end{array}$

As

$\begin{array}[]{l}(\mu_{2}(a)+\mu_{1}(a))(\mu_{2}(a)-\mu_{1}(a)=\mu_{2}(a)^{2}-\mu_{1}(a)^{2}=\frac{(2a-1)(2-a)}{18a^{2}}>0\end{array}$
(7.16)

for $1<2a\leq 2$, we see that $\mu_{2}(a)+\mu_{1}(a)$ and
$\mu_{2}(a)-\mu_{1}(a)$ do not change sign in $1<2a\leq 2$. But, since for
$a=1$ we have

$\mu_{2}(1)=\frac{1}{c}-\frac{2}{3}=\frac{3-c}{3c}>0,\quad\mu_{1}(1)=\frac{1}{c}\sqrt{\frac{3c^{2}-8c+6}{6}}>0,$

we see that $\inf_{1<2a\leq 2}[\mu_{2}(a)+\mu_{1}(a)]>0$. Moreover
$\mu_{1}(\frac{1}{2})=\frac{1}{c}-1$ and
$\mu_{2}(\frac{1}{2})=\sqrt{\frac{c^{2}-2c+1}{c^{2}}}=\frac{1-c}{c}>0$ and
that $\mu_{1}(\frac{1}{2})+\mu_{2}(\frac{1}{2})=2\frac{1-c}{c}>0$. So

$\inf_{1\leq 2a\leq 2}[\mu_{2}(a)+\mu_{1}(a)]>0$

and then

$\frac{1}{u_{-}(t,y)}\geq\frac{1}{2}\inf_{y\in\mathbb{R}}\\{\mu_{2}(a_{-})+\mu_{1}(a_{-})\\}e^{-t\inf_{y\in\mathbb{R}}a_{-}}.$

As a consequence $u_{-}(t,y)>0$ for $t\leq 0$ and

$\lim_{t\to-\infty}u_{-}(t,y)=0$

uniformly with respect to $y\in\mathbb{R}$. From (7.12) we also see that
$\dot{u}_{-}(t,y)>0$ and $\displaystyle\lim_{t\to-\infty}\dot{u}_{-}(t,y)=0$
uniformly with respect to $y$. In particular, for $t\leq 0$,
$u_{-}(t,y),\dot{u}_{-}(t,y)$ are bounded, positive, functions and
$u_{-}(t,y)$ is strictly increasing from $0$ to $u_{-}(0,y)=c$.

Now, we look for a (strictly) increasing solution of the second equation in
(7.9) on $t\geq 0$. To this end we observe that, $u_{+}(t,y)$ is a (strictly)
increasing solution, on $t\geq 0$, of the second equation in (7.9) such that

$u_{+}(0,y)=c,\quad\lim_{t\to\infty}u_{+}(t,y)=1$

if and only if $v(t)=1-u_{+}(-t)$ is a (strictly) increasing solution, on
$t\leq 0$, of

$\ddot{v}+v(v-(1-a_{+}))(1-v)=0$

such that

$v(0)=1-c,\quad\lim_{t\to-\infty}v(t)=0.$

Since $0<a_{+}(y)\leq\frac{1}{2}$ is equivalent to $\frac{1}{2}\leq
1-a_{+}(y)<1$, from the previous part we conclude that the second equation in
(7.9) has a unique strictly increasing solution $u_{+}(t,y)$ such that

$u_{+}(0,y)=c,\quad\lim_{t\to\infty}u_{+}(t,y)=1$

the limit being uniform with respect to $y\in\mathbb{R}$.

Then for any $y\in\mathbb{R}$, equation (7.9) with $\varepsilon=0$ has a pair
of solutions $(u_{-}(t,y),\dot{u}_{-}(t,y))$, $(u_{+}(t,y),\dot{u}_{+}(t,y))$
defined for $t\leq 0$ and $t\geq 0$ resp., such that $u_{\pm}(t)$ are
increasing in their interval of definition and

$\begin{array}[]{l}0<u_{-}(t,y)<c,\;\hbox{for $t<0$ and $u_{-}(0,y)=c$}\\\
c<u_{+}(t,y)<1,\;\hbox{for $t>0$ and $u_{+}(0,y)=c$}\\\
\displaystyle\lim_{t\to-\infty}u_{-}(t,y)=\lim_{t\to\infty}1-u_{+}(t,y)=0\\\
\displaystyle\lim_{t\to-\infty}\dot{u}_{-}(t,y)=\lim_{t\to\infty}\dot{u}_{+}(t,y)=0\end{array}$

uniformly with respect to $y$.

Now, since
$\displaystyle\lim_{t\to-\infty}(u_{-}(t,y),\dot{u}_{-}(t,y))=(0,0)$ and
$\displaystyle\lim_{t\to\infty}(u_{+}(t,y),\dot{u}_{+}(t,y))=(1,0)$ we get:

$\begin{array}[]{l}\displaystyle\frac{1}{2}\dot{u}_{-}(0,y)^{2}=\int_{0}^{c}u(u-a_{-}(y))(u-1)du\\\
\displaystyle\frac{1}{2}\dot{u}_{+}(0,y)^{2}=\int_{c}^{1}u(u-a_{+}(y))(1-u)du.\end{array}$
(7.17)

So equation (7.9), with $\varepsilon=0$, has a $C^{1}$ solution heteroclinic
to the fixed points $w_{-}(y)$ and $w_{+}(y)$ if and only if

$\int_{0}^{c}u(u-a_{-}(y))(u-1)du+\int_{c}^{1}u(u-a_{+}(y))(u-1)du=0.$

It is easy to check that

$\int_{0}^{c}u(u-a_{-}(y))(u-1)du+\int_{c}^{1}u(u-a_{+}(y))(u-1)du=\frac{1}{6}D(y).$
(7.18)

and hence equation (7.9) has a $C^{1}$ solution heteroclinic to the fixed
points $w_{-}(y)$ and $w_{+}(y)$ if and only if

$D(y)=0.$ (7.19)

From (7.10) we see that equation (7.19) has a unique solution only for
$y=y_{0}$.

Then equation (7.11)0 has a unique $C^{1}$ solution $u_{0}(t)=u(t,y_{0})$
asymptotic to $u=0$ as $t\to-\infty$ and to $u=1$ as $t\to\infty$ only for
$y=y_{0}$ and this solution breaks when $y\neq y_{0}$ (in the sense that it is
no longer $C^{1}$). Since $f_{+}(x,y)-f_{-}(x,y)$ is tangent to the manifold
$h(x,y)=c$ and

$f_{y}(u_{0}(t),\dot{u}_{0}(t),y)=\left\\{\begin{array}[]{ll}\begin{pmatrix}0\\\
-a_{-}^{\prime}(y_{0})u_{0}(t)(u_{0}(t)-1)\end{pmatrix}&\hbox{if $t\leq 0$}\\\
\\\ \begin{pmatrix}0\\\
-a_{+}^{\prime}(y_{0})u_{0}(t)(u_{0}(t)-1)\end{pmatrix}&\hbox{if
$t>0$}\end{array}\right.$

the Melnikov function reads:

$\begin{array}[]{l}\displaystyle-\int_{-\infty}^{0}a_{-}^{\prime}(y_{0})\dot{u}_{0}(t)u_{0}(t)(u_{0}(t)-1)dt-\int_{0}^{\infty}a_{+}^{\prime}(y_{0})\dot{u}_{0}(t)u_{0}(t)(u_{0}(t)-1)dt\\\
=\displaystyle-a_{-}^{\prime}(y_{0})\int_{0}^{c}u(u-1)du-
a_{+}^{\prime}(y_{0})\int_{c}^{1}u(u-1)du=\frac{1}{6}D^{\prime}(y_{0})\neq
0.\end{array}$

From Theorem 6.2 the existence follows of $0<\rho<1$ and $\varepsilon_{0}>0$
such that for $0<\varepsilon<\varepsilon_{0}$ there exist a $C^{r}$-function
$y_{0}(\varepsilon)$ with $\lim_{\varepsilon\to 0}y_{0}(\varepsilon)=y_{0}$
and a solution $\tilde{u}(t,\varepsilon)$ of (7.11) such that

$\begin{array}[]{l}\sup_{t\in\mathbb{R}}|\tilde{u}(t,\varepsilon)-u(t,\varepsilon
t+y_{0}(\varepsilon))|+\sup_{t\in\mathbb{R}}|\dot{\tilde{u}}(t,\varepsilon)-\dot{u}(t,\varepsilon
t+y_{0}(\varepsilon))|<\rho,\\\ \lim_{\varepsilon\to
0}\\{\sup_{t\in\mathbb{R}}|\tilde{u}(t,\varepsilon)-u(t,\varepsilon
t+y_{0}(\varepsilon))|+\sup_{t\in\mathbb{R}}|\dot{\tilde{u}}(t,\varepsilon)-\dot{u}(t,\varepsilon
t+y_{0}(\varepsilon))|\\}=0.\end{array}$ (7.20)

∎

######  Remark 7.6.

Suppose that $a_{-}(y)-a_{+}(y)=\frac{1}{2}$ with $0<a_{+}(y)<\frac{1}{2}$.
Then (7.10) reads:

$\begin{array}[]{l}2a_{+}(y_{0})=2c^{3}-3c^{2}+1\\\ a^{\prime}_{+}(y_{0})\neq
0.\end{array}$

Note that $0<2c^{3}-3c^{2}+1<1$ for $0<c<1$.

##  8 Concluding remark

The assumption $0<a_{+}(y)\leq\frac{1}{2}\leq a_{-}(y)<1$ can be slightly
weakened. Indeed, suppose that

$0<a_{\min}=\min\\{a_{-}(y):y\in\mathbb{R}\\}\leq a_{-}(y)<1$

where $a_{\min}<\frac{1}{2}$.

By uniqueness of analytical continuation, the function defined in (7.14) is a
solution of the first equation in (7.5) (with $\varepsilon=0$) for any value
of $a_{-}=a_{-}(y)$ for which $\mu_{2}(a_{-})$ has a meaning, that is such
that

$3c^{2}-4(a_{-}+1)c+6a_{-}\geq 0$ (8.1)

Moreover to prove that $u_{-}(t,y)\to 0$ as $t\to-\infty$, uniformly with
respect to $y$, following the same argument of Proposition 7.5, we need that
$\mu_{1}(a)+\mu_{2}(a)>0$ that is

$3a-(a+1)c+3a\sqrt{\frac{3c^{2}-4(a+1)c+6a}{6a}}>0.$ (8.2)

We prove that (8.1) and (8.2) hold if and only if

$c<\frac{1}{3}\left(2(a_{\min}+1)-\sqrt{4a_{\min}^{2}-10a_{\min}+4}\right).$
(8.3)

As the function of $a$: $3c^{2}-4(a+1)c+6a$ is increasing, for $2c<3$, (8.1)
holds if and only if

$3c^{2}-4(a_{\min}+1)c+6a_{\min}\geq 0$

that is

$\hbox{either}\;c\leq\frac{2(a_{\min}+1)-\sqrt{2\Delta}}{3}\;\hbox{ or
}\;c\geq\frac{2(a_{\min}+1)+\sqrt{2\Delta}}{3}$

where $\Delta=2a_{\min}^{2}-5a_{\min}+2$. However it is easy to check that,
for $a_{\min}<\frac{1}{2}$, $\frac{2(a_{\min}+1)+\sqrt{2\Delta}}{3}>1$. Hence
(8.1) is equivalent to

$c\leq\frac{2(a_{\min}+1)-\sqrt{2\Delta}}{3}$ (8.4)

Next, the function of $a$

$\mu(a):=(a+1)c-\sqrt{\frac{3a}{2}}\sqrt{3c^{2}-4(a+1)c+6a}-3a$

is convex and its values at $a=a_{\min}<\frac{1}{2}$ and $a=\frac{1}{2}$ are

$\begin{array}[]{l}\mu(a_{\min})=(a_{\min}+1)c-\sqrt{\frac{3a_{\min}}{2}}\sqrt{3c^{2}-4(a_{\min}+1)c+6a_{\min}}-3a_{\min}\\\
\mu(\frac{1}{2})=\frac{3}{2}c-\sqrt{\frac{3}{4}}\sqrt{3c^{2}-6c+3}=\frac{3}{2}(c-|c-1|)=\frac{3}{2}(2c-1)-3=3c-\frac{9}{2}<0\end{array}$

(since $c<1$). Condition $\mu(a_{\min})<0$ for $a_{\min}\leq\frac{1}{2}$ is
equivalent to: