Abstract

A fractional generalization of the gradient system to the Birkhoff mechanics, that is, a new fractional gradient representation of the Birkhoff system is investigated, in this paper. The definitions of the fractional gradient system are generalized to Birkhoff mechanics. Based on the definition, a general condition that a Birkhoff system can be a fractional gradient system is derived, the former studies for fractional gradient representation of the Birkhoff system are special cases of this paper. As applications of the results, the Birkhoff equations and fractional gradient expression of several classical nonlinear models are derived, such as the Hénon–Heiles equation, the Duffing oscillator model, and the Hojman–Urrutia equations. The results indicate, different from the former studies, that only gave the second-order gradient (integer order) expression for the Birkhoff system, an arbitrary fractional order gradient representation, exists for the Birkhoff system. The fractional potential function obtained from the general condition can determine the stationary states of these models in Birkhoffian expression.

1. Introduction

Fractional integrals and derivatives are increasingly important in the modeling of science and engineering problems because they are more suitable for describing complex phenomena in science and engineering. The discussion of noninteger order derivatives dates back to Leibniz and L’Hopital in 1695, but until 1974 the first monography [1] written by Oldham and Spanier about fractional integrals and derivatives was published. Since then, fractional calculus has been applied in many fields such as physics, mechanics, etc. [215].

A gradient system is an important dynamical system. The characteristics of a gradient system are suitable, especially to determine the stability of a system by using the Lyapunov function [1620]. If a mechanical system can be transformed into a gradient system, the stability of the mechanical system can be studied through the characteristics of the gradient system. Many researchers have performed the gradient representations of constrained mechanical systems such as the Lagrange system, the Hamilton system, and the Birkhoff system and studied the stability of their solutions [2127]. Tarasov [2830] generalized the Hamilton system and the Lagrange system into fractional gradient expressions. Mei et al. [31, 32] studied the fractional gradient representation of the Birkhoff system, Chen et al. [33] gave a fractional gradient representation of Poincaré equations. However, their studies only provide a special condition with second order and only can provide a second order gradient system, which limits its application in constrained mechanical systems. In this paper, our goal is to generalize the fractional gradient representation of the Birkhoff system to any order.

The Birkhoff system is an important dynamical system; it is an extension of the Hamilton system. Birkhoff mechanics has important applications in hadronic physics, statistic mechanics, space mechanics, biophysics, and engineering [34]. Birkhoff’s mechanical system is a more general, constrained mechanical system [3537]. Conservative and nonconservative systems and holonomic and nonholonomic systems can all be expressed by Birkhoffian formulations, which have a generalized symplectic structure [38]. So, generalizing the fractional gradient system into the Birkhoff system has general significance. Since former studies about the fractional gradient representation of the Birkhoff system only gave its second-order expression. In this paper, we will show the order of the fractional gradient Birkhoff system is not only with but also can be other orders (integer or noninteger order), and give a general condition for the Birkhoff system transforming into a fractional gradient Birkhoff system.

The structure of this paper is as follows: In Section 2, we will fix the notations of fractional derivative and fractional exterior derivative and give a brief review of the definition of the fractional gradient system and its formulation. In Section 3, we will discuss under what conditions a Birkhoff system can transform into a fractional gradient system. In Section 4, we will apply the results that we got in Section 3 to several examples, such as the Hénon–Heiles equation, the Duffing oscillator model, and the Hojman–Urrutia equations. Section 5 is the conclusion.

2. Review of the Fractional Gradient System

Fractional derivatives have different definitions. In this paper, we will adopt the Caputo and Riemann–Liouville fractional derivative definitions. Based on the definition of the Riemann–Liouville fractional derivative, we have the definition of in the Riemann–Liouville fractional derivative sense.where is the first whole number greater than . The initial point of a fractional derivative is set to be zero. The Riemann–Liouville fractional derivative has some disadvantages in physical applications, such as the hyper-singular improper integral, where the order of singularity is higher than the dimension, and the nonzero of the fractional derivative of constants, which would make dissipation not vanish for a system in equilibrium. The Caputo fractional derivatives can choose the same initial value as the integer derivatives, which are more suitable for applying to physical problems. The Caputo fractional derivative has the following definition:where

When , and , we have

We can introduce fractional exterior derivatives as [2931, 39]

For the convenience of expression, we use to denote order fractional derivative, and only point out whether it is the Riemann–Liouville fractional derivative or the Caputo fractional derivative as calculation.

The differential equations of a gradient system have the form aswhere , and is a potential function.

Following, we will review the definition and conditions of the fractional gradient system proposed by Tarasov [2830].

Definition 1. For a dynamical system,If a fractional differential 1-form,This is an exact form, that is, , where is a continuous differential function, then system (7) is a fractional gradient system. Thus, we have .

It is known that to be exact, it is a sufficient condition to be closed; however, the converse statement is not necessarily true. From the Poincaré theorem that is verified to be true for fractional exterior derivatives, it can be deduced that any smooth 1-form (8) that is closed is also exact on a contractible open subset of [28]. We have the following proposition.

Proposition 1. If smooth functions on a contractible open subset of satisfy the relationsthen dynamical system (7) is a fractional gradient system with formwhere is a continuous differential function and . Because the Riemann–Liouville fractional derivative of a constant is not necessary to be zero, so, cannot define the stationary state of a fractional system (10), where is constant. The stationary state of the fractional gradient system with the Riemann–Liouville fractional derivative can be defined bywhere are constants and is the first whole number greater than or equal to .
If the fractional exact 1-form 8 equals to zero that is in the Caputo derivative, we can get , which defines the stationary state of the fractional gradient system (10), where is constant.

3. Fractional Gradient Representation of the Birkhoff System

An automatic Birkhoff system has formwhere is a Birkhoffian, andis called the Birkhoffian tensor and are Birkhoff’s functions. If the system is nonsingularity, that is , the automatic Birkhoff system (12) can be expressed aswhere

Definition 2. For a Birkhoff system (12) or (14), if the fractional differential 1-formis an exact form, that is , then it is a fractional gradient expression of the Birkhoff system.
Make exterior derivative to fractional differential 1-form (16), we haveIf is closed, we haveSo, we have the following proposition.

Proposition 2. If equation (14) on a contractible open subset of satisfies condition (18), it is an order fractional gradient system.

Remark 1. When , condition (18) degenerates to an integer derivative case, which gives the condition for a Birkhoff system transforming into a classical gradient expression [21].

When (including integer and noninteger numbers), we call it an order fractional gradient expression of the Birkhoff system.

Remark 2. When , condition (18) gives the condition for a Birkhoff system transforming into a second-order fractional gradient system [3133].

We point out that cannot be a fractional gradient system for the Birkhoff system, it may be a fractional gradient system for the fractional Birkhoff system, which question we will discuss in future works. From (10), we have

According to Birkhoff 1-form (19), we can get the expression of the potential function .where is the fractional Riemann–Liouville integral symbol with definition

Remark 3. Here can be any order, but we should point out that only (20) makes the Birkhoff 1-form is exact, and (19) is a fractional order gradient expression of the Birkhoff system.

4. Results

4.1. Application A: Hénon-Heiles Equation as a Birkhoff System and Fractional Gradient System

In this section, we will give the Birkhoff expression of the Hénon–Heiles equation and prove it can be a fractional gradient system with the Riemann–Liouville or the Caputo fractional derivatives, respectively, in the Birkhoff mechanics frame.

The Hénon–Heiles equation is a sort of nonlinear, nonintegrable Hamilton system, which has broad applications in physics and applied mathematics [10, 40, 41]. The Hénon–Heiles equation’s form is

Let Birkhoff variables be ; a Birkhoff expression of H non-Heiles problem is

Therefore, we can get

The Birkhoff equations can be expressed as

From (16), we can get a fractional differential 1-form of Birkhoff (26).

Its fractional exterior derivative is

Next, we will discuss its fractional gradient expression by the Caputo and Riemann–Liouville fractional derivatives, respectively.

4.2. Fractional Gradient Representation of the Birkhoff Equation (27) with the Caputo Fractional Derivative

Firstly, we study its fractional gradient system using the Caputo fractional derivative. As we take the Caputo fractional derivative, that is, in (28), when , (29) becomes

So , then (26) cannot be expressed as a fractional gradient system.

When , we can get

So , then (26) cannot be expressed a fractional gradient system.

When with , , so (26) can be expressed as an order Caputo type fractional gradient system, where . From (19), we can get

According to (20), we can get the potential function aswhere . We can directly calculate that is not the Birkhoff system (27). In order to make the potential function V (33) is the potential function of the Birkhoff system (27), we have to modify (33) as

Here, is the first integer greater than or equal to . However, it does not guarantee the system is a fractional gradient system, so considering , we have , the Birkhoff system (27) can be expressed as an order Caputo fractional gradient system.

The stationary states of the Birkhoff systems are defined by equationwhere is constant.

4.3. Fractional Gradient Representation of the Birkhoff Equation (27) with the Riemann–Liouville Fractional Derivative

Secondly, we study the fractional gradient representation of the Birkhoff system (27) using the Riemann–Liouville fractional derivative. When we take the Riemann–Liouville fractional derivative that is , (29) becomes

When , we have

So, Birkhoff system (27) is not a gradient system.

When , we have

So (27) is not an order fractional gradient system.

When , , so (27) is an order fractional gradient system.

From (19), we can get

The potential function can be written aswith .

4.4. Conclusion

In this section, we give the Birkhoff expression of the Hénon-Heiles equation and give its fractional gradient representation using the Caputo fractional derivative and the Riemann–Liouville fractional derivative, respectively. We can conclude that from this example, the Hénon–Heiles equation in the Birkhoff mechanics frame can be expressed as an (noninteger) order fractional gradient system using the Caputo fractional derivative, or as a third-order fractional gradient system using the Riemann–Liouville fractional derivative. However, it cannot be a second-order gradient system, so our results are more general than those in [3133]. We also conclude that a fractional gradient system is not a classical gradient system in general. From the calculus, we find that using the Riemann–Liouville fractional derivative can only give an integer order gradient system, but using the Caputo fractional derivative can give an integer or noninteger order gradient system, so the Caputo fractional derivative is more suitable to study fractional gradient systems than the Riemann–Lioville fractional derivative.

5. Application B: The Birkhoff Equation of Duffing Oscillator Models and Its Fractional Gradient System

In this section, we will give the Birkhoff equation expression of the Duffing oscillator models and prove it can be a fractional gradient system with the Riemann–Liouville or the Caputo fractional derivatives, respectively, in the Birkhoff mechanics frame.

Duffing oscillator models and their equation are widely discussed in nonlinear science: where is the nonlinear stiffness coefficient of the Duffing oscillator. Take Birkhoff’s variables as

We can get Birkhoff expression of the Duffing oscillator equationwith Birkhoff equations

From (16), we can get a fractional differential 1-form of Birkhoff equation (43).

Its fractional exterior derivative is

Next, we will discuss its fractional gradient expression using the Caputo and the Riemann–Liouville fractional derivatives, respectively.

5.1. Fractional Gradient Representation of Birkhoff System (43) with the Caputo Fractional Derivative

Firstly, we study its fractional gradient system using the Caputo fractional derivative. As we take the Caputo fractional derivative, that is, in (43), When , (45) becomes

Then , so the Duffing oscillator system expressed by the Birkhoff equation (43) is not a fractional gradient system.

When , we can get

Then , so the Duffing oscillator system expressed by the Birkhoff equation (43) is not a fractional gradient system.

When , we can get

Then , so the Duffing oscillator system expressed by the Birkhoff equation (43) is not a fractional gradient system.

When with , , so equation (43) can be expressed as an order Caputo type fractional gradient system, where . From (19), we can get

According to (20) and condition (18), we can get the potential function as

The stationary state of system 43 in the Caputo derivative is defined by equation (35).

5.2. Fractional Gradient Representation of Birkhoff System (43) with the Riemann–Liouville Fractional Derivative

Secondly, we study the fractional gradient representation of the Birkhoff system (43) using the Riemann–Liouville fractional derivative. When we take the Riemann–Liouville fractional derivative, that is , (45) becomes

When , we have

Then , so (43) is not a fractional gradient system.

When , we have

Then , so (43) is not a fractional gradient system.

When , we have

Then , so (43) is not a fractional gradient system.

When , , so (43) is an order fractional gradient system. From equation (19), we can get

The potential function is

5.3. Conclusion

In this section, we give the Birkhoff equation expression of the Duffing oscillator problem and give its fractional gradient representation using the Caputo fractional derivative and the Riemann–Liouville fractional derivative, respectively. We conclude that, from this example, the Duffing oscillator model in the Birkhoff mechanics frame can be expressed as a (noninteger) order fractional gradient system using the Caputo fractional derivative or a fourth order fractional gradient system using the Riemann–Liouville fractional derivative. However, it cannot be a second-order gradient system, so our results are more general than those in [3133]. We also conclude that a fractional gradient system is not a classical gradient system in general. From the calculus, we find that using the Riemann–Liouville fractional derivative can only give an integer order gradient system, but using the Caputo fractional derivative can give an integer or noninteger order gradient system, so the Caputo fractional derivative is more suitable to study fractional gradient systems than the Riemann–Liouville fractional derivative.

6. Application C: Hojman-Urrutia Equations as a Birkhoff System and Fractional Gradient System

In this section, we will give the Birkhoff equation expression of the Hojman–Urrutia equations, and prove it can be a fractional gradient system with the Riemann–Liouville or Caputo fractional derivatives, respectively, in the Birkhoff mechanics frame.

The Hojman–Urrutia equations have an important status in the development of inverse problems in Lagrangian mechanics and Birkhoffian mechanics [36]. Their equations are written in the following form:

Let ; we can get one types of its Birkhoff expression.

So we can getandwith Birkhoff equations

From (16), we can get a fractional differential 1-form of Birkhoff (57).

Its fractional exterior derivative is

6.1. Fractional Gradient Representation of Birkhoff Systme (60) with the Caputo Fractional Derivative

Next, we will discuss its fractional gradient expression by the Caputo and the Riemann–Liouville fractional derivative, respectively.

Firstly, we study its fractional gradient system using the Caputo fractional derivative. As we take the Caputo fractional derivative, that is, in (59), When , (62) becomes

Then, , so the Hojman–Urrutia equation expressed by Birkhoff (57) is not a fractional gradient system.

When with , we can get , so (57) can be expressed as an order Caputo type fractional gradient system, where . From (19), we can get

According to (20) and condition (18), we can get the potential function as

Which can define the stationary state of system (60) in the sense of the Caputo fractional derivatives in equation (35).

6.2. Fractional Gradient Representation of Birkhoff System (60) with the Riemann–Liouville Fractional Derivative

Secondly, we study fractional gradient representation of the Birkhoff system (60) using the Riemann–Liouville fractional derivative. When we take the Riemann–Liouville fractional derivative that is , (62) becomes

When , we get

Then , so (60) is not a fractional gradient system.

When , , so (60) is an order fractional gradient system with potential function with . From equation (19), we can get

The potential function is

6.3. Conclusion

In this section, we give the Birkhoff expression of the Hojman–Urrutia equation and give its fractional gradient representation using the Caputo fractional derivative and the Riemann–Liouville fractional derivative, respectively. We can conclude that from this example, the Hojman–Urrutia equation in the Birkhoff mechanics frame can be expressed as an (noninteger) order fractional gradient system using the Caputo fractional derivative or a second-order fractional gradient system using the Riemann–Liouville fractional derivative, so our results are more general than those in [3133]. We also conclude that a fractional gradient system is not a classical gradient system in general. From the calculus, we find that using the Riemann–Liouville fractional derivative can only give an integer order gradient system, but using the Caputo fractional derivative can give an integer or noninteger order gradient system, so the Caputo fractional derivative is more suitable to study fractional gradient systems than the Riemann–Liouville fractional derivative.

7. Conclusion

In this paper, we give a definition and a general condition of a Birkhoff system to be a fractional gradient system with the Caputo fractional derivative and the Riemann–Liouville fractional derivative, respectively. When , the definition and condition degenerate into those of the classical gradient system, and when , our results goes back to those of the second order gradient representation [3133], so the classical gradient system and a second-order gradient system can be considered as special cases of the fractional gradient representation of the Birkhoff system obtained in this paper.

As applications, we apply the results we get in this paper to the Hénon–Heiles equation, the Duffing oscillator model, and the Hojman–Urrutia equations. The Birkhoff equation expression and fractional gradient expression of these classical equations are given. It shows that the orders of the fractional gradient Birkhoff system are not only with but can also be other integers or noninteger orders in this paper, so our results are more general.

The limitation of this paper is that we only discuss the fractional gradient representation of an autonomous Birkhoff system. The more general cases of the fractional gradient representation of nonautonomous Birkhoff systems and the fractional gradient representation of the fractional Birkhoff system [42] will be considered in future work.

Abbreviations

:Symbol of Riemann–Liouville fractional derivative
:Symbol of Caputo fractional derivative
:Riemann–Liouville fractional integral symbol
:Symbol of Gamma function
:Symbol of fractional exterior derivatives
:Symbol of Birkhoffian
:Symbol of Birkhoffian tensor
:Symbols of Birkhoff’s functions.

Data Availability

All data generated or analyzed during this study are included in this article.

Conflicts of Interest

The author declares that there are no conflicts of interest.

Acknowledgments

This work was partly supported by the National Natural Science Foundation of China (Grant Nos. 12272148, 11772141).