Remarks on the Generalized Fractional Laplacian Operator

The fractional Laplacian, also known as the Riesz fractional derivative operator, describes an unusual diffusion process due to random displacements executed by jumpers that are able to walk to neighbouring or nearby sites, as well as perform excursions to remote sites by way of Lévy flights. The fractional Laplacian has many applications in the boundary behaviours of solutions to differential equations. The goal of this paper is to investigate the half-order Laplacian operator (−∆) 1 2 in the distributional sense, based on the generalized convolution and Temple’s delta sequence. Several interesting examples related to the fractional Laplacian operator of order 1/2 are presented with applications to differential equations, some of which cannot be obtained in the classical sense by the standard definition of the fractional Laplacian via Fourier transform.


Introduction
Let s ∈ (0, 1) and ∆ = ∂ 2 /∂x 2  1 The fractional Laplacian of a function u: R n → R is defined as: where P.V. stands for the Cauchy principal value, and the constant C n,s is given by: s, On the other hand, the fractional Laplacian in R n can be written by the Fourier transform: (−∆) s u(x) = 1 (2π) n R n |ζ| 2s (u, e −iζx ) L 2 e iζx dζ = F −1 {|ζ| 2s F (u)(ζ)}(x) where: Hence, the fractional Laplacian is really a pseudo-differential operator with symbol |ζ| 2s .Let L p be the Lebesgue space with p ∈ [1, ∞), B 0 be the space of continuous functions vanishing at infinity, and B bu be the space of bounded uniformly-continuous functions.M.Kwaśnicki recently presented ten equivalent definitions in [1] for defining (−∆) s over these three spaces, including the above Fourier definition.
There have been many studies, including numerical analysis approaches, on the fractional Laplacian with applications in solving certain differential equations on bounded domains and in the theory of stochastic processes and anomalous diffusion [2][3][4][5][6][7][8][9][10].For example, the work in [11] used the fractional Laplacian for linear and nonlinear lossy media, as well as studying a linear integro-differential equation wave model.The work in [12] studied a finite difference method of solving parabolic partial integro-differential equations, with possibly singular kernels.These arise in option pricing theory when the random evolution of the underlying asset is driven by a Lévy process related to the fractional Laplacian or, more generally, a time-inhomogeneous jump-diffusion process.Using the fractional Laplacian operator, Araci et al. [13] investigated the following q-difference boundary value problem: with the boundary conditions: They proved the existence and uniqueness of a positive and nondecreasing solution for this problem by means of a fixed point theorem involving partially-ordered sets.
Let Ω ⊂ R n denote a bounded, open domain.For u(x): Ω → R, D'Elia and Gunzburger [14] investigated the action of the nonlocal diffusion operator L on the function u(x) as: where the volume of Ω is non-zero and the kernel γ(x, y): Ω × Ω → R is a non-negative symmetric mapping, as well as the nonlocal, steady-state diffusion equation: An example of γ(x, y) is given by: with σ(x, y) bounded from above and below by positive constants.This nonlocal diffusion operator has, as special cases, the fractional Laplacian and fractional differential operators that arise in several applications.The corresponding evolution model was further studied in [15].Recently, Hu et al. [16] studied the following high-dimensional Caputo-type parabolic equation with the fractional Laplacian by using the finite difference method: where α ∈ (0, 1), s ∈ (0, 1) and Ω ⊂ R n .In particular, this involves the half-order Laplacian operator (−∆) On the other hand, the Riesz fractional derivative is generally given as: where 0 < α < 2 and α = 1, and: Note that in space fractional quantum mechanics, the α = 2 case corresponds to the Schrödinger equation for a massive non-relativistic particle, while the α = 1 case needs to be examined carefully, both on physical and mathematical grounds, since Equation ( 2) is undefined for α = 1.
It follows from [10,[17][18][19][20][21][22] that: Applying two identities: we get: dζ, 1 < α < 2 by removing the term: due to the meaning of the finite part integral as the integral is divergent and the finite part: With the same argument, we come to: Hence, we have another representation of the Riesz fractional derivative from Equation (2): Similarly, we can claim that this representation still holds for the entire range 0 < α ≤ 2 [23].In particular, based on the formula: Clearly, Equation (1) becomes: for s = 1/2 and n = 1.Therefore, investigations of the half-order Laplacian operator (−∆) 2 on R are equivalent to studies of the first-order Riesz derivative.We can define −(−∆) 1 2 u as the Riesz derivative RZ D α x u(x) in the case of α = 1, which is undefined in Equation ( 2) in the classical sense.The aim of this work is to study the operator (−∆) 1 2 on R in distribution explicitly and implicitly, using a particular delta sequence and the generalized convolution.We also present several interesting examples, such as (−∆) , which are undefined in the classical sense.At the end of this work, we describe applications of such studies to solving the differential equations involving the half-order Laplacian operator.

The Explicit Approach to (−∆) 1/2 u
In order to extend the fractional Laplacian (−∆) 1/2 distributionally, we briefly introduce the following basic concepts of distributions.Let D(R) be the Schwartz space (testing function space) [24,25] of infinitely-differentiable functions with compact support in R and D (R) the (dual) space of distributions defined on D(R).
and only if these functions vanish outside a certain fixed bounded set and converge to zero uniformly together with their derivatives of any order.
The functional δ (n) (x − x 0 ) is defined as: where φ ∈ D(R).Clearly, δ (n) (x − x 0 ) is a linear and continuous functional on D(R), and hence, Define the unit step function θ(x) as: Then, The distributional derivative of f , denoted by f or d f /dx, is defined as: Clearly, f ∈ D (R) and every distribution has a derivative.As an example, we are going to show that θ (x) = δ(x) distributionally, although θ(x) is not defined at x = 0 in the classical sense.Indeed, which claims: It can be shown that the ordinary rules of differentiation also apply to distributions.For instance, the derivative of a sum is the sum of the derivatives, and a constant can be commuted with the derivative operator.Definition 1.Let f and g be distributions in D (R) satisfying either of the following conditions: (a) either f or g has bounded support (set of all essential points).(b) the supports of f and g are bounded on the same side (either on the left or right).
Then, the convolution f * g is defined by the equation: for φ ∈ D. Clearly, we have: It follows from the definition above that: for any distribution f ∈ D (R).
Let δ n (x) = nρ(nx) be Temple's δ-sequence for n = 1, 2, • • • , where ρ(x) is a fixed, infinitely-differentiable function on R, having the following properties [26,27]: An example of such a ρ(x) function is given as: where: This delta-sequence plays an important role in defining products of distributions [28,29].Let f be a continuous function on R.Then, uniformly converges to f on any compact subset of R. Indeed, we assume that f is continuous on R, and L is any compact subset of R.Then, f is uniformly continuous on L, and for all > 0, there exists for all x ∈ L and |t| < δ 1 .This implies that: for all x ∈ L and 1/n < δ 1 by noting that: by noting that: lim in D(R).
In order to study the half-order Laplacian operator in the distribution, we introduce an infinitely-differentiable function τ(x) satisfying the following conditions: Define: From Equation (3), we have: This suggests the following explicit definition for defining (−∆) . This explicit definition directly evaluates the half-order fractional Laplacian of u(x) as a function of x, without relating to any testing function in the Schwartz space.Definition 2. Let u ∈ D (R) and u n = u * δ n = (u (t), δ n (x − t)) for n = 1, 2, • • • .We define the half-order Laplacian operator (−∆) if it exists.
Clearly, the integral: is well defined as φ m (t) has a bounded support and: Proof.Assuming x > 0, we choose a large n such that 1/n < x.This infers that δ n (x + t) = 0, and: from Definition 2. Making the substitution w = x − t, we get: by noting that δ n (w) is a delta sequence, and: x − w is a testing function of w if w < x, and: x − w is identical to zero if w ≥ x as δ n (w) = 0. Choosing m such that x < m, we derive that: If x < 0, we set y = −x and: which implies that: Finally, we have for x = 0 by making the variable change u = nt: by noting that the term φ m (u/n) = 1, and: is well defined as ρ(u) is an even function.This completes the proof of Theorem 1.
We should note that Theorem 2 cannot be derived by the standard definition of the fractional Laplacian via Fourier transform, and: holds in the distributional sense and: although θ(x) is discontinuous at x = 0 in the classical sense.(−∆) (−∆) (−∆) where a and b are arbitrary constants.
Proof.Clearly, sin x is an infinitely-differentiable function on R.This claims that: uniformly converges to cos x on any compact subset of R. Therefore, (−∆) Similarly, (−∆) Finally, This completes the proof of Theorem 3.
We would like to mention that Theorem 3 extends several classical results obtained in [30] to distributions by Definition 2 with a new approach.(−∆) Proof.Clearly, (arctan x) = 1 1 + x 2 is continuous on R, and uniformly converges to 1/(1 + x 2 ) on any compact subset of R. Therefore, (−∆) by noting that the integral: is well defined for every point x ∈ R. It remains to show that: for all x ∈ R. First, we note that: is even with respect to z and has two singular points z 1 = x + i and z 2 = i − x in the upper half-plane.Clearly, we have for x = 0 that: . By Cauchy's residue theorem, we get: This implies that: for all nonzero x.Furthermore, we derive for x = 0 that: π by using the identity [31]: This completes the proof of Theorem 4.

3.
The Implicit Approach to (−∆) It seems infeasible to calculate directly the fractional Laplacian operator of some functions or distributions by Definition 2. For example, (−∆) where: uniformly converges to e x as it is continuous on R. Clearly, the right-hand side of the above integral is divergent as: In this section, we are going to provide another definition for dealing with (−∆) 1 2 u(x) efficiently, based on a testing function with compact support.This definition is implicit and only used to define the meaning of: ((−∆) rather than finding an explicit function of x.It clearly makes sense in the distribution as we regard (−∆) 1 2 u(x) as a functional (not a function) on the Schwartz testing space D(R).We must point out that this implicit-definition, using a different generalization, is independent of the explicit one provided in Section 2. Definition 3. Let u ∈ D (R) and u n = u * δ n = (u (t), δ n (x − t)) for n = 1, 2, • • • .We define the half-order Laplacian operator (−∆) if it exists.
Clearly, the integral: is well defined, as φ(t) has a bounded support and: It follows from Definition 3 that: by noting that: is well defined and is a number for each φ, since e x φ(x) is also a testing function in the Schwartz space.Furthermore, ((−∆) Theorem 5. Let u(x) be an infinitely-differentiable function satisfying: Then, ((−∆) Proof.Clearly, and: uniformly converges to u(x) on any compact subset of R. By Definition 3, ((−∆) We should note that the term: is well defined for every point x ∈ R, and it is indeed not Cauchy's principal value as: This completes the proof of Theorem 5.
It follows from Theorem 5 that: ((−∆) which implies that: (−∆)  Remark 1.To end this section, we must point out that if we replace φ(t) by φ m (t) used in Section 2 and add the limit, then the two operators (−∆)

Conclusions
This paper introduced two independent definitions for defining the fractional Laplacian of the half-order (−∆)

1 2
when s = 1/2.The convergence and error estimate of the established finite difference scheme are shown with several examples.

1 2
are identical for some functions.For instance,

1 2 1 2
in the distribution, both explicitly and implicitly.We demonstrate several examples, such as (−∆) 1 2 δ (m) (x) and (−∆) arctan x, some of which are undefined in the classical sense.The results obtained have potential applications in solving the differential equations involving the half-order Laplacian operator.For example, the differential equation: : u(x) = arctan x + ax + b on any non-empty subset of R, and the differential equation: has a solution: u(x) = πθ(x) + ax + b where a and b are arbitrary constants.