## Lie symmetries of fundamental solutions of one (2+1)-dimensional ultra-parabolic Fokker-Planck-Kolmogorov equation

Sergii Kovalenko, 1 Valeriy Stogniy 2 and Maksym Tertychnyi 3

- 1
- Department of Physics, Faculty of Oil, Gas and Nature Engineering, Poltava National Technical University, Poltava, Ukraine

(e-mail: kovalenko@imath.kiev.ua)

2 Department of Mathematical Physics, Faculty of Physics and Mathematics, National Technical University of Ukraine 'Kyiv Polytechnic Institute', Kyiv, Ukraine (e-mail: valeriy stogniy@mail.ru)

- 3 Department of Mathematics and Statistics, Faculty of Science, University of Calgary, Calgary, Canada

(e-mail: maksym.tertychnyi@gmail.com)

## Abstract

A (2+1)-dimensional linear ultra-parabolic Fokker-Planck-Kolmogorov equation is investigated from the group-theoretical point of view. By using the Berest-Aksenov approach, an algebra of invariance of fundamental solutions of the equation is found. A fundamental solution of the equation under study is computed in an explicit form as a weak invariant solution.

Keywords: Fokker-Planck-Kolmogorov equation, fundamental solution, Lie symmetries. 2010 Mathematical Subject Classification: 22E70, 35K70, 35Q84.

## 1 Introduction

The idea that a solution of any well-posed boundary problem for a particular linear partial differential equation (PDE) can be reduced to the construction of some solution of a special type (known as a fundamental solution) is one of the most powerful approaches in a classical mathematical physics. The method of integral transformations is a powerful and well-developed tool for the construction of the solutions of such a type. Unfortunately, this method is efficient only for linear PDE's with constant(or, at least, analytic) coefficients. Moreover, in some cases it is very difficult to investigate qualitative properties of fundamental solutions, because they can not be presented in an explicit form, but only in the form of inverse integral transformations. As a result, the development of more direct methods for the construction of fundamental solutions is an important problem of a modern mathematical physics. Especially, it is urgent in the case of PDE's with alternating coefficients. One of such modern methods is a classical Lie method for investigation of symmetrical properties of PDE's.

Group analysis of differential equations is a mathematical theory with an area of interest in the symmetry properties of differential equations. Basics of group analysis are in fundamental works of S. Lie and his scholars (see, for instance, [1]). Namely, Lie developed and was the first to use the tools of symmetry reduction. The main idea of this method is to search for the solution of the equation under consideration in the form of a special substitution (ansatz), that reduces the given equation to the differential equation with less number of independent variables.

Further development of group-theoretical methods of differential equations is mainly connected with works of the following mathematicians: G. Bikhoff [2], L. I. Sedov [3], A. J. A. Morgan [4], L. V. Ovsiannikov [5], N. H. Ibragimov [6-8], P. J. Olver [9,10], W. I. Fushchich [11-13] and others. At present, symmetry properties of many well-known equations of mechanics, gas dynamics, quantum physics, etc were investigated. Detailed analysis of results of research of symmetry properties of a wide range of linear and non-linear differential equations can be found in the monographs [5,6,11-16].

It is well-known, that as a rule fundamental solutions of linear PDE's are invariant with respect to the transformations, admitted by the prescribed equation [7, 17-25]. In particular, fundamental solutions of classical equations of mathematical physics such as the Laplace equa-

tion, the Heat equation, the Wave equation, etc have this property (see, for example, [18]). It means, that if a linear PDE (especially in case of an equation with alternating coefficients) has non-trivial symmetry properties, then we can use group-theoretical methods to construct its fundamental solution.

The object of our investigation is the linear (2+1)-dimensional ultra-parabolic FokkerPlanck-Kolmogorov equation

$$u _ { t } - u _ { x x } + x u _ { y } = 0, \quad ( x, y, t ) \in \mathbb { R } ^ { 3 },$$

where u = u t, x, y ( ), u t = ∂u ∂t , u y = ∂u ∂y , u xx = ∂ 2 u ∂x 2 .

Eq. (1) is the simplest and lowest dimensional version of the following (2 n +1)-dimensional linear ultra-parabolic equation

$$u _ { t } - \sum _ { j, l = 1 } ^ { n } ( k _ { j l } ( t, x, y ) u ) _ { x _ { j } x _ { l } } + \sum _ { j = 1 } ^ { n } ( f _ { j } ( t, x, y ) u ) _ { x _ { j } } + \sum _ { j = 1 } ^ { n } x _ { j } u _ { y _ { j } } = 0, \quad ( x, y, t ) \in \mathbb { R } ^ { 2 n + 1 }. \quad ( 2 ) \\ \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \end{$$

A.N. Kolmogorov introduced (2) in 1934 to describe the probability density of a system with 2 n degrees of freedom [26]. Here, the 2 n -dimensional space is the phase space, x = ( x , . . . , x 1 n ) is the velocity and y = ( y , . . . , y 1 n ) is the position of the system. It should be also stressed that Eq. (1) arises in mathematical finance in some generalization of the Black-Scholes model (see, for instance, [27, 28]).

Research of symmetry properties of Eq. (1) was commenced in work [29], where several point transformations of symmetries were found. But complete group analysis for this equation was not performed.

Maximal invariance algebra of Eq. (1) in Lie sense was calculated in work [30]. In addition, in this article an optimal system of sub-algebras was found for the calculated invariance algebra. Using two-dimensional sub-algebras, it was also performed the symmetry reduction and constructed several explicit invariant solutions of Eq. (1).

In this article, we continue research of symmetry properties of Eq. (1). The goal of our work is to find an invariance algebra of fundamental solutions of this equation and to construct in an explicit form the fundamental solution of Eq. (1) using already found symmetry algebra.

## 2 Symmetries of fundamental solutions of linear PDE's

Consider linear homogeneous PDE of p -th order with m independent variables

$$L u \equiv \sum _ { | \alpha | = 0 } ^ { p } A _ { \alpha } ( x ) D ^ { \alpha } u = 0, \quad x \in \mathbb { R } ^ { m }. \\. \. \. \. \. \. \. \. \. \. \. \. \. \.$$

In (3) we use standard notations: x = ( x , . . . , x 1 m ), α = ( α , . . . , α 1 m ) is a multi-index with integer non-negative components, | α | = α 1 + . . . + α m ;

$$D ^ { \alpha } \equiv \left ( \frac { \partial } { \partial x ^ { 1 } } \right ) ^ { \alpha _ { 1 } } \dots \left ( \frac { \partial } { \partial x ^ { m } } \right ) ^ { \alpha _ { m } } ;$$

A α ( x ) are some smooth functions of the variable x .

Fundamental solution of Eq. (3) is a function u x, x ( 0 ) (namely, generalized), that yields the following equation

$$L u = \delta ( x - x _ { 0 } ),$$

where δ x ( -x 0 ) is the Dirac delta function.

Standard methods to find fundamental solutions of linear PDE's are the method of integral transformations (especially, in case of the equations with constant coefficients), the Green's functions method, etc. [31,32]. Here we consider an algorithm to find fundamental solutions of linear homogeneous PDE's by using symmetry groups of this equation.

Remind, that a non-degenerate local substitution of the variables x, u

$$\bar { x } ^ { i } = f ^ { i } ( x, u, a ), \ \bar { u } = g ( x, u, a ), \quad i = 1, \dots, m,$$

depending on a continuous parameter a is called a symmetry transformation of Eq. (3), if this equation does not change its form with respect to the new variables ¯ and ¯. x u A set G of all such transformations forms a Lie group (local, more precisely), that is called a symmetry group (or an acceptable group) of Eq. (3).

According to the Lie theory, the construction of the symmetry group G of Eq. (3) is equivalent to finding its infinitesimal transformations:

$$\bar { x } ^ { i } \approx x ^ { i } + a \cdot \xi ^ { i } ( x, u ), \ \bar { u } \approx u + a \cdot \eta ( x, u ), \ \ i = 1, \dots, m,$$

where x i ( x, u ) and η x, u ( ) are some smooth functions.

Linear differential operator of the first order

$$X = \sum _ { i = 1 } ^ { m } \xi ^ { i } ( x, u ) \frac { \partial } { \partial x ^ { i } } + \eta ( x, u ) \frac { \partial } { \partial u } \\ \prod _ { \substack { \text{immatrix} \text{convert} } } \, \text{for the amount } \mathcal { T } \text{ to one of $Y$ is also called a minimum} }.$$

is called an infinitesimal operator of the group G . The operator X is also called a symmetry operator of Eq. (3).

The group transformations (5) corresponding to the infinitesimal transformations (6) with the operator (7) are found using so called Lie equations:

$$\frac { d \bar { x } ^ { i } } { d a } = \xi ^ { i } ( \bar { x }, \bar { u } ), \ \frac { d \bar { u } } { d a } = \eta ( \bar { x }, \bar { u } ), \ \ i = 1, \dots, m$$

with the initial conditions

$$\bar { x } ^ { i } | _ { a = 0 } = x ^ { i }, \ \bar { u } | _ { a = 0 } = u.$$

So called infinitesimal criterion of invariance plays a fundamental role in the symmetry analysis of differential equation, that is in case of Eq. (3) can be formulated in the following form.

Theorem 1. The infinitesimal operator (7) is a symmetry operator of Eq. (3) , if and only if there exists such a function λ = λ x ( ) , that yields the following identities:

$$X ( L u ) \equiv \lambda ( x ) \cdot L u$$

for any function u = u x ( ) from the domain of Eq. (3) .

In Eq. (8), X p is the prolongation of p -th order of the infinitesimal operator (7), that is calculated by the well-known formula [5,9,13,16]:

$$X = X + \sum _ { k = 1 } ^ { p } \sum _ { i _ { 1 }, \dots, i _ { k } = 1 } ^ { m } \varphi ^ { i _ { 1 } \cdots i _ { k } } \frac { \partial } { \partial u _ { i _ { 1 } \cdots i _ { k } } },$$

where

$$\varphi ^ { i _ { 1 } \cdots i _ { k } } = D _ { i _ { 1 } } \cdots D _ { i _ { k } } \left ( \eta - \sum _ { i = 1 } ^ { m } \xi ^ { i } u _ { i } \right ) + \sum _ { i = 1 } ^ { m } \xi ^ { i } u _ { i _ { 1 } \cdots i _ { k } i }, \\ i _ { 1 }, \dots, i _ { k } = 1, \dots, m, \ \ k = 1, \dots, p.$$

Here we denote by D i the operator of total differentiation with respect to the variable x i :

$$D _ { i } = \frac { \partial } { \partial x ^ { i } } + u _ { i } \frac { \partial } { \partial u } + \sum _ { j = 1 } ^ { \infty } \sum _ { i _ { 1 }, \dots, i _ { j } = 1 } ^ { m } u _ { i _ { 1 } \cdots i _ { j } i } \frac { \partial } { \partial u _ { i _ { 1 } \cdots i _ { j } } }, \quad i = 1, \dots, m.$$

It was shown in [33], that the linear homogeneous PDE (3) given additional conditions p ≥ 2 and m ≥ 2 admits symmetry operators only of such a form

$$X = \sum _ { i = 1 } ^ { m } \xi ^ { i } ( x ) \frac { \partial } { \partial x ^ { i } } + \left ( \alpha ( x ) u + \beta ( x ) \right ) \frac { \partial } { \partial u }. \\ \inf i n t i o & \text{internal operators of the form } \mathbb { m } ( \mathbb { m } \text{ }the maximal invariance also for } \mathbb { e }$$

In the class of infinitesimal operators of the form (9), the maximal invariance algebra of Eq.

- (3) as a vector space is a direct sum of two sub-algebras: the sub-algebra, that consists of the operators of the form

and the infinite-dimensional subalgebra, that is generated by the operators

$$X & = \sum _ { i = 1 } ^ { m } \xi ^ { i } ( x ) \frac { \partial } { \partial x ^ { i } } + \alpha ( x ) u \frac { \partial } { \partial u }, \\ \text{nsional subalgebra, that is generated by the operators}$$

$$X = \beta ( x ) \frac { \partial } { \partial u },$$

where β = β x ( ) is an arbitrary smooth solution of Eq. (3).

It is clear, that the infinitesimal operators (11) are symmetry operators of Eq. (4). Hence, in the sequel we will consider only operators of the form given by (10).

Constructive method to find symmetries of the form (10) of linear inhomogeneous PDE's with δ -function in a right-hand side was proposed in works [17,20]. In the same articles, it was introduced an algorithm for construction of invariant fundamental solutions of the equations of the form (3).

Main result of these works is in the following statement.

Theorem 2. The Lie algebra of symmetry operators of the form (10) of Eq. (4) is a subalgebra of the Lie algebra of symmetry operators of Eq. (3) , which is defined by the following conditions:

$$\xi ^ { i } ( x _ { 0 } ) = 0,$$

where i = 1 , . . . , m .

$$\lambda ( x _ { 0 } ) + \sum _ { i = 1 } ^ { m } \frac { \partial \xi ^ { i } ( x _ { 0 } ) } { \partial x ^ { i } } = 0,$$

Formulate the algorithm to find fundamental solutions of a linear PDE using properties of its invariance algebra:

- 1. find a general form of symmetry operator of Eq. (3) and a relevant function λ x ( ), that yields identities (8);

- 2. using conditions (12) and (13), find a Lie algebra of symmetry operators of Eq. (4);
- 3. construct invariant fundamental solutions of Eq. (3) using symmetry operators of Eq. (4).

Remark 1. Formulated algorithm for construction of fundamental solutions using symmetries of linear PDE's with δ -function in a right-hand side is especially efficient for multi-dimensional linear equations and for the equations with alternating coefficients in case, if they allow rather wide invariance algebras.

Remark 2. To find generalized invariant fundamental solutions, it is necessary to solve reduced equations (these equations are written in invariants of corresponding transformation groups) on the set of generalized functions (see, for example, [19,20,34]).

## 3 Symmetries of fundamental solutions of Eq. (1)

Apply the method described in the previous section to the object of our research, in other words to the linear Fokker-Planck-Kolmogorov equation (1).

We define fundamental solution of Eq. (1) as a generalized function u = u t, x, y, t ( 0 , x 0 , y 0 ), that depends on t 0 , x 0 , y 0 as parameters and yields the equation

$$u _ { t } - u _ { x x } + x u _ { y } = \delta ( t - t _ { 0 }, x - x _ { 0 }, y - y _ { 0 } ),$$

given the additional condition u | t&lt;t 0 = 0.

In work [30], it was shown that the maximal invariance algebra of Eq. (1) is generated by the following infinitesimal operators:

$$X _ { 1 } = \partial _ { x } + t \partial _ { y }, \ X _ { 2 } = 2 t \partial _ { t } + x \partial _ { x } + 3 y \partial _ { y } - 2 u \partial _ { u },$$

$$X _ { 3 } = t ^ { 2 } \partial _ { t } + ( t x + 3 y ) \partial _ { x } + 3 t y \partial _ { y } - ( 2 t + x ^ { 2 } ) u \partial _ { u },$$

$$X _ { 4 } = 3 t ^ { 2 } \partial _ { x } + t ^ { 3 } \partial _ { y } + 3 ( y - t x ) u \partial _ { u }, \ X _ { 5 } = 2 t \partial _ { x } + t ^ { 2 } \partial _ { y } - x u \partial _ { u },$$

$$X _ { 6 } = \partial _ { t }, \ X _ { 7 } = \partial _ { y }, \ X _ { 8 } = u \partial _ { u }, \ X _ { \infty } = \beta ( t, x, y ) \partial _ { u }.$$

In the last operator, the function β = β t, x, y ( ) is an arbitrary smooth solution of the equation under study.

Theorem 3. Eq. (14) admits an infinite-dimensional Lie algebra of symmetry operators with the following basis of finite-dimensional part:

$$\begin{array} { c c c c c } \cdot & \cdot & \cdot & \cdot \\ & Y _ { 1 } = 2 ( t - t _ { 0 } ) \partial _ { t } + ( x - x _ { 0 } ) \partial _ { x } - ( x _ { 0 } ( t - t _ { 0 } ) - 3 ( y - y _ { 0 } ) ) \partial _ { y } - 4 u \partial _ { u }, \\ & Y _ { 2 } = ( t ^ { 2 } - t _ { 0 } ^ { 2 } ) \partial _ { t } + ( ( t x + 3 y ) - ( t _ { 0 } x _ { 0 } + 3 y _ { 0 } ) ) \partial _ { x } + ( 3 ( y - y _ { 0 } ) t - t _ { 0 } x _ { 0 } ( t - t _ { 0 } ) ) \partial _ { y } - \\ & & - ( 2 ( t - t _ { 0 } ) + x ^ { 2 } - x _ { 0 } ^ { 2 } ) u \partial _ { u }, \\ & Y _ { 3 } = 3 ( t ^ { 2 } - t _ { 0 } ^ { 2 } ) \partial _ { x } + ( t ^ { 3 } - 3 t _ { 0 } ^ { 2 } t + 2 t _ { 0 } ^ { 3 } ) \partial _ { y } - 3 ( t x - y - ( t _ { 0 } x _ { 0 } - y _ { 0 } ) ) u \partial _ { u }, \\ & Y _ { 4 } = 2 ( t - t _ { 0 } ) \partial _ { x } + ( t - t _ { 0 } ) ^ { 2 } \partial _ { y } - ( x - x _ { 0 } ) u \partial _ { u }. \\ \prod _ { \infty } \tilde { \ } W i t a \ t h a \ a n o n a l \ t h a \ a n o n a l \ t h a \ a n o n a l \ t h a \ a n o n a l \ t h a \ a n o n a l \ t h a \ a n o n a l \ t h a \ a n o n a l \ t h a \ a n o n a l \ t h a \ a n o n a l \ t h a \ a n o n a l \ t h a \ a n o n a l \ t h a \ a n o n a l \ t h a \ a n o n a l \ t h a \ a n o n a l \ t h a \ a n o n a l \ t h a \ a n o n a l \ t h a \ a n o n a l \ t h a \ a n o n a l \ t h a \ a n o n a l \ a n o n a l \ t h a \ a n o n a l \ t h a \ a n o n a l \ t h a \ a n o n a l \ t h a \ a n o n a l \ t h a \ a n o n a l \ t h a \ a n o n a l \ t h a \ a n o n a l \ t h a \ a n o n a l \ t h a \ a n o n a l \ t h a \ a n o n a l \ t h a \ a n o n a l \ t h a \ a n o n a l \ t h a \ a n o n a l \ t h a \ a n o n a l \ t h a \ a n o n a l \ t h a \ a n o n a l \ t h a \ a n o n a l \ t h a \ t h a \ a n o n a l \ t h a \ a n o n a l \ t h a \ a n o n a l \ t h a \ a n o n a l \ t h a \ a n o n a l \ t h a \ a n o n a l \ t h a \ a n o n a l \ t h a \ a n o n a l \ t h a \ a n o n a l \ t h a \ a n o n a l \ t h a \ a n o n a l \ t h a \ a n o n a l \ t h a \ a n o n a l \ t h a \ a n o n a l \ t h a \ a n o n a l \ t h a \ a n o n a l \ t h a \ a n o n a l \ t h a \ a n o n a l$$

Proof. Write the general form of symmetry operator of the finite-dimensional part of invariance algebra of Eq. (1)

or in a more detailed form

$$X = \sum _ { i = 1 } ^ { 8 } a _ { i } X _ { i },$$

$$X = & ( 2 a _ { 2 } t + a _ { 3 } t ^ { 2 } + a _ { 6 } ) \partial _ { t } + ( a _ { 1 } + a _ { 2 } x + a _ { 3 } ( t x + 3 y ) + 3 a _ { 4 } t ^ { 2 } + 2 a _ { 5 } t ) \partial _ { x } + \\ & + ( a _ { 1 } t + 3 a _ { 2 } y + 3 a _ { 3 } t y + a _ { 4 } t ^ { 3 } + a _ { 5 } t ^ { 2 } + a _ { 7 } ) \partial _ { y } + \\ & + ( - 2 a _ { 2 } - a _ { 3 } ( 2 t + x ^ { 2 } ) + 3 a _ { 4 } ( y - t x ) - a _ { 5 } x + a _ { 8 } ) u \partial _ { u },$$

where a i ( i = 1 , . . . , 8) are any real constants.

Substituting the infinitesimal operator (15) into Eq. (8), where we put Lu ≡ u t -u xx + xu y , p = 2, we find the function λ = λ t, x, y ( ), that corresponds to this operator:

$$\lambda ( t, x, y ) = - 4 a _ { 2 } - a _ { 3 } ( 4 t + x ^ { 2 } ) + 3 a _ { 4 } ( y - t x ) - a _ { 5 } x + a _ { 8 }.$$

Substituting (15) and (16) into Eqs. (12) and (13) (see, Th. 1), we obtain the following equalities:

$$& 2 a _ { 2 } t _ { 0 } + a _ { 3 } t _ { 0 } ^ { 2 } + a _ { 6 } = 0 ; \\ & a _ { 1 } + a _ { 2 } x _ { 0 } + a _ { 3 } ( t _ { 0 } x _ { 0 } + 3 y _ { 0 } ) + 3 a _ { 4 } t _ { 0 } ^ { 2 } + 2 a _ { 5 } t _ { 0 } = 0 ; \\ & a _ { 1 } t _ { 0 } + 3 a _ { 2 } y _ { 0 } + 3 a _ { 3 } t _ { 0 } y _ { 0 } + a _ { 4 } t _ { 0 } ^ { 3 } + a _ { 5 } t _ { 0 } ^ { 2 } + a _ { 7 } = 0 ; \\ & 2 a _ { 2 } - a _ { 3 } ( 2 t _ { 0 } - x _ { 0 } ^ { 2 } ) + 3 a _ { 4 } ( y _ { 0 } - t _ { 0 } x _ { 0 } ) - a _ { 5 } x _ { 0 } + a _ { 8 } = 0,$$

from which we easily derive

$$a _ { 1 } = - a _ { 2 } x _ { 0 } - a _ { 3 } ( t _ { 0 } x _ { 0 } + 3 y _ { 0 } ) - 3 a _ { 4 } t _ { 0 } ^ { 2 } - 2 a _ { 5 } t _ { 0 } ;$$

$$a _ { 6 } = - 2 a _ { 2 } t _ { 0 } - a _ { 3 } t _ { 0 } ^ { 2 } ;$$

$$a _ { 7 } = a _ { 2 } ( t _ { 0 } x _ { 0 } - 3 y _ { 0 } ) + a _ { 3 } t _ { 0 } ^ { 2 } x _ { 0 } + 2 a _ { 4 } t _ { 0 } ^ { 3 } + a _ { 5 } t _ { 0 } ^ { 2 } ;$$

$$a _ { 8 } = - 2 a _ { 2 } + a _ { 3 } ( 2 t _ { 0 } - x _ { 0 } ^ { 2 } ) - 3 a _ { 4 } ( y _ { 0 } - t _ { 0 } x _ { 0 } ) + a _ { 5 } x _ { 0 }.$$

$$& - a _ { 3 } ( t _ { 0 } x _ { 0 } + 3 y _ { 0 } ) - 3 a _ { 4 } t _ { 0 } ^ { 2 } - 2 a _ { 5 } t _ { 0 } ; \\ & - a _ { 3 } t _ { 0 } ^ { 2 } ; \\ & - 3 y _ { 0 } ) + a _ { 3 } t _ { 0 } ^ { 2 } x _ { 0 } + 2 a _ { 4 } t _ { 0 } ^ { 3 } + a _ { 5 } t _ { 0 } ^ { 2 } ; \\ & \cdot a _ { 3 } ( 2 t _ { 0 } - x _ { 0 } ^ { 2 } ) - 3 a _ { 4 } ( y _ { 0 } - t _ { 0 } x _ { 0 } ) + a _ { 5 } x _ { 0 }. \\ & \ n t c \ a. \ a. \ a. \ \text{ and } a. \ \text{br the calculated over}.$$

Substituting in (15) the constants a , a 1 6 , a 7 , and a 8 by the calculated expressions and having split with respect to the independent constants a , a 2 3 , a 4 , a 5 , we obtain that the finitedimensional part of invariance algebra of Eq. (14) is four-dimensional and generated by the following operators:

$$\text{nowing operators} \colon \\ Y _ { 1 } & = 2 ( t - t _ { 0 } ) \partial _ { t } + ( x - x _ { 0 } ) \partial _ { x } - ( x _ { 0 } ( t - t _ { 0 } ) - 3 ( y - y _ { 0 } ) ) \partial _ { y } - 4 u \partial _ { u }, \\ Y _ { 2 } & = ( t ^ { 2 } - t _ { 0 } ^ { 2 } ) \partial _ { t } + ( ( t x + 3 y ) - ( t _ { 0 } x _ { 0 } + 3 y _ { 0 } ) ) \partial _ { x } + ( 3 ( y - y _ { 0 } ) t - t _ { 0 } x _ { 0 } ( t - t _ { 0 } ) ) \partial _ { y } - \\ & \quad - ( 2 ( t - t _ { 0 } ) + x ^ { 2 } - x _ { 0 } ^ { 2 } ) u \partial _ { u }, \\ Y _ { 3 } & = 3 ( t ^ { 2 } - t _ { 0 } ^ { 2 } ) \partial _ { x } + ( t ^ { 3 } - 3 t _ { 0 } ^ { 2 } t + 2 t _ { 0 } ^ { 3 } ) \partial _ { y } - 3 ( t x - y - ( t _ { 0 } x _ { 0 } - y _ { 0 } ) ) u \partial _ { u }, \\ Y _ { 4 } & = 2 ( t - t _ { 0 } ) \partial _ { x } + ( t - t _ { 0 } ) ^ { 2 } \partial _ { y } - ( x - x _ { 0 } ) u \partial _ { u }. \\ \text{Invariance property of Eq. } ( 14 ) \text{ with respect to the operators of the form } \beta ( t, x, y ) \partial _ { u }, \text{ where}$$

Invariance property of Eq. (14) with respect to the operators of the form β t, x, y ( ) ∂ u , where β t, x, y ( ) is an arbitrary smooth solution of Eq. (1) is straightforward.

$$\text{The proof is now completed.}$$

Show how the results of Th. 3 can be applied to the construction of invariant fundamental solutions of Eq. (1). Use for this, for example, the two-dimensional algebra of operators 〈 Y , Y 1 4 〉 . Classical invariants corresponding to these operators can be derived from the system of equations

$$Y _ { 1 } I = 0, \ Y _ { 4 } I = 0,$$

where I = I ( t, x, y, u ), and the equations are solved in a classical sense. We obtain after related calculations:

$$I _ { 1 } = ( t - t _ { 0 } ) ^ { 2 } \exp \left [ \frac { ( x - x _ { 0 } ) ^ { 2 } } { 4 ( t - t _ { 0 } ) } \right ] u, \ I _ { 2 } = \frac { ( t - t _ { 0 } ) ( x + x _ { 0 } ) - 2 ( y - y _ { 0 } ) } { ( t - t _ { 0 } ) ^ { 3 / 2 } }. \\ \intertext { a n d i n d } \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \ \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \cdots \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \pm \, \, \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \. \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \end{cases}$$

Classical invariant solution is found from the equality I 1 = ϕ I ( 2 ), from which we obtain a substitution (ansatz)

$$u = \frac { 1 } { ( t - t _ { 0 } ) ^ { 2 } } \exp \left [ - \frac { ( x - x _ { 0 } ) ^ { 2 } } { 4 ( t - t _ { 0 } ) } \right ] \varphi ( \omega ), \ \omega = I _ { 2 },$$

that reduces the equation under study (1) to the following ordinary differential equation:

$$\frac { d ^ { 2 } \varphi } { d \omega ^ { 2 } } + \frac { 3 } { 2 } \, \omega \frac { d \varphi } { d \omega } + \frac { 3 } { 2 } \, \varphi = 0.$$

It is well-known, that a particular solution of this equation is the following function [35, P. 216]:

$$\varphi = C \exp \left [ - \frac { 3 } { 4 } \, \omega ^ { 2 } \right ]. \quad \quad \quad$$

Substituting (18) to (17), we derive a classical invariant solution of Eq. (1) corresponding to the operators Y 1 and Y 4 :

$$u = \frac { C } { ( t - t _ { 0 } ) ^ { 2 } } \, \exp \left [ - \frac { ( x - x _ { 0$$

Substituting (19) to Eq. (14), it is easy to show that Lu t, x, y ( ) = 0. As a result, the solution (19) is not a fundamental solution of the Fokker-Planck-Kolmogorov equation (1). To construct its weak invariant solution we use Statement 1 from work [36], in other words we search for weak invariant solutions in the form

$$u = \frac { h ( t, x, y ) } { ( t - t _ { 0 } ) ^ { 2 } } \, \exp \left [ - \frac {$$

where h = h t, x, y ( ) ∈ D ′ ( R 3 ) (denote by D ′ = D ′ ( R 3 ) the space of generalized functions). The equations Y u 1 = 0 and Y u 4 = 0 give us correspondingly:

$$2 ( t - t _ { 0 } ) h _ { t } + ( x - x _ { 0 } ) h _ { x } - ( x _ { 0 } ( t -$$

$$2 ( t - t _ { 0 } ) h _ { x } + ( t - t _ { 0 } ) ^ { 2 } h _ { y } = 0.$$

$$x - x _ { 0 } ) h _ { x } & - ( x _ { 0 } ( t - t _ { 0 } ) - 3 ( y - y _ { 0 } ) ) h _ { y } = 0, \\ t - t _ { 0 } ) ^ { 2 } h _ { y } & = 0.$$

It is easy to see that the generalized function h t, x, y ( ) = C θ t 1 ( -t 0 ) + C 0 is a solution of these equations (here, θ t ( -t 0 ) is the Heaviside step function). We obtain using this fact:

$$u = \frac { C _ { 1 } \theta ( t - t _ { 0 } ) + C _ { 0 } } { ( t - t _ { 0 } ) ^ { 2 } } \, \exp \left [ - \frac { ( x - x _ { 0 } ) ^ { 2 } } { 4 ( t - t _ { 0 } ) } - \frac { 3 } { ( t - t _ { 0 } ) ^ { 3 } } \left ( y - y _ { 0 } - ( t - t _ { 0 } ) \frac { x + x _ { 0 } } { 2 } \right ) ^ { 2 } \right ]. \quad \ \ ( 2 0 )$$

Substituting (20) to Eq. (14), we find the constant C 1 = √ 3 2 π .

Hence, the fundamental solution of the Fokker-Planck-Kolmogorov equation (1)

$$u & = \frac { \frac { 3 } { 2 \pi } \, \theta ( t - t _ { 0 } ) + C _ { 0 } } { ( t - t _ { 0 } ) ^ { 2 } } \, \exp \left [ - \frac { ( x - x _ { 0 } ) ^ { 2 } } { 4 ( t - t _ { 0 } ) } - \frac { 3 } { ( t - t _ { 0 } ) ^ { 3 } } \left ( y - y _ { 0 } - ( t - t _ { 0 } ) \frac { x + x _ { 0 } } { 2 } \right ) ^ { 2 } \right ] \quad \ \ ( 2 1 ) \\. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \.. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \.$$

is found as a weak invariant solution with respect to a two-dimensional algebra 〈 Y , Y 1 4 〉 of point symmetries of Eq. (14). In formula (21) we can put C 0 = 0, because a fundamental solution is defined up to the addition of any solution of a homogeneous equation.

Remark 3. The fundamental solution (21) was found by A. N. Kolmogorov [37] without using the methods of symmetry analysis of differential equations. Our calculations give grouptheoretical background for this solution and justify once again an empiric observation, that fundamental solutions of linear PDE's should be searched among the invariant solutions.

Remark 4. Using a similar scheme, as it was performed for Eq. (1) , we can investigate the symmetry properties of fundamental solutions of other many-dimensional and more complicated Kolmogorov type equations (2) . In this article, we restricted our attention to Eq. (1) in order not to lose the main idea of Berest-Aksenov method of construction of fundamental solutions of linear homogeneous PDE's using its group properties under complicated technical calculations.

## 4 Conclusions

In this article, using the Berest-Aksenov method [17,20] we found an invariance algebra of fundamental solutions of linear Fokker-Planck-Kolmogorov equation (1). Its operators were used to construct invariant fundamental solutions of this equation. It was shown that the fundamental solution (21) of Eq. (1), found by A. N. Kolmogorov is a weak invariant fundamental solution.

## References

- [1] S. Lie: Theory of Transformation Groups , vol. 1-3, B. G. Teubner, Leipzig 1888-93 (in German).
- [2] G. Birkhoff: Hydrodynamics: A Study in Logic, Fact and similitude , 1st ed., Princeton University Press, Princeton, N.J. 1950.
- [3] L. I. Sedov: Similarity and Dimensional Methods, 4th ed., Academic Press, New York 1959.
- [4] A.J. A. Morgan: The reduction by one of the number of independent variables in some systems of partial differential equations, Quart. J. Math. Oxford 3 (12) (1952), 250-259.

- [5] L. V. Ovsiannikov: Group Analysis of Differential Equations , Academic Press, New York 1982.
- [6] N.H. Ibragimov: Transformation Groups Applied to Mathematical Physics , D. Reidel, Dordrecht 1985.
- [7] N.Kn. Ibragimov: Group analysis of ordinary differential equations and the invariance principle in mathematical physics (for the 150th anniversary of Sophus Lie), Russ. Math. Surv. 47 (4) (1992), 89-156.
- [8] Yu.N. Grigoriev, NH. Ibragimov, VF. Kovalev and S. V. Meleshko: Symmetries of Integro-Differential Equations: With Applicatons in Mechanics and Plasma Physics , Springer, Dordrecht 2010.
- [9] P. J. Olver: Applications of Lie Groups to Differential Equations , 1st ed., Springer-Verlag, New York 1986.
- [10] P. J. Olver: Equivalence, Invariants, and Symmetry , Cambridge University Press, Cambridge 1995.
- [11] W.I. Fushchich and A. G. Nikitin: Symmetries of Maxwell's Equations , D. Reidel, Dordrecht 1987.
- [12] W.I. Fushchich and A. G. Nikitin: Symmetries of Equations of Quantum Mechanics , Allerton Press Inc., New York 1994.
- [13] W.I. Fushchich, W. M. Shtelen and N. I. Serov: Symmetry analysis and Exact Solutions of Equations of Nonlinear Mathematical Physics , Kluwer, Dordrecht 1993.
- [14] V.K. Andreev, O.V. Kaptsov, V.V. Pukhnachov and A.A. Rodionov: Applications of Group-Theoretical Methods in Hydrodynamics , Kluwer, Dordrecht 1998.
- [15] W. Miller (Jr.): Symmetry and Separation of Variables , Addison-Wesley, Reading, Mass. 1977.
- [16] V.I. Lahno, S. V. Spichak and V. I. Stogniy: Symmetry Analysis of Evolution Type equations , RCD, Moskow-Izhevsk 2004.

- [17] A.V. Aksenov: Symmetries of linear partial differential equations and fundamental solutions, Dokl. Math. 51 (3) (1995), 329-331.
- [18] A.V. Aksenov: Symmetries of fundamental solutions of linear partial differential equations, in: Symmetries of Differential Equations (G. N. Yakovenko, ed.), pp. 6-35, MFTI, Moskow 2009.
- [19] Yu.Yu. Berest: Construction of fundamental solutions for Huygens' equations as invariant solutions, Sov. Math. Dokl. 43 (2) (1991), 496-499.
- [20] Yu.Yu. Berest: Group analysis of linear differential equations in distributions and the construction of fundamental solutions, Differ. Equations 29 (11) (1993), 1700-1711.
- [21] Yu. Berest and Yu. Molchanov: Fundamental solutions for partial differential equations with reflection group invariance, J. Math. Phys. 36 (8) (1995), 4324-4339.
- [22] R.K. Gazizov and N.H. Ibragimov: Lie symmetry analysis of differential equations in finance, Nonlinear Dynam. 17 (4) (1998), 387-407.
- [23] M. Craddock and E. Platen: Symmetry group methods for fundamental solutions, J. Differential Equations 207 (2) (2004), 285-302.
- [24] M. Craddock: Fundamental solutions, transition densities and the integration of Lie symmetries, J. Differential Equations 246 (6) (2009), 2538-2560.
- [25] M. Craddock: Lie symmetry methods for multi-dimensional parabolic PDEs and diffusions, J. Differential Equations 252 (1) (2012), 56-90.
- [26] A.N. Kolmogoroff: Zuf¨llige Bewegungen (Zur Theorie der Brownschen Bewegung), a Ann. Math. 35 (1) (1934) 116-117.
- [27] E. Barucci, S. Polidoro and V. Vespri: Some results on partial differential equations and Asian options, Math. Models Methods Appl. Sci. 11 (3) (2001), 475-497.
- [28] A. Pascussi: Kolmogorov equations in physics and in finance, Progres in Nonlinear Differential Equations and their Applications 63 (2005), 313-324.

- [29] E. Lanconelli and S. Polidoro: On a class of hypoelliptic evolution operators, Rend. Semin. Mat. Torino 52 (1) (1994), 29-63
- [30] S. V. Spichak, V. I. Stogniy and I. M. Kopas: Symmetry properties and exact solutions of the linear Kolmogorov equation, Research Bulletin of NTUU 'Kyiv Politechnic Institute' 2011 (4) (2011), 93-97 (in Ukrainian).
- [31] A. Friedman: Partial Differential Equations of Parabolic Type , Prentice-Hall, Englewood Cliffs, N.J. 1964.
- [32] L. C. Evans: Partial Differential Equations , 2sd ed., American Mathematical Society, Providence 2010.
- [33] G. Bluman: Simplifying the form of Lie groups admitted by a given differential equation, J. Math. Anal. Appl. 145 (1) (1990), 52-62.
- [34] N. Dapi´ c, M. Kunzinger and S. Pilipovi´ c: Symmetry group analysis of weak solutions, Proc. London Math. Soc. 84 (3) (2002), 686-710.
- [35] A.D. Polyanin and V. F. Zaitsev: Handbook of Exact Solutions for Ordinary Differential Equations , 2sd ed., Chapmam &amp; Hall/CRC, Boca Raton 2002.
- [36] Yu.Yu. Berest: Weak invariants of local trasformation groups, Differ. Equations 29 (10) (1993), 1561-1567.
- [37] A.N. Kolmogoroff: Zur theorie der steigen Zuf¨lligen prozesse, a Math. Ann. 108 (1) (1933), 149-160.
Sergii Kovalenko Valeriy Stogniy Maksym Tertychnyi
Lie symmetries of fundamental solutions of one (2+1)-dimensional ultra-parabolic Fokker--Planck--Kolmogorov equation
A (2+1)-dimensional linear ultra-parabolic Fokker--Planck--Kolmogorov
equation is investigated from the group-theoretical point of view. By using the
Berest--Aksenov approach, an algebra of invariance of fundamental solutions of
the equation is found. A fundamental solution of the equation under study is
computed in an explicit form as a weak invariant solution.