# Spatial Rotation of the Fractional Derivative in Two-Dimensional Space

Authors: Ehab Malkawi
Journal: Advances in Mathematical Physics
 (2015)

---

## Abstract

The transformations of the partial fractional derivatives under spatial rotation inR
2 are derived for the Riemann-Liouville and Caputo definitions. These transformation properties link the observation of physical quantities, expressed through fractional derivatives, with respect to different coordinate systems (observers). It is the hope that such understanding could shed light on the physical interpretation of fractional derivatives. Also it is necessary to be able to construct interaction terms that are invariant with respect to equivalent observers.


---

## Body

## 1. Introduction

Fractional calculus deals with differentiation and integration to arbitrary real or complex orders. The idea is as old as the integer-order calculus. Extensive mathematical discussion of fractional calculus can be found in [1–5] and references therein. The techniques of fractional calculus have been applied to wide range of fields, such as physics, engineering, chemistry, biology, economics, control theory, signal image processing, and groundwater problems.Physical applications of fractional calculus span a wide range of topics and problems (for a review see [6–11] and references therein). Generalizing fractional calculus to several variables, multidimensional space, and generalization of fractional vector calculus has been reported [12–19]. Also, progress has been reported on generalization of Lagrangian and Hamiltonian systems [20–25].Despite the significant progress in applying fractional calculus to a wide range of physical problems, there is still a lack of satisfactory geometric and physical interpretation of fractional calculus, in comparison with the simple interpretations of their integer-order counterparts (see [26] and references therein). Effort has been devoted to relate fractional calculus and fractal geometry [27–31]. A different approach to geometric interpretation of fractional calculus is based on the idea of the contact of αth-order [32]. However, a satisfactory interpretation is still missing [26].A study of the symmetry of physical systems described by fractional calculus and investigating the transformation properties of fractional derivative operators under specific groups can shed light on the nature of such operators. For example, studying the transformation of the fractional operators under point transformation (or change of variables) can be found in [33–36], even though the class of variables substitution preserving the form of the fractional derivative is found to be very narrow. Generalizing the definition of fractional derivative of a function with respect to another function can be used to study general transformation of point transformation [37].The question that we are trying to address in this work is basically the following: given a quantity that is described by a fractional derivative of some function in a specific coordinate system(
x
,
y
), how can this quantity be expressed in a rotated coordinate system (
x
′
,
y
′
)? Thus, we try to relate the same fractional quantity in two equivalent coordinate systems. In this work, we only consider the Cartesian space coordinates in two dimensions, (
x
,
y
), and adhere to the notion of time as invariant in all coordinate systems (nonrelativistic). In other words, we are looking into the transformation of fractional derivative operators and their effect on physical quantities under the rotation group S
O
(
2
). We are not interested in constructing the irreducible representations of S
O
(
2
), as done in the standard integer-order differentiation case [38–41]; rather we focus on the orthonormal basis of the two-dimensional Cartesian space of R
2. Intuitively, quantities expressed through fractional derivative operators are expected to behave differently from scalar, vector, and tensor quantities.In this work we investigate the transformation properties of the fractional derivative and its action on an invariant scalar field, under space rotation in two dimensions, and using both the Riemann-Liouville and Caputo definitions. We compare between the Riemann-Liouville and Caputo definitions in relating the same physical quantity under rotation. We check if there is a major difference in the transformation properties between the two definitions. In [42], the transformation properties of the Riemann-Liouville fractional derivative of a scalar field under infinitesimal transformation of S
O
(
2
) are derived. The current work is more general and can be applied to any function and is not bounded to infinitesimal transformation. Also, we include the case of the Caputo definition of fractional derivative. Their special result agrees with our general results. In Section 2, we give a brief introduction to fractional calculus to lay out our notation. In Section 3, we lay out the S
O
(
2
) transformation properties of fields and their derivatives expressed in Cartesian coordinates. In Section 4, we apply our method to the Riemann-Liouville and Caputo fractional derivatives and lay out their corresponding transformations under S
O
(
2
). Finally, in Section 5 we give a brief discussion of our results.


## 2. Fractional Calculus

For mathematical properties of fractional derivatives and integrals one can consult [1–5] and the references therein. In this section we lay out the notation used in the rest of this work. We consider real analytic functions on R
2 and only discuss the Riemann-Liouville and Caputo definitions of the fractional derivative. Let f
(
x
,
y
) to be a real analytic function in a specific domain in the Euclidian space R
2; f
:
R
2
→
R. The x-partial fractional derivative of order α (keeping y constant), where a is the lower limit of x, is written as D
x
α
a
f
(
x
,
y
). Similarly the y-partial fractional derivative of order α (keeping x constant), where b is the lower limit of y, is written as D
y
α
b
f
(
x
,
y
). Since f
(
x
,
y
) is analytic, then the partial fractional derivatives are assumed to commute; that is, D
x
α
a
,
D
y
α
b
=
0.Definition 1.
Cauchy’s repeated integration formula of thenth-order integration of the function f
(
x
,
y
) along x, keeping y constant, can be written as(1)
D
x
-
n
a
f
x
,
y
=
∫
a
x
d
x
n
-
1
∫
a
x
n
-
1
d
x
n
-
2
⋯
∫
a
x
1
f
x
0
,
y
d
x
0
=
1
n
-
1
!
∫
a
x
f
u
,
y
d
u
x
-
u
1
-
n
.A similar relation for the nth-order integration of the function f
(
x
,
y
) along y, keeping x constant,(2)
D
y
-
n
b
f
x
,
y
=
v
∫
b
y
d
x
n
-
1
∫
b
y
n
-
1
d
y
n
-
2
⋯
∫
b
y
1
f
x
,
y
0
d
y
0
=
1
n
-
1
!
∫
b
y
f
x
,
u
d
u
y
-
u
1
-
n
.Definition 2.
The Riemann-Liouville fractional integration of orderα
<
0 and along x, keeping y constant, is defined as (3)
D
x
α
a
f
x
,
y
=
1
Γ
-
α
∫
a
x
f
u
,
y
d
u
x
-
u
1
+
α
,and along y, keeping x constant, is (4)
D
y
α
b
f
x
,
y
=
1
Γ
-
α
∫
b
y
f
x
,
u
d
u
y
-
u
1
+
α
,where Γ
(
·
) is the Gamma function.We considera
=
b
=
0 and drop them from the notation henceforth.Definition 3.
The Riemann-Liouville partial fractional derivatives of the orderα
>
0, where n
-
1
<
α
<
n and n
∈
N, are defined as (5)
D
x
α
f
x
,
y
=
D
x
n
D
x
α
-
n
f
x
,
y
=
1
Γ
n
-
α
∫
0
x
f
u
,
y
d
u
x
-
u
α
-
n
+
1
,
D
y
α
f
x
,
y
=
D
y
n
D
y
α
-
n
f
x
,
y
=
1
Γ
n
-
α
∫
0
y
f
x
,
u
d
u
y
-
u
α
-
n
+
1
,where to unify notation we used D
x
n
≡
d
n
/
d
x
n and D
y
n
≡
d
n
/
d
x
n.The Leibniz rule applied to the Riemann-Liouville partial fractional derivatives can be written as [1–5] (6)
D
x
α
f
x
,
y
g
x
,
y
=
∑
k
=
0
∞
Γ
α
+
1
Γ
k
+
1
Γ
α
-
k
+
1
D
x
α
-
k
f
x
,
y
D
x
k
g
x
,
y
,
(7)
D
y
α
f
x
,
y
g
x
,
y
=
∑
k
=
0
∞
Γ
α
+
1
Γ
k
+
1
Γ
α
-
k
+
1
D
y
α
-
k
f
x
,
y
D
y
k
g
x
,
y
.Definition 4.
The Caputo partial fractional derivatives of orderα
>
0, where n
-
1
<
α
<
n and n
∈
N, are defined as (8)
D
x
α
C
f
x
,
y
≡
1
Γ
n
-
α
∫
0
x
D
u
n
f
u
,
y
d
u
x
-
u
α
-
n
+
1
,
D
y
α
C
f
x
,
y
≡
1
Γ
n
-
α
∫
0
y
D
u
n
f
x
,
u
d
u
y
-
u
α
-
n
+
1
.The Riemann-Liouville and Caputo definitions of the fractional derivative are related [1–5] (9)
D
x
α
C
f
x
,
y
=
D
x
α
f
x
,
y
-
∑
k
=
0
n
-
1
x
k
-
α
Γ
k
+
1
-
α
D
x
k
f
x
,
y
x
=
0
,
D
y
α
C
f
x
,
y
=
D
y
α
f
x
,
y
-
∑
k
=
0
n
-
1
y
k
-
α
Γ
k
+
1
-
α
D
y
k
f
x
,
y
y
=
0
.As an example, for the functionf
(
x
,
y
)
=
x
β
y
λ, where β
,
λ
>
-
1, we find (10)
D
x
α
x
β
y
λ
=
Γ
β
+
1
Γ
β
-
α
+
1
x
β
-
α
y
λ
,
D
y
α
x
β
y
λ
=
Γ
λ
+
1
Γ
λ
-
α
+
1
x
β
y
λ
-
α
,
D
x
α
D
y
α
x
β
y
λ
=
Γ
β
+
1
Γ
λ
+
1
Γ
β
-
α
+
1
Γ
λ
-
α
+
1
x
β
-
α
y
λ
-
α
.The Caputo partial fractional derivatives give the same result for x
β
y
λ, except for β
,
λ
∈
N
<
n, where in this case the Caputo partial fractional derivatives vanish, as expected from (8).


## 3. Transformation of Functions and Integer-Order Derivatives under Spatial Rotation

We consider the Euclidian spaceR
2 and work with orthonormal basis of the Cartesian coordinate system. The distinction between covariant and contravariant quantities disappears and we drop this notation in our discussion. As discussed before, we are not interested in discussing the irreducible representations of S
O
(
2
) [38–41]; instead we keep working with the Cartesian coordinates (
x
,
y
) since we want to express the partial fractional derivatives in the same way. Consider the two Cartesian coordinate systems x
y and x
′
y
′ sharing the same origin (
0,0
). The (
x
′
,
y
′
) coordinates are related to the (
x
,
y
) coordinates though a counterclockwise rotation by an angle ϕ about their identical axis of rotation z
′
=
z. For ϕ
=
0 the two systems of coordinates are overlapping. In this work we consider the case of passive transformation of the coordinate system; thus we consider the observation of the same quantity from two different coordinate frames, without actually rotating physical quantities. The two coordinate systems are related by the transformation matrix (11)
R
ϕ
=
cos
⁡
ϕ
sin
⁡
ϕ
-
sin
⁡
ϕ
cos
⁡
ϕ
.For the active transformation case the angle ϕ would be replaced by (
-
ϕ
). The coordinates (
x
′
,
y
′
) are related to the coordinates (
x
,
y
) through (12)
x
′
y
′
=
cos
⁡
ϕ
sin
⁡
ϕ
-
sin
⁡
ϕ
cos
⁡
ϕ
x
y
.We present some of the properties of the functional transformations of the Lie group as applied to functionf
(
x
,
y
). The generator of the S
O
(
2
) rotation (angular momentum), around the z-axis, L is written as (13)
L
=
y
D
x
1
-
x
D
y
1
.The group elements of the rotation group are written, in the exponential form, as e
ϕ
L, where the exponential is defined as an infinite sum of differential operators, explicitly (14)
e
ϕ
L
≡
∑
n
=
0
∞
ϕ
n
n
!
L
n
,where ϕ is the angle of rotation. Since L
x
=
y and L
y
=
-
x, it is easy to derive the action of the element e
ϕ
L on the coordinates x and y, explicitly (15)
e
ϕ
L
x
=
cos
⁡
ϕ
x
+
sin
⁡
ϕ
y
=
x
′
,
e
ϕ
L
y
=
-
sin
⁡
ϕ
x
+
cos
⁡
ϕ
y
=
y
′
.Thus we retrieve the standard result in (12).Remark 5.
In general, forn
,
m
∈
Z
+,(16)
x
′
n
y
′
m
=
e
ϕ
L
x
n
e
ϕ
L
y
m
=
e
ϕ
L
x
n
y
m
.This result emerges from proving the following two lemmas.Lemma 6.
Given any two functionsA
(
x
,
y
) and B
(
x
,
y
) that are analytic (of Class C
∞) in their common domain, then(17)
L
n
A
B
=
∑
k
=
0
n
n
!
k
!
n
-
k
!
L
k
A
L
n
-
k
B
,where n
∈
Z
+. This is just the general Leibniz rule applied to the differential operator L and the proof is similar using mathematical induction. The second important lemma is given as follows.Lemma 7.
Given any two functionsA
(
x
,
y
) and B
(
x
,
y
) that are analytic (of Class C
∞) in their common domain, then(18)
e
ϕ
L
A
x
,
y
B
x
,
y
=
e
ϕ
L
A
x
,
y
e
ϕ
L
B
x
,
y
.Proof.
We first write(19)
e
ϕ
L
A
x
,
y
B
x
,
y
=
∑
n
=
0
∞
ϕ
n
n
!
L
n
A
x
,
y
B
x
,
y
=
∑
n
=
0
ϕ
n
n
!
∑
k
=
0
n
n
!
k
!
n
-
k
!
L
k
A
L
n
-
k
B
.Noting that 1
/
(
n
-
k
)
! vanishes for k
>
n we rewrite the above results as (20)
e
ϕ
L
A
x
,
y
B
x
,
y
=
∑
n
=
0
∞
∑
k
=
0
∞
ϕ
n
k
!
n
-
k
!
L
k
A
L
n
-
k
B
.Shifting the order of the sum and then shifting the index n
→
n
+
k, we write (21)
e
ϕ
L
A
x
,
y
B
x
,
y
=
∑
k
=
0
∞
∑
n
=
-
k
∞
ϕ
n
+
k
k
!
n
!
L
k
A
L
n
B
.Noting that 1
/
n
! vanishes for n
=
-
k
,
-
k
+
1
,
…
,
-
1, we reach the final conclusion (22)
e
ϕ
L
A
x
,
y
B
x
,
y
=
∑
k
=
0
∞
ϕ
k
k
!
L
k
A
∑
n
=
0
∞
ϕ
n
n
!
L
n
B
=
e
ϕ
L
A
x
,
y
e
ϕ
L
B
x
,
y
.Thus, (16) follows.
Based on (16) we can discuss the transformation of functions under rotation. A function Ψ
(
x
,
y
), which is analytic at (
0,0
), can be always expanded as(23)
Ψ
x
,
y
=
∑
n
,
m
a
n
m
n
!
m
!
x
n
y
m
.In the x
′
y
′ coordinates frame, the transformed function Ψ
′
(
x
′
,
y
′
) can also be written as (24)
Ψ
′
x
′
,
y
′
=
∑
n
,
m
a
n
m
′
n
!
m
!
x
′
n
y
′
m
=
e
ϕ
L
∑
n
,
m
a
n
m
′
n
!
m
!
x
n
y
m
.If the constants a
n
m are independent of the coordinates, then we develop the following important result.Remark 8.
A functionΨ
(
x
,
y
) that is analytic in a specific domain, where its dependence on the rotation angle ϕ comes only through its arguments (
x
,
y
), is transformed as(25)
Ψ
′
x
′
,
y
′
=
e
ϕ
L
Ψ
x
,
y
.
A scalar functionΨ
(
x
,
y
) that is invariant under spatial rotation must transform as Ψ
′
(
x
′
,
y
′
)
=
Ψ
(
x
,
y
); thus (26)
L
Ψ
x
,
y
=
x
D
y
1
Ψ
-
y
D
x
1
Ψ
=
0
,assuming the dependence on the rotation angle ϕ occurs only through the argument (
x
,
y
).Next we investigate the transformation of derivatives; one notes the following:(27)
D
x
′
1
x
′
n
=
n
x
′
n
-
1
=
n
e
ϕ
L
x
n
-
1
=
e
ϕ
L
D
x
1
x
n
=
e
ϕ
L
D
x
1
e
-
ϕ
L
e
ϕ
L
x
n
=
e
ϕ
L
D
x
1
e
-
ϕ
L
x
′
n
.Thus we derive the important result (28)
D
x
′
1
=
e
ϕ
L
D
x
1
e
-
ϕ
L
,and similarly (29)
D
y
′
1
=
e
ϕ
L
D
y
1
e
-
ϕ
L
.The following commutation relations are easy to derive: (30)
D
x
1
,
L
=
-
D
y
1
,
D
y
1
,
L
=
D
x
1
.Transformation of higher derivatives can be derived: for example,(31)
D
x
′
2
=
e
ϕ
L
D
x
1
e
-
ϕ
L
e
ϕ
L
D
x
1
e
-
ϕ
L
=
e
ϕ
L
D
x
2
e
-
ϕ
L
.Remark 9.
The partial derivatives transform under spatial rotation as(32)
D
x
′
n
=
e
ϕ
L
D
x
n
e
-
ϕ
L
,
D
y
′
n
=
e
ϕ
L
D
y
n
e
-
ϕ
L
,
D
x
′
n
D
y
′
m
=
e
ϕ
L
D
x
n
D
y
m
e
-
ϕ
L
,where n
,
m
∈
N.
As an illustration, one can show that(33)
D
x
2
,
L
=
-
2
D
x
1
D
y
1
,
D
y
2
,
L
=
+
2
D
x
1
D
y
1
.Thus, (34)
D
x
2
+
D
y
2
,
L
=
0
,indicating that the Laplace operator D
x
2
+
D
y
2 is a scalar operator, as expected.
An important remark is to check the transformation of the operatorL itself under spatial rotation. Consider(35)
L
′
=
e
ϕ
L
L
e
-
ϕ
L
=
L
.Thus the operator L is a scalar itself, a result that can be verified easily. This is an important point as it allows for the definition of inverse transformation between the x
y and x
′
y
′ coordinate frames.


## 4. Transformation of the Fractional Derivative under Spatial Rotation

In this section we generalize the results of the previous section to fractional calculus, in particular to the Riemann-Liouville and Caputo definitions of the fractional derivatives. We consider only real functions and real values of the fractional orderα
>
0. This restriction limits the allowed domain space (
x
,
y
). This is clear if we consider the transformation of x
β and y
β, where β
∈
R. Consider(36)
e
ϕ
L
x
β
=
cos
⁡
ϕ
x
+
sin
⁡
ϕ
y
β
=
x
′
β
,
e
ϕ
L
y
β
=
-
sin
⁡
ϕ
x
+
cos
⁡
ϕ
y
β
=
y
′
β
.Therefore, we require x
,
y
>
0 and x
′
,
y
′
>
0, which is equivalent to only consider the allowed domain: (
x
,
y
)
:
x
,
y
>
0
∩
y
/
x
>
tan
⁡
ϕ. Despite the restriction on the domain, one can still study the transformation of the fractional derivatives and deduce their properties. An extension to the whole domain can be pursued in the complex plane once a clear understanding of such transformation is acquired in the restricted domain.Based on the results of the previous section, we have the following.Remark 10.
In the domainx
,
y
>
0
∩
y
/
x
>
tan
⁡
ϕ, the quantity x
β
y
λ transforms as (37)
x
′
β
y
′
λ
=
e
ϕ
L
x
β
e
ϕ
L
y
λ
=
e
ϕ
L
x
β
y
λ
,where β
,
λ
∈
R. Next we look at the transformation of the partial fractional derivatives. Using the Riemann-Liouville and Caputo definition of the fractional derivative we note the following: (38)
D
x
′
α
x
′
β
=
Γ
β
+
1
Γ
β
-
α
+
1
x
′
β
-
α
=
Γ
β
+
1
Γ
β
-
α
+
1
e
ϕ
L
x
β
-
α
=
e
ϕ
L
D
x
α
x
β
=
e
ϕ
L
D
x
α
e
-
ϕ
L
e
ϕ
L
x
β
=
e
ϕ
L
D
x
α
e
-
ϕ
L
x
′
β
.Remark 11.
The partial fractional derivatives according to the Riemann-Liouville and Caputo definitions transform as(39)
D
x
′
α
=
e
ϕ
L
D
x
α
e
-
ϕ
L
,
D
y
′
α
=
e
ϕ
L
D
y
α
e
-
ϕ
L
.
Considering an infinitesimal angle transformation,ϕ
≪
1, and keeping only linear terms of ϕ, we write (40)
D
x
′
α
≈
D
x
α
+
ϕ
L
,
D
x
α
,
D
y
′
α
≈
D
y
α
+
ϕ
L
,
D
y
α
.Remark 12.
Given an analytic functionf
(
x
,
y
), its partial fractional derivatives transform as (41)
D
x
′
α
f
′
x
′
,
y
′
=
e
ϕ
L
D
x
α
f
x
,
y
,
D
y
′
α
f
′
x
′
,
y
′
=
e
ϕ
L
D
y
α
f
x
,
y
.Considering the infinitesimal angle approximation, ϕ
≪
1, and keeping only linear terms of ϕ, we write (42)
D
x
′
α
f
′
x
′
,
y
′
≈
D
x
α
f
x
,
y
+
ϕ
L
D
x
α
f
x
,
y
,
(43)
D
y
′
α
f
′
x
′
,
y
′
≈
D
y
α
f
x
,
y
+
ϕ
L
D
x
α
f
x
,
y
.Thus we derived the transformation of the partial fractional derivatives according to the Riemann-Liouville and Caputo definitions. Next we consider the transformation of the partial fractional derivatives as applied to a scalar field Ψ
(
x
,
y
) and using both Riemann-Liouville and Caputo definitions.


### 4.1. Transformation of the Riemann-Liouville Partial Fractional Derivatives of a Scalar Field

Consider a scalar functionΨ
(
x
,
y
) that is invariant under spatial rotation; that is, L
Ψ
(
x
,
y
)
=
0, and then according to (42) (44)
D
x
′
α
Ψ
′
x
′
,
y
′
≈
D
x
α
Ψ
x
,
y
+
ϕ
L
D
x
α
Ψ
x
,
y
,where (45)
L
D
x
α
Ψ
′
x
′
,
y
′
=
y
D
x
1
D
x
α
Ψ
x
,
y
-
x
D
y
1
D
x
α
Ψ
x
,
y
.Since y and D
y
1 commute with D
x
α we can write the above equation as (46)
L
D
x
α
Ψ
′
x
′
,
y
′
=
D
x
1
D
x
α
y
Ψ
-
x
D
x
α
D
y
1
Ψ
.Rewriting the above equation using commutation relations, we find (47)
L
D
x
α
Ψ
′
x
′
,
y
′
=
D
x
1
,
D
x
α
y
Ψ
-
x
,
D
x
α
D
y
1
Ψ
+
D
x
α
y
D
x
1
-
x
D
y
1
Ψ
.The last term vanishes since L
Ψ
=
0. Thus we arrive at(48)
L
D
x
α
Ψ
x
,
y
=
D
x
1
,
D
x
α
y
Ψ
-
x
,
D
x
α
D
y
1
Ψ
.Note that, from the properties of the fractional derivative, D
x
m
D
x
α
=
D
x
m
+
α, while D
x
1
,
D
x
α
≠
0 except for α
∈
Z
+. We use the known relation [1–5] (49)
D
x
α
D
x
1
Ψ
x
,
y
=
D
x
α
+
1
Ψ
x
,
y
-
x
-
α
-
1
Ψ
0
,
y
Γ
-
α
.Thus, we conclude (50)
D
x
1
,
D
x
α
y
Ψ
x
,
y
=
+
y
Ψ
0
,
y
x
-
α
-
1
Γ
-
α
.Using Leibniz rule in (6), substituting f
(
x
)
=
Ψ
(
x
,
y
) and g
(
x
)
=
x, one can show that (51)
x
,
D
x
α
D
y
1
Ψ
x
,
y
=
-
α
D
x
α
-
1
D
y
1
Ψ
x
,
y
.Thus we conclude that (52)
D
x
′
α
Ψ
′
x
′
,
y
′
≈
D
x
α
Ψ
x
,
y
+
ϕ
α
D
x
α
-
1
D
y
1
Ψ
+
y
Ψ
0
,
y
x
-
α
-
1
Γ
-
α
.A similar relation holds for y
(53)
D
y
′
α
Ψ
′
x
′
,
y
′
≈
D
y
α
Ψ
x
,
y
-
ϕ
α
D
y
α
-
1
D
x
1
Ψ
+
x
Ψ
x
,
0
y
-
α
-
1
Γ
-
α
.Note that, for α
∈
Z
+, the last term in the above result vanishes. Also for 0
<
α
<
1, the term D
y
α
-
1 represents a fractional integral.Remark 13.
Forα
=
1 we recover the standard transformation of the gradient operator: (54)
D
x
′
1
Ψ
′
x
′
,
y
′
=
D
x
1
Ψ
x
,
y
+
ϕ
D
y
1
Ψ
x
,
y
,
D
y
′
1
Ψ
′
x
′
,
y
′
=
D
y
1
Ψ
x
,
y
-
ϕ
D
x
1
Ψ
x
,
y
.Remark 14.
Forα
=
2 we recover the known scalar invariance of the Laplacian operator(55)
D
x
′
2
+
D
y
′
2
Ψ
′
x
′
,
y
′
=
D
x
2
+
D
y
2
Ψ
x
,
y
.


### 4.2. Transformation of the Caputo Partial Fractional Derivatives of a Scalar Field

To derive the transformation of the Caputo partial fractional derivatives we consider0
<
α
<
1 and use the relation between the Riemann-Liouville and Caputo definitions, as expressed in (9). Also we use the results of (50) and (51). One can show that (56)
x
,
D
x
α
C
D
y
1
Ψ
=
x
,
D
x
α
D
y
1
Ψ
-
x
1
-
α
D
y
1
Ψ
0
,
y
Γ
1
-
α
,
D
x
1
,
D
x
α
C
y
Ψ
=
D
x
1
,
D
x
α
y
Ψ
-
x
-
α
-
1
y
Ψ
0
,
y
Γ
-
α
.Therefore, the Caputo partial fractional derivative transforms, under infinitesimal angle ϕ
≪
1, as (57)
D
x
′
α
C
Ψ
′
x
′
,
y
′
≈
D
x
α
C
Ψ
x
,
y
+
ϕ
α
D
x
α
-
1
D
y
1
Ψ
+
x
1
-
α
D
y
1
Ψ
0
,
y
Γ
1
-
α
,
(58)
D
y
′
α
Ψ
′
x
′
,
y
′
≈
D
y
α
C
Ψ
x
,
y
-
ϕ
α
D
y
α
-
1
D
x
1
Ψ
+
y
1
-
α
D
x
1
Ψ
x
,
0
Γ
1
-
α
.Note that α
-
1
<
0 and thus D
x
α
-
1 and D
y
α
-
1 represent fractional integrals.


## 4.1. Transformation of the Riemann-Liouville Partial Fractional Derivatives of a Scalar Field

Consider a scalar functionΨ
(
x
,
y
) that is invariant under spatial rotation; that is, L
Ψ
(
x
,
y
)
=
0, and then according to (42) (44)
D
x
′
α
Ψ
′
x
′
,
y
′
≈
D
x
α
Ψ
x
,
y
+
ϕ
L
D
x
α
Ψ
x
,
y
,where (45)
L
D
x
α
Ψ
′
x
′
,
y
′
=
y
D
x
1
D
x
α
Ψ
x
,
y
-
x
D
y
1
D
x
α
Ψ
x
,
y
.Since y and D
y
1 commute with D
x
α we can write the above equation as (46)
L
D
x
α
Ψ
′
x
′
,
y
′
=
D
x
1
D
x
α
y
Ψ
-
x
D
x
α
D
y
1
Ψ
.Rewriting the above equation using commutation relations, we find (47)
L
D
x
α
Ψ
′
x
′
,
y
′
=
D
x
1
,
D
x
α
y
Ψ
-
x
,
D
x
α
D
y
1
Ψ
+
D
x
α
y
D
x
1
-
x
D
y
1
Ψ
.The last term vanishes since L
Ψ
=
0. Thus we arrive at(48)
L
D
x
α
Ψ
x
,
y
=
D
x
1
,
D
x
α
y
Ψ
-
x
,
D
x
α
D
y
1
Ψ
.Note that, from the properties of the fractional derivative, D
x
m
D
x
α
=
D
x
m
+
α, while D
x
1
,
D
x
α
≠
0 except for α
∈
Z
+. We use the known relation [1–5] (49)
D
x
α
D
x
1
Ψ
x
,
y
=
D
x
α
+
1
Ψ
x
,
y
-
x
-
α
-
1
Ψ
0
,
y
Γ
-
α
.Thus, we conclude (50)
D
x
1
,
D
x
α
y
Ψ
x
,
y
=
+
y
Ψ
0
,
y
x
-
α
-
1
Γ
-
α
.Using Leibniz rule in (6), substituting f
(
x
)
=
Ψ
(
x
,
y
) and g
(
x
)
=
x, one can show that (51)
x
,
D
x
α
D
y
1
Ψ
x
,
y
=
-
α
D
x
α
-
1
D
y
1
Ψ
x
,
y
.Thus we conclude that (52)
D
x
′
α
Ψ
′
x
′
,
y
′
≈
D
x
α
Ψ
x
,
y
+
ϕ
α
D
x
α
-
1
D
y
1
Ψ
+
y
Ψ
0
,
y
x
-
α
-
1
Γ
-
α
.A similar relation holds for y
(53)
D
y
′
α
Ψ
′
x
′
,
y
′
≈
D
y
α
Ψ
x
,
y
-
ϕ
α
D
y
α
-
1
D
x
1
Ψ
+
x
Ψ
x
,
0
y
-
α
-
1
Γ
-
α
.Note that, for α
∈
Z
+, the last term in the above result vanishes. Also for 0
<
α
<
1, the term D
y
α
-
1 represents a fractional integral.Remark 13.
Forα
=
1 we recover the standard transformation of the gradient operator: (54)
D
x
′
1
Ψ
′
x
′
,
y
′
=
D
x
1
Ψ
x
,
y
+
ϕ
D
y
1
Ψ
x
,
y
,
D
y
′
1
Ψ
′
x
′
,
y
′
=
D
y
1
Ψ
x
,
y
-
ϕ
D
x
1
Ψ
x
,
y
.Remark 14.
Forα
=
2 we recover the known scalar invariance of the Laplacian operator(55)
D
x
′
2
+
D
y
′
2
Ψ
′
x
′
,
y
′
=
D
x
2
+
D
y
2
Ψ
x
,
y
.


## 4.2. Transformation of the Caputo Partial Fractional Derivatives of a Scalar Field

To derive the transformation of the Caputo partial fractional derivatives we consider0
<
α
<
1 and use the relation between the Riemann-Liouville and Caputo definitions, as expressed in (9). Also we use the results of (50) and (51). One can show that (56)
x
,
D
x
α
C
D
y
1
Ψ
=
x
,
D
x
α
D
y
1
Ψ
-
x
1
-
α
D
y
1
Ψ
0
,
y
Γ
1
-
α
,
D
x
1
,
D
x
α
C
y
Ψ
=
D
x
1
,
D
x
α
y
Ψ
-
x
-
α
-
1
y
Ψ
0
,
y
Γ
-
α
.Therefore, the Caputo partial fractional derivative transforms, under infinitesimal angle ϕ
≪
1, as (57)
D
x
′
α
C
Ψ
′
x
′
,
y
′
≈
D
x
α
C
Ψ
x
,
y
+
ϕ
α
D
x
α
-
1
D
y
1
Ψ
+
x
1
-
α
D
y
1
Ψ
0
,
y
Γ
1
-
α
,
(58)
D
y
′
α
Ψ
′
x
′
,
y
′
≈
D
y
α
C
Ψ
x
,
y
-
ϕ
α
D
y
α
-
1
D
x
1
Ψ
+
y
1
-
α
D
x
1
Ψ
x
,
0
Γ
1
-
α
.Note that α
-
1
<
0 and thus D
x
α
-
1 and D
y
α
-
1 represent fractional integrals.


## 5. Discussion and Conclusions

In Section4, we derived the transformation of the Riemann-Liouville and Caputo partial fractional derivatives under spatial rotation. We also derived the transformation of the Riemann-Liouville and Caputo partial fractional derivatives of a scalar function. We observe no fundamental difference between the transformations of the Riemann-Liouville and Caputo partial fractional derivatives. The fact that the Caputo fractional derivative of a constant vanishes does not lead to a significant difference in its transformation. The difference we observe is that, in the Riemann-Liouville transformation, there is a term that depends on the value of the function at the boundary, while, for the Caputo transformation, a similar term depends on the normal derivative of the function at the boundary. We observe that the Riemann-Liouville and Caputo partial fractional derivatives of a scalar function do not transform as the components of a gradient vector, as in the case of the integer-order, α
=
1, partial derivatives. Thus they are classified to have a new form of transformation, different from the scalar, vector, and tensor quantities. They represent new entities that deserve to be studied thoroughly in future work.Next we provide few examples of scalar functions. Consider the constant scalar functionf
(
x
,
y
)
=
1. According to (52) and (53), the Riemann-Liouville transformations of the fractional derivatives are, for infinitesimal ϕ
≪
1, (59)
D
x
′
α
1
≈
D
x
α
1
+
ϕ
y
x
-
α
-
1
Γ
-
α
,
D
y
′
α
1
≈
D
y
α
1
-
ϕ
x
y
-
α
-
1
Γ
-
α
.One can easily check that the result is as expected, for example, (60)
D
x
′
α
1
=
x
′
-
α
Γ
-
α
≈
x
+
ϕ
y
-
α
Γ
-
α
≈
D
x
α
1
+
ϕ
y
x
-
α
-
1
Γ
-
α
,as required. It is easy to show that the following quantities are invariant: (61)
x
α
D
x
α
1
,
y
α
D
y
α
1
,
y
-
α
D
x
α
1
-
x
-
α
D
y
α
1
,
x
2
+
α
D
x
α
1
+
y
2
+
α
D
y
α
1
.One can use (
d
x
)
α and (
d
y
)
α instead of x
α and y
α in the above invariant forms. The Caputo partial fractional derivatives of a constant field vanish.Consider the scalar functionΨ
(
x
,
y
)
=
x
2
+
y
2. The Riemann-Liouville and Caputo x-partial fractional derivatives are, respectively, (62)
D
x
α
Ψ
x
,
y
=
2
Γ
3
-
α
x
2
-
α
+
y
2
x
-
α
Γ
1
-
α
,
D
x
α
C
Ψ
x
,
y
=
2
Γ
3
-
α
x
2
-
α
.Again it is easy to check that the transformations of the partial fractional derivatives, for infinitesimal angle ϕ
≪
1, agree with (52) and (57) for both the Riemann-Liouville and Caputo fractional derivatives. A similar transformation holds for D
y
α and D
y
α
C with x and y interchanged and ϕ
→
-
ϕ. One can easily check that the quantities (63)
x
α
D
y
α
Ψ
x
,
y
+
y
α
D
x
α
Ψ
x
,
y
,
x
α
y
α
D
y
α
D
y
α
Ψ
x
,
y
,
(64)
x
α
D
y
α
C
Ψ
x
,
y
+
y
α
D
x
α
C
Ψ
x
,
yare invariant under spatial rotation. Again, one can use (
d
x
)
α and (
d
y
)
α instead of x
α and y
α in the above invariant forms.For the functionΨ
(
x
,
y
)
=
(
x
2
+
y
2
)
2, one can check the transformations correctly. However, the quantities in (63) are not invariant separately. Nevertheless, it is possible to find a linear combination of the two quantities to form an invariant quantity: namely, (65)
x
α
D
y
α
Ψ
x
,
y
+
y
α
D
x
α
Ψ
x
,
y
+
A
x
α
y
α
D
y
α
D
y
α
Ψ
x
,
y
,where A is some constant to make the combination invariant. Again, one can use (
d
x
)
α and (
d
y
)
α instead of x
α and y
α in the above invariant forms.In [42], the authors derived the infinitesimal transformation of the partial fractional derivatives of a scalar field, by using the function’s Taylor expansion and then applying the Riemann-Liouville fractional derivative to each term of the expansion. Their result agrees with our general result.



---

*Source: 719173-2015-07-27.xml*