Next Article in Journal
Euclidean Space Controllability Conditions for Singularly Perturbed Linear Systems with Multiple State and Control Delays
Previous Article in Journal
Fixed Point Theorems for Geraghty Contraction Type Mappings in b-Metric Spaces and Applications
Previous Article in Special Issue
Harrod–Domar Growth Model with Memory and Distributed Lag
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Operator Ordering and Solution of Pseudo-Evolutionary Equations

1
Institut de Recherche en Informatique Fondamentale (IRIF), Université Paris-Diderot, F-75013 Paris, France
2
ENEA—Frascati Research Center, Via Enrico Fermi 45, 00044 Rome, Italy
3
H. Niewodniczański Institute of Nuclear Physics, Polish Academy of Science, 31-342 Kraków, Poland
*
Author to whom correspondence should be addressed.
Axioms 2019, 8(1), 35; https://doi.org/10.3390/axioms8010035
Submission received: 1 February 2019 / Revised: 8 March 2019 / Accepted: 13 March 2019 / Published: 16 March 2019
(This article belongs to the Special Issue Fractional Differential Equations)

Abstract

:
The solution of pseudo initial value differential equations, either ordinary or partial (including those of fractional nature), requires the development of adequate analytical methods, complementing those well established in the ordinary differential equation setting. A combination of techniques, involving procedures of umbral and of operational nature, has been demonstrated to be a very promising tool in order to approach within a unifying context non-canonical evolution problems. This article covers the extension of this approach to the solution of pseudo-evolutionary equations. We will comment on the explicit formulation of the necessary techniques, which are based on certain time- and operator ordering tools. We will in particular demonstrate how Volterra-Neumann expansions, Feynman-Dyson series and other popular tools can be profitably extended to obtain solutions for fractional differential equations. We apply the method to several examples, in which fractional calculus and a certain umbral image calculus play a role of central importance.

1. Introduction

We discuss in this paper the notion of so-called pseudo-evolutionary differential equations, which we define as equations of the Cauchy type,
D ^ t F ( x , t ) = α O ^ x F ( x , t ) F ( x , 0 ) = f ( x ) .
Here, D ^ t is an operator playing the role of the derivative with respect to the variable t, α is a constant, O ^ x is an operator of differential or pseudo-differential nature (to be specified later in our discussion), and where f ( x ) is the initial condition.
The development of suitable solution methods, mimicking those exploited for the ordinary Cauchy problem setting, requires the following steps:
  • Finding an eigenfunction of the D ^ t operator, such that
    D ^ t E ( λ t ) = λ E ( λ t ) .
  • Constructing a pseudo evolution operator (PEO) as
    U ^ ( t ) = E ( α t O ^ x ) ,
    which leads to the solution of Equation (1) via
    F ( x , t ) = E ( α t O ^ x ) f ( x ) .
  • Establishing rules that permit the explicit evaluation of the action of the PEO U ^ ( t ) on the initial function f ( x ) in the formal solution Equation (4).
Consequently, the precise properties of the eigenfunctions E that arise in the first step of this procedure are of crucial importance in terms of developing explicit calculation techniques. It is evident that the strategy we have in mind is that of treating E as a kind of exponential function (the eigenfunction of the ordinary derivative operator)—in this way, we can recover most of the techniques associated with operator disentanglement and time-ordering (whenever necessary).
A first study in this direction has been addressed by D. Babusci and one of the present authors in [1], which has later been specialized in several papers for situations in which D ^ t is a fractional [2] or Laguerre [3] derivative. In elaborating this program, we have been faced with some difficulties associated with the fact that the corresponding eigenfunctions E lack the semi-group property, namely
E ( x + y ) E ( x ) E ( y ) = E ( y ) E ( x ) .
This is an additional source of difficulties also when dealing with operator-ordering in pseudo-evolutionary problems. If x ^ , y ^ are non-commuting operators, such that
x ^ , y ^ 0 ,
it is well known that even for the exponential function one finds that [4]
e x ^ + y ^ e x ^ e y ^ e y ^ e x ^ .
Instead, one finds the following types of expansions [5,6]
e x ^ + y ^ = e x ^ e y ^ e f ^ 1 e f ^ n = e y ^ e x ^ e g ^ 1 e g ^ n ,
where f ^ n , g ^ n are expressed in terms of chains of commutators involving x ^ , y ^ , as it happens e.g., for the Zassenhaus expansion [7]. It is evident that in the special case x ^ , y ^ = 0 , the f ^ n and g ^ n operator functions vanish, and the semi-group property is restored.
Time-ordering is a further element of complication, which may arise in evolutionary problems (ordinary or pseudo) whenever the operator O ^ x on the rhs of Equation (1) is explicitly time-dependent and does not commute with itself at different times.
In this paper, we reconsider the operator-ordering problem for pseudo-evolutionary problems in more general terms than those considered in [1,2,3], and we address the time-ordering problem by discussing the possibility of adapting for this purpose the use of expansions such as the Volterra-Neumann or Feynman-Dyson series [8]. The paper is organized as follows: In Section 2, we summarize and extend the results contained in [2,3]. In Section 3, we provide a general view on the problem of time-ordering for pseudo evolution equations, while Section 4 contains several specific examples and final comments.

2. Laguerre Derivative, Laguerre Exponential and Operator-Ordering

The Laguerre transform and the associated operational calculus have played a crucial role in the theory of operational calculus [9,10]. They have offered elements of key importance within the context of the monomiality theory [11] and for the study of integro-differential equations of Volterra type [12]. We will outline the procedure allowing the merging of ordering procedures and umbral image type methods, using as a reference example the so-called “Laguerre-calculus” along the lines of [1].
We introduce the forthcoming discussion by going back to the paradigmatic strategy sketched in the introductory section and, accordingly, fix the following specific steps:
  • We specialize the operators in Equation (1) to
    D ^ t = l t = t t t , O ^ x = x ,
    where l t is the Laguerre derivative [3,10,11,13].
  • The eigenfunction of the Laguerre derivative operator is the Bessel-like function l e ( x ) [14],
    l e ( x ) = r 0 x r ( r ! ) 2 ,
    which satisfies D ^ t l e ( λ t ) = λ l e ( λ t ) .
  • In view of explicit computations, it will prove advantageous to express l e ( x ) via an umbral image [15,16] (where we refer to Appendix A for the explicit definition of the full formalism)
    l e ( x ) = I ^ v e v x .
    Here, v is a formal variable, and I ^ a formal integration operator, which acts according to
    I ^ v α : = 1 Γ ( α ) ( α C ) .
We can therefore write the solution of our problem as
F ( x , t ) = I ^ v e α v t x f ( x ) = I ^ v f ( x + α v t ) ,
where we have just adopted the properties of the exponential shift operator (i.e., Taylor’s formula). This illustrates one particularly simple scenario in which the approach sketched in the introduction may be explicitly carried out.
The next example addresses the problem of operator-ordering. We consider a Laguerre-type evolution problem (with x ^ the operator of multiplication by x, i.e., x ^ ( x n ) : = x n + 1 ,
O ^ x = ( α x ^ β x ) ,
in which the novelty and the difficulty stems from the fact that it consists of the sum of two terms that are not commuting with each other (due to [ x , x ^ ] = 1 0 ). The solution of our problem can be cast in a first step into the form
F ( x , t ) = I ^ v e v t ( α x ^ β x ) f ( x ) .
However, the PEO in this expression cannot be disentangled into the product of two exponentials, because the operators in the argument of the exponential do not commute. We thus proceed as follows:
  • We define the auxiliary operators
    X ^ : = α v t x ^ , Y ^ : = β v t x .
  • Applying the Weyl disentanglement rule (taking advantage of the fact that [ [ X ^ , Y ^ ] , X ^ ] = [ [ X ^ , Y ^ ] , Y ^ ] = 0 ), we find that
    e X ^ + Y ^ = e 1 2 [ X ^ , Y ^ ] e X ^ e Y ^ .
  • We then eventually arrive at the closed-form expression
    F ( x , t ) = I ^ v e ( v t ) 2 2 α β e v t α x ^ e v t β x f ( x ) = I ^ v e ( v t ) 2 2 α β e v t α x ^ f ( x + v β t ) .
It important to emphasize that the operational ordering as performed above thus brings into play a term depending on the square of the formal variable v (which commutes with the differential operators x ^ and x ). Assuming for simplicity f ( x ) = 1 , we find that
F ( x , t ) = I ^ v e ( v t ) 2 2 α β e v t α x .
Finally, after Taylor-expanding the exponential and evaluating the action of the formal integral operator I ^ according to Equation (12), we find the closed-form expression
F ( x , t ) = I ^ v n 0 ( v t α ) n n ! e ( v t ) 2 2 α β = n 0 ( α t ) n n ! l e ( 2 ) α β t 2 2 ,
where l e n ( m ) ( x ) is the Bessel-like function defined as
l e n ( m ) ( x ) : = I ^ v n + 1 e x v m = r 0 x r r ! Γ ( m r + n + 1 ) .
The example we have discussed is sufficient to demonstrate that the umbral image formalism naturally yields the solution of evolution problems involving a Laguerre derivative and non-commuting operators.
In order to stress the generality and the flexibility of our method, we consider the further example concerning the Schrödinger-type equation
i l t Ψ ( x , t ) = O ^ x Ψ ( x , t ) , O ^ x = α x + β 2 x 2 .
We will derive the PEO associated with Equation (22) in complete analogy to the procedure discussed in the previous example (see also Appendix B for an alternative approach using a Berry-type formula). As a preparatory step, let us recall for the readers’ convenience the Zassenhaus formula in its “right- and left-oriented” forms (see e.g., [7]), whereby for a formal variable λ and for two composable linear operators X and Y one has
e λ ( X + Y ) = e λ X e λ Y e λ 2 C 2 ( X , Y ) e λ 3 C 3 ( X , Y ) e λ 4 C 4 ( X , Y ) e λ ( X + Y ) = e λ 4 C ^ 4 ( X , Y ) e λ 3 C ^ 3 ( X , Y ) e λ 2 C ^ 2 ( X , Y ) e λ Y e λ X C 2 ( X , Y ) = 1 2 [ Y , X ] , C 3 ( X , Y ) = 1 3 [ C 2 ( X , Y ) , X + 2 Y ] C n ( X , Y ) = 1 n ! d n d λ n e λ n 1 C n 1 ( X , Y ) e λ 3 C 3 ( X , Y ) e λ 2 C 2 ( X , Y ) e λ Y e λ X e λ ( X + Y ) | λ 0 ( n 3 ) C ^ n ( X , Y ) = ( 1 ) n + 1 C n ( X , Y ) ( n 2 ) .
Introducing the auxiliary operators
A ^ : = λ x 2 , B ^ : = κ x ^ , λ : = i 2 β v t , κ : = i α v t ,
we find the commutation relations
A ^ , B ^ = 2 κ λ x , [ A ^ , [ A ^ , B ^ ] ] = 0 , [ B ^ , [ A ^ , B ^ ] ] = 2 κ 2 λ ,
with all higher nested commutators vanishing. Thus invoking the “left-oriented” form of the formula (setting X = A ^ , Y = B ^ and λ = 1 ) and using the commutation relations Equation (25), we obtain
e A ^ + B ^ = e 1 3 ! ( [ A ^ , [ A ^ , B ^ ] ] + 2 [ B ^ , [ A ^ , B ^ ] ] ) e 1 2 [ A ^ , B ^ ] e B ^ e A ^ = ( * ) e 2 3 κ 2 λ e κ λ x e κ x ^ e λ x 2 = e 2 3 κ 2 λ e κ ( x ^ + κ λ ) e κ λ x e λ x 2 = e 1 3 κ 2 λ e κ x ^ e κ λ x e λ x 2
Here, in the step marked ( * ) , we have taken advantage of the Crofton-Glaisher identity [17] (see also [11,18] Equation (I.3.17)), whereby for a formal power series f ( x ) and for an integer-valued parameter m 1 ,
e λ x m f ( x ^ ) = f x ^ + m λ x m 1 e λ x m
Note that we provide this identity in “operational form”, i.e., by using the formal multiplication operator x ^ , this expression is also valid when part of larger expressions.
Coincidentally, this identity also permits us to compute the action of the operational expression computed in Equation (26) on an initial condition Ψ ( x , 0 ) = φ ( x ) , resulting in
e 1 3 κ 2 λ e κ x ^ e κ λ x e λ x 2 φ ( x ) = e 1 3 κ 2 λ e κ x ^ φ ( x ^ + κ λ + 2 λ x ) 1 .
Here and throughout this paper, in expressions such as φ ( x ^ + κ λ + 2 λ x ) 1 , the occurrence of the symbol “1” entails that the expression is to be evaluated by expanding φ ( x ^ + κ λ + 2 λ x ) into normal-ordered form (i.e., into a series in the normal-ordered monomials x ^ r x s for r , s 0 ), followed by acting on 1 (which due to x s 1 = 0 for s > 0 in effect amounts to dropping all terms of the expansion involving non-zero powers of x ).
Combining this result with the explicit formula for the eigenfunctions l e ( x ) as provided in Equation (11), we thus finally arrive at the explicit solution for Ψ ( x , t ) evolving according to the pseudo-evolution Equation (22) with initial condition Ψ ( x , 0 ) = φ ( x ) :
Ψ ( x , t ) = I ^ v e 1 6 ( i v t ) 3 α 2 β e ( i v t ) α x ^ φ ( x ^ + α β 2 ( i v t ) 2 + ( i v t ) β x ) 1 .
Specializing for simplicity to the case φ ( x ) = 1 , Equation (29) evaluates to
Ψ ( x , t ) = I ^ v e ( i v t ) α x + 1 6 ( i v t ) 3 α 2 β .
Coincidentally, the expression obtained in the last step has an interesting formal meaning: consider the third order Hermite polynomials H n ( 3 ) ( x , y ) , which are defined as
H n ( 3 ) ( x , y ) : = e y x 3 x n = n ! r = 0 n 3 x n 3 r y r ( n 3 r ) ! r ! ,
and whose exponential generating function (EGF) reads
H ( 3 ) ( t ; x , y ) : = n 0 t n n ! H n ( 3 ) ( x , y ) = e y x 3 e t x = e t x + t 3 y .
One may in fact derive the explicit formula for this generating function directly via use of the Crofton-Glaisher identity Equation (27) combined with the results of Equation (26).
We thus recognize the occurrence of the above EGF as a term in Equation (30), which allows us to express Ψ ( x , t ) in the alternative form
Ψ ( x , t ) = I ^ v n 0 ( i v t ) n n ! H n ( 3 ) α x , α 2 β 6 = n 0 ( i t ) n ( n ! ) 2 H n ( 3 ) α x , α 2 β 6 .
The examples of this introductory section have shown that a judicious combination of our suggested pseudo evolution operator (PEO) method with various elements from the theory of generalized functions and umbral image type techniques results in a toolset that allows to deal with non-standard forms of partial differential equations efficiently.
Before closing this section, it is worth commenting on the role played by the concepts associated with the semi-group property (or, rather, lack thereof) of the Laguerre exponential. We follow the point of view of [2], where these problems have been systematically investigated. As may be verified via an explicit calculation, one finds that the Laguerre exponential does not satisfy the semi-group property e x + y = e x e y (for x , y commuting variables) of the ordinary exponential function, whence
l e ( x + y ) l e ( x ) l e ( y ) ,
but rather satisfies
l e ( x ) l e ( y ) = r 0 x r r ! 2 s 0 y s s ! 2 = n 0 1 n ! 2 x l y n .
Here, the symbol l denotes the composition rule
x l y n : = s = 0 n n s 2 x n s y s ,
thus yielding the so-called “Laguerre Newton binomial”. According to the previous identities we may redefine the semi-group property for the Laguerre exponential as
l e ( x ) l e ( y ) = l e ( x l y ) .
The recently introduced reformulation [16] of the umbral calculus framework in terms of umbral image type techniques permits to understand the calculations that lead to Equation (37) in a very direct manner: taking advantage of the identity (see Appendix A for further details)
I ^ u α v α = Γ ( α ) Γ ( α ) = 1 ( α C { 0 , 1 , 2 , } ) ,
we may compute a “Laguerre Newton binomial” type result as follows:
l e ( x ) l e ( y ) = I ^ v 1 v 2 e v 1 x + v 2 y = I ^ u v 1 v 2 v 3 e v 3 ( u ( v 1 x + v 2 y ) ) = I ^ u v 1 v 2 l e ( u ( v 1 x + v 2 y ) ) .
We thus indeed find that
I ^ u v 1 v 2 ( u ( v 1 x + v 2 y ) ) n = r = 0 n n r n ! x r y n r r ! ( n r ) ! = x l y n .
This form of “umbral image reshaping” will prove particularly useful when considering more complicated types of special functions in the sequel.

3. Pseudo-Evolutive Problems and Matrix Calculus

In this section, we will demonstrate an extension of the previously introduced pseudo-evolution equation techniques to a form of matrix calculus. The problem we wish to address is the search of a solution for equations of the type
l t ( x ) Y ̲ ( t ) = M ^ Y ̲ ( t ) , Y ̲ ( 0 ) = Y ̲ 0 ,
with l t = t t t the Laguerre-type time-derivative, and where M ^ and Y ̲ ( t ) , Y ̲ 0 denote an n × n matrix and n-element column vectors, respectively. We specialize our discussion to the case of a non-singular 2 × 2 matrix with eigenvalues λ ± . Following the paradigm of the PEO method introduced in Section 2, one may obtain a solution of Equation (41) in the form
Y ̲ ( t ) = U ^ ( t ) Y ̲ 0 , U ^ ( t ) = l e ( M ^ t ) = n = 0 ( M ^ t ) n ( n ! ) 2 .
By application of the Cayley-Hamilton theorem, we thus obtain
l e ( M ^ t ) = 1 λ + λ ( λ + 1 ^ M ^ ) l e ( λ t ) ( λ 1 ^ M ^ ) l e ( λ + t ) ) ,
where 1 ^ denotes the 2 × 2 unit matrix. Let us now consider for illustration a matrix with zero diagonal entries and imaginary eigenvalues, namely
M ^ = 0 α β 0 .
According to Equation (43), the corresponding PEO U ^ ( t ) can be reduced to the “pseudo rotation matrix”
l e ( M ^ t ) = l c ( α β t ) α β l s ( α β t ) β α l s ( α β t ) l c ( α β t ) ,
where l c ( x ) and l s ( x ) denote the Laguerre (co-)sine functions [14] defined as
l c ( x ) = l e ( i x ) + l e ( i x ) 2 , l s ( x ) = l e ( i x ) l e ( i x ) 2 i ( i 2 = 1 ) .
The relevant geometrical meaning differs from that of the ordinary circular functions and is illustrated in Figure 1, where we have plotted l s ( x ) against l c ( x ) in the region of the first negative and positive zeros of the Laguerre sine function (cf. Figure 2).
The functions in Equation (46) can be recognized as Bessel type functions (in particular as Kelvin ber, bei functions), and they satisfy the differential equations
l t 2 l c ( ω t ) = ω 2 l c ( ω t ) , l t 2 l s ( ω t ) = ω 2 l s ( ω t ) .
We also note that the Laguerre derivative satisfies the identity
l t α = t α t α t α ( α R ) .
Our next example illustrates a possible generalization of the PEO methods in a different direction, namely in the form of fractional evolutive problems such as
t μ Y ̲ ( t ) = M ^ Y ̲ ( t ) + t μ Γ ( 1 μ ) Y ̲ 0 , Y ̲ ( 0 ) = Y ̲ 0 ( 0 < μ < 1 ) .
Following the paradigm of the PEO method, the first step in solving Equation (49) consists in finding an eigenfunction of the fractional differential operator D ^ t = t μ . To this end, we recall the definition of the so-called Mittag-Leffler function (ML-f) [19] E α , β ( x ) (for α , β C ),
E α , β ( x ) = r 0 x r Γ ( α r + β ) = I ^ v β 1 x v α .
Here, we have yet again taken advantage of the formal integral operator I ^ as introduced in [16] (compare Equation (12)). An alternative useful expression for the ML-f E α , β ( x ) may be obtained via taking a suitable Laplace transform of Equation (50), whence
E α , β ( x ) = I ^ v β 0 e s e v α x d s .
For example, the alternative form Equation (51) for the ML-f allows to derive (via a suitable umbral image reshaping) a multiplicative law in a straightforward fashion:
E α , β ( x ) E α , β ( y ) = I ^ ( v 1 v 2 ) β 0 e s e v 1 α x + v 2 α y d s = I ^ v 3 β ( u v 1 v 2 ) β 0 e s e v 3 α ( u α ( v 1 α x + v 2 α y ) ) d s = I ^ ( u v 1 v 2 ) β E α , β ( u α ( v 1 α x + v 2 α y ) ) E α , β ( x E α , β y ) .
Here, in the second step we have introduced two additional formal integration variables u and v 3 and took advantage of the identity Equation (38) in order to suitably “reshape” the umbral image type expression without changing its evaluation result. This then permits use of the formal variable v 3 in order to realize the defining equation for the ML-f according to Equation (50), resulting in the ML-f at modified argument depending on the remaining formal integration variables as presented in the third step above. We thus conclude that the “Mittag-Leffler binomial” law should read
( x E α , β y ) n = I ^ ( u v 1 v 2 ) β ( u α ( v 1 α x + v 2 α y ) n = r = 0 n n r I ^ u α n + β v 1 α r + β x r v 2 α ( n r ) + β y ( n r ) = r = 0 n n r Γ ( n α + β ) x r y n r Γ ( α r + β ) Γ ( α ( n r ) + β ) .
Back to the fractional pseudo evolution problem described in Equation (49), note that the Mittag-Leffler function E α , β ( x ) may be utilized to construct a “pseudo eigenfunction” of the fractional time-derivative operator D ^ t = t μ as follows:
t μ D ^ t E μ , 1 ( M ^ t μ ) = M ^ E μ , 1 ( M ^ t μ ) + t μ Γ ( 1 μ ) .
Here, the last term in Equation (54) arises due to the action of the fractional derivative in the sense of Riemann-Liouville onto the constant term 1 of E α , β ( M ^ t μ ) , i.e., it is the contribution t μ 1 = t μ / Γ ( 1 μ ) (compare [2]). This permits us to determine the solution of the fractional pseudo-evolution Equation (49) in closed form as
Y ̲ ( t ) = E μ , 1 ( M ^ t μ ) Y ̲ 0 .
As a final example of an interesting fractional pseudo evolution problem, which in a sense combines the technique of the previous example with the one presented in the beginning of this section, consider
t μ F ( x , t ) = α O ^ x F ( x , t ) + t μ Γ ( 1 μ ) f ( x ) , F ( x , 0 ) = f ( x ) ( 0 < μ < 1 ) .
Following the previously presented strategy, we find an explicit solution of Equation (56) in the form
F ( x , t ) = E μ , 1 ( α t μ O ^ x ) f ( x ) .
Let us then specialize this result to the case of α = 1 , and with a differential operator O ^ x as in Equation (22),
O ^ x = α x + β 2 x 2 .
Taking advantage of the form of the Mittag-Leffler function E α , β ( x ) as presented in Equation (51),
E μ , 1 ( O ^ x t μ ) = I ^ v 0 e s e ( v t ) μ O ^ x d s ,
and by performing an analysis based on the Zassenhaus Formula (23) and the Crofton-Glaisher identity Equation (27), we obtain the general formula (for F ( x , 0 ) = f ( x ) ):
F ( x , t ) = E μ , 1 ( α t μ O ^ x ) f ( x ) = I ^ v 0 e s e 1 6 α 2 β ( v t ) 3 μ e α ( v t ) μ x ^ f x ^ + α β 2 ( v t ) μ β ( v t ) μ x 1
In particular, specializing further to the case F ( x , 0 ) = f ( x ) = 1 , the result of Equation (59) may be evaluated by using the third-order Hermite polynomial exponential generating function Formula (32), namely
F ( x , t ) = E μ , 1 ( α t μ O ^ x ) 1 = I ^ v 0 e s e α ( v t ) μ x 1 6 α 2 β ( v t ) 3 μ = e α 2 β 6 z 3 I ^ v 0 e s e ( v t ) μ z ( z = α x ) = e α 2 β 6 z 3 E μ , 1 ( z t μ ) .
Here, the second step follows by invoking the generating function formula given in Equation (32), while the third step is a consequence of Equation (51). Recalling both the definition of the third-order Hermite polynomials H n ( 3 ) ( x , y ) as given in Equation (31) and of the Mittag-Leffler function as given in Equation (50), we may reformulate the above result in the more explicit form
F ( x , t ) = r 0 t μ r Γ ( μ r + 1 ) H r ( 3 ) α x , α 2 β 6 .

4. Time-Ordering and Concluding Comments

In this section we touch upon on the problem of time-dependent pseudo-evolutive equations. We will then take a cursory look at the possibility of extending our techniques as introduced thus far to time-ordering problems, which as will become apparent is a nontrivial challenge even for the standard Cauchy problems.
To illustrate the difficulties we are going to meet, we first consider the non-homogeneous first order Laguerre differential equation
l t Y ( t ) = f ( t ) , Y ( 0 ) = Y 0 ,
where f ( t ) is a time-dependent function. The formal solution of this equation reads
Y ( t ) = l t 1 f ( t ) + Y 0 = t 1 t 1 t 1 f ( t ) + Y 0 ,
with the second identity obtained using Equation (48). We may then rewrite Equation (63) into the form
Y ( t ) = 0 t d t 1 t 1 0 t 1 f ( t 2 ) d t 2 + Y 0 .
The Laguerre integration can be carried out straightforwardly. For instance, if the integrand function is expandable as a (summable) series f ( t ) = n 0 t n n ! a n , we obtain
Y ( t ) = n 0 a n n ! ( n + 1 ) 2 t n + 1 + Y 0 ,
which is valid provided that we can exchange summation and integral sign, and given suitable convergence properties.
The problem becomes more difficult if we consider the equation
l t Y ( t ) = f ( t ) Y ( t ) , Y ( 0 ) = Y 0 ,
whose solution is obtained through a judicious application of the recipes we have discussed in the previous sections. We may indeed use the corresponding equation for the ordinary derivative, namely
y ( t ) = y 0 e 0 t f ( t ) d t ,
by replacing the exponential by its Laguerre counterpart, and the relevant argument by a suitable integration of the function f ( t ) . To better illustrate this technique, let us proceed by transforming Equation (66) into an integral equation,
Y ( t ) = l t 1 f ( t ) Y ( t ) = 0 t d t 1 t 1 0 t 1 f ( t 2 ) Y ( t 2 ) d t 2 + Y 0 .
We then eventually apply a Volterra-Neumann expansion, defined as
Y ( t ) = n 0 Y n ( t ) Y 0 ( t ) = Y 0 Y n + 1 ( t ) = 0 t d t 1 t 1 0 t 1 f ( t 2 ) Y n ( t 2 ) d t 2 .
It is worth noting that the inclusion of a non-homogeneous term l t Y ( t ) = f ( t ) Y ( t ) + g ( t ) does not introduce any further conceptual complication, but leads to additional inessential computational details, whence we omit the discussion of this more general case for brevity.
We may verify the correctness of the procedure by considering the example with f ( t ) = t in Equation (66), resulting in the fractional differential equation
1 t t t t Y ( t ) = Y ( t ) , Y ( 0 ) = Y 0 .
Noting that
1 t t t t = t 2 2 t 2 2 t 2 2 ,
we may conclude that the solution of Equation (70) with the initial condition Y ( 0 ) = 1 is just the Bessel function
J 0 ( t ) = l e t 2 2 .
We may then verify that the same result may be obtained by summing the series in Equation (69) directly, noting that
0 t d t 1 t 1 0 t 1 t 2 d t 2 = t 2 2 ,
which then entails that performing the Volterra-Neumann expansion Equation (69) indeed evaluates to
Y ( t ) = I ^ v e v t 2 2 = l e t 2 2 .
An entirely analogous computation permits to derive the solution of Equation (66) with f ( t ) = t m and with initial condition Y ( 0 ) = Y 0 = 1 , which reads
Y ( t ) = l e t m + 1 ( m + 1 ) 2 .
Unfortunately, the above procedures become considerably more complicated if f ( t ) is not just a monic function such as f ( t ) = t m as above. This is in fact a direct consequence of the lack of the semi-group property of the Laguerre exponential (see Equation (34)). If for instance f ( t ) = cos ( t ) , the solution of Equation (66) with initial condition Y ( 0 ) = Y 0 = 1 indeed becomes rather intricate:
Y ( t ) = n 0 Y n ( t ) Y n ( t ) = r = 0 ( 1 ) r t 2 r + n ( 2 r + n ) 2 n a r ( n 1 ) , Y 0 ( t ) = Y 0 = 1 n a r = k = 0 r n 1 a k ( 2 k + n 1 ) 2 [ 2 ( r k ) ] ! ( n 2 ) , 1 a r = 1 ( 2 r ) ! .
Our formalism remains applicable when passing to the setting of fractional derivatives. Consider for illustration the fractional evolution equation
t α Y ( t ) = f ( t ) Y ( t ) + Y 0 t α Γ ( 1 α ) , Y ( 0 ) = Y 0 ,
which may be transformed into integral form via noting that according to the definition of the fractional derivative in the sense of Riemann-Liouville, one finds that t α 1 = t α / Γ ( 1 α ) , and thus
Y ( t ) = t α f ( t ) Y ( t ) + Y 0 .
The use of the Riemann-Liouville integral in order to evaluate the action of t α yields [20]
Y ( t ) = 1 Γ ( α ) 0 t f ( τ ) Y ( τ ) ( t τ ) α 1 d τ + Y 0 ,
and the coefficients Y n ( t ) of the Volterra-Neumann expansion Equation (69) consequently satisfy the following recursion equation (with Y 0 ( τ ) = Y 0 as before):
Y n + 1 ( t ) = 1 Γ ( α ) 0 t f ( τ ) Y n ( τ ) ( t τ ) α 1 d τ ( n 0 ) .
Specializing for illustration to the case of f ( t ) = t and Y 0 = 1 , the expansion terms evaluate to
Y n ( t ) = t α + 1 Γ ( α ) n k = 0 n 1 B k ( α + 1 ) + 2 , α ,
where B ( x , y ) denotes the Euler Beta function.
The notions we have developed so far are in fact a necessary prerequisite for the development of the concepts associated with time-ordering. Let us thus pass to an illustrative first problem requiring explicit time-ordering, in the form of the matrix equation
t α Y ̲ ( t ) = M ^ ( t ) Y ̲ ( t ) + t α Γ ( 1 α ) Y ̲ 0 , Y ̲ ( 0 ) = Y ̲ 0 .
Here, Y ̲ ( t ) and Y ̲ 0 denote n-column vectors, while M ^ ( t ) denotes a non-singular time-dependent n × n matrix, which is in general assumed to be non-commuting with itself at different times (i.e., [ M ^ ( t ) , M ^ ( t ) ] 0 for t t ). It is evident that also in this case the most appropriate treatment is a series expansion, but ordering criteria such as those inherent in the well-known Dyson expansion are necessary. We will therefore write the formal solution of Equation (82) as
Y ̲ ( t ) = 1 Γ ( α ) 0 t M ^ ( τ ) Y ̲ ( τ ) ( t τ ) α 1 d τ + Y ̲ 0 .
The corresponding Volterra-Neumann series reads
Y ̲ ( t ) = U ^ ( t ) Y ̲ 0 ,
with a matrix-valued evolution operator U ^ ( t ) defined as
U ^ ( t ) = 1 ^ + 1 Γ ( α ) 0 t M ^ ( t 1 ) ( t t 1 ) α 1 d t 2 + 1 Γ ( α ) 2 0 t M ^ ( t 2 ) ( t t 2 ) α 1 0 t 2 M ^ ( t 1 ) ( t t 1 ) α 1 d t 1 d t 2 + 1 Γ ( α ) 3 0 t M ^ ( t 3 ) ( t t 3 ) α 1 0 t 3 M ^ ( t 2 ) ( t t 2 ) α 1 0 t 2 M ^ ( t 1 ) ( t t 1 ) α 1 d t 1 d t 2 d t 3 +
The series in Equation (85) has been obtained by translating to the fractional integration setting the usual expression given e.g., in [21]. The derivation of the corresponding Feynman-Dyson series along with the associated diagrammatic interpretation will be discussed in a forthcoming investigation.
In this paper, we have demonstrated that the operator and time-ordering techniques familiar from the setting of ordinary differential calculus may be exploited for analyzing fractional and for Laguerre operators as well. The methods we have highlighted are based on a suitable interpretation of operators and functions in umbral image form. The price to be paid is the demand for a certain level of abstraction allowing the search of a common thread yielding the pathway to generalized ordered formulae of Feynman-Dyson type.

Author Contributions

Conceptualization, G.D.; methodology, N.B., G.D.; validation, N.B., G.D.; formal analysis, N.B., G.D., A.L.; writing—original draft preparation, N.B., G.D., A.L.; writing— review and editing, N.B., A.L.

Funding

The work of N.B. was supported by a H2020 Marie Skłodowska-Curie Actions Individual Fellowship (grant # 753750—RaSiR). A.L. was supported by the NCN research project OPUS 12 no. UMO-2016/23/B/ST3/01714.

Acknowledgments

N.B. would like to thank the LPTMC (Paris 06) and ENEA Frascati for warm hospitality.

Conflicts of Interest

The authors declare no conflicts of interest. The founding sponsors had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, and in the decision to publish the results.

Appendix A. The Umbral Image Type Technique

For the readers’ convenience, we briefly recall the central definition of the umbral image type technique as introduced in [16]:
Definition A1
(Definition 2 of [16]). Let A = { λ } U V X be an alphabet of formal variables, where ⊎ denotes the operation of disjoint union, and where { λ } , U , V and X are four (disjoint) alphabets of auxiliary formal variables. We will typically employ notations such as X = { x , y , x 1 , x 2 , } , where we make use of the indexed variable notations in case of many variables for convenience. Let furthermore A = A { λ } .
We define a formal integration operator I ^ via specifying first its domain d o m ( I ^ ) as
d o m ( I ^ ) : = { S C [ G C ( A ) ] [ [ λ ] ] f o r   a l l A α s u p p ( S ) : r a n g e ( α | U ) C Z 0 } ,
whence elements of d o m ( I ^ ) are formal power series in λ with coefficients that are generalized polynomials over the alphabet A (where the extension to formal power series requires a suitable notion of summability, see below). Here, the notation r a n g e ( α | U ) C Z 0 entails that functions α : A C are required to not take zero or negative integer values when evaluated on elements of U . Then for some monomial A α d o m ( I ^ ) , which reads more explicitly (recall that by definition of C ( A ) , we have that α ( u i ) 0 and α ( v j ) 0 for only finitely many indices i and j)
A α = u 1 α ( u 1 ) u 2 α ( u 2 ) v 1 α ( v 1 ) v 2 α ( v 2 ) ,
the action of I ^ on A α is defined as
I ^ A α : = u i U ( α ) 0 d u i e u i u i v j V ( α ) 1 2 π i γ d v j e v j A α ˜ U ( α ) : = s u p p ( α ) U , V ( α ) : = s u p p ( α ) V α ˜ ( a ) : = α ( a ) if a A V α ( a ) if a V .
We extend I ^ by linearity to finite sums. For infinite sums, this requires an appropriate notion of convergence. A series i I c i A α i will be in d o m ( I ^ ) (the domain of I ^ ) if the family c i I ^ ( A α i ) i I is summable in the target (in the sense of discrete summability or compact convergence for entire functions).
Intuitively, the above definition allows to define a type of transformation of formal power series in which monomials in the auxiliary variable alphabet U = { u , u , u 1 , u 2 , } yield occurrences of Gamma functions in the numerators (i.e., I ^ u α = Γ ( α ) for α C { 0 , 1 , 2 , } ), while monomials in the formal variables of the alphabet V = { v , v , v 1 , v 2 , } yield occurrences of reciprocal Gamma functions (i.e., I ^ v β = 1 / Γ ( β ) for β C ). The only “rule” in translating an expressions involving Gamma and reciprocal Gamma functions into the “umbral image form” via using the formal integration operator I ^ is that if one wishes to express the product of two expressions via the operator, we have that
I ^ S 1 I ^ S 2 = I ^ S 1 S 2
if and only if the support of the formal power series S 1 and S 2 is disjoint (i.e., if S 1 and S 2 do not share any of the formal variables; cf. Lemma 1 of [16]). It is this “rule” that replaces a lot of more heuristic concepts in the traditional umbral calculus literature. To provide a quick application example, we present the following short list of expressions for illustration (where the exponents of the formal variables of the alphabet U are constrained to only take complex values not equal to 0 or a negative integer):
I ^ u α = Γ ( α ) , I ^ v β = 1 Γ ( β ) I ^ u α + n v α = Γ ( α + n ) Γ ( α ) = ( α ) n , I ^ u β v β + n = Γ ( β ) Γ ( β + n ) = 1 ( β ) n
A more complex set of examples is provided by the following expression quoted from [16] for the generalized hypergeometric functions, illustrating further the utility of the umbral image type formalism:
p F q α β ; z p F q α i 1 i p β j 1 j q ; z : = n 0 z n n ! ( α 1 ) n ( α p ) n ( β 1 ) n ( β q ) n = I ^ i = 1 p ( u i v i ) α i j = 1 q ( u j + p v j + p ) β j e z u 1 u p v p + 1 v p + q .

Appendix B. An Alternative to the Zassenhaus Formula

In Section 2, we have exploited the Zassenhaus formula to derive the solution of Laguerre or fractional Schrödinger equations, in which the operators appearing in the arguments of the respective pseudo evolution operators may be reduced to exponentials of the type
e X ^ + Y ^ = e X ^ e Y ^ e 1 2 [ X ^ , Y ^ ] e 1 3 [ Y ^ , [ X ^ , Y ^ ] ] + 1 6 [ X ^ , [ X ^ , Y ^ ] ] .
An alternative to the use of Equation (A6) is provided by the following Berry-type rule [22]
e X ^ + Y ^ = e m 2 12 m 2 X ^ 1 2 + X ^ e Y ^ ,
which is valid if
[ X ^ , Y ^ ] = m X ^ 1 2 .
To verify the correctness of Equation (A7), let us consider the following realization of the operators X ^ and Y ^ :
X ^ = α x 2 and Y ^ = β x ^
By applying the identity Equation (A7), noting that here indeed [ X ^ , Y ^ ] = m X ^ 1 2 (with m = 2 α β ), we obtain
e X ^ + Y ^ = e m 2 12 m 2 X ^ 1 2 + X ^ e Y ^ = e 1 3 α β 2 α β x + α x 2 e β x .
This is thus consistent with the result calculated via the Zassenhaus-type Formula (A6). A detailed proof of the Berry-type identity Equation (A7) independent of the Zassenhaus expansion may also be found in [18,23].

References

  1. Babusci, D.; Dattoli, G. Umbral methods and operator ordering. arXiv, 2011; arXiv:1112.1570. [Google Scholar]
  2. Dattoli, G.; Gorska, K.; Horzela, A.; Licciardi, S.; Pidatella, R.M. Comments on the properties of Mittag-Leffler function. Eur. Phys. J. Spec. Top. 2017, 226, 3427–3443. [Google Scholar] [CrossRef]
  3. Dattoli, G. Hermite-Bessel and Laguerre-Bessel functions: A by-product of the monomiality principle. In Proceedings of the Workshop on Advanced Special Functions and Applications; Aracne: Rome, Italy, 1999. [Google Scholar]
  4. Weyl, H. The Theory of Groups and Quantum Mechanics; Courier Corporation: North Chelmsford, MA, USA, 1950. [Google Scholar]
  5. Wilcox, R.M. Exponential operators and parameter differentiation in quantum physics. J. Math. Phys. 1967, 8, 962–982. [Google Scholar] [CrossRef]
  6. Dattoli, G.; Gallardo, J.C.; Torre, A. An algebraic view to the operatorial ordering and its applications to optics. Riv. Nuovo Cimento Ser. 3 1988, 11, 1–79. [Google Scholar] [CrossRef]
  7. Casas, F.; Murua, A.; Nadinic, M. Efficient computation of the Zassenhaus formula. Comput. Phys. Commun. 2012, 183, 2386–2391. [Google Scholar] [CrossRef]
  8. Gill, T.L. The Feynman-Dyson view. J. Phys. Conf. Ser. 2017, 845, 012023. [Google Scholar] [CrossRef]
  9. Ditkin, V.A.; Prudnikov, A.P. Operational calculus. Itogi Nauk. Tekhnik. Ser. Mat. Anal. 1967, 4, 7–82. [Google Scholar]
  10. Penson, K.A.; Blasiak, P.; Horzela, A.; Duchamp, G.H.E.; Solomon, A.I. Laguerre-type derivatives: Dobiński relations and combinatorial identities. J. Math. Phys. 2009, 50, 083512. [Google Scholar] [CrossRef]
  11. Dattoli, G. Generalized polynomials, operational identities and their applications. J. Comput. Appl. Math. 2000, 118, 111–123. [Google Scholar] [CrossRef]
  12. Dattoli, G.; Srivastava, H.M.; Zhukovsky, K. A new family of integral transforms and their applications. Integr. Trans. Spec. Funct. 2006, 17, 31–37. [Google Scholar] [CrossRef]
  13. Dattoli, G.; He, M.; Ricci, P. Eigenfunctions of laguerre-type operators and generalized evolution problems. Math. Comput. Modell. 2005, 42, 1263–1268. [Google Scholar] [CrossRef]
  14. Dattoli, G.; Licciardi, S.; Pidatella, R. Theory of generalized trigonometric functions: From Laguerre to Airy forms. J. Math. Anal. Appl. 2018, 468, 103–115. [Google Scholar] [CrossRef]
  15. Babusci, D.; Dattoli, G.; Górska, K.; Penson, K. The spherical Bessel and Struve functions and operational methods. Appl. Math. Comput. 2014, 238, 1–6. [Google Scholar] [CrossRef]
  16. Behr, N.; Dattoli, G.; Duchamp, G.H.; Penson, S. Operational methods in the study of Sobolev-Jacobi polynomials. Mathematics 2019, 7, 124. [Google Scholar] [CrossRef]
  17. Crofton, M.W. Theorems in the calculus of operations. Q. J. Math. 1879, 16, 323–352. [Google Scholar]
  18. Dattoli, G.; Ottaviani, P.L.; Torre, A.; Vázquez, L. Evolution operator equations: Integration with algebraic and finitedifference methods. Applications to physical problems in classical and quantum mechanics and quantum field theory. Riv. Nuovo Cimento 1997, 20, 3–133. [Google Scholar] [CrossRef]
  19. Mittag-Leffler, G. Sur la nouvelle fonction Eα (x). CR Acad. Sci. Paris 1903, 137, 554–558. [Google Scholar]
  20. Diethelm, K. The Analysis of Fractional Differential Equations; Springer: Berlin/Heidelberg, Germany, 2010. [Google Scholar] [CrossRef]
  21. Louisell, W.H. Quantum Statistical Properties of Radiation; Wiley: New York, NY, USA, 1973; Volume 7. [Google Scholar]
  22. Berry, M.V. The Diffraction of Light by Ultrasound; Academic Press: Cambridge, MA, USA, 1966. [Google Scholar]
  23. Babusci, D.; Dattoli, G.; Del Franco, M. Lectures on Mathematical Methods for Physics; Thecnical Report; ENEA: Rome, Italy, 2010. [Google Scholar]
Figure 1. Polar plot of the fundamental l-trigonometric relation.
Figure 1. Polar plot of the fundamental l-trigonometric relation.
Axioms 08 00035 g001
Figure 2. Plot of the l s ( x ) function.
Figure 2. Plot of the l s ( x ) function.
Axioms 08 00035 g002

Share and Cite

MDPI and ACS Style

Behr, N.; Dattoli, G.; Lattanzi, A. Operator Ordering and Solution of Pseudo-Evolutionary Equations. Axioms 2019, 8, 35. https://doi.org/10.3390/axioms8010035

AMA Style

Behr N, Dattoli G, Lattanzi A. Operator Ordering and Solution of Pseudo-Evolutionary Equations. Axioms. 2019; 8(1):35. https://doi.org/10.3390/axioms8010035

Chicago/Turabian Style

Behr, Nicolas, Giuseppe Dattoli, and Ambra Lattanzi. 2019. "Operator Ordering and Solution of Pseudo-Evolutionary Equations" Axioms 8, no. 1: 35. https://doi.org/10.3390/axioms8010035

APA Style

Behr, N., Dattoli, G., & Lattanzi, A. (2019). Operator Ordering and Solution of Pseudo-Evolutionary Equations. Axioms, 8(1), 35. https://doi.org/10.3390/axioms8010035

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop