An Introduction to Hyperholomorphic Spectral Theories and Fractional Powers of Vector Operators

The aim of this paper is to give an overview of the spectral theories associated with the notions of holomorphicity in dimension greater than one. A first natural extension is the theory of several complex variables whose Cauchy formula is used to define the holomorphic functional calculus for n-tuples of operators (A1,…,An)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$(A_1,\ldots ,A_n)$$\end{document}. A second way is to consider hyperholomorphic functions of quaternionic or paravector variables. In this case, by the Fueter-Sce-Qian mapping theorem, we have two different notions of hyperholomorphic functions that are called slice hyperholomorphic functions and monogenic functions. Slice hyperholomorphic functions generate the spectral theory based on the S-spectrum while monogenic functions induce the spectral theory based on the monogenic spectrum. There is also an interesting relation between the two hyperholomorphic spectral theories via the F-functional calculus. The two hyperholomorphic spectral theories have different and complementary applications. We finally discuss how to define the fractional Fourier’s law for nonhomogeneous materials using the spectral theory on the S-spectrum.


Introduction
The problem to define functions of an operator A or of an n-tuple of operators (A 1 , . . . , A n ) is very important both in mathematics and in physics and has been investigated with different methods starting from the beginning of the The first mathematicians who understood the importance of hypercomplex analysis to define functions of noncommuting operators on Banach spaces were A. McIntosh and his collaborators, staring from preliminary results in [74]. Using the theory of monogenic functions they developed the monogenic functional calculus and several of its applications, see [73]. The S-functional calculus, and in general the spectral theory on the S-spectrum, started its development only in 2006 when F. Colombo and I. Sabadini discovered the S-spectrum. The discovery of the S-spectrum and of the S-functional calculus is well explained in the introduction of the book [24] with a complete list of the references and it is also described how hypercomplex analysis methods were used to identify the appropriate notion of quaternionic spectrum whose existence was suggested by quaternionic quantum mechanics.
If we denote by B(V ) the Banach space of all bounded right linear operators acting on a two sided quaternionic Banach space V then the appropriate notion of quaternionic spectrum, the S-spectrum, is defined in a very counterintuitive way because it involves the square of the quaternionic linear operator T and it is defined as: The S-spectrum for quaternionic operators can be naturally defined also for paravector operators when we work in a Clifford algebra, see [45] and the book [48].
The problem of defining the quaternionic spectrum for the quaternionic spectral theorem has been an open problem for long time even though several attempts have been done by several authors in the past decades, see e.g. [85,87], however the correct definition of spectrum was not specified. Finally in [2] the spectral theorem on the S-spectrum was proved for both bounded and unbounded normal operators on a quaternionic Hilbert space.
The main problems with the quaternionic notion of spectrum can be described with the following considerations related to bounded linear operators just for the sake of simplicity. Let T : V → V be a right linear bounded quaternionic operator acting on a two sided quaternionic Banach space V . If we readapt the notion of spectrum for a complex linear operator to the quaternionic setting we obtain two different notions of spectra because of the noncommutativity of the quaternions. The left spectrum σ L (T ) of T is defined as where the notation sI in B(V ) means that (sI)(v) = sv. The right spectrum σ R (T ) of T is associated with the right eigenvalue problem, i.e., the search of those quaternions s such that there exists a nonzero vector v ∈ V satisfying T (v) = vs.
In both spectral problems it is unclear how to associate to the spectrum a resolvent operator with the property of being an hyperholomorphic function operator-valued. In fact, for the left spectrum σ L (T ) it is not clear what notion of hyperholomorphicity is associated to the map s → (sI − T ) −1 , for s ∈ H \ σ L (T ) and for the right spectrum it is even more weird because the operator Is − T (where Is means (Is)(v) = vs) is not linear, so it is not clear which operator is the candidate to be the resolvent operator.
Remark 1.1. One of the main motivations that suggested the existence of the S-spectrum is the paper [13] by G. Birkhoff and J. von Neumann, where they showed that quantum mechanics can be formulated also on quaternionic numbers. Since that time, several papers and books treated this topic, however it is interesting, and somewhat surprising, that an appropriate notion of spectrum for quaternionic linear operators was not present in the literature. Moreover, in quaternionic quantum mechanics the right spectrum σ R (T ) is the most useful notion of spectrum to study the bounded states of a quantum systems. Before 2006 only in one case the quaternionic spectral theorem was proved specifying the spectrum and it is the case of quaternionic normal matrices, see [53], where the right spectrum σ R (T ) has been used. Now we recall some research directions and applications of the hyperholomorphic function theories and related spectral theories.
The first step of FSQ-construction generates slice hyperholomorphic functions and the spectral theory of the S-spectrum, has the following research directions: • The foundation of the quaternionic spectral theory on the S-spectrum are organized in the books [23,24], and for paravector operators see [48]. • The mathematical tools for quaternionic quantum mechanics is the spectral theorem based on the S-spectrum [2,57]. • Quaternionic evolution operators, Phillips functional calculus, H ∞functional calculus, see [23]. • Quaternionic approximation [56].
• The theory of function spaces of slice hyperholomorphic functions [6]. • New classes of fractional diffusion problems based on fractional powers of quaternionic linear operators, see the book [23] and the more recent contributions [17][18][19]31,32].
In the last section of this paper we explain how to treat fractional diffusion problems using the quaternionic spectral theory on the S-spectrum and we show some of the results on fractional Fourier's law for nonhomogeneous materials recently obtained. An example of problems that we can treat is the following.
We warn the reader that in this paper, with an abuse of notation, we use the symbol x for both the coordinates of a point (x 1 , x 2 , . . . , x n ) ∈ R n or for the vector part of a quaternion or for the imaginary part of a paravector in a Clifford algebra.
Let Ω be a bounded or an unbounded domain in R 3 and let τ > 0 and denote by v the temperature of the material contained in Ω. Let x = (x 1 , x 2 , x 3 ) ∈ Ω and consider the evolution problem where f is a given datum and the heat flux for the nonhomogeneous material contained in Ω, is given by the vector differential operator: We determine the conditions on the coefficients a, b, c : Ω → R under which the operator T (x) generates the fractional powers P α (T (x)) of T (x), for α ∈ (0, 1). The vector part of P α (T (x)) of T (x) is defined to be the nonlocal Fourier's law associated with T (x). The second step of FSQ-construction generates Fueter or monogenic functions and the spectral theory on the monogenic spectrum. We highlight some references for the research directions in this area: • Monogenic spectral theory and applications [73]. Here one can also find the relations of the monogenic functional calculus with the Taylor functional calculus and the Weyl functional calculus see also some of the original contributions [68][69][70]72,79]. • Harmonic analysis in higher dimension, singular integrals and Fourier transform see the recent book [81]. • Algebraic Analysis of Dirac systems [44].
Prior to the recent developments on slice hyperholomorphic functions, this function theory was simply seen as an intermediate step in the Fueter-Sce-Qian construction. The literature on hyperholomorphic function theories and related spectral theories is nowadays very large. For the function theory  [6,46,48,56,60], while for the spectral theory on the S-spectrum we mention the books [5,23,24,48]. For the Fueter and monogenic function theory and related topics see the books [15,44,52,65,66,73,82].

Spectral Theory in the Complex Setting
In this section we discuss what is a functional calculus of a single operator on a Banach space and also for the case of several operators. When we consider a closed linear operator A with domain D(A) ⊂ X, where X is a Banach space, the resolvent set ρ(A) of A is defined as where B(X) in the space of all bounded linear operators on X and the spectrum of A is the set σ(A) = C\ρ(A) and for λ ∈ ρ(A) the map λ → (λI−A) −1 is called the resolvent operator. We start with the following intuitive definition of what is a functional calculus.
A functional calculus for a closed linear operator A on complex Banach space X is a mathematical technique that allows to construct in a meaningful way an operator f (A) for any function f in a certain class of functions F defined on sets that contain the spectrum σ(A) of A.
The formulation in a meaningful way usually means that the functional calculus is compatible with formally plugging A into the function f , whenever this is possible. That is, whenever f (z) can be expressed by a formula so that formally replacing z by the operator A yields an expression that is meaningful, then f (A) should correspond to this expression. Several examples below illustrate this idea: (a) For any polynomial p(z) = a n z n + · · · + a 1 z + a 0 ∈ F with a ∈ C, the operator p(A) should be given by p(A) = a n A n + · · · + a 1 A + a 0 I.
the functional calculus is compatible with the resolvent operator at λ. That is, we have R λ (A) = (λI − A) −1 . (c) For any rational function r(z) = p(z)/q(z) ∈ F with polynomials p(z) = a n z n + · · · + a 1 z + a 0 and q(z) = b n z n + · · · + b 1 z + b 0 with a , b ∈ C, the operator r(A) should be given by where p(A) = a n A n + · · · + a 1 A + a 0 I and (d) If A is the infinitesimal generator of a strongly continuous group U A (t), t ≥ 0 and exp(tz) ∈ F, then exp(tA) = U A (t).
Of course in the case of unbounded operators one has to pay attention to the domain of the operators. Usually the class F constitutes an algebra, often even a Banach algebra, and the meaningfulness of the functional calculus as described above follows from the compatibility of the functional calculus with the algebraic operation. Precisely, a functional calculus usually satisfies several (or all) of the following conditions: (I) The functional calculus is an algebra homomorphism that is (af + g)(T ) = af (T ) + g(T ) for all f, g ∈ F and all a ∈ C. (II) For f and g ∈ F such that fg ∈ F we expect We restrict ourselves to the case F consists of functions that are at least continuous and to the case that the topology on F is coarses than the topology of locally uniform convergence. In particular, convergence in F implies pointwise convergence.
For the measurable functional calculus the main statements of this section hold true, but their justification is based on different arguments.
There are two main methods for defining a functional calculus. There is also another concept behind the notion of functional calculus: the operator f (A) should be defined by letting f act on the spectral values of A. In particular, this means that This idea is usually not so much emphasized in the complex setting when one introduces and explains the concept of a functional calculus and we shall see here in the sequel the reason way this happens. However, it is this relation that explains why functional calculi are the fundamental techniques for investigating linear operators. If it does not hold, then a functional calculus does not provide any information about the operator even though it generates functions of operators. Interestingly enough a deep difference between the theory of complex and the theory of quaternionic linear operators (or more in general for hyperholomorphic spectral theories) is revealed here so that in the latter, the relation (1) needs to be addressed explicitly.
Let us start our considerations by justifying the importance of the relation (1). We therefore recall the easiest result that is shown by an application of a functional calculus. It is obvious that even for the above, very easy, and fundamental result, the fact that the polynomial functional calculus satisfies the relation (1) is crucial.
The fact that the relation (1) trivially holds true for any known functional calculus in the complex setting is shown in the next two results, and this is the reason for which it is not explicitly mentioned.  (1).
We now recall that the spectral theorem works as a functional calculus. In the finite dimensional case when we pick an n × n matrix A = (a i,j ) for i, j = 1, . . . n of complex numbers such that a i,j = a j,i for all i, j = 1, . . . , n. The spectrum σ(A) consists of eigenvalues of A, that is, complex numbers λ for which the equation Av = λv has a nonzero vector v ∈ C n as a solution. The hermitian matrix A has a unique decomposition as a finite sum where E λj ,A is the orthogonal projection onto the eigenspace of the eigenvalue λ j . In the case of bounded selfadjoint (or more in general normal operators) operators A acting in Hilbert space the spectral theorem is the most important tool for the complete description of such operators, in fact we have with respect to a spectral measure E λ;A associated with A. From the spectral theorem we can define f (A) by is valid for the joint spectral measure E λ;A1,...,An associated with A 1 , . . . , A n . The joint spectrum of A 1 , . . . , A n in R n is the support of E λ;A1,...,An and f : σ(A 1 , . . . , A n ) → C is any bounded Borel measurable function. The theorem holds more in general for unbounded normal operators but one has to pay attention to the definition of commutativity in this case, see the book [84].
In the case we work in a Banach space the most natural way to define functions of bounded (and also of unbounded) operators is the Riesz-Dunford functional calculus which holds for all holomorphic functions f defined in a neighborhood of σ(A) in the complex plane. The simple closed contour C surrounds σ(A) and is contained in the domain of the function f . There is a natural generalization of Riesz-Dunford functional calculus for n-tuples of bounded operators where dλ = dλ 1 · · · dλ n and f is any holomorphic function in a neighborhood of σ(A 1 ) × · · · × σ(A n ) in C n . For each j = 1, . . . , n the simple closed contour C j surrounds σ(A j ) and C 1 × · · · × C n is contained in the domain of f in C n . Also when the operators A 1 , . . . , A n do not commute with each other, the functional calculus makes sense with any change in the operator ordering of the function . For noncommuting operators, the results are in general more complicated and we will not enter into the details here (see, e.g., the book [73]). The material discussed in this section can be found is several classical books, such as [86].
Some remarks in view of the hyperholomorphic spectral theories.
(I) In operator theory on the S-spectrum, the statements of Theorems 2.4,2.5 still hold true but just for a subclass of functions. The conditions that the functions in the dense subspace F 0 resp. the kernel K(λ, A) satisfy (1) is not true in this setting. It is neither satisfied by the S-resolvent operator, nor by the F -resolvent operator, nor by slice hyperholomorphic rational functions with non-real coefficients. In particular it is not satisfied, whenever the leftlinear structure of the space has a prominent role.
(II) In general for the definitions of hyperholomorphic functional calculi (S-functional calculus, the F -functional calculus, the monogenic functional calculus) the product rule, the composition rule or the spectral mapping theorem do not hold. One needs additional arguments to show that functional calculi based on the S-spectrum satisfy (1) at least for a subclass of functions namely, the class of intrinsic functions and for these functions, the problems mentioned before do not occur.

Spectral Theories in the Hyperholomorphic Setting
At the beginning of the last century several authors started the study of hyperholomorphic functions and the most popular classes of functions are nowadays called Fueter (or Cauchy-Fueter) regular functions in the case of the quaternions and monogenic functions (or functions in the kernel of the Dirac operator) for Clifford algebra setting. The second class of hyperholomorphic functions have been developed more recently, just at the beginning of this century, and different definitions are possible even though they are not totally equivalent.
In the following we will discuss mainly the implications of the Fueter-Sce-Qian construction in the Clifford setting, the quaternionic setting is similar and we will use it for the fractional powers of vector operators in the last section of this paper.
Let R n be the real Clifford algebra over n imaginary units e 1 , . . . , e n satisfying the relations e e m + e m e = 0, = m, e 2 = −1. An element in the Clifford algebra will be denoted by x j e j ∈ R n called paravector and the real part x 0 of x will also be denoted by Re(x). The imaginary part of x is defined by Im(x) = x 1 e 1 + · · · + x n e n and for the sake of simplicity we also use the notation x for Im(x). The conjugate of x is denoted by x = x 0 − Im(x) and the Euclidean modulus of x is given by |x| 2 = x 2 0 + · · · + x 2 n . The sphere of purely imaginary paravectors with modulus 1, is defined by S = {x = e 1 x 1 + · · · + e n x n | x 2 1 + · · · + x 2 n = 1}. The element I ∈ S are such that I 2 = −1 so we will denote the complex space with imaginary unit I by C I . For this reason the elements of S are also called we can associate to it the sphere defined by The set of quaternions will be denoted by H and the above definitions adapt to this setting in a natural way.
We say that U is axially symmetric if, for every u + Iv ∈ U , all the elements u + Jv for J ∈ S are contained in U .
For operator theory the most appropriate definition of slice hyperholomorphic functions is the one that comes from the Fueter-Sce-Qian mapping theorem because it allows to define functions on axially symmetric open sets.

Definition 3.2. Let U ⊆ R n+1 be an axially symmetric open set and let
and if f 0 and f 1 satisfy the Cauchy-Riemann system The above definition adapts naturally to the quaternionic setting. Since we will restrict just to left slice hyperholomorphic function on U we introduce the symbol SH L (U ) to denote them. The subset of intrinsic functions consist of those slice hyperholomophic functions such that f 0 , f 1 are real-valued and is denoted by N (U ). We recall that right slice hyperholomorphic functions are of the form where f 0 , f 1 satisfy the above conditions.
where D is the Dirac operator defined by The definition of slice hyperholomorphic functions and of monogenic functions can be seen as two steps in the Fueter-Sce-Qian constructions to extend holomorphic functions to dimension greater than one for the vectorvalued functions (quaternionic or Clifford valued-functions).
In fact, starting from holomorphic functions, R. Fueter in 1935, see [55], showed an interesting way to generate Cauchy-Fueter regular functions. More then 20 years later in 1957 M. Sce, see [83], extended this result in a very pioneering and general way that includes Clifford algebras, see the English translation of his works in hypercomplex analysis with commentaries collected in the recent book [33].
In the original construction of R. Fueter the holomorphic functions are defined on open sets of the complex upper half plane. This condition can be relaxed by taking function defined in a set D ⊆ C, symmetric with respect to the real axis such that namely if g 0 and g 1 are, respectively, even and odd functions in the variable v. Additionally the pair (g 0 , g 1 ) satisfies the Cauchy-Riemann system. The above remark holds also for M. Sce's theorem that we state in the following for Clifford algebras.

be a holomorphic function defined in a domain (open and connected) D in the upper-half complex plane and let
be the open set induced by D in R n+1 . The following map takes the holomorphic functionsf (z) and induces the Clifford-valued function f (x). Then the function where T F S2 := Δ n−1 2 n+1 and Δ n+1 is the Laplacian in n + 1 dimensions, is in the kernel of the Dirac operator, i.e., The case where the operator Δ n−1 2 n+1 has a fractional index has been treated by T. Qian in [78]. Observe that for the Fueter's theorem the operator T F S2 is equal to the Laplacian Δ in 4 dimensions. Further developments can be found in [75][76][77] see also the survey [80].
where T F S1 denotes the first linear operator of the Fueter-Sce construction and T F S2 the second one. The Fueter-Sce mapping theorem induces two spectral theories according to the two classes of hyperholomorphic functions it generates.
Recently also the problem of construction the inversion of the maps that appear in the Fueter-Sce-Qian extension has been treated, we mention the papers [11,12,34,42,43], while a different method to connect slise monogenic and monogenic functions is via the Radon and dual Radon transform, see [30].
Remark 3.5. The theory of slice hyperholomorphic functions was somewhat abandoned until 2006 when G. Gentili and D. C. Struppa (inspired by C. G. Cullen [51]) introduced in [59] the notion of slice regular functions for the quaternions. Further developments of the theory of slice regular functions were discussed also in [28] and the above definition was extended by F. Colombo, I. Sabadini and D.C. Struppa, in [47], (see also [35,49,50]) to the Clifford algebra setting. Slice regular functions as defined in [59] and their generalization to the Clifford algebra as in [47], called slice monogenic functions, possess good properties on specific open sets that are called axially symmetric slice domains. When it is not necessary to distinguish between the quaternionic case and the Clifford algebra case we call these functions slice hyperholomorphic. The extension of the notion of a slice hyperholomorphic function on real alternative algebras can be found in [62].
Remark 3.6. It is also possible to define slice hyperholomorphic functions, as functions in the kernel of the first order linear differential operator (introduced in [29]) While, a fourth way to introduce slice hyperholomorphicity, done in 1998 by G. Laville and I. Ramadanoff in the paper [67], is inspired by the Fueter-Sce-Qian mapping theorem. They introduce the so called Holomorphic Cliffordian functions defined by the differential equation DΔ m f = 0 over R 2m+1 , where D is the Dirac operator. Observe that the definition via the global operator G requires less regularity of the functions than the definition in [67].
We now recall the hyperholomorphic Cauchy formulas that are the heart of the hyperholomorphic spectral theories. It is important to remark that the hypotheses of the following Cauchy formula are related to the Definition 3.2 of slice hyperholomorphic functions.
where the slice hyperholomorphic Cauchy kernel is given by and the value of the integral (2) depends neither on U nor on the imaginary unit I ∈ S. Moreover, the Cauchy kernel can be written in two ways as follows: Remark 3.8. The notion of slice regularity was introduced in the paper [59], besides the definition, the authors treated power series centered at the origin and some consequences. Without any tools the Cauchy formula with slice hyperholomorphic kernel and the representation formula were originally determined via the following elementary considerations. To determine the slice hyperholomorphic Cauchy kernel S −1 L (s, q) we observe that from the definition of slice regularity its expansion is true when q and s belong to the same complex plane C I , for I ∈ S. Then we ask ourselves what is the closed form of the series in the case q and s do not belong to the same complex plane C I observing that is true also when s and q do not belong to the same complex plane C I . In the quaternionic case it was observed that the inverse S of S −1 L (s, q) is the non trivial solution of the quaternionic equation which easily follows from (5). The unknown S was determined using the Niven's Algorithm as it is shown in the historical Note 4.18.3 in the book [48] and it gives Taking the inverse of S(s, q) we have the Cauchy kernel defined in (3) for the quaternions. This strategy to determine the Cauchy kernel shows that in the Clifford setting the Cauchy kernel remains the same if we consider ..+x n e n and s = x 0 +s 1 e 1 +...+s n e n are paravectors. Moreover, observe that a direct computation of the integral by computing the residues of the singularities of the kernel S −1 L (s, x) in the complex plane C I , gives: choosing any J ∈ S for all x = u + Iv ∈ U . From here one can clearly see the existence of the structure formula (or representation formula) for slice monogenic functions, see for example [35] (or [36]): The quaternionic setting is just a particular case and from these observations started a full development of the theory of slice hyperholomorphic functions.
The second ingredient for our discussion below is the Cauchy formula for monogenic functions.
for every x in U , where the monogenic Cauchy kernel is given by Before to introduce the basic fact on the hyperholomorphic spectral theories we need some important considerations.
(I) Holomorphic functions of one complex variable and harmonic analysis are strongly connected since the Cauchy-Riemann operator factorizes the Laplace operator. The holomorphic functional calculus and the spectral theorem are based on the same notion of spectrum.
(II) In order to restore the analogy with the holomorphic functional calculus and the spectral theorem in the complex setting we have to replace the classical spectrum with the S-spectrum. In fact the S-functional calculus The monogenic functional calculus, based on the monogenic spectrum, has applications in harmonic analysis and in other related fields. Let us consider a Banach space V over R with norm · . It is possible to endow V with an operation of multiplication by elements of R n which gives a two-sided module over R n and by V n we indicate the two-sided Banach module over R n given by V ⊗ R n .
We start with the definition of a functional calculus for (n + 1)-tuples of not necessarily commuting operators using slice hyperholomorphic functions. So we consider the paravector operator  even though the components of T do not commute. From this relation we justify the definition of the S-resolvent operator and of the S-spectrum. The quaternionic setting is just a particular case.
We have the following definition. Definition 3.11. (S-spectrum) Let T ∈ B(V n ) be a paravector operator. We define the S-spectrum σ S (T ) of T as: where I denotes the identity operator. The S-resolvent set of T is defined as A similar definition can be given for the right resolvent operator. A crucial result for the definition of the S-functional calculus is that depends neither on U nor on the imaginary unit I ∈ S, so the S-functional calculus turns out to be well defined.
Definition 3.14. (S-functional calculus) Let T ∈ B(V n ) and let U ⊂ R n+1 be as above. We set ds I = −Ids and we define the S-functional calculus as for f ∈ SH L σS (T ) . Observe that the definition of the S-functional calculus is very natural for non commuting operators in noncommutative spectral theory. The heart of the general version of the S-functional calculus can be found in the original papers [1,36,39,40] and its commutative version [38].
Warning. In the monogenic setting the natural functional calculus is for vector operators that is when we set T 0 = 0 in the paravector operator T = T 0 + n j=1 e j T j . The reason will be clear in the sequel, but to point out this fact we use the symbol A = (A 1 , . . . , A n ) or A = n j=1 e j A j instead of (T 1 , . . . , T n ) or T = n j=1 e j T j . Using the Cauchy integral formula for monogenic functions, we establish the monogenic functional calculus for the n-tuple A = (A 1 , . . . , A n ) of bounded linear operators on a Banach space X by substituting the n-tuple A for the vector x ∈ R n .
In the following for the monogenic functional calculus we limit ourselves to the most simple case when n is odd and the n-tuple A = (A 1 , . . . , A n ) of bounded linear operators commute among themselves. Such restrictions can be removed but one needs to do further considerations.

Remark 3.15.
If n is odd, A is a commutative n-tuple, that is, A j A k = A k A j for j, k = 1, . . . , n, and each operator A j has real spectrum σ(A j ) ⊂ R for j = 1, . . . , n, then for suitable ω ∈ R n+1 , the expression makes sense as an element of B(V n ) and it is called the monogenic resolvent. for ω = ω 0 + n j=1 ω j . Observe that the operator is called the monogenic spectrum.
where G ω (A) is the monogenic resolvent operator (9), η(ω) is the outer unit normal to ∂U and dS(ω) is the scalar element of surface area on ∂U .
In the case m = 2, 4, 6, ... the operator |ωI −A| −m needs to be defined in a suitable way. The direct formulation employs Taylor's functional calculus, but by using the plane wave decomposition of the Cauchy kernel, the case of even n and noncommuting operators can be treated simultaneously. For an n-tuple (A 1 , . . . , A n ) of commuting bounded linear operators on a Banach space V with real spectra, the nonempty compact subset γ(A) of R n coincides with Taylor's joint spectrum defined in terms of the Koszul complex.
The Cauchy formula of slice hyperholomorphic functions allows to define the notion of S-spectrum, while the Cauchy formula for monogenic functions induces the notion of monogenic spectrum, as illustrated by the diagram: In the above diagram we have replaced the set of intrinsic functions N by the larger set of slice hyperholomorphic functions SH. This is clearly possible because the map T F S2 is the Laplace operator or its powers. We finally recall that the quaternionic spectral theorem is based on the S-spectrum and not on the monogenic spectrum. In 2015 (and published in 2016) the quaternionic spectral theorem for quaternionic normal operators was finally proved, see [2] (the case of unitary operators is treated in [3] and the case of compact normal operators is in [61]). Later on perturbation results of quaternionic normal operators were proved in [16]. Beyond the spectral theorem there are more recent developments in the direction of the characteristic operator functions, see [5] and the theory of quaternionic spectral operators was developed in [58].
Finally, we wish to give an idea of the structure of the quaternionic spectral theorem. For a complete treatment see [24]. If T ∈ B(H) is a bounded normal quaternionic linear operator, on a quaternionic Hilbert space H, then there exist three quaternionic linear operators A, J, B such that T = A + JB, where A is self-adjoint and B is positive, J is an anti self-adjoint partial isometry (called imaginary operator). Moreover, A, B and J mutually commute.
There exists a unique spectral measure E I on σ S (T ) ∩ C + I so that for any slice continuous intrinsic function f = f 0 + f 1 I we have: where x, y ∈ H. This theorem extends to the case of unbounded operators as well and holds true for a larger class of functions that are not necessarily continuous. The continuous functional calculus is treated in [63].

Intersection of the Hyperholomorphic Spectral Theories
Now we formulate the Fueter-Sce-Qian theorem in integral form and we use it to define the F -functional calculus. This gives a version of the monogenic functional calculus for n-tuples of commuting operators but it is based on the S-spectrum instead of the monogenic spectrum. This calculus was introduced in [41] and further investigated in [21,37]. It is important to recall that the monogenic functional calculus is defined for n-tuples of operators A j , j = 1, . . . , n that have real spectrum considered as operators A j : V → V on the real Banach space V . The F -functional calculus has advantages and disadvantages with respect to the monogenic functional calculus. Precisely, the F -functional calculus allows to consider a much larger class of operators because it does not require that the spectrum of the operators A j , j = 1, . . . , n has to be real. Moreover, this calculus allows to consider paravector operators and not only vector operators as the monogenic functional calculus imposes. On the other hand, from the hyperholomorphic functions point of view the F -functional calculus is less general with respect to the monogenic functional calculus because it works for the subset of monogenic function given byM We now show how the Fueter-Sce mapping theorem provides an alternative way to define the functional calculus for monogenic functions. The main idea is to apply the Fueter-Sce operator T F S2 to the slice hyperholomorphic Cauchy kernel as illustrated by the diagram: Observe that in the above diagram the arrow from the space of axially monogenic function AM (U ) is missing because the F -functional calculus is deduced from the slice hyperholomorphic Cauchy formula.
This method generates an integral transform, called the Fueter-Sce mapping theorem in integral form, that allows to define the so called F -functional calculus. This calculus uses slice hyperholomorphic functions and the commutative version of the S-spectrum and now we show how it works. We point out that the operator T F S2 has a kernel and one has to pay attention to this fact with the definition of the F -functional calculus, more details are given in [24]. Now observe that one can apply the powers of the Laplace operators to both sides of (2) so that we have In general, it is not easy to compute Δ h f and when we apply Δ h to the Cauchy kernel written in the form (3), we do not get a simple formula. However, S −1 L (s, x) can be written in two equivalent ways as follows. If we use the second expression for the Cauchy kernel (see formula (4)) we find a very simple expression for Δ h S −1 L (s, x). In fact, we have:  (x 0 , x 1 , ..., x n ). Then, for h ≥ 1, we have: where C n,h : The function Δ h S −1 (s, x) is slice hyperholomorphic in s for any h ∈ N but is monogenic in x if and only if h = (n + 1)/2, namely if and only if h equals the Sce's exponent. We define the kernel which can be used to obtain the Fueter-Sce mapping theorem in integral form.
is monogenic and it admits the integral representation where the integral depends neither on U nor on the imaginary unit I ∈ S.
In the sequel, we will consider bounded paravector operators T , with commuting components T ∈ B(V ) for = 0, 1, . . . , n. Such subset of B(V n ) will be denoted by BC 0,1 (V n ). BC 0,1 (V n ) The F -functional calculus is based on the commutative version of the S-spectrum given by where the operator T is defined by We observe that for historical reasons the commutative version of the Sspectrum is sometimes called F -spectrum because it is used for the F -functional calculus. So we define the F -resolvent operators.
where the constants γ n are given in (12).

Definition 4.5. (T heF -functional calculus for bounded operators)
Let n be an odd number, let T = T 0 + T 1 e 1 + · · · + T n e n ∈ BC 0,1 (V n ) and set ds I = ds/I, for I ∈ S. Let SH L σS (T ) and U be as in Definition 3. 13. We definȇ The definition of the F -functional calculus is well posed since the integral in (15) depends neither on U and nor on the imaginary unit I ∈ S.
We conclude this section with several remarks on the hyperholomorphic functional calculi to stress the difference with respect to the complex case.
(I) The product rule holds for the S-functional calculus but only in the case one of the two functions is intrinsic function. For the monogenic functional calculus the product rule does not hold. This is due to the fact that the product of two monogenic functions is not monogenic. For the Ffunctional calculus the product rule does not hold.
(II) Regarding the compatibility with polynomials we have that the S-functional calculus and the monogenic functional calculus are compatible with slice hyperholomorphic polynomials and with monogenic polynomials, respectively. For the F -functional calculus the compatibility with polynomials holds if we considerP (q) = ΔP (q) whereP (q) is a monogenic (or Fueter) and P is a slice monogenic polynomial q → T ⇒P (q) ⇒P (T ) (III) The spectral properties of the operator T can be deduced by the S-functional calculus and the quaternionic spectral theorem for which when we use intrinsic functions.

The S-Spectrum Approach to Fractional Diffusion Problems
An important extension of the S-functional calculus to unbounded sectorial operators is the H ∞ -functional calculus which is one of the ways to define functions of unbounded operators. The H ∞ -functional calculus has been used to define fractional powers of paravector operators and of quaternionic linear operators that define fractional Fourier laws for nonhomogeneous material in the theory of heat propagation. For the original contributions on fractional powers of vector operators and of quaternionic operators and of the H ∞functional calculus based on the S-spectrum see [7,20,22]. For a systematic and recent treatment of quaternionic spectral theory on the S-spectrum and the fractional diffusion problems based on techniques on the S-spectrum see the books [23,24] published in 2019. Moreover, in the monograph [48], published 2011, one can find also the foundations of the spectral theory on the S-spectrum for n-tuples of noncommuting operators. The theory on the fractional powers of quaternionic operators has been recently applied to physical problems and in particular to generate the fractional Fourier law for the heat equation that is collected in the papers [17][18][19]31,32], here we give an overview of some of our results.
We denote by x := (x 1 , x 2 , x 3 ) a generic point in R 3 . Let Ω ⊂ R 3 bounded or unbounded domain (with C 1 boundary), the heat equation for nonhomogeneous materials with the associated initial-boundary conditions descibes the evolution of the heat. Precisely, we determine v : where f is a given datum and where we suppose that the coefficients a 1 , a 2 , a 3 : Ω ⊂ R 3 → R of T belong to C 1 (Ω) and they are not necessarily constant. We also consider the heat equation for nonhomogeneous materials with Robin boundary conditions, that consists in finding v : where n = (n 1 , n 2 , n 3 ) is the outward unit normal vector to ∂Ω, and b : ∂Ω → R is a given continuous function. From the physical point of view, if we call q(x, t) the flux of the quantity described by v(x, t) at the instant t, the Fourier's law states that q(x, t) = T (x)(v(x, t)). The simpler case is when we consider Ω = R 3 and the homogeneous diffusion problem is the consequence of Fourier's law and of the conservation of the energy In this case T is reduced to the negative gradient operator T = −∇ and observing that div •∇ = Δ the fractional diffusion model is obtained by replacing in the heat equation the Laplace operator by its fractional powers The fractional versions of the evolution equation in R 3 is given by We observe that the fractional diffusion problem modifies both the Fourier's law and the conservation of the energy. Using the quaternionic functional calculus we are able to define the fractional powers of vector operators, such as ∇ or T , in a bounded or unbounded domain Ω of R 3 . Denoting, just for the moment, by T α or ∇ α these fractional operators, we can define the fractional diffusion problem (20) in the divergence form The boundary conditions to associate with the fractional evolution problem are a very delicate issue and will not discussed here. We just mention that the most natural boundary conditions are v = 0 at infinity in the case Ω = R 3 .
Remark 5.1. The boundary condition that we have to assume to generate the fractional powers of T are given by where a : ∂Ω → R is a given continuous function. These Robin-like boundary condition differs from the boundary condition in (19) by a power two on the coefficients a 's. This is due to the fact that in (19) the boundary condition rise from a physical condition on the flux through the boundary, instead the boundary condition in (22) naturally comes from the definition of T α . In any case, the two boundary conditions are related when T has coefficients that become constant on the boundary ∂Ω. Indeed, suppose that there exists a constant μ such that the functions a 1 , a 2 , a 3 satisfy the conditions and the coefficients a and b are such that Then the relation For, using (23) and (24), we have This approach has several advantages. • It generates the fractional Fourier law from the Fourier law using the boundary conditions of the problem and without modifying the conservation of energy law. • We can define the fractional heat equation for nonhomogeneous materials.
• The fractional differential equation remains in the divergence form, so the definition of a weak solution in obtained in a simple way. • It turns out that the approach through the quaternionic functional calculus for defining the fractional heat equations is consistent with the classical one. Indeed, we have that for any α ∈ (0, 1) Now we present how the quaternions can be used to describe the vector operators. Let e , for = 1, 2, 3, be an orthogonal basis for the quaternions H. We identify the vector operator T , described in (18), with the quaternionic gradient operator with non constant coefficients ⎛ ⎜ ⎝ where the components T , = 1, 2, 3, are defined by T := a (x)∂ x , x ∈ Ω. From the physical point of view the operator T , defined in (26), can represent the Fourier law for nonhomogeneous materials, but it can represent also different physical laws. Our goal is to generate the fractional powers of T , that we denote with P α (T ) for α ∈ (0, 1), when the operators T , for = 1, 2, 3 do not commute among themselves.
Remark 5.2. The notation P α (T ), for the fractional powers of T , is more precise with respect to the formal notation T α in a sense that will be clear just in the following with the precise definition. In simple words P α (T ) is defined by a projector when there are points of the S-spectrum that belong to the negative real line (see for example (30)). The formal notation T α is used in (21) and (22) (see (25)). The operator Q s (T ) −1 is called the pseudo S-resolvent operator. For s ∈ ρ S (T ), the left S-resolvent operator is defined as (28) and the right S-resolvent operator is given by The fractional powers of T , denoted by P α (T ), are defined as follows: for any I ∈ S, for α ∈ (0, 1) and v ∈ dom(T ) we set or where ds j = ds/I. These formulas are a consequence of the quaternionic version of the H ∞ -functional calculus based on the S-spectrum, see the book [23] for more details. For the generation of the fractional powers P α (T ) a crucial assumption on the S-resolvent operators is that, for s ∈ H \ {0} with Re(s) = 0, the estimates hold with a constant Θ > 0 that does not depend on the quaternion s. It is important to observe that the conditions (32) assure that the integrals (30) and (31) are convergent and so the fractional powers are well defined. For the definition of the fractional powers of the operator T we can use equivalently the integral representation in (30) or the one in (31). Moreover, they correspond to a modified version of Balakrishnan's formula that takes only spectral points with positive real part into account.
We want to apply the previous theory to the case V := L 2 (Ω, H) and T ∈ K(L 2 (Ω, H)) defined as in (26) (dom(T ) ⊂ L 2 (Ω, H) is a densely subset). A crucial problem is to determine the conditions on the coefficients a 1 , a 2 , a 3 : Ω ⊂ R 3 → R such that (30) and (31) are convergent. This problem is split into two problems: • The first is to find appropriate conditions for the coefficients a i 's such that the purely imaginary quaternions are in the S-resolvent set ρ S (T ) (i.e. Q s (T ) : dom(T 2 ) → L 2 (Ω, H) is invertible and bounded). This is a necessary condition, see formulas (30) and (31). Then, observe that in the quaternionic case the map s → s α , for α ∈ (0, 1) is not defined for s ∈ (−∞, 0). For this reason it is of great importance to assume the condition Re(s) ≥ 0 that avoids the half real line (−∞, 0]. • The second crucial fact is to determine the conditions on the coefficients a i 's such that the estimate (32) for the S-resolvent operator of T holds true.
Both these problems are solved by considering the following approach. According to the initial condition of the boundary-value problems we invert the operator Q s (T ) on the space H 1 0 (Ω, H), when we consider the Dirichlet boundary condition, and on the space when we consider the Robin boundary condition. The invertibility of Q s (T ) is thus reduced to solve in a weak sense the following two partial differential equations: given F ∈ L 2 (Ω, H) and s ∈ H \ {0} such that Re(s) = 0 and To solve in the weak sense (33) (resp. (34)) means that for any F ∈ L 2 (Ω, H) we have to find u F ∈ H 1 0 (Ω, H) (resp. u F ∈ H(Ω, H)) such that: for where the angle-brackets means that Q s (T ) is applied to u F in the sense of distribution. If we solve ( (33): or the Robin boundary condition for (34) In [17][18][19]31,32] we found suitable conditions for the coefficients a i 's such that the two sesquilinear forms b s,1 and b s,2 are coercive and continuous when Ω is bounded or unbounded. In conclusion by the Lax-Milgram lemma we obtain the solvability of the equation (35) (i.e. the invertibility of Q s (t)) and the estimate (32) for the S-resolvent operator. Thus the problem of the convergence of (30) and (31) is solved. In the next three paragraphs we summarize the conditions we found on the coefficients of T to obtain the convergence of (30) and (31) .
The fractional Fourier's law in the problem (21) with Ω bounded. In the paper [19] it was considered the commutative Fourier's law T com , that is an operator of the form where the real operators a 1 (x)∂ x1 , a 2 (x 2 )∂ x2 and a 3 (x 3 )∂ x3 commute among themselves. It has been shown that if the coefficients a : Ω → R, for = 1, 2, 3 belong to C 1 (Ω, R) and if a , for = 1, 2, 3 are suitably large and their derivative are suitably small then the integrals in (30) and (31) are convergent.
In [31] we replace the commutative Fourier's law T com by the more general Fourier's law where now the real operators a 1 (x)∂ x1 , a 2 (x)∂ x2 e 2 and a 3 (x)∂ x3 do not commute among themselves. In this case the conditions for the existence of the fractional powers are more complicated.
The main result is summarized in the following theorem (see for more details Theorems 4.1, 4.4 and 4.5 in [31]).
The sesquilinear form b s,1 (u, v), defined in (36), associated with the invertibility of the operator Q s (T ) := T 2 − 2s 0 T + |s| 2 I, with homogeneous Dirichlet boundary conditions, has to be considered with care. The next Remark 5.5 summarizes several facts concerning b s,1 (u, v).  (IV) If we consider the quadratic form in dimension n = 3, that is when Ω is a C 1 bounded set in R 3 and a 1 = 0, a 2 = 0, a 3 = 0, then the quadratic form is not coercive because of a 3 = 0. It seems that this case cannot be treated using Lax-Milgram Lemma, but a suitable method for degenerate equations has to be used.
(V) In [31] we proved, under more restrictive hypothesis on the coefficients a i 's, that the sesquilinear form b s,1 (u, v) is continuous and coercive in H 1 0 (Ω, H).
(VI) From the physical point of view the case for a 1 = 0, a 2 = 0, a 3 = 0, in dimension n = 3 is the case in which the conductivity is the direction z goes to zero.
(VII) The proofs for the continuity and coercivity are similar in any dimension, the estimate for the S-resolvent operators have some differences according to the fact that we work in Y or in X .
The fractional Fourier's law in the problem (22) with Ω bounded. Regarding the initial boundary value problem of Robin-type, the conditions on the coefficients a i 's depend on two positive constants which appear in the following two inequalities: • for all u ∈ H 1 (Ω, R) the following inequality holds: where C ∂Ω does not depend on u • for all u ∈ H 1 (Ω, R) the following inequality holds: where C P does not depend on u.
Theorem 5.6. Let Ω be a bounded domain in R 3 with boundary ∂Ω of class C 1 . Assume that a ∈ C 0 (∂Ω, R) and let T be the operator defined in (18) with coefficients a 1 , a 2 , a 3 ∈ C 1 (Ω, R). Define the following constants: where · ∞ denotes the sup norm and C ∂Ω is the constant in (42). Moreover, assume that C T − C T C P − K a,Ω 1 + C 2 P > 0 and C T > 0, (44) where C P is the constant in (43). Then for any α ∈ (0, 1) and for any v ∈ dom(T ), the integrals (30) and (31) converge absolutely. Remark 5.7. In [17] we treated the case of the operator T with commutative coefficients and of Ω bounded with Robin-type boundary condition. The fractional Fourier's law in the problem (21) with Ω unbounded. Regarding the initial boundary value problem of Dirichlet-type for the unbounded domains, we obtained in [18] the following result (see Theorems 4.8, 5.1 and 5.2 in [18]). Theorem 5.8. Let Ω be an unbounded domain in R 3 with boundary ∂Ω of class C 1 . Let T be the operator defined in (18) with coefficients a 1 , a 2 , a 3 ∈ C 1 (Ω, R) ∩ L ∞ (Ω, R). Suppose that and Then for any α ∈ (0, 1) and for any v ∈ dom(T ), the integrals (30) and (31) converge absolutely.

Concluding Remarks
During the pandemic period several efforts have been done by D.P Kimsey and I. Sabadini with the authors for further developments of the spectral theory on the S-spectrum and its applications. During the year 2020 part of this material has been collected in some papers. Precisely in [26] it was proved the spectral theorem for normal operators on a Clifford module, this paper has motivated the theory of slice monogenic functions of a Clifford variable introduced in [27] and the universality property of the S-functional calculus considered in [25].