## Anais da Academia Brasileira de Ciências

##
*Print version* ISSN 0001-3765*On-line version* ISSN 1678-2690

### An. Acad. Bras. Ciênc. vol.73 no.4 Rio de Janeiro Dec. 2001

#### http://dx.doi.org/10.1590/S0001-37652001000400002

**On Huygens' principle for Dirac operators associated to electromagnetic fields**

**FABIO A.C.C. CHALUB***

IMPA, Estrada Dona Castorina, 110 - 22460-320 Rio de Janeiro, RJ, Brazil

* Manuscript received on July 20, 2001; accepted for publication on August 22, 2001;*

* presented by *MANFREDO DO CARMO

**ABSTRACT**

We study the behavior of massless Dirac particles, i.e., solutions of the Dirac equation with *m* = 0 in the presence of an electromagnetic field. Our main result (Theorem 1) is that for purely real or imaginary fields any Huygens type (in Hadamard's sense) Dirac operators is equivalent to the free Dirac operator, equivalence given by changes of variables and multiplication (right and left) by nonzero functions.

**Key words:** Huygens' principle, Dirac operators, electromagnetic fields, Hadamard's problem.

**1. INTRODUCTION**

By 1690 two different theories were disputing the nature of light: for Newton (Newton, 1979), light was composed of particles propagating along straight lines; whereas Huygens believed that light was a mechanical wave propagating in the æther, the hypothetical medium that fulfilled all the space (Shapiro 1980).

Huygens' theory was particularly useful to explain diffraction ("the strange refraction of the Iceland crystal"), but had some remarkable difficulties to explain the apparent propagation of the light along straight lines. To overcome this difficulty, he wrote his masterpiece "Treatise on Light" (Huygens 1920) where he formulated the principle that bears his name.

In more modern terms (in a formulation essentially due to Fresnel (Fresnel 1816) -- for a historical overview, see (Born and Wolf 1980)) we say that each point in the wave front acts as a secondary source of wave and all these waves interfere to produce the new wave front. See Figure 1.

Despite the large use of the term "Huygens' principle" in the physical community, its precise mathematical formulation is seldom made clear. In 1923, the French mathematician J. Hadamard, in a series of lectures at Yale University, formulated in mathematical terms three different meanings of "Huygens' principle" he found in the literature of his time (Hadamard 1953). This work is concerned with the second definition, called "Hadamard's minor premise."

We define the forward light cone with vertex in *y* as

*y*.

We say that (and, by extension, its fundamental solution ) possesses Huygens' property (or obeys Huygens' principle) if, for every *y*, is supported in the forward light cone with vertex in *y*.

This definition is called "Hadamard's minor premise." See Figure 1.

As Huygens' principle was originally formulated for light propagation, a physical phenomenon modeled by the wave equation, it is useful to discuss Huygens' principle in Hadamard's minor premise sense in this case first.

The wave equation is a partial differential equation in *n* + 1 variables, (*x*^{0}, *x*^{1},^{ ... }, *x*^{n}), where *x*^{0} is taken as * time*, and the other variables are the * space* variables. Explicitly, it is an equation for a function *u* such that

*n*space variables. With these definitions we take the velocity of propagation

*c*equals 1. If we do not want to fix the value of

*c*, we change in the above and the following expression the value of

*x*

^{0}by

*cx*

^{0}.

It is a well known fact that the wave operator possesses Huygens' principle for odd *n* greater than or equal to 3 and does not posses it otherwise. See (Folland 1997).

* Hadamard's problem* consists in classifying all second order hyperbolic operators which obey Huygens' principle, up to trivial relations. The latter consist in successive combinations of changes of variables, and left or right multiplication of the operator by a never-singular smooth function. This problem has not been fully solved, but progress has been made recently. An extensive and up-to-date review can be found in Berest's work (Berest 1998). Previous work on the subject of Huygens' principle for Dirac operators are (Chalub and Zubelli 2000, Chalub and Zubelli 2001a).

In this work we will study the validity of Huygens' principle for Dirac operators in the presence of an electromagnetic field. In Section 2, we briefly introduce Clifford Algebras, that are instrumental to the definition of Dirac Operators. In Section 3, we present the Hadamard expansion and re-define Huygens' principle in terms of its truncation. In Section 4, we develop some comments about the physical nature of Huygens' principle, in particular in the presence of electromagnetic field and prove the main theorem of this article, namely, that every Dirac operator in the presence of a real or purely imaginary field with Huygens' property is trivially equivalent to the free Dirac operator (Theorem 1).

**2. CLIFFORD ALGEBRAS AND DIRAC OPERATORS**

Dirac operators appear in Mathematical Physics in the framework of the construction of a relativistic electron theory. In a celebrated paper, Dirac (Dirac 1928) defined a Lorentz-covariant Hamiltonian (and consequently first-order in space variables, as Schrödinger equation is first order in time) whose square furnishes the Klein-Gordon equation.

Let *g*^{mn} = *diag*[1, - 1,^{ ... }, - 1] denote the Minkowski tensor. Associated to *g*^{mn }we can construct a Clifford Algebra. It is an associative algebra (with identity *I*) over the reals generated by all linear combinations of the form , where the matrices {}, obey the relation + = 2*g*^{mn} *I*. The matrices {}, which are called Dirac matrices, are linearly independent (Gilbert and Murray 1991, Marchuk 1998).

Dirac operators are defined by:

= + *v* ,

*v*is of electromagnetic (or

*vector*) form,

*v*=

*a*

_{m}.

**3. HUYGENS' PRINCIPLE AND HADAMARD EXPANSIONS **

In this section we shall look for fundamental solutions of Dirac operators. A * fundamental solution* for is a solution of where denotes Dirac-delta distribution supported at an arbitrary point *y* in space-time.

Let (*x*, *y*) = denote the Minkowski distance between two points *x* and *y*.

The first important concept to be recalled is the so-called Riesz kernel. It was introduced by M. Riesz in order to unify the treatment of elliptic and hyperbolic problems. They are given by distributions in '(^{n + 1}) defined first for **Re**(a) sufficiently large by

*N*(a) is given by Then, is extended analytically for all values of a in the complex plane by means of = , which is satisfied for all a with sufficiently large real part. It is not hard to check that = . See (Folland 1997).

We shall say that a fundamental solution obeys * Huygens' principle* if *supp* ^{+}(*y*) *y* ^{n + 1}.

Adapting Hadamard's seminal idea from the wave operator case to the one at hand, we look for series expansion for the fundamental solution

where *s*_{m} is a matrix coefficient, , which we call Dirac kernels, and are Riesz kernels. It can be shown that satisfies Huygens' property for a = 0, - 2, - 4,^{ ... }; while does it for a = - 2, - 4,^{ ... }.

If is a fundamental solution, using the properties of Riesz and Dirac kernels, we have the ** Hadamard's recursion** for the coefficients *s*_{m}

In all these cases *s*_{m} should be smooth (in particular when *y*^{m}*x*^{m}) and we impose normalization *s*_{0}(*y*, *y*) = 1.

We solve equation (1)

An important property of this series is that it is unique, provided we require*s*

_{m}to be smooth. The uniqueness for (2) is trivial. In the case of (3) it should be proved. Let

*s*

^{1}and

*s*

^{2}be two different solutions of (3). So, the difference

*s*

^{1}-

*s*

^{2}=

*S*

_{2m}is smooth and satisfies

We find that

with satisfying

The only smooth solution in the vicinity of the vertex *y* is given by = 0, then *S*_{2m} = 0.

Bearing in mind this fact it is easy to prove the equivalence between Huygens' principle and the truncation of the Hadamard series:

REMARK 1. A Dirac operator in the form possesses Huygens' property in odd spatial dimension *n* its Hadamard series truncates at *n*, i.e, *s*_{m} = 0 for *m**n*.

PROOF. Proving this remark is quite simple. In the direction, it follows from the properties of Dirac and Riesz kernels. In the other direction, it is a consequence of the uniqueness of the Hadamard expansion.

**4. HUYGENS' PRINCIPLE IN AN ELECTROMAGNETIC FIELD**

In this section we shall consider potentials of the form .

Potentials of the form are called * electromagnetic potentials* because they correspond to including an electromagnetic field in the physical situation under consideration, except for a multiplication by the complex unit *i*. If the physical potential is given by the "four-vector" *A*^{m} the new momentum is given by , so the "physical case", in the expression of *v*, is given by purely imaginary components.

If we write the potential as a 1-form = , the electromagnetic field is given by its exterior derivative

The 2-form is often described by its tensor components,Electromagnetic fields obey Maxwell's equation, which in vacuum are written

In these equations *j* is the density of electrical charge and current, and the notation means that we sum over all the permutation of the indexes in brackets, taking in consideration the parity of the permutation (For example: *B*_{[mn]} = *B*_{mn} – *B*_{nm}, *B*_{[m}*C*_{n]} = *B*_{m}*C*_{n} – *B*_{n}*C*_{m} and *B*_{[m}*C*_{nr]} = *B*_{m}*C*_{nr} + *B*_{n}*C*_{rm} + *B*_{r}*C*_{mn} – *B*_{m}*C*_{rn} – *B*_{n}*C*_{mr} – *B*_{r}*C*_{nm}.). The second equation is automatically obeyed, as long as *F*_{mn} is obtained from a potential. Inside matter the second Maxwell equation remains valid (which implies, from Poincaré's lemma that, locally, the electromagnetic field can be obtained from a potential).

In classical physics the field at each point is observable, while the potentials are considered just a technique for calculating such field, having no separate reality (This changes dramatically in Quantum Mechanics, as can be seen by the Aharonov-Bohm effect. See, for example (Sakurai 1967)). This means that adding an exact form *d* to the potential should not change the field (which is obvious from the fact that *d*^{2}=0). In terms of its components, this means to change from *a*_{m} to . As the Huygens' property is a physically sound problem, the validity or not of the Huygens' principle cannot depend on the particular choice of the electromagnetic gauge (to understand the idea that Huygens' principle is physically observable we should understand the meaning of , the solution of Dirac equation. Following Section 1.7.3 of (Thaller 1992), we say that for every Borel-measurable set ^{3} the probability of finding the particle in is given by , *E*(), where *E*() is a projection operator in and the inner product should be interpreted as Hilbert space inner product where the Dirac operator is acting and the solution lives. In simpler words, we have a non-zero probability of finding the particle only on the support of the wave-function.).

The change of the gauge is a particular kind of trivial transformation, as defined in the introduction. We see immediately

so the change of the gauge does not change the validity of the Huygens' principle.

Thus, the invariance of Huygens' principle by change of gauge can be concluded by two distinct reasonings, coming from two different points of view, historically unrelated. We point out that, as far as we know, the fact that the change of electromagnetic gauge is a particular case of trivial transformations as defined by Hadamard was never mentioned in the literature.

After these remarks we are ready to state:

THEOREM 1. * If a Dirac operator obeys Huygens' property for real or purely imaginary potentials a^{m} in 3+1 dimensions, then it is equivalent to the free operator .*

PROOF. Let us define

where . The first remark is that

which shows that and are gauge equivalent. Explicitly, the change of gauge is given bywhere

From now on, the calculations will be made in the new gauge.

The functions *f*_{m} can be re-written as

(4) |

and, so, by the anti-symmetry of the tensor *F*_{mn} we have

(5) |

Since Hadamard's recursion stops at the second term, Hadamard's recursion becomes

(6) |

The first equation can be easily solved to get *s*_{0} = 1.

We shall study compatibility conditions between the last two equations. We start applying Dirac operator to the first:

With the second condition this can be written

Applying ¶_{n} to equation (5) we find that (*x*^{m} - *y*^{m})¶_{n }*f*_{m} = - *f*_{n} which allows one to write the equation

Finally we write

(7) |

The RHS of last equation is

For *x* = *y*, it is immediate to see that *f* (*x*, *y*) = 0. So, the only non-zero term in *x* = *y* in the RHS is *f*. Actually,

Letting *y**x*

As the LHS of equation (7) is zero at (*x*, *x*), we conclude that

(8) |

is a necessary condition for the operator to be a Huygens one. This can be interpreted as the absence of sources.

To get a second necessary condition, we apply

to the LHS of equation (7) and evaluate it at (*x*,

*x*). We get

Taking into account equations (4) and (8) the above system can be re-written as

where (*F*^{2})_{mn} indicates the (m, n) component of the tensor given by the square of the matrix *F*, .

The RHS of equation (7) gives

and *u*_{2}(*x*, *x*) can be easily obtained from equation (6):

We gather all this information in a single equation, written in the Clifford base as

where, as before means that we sum over all the permutation of m, n and r taking in consideration the parity of the permutation.

From the independence of the generators of the Clifford Algebra, we conclude that

Let us write (*F*^{2}) componentwise:

where *E*_{x}, *E*_{y}, *E*_{z}, *B*_{x}, *B*_{y} and *B*_{z} are the *x*, *y* and *z* components of the electric and magnetic fields, respectively.

From the condition *F*^{2} = 0 and the hypothesis that all the potentials are either real or purely imaginary (so are the fields), we conclude that *E*_{x} = *E*_{y} = *E*_{z} = *B*_{x} = *B*_{y} = *B*_{z} = 0, so *F*_{mn} = 0 and is gauge equivalent to the zero field.

REMARK 2. If we relax the hypothesis that the potential is real or purely imaginary, then we can find a large class of Huygens potentials. One family of examples is given by

where e(v) is any function such that e(v) = e*(– v) and * denotes complex conjugation.

**5. CONCLUSIONS**

Huygens' principle has been subject of extensive investigation in recent times. For reviews in two different directions, see (Berest 1998) and (Belger et al. 1997). One of its most striking characteristics is its relations to integrability (see (Berest and Veselov 1994)). In the wave operator case, potentials of Huygens' type are, after a suitable change of variables, solutions of the Korteweg-de Vries equation

*u*_{t} = *u*_{xxx} - 6*uu*_{x} .

For Dirac operators, the relation between Huygens' principle and integrability involves the AKNS (after Ablowitz, Kaup, Newell and Segur) hierarchy, which in turn is an example of integrable matrix systems. This is proved in (Chalub and Zubelli 2001b).

A natural follow up of the present work is to relate Huygens' principle for Dirac operators associated to the electromagnetic field and integrable systems.

**ACKNOWLEDGMENTS**

I thank Prof. Jorge P. Zubelli for introducing me on this subject and for many suggestions. This work was supported by the Brazilian Agencies CNPq and CAPES.

**RESUMO**

Neste artigo estudamos o comportamento das soluções da equação de Dirac sem massa ("massless Dirac particles") na presença de um campo eletromagnético. Nosso resultado (Teorema 1) indica que para campos reais ou puramente imaginários todo operador de Dirac que obedece ao princípio de Huygens é equivalente ao operador trivial, equivalência dada por mudanças de variáveis e multiplicação, à direita ou esquerda, por funções não singulares.

**Palavras-chave:** princípio de Huygens, operadores de Dirac, campos eletromagnéticos, problema de Hadamard.

**REFERENCES**

BELGER M, SCHIMMING R AND WÜNSCH V. 1997. A survey on Huygens' principle. Z Anal Anwendungen, 16(1): 9-36. Dedicated to the memory of Paul Günther. [ Links ]

BEREST Y AND VESELOV AP. 1994. The Huygens principle and integrability. Uspekhi Mat Nauk, 49(6(300)): 7-78. [ Links ]

BEREST Y. 1998. Hierarchies of Huygens' operators and Hadamard's conjecture. Acta Appl Math, 53(2): 125-185. [ Links ]

BORN M AND WOLF E. 1980. Eletromagnetic Theory of Propagation, Interference and Diffraction of Light. Pergamon Press, sixth edition. [ Links ]

CHALUB FACC AND ZUBELLI JP. 2000. Integrable systems, Huygens' principle, and Dirac operators. In Proceedings of the Workshop on Nonlinearity, Integrability and All That: Twenty Years after NEEDS '79 (Gallipoli, 1999), pages 89-96, River Edge, NJ. World Sci. Publishing. [ Links ]

CHALUB FACC AND ZUBELLI JP. 2001a. On Huygens' Principle for Dirac Operators and Nonlinear Evolution Equations. J Nonlin Math Phys, 8, Supplement: 62-68. [ Links ]

CHALUB FACC AND ZUBELLI JP. 2001b. Huygens' Principle, Dirac Operators and Rational Solutions of the AKNS Hierarchy. Submitted.

DIRAC PAM. 1928. The Quantum Theory of the Electron. Proc Roy Soc of London A, 117: 610-624. [ Links ]

FRESNEL A. 1816. Ann Chim et Phys, (2), 1, 239. in Oeuvres, Vol 1, 89, 129.

FOLLAND GB. 1997. Fundamental solutions for the wave operator. Exposition Math, 15(1): 25-52. [ Links ]

GILBERT JE AND MURRAY MAM. 1991. Clifford algebras and Dirac operators in harmonic analysis. Cambridge University Press, Cambridge. [ Links ]

HADAMARD J. 1953. Lectures on Cauchy's problem in linear partial differential equations. Dover Publications, New York. [ Links ]

HUYGENS C. 1920. Traité de La Lumiére, où sont expliquées les causes de ce qui luy-arrive dans la reflexion, et dans la refraction et particulièrement dans l'etrange refraction du Cristal D'Islande. Societé Hollandaise des Sciences, 1920. Reprint of the 1690 original, Librarie Marchand. [ Links ]

MARCHUK NG. 1998. Dirac g-equation, classical gauge fields and Clifford algebra. Adv Appl Clifford Algebras, 8(1): 181-225. [ Links ]

NEWTON I. 1979. Opticks. A treatise of the Reflexion, Inflexion and Colours of Light. Dover Publications, New York. [ Links ]

SAKURAI JJ. 1967. Advanced Quantum Mechanics. Addison Wesley. [ Links ]

SHAPIRO A. 1980. Huygens' kinematic theory of light. In Studies on Christiaan Huygens. Swets and Zeitlinger B.V., Lisse.

THALLER B. 1992. The Dirac equation. Springer-Verlag, Berlin. [ Links ]

* Present address :

Dept. of Mathematics, University of Vienna

Boltzmangasse 9 - A-1090, Wien, Austria

E-mail: fabio@impa.br