An Introduction to Space-Time Exterior Calculus

The basic concepts of exterior calculus for space-time multivectors are presented: interior and exterior products, interior and exterior derivatives, oriented integrals over hypersurfaces, circulation and flux of multivector fields. Two Stokes theorems relating the exterior and interior derivatives with circulation and flux respectively are derived. As an application, it is shown how the exterior-calculus space-time formulation of the electromagnetic Maxwell equations and Lorentz force recovers the standard vector-calculus formulations, in both differential and integral forms.


Introduction
Vector calculus has, since its introduction by J. W. Gibbs [1] and Heaviside, been the tool of choice to represent many physical phenomena. In mechanics, hydrodynamics and electromagnetism, quantities such as forces, velocities and currents are modeled as vector fields in space, while flux, circulation, divergence or curl describe operations on the vector fields themselves.
With relativity theory, it was observed that space and time are not independent but just coordinates in spacetime [2] (pp. 111-120). Tensors like the Faraday tensor in electromagnetism were quickly adopted as a natural representation of fields in space-time [3] (pp. 135-144). In parallel, mathematicians such as Cartan generalized the fundamental theorems of vector calculus, i.e., Gauss, Green, and Stokes, by means of differential forms [4]. Later on, differential forms were used in Hamiltonian mechanics, e. g. to calculate trajectories as vector field integrals [5] (pp. [194][195][196][197][198].
A third extension of vector calculus is given by geometric and Clifford algebras [6], where vectors are replaced by multivectors and operations such as the cross and the dot products subsumed in the geometric product. However, the absence of an explicit formula for the geometric product hinders its widespread use. An alternative would have been the exterior algebra developed by Grassmann which nevertheless has received little attention in the literature [7]. An early work in this direction was Sommerfeld's presentation of electromagnetism in terms of six-vectors [8].
We present a generalization of vector calculus to exterior algebra and calculus. The basic notions of space-time exterior algebra, introduced in Section 2, are extended to exterior calculus in Section 3 and applied to rederive the equations of electromagnetism in Section 4. In contrast to geometric algebra, our interior and exterior products admit explicit formulations, thereby merging the simplicity and intuitiveness of standard vector calculus with the power of tensors and differential forms.

Exterior Algebra
Vector calculus is constructed around the vector space R 3 , where every point is represented by three spatial coordinates. In relativity theory the underlying vector space is R 1+3 and time is treated as a coordinate in the same footing as the three spatial dimensions. We build our theory in space-time with k time dimensions and n space dimensions. The number of space-time dimensions is thus k + n and we may refer to a (k, n)-or (k + n)space-time, R k+n . We adopt the convention that the first k indices, i.e., i = 0, . . . , k − 1, correspond to time components and the indices i = k, . . . , k + n − 1 represent space components and both k and n are non-negative integers. A point or position in this space-time is denoted by x, with components {x i } k+n−1 i=0 in the canonical basis (1) Given two arbitrary canonical basis vectors e i and e j , then their dot product in space-time is (2) For convenience, we define the symbol ∆ ij = e i · e j as the metric diagonal tensor in Minkowski space-time [2] (pp. 118-120), such that time unit vectors e i have negative norm ∆ ii = −1, whereas space unit vectors e i have positive norm ∆ ii = +1. The dot product of two vectors x and y is the extension by linearity of the product in Equation (2), namely

Grade, Multivectors, and Exterior Product
In addition to the (k + n)-dimensional vector space R k+n with canonical basis vectors e i , there exist other natural vector spaces indexed by ordered lists I = (i 1 , . . . , i m ) of m non-identical space and time indices for every m = 0, . . . , k + n. As there are k+n m such lists, the dimension of this vector space is k+n m . We shall refer to m as grade and to these vectors as multivectors or grade-m vectors if we wish to be more specific. A general multivector can be written as v = where the summation extends to all possible ordered lists with m indices. If m = 0, the list is empty and the corresponding vector space is R. The direct sum of these vector spaces for all m is a larger vector space of dimension k+n m=0 k+n m = 2 k+n , the exterior algebra. In tensor algebra, multivectors correspond to antisymmetric tensors of rank m. In this paper, we study vector fields v(x), namely multivector-valued functions v varying over the space-time position x.
The basis vectors for any grade m may be constructed from the canonical basis vectors e i by means of the exterior product (also known as wedge product), an operation denoted by ∧ [9] (p. 2). We identify the vector e I for the ordered list I = (i 1 , i 2 , . . . , i m ) with the exterior product of e i1 , e i2 , . . . , e im : In general, we may compute the exterior product as follows. Let two basis vectors e I and e J have grades m = |I| and m = |J|, where |I| and |J| are the lengths of the respective index lists. Let (I, J) = {i 1 , . . . , i m , j 1 , . . . , j m } denote the concatenation of I and J, let σ(I, J) denote the signature of the permutation sorting the elements of this concatenated list of m + m indices, and let ε(I, J) denote the resulting sorted list, which we also denote by I + J. Then, the exterior product e I of e J is defined as e I ∧ e J = σ(I, J)e ε(I,J) .
The exterior product of vectors v and w is the bilinear extension of the product in Equation (6), Since permutations with repeated indices have zero signature, the exterior product is zero if m + m > k + n or more generally if both vectors have at least one index in common. Therefore, the exterior product is either zero or a vector of grade m + m . Further, the exterior product is a skew-commutative operation, as we can also write Equation (6) as e I ∧ e J = (−1) |I||J| e J ∧ e I . At this point, we define the dot product · for arbitrary grade-m basis vectors e I and e J as where I and J are the ordered lists I = (i 1 , i 2 , . . . , i m ) and J = (j 1 , j 2 , . . . , j m ). As before, we extend this operation to arbitrary grade-m vectors by linearity. Finally, we define the complement of a multivector. For a unit vector e I with grade m, its Grassmann or Hodge complement [10] (pp. 361-364), denoted by e H I , is the unit (k + n − m)-vector where I c is the complement of the list I, namely the ordered sequence of indices not included in I. As before, σ(I, I c ) is the signature of the permutation sorting the elements of the concatenated list (I, I c ) containing all space-time indices. In other words e I c is the basis vector of grade k + n − m whose indices are in the complement of I. In addition, we define the inverse complement transformation as We extend the complement and its inverse to general vectors in the space-time algebra by linearity.

Interior Products
While the exterior product of two multivectors is an operation that outputs a multivector whose grade is the addition of the input grades, the dot product takes two multivectors of identical grade and subtracts their grades, yielding a zero-grade multivector, i.e., a scalar. We say that the exterior product raises the grade while the dot product lowers the grade. In this section, we define the left and right interior products of two multivectors as operations that lower the grade and output a multivector whose grade is the difference of the input multivector grades.
As always, we start by defining the operation for the canonical basis vectors. Let e I and e J be two basis vectors of respective grades |I| and |J|. The left interior product, denoted by , is defined as e I e J = ∆ I,I σ ε(I, J c ) c , I e ε(I,J c ) c .
If I is not a subset of J, that is when there are elements in I not present in J, e. g. for |I| > |J|, the signature of the permutation sorting the concatenated list ε(I, J c ) c , I is zero as there are repeated indices in the list to be sorted, and the left interior product is zero. Otherwise, if I is a subset of J, the permutation rearranges the indices in J in such a way that the last |I| positions coincide with I and ε(I, J c ) c represents the first |J| − |I| elements in the rearranged sequence, that is ε(I, J c ) c = J \ I. The right interior product, denoted by , of two basis vectors e I and e J is defined as As with the left interior product, if J is a subset of I, ε(I c , J) c = I \ J then the permutation rearranges the indices in I so that the first |J| positions coincide with J, otherwise the right interior product is zero.
In general, we have that e I e J = e J e I (−1) |I|(|J|−|I|) , as verified in Appendix A.1. We note that these interior products are not commutative, unless either |J| − |I| or |I| is an even number, e. g. when |I| = |J|, in which case both interior products coincide with the dot product of the two vectors. The interior products may therefore be seen as generalizations of the dot product.
As with the dot and the exterior products, the value of the interior products does not depend on the choice of basis and we may thus compute the left interior product of two vectors v and w as v w = and a similar expression holds for the right interior product v w. Both are grade-lowering operations, as the left (resp. right) interior product is either zero or a multivector of grade m − m (resp. m − m ).
The interior products are not independent operations from the exterior product, as they can be expressed in terms of the latter, the Hodge complement and its inverse (proved in Appendix A.2): e I e J = e H −1 If u and v are 1-vectors and w is an r-vector, then we have the following expression as proved in Appendix A.3. This expression can be seen as a generalization of the vectorial expression in the vector space R 3 , i.e., a k = 0, n = 3 space-time. This fact is built of the realization that the cross product between two vectors v and w can be expressed in the following alternative ways Whenever it holds that I ⊆ J, the interior and exterior products are related by the following: Having introduced the basic notions of space-time exterior algebra, the next section focuses on operations with elements in the exterior algebra, namely integrals and derivatives of vector fields.

Oriented Integrals
Integrals are, together with derivatives, the fundamental mathematical objects of calculus. For example, operations on vectors fields lying in exterior algebra such as the flux and the circulation are expressed in terms of integrals over high-dimensional geometric objects. The integral of an m-graded vector field v over a hypersurface V m of the same dimension, denoted as is the limit of the Riemann sums for the dot product d m x · v over points in the hypersurface, where d m x is an m-dimensional infinitesimal vector element. For any = 0, . . . , k + n, the infinitesimal vector element d x is given by the sum of all possible differentials for -dimensional hypersurfaces in a (k, n) space-time, and is represented in the canonical basis as where for a given list I = (i 1 , . . . , i ) each differential is given by dx I = dx i1 · · · dx i . As in traditional calculus, the integral in Equation (21) exhibits coordinate invariance, while the integrand d m x · v is regarded as an oriented object. Orientation is well defined for integrals along a curve from one point to another, or integrals over a surface oriented at the direction of the normal to the surface. Switching the extreme points of the curve, or taking the opposite direction of the normal would induce a change of sign in the line and surface integrals. In our generalization of vector calculus, a positive orientation is implicit in the ordering of the canonical basis. The skew-symmetry property of the exterior product Equation (6) may introduce sign changes to compensate an eventual change of orientation after changes of coordinates such as permutations of the space-time components.
For a given hypersurface V m , a convenient transformation for solving the integral in Equation (21) is one such that, at a given point x in the hypersurface, the infinitesimal vector element d m x has one component that is tangent to the hypersurface at that point. Let e be a unit m-graded vector parallel to V m at point x, and let e 0 , . . . , e k+n−1 form an orthonormal basis of R k+n such that e = e k+n−m ∧ · · · ∧ e k+n−1 for the given point x in V m . This change of coordinates from the canonical basis to the new basis is described by a unitary matrix U , dependent on x, and that satisfies e 0 ∧ · · · ∧ e k+n−1 = det(U ) e 0 ∧ · · · ∧ e k+n−1 . (23) Being a unitary matrix, the determinant of U is ±1. Assuming an orientation-preserving change of coordinates, that is det(U ) = 1, the infinitesimal vector element in Equation (22) for = m can be expressed as where ⊥= {0, . . . , k + n − m − 1} is the set of indices for the unit vectors in the new basis orthogonal to V m . Since all elements in the summation in Equation (24) have at least one differential element lying outside the integration hypersurface, their integrals vanish and therefore In analogy to e , a multivector of grade m, we define a unit (k + n − m)-grade vector e ⊥ normal to V m at point x such that e ⊥ ∧ e = e 0 ∧ · · · ∧ e k+n−1 . From Equation (10), we see that one such normal multivector with the correct orientation is For the common spaces considered in vector calculus, R 2 and R 3 , and according to Equation (23), orientationpreserving changes of coordinates must respectively satisfy e ⊥ ∧ e = e 0 ∧ e 1 and e ⊥ ∧ e = e 0 ∧ e 1 ∧ e 2 , where e ⊥ is the basis element normal to V m . These two equalities turn out to describe the counterclockwise (resp. right-hand rule) orientation when e ⊥ conventionally points outside an integration path for R 2 (resp. a surface for R 3 ) [5] (pp. 184-185). Building on the concepts and operations of circulation and flux in vector calculus, the right and left interior products lead to general definitions of circulation and flux of multivector fields in exterior algebra along and across hypersurfaces of arbitrary number of dimensions.

Circulation and Flux of Multivector Fields
Definition 1. The circulation of a vector field v(x) of grade m along an -dimensional hypersurface V , denoted by C(v, V ), is given by Expressing the vector field in the canonical basis and using the definition of d x in Equation (22), the circulation can be specified in some cases of interest. For = m, the circulation reads For instance, for = m = 1 and R n , this formula recovers the definition the circulation of a vector field along a closed path with the appropriate orientation. Alternatively, using Equation (25), we note that v is integrated along the direction of e , tangential to the hypersurface, in an orientation-preserving change of coordinates, that is Intuitively, the circulation Equation (27)  Definition 2. The flux of a vector field v(x) of grade m across an -dimensional hypersurface V , denoted by F(v, V ), is given by Expressing both v and d x in the canonical basis, and using the inverse Hodge operation in Equation (10), the flux in the special case of = k + n − m can be written as As an example in R 3 , the flux of a vector field v through a surface V 2 reads The right-hand side of Equation (32) is a conventional surface integral, upon the identification of I,i / ∈I dx I σ(i, I)e i as an infinitesimal surface element dS.
Alternatively, using the analogous of Equation (25) for the differential vector element d x H −1 , the equivalent to Equation (29) for the flux is This equation implies that v is integrated along a normal component to the hypersurface since e H −1 is a multivector of grade k + n − orthogonal to V . Intuitively, the flux Equation (30) measures the magnitude of the multivector field crossing the hypersurface. In general, the flux is a vector of grade (m + − n − k) if ≥ k + n − m and zero otherwise. For instance, if = k + n, the flux of v over an (k + n)-dimensional hypersurface V k+n gives the integral of v over V k+n , an extension of the volume integral to R k+n , where we used the relation 1 H = e i1,...,i k+n , implying that d k+n x H −1 = dx i1,··· ,i k+n , and that 1 v = v.

Exterior and Interior Derivatives
In vector calculus, extensive use is made of the nabla operator ∇, a vector operator that takes partial space derivatives. For instance, operations such as gradient, divergence or curl are expressed in terms of this operator. In our case, we need the generalization to (k, n) space-time to the differential vector operator ∂, defined as For a given vector field v of grade m, we define the exterior derivative of v as ∂ ∧ v, namely The grade of the exterior derivative of v is m + 1, unless m = k + n, in which case the exterior derivative is zero, as can be deduced from the fact that all signatures are zero.
In addition, we define the interior derivative of v as ∂ v, namely The grade of the interior derivative of v is m − 1, unless m = 0, in which case the interior derivative is zero, as implied by the fact that the grade of ∂ is larger than the grade of v. Using Equation (16) with u = ∂ and assuming that v and w are 1-vectors, we obtain a generalization of Leibniz's product rule The formulas for the exterior and interior derivatives allow us express some common expressions in vector calculus. For a scalar function φ, its gradient is given by its exterior derivative ∇φ = ∂ ∧ φ, while for a vector field v, its divergence ∇ · v is given by its interior derivative ∇ · v = ∂ v. From Equation (16) we further observe that for a scalar function φ we recover the relation In addition, for a vector fields v in R 3 , taking into account Equation (18) then the curl can be variously expressed as This formula allows us to write the curl of a vector field ∇ × v in terms of the exterior and interior products and the Hodge complement, while generalizing both the cross product and the curl to grade-m vector fields in space-time algebras with different dimensions. Moreover, from Equation (16) we can recover for r = 1 the well-known formula for the curl of the curl of a vector, It is easy to verify that the exterior derivative of an exterior derivative is zero, as is the interior derivative of an interior derivative, that is for any vector field v, we have that In regard to the vector space R 3 , and using Equation (18), these expressions imply the well-known facts that the curl of the gradient and the divergence of the curl are zero:

Stokes Theorem for the Circulation
In vector calculus in R 3 , the Kelvin-Stokes theorem for the circulation of a vector field v of grade 1 along the boundary ∂V 2 of a bidimensional surface V 2 relates its value to that of the surface integral of the curl of the vector field over the surface itself. In the notation used in the previous section, the surface integral is the flux of the curl of the vector field across the surface and this theorem reads Taking into account the identity ∇ × v = (∇ ∧ v) H −1 in Equation (40), we rewrite the right-hand side in Equation (46) as where we used that u · w = u H −1 · w H −1 = u H · w H for vectors u, w. The flux of the curl of the vector field across a surface is also the circulation of the exterior derivative of the vector field along that surface. The generalized Stokes theorem for differential forms [4] (p. 80) allows us to extend the Kelvin-Stokes theorem to multivectors of any grade m as we do in the following theorem.
Theorem 1. The circulation of a grade-m vector field v along the boundary ∂V of an -dimensional hypersurface V is equal to the circulation of the exterior derivative of v along V : (48) ¶ As hinted at above, the role of the vector curl in the right-hand side of Equation (46) is played by the exterior derivative in this generalized theorem.
Proof. We start by stating the generalized Stokes Theorem for differential forms [4] (pp. 80) where ω is a differential form and dω its exterior derivative, represented by the operator Expressing the circulations in Equation (48) by means of the integrals in Equation (27), we obtain In the integral in the left-hand side of Equation (51), the integrand is a differential form ω = d −1 x v. After expanding the interior product using the definitions of d −1 x and v we obtain We next write down the integrand in the right-hand side of Equation (51) dx K e K Im j ∈I ∆ j,j ∂ j v I σ(j, I) e ε(j,I) = K ,Im : ε(j,I)⊆K j ∈I ∆ j,j ∂ j v I dx K σ(j, I)∆ ε(j,I),ε(j,I) σ(ε(j, I), ε(K c , ε(j, I)) c ) e ε(K c ,ε(j,I)) c = K ,Im : ε(j,I)⊆K +1 j ∈I and verify that it coincides with exterior derivative in Equation (53). As the set of m indices I m is included in the sets J −1 or K in Equation (53) or (54), we may write K = ε(J −1 , j) for some j / ∈ J −1 . Then, we obtain the following chain of equalities for the basis elements in Equations (53) and (54): Therefore, and using that ε(J c \ {j}, ε(j, I)) c = J \ I, we can write Equation ( To prove Equation (57) we exploit that the σ are permutation signatures and that the signature of the composition of permutations is the product of the respective signatures. We proceed with the help of a visual aid in Figure 1, which depicts the identity between two different ways of sorting the concatenated list (j, I, J \ I). On the left column we first sort the list (I, J \ I) to obtain J and then sort the list (j, J). On the right column, we first sort the list (j, I) and then the list (j + I, J \ I). This proves Equation (57) and the theorem. Finally, we note that, had we defined the circulation with the left interior product, we would have got an incompatible relation in Equation (57), which could not be solved.

Stokes Theorem for the Flux
In vector calculus in R 3 , the Gauss theorem relates the volume integral of the divergence of a vector field v over a region V 3 to the surface integral of the vector field over the region boundary ∂V 3 . In the notation used in previous sections, and taking into account that both the surface integral and the volume integral can be expressed as fluxes for R 3 , this theorem reads Making use of the identity ∇ · v = ∇ v, we can rewrite the right-hand side in Equation (58) as In other words, the Gauss theorem relates the flux of the interior derivative of a vector field v across a region V 3 to the flux of the vector field itself across the region boundary ∂V 3 . The generalized Stokes theorem for differential forms allows us to extend the Gauss theorem to multivectors of any grade m as we do in the following theorem.
Theorem 2. The flux of a grade-m vector field v across the boundary ∂V of an -dimensional hypersurface V is equal to the flux of the interior derivative of v across V : (60) ¶

Proof. Expressing the fluxes in Equation (60) by means of the integrals in Equation (30), we obtain
As in the proof of Theorem 1, we apply the Stokes theorem for differential forms in Equation (49) upon the identifications ω with d −1 x H −1 v and dω with d x H −1 (∂ v). First, for ω, we get Now, taking the exterior derivative of Equation (62), we obtain dω = This quantity should be equal to d x H −1 (∂ v) in the right-hand side of Equation (61), which we expand as We first consider the sets in the summations in the alternative expressions for dω, Equations (63) and (64). Since J c contains j and is a subset of I, but K c does not contain j and is also a subset of I (with j ∈ I), then we can assert that K = J ∪ {j} so that the conditions in the summations are equivalent. The basis elements coincide and so do the differentials and derivatives, and it remains to verify the identity With the definition L = I \ J c , and expressed in terms of j, J, and L, this condition gives Multiplying both sides of the equation by σ(J c , J), σ(J c \ j, J + j) and σ(J c \ j, j), and taking into account that the square of a signature is +1, we obtain We start by simplifying Equation (67) by noting that with help of the visual aid in Figure 2. The permutations on the left column first merge (J c \ j) with j and then the resulting J c with L. Similarly, on the right column, we start with L, (J c \ j) and {j}, then concatenate (L, J c \ j) and then add j, getting the same result as the left column. Therefore, we have reduced Equation (67) to the simpler form which we prove with the aid depicted in Figure 3. On the left column, j and J are first merged and then the concatenation (J c \ j, J + j) gives the sorted ε(J c , J). On the right column, after sorting (J c \ j) with j, merging it with J leads to the same final sequence. density vector J. In the appropriate units, the bivector field F can be decomposed as F = F E + F B , where F E contains the electric-field E time-space components and F B contains the space-space components for the magnetic field B. Similarly, the current density depends on the charge density ρ and the spatial current density j. More specifically, Here the Hodge complement acts only on the space components, and B H = B H −1 . The bivector field F is closely related to the Faraday tensor, a rank-2 antisymmetric tensor. Maxwell equations, in their differential form, constrain the divergence of the electric and the magnetic field, Equations (72) and (73), respectively, and the curl of E and B, namely Equations (74) and (75) [11] (p. 4-1).
We refer to Equations (73) and (74) as homogeneous Maxwell equations and to Equations (72) and (75) as inhomogeneous Maxwell equations, as they include the fields and the sources given by charge and current densities. In exterior-calculus notation, both pairs of equations can be combined into simple multivector equations, where ∂ is the differential operator ∂ = −∂ 0 e 0 + ∇ for k = 1 and n = 3. As a consistency check, note that the wedge product raises the grade of F, and the zero in Equation (76) is the zero trivector; also, as the left interior product lowers the grade of F, both sides of Equation (77) relate space-time vectors. Next to Maxwell equations, the Lorentz force density f f f characterizes, after integrating over the appropriate region, the force exerted by the electromagnetic field upon a system of charges described by the charge and current densities ρ and j [11] (pp. 13-1-13-3), In relativistic form, the Lorentz force density becomes a four-dimensional vector f [2] (pp. 153-157). The time component of this vector is j · E, the power dissipated per unit of volume, or after integrating over the appropriate region, the rate of work being done on the charges by the fields. In exterior-calculus notation, the Lorentz force density vector can be computed as a left interior product, namely

Equivalence of the Lorentz Force Density
In this section, we prove that Equation (79) indeed recovers the relativistic Lorentz force density by verifying that its components in both vector-calculus and exterior-calculus coincide. From the definitions of J and F, and using the distributive property of the interior product, we get Some straightforward calculations give e 0 (e 0 ∧ E) = E, e 0 B H = 0, and j (e 0 ∧ E) = e 0 j · E. In addition, the formula for the left interior product in Equation (18) gives j B H −1 = (j ∧ B) H −1 = j × B, where the cross product is only valid for three dimensions. With these calculations, we obtain namely, a time-component j · E and a spatial component equal to the Lorentz force density ρE + j × B.

Equivalence of the Differential Form of Maxwell Equations
In this section, we prove that Equation (76) indeed recovers the homogeneous Maxwell equations and that Equation (77) recovers the inhomogeneous Maxwell equations. First, we observe that the exterior derivative ∂ ∧ F gives a trivector with 4 components, while the homogeneous Maxwell equations are a scalar, Equation (73), and a vector, Equation (74). We shall verify that the scalar equation turns out to be given by the trivector component e 123 of ∂ ∧ F, while the vector equation is given by the trivector components e 012 , e 013 , and e 023 of the exterior derivative.
We evaluate the exterior derivative ∂ ∧ F using the decomposition of F in Equation (71), We move on to the inhomogeneous Maxwell equations. We compute the interior derivative ∂ F, We can verify that the scalar equation turns out to be given by the vector component e 0 of ∂ F, while the spatial vector equation is given by the spatial vector components e 1 , e 2 , and e 3 of ∂ F. Indeed, if we match this expression with the current density vector J, then the time component e 0 of ∂ F gives Equation (72). Selecting the space components of ∂ F, the differential equation is which, using the relation ∇ B H = ∇ × B can be written as Equation (75).

Equivalence of the Integral Form of Maxwell Equations
After studying the exterior-calculus differential formulation of Maxwell equations, we recover the standard integral formulation. Applying the Stokes Theorem 1 to Equation (76), we find that the circulation of the bivector field F along the boundary of any three-dimensional space-time volume V 3 is zero: At this point, Equation (86) is a scalar equation and we obtain the pair of homogeneous Maxwell equations by considering two different hypersurfaces V 3 . First, let the domain V 3 = V contain only spatial coordinates. There are no tangential components to V with time indices and the contribution of F E to the circulation of F over ∂V 3 in Equation (86) is zero, i.e., Using that u · w = u H −1 · w H −1 for any vectors u, w, and therefore d 2 and the definition F B = B H , the integral in the right-hand side of Equation (87) becomes where we used Equation (32) to write the last surface integral. Substituting Equation (88) back into Equation (86) gives the Gauss law for the magnetic field [11] (pp. 1-5-1-9). Let now V 3 be a time-space domain (t 0 , t 1 ) × S, where S is a two-dimensional spatial surface. With no real loss of generality we assume that S lies on the e 1 ∧ e 2 plane. Its boundary ∂V 3 is the union of the sets (t 0 , t 1 ) × ∂S, t 0 × S and t 1 × S. For the first set, we choose e ⊥ as the vector normal to ∂S pointing outwards on the plane defined by S and e = e 0 ∧ e ∂S , where e ∂S is a vector tangent to ∂S with a counterclockwise orientation, so that e ⊥ ∧ e = −e 012 . Further, since e is a time-space bivector, the contribution of F B to the circulation of F over this first set in Equation (86) is zero, and (t0,t1)×∂S Writing the differential vector as d 2 x = dt dxe 0x , parameterizing the line integral over the boundary ∂S by the variable x with unit vector e 0x , and using that e 0x · F E = −e x · E and therefore dx e 0x · F E = − dx · E, the integral of the right-hand side of Equation (89) becomes For the second and third sets the normal vector to the integration surface pointing outwards are e ⊥ = −e 0 and e ⊥ = e 0 respectively. Since e is a space-space bivector in both cases, then the contribution of F E to the circulation is zero. We express the circulations of F B as fluxes of B and surface integrals as done in Equation (88).Using these observations the integral for the circulation of F over these two sets in Equation (86) is given by Combining Equations (90) and (91) in Equation (86) we recover the integral over time of the so called Faraday law [11] (pp. 17-1-17-2). Equivalently, taking the time derivative recovers the usual Faraday law, namely In regard to the inhomogeneous Maxwell equations, applying the Stokes Theorem 2 to Equation (77), we find that the flux of the bivector field F across the boundary of any three-dimensional space-time volume is equal to the flux of the current density J across the three-dimensional space-time volume: As with the homogeneous Maxwell equations, the scalar Equation (93) yields the inhomogeneous Maxwell equations by considering two different hypersurfaces V 3 . First, let the integration domain V 3 be a spatial volume V . Since there are no normal components to V with space indices only, the contribution of F B to the flux is zero so that Equation (93) becomes From the definition of inverse Hodge complement in Equation (10), we write the differential vectors Plugging these expressions in Equation (94), using the definitions of F E and J, and computing the dot products on both sides of the equality, we obtain that Equation (94) simplifies as In Equation (98) we used that σ(0i, I) = σ(i, I) and in Equation (99) we used that Since the Hodge complement is over space, the result is a surface integral with positive orientation as in Equation (32). We recovered in Equation (99) the Gauss law for the electric field [11] (pp. 4-7-4-9).
For V 3 = (t 0 , t 1 ) × S where S is a two-dimensional surface lying on the e 1 ∧ e 2 plane, the boundary ∂V 3 is the union of the sets (t 0 , t 1 ) × ∂S, t 0 × S and t 1 × S. For the first set, since d 2 x H −1 has no time components, the contribution of F E to this set is zero, that is As in the homogeneous case, we choose e ⊥ as the vector normal to ∂S pointing outwards on the plane defined by S and e = e 0 ∧ e ∂S , where e ∂S is a vector tangent to ∂S with a counterclockwise orientation, such that e ⊥ ∧ e = −e 012 introduces a change of sign. Expressing d 2 x H −1 and F B in the canonical basis, defining I = (0, i) so that I c contains only space indexes, and using that e I c · e i c = 1 and σ(I c , I) = σ(I c , 0, i) = σ(0, I c , i) = σ(I c , i), we obtain that Equation (100) simplifies to For the second and third sets we respectively choose e ⊥ = −e 0 and e ⊥ = e 0 pointing outside V 3 , implying that the contribution of F B is zero for this set as the inverse Hodge complement of e ⊥ is a space vector. Expressing d 2 x H −1 in Equation (95) and using similar steps as in Equations (97)-(99), the left-hand side of Equation (93) over these two sets is given by Finally, for the right-hand side of Equation (93), we choose e ⊥ as the vector normal to V 3 pointing outside. Since e = e 012 implies that e ⊥ ∧ e = −e 0123 , we obtain that We have thusly recovered the integral form of the Ampere-Maxwell equation [11] (p. 18-1-18-4) integrated over the time interval (t 0 , t 1 ) by combining Equations (101)-(103) into Equation (93), that is

Summary
In this paper, we aimed at showing how exterior calculus provides a tool merging the simplicity and intuitiveness of standard vector calculus with the power of tensors and differential forms. Set in the context of a general space-time algebra with multiple space and time components, we provided the basic concepts of exterior algebra and calculus, such as multivectors, wedge product and interior products, with a distinction between left and right products, Hodge complement, and exterior and interior derivatives. While a space-time with multiple time coordinates leads to several issues from the physical point of view [12], we did not deal with these problems as this paper focuses on the mathematical constructions. We also defined oriented integrals, with two important examples being the flux and circulation of grade m-vector fields as integrals of the normal and tangent components of the field to a hypersurface respectively. These operations extend the standard circulation of a vector field as a line integral and the flux of a vector field as a surface integral in three dimensions to any number of dimensions and any vector grade.
Armed with the theory of differential forms, we proved two exterior-calculus Stokes theorems, one for the circulation and one for the flux, that generalize the Kelvin-Stokes, Gauss and Green theorems. We saw that the flux of the curl of a vector field in three dimensions across a surface is also the circulation of the exterior derivative of the vector field along that surface. In exterior calculus, these Stokes theorems hold for any number of dimensions and any vector grade and are simply expressed in terms of the exterior and interior derivatives for the circulation and flux respectively.
As an application of our tools, we showed how to recover the classical laws of electromagnetism, Maxwell equations and Lorentz force, from a exterior-calculus formalism in relativistic space-time with one temporal and three spatial dimensions. The electromagnetic field is described by a bivector field with six components, closely related to Faraday's antisymmetric tensor, containing both electric and magnetic fields. The differential form of Maxwell equations relates the exterior derivative of the bivector field with the zero trivector and the interior derivative of the field with the current density vector. In the integral form, these equations correspond to the statements that the circulation of the bivector field along the boundary of any three-dimensional space-time volume is zero, and that the flux of the bivector field across the boundary of any three-dimensional space-time volume is equal to the flux of the current density across the same space-time volume.

A Proofs of product identities
In this appendix, we verify the relations about interior products introduced in Section 2.

A.1 Relation between left and right interior products
We now prove the formula e I e J = e J e I (−1) |I|(|J|−|I|) , relating left and right interior products. For two lists I and J, we have e I e J = ∆ I,I σ J \ I, I e J\I , e J e I = ∆ I,I σ I, J \ I e J\I , where we assumed that I ⊆ J with no loss of generality and used that ε(I, J c ) c = J \ I in this case. The only difference between the expressions lies in the signatures, that are related by setting A = J \ I and B = I in the following lemma.