Next Article in Journal
On r-Noncommuting Graph of Finite Rings
Previous Article in Journal
Fixed Point Results for Frum-Ketkov Type Contractions in b-Metric Spaces
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Two Forms of An Inverse Operator to the Generalized Bessel Potential

1
Department of Mathematical Analysis, Algebra and Geometry, Chechen State University, Sheripov Str. 32, 364024 Grozny, Russia
2
Department of Mathematics, Physics, and Chemistry, Beuth Technical University of Applied Sciences Berlin, Luxemburger Str. 10, 13353 Berlin, Germany
3
Department of Mathematical and Applied Analysis, Voronezh State University, Universitetskaya pl., 1, 394018 Voronezh, Russia
4
Department of Applied Mathematics and Computer Modeling, Belgorod State National Research University (BelGU), Pobedy Street, 85, 308015 Belgorod, Russia
*
Author to whom correspondence should be addressed.
Axioms 2021, 10(3), 232; https://doi.org/10.3390/axioms10030232
Submission received: 28 July 2021 / Revised: 11 September 2021 / Accepted: 14 September 2021 / Published: 18 September 2021

Abstract

:
In this paper, we treat a convolution-type operator called the generalized Bessel potential. Our main result is the derivation of two different forms of its inversion. The first inversion is provided in terms of an approximative inverse operator using the method of an improving multiplier. The second one employs the regularization technique for the divergent integrals in the form of the appropriate segments of the Taylor–Delsarte series.

1. Introduction

The potential theory originates from the theory of electrostatic and gravitational potentials and the Laplace, wave, Helmholtz, and Poisson equations. The famous Riesz potentials are known to be the realizations of the real negative powers of the Laplace and wave operators. In the meantime, a lot of attention in the potential theory is given to the Bessel potential
G α f ( x ) = R n G α ( x y ) f ( y ) d y , α > 0 ,
where
G α ( x ) = 2 2 n α 2 π n 2 Γ α 2 K n α 2 ( | x | ) | x | n α 2
and K ν stands for the modified Bessel functions of the second kind. The operator G α can be interpreted as a realization of the real negative powers of the operator ( I Δ ) . In particular, the Bessel potential was treated by N. Aronzajn and K.T. Smith in [1] and by Calderon in [2]. It is worth mentioning that the Bessel potentials generate some special spaces of functions with a fractional smoothness α . They are very useful for the investigation of some fractional PDEs of the elliptic type. These spaces are an extension of the Sobolev spaces L p m ( R n ) to the case of the fractional order α . They are often referred to as Sobolev spaces of fractional order. The first results regarding the spaces of Bessel potentials were obtained by I. Stein in [3] in the case 0 < α < 2 and by Lizorkin in [4] in the general case. The inversion of the Bessel potentials was first derived by V.A. Nogin in [5,6] using the hypersingular integrals. Later on, V.S. Guliev and Z.V. Safarov studied the Bessel potentials generated by the Bessel differential operators in [7]. In particular, they proved the boundedness of these potentials in the weighted Lebesgue spaces and derived some embedding theorems in B k , n –Sobolev–Liouville spaces. In [8], the Young inequality for the B-convolution operators in the B-Bessel potential spaces was proved, and the Bessel potentials were characterized in terms of the B–Lizorkin–Triebel spaces. The optimal embedding of spaces of Bessel-type potentials were obtained in [9,10,11].
An inversion formula for the Bessel potential G α was derived in [12] using a special type of the weighted wavelet transforms. Further investigation of the hypersingular integrals associated to G α was performed in [13]. In [14], the fractional powers I k = 1 n 2 x k 2 + 2 ν x n x k α of the operator I k = 1 n 2 x k 2 + 2 ν x n x k , ν > 0 were introduced, and their approximation properties were studied as α tends to 0 + .
In this paper, we first construct the positive fractional powers of the operator ( I Δ γ ) , where Δ γ = k = 1 n ( B γ k ) x k is the Laplace–Bessel operator with ( B γ k ) x k = 2 x k 2 + γ k x k x k . The negative fractional powers of this operator are the generalized Bessel potentials considered in [15]. The most important difference between the generalized Riesz potential ( Δ γ ) α 2 , α > 0 , and the generalized Bessel potential ( I Δ γ ) α 2 is that the kernel of the generalized Bessel potential has an exponential decay at infinity.
It is worth mentioning that the integral equations involving the Bessel potentials and their generalizations are actively employed, say, for the investigation of the corresponding PDEs, such as the static Shrödinger equations with the critical and the super-critical exponents.
The generalized Bessel potential or Bessel B-potential is defined by the formula [15]
( G γ α φ ) ( x ) = ( G α γ φ ) γ = R + n G α γ ( y ) ( γ T x y φ ( x ) ) y γ d y , α > 0
with the kernel
G α γ ( x ) = 2 n | γ | α 2 + 1 | x | n + | γ | α 2 Γ α 2 i = 1 n Γ γ i + 1 2 K n + | γ | α 2 ( | x | )
and the generalized translation γ T x y defined by the Formula (17) (see the next section). Equation (1) realizes the negative fractional power ( I Δ γ ) α / 2 , α > 0 , where I is the identity operator, and Δ γ is the Laplace–Bessel operator
Δ γ = ( γ ) x = k = 1 n ( B γ k ) x k
and
( B γ ) t = 2 t 2 + γ t t = 1 t γ t t γ t , t > 0 , γ R
is the Bessel operator. The results derived in [15] allow us to represent the generalized Bessel potential in the form
G γ α φ = F γ 1 ( 1 + | x | 2 ) α 2 F γ φ , α > 0 ,
where F γ stands for the multidimensional Hankel transform defined by the Equation (23) (see the next section).
In this paper, we construct two different explicit forms of the operator inverse to Equation (1) in the weighted Lebesgue spaces, or, in other words, we get two representations of the positive fractional powers ( I Δ γ ) α / 2 , α > 0 . The operator ( I Δ γ ) α / 2 , α > 0 is first constructed by the approximative inverse operator (AIO) method, and then, another realization is obtained using the so-called B-hypersingular integral with a stabilizing characteristic.

2. Preliminary Results

In what follows, by R + n , we denote the open orthant
R + n = { x = ( x 1 , , x n ) R n , x 1 > 0 , , x n > 0 } .
and by R ¯ + n its closure:
R ¯ + n = { x = ( x 1 , , x n ) R n , x 1 0 , , x n 0 } .
By γ = ( γ 1 , . . . , γ n ) , we denote a multi-index that consists of some positive real numbers γ i , i = 1 , . . . , n , and | γ | = γ 1 + + γ n .
A part of a sphere of radius r with the center at the origin that belongs to R + n will be denoted by S r + ( n ) :
S r + ( n ) = { x R ¯ + n : | x | = r } { x R ¯ + n : x i = 0 , | x | r , i = 1 , . . . , n } .
For our derivations, we employ the following two formulas (see, e.g., the formulas 1.107 and 3.140 in [16]):
S 1 + ( n ) σ γ d S = i = 1 n Γ γ i + 1 2 2 n 1 Γ n + | γ | 2 , σ γ = i = 1 n σ i γ i ,
S 1 + ( n ) j γ ( x , r σ ) σ γ d S = i = 1 n Γ γ i + 1 2 2 n 1 Γ n + | γ | 2 j n + | γ | 2 1 ( r | x | ) ,
where
j γ ( x ; ξ ) = i = 1 n j γ i 1 2 ( x i ξ i )
and j ν stands for the normalized Bessel function of the first kind
j ν ( x ) = 2 ν Γ ( ν + 1 ) x ν J ν ( x ) ,
where J ν is the Bessel function of the first kind.
Let Ω be finite or infinite open set in R n symmetric with respect to each hyperplane x i = 0 , i = 1 , . . . , n , Ω + = Ω R + n , and Ω ¯ + = Ω R ¯ + n , where R ¯ + n = { x = ( x 1 , , x n ) R n , x 1 0 , , x n 0 } . In this paper, we employ the class C m ( Ω + ) of functions that are m-times differentiable on Ω + . By C m ( Ω ¯ + ) , we denote the subset of C m ( Ω + ) that contains the function with all continuous derivatives with respect to x i for any i = 1 , . . . , n up to x i = 0 . The class C e v m ( Ω ¯ + ) of functions consists of all functions from C m ( Ω ¯ + ) such that 2 k + 1 f x i 2 k + 1 | x i = 0 = 0 for all non-negative integers k m 1 2 (see [17], p. 21). In what follows, we denote the space C e v m ( R ¯ + n ) just by C e v m and set
C e v ( Ω ¯ + ) = C e v m ( Ω ¯ + ) ,
where the intersection is taken over all m N 0 and all C e v ( R ¯ + ) = C e v .
Our results in the next sections will be formulated for the functions from some spaces of generalized functions. Now we proceed with a description of these spaces. As a space of the basic (or test) functions, the following subspace of the space of rapidly decreasing functions will be used:
S e v = f C e v : sup x R + n x α D β f ( x ) < α , β Z + n ,
where α = ( α 1 , . . . , α n ) , β = ( β 1 , . . . , β n ) , α 1 , . . . , α n , β 1 , . . . , β n are integer non-negative numbers, x α = x 1 α 1 x 2 α 2 x n α n , D β = D x 1 β 1 . . . D x n β n , D x j = x j .
As usual, the space of the weighted generalized functions S e v ( R + n ) = S e v is defined as a class of continuous linear functionals that map the test functions from f S e v into the set of real numbers. Each function u such that u ( x ) x γ L 1 , l o c , x γ = x 1 γ 1 . . . x n γ n will be identified with the functional u S e v ( R + n ) = S e v that acts according to the formula
( u , f ) γ = R + n u ( x ) f ( x ) x γ d x , f S e v .
The generalized functions u S e v defined by Equation (11) with u ( x ) x γ L 1 , l o c will be called regular weighted generalized functions. All other generalized functions u S e v will be called singular weighted generalized functions.
By L p γ ( R + n ) = L p γ , 1 p < , we denote the space of all functions measurable in R + n even with respect to each variable x i , i = 1 , . . . , n such that
R + n | f ( x ) | p x γ d x < ,
where x γ = i = 1 n x i γ i . For a real number p 1 , the L p γ -norm of f is defined by
| | f | | L p γ ( R + n ) = | | f | | p , γ = R + n | f ( x ) | p x γ d x 1 / p .
It is known (see [17]) that L p γ is a Banach space.
The modified Bessel functions I α ( x ) and K α ( x ) (occasionally called the hyperbolic Bessel functions) of the first and second kind, respectively, are defined as follows (see [18,19]):
I α ( x ) = i α J α ( i x ) = m = 0 1 m ! Γ ( m + α + 1 ) x 2 2 m + α ,
K α ( x ) = π 2 I α ( x ) I α ( x ) sin ( α π ) ,
where α is non-integer. For integer values of the parameter α , the functions are defined by taking the limit in the expressions presented above. Obviously, the function K α ( x ) is an even function, i.e., the relation
K α ( x ) = K α ( x )
holds valid. For r 0 ( 0 < | r | ν + 1 ), the known asymptotic formula for the function K α ( x ) is as follows [20]:
K ν ( r ) ln r 2 ϑ , if ν = 0 , Γ ( ν ) 2 1 ν r ν , if ν > 0 ,
where
ϑ = lim n ln n + k = 1 n 1 k = 1 1 x + 1 x d x
is the Euler–Mascheroni constant.
For r + , we have another known asymptotic formula:
K ν ( r ) π 2 r e r .
Some of our results will be formulated in terms of the generalized translation defined as follows (see [21]):
γ T x y f ( x ) = C ( γ ) 0 π f ( x 2 + y 2 2 x y cos φ ) sin γ 1 φ d φ , γ > 0 ,
C ( γ ) = 0 π sin γ 1 φ d φ 1 = Γ γ + 1 2 π Γ γ 2 , γ > 0 .
For γ = 0 , the generalized translation γ T x y is reduced to the following simple translation:
0 T x y = T x y f ( x ) = f ( x + y ) + f ( x y ) 2 .
The multidimensional generalized translation is defined by the equality
( γ T x y f ) ( x ) = γ T x y f ( x ) = ( γ 1 T x 1 y 1 . . . γ n T x n y n f ) ( x ) ,
where each of the one-dimensional generalized translations γ i T x i y i , i = 1 , . . . , n acts according to Equation (17) for γ > 0 and according to Equation (18) for γ = 0 .
Another object useful for our derivations is a generalized convolution generated by a multidimensional generalized translation γ T x y . It is defined as follows:
( f g ) γ ( x ) = ( f g ) γ = R + n f ( y ) ( γ T x y g ) ( x ) y γ d y .
It is worth mentioning that the operators γ T x h can be expanded in terms of the powers of the Bessel operator by the following formula from [21] (compare to the well-known representation of the shift operator in terms of powers of the differentiation operator D):
γ T x h u ( x ) = k = 0 m φ k ( h ) B γ k u ( x ) + r m ( x , h ) ,
where
φ k ( h ) = 1 k ! Γ γ + 1 2 Γ γ + 1 2 + k h 2 2 k ,
r m ( x , h ) = φ m + 1 ( h ) B γ m + 1 u ( ξ ) | ξ = x + θ h , 1 < θ < 1 .
In particular, we obtain
φ 0 ( y ) = 1 , φ 1 ( h ) = Γ γ + 1 2 Γ γ + 1 2 + 1 h 2 2 = h 2 2 ( γ + 1 ) .
Equation (21) is called the Taylor–Delsarte formula. Letting m go to infinity in Equation (21), we arrive at the the Taylor–Delsarte series
γ T x h u ( x ) = k = 0 φ k ( h ) B γ k u ( x ) .
The Taylor–Delsarte Formula (21) can be easily generalized to the case of functions depending on several variables. In what follows, by
( P y l 1 f ) ( x ) = | 2 m | l 1 ζ m γ ( B x m f ) ( x ) y 2 m
we denote a segment of the multidimensional Taylor–Delsarte series, where m = ( m 1 , . . . , m n ) is a multi-index containing the non-negative integers, y 2 m = y 1 2 m 1 . . . y n 2 m n , ( B x m f ) ( x ) = B γ 1 m 1 . . . B γ n m n f ( x 1 , . . . , x n ) , ζ m γ = 1 2 2 | m | m ! i = 1 n Γ γ i + 1 2 Γ m i + γ i + 1 2 , and m ! = m 1 ! · . . . · m n ! .
In this paper, we also employ the multidimensional Hankel transform with the kernel defined by Equation (9) acting on the functions f L 1 γ ( R + n ) . It is provided by the formula
F γ [ f ] ( ξ ) = F γ [ f ( x ) ] ( ξ ) = f ^ ( ξ ) = R + n f ( x ) j γ ( x ; ξ ) x γ d x ,
where R + n is defined by Equation (6), the kernel j γ is given by Equation (9), and x γ = i = 1 n x i γ i . Let f L 1 γ ( R + ) be a function of bounded variation in a neighborhood of a point x of its continuity. Then, for γ > 0 , the following inversion formula for the multidimensional Hankel transform holds valid:
F γ 1 [ f ^ ( ξ ) ] ( x ) = f ( x ) = 2 n | γ | j = 1 n Γ 2 γ j + 1 2 R + n j γ ( x , ξ ) f ^ ( ξ ) ξ γ d ξ .
For the multidimensional Hankel transform and the following convolution operator with respect to the generalized convolution
A f = ( T f ) γ , f S e v ,
the convolution theorem
F γ [ A f ] = F γ [ ( T f ) γ ] = F γ [ T ] · F γ [ f ]
holds valid (see [16], p. 156, formula 3.176).

3. Some Useful Properties of the General Poisson Kernel

For our construction of an inverse operator to the generalized Bessel potential in Equation (1) by the method of approximative inverse operators, we employ the generalized convolution with the general Poisson kernel in the form
P γ ( x , ε ) = 2 n Γ n + | γ | + 1 2 π j = 1 n Γ γ j + 1 2 ε ( ε 2 + | x | 2 ) n + | γ | + 1 2 , ε > 0 .
In this section, we present some properties of this kernel needed for the further discussions and start with the following result:
Lemma 1.
The image of the function e ε | x | , ε > 0 by the Hankel transform defined by Equation (23) can be represented in the form
F γ [ e ε | x | ] ( ξ ) = 2 | γ | ε i = 1 n Γ γ i + 1 2 Γ n + | γ | + 1 2 π ( ε 2 + | ξ | 2 ) n + | γ | + 1 2 .
Proof. 
By definition and employing a simple variables substitution, we first get the representation
F γ [ e ε | x | ] ( ξ ) = R + n e ε | x | j γ ( x ; ξ ) x γ d x = { x = ρ σ } =
= 0 e ε ρ ρ n + | γ | 1 d ρ S 1 + ( n ) j γ ( ρ σ ; ξ ) σ γ d S .
Then, applying Equation (8), we obtain
F γ [ e ε | x | ] ( ξ ) = i = 1 n Γ γ i + 1 2 2 n 1 Γ n + | γ | 2 0 e ε ρ j n + | γ | 2 1 ( ρ | ξ | ) ρ n + | γ | 1 d ρ =
= i = 1 n Γ γ i + 1 2 2 n | γ | 2 | ξ | n + | γ | 2 1 0 e ε ρ J n + | γ | 2 1 ( ρ | ξ | ) ρ n + | γ | 2 d ρ .
Now, using the formula 2.12.8.4 from [22], p. 164 in the form
0 x ν + 2 e p x J ν ( c x ) d x = 2 p ( 2 c ) ν Γ ν + 3 2 π ( p 2 + c 2 ) ν + 3 2 , R e ν > 1 ,
we arrive at the representation
0 e ε ρ J n + | γ | 2 1 ( ρ | ξ | ) ρ n + | γ | 2 d ρ = 2 ε ( 2 | ξ | ) n + | γ | 2 1 Γ n + | γ | + 1 2 π ( ε 2 + | ξ | 2 ) n + | γ | + 1 2 .
Finally, we get the desired result:
F γ [ e ε | x | ] ( ξ ) = i = 1 n Γ γ i + 1 2 2 n | γ | 2 | ξ | n + | γ | 2 1 2 ε ( 2 | ξ | ) n + | γ | 2 1 Γ n + | γ | + 1 2 π ( ε 2 + | ξ | 2 ) n + | γ | + 1 2 =
= 2 | γ | ε i = 1 n Γ γ i + 1 2 Γ n + | γ | + 1 2 π ( ε 2 + | ξ | 2 ) n + | γ | + 1 2 .
Now we deduce some properties of the general Poisson kernel needed for the further discussions.
Lemma 2.
The general Poisson kernel P γ ( x , ε ) possesses the following properties:
1. 
F γ [ P γ ( x , ε ) ] ( ξ ) = e ε | ξ | ,
2. 
R + n P γ ( x , ε ) x γ d x = R + n P γ ( x , 1 ) x γ d x = 1 ,
3. 
P γ ( x , ε ) L p γ , 1 p .
Proof. 
1. This property immediately follows from Lemma 1.
2. By the variables substitution x = ε y , we first obtain
R + n P γ ( x , ε ) x γ d x = 2 n ε Γ n + | γ | + 1 2 π j = 1 n Γ γ j + 1 2 R + n x γ d x ( ε 2 + | x | 2 ) n + | γ | + 1 2 = { x = ε y } =
= 2 n Γ n + | γ | + 1 2 π j = 1 n Γ γ j + 1 2 R + n y γ d y ( 1 + | y | 2 ) n + | γ | + 1 2 = R + n P γ ( x , 1 ) x γ d x .
Now we show that R + n P γ ( x , 1 ) x γ d x = 1 . Indeed, first, we get the representation
R + n y γ d y ( 1 + | y | 2 ) n + | γ | + 1 2 = { y = ρ σ } = 0 ρ n + | γ | 1 d ρ ( 1 + ρ 2 ) n + | γ | + 1 2 S 1 + ( n ) σ γ d S .
To simplify the integral at the right-hand side of the last relation, we employ Equation (7) and obtain
R + n y γ d y ( 1 + | y | 2 ) n + | γ | + 1 2 = i = 1 n Γ γ i + 1 2 2 n 1 Γ n + | γ | 2 0 ρ n + | γ | 1 d ρ ( 1 + ρ 2 ) n + | γ | + 1 2 = { ρ 2 = r } =
= i = 1 n Γ γ i + 1 2 2 n Γ n + | γ | 2 0 r n + | γ | 2 1 ( 1 + r ) n + | γ | + 1 2 d r .
Now we use the formula 2.2.5.24 from [23], p. 239 in the form
0 x α 1 ( x + z ) β d x = z α β B ( α , β α ) , 0 < R e α < R e β ,
and obtain the relation
2 0 ρ n + | γ | 1 d ρ ( 1 + ρ 2 ) n + | γ | + 1 2 = 0 r n + | γ | 2 1 ( 1 + r ) n + | γ | + 1 2 d r = π Γ n + | γ | 2 Γ n + | γ | + 1 2 .
Thus, the formula
R + n y γ d y ( 1 + | y | 2 ) n + | γ | + 1 2 = i = 1 n Γ γ i + 1 2 2 n Γ n + | γ | 2 π Γ n + | γ | 2 Γ n + | γ | + 1 2 = π i = 1 n Γ γ i + 1 2 2 n Γ n + | γ | + 1 2
holds valid. Finally, we get the desired property of the general Poisson kernel:
R + n P γ ( x , 1 ) x γ d x = 2 n Γ n + | γ | + 1 2 π j = 1 n Γ γ j + 1 2 R + n y γ d y ( 1 + | y | 2 ) n + | γ | + 1 2 = 2 n Γ n + | γ | + 1 2 π j = 1 n Γ γ j + 1 2 π i = 1 n Γ γ i + 1 2 2 n Γ n + | γ | + 1 2 = 1 .
3. To show that P γ ( x , ε ) L p γ , 1 p , we start with the representation
R + n x γ d x ( ε 2 + | x | 2 ) p n + | γ | + 1 2 = ε ( n + | γ | ) ( 1 p ) p R + n x γ d x ( | x | 2 + 1 ) p n + | γ | + 1 2 =
= { x = ρ σ , | x | = ρ } = ε ( n + | γ | ) ( 1 p ) p 0 ρ n + | γ | 1 d ρ ( ρ 2 + 1 ) p n + | γ | + 1 2 S 1 + ( n ) σ γ d S .
For 1 p < , we apply Equations (7) and (27) and obtain the following norm estimate:
| | P γ ( x , ε ) | | p , γ = ε ( n + | γ | ) ( 1 p ) p π Γ n + | γ | 2 2 Γ n + | γ | + 1 2 i = 1 n Γ γ i + 1 2 2 n 1 Γ n + | γ | 2 1 p =
= ε ( n + | γ | ) ( 1 p ) p π i = 1 n Γ γ i + 1 2 2 n Γ n + | γ | + 1 2 1 p < .
In the case p = + , the inequality | | P γ ( x , ε ) | | , γ < is justified by passing to the limit p in the last formula and using the relation [16]
| | f | | , γ = lim p | | f | | p , γ , f L γ .

4. Inversion of the Bessel B-Potential by the Method of an Approximative Inverse Operator

In this section, we derive an inverse operator to the generalized Bessel potential by applying the method of an approximative inverse operator.
The approximative inverse operator (AIO) method presented in [24] is one of the possible ways for the derivation of an inversion to certain convolution operators. The idea of the AIO method is that the inverse operator is constructed as the limit of a certain sequence of the approximative convolution operators with integrable kernels.
To apply this method for inversion of the Bessel B-potential, we first introduce the function
ω α , γ ( | x | ) = 2 n | γ | α 2 + 1 Γ α 2 i = 1 n Γ γ i + 1 2 | x | n + | γ | α 2 K n + | γ | α 2 ( | x | ) = = 2 n α + 2 Γ α 2 i = 1 n Γ γ i + 1 2 0 t n + | γ | α 2 1 e t | x | 2 4 t d t .
Then the generalized Bessel potential in Equation (1) can be represented as a convolution operator with respect to the generalized convolution defined by Equation (20):
G γ α φ = ω α , γ ( | x | ) | x | n + | γ | α φ γ , α > 0 .
Now we consider the convolution operators
( G γ , ε α ) 1 φ = ( g γ , ε α φ ) γ
with the kernel
g γ , ε α ( x ) = F γ 1 ( 1 + | ξ | 2 ) α / 2 · e ε | ξ | ( x ) .
To obtain a suitable integral representation of the kernel g γ , ε α , we change to the spherical coordinates ξ = r σ , σ R + n in the integral at the right-hand side of the previous relation and then apply Equation (8):
g γ , ε α ( x ) = 2 n | γ | j = 1 n Γ 2 γ j + 1 2 R + n ( 1 + | ξ | 2 ) α / 2 · e ε | ξ | j γ ( x , ξ ) ξ γ d ξ =
= 2 n | γ | j = 1 n Γ 2 γ j + 1 2 0 ( 1 + r 2 ) α / 2 e ε r r n + | γ | 1 d r S 1 + ( n ) j γ ( x , r σ ) σ γ d S =
= 2 1 | γ | j = 1 n Γ γ j + 1 2 Γ n + | γ | 2 0 j n + | γ | 2 1 ( r | x | ) ( 1 + r 2 ) α / 2 e ε r r n + | γ | 1 d r .
Now we introduce the notion of a B-multiplier in L p γ . Let M S e v . The weighted generalized function M is called a B-multiplier in L p γ if for f S e v , the generalized convolution ( F γ 1 M f ) γ belongs to the space L p γ and the inequality
sup | | f | | p , γ = 1 | | ( F γ 1 M f ) γ | | p , γ < +
holds true. The linear space of all B-multipliers is denoted by M p , γ = M p , γ ( R + n ) . On the space M p , γ , the left-hand side of the inequality in Equation (30) defines a norm that is used in the further discussion.
Lemma 3.
Let f S e v . The operator
( G γ , ε α ) 1 φ = ( g γ , ε α φ ) γ = R + n g γ , ε α ( t ) ( γ T x t φ ( x ) ) t γ d t
is bounded in L p γ , 1 < p < + .
Proof. 
Let ε > 0 be a fixed number. The function ( 1 + | ξ | 2 ) α / 2 · e ε | ξ | is a B-multiplier for any p , 1 < p < + . Since the kernel function g γ , ε α ( x ) = F γ 1 ( 1 + | ξ | 2 ) α / 2 · e ε | ξ | ( x ) is represented as an operator generated by the B-multiplier ( 1 + | ξ | 2 ) α / 2 · e ε | ξ | in L p γ , the inclusion g γ , ε α ( x ) L p γ holds valid.
Because the operator
( G γ , ε α ) 1 φ ( x ) = F γ 1 [ ( 1 + | ξ | 2 ) α / 2 · e ε | ξ | ] ( x ) φ ( x ) γ
is a generalized convolution with a B-multiplier, its image also belongs to L p γ . □
Lemma 4.
Let φ S e v . Then the formula
( ( G γ α ) ε 1 G γ α φ ) ( x ) = ( P γ , ε f ) ( x ) ,
holds valid, where ( P γ , ε f ) ( x ) is the generalized convolution with the Poisson kernel
( P γ , ε f ) ( x ) = ( f ( x ) P γ ( x , ε ) ) γ .
Proof. 
By the convolution theorem, we get the relation
F γ [ G γ α φ ] ( ξ ) = ( 1 + | ξ | 2 ) α / 2 F γ φ .
Let us denote G γ α φ by ψ . Then we have
F γ ( ( G γ α ) ε 1 ψ ) ( x ) = F γ F γ 1 [ ( 1 + | ξ | 2 ) α / 2 · e ε | ξ | ] ( x ) ψ ( x ) γ
= ( 1 + | ξ | 2 ) α / 2 · e ε | ξ | ( x ) · F γ ψ = ( 1 + | ξ | 2 ) α / 2 · e ε | ξ | ( x ) · ( 1 + | ξ | 2 ) α / 2 F γ φ = e ε | ξ | · F γ φ
that leads to the representation
( ( G γ α ) ε 1 G γ α φ ) ( x ) = F γ 1 e ε | ξ | · F γ φ ( x ) = F γ 1 ( e ε | ξ | ) ( x ) F γ φ γ .
Applying Equation (26), we arrive at the relation
( F γ 1 e ε | ξ | ) ( x ) = 2 n | γ | j = 1 n Γ 2 γ j + 1 2 2 | γ | ε i = 1 n Γ γ i + 1 2 Γ n + | γ | + 1 2 π ( ε 2 + | x | 2 ) n + | γ | + 1 2 = = 2 n Γ n + | γ | + 1 2 π j = 1 n Γ γ j + 1 2 ε ( ε 2 + | x | 2 ) n + | γ | + 1 2 = P γ ( x , ε ) ,
where P γ ( x , ε ) is the generalized Poisson kernel in Equation (25). By Lemma 2, the inclusion P γ ( x , ε ) L p γ holds valid. □
Theorem 1.
Let φ L p γ , 1 p + and
( G γ α ) 1 φ ( x ) = lim ε 0 L p γ ( G γ , ε α ) 1 φ ( x ) = F γ 1 [ ( 1 + | ξ | 2 ) α / 2 · e ε | ξ | ] ( x ) φ ( x ) γ .
Then the relation
( G γ α ) 1 G γ α φ ( x ) = φ ( x )
holds valid, where the limit is understood in the sense of the norm in L p γ .
Proof. 
According to Lemma 4, we have to prove that
| | ( P γ , ε f ) ( x ) f ( x ) | | p , γ 0 as ε 0 .
Taking into account the 2nd property from Lemma 2, we get the representation
( f ( x ) P γ ( x , ε ) ) γ f ( x ) = R + n [ γ T x y f ( x ) f ( y ) ] P γ ( y , ε ) y γ d y .
Applying the generalized Minkowski inequality [16] to the integral at the right-hand side of the last formula, we have
| | ( f ( x ) P γ ( x , ε ) ) γ f ( x ) | | p , γ R + n R + n [ γ T x y f ( x ) f ( x ) ] p x γ d x 1 p | P γ ( y , ε ) | y γ d y = { y = ε t }
= R + n R + n [ γ T x ε t f ( x ) f ( x ) ] p x γ d x 1 p | P γ ( t , 1 ) | t γ d t .
Lemma 3.6 on p. 166 from [25] ensures that
| | γ T x ε t f ( x ) f ( x ) | | p , γ c | | f ( x ) | | p , γ
for f L p γ . The Proposition 4.1, p. 182 from [26] and p. 50 from [27], leads to the formula
lim ε 0 R + n [ γ T x ε t f ( x ) f ( x ) ] p x γ d x 1 p = 0 .
By the Lebesgue majorized convergence theorem, the integral in Equation (33) tends to zero as ε 0 since the integrand is majorized by the integrable function c | | f | | p , γ | P γ ( t , 1 ) | t γ . □

5. Inversion of the Bessel B-Potential by the Hadamard Finite Part Regularization Technique

In this section, we construct another form of an inverse operator to the Bessel B-potential using the concept of the “finite part” of a singular integral introduced by Hadamard.
Let a function f = f ( x ) be integrable in a layer ε < | x | < A for any 0 < ε , ε < A < and the representation
ε < | x | < A f ( x ) d x = k = 1 N a k ε λ k + b ln 1 ε + J ε
hold valid, where a k , b, and λ k are some constant positive numbers independent of A. If the limit lim ε 0 J ε exists, then it is called the Hadamard finite part of the singular integral of the function f. The function f = f ( x ) is said to possess the Hadamard property at the origin. The standard notation for the finite part of the Hadamard singular integral is as follows
p . f . | x | < A f ( x ) d x = lim ε 0 J ε .
In the case b = 0 in Equation (34), the function f = f ( x ) is said to possess the non-logarithmic-type Hadamard property at the origin.
If a function f = f ( x ) is integrable on the set R + n \ { | x | > ε } and possesses the Hadamard property at the origin, then its finite part is defined as follows:
p . f . R + n f ( x ) d x = p . f . | x | < A f ( x ) d x + | x | > A f ( x ) d x .
We start by mentioning that the asymptotic properties of the modified Bessel function K ν ensure that the kernel function ω α , γ ( | x | ) exponentially decays at infinity (see Equation (16)) and goes to a constant at the origin (see Equation (15)):
lim | x | ω α , γ ( | x | ) = 2 n | γ | + α 2 + 1 Γ α 2 i = 1 n Γ γ i + 1 2 lim | x | | x | n + | γ | + α 2 K n + | γ | + α 2 ( | x | ) =
= 2 n | γ | + α 2 + 1 Γ α 2 i = 1 n Γ γ i + 1 2 π 2 lim | x | | x | n + | γ | + α 2 1 2 e | x | = 0 ,
lim | x | 0 ω α , γ ( | x | ) = 2 n | γ | + α 2 + 1 Γ α 2 i = 1 n Γ γ i + 1 2 lim | x | 0 | x | n + | γ | + α 2 K n + | γ | + α 2 ( | x | ) =
= 2 n | γ | + α 2 + 1 Γ α 2 i = 1 n Γ γ i + 1 2 lim | x | 0 | x | n + | γ | + α 2 Γ n + | γ | + α 2 2 1 n + | γ | + α 2 | x | n + | γ | + α 2 = 2 n + α Γ n + | γ | + α 2 Γ α 2 i = 1 n Γ γ i + 1 2 .
Therefore,
lim | x | ω α , γ ( | x | ) = 0 , lim | x | 0 ω α , γ ( | x | ) = 2 n + α Γ n + | γ | + α 2 Γ α 2 i = 1 n Γ γ i + 1 2 .
Let as note that the function Γ α 2 in the denominator of ω α , γ ( 0 ) has the simple poles at the points α = 0 , 2 , 4 , 6 , . . . , and therefore, the weighted generalized function 1 | x | n + | γ | α has the simple poles at the same points.
As already mentioned in the previous section, the operator G γ α can be represented in terms of the generalized convolution, as in Equation (20)
G γ α φ = ω α , γ ( | x | ) | x | n + | γ | α φ γ , α > 0 .
Then we can apply the convolution theorem for the multidimensional Hankel transform (Equation (24)) and obtain the relation
F γ [ G γ α φ ] = F γ ω α , γ ( | x | ) | x | n + | γ | α · F γ [ φ ] .
Furthermore, the inclusion ω α , γ ( | x | ) | x | n + | γ | α L γ 1 ( R + n ) and the formula
F γ ω α , γ ( | x | ) | x | n + | γ | α ( ξ ) = ( 1 + | ξ | 2 ) α / 2
hold valid, and we arrive at the representation
F γ [ G γ α φ ] = ( 1 + | ξ | 2 ) α / 2 · F γ [ φ ] .
Because of the semi-group property for the generalized Bessel potential [15], its inverse operator can be constructed in the form of a regularization of the generalized convolution
( G γ α ) 1 φ = ω α , γ ( | x | ) | x | n + | γ | + α φ γ , α > 0 .
The inverse operator ( G γ α ) 1 can be formally represented in terms of the multidimensional Hankel transform as follows:
F γ ( G γ α ) 1 φ = F γ ω α , γ ( | x | ) | x | n + | γ | + α · F γ [ φ ] , α > 0 .
However, Equation (37) shows that the integrand of the integral in Equation (39) has a singularity at the origin with the order greater than the dimension n of R + n , and therefore, the integral diverges. Thus, this integral cannot be used in its present form and should be suitably regularized. In what follows, we present some results regarding the regularization of this integral by employing the method of the Hadamard finite part of the singular integrals and construct the inverse operator ( G γ α ) 1 in the following form
( G γ α ) 1 φ = p . f . ω α , γ ( | x | ) | x | n + | γ | + α f γ .
For further derivations, some advanced results regarding asymptotic behavior of the modified Bessel function K ν are needed. It is known ([28], p. 273) that the inclusion r ν K ν ( r ) C ( 0 , ) , ν R holds valid. Moreover, we have the estimates
d m d r m r ν K ν ( r ) c · r ν 1 2 e r
for large values of r and the estimates
d m d r m r ν K ν ( r ) c r θ i f   m < 2 ν , ν > 0 ; r 2 ν m i f   m > 2 ν ,
for small values of r, where θ = 1 if m is odd, and θ = 0 if m is even, and ν is not an integer. If ν Z , then
d m d r m r ν K ν ( r ) c 1 i f   m < 2 ν , ν > 0 ; r 2 ν m ln 1 r i f   m > 2 ν .
For m < 2 ν , the formula
d m d r m r ν K ν ( r ) = ( 1 ) m 2 2 ν m 2 1 ( m 1 ) ! ! Γ ν m 2 , m = 2 , 4 , 6 , . . . < 2 ν
holds valid.
Equations (40)–(42) lead to the following statement:
Lemma 5.
For m < n + | γ | + α , the inclusion ω α , γ ( r ) C m ( [ 0 , ] ) holds valid as well as the formulas
ω α , γ ( m ) ( 0 ) = 0 , m = 1 , 3 , 5 , . . . < n + | γ | + α ,
-4.6cm0cm
ω α , γ ( m ) ( 0 ) = ( 1 ) m 2 2 n + α m 2 ( m 1 ) ! ! Γ n + | γ | + α 2 m 2 Γ α 2 i = 1 n Γ γ i + 1 2 , m = 2 , 4 , 6 , . . . < n + | γ | + α .
In the next theorem, a regularization of the generalized convolution ω α , γ ( | x | ) | x | n + | γ | + α φ γ is introduced using the appropriate segment of the Taylor–Delsarte series of the function φ .
Theorem 2.
Let φ ( x ) S e v ( R + n ) and α > 0 . Then the representation
p . f . ω α , γ ( | x | ) | x | n + | γ | + α φ γ = R + n γ T x y φ ( x ) ( P y l 1 φ ) ( x ) | y | n + | γ | + α ω α , γ ( | y | ) y γ d y + + | 2 m | l 1 ζ m γ ( B x m φ ) ( x ) 2 2 | m | i = 1 n Γ γ i + 1 2 + m i Γ | m | α 2 Γ α 2 i = 1 n Γ γ i + 1 2
holds valid, where α < l N , ( P y l 1 φ ) ( x ) = | 2 m | l 1 ζ m γ ( B x m φ ) ( x ) y 2 m is a segment of the Taylor–Delsarte series, m = ( m 1 , . . . , m n ) is the multi-index of non-negative integers, y 2 m = y 1 2 m 1 . . . y n 2 m n , ( B x m φ ) ( x ) = B γ 1 m 1 . . . B γ n m n φ ( x 1 , . . . , x n ) , ζ m γ = 1 2 2 | m | m ! i = 1 n Γ γ i + 1 2 Γ m i + γ i + 1 2 , and m ! = m 1 ! · . . . · m n ! .
Proof. 
We start with the representation
p . f . ω α , γ ( | x | ) | x | n + | γ | + α φ γ = p . f . R + n ω α , γ ( | y | ) γ T x y φ ( x ) | y | n + | γ | + α y γ d y =
= R + n γ T x y φ ( x ) ( P y l 1 φ ) ( x ) | y | n + | γ | + α ω α , γ ( | y | ) y γ d y + p . f . R + n ( P y l 1 φ ) ( x ) | y | n + | γ | + α ω α , γ ( | y | ) y γ d y ,
where
( P y l 1 φ ) ( x ) = | 2 m | l 1 ζ m γ ( B x m φ ) ( x ) y 2 m
is a segment of the Taylor–Delsarte series, l > α . Then we have the relation
p . f . R + n ( P y l 1 φ ) ( x ) | y | n + | γ | + α ω α , γ ( | y | ) y γ d y =
= | 2 m | l 1 ζ m γ ( B x m φ ) ( x ) p . f . R + n y 2 m | y | n + | γ | + α ω α , γ ( | y | ) y γ d y .
In order to derive a suitable representation for p . f . R + n y 2 m | y | n + | γ | + α ω α , γ ( | y | ) y γ d y , let us consider an integral with an addition parameter β
I α , β , γ = R + n y 2 m | y | n + | γ | + β ω α , γ ( | y | ) y γ d y , β C ,
and suppose that α > 0 is fixed.
Passing to the spherical coordinates in the integral I α , β , γ , we first obtain the representation
I α , β , γ = R + n y 2 m | y | n + | γ | + β ω α , γ ( | y | ) y γ d y = { y = r σ } = 0 r 2 | m | β 1 ω α , γ ( r ) d r S 1 + ( n ) σ γ + 2 m d S .
Using the formula 1.107, p. 49 from [16], we then get
S 1 + ( n ) σ 2 m + γ d S = i = 1 n Γ γ i + 1 2 + m i 2 n 1 Γ n + | γ | 2 + | m | = I m γ
and
I α , β , γ = I m γ 0 r 2 | m | β 1 ω α , γ ( r ) d r =
= I m γ 2 n | γ | + α 2 + 1 Γ α 2 i = 1 n Γ γ i + 1 2 0 r 2 | m | β 1 + n + | γ | + α 2 K n + | γ | + α 2 ( r ) d r .
The next step is applying the formula 2.16.2 from [22] in the form
0 x a 1 K ν ( c x ) d x = 2 a 2 c a Γ a + ν 2 Γ a ν 2 , ( a ) > | ( ν ) | , c > 0 .
In our case, the parameters are as follows: a = 2 | m | β + n + | γ | + α 2 , ν = n + | γ | + α 2 > 0 , and the convergence condition takes the form of 2 | m | ( β ) + n + | γ | + α 2 > n + | γ | + α 2 or ( β ) < 0 since | m | = 0 , 1 , 2 , . . . . Thus, we can calculate the integral I α , β , γ under the condition ( β ) < 0 :
I α , β , γ = I m γ 2 n + 2 | m | + α β 1 Γ α 2 i = 1 n Γ γ i + 1 2 Γ | m | + n + | γ | + α β 2 Γ | m | β 2 .
The formula above can be analytically extended with respect to β to the half-plane ( β ) > 0 . Setting β = α in this extended formula, we arrive at the representation
p . f . R + n y 2 m | y | n + | γ | + α ω α , γ ( | y | ) y γ d y = I m γ 2 n + 2 | m | 1 Γ | m | + n + | γ | 2 Γ | m | α 2 Γ α 2 i = 1 n Γ γ i + 1 2 .
Substituting the value of the integral I m γ into the formula
p . f . ω α , γ ( | x | ) | x | n + | γ | + α φ γ = R + n γ T x y φ ( x ) ( P y l 1 φ ) ( x ) | y | n + | γ | + α ω α , γ ( | y | ) y γ d y +
+ | 2 m | l 1 ζ m γ ( B x m φ ) ( x ) I m γ 2 n + 2 | m | 1 Γ | m | + n + | γ | 2 Γ | m | α 2 Γ α 2 i = 1 n Γ γ i + 1 2
we arrive at the final result in Equation (45). □
Now we proceed with a construction of an inversion of the Bessel B-potential by the Hadamard finite part regularization technique. For φ S e v ( R + n ) and α > 0 , we define the operator
( I Δ γ ) α 2 φ = R + n γ T x y φ ( x ) ( P y l 1 φ ) ( x ) | y | n + | γ | + α ω α , γ ( | y | ) y γ d y + + | 2 m | l 1 ζ m γ ( B x m φ ) ( x ) 2 2 | m | i = 1 n Γ γ i + 1 2 + m i Γ | m | α 2 Γ α 2 i = 1 n Γ γ i + 1 2
that turns to be inverse to the Bessel B-potential G γ α φ .
Theorem 3.
On the space S e v , the operator ( I Δ γ ) α 2 , ( α ) > 0 defined by Equation (46) is the right- and left-inverse to the Bessel B-potential, i.e., the formulas
( I Δ γ ) α 2 G γ α φ = G γ α ( I Δ γ ) α 2 φ = φ
hold valid for any φ S e v .
Proof. 
The main idea of the proof is first to consider the generalized convolution of the weighted functional F γ 1 ( 1 + | x | 2 ) α / 2 with a function f from S e v . Then an analytic continuation of the weighted functional F γ 1 ( 1 + | x | 2 ) α / 2 with respect to α is constructed. It turns out that the generalized convolution of this analytic continuation with a function f from S e v coincides with the representation in Equation (46). The final part of the proof involves the application of the inverse multidimensional Hankel transform.
Now we consider the steps mentioned above in detail and start by considering the inverse Hankel integral transform F γ 1 ( 1 + | x | 2 ) α / 2 in the sense of the weighted distributions from the space S e v . On the space S e v , the operator inverse to G γ α has the form
( G γ α ) 1 φ =   ( F γ 1 ( 1 + | x | 2 ) α / 2 φ ) γ = ( F γ 1 ( 1 + | x | 2 ) α / 2 , γ T x y φ ) γ , φ S e v .
For φ S e v , the generalized convolution from the right-hand side of the last formula is analytical with respect to the parameter α on the whole complex plane.
Let us introduce a parameter β C that is not an integer number. It is known (see [21]) that the inclusion f S e v holds valid for f = γ T x y φ and φ S e v . Thus, for ( n + | γ | ) < ( β ) < l , the weighted functional ω β , γ ( | x | ) | x | n + | γ | + β , f γ , f S e v can be represented in terms of an appropriate segment of the Taylor–Delsarte series:
ω β , γ ( | x | ) | x | n + | γ | + β , f γ =
= R + n f ( y ) ( P y l 1 f ) ( 0 ) | y | n + | γ | + β ω β , γ ( | y | ) y γ d y + | 2 m | l 1 ζ m γ ( B x m f ) ( 0 ) I m γ 2 n + 2 | m | 1 Γ | m | + n + | γ | 2 Γ | m | β 2 Γ β 2 i = 1 n Γ γ i + 1 2 ,
where ( P y l 1 f ) ( 0 ) = | 2 m | l 1 ζ m γ ( B x m f ) ( 0 ) y 2 m . The integral at the right-hand side of Equation (48) is analytical with respect to the parameter β in the domain ( n + | γ | ) < ( β ) < l since the expression ( f ( y ) ( P y l 1 f ) ( 0 ) ) compensates the singularity of the denominator at the origin as in the case of the B-hyper-singular integral with a homogeneous characteristic [29]. Therefore, the weighed functional in Equation (48) is analytic with respect to β in the domain ( n + | γ | ) < ( β ) < l . Thus, under the condition ( n + | γ | ) < ( β ) < 0 , we arrive at the representation
ω β , γ ( | x | ) | x | n + | γ | + β , f γ = R + n ω β , γ ( | y | ) | y | n + | γ | + β f ( y ) y γ d y = ( F γ 1 ( 1 + | x | 2 ) β / 2 , f ) γ , f S e v .
In its turn, the weighed functional on the right-hand side of Equation (49) is analytic with respect to β in the whole complex plane because it can be represented in the form
( F γ 1 ( 1 + | x | 2 ) β / 2 , f ) γ = 2 n | γ | j = 1 n Γ 2 γ j + 1 2 ( ( 1 + | x | 2 ) β / 2 , F γ f ) γ
and because all weighted functions ( 1 + | x | 2 ) β / 2 are analytic with respect to β (see Chapter 4 in [16]). Therefore, the representation in Equation (49) also keeps its validity in the case β = α , 0 < ( α ) < l . Using the relation f = γ T x y φ , we rewrite the representation from Equation (48) for the function φ as follows:
( G γ α ) 1 φ = ( F γ 1 ( 1 + | x | 2 ) α / 2 φ ) γ = ( I Δ γ ) α 2 φ ,
where ( I Δ γ ) α 2 is defined as in Equation (46). Equations (47) now follow from their representations as the images of the Hankel integral transform. To derive the first formula, we first get the representation
F γ ( I Δ γ ) α 2 G γ α φ = F γ ( F γ 1 ( 1 + | x | 2 ) α / 2 G γ α φ ) γ = ( 1 + | x | 2 ) α / 2 · F γ G γ α φ = ( 1 + | x | 2 ) α / 2 · ( 1 + | ξ | 2 ) α / 2 F γ φ = F γ φ .
Applying the inverse Hankel transform to the last formula, we immediately arrive at the desired relation ( I Δ γ ) α 2 G γ α φ = φ . The second formula from Equation (47), i.e., G γ α ( I Δ γ ) α 2 φ = φ is proved by employing the same arguments. □
Finally, let us consider a particular case of the inversion of the Bessel B-potential in the form of Equation (46) under the condition 0 < α < 1 . In this case, for l = 1 , we have the formula ( P y 0 f ) ( x ) = f ( x ) and the inverse operator in Equation (46) takes a simple form
( I Δ γ ) α 2 φ = R + n γ T x y f ( x ) f ( x ) | y | n + | γ | + α ω α , γ ( | y | ) y γ d y + f ( x ) .

6. Conclusions

In this paper, we considered a convolution-type operator called the generalized Bessel potential. This is a very important object in the potential theory with several applications, including analytical treatment of the non-homogeneous equations of potential type, theory of functional spaces, and functional completion. Thus, a problem of the derivation of its inversion is worth investigation. In this paper, we presented two different forms of the operators inverse to the generalized Bessel potential. The first form was derived in terms of an approximative inverse operator using the method of an improving multiplier. The second form employs the regularization technique for the divergent integrals in terms of the appropriate segments of the Taylor–Delsarte series. It is worth mentioning that both inversions can be interpreted as new fractional differentiation operators of arbitrary positive order. Further research problems in this direction include detailed investigation of these new operators on different spaces of functions as well as the analysis of the partial differential equations of fractional order with the fractional derivatives in form of the inversions of the generalized Bessel potential.

Author Contributions

Conceptualization: Y.L. and E.S.; methodology: Y.L. and E.S.; formal analysis: A.D. and Y.L.; investigation: Y.L., A.D., and E.S.; resources: Y.L. and E.S.; writing—original draft preparation: Y.L., A.D., and E.S.; writing—review and editing: Y.L.; visualization: Y.L., A.D., and E.S.; supervision: Y.L.; project administration: Y.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

G α the generalized Bessel potentialEquation (1)
G α γ the generalized Bessel potential kernelEquation (2)
Δ γ the Laplace–Bessel operatorEquation (3)
B γ the Bessel operatorEquation (4)
j γ the multidimensional normalized Bessel function of the first kindEquation (9)
j ν the normalized Bessel function of the first kindEquation (10)
γ T x y the one-dimensional generalised translationEquation (17)
γ T x y the multidimensional generalized translationEquation (19)
( · · ) γ the generalized convolutionEquation (20)
F γ the multidimensional Hankel transformEquation (23)
P γ the general Poisson kernelEquation (25)
g γ , ε α the regularisation kernelEquation (29)
P γ , ε the generalized convolution with the Poisson kernelEquation (31)
p . f . the finite part of the Hadamard singular integralEquation (35)

References

  1. Aronszajn, N.; Smith, K.T. Theory of Bessel potentials. I. Ibid. 1961, 11, 365–475. [Google Scholar] [CrossRef]
  2. Calderon, A.P. Lebesgue spaces of differentiable functions and distributions. In Proceedings of the Symposia in Pure Mathematics 4, Berkeley, CA, USA, 21–22 April 1960; Volume 4, pp. 33–50. [Google Scholar]
  3. Stein, E.M. The characterization of functions arising as potentials. I. Bull. Amer. Math. Soc. 1961, 67, 102–104. [Google Scholar] [CrossRef] [Green Version]
  4. Lizorkin, P.I. Characterization of the spaces Lpr(Rn) in terms of difference singular integrals. Russian Mat. Sb. 1970, 81, 79–91. [Google Scholar]
  5. Nogin, V.A. Inversion of Bessel potentials. Differ. Equat. 1982, 18, 1407–1411. [Google Scholar]
  6. Nogin, V.A. Inversion of Bessel potentials by means of hypersingular integrals. Soviet Math. 1985, 29, 73–83. [Google Scholar]
  7. Guliev, V.S.; Safarov, Z.V. Bk,n-Bessel potentials and certain imbedding theorems in Bk,n-Sobolev–Liouville spaces. Proc. Inst. Math. Mech. Natl. Acad. Sci. Azerb. 2001, 15, 68–80. [Google Scholar]
  8. Guliev, V.S.; Serbetci, A.; Akbulut, A.; Mammadov, Y.Y. Nikolskii–Besov and Lizorkin–Triebel spaces constructed on the base of the multidimensional Fourier–Bessel transform. Eurasian Math. J. 2011, 2, 42–66. [Google Scholar]
  9. Gol’dman, M.L. Integral properties of generalized Bessel potentials. DAN 2007, 414, 159–164. [Google Scholar] [CrossRef]
  10. Gol’dman, M.L. Rearrangement-invriant spans for generalized Bessel and Riesz Potentials. DAN 2008, 423, 14–18. [Google Scholar]
  11. Gol’dman, M.L. The cone of rearrangements for generalized Bessel Potentials. Tr. Mat. Inst. Steklova 2008, 260, 151–163. [Google Scholar] [CrossRef]
  12. Aliev, I.A.; Eryigit, M. Inversion of Bessel potentials with the aid of weighted wavelet transforms. Math. Nachr. 2002, 242, 27–37. [Google Scholar] [CrossRef]
  13. Aliev, I.A.; Eryigit, M. On a rate of convergence of truncated hypersingular integrals associated to Riesz and Bessel potentials. J. Math. Anal. Appl. 2013, 406, 352–359. [Google Scholar] [CrossRef]
  14. Aliev, I.A.; Gadjiev, A.D.; Aral, A. On approximation properties of a family of linear operators at critical value of parameter. J. Approx. Theory 2006, 138, 242–253. [Google Scholar] [CrossRef] [Green Version]
  15. Shishkina, E.; Ekincioǧlu, I.; Keskin, C. Generalized Bessel potential and its application to non-homogeneous singular screened Poisson equation. Integral Transform. Spec. Funct. 2020, 1–16. [Google Scholar] [CrossRef]
  16. Shishkina, E.L.; Sitnik, S.M. Transmutations, Singular and Fractional Differential Equations with Applications to Mathematical Physics; Elsevier: Amsterdam, The Netherlands, 2020. [Google Scholar]
  17. Kipriyanov, I.A. Singular Elliptic Boundary Value Problems; Nauka: Moscow, Russia, 1997. [Google Scholar]
  18. Erdelyi, A.; Magnus, W.; Oberhettinger, F.; Tricomi, F.G. Higher Transcendental Functions; McGraw-Hill Book Co.: New York, NY, USA, 1953; Volume 2. [Google Scholar]
  19. Watson, G.N. A Treatise on the Theory of Bessel Functions; University Press: Cambridge, UK, 1922. [Google Scholar]
  20. Abramowitz, M.; Stegun, I.A. Handbook of Mathematical Functions with Formulas, Graphs and Mathematical Tables; Dover Publ., Inc.: New York, NY, USA, 1972. [Google Scholar]
  21. Levitan, B.M. Expansion in Fourier Series and Integrals with Bessel Functions. Uspekhi Mat. Nauk 1951, 6, 102–143. [Google Scholar]
  22. Prudnikov, A.P.; Brychkov, Y.A.; Marichev, O.I. Integrals and Series, Vol. 2, Special Functions; Gordon & Breach Sci. Publ.: New York, NY, USA, 1990. [Google Scholar]
  23. Prudnikov, A.P.; Brychkov, Y.A.; Marichev, O.I. Integrals and Series, Vol. 1, Elementary Functions; Gordon & Breach Sci. Publ.: New York, NY, USA, 1992. [Google Scholar]
  24. Samko, S.G. A new approach to the inversion of the Riesz potential operator. Fract. Calc. Appl. Anal. 1998, 1, 225–245. [Google Scholar]
  25. Platonov, S.S. Bessel generalized translations and some problems of approximation theory for functions on the half-line. Sibirsk. Mat. Zh. 2009, 50, 154–174. [Google Scholar] [CrossRef]
  26. Platonov, S.S. Bessel harmonic analysis and approximation of functions on a semiaxis. Izv. RAN. Ser. Mat. 2007, 71, 149–196. [Google Scholar]
  27. Platonov, S.S. Generalized Bessel translations and some inverse theorems of the theory of approximation of functions on the half-line. Tr. PetrGU. Ser. Mat. 2007, 14, 44–57. [Google Scholar]
  28. Samko, S.G. Hypersingular Integrals and Their Applications; Taylor & Francis: London, UK, 2002. [Google Scholar]
  29. Lyakhov, L.N.; Shishkina, E.L. General B-Hypersingular Integrals with Homogeneous Characteristic. DAN 2007, 75, 39–43. [Google Scholar] [CrossRef]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Dzhabrailov, A.; Luchko, Y.; Shishkina, E. Two Forms of An Inverse Operator to the Generalized Bessel Potential. Axioms 2021, 10, 232. https://doi.org/10.3390/axioms10030232

AMA Style

Dzhabrailov A, Luchko Y, Shishkina E. Two Forms of An Inverse Operator to the Generalized Bessel Potential. Axioms. 2021; 10(3):232. https://doi.org/10.3390/axioms10030232

Chicago/Turabian Style

Dzhabrailov, Akhmed, Yuri Luchko, and Elina Shishkina. 2021. "Two Forms of An Inverse Operator to the Generalized Bessel Potential" Axioms 10, no. 3: 232. https://doi.org/10.3390/axioms10030232

APA Style

Dzhabrailov, A., Luchko, Y., & Shishkina, E. (2021). Two Forms of An Inverse Operator to the Generalized Bessel Potential. Axioms, 10(3), 232. https://doi.org/10.3390/axioms10030232

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