Research

# Homoclinic and heteroclinic solutions for a class of second-order non-autonomous ordinary differential equations: multiplicity results for stepwise potentials

Elisa Ellero and Fabio Zanolin*

Author Affiliations

Department of Mathematics and Computer Science, University of Udine, Via delle Scienze 206, Udine, 33100, Italy

For all author emails, please log on.

Boundary Value Problems 2013, 2013:167  doi:10.1186/1687-2770-2013-167

 Received: 30 November 2012 Accepted: 11 June 2013 Published: 15 July 2013

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

### Abstract

We prove some multiplicity results for a class of one-dimensional nonlinear Schrödinger-type equations of the form

x + 2 k x 2 g ( t ) x 3 = 0 ,

where k > 0 and the weight g ( t ) is a positive stepwise function. Instead of the cubic term, more general nonlinearities can be considered as well.

MSC: 34C37, 34B15.

##### Keywords:
homoclinic solutions; heteroclinic solutions; multiplicity; nonlinear Schrödinger equation; stepwise potential; topological methods

### Dedication

Dedicated to Professor Jean Mawhin

### 1 Introduction

This paper deals with the study of homoclinic and heteroclinic solutions for a class of nonlinear second-order differential equations of the form

x + σ x a ( t ) h ( x ) = 0 , (1.1)

where σ > 0 is a fixed (positive) coefficient and a ( ) : R R 0 + : = ] 0 , + ) is a bounded weight function. For the nonlinear term, which we split as

h ( x ) = x f ( x ) ,

we suppose that

(∗) f : R R + : = [ 0 , + ) is a locally Lipschitz function which is even, strictly increasing on [ 0 , + ) such that f ( 0 ) = 0 and f ( + ) = + .

Examples of functions f which are suitable for our considerations are, for instance, f ( s ) = s 2 as in [1], or f ( s ) = s 2 + b s 4 ( b > 0 ) as in [2]. In both these cases, the function f is even. In Section 3 we briefly describe how this restriction could be avoided.

Equations of the form (1.1) naturally arise in the search of particular solutions for some classes of nonlinear Schrödinger equations (NLSE) with inhomogeneous nonlinearities. Typical examples are related to the NLSE

i ħ Ψ t = ħ 2 2 m Δ Ψ + V 1 ( x ) Ψ + V 2 ( x ) | Ψ | p 1 Ψ , x R N , (1.2)

where ħ is the Planck constant, m is the particle’s mass, i is the imaginary unit, Δ is the Laplace operator and p > 1 . In literature, solutions of the form

Ψ ( t , x ) = u ( x ) exp ( i λ ħ 1 t ) ,

where λ R and u ( x ) is a real-valued function, are called stationary waves. Their search leads to the equation

ħ 2 2 m Δ u + λ u + V 1 ( x ) u V 2 ( x ) | u | p 1 u = 0 , x R N (1.3)

(see [3]). Usually, the Planck constant and the mass are omitted in (1.3) after rescaling.

In many significant models of NLSE, one-dimensional waves are considered. They are studied, for example, in nonlinear optics, in the theory of ocean rogue waves and for Bose-Einstein condensates (just to mention a few cases). For instance, in [1], Belmonte-Beitia and Torres analyzed the one-dimensional equation

1 2 u + λ u + g ( x ) u 3 = 0 , x R , (1.4)

which is a particular case of (1.3) with N = 1 , p = 3 and V 1 = 0 . In equation (1.4) the nonlinearity is called inhomogeneous due to the presence of a non-constant weight V 2 ( x ) = g ( x ) , which in [1] is assumed to be positive. The sign condition on g ( x ) implies that nontrivial bounded solutions of (1.4) can exist only for λ < 0 (see [[1], Theorem 1]). For this reason, we prefer to set λ = k with k > 0 .

In order to study equation (1.4) or its variants, we are going to follow a dynamical system approach, hence we choose to treat the independent variable (which in the applications has a spatial connotation) as a time variable, via the substitution x t . Similarly, for the dependent variable, we make the substitution u x . In this way, equation (1.4) reads as

x + 2 k x 2 g ( t ) x 3 = 0 , (1.5)

which belongs to the same class of (1.1).

In spite of the apparent simplicity of equation (1.5), a throughout study of its solutions may be a rather difficult task for a general nonconstant weight function g ( t ) . Looking for homoclinic and heteroclinic solutions of (1.5), it will be natural to focus on the behavior of g ( t ) at ±∞. In similar situations, various authors have confined their study to the case in which g ( t ) is asymptotically constant [4-6], or eventually constant [7,8]. These assumptions are also justified by the analysis of some physical underlying models, in which a layered structure is present. With this respect, see the introduction in [8], where different eventually autonomous cases are listed for related NLSEs arising in nonlinear optics. Examples in which the nonlinear term presents a piecewise constant weight function have been studied also in biological and chemical models. In particular, these situations occur in the theory of wave propagation for reaction-diffusion systems; see, for instance, [9-11]. In the context of equation (1.5), examples in which g ( t ) is a piecewise constant function have been considered as well (see [12]).

The aim of the present paper is to provide multiplicity results regarding homoclinic and heteroclinic solutions for equation (1.1) under particular assumptions on the weight function a ( t ) . Actually, we suppose that a ( t ) is an eventually constant piecewise function with only two steps.

Our approach combines phase-plane analysis with time-mapping estimates. As in [5,7,8,10,11,13], the solutions are obtained by connecting the unstable and stable manifolds of the equilibrium points of the asymptotically autonomous equations. Such connections are performed by means of orbit paths of an intermediate equation, which represents the behavior of the system during a suitable interval of transition [ t , t + ] between the asymptotic states. Multiple connecting solutions arise when such interval length is sufficiently large. Lower estimates for t + t will be provided in terms of time mappings, which can be expressed by Abelian integrals. This kind of approach is also reminiscent of some topological methods for the study of Sturm-Liouville boundary value problems. Indeed, a solution that satisfies the Sturm-Liouville boundary conditions can be interpreted as a trajectory in the phase-plane that connects two lines (see [14-17]). Generalized Sturm-Liouville solutions, which connect the graphs of two functions or given planar continua, have been considered as well (see [18-20]). Often these problems can be settled in the framework of the theory of ODEs with nonlinear boundary conditions (see [21,22]).

In order to make our approach more transparent, we are going to perform our analysis for equation (1.5). This choice is motivated by the sake of avoiding unnecessary technicalities. Our arguments can be modified in a straightforward manner in the case of more general equation (1.1), with f satisfying (∗) (see Section 3). Homoclinic and heteroclinic orbits can be interpreted as solutions for some boundary value problems on unbounded intervals. In the last section we also outline possible applications of our approach to boundary value problems on a compact interval (like the Sturm-Liouville one).

Besides his manifold achievements in different areas of mathematics, Professor Jean Mawhin is one of the pioneers in the study of topological methods for nonlinear boundary value problems. It is a pleasure and an honor to have the possibility to dedicate our work to his important contributions in this area.

### 2 Homoclinic and heteroclinic solutions: multiplicity results

#### 2.1 General setting

We consider the second-order nonlinear differential equation

x + 2 k x 2 g ( t ) x 3 = 0 (2.1)

and its equivalent system in the phase-plane

{ x = y , y = 2 k x + 2 g ( t ) x 3 , (2.2)

where k > 0 is a fixed coefficient and g : R R 0 + : = ] 0 , + ) is a bounded measurable function. Solutions of (2.1) are meant in the generalized (i.e., Carathéodory) sense (see [23]). Actually, in our results a step function g ( t ) with only two jumps is considered and the solutions are piecewise smooth.

In the particular case of a constant coefficient g ( t ) μ > 0 , an elementary analysis of the system

{ x = y , y = 2 k x + 2 μ x 3 (2.3)

shows that the associated phase portrait presents a center O = ( 0 , 0 ) and two saddle points P ± : = ( ± ( k / μ ) 1 2 , 0 ) . These saddles are connected by two heteroclinic solutions (one connecting P to P + in the upper half-plane and a symmetric one from P + to P in the lower half-plane). The heteroclinic orbits are described by

y 2 = μ ( x 2 ( k / μ ) ) 2 for  ( k / μ ) 1 2 < x < ( k / μ ) 1 2 . (2.4)

In [1] Belmonte-Beitia and Torres supposed that the weight function g ( t ) is even andT-periodic and proved the existence of a heteroclinic solution connecting two periodic stationary states. Such a result generalizes to the case of periodic coefficients the situation described above for the autonomous system.

As already observed in the introduction, equation (2.1) has been already studied by various authors for its relevance in many applicative models. The aim of the present paper is to provide multiplicity results for heteroclinic and also homoclinic solutions of equation (2.1). A possible way to obtain this goal is to assume that the weight function g ( t ) has a different behavior at infinity and in some intermediate time interval. In this setting, a first step consists of analyzing the situation in which g ( t ) is a stepwise function assuming only two values. We believe that the study of such a simplified case may lead to more general considerations, in which, for instance, g ( t ) is eventually periodic. In this way, multiplicity results for Belmonte-Beitia and Torres’s model can be obtained (see Section 3 for a brief overview concerning possible applications of our approach). Equation (2.1) with a periodic stepwise coefficient g ( t ) has recently been studied in [24] in the context of chaotic-like dynamics.

As a tool, in our paper, we use the Poincaré map associated to system (2.2). Given a fixed time interval [ t 0 , t 1 ] and a point z R 2 , we indicate by ζ ( , t 0 , z ) the solution ζ ( t ) = ( x ( t ) , y ( t ) ) of (2.2) satisfying the initial condition ζ ( t 0 ) = z . The Poincaré map on [ t 0 , t 1 ] is defined as

Φ = Φ t 0 t 1 : z ζ ( t 1 , t 0 , z ) .

From the fundamental theory of differential equations it follows that the domain D Φ of Φ is an open subset of the plane and Φ is an orientation preserving homeomorphism of D Φ onto its image Φ ( D Φ ) .

#### 2.2 Analysis of the equation

Let us consider system (2.2) for the stepwise weight function g ( t ) , defined as

g ( t ) = { μ 1 for  t t  or  t t + , μ 0 for  t < t < t + , (2.5)

where

t < t + and 0 < μ 0 < μ 1 .

System (2.2) can be seen as the superposition of the autonomous systems

{ x = y , y = 2 k x + 2 μ 1 x 3 (2.6)

and

{ x = y , y = 2 k x + 2 μ 0 x 3 . (2.7)

The first one describes the asymptotic behavior of the equation. We are going to use the orbits of the second system for connecting unstable/stable manifolds of (2.6) during the time interval [ t , t + ] . Figure 1 shows the superposition of the phase portraits of the two systems.

Figure 1. In the present figure we have considered the superposition of the phase portraits of systems (2.6) (in darker color) and (2.7) (in lighter color) with k = 1 , μ 1 = 4 , μ 0 = 1 / 4 . For graphical reasons, a slightly different x- and y-scaling has been used.

First of all, we briefly analyze the structure of system (2.6). Its equilibrium points are the origin (which is a center) and two saddle points P ± = ( ± ξ 1 , 0 ) , where

ξ 1 : = ( k / μ 1 ) 1 2 .

The stable and unstable manifolds of P ± are illustrated in Figure 2. The sets I ( P ) L ( P + , P ) and O ( P ) L ( P , P + ) are, respectively, the stable and the unstable manifolds of  P . Symmetrically, I ( P + ) L ( P , P + ) and O ( P + ) L ( P + , P ) are the stable and the unstable manifolds of P + .

Figure 2. The stable and unstable manifolds of P and P + . The arrows show the direction of the flow along the orbits.

In the case of the cubic nonlinearity, we find

I ( P ) = { ( x , y ) : y = | μ 1 x 2 k | / μ 1 , x < ( k / μ 1 ) 1 2 } , O ( P ) = { ( x , y ) : y = | μ 1 x 2 k | / μ 1 , x < ( k / μ 1 ) 1 2 } , I ( P + ) = { ( x , y ) : y = | μ 1 x 2 k | / μ 1 , x > ( k / μ 1 ) 1 2 } , O ( P + ) = { ( x , y ) : y = | μ 1 x 2 k | / μ 1 , x > ( k / μ 1 ) 1 2 } , L ( P , P + ) = { ( x , y ) : y = | μ 1 x 2 k | / μ 1 , ( k / μ 1 ) 1 2 < x < ( k / μ 1 ) 1 2 } , L ( P + , P ) = { ( x , y ) : y = | μ 1 x 2 k | / μ 1 , ( k / μ 1 ) 1 2 < x < ( k / μ 1 ) 1 2 } .

The solutions of system (2.2) with g ( t ) defined as in (2.5) will be obtained by connecting suitably chosen parts of stable and unstable manifolds of (2.6) with trajectories of (2.7). More in detail, we proceed as follows.

Due to the Hamiltonian nature of the equation under consideration, we can consider the ‘energy’ level lines of system (2.7), given by

E ( x , y ) : = y 2 + 2 k x 2 μ 0 x 4 = c . (2.8)

We also set

Γ c : = { ( x , y ) : ξ 0 x ξ 0 , E ( x , y ) = c } ,

where

ξ 0 : = ( k / μ 0 ) 1 2 .

The curve Γ c is the part of the line at level c R in the strip [ ξ 0 , ξ 0 ] × R . For every c such that

0 < c < 0 : = k 2 μ 0 ,

the set Γ c is a periodic orbit of (2.7) run in the clockwise sense, which intersects the y-axis at the points ( 0 , ± c ) . The (minimal) period τ c of Γ c can be expressed by the following time-mapping formula:

τ c = 2 s 0 s 1 d x c 2 F ( x ) , (2.9)

where F ( x ) is the potential associated to the equation and s 0 , s 1 are such that

F ( s 0 ) = F ( s 1 ) = c / 2 with  ξ 0 < s 0 < 0 < s 1 < ξ 0 .

In our case,

F ( x ) = k x 2 μ 0 2 x 4

and by symmetry s 0 = s 1 . More in detail, we can express the period τ c by means of an elliptic integral in the following way:

τ c = 4 μ 0 0 s 1 d x ( a 2 x 2 ) ( b 2 x 2 ) , (2.10)

for a : = s 1 , b : = ( 2 k 2 μ 0 2 s 1 2 ) 1 / 2 and 0 < s 1 < ξ 0 , 0 < c = 2 k s 1 2 μ 0 s 1 4 < 0 . Notice that 0 < a < b in (2.10). Therefore

τ c = 4 b μ 0 sn 1 ( 1 , a / b ) ,

where sn ( , k ) is the Jacobi elliptic sine function of modulus k : = a / b (see [25,26]).

The level line Γ 0 contains the saddle points Q ± : = ( ± ξ 0 , 0 ) of system (2.7) and their heteroclinic connections L ( Q , Q + ) in the upper half-plane and L ( Q + , Q ) in the lower half-plane. The energy level

1 : = k 2 μ 1

corresponds to the closed curve Γ 1 , which is tangent to L ( P , P + ) and L ( P + , P ) . Moreover, for

: = k 2 μ 1 ( 2 μ 0 μ 1 ) ,

the energy level line Γ contains the saddle points P ± of system (2.6).

Using the parameters 0 , 1 and , we define the regions

E : = { ( x , y ) : ξ 0 x ξ 0 , 1 E ( x , y ) } , F : = { ( x , y ) : ξ 0 x ξ 0 , < E ( x , y ) < 0 } .

Both ℰ and ℱ are invariant sets for (2.7); they are filled by periodic orbits. In the sequel, the orbits in the region ℰ will be called internal, while those in ℱ will be called external. The choice of these names is made in order to distinguish the trajectories with respect to  Γ , which contains the saddle points P ± of system (2.6). (See Figure 3.)

Figure 3. The energy level lines Γ 0 , Γ 1 and Γ for system (2.7) superimposed over the stable and unstable manifolds of P ± for system (2.7). The colored regions are the invariant sets ℰ (the darker one) and ℱ (lighter). The figure has been plotted with the parameters k = 1 , μ 0 = 1 / 2 , μ 1 = 2 . For graphical reasons, a slightly different x- and y-scaling has been used.

In order to obtain homoclinic or heteroclinic solutions for system (2.2) with g ( t ) defined as in (2.5), we connect the unstable manifold of the point A { P , P + } to the stable manifold of B { P , P + } via an orbit path of Γ c for c [ 1 , 0 [ . Actually, this is the only way to get the desired solutions. If c , such connection will follow a trajectory ( x ( t ) , y ( t ) ) with ξ 0 x ( t ) ξ 0 , so it will lie in the invariant annular region ℰ bounded by Γ 1 and Γ . On the contrary, for c > the connection lies in the set ℱ between Γ and Γ 0 .

To provide more explicit details, we consider the case A = P + . A solution of (2.2), which is homoclinic to P + , can be produced in two ways. One consists in connecting O ( P + ) to I ( P + ) , using an orbit path of Γ c for c ] , 0 [ . Another possibility is given by connecting a point of L ( P + , P ) to a point of L ( P , P + ) via an orbit path of Γ c for c [ 1 , [ . For c = , we can obtain our solution if and only if t + t is an integer multiple of the fundamental period τ of the orbit Γ . Indeed, such a solution will be constantly equal to P + for t ] t , t + [ , and for t [ t , t + ] it coincides with the periodic solution of (2.7) such that x ( t ) = ξ 1 and y ( t ) = 0 . This solution makes ( t + t ) / τ turns around the origin for t [ t , t + ] .

Similar considerations can be developed with respect to heteroclinic solutions. For instance, a heteroclinic orbit from P + to P can be obtained as follows: by connecting a point of O ( P + ) with one of I ( P ) using an orbit path of Γ c for c ] , 0 [ , or by a connection of two different points of L ( P + , P ) via an orbit path of Γ c for c [ 1 , [ . For c = , we can obtain our solution if and only if t + t is an odd multiple of τ / 2 , where τ is the period of the orbit Γ . Indeed, such a solution will be constantly equal to P + for t ( , t ] and constantly equal to P for t [ t + , + ) . Moreover, it coincides for t [ t , t + ] with the solution of (2.7) such that x ( t ) = ξ 1 and y ( t ) = 0 . This solution makes 2 ( t + t ) / τ half-turns around the origin for t [ t , t + ] .

Figure 4, although not exhaustive of all the conceivable cases, summarizes several different possibilities. Indeed, a trajectory homoclinic to P + can be obtained as follows: move by system (2.6) along O ( P + ) from P + to the intersection point of O ( P + ) with a closed curve external to P + , namely δ (for t ( , t ] ). This point is put in evidence with a small black circle. Next, follow δ (by system (2.7)) in order to reach the intersection of δ with I ( P + ) (for t [ t , t + ] ). Such an intersection point is indicated by a grey square. Finally, switch to system (2.6) and move toward P + along I ( P + ) (for t [ t + , + ) ). With a similar procedure, we can obtain a heteroclinic trajectory from P + to P . In fact, we can move by system (2.6) along O ( P + ) from P + to the intersection point of O ( P + ) with a closed curve external to P + , namely γ (for t ( , t ] ). Next, follow γ (by system (2.7)) to the intersection with I ( P ) (for t [ t , t + ] ). Finally, switch to system (2.6) and move toward P along I ( P ) (for t [ t + , + ) ).

Figure 4. In the present figure we have considered the superposition of (2.6) and (2.7) for k = 1 , μ 1 = 4 , μ 0 = 1 / 4 . The closed curves α, β, γ, δ represent four different level lines of system (2.7): the curves α and β in ℰ correspond to a level c ] 1 , [ , while γ and δ in ℱ refer to c ] , 0 [ . In order to describe a connection from an unstable manifold to a stable one of system (2.6) via an orbit path of (2.7), we have marked with a black circle possible starting points and with a grey square some available end points on the same level lines.

In general, the connections through the trajectories of system (2.7) will either involve only an arc of the closed curve Γ c , or they will require to perform a certain number of winds around the origin, depending on the available time. If the length of the time interval [ t , t + ] is small, only the first case is possible. However, if t + t is large enough, more choices arise. For instance, looking at the orbit δ in Figure 4, we could make a certain number of loops from the black circle before reaching the grey square on the same orbit. Similar considerations apply to the heteroclinic trajectory described above with reference to γ.

Until now we have described the ‘external connections’, namely those which lie in the region ℱ. Further possibilities appear if we consider ‘internal connections’ by means of orbit paths contained in ℰ. For example, if we look at the orbit β in Figure 4, we can obtain trajectories homoclinic to P + by choosing one of the two black circles as a starting point from L ( P + , P ) and one of the two grey squares as an end point on L ( P , P + ) . Loops along β will be permitted if the time is sufficiently large. The construction of heteroclinic orbits from P + to P via the line α follows a similar procedure.

#### 2.3 Study of the Poincaré map

Another possible point of view to describe the previous construction of homoclinic and heteroclinic solutions consists in considering the Poincaré map associated to system (2.7) for the time interval [ t , t + ] . Due to the autonomous nature of this system, we have that Φ t t + = Φ 0 t + t , denoted simply by Φ when no confusion may occur. Moreover, the region

W : = { ( x , y ) : ξ 0 x ξ 0 , E ( x , y ) 0 }

is a compact invariant set contained in the domain of Φ. Let us denote by Y : = O ( P + ) ¯ F ¯ the part of the unstable manifold O ( P + ) between P + and Γ 0 , including the extreme points. Hence the external connections can be precisely described by looking at the intersection points of Φ ( Y ) with I ( P + ) (for the homoclinics) and with I ( P ) (for the heteroclinics), as illustrated in Figure 5.

Figure 5. In the present figure we show the transformation ofunder the Poincaré mapΦfor t + t = 5.5 . The parameters used are k = 1 , μ 1 = 2 , μ 0 = 1 / 2 . For graphical reasons a slightly different x- and y-scaling has been used.

With reference to Figure 5, we observe that the intersection point of O ( P + ) with Γ 0 (indicated by a small circle) is moved by Φ along Γ 0 to a point very close to Q + . On the other hand, during the time interval [ t , t + ] , the point P + makes a little more than a complete turn around the origin; the final position Φ ( P + ) is indicated by a cross. The arc Φ ( Y ) is a spiral-like curve with one end near Q + and the other on Γ . For the time length t + t = 5.5 , considered in Figure 5, there are precisely two intersections of Φ ( Y ) with I ( P + ) , corresponding to external homoclinic solutions. Meanwhile, we have also an external heteroclinic solution due to the intersection of Φ ( Y ) with I ( P ) . All these three intersection points have been indicated by a grey square (following the same convention used in Figure 4).

In order to describe further these solutions, let us consider their energy at the time t + that we denote by k 1 , k 2 , k 3 , with < k 3 < k 2 < k 1 < 0 . Let us define the points H i as { H i } = O ( P + ) Γ k i for i = 1 , 2 , 3 . The solution ζ ( t ) = ( x ( t ) , y ( t ) ) of (2.2) with ζ ( t ) = H 1 or ζ ( t ) = H 3 has the following behavior: ζ ( t ) P + for t ± , with x ( t ) > ξ 1 and convex for t ] t , t + [ . Moreover, x ( t ) is increasing on ( , t ] and decreasing on [ t + , + ) . On the interval [ t , t + ] , if i = 1 , the solution x ( t ) is concave, while if i = 3 , the solution has two maxima and one (negative) minimum separated by two simple zeros. The solution of system (2.2) with ζ ( t ) = H 2 is such that ζ ( t ) P + for t with x ( t ) > ξ 1 and ζ ( t ) P for t + with x ( t ) < ξ 1 ; x ( t ) is convex for t ] t , t + [ . Moreover, the solution has one maximum and one minimum separated by one simple zero.

Even if the time interval length t + t is small, we always have at least one intersection of Φ ( Y ) with I ( P + ) and thus a homoclinic solution. On the other hand, if t + t grows, the spiral curve Φ ( Y ) will wind more times around the origin, hence more homoclinic/heteroclinic solutions will appear.

If we look for the internal connections (made by orbit paths lying in the region ℰ), we proceed as follows.

Let us set L : = L ( P + , P ) ¯ and denote by L : = L ( P , P + ) ¯ its symmetric part with respect to the x-axis. Figure 6 illustrates the involved geometry from the point of view of the Poincaré map. We observe that all the points of ℒ are contained in ℰ, hence they are periodic points (of different periods). As before, during the time interval [ t , t + ] , the point P + moves clockwise around the origin along Γ ; the final position Φ ( P + ) Γ is indicated by a cross. Symmetrically, the point P performs exactly the same angle around the origin, reaching the final position Φ ( P ) Γ , indicated by a small circle. All the other points of L ( P + , P ) move, under the action of Φ, on the energy level lines Γ c with 1 c < . The arc Φ ( L ) is a curve contained in ℰ, connecting the two antipodal points Φ ( P ± ) Γ and leaning on one point of Γ 1 . This tangent point is the image through Φ of the point R : = ( 0 , 1 ) . Such a property of Φ ( L ) implies that Φ ( L ) intersects both ℒ and L (at different points). These two points have been indicated by a grey square in the figure. Accordingly, we always find at least an internal homoclinic solution and an internal heteroclinic one.

Figure 6. In the present figure we show the transformation of L : = L ( P + , P ) ¯ under the Poincaré mapΦfor t + t = 5.5 , the same time considered in Figure5. The parameters used are k = 1 , μ 1 = 2 , μ 0 = 1 / 2 , as before. The point indicated by a cross is precisely the same point Φ ( P + ) of Figure 5. Notice that for this time interval the point makes a little more than a complete turn around the origin. For graphical reasons, a slightly different x- and y-scaling has been used.

Summarizing the above information, we conclude that, for every time interval [ t , t + ] , we have at least two solutions which are homoclinic to P + (one external and one internal) and one internal heteroclinic solution from P + to P .

For the external connections, if the time interval length t + t grows, the situation can be summarized as follows. The number of external homoclinic/heteroclinic solutions increases, depending on the number of winds around the origin of the curve Φ ( Y ) . Indeed, the end point of such a curve on Γ 0 cannot go beyond Q + , while the other end point Φ ( P + ) is free to move on the periodic orbit Γ . If we denote by η ext and ν ext the number of external homoclinic and heteroclinic solutions respectively, it holds that 0 η ext ν ext 1 . From the above discussion, we can conclude that if, for some nonnegative integer n, we have t + t > n τ (where τ is the period of Γ ), then η ext n + 1 . An analogous lower bound can be provided for ν ext . For a formal proof, see Theorem 2.1.

When t + t grows, the situation for the internal homoclinic/heteroclinic solutions is more intriguing, as illustrated in Figure 7 and Figure 8. A twist effect depending on the different periods of Γ 1 and Γ produces a double spiral-like curve Φ ( L ) when t + t is sufficiently large.

Figure 7. In the present figure we show the transformation ofunder the Poincaré mapΦfor t + t = 50 . The parameters used are k = 1 , μ 1 = 2 , μ 0 = 1 / 2 , as before. For graphical reasons, a slightly different x- and y-scaling has been used.

Figure 8. The present figure shows the evolution ofunder the Poincaré mapΦfor time intervals of different lengths. The parameters used are k = 1 , μ 1 = 2 , μ 0 = 1 / 2 , as before, while for t + t , we have considered the cases 35, 60, 100, 150 (from the left to the right). For graphical reasons, a slightly different x- and y-scaling has been used.

The fact that Φ ( L ) looks like a double spiral depends on the different velocities of the points P ± and R = ( 0 , 1 ) (which is the intersection of ℒ with the negative y-axis). In fact, since τ 1 < τ , the points on Γ 1 move faster than those of Γ . By construction, P ± Γ , while R Γ 1 . Therefore, when the time gap t + t is sufficiently large, the number of turns of Φ ( R ) around the origin will exceed the number of turns of Φ ( P ± ) . Accordingly, the image through Φ of the right part of ℒ (the sub-arc of ℒ connecting P + to R) is a spiral that winds a certain number of times in the clockwise sense around the origin. Indeed, it connects the slower point Φ ( P + ) to the faster one Φ ( R ) inside the region ℰ, as illustrated in Figure 9. The number of half-turns of this spiral depends on a ‘rotational gap’ that, for our purposes, it will be convenient to define as

w ( t + t ) : = 2 ( t + t ) τ 1 2 ( t + t ) τ 2 . (2.11)

Similarly, the image through Φ of the left part of ℒ (the sub-arc of ℒ connecting R to  P ) is a spiral-like curve winding a certain number of times in the counterclockwise sense around the origin. Indeed, it connects the points Φ ( R ) and Φ ( P ) inside the region ℰ (the first point moves at a faster speed than the second one). Again, the number of half-turns of this second curve will depend on w ( t + t ) . In the end, gluing together the two spiral-like curves, we conclude that, if t + t is sufficiently large, the arc Φ ( L ) will appear like a double spiral with a central ‘hook’. To better describe the structure of Φ ( L ) , we should observe that the ‘tip of the hook’ is not Φ ( R ) , but it is the image through Φ of a point (close to R) on the left part of L { R } (we owe this remark to the referee). For a formal proof based on the argument outlined above, see Theorem 2.2.

Figure 9. The present figure explains the twist effect due to the different periods of Γ 1 and Γ . We have considered the evolution of the part of ℒ with x 0 connecting P + to the point R for two different time gaps: t + t = 5.5 and t + t = 50 . The corresponding two images Φ ( P + ) are indicated by a cross, while the images Φ ( R ) are marked by a black dot. The parameters used are k = 1 , μ 1 = 2 , μ 0 = 1 / 2 , as before. Notice that the two images Φ ( P + ) , although close to each other in the figure, have made a different number of turns around the origin. For graphical reasons, a slightly different x- and y-scaling has been used.

In preparation for this theorem, we introduce the following notation. Given a positive real number S and an energy level [ 1 , ] , we denote by n ̲ ( S , ) : = 2 S τ the lower integer part of S / ( τ / 2 ) , and by n ¯ ( S , ) : = 2 S τ the upper integer part of S / ( τ / 2 ) . By definition, during a time interval of length S, a point on the orbit Γ makes more that n ¯ ( S , ) 1 half-turns around the origin, but less that n ̲ ( S , ) + 1 half-turns. Thus the rotational gap defined above can be written as w ( T ) = n ¯ ( T , 1 ) n ̲ ( T , ) 2 for T = t + t . By definition, in the open interval ] 2 T / τ , 2 T / τ 1 [ , the smallest integer is n ̲ ( T , ) + 1 , while the largest integer is n ¯ ( T , ) 1 . Hence, the open interval ] 2 T / τ , 2 T / τ 1 [ contains w ( T ) + 1 positive integers.

#### 2.4 Conclusion

After the preliminary study of the previous subsections, we are now in a position to express our results for the equation

as statements with a formal proof. Recall that the weight function g ( t ) is defined as follows:

g ( t ) = { μ 1 for  t t  or  t t + , μ 0 for  t < t < t + ,

where

t < t + and 0 < μ 0 < μ 1 .

We also set T : = t + t and consider the Poincaré map Φ 0 T associated to the system

In order to make the exposition more clear, we have decided to consider separately the cases of external and internal connections.

Theorem 2.1 (External connections)

Under the above assumptions, the following results hold for solutions which satisfy the energy condition

x ( t ) 2 + 2 k x ( t ) 2 μ 0 x ( t ) 4 ] , 0 [ , t R :

there always exist a (positive) solution homoclinic to P + and a (negative) solution homoclinic to P ;

if t + t > τ / 2 , there always exist a heteroclinic solution from P + to P and a heteroclinic solution from P to P + , both with exactly one zero;

if t + t > n τ (for some integer n 1 ) then, for each integerjwith 1 j n , there exist at least one solution homoclinic to P + and one solution homoclinic to P , both with exactly 2jzeros;

if t + t > ( n + 1 ) τ / 2 (for some integer n 1 ) then, for each integerjwith 1 j n , there exist at least one heteroclinic solution from P + to P and one from P to P + , both with exactly 2 j + 1 zeros.

Proof As a first step, we focus our attention on the search of solutions homoclinic to the point P + . The case of homoclinics to P is analogous, thus it will be omitted.

Recalling that T = t + t , our goal is to find points on O ( P + ) which are moved by the Poincaré map Φ 0 T on I ( P + ) .

We introduce in the phase-plane for (2.7), a system of polar coordinates ( θ , ρ ) , with center in the origin. The initial points in the arc (which is the closure of the intersection of O ( P + ) with the external region ℱ) are parameterized as γ ( s ) : = ( ρ ( s ) cos θ ( s ) , ρ ( s ) sin θ ( s ) ) , for s [ 0 , 1 ] , with γ ( 0 ) = P + , γ ( 1 ) Γ 0 and E ( γ ( s ) ) strictly increasing with s. The target set is the symmetric of with respect to the x-axis. It will be denoted by Y and parameterized by reversing the angle θ ( s ) .

Using the same polar coordinates to represent the solutions of (2.7), we can express the final points Φ 0 T ( γ ( s ) ) as

Φ 0 T ( γ ( s ) ) = ( ρ ( T , s ) cos θ ( T , s ) , ρ ( T , s ) sin θ ( T , s ) ) ,

with

ρ ( 0 , s ) = ρ ( s ) , θ ( 0 , s ) = θ ( s ) .

We notice that

( x ( t ) , y ( t ) ) = ( ρ ( t , s ) cos θ ( t , s ) , ρ ( t , s ) sin θ ( t , s ) )

is the solution of (2.7) (at the time t), which departed from the point γ ( s ) at the time t = 0 . Observe that, for every s [ 0 , 1 ] , the map

t θ ( t , s )

is strictly decreasing on [ 0 , + ) (this is an equivalent way to express the fact that the solutions turn around the origin in the clockwise sense).

With these positions, we obtain an external solution homoclinic to P + if and only if there exists s ˆ ] 0 , 1 [ such that Φ 0 T ( γ ( s ˆ ) ) Y . This happens if and only if

ρ ( T , s ˆ ) = ρ ( s ˆ ) and θ ( T , s ˆ ) = θ ( s ˆ ) 2 j π , (2.12)

for some nonnegative integer j. In this case, if we denote by x ˆ ( ) the corresponding homoclinic solution of (2.1) such that ( x ˆ ( t ) , x ˆ ( t ) ) = γ ( s ˆ ) , we have that E ( x ˆ ( t ) , x ˆ ( t ) ) ] , 0 [ for all t R . Moreover, if j = 0 , x ˆ ( t ) > ξ 1 for all t [ t , t + ] , while x ˆ ( t ) has precisely 2j zeros in ] t , t + [ if j 1 .

Using the fact that the energy level lines of (2.7) in the region W { O } = { ( x , y ) : ξ 0 x ξ 0 , 0 < E ( x , y ) 0 } are strictly star-shaped with respect to the origin and symmetric with respect to the x-axis, we find that (2.12) holds (for some nonnegative integer j) if and only if

θ ( T , s ˆ ) = θ ( s ˆ ) 2 j π (2.13)

is satisfied (for the same j).

By virtue of (2.13), we can refer only to the angular coordinates, hence we will obtain solutions as follows. Since the point Φ 0 T ( γ ( 1 ) ) lies on Γ 0 (which is the trajectory of a heteroclinic solution of (2.7)), it can never reach the point Q + . Accordingly, the angle θ ( T , 1 ) satisfies

θ ( T , 1 ) > 0 , T 0 .

On the other hand, as γ ( 0 ) = P + , we have that θ ( 0 , 0 ) = θ ( 0 ) = 0 , therefore

θ ( T , 0 ) < 0 , T > 0 .

As a consequence, θ ( T , 0 ) + θ ( 0 ) < 0 < θ ( T , 1 ) + θ ( 1 ) and the intermediate value theorem ensures the existence of an s ˆ such that θ ( T , s ˆ ) = θ ( s ˆ ) . In this way we have found a positive homoclinic solution to P + , independently of the length T of the time interval [ t , t + ] .

Suppose now that T > n τ for some positive integer n. In this case,

θ ( T , 0 ) < 2 n π , T > 0 .

As a consequence, for every integer j with 1 j n , we have θ ( T , 0 ) + θ ( 0 ) < 2 j π < 0 < θ ( T , 1 ) + θ ( 1 ) , and again the intermediate value theorem ensures the existence of an s ˆ such that θ ( T , s ˆ ) = θ ( s ˆ ) 2 j π . This ends the proof for homoclinic solutions.

For the search of heteroclinic solutions from P + to P , we follow a similar procedure, choosing as a target the set Y , which is the symmetric of with respect to the y-axis (here we exploit the oddness of the nonlinear term x 3 ). In terms of polar coordinates, the desired solutions will be obtained if and only if there exists s ˆ ] 0 , 1 [ such that

θ ( T , s ˆ ) = θ ( s ˆ ) ( 2 j 1 ) π (2.14)

for some integer j 1 . Then the proof can be concluded as above, via the intermediate value theorem. We observe that the case of heteroclinic solutions from P to P + is analogous, thus it will be omitted. □

The next theorem deals with the internal connections. For this result, it is useful to recall the rotational gap defined in (2.11).

Theorem 2.2 (Internal connections)

Under the above assumptions, the following results hold for solutions which satisfy the energy condition

x ( t ) 2 + 2 k x ( t ) 2 μ 0 x ( t ) 4 [ 1 , ] , t R :

there always exist a solution homoclinic to P + and a solution homoclinic to P ;

there always exist a heteroclinic solution from P + to P and a heteroclinic solution from P to P + ;

if t + t > τ 1 τ τ τ 1 , there exist at least 2 w ( t + t ) + 1 solutions homoclinic to P + and 2 w ( t + t ) + 1 solutions homoclinic to P ;

if t + t > τ 1 τ τ τ 1 , there exist at least 2 w ( t + t ) + 1 heteroclinic solutions from P + to P and 2 w ( t + t ) + 1 from P to P + .

Proof As the first step, we focus our attention on the search of solutions homoclinic to the point P + . The case of homoclinics to P is analogous, thus it will be omitted.

Our goal is to find points on L = L ( P + , P ) ¯ which are moved by the Poincaré map Φ 0 T on L = L ( P , P + ) ¯ .

As a preliminary remark, we note that Φ 0 T ( L ) is a simple arc contained in ℰ, connecting its end points Φ 0 T ( P + ) and Φ 0 T ( P ) on Γ through the ‘intermediate’ point Φ 0 T ( R ) on Γ 1 . Observe also that Φ 0 T ( P + ) and Φ 0 T ( P ) are antipodal. In fact, P and P + are antipodal and Φ 0 T is an odd map (indeed, ( x ( t ) , y ( t ) ) is a solution of (2.7) if and only if ( x ( t ) , y ( t ) ) is a solution of the same equation). Now, if the trivial situation { Φ 0 T ( P ) , Φ 0 T ( P + ) } = { P , P + } occurs, we have Φ 0 T ( L ) L = { P , P + } . Otherwise, we find that Φ 0 T ( P ) and Φ 0 T ( P + ) belong to the two different components of E L . Since Φ 0 T ( P ) and Φ 0 T ( P + ) are the end points of the arc Φ 0 T ( L ) E , by an elementary connectivity argument, we conclude that Φ 0 T ( L ) L . This proves the first assertion of the theorem.

For the search of heteroclinic solutions from P + to P , we adopt a similar procedure, choosing as a target the set ℒ itself. The fact that, for every T > 0 , there is always at least one intersection point between Φ 0 T ( L ) and ℒ follows by the same argument developed for homoclinic solutions. Namely, if the trivial situation { Φ 0 T ( P ) , Φ 0 T ( P + ) } = { P , P + } occurs, we have Φ 0 T ( L ) L = { P , P + } . Otherwise, we find that Φ 0 T ( P ) and Φ 0 T ( P + ) belong to the two different components of E L . Therefore, it follows that Φ 0 T ( L ) L .

We study now the problem of multiplicity of solutions. As before, we consider at first the case of solutions homoclinic to P + .

We introduce in the inner region ℰ of a phase-plane for (2.7) a system of modified polar coordinates ( θ , E ) with center in the origin. In this system, every point is determined by its angular coordinate θ and its energy E defined in (2.8). Since the energy level lines in ℰ are strictly star-shaped with respect to the origin, we obtain a coordinate system equivalent to the polar one.

Observe that the set L L intersects every half-line from the origin exactly in one point. Accordingly, we can parameterize those points using the angular coordinate. With this convention, the initial points in the arc ℒ are parameterized as γ ( s ) for s [ π , 2 π ] (s is the angle) with γ ( π ) = P , γ ( 2 π ) = P + and γ ( 3 π / 2 ) = R . As a consequence, E ( γ ( s ) ) is strictly decreasing for s [ π , 3 π / 2 ] and strictly increasing in [ 3 π / 2 , 2 π ] . The target set L (the symmetric of ℒ with respect to the x-axis) is parameterized with the angle s [ 0 , π ] .

For the specific case of system (2.6), an analytic expression for ℒ as a graph in the ( θ , E ) -plane, is given as follows:

E ( θ ) : = E ( ρ ( θ ) cos θ , ρ ( θ ) sin θ ) = ρ 2 ( θ ) sin 2 θ + 2 k ρ 2 ( θ ) cos 2 θ μ 0 ρ 4 ( θ ) cos 4 θ , (2.15)

with

ρ ( θ ) = ( 2 k / μ 1 ) ( sin 2 θ + 4 k cos 2 θ ) 1 / 2 sin θ , θ [ π , 2 π ] .

Hence, the natural parametrization of ℒ in the ( θ , E ) -plane is given by γ ( s ) = ( s , E ( s ) ) for s = θ [ π , 2 π ] .

Using the same modified polar coordinates to represent the solutions of (2.7), we can express the final points Φ 0 T ( L ) in the ( θ , E ) -plane by means of their angular coordinate θ ( t , s ) and energy E ( t , s ) . As in the previous proof, we observe that, for every s [ π , 2 π ] , the map t θ ( t , s ) is strictly decreasing on [ 0 , + ) , while the energy is constant with respect to t. Observe also that in the new coordinates ( θ , E ) , system (2.7) becomes

{ θ = sin 2 θ 2 k cos 2 θ + 2 μ 0 Λ ( θ , E ) cos 4 θ , E = 0 , (2.16)

with

Λ ( θ , E ) : = 2 E sin 2 θ + 2 k cos 2 θ + ( ( sin 2 θ + 2 k cos 2 θ ) 2 4 μ 0 E cos 4 θ ) 1 / 2 .

In the ( θ , E ) -plane, the points of ℒ move with a negative angular speed along the lines of constant energy. Thus, under the action of the flow associated to (2.16), the points of ℒ shift from the right to the left in the strip

E ( , 2 π ] : = { ( θ , E ) : θ 2 π , 1 E } .

As we have explained before, we obtain an internal solution homoclinic to P + whenever there is a point z L such that Φ 0 T ( z ) L , with T = t + t . In the ( θ , E ) -plane, this target set is expressed as the union of the graphs

Ξ : = j = 0 { ( θ , ϒ ( θ ) ) : 2 j π θ π 2 j π } ,

where

ϒ ( θ ) : = E ( ϱ ( θ ) cos θ , ϱ ( θ ) sin θ ) = ϱ 2 ( θ ) sin 2 θ + 2 k ϱ 2 ( θ ) cos 2 θ μ 0 ϱ 4 ( θ ) cos 4 θ ,

with

ϱ ( θ ) = ( 2 k / μ 1 ) ( sin 2 θ + 4 k cos 2 θ ) 1 / 2 + sin θ

(see Figures 10, 11 for a graphical description).

Figure 10. The initial setand the target set L represented in the ( θ , E ) -plane. The present figure is drawn for the parameters k = 1 , μ 1 = 2 , μ 0 = 1 / 2 . The internal region ℰ corresponds to the strip R × [ 1 , ] . Since we are interested in the evolution of the set ℒ through the Poincaré map, we consider only the half-strip E ( , 2 π ] = ( , 2 π ] × [ 1 , ] . For the flow associated to (2.16), all the points of ℒ move from the right to the left on lines parallel to the θ-axis. The point R (as well as the points on the line E = 1 ) moves faster than the points P ± which are on the line E = . During all the evolution, the distance between the images of P and P + remains constantly equal to π.

Figure 11. The present figure describes the same situation of Figure7in the angle-energy coordinates (with the same choice of the coefficientsk, μ 1 , μ 0 and time gap T = t + t = 50 ).

In this setting, we obtain an internal solution homoclinic to P + if and only if there exists s ˆ [ π , 2 π ] such that Φ ( γ ( s ˆ ) ) Ξ , where we have denoted by Φ the Poincaré map for (2.16) in the time interval [ 0 , T ] . Of course Φ is exactly Φ 0 T (which was the Poincaré map associated to (2.7)) in the new ( θ , E ) -coordinates. Using the parameterized curve γ ( s ) for the initial points, we can express Φ as

Φ ( γ ( s ) ) = ( θ ( T , s ) , E ( s ) ) for  s [ π , 2 π ]

(recall that γ ( s ) = ( s , E ( s ) ) for s = θ [ π , 2 π ] ).

The points P and P + are antipodal, lie on the same energy line Γ and move with the same angular speed. In the ( θ , E ) -plane, P and P + are expressed by P = ( π , ) , P + = ( 2 π , ) . Therefore,

θ ( t , π ) θ ( t , 2 π ) = π , t 0 . (2.17)

Moreover, for every integer j 0 , we have that

θ ( t , π ) π j π if and only if t j τ 2

(this follows from the fact that the orbit Γ has a period τ and is symmetric with respect to the x-axis).

In the same plane, the point R is indicated by R = ( 3 π / 2 , 1 ) . Using the fact that the orbit Γ 1 has a period τ 1 and is symmetric with respect to the y-axis, we find that, for every integer i 0 ,

θ ( t , 3 π / 2 ) π 2 2 i π if and only if t τ 1 2 + i τ 1 .

Suppose now that

2 T τ 1 2 T τ > 2 . (2.18)

Then, the open interval ] 2 T / τ , 2 T / τ 1 [ contains at least one odd integer.

Let m o ( T ) be the largest odd integer contained in ] 2 T / τ , 2 T / τ 1 [ . From T > m o ( T ) τ 1 / 2 , we have that

θ ( T , 3 π / 2 ) < π 2 ( m o ( T ) 1 ) π = 3 π 2 m o ( T ) π . (2.19)

Suppose also that j is a positive integer with 2 T / τ < j m o ( T ) . In this case, from T < j τ / 2 , we have that

θ ( T , π ) > π j π and θ ( T , 2 π ) > 2 π j π . (2.20)

Let

L left : = { ( θ , E ( θ ) ) : π θ 3 π / 2 }

and

L right : = { ( θ , E ( θ ) ) : 3 π / 2 θ 2 π } ,

with L = L left L right be the left and the right parts of ℒ in the ( θ , E ) -plane. Similarly, we define

L left : = { ( θ , ϒ ( θ ) ) : 0 θ π / 2 }

and

L right : = { ( θ , ϒ ( θ ) ) : π / 2 θ π } .

By definition, for any nonnegative integer i, the set L left ( 2 i π , 0 ) is a simple arc connecting the points ( 2 i π , ) and ( 2 i π + ( π / 2 ) , 1 ) in the rectangle [ 2 i π , 2 i π + ( π / 2 ) ] × [ 1 , ] and, similarly, L right ( 2 i π , 0 ) is a simple arc connecting the points ( 2 i π + ( π / 2 ) , 1 ) and ( 2 i π + π , ) in the rectangle [ 2 i π + ( π / 2 ) , 2 i π + π ] × [ 1 , ] . On the other hand, the set Φ ( L left ) is a simple arc connecting the points ( θ ( T , π ) , ) and ( θ ( T , 3 π / 2 ) , 1 ) in E ( , 2 π ] and, similarly, Φ ( L right ) is a simple arc connecting the points ( θ ( T , 3 π / 2 ) , 1 ) and ( θ ( T , 2 π ) , ) = ( π + θ ( T , π ) , ) in the same strip.

Suppose now that j = m o ( T ) = 2 k + 1 (in this case, j is odd). In such a situation, from (2.19) and (2.20) we have that

θ ( T , 3 π / 2 ) < π 2 2 k π , θ ( T , π ) > 2 k π , θ ( T , π ) > π 2 k π ,

for 2 k = j 1 an even integer. Hence the points ( θ ( T , 3 π / 2 ) , 1 ) and ( θ ( T , π ) , ) are separated by the arc L left ( 2 k π , 0 ) , while the points ( θ ( T , 3 π / 2 ) , 1 ) and ( θ ( T , 2 π ) , ) are separated by L left ( 2 k π , 0 ) and L right ( 2 k π , 0 ) . As a consequence, Φ ( L ) Ξ contains at least three points, precisely the nonempty intersections of Φ ( L left ) ( L left ( 2 k π , 0 ) ) , Φ ( L right ) ( L left ( 2 k π , 0 ) ) and Φ ( L right ) ( L right ( 2 k π , 0 ) ) .

Suppose that j = m o ( T ) 1 (in this case, j is even). In such a situation, from (2.19) and (2.20) we have that

θ ( T , 3 π / 2 ) < π 2 2 k π , θ ( T , π ) > π 2 k π , θ ( T , π ) > 2 π 2 k π

for 2 k = m o ( T ) 1 an even integer. Hence the points ( θ ( T , 3 π / 2 ) , 1 ) and ( θ ( T , π ) , ) are separated by the arcs L left ( 2 k π , 0 ) and L right ( 2 k π , 0 ) . The points ( θ ( T , 3 π / 2 ) , 1 ) and ( θ ( T , 2 π ) , ) are separated by L left ( 2 k π , 0 ) , L right ( 2 k π , 0 ) and L left ( ( 2 k 1 ) π , 0 ) . As a consequence, Φ ( L ) Ξ contains at least five points: two points coming from the intersections of Φ ( L left ) with the target set Ξ and three from the intersections of Φ ( L right ) .

Proceeding by induction, with the same argument, the following result is obtained.

Claim 1 Φ ( L ) Ξ has at least 2 n # + 1 solutions, where n # is the number of integers less than or equal to m o ( T ) which are contained in the open interval ] 2 T / τ , 2 T / τ 1 [ .

A completely symmetric argument leads to the same multiplicity result for solutions which are homoclinic to P .

At last, we look for a multiplicity result for heteroclinic solutions from P + to P . We follow step by step the argument described in the part of the proof devoted to the search of multiple internal homoclinic solutions and, therefore, we transform our equation in the ( θ , E ) -coordinates. As we have explained before, we obtain an internal heteroclinic (from P + to P ) whenever there is a point z L such that Φ 0 T ( z ) L , with T = t + t . In the ( θ , E ) -plane, this target set is expressed as the union of the graphs

Ξ : = j = 0 { ( θ , E ( θ ) ) : π 2 j π θ 2 π 2 j π } ,

where E ( θ ) is the function defined in (2.15).

As before, we assume the validity of condition (2.18). Then the open interval ] 2 T / τ , 2 T / τ 1 [ contains at least one even integer.

Let m e ( T ) be the largest even integer contained in ] 2 T / τ , 2 T / τ 1 [ . From T > m e ( T ) τ 1 / 2 , we have that

θ ( T , 3 π / 2 ) < 3 π 2 m e ( T ) π .

For any positive integer j with 2 T / τ < j m e ( T ) , we have that (2.20) holds too. At this point, we have simply to repeat (with obvious changes due to the fact that the target set is Ξ which is a shift of Ξ by ( 2 p i , 0 ) ) the argument presented above for the case of homoclinic connections and obtain the following conclusion.

Claim 2 Φ ( L ) Ξ has at least 2 n # + 1 solutions, where n # is the number of integers less than or equal to m e ( T ) which are contained in the open interval ] 2 T / τ , 2 T / τ 1 [ .

A completely symmetric argument leads to a multiplicity result for heteroclinic solutions from P to P + .

As final step, we have just to show how the condition on the rotational gap allows us to achieve the conclusion from Claim 1 and Claim 2 that we have obtained along the proof. Now, it is sufficient to observe that t + t = T > τ 1 τ τ τ 1 implies (2.18) and, therefore, the open interval ] 2 T / τ , 2 T / τ 1 [ contains at least two integers. More precisely, by the definition of w ( T ) , the interval ] 2 T / τ , 2 T / τ 1 [ contains at least w ( T ) + 1 integers and, therefore,

n # , n # w ( T ) .

This concludes the proof of the theorem.  □

### 3 Remarks and related results

We end the paper with a list of remarks about possible variants and extensions of the main results obtained for equation

1. A problem which naturally arises from the analysis that we have performed concerns what happens if, for a stepwise weight function satisfying (2.5), we suppose that 0 < μ 1 < μ 0 . Repeating the preliminary phase-plane analysis of Section 2.2, one can easily check that for any gap T = t + t > 0 , there always exist a solution homoclinic to P + and another homoclinic to P , as well as a heteroclinic from P + to P and another one from P to P + . However, in general, if 0 < μ 1 < μ 0 , one cannot obtain multiplicity results like those achieved in Section 2 without some extra assumptions on h ( x ) . This is the reason for which, in the study of equation (2.1), we have considered only the case 0 < μ 0 < μ 1 . For different examples on related equations in which a weight coefficient can be above or below its limits at infinity, see, for instance, [4-6].

On the other hand, if we assume that the weight function changes its sign, with μ 0 < 0 < μ 1 , some interesting multiplicity results could be produced. In fact the phase portrait of (1.1) in the time interval [ t , t + ] shows a global center, hence if t + t is large, a lot of connections between the unstable and stable manifolds of the points P ± can be obtained. For the sake of conciseness, we omit the study of this latter situation, which is beyond the goal of the present paper.

2. We notice that the same argument of the proofs applies to an equation of the form

for σ > 0 , h ( x ) = x f ( x ) , and with f ( x ) satisfying condition (∗), if we assume that a ( t ) is a stepwise weight function, playing the same role of g ( t ) in (2.1). In such a case, the saddle points become P ± = ( ± ξ 1 , 0 ) with ξ 1 : = f 1 ( σ / μ 1 ) , where f 1 is the inverse of f restricted to [ 0 , + ) . In this way, we can apply our results to nonlinearities like those considered in [2,27]. More precisely, Theorem 2.1 holds without any further assumption, while for Theorem 2.2 we need to require a gap between the periods of the orbits Γ 1 and Γ . In order to obtain this gap, we can apply (for instance) some results ensuring the monotonicity of the time-map (like [[28], Theorem A]).

3. Due to the special form of the weight coefficient, it is standard to verify (via a simple rescaling procedure) that (1.1) is equivalent to an equation of the form

x + b ( t ) x a h ( x ) = 0 , (3.1)

where a > 0 and b ( t ) is a stepwise coefficient. In this manner, one can deal with some nonlinear Schrödinger-type equations related to the case of potential wells or potential walls [29,30].

4. The approach used in the proofs, based on the properties of the Poincaré map, guarantees that our results are stable with respect to small perturbations. More precisely, fixed a suitable length T = t + t for the time interval, Theorems 2.1 and 2.2 provide a lower bound for the number of solutions. We can state that the same lower bound persists for a small perturbation of the coefficient in the L 1 -norm on [ t , t + ] . Therefore, the assumption that the weight a ( t ) in (1.1) or b ( t ) in (3.1) are stepwise functions can be slightly relaxed, so we can ‘smooth’ them.

5. With reference to equations (1.1) or (3.1) with stepwise coefficients, we observe that our approach can be adapted to boundary value problems on a compact interval [ 0 , T ] like, e.g., the Dirichlet (two-point) or the Neumann problem. In these cases, we have to find solutions connecting given lines which depart from the origin. For the sake of conciseness, we cannot describe the most general situation, but we just outline a possible application for the Neumann problem

{ x + 2 k x 2 g (