Application of q-Calculus to the Solution of Partial q-Differential Equations

Abstract

We introduce the concept of q-calculus in quantum geometry. This involves the q-differential and q-integral operators. With these, we study the basic rules governing q-calculus as compared with the classical Newton-Leibnitz calculus, and obtain some important results. We introduce the reduced q-differential transform method (RqDTM) for solving partial q-differential equations. The solution is computed in the form of a convergent power series with easily computable coefficients. With the help of some test examples, we discover the effectiveness and performance of the proposed method and employing MathCAD 14 software for computation. It turns out that when q = 1, the solution coincides with that for the classical version of the given initial value problem. The results demonstrate that the RqDTM approach is quite efficient and convenient.

Share and Cite:

Sadik, M. and Orie, B. (2021) Application of q-Calculus to the Solution of Partial q-Differential Equations. Applied Mathematics, 12, 669-678. doi: 10.4236/am.2021.128047.

1. Introduction

Every physical theory is formulated in terms of mathematical objects. It is thus necessary to establish a set of rules to map physical concepts and objects into mathematical objects that we use to represent them. Often times this mapping is evident, as in classical mechanics, while for other theories, such as quantum mechanics, the mathematical objects are not so intuitive. The background of our study is based on the dynamics of closed quantum systems, with intrinsic components such as states, observables, measurements and evolution. Quantum geometry on the other hand can be traced to the early days of quantum mechanics. Specifically, if we consider Heisenberg’s commutation relations [1] [2]

[ x , p ] = i (1)

it becomes obvious that the geometry of classical phase space is completely lost. When coordinates such as x (position) and p (momentum) on a phase space cease to commute, then there can be no such space! Moreover, we discover that this operator algebra [3] [4] [5] forms some kind of noncommutative geometric space. This is in contrast to algebraic geometry [6] [7], which is built on a correspondence betweenspaces and commutative algebras. This correspondence in particular associates with any given space, the algebra of functions on it, and geometric notions are then expressed in a purely algebraic format. This principle turns out to be the most logical starting point for a generalized geometry such as quantum geometry. While for algebraic geometry, the spaces are affine schemes, a correspondence that is closer to differential geometry is given by the Gelfand-Naimark theorem [8]. In this case, the spaces are topological spaces and the algebras are commutative C*-algebras.

In 2014, Maliki et al. [9], discussed the notion of q-deformed calculus in quantum geometry. Here they showed that the mathematical study of noncommutative geometry is intimately related to the so-called q-calculus, which is a generalization of the Newton-Leibnitz classical calculus.

In this work, we summarize some important q-calculus results which will enable us to study non-commutative differential equations, specifically we shall employ the reduced q-differential transform method (RqDTM) to solve partial q-differential equations.

1.1. The q-Differential Operator

For 1 q , we define the q-differential operator D q as;

D q f ( x ) = f ( x ) f ( q x ) ( 1 q ) x (2)

Note that D q d / d x D as q 1 . We make the following remarks.

1) At the beginning of the 20th century, F.H. Jackson [10] studied this modified derivative and many of its consequences.

2) This q-derivative can be applied to function not containing 0 in their domain of definition. Then it reduces to the ordinary derivative when q goes to 1.

3) One can easily check that the q-derivative operator is linear, i.e.,

(a) D q ( f + g ) = D q f + D q g (b) D q ( λ f ) = λ D q f (3)

1.2. Basic Notions of q-Calculus

The mathematical study of non-commutative geometry is intimately related to the so-called q-calculus. We begin with the q-differentials of some standard functions, then we give a brief introduction to q-numbers, q-factorials and q-integrals.

1.2.1. q-Derivative of Some Standard Functions

Following the procedure for computing the q-derivative from first principles, we now obtain important results for the q-derivative of the following standard functions such sin x , and e x .

1.2.2. q-Derivative of the Function f ( x ) = s i n x

By definition, we have

D q ( sin x ) = sin x sin q x ( 1 q ) x = ( x 1 3 ! x 3 + 1 5 ! x 5 1 7 ! x 7 + ) ( q x 1 3 ! q 3 x 3 + 1 5 ! q 5 x 5 1 7 ! q 7 x 7 + ) ( 1 q ) x = 1 1 3 ! x 2 ( 1 q 3 ) ( 1 q ) + 1 5 ! x 4 ( 1 q 5 ) ( 1 q ) 1 7 ! x 6 ( 1 q 7 ) ( 1 q ) + (4)

**Note: for typesetting use Dq(sinx) below, similarly for D(sinx).

Using the fact that 1 q r 1 q = 1 + q + q 2 + + q r 1

D q sin x = 1 1 3 ! x 2 ( 1 + q + q 2 ) + 1 5 ! x 4 ( 1 + q + q 2 + q 3 + q 4 ) 1 7 ! x 6 ( 1 + q + q 2 + q 3 + q 4 + q 5 + q 6 ) +

When q = 1 , we obtain the classical derivative of sin x function, i.e.

D sin x = 1 1 3 ! x 2 ( 3 ) + 1 5 ! x 4 ( 4 ) 1 7 ! x 6 ( 6 ) + = 1 1 2 ! x 2 + 1 4 ! x 4 1 6 ! x 6 + = cos x

1.2.3. q-Derivative of the Exponential Function f ( x ) = e x

We have by definition;

D q e x = e x e q x ( 1 q ) x = e x ( 1 e x ( 1 q ) ) ( 1 q ) x (5)

Since e u = n = 0 u n n ! , letting u = ( 1 q ) x we have;

e x ( 1 q ) = n = 0 ( x ( 1 q ) ) n n ! = n = 0 ( 1 ) n ( 1 q ) n x n n ! = 1 ( 1 q ) x + ( 1 q ) 2 x 2 2 ! + + ( 1 ) n ( 1 q ) n x n n ! +

1 e x ( 1 q ) = ( 1 q ) x ( 1 q ) 2 x 2 2 ! + ( 1 q ) 3 x 3 3 ! + ( 1 ) n + 1 ( 1 q ) n x n n ! +

D q e x = e x ( 1 e x ( 1 q ) ) ( 1 q ) x = e x [ 1 ( 1 q ) x 2 ! + ( 1 q ) 2 x 2 3 ! + ( 1 ) n ( 1 q ) n 1 x n 1 n ! + ] (6)

From the above expression it is clear that when q = 1 , we obtain the classical exponential function. We remark that the q-derivative of other standard functions can be obtained similarly.

1.2.4. q-Numbers and q-Factorials

We adopt basically the notations in [11]. Thus the set of positive integers is denoted by . Furthermore, throughout this paper K denotes a field of characteristic 0 and K ( q ) denotes the field of rational functions in one parameter q over K ( q ) . K ( q ) is our base field in the q-deformed setting, while K is the base field in the classical setting. We define respectively the q-integers, q-factorials and q-binomials as follows:

1) n q : = i = 0 n 1 q i = 1 q n 1 q (7)

2) n q ! : = 1 q 2 q n q ; 0 q ! = 1 (8)

3) { n m } q : = n q ! m q ! n m q ! , n , m { 0 } , n m (9)

For example, given f ( x ) = x n then D q f ( x ) = n x n 1 .

1.2.5. Remark

The properties of the q-integers, q-factorials and q-binomials and their proofs are presented in [9]. We also have the following important results on the q-differential operator. For 1 q , and with D q as defined previously;

1) D q f ( x ) = r = 0 ( q 1 ) r ( r + 1 ) ! x r d r + 1 d x r + 1 f ( x ) (10)

2) D q n f ( x ) = x n q n ( n 1 ) / 2 ( q 1 ) n j = 0 n { n j } q ( 1 ) j q j ( j 1 ) / 2 f ( q n j x ) (11)

3) D q { u ( x ) v ( x ) } = v ( x ) D q u ( x ) + u ( q x ) D q v ( x ) (12)

4) D q { u ( x ) / v ( x ) } = v ( q x ) D q u ( x ) u ( q x ) D q v ( x ) v ( x ) v ( q x ) (13)

1.2.6. The q-Integral Operator

A function F ( x ) is a q-antiderivative of f ( x ) if D q F ( x ) = f ( x ) . It is denoted by F ( x ) = f ( x ) d q x and called the Jackson integral [10]. We make the following remarks.

· Similar to classical integral calculus, any given function has multiple q-antiderivatives.

· Though the q-antiderivative of a function might not be unique, it can prove that if q ( 0 , 1 ) , a function has up to an additive constant oneq-antiderivative which is continuous at x = 0 [12].

· In [8] it is shown that;

F ( x ) = f ( x ) d q x = ( 1 q ) x k = 0 q k f ( q k x ) (14)

1.2.7. Properties of the q-Integral

1) The Jackson integral gives a q-antiderivative F ( x ) which is continuous at x = 0 , and is unique up to additive constant.

2) Generally, given g ( x ) is another function and D q g ( x ) denotes its q-derivative, we have formally

f ( x ) D q g ( x ) d q x = ( 1 q ) x r = 0 q r f ( q r x ) D q g ( q r x ) = ( 1 q ) x r = 0 q r f ( q r x ) g ( q r x ) g ( q r + 1 x ) ( 1 q ) q r x

f ( x ) d q g ( x ) = r = 0 f ( q r x ) ( g ( q r x ) g ( q r + 1 x ) ) (15)

giving a q-analogue of the Riemann–Stieltjes integral [12]. As an example of the foregoing, let f ( x ) = x n , n + . We have

x n d q x = ( 1 q ) x j = 0 q j q j n x n = ( 1 q ) x n + 1 j = 0 q j ( n + 1 ) = ( 1 q ) 1 q n + 1 x n + 1 = x n + 1 n + 1 q (16)

3) The integration by parts formula of Newton-Leibnitz calculus is interpreted in the present non-commutative context as;

f 2 ( x ) D q f 1 ( x ) d q x = f 1 ( x ) f 2 ( x ) f 1 ( q x ) D q f 2 ( x ) d q x (17)

1.3. Partial q-Derivative of a Multivariable Function

We define the partial q-derivative of a multivariable real continuous function f ( x 1 , x 2 , , x n ) with respect to a variable x i by;

D q , x i f ( x ) = ( ε q , i f ) ( x ) f ( x ) ( q 1 ) x i , x i 0 , q ( 0 , 1 ) (18)

[ D q , x i f ( x ) ] x i = 0 = lim x i 0 D q , x i f ( x ) (19)

where x = ( x 1 , x 2 , , x n ) and, ( ε q , i f ) ( x ) = f ( x 1 , x 2 , , q x i , , x n ) .

We adopt subsequently the identity D q , x k q k q x k for the kth order q-derivative with respect to x k .

Our objective in this section is to solve partial q-differential equations using a novel q-differential transform method introduced in [13].

1) Reduced q-Differential Transform Method (RqDTM)

Given that all q-differentials of u ( x , t ) exist in some neighborhood of t = a , then let

U k ( x ) = 1 k q ! [ q k k t k u ( x , t ) ] t = a (20)

where the t-dimensional spectrum function U k ( x ) is the transformed function. Subsequently, the lowercase u ( x , t ) represents the original function while the uppercase U k ( x ) stands for the transformed function. We have the following important definition.

2) Definition

The q-differential inverse transform of U k ( x ) is defined by;

u ( x , t ) = k = 0 U k ( x ) ( t a ) ( k ) (21)

Substituting Equations (20) in (21) we obtain

u ( x , t ) = k = 0 1 k q ! [ q k k t k u ( x , t ) ] t = a ( t a ) ( k ) (22)

In the subsequent theorems, we set a = 0 so that ( t a ) ( k ) = ( t 0 ) ( k ) = t k . From the linearity of the q-derivative, we can establish the following fact; given w ( x , t ) = α u ( x , t ) ± v ( x , t ) then W k ( x ) = α U k ( x ) ± V k ( x ) , α being a constant. We have the following important theorem.

3) Theorem: Given w ( x , t ) = x m t n then W k ( x ) = x m δ ( k n ) where

δ ( k ) = { 1 , k = 0 0 , k 0 (23)

Proof: By definition (20), we have;

W k ( x ) = 1 k q ! [ q k ( x m t n ) q t k u ( x , t ) ] t = 0 = x m k q ! [ q k ( x m t n ) q t k u ( x , t ) ] t = 0 = { x m k q ! k q ! = x m , k = n x m n q ( n 1 ) q ( n k + 1 ) q k q ! t n k | t = 0 , k 0 x m 0 = 0 , k > n = x m δ ( k n ) (24)

4) Theorem: Given w ( x , t ) = q q x u ( x , t ) then W k ( x ) = q q x U ( x ) .

Proof

W k ( x ) = 1 k q ! [ q k q t k ( q q x u ( x , t ) ) ] = 1 k q ! [ q q x ( q k q t k u ( x , t ) ) ] t = 0 = q q x [ 1 k q ! q k q t k u ( x , t ) ] t = 0 = q q x U k ( x ) (25)

5) Theorem: Given w ( x , t ) = q r q x u ( x , t ) then

W k ( x ) = k + 1 q k + 2 q k + r q U k + r ( x ) (26)

Proof

W k ( x ) = 1 k q ! [ q k q t k ( q r q t r u ( x , t ) ) ] = k + r q ! k q ! 1 k + r q ! ( q k + r q t k + r u ( x , t ) ) t = 0 = k + 1 q k + 2 q k + r q U k + r ( x )

6) Example

q q t u ( x , t ) = u 2 ( x , t ) + q q x u ( x , t ) , u ( x , 0 ) = 1 + 3 x (27)

Taking the RqDTM of the given partial q-differential equation, we have

k + 1 q U k + 1 ( x ) = n = 0 k U k n ( x ) U k ( x ) + q q x U k ( x ) , k = 0 , 1 , 2 , (28)

The initial condition becomes

U 0 ( x ) = u ( x , 0 ) = 1 + 3 x .

Starting with k = 0 , the values of U k ( x ) are computed successively as follows;

1 q U 1 ( x ) = U 0 ( x ) U 0 ( x ) + q q x U 0 ( x ) = ( 1 + 3 x ) 2 + q q x ( 1 + 3 x ) = ( 1 + 3 x ) 2 + 3 U 1 ( x ) = 4 + 6 x + 9 x 2 (29)

When k = 1 , we have

2 q U 2 ( x ) = 2 U 0 ( x ) U 1 ( x ) + q q x U 1 ( x ) = 2 ( 1 + 3 x ) ( 4 + 6 x + 9 x 2 ) + q q x ( 4 + 6 x + 9 x 2 ) = 2 ( 4 + 6 x + 9 x 2 + 12 x + 18 x 2 + 27 x 3 ) + ( 6 + 9 x ( 1 + q ) )

U 2 ( x ) = 14 + 9 ( 5 + q ) x + 54 x 2 + 54 x 3 1 + q (30)

Following the same procedure, it is easy to compute an expression for k = 2 .

Formally, we have the required solution of the partial q-differential equation to be

u ( x , t ) = 1 + 3 x + ( 4 + 6 x + 9 x 2 ) t + ( 14 + 9 ( 5 + q ) x + 54 x 2 + 54 x 3 1 + q ) t 2 + (31)

Let us now consider the classical version of the given partial q-differential equation, namely

t u ( x , t ) = u 2 ( x , t ) + x u ( x , t ) , u ( x , 0 ) = 1 + 3 x (32)

The above is a first order quasilinear partial differential equation easily solved by the method of characteristics. The associated auxiliary equations are;

d t 1 = d x 1 = d u u 2 (33)

These provide us with two possible integrals given by;

t + x = c 1 and t + 1 u = c 2 (34)

c 1 , c 2 being arbitrary constants of integration. Using the initial condition, we have at t = 0 , u = 1 + 3 x . Hence, x = c 1 and 1 1 + 3 x = c 2 . It then follows that 1 1 + 3 c 1 = c 2 .

Hence the required solution is

t + 1 u = 1 1 + 3 ( t + x ) or u ( t , x ) = 1 + 3 ( t + x ) 1 t ( 1 + 3 x ) 3 t 2 (35)

Using MathCAD14, a numerical algebra software [14] to expand the expression for u in terms of powers of t, we obtain

u ( t , x ) = ( 1 + 3 x ) + ( 4 + 6 x + 9 x 2 ) t + ( 7 + 27 x 2 + 27 x 3 ) t 2 + (36)

We make the interesting observation that when we set q = 1 in (31) we obtain exactly the solution to the classical PDE.

7) Example. As an example of a second order partial q-differential equation we consider the q-diffusion Cauchy problem

q q t u ( x , t ) = q 2 q x 2 u ( x , t ) , u ( x , 0 ) = e q x (37)

Again we employ the RqDTM transform. Here e q x in the initial data, is the q-exponential function defined by e q x = i = 0 x i i q ! , with q-derivative e q x . The RqDTM transform of the q-diffusion equation, gives the following recursion;

k + 1 q U k + 1 ( x ) = q 2 q x 2 U k ( x ) , k = 0 , 1 , 2 , (38)

The initial data given is then written

U 0 ( x ) = u ( x , 0 ) = e q x (39)

Now, substituting (38) into (37), we obtain the following U k ( x ) values successively

U 1 ( x ) = 1 1 q ! e q x , U 2 ( x ) = 1 2 q ! U 1 ( x ) = 1 1 q ! 2 q ! e q x = 1 2 q ! e q x (40)

U 3 ( x ) = 1 3 q ! e q x , , U k ( x ) = 1 k q ! e q x (41)

In view of (28), the differential inverse transform of U k ( x ) gives

u ( x , t ) = k = 0 U k t k = k = 0 1 k q ! e q x t k = e q x e q t (42)

which is the analytic solution of the problem (36). Now consider the classical diffusion equation is written

t u ( x , t ) = 2 x 2 u ( x , t ) , u ( x , 0 ) = e x (43)

By the method of separation of variables, it is easy to show that the solution is given by

u ( x , t ) = X ( x ) T ( t ) = e x e t = e x + t (44)

Comparing with the q-solution u ( x , t ) = e q x e q t , we see that the two solutions are in perfect agreement when q = 1, and satisfying their appropriate initial conditions.

1.4. Conclusion

In this research article, we introduced the concept of q-calculus in quantum geometry. This involves the study of the basic rules governing q-calculus as compared with the classical Newton-Leibnitz calculus. Our main objective is to employ the results obtained to solve partial q-differential equations. To this end we introduced the reduced q-differential transform method which provides the solution in the form of a convergent power series with easily computable components. With the help of a few examples, we were able to show that the proposed iteration technique is very effective and convenient. It turns out that when q = 1, the solution coincides with the classical version of the given initial value problem. In conclusion, q-calculus is a non-commutative calculus that generalizes the Newton-Leibnitz classical calculus. The Reduced Differential Transform method for solving differential equations was introduced and extended in this work to solve partial q-differential equations, which represent some form of dynamics in non-commutative spaces.

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

References

[1] Laragno, A. and Gervino, G. (2009) Quantum Mechanics in q-Deformed Calculus. Journal of Physics Conference Series, 174, 012071.
https://doi.org/10.1088/1742-6596/174/1/012071
[2] Connes, A. (1994) Noncommutative Geometry. Academic Press, New York.
[3] Connes, A. (1986) Non-Commutative Differential Geometry. Publications Mathématiques de l’Institut des Hautes Études Scientifiques, 62, 41-144.
https://doi.org/10.1007/BF02698807
[4] Connes, A. (1985) Noncommutative Differential Geometry. Publications Mathématiques de l’I.H.É.S, No. 62, 257-360.
[5] Connes, A. (1994) Noncommutative Geometry. Academic Press, London, 126.
[6] Majid, S. (1990) Physics for Algebraists: Non-Commutative and Non-Commutative Hopf Algebras by a Bicrossproduct Construction. Journal of Algebra, 130, 17-64.
[7] Drinfeld, G. (1986) Quantum Groups. In: Gleason, A., Ed., Proceedings of the ICM 1986, AMS, Rhode Island, 798-820.
[8] Brown, L.G., Douglas, R.G. and Filmore, P.G. (1977) Extensions of C*-Algebras and K-Homology. Annals of Mathematics, 105, 265-324.
https://doi.org/10.2307/1970999
[9] Maliki, O.S. and Ugwu, E.I. (2014) Onq-Deformed Calculus in Quantum Geometry. Applied Mathematics, 5, 1586-1593.
https://doi.org/10.4236/am.2014.510151
[10] Jackson, H.F. (1910) q-Difference Equations. American Journal of Mathematics, 32, 305-314.
https://doi.org/10.2307/2370183
[11] Connes, A. (1980) C*-algèbres et géométrie différentielle. Comptes Rendus de l’Académie des Sciences, 290, 599-604.
[12] Edward, N. (1974) Notes on Non-Commutative Integration. Journal of Functional Analysis, 15, 103-116.
https://doi.org/10.1016/0022-1236(74)90014-7
[13] Jafari, H., Haghbin, A., Hesam, S. and Baleanu, D. (2014) Solving Partial q-Differential Equations within Reduced q-Differential Transformation Method. Romanian Journal of Physics, 59, 399-407.
[14] Mathcad Version 14 (2007) PTC (Parametric Technology Corporation) Softwar Products.
http://communications@ptc.com

Copyright © 2024 by authors and Scientific Research Publishing Inc.

Creative Commons License

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.