On the General Solutions of Some Non-Homogeneous Div-Curl Systems with Riemann–Liouville and Caputo Fractional Derivatives

: In this work, we investigate analytically the solutions of a nonlinear div-curl system with fractional derivatives of the Riemann–Liouville or Caputo types. To this end, the fractional-order vector operators of divergence, curl and gradient are identiﬁed as components of the fractional Dirac operator in quaternionic form. As one of the most important results of this manuscript, we derive general solutions of some non-homogeneous div-curl systems that consider the presence of fractional-order derivatives of the Riemann–Liouville or Caputo types. A fractional analogous to the Teodorescu transform is presented in this work, and we employ some properties of its component operators, developed in this work to establish a generalization of the Helmholtz decomposition theorem in fractional space. Additionally, right inverses of the fractional-order curl, divergence and gradient vector operators are obtained using Riemann–Liouville and Caputo fractional operators. Finally, some consequences of these results are provided as applications at the end of this work.


Introduction
From an analytical point of view, the literature offers a wide range of reports that focus on the extension of integer-order methods and results for the fractional case. From a more particular point of view, the fractional generalization of the classical vector calculus operators (that is, the gradient, divergence, curl and Laplacian operators) has been also an active topic of research, which has been developed from different approaches. Some of the first attempts to extend these operators to the fractional scenario are described in [1,2] using the Nishimoto fractional derivative. These operators were used later on in [3] to provide a physical interpretation for the fractional advection-dispersion equation for flow in heterogeneous porous media. In 2008, Vasily E. Tarasov described different approaches to formulate a fractional form of vector calculus with physical applications in [4] (see also references therein). More recently, a new generalization of the Helmholtz decomposition theorem for both fractional time and space was proposed in [5,6] using the discrete Grünwald-Letnikov fractional derivative. Another related work is [7], where the authors investigate the dynamic creation of fractionalized half-quantum vortices in Bose-Einstein condensates of sodium atoms.
In this work, we consider fractional derivatives of the Riemann-Liouville and the Caputo types and provide extensions of the definitions of the main differential operators from vector calculus using these fractional operators. In such a way, we present fractional forms of the divergence, the rotational and the gradient operators. Moreover, we also consider generalized forms of the Dirac and the Laplace operators using fractional derivatives. Our goal in this work is to extend quaternionic analysis to consider fractional forms of the classical differential operators. Analogues of the properties satisfied by the classical operators will be mathematically established in this work. For instance, we will show that RL D α a + [w] = − RL div α a + w + RL grad α a + w 0 + RL curl α a + w, (1) C D α a + [w] = − C div α a + w + C grad α a + w 0 + C curl α a + w.
In other words, when we apply the fractional Dirac operator to a quaternion valued function w = w 0 + w, this expression can be decomposed in terms of a fractional divergence, a fractional gradient and a fractional rotational operator. Based on this, we will also provide explicit expressions of general weak solutions for some fractional forms of the div-curl system, considering various analytical hypotheses. More precisely, we will prove that if we restrict ourselves to the class of functions with fractional divergence zero and whose Riemann-Liouville fractional integral has zero normal trace, then the fractional Teodorescu transform represents a solution of the above-mentioned div-curl systems (see Theorems 2 and 3).
This manuscript is organized as follows. In Section 2, we recall some important definitions from the literature, including those of the left and right Riemann-Liouville fractional derivatives, the left and right Caputo fractional derivatives and the two-parameter Mittag-Leffler function. A useful characterization of the functions with summable fractional derivatives is also recalled. In Section 2.2, the Riemann-Liouville and the Caputo Dirac and Laplace operators are introduced. Moreover, some fundamental solutions for the fractional Dirac operators are recalled in Section 2.3. In Section 3, we introduce fractional extensions of the divergence, rotational and gradient differential operators. Some properties among these operators are established, and a useful factorization theorem for the fractional Laplace operators is proven. It is worth mentioning that this factorization is not new; however, we were able to derive it only using the identities preserving the fractional gradient, divergence and rotational operators. Among the most important results, we establish that the fractional Teodorescu transform is a right inverse of the fractional Dirac operator under suitable analytical conditions, and we prove a fractional form of the Divergence Theorem. Section 4 is devoted to establishing the existence of weak solutions for Riemann-Liouville and Caputo fractional div-curl systems. The explicit form of the operators involved in the solution, as well as some of their properties, allow the solution to be reexpressed as the sum of the fractional gradient of a scalar potential plus a fractional curl of a vector potential; we can say that our solutions preserve a Helmholtz-type decomposition (see Propositions 6 and 7). As a consequence, right inverses of the fractional rotational and divergence operators are provided in a subclass of the fractional divergence-free vector fields. In turn, Section 5 provides some consequences of the factorization results proven in Section 2 to the construction of fractional hyper-conjugate pairs. A theorem providing necessary and sufficient conditions for the existence of Caputo fractional hyper-conjugate pairs is proven in this stage, along with a result of the existence of a right inverse for the fractional gradient. Finally, this manuscript closes with a section of concluding remarks.

Fractional Calculus
The present section is devoted to recalling some useful definitions from fractional calculus. Throughout, we assume that a, b, α ∈ R satisfy α > 0 and a < b. Meanwhile, we suppose that f : R → R is a sufficiently smooth function, with the property that f is identically equal to zero outside of the interval [a, b].
Definition 1. The left and right Riemann-Liouville fractional integrals of f of order α with respect to the interval [a, b] (whenever they exist) are the functions I α a + [ f ] and I α b − [ f ], defined, respectively, by (see [8]) Let m = [α] + 1 ∈ Z. The left Riemann-Liouville and the left Caputo fractional derivatives of order α with respect to the interval [a, b] are, respectively, defined as follows: Finally, we define the right Riemann-Liouville and the right Caputo fractional derivatives of order α with respect to the point a, respectively, as the functions For the sake of convenience, we will employ the notation D α a ± when we present properties satisfied by both fractional derivatives RL D α a ± and C D α a ± . Throughout, we let I α a + (L 1 ) denote the class of all functions f that are represented by the fractional integral (3) of some integrable function, i.e., f = I α a + [g], for some g ∈ L 1 (a, b). Using this notation, the following result provides a characterization of these functions. Definition 2. If (I m−α a + [ f ]) (k) (a) = 0, for each k ∈ {0, 1 . . . , m − 1}, then it follows that f (k) (a) = 0 holds, for each k ∈ {0, . . . , m − 1} (see [8,9]). In light of this fact, we will say that f has a summable fractional derivative D α In the following discussion, suppose that α > 0, f admits a summable fractional derivative of order α > 0 on [a, b], and let m = [α] + 1. Then, the following composition rules are satisfied: On the other hand, we know that both fractional derivatives RL D α a + and C D α a + satisfy the one-sided invertibility property D α a + I α a + [ f ] = f . It is worth noting that this identity is a particular case of the property D α It is important to recall also that the following semi-group property for the composition of fractional derivatives is not generally satisfied: However, if f (k) (a) = 0 for k = 0, 1, . . . , max{[α] + 1, [β] + 1} − 1, then (11) does hold; see Section 2.2.6 [8]. An analogous condition for the semi-group property in the context of the Riemann-Liouville derivative is found in Section 2.3.6 [8]. Finally, the following relation between the Riemann-Liouville and Caputo fractional derivatives is valid: Definition 3 (Gorenflo et al. [10]). Let µ, ν ∈ R be such that µ > 0. We define the twoparameter Mittag-Leffler function E µ,ν : C → C with parameters µ and ν in terms of the following power series: , ∀z ∈ C.

Fractional Quaternionic Analysis
In this section, we will mention some recent results in fractional Clifford analysis. The Dirac operator in Clifford analysis, also known as the Moisil-Teodorescu differential operator, represents the cornerstone of the analysis in higher dimensions. A remarkable number of systems of differential equations have been analyzed using this operator or a perturbation of it, and the monographs [11][12][13] of the authors Gürlebeck and Sprößig contain many examples of the applications that have been made over the years. See also [14], where the authors introduced the fractional Dirac operator with Caputo derivatives as well as the basic tool of a fractional function theory in more dimensions.
More precisely, this section is devoted to the collection of some recent results of the authors Ferreira et al. [15][16][17], by whom fundamental solutions of the fractional Laplacian were found, where the derivatives are of Riemann-Liouville and Caputo types, as well as of the fractional Dirac operators.
For the remainder, let a i , b i ∈ R satisfy a i < b i , for each i = 1, 2, 3. We will suppose that is a bounded open rectangular domain in R 3 , and let α = (α 1 , α 2 , α 3 ), with α i ∈ (0, 1], for all i = 1, 2, 3. Definition 4 (Ferreira et al. [15][16][17]). The fractional Riemann-Liouville and fractional Caputo Dirac operators are represented by RL D α a + and C D α a + , respectively, and they are defined as Here, RL ∂ 1+α i 2 x i ,a i and C ∂ 1+α i 2 x i ,a i are, respectively, the Riemann-Liouville and the Caputo fractional derivative operators of order (1 + α i )/2 with respect to the variable x i ∈ (a i , b i ), for each i = 1, 2, 3. We define the fractional Laplace operators RL ∆ α a + and RL ∆ α a + , respectively, by where RL ∂ 1+α i x i ,a i and C ∂ 1+α i x i ,a i are, respectively, the Riemann-Liouville and Caputo fractional derivatives of order 1 + α i with respect to the variable x i ∈ (a i , b i ), for each i = 1, 2, 3.

Fundamental Solutions
The purpose of this subsection is to determine fundamental solutions for the fractional Dirac operator and use their properties in the investigation of the solutions of fractional div-curl systems. Beforehand, it is worth recalling that a family of fundamental solutions for the fractional Laplace operators RL ∆ α a + and C ∆ α a + , and a family of fundamental solutions for the fractional Dirac operators RL D α a + and C D α a + , were obtained in [15,16] for the Riemann-Liouville and Caputo case, respectively. In the case of Riemann-Liouville fractional operators, the authors employed some properties of the Mittag-Leffler function and the Laplace transform in two dimensions.
We will begin this section by recalling some relevant results derived in [15]. To this end, let u be an eigenfunction of the fractional Laplace operator, i.e., suppose that RL ∆ α a + u = λu for some λ ∈ C, and assume that u( x) admits a summable fractional derivative RL ∂ 1+α 1 2 x 1 ,a + 1 in the variable x 1 , and that it belongs to I In what follows, we will consider the following integral and differential conditions of Cauchy type: Lemma 1 (Ferreira et al. [15,17]). A family of eigenfunctions of the fractional Laplace operator RL ∆ α a + is given by the function Meanwhile, a family of fundamental solutions of the fractional Dirac operator RL D α a + is obtained by considering λ ≡ 0 in (19). More precisely, this family of solutions is given by where u 0 is a fundamental solution of RL ∆ α a + , i.e., Here, f 0 and f 1 satisfy the conditions (18).
Let v be an eigenfunction of the fractional Laplace operator, i.e., suppose that C ∆ α a + v = λv, for some λ ∈ C. Assume that v( x) admits a summable fractional derivative RL ∂ 1+α 1 2 x 1 ,a + 1 in the variable x 1 , and that it belongs to I In what follows, we will consider the following Cauchy conditions: As a consequence, g 0 (a 2 , Lemma 2 (Ferreira et al. [16,17]). A family of eigenfunctions for the fractional Laplace operator C ∆ α a + is given by the function Meanwhile, a family of fundamental solutions of the fractional Dirac operator C D α a + is obtained by considering λ ≡ 0 in (23). More precisely, this family of solutions is given by where g 0 (x 2 , x 3 ) and g 1 (x 2 , x 3 ) satisfy (22).

Fractional Vector Calculus
For the remainder of this section, we will study the fractional divergence, gradient and rotational operators as parts of a decomposition of the fractional Dirac operator in three dimensions. More precisely, recall that if w = w 0 + w is a quaternionic-valued function, then the following decomposition in quaternionic form is satisfied: Here, D = ∑ 3 i=1 e i ∂ i is the classical Dirac operator, which is also called the Moisil-Teodorescu differential operator. For more details about quaternionic analysis, see [11,18,19]. Our goal in this section is to provide an extension of this decomposition (27) using fractional operators of the Riemann-Liouville and Caputo types.
Let w = w 0 + ∑ 3 i=1 e i w i be a quaternionic-valued function in AC(Ω), whose scalar part is denoted by Sc[w] = w 0 and its vector part by Vec[w] = w = ∑ 3 i=1 e i w i . Then, the action of the operator RL D α a + on w reduces to The above decomposition of Equation (27) originates a fractional version of the classical divergence, rotational and gradient differential operators from vector calculus. These operators are, respectively, the scalar component of (27), the vector term acting over w and the vector term of the equation acting over w 0 . These facts are the motivation to analyze the following fractional differential operators.

Definition 5.
Define the fractional divergence, curl and gradient operators in the Riemann-Liouville sense by It is important to point out that the fractional operators (28)-(30) reduce, respectively, to the classical div, curl, and grad operators from vector calculus when α i = 1, for each i = 1, 2, 3. Moreover, if α * > 0 and α i = α * , for each i = 1, 2, 3, then the above fractional operators coincide with the divergence, curl and gradient operators defined in [1,2] up to a constant factor. See also [4] and references therein for a historic account of the efforts to formulate a fractional form of vector calculus. Unlike the classical vector calculus operators, these fractional operators are non-local. Consequently, the fractional divergence, curl and gradient depend on the domain Ω.
Notice now that (27) can be rewritten as the following decomposition Since the specific form of the Riemann-Liouville fractional derivative does not affect the above decomposition, we analogously obtain the following decomposition in terms of Caputo fractional derivatives: Here, the operators C div α a + , C curl α a + and C grad α a + are defined as in (28)-(30), respectively, using Caputo fractional derivatives instead of Riemann-Liouville operators.
We define now a class of functions in AC 1 (Ω) where we can apply the semi-group property (11).
. Moreover, the identities (i)-(iv) also hold for the Caputo fractional operators.
Proof. The results readily follow from the identities which are trivially satisfied in Z α a + (Ω). The identities with Caputo fractional operators are established in a similar fashion. Similar identities using Caputo fractional derivatives were proven in [4] when all the orders of the fractional derivatives are equal, i.e., when there exists α * > 0 such that α i = α * , for each i = 1, 2, 3. On the other hand, a direct consequence of Proposition 1 is that the fractional Dirac operator RL D α a + factorizes the fractional Laplace operator RL ∆ α a + . A more general factorization for functions taking values in Cl 0,3 was proven in Section 4 [16]. In light of these remarks, the following result is a direct consequence of Proposition 1.
a + (Ω), then the following factorizations of the fractional Laplace operators are satisfied: We turn our attention now to the right Caputo fractional Dirac operator, which is given by the expression The following fractional Stokes formula was proven in Theorem 10 [17]: Here, we require that f , h ∈ AC 1 (Ω) ∩ AC(Ω) and I α a It is worth noting here that the operator C D α b − acts on the right, while RL D α a + acts on the left. Intuitively, the last formula shows that the left Riemann-Liouville and right Caputo fractional Dirac operators act by 'intertwining' to obtain the fractional analogue of the Stokes formula.
The following result is a fractional form of the well-known Divergence Theorem.
Before closing this section, it is natural to compare qualitatively the results obtained in traditional vector calculus against those in the fractional case. In classical vector calculus, the following product rules are satisfied: On the other hand, in the fully fractional case considered in Proposition 3, when we restrict f to the class of functions Z α a + (Ω), the first part of these identities is also satisfied, except that the second terms on the right-hand sides of (46) and (47) are not present anymore. Notice that it is not difficult to construct a family of functions belonging to Z α a + (Ω), for instance f (x) = (x 1 − a 1 ) γ 1 (x 2 − a 2 ) γ 2 (x 3 − a 3 ) γ 3 g(x) for all g ∈ AC 1 (Ω) and γ i ≥ 0 for all i = 1, 2, 3.

Properties of the Fractional Teodorescu Transform
As a derivation of the fractional Borel-Pompeiu formula [17], the authors defined the Caputo-type Teodorescu transform in a very similar way to the following definition, the difference being that the kernel is now a fundamental solution of the fractional Dirac operator defined in (20). Definition 7. Let x = (x 1 , x 2 , x 3 ) with x i > a i , for all i = 1, 2, 3. Define the Riemann-Liouville and Caputo fractional Teodorescu transform by Here, we follow the nomenclature and conventions of Lemma 1. Moreover, the derivatives RL E α a + = − RL D α a + [u 0 ] and C E α a + = − C D α a + [v 0 ] that appear in the kernel of (48) and (49), respecrively, are with respect to the variable x.
In the following, it will be convenient to employ the translation operator, which is defined by T z f ( y) = f ( y + z). Similarly, we will use the reflection operator given by R y f ( y) = f (− y). An important relation between D α a + and T θ is where the derivative is taken with respect to the variable x. The proof is analogous to Theorem 11 [17], which is for the Caputo case, but we will give it here for completeness.
Proof. Observe firstly that the fundamental solution RL E α a + ( x) of the fractional Dirac operator RL D α a + is defined only for x i > a i , for each i = 1, 2, 3. Moreover, it satisfies the iden- . In the following, the derivatives RL D α y + and RL D α a + are with respect to the variable x. Using (50) with θ = a − y yields which we wished to prove.
In the following, we will employ a key decomposition of the classical Teodorescu operator used in [20][21][22][23] for different kinds of bounded or unbounded domains in R 3 . For the fractional version, we denote the component operators of the fractional Teodorescu transform as follows: The first term on the right-hand side of Equation (52) is the scalar part, while the last two summands represent the vector part, and it has been split into two components for the sake of convenience. These three terms are given by We will see later in Corollary 2 how some of these component operators themselves represent right inverse operators of the fractional divergence and rotational operators, under certain conditions. Moreover, RL T α a + , as a good generalization of the classical Teodorescu operator in quaternionic analysis, preserves many of its properties. To this end, we use the above decomposition (52), RL , to see necessary and sufficient conditions to guarantee that both its scalar part and its vector part belong to the kernel of the fractional Laplacian RL ∆ α a + . In order to apply the factorization of Corollary 1 to the fractional Teodorescu transform, we will need to guarantee that the condition RL T α a + [ g] ∈ Z α a + (Ω) can be satisfied.
, belongs to the kernel of RL ∆ α a + if and only if RL curl α a + g = 0. Moreover, the statements (i) and (ii) also hold for the Caputo fractional Teodorescu transform.
Proof. Using the factorization in Corollary 1, Proposition 4 and the decomposition (31), it readily follows that Taking its scalar part or vector part, respectively, we obtain the desired result.

Riemann-Liouville System
In the following, we will study a fractional form of the classical div-curl system and construct its solution. More precisely, we fix the domain Ω, and consider the fractional system where g 0 ∈ L p (Ω, R) and g ∈ L p (Ω, R 3 ), for some 1 < p < 2/(1 − α * ). Notice that if the solution is such that w ∈ Z α a + (Ω), then g satisfies RL div α a + g = 0, i.e., g is a 'fractional divergence-free' vector field. Let f 0 ∈ AC(Ω). The following relations will be fundamental in the sequel: Here, the derivatives are taken with respect to the variable y. In the sequel and for the sake of convenience, we will employ Ker to denote the kernel of operators. As in the previous section, we will let Ω = Π 3 i=1 (a i , b i ) be a bounded open rectangular domain in R 3 , and assume that α = (α 1 , α 2 , α 3 ), with α i ∈ (0, 1), for each i = 1, 2, 3. Using this notation, we have the following result.
= 0 and the normal trace of I α a + [ g] vanishes, then a general weak solution of the fractional Riemann-Liouville div-curl system (57) is given by where h ∈ Ker( RL ∆ α a + ) ∩ Z α a + (Ω) is an arbitrary scalar function.
Proof. By Proposition 4 and the decomposition (31), the fractional Teodorescu transform RL T α a + [−g 0 + g] is a quaternionic solution of the system (57). To obtain a pure-vector solution, note that the decomposition (52) yields Taking h( y) = T x+a R y u 0 ( y) and f ( y) = g( y) in the Stokes formula (38) (all derivatives with respect to the variable y) and letting u 0 be a fundamental solution of RL ∆ α a + given in Lemma 1, we obtain Calculate the scalar part of (61) and use the hypotheses −Sc D α By differentiating now under the integral sign and using the traditional Leibniz' rule, we readily obtain the following fundamental relation: Here, the second sub-index indicates whether we are taking derivatives with respect to the variable x or y. Recall that RL E α a + = − RL grad α a + , x [u 0 ] with respect to the variable x. On the other hand, due to u 0 ∈ Z α a + (Ω) and relations (12), (63) and (58), we obtain However, we know that RL E α a + ( x + a − y) is only defined when x i > y i , for all i = 1, 2, 3. As a consequence, we can readily replace the operator C grad α b − by C grad α x − in (62). Finally, by employing (64) and (62), we readily reach This means that RL T α a is purely vectorial and, moreover, Setting the scalar and vector parts equal to each other, we obtain that RL T α a is a solution of the fractional div-curl system (57). Finally, the fact that the solution is not unique is a consequence of the identities (ii) and (iii) of Proposition 1.
In the limit α i → 1 − , the fractional div-curl system (57) reduces to the well-known integer-order system from vector calculus, and the hypothesis RL div α a + [ g] = 0 reduces to the evident requirement that g be a divergence-free vector field. Moreover, I This means that it is sufficient to require that g has zero normal trace. Taking g 0 ≡ 0 or g ≡ 0 in (57), we readily obtain Corollary 2. Under the same assumptions of Theorem 2, RL − → T α 2,a + is a right inverse operator of RL curl α a + in the class of functions considered in Theorem 2. Meanwhile, − RL − → T α 1,a + is always a right inverse operator of RL div α a + in L p (Ω).
Another important analogy with the classical vector calculus is that the solution of the non-linear fractional system (57) also admits a Helmholtz-type fractional decomposition as follows.

Proposition 6.
Under the same assumptions of Theorem 2, the general weak solution of the fractional Riemann-Liouville div-curl system (57) admits a fractional Helmholtz decomposition as follows where the scalar potential ϕ 0 and the vector potential ϕ are given by Proof. Due to RL E α a + = − RL grad α a + [u 0 ] and by relation (58), we can write Finally, since the fractional gradient involved in the above re-expressions of the component operators RL − → T α 1,a + and RL − → T α 2,a + is taken with respect to the variable x, we readily obtain (67).

Caputo System
For the corresponding fractional div-curl system in the sense of Caputo derivatives, we will follow the same approach as that used with the Riemann-Liouville div-curl system (57). Let us consider the system C div α a + w = g 0 , C curl α a + w = g, where g 0 ∈ L p (Ω, R) and g ∈ L p (Ω, R 3 ), for some 1 < p < 2/(1 − p * ). The cornerstone in our analysis will be again the fractional Teodorescu transform associated with Caputo derivatives, which means that its kernel is a fundamental solution of the fractional Dirac operator C D α a + . As in the case of the Riemann-Liouville fractional Teodorescu operator, we define the decomposition where C T α 0,a + , C − → T α 1,a + and C − → T α 2,a + are given by (53), (54) and (55), respectively, but using now the Caputo kernel C E α a + in the integrand, instead of the Riemann-Liouville kernel RL E α a + . Analogously to the proof of Theorem 2, it follows that C T α a + [−g 0 + g] is a quaternionic solution of the fractional div-curl system (68). This is a direct consequence of the fact that C T α a + is a right inverse of C D α a + in L p (see [17], Theorem 11). Apply the decomposition (69) to the quaternion-valued function g = g 0 + g, where g 0 and g are the known data provided by the fractional div-curl system (68). In this way, we notice that To obtain a purely vectorial solution, we will impose suitable conditions over g in order to guarantee that the scalar part of C T α a + [−g 0 + g] vanishes in Ω, i.e., that C T α 0,a + [ g] ≡ 0 is satisfied in Ω. We will see that the fractional divergence-free functions of the Riemann-Liouville type, whose Riemann-Liouville fractional integral has zero normal trace, belong to the kernel of the operator C T α 0,a + , in the same way as seen in Theorem 2 for the operator RL T α 0,a + . Nevertheless, the upcoming proof is relatively more straightforward.
Proof. It suffices to prove that C T α 0,a + [ g] = 0 in Ω, in light of the previous discussion. Take h( y) = T x+a R y v 0 ( y) and f ( y) = g( y) in the Stokes formula (38), and let v 0 be the fundamental solution of C ∆ α a + (25). It follows that Taking now the scalar part of (71) and using the hypotheses, we readily obtain As a consequence of (58), we . However, we know that C E α a + ( x + a − y) is only defined for x i > y i , for all i = 1, 2, 3. This implies that we can replace the operator C grad α b − with C grad α x − in (72). It is easy to see then that whence the conclusion readily follows.

Corollary 3.
Under the same assumptions of Theorem 3, C − → T α 2,a + is a right inverse operator of C curl α a + in the class of functions considered in Theorem 3. Moreover, − C − → T α 1,a + is a right inverse operator of C div α a + in L p (Ω).
Analogously to the Riemann-Liouville case, the fractional Caputo div-curl system (68) also admits a fractional Helmholtz decomposition, but now the potential is in terms of v 0 defined in (25), which is a fundamental solution of the fractional Laplace operator of Caputo type.

Proposition 7.
Under the same assumptions of Theorem 3, the general weak solution of the fractional Caputo div-curl system (68) admits a fractional Helmholtz decomposition as follows where the scalar potential ψ 0 and the vector potential ψ are given by Proof. The proof is analogous to that of Proposition 6 for the Riemann-Liouville case.

Application
Let Ω = Π 3 i=1 (a i , b i ) be as in the previous sections, and assume 0 < α i < 1, for all i = 1, 2, 3. The present section provides some consequences of the factorization provided by Corollary 1 to the construction of fractional hyper-conjugate pairs. In addition, we will give an explicit expression of a right inverse of the fractional gradient of Caputo type considering different derivative orders. Definition 8. Let w = w 0 + w ∈ AC(Ω). We say that (w 0 , w) is a Riemann-Liouville fractional hyper-conjugate pair if w ∈ Ker( RL D α a + ). Analogously, (w 0 , w) is a Caputo fractional hyperconjugate pair if w ∈ Ker( C D α a + ).
The following result is a straightforward consequence of the factorization of the fractional Laplace operator in the class Z α a + (Ω) provided by Corollary 1. For this reason, we omit the proof.
By Definition 8, it is obvious that (w 0 , w) forms a Riemann-Liouville fractional hyperconjugate pair if and only if the following fractional div-curl system is satisfied: Similarly, (w 0 , w) is a Caputo fractional hyper-conjugate pair if and only if The above systems (75) and (76) can be considered fractional generalizations of the Moisil-Teodorescu system studied for the first time in [24].
Let us define the following integral operator in terms of the Riemann-Liouville frac- As the following result shows, it turns out that A α a + behaves as a right-inverse operator of C grad α a + in the class of functions satisfying C curl α a + f = 0. For this reason, A α a + is called the fractional anti-gradient operator.
Proof. Using the characterization of Caputo fractional hyper-conjugate pairs given under Corollary 4 and differentiating under the integral sign, we readily obtain Now, by hypothesis C curl α a + f = 0 or, equivalently, the following identities are satisfied: Substituting (79) into (78) and using the composition rule (10), it follows that Analogously, one can establish that C ∂ 1+α i 2 , for i = 2, 3. The conclusion readily follows now.
The next proposition shows that the fractional anti-gradient operator A α a + allows us to construct a Caputo fractional hyper-conjugate pair w 0 when w ∈ Ker( C ∆ α a + ) is known beforehand. The proposition is clearly a generalization of [20], Proposition 2.1. Proposition 9. Let w ∈ Ker( C ∆ α a + ) ∩ Z α a + (Ω). A necessary and sufficient condition for the existence of a Caputo fractional hyper-conjugate pair of w is that C div α a + w = 0. In that case, there is w 0 such that w = w 0 + w ∈ Ker( C D α a + ).
Proof. The necessity is clear due to the characterization of Caputo fractional hyperconjugate pairs provided by (76). Suppose now that C div α a + w = 0. By Proposition 1(iv), it follows that C curl α a + C curl α a + w = 0. On the other hand, Proposition 8 ensures that C grad α a + A α a + [ C curl α a + w] = C curl α a + w.
The conclusion of this result follows from (76) if we let w 0 = −A α a + [ C curl α a + w].

Conclusions
In this work, we extend some results from vector calculus to the fractional case-for instance, the space fractional Helmholtz Decomposition Theorem provided by Propositions 6 and 7.
The key tools used are the decompositions of the fractional Teodorescu transform in the Riemann-Liouville case (52) and in the Caputo case (69) as well as various properties associated with these fractional operators, which are thoroughly established in this manuscript. To this end, we consider fractional derivatives in the senses of Riemann-Liouville and Caputo, and we analyze fractional forms of various integer-order differential operators, including the divergence, the rotational, the gradient, the Dirac and the Laplace operators. As the most important result, we prove an existence theorem for the solutions of a div-curl system, considering fractional differential operators of the Riemann-Liouville and Caputo types (see Theorems 2 and 3). Other important generalizations of well-known theorems from vector calculus are proven in this way. More precisely, we present fractional versions of the classical Divergence and Stokes Theorems for vector fields (see Proposition 2). Furthermore, we focus on the construction of fractional hyper-conjugate pairs, which represent a fractional generalization of the well-known Moisil-Teodorescu system in quaternionic analysis. Finally, we note that we are also able to provide an explicit expression for an inverse of the fractional gradient operator when we restrict ourselves to vector fields whose fractional rotational is zero, when we consider fractional derivatives of the Caputo type (see Proposition 8).

Data Availability Statement:
No new data were created or analyzed in this study. Data sharing is not applicable to this article.