Medical, Pharma, Engineering, Science, Technology and Business

**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.

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]:

(1)

where is the complex susceptibility, *f*_{*} is the complex permittivity of the dielectric, the real value *f*_{0} 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

(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 given by:

(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:

(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.

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

(5)

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

(6)

where 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:

(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]:

(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 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):

(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

(10)

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

(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:

(12)

(13)

(14)

(15)

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

(16)

(17)

(18)

(19)

with 0 < α ≤ 1, 0 < β ≤ 1 and where 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

(20)

and

(21)

where Re(a)> 0, Re(b)> 0, Re(c)> 0 and 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:

(22)

(23)

(24)

(25)

Expression denotes the Riemann-Liouville fractional derivative, defined by

(26)

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

(27)

(28)

(29)

(30)

The graphic representation of solutions (27)-(30) can be seen in **Figures 1**-**4**, respectively^{1}.

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:

(31)

with u_{0} 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

(33)

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

(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 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

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

Still supposing β = 0 one finds the solution u(t) = E_{1,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 replaces the derivatives in the relaxation equation.

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

(35)

In this equation, operator is the Riemann-Liouville fractional derivative, defined in Equation (26). Imposing the normalization condition 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:

(36)

Applying the Laplace transform to it, we obtain

(37)

Then, isolating we get

(38)

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

(39)

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

(40)

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

hence emerges the C-C fractional relaxation function:

(41)

(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

(43)

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

(44)

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

(45)

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

(46)

which results in

(47)

or, isolating

(48)

Thus, using the inverse Laplace transform we finally get

(49)

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

(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:

(51)

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

(52)

we obtain

(53)

or, equivalently,

(54)

From this expression it follows that

(55)

Applying the inverse Laplace transform, we obtain

(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:

(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.

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.

^{1}In all graphs we have σ = 1.

- Tenreiro MJ,KiryakovaV,Mainardi F (2010) A poster about the old history of fractional calculus.FractCalcAppl Anal13: 447-454.
- Hilfer R (2000) Applications of Fractional Calculus in Physics. World Scientific Publishing Co. Pte. Ltd, Singapore.
- Kiryakova V (1994) Generalized Fractional Calculus and Application. Longman Scientific and Technical, Harlow, New York.
- Mainardi F (2010) Fractional Calculus and Waves in Linear Viscoelasticity. Imperial College Press, London.
- Haubold HJ, MathaiAM,Saxena RK(2011) Mittag-Leffler functions and their applications. J ApplMathl: 1-51.
- Oustaloup A (1991) CRONE: CommandeRobusted’Ordre NonEntier. Hermes, Paris.
- Pan I, Das S (2013) Intelligent Fractional Order Systems and Control: An Introduction, Springer Publishing CompanyINC, Beijing.
- Podlubny I (1999) Fractional Differential Equations. Academic Press, San Diego.
- Uchaikin VV (2013) Fractional Derivatives for Physicists an engineers. Nonlinear Physical Science, Beijing.
- UchaikinVV, Sibatov R (2013) Fractional Kinetics in Solids: Anomalous Charge Transport in Semiconductors. Dielectrics and Nanosystems, World Scientific Publishing Co. Pte Ltd, Russia.
- Capelas de Oliveira E, MainardiF, VazJr J (2011) Models based on Mittag-Leffler functions for anomalous relaxation in dielectrics. The European Physical Journal 193: 261-171.
- Hanyga,Seredynska M (2008)On mathematical framework for the constitutive equations of anisotropic dielectric relaxation. J Stat Phys131: 269-303.
- Hilfer H (2002) Analytical representations for relaxation functions of glasses. J Non-Cryst Solids 305: 122-126.
- Garra R, Giusti A, MainardiR, PagniniG(2014) Fractional relaxation with time-varying coefficient FractCalcAppl Anal 17: 424-439.
- Khamzin RR,NigmatullinR, Popov II (2014) Justification of the empirical laws of the anomalous dielectric relaxation in the framework of the memory function formalism. FractCalcAppl Anal 17: 247-258.
- Gorenflo R, Kilbas AA, MainardiF, Rogosin SV (2014)Mittag-Leffler Functions, Related Topics and Applications, Springer in Monographs in Mathematics. Springer-Verlag, Berlin- Heidelberg.
- MainardiF, Garrappa R (2015)On complete monotonicity of the Prabhakar function and non-Debye relaxation in dielectrics. J ComputPhys293: 70-80.
- Gross B, Course A (1981) On the Theory of Dielectric Relaxation and Topics on the Linear Re- sponse Theory, IFQSC-USP, Sa˜o Carlos.
- BottcherCJF,Bordewijk P (1978) Theory of Electric Polarization, Elsevier Sci. Publ, Amsterdam.
- Debye PJW (1929) Polar Molecules, Dover, New York.
- Cole S, Cole RH (1941) Dispersion and absorption in dielectrics I Alternating current characteristics. J ChemPhys 9: 341-351.
- Cole S, Cole RH (1942) Dispersion and absorption in dielectrics II Direct current characteristics. J ChemPhys10: 98-105.
- Davidson DW, Cole RH (1950) Dielectric relaxation in glycerine. J ChemPhys18: L1417.
- Davidson DW, Cole RH (1951) Dielectric relaxation in glycerol, propylene glycol, and n-propanol. J ChemPhys19: 1484-1490.
- Havriliak Jr. S, Negami S (1966) A complex plane analysis of a-dispersions in some polymer systems. J Polymer Sci14: 99-117.
- Havriliak Jr. S, Negami S (1967)A complex plane representation of dielectric and mechanical relaxation processes in some polymers. JPolymer Sci 8: 161-210.
- Jurlewicz, Weron K (2002) Relaxation of dynamically correlated clusters. J Non-Cryst Solids 305: 112-121.
- Jurlewicz A, Weron K, Teuerle M (2008) Generalized Mittag-Leffler relaxation: clustering-jump continuous-time random walk approach. See comment in PubMed Commons below Phys Rev E Stat Nonlin Soft Matter Phys 78: 011103.
- Stanislavsky K, WeronJ (2010) Trzmiel, Subordination model of anomalous diffusion leading to the two-power-law relaxation responses.EurophysLett 91: 40003/1.
- Szabat K, Weron, Hetman P (2007) Heavy-tail properties of relaxation time distributions underlying the Havriliak-Negami and the Kohlrausch-Williams-Watts relaxation patterns. JNon-Cryst Solids 353: 4601-4607.
- Hilfer R (2002) H-function representations for stretched exponential relaxation and non-Debye susceptibilities in glassy systems. Phys Rev E Stat Nonlin Soft Matter Phys 65: 061510.
- Hanyga, Seredynska M (2008)On a mathematical framework for the constitutive equations of anisotropic dielectric relaxation. J Stat Phys 131: 269-303.
- Frohlich H (1958) Theory of Dielectrics. Oxford University Press, London.
- Manning F, Bell ME (1940) Electrical conduction and related phenomena in solid di- electrics. Rev Mod Phys12: 215-257.
- Williams G (1972) Use of the dipole correlation function in dielectric relaxation. J Chem Rev72: 55-69.
- Mori H (1965) A continued-fraction representation of the time-correlation functions.ProgTheorPhys34: 399-416.
- Zwanzig R (1961) Lectures in Theoretical Physics.Interscience, New York.
- Boon P, Yip S (1980) Molecular Hydrodynamics. Dover, New York.
- Capelas de Oliveira E, MainardiF,Vaz Jr. J(2014) Fractional models of anomalous relaxation based on the Kilbas and Saigo function. New Trends in Fluid and Solid Mechanical Models 49: 2049-2060.
- Williams G, Watts DC (1970) Non-symmetrical dielectric relaxation behaviourarisisng from a simple empirical decay function. Trans Faraday Soc66: 80-85.
- Anderssen RS, Husain SA, Loy RJ (2004) The Kohlrausch function: properties and applications.Anziam J 45E: C800-C816.
- Kohlrausch R (1854) Theorie des elektrischenruckstandes in der leidenerflasche.Pogg Ann PhysChem91: 179-214.
- KilbasA, Saigo M(1995) Fractional integral and derivatives of Mittag-Leffler type function. DoklAkadNaukBelarusi 39: 22-26.
- KilbasA, Saigo M (1995) Solution of Abel integral equations of the second kind and of differential equation of fractional order.DoklAkadNaukBelarusi 39: 29-34.
- Ester CFAR (2015)On the complete monotonicity of fractional relaxation functions. To be submitted.

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

- Algebraic Geometry
- Analytical Geometry
- Axioms
- Behaviometrics
- Big Data Analytics
- Binary and Non-normal Continuous Data
- Binomial Regression
- Biometrics
- Biostatistics methods
- Clinical Trail
- Complex Analysis
- Cross-Covariance and Cross-Correlation
- Differential Equations
- Fourier Analysis
- Genetic Linkage
- Hamilton Mechanics
- Hypothesis Testing
- Integration
- Large-scale Survey Data
- Matrix
- Microarray Studies
- Multivariate-Normal Model
- Noether's theorem
- Non rigid Image Registration
- Physical Mathematics
- Quantum Mechanics
- Quantum electrodynamics
- Regressions
- Relativity
- Riemannian Geometry
- Robust Method
- Soft biometrics
- Spatial Gaussian Markov Random Fields
- Statistical Methods
- Theoretical Physics
- Theory of Mathematical Modeling
- Topology
- mirror symmetry
- vector bundle

- Total views:
**12823** - [From(publication date):

November-2015 - Jan 23, 2018] - Breakdown by view type
- HTML page views :
**8848** - PDF downloads :
**3975**

Peer Reviewed Journals

International Conferences
2018-19