_{1}

^{*}

This paper proposes a method for constructing partial differential equation (PDE) systems with chaotic solitons by using truncated normal forms of an ordinary differential equation (ODE). The construction is based mainly on the fact that the existence of a soliton in a PDE system is equal to that of a homoclinic orbit in a related ODE system, and that chaos of ?i’lnikov homoclinic type in the ODE imply that the soliton in the PDE changes its profile chaotically along propagation direction. It is guaranteed that the constructed systems can self-generate chaotic solitons without any external perturbation but with constrained wave velocities in a rigorously mathematical sense.

Chaotic solitons have been a subject of many theoretical papers over the last decades of years. They have been experimentally observed in several physical systems [1- 4]. Despite this large amount of effort, many key issues remain open. The central question addressed in this paper is: How can one construct a partial differential equation (PDE) system that self-generates a chaotic solitary-wave pulse that exists in a rigorous sense? There are few rigorous results that addresses this question. In most cases, rigorous results for the generation of chaotic solitons are obtained through perturbations to a known system, e.g., Schrödinger equation and Ginzburg-Landau equation [5- 8], that can generate solitons. Unfortunately, this method is not regular since, in practice, it can be applied only to particular examples. Another rigorous approach is that of C. L. Zheng [

In this paper, we propose a method for constructing a PDE system self-generating chaotic solitons. The construction is based on the truncated normal forms of an autonomous ODE and on the relationship between ODEs and PDEs when solitary waves are under consideration. The constructed system can self-generate a chaotic solitary-wave pulse without any external perturbation. Moreover, the existence of the solitary-wave pulse in the PDE system is ensured by that of one homoclinic orbit of an ODE system whereas its chaotic behavior (i.e., the soliton changes its profile chaotically along propagation direction) results from the chaos property of the ODE system which is guaranteed by Ši’lnikov’s homoclinic theorem [

Since finding a soliton in a PDE system can be changed into proving the existence of a homoclinic orbit in a related ODE system, in this section we will introduce a class of autonomous ODE systems as candidates for the construction of PDE systems that can self-generate chaotic solitons. The introduction of the candidate ODE systems is based on amplitude equations which are derived from normal form theory. Then, we discuss the existence of homoclinic orbits in the ODE systems. We do not prove the existence of the homoclinic orbit, and instead obtain some so-called principle homoclinic orbits by applying perturbation theory. Furthermore, we show that such homoclinic orbits are of Ši’lnikov’s type, implying that the ODE systems have chaos of Ši’lnikov’s type.

We begin by a general theory about the standard form of nonlinear ODE systems. In the analysis of the stability of static states of autonomous nonlinear ODE systems, linear theory produces normal modes that vary like, where t is time, k is an integer, and is a root of the characteristic value equation of the linear theory. Consider an autonomous nonlinear ODE system having a polycritical state in which for l roots and with for m pairs of complex conjugate roots. Suppose that, for this state, is bounded away from zero for all the other roots of the characteristic equation. By using normal form theory, we can obtain amplitude equations of the form:

where with is the amplitude vector, is a strictly nonlinear vector-valued function, and is a set of parameters. We restrict ourselves to cases where the number of parameters in the characteristic polynomial for is equal to n (called criticality). In this case, the characteristic equation (called the critical polynomial) for matrix may be expressed as

We are interested here in the leading-order terms in, e.g., for, they are of order. These leading terms in are dictated by. For the critical condition, in the case of, the characteristic equation has zero as a double root. We refer to this case as when is an irreducible matrix. Similarly, we have other (see [

where the are constants which can be analytically calculated [

Let be a small positive number and introduce the scaled variables

where. If we arrange matters so that the are of order unity, the values of are near to the critical condition. That is, we position the system in that part of parameter space where, for,. Then Equation (3) becomes

If we keep only the leading order terms in g, we may reduce Equation (5) to

where and we are ignoring terms. Once we are within the region of parameter space delimited by, we may use Equation (6) as the amplitude equation for the configuration, for any n, with the error of order.

We are interested especially in the case of. In this case, we have the following truncated normal form as the amplitude equation for (see [13,14]):

where and are constants. Equation (7) has two singular points and. We look for solutions of the form. If is small, , and we obtain

The parameter is a measure of the dissipation, and assumed to be positive. A change in the sign of essentially interchanges the rules of the two singular points, so we need consider only. When and, we can obtain some principle homoclinic orbits by applying perturbation theory [

A particular example is the asymptotic normal form for which is of the form:

where X has been scaled to give a convenient coefficient to. Based on such an ODE and inspired by the relation between ODE and PDE when the travelling waves are considered, we construct the following model:

where a and are positive constants. Assume that

Equation (10) has travelling waves of the standard form, where and c represents travelling wave speed (assume c > 0). Then, U satisfies Equation (9) with. Similarly, we can construct the classical Korteweg-de Vries (KdV) equation which has been applied to the optical fibers communication field if we set, , and, , based on the ODE (8). Now, we seek a solitary wave solution, that is, we find the solution of the following equation subjective to the constrain aswhere the subscript stands for derivatives of function U with respect to the argument. For this, we rewrite Equation (9) into an ODE

with the condition, and as, where. Note that finding a soliton solution of Equation (10) satisfying the condition is currently changed into finding a homoclinic orbit of Equation (11) based at the equilibrium point.

Fourth-order Rounge-Kutta method is used to solve the system (11) with time step size equal to 0.01, the absolute and relative error equal to. Numerical simulation shows that the system (11) is chaotic when and, as shown in

We can verify analytically that Equation (10) indeed has a Ši’lnikov homoclinic orbit and further chaos of Ši’lnikov type from a theorem of Ši’lnikov [

In addition, the characteristic polynomial of Jacobian matrix of Equation (10) evaluated at is which has three roots: one negative and a conjugate pair of complex with the positive real part for and. A rigorously theoretical proof for the existence of chaos of Ši’lnikov type may be referred to a standard treatment in [

Based on normal form theory of ODEs, we have proposed a general method of constructing systems that selfgenerate chaotic solitons. By adjusting some system parameters in truncated normal forms of ODEs (e.g., Equation (8)), one can construct many systems with chaotic solitons. Note that the constructed system (10) is conservative. It seems the first example that self-generates chaotic solitons.