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,
Here, is an operator playing the role of the derivative with respect to the variable t, is a constant, is an operator of differential or pseudo-differential nature (to be specified later in our discussion), and where 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
operator, such that
Constructing a
pseudo evolution operator (PEO) as
which leads to the solution of Equation (
1) via
Establishing rules that permit the explicit evaluation of the action of the PEO
on the initial function
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
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
This is an additional source of difficulties also when dealing with operator-ordering in pseudo-evolutionary problems. If
are non-commuting operators, such that
it is well known that even for the exponential function one finds that [
4]
Instead, one finds the following types of expansions [
5,
6]
where
are expressed in terms of chains of commutators involving
, as it happens e.g., for the Zassenhaus expansion [
7]. It is evident that in the special case
, the
and
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
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
where
is the
Laguerre derivative [
3,
10,
11,
13].
The eigenfunction of the Laguerre derivative operator is the Bessel-like function
[
14],
which satisfies
.
In view of explicit computations, it will prove advantageous to express
via an
umbral image [
15,
16] (where we refer to
Appendix A for the explicit definition of the full formalism)
Here,
v is a formal variable, and
a formal integration operator, which acts according to
We can therefore write the solution of our problem as
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
the operator of multiplication by
x, i.e.,
,
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
). The solution of our problem can be cast in a first step into the form
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
Applying the Weyl disentanglement rule (taking advantage of the fact that
), we find that
We then eventually arrive at the closed-form expression
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
and
). Assuming for simplicity
, we find that
Finally, after Taylor-expanding the exponential and evaluating the action of the formal integral operator
according to Equation (
12), we find the closed-form expression
where
is the Bessel-like function defined as
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
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
Introducing the auxiliary operators
we find the commutation relations
with all higher nested commutators vanishing. Thus invoking the “left-oriented” form of the formula (setting
,
and
) and using the commutation relations Equation (
25), we obtain
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
and for an integer-valued parameter
,
Note that we provide this identity in “operational form”, i.e., by using the formal multiplication operator
, 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
, resulting in
Here and throughout this paper, in expressions such as
, the occurrence of the symbol “1” entails that the expression is to be evaluated by expanding
into normal-ordered form (i.e., into a series in the normal-ordered monomials
for
), followed by acting on 1 (which due to
for
in effect amounts to dropping all terms of the expansion involving non-zero powers of
).
Combining this result with the explicit formula for the eigenfunctions
as provided in Equation (
11), we thus finally arrive at the explicit solution for
evolving according to the pseudo-evolution Equation (
22) with initial condition
:
Specializing for simplicity to the case
, Equation (
29) evaluates to
Coincidentally, the expression obtained in the last step has an interesting formal meaning: consider the
third order Hermite polynomials , which are defined as
and whose exponential generating function (EGF) reads
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
in the alternative form
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
(for
commuting variables) of the ordinary exponential function, whence
but rather satisfies
Here, the symbol
denotes the
composition rule
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
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)
we may compute a
“Laguerre Newton binomial” type result as follows:
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
with
the Laguerre-type time-derivative, and where
and
denote an
matrix and
n-element column vectors, respectively. We specialize our discussion to the case of a non-singular
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
By application of the Cayley-Hamilton theorem, we thus obtain
where
denotes the
unit matrix. Let us now consider for illustration a matrix with zero diagonal entries and imaginary eigenvalues, namely
According to Equation (
43), the corresponding PEO
can be reduced to the “pseudo rotation matrix”
where
and
denote the
Laguerre (co-)sine functions [
14] defined as
The relevant geometrical meaning differs from that of the ordinary circular functions and is illustrated in
Figure 1, where we have plotted
against
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
We also note that the Laguerre derivative satisfies the identity
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
Following the paradigm of the PEO method, the first step in solving Equation (
49) consists in finding an eigenfunction of the fractional differential operator
. To this end, we recall the definition of the so-called
Mittag-Leffler function (ML-f) [
19]
(for
),
Here, we have yet again taken advantage of the formal integral operator
as introduced in [
16] (compare Equation (
12)). An alternative useful expression for the ML-f
may be obtained via taking a suitable Laplace transform of Equation (
50), whence
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:
Here, in the second step we have introduced two additional formal integration variables
u and
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
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
Back to the fractional pseudo evolution problem described in Equation (
49), note that the Mittag-Leffler function
may be utilized to construct a “pseudo eigenfunction” of the fractional time-derivative operator
as follows:
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
, i.e., it is the contribution
(compare [
2]). This permits us to determine the solution of the fractional pseudo-evolution Equation (
49) in closed form as
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
Following the previously presented strategy, we find an explicit solution of Equation (
56) in the form
Let us then specialize this result to the case of
, and with a differential operator
as in Equation (
22),
Taking advantage of the form of the Mittag-Leffler function
as presented in Equation (
51),
and by performing an analysis based on the Zassenhaus Formula (
23) and the Crofton-Glaisher identity Equation (
27), we obtain the general formula (for
):
In particular, specializing further to the case
, the result of Equation (
59) may be evaluated by using the third-order Hermite polynomial exponential generating function Formula (
32), namely
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
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
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
where
is a time-dependent function. The formal solution of this equation reads
with the second identity obtained using Equation (
48). We may then rewrite Equation (
63) into the form
The Laguerre integration can be carried out straightforwardly. For instance, if the integrand function is expandable as a (summable) series
, we obtain
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
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
by replacing the exponential by its Laguerre counterpart, and the relevant argument by a suitable integration of the function
. To better illustrate this technique, let us proceed by transforming Equation (
66) into an integral equation,
We then eventually apply a
Volterra-Neumann expansion, defined as
It is worth noting that the inclusion of a non-homogeneous term 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
in Equation (
66), resulting in the fractional differential equation
Noting that
we may conclude that the solution of Equation (
70) with the initial condition
is just the Bessel function
We may then verify that the same result may be obtained by summing the series in Equation (
69) directly, noting that
which then entails that performing the Volterra-Neumann expansion Equation (
69) indeed evaluates to
An entirely analogous computation permits to derive the solution of Equation (
66) with
and with initial condition
, which reads
Unfortunately, the above procedures become considerably more complicated if
is not just a monic function such as
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
, the solution of Equation (
66) with initial condition
indeed becomes rather intricate:
Our formalism remains applicable when passing to the setting of fractional derivatives. Consider for illustration the fractional evolution equation
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
, and thus
The use of the Riemann-Liouville integral in order to evaluate the action of
yields [
20]
and the coefficients
of the Volterra-Neumann expansion Equation (
69) consequently satisfy the following recursion equation (with
as before):
Specializing for illustration to the case of
and
, the expansion terms evaluate to
where
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
Here,
and
denote
n-column vectors, while
denotes a non-singular time-dependent
matrix, which is in general assumed to be non-commuting with itself at different times (i.e.,
for
). 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
The corresponding Volterra-Neumann series reads
with a matrix-valued evolution operator
defined as
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.