On Using Curvature to Demonstrate Stability

A new approach for demonstrating the global stability of ordinary differential equations is given. It is shown that if the curvature of solutions is bounded on some set, then any nonconstant orbits that remain in the set, must contain points that lie some minimum distance apart from each other. This is used to establish a negative-criterion for periodic orbits. This is extended to give a method of proving an equilibrium to be globally stable. The approach can also be used to rule out the sudden appearance of large-amplitude periodic orbits.


Introduction
A key issue in the analysis of a system of ordinary differential equations is to determine the long-term dynamics. Generally, it is desirable, if possible, to know whether a system settles down to an equilibrium, approaches a limit cycle, or exhibits some other long-term behavior. In this paper, we give a method that uses the curvature of solutions to help show that a system approaches an equilibrium.
A nice feature of considering the curvature is that the calculation only involves local information. Furthermore, it is not necessary to calculate solutions to the differential equations since the curvature calculation uses the derivative of a solution, and this is given by the vector field that defines the flow.
On any set for which the curvature is bounded, there is a minimum length of a closed orbit not including equilibria . Also, in such a set, any closed orbit must contain points that are a certain distance away from each other. The smaller the curvature, the farther apart certain points on the orbit must be from each other. This idea is used to formulate a negative criterion for periodic orbits. The method is then extended to give a criterion that shows that all forward trajectories leave the set.

D i fferential Equations and Nonlinear Mechanics
Unfortunately, if the eigenvalues of the Jacobian at an equilibrium are not all equal, then the curvature near the equilibrium is unbounded. This makes it necessary to deal with the behavior near an equilibrium separately. For instance, one may use a Lyapunov function near an equilibrium to show that a certain neighborhood is contained in the equilibrium's basin of attraction. Then, analysis of the curvature can be used to show that all other solutions enter this set by showing that they leave the complement , making the stability global.
This approach can also be useful in ruling out nonlocal bifurcations as parameters are varied. It is often possible to show that an equilibrium is locally asymptotically stable for all parameter values, and that it is globally stable for a particular set of parameter values. As parameters are varied, the continued local stability precludes the possibility of local bifurcations such as Hopf. Nonlocal bifurcations, however, may still occur. The method presented here can be used to preclude the existence and sudden appearance of largeamplitude periodic orbits and, therefore, help rule out nonlocal bifurcations.
We now give a brief outline of the paper. In Section 2, the curvature of parametric curves is discussed. This is applied to periodic solutions of differential equations in Section 3. In Section 4, results are given for omega limit points. Section 5 stability gives stability theorems. A discussion of how to apply the results is included in Section 6.

Curvature Calculations
Consider a curve Γ given by the function r : R→R n , where r is C 2 and r is nonzero for all t ∈ R.
In this paper, we use the convention that represents differentiation with respect to time. Let s be the arclength along Γ from some reference point, that is, s t t t 0 r τ dτ. We use the convention that when the arguement of a function is changed to s, the function is considered to be reparameterized in terms of arc length. Note that ds/dt r t , and so d/ds d/dt dt/ds 1/ r t d/dt . For each t, the unit tangent at the point r t is

2.2
The curvature 1 of Γ is defined to be K dT/ds and is given by

2.4
Let T 0 be a unit vector in R n . Let θ t be the angle between T t and T 0 : Proof. Noting that cos θ T 0 ·T , we have − sin θ dθ/ds T 0 ·dT/ds. Letting β be the angle between dT/ds and T 0 , then, T 0 ·dT/ds dT/ds cos β . Thus,

2.6
We now show that |cos β / sin θ | is at most one. Since T s is a unit vector for all s, it follows that dT/ds is orthogonal to T . Since T 0 is a unit vector separated from T by angle θ, the component of T 0 normal to T has magnitude sin θ . Thus, the component of T 0 in the direction of dT/ds has magnitude less than or equal to sin θ . On the other hand, this magnitude is given by |cos β |, showing that |cos β / sin θ | ≤ | sin θ / sin θ | 1. Combining with 2.6 yields the proposition. Now, suppose Γ is a simple closed curve given by the C 2 ω-periodic function r : R→R n , where r is nonvanishing. Let x 0 ∈ Γ. Choose t 0 ∈ R such that r t 0 x 0 , and let T 0 T t 0 , that is, T 0 is the unit tangent to Γ at x 0 .
Since r is ω-periodic, it follows that 0 r t ω − r t t ω t r τ dτ for all t. Hence, 0 t ω t r τ ·T 0 dτ. Combining this with the fact that r t ·T 0 is ω-periodic, continuous and positive at t t 0 , we see that r t ·T 0 oscillates in sign and, therefore, must have at least two zeros in every interval of length ω. Let t − , t be the maximal interval containing t 0 such that r t ·T 0 is positive. Note that θ t − θ t π/2. Let d x, y be the Euclidean distance between points x, y ∈ R n . We now work toward finding a lower bound on the distance between the points r t 0 and r t :

4 D i fferential Equations and Nonlinear Mechanics
Let K Γ be the least upper bound of K on Γ. Then, recalling that ds/dt r t gives This leads to the following result.
Theorem 2.2. Let Γ be a simple closed curve and let K be an upper bound for the curvature on Γ. Then, given a point By not requiring that the curve Γ be closed, we obtain the following result, which will be useful in Section 4. Theorem 2.3. Let K be an upper bound for the curvature on a curve Γ which is given by r : R→R n , with r nonvanishing. Suppose further that there exist t − , t 0 , t ∈ R, with t − < t 0 < t such that the tangent at x 0 r t 0 is orthogonal to the tangents at x − r t − and x r t , and r t ·r

Periodic solutions of differential equations
Consider the differential equation where f : R n →R n is C 1 . We denote by x t; x 0 the solution to 3.1 which passes through x 0 at time t 0. Suppose x t; x 0 is a nonconstant solution to 3.1 . Noting that x d x /dt ∂ x /∂x dx/dt ∂f/∂x f x , it follows from 2.4 that the curvature at x is Thus, as long as the point x is not an equilibrium of 3.1 , the curvature at x of the solution through x can be precisely calculated without any explicit knowledge of the solution. We now relate Theorem 2.2 to periodic solutions of ordinary differential equations.
Then, any periodic solution to 3.1 which lies entirely in Δ has maximum curvature less than or equal to K Δ .

Omega limit points
Theorem 4.1. Suppose Δ ⊆ R n has the negative property, and that ϕ t is a forward trajectory of 3.1 which is contained entirely in Δ. Then, the omega limit set of ϕ consists entirely of equilibria which lie in the closure of Δ, but not in Δ.
Proof. Suppose x is an omega limit point of ϕ t , where ϕ t ∈ Δ for all t ≥ 0, and that f x / 0. Then, there is an increasing sequence of times {t n } which limits to ∞ such that lim n→∞ ϕ t n x. Then

4.1
Fix , δ > 0 such that f y ·f z > 2 if y, z ∈ N δ , the open ball centred at x with radius δ. Choose τ > 0 so that f x t; y ·f z > if y, z ∈ N δ and |t| < τ.
Taking the dot product of each side of 4.1 with a vector v ∈ f N δ gives For sufficiently large n, we have ϕ t n ∈ N δ and, therefore, f ϕ t ·v > for |t − t n | < τ. Thus, Hence, the integrand must oscillate between positive and negative for large enough t, in order that 4.2 be satisfied. Thus, for large enough n, there exist t − n and t n with t − n < t n < t n , such that f ϕ t − n and f ϕ t n are orthogonal to v. Furthermore, t − n and t n can be chosen so that f ϕ t ·v is positive on the interval with these endpoints.
The only restriction on v, thus far, is that v ∈ f N δ . In order for n to be sufficiently large for t − n and t n to be defined, it is only necessary that ϕ t n ∈ N δ . For large n, a change in the choice of v will change the values of t − n and t n but will not change the fact that they exist. Now, fix n to be sufficiently large that ϕ t n ∈ N δ , and choose v f ϕ t n .
Applying Theorem 2.3, with This contradicts the negative property, and so if x is an omega limit point, then we must have f x 0. Since Δ satisfies the negative property, there are no equilibria in Δ and, therefore, x is not an element of Δ. However, since x is an omega limit point of ϕ, it must be an element of the closure of Δ, completing the proof. Proof. Since Δ is compact, any forward trajectory contained entirely in Δ has a nonempty omega limit set Ω ⊆ Δ. By Theorem 4.1, Ω consists entirely of equilibria, but by the negative property Δ contains no equilibria, so Δ cannot contain an entire forward trajectory.

Stability theorems
Suppose D ⊆ R n is positively invariant under 3.1 . If Δ ⊆ D is compact and has the negative property, then by Corollary 4.2 there are no forward trajectories contained entirely in Δ. If it is known that all solutions intersecting D \ Δ limit to an equilibrium, then we can conclude that all solutions in D limit to an equilibrium giving the following result. where α : R n →R n is a diffeomorphism. Then, x ∂α/∂y y , and so

5.2
If the set α −1 Δ is compact and has the negative property for 5.2 , then Corollary 4.2 implies solutions of 5.2 leave α −1 Δ , which then implies solutions of 3.1 leave Δ. This allow Theorem 5.1 to be generalized as follows.
Theorem 5.3. Suppose D is positively invariant under 3.1 , Δ ⊆ D is compact and that all solutions in D \ Δ limit to an equilibrium. If there exists a diffeomorphism α : R n →R n such that α −1 Δ has the negative property for 5.2 , then all solutions in D limit to an equilibrium.

Discussion
If it is known that the curvature of a vector field is bounded on a set, then it follows that the rate at which trajectories turn is bounded. A consequence of this is that a closed orbit has a minimum size. In particular, there exist points on the orbit which are separated by a distance at least twice as large as the reciprocal of the bound on the curvature. This approach can be used to rule out periodic trajectories of a differential equation and to show that solutions which remain in a given set must limit to sets consisting of equilibria.
It may appear that this alone could be used to show global stability of an equilibrium in an invariant compact set, but there is a problem. Near an equilibrium the curvature is usually unbounded. This is true whenever the linearization at the equilibrium has at least two distinct eigenvalues.
Thus, in determining the global stability of an equilibrium, it would be necessary to deal with the behavior near the equilibrium in some other manner. For instance, one could use the linearization to obtain a quadratic form that could be used to show that a given neighborhood limits to the equilibrium. Then, by showing that the remainder of the compact set has the negative property, it would follow that the equilibrium is globally stable.