Hindawi Publishing Corporation ξ e Scientiο¬c World Journal Volume 2014, Article ID 454865, 8 pages http://dx.doi.org/10.1155/2014/454865
Research Article Solution of Some Types of Differential Equations: Operational Calculus and Inverse Differential Operators K. Zhukovsky Faculty of Physics, M. V. Lomonosov Moscow State University, Leninskie Gory, Moscow 119899, Russia Correspondence should be addressed to K. Zhukovsky;
[email protected] Received 23 January 2014; Accepted 23 February 2014; Published 30 April 2014 Academic Editors: D. Baleanu and H. Jafari Copyright Β© 2014 K. Zhukovsky. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. We present a general method of operational nature to analyze and obtain solutions for a variety of equations of mathematical physics and related mathematical problems. We construct inverse differential operators and produce operational identities, involving inverse derivatives and families of generalised orthogonal polynomials, such as Hermite and Laguerre polynomial families. We develop the methodology of inverse and exponential operators, employing them for the study of partial differential equations. Advantages of the operational technique, combined with the use of integral transforms, generating functions with exponentials and their integrals, for solving a wide class of partial derivative equations, related to heat, wave, and transport problems, are demonstrated.
1. Introduction Most of physical systems can be described by appropriate sets of differential equations, which are well suited as models for systems. Hence, understanding of differential equations and finding its solutions are of primary importance for pure mathematics as for physics. With rapidly developing computer methods for the solutions of equations, the question of understanding of the obtained solutions and their application to real physical situations remains opened for analytical study. There are few types of differential equations, allowing explicit and straightforward analytical solutions. It is common knowledge that expansion into series of Hermite, Laguerre, and other relevant polynomials [1] is useful when solving many physical problems (see, e.g., [2, 3]). Generalised forms of these polynomials exist with many variables and indices [4, 5]. In what follows, we develop an analytical method to obtain solutions for various types of partial differential equations on the base of operational identities, employing expansions in series of Hermite, Laguerre polynomials, and their modified forms [1, 6]. The key for building these solutions will be an operational approach and development of the formalism of inverse functions and inverse differential operators, already touched in [7, 8]. We will demonstrate in what follows that when used properly and combined,
in particular, with integral transforms, such an approach leads to elegant analytical solutions with transparent physical meaning without particularly cumbersome calculations.
2. Inverse Derivative For a common differential operator π· = π/ππ₯ we can define an inverse derivative, such that upon the action on a function π(π₯) it gives another function πΉ(π₯): π·β1 π (π₯) = πΉ (π₯) ,
(1)
whose derivative is πΉσΈ (π₯) = π(π₯). Evidently, the inverse derivative π·β1 is executed by an integral operator being the inverse of differential operator, acting on π(π₯), and its general form is β« π(π₯) = πΉ(π₯) + πΆ, where πΆ is the constant of integration. The action of its πth order π·π₯βπ π (π₯) =
π₯ 1 β« (π₯ β π)πβ1 π (π) ππ, (π β 1)! 0
(π β π = {1, 2, 3, . . .})
(2)
2
The Scientific World Journal
can be complemented with the definition of the action of zero order derivative as follows: π·π₯0 π (π₯)
= π (π₯) ,
(3)
Laguerre polynomials of two variables [4] πΏ π (π₯, π¦) = exp (βπ¦
so that evidently π·π₯βπ 1 =
π₯π , π!
(π β π0 = π βͺ {0}) .
(4)
are related to the following operator ππ₯ π₯ππ₯ [10]: πΏ π·π₯
In what follows we will appeal to various modifications of the following equation: ]
(π½2 β π·2 ) πΉ (π₯) = π (π₯) .
(5)
Thus, it is important to construct the particular integral πΉ(π₯) with the help of the following operational identity (see, e.g., [6]): β]
πΜ
β 1 = ππ‘) π‘]β1 ππ‘, β« exp (βΜ Ξ (]) 0
π π π , π₯ = ππ₯ ππ₯ ππ·π₯β1
[ πΏ π·π₯ , π·π₯β1 ] = β1,
([π΄, π΅] = π΄π΅ β π΅π΄) .
ππ¦ πΏ π (π₯, π¦) = β (ππ₯ π₯ππ₯ ) πΏ π (π₯, π¦) ,
(π½ β π· ) π (π₯)
(β1)π π¦πβπ π₯π , 2 π=0 (π β π)!(π!)
(7)
2
(8)
Exponential operator mentioned above is closely related to Hermite orthogonal polynomials: [π/2]
π»π (π₯, π¦) = π! β
π=0
π₯πβ2π π¦π , (π β 2π)!π!
(9)
as demonstrated in [4, 5, 10] by operational relations: π»π (π₯, π¦) = exp (π¦
π2 ) π₯π . ππ₯2
(10)
Moreover, the following generating function for Hermite polynomials exists: β
π‘π π»π (π₯, π¦) . π=0 π!
exp (π₯π‘ + π¦π‘2 ) = β
(11)
πΏ π (π₯, 0) =
(βπ₯)π . π!
(18)
In the following sections, we will investigate the possibilities to solve some partial differential equations, involving the differential operators studied above. Now we just note how the technique of inverse operator, applied for derivatives of various orders and their combinations and combined with integral transforms, allows for easy and straightforward solutions of various types of differential equations.
3. Diffusion Type, Heat Propagation Type Problems, and Inverse Derivatives Heat propagation and diffusion type problems play a key role in the theory of partial differential equations. Combination of exponential operator technique and inverse derivative together with the operational identities of the previous section is useful for the solution of a broad spectrum of partial differential equations, related to heat and diffusion processes. Some of them have been already studied by operational method (see, e.g., [10, 12]). Below we will focus on the generalities of the solution of the following problem: π Μ {πΉ (π₯, π¦)} , πΉ (π₯, π¦) = (πΜ + π) ππ¦
(19)
with initial conditions
Note also an easy to prove and useful relation [11]: π§π π»π (π₯, π¦) = π»π (π₯π§, π¦π§2 ) .
(17)
with proper initial conditions:
We will also make explicit use of the generalized form of the Glaisher operational rule [9]; action of the operator πΜ = exp(π‘π·π₯2 ) on the function π(π₯) = exp(βπ₯2 ) yields
1 π₯2 = exp (β ). 1 + 4π¦ β1 + 4π¦
(16)
π
πΏ π (π₯, π¦) = π!β
Μ (π₯) = exp (π¦ π ) exp (βπ₯2 ) ππ ππ₯2
(15)
where
2 β]
β 1 β« exp (βπ½2 π‘) π‘]β1 exp (π‘π·2 ) π (π₯) ππ‘. Ξ (]) 0
(14)
They represent solutions of the following partial differential equation with proper initial conditions:
(6)
For the operator πΜ = π½2 β π·2 we have
=
=
sometimes called Laguerre derivative πΏ π·π₯ . Note their noncommutative relation with the inverse derivative operator:
min {Re (π) , Re (])} > 0.
2
π (β1)π π¦πβπ π₯π π π (βπ₯)π π₯ ) = π!β 2 ππ₯ ππ₯ π! π=0 (π β π)!(π!) (13)
(12)
πΉ (π₯, 0) = π (π₯) .
(20)
The Scientific World Journal
3
We will employ operational approach, combined with integral transforms and exponential operator technique. Formal solution for our generic formulation reads as follows: Μ + π)) Μ π (π₯) . πΉ (π₯, π¦) = exp (π¦ (π
(21)
Μ may not We would like to underline that operators πΜ and π commute, so, dependently on the value of their commutator, we will obtain different sequences of operators in (23), Μ [13, 14]. In the simplest case, when disentangling πΜ and π Μ and πΜ are multiplication and differentiation operators π operators, respectively, they can be easily disentangled in the exponent with account for Μ π] Μ = 1, [π,
(22)
which yields
π¦2 Μ exp (π¦π) Μ π (π₯) . ) exp (π¦π) 2
(23)
Explicit expressions for the action of these operators on the initial condition function π(π₯) can be obtained by integral transforms, series expansions, and operational technique to be evaluated in every special case. This formulation, simple in its essence, nevertheless has wide application and allows us to frame either some of integrodifferential equations in this scheme. An elegant and interesting example is given by the following equation: π₯ π π π πΉ (π₯, π‘) = β π₯ πΉ (π₯, π‘) + β« πΉ (π, π‘) ππ, ππ‘ ππ₯ ππ₯ 0
(24)
πΉ (π₯, 0) = π (π₯) .
π ) π (π₯) , ππ·π₯β1
Μ + π΅) Μ π (π₯) , πΉ (π₯, π‘) = exp (π΄
(25)
where operators π΅Μ = π‘π·π₯β1
(26)
π π π π (π₯, π‘) = β π₯ π (π₯, π‘) , ππ‘ ππ₯ ππ₯
π (π₯, 0) = π (π₯) .
(30)
Relevant studies were performed in [10, 12]. The initial condition function π(π₯, 0) = π(π₯) can be written as follows: (31)
and the image π(π₯) is explicitly given by the following integral: β
π (π₯) = β« exp (βπ) π (π₯π) ππ, 0
(32)
which is supposed to converge. Then the result of the Laguerre diffusion (29) appears in the form of the translation of the image function π: π (π₯, π¦) = π (π·π₯β1 β π‘) 1.
(33)
Consequently, we have to apply the exponential operator exp(π‘π·π₯β1 ), which can be expanded in series: (34)
The simplest example of the initial function π(π₯) = exp(βπ₯) demonstrates the technique sketched above, resulting in π (π₯) =
1 (1 + π₯)
when |π₯| < 1,
(35)
and the Laguerre diffusion contribution (29) produces the following function:
do not commute: Μ π΅] Μ = βπ‘2 . [π΄,
(29)
where π·π₯β1 is defined in (2). Equation (29) represents, in fact, the diffusion process and it is the solution of the following initial value problem:
π βπ π‘2 β π‘ π·π₯ πΉ (π₯, π‘) = exp ( ) β π (π₯, π‘) . 2 π=0 π!
Its formal solution reads
Μ = βπ‘πΏ π·π₯ = βπ‘ π , π΄ ππ·π₯β1
π (π₯, π‘) = exp (βπ‘
π (π₯) = π (π·π₯β1 ) 1,
Μ + π)) Μ π (π₯) πΉ (π₯, π¦) = exp (π¦ (π = exp (β
depend on the explicit form of this function. In the most general case of π(π₯) we may take advantage of the inverse derivative technique. First, consider the action of the operator exp(βπ/(ππ·π₯β1 )) on π(π₯):
(27)
Evidently, operators in the exponential disentangle: Μ + π΅) Μ π (π₯) πΉ (π₯, π‘) = exp (π΄ (28) π π‘2 = exp ( ) exp (π‘π·π₯β1 ) exp (βπ‘ β1 ) π (π₯) . 2 ππ·π₯ Thus, we have obtained the solution of the integrodifferential equation (24) as a sequence of exponential operators, transforming the initial condition π(π₯). Our further steps
π (π₯, π‘) =
1 1 exp (β π₯) . 1βπ‘ 1βπ‘
(36)
The inverse derivative action on the exponent reads π·π₯βπ exp(πΌπ₯) = exp(πΌπ₯)/πΌπ and eventually we obtain π βπ π‘2 1 β π‘ π·π₯ 1 πΉ (π₯, π‘) = exp ( ) exp (β π₯) β 2 1 β π‘ π=0 π! 1βπ‘
1 π₯ π‘2 = exp (β ) exp (βπ‘ (1 β π‘) + ) . 1βπ‘ 1βπ‘ 2
(37)
4
The Scientific World Journal
In conclusion of the present chapter we consider the example of the solution of a heat propagation type equation by operational method, involving the inverse derivative operator and exponential operator technique. We recall that the common heat equation with initial condition problem π π (π₯, π‘) = ππ₯2 π (π₯, π‘) , ππ‘
π (π₯, 0) = π (π₯)
exp (πΜ2 ) =
exp (ππ₯ππ₯ ) π (π₯) = π (ππ π₯) ,
(39)
π΄ (π₯, π‘) =
π (π₯, 0) = π (π₯) .
(40)
The Laguerre heat-type propagation problem (40) possesses the following solution (see the Appendix): β 2 βππ₯ π₯2 1 , β ) π (π) ππ. β« πβπ /4π‘ π»πΆ0 ( 2π‘ 4π‘ 2βππ‘ ββ (41)
4. Operational Approach and Other Types of Differential Equations Operational approach to solution of partial differential equations, demonstrated on the examples of diffusion-like and heat-like equations with ππ₯ π₯ππ₯ derivatives, can be further extended to other equation types. Consider the following example of a rather complicated differential equation: π2 π 1 π π΄ (π₯, π‘) = π₯2 2 π΄ (π₯, π‘) + ππ₯ π΄ (π₯, π‘) π ππ‘ ππ₯ ππ₯ β ππ΄ (π₯, π‘) ,
(47)
(42)
π (π₯) = π΄ (π₯, 0) ,
where π, π, and π are some arbitrary constant coefficients and function π(π₯) = π΄(π₯, π‘ = 0) is the initial condition. By introducing operator π· = π₯ππ₯ and distinguishing the perfect square, we rewrite (42): 1 π π 2 π΄ (π₯, π‘) = ((π· + ) β π) π΄ (π₯, π‘) , ππ π‘ 2
(43)
where
exp (βπππ‘) βπ Γβ«
In complete analogy with the above statements, we can solve the heat-type equation with differential operator πΏ π·π₯ (14) with, for example, the following initial condition:
π (π₯, π‘) =
(46)
we obtain the following compact expression for π΄(π₯, π‘):
2 β 1 (π₯ β π) ) π (π) ππ. β« exp (β 4π‘ 2βππ‘ ββ
π π (π₯, π‘) = πΏ π·π₯2 π (π₯, π‘) , ππ‘
1 β Μ ππ β« exp (βπ2 + 2ππ) βπ ββ
and applying exp(ππ·) according to
(38)
can be solved by Gauss transforms: π (π₯, π‘) =
Now making use of the operational identity
β
ββ
exp [βπ2 + ππΌ
π ] π (π₯ exp (ππΌ)) ππ, 2π (48)
where πΌ = πΌ (π‘) = 2βππ‘.
(49)
Consider two simple examples of initial condition functions. The first one is π (π₯) = π₯π .
(50)
Then we immediately obtain the solution of (42) as follows: π΄ (π₯, π‘) = π₯π exp {ππ‘ (π2 + ππ β π)} .
(51)
The second example is given by the following initial condition function: π (π₯) = ln π₯.
(52)
Trivial computations yield the following solution: π΄ (π₯, π‘) = (ln π₯ + ππ‘π) exp (βππ‘π) .
(53)
Thus, operational technique, combined with integral transforms, operational identities, and extended forms of orthogonal polynomials, represents powerful tool for finding solutions of various classes of differential equations and initial value problems. Note that within the framework of inverse differential operators, developed and described above, the usage of the evolution operator method opens new possibilities, which we will elucidate in what follows. Let us consider the following generalization of the heat equation: ππ‘ πΉ (π₯, π‘) = πΌππ₯2 πΉ (π₯, π‘) + π½π₯πΉ (π₯, π‘) ,
(54)
with the initial condition: 2
π π=π+( ) . 2
(44)
Thus, the following exponential solution for (42) appears: π΄ (π₯, π‘) = exp {ππ‘ ((π· +
π 2 ) β π)} π (π₯) . 2
(45)
πΉ (π₯, 0) = π (π₯) .
(55)
The evolution type equation (54) contains linear coordinate term in addition to the second order derivative. Its formal Μ solution can be written via the evolution operator π: Μ (π₯) , πΉ (π₯, π‘) = ππ
(56)
The Scientific World Journal
5
where
and we write Μ = exp (π΄ Μ + π΅) Μ , π
Μ = πΌπ‘π2 , π΄ π₯
π΅Μ = π½π‘π₯.
Μ in (56) is the The exponential of the evolution operator π sum of two noncommuting operators and it can be written as the ordered product of two exponential operators. Indeed, Μ and π΅Μ has the following nonzero value: the commutator of π΄ Μ π΅] Μ = ππ΄ Μ1/2 = 2πΌπ½π‘2 ππ₯ , [π΄,
π = 2π½π‘3/2 πΌ1/2 .
(58)
Then we can apply the disentanglement operational identity: Μ Μ
Μ1/2 +π΄ Μ π΅Μ
2
ππ΄+π΅ = π(π /12)β(π/2)π΄
π ,
(59)
π
ππ2 ππ₯ 2ππππ₯ πππ₯2
π π (π₯) = π
π π
π
π (π₯) ,
(60)
where π, π are constant parameters. With their help, we obtain Μ for (54): the evolution operator π Μ = πΞ¦(π₯,π‘;π½) Ξ Μ π. Μ π
(61)
Note that we have factorized two commuting operators: the operator of translation in space 2
Μ = ππΌπ½π‘ ππ₯ Ξ
(62)
Thus, (54) with initial condition (55) has the following explicit solution: πΉ (π₯, π‘) = πΞ¦(π₯,π‘;πΌ,π½)
β 2 2 1 πβ(π₯+πΌπ½π‘ βπ) /4π‘πΌ π (π) ππ, β« 2βππΌπ‘ ββ
2
1 Ξ¦ (π₯, π‘; πΌ, π½) = πΌπ½2 π‘3 + π½π‘π₯. 3
(64)
Μ on the initial condition function π(π₯) yields The action of π the following solution for our problem: Ξ¦(π₯,π‘;πΌ,π½) Μ Μ
πΉ (π₯, π‘) = π
Ξππ (π₯) .
2 2
πΉ(π₯, π‘)|π(π₯)=exp(βπ₯2 )
(π₯ + π½π‘ ) exp (Ξ¦ (π₯, π‘)) = exp (β ), β1 + 4π‘ 1 + 4π‘ π‘ β₯ 0. (71)
ππ ππ₯ π (π₯) = π (π₯ + π )
(66)
and that the action of πΜ on the function π(π₯) yields the solution of the ordinary heat equation, through GaussWeierstrass transform (A.1). Accordingly, we denote πΌπ‘ππ₯2
Μ (π₯) β‘ π π (π₯, π‘) β‘ ππ
π (π₯) ,
The above result is exactly the generalization of Gleisher rule (8), considered earlier in the context of the heat equation. Thus, (71) is the solution of the ordinary heat equation with π½ = 0, when the initial function is Gaussian. Another interesting example of solving (54) appears when the initial function π(π₯) allows the expansion in the following series: π (π₯) = β ππ π₯π ππΎπ₯ .
(65)
Thus, we conclude from the form of (65) that the problem (54) with the initial condition (55) can be solved by the Μ (62) and πΜ consequent application of commuting operators Ξ Ξ¦(π₯,π‘;πΌ,π½) (63) to π(π₯), apart from the factor π . Now the explicit Μ form of the solution (65) can be obtained by recalling that Ξ acts as a translation operator
(67)
(70)
Moreover, this is a general observation for this type of equation, valid for any function π(π₯) (provided the integral converges). In other words, we have obtained the solution of Fokker plank equation as a consequent action of operator πΜ Μ of translation on the initial of heat diffusion and operator Ξ condition function. Note that π(π₯, π‘) is the solution of the heat equation, representing a natural propagation phenomenon. Μ and the The effect, produced by the translation operator Ξ Μ is best illustrated with the example of Gaussian operator π, π(π₯) = exp(βπ₯2 ) evolution, when (69) becomes
(63)
Μ is which is, in fact, operator πΜ = exp(π‘π·π₯2 ). The phase in π written as follows:
(69)
provided that the integral converges. Summarizing the above outlined procedure, we conclude that a solution for (54) consists in finding a Gauss transformed function π with a shifted argument:
and operator πΜ = ππΌπ‘ππ₯ ,
(68)
πΉ (π₯, π‘) = πΞ¦(π₯,π‘;πΌ,π½) π (π₯ + πΌπ½π‘2 , π‘) .
and the following chain rule: πππ₯2 ππ₯
Μ (π₯, π‘) = π (π₯ + πΌπ½π‘2 , π‘) . Ξπ
(57)
(72)
π
In this case, we refer to the identity 2
exp (π¦π·π₯2 ) π₯π ππΎ π₯ = π(πΎ π₯+πΎ π¦) π»π (π₯ + 2πΎπ¦, π¦) ,
(73)
which arises from exp (π¦
ππ ππβ1 ) π = π (π₯ + ππ¦ )1 (π₯) ππ₯π ππ₯πβ1
(74)
with account for the generating function of Hermite polynomials (11). Then the action (67) of operator πΜ on the initial function (72) produces π (π₯, π‘) = β ππ ππ (π₯, π‘) . π
πΎ(π₯+πΎπ)
ππ (π₯, π‘) = π
π»π (π₯ + 2πΎπ, π) ,
(75) π = πΌπ‘.
6
The Scientific World Journal
Μ shifts the arguIn addition, the translation, operated by Ξ, ment: π(π₯, π‘) β π(π₯ + ππ, π‘), π = π½π‘. Thus, we obtain the solution in a form of Hermite polynomial as follows: πΉ (π₯, π‘) = β ππ πΞ¦+Ξ¦1 π»π (π₯ + π
πΌπΎ2 π½2 π½ (2π‘ + π‘2 2 ) , πΌπ‘) , π½ πΎ πΎ (76)
where Ξ¦ is defined by (64) and Ξ¦1 = πΎ (π₯ +
πΌπΎ2 π½ π½ 2 (π‘ + (π‘ ) )) . π½ πΎ πΎ
(77)
It is now evident that, for a short time, when π‘ βͺ πΎ/π½, the solution will be spanning in space off the initial function, modulated by π»π (π₯, πΌπ‘) and exponentially depending on time: πΉ(π₯, π‘)|π‘βͺπΎ/π½ β β ππ exp (πΌπΎ2 π‘ + πΎπ₯) π»π (π₯, πΌπ‘) , π
2
πΉ(π₯, π‘)|π‘βͺπΎ/π½,π‘βͺ1 β β ππ ππΌπΎ π‘ ππΎπ₯ π₯π .
(78)
Note that for extended times π‘ β« πΎ/π½ we have dominance of time variable: Ξ¦ β« Ξ¦1 , and the solution asymptotically behaves as exp(π½π‘π₯), while π₯ plays minor role in π»π (π₯+2πΎπΌπ‘+ πΌπ½π‘2 , πΌπ‘). The same operational technique as employed for the treatise of (54) can be easily adopted for the solution of the SchrΒ¨odinger equation: β2 2 π Ξ¨ (π₯, π‘) + πΉπ₯Ξ¨ (π₯, π‘) , 2π π₯
(79)
Ξ¨ (π₯, 0) = π (π₯) , where πΉ is constant (πΉ has the dimension of force). Indeed, rescaling variables in (79), we obtain the form of equation, similar to (54): πππ Ξ¨ (π₯, π) =
βππ₯2 Ξ¨ (π₯, π)
+ ππ₯Ξ¨ (π₯, π) ,
(80)
2πΉπ . β2
(81)
where π=
βπ‘ , 2π
π=
Following the operational methodology, developed for (54), we write the following solution of (80) in the form (56): Μ (π₯) , Ξ¨ (π₯, π) = ππ
(82)
where Μ = exp (πππ») Μ , π
Μ = π2 β ππ₯. π» π₯
Μ (π₯, π) = π (π₯ + ππ2 , π) , Ξπ
πΜ = exp (ππππ₯2 ) ,
Μ (π₯) = π (π₯, ππ) ππ
(86)
where Ξ¦ is defined by (64). The integral form of the solution then is written as follows: Ξ¨ (π₯, π‘) = exp (βπΞ¦ (π₯, π; π))
Γβ«
β
ββ
1 2βπππ (87)
2
exp (β
(π₯ + ππ2 β π) 4ππ
) π (π) ππ.
Again, as well as in (70), without any assumption on the nature of the initial condition function π(π₯) of the SchrΒ¨odinger equation, its solution (88)
where π(π₯, π‘) is given by (85) and is expressed in terms of the function of two variables, obtained by the consequent application of the heat propagation and translation operators (85) and (84) toπ(π₯). So far we have demonstrated as Gauss-Weierstrass transform (A.1) describe the action of πΜ on the initial probability amplitude and how the shift (84) finally yields the explicit form of the solution of SchrΒ¨odinger equation. It means that the result of the action of evolution operator (83) on π(π₯) is the product of combined action Μ Μ and heat propagation operator π, of translation operator Ξ representing the evolution operator of the free particle. Now let us consider another Fokker-Plank type equation, that is, the following example: ππ‘ πΉ (π₯, π‘) = πΌππ₯2 πΉ (π₯, π‘) + π½π₯ππ₯ πΉ (π₯, π‘) ,
(89)
with initial condition (55). Proceeding along the above outlined scheme of the solution of (54), we write the solution of (89) in general form (56): Μ (π₯) , πΉ (π₯, π‘) = ππ
Μ = exp (π΄ Μ + π΅) Μ , π
(90)
Μ and π΅Μ are defined as follows: where operators π΄ Μ = π‘πΌπ2 , π΄ π₯
π΅Μ = π‘π½π₯ππ₯ .
(91)
Μ and π΅Μ evidently do not commute: Quantities π΄ (83)
Then, on account of the substitution π‘ β ππ, π½ β βπ, operators Μ = exp (ππ2 π ) , Ξ π₯
Μ ππ Μ (π₯) , Ξ¨ (π₯, π‘) = exp (βπΞ¦ (π₯, π; π)) Ξ
Ξ¨ (π₯, π) = exp (βπΞ¦ (π₯, π; π)) π (π₯ + ππ2 , ππ) ,
π
πβππ‘ Ξ¨ (π₯, π‘) = β
arise. Thus, the solution of the SchrΒ¨odinger equation is a result of consequent action of the operator πΜ and further action of Μ on the initial condition function: Ξ
(84) (85)
Μ Μ π΅] Μ = 2πΌπ½π‘2 π2 = ππ΄, [π΄, π₯
π = 2π‘π½.
(92)
It allows disentanglement of the operators in the exponential according to the following rule: Μ + π΅) Μ = exp ( exp (π΄
1 β exp (βπ) Μ Μ . π΄) exp (π΅) π
(93)
The Scientific World Journal
7
Μ action on π(π₯), given below: Thus, the evolution operator π Μ (π₯) = exp (ππ2 ) exp (π½π‘π₯ππ₯ ) π (π₯) ππ π₯ = exp (πππ₯2 ) π (ππ½π‘ π₯) ,
(94)
simply reduces to a Gauss transforms, where the parameter π = π(π‘) reads as follows: (1 β exp (β2π½π‘)) πΌ . π½
π (π‘) =
(95)
Upon the trivial change of variables, we obtain Μ (π₯) = ππ Μ (π¦) , ππ
π¦ = π₯ exp (π) ,
π πΜ = exp ( ππ¦2 ) , 2 (96)
where π (π‘) =
πΌ 2π½π‘ (π β 1) π½
(97)
and, eventually, we end up with the following simple solution of (89): πΉ (π₯, π‘) =
β 2 1 β« πβ(exp(π½π‘)π₯βπ) /2π π (π) ππ. β2ππ ββ
(98)
The same example of the initial Gaussian function π(π₯) = exp(βπ₯2 ) as in the case of (54) yields (compare with (71)) πΉ(π₯, π‘)|π(π₯)=exp(βπ₯2 ) =
π2π½π‘ π₯2 1 ), exp (β 1 + 2π (π‘) β1 + 2π (π‘) (99)
where π is defined in (97). Note that we can meet the following modified form of Fokker-Plank type equation (89): ππ‘ πΉ (π₯, π‘) = πΌππ₯2 πΉ (π₯, π‘) + π½ππ₯ π₯πΉ (π₯, π‘) ,
=
Operational method is fast and universal mathematical tool for obtaining solutions of differential equations. Combination of operational method, integral transforms, and theory of special functions together with orthogonal polynomials closely related to them provides a powerful analytical instrument for solving a wide spectrum of differential equations and relevant physical problems. The technique of inverse operator, applied for derivatives of various orders and combined with integral transforms, allows for easy and straightforward solutions of various types of differential equations. With operational approach, we developed the methodology of inverse differential operators and derived a number of operational identities with them. We have demonstrated that using the technique of inverse derivatives and inverse differential operators, combined with exponential operator, integral transforms, and special functions, we can make significant progress in solution of various mathematical problems and relevant physical applications, described by differential equations.
Appendix In complete analogy with the heat equation solution by Gauss-Weierstrass transform [16]: exp (πππ₯2 ) π (π₯) =
2 β 1 (π₯ β π) ) π (π) ππ, β« exp (β 4π 2βππ ββ (A.1)
accounting for noncommutative relation for operators of inverse derivative π·π₯β1 and πΏ π·π₯ , defined through the operational relation (14) and accounting for (31) and (32), we write the solution of (40) in the following form: 2
(100)
in problems, related to propagation of electron beams in accelerators. Its solution arises from (99) immediately and differs from it just by a factor exp(π½π‘), as written below for a Gaussian π(π₯): πΉ(π₯, π‘)|π(π₯)=exp(βπ₯2 ) =
5. Conclusions
π2π½π‘ π₯2 ππ½π‘ exp (β ) 1 + 2π (π‘) β1 + 2π (π‘)
Μβ1 β π) β (π· 1 π₯ ) 1π (π) ππ, (A.2) π (π₯, π‘) = β« exp (β 4π‘ 2βππ‘ ββ where π is the image (32) of the initial condition function π. The kernel of the integral in the above formula can be expanded into series of two-variable Hermite polynomials π»π (π₯, π¦): 2
exp (β
π₯2 1 ), exp (β π (π‘) βπ (π‘)
π½ β2π½π‘ πΌ π (π‘) = 2 (1 β πβ2π½π‘ + π ). π½ 2πΌ
2
Μβ1 β π) (π· π₯ 4π‘
)1
Μβπ π· π 1 π₯ π»π ( , β ) 2π‘ 4π‘ π=0 π! β
= πβπ /4π‘ β (101)
However, differently from the solution of (54), where we had consequent transforms of the initial condition function π(π₯) by operators of translation and heat diffusion πΜ (63) (see also [15]), here we have just the action of πΜ alone with much more complicated dependence of the solution on time.
2
(A.3)
β
π₯π π 1 π» ( ,β ). 2 π 2π‘ 4π‘ π=0 (π!)
= πβπ /4π‘ β
Taking into account formula (12) for π₯π π»π (π/2π‘, β1/4π‘) in the operational identity above, which can be viewed as a generating function in terms of inverse derivative for
8
The Scientific World Journal
π»π (π₯, π¦), we obtain the following expansion for the kernel of the integral: exp (β
Μβ1 β π) (π· π₯ 4π‘
2
β
1 π₯π π₯2 π» ( ,β ), 2 π 2π‘ 4π‘ π=0 (π!) (A.4)
) 1 = πβπ /4π‘ β
where the series of Hermite polynomials of two arguments π»π (π₯, π¦) can be expressed in terms of Hermite-BesselTricomi functions π»πΆπ (π₯, π¦)βgeneralization of BesselTricomi functionsβand related to Bessel-Wright functions and to common Bessel functions [17]: β
π»π (βπ₯, π¦) , π=0 π! (π + π)!
π» πΆπ (π₯, π¦) = β
π β π0 .
(A.5)
In particular, for π = 0, we immediately find our series: π» πΆ0
β
(π₯, π¦) = β π=0
π»π (βπ₯, π¦) (π!)2
,
(A.6)
which obviously leads to the solution of the heat-type propagation problem (40) by the following appropriate Gauss transform: π (π₯, π‘) =
β 2 βππ₯ π₯2 1 , β ) π (π) ππ. β« πβπ /4π‘ π»πΆ0 ( 2π‘ 4π‘ 2βππ‘ ββ (A.7)
Conflict of Interests The author declares that there is no conflict of interests regarding the publication of this paper.
References [1] A. App`el and J. KampΒ΄e de FΒ΄eriet, Fonctions HypergΒ΄eomΒ΄etriques et HypersphΒ΄eriques: PolynΛomes dβHermite, Gauthier-Villars, Paris, France, 1926. [2] D. T. Haimo and C. Markett, βA representation theory for solutions of a higher order heat equation, I,β Journal of Mathematical Analysis and Applications, vol. 168, no. 1, pp. 89β107, 1992. [3] D. T. Haimo and C. Markett, βA representation theory for solutions of a higher order heat equation, II,β Journal of Mathematical Analysis and Applications, vol. 168, no. 2, pp. 289β 305, 1992. [4] G. Dattoli, βGeneralized polynomials, operational identities and their applications,β Journal of Computational and Applied Mathematics, vol. 118, no. 1-2, pp. 111β123, 2000. [5] G. Dattoli, H. M. Srivastava, and K. Zhukovsky, βOrthogonality properties of the Hermite and related polynomials,β Journal of Computational and Applied Mathematics, vol. 182, no. 1, pp. 165β 172, 2005. [6] A. ErdΒ΄elyi, W. Magnus, F. Oberhettinger, and F. G. Tricomi, Higher Transcendental Functions, vol. 2, McGraw-Hill Book Company, New York, NY, USA, 1953. [7] G. Dattoli, B. Germano, M. R. Martinelli, and P. E. Ricci, βNegative derivatives and special functions,β Applied Mathematics and Computation, vol. 217, no. 8, pp. 3924β3928, 2010.
[8] G. Dattoli, M. Migliorati, and S. Khan, βSolutions of integrodifferential equations and operational methods,β Applied Mathematics and Computation, vol. 186, no. 1, pp. 302β308, 2007. [9] H. M. Srivastava and H. L. Manocha, A Treatise on Generating Functions, Mathematics and its Applications (Ellis Horwood Ltd), Halsted Press, John Wiley and Sons, Chichester, UK, 1984. [10] G. Dattoli, H. M. Srivastava, and K. Zhukovsky, βA new family of integral transforms and their applications,β Integral Transforms and Special Functions, vol. 17, no. 1, pp. 31β37, 2006. [11] H. W. Gould and A. T. Hopper, βOperational formulas connected with two generalizations of Hermite polynomials,β Duke Mathematical Journal, vol. 29, no. 1, pp. 1β174, 1962. [12] G. Dattoli, H. M. Srivastava, and K. Zhukovsky, βOperational methods and differential equations with applications to initialvalue problems,β Applied Mathematics and Computation, vol. 184, no. 2, pp. 979β1001, 2007. [13] G. Dattoli and K. V. Zhukovsky, βApp`el polynomial series expansion,β International Mathematical Forum, vol. 5, no. 14, pp. 649β6662, 2010. [14] G. Dattoli, M. Migliorati, and K. Zhukovsky, βSummation formulae and stirling numbers,β International Mathematical Forum, vol. 4, no. 41, pp. 2017β22040, 2009. [15] J. E. Avron and I. W. Herbst, βSpectral and scattering theory of SchrΒ¨odinger operators related to the stark effect,β Communications in Mathematical Physics, vol. 52, no. 3, pp. 239β254, 1977. [16] K. B. Wolf, Integral Transforms in Science and Engineering, Plenum Press, New York, NY, USA, 1979. [17] G. N. Watson, A Treatise on the Theory of Bessel Functions, Cambridge University Press, Cambridge, UK, 2nd edition, 1944.