Dynamics of the Morse Oscillator: Analytical Expressions for Trajectories, Action-Angle Variables, and Chaotic Dynamics

TL;DR

This paper derives explicit analytical expressions for the trajectories, action-angle variables, and chaotic dynamics of the Morse oscillator, enhancing understanding of its behavior under perturbations.

Contribution

It provides new explicit formulas for Morse oscillator trajectories and establishes conditions for chaos using a Melnikov approach.

Findings

Explicit formulas for trajectories and periods

Conditions for chaotic dynamics under perturbation

Application of Melnikov method to Morse oscillator

Abstract

We consider the one degree-of-freedom Hamiltonian system defined by the Morse potential energy function (the "Morse oscillator"). We use the geometry of the level sets to construct explicit expressions for the trajectories as a function of time, their period for the bounded trajectories, and action-angle variables. We use these trajectories to prove sufficient conditions for chaotic dynamics, in the sense of Smale horseshoes, for the time-periodically perturbed Morse oscillator using a Melnikov type approach.

Equations

$$V(q)=D{\left(1-e^{-\alpha q}\right)}^2,$$

$$H(q,p)=\frac{p^2}{2m}+D{\left(1-e^{-\alpha q}\right)}^2,(q,p)\in {\mathbb{R}}^2,$$

$$\dot{q}=\frac{\partial H}{\partial p}$$

$$\dot{p}=-\frac{\partial H}{\partial q}$$

$$\left\{(q,p)\in {\mathbb{R}}^2|H(q,p)=h=\text{mbox}constant\right\}.$$

$$(q,p)=(\infty ,0),(0,0).$$

$$J=\left(\begin{array}[]{cc}0&\frac{1}{m}\\ -2D\alpha^{2}\left(-e^{-\alpha q}+2e^{-2\alpha q}\right)&0\end{array}\right).$$

$$\pm \sqrt{\det \mathrm{J}}.$$

$$(q,p)=(0,0)\Rightarrow \det \mathrm{J}=-\frac{2D{\alpha }^2}{\mathrm{m}},$$

$$\pm i\sqrt{\frac{2D}{m}}\alpha ,$$

$$(q,p)=(\infty ,0)\Rightarrow \det \mathrm{J}=0,$$

$$(q,p)=(\infty ,0)\Rightarrow H=D>0.$$

$$(q,p)=(0,0)\Rightarrow H=0.$$

$$h=\frac{p^2}{2m}+D{\left(1-e^{-\alpha q}\right)}^2,0<h<D,$$

$$\frac{p^2}{2m}=h-D{\left(1-e^{-\alpha q}\right)}^2.$$

$$h=D{\left(1-e^{-\alpha q}\right)}^2$$

$$0\le \sqrt{\frac{h}{D}}\le 1,$$

$$1\le 1+\sqrt{\frac{h}{D}}\le 2,$$

$$0\le 1-\sqrt{\frac{h}{D}}\le 1,$$

$$1-e^{-\alpha q}=\sqrt{\frac{h}{D}}.$$

$$q_{+}=-\frac{1}{\alpha }\log \left(1-\sqrt{\frac{h}{D}}\right)>0.$$

$$1-e^{-\alpha q}=-\sqrt{\frac{h}{D}},$$

$$q_{-}=-\frac{1}{\alpha }\log \left(1+\sqrt{\frac{h}{D}}\right)<0.$$

$$D=\frac{p^2}{2m}+D{\left(1-e^{-\alpha q}\right)}^2,$$

$$\frac{dq}{dt}=\pm \sqrt{\frac{2}{m}}\sqrt{h-D{\left(1-e^{-\alpha q}\right)}^2},$$

$$\frac{dq}{\sqrt{h-D{\left(1-e^{-\alpha q}\right)}^2}}=\pm \sqrt{\frac{2}{m}}dt.$$

$$T(h)$$

$$u=e^{-\alpha q}.$$

$$T(h)=\frac{\pi \sqrt{2m}}{\alpha \sqrt{D-h}}$$

$${\int }_{q_{+}}^q\frac{d{q}^{\prime }}{\sqrt{h-D{\left(1-e^{-\alpha q^{\prime }}\right)}^2}}=\sqrt{\frac{2}{m}}t.$$

$$q(t)=\frac{1}{\alpha }\log \frac{\sqrt{Dh}\cos \left(\sqrt{\frac{2(D-h)}{m}}\alpha t\right)+D}{D-h}.$$

Full text

*Keywords: * Morse oscillator, action-angle variables, homoclinic orbit, Melnikov function, chaos.

1 Introduction

The potential energy function derived by P. M. Morse is truly a ”workhorse” potential energy function in theoretical chemistry [Mor29]. Originally it was devised to describe the intermolecular force between the two atoms in a diatomic molecule. It has the functional form:

[TABLE]

where $q$ represents the distance between the two atoms, $D>0$ represents the depth of the potential well (defined relative to the dissociated atoms), and $\alpha>0$ controls the width of the potential well ($\alpha$ small corresponds to a ”wide” well, $\alpha$ large corresponds to a narrow well).

The Morse potential defines a one degree-of-freedom Hamiltonian system, i.e. the phase space is two dimensional described by coordinates $(q,p)$, where $p$ is the momentum conjugate to the position variable $q$. The Hamiltonian has the form of the sum of the kinetic energy and the potential energy (the Morse potential). The Hamiltonian system is integrable, and all trajectories lie on the level sets of the Hamiltonian function. The level sets can be used to derive integrals that give (time) parametrizations of trajectories. Regions of closed (bounded) trajectories can be used to construct special coordinates–action-angle coordinates, where the angle denotes a particular location on the closed level set and the action is the area enclosed by a closed level set (divided by $2\pi$). The transformation to action-angle coordinates for integrable Hamiltonan systems is a standard topic in good classical mechanics textbooks, see, e.g. [LL60, Arn13]. The transformation preserves the Hamiltonian nature of the system, i.e. it is a canonical transformation, and therefore the standard approach to constructing such transformations is through the use of generating functions. However, for one degree-of-freedom time-independent Hamiltonian systems (such as the one described by the Morse potential) there is a simpler approach to generating the action-angle transformation that uses the geometry of the closed level set of the Hamiltonian function and the explicit (time) parametrization of the trajectories that can be obtained (in principle, if the necessary integrals can be explicitly computed) for one degree-of-freedom Hamiltonian systems. The approach is inspired by the seminal paper of Melnikov [Mel63]. This approach was developed in detail in [Wig90] (the 1990 edition, not the 2003 edition) and is also described in [MW94]. This is the approach that we will follow here.

Action-angle variables are important in Hamiltonian mechanics for a number of reasons. From the point of view of classical mechanics they are the coordinate system used for the development of the Kolmogorov-Arnold-Moser and Nekhoroshev theorems [Dum14]. They also play a central role in the quantization of classical Hamiltonian systems (in fact, the action and the constant $\hbar$ have the same units) [Sto05, Kel58, KR60, Kel85].

This paper is outlined as follows. In Section 2 we describe the Hamiltonian system described by the Morse potential (1). In Section 2.1 we describe the equilibria and determine their stability properties in the linear approximation. In Section 3 we discuss the geometry of the region of bounded motion, i.e. the region of periodic orbits bounded by a homoclinic orbit (”separatrix”). In Section 4 we compute the period of the periodic orbits and show how it depends on the energy and other system parameters. In Section 5 we compute explicit expressions for the trajectories in the region of bounded motion. In Section 5.1 we compute an explicit expression for the homoclinic orbit. In Section 6 we compute the action-angle variables. In Section 7 we compute explicit expressions for the trajectories in the region of unbounded motion. In section 8 we apply our results to exploring the nature of chaotic dynamics in the time periodically forced Morse oscillator and in Section 9 we discuss our conclusions.

2 The Hamiltonian

The dynamical system defined by the Morse potential is Hamiltonian, with Hamiltonian function given by:

[TABLE]

and Hamilton’s equations defined by:

[TABLE]

The level sets of the Hamiltonian function have the form:

[TABLE]

They are (in general) one dimensional curves that are invariant under the Hamiltonian dynamics, i.e. the Hamiltonian vector field is tangent to the level sets. Since the trajectories lie on these curves, we can use the form of the level curves to obtain parametrizations of certain trajectories, as we will demonstrate.

2.1 Equilibria and their Linearized Stability

It is straightforward to verify the the following two points are equilibria for Hamilton’s equations:

[TABLE]

Next we check their linearized stability properties. The Jacobian matrix, denoted $J$, of Hamilton’s equation is given by:

[TABLE]

The eigenvalues of $J$ are given by:

[TABLE]

Hence for the two equilibria we have:

[TABLE]

with corresponding eigenvalues:

[TABLE]

and

[TABLE]

where both eigenvalues are zero.

The equilibrium $(q,p)=(0,0)$ is stable (”elliptic” in the Hamiltonian dynamics terminology) and the $(q,p)=(\infty,0)$ is unstable (a ”parabolic” saddle point in the Hamiltonian dynamics terminology).

3 The Region of Bounded Motion: Periodic Orbits

Using the Hamiltonian (2) it is straightforward to verify that the equilibria have the following energies:

[TABLE]

and

[TABLE]

Trajectories with energies larger than $D$ have unbounded motion in $q$. Trajectories having energies $h$ satisfying $0<h<D$ correspond to periodic motions. The level sets of these periodic orbits are given by:

[TABLE]

and surround the stable equilibrium point $(q,p)=(0,0)$ as shown in figure 1. The periodic orbits intersect the $q$ axis at two distinct points, $q_{+}>0$ and $q_{-}<0$, which are referred to as turning points. These turning points are computed as follows.

Rewriting (13) gives:

[TABLE]

The turning points are obtained from (14) by setting $p=0$:

[TABLE]

Note that we have:

[TABLE]

from which we obtain the following relations:

[TABLE]

Taking the positive root of (15) gives:

[TABLE]

from which we obtain the positive turning point:

[TABLE]

Taking the negative root of (15) gives:

[TABLE]

from which we obtain the negative turning point:

[TABLE]

The level curve with energy equal to the dissociation energy $h=D$ has the form:

[TABLE]

and is a separatrix connecting the (parabolic) saddle point. In the terminology of Hamiltonian dynamics it is a homoclinic orbit. It separates bounded from unbounded motion as illustrated in figure 1.

4 Calculation of the Period of a Periodic Orbit

In this section we calculate the period of the periodic orbits. From (3) we have $\dot{q}=\frac{p}{m}$. Using this expression, and the expression for the level set of the Hamiltonian defining a periodic orbit given in (13), we have

[TABLE]

or

[TABLE]

We denote the period of a periodic orbit corresponding to the level set with energy value $h$ by $T(h)$. We can obtain the period by integrating $dt$ around this level set. Using (25), this becomes:

[TABLE]

Computation of this integral is facilitated by the substitution:

[TABLE]

After computing the integral using integral 2.266 in [GR80] we obtain:

[TABLE]

There are two limits in which this expression can be checked with respect to previously obtained results. First, we note that for $h=D$ (i.e. the energy of the homoclinic orbit, or ”separatrix”) $T(D)=\infty$, which is what we expect for the ”period” of a separatrix.

Second, we consider $h=0$, which is the energy of the elliptic equilibrium point. In this case we have $T(0)=\frac{\pi\sqrt{2m}}{\alpha\sqrt{D}}$, which is $2\pi$ divided by the imaginary part of the magnitude of eigenvalue of the Jacobian evaluated at the stable equilibrium point, as we expect.

5 Expressions for $q(t)$ and $p(t)$, $0<h<D$

In this section we derive an expression for $q(t)$. Differentiating the expression for $q(t)$ will give the expression for $p(t)$ through the relation $\dot{q}=\frac{p}{m}$. Using (25), we have

[TABLE]

Choosing the lower limit if the integral to be $q_{+}$ is arbitrary, but it is equivalent to the choice of an initial condition. After computing this integral, we obtain:

[TABLE]

It is straightforward to check that the period of (30) is (28).

5.1 The Homoclinic Orbit

As noted earlier, the homoclinic orbit, corresponding to $h=D$, is given by the level set (23). Hence, the integral expression for the homoclinic orbits is obtained from (25) by setting $h=D$. Computing the integral gives:

[TABLE]

It is a simple matter to check that $\lim_{t\rightarrow\pm\infty}q(t)=\infty$. Subsequently we obtain $p_{0}(t)$ from $\dot{q}=\frac{p}{m}$ as

[TABLE]

6 Expressions for Action and the Angle, $0<h<D$

In this section we compute the action-angle representation of the orbits in the bounded region following [Mel63, Wig90, MW94].

We consider a level set defined by the Hamiltonian (4), for $0<h<D$. i.e. we consider a periodic orbit with period $T(h)$. Choosing an arbitrary reference point on the periodic orbit, the angular displacement of a trajectory starting from this reference position after time $t$ is given by:

[TABLE]

Using the substitution (27) and integral 2.266 in [GR80], we have:

[TABLE]

The action associated with this periodic orbit is the area that it encloses (in phase space) divided by $2\pi$:

[TABLE]

Recalling (25)

[TABLE]

we obtain:

[TABLE]

and therefore

[TABLE]

Using the substitution (27) and integral 2.267 in [GR80], we have:

[TABLE]

7 Expressions for $q(t)$ and $p(t)$, $0<D<h$

Increasing the total energy $h$ to the value of $D$ and beyond results in unbounded motion. Trajectories retain the turning point $q_{-}$, while $q_{+}$ becomes infinite. The expression for $q_{-}$ is identical to low energies and is obtained from (2) by setting $p=0$ and solving for $q$. Recall from (22) that $q_{-}=-\frac{1}{\alpha}\log\left(1+\sqrt{\frac{h}{D}}\right)$.

For unbounded trajectories it is not possible to define a (finite) period, but we can obtain an expression for $t$ as a function of position. It can be derived by integrating (25) from $q_{-}$ to an arbitrary position $q$ by using the substitution (27) as follows:

[TABLE]

We obtain an explicit solution $q=q(t)$ by inverting (40).

[TABLE]

where $\beta=\alpha\sqrt{\frac{2(h-D)}{m}}$. Differentiating (41) and using the relation $\dot{q}=\frac{p}{m}$ yields the expression of $p(t)$.

8 Application: Chaos in the Periodically Forced Morse Oscillator

Once explicit solutions are known, it is insightful to consider how they change when the vector field is changed, slightly. This is considered in the field of perturbation theory. The subharmonic and homoclinic Melnikov [Wig90] are examples of global perturbation methods that consider the effect of perturbations on the entire unperturbed integrable system, such as the Morse oscillator that we have considered. The type of perturbation that we will consider is a time periodic excitation of the Morse oscillator. The subharmonic and homoclinic Melnikov methods for this system have been previously considered in [GHZ03], however, with respect to chaotic dynamics, important technical details for this periodically perturbed Morse oscillator were not considered. We will describe these in detail here, as well as consider the nature of the effect of the parameters of the Morse potential, and the periodic excitation, on chaotic dynamics.

The homoclinic Melnikov function is a well-known and popular method for proving the existence of chaos, in the sense of Smale horseshoes, in time periodically perturbed one degree-of-freedom Hamiltonian systems. The method allows one to prove the existence of transverse homoclinic periodic orbits to a hyperbolic periodic orbit. Such orbits admit a Smale horseshoe construction. There are two issues with the application of the standard homoclinic Melnikov function to the time periodically perturbed Morse oscillator, and both issues arise from the fact that the saddle point at infinity is parabolic, not hyperbolic. One issue is the fact that the standard derivation of the homoclinic Melnikov function uses hyperbolicity of the fixed point to show that certain terms in the Melnikov function vanish. The other issue involves the construction of the Smale horseshoe map for orbits homoclinic to a parabolic point. Both issue are considered for the perturbed Morse oscillator in [BW92] where it is shown that the growth and decay properties for the parabolic point in the Morse oscillator are sufficient for the standard Melnikov function to be valid as well as the Smale horseshoe map construction to be valid. Hence, the ”Melnikov approach” is sufficient for determining the existence of chaos for the time-periodically perturbed Morse oscillator.

We consider the following time-periodic perturbation of the Morse oscillator:

[TABLE]

where $\varepsilon>0$ is the magnitude of the perturbation and $\omega>0$ the frequency of the perturbation. The unperturbed system (3) is obtained by setting $\varepsilon=0$. We formulate the homoclinic Melnikov function using expressions (31), (32) for $q_{0}(t)$, $p_{0}(t)$, the homoclinic orbit of the equilibrium point $(q,p)=(\infty,0)$, as follows:

[TABLE]

where $DH$ is the gradient of the unperturbed Hamiltonian (2), $t_{0}$ defines the point $(q_{0}(t_{0}),p_{0}(t_{0}))$ at which the Melnikov funtion $M$ is evaluated and $\phi_{0}$ is the time at which $M$ is evaluated. We solve (43) using integral 3.723 in [GR80].

[TABLE]

Clearly, this function has ”simple zeros” indicating the existence of homoclinic orbits, and the associated Smale horseshoe type chaos.

The magnitude of the Melnikov function is a measure of the intensity of the chaos. Since the exponent in the exponential term is negative,

[TABLE]

The exponent is decreased by decreasing $\omega$ and $m$, and increasing $\alpha$ and $D$. Note that $M$ is monotonically increasing in dissociation energy $D$, yet the homoclinic orbit ceases to exist for $D=\infty$.

The dominant term in the Melnikov functions is

[TABLE]

To reach the highest intensity of chaos in the system, $m$ has to be large and $\alpha$ small, whereby the mitigating effect of the exponential term can be countered by choosing $D>\frac{m}{\alpha^{2}}$. Then for any $\omega$ we can find $(t_{0},\phi_{0})$, such that $\sin(\omega t_{0}+\phi_{0})=1$ and

[TABLE]

In figure 2 we show a comparison of different intensities of chaos using Lagrangian descriptors, a method introduced by [JMM09] and proven to display invariant manifolds in two dimensional area-preserving maps by [LBWM15].

9 Conclusion

In this paper we have considered the one degree-of-freedom Hamiltonian system defined by the Morse potential energy function (the ”Morse oscillator”). We used the geometry of the level sets to construct explicit expressions for the trajectories as a function of time, their period (for the bounded trajectories, and action-angle variables. We used these trajectories to prove sufficient conditions for chaotic dynamics, in the sense of Smale horseshoes, for the time-periodcally perturbed Morse oscillator using a Melnikov type approach. Using this approach we were able to determine the influence of the parameters of the Morse potential on the chaotic dynamics.

Acknowledgments

We acknowledge the support of EPSRC Grant no. EP/P021123/1 and Office of Naval Research (Grant No. N00014-01-1-0769).