Arxiv:2009.10356V2 [Gr-Qc] 24 Sep 2020 Aiirt Vr Tdn.I Resae Hs Equations By: These Units Space, SI Free in in Given Are Student

Home , Pseudovector

arXiv:2009.10356v2 [gr-qc] 24 Sep 2020 aiirt vr tdn.I resae hs equations by: these units space, SI free in In given are student. form every to compact to more calculus the familiar vector into developed equations Hertz Maxwell’s rewrite H. magnetism and Heaviside and O. electricity his on in J.C. algebra tise Although quaternions used vector. originally three-dimensional Maxwell of use of extensive concept an the common: in something keep di- still (EM) from remote objects intangible experience. of accessible favor famil- rectly in the mass, up (acceleration, force...) give mechanics Newtonian to is of concepts requires what introduced iar field understanding field, Actually, of electromagnetic concept Faraday. the the M. in by lies firstly answer the gests steeeti u density, flux electric the is etrfil vnsiga nnt)i nqeadwell- specified are and divergence unique its and is curl infinity) its provided at defined (vanishing any field that vector states which theorem, decomposition Helmholtz density. flux magnetic the is 2 n 3 are (3) equations and nonhomogeneous (2) the in source-terms the Here, h lcrcfe urn density. current free electric the osdrddffiut na sa entitled essay an understand? to In hard so theory Maxwell’s difficult? considered INTRODUCTION I. rlRrld enmuo 27-0,Rcf,P,Brazil PE, Recife, 52171-900, Pernambuco, de Rural eral b) a) c the electrodynamics: of understanding student Improving lcrncmi:[email protected] mail: Electronic emnn drs:Dpraet eF´ısica, Universidade de Departamento address: Permanent e ehnclmdl n lsia electromagnetism classical and models mechanical Yet h rtfu qain r eial eas of because desirable are equations four first The h sudrrdaeeetoants ocommonly so electromagnetism undergraduate is Why oiae h hsc ece obigtesbeto differential situ po of nontrivial the subject in illustrate the electromagnetism bring to to to with used teacher due are dealing physics vacuum the effect in motivates classical Casimir forms kn the ca the differential previous of that in no of insights polarization vacuum that relevant effective quantum in the The self-contained the of is some through. article forth The it bring d E&M. seldom to to is article approach electromagnetism this in forms of differential aim of role illuminating The 701 UMR Théoriques, CNRS Chimie France. Vandœuvre-les-Nancy, et Physique de Laboratoire .Fumeron, S. ρ h lcrcfe hredniy and density, charge free electric the , ulH curl ulE curl 1, a) 2,3 .Berche, B. − + div div H ∂ ∂ ∂t ∂t D B D D stemgei edand field magnetic the is B B = = (1) 0 = = (4) = 0 = E j ρ µ ε 1 steeeti field, electric the is 0 0 17) ..Gibbs, J.W. (1873), n .Moraes F. and E H , 1 .Dsnsug- Dyson F. h is Why 1, Trea- b) Fed- (3) (2) (6) (5) D B j , 1 lascle eeti ed n only and field’ ‘electric called always oc’(lcrclegneig:telte he suggest three example for latter the ‘magnetising even that engineering): or in- (electrical field strength’ ‘magnetic force’ field also but ‘magnetic here), tensity’, choice (our field’ netic ato h edhsals iprt terminology: disparate less a field’. has ‘magnetic field the the some- as of textbooks it part EM to but refer here) times choice (our ‘magnetic law), density’ Faraday’s flux to (ought field’ induction netic atasto uhcoc sta tepaie h parallel the emphasizes it that impor- is between (an choice density’ such flux of ’electric asset of tant term less-used the for n awl’ qain nmte sefreapethe example for (see matter in review equations Maxwell’s ing hand). other the nuto’ hr gi,telte w neln h idea the underline two latter ’electric of the or again, field’ There displacement induction’. ’electric density’, flux tric iua,tesotigo oretrsgvsrs onew to rise namely gives fields, source-terms effective of macroscopic par- smoothing in the and quantities ticular, macroscopic retrieve proce- to averaging allows Lorentz-Rosenfeld dure a Performing fields space. microscopic varying rapidly produce e medium a in n ad eedn ntecontext, the on the On depending terminology. hand, the one to attention paying by obtained close equations. to of order known system in required equations, the are two the relations, time last constitutive from partial the The as the only fields). from the not also of but supplied derivatives currents, are and charges latter free these that (note rmrl nmte.Teeoe nvcu constitutive vacuum in relevant Therefore, excitations) relations, matter. (or in quantities primarily auxiliary as appear oeaiehwtefilstasomudra under transform fields transformation the coordinate how general examine to multiplicative a to up equal constant. simply physical are they as fields dant { aoinmti,Pnr n Ward and Pendry matrix, Jacobian u ( s o ieeta forms differential for ase x l hsmycm rmawdsra a o build- for way widespread a from come may this All oeat ftemsocpin nonee nE is EM in encountered misconceptions the of foretaste A atclryisgtu xriei then is exercise insightful particularly A ( D ,y z y, x, t , H and ) en naxlayqatt n hrfr,w opted we therefore, and quantity auxiliary an being 7 :a ugse yLrnz lcrn n nuclei and electrons Lorentz, by suggested as ): snttefnaetlmgei uniy(see quantity magnetic fundamental the not is E ) v , B and b 4 ( ( .O h te hand, other the On ). om oteclassroom. the to forms and ,y z y, x, ,Uiested orie 54506, Lorraine, Université de 9, x t , H ,ta byMxelseutosi free in equations Maxwell’s obey that ), H nteoehn,between hand, one the on ) elan rmadffrnilforms differential a from learnt be n weg ffrsi eddt follow to needed is forms of owledge w , nfr rvttoa edadof and field gravitational uniform a (resp. tos ehp htti article this that hope We ations. e n aieso interpretation of easiness and wer susdi h lsro.I sthe is It classroom. the in iscussed ( ,y z y, x, D ) } and eoigby Denoting . E D B emt eredun- be to seem ) D and skona ‘mag- as known is H skona ’elec- as known is 8 5,6 scle ‘mag- called is { H ,y z y, x, hwdthat showed h electric The D hc thus which , and → } J B E the on is Maxwell’s equations remain invariant provided the old exterior-calculus-based electrodynamics course.11 Under- and the new fields (primed quantities) are related by: graduate students are generally familiar with exact (or total) differentials, as those abound in thermodynamics ′ −T ′ −T E = J E H = J H (7) courses. Given a Cartesian coordinate system on the Eu- J J clidean space R3, the total differential of any scalar func- D′ = D B′ = B (8) detJ detJ tion F writes as (Einstein’s summation convention on repeated indices is used) D and E (resp. B and H) do not transform in the same way and hence, they are not redundant although mathe- ∂F ∂F ∂F ∂F dF = dx + dy + dz = dxa. (9) matically, they are described by the same kind of object: ∂x ∂y ∂z ∂xa vectors. Another intriguing point is that E and H obey a sim- When changing the function, only the partial derivatives ilar transformation law (as D and B do), suggesting the will change: this means that anytime an exact differential pair share a common mathematical nature. Sometimes, is computed, one is working in a vector space whose basis an additional clarification is introduced based on mirror elements are the dxa,a =1..3 = dx, dy, dz . { } { } symmetries, dividing the fields into polar (or true) vec- More generally, any object (not necessarily a total dif- tors such as E and D, and axial vectors (or pseudovec- ferential) that is written as a linear combination of dxj is tors) such as B and H.3 called a 1-form and it belongs to the 1-form vector space However, it is obviously not enough to enlighten the (or cotangent space) denoted by Λ1 R3 . Once given the form of all (7)-(8) and as remarked by Kitano,9 “in spite 1-form vector space, more general objects can be built of the simple appearance, the constitutive relations, even in a natural way: this is the idea underlying the more for the case of vacuum, are the non-trivial part of the EM general concept of tensors (technically, a (m,p)-tensor is theory.”. simply a multilinear map acting on a collection of m 1- The fact that D and E (resp. B and H) are not redun- forms and p vectors to produce a real number - for more 30 dant in free space is largely unknown and it raises at least details, see for example ). Taking the antisymmetrized tensor product (denoted for short by ) for each pair of two questions: 1) What is the true nature of their con- ∧ nection? 2) Are there pedagogical examples that could 1-form basis elements gives help illustrate their different physical contents? In this dxa dxb dxb dxa = dxa dxb. (10) paper, we use exterior calculus to answer these two ques- ⊗ − ⊗ ∧ tions and in doing so, all the points raised in the previous Here, the regular tensor product is defined as the or- paragraphs. dered product of pairs of 1-forms⊗ (and of vectors) and it Exterior calculus originates from the pioneering works is associative. of Grassmann and Cartan and it is concerned with Now, it appears that one can generate only a finite the properties of mathematical objects called differential number of non-zero terms, which are a linearly inde- forms. These latter have raised considerable attention pendent and spanning subset of a new vector space: because of their many applications in physics, as testified Λ2 R3 , the space of 2-forms. For example, in Carte- by the number of articles9–17 and books18–27 dedicated to sian coordinates, the 2-form basis written in the right differential forms, but they are usually not introduced at cyclic order is the set: dy dz, dz dx, dx dy an undergraduate level and as far as we know, only more (no element dxa is repeated as∧ (10) would∧ return∧ 0). advanced textbooks make abundant use of this formalism 19,25,28,29 That process can be iterated for p-uples basis elements (see e.g. ). a1 ap dx .. dx ,a1 = .. = ap and generates forms of In the first section of this work, the basic ideas of ex- {degree∧p ∧(or p-forms)6 that6 } belong to a vector space terior calculus are introduced in an intuitive way. Only a Λp R3 of dimension Cp = 3 !/(p !(3 p) !). By con- very basic knowledge of linear algebra will be used to de- 3 − struction, a p-form is a completely antisymmetric (0,p) fine differential forms along with the different operations tensor. Concretely, for Cartesian coordinates in R3, the which they allow to perform. A particular emphasis will general expression for a form of degree be put on the assets of this formalism from the standpoint of teaching EM in introductory university courses. Then, 0 is f(x,y,z) we illustrate and discuss some key differences between D 1 is f1(x,y,z)dx + f2(x,y,z)dy + f3(x,y,z)dz and E by working out two examples. 2 is g (x,y,z)dy dz + g (x,y,z)dz dx 1 ∧ 2 ∧ +g (x,y,z)dx dy 3 ∧ II. A PRIMER ON EXTERIOR CALCULUS 3 is g(x,y,z)dx dy dz ∧ ∧ As seen above, the cotangent space is a vector space A. What are differential forms? associated to 1-forms. In a similar fashion one associates ordinary vectors to a ‘tangent space’ as follows. A vec- Introducing differential forms in a handy but yet ac- tor can be seen as an operator which, when acting on a curate way is probably the most challenging part of an function, gives its directional derivative along the vector

2 direction. That is, taking the directional derivative of a xk are points evenly spaced in the interval [a,b]. This function F along the vector v at the point x one gets suggests that “dx” is in practice a shorthand notation for an infinitely thin width, which it is not. As explained a 31 DvF (x)= v F (x)= v ∂aF (x). (11) by Spivak, “Classical differential geometers (and clas- · ∇ sical analysts) did not hesitate to talk about ‘infinitely We see that there is a one-to-one correspondence between small’ changes dxi of the coordinates xi, just as Leibnitz vectors and derivations such that the vector v may be had. No one wanted to admit that this was nonsense, represented by its components in the basis ∂a . It fol- { } because true results were obtained when these infinitely lows that the linear space spanned by this basis is a vector small quantities were divided into each other”. A rig- space called the tangent space at the point x R3. Like ∈ orous connection between infinitesimal elements and 1- in the case of the p-forms, the wedge product between forms belongs to the realm of non-standard analysis, a vectors can be used to generate higher order vectors, or branch of mathematics developed in the 60s by Abraham p-vectors. Robinson.32 The spaces of p-forms and of p-vectors can be defined on base spaces that are more general than R3. For any To sum up, the differential form dx is rightly under- manifold of dimension n (that is a smooth hypersur- stood as a basis element of the cotangent space, but not face thatM locally looks like Rn), once given a coordinate as a tiny change: this is why one favors the introduc- system xa,a =1..n , one gets the vector coordinate ba- tion of differential forms from the algebraic standpoint { a } a instead of the more widespread analysis standpoint. sis ea = ∂/∂x = ∂a and a 1-form basis dx . These { } { } two sets are said dual, in the sense that when they act A popular way to gain insight on 1, 2 and 3-forms in upon each other, it returns a real number as prescribed R3 is to use graphical representations known as Schouten by pictograms.33 The main idea is to plot a p-form field ω as a foliation built from its kernel Kerω. Recalling that dif- a a a dx (∂b)= ∂b (dx )= δb . (12) ferential forms act on vectors returning real numbers (see Eq. (12)), the kernel of a form α is the set of all vectors In other words, a 1-form is a linear map from a vector v such that α(v) = 0. A foliation is a family of hypersur- R space onto and more generally, a p-form is a multilinear faces of dimension 3 p, filling R3 with no overlap. One R map acting on a collection of p-vectors onto . Note that can define it from the− kernel if and only if Kerω is an inte- in mathematics, there are other notions of duality, such grable plane field (this is known as Pfaff’s problem and a as Hodge duality which will be introduced in section II C. general criterion involving the exterior derivative d intro- Next to that algebraic standpoint on forms, another duced in the next paragraph is ω dω = 0). For 1-forms, way to grasp forms is from integral calculus. Indeed, dif- integrability means that there is∧ always a 2-surface that ferential forms “are the things which occur under integral is locally tangent to the plane field defined by Ker ω. To 22 signs”: a form of degree p is simply an object that one be more concrete, let us focus on the 1-form ω = dr in integrates p-times to get a scalar. Stated otherwise, “a cylindrical coordinates r,θ,z : its kernel is simply ob- differential form is simply this: an integrand”.26 A par- { } tained at each point from (12) dr(∂θ)=0, dr(∂z) = 0. ticular care must be taken when introducing forms that Hence, the plane field generated by vectors ∂θ, ∂z is way. the family of planes, tangent to 2-surfaces r ={Cst. } First, exterior calculus relies strictly speaking on Lebesgue integrals, not on Riemann integrals. Although it might look appealing, this visual interpre- Second, it is customary in calculus to drop the wedge tation may help for a first initiation to the language of sign when computing double (multiple) integrals such as forms, but this rarely goes beyond an introductory level: f(x, y)dxdy, where the domain of integration is de- as discussed in Ref. 26, there are many simple forms R whose kernels are not integrable and present “contact finedRR as R = [x1, x2] [y1,y2]. But then, the Fubini’s exchange-of-integration-order× formula may look contra- structures” (for example xdy + dz). Hence, we will not dictory with the antisymmetry property of the 2-form rely on that point of view in the reminder of this work. dx dy. Although subtle, the disagreement is nothing∧ but apparent: indeed, switching dx with dy also involves a reparametrization of the integration domain p : R R′ = [y ,y ] [x , x ]. As for any p-form u → 1 2 × 1 2 u = u if the map p reverses the orientation B. Wedge product, exterior derivative and R − p(R) RRof the integrationRR domain (conversely, there is a plus sign vector calculus identities if the map preserves orientation), then the antisymmetry feature of forms is retained in the customary calculus formulas. Relation (10) introduces an algebraic operation be- Third, introductory calculus courses define integrals as tween 1-forms: the wedge (or exterior) product, denoted by the symbol . More generally, given two forms u, v the limits of Riemann sums when the partition becomes ∧ b (of degree m), w (of degree p) and z (of degree q), it can infinitely fine, f(x)dx = lim f(xk)∆x, where the a ∆x→0 be shown easily that the wedge product obeys the three R Pk

3 following properties: necessarily closed (its exterior derivative vanishes). Con- versely, Poincaré’s lemma states that a closed p-form is u w = ( 1)mpw u (supercommutativity) (13) ∧ − ∧ always locally exact but generally not globally, depending u (w z) = (u w) z (associativity) (14) on the topology of the base space (this is the object ∧ ∧ ∧ ∧ M (αu + βv) w = αu w + βv w (linearity) (15) of the de Rham cohomology). ∧ ∧ ∧ The exterior derivative of terms such as u w is a form of degree m + p. As remarked previously, g x1, x2... dxa1 .. dxap is obtained using both ∧ a a ∧ ∧ supercommutativity implies that dx dx = 0: hence, the Leibniz formula and the nilpotence property in (18) for a base space of dimension n, no form∧ of degree p>n can exist. The set of all Λp ( ) ,p =0..n equipped d g x1, x2... dxa1 .. dxap = (dg) dxa1 .. dxap ∧ ∧ ∧ ∧ ∧ with the exterior product{ definesM an algebra} known as (21) the Grassmann algebra. In practice, the exterior derivative of a p-form simply Although the meaning of the wedge product may be consists in computing the sum of the partial derivatives obscure, it is still a familiar object to students in the of its components and then in discarding terms within usual Euclidean space. which a same dxa is repeated. When u is a 0-form (a scalar) and w is a p-form, the It is interesting to point out the coherence of the no- wedge is simply the scalar multiplication and it will be tations employed heretofore: for the 0-form u = x1, (21) omitted. leads to du = d(x1)= dx1. When computing the wedge product of two 1-forms The exterior derivative must not be confused with ei- a a u = uadx and v = vadx , it comes ther the Lie derivative (which computes variations with respect to a vector field and is connected to d via the u v =(uyvz uzvy) dy dz + (uzvx uxvz) dz dx ∧ − ∧ − ∧ Cartan formula) or with the covariant derivative (which + (uxvy uyvx) dx dy. (16) is dedicated to parallel transport and requires a connec- − ∧ tion on the manifold). The components of u v look like the components ob- ∧ Similarly to the wedge product, some insight on d can tained from the usual cross product between vectors be gained by expressing (19)-(20) with forms in R3. For (uxex + uyey + uzez) (vxex + vyey + vzez) (the exact × example in Cartesian coordinates, the exterior derivative forms-vectors connection and its inferences will be refined of a 0-form F (x,y,z) gives (9), but now with a clarified later, see also Ref. 12). Roughly speaking, when applied meaning for the dxa. Hence, the exterior derivative of a to 1-forms in R3, has the meaning of a cross-product. ∧ 0-form returns a 1-form with components similar to the What about the wedge product between a 1-form and a ordinary gradient operator. Likewise, the exterior deriva- 2-form? A straightforward calculation shows that for a 1- tive of a 1-form u gives after straightforward algebra form u and a 2-form w = wxdy dz+wydz dx+wzdx dy, ∧ ∧ ∧ the wedge product yields du =(∂xuy ∂yux) dx dy + (∂yuz ∂zuy) dy dz − ∧ − ∧ + (∂zux ∂xuz) dz dx. (22) u w = (uxwx + uywy + uzwz) dx dy dz. (17) − ∧ ∧ ∧ ∧ One recognizes a 2-form whose components correspond to This time, behaves like the ordinary dot product be- the ordinary curl operator of a vector field with appro- ∧ tween the components of u and w. priate components (discussed later). Finally, applying The wedge product concatenates two forms to produce the exterior derivative to a 2-form w leads to a 3-form a new form of higher degree. Another possibility to in- with components analogous to the ordinary divergence crease the degree of a form is to take the exterior deriva- operator tive of a form. The exterior derivative of a p-form is a linear mapping from Λp ( ) to Λp+1 ( ) defined for- dw = (∂xwx + ∂ywy + ∂zwz) dx dy dz. (23) M M ∧ ∧ mally in n dimensions as To sum up, grad, curl and div are just a unique operator in disguise, d, applied to forms of different degrees. ∂ a dω a dx ω a =1..n (18) But there is more: the exterior derivative not only unifies ≡ ∂x ∧ grad, curl and div, but it also generalizes them to arbitrary dimensions (this is responsible for Burke’s thought- It generalizes the notion of total differential of a scalar 34 function (9). Unlike the ordinary derivative, it is a di- provoking statement “Div, Grad, Curl are Dead”. ) mensionless operator. Besides, it obeys the two following As a sidenote, it is worth noting that the inverse op- properties eration to derivation, integration, also exists in exterior calculus, but there are some subtleties compared to or- d (u v) = (du) v + ( 1)p u dv (Leibniz) (19) dinary integration familiar to students. There is indeed ∧ ∧ − ∧ d (du) = 0 (Nilpotence)(20) a natural way to integrate a p-form over a p-dimensional submanifold (note that the degree of the form has to be with p the degree of u. the same as the dimension of the domain of integration), From the nilpotence property, if a p-form ω is exact which does not rely on lengths, angles or scalar prod- (i.e. there is a (p 1)-form u such that ω = du), then it is ucts and is hence purely topological. Such integral is not −

4 defined in the sense of Riemann but it is properly per- the displacement vector ∆s = (∆x, ∆y) has a length formed with respect to Lebesgue measure.35 In practice, squared given by Pythagoras theorem ∆s2 = ∆x2 +∆y2, this means that integration of a 1-form u on a path so that g11 = 1,g22 = 1,g12 = g21 = 0. However, on C the surface of a sphere (of radius R), an arc of length 2 2 2 2 2 2 u = uxdx + uydy + uzdz squared is given by ∆s = R ∆θ + R sin θ∆φ so that 2 2 2 ZC ZC g11 = R ,g22 = R sin θ,g12 = g21 = 0. From these examples, one understands that the components of g are can not be interpreted as the usual circulation of a vector related in particular to the curvature of the manifold, a field u along the curve C property which remains true in higher dimensions. The metric is real and symmetric, hence the inverse metric u ds cb b . exists and obeys gacg = δa . This provides a one-to- ZC one correspondence between the components of a 1-form as the scalar product requires the knowledge of the geom- v and those of a vector v: etry of (technically, the metric). But there is worse: M b there is no consistent way to define the flux of u, as that va = gabv . (26) would involve integration of a 1-form on a domain of dimension 2. As can be guessed, the degree of u must Locally, there is always an orthonormal coordinate ba- be increased beforehand, as expressed by the generalized sis such that gab reduces to a diagonal matrix having only Stokes theorem:13 +1’s or 1’s on its diagonal: the signature of the metric is the couple− (s,n s) where s is the number of -1’s. From the standpoint of− general relativity, the metric tensor is u = du. (24) Lorentzian (e.g. with signature=(1,3)) and it is obtained 18,30,37,38 Z∂D ZD as a solution of Einstein’s field equations. The inverse metric gab defines the inner product , between two p-forms according to: h i C. Hodge star operator p = 1 : dxa, dxb = gab (27) Another ingredient is necessary to reformulate EM p> 1 : dxa1 .. dxap , dxb1 .. dxbp with differential forms. To apprehend this, consider the ∧ ∧ ∧ ∧ ga1b1 ga1bp example of a vector field whose counterpart is a 1-form ···  . . .  field. As we just saw, the exterior derivation only allows = . .. . . (28) us to build a curl-like equation, leaving the divergence gapb1 gapbp  part undefined. This violates Helmholtz theorem (what-  ···  36 ever its time-dependent form may be ), meaning that The inner product between 1-forms is the equivalent of the field is not defined unambiguously. Another vari- the dot product for vectors. ation on the same theme: in virtue of the nilpotence Finally, having at our disposal the metric and the inner property (20), having only and d at our disposal leaves ∧ product, the Hodge dual operator (or star operator ⋆) is no room for equations involving a Laplacian. As we shall an invertible linear map between v Λp ( ) and ⋆ v see, these two problems are solved by the Hodge dual- Λn−p ( ) such that21 ∈ M ∈ ity, which consists in a mapping between p-forms and M (n p)-forms (remember that n is the dimension of the u (⋆ v)= u, v detg dx1 .. dxn (29) manifold).− Beforehand, one must introduce two objects: ab ∧ h i p| | ∧ ∧ the metric tensor and the inner product. (here u is of the same degree as v). So far, all mathematical objects that have been intro- An equivalent and more explicit definition is given by39 duced are purely topological, i.e. they are relying neither on angles nor on distances. These latter can be computed 1 µ1...µp µp+1 µn ⋆ v = v detgab ǫµ ...µ dx ... dx , only when the base space is equipped with a metric struc- (n p)! | | 1 n ∧ ∧ ture ( is then called a pseudo-Riemannian manifold). − p M (30) The metric consists in a symmetric (0, 2) tensor denoted where by g, which defines the scalar product between two vectors: µ1 µp v = vµ ...µ dx ... dx . (31) 1 p ∧ ∧

a b and the totally anti-symmetric Levi-Civita symbol is de- g (u, v)= gab dx (u) dx (v) . (25) ⊗ ﬁned as From (12), its appears that the components of g are ǫ = +1 if µ , ..µ is an even permutation of 1,..n simply the scalar products between the coordinate ba- µ1...µn 1 n = 1 if µ , ..µn is an odd permutation of 1,..n sis elements gab = g (∂a, ∂b). What is the physical con- − 1 tent of these components? In the ﬂat Eulidean plane, = 0 otherwise. (32)

5 Contrary to and d, the Hodge dual is a metric- This is the main result of the Hodge decomposition the- dependent operator∧ and it obeys the properties: orem, that generalizes the Helmholtz theorem.41 More- − over, the Hodge operator connects in a precise manner ⋆ (⋆ v) = ( 1)s+p(n p) v (33) − , d and their counterparts in vector analysis. The case u (⋆ v)= v (⋆u) (34) of∧ 1-forms has already been settled: they are translated ∧ ∧ u (⋆u)=0 u =0. (35) unambiguously into vectors as prescribed by (26). Now, ∧ ⇒ looking back at (17), we remarked that when dealing with To understand what (29) really does, let us illustrate a 1-form u and a 2-form v, u v behaves like the ordinary its action on the p-form basis for the flat Euclidean space dot product: but as the right∧ hand side is a 3-form and in spherical coordinates r, θ, φ : not a 0-form (a scalar), it is in fact ⋆(u v) (and not { } ∧ 1 simply u v) that actually returns a 0-form and corre- ⋆ (dθ dφ)= dr, (36) sponds to∧ the ordinary dot product. A similar statement ∧ r2 sin θ holds for divergence which corresponds to ⋆ d of a 2-form ⋆ (dr dθ) = sin θdφ (37) ∧ (and not simply d) as appears from (23). More generally, 1 ⋆ (dφ dr)= dθ, (38) taking the Hodge dual translates unambiguously 3-forms ∧ sin θ into scalars. The case of 2-forms may look trickier at 1 ⋆ (dr dθ dφ)= . (39) first glance, as they can be translated either as bivectors ∧ ∧ r2 sin θ (using the metric two times to raise each index) or as These relations are of course completed by means of vectors (taking the Hodge star and then using the metric (33). When the metric is diagonal, they are obtained on the 1-form obtained). However, it can be shown after for any curvilinear coordinates by following the simple some algebra that, for a diagonal metric in 3 dimensions, recipe: respecting the right cyclic order, the action of these two procedures lead to the same components (this the Hodge dual on a p-uple returns the missing (n is no longer true in 4 dimensions). Hence, as summarized p)-uple, each dxa being multiplied by its corresponding− in Ref. 11, for u, v Λ1 , u v identifies with ⋆(u v) ∈ M × ∧ Lamécoefficient (the square root of the metric √gaa (no and curl u with ⋆ du. index summation)). Although ⋆ v appears as some kind of orthogonal complement of the p-form v, we will see later that in EM it can be pictured more accurately as a III. OUTCOMES rotation between the electric and magnetic parts of the field. A. Unification of vector analysis Having the Hodge operator at our disposal, one is now in position to address the questions asked at the begin- As seen previously, depending on the degree of the ning of this section. Symmetrically to the exterior deriva- forms between which it is applied, and up to a Hodge tive which increases the degree of forms, one can now de- dual, the wedge product unifies scalar multiplication, fine a derivation operation that lowers the degree of the cross product and scalar product whereas exterior deriva- forms: this is the coderivative δ defined as tive unifies grad, curl and div operators. Bearing that in s+p(n−p) mind, let us investigate on a few examples how (19)-(20) δv = ( 1) ⋆ d ⋆ v. (40) − are translated for forms of different degrees in R3. Be- The coderivative maps Λp ( ) Λp−1 ( ), in a similar ginning with the nilpotence property applied to a 0-form fashion as a divergence operatorM → whichM returns a scalar F , dF stands for the gradient and as it is a 1-form, its from a higher-degree object, vector or pseudovector. It exterior derivative gives the curl, so that in the end one obeys the two following properties:40 retrieves the well-known formula δ(δu) = 0 (Nilpotence) (41) ⋆ d (dF )=0 curl (grad F )=0. (44) → du (⋆ v)= u ⋆(δv)+ d (u ⋆v ) (42) ∧ ∧ ∧ For a 1-form u, the same approach leads to That is the missing piece required to construct second derivatives that do not identically vanish. Indeed, the ⋆ d (du)=0 div (curl u)=0. (45) anticommutator between d and δ, known as the Laplace- → de Rham operator, generalizes the usual Laplace operator Considering the Leibniz formula, in the case of two 0- for vectors in R3 to p-forms on an arbitrary base space forms F and G, it comes that : 0 M d (F G) = (dF ) G + ( 1) F dG ∧ ∧ − ∧ ∆v = (δd + dδ) v. (43) grad (F G)= G grad F + F grad G. (46) → It also provides a way to close the set of equations re- For F a 0-form and u a 1-form, one obtains quired to define the field: instead of the curl and diver- ⋆ d (F u)= ⋆ (dF u) + ( 1)0 F (⋆ du) gence required for vector fields, form fields are unambigu- ∧ ∧ − ously defined by their derivative and their coderivative. curl (F u)= grad F u + F curl u. (47) → ×

6 For two 1-forms, this now gives , it means that one can always count (at least classi- cally)V the number of elementary charges localized within ⋆ d (u v)= ⋆ (du v) + ( 1)1 ⋆ (u dv) ∧ ∧ − ∧ and then obtain the total charge Q. This deﬁnes a div (u v) = (rot u) .v u. (rot v) . (48) 3-formV electric charge density ρ such that → × − As can be understood from these examples, the various identities of vector analysis simply come down to two Q = ρ. (52) ZV short formulas, (19) and (20). Other identities involving the mixed product and triple vector product can also be This relation is purely topological, as it does not rely on derived from the antisymmetry properties of the wedge lengths or angles. As ρ is a 3-form in R3, then necessarily product.11 Finally, all formulas involving integration of scalar and vector ﬁelds can also be retrieved from (24). dρ =0. (53) Indeed, when applied to a 0-form F (x), it simply gives the fundamental theorem of calculus This equation has no analog in the usual Maxwell-

∂F i Heaviside formulation of EM and comes from the fact dF = F i dx = F (b) F (a) . that electric charge is quantifiable. A first outcome arises Z[a,b] Z∂[a,b] → Z[a,b] dx − in virtue of the Poincarélemma: as ρ is a closed form, it Once a metric is given, one can draw connections be- is locally exact so that there exist a potential 2-form D tween the usual integral identities and the generalized such that Stokes theorem: the Stokes formula is obtained from (24) applied to a 1-form, whereas the Gauss formula comes ρ = dD. (54) from (24) applied to a 2-form. The Green identities are obtained by integrating (42) and then using the Stokes One recognizes here the exterior calculus version of the theorem40 Maxwell-Gauss equation (2). Hence, electric flux density translates into a 2-form D. dF (⋆dG)+ F (⋆∆G)= F ⋆dG (49) A second outcome appears when taking the time ZD ∧ ZD ∧ Z∂D ∧ derivative of ρ and using Schwarz’s theorem: as d(∂tρ)= F (⋆∆G) G(⋆∆F )= F (⋆dG) G(⋆dF ) (50) 0, in virtue of the Poincarélemma again, ∂tρ is exact ZD − Z∂D − so that there exists in any contractible domain a 2-form for F and G two 0-forms. Integration of the Leibniz potential j such that formula also provides an interesting result: ∂tρ = d( j). (55) − (du) v = u v ( 1)p u dv (51) ZD ∧ Z∂D ∧ − − ZD ∧ One recognizes here the continuity equation for the electric charge, which is usually built as a spin-off of (2) and where (24) was used to simplify the second integral: this (3). Note that as suggested by the usual continuity equa- expression generalizes integration by parts at arbitrary tion, the term in div j confirms that the electric current dimensions. density is indeed a 2-form. But there is more: combining (54) and (55), it comes that d (j + ∂tD) = 0, so that once again from the Poincarélemma, there exists a 1-form po- B. “Forms illuminate EM...” tential H such that

In this section, we reformulate electromagnetic theory j + ∂tD = dH. (56) in the language of differential forms and in so doing we discuss its axiomatics. Indeed, there are some redundan- This is the Maxwell-Ampère equation (3). cies within the set of equations, as can be seen from the To summarize, D and H are potentials associated to Maxwell-Thomson equation (4). In virtue of Schwarz’s sources ρ and j and they are introduced from topological theorem, div B = 0 div ∂tB = 0, so that ∂tB is a equations, which result from the fact that electric charges ⇒ purely solenoidal field. Hence, it expresses as the curl of a are quantifiable. Hence, they are relevant not only in “potential”, which is nothing else than Maxwell-Faraday matter, but also in free space and they “are microphysical equation (1). In order to clarify the foundations and the quantities of the same type likewise - in contrast to what structure of EM, we therefore proceed by trying to use the is stated is most textbooks”.42 minimal set of mathematical and physical objects at each The last set of equations is based on the absence of step: in particular, following the prescription of Obukhov magnetic monopoles, which can also be expressed from a and Hehl,33 we will try to postpone the use of the metric counting procedure on the 3D compact domain (also (which originates from another interaction, gravitation) coming from the idea that charges are quantifiable):V to as late as possible. Let us begin with a very general statement: electric charge is a quantifiable physical prop- QM =0= ρM , (57) erty. In practice, given a 3-dimensional compact domain ZV

7 so that the 3-form magnetic charge density identically for E and D are different because, fundamentally, these vanishes, ρM = 0. This is the physical content of the correspond to two different mathematical objects: forms Maxwell-Thomson equation (4) of different degrees. Now, considering (54), (56), (58) and (59), it appears dB =0, (58) that the number of unknowns exceeds the number of where B is the 2-form magnetic flux density. Hence, from equations. The role devoted to the constitutive relation- the Poincarélemma, once again deriving (58) with re- ships is precisely to fill this gap by providing additional spect to time and using the Schwarz’s theorem, there connections between the different parts of the electro- exists a 1-form potential E, the electric field, such that magnetic field (as the purpose of the present work is mainly pedagogical, we will not discuss the more tech- ∂tB = d( E). (59) nical cases of non-linear media and magneto-electric ef- − fects). As in 3D 1-forms and 2-forms are naturally This is the exterior algebra version of the Maxwell- mapped into each other by the Hodge dual operator, the Faraday equation (1) and its similarity with (55) high- only admissible translation of (5)-(6) is: lights its status of continuity equation for the magnetic flux. As suggested by the usual form of the Maxwell- D = ε0 ⋆ E, (61) Faraday equation, the term in curl E confirms that the B = µ0 ⋆ H. (62) electric field is indeed a 1-form. Moreover, forms and their associated vector fields do The constitutive relations are metric-coupled dual re- not share the same dimension. Combining (52) with (54) lations needed to close Maxwell’s equations. This and using the Stokes theorem gives the Gauss theorem highlights one extremely important point: as the four Maxwell’s equations can be built without the metric, it is the main purpose of constitutive relations to re- Q = D. (60) Z∂V veal how the EM field couples to geometry. Moreover, in free space, it is customary from the standpoint of A quick dimensional analysis reveals that the dimen- vector analysis to proceed to Gauss field identification sion of D is that of a charge [D] = C (but its com- 2 E = D, H = B (up to dimensional multiplicative con- ponents recover the regular dimension C/m ), [j] = C/s stants). As mentioned in section II.C, when the metric [H] = C/s,[B] = flux = J.s/C and [E] = J/C. We are is diagonal, this identity applies but anytime the met- now in position to address the problem of the transfor- ric has off-diagonal terms, this identification fails (see for mation laws (7)-(8) mentioned in the introduction. In- instance the mixing of EM components in the Sagnac deed, consider a general coordinate change x,y,z effect.43) Writing constitutive relations from exterior al- u (x,y,z) , v (x,y,z) , w (x,y,z) , one is looking{ for} the → { } gebra clarifies why Gauss identification may not apply: components of E in the new coordinate system. Straight- the action of Hodge operator - and consequently the con- forward algebra gives stitutive properties of vacuum - results from the geomet- a ∂x ∂x ∂x ric properties of space, so that anytime these latter are E = Eadx = Ex du + dv + dw + ... ∂u ∂v ∂w non-trivial, the links between the different fields are non- trivial as well (in 4D, vacuum constitutive properties are ∂x ∂y ∂z 44 = Ex + Ey + Ez du + ... only related to the conformal part of the metric ). ∂u ∂u ∂u In other words, vacuum has to be considered as a full- −T −T −T fledged medium with its own optical properties, as testi- = J 11 Ex + J 12 Ey + J 13 Ez du + ..., fied by the terminology of ε0, i.e. the permittivity of free which agrees with (7). On the other hand, for the 2-form space. Like in any material medium, the electric part of D, it comes that the field is determined by both E and D (resp. by H and

D= Dxdy dz + Dydz dx + Dzdx dy B for the magnetic part) that encompass information of ∧ ∧ ∧ different kinds and are not redundant quantities even in ∂y ∂z a b free space. This is in fact very close to Maxwell’s original = Dx a du b du + ... ∂u ∧ ∂u view on the field, as expressed explicitly in his Treatise: ∂y ∂z ∂z ∂y “We are thus led to consider two different quantities, the = Dx dv dw + ... + ... ∂v ∂w − ∂v ∂w ∧ magnetic force and the magnetic induction, both of which −1 −1 are supposed to be observed in a space from which the = Dx Com J dv dw + Com J dw du 11 ∧ 12 ∧ magnetic matter is removed.”.45 + Com J −1 du dv... + ... 13 ∧ J J J = Dx + Dy + Dz C. “... and EM illuminates forms.”18 det J11 det J21 det J31 dv dw + ..., ∧ The formula defining the Hodge dual makes its mean- which agrees with (8). The notation Com M refers to ing rather abstruse. Although ⋆u might appear as some the matrix of cofactors of M. The transformation laws kind of orthogonal complement of the p-form u, a deeper

8 way to understand Hodge star is related to an internal IV. VACUUM AS A POLARIZABLE symmetry of Maxwell’s equations, known as S-duality. MEDIUM In the absence of sources, these latter are left invariant when performing “cross-rotations” between the electric A. Classical vacuum: electrostatics in a 3,6 and magnetic parts of the field: uniform gravitational field E′ = E cos α B sin α (63) − B′ = E sin α + B cos α (64) As it is well known, Newton’s law of gravitation as r where for simplicity we use units for which c = 1. Fol- given by the potential Φ( ), may be obtained from any lowing a remark by Weber,46 this duality rotation can be metric theory of gravitation which, in the weak field limit, expressed in a more compact form from the Weber vector turns into F = E + iB (cf. Ref. 47): ds2 = (1 + 2Φ) dt2 + dℓ2, (75) F′ = eiα F = (cos α + i sin α) (E + iB) (65) − = (E cos α B sin α)+ i (E sin α + B cos α) . (66) where dℓ represents the space part of the metric. − Whittaker51 took the limit of Schwarzschild’s metric, α is a mixing angle related to the relative proportions of electric and magnetic fields. 2M 1 A direct connection with the Hodge operator can be ds2 = 1 dt2+ dr2+r2 dθ2 + sin2 θdϕ2 , − − r 1 2M made but at the cost of working in four dimensions. In − r this case, the electric 1-form E and the magnetic 2-form (76) B are gathered within the Faraday 2-form at a large distance from the gravitating centre such that the Newtonian potential gz appears in the time part of F = E dt + B. (67) − ∧ the metric. The result is Taking the 4D-exterior derivative and using (59) gives: 2 2 2 2 2 dz d F = dE dt + dB + ∂tB dt ds = (1+2gz) dt + dx + dy + , (77) 4 ∧ ∧ − 1+2gz = ( ∂tB) dt + ∂tB dt =0. (68) − ∧ ∧ Similarly, one defines the Maxwell 2-form as which thus corresponds to an attractive uniform Newto- nian gravitational field parallel to the z-axis. G = ⋆ F = D H dt (69) 4 − ∧ Considering an electric field weak enough such that its and this latter obeys the source equation: contribution to the gravitational field of this spacetime is negligible, Whittaker derived Laplace’s equation for the d G = dD + ∂tD dt dH dt 4 ∧ − ∧ scalar electrostatic potential in the absence of sources: = ρ + ∂tD dt (j + dD) dt ∧ − ∧ 2 2 2 = ρ j dt = J. (70) ∂ V ∂ V ∂ V − ∧ 2 + 2 +(1+2gz) 2 =0. (78) When applied on a 2-form in a 4D Lorentzian metric, ∂x ∂y ∂z 2 (33) becomes ⋆4 = 1, in a similar fashion to the imag- Now, one considers a uniformly charged infinite x-y inary unit i2 = 1.− Hence, one constructs the following − plane. Denoting by σ the uniform areal charge den- operator on a 2-form as sity, the corresponding volume charge density 3-form is exp(⋆4 α) = cos α + ⋆4 sin α. (71) ρ = σδ(z) dx dy dz. It can be easily seen that the displacement field∧ 2-form∧ Like the complex number exp(iθ), it obeys e⋆4 αe⋆4 β = e⋆4 (α+β) = e⋆4 βe⋆4 α (72) 1 π D = [Θ(z) Θ( z)]σ dx dy (79) ⋆4 e 2 = ⋆4. (73) 2 − − ∧ In Cartesian coordinates, applying the dual rotation obeys the Maxwell-Gauss law dD = ρ and the expected operator (71) to the Faraday 2-form leads to planar symmetry D( z) = D(z). Here Θ(z) is the − − ⋆4 α Heaviside step function. e F = cos α F + sin α⋆4 F Now, the 2-form D is related to E via D = ε0 ⋆4 (E dt) = cos α (Exdx dt + .. + Bxdy dz + ..) ∧ ∧ ∧ and symmetry tells us that the electric field 1-form writes + sin α⋆4 (Exdx dt + .. + Bxdy dz + ..) E = E˜ dz. We use (30) to compute ⋆ (dz dt)= dx dy. ∧ ∧ 4 ∧ ∧ = (Ex cos α Bx sin α) dx dt + .. Note that detgab = 1 for the metric (77). We then have − ∧ | | + (Ex sin α + Bx cos α) dy dz + .. (74) ∧ D = ǫ0Edx˜ dy (80) A comparison with (66) shows that the Hodge operator ∧ acts as a dual rotation between electric 2-forms and mag- and finally that netic 2-forms. In the next section, we explore how the Hodge dual applies on the EM field in the presence of 1 E = [Θ(z) Θ( z)]σ dz. (81) non trivial vacuua. 2ε0 − −

9 ∂2V On the other hand, Eq. (78) becomes (1 + 2gz) ∂z2 = with no summation on the j index. Ought to the sym- ∂2V metries of the problem, only the z-component of the field 0, or ∂z2 = 0, which gives the usual g-independent 0- form potential is relevant and therefore using (86) leads to 2 σ 11π 2 2 V (x,y,z)= z . (82) ε =1 4 4 α + O(α ) (89) −2ǫ0 | | − 8100 a me With E = dV one confirms the electric field 1-form Although the second term is probably far too small to given by Eq.− (81). be detectable, quantum vacuum behaves in principle as Comparing the expressions (79) and (81) reveals that a polarized medium, as is well-known in quantum elec- Dxy = ε0Ez, similarly to what happens in flat spacetime: trodynamics. As there are no gradient of the dielectric hence, there is no vacuum polarization due to the uniform constant, the refractive index of the Casimir vacuum is a gravitational field. constant given by the square root of (89). Note that as emphasized in Ref. 48, the dielectric constant being the same along the three axes, no birefringence is expected. B. Quantum vacuum: wave propagation in a Casimir vacuum V. CONCLUDING REMARKS The Casimir effect was first discovered by Hendrik Casimir in 1948 and it corresponds to a macroscopic force In the wake of Ref. 49, we hold exterior calculus to be due the quantum vacuum fluctuations. That force ap- a precious pedagogical tool for teaching electrodynam- pears when two parallel plates are separated from one ics. Although students are more comfortable with vector another by a vacuum layer of thickness a and it is gener- analysis, picturing the EM field as three-dimensional vec- ally attractive. Propagation of an electromagnetic field tors turns out to be misleading in some cases. Conversely, in the so-called Casimir vacuum can be adequately de- picturing the EM field as differential forms is probably scribed as experiencing the following geometry:48 less intuitive at first sight, but it leads to a deeper understanding of EM. d2(z) The main concepts of exterior calculus were introduced gµν = ηµν nµnν (83) − d1(z)+ d2(z) from undergraduate linear algebra and calculus. Con- trary to a widespread misconception, differential forms where ηµν = ( 1, +1, +1, +1) is the flat Minkowski met- − are not tiny quantities, they are dual vectors. The many ric and nµ = (0, 0, 0, 1) is the unit vector orthogonal to formulas used in vector analysis can be concisely written 3 the plates. The variable z is the coordinate x , which in terms of the wedge product and the exterior deriva- is measured in the direction perpendicular to the plates. tive. Exterior calculus also helps clarifying the axiomat- Functions d1(z) and d2(z) are determined from quantum ics of EM. From two postulates on electric charge (Eq. electrodynamics and a perturbative expansion up to or- (52)) and magnetic charge (Eq. (57)) or equivalently 2 48 der α gives: (58), the whole set of Maxwell’s equations can be re- covered without relying on the background geometry: in 1 2 d1(z)= + O(α ) (84) other words, their foundation is purely topological and −16π does not require any metric structure. Metric properties 11πα2 of the manifold appear on the other hand in the constitu- d2(z)= 4 4 (85) −64800 a me tive relations. D and E (resp. H and B) must indeed be with me the electron mass. In Cartesian coordinates, the understood as forms of different degrees that carry differ- Casimir vacuum metric can thus be approximated by ent information. They are connected by a duality relation that requires the knowledge of the background geometry, 100 0 as the Hodge star operator is explicitly depending on the −010 0 metric. g =   (86) µν 001 0 According to Sir W.L. Bragg, “the important thing  2 about science is not so much to obtain new facts as  0 001 Cα  − to discover new ways of thinking about them” (cited 2 4 4 where C = 11π /4050 a me. In the absence of magnetic in Ref. 50). For such a purpose, the relevance of ex- sources, Eq. (69) becomes D = ⋆4F . Now, using (30) we terior calculus goes far beyond the realm of electrody- have for the components of D namics. As a matter of fact, this formalism is likely to enlighten conceptual issues arising in thermodynamics, 1 1 D = F µν g ǫ = gµαgνβF g ǫ . (87) fluid mechanics...22 Forthcoming developments may also ρσ 2 µνρσ 2 αβ µνρσ p| | p| | involve the Yang-Mills theories: these are generalizations This leads to of Maxwell’s equations but built from larger symmetry groups (SU(N) instead of U(1)) and they require addi- j 00 jj D = √ g g g Ej (88) tional ingredients taken from exterior calculus. − | |

10 ACKNOWLEDGMENTS 17K.F. Warnick and P.H. Russer, “Differential forms and electromagnetic field theory”, Progress In Electromag- netics Research 148, 83–112 (2014) The authors thank T. Bertels for insightful discus- 18 sions. FM wishes to acknowledge the support of C.W. Misner, K.S. Thorne and J.A. LPCT/U.Lorraine and CNPq. Wheeler, Gravitation, (Freeman, San-Francisco, 1973) 19W. Thirring, A Course in Mathematical Physics 2: Classical Field Theory, (Springer-Verlag New York, 1978) 20 REFERENCES W.L. Burke, Applied differential geometry, (Cambridge University Press, Cambridge, 1985) 21 1 J. Baez and J.P. Muniain, Gauge fields, knots and F.J. Dyson, Why is Maxwell’s theory so hard to gravity, (World Scientific Publishing Company, Singa- understand?, (James Clerk Maxwell commemorative pore, 1994) booklet: James Clerk Maxwell Foundation, Edin- 22H. Flanders, Differential Forms with Applications to burgh, 1999) 2 the Physical Sciences, (Dover Publication Inc., New- M. Born and E. Wolf, Principles of Optics, (Cambridge York, 1989) University Press, Cambridge, 1999) 23 3 I.V. Lindell, Differential forms in J.D. Jackson, Classical Electrodynamics, (John Wiley electromagnetics, (John Wiley & Sons, 2004) & Sons, New-York, 1999) 24T. Frankel, The geometry of physics: an 4J. Gratus, P. Kinsler and M.W. McCall, “Maxwell’s introduction, (Cambridge University Press, Cam- (D, H) excitation fields: lessons from permanent mag- bridge, 2011) nets”, European Journal of Physics 40, 025203 (2019) 25 5 B. Felsager, Geometry, Particles and Fields, (Springer, D.J. Griffith, Introduction to Electrodynamics, (Cam- New-York, 1998) bridge University Press, Cambridge, 2017) 26 6 D. Bachman, A geometric approach to differential A. Zangwill, Modern electrodynamics, (Cambridge forms, (Birkhäuser, Boston, 2012) University Press, Cambridge, 2013) 27 7 T. Dray, Differential forms and the geometry of general J. Van Kranendonk and J.E. Sipe, “Foundations of the relativity, (CRC Press, Boca Raton, 2014) macroscopic electromagnetic theory of dielectric me- 28A. Zee, Quantum Field Theory in a Nutshell, (Prince- 15 dia”, Progress in Optics , 245–350 (1977) ton University Press, Princeton, 2003) 8A.J. Ward and JB. Pendry, “Refraction and geome- 29M. Göckeler and T. Schücker, Differential Geometry, try in Maxwell’s equations”, Journal of Modern Op- gauge theories, and gravity, (Cambridge University tics 43, 773–793 (1996) 9 Press, Cambridge, 1987) M. Kitano, Reformulation of electromagnetism with 30B. Schutz, Geometrical methods of mathematical differential forms, (Trends in ElectromagnetismFrom physics, (Cambridge University Press, Cam- Fundamentals to Applications; V. Barsan and R.P. bridge, 1980) Lungu Eds. IntechOpen, 2012) 31 10 M.D. Spivak, A comprehensive introduction to G.A. Deschamps, “Electromagnetics and differential differential geometry, (Publish or perish Inc., Hous- forms”, Proceedings of the IEEE 69, 676–696 (1981) 11 ton, 1970) N. Schleifer, “Differential forms as a basis for vector 32A. Robinson, Non-standard analysis, (Princeton land- analysiswith applications to electrodynamics”, Ameri- marks in mathematics and physics, Princeton Univer- can Journal of Physics 51, 1139–1145 (1983) 12 sity Press, Princeton, 1996) Y. Hoshino, “An elementary application of exterior dif- 33F.W. Hehl and Y.N. Obukhov, Foundations of classical ferential forms in quantum mechanics”, American Jour- electrodynamics: Charge, flux, and metric, (Birhäuser, 46 nal of Physics , 1148-1150 (1978) Boston, 2003) 13H. Amar, “Elementary application of differential 34W.L. Burke, Div, grad, curl are dead, 1995 forms to electromagnetism”, American Journal of https://people.ucsc.edu/ rmont/papers/Burke_DivGradCurl.pdf Physics 48, 252-252 (1980) ~ 14 D. Baldomir, “Differential forms and electromagnetism 35Y. Choquet-Bruhat and C. DeWitt-Morette, with M. in 3-dimensional Euclidean space R3”, IEE Proceed- Dillard-Bleick, Analysis, manifolds and physics, (North ings A-Physical Science, Measurement and Instrumen- Holland, Amsterdam, 1982) tation, Management and Education-Reviews 133, 139– 36R. Heras, “The Helmholtz theorem and retarded 143 (1986) fields”, European Journal of Physics 37, 065204 (2016) 15A. Bossavit, “Differential forms and the computation of 37S. Carroll, Spacetime and geometry: An introduction fields and forces in electromagnetism”, European Jour- to general relativity, (Pearson, Harlow, 2003) nal of Mechanics B 10, 474–488 (1991) 38 16 R.H. Price, “General relativity primer”, American F.L. Teixeira and W.C. Chew, “Differential forms, met- Journal of Physics 50, 300–329 (1982) rics, and the reflectionless absorption of electromag- 39M. Gasperini, Theory of Gravitational netic waves”, Journal of Electromagnetic Waves and interactions, (Springer, Milan, 2013) Applications 13, 665–686 (1999)

11 40M. Fecko, Differential geometry and Lie groups 47I. Bialynicki-Birula and Z. Bialynicka-Birula, “The for physicists, (Cambridge University Press, Cam- role of the Riemann–Silberstein vector in classical bridge, 2006) and quantum theories of electromagnetism”, Jour- 41D.H. Kobe, “Helmholtzs theorem revisited”, American nal of Physics A: Mathematical and Theoreti- Journal of Physics 54, 552-554 (1986) cal 46, 053001 (2013) 42F.W. Hehl and Y.N. Obukhov, A gentle introduction 48S. Liberati, S. Sonego S and M. Visser, “Scharn- to the foundations of classical electrodynamics: The horst effect at oblique incidence”, Physical Review meaning of the excitations (D, H) and the field D 63, 085003 (2001) strengths (E, B), arXiv preprint physics/0005084, 2000 49K.F. Warnick, R.H. Selfridge and D.V. Arnold, “Teach- 43E.J. Post, “Sagnac effect”, Reviews of Modern ing electromagnetic field theory using differential Physics 39, 475 (1967) forms”, IEEE Transactions on education 40, 53– 44V. Guillemin and S. Sternberg, Symplectic techniques 68 (1997) in physics, (Cambridge university Press, Cam- 50H.J. Eysenck, Genius: The natural history of bridge, 1990) creativity, (Cambridge University Press, Cam- 45J.C. Maxwell, A treatise on electricity and bridge, 1995) magnetism, (Clarendon press, Oxford, 1881) 51E. S. Whittaker, “On electric phenomena in gravita- 46H. Weber, Die partiellen Differential-Gleichungen der tional fields”, Proceedings of the Royal Society of Lon- mathematischen Physik nach Riemann’s Vorlesungen don. Series A 116, 720–735 (1927) bearb, (F. Vieweg & Sohn, 1919)