On exact discretization of the cubic-quintic Duffing oscillator

On exact discretization of the cubic-quintic Duffing oscillator


Application of intersection theory to construction of -point finite-difference equations associated with classical integrable systems is discussed. As an example, we present a few exact discretizations of one-dimensional cubic and quintic Duffing oscillators sharing form of Hamiltonian and canonical Poisson bracket up to the integer scaling factor.

1 Introduction

A completely integrable system on symplectic manifold with form of dimension is defined by smooth functions in involution

with the independent differentials at each cotangent space , .

If is a regular value of , then the corresponding level is a smooth -dimensional Lagrangian submanifold of . Geometrically this means that locally around the regular value the map collecting the integrals of motion is a Lagrangian fibration, i.e. it is locally trivial and the fibers are Lagrangian submanifolds.

Let us consider a finite-difference equation


relating points of submanifold . Ordinary finite difference equations of this type can be viewed as a dynamical -point map, see [15].

Because all points in (1.1) belong to the given Lagrangian submanifold we may suppose that the corresponding map preserves functions and symplectic form up to a scaling factor. Finite-difference equation sharing integrals of motion with the continuous time system and the symplectic structure is the so-called exact discretization of integrable systems. Nowadays, refactorization in the Poisson-Lie groups is viewed as one of the most universal constructions of finite difference equations (1.1), see discussion in [3, 4, 15, 18, 19, 27] and references within.

The idea is to identify points in (1.1) with intersection points of with auxiliary curve . If and are algebraic, then we can consider the standard equation for their intersection divisor

as the finite-difference equation (1.1) for the corresponding completely integrable system. Here is the intersection divisor of two algebraic varieties and is a suitable equivalence relation [6, 11, 13, 16]. Our main objective is to study properties of such -point finite-difference equations for different integrable systems [32, 33, 34, 35, 36]. In this paper we restrict ourselves by consideration of cubic and quintic nonlinear Duffing oscillators in order to clarify our view point on relations between the exact discretizations and the intersection divisors.

Thus, we consider integrable systems on two-dimensional plane with a pair coordinates and symplectic form . Because any smooth curve on the plane is a Lagrangian submanifold, we can directly apply classical intersection theory [2, 12] to exact discretization of one-dimensional Hamiltonian systems with the algebraic Hamilton function. Below we consider Hamiltonians associated with hyperelliptic curves on the projective plane defined by equation

at and various intersections of with the line, quadric and cubic on the plane defined by


From now on and mean coordinates on the projective plane, whereas and are coordinates on the phase space .

2 Cubic oscillator

The Hamilton function


and canonical Poisson bracket determine Hamiltonian equations


and equation of motion


for the generalized oscillator with the cubic nonlinearity [22].

At this integrable system is called a cubic Duffing oscillator without forcing. Duffing oscillators have received remarkable attention in recent decades due to the variety of their engineering applications. For instance magneto-elastic mechanical systems, large amplitude oscillations of centrifugal governor systems, nonlinear vibration of beams, plates and fluid flow induced vibration, seismic waves before earthquake, ecology or cancer dynamics, financial fluctuations and so on are modeled by the nonlinear Duffing equations.

In the numerical integration of nonlinear differential equations, discretization of the nonlinear terms poses extra ambiguity in reducing the differential equation to a discrete difference equation. For instance, in the framework of the standard-like discretization differential equation


can be transformed to the finite-difference equation

where is a discrete time interval [23, 24]. This equation may be reduced to the expression with the mapping function

or to the area preserving map on the plane

where is a rational control function [20, 21]. This integrable map admits the invariant integral

see details in [20, 21, 23, 24, 25, 26, 27], but it is not exact discretization of the Duffing oscillator, i.e. trajectories of the discrete flow do not coincide with the trajectories of the continuous flow.

2.1 Exact discretization and intersection divisors

In order to get exact discretization sharing integrals of motion with the continuous time system and the Poisson bracket we can start with the well-known analytical solutions of the Duffing equation (2.4), which are expressed via Jacobi elliptic functions.

Indeed, let us consider the equation (2.4) with initial condition

For and periodic solution is

For and periodic solution reads as

For and periodic solution has the form

Here and are the Jacobi elliptic functions. Discussion of the non-periodic solution can be also found in [23, 24].

Following [3] we can construct exact discretizations of the Duffing equation using these explicit solutions and well-known addition theorems for Jacobi elliptic functions, for instance

However, it is more easy and convenient to apply standard algorithms of the intersection theory for this purpose.

In order to apply the intersection theory to the exact discretization of the cubic oscillator we put and consider the corresponding level curve on the projective plane defined by equation


where . Any partial solution and of the Hamiltonian equations (2.2) at is a point of with abscissa and ordinate . It allows us to study relations between points of instead of relations between solutions of the differential equations (2.2).

Let be a smooth nonsingular algebraic curve on a projective plane. Prime divisors are rational points on denoted and is a point at infinity. Divisor

is a formal sum of prime divisors, and the degree of divisor is a sum deg of multiplicities of points in support of the divisor. Group of divisors is an additive Abelian group under the formal addition rule

Two divisors are linearly equivalent

if their difference is principal divisor

i.e. divisor of rational function on .

Intersection divisor of with some auxiliary smooth nonsingular plane curve

is equal to zero with respect to the linear equivalence of divisors. It allows us to identify intersection divisor with some finite-difference equation (1.1)


Here we divide intersection divisor in two parts

where prime divisors are parameters of discretization implicitly depending on .

2.2 Examples of the intersections

Let us consider the intersection of plane curve (2.5) with a parabola

and the corresponding intersection divisor of degree four, see [2], p.113 or [12], p.166. Following Abel we substitute into (2.5) and obtain the so-called Abel polynomial

Divisor of this polynomial on coincides with , i.e. roots of this polynomial are abscissas of intersection points and forming support of the intersection divisor .

At one of the intersection points is , see examples in Figure 1.

a) b)
Figure 1: Interection of curve (2.5) with parabola

In this case polynomial is equal to

Equating coefficients of one gets relation between abscissas of the remaining rational points and in support of the intersection divisor


If as in Figure 1a, we can define parabola using the Lagrange interpolation by any pair of points , or . For instance, taking the following pair of points one gets

which allows us to determine as functions on and . Substituting coefficients of into the equation (2.6) we obtain an explicit expression for abscissa as a function of coordinates and


If we have a double intersection point, for instance as in Figure 1b, then


where function is defined by due to the Hermite interpolation

In modern terms, we consider two partitions of the intersection divisor

Using brackets we separate a part of the intersection divisor which is necessary for polynomial interpolation of auxiliary curve . Because

these partitions can be rewritten as addition and doubling of prime divisors

where we use standard hyperelliptic inversion , see Figure 1.

At support of the intersection divisor consists of four rational points up to multiplicity. Let us consider the following partitions of this divisor

see Figure 2. In the first case parabola is defined by the Lagrange interpolation using three ordinary points and . In the second and third cases parabola is defined by the Hermite interpolation using either double and ordinary points or one triple point , respectively.

a) b) c)
Figure 2: Intersection of curve (2.5) with parabola

In the first case abscissa of the fourth intersection point is


where function is defined using coefficients of quadratic polynomial


In the second case expression for the abscissa looks like


Here function is defined via coefficients of the same polynomial and Hermite interpolation formulae

In the third case, when we consider tripling the prime divisor on

second abscissa is equal to


where function is defined via coefficients of the polynomial


At we have the intersection divisor of with line , which can be represented in the following form

It means that line is interpolated by two points and

whereas abscissas of remaining two points and are the roots of polynomial

Thus, are algebraic functions on coordinates and




In the generic cases, using intersection divisors of plane curve with auxiliary curves

we can describe multiplication of the prime divisor on any integer , which is a key ingredient of the modern elliptic curve cryptography, and other configurations of the prime divisors entering into the intersection divisor.

All the relations between abscissas (2.6-2.14) are well known, here we only repeat the fairly simple calculations based on Abel’s ideas and their geometric interpretation proposed by Clebsch, see the historical comments in [16]. The modern intersection theory gives a common language for the compact description of these partial cases of intersections [6, 11], whereas modern cryptography equips us with the effective algorithms for such computations [14].

2.3 Examples of finite-difference equations

Our aim is to interpret well-studied relations between prime divisors as the finite-difference equations (1.1) realizing various exact discretizations of the given Hamiltonian system and to study the properties of the corresponding discrete maps. For this purpose, we will identify partial solutions and of the Hamilton equations (2.2)

at with a prime divisor , where and .

For instance, substituting

in (2.6) and one gets finite-difference equations


where is the rational function on variables and

We can directly verify the following properties of the corresponding discrete mapping.

Proposition 1

Relations (2.16 ) determine 3-point mapping

preserving the form of Hamiltonain (2.1) and Poisson bracket, i.e. from , and (2.16) will follow that  .

In order to get an iterative system of finite-difference equations we identify a part of abscissas of the intersection points with arbitrary numbers

In this case finite-difference equations (2.6) implicitly depend on the independent variable via parameters of discretization . For instance, addition of prime divisors


determines the following iterative system of 2-point invertible mappings


where is given by (2.7) and

Here are arbitrary numbers, whereas the corresponding ordinates

are the functions on the phase space . We have to use this fact to calculate Poisson bracket between variables and (2.17), obtained from variables and .

Proposition 2

Relations (2.17) determine iterative system of 2-point invertible mappings

preserving the form of Hamilton function


and Poisson bracket, i.e. from and (2.17) will follow that .

The proof is a straightforward calculation.


in (2.8) and (2.12) one gets two other iterative systems of 2-point mappings

associated with multiplication of prime divisor on integer at . For the cubic Duffing oscillator at we present these mapping explicitly

Proposition 3

Relations (2.19) define two iterative systems of 2-point maps

which can be considered as the counterparts of usual geometric progression. These 2-points maps are canonical transformations of valence preserving the form of Hamiltonian (2.18), i.e. from and (2.19) will follow that .

The proof is a straightforward calculation.

Let us now take intersection divisor

see Fig.2a. If we identify coordinates of all the intersection points with the partial solutions of the Hamilton equations, relations (2.9) and define 4-point mapping

which has the standard properties.

Proposition 4

Discrete map

where and are given by (2.9,2.10), preserves the form of Hamiltonian and original Poisson bracket.

The proof is a straightforward calculation.

In order to get iterative systems of finite-difference equations we identify one of the intersection points with parameter of discretization. For instance, we can substitute