alexa Relaxation Equations: Fractional Models | OMICS International
ISSN: 2090-0902
Journal of Physical Mathematics
Make the best use of Scientific Research and information from our 700+ peer reviewed, Open Access Journals that operates with the help of 50,000+ Editorial Board Members and esteemed reviewers and 1000+ Scientific associations in Medical, Clinical, Pharmaceutical, Engineering, Technology and Management Fields.
Meet Inspiring Speakers and Experts at our 3000+ Global Conferenceseries Events with over 600+ Conferences, 1200+ Symposiums and 1200+ Workshops on
Medical, Pharma, Engineering, Science, Technology and Business

Relaxation Equations: Fractional Models

Capelas de Oliveira E* and Ester CFA Rosa

Department of Applied Mathematics, IMECC - University of Campinas, Brazil

*Corresponding Author:
Capelas de Oliveira E
Department of Applied Mathematics
IMECC - University of Campinas, Brazil
Tel: +55 19 3521-7000
E-mail: [email protected]

Received Date: September 01, 2015; Accepted Date: September 02, 2015; Published Date: September 14, 2015

Citation: Oliveira CD, Rosa CFA (2015) Relaxation Equations: Fractional Models. J Phys Math 6:146. doi:10.4172/2090-0902.1000146

Copyright: © 2015 Oliveira CD, et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Visit for more related articles at Journal of Physical Mathematics


The relaxation functions introduced empirically by Debye, Cole-Cole, Cole-Davidson and Havriliak-Negami are, each of them, solutions to their respective kinetic equations. In this work, we propose a generalization of such equations by introducing a fractional differential operator written in terms of the Riemann-Liouville fractional derivative of order γ, 0<γ≤1. In order to solve the generalized equations, the Laplace transform methodology is introduced and the corresponding solutions are then presented, in terms of Mittag-Leffler functions. In the case in which the derivative’s order is γ = 1, the traditional relaxation functions are recovered. Finally, we presented some 2D graphs of these function.


Mittag-Leffler functions; Laplace transform; Riemann- Liouville derivative; Fractional differential equations; Dielectric relaxation; Complex susceptibility; Relaxation function; Response function; Debye; Cole-Cole; Davidson-Cole; Havriliak-Negami


Fractional calculus has greatly developed during the last years and has established itself as a generalization of classical differential and integral calculus. Over the last three decades the investigation on this subject has intensified and applications were found in many fields of science such as physics, mathematics, chemistry, financial systems, economy and engineering [1]. As examples we may mention the following applications: the fractional control of engineering systems; variational calculus and optimum control for fractional dynamic systems; numerical and analytical tools and techniques; fundamental explorations of mechanical, electric and thermal constitutive relations; investigation of properties of several engineering materials such as viscoelastic polymers, foams, gels and animal tissues; diffusion phenomena; bioengineering and biomedical applications; thermal modeling for engineering systems such as breaks and tool-machines and, finally, imaging and signal processing [2-10].

The mathematical investigation of relaxation processes in dielectrics has been conducted with the use of fractional calculus tools [11-14]. The properties of dielectric materials are usually represented by the empirical susceptibility functions as they have been initially modeled by Debye (D), then by the Cole brothers (C-C) and, later, in the works of Cole-Davidson (C-D) and Havriliak-Negami (H-N).

With the help of the concept of memory developed in the formalism of Mori-Zwanzig, it is possible to deduce the kinetic equations for the relaxation processes described by those empirical functions [15]. The equations obtained are expressed in terms of the fractional derivative of Riemann-Liouville and their solutions are given in terms of Mittag- Leffler functions [16].

The Mittag-Leffler functions of a real variable t, with one, two and three parameters, are convenient to study anomalous processes such as those occurring in dielectrics, including the aforementioned classical cases C-C, C-D and H-N. It is crucial to emphasize that these functions are completely monotonic for t > 0 [11,17].

In this work we study fractional kinetic equations with a derivative of order γ, 0 < γ ≤ 1, proposed as generalizations of the classic kinetic equations associated with relaxation processes in dielectrics, in such a way that the traditional kinetic equations become particular cases of the fractional equations when γ = 1.

This work is organized as follows: in the first section a brief theoretical introduction is presented of the four empirical expressions of complex susceptibility (in the frequency domain) related to Debye’s model and to the anomalous relaxation processes C-C, C-D e H-N. On the basis of these expressions of complex dielectric permittivity, we present in the second section the construction of the kinetic equations for the Debye and the anomalous relaxation functions. The construction of these kinetic equations is carried with the introduction of an integral memory function. At the end of the second section the equations are solved and the relaxation functions (in the time domain) are written in terms of Mittag-Leffler functions; we also show their graphic representations. In the third section, the previously mentioned fractional kinetic equations are solved and the graphs of such solutions are shown. Finally, the fourth section, dedicated to the conclusions, presents a comparative analysis of the results obtained and to the discussion of outline a proposal for future study.

Classic Empirical Models

The most important phenomenon associated with a dielectric material is its polarization, which consists in the change of the distribution of its molecular and atomic charges when it is subjected to the action of an electric field. Thus, when an electric field is applied to a dialectric it produces a very small electric current called dielectric loss, and its constituent particles, ions or molecules suffer small dislocations or rearrangements, thus altering their equilibrium positions [18]. These molecular parts do not leave nor reach their state of equilibrium instantaneously: a variable amount of time is necessary for this change of positions to take place. This time lapse necessary for the material to respond to the electric field applied is called relaxation time. Thus, when submitted to an electric excitation, the dielectric will respond to this action in an attempt to reestablish its equilibrium during and after the electric stimulus. The polarization of each dielectric material depends on the nature of its molecular and atomic chemical bonds, and there is presently no universal model which can explain the polarization phenomenon in all materials [19]. Debye’s response function was the first theoretical model for the dielectric behavior of some substances [20]. However, due to its limitations, this model is incapable of describing in detail the dielectric response of a large number of solids and liquids.

Therefore, in the years following the appearance of Debye’s theory, several other response functions were proposed in the literature to serve as models for describing the dielectric relax-ation of many materials. Among them, we here focus our attention on the model proposed by the Cole brothers [21,22], and the one by Cole-Davidson [23,24], both of which emerged in attempts to adjust the response function to the experimental behavior of some dielectric materials. There is also the model of H-N [25,26] which can be considered a generalization of the latter models. These last three models, called anomalous models, are the most relevant in the literature; however, there are other models which approach the phenomena from a different perspective, as the models by Weron and collaborators [13,27-32].

During the last years, some of these authors and others not mentioned above have addressed the anomalous relaxation processes of types C-C, C-D and H-N with the help of fractional calculus, for example, by associating the response functions to the Mittag-Leffler functions (which turn out to be related with the complex susceptibility s = iω through the Laplace transform) [11]; by modeling relaxation processes through equations with fractional derivatives, whose solutions are also given in terms of Mittag-Leffler functions; and establishing other theoretical connections discovered through the use of fractional calculus tools.

In 1929, Debye conceived a simple model for the relaxation process in which he supposed a unique relaxation time for all molecules, obtaining the following expression [18]:

equation (1)

where equation is the complex susceptibility, f* is the complex permittivity of the dielectric, the real value f0 is the low frequency dielectric constant, the real value f is the high frequency dielectric constant and σ is the constant associated to the dipole’s characteristic relaxation time.

Cole and collaborators [21,22] formulated a more complete model based on experimental data on dielectrics. This new formulation contains a parameter α which can assume values between 0 and 1 and is given by

equation (2)

In their studies of the dielectric relaxation process, C-D [24] proposed a modified equation for the dielectric liquids glycerol and propylene-glycol. They also studied n-propanol and confirmed that, in its case, the process is described by Debye’s equation. The difference from the former two liquids lies in the broadening of the dispersion range under higher frequencies [21].

The works of C-D [23,24] allowed for an almost complete determination of the dielectric properties and a more accurate quantitative description. The empirical expression for the complex susceptibility equation given by:

equation (3)

where 0 < β ≤ 1.

In their works, H-N [25,26] studied the complex dielectric behavior of twenty-one polymers and noticed that they had approximately the same form. They then arrived at an empirical expression which generalizes the dispersion models of Debye, C-C and C-D. This is the expression representing the relaxation process:

equation (4)

It is important to observe that for β = 1 we have the C-C model, while α = 1 takes to C-D model and, finally, with α = β = 1 we recover the first model proposed by Debye.

Kinetic Equations

The properties of dielectric materials are usually described by two constants ε' and ε'' which are called dielectric constants (or loss factors). They can be combined in a complex dielectric constant given by

equation (5)

This constant, known as complex dielectric permittivity, is given by the following superposition relation [33]:

equation (6)

where equation is the Laplace transform of f(t) in variable iω. We have that f(t) is the normalized polarization decay function when a macroscopic electric field is removed from its medium. Function f(t) contains only the contributions from the relaxation process and we have chosen f(0) = 1 [34].

In the case of linear approximation response, the polarization changes caused by thermal motion are the same as for the macroscopic function dipole relaxation induced by the electric field [35]. Therefore, the laws governing the dipole correlation function φ(t) are directly related to the kinetic properties and macroscopic structures of the dielectric system, represented by function f(t). Thus, it is possible to equate the relaxation function f(t) to the macroscopic dipole correlation function φ(t) as follows:

equation (7)

where M(t) is the fluctuating macroscopic dipole moment.

In the projection operator formalism developed by Mori [36] and Zwanzig [37], function φ(t) is called temporal correlation function, as the dipole correlation function defined above is a specific case of the temporal correlation function. Thus, function φ(t), now called temporal correlation function, has the following form in the approximation context [38]:

equation (8)

This is an integro-differential equation which takes into account the effects of memory. Therefore, introducing the concept of an integral memory function given by equation and using the fact that the relaxation function f(t) also satisfies Equation (8), it is possible to obtain the following relation involving function f(t):

equation (9)

where * denotes a convolution product.

The relations given by Equations (9) and (6) can now be used to calculate the integral memory function M(t). Applying the Laplace transform to Equation (9) we obtain

equation (10)

where equation is the Laplace transform of M(t). Substituting Equation (10) into Equation (6), we obtain the following relation:

equation (11)

Comparing the relation expressed by Equation (11) with the classical empirical laws (1)-(4) (where s = iω), it is possible to obtain the corresponding memory functions:

equation (12)

equation (13)

equation (14)

equation (15)

Applying the inverse Laplace transform to functions (12)-(15) we find that the memory functions in the time variable are respectively given by:

equation (16)

equation (17)

equation (18)

equation (19)

with 0 < α ≤ 1, 0 < β ≤ 1 and where equation is the Mittag-Leffler function with three parameters. This function contains as particular cases the Mittag-Leffler with two parameters (c = 1) and one parameter (b = c = 1).

The Mittag-Leffler functions with two and three parameters are defined, respectively, by

equation (20)


equation (21)

where Re(a)> 0, Re(b)> 0, Re(c)> 0 and equation is the Pochhammer symbol. Substituting the memory functions given by Equations (16)-(19) into Equation (9), it is possible to obtain the kinetic equations associated with the models of Debye, C-C, C-D and H-N, respectively:

equation (22)

equation (23)

equation (24)

equation (25)

Expression equation denotes the Riemann-Liouville fractional derivative, defined by

equation (26)

The solutions of kinetic equations (22)-(25) are given respectively by:

equation (27)

equation (28)

equation (29)

equation (30)

The graphic representation of solutions (27)-(30) can be seen in Figures 1-4, respectively1.


Figure 1: Debye function.


Figure 2: Cole-Cole function.


Figure 3: Cole-Davidson function.


Figure 4: Havriliak-Negami function.

Fractional Kinetic Equations

As mentioned above, the generalization of the laws describing anomalous relaxation phenomena in dielectrics has been the object of studies by several important researchers. Here we turn our attention to a recent study [39] in which the following initial value problem has been considered:

equation (31)

with u0 and λ positive constants. Parameters α and β are subject to the conditions

0 < α ≤ 1 and − α < β ≤ 1 − . (32)

Fractional equation (31) appears as a generalization of the wellknown model by Kohlrausch (Kohlrausch-Willians-Watts) based on the “stretched” exponential function [40-42] and its solution is represented in terms of a function introduced by Kilbas and Saigo [43,44], which is given by the series

equation (33)

with α, m, l ∈equation, α > 0, m > 0 and α(im + l) ≠ −1, −2, −3, . . . Without loss of generality, we can take λ = 1 in problem (31); the solution is then

equation (34)

The conditions given in Equation (32) guarantee the existence and complete monotonicity of solution u(t) for t ≥ 0. We now mention some particular cases. For β = 0 we have the particular case equation where Eα(•) is the classical Mittag- Leffler function [16]. The second particular case is associated with the well- known “stretched” exponential function and is given by

equation where α = 1 and −1 < β ≤ 0.

Still supposing β = 0 one finds the solution u(t) = E1,1,0(−t) = exp(−t), which is the exponential solution to the relaxation equation. With an analogous generalization purpose, Garra et al. [14], using the general theory of the hyper-Bessel operator [3], solved a generalized relaxation equation which was called fractional differential equation with time-variant coefficient.

In that equation, operator equation replaces the derivatives in the relaxation equation.

Here, the aim of the paper, we will consider the following fractional differential equation:

equation (35)

In this equation, operator equation is the Riemann-Liouville fractional derivative, defined in Equation (26). Imposing the normalization condition equation it is possible to obtain the fractional kinetic equations with the help of the memory functions explained in the previous section. First, memory function Equation (16) is substituted into Equation (35) in order to obtain:

equation (36)

Applying the Laplace transform to it, we obtain

equation (37)

Then, isolating equation we get

equation (38)

Calculating its inverse Laplace transform, Debye’s fractional relaxation function turns out to be

equation (39)

For γ = 1 we recover Debye’s exponential solution equation In the same way, substituting Equation (17) into Equation (35) we have

equation (40)

Again, applying Laplace transform to Equation (40) we obtain the following expression:

hence emerges the C-C fractional relaxation function:

equation (41)

equation (42)

For α = 1 we recover Debye’s fractional function Equation (39). For γ = 1 one has the C-C model given by Equation (28).

In order to generalize the C-D model, Equation (35) can be substituted into Equation (18), in order to obtain

equation (43)

Applying the Laplace transform to the last expression, we obtain:

equation (44)

Expanding the exponential function and the Mittag-Leffler function with two parameters in power series, we get

equation (45)

Multiplying and dividing by Γ(β + j + βk) and rearranging factors, this gives:

equation (46)

which results in

equation (47)

or, isolating equation

equation (48)

Thus, using the inverse Laplace transform we finally get

equation (49)

Now, taking β = 1 in Equation (49) and using the following relation involving the Mittag-Leffler function with two parameters:

equation (50)

we recover the function of Debye’s fractional model, given by Equation (39). For γ = 1 we recover the C-D function, given by Equation (29).

Finally, we substitute H-N’s memory function, given by Equation (19), in Equation (35), we have:

equation (51)

Applying the Laplace transform and using the following equality given [11],

equation (52)

we obtain

equation (53)

or, equivalently,

equation (54)

From this expression it follows that

equation (55)

Applying the inverse Laplace transform, we obtain

equation (56)

If α = 1 in Equation (56), we recover the function of the fractional C-D model Equation (49). If we put β = 1 in Equation (56), it is then possible to recover the function of the fractional C-C model given by Equation (42). Finally, for α = β = 1 we have function of the fractional Debye model, Equation (39). The H-N model Equation (30) is recovered when γ = 1.

Therefore, we have just shown that the models by Debye, C-C, C-D and H-N given by Equations (27), (28), (29) and (30) can be considered particular cases of the fractional models given by Equations (39), (42), (49) and (56), respectively.

The complex dielectric permittivity (in variable s) is given in the fractional case by the following superposition relation:

equation (57)

With this definition, we recover in the four fractional relaxation function models the corresponding classical empirical responses described by Equations (1), (2), (3) and (4). To see this we need only to substitute Equations (38), (41), (48) and (55) into Equation (57). Figures 5-8 graphically represent the solutions given by Equations (39), (42), (49) and (56), respectively.


Figure 5: Fractional Debye function.


Figure 6: Fractional Cole-Cole functions with α = 0.5.


Figure 7: Fractional Cole-Davidson function with β = 0.5.


Figure 8: Fractional Havriliak-Negami functions with α = 0.5.


This study has shown that certain fractional kinetic equations with a parameter 0 < γ ≤ 1, whose solutions are given in terms of Mittag- Leffler functions, can be used as generalized models of the classical relaxation kinetic equations. The classical solutions turned out to correspond to the particular case γ = 1.

In all the cases discussed, the results of the corresponding classical models (response functions) were recovered, as well as the empirical dielectric complex functions (in variable s).

The graphs of the fractional relaxation functions and of the solutions of the fractional kinetic equations were constructed. The analysis of fractional relaxation function graphs and their algebraic expressions given by Equations (39), (42), (49) and (56) has shown that the function’s signals and their complete monotonicity depend on the parameters involved. These aspects deserve a more careful analysis and will be the subject of a future work [45].


We are indebted to Dr. J. Emlio Maiorino for several and usefull discussions. (ECFAR) thanks CAPES for a research grant.

1In all graphs we have σ = 1.


Select your language of interest to view the total content in your interested language
Post your comment

Share This Article

Article Usage

  • Total views: 12823
  • [From(publication date):
    November-2015 - Jan 23, 2018]
  • Breakdown by view type
  • HTML page views : 8848
  • PDF downloads : 3975

Post your comment

captcha   Reload  Can't read the image? click here to refresh

Peer Reviewed Journals
Make the best use of Scientific Research and information from our 700 + peer reviewed, Open Access Journals
International Conferences 2018-19
Meet Inspiring Speakers and Experts at our 3000+ Global Annual Meetings

Contact Us

Agri & Aquaculture Journals

Dr. Krish

[email protected]

1-702-714-7001Extn: 9040

Biochemistry Journals

Datta A

[email protected]

1-702-714-7001Extn: 9037

Business & Management Journals


[email protected]

1-702-714-7001Extn: 9042

Chemistry Journals

Gabriel Shaw

[email protected]

1-702-714-7001Extn: 9040

Clinical Journals

Datta A

[email protected]

1-702-714-7001Extn: 9037

Engineering Journals

James Franklin

[email protected]

1-702-714-7001Extn: 9042

Food & Nutrition Journals

Katie Wilson

[email protected]

1-702-714-7001Extn: 9042

General Science

Andrea Jason

[email protected]

1-702-714-7001Extn: 9043

Genetics & Molecular Biology Journals

Anna Melissa

[email protected]

1-702-714-7001Extn: 9006

Immunology & Microbiology Journals

David Gorantl

[email protected]

1-702-714-7001Extn: 9014

Materials Science Journals

Rachle Green

[email protected]

1-702-714-7001Extn: 9039

Nursing & Health Care Journals

Stephanie Skinner

[email protected]

1-702-714-7001Extn: 9039

Medical Journals

Nimmi Anna

[email protected]

1-702-714-7001Extn: 9038

Neuroscience & Psychology Journals

Nathan T

[email protected]

1-702-714-7001Extn: 9041

Pharmaceutical Sciences Journals

Ann Jose

[email protected]

1-702-714-7001Extn: 9007

Social & Political Science Journals

Steve Harry

[email protected]

1-702-714-7001Extn: 9042

© 2008- 2018 OMICS International - Open Access Publisher. Best viewed in Mozilla Firefox | Google Chrome | Above IE 7.0 version