Basic Properties of Solutions of Boundary Value Problem for the Dirac-Wave Operator

We study local and global properties of solutions of the boundary value problem for the Dirac wave operator. This includes explicit formulas, domain of dependence, range of influence, finite propagation speed and energy estimates. We also introduce the space of Dirac wave homogeneous polynomials and the corresponding Dirac-wave conjugate functions.


Introduction
As it is well know from a basic course in PDE, solutions of the boundary value problem (BVP) for the wave operator or d'Alembertian : = ∂ 2 t − = 0 have as main local properties explicit formulas, domain of dependence, range of influence, finite propagation speed and as global properties energy estimates among others.On the other hand, we also recall that the fundamental solution of the wave operator is not C ∞ and this has important implication for the regularity of their solutions.In particular, a singularity at a point propagates along the characteristic cone and focuses at its vertex.
The Dirac-wave operator is a first order linear operator with coefficients in a clifford algebra which factor the wave operator.We will present analogues properties as the ones mentioned above for the wave operator for solutions of the BVP for the Diracwave operator.This is a first step to have a kind of complex analysis for the wave operator.We also will present the corresponding wave conjugate functions and its Hilbert transform and the space of homogeneous polynomials which satisfy the Diracwave operator.
Previous works from David Eelbode and Frank Sommen compute the fundamental solution for the Dirac-wave operator on hyperbolic spaces [1-3] see also Rosen [4] Chapter 9.As we said we will be interested in local and global properties of solutions.We point out that the behavior of solutions of the hyperbolic operator are rather different of solutions of the elliptic operator and the same will applies to solutions of the Dirac operator and the Dirac-wave operator.
As an introduction and for future reference let us recall the BVP for the wave operator and its solution.We will follow [5][6][7].Recall that the BVP for the wave operator is to find u ∈ C 2 (0, ∞) × R n such that given the initial value φ(x) and the initial velocity ψ(x).
We state the solution to the BVP in dimension three first and comment about the other cases below.
) and satisfy the wave operator:

Moreover, for any x
Note that the solution to the BVP u is C k -continuous up to the boundary.Given k ≥ 2 and any ϕ(x) ∈ C k+1 (R 3 ) it is defined its Riemann transform on (0, ∞) × R 3 as (1.4) We will need that following identity which is easy to see It follows that Kirchhoff formula could be written as Also note that Rψ solves the BVP with zero initial values and with initial velocity ψ and ∂ t (Rφ) solves the BVP with initial values φ and with zero initial velocity.Now a few comments on the other dimensions.When n > 3 is odd the explicit formula for the solution of the BVP need φ with k ≥ 2 and the corresponding Riemann transform of ϕ is given by where and ω n is the surface area of the unit sphere in R n .Then the solution of the BVP is given in an analogous way as in (1.6).The solution for the case when n is even is derived from the solution for n + 1 by the method of decent.The corresponding transform is given by dy . (1.8) Again the solution of the BVP in the even case is given in an analogous way as in (1.6).This explicit formulas (1.7), (1.8) demonstrate the necessity of making more and more smoothness assumption on the data to ensure the existence of a C 2 -solution.The key point here is we can work directly with C ∞ -functions as data.For the distributional case we recommended [8].The layout of this note is as follows.In Sect. 2 we introduce the Clifford algebra and the Dirac-wave operator and some of their properties.We also introduce the Teodorescu transform that we will use later.In Sect. 3 we study the BVP for the Dirac-wave operator.We give several equivalences of their solutions.Here we see the domain of dependences and the range of influences of the wave-monogenic functions.As an application we find a basic for the space of Dirac-wave homogeneous polynomials.In Sect. 4 we study the energy estimates for the Dirac-wave operator.In particular, we see here uniqueness for the Dirac wave equation in a bounded set and finite propagation speed.In the last Sect.5, we see how to extend a set of wave functions to a monogenic-wave function and from this we define the corresponding Dirac-wave conjugate function.It turns out that the Teodorescu transform corresponds to the wave Hilbert transform.We also decompose any solution of the wave BVP with values in the Clifford algebra as a sum of a Dirac-wave BVP solution and an anti Dirac-wave BVP solution.

Clifford Algebras and Dirac Operators
Here we collect basic properties of the Clifford algebra with will be working with.For more background material and further general references on Clifford algebras and related matters, the interested reader is referred to the monographs [9][10][11].See also [4].
We start with the Minkowski space R 1,n , this is R n+1 = (t, x 1 , x 2 , . . ., x n ) with the quadratic form and with the standard orthonormal basis { , e 1 , e 2 , . . ., e n }.The Clifford algebra C (1, n) with signature (1, n) is the minimum enlargement of R 1,n to a unitary algebra.This implies that the following relations holds = 1, e j e j = −1, e j = −e j , j = 1, . . ., n, e i e j = −e j e i i = j.
(2.2) Here e I stands for the product ) and e 0 : = , e ∅ : = 1 is the multiplicative unit.Also, indicates that the sum is performed only over strictly increasing multi-indices, i.e., l-tuples We will use properties of the standard Clifford algebra C (0, n).In particular, any The Clifford conjugacion on C (1, n), denoted by 'bar', is defined as the extension on the conjugacion on C (0, n) such that = − and u = u for any u ∈ C (1, n).
More specifically, recall that on C (0, n) e I e I = e I e I = 1 for any multi-index I , One can also verify without difficulty that Next, for each ∈ {0, 1, . . ., n + 1} consider the projection map onto the -homogeneous part of u, i.e., and denote by the range of (2.9) We will use the subalgebra of even homogeneous element of C (1, n) which we denote by (2.12) Call u monogenic in if Du = 0 in .The fundamental solution of D is given by (2.13) We will use, only for functions of compact support, the Teodorescu transform that we recall.For a function ϕ(x) defined in R n with values in C (1, n) of compact support the Teodorescu transform of ϕ(x) is defined as The main property of T that we will use [11] Theorem 8.2 is that Note that the property above is true for C (0, n)-valued functions but it is easy to see that it extends to C (1, n)-valued functions.Now the Dirac-wave operator or Hyperbolic-Dirac operator associated with R 1,n , is given by where T ∈ (0, ∞]. Call u wave-monogenic in (0, T ) × if D w u = 0 in (0, T ) × .We will always act D w on the left.There is a parallel theory where D w acts on the right.This follows from the observation that D w u = 0 if and only if uD w = 0 (D w acting on the right).
One of the most fundamental properties of Dirac operators introduced above is that they can be thought of as square-roots of familiar second-order differential operators.More precisely, D and D w satisfy Another operator of interest will be D w : = ∂ t + D, which also satisfy D w 2 = .
Moreover as a consecuence of D = − D it is easy to see that If D w (u) = 0 we call u anti wave-monogenic.
We will be interesting in study the space of solutions of the Dirac-wave operator in R + × R n , namely (2.20) A few remarks that can be verify without difficulty are in order here.

Remark 2.1
The following properties hold (i) We requiere wave-monogenic functions to be at least C 2 in order to satisfy the wave operator.(ii) The space ) is a scalar-valued solution of the wave operator then D w u(t, x) is a wave-monogenic function.(iv) For any φ : x) then D w u(t, x) = 0 if and only if u 1 and u 2 satisfy the system (2.21) As a key example we present, when n=3, that the whole Maxwell system could be seen as a Dirac-wave function: Let E(t, x), B(t, x) be the electric field and magnetic field, in a vacuum, the Maxwell system is given by the following four equations: Here ρ is the charge density and J the electric current.Units are chosen so that the speed of light is 1.See for example [12] page 165.We could see the Maxwell system as a nonhomogeneous solution of D w in the following way: If (2.23) and (2.24) A simply computation shows that Maxwell system (2.22) is equivalent to D w u = Ĵ .In particular, if there is no charge or current then (2.22) is equivalent to D w u = 0.

Boundary Value Problem for the Dirac-Wave Operator
Here we discuss the initial-value problem for the Dirac-wave operator.We will work with data functions of class at least C 2 and in fact we could assume C ∞ .We will be concerned with local properties of solutions of the initial value problem for the Dirac-wave operator.Given φ : R n → C (1, n) a sufficiently smooth function.Consider the initial value problem for the Dirac-wave equation: Recall that we want our solution to be C 2 -continuous up to the boundary.From the remark (2.1)(v) we see that is it enough to take φ(x) with values in C + (1, n) and u will take values also in The case when n = 1, the d'Alembert formula is simply but nevertheless gives us intuition about what to expect in the case when n > 1.Consider φ(x) = f (x) + e 1 g(x) and u(t, x) = u 0 (t, x) + e 1 u 01 (t, x) to be C 1 and C 2 -function in C + (1, 1) respectively.Then it is straightforward to see that the following are equivalent.
• u satisfy the wave Cauchy-Riemann system and u(0, x) = f (x) + e 1 g(x).• Each component of u satisfy the BVP for u = 0 • u(t, x) is given by any of the formulas All the standard properties (slightly better) really follows.For example, the domain of dependence of u(t, x) are the endpoints of the interval [x − t, x + t] and the range of influence of an initial value φ(x 0 ) is the boundary of the wedge-shaped region {(t, x) : For most part of this section we shall entirely restrict our attention to the physically most relevant case, namely n = 3.All other cases are done in a analogous way using with n is even (1.8) and when n > 3 is odd (1.7).Nevertheless we emphasize the necessity of making more and more smoothness assumptions on the initial data to ensure the existence of a C 2 -solution.
First note that the Riemann transform act naturally in Clifford valued functions componentwise.i.e. if φ : R 3 → C (1, 3), then (3.5) In the next Theorem we find the solution to (3.1) and several equivalencies.(ii) u(0, x) = φ(x) and u satisfy in (0, ∞) × R 3 the 8x8 first order differential system (iii) For each multi-index I , u I satisfy the BVP u I = 0 ] This is just to write down u as in (2.10) and the equation We know that u = 0 taking components give us that for each multi-index I , u I = 0 and u I (0, x) = φ I (x).Note that D w u = 0 is equivalent to ∂ t u(t, x) = Du(t, x) taking the limit as t → 0 give the condition on the velocity.
Several consequences of this Theorem follows.We collect the first in a form of a remark.

Remark 3.2
The following properties follows for Dirac-wave solutions of (3.1) (i) Domain of dependence: u(t, x) depends only on the values of φ(y) and Dφ(y) for y ∈ ∂ B(x, t).Specifically, for each multi-index I , u I (t, x) depends only on the values of φ I (y) and ( Dφ(y)) I for y ∈ ∂ B(x, t).(ii) Range of influence: φ(y) influence u(t, x) only on the surface x − y = t.
Specifically, for each multi-index I , φ I (y) influence the components u I (t, x), u I 1 (t, x), u I 2 (t, x), u I 3 (t, x), only on the surface x − y = t, where I j is such that e j e I j = ±e I for j = 1, 2, 3.By way of example, φ 01 (y) influence the components u 01 (t, x), u 0 (t, x), u 12 (t, x), u 13 (t, x) on the surface x − y = t.
In the next corollary we will see conditions on which the solution to the BVP is of a particular form.Also we note that if the initial data has only certain components the solution u will have more components in general.
Corollary 3.3 Suppose the initial data φ for the BVP (3.1) is given as

.12)
This means all the other components are zero.Assume that φ is zero at infinity, consider the vector field F = (F i ), if u is the solution to (3.1) then • If curl F = 0 there exist a scalar valued wave-solution such that u(t, x) = (D w (t, x)) .• if curl F is not zero u has more that 4 components not zero.
• if Div F = 0 and φ 0 ≡ 0 then u has values in 2 and they satisfy a Maxwell system.

Proof
The proof is straightforward.
Here we consider n arbitrary again.We will give an application of the results above.For any k integer let P k (R n , C (1, n)) denote the space of homogeneous polynomials of degree k with values in C (1, n).We will denote by W k (R × R n ) the space of homogeneous wave-monogenic polynomials of degree k with values in C (1, n).We will find a nice basis for W k (R × R n ).
For each multiindex α Then it is well know that a basis of P k (R n , R) is given by this set of p α .Moving on, let P α (t, x) be the solution of the BVP (3.1) with initial data p α (x).We find an expression of P α as follows. Let then clearly W i are wave-monogenic polynomials of degree 1.We will follow the construction of the Fueter polynomials as in [11].For α as above, let . .W j k where the sequence of indices j 1 , j 2 , . . ., j k has the first k 1 indices equal 1, the next k 2 indices equal 2 and so on.Next we define where Perm(S m ) denote the permutation group of m elements.This is a symmetrization process.We see from [11] Theorem 6.2 (who proof applies to this case with minors changes) that where a α ∈ C (1, n).

Energy Estimates for the Dirac-Wave Operator
In this section we will see energy estimates for the Dirac-wave operator.Here n is arbitrary again.Note that we only use that u(x, t) is C 1 in this section.In the next result we show conservation of the energy.The proof use the Clifford algebra structure in a subtle way.
Theorem 4.1 Suppose that the initial data φ is a function of compact support and let u(x, t) be the solution of (3.1).Then ∀t ≥ 0 Proof We decompose u as in (2.4).Then we have that here we recall that a 0 is the real part of a. Therefore thanks to (2.21) and the fact that for any function f , D f = − f D where D acts on the right in the right side of the equation.
We will use the following version of the diverge Theorem Here η is the unit outward normal.In fact, we will only need the real component of the equation.Then we have thanks to (4.3) and (4.4) that Here we had use that u is a function of compact support.for some function g.
Proof if u and û are solutions of (4.6) let w = u − û then we proceed to compute the energy of w as above.We have no boundary terms.Therefore E(t) ≡ 0 and w ≡ 0.  Hence E (t) ≤ 0 for 0 ≤ t ≤ t 0 and E(0) = 0 implies that E(t) ≡ 0. Therefore u(t, x) ≡ 0 within C.

Dirac Wave Conjugates Functions
In this section we define the Dirac-wave conjugate to a wave function when n = 3. Again the other cases are analogous.
For starter given a scalar wave function u 0 we would like to find a Clifford-valued function u with real part u 0 which satisfaces the Dirac-wave operator.This is a wavemonogenic extension of u 0 .The Dirac-wave conjugate function will be then u − u 0 .Also from this, as it is done in the complex analysis case, we would like to define a Hilbert-wave transform as the limit as t go to zero of the Dirac-wave conjugate function.Actually more is true.We can give four wave functions u i (t, x) each one

5 )
with a, b, c, d ∈ C (0, n).If we now define a complex conjugation on C (1, n) by setting u c : = I u c I e I for each u = I u I e I , where z c denotes the usual complex conjugation of z ∈ C then this semi-inner structure induces the semi-norm and we have that on C (1, n), e I e I = e I e I = −1.
Here [m]  denotes the even number small o equal to m).As an example and for future references we have that any u ∈ C + (1, 3) can be written as u = u I e I = u 0 + u 01 e 1 + u 02 e 2 + u 03 e 3 +u 23 e 23 + u 13 e 13 + u 12 e 12 + u 0123 e 123 .(2.10) Moving on, recall that if we let be an open set in R n .Then the classical (homogeneous) Dirac operator associated with R n is given by D: = n j=1 e j ∂ j .(2.11)This acts on some u ∈ C 1 ( , C (0, n)) (or in C (1, n)) according to Du: = n j=1 e j ∂ j u.