Abstract

We propose a peculiar fractional oscillator. By assuming that the motion takes place in a complex media where the level of fractionality is low, we find that the time rate of change of the energy of this system has an oscillatory behavior.

1. Introduction

In complex media such as glasses, liquid crystals, polymers, and biopolymers, the dynamical variable of interest often obeys fractional differential equations [16].

For instance, the mean squared displacement of a Brownian particle is given by ; this linear dependence on time is referred to as normal diffusion. In complex media this kind of behavior is often violated, leading to anomalous diffusion. For a subdiffusive process , with . For this process, fractional dynamic equations emerge naturally in the physical concept of continuous time random walks [6, 7].

Fractional differential equations have many applications in applied science and engineering [811]. Fractional differential equations have been investigated in pure sciences, such as pure mathematics [12].

As a fractional generalization of the oscillation phenomena, one can consider [13]

The case with corresponds to attenuated oscillation phenomenon [1419]. With the initial conditions and , the solution is , the so-called Mittag-Leffler function.

The case with corresponds to amplified oscillation phenomenon. With the initial conditions , , and , the solution is . In [14] it has been discarded on the ground that it signifies the instability of the system. In a recent study [11] it has been proven that the Mittag-Leffler function of this order is a Nussbaum function. Hence it may have applications in the control theory of electrical engineering. In another study [12] general equations of the type have been considered; however in this work the emphasis is on the existence and uniqueness of solutions.

There are various representations of the Mittag-Leffler function:(i)series representation [13],(ii)integral representation [14],(iii)approximate representation for a medium with low level of fractionality [17, 19, 20].

In this work we use the approximate representation to probe fractional differential equations of order .

The plan of this paper is as follows.

In Section 2 we describe this approximate representation. In Section 3 we obtain the solution for (1) and we compare it with the exact results of [11, 14].

In Section 4 we propose a new fractional oscillator of the form where . We obtain the solution for this damped oscillator, and we discuss the time rate of change of this oscillator. Finally in Section 5 we present our conclusions.

2. The Perturbation Scheme

There exists a multitude of definitions [15] like Riemann-Lioville, Weyl, Reisz, and Caputo for the fractional derivative. However, the fractional operator is uniquely defined. For instance, Riemann-Liouville and Caputo formulations emerge by choice (how to include the initial values) and are fully equivalent.

In this paper we only use the Caputo fractional derivative as it is easier to apply the initial conditions in this type. The left(forward) Caputo fractional derivative for is defined by where is an integer number and denotes the th derivative of the function .

To construct a solution for a process described by an equation with Caputo fractional derivative, one needs the initial conditions that can be written as

In this work we seek causal solution. Hence we require the condition for [13]. At the time the solution is determined by initial condition equation (5). And for it is obtained from the fractional differential equation describing the process under consideration. Now if the order of the fractional derivative is close to a positive integer, namely, , with small positive [17], then we will have

As noted in [17] this representation is not convenient in the limit , since

However, in this case from integration by parts we find

Now from the expansion where and is the Euler constant. Hence we obtain where is

In order to get a good approximation from (6) we find the following conditions:(i) , (ii) .

Next we generalize this expansion for the case, where , with small positive . We note that in this case ; therefore So

By expansion in terms of the parameter we get

But the process described by fractional derivative has one extra initial condition more than the process described by , namely, the condition . However, if we assume , then our expansion and hence the solutions will have a continuous dependence on the parameter in transition from to . In other words if we set , then we get

Now we consider (1) with and the initial conditions which are specified by

The perturbation expansion of the dynamical variable is given by where is the solution of the integer case, and it is described by and it is subject to

The term is a correction term, and it designates the deviation from the integer case. By substituting (15), (17) in (1) we obtain subject to

By assuming these set of initial conditions we obtain the correct result in the limit of .

3. Amplified Oscillations

The phenomena of amplified oscillations correspond to the case , in (1). It has been discussed in [11, 12, 14].

In control theory, a branch of electrical engineering, one utilizes a closed feedback loop to guarantee the stability of the system. In [11] it is proved that the solution of (1) is a Nussbaum function. Hence when it is used as a part of a general system, it can guarantee the asymptotic stability of this larger system. But this discussion is restricted to a specific field of applied science, namely, control theory.

In [12] the general properties of (2) such as uniqueness and existence have been investigated. But it does not consider the behavior of the solution.

In [14] (1) is briefly discussed; however it was discarded on the grounds that it signals the instability of the system under considerations.

Our strategy to overcome this difficulty is(1)to consider the media with low-level fractionality,(2)to add a damping term to remove this instability.

In this section we investigate the limit, where , so

For the normalized initial conditions , , and one obtains [19] where

Here and are sine and cosine integral functions, respectively [21],

For this process we define the momenta by

From (26) we find [19] where

Hence the total energy of this oscillator is

It is possible to calculate the time rate of change of the energy. The result is

Therefore, the time rate of change of the energy is always positive. Hence the total energy of the system is a monotonous increasing function of time, which means that the system absorbs energy from the environment. For instance, if we take , then we have about one percent increase in the value of in a year. For we have about one percent increase in the value of in a day. These values are chosen in a way that the condition is satisfied [17].

Thus, our study suggests an absorption interpretation for fractional derivative of order for the solution of (1), to be contrasted with the common dissipation interpretation for fractional derivative of order .

3.1. Comparison with the Exact Result

With our choice of the initial condition, the exact result is where the one parameter Mittag-Leffler function is where

The function is given by

Using and the condition , we obtain

By substituting (36) into (33) we recover the result expressed in (23).

4. An Especial Fractional Oscillator

In this section we investigate the damped fractional oscillator described by (3). For simplicity we consider a case where the damping coefficient is . In the media with low-level fractionality . By inserting (17) in (3) we get

With the initial conditions and , we will have The perturbation term satisfies with the initial conditions and . By assuming these initial conditions we obtain the correct result in the limit of . From (39) and after some calculations we obtain

For this process we define the momenta by

From (41) we find where

Hence the total energy of this oscillator is

It is possible to calculate the time rate of change of the energy. The result is

The first term and the second term in the right hand side of (45) are due to the first and the second term of (3), respectively. The first term is positive and denotes absorption of energy from the environment, and the second term is always negative and shows the dissipation of energy from the system to the environment. For large values of we have . Hence from (45)

We see that the time average of time rate of change of energy of this system is zero; namely,

5. Conclusions

Previous studies of fractional oscillation were limited to the case of ; see [1319]. The solutions portray damped motion. At the limit of we have the algebraic decay of the solutions. These solutions are of relevance in vibration of mechanical systems or oscillations of electrical network where the energy is dissipated in the form of heat. Hence in this domain of the values of the parameter , fractional derivative is a convenient tool for modeling of damping.

The results of Section 3 indicate that we may use fractional derivative to model amplification process as well.

In Section 4 we introduced a damped oscillator for a case where the order of fractional derivative is . The model presented in Section 4 has not been considered for general values of . It will be of interest to go beyond the small epsilon expansion. We conjecture that analytical solution for this model exists. We plan to report on these and other related issues in the future.

Acknowledgment

The author is grateful to the reviewers for their useful suggestions and comments.