# Fractional-calculus diffusion equation

- Abdul-Wali MS Ajlouni
^{1}Email author and - Hussam A Al-Rabai'ah
^{2}

**4**:3

**DOI: **10.1186/1753-4631-4-3

© Ajlouni and Al-Rabai'ah; licensee BioMed Central Ltd. 2010

**Received: **2 November 2009

**Accepted: **21 May 2010

**Published: **21 May 2010

## Abstract

### Background

Sequel to the work on the quantization of nonconservative systems using fractional calculus and quantization of a system with Brownian motion, which aims to consider the dissipation effects in quantum-mechanical description of microscale systems.

### Results

The canonical quantization of a system represented classically by one-dimensional Fick's law, and the diffusion equation is carried out according to the Dirac method. A suitable Lagrangian, and Hamiltonian, describing the diffusive system, are constructed and the Hamiltonian is transformed to Schrodinger's equation which is solved. An application regarding implementation of the developed mathematical method to the analysis of diffusion, osmosis, which is a biological application of the diffusion process, is carried out. Schrödinger's equation is solved.

### Conclusions

The plot of the probability function represents clearly the dissipative and drift forces and hence the osmosis, which agrees totally with the macro-scale view, or the classical-version osmosis.

## Introduction

In this paper we aim to consider the dissipation effects, appeared in the will-known diffusion process, quantum-mechanically depending on the procedure of the quantization of nonconservative systems using fractional calculus [1–4], which was also applied on the related phenomenon, the Brownian motion [5].

Most of the natural laws of physics, such as Maxwell's equations, Newton's laws of motion, and Schrödinger equation, are stated, or can be, in terms of partial deffirential equations (PDEs), that is, these laws describe physical phenomena by relating space and time derivatives. Diffusion equation, or heat flow, is one of the most important PDEs in physical sciences. The basic process in the diffusion phenomenon is the flow of the fluid from a region of higher density to one of lower density [6].

The tendency of a statistical ensemble to achieve thermodynamic equilibrium with a uniform distribution of states for its constituent subsystems does not have to be monotonic in time. In general, equilibration takes place in stages and is characterized by several stochastization times with vastly different orders of magnitude. Thermodynamic equilibrium has no absolute meaning and depends on the time scale over which a given process is analyzed.

In a diffusion process or chemical reaction, Fick's law provides a linear relationship between the flux of molecules and the chemical potential difference. Likewise, a direct proportionality exists between the heat flux and the temperature difference in a thermally conducting slab, as expressed by Fourier's law. Diffusion of gases between air in the lungs and blood proceeds in the direction from high to low concentration, and the rate of diffusion is greatest when the difference in concentration is greatest. Diffusion obeys Fick's law, but the actual rate of exchange is greatly affected by hemoglobin in the blood.

## Diffusion equation

*diffusion coefficient*,

*D*, of the material, given in general by the

*Stokes-Einstein equation*:

where *k*_{
B
}is the Boltzmann constant, *T* is the absolute temperature in K, and *f* is a frictional coefficient. The diffusion coefficient, D, is defined as the net flow of particles per unit time across an imaginary plane of unit area lying at right angles to the concentration gradient, that gradient also having unit strength.

*f*is the frictional coefficient, (

*dv/dt*) is the acceleration and m is the mass of the molecule. In the case of spherical particles, the translational frictional force

*f*is proportional to the fluid viscosity,η, and radius r of the particle. Thus the coefficient of friction for spherical particles, known as Stokes law, is

*f*comes into effect when a molecule moves through a medium. The movement of the molecule could be either diffusion or sedimentation and the driving force, F, can be the concentration gradient, the force of gravity or the centrifugal force. According to Fick's law, the rate of diffusion across a boundary (

*dn/dt*: the number of molecules which pass through a cross section

*A*in unit time) for a single solute component diffusing in a system at constant temperature and pressure is given by [7]:

where
is the concentration gradient. A concentration gradient implies that the concentration of the molecules (i.e. the solute) varies with distance *r* in one dimension.

*A*

_{1},

*A*

_{2},

*v*

_{1}, and

*v*

_{2}, from the equation of continuity, we have:

*C*

_{1}to a lower value of

*C*

_{2}over a short length (

*d*), then the mass (

*m*) of the solute diffusing down the pipe in time (

*t*) is

where C_{0} is the total solute concentration difference across the boundary.

## Quantization of diffusion process

*V'*(

*x*) by -

*δ*(

*x'*-

*x*) and the entire potential,

*V*(

*x*), will be written as

*F*(

*x'*) is obtained directly

Thus, we obtain a definition of the random force, that agrees with our assumption and with the fact that ⟨*F(x)*⟩ = 0 of Eq.(2).

which is the classical equation of motion of the diffused particle, Eq.(2).

Here, *p*_{0} and *p*_{1/2} are the canonical conjugate momenta to *q*_{0} and *q*_{1/2} respectively.

where *q*_{0} = *x* and *q*_{1/2} = *y*.

*u*. This leads to

*H*_{
n
}being Hermite polynomials.

*y*=

*y'*of Eq.(39) will be ignored since the impulse potential effects will be considered in the part of Eq. (37) which will be written as

where *H*(*u* - *u'*) Heaviside step function.

*q*

_{ s(i) }, Ψ is expressed as

## Application: osmosis

Osmosis is a physical phenomenon that has been extensively studied by scientists in various disciplines of science and engineering. Early researchers studied the mechanism of osmosis through natural materials, and from the 1960s, special attention has been given to osmosis through synthetic materials. Following the progress in membrane science in the last few decades, especially for reverse osmosis applications, the interests in engineered applications of osmosis has been spurred. Osmosis, or as it is currently referred to as forward osmosis, has new applications in separation processes for wastewater treatment, food processing, and seawater/brackish water desalination. Other unique areas of forward osmosis research include pressure-retarded osmosis for generation of electricity from saline and fresh water and implantable osmotic pumps for controlled drug release. This paper provides the state-of-the-art of the physical principles and applications of forward osmosis as well as their strengths and limitations.

*Osmosis*is usually defined as the transport of molecules in a fluid through a semipermeable membrane due to an imbalance in its concentration on either side of the membrane [6]. Osmosis may be by diffusion, but it may also be a bulk flow through pores in a membrane. If a plant cell is put into a concentrated solution of sugar, for example, Fig. 1, the pressure on the right is then, in either case, water moves from a region of high concentration to a region of low concentration greater than the pressure on the left by an amount

*hρg*, where

*ρ*is the density of the liquid on the right and is called the relative osmotic pressure. The general formula for the osmotic pressure

*P*of a solution containing

*n*moles of solute per volume

*V*of solvent is [6]

The net osmotic pressure exerted on a semipermeable membrane separating the two compartments is thus the difference between the osmotic pressures of both compartments.

Thus, we obtain a definition of the new random force, but with a complete disagreement with our assumption with the fact that ⟨*F(x)*⟩ = 0 of Eq.(2). This disagreement appears from the fact that pressure exerted a net drift force *F*_{
p
}on the particles in the direction of osmosis (Fig. 1).

## Discussion

Fraction calculus is very helpful expressing the dissipation, as well as in quantizing nonconservative systems associated with many important physical problems: either where the ordinary quantum-mechanical treatment leads to an incomplete description, such as the energy loss by charged particles when passing through matter; or where it leads to complex nonlinear equations, such as Brownian motion, and diffusion.

By using fractional calculus in physical problems it is possible to create a whole mechanical description of nonconservative systems, including Lagrangian and Hamiltonian mechanics, canonical transformations, Hamilton-Jacobi theory, and quantum mechanics. In this paper, an important physical nonconservative system, which is the diffusion process is treated quantum-mechanically, for the first time using fractional-calculus.

^{2}of Eq. (51), connected to of the "osmosisized" particle including the drift, and frictional forces; the osmosis process is manifested very clearly, where the confinement of these particles to one region of space gradually leads to a situation in which the particles uniformly fills all the available space in the high concentration, where the new Heaviside step function,

*H*

_{ p }, is modified due to drift force to show non-step behavior; which agrees totally with the macro-scale view, or the classical-version osmosis.

## Conclusion

Diffusion and the diffusion equation are central topics in both Physics and Mathematics, and their ranges of applicability span from astrophysical dynamics to the diffusion of particles governed by Schrödinger's equation.

The quantization of a system with diffusion process has been carried out according to the theory proposed recently [1, 2]. A potential, and a Hamiltonian, corresponding to the random force, and dissipative force, were constructed. The relevant Schrödinger's equation has then been set and solved. The classical equation of motion, connected to diffused particle, could be obtained easily from the fractional Lagrangian. The random and frictional forces were plotted; the diffusion process manifested very clearly. The next step could be to study problems such as the correlation functions, transport equation, chemical potential, entropy, etc., on a quantum-mechanical basis.

An application of the developed mathematical method to the analysis of diffusion in a biological medium, osmosis, is carried out. Schrödinger's equation is solved. The plot of the probability function represents clearly the dissipative and drift forces and hence the osmosis, which shows the same macro-scale view of the osmosis.

## Appendix: Fractional calculus

where *J*^{α} represents the fractional integral operator of order *α*, and R^{+} represents the set of positive real numbers.

*m*such that

*m*-1 <

*α*≤

*m*the fractional derivative of order

*α*> 0 may be defined as

*D*

^{ a }being the fractional deferential operator of order

*a*Equation (2) may be rewritten using Eq. (1) as follows:

Here, we formulate the problem in terms of the left fractional derivative the left Riemann-Liouville fractional derivatives, which are defined in Eqs. (A.1, A.2). Most of the left fractional operations also hold for the right ones. For the left operations *f*(*t*) must vanish for *t* < *a* while *f*(*t*) = 0 for *t* > *b* for the right operation. Thus, the left operations are causal. Conversely, the right operations are anti-causal [16]. From the physical point of view, when we differentiate with respect to time, the right differentiation represents an operation performed on the future state of the process *f*(*t*) [17].

Fractional integral and differential operators have the following properties [14, 15]:

*I*, the identity operator:

but the inverse application of the two operators is not necessarily true.

*n*> 0,

*J*

^{ n }and

*D*

^{ n }are linear operators, i.e.,

*c*,

*J*

^{ n }and

*D*

^{ n }are homogeneous operators, i.e.,

*α*,

*β*> 0,

*J*

^{ n }obeys the additive index law, but not necessarily

*D*

^{ n }, i.e.,

*t*-

*a*)

^{ β }, which are given by

For *α* = 1/2 this equation is called semi-derivative; for *α* = - 1/2 it is called semi-integral.

## Authors' information

**Abdul-Wali Ajlouni (Jordan)** is the chair of the Applied Physics Department at Tafila Technical University, Jordan. Prior to this post, he served as an Associate Dean of the Faculty of Science, and as an assistant Professor. Before joining TTU, Dr. Ajlouni was a researcher on radiation safety for the Ministry of Energy. He has received training on radiation protection and radiological emergency preparedness in Jordan, Turkey, the United States, Czechoslovakia, Egypt, Iran, and Syria. In addition, he has attended numerous conferences on nuclear energy, and Applied mathematics around the world. Dr. Ajlouni has published extensively on the two subjects, with titles such as: "Mathematical Model for Dispersion of Nuclear Pollutants in Jordan Atmosphere," and "Quantization of Brownian Motion." Dr. Ajlouni received a SCOPUS Scientific Research and the International Conference on Mathematics and Information Security (ICAMIS2009) awards in 2009. He holds a B.Sc. in Physics from Yarmouk University, Jordan, an M.S. in Nuclear Physics from Mustansyria University, Iraq, and a Ph.D. in Mathematical/Radiation Physics from the University of Jordan.

## Declarations

## Authors’ Affiliations

## References

- Ajlouni A-W: Ph.D. thesis: Quantization of Nonconservative Systems. 2004, Jordan University, Amman, JordanGoogle Scholar
- Ajlouni A-W, Invited Talk: Fractional Calculus Tools in Physics. Applied Mathematics & Information Sciences. 2010, Presented in the International Conference on Mathematics and Information Security (ICAMIS2009) Nov. 13th – 15th , 2009, Sohag University, Egypt.Google Scholar
- Rabei E, Ajlouni A-W, Ghassib H: Quantization with Fractional Calculus. 9th WSEAS International Conference On Applied Mathematics (MATH '06), MATH, TELE-INFO and SIP '06, Istanbul, Turkey, May 27-29. 2006Google Scholar
- Rabei E, Ajlouni A-W, Ghassib H: Quantization of Nonconservative Systems Using Fractional Calculus. Wseas Transactions on Mathematics. 2006, 5: 853-864.MathSciNetGoogle Scholar
- Rabei E, Ajlouni A-W, Ghassib H: Quantization of Brownian Motion. Int J of Theoretical Physics. 2006, 45: 1619-1629.View ArticleADSMathSciNetMATHGoogle Scholar
- Tuszynski J, Kurzynski M: Introduction to Molecular Biophysics. 2003, CRC Press, FloridaView ArticleGoogle Scholar
- Vasantha Pattabhi V, Gautham N: Biophysics. 2002, Kluwer Academic Publishers,New York, USAGoogle Scholar
- Rabei EM, Al-halholy T, Rousan A: Potentials of Arbitrary Forces with Fractional Derivatives. International Journal of Modern Physics A. 2004, 19: 3083-3089. 10.1142/S0217751X04019408.View ArticleADSMathSciNetMATHGoogle Scholar
- Duffy D: Green's Functions with Applications. 2001, Chapman & Hall/CRC, New York, 1View ArticleMATHGoogle Scholar
- Riewe F: Nonconservative Lagrangian and Hamiltonian Mechanics. Physical Review E. 1996, 53: 1890-1898. 10.1103/PhysRevE.53.1890.View ArticleADSMathSciNetGoogle Scholar
- Riewe F: Mechanics with Fractional Derivatives. Physical Review E. 1997, 55: 3581-3592. 10.1103/PhysRevE.55.3581.View ArticleADSMathSciNetGoogle Scholar
- Griffiths DJ: Introduction to quantum Mechanics. 1995, Prentice Hall. New JerseyMATHGoogle Scholar
- Merzbacher E: Quantum Mechanics. 1970, Wiley, 2Google Scholar
- Oldham B, Spanier J: The Fractional Calculus. 1974, Academic Press, NewYorkMATHGoogle Scholar
- Carpintri A, Mainardi F: Fractals and Fractional Calculus in Continuum Mechanics. 1997, Springer, New YorkView ArticleGoogle Scholar
- Dreisigmeyer DW, Young PM: Nonconservative Lagrangian mechanics: a generalized function approach. J Phys A: Math Gen. 2003, 36: 8297-8310. 10.1088/0305-4470/36/30/307.View ArticleADSMathSciNetMATHGoogle Scholar
- Agrawal OP: Formulation of Euler-Lagrange equations for fractional variational problems. Journal of Mathematical Analysis and Applications. 2002, 272: 368-379. 10.1016/S0022-247X(02)00180-4.View ArticleMathSciNetMATHGoogle Scholar

## Copyright

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.