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
where
and
stands for the modified Bessel functions of the second kind. The operator
can be interpreted as a realization of the real negative powers of the operator
. 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
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
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
–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
was derived in [
12] using a special type of the weighted wavelet transforms. Further investigation of the hypersingular integrals associated to
was performed in [
13]. In [
14], the fractional powers
of the operator
,
were introduced, and their approximation properties were studied as
tends to
.
In this paper, we first construct the positive fractional powers of the operator
, where
is the Laplace–Bessel operator with
. The negative fractional powers of this operator are the generalized Bessel potentials considered in [
15]. The most important difference between the generalized Riesz potential
,
, and the generalized Bessel potential
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]
with the kernel
and the generalized translation
defined by the Formula (
17) (see the next section). Equation (
1) realizes the negative fractional power
,
, where
I is the identity operator, and
is the Laplace–Bessel operator
and
is the Bessel operator. The results derived in [
15] allow us to represent the generalized Bessel potential in the form
where
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
,
. The operator
,
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
, we denote the open orthant
and by
its closure:
By
, we denote a multi-index that consists of some positive real numbers
,
and
.
A part of a sphere of radius
r with the center at the origin that belongs to
will be denoted by
:
For our derivations, we employ the following two formulas (see, e.g., the formulas 1.107 and 3.140 in [
16]):
where
and
stands for the normalized Bessel function of the first kind
where
is the Bessel function of the first kind.
Let
be finite or infinite open set in
symmetric with respect to each hyperplane
,
,
, and
, where
In this paper, we employ the class
of functions that are
m-times differentiable on
. By
, we denote the subset of
that contains the function with all continuous derivatives with respect to
for any
up to
. The class
of functions consists of all functions from
such that
for all non-negative integers
(see [
17], p. 21). In what follows, we denote the space
just by
and set
where the intersection is taken over all
and all
.
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:
where
,
,
are integer non-negative numbers,
,
,
.
As usual, the space of the weighted generalized functions
is defined as a class of continuous linear functionals that map the test functions from
into the set of real numbers. Each function
u such that
,
will be identified with the functional
that acts according to the formula
The generalized functions
defined by Equation (
11) with
will be called regular weighted generalized functions. All other generalized functions
will be called singular weighted generalized functions.
By
,
, we denote the space of all functions measurable in
even with respect to each variable
,
such that
where
For a real number
, the
-norm of
f is defined by
It is known (see [
17]) that
is a Banach space.
The modified Bessel functions
and
(occasionally called the hyperbolic Bessel functions) of the first and second kind, respectively, are defined as follows (see [
18,
19]):
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
is an even function, i.e., the relation
holds valid. For
(
), the known asymptotic formula for the function
is as follows [
20]:
where
is the Euler–Mascheroni constant.
For
, we have another known asymptotic formula:
Some of our results will be formulated in terms of the generalized translation defined as follows (see [
21]):
For
, the generalized translation
is reduced to the following simple translation:
The multidimensional generalized translation is defined by the equality
where each of the one-dimensional generalized translations
acts according to Equation (
17) for
and according to Equation (
18) for
.
Another object useful for our derivations is a generalized convolution generated by a multidimensional generalized translation
. It is defined as follows:
It is worth mentioning that the operators
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):
where
In particular, we obtain
Equation (
21) is called the Taylor–Delsarte formula. Letting
m go to infinity in Equation (
21), we arrive at the the Taylor–Delsarte series
The Taylor–Delsarte Formula (
21) can be easily generalized to the case of functions depending on several variables. In what follows, by
we denote a segment of the multidimensional Taylor–Delsarte series, where
is a multi-index containing the non-negative integers,
,
,
and
.
In this paper, we also employ the multidimensional Hankel transform with the kernel defined by Equation (
9) acting on the functions
. It is provided by the formula
where
is defined by Equation (
6), the kernel
is given by Equation (
9), and
. Let
be a function of bounded variation in a neighborhood of a point
x of its continuity. Then, for
, the following inversion formula for the multidimensional Hankel transform holds valid:
For the multidimensional Hankel transform and the following convolution operator with respect to the generalized convolution
the convolution theorem
holds valid (see [
16], p. 156, formula 3.176).
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
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):
Now we consider the convolution operators
with the kernel
To obtain a suitable integral representation of the kernel
, we change to the spherical coordinates
,
in the integral at the right-hand side of the previous relation and then apply Equation (
8):
Now we introduce the notion of a B-multiplier in
. Let
. The weighted generalized function
M is called a B-multiplier in
if for
, the generalized convolution
belongs to the space
and the inequality
holds true. The linear space of all B-multipliers is denoted by
. On the space
, the left-hand side of the inequality in Equation (
30) defines a norm that is used in the further discussion.
Lemma 3. Let . The operatoris bounded in , . Proof. Let be a fixed number. The function is a B-multiplier for any . Since the kernel function is represented as an operator generated by the B-multiplier in , the inclusion holds valid.
Because the operator
is a generalized convolution with a B-multiplier, its image also belongs to
. □
Lemma 4. Let . Then the formulaholds valid, where is the generalized convolution with the Poisson kernel Proof. By the convolution theorem, we get the relation
Let us denote
by
. Then we have
that leads to the representation
Applying Equation (
26), we arrive at the relation
where
is the generalized Poisson kernel in Equation (
25). By Lemma 2, the inclusion
holds valid. □
Theorem 1. Let , andThen the relationholds valid, where the limit is understood in the sense of the norm in . Proof. According to Lemma 4, we have to prove that
Taking into account the 2nd property from Lemma 2, we get the representation
Applying the generalized Minkowski inequality [
16] to the integral at the right-hand side of the last formula, we have
Lemma 3.6 on p. 166 from [
25] ensures that
for
. The Proposition 4.1, p. 182 from [
26] and p. 50 from [
27], leads to the formula
By the Lebesgue majorized convergence theorem, the integral in Equation (
33) tends to zero as
since the integrand is majorized by the integrable function
. □
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
be integrable in a layer
for any
,
and the representation
hold valid, where
,
b, and
are some constant positive numbers independent of
A. If the limit
exists, then it is called the Hadamard finite part of the singular integral of the function
f. The function
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
In the case
in Equation (
34), the function
) is said to possess the non-logarithmic-type Hadamard property at the origin.
If a function
is integrable on the set
and possesses the Hadamard property at the origin, then its finite part is defined as follows:
We start by mentioning that the asymptotic properties of the modified Bessel function
ensure that the kernel function
exponentially decays at infinity (see Equation (
16)) and goes to a constant at the origin (see Equation (
15)):
Therefore,
Let as note that the function
in the denominator of
has the simple poles at the points
, and therefore, the weighted generalized function
has the simple poles at the same points.
As already mentioned in the previous section, the operator
can be represented in terms of the generalized convolution, as in Equation (
20)
Then we can apply the convolution theorem for the multidimensional Hankel transform (Equation (
24)) and obtain the relation
Furthermore, the inclusion
and the formula
hold valid, and we arrive at the representation
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
The inverse operator
can be formally represented in terms of the multidimensional Hankel transform as follows:
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
, 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
in the following form
For further derivations, some advanced results regarding asymptotic behavior of the modified Bessel function
are needed. It is known ([
28], p. 273) that the inclusion
,
holds valid. Moreover, we have the estimates
for large values of
r and the estimates
for small values of
r, where
if
m is odd, and
if
m is even, and
is not an integer. If
, then
For
, the formula
holds valid.
Equations (
40)–(
42) lead to the following statement:
Lemma 5. For , the inclusion holds valid as well as the formulas In the next theorem, a regularization of the generalized convolution is introduced using the appropriate segment of the Taylor–Delsarte series of the function .
Theorem 2. Let and . Then the representationholds valid, where , is a segment of the Taylor–Delsarte series, is the multi-index of non-negative integers, , , and . Proof. We start with the representation
where
is a segment of the Taylor–Delsarte series,
. Then we have the relation
In order to derive a suitable representation for
, let us consider an integral with an addition parameter
and suppose that
is fixed.
Passing to the spherical coordinates in the integral
, we first obtain the representation
Using the formula 1.107, p. 49 from [
16], we then get
and
The next step is applying the formula 2.16.2 from [
22] in the form
In our case, the parameters are as follows:
,
, and the convergence condition takes the form of
or
since
. Thus, we can calculate the integral
under the condition
:
The formula above can be analytically extended with respect to
to the half-plane
. Setting
in this extended formula, we arrive at the representation
Substituting the value of the integral
into the formula
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
and
, we define the operator
that turns to be inverse to the Bessel B-potential
.
Theorem 3. On the space , the operator , defined by Equation (46) is the right- and left-inverse to the Bessel B-potential, i.e., the formulashold valid for any . Proof. The main idea of the proof is first to consider the generalized convolution of the weighted functional
with a function
f from
. Then an analytic continuation of the weighted functional
with respect to
is constructed. It turns out that the generalized convolution of this analytic continuation with a function
f from
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
in the sense of the weighted distributions from the space
. On the space
, the operator inverse to
has the form
For
, 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
that is not an integer number. It is known (see [
21]) that the inclusion
holds valid for
and
. Thus, for
, the weighted functional
,
can be represented in terms of an appropriate segment of the Taylor–Delsarte series:
where
. The integral at the right-hand side of Equation (
48) is analytical with respect to the parameter
in the domain
since the expression
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
. Thus, under the condition
, we arrive at the representation
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
and because all weighted functions
are analytic with respect to
(see Chapter 4 in [
16]). Therefore, the representation in Equation (
49) also keeps its validity in the case
,
. Using the relation
, we rewrite the representation from Equation (
48) for the function
as follows:
where
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
Applying the inverse Hankel transform to the last formula, we immediately arrive at the desired relation
. The second formula from Equation (
47), i.e.,
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
. In this case, for
, we have the formula
and the inverse operator in Equation (
46) takes a simple form