1

Electromagnetism in Curved Spacetimes: Coupling of the Doppler and Gravitational Redshifts

Cameron R. D. Bunney, Gabriele Gradoni, Member, IEEE

Abstract—We present the basic prerequisites of electromag- tensor, we are able to simplify Einstein’s field equations to netism in flat spacetime and provide the description of elec- the Maxwell-Einstein field equations. We then discuss the tromagnetism in terms of the Faraday tensor. We generalise Lorenz gauge and the canonical wave equation in this theory, electromagnetic theory to a general relativistic setting, introduc- ing the Einstein field equations to describe the propagation of yielding the de Rham wave equation. We provide a solution to electromagnetic radiation in curved space-time. We investigate this, as in [14], through the geometrical optics approximation. gravitational redshift and derive formulae for the combined The discussion again moves to discuss electrodynamics in effect of gravitational redshift and the Doppler shift in curved curved spacetimes. Key concepts such as the four-velocity spacetime. and proper time are explained, leading to the key dynamical Index Terms—Dipole Antennas, Doppler Effect, Electromag- equation of motion. Changing from Euclidean space to a curved netic Propagation, Formal Concept Analysis, Lorentz Covariance spacetime, the geodesic equation is found in an approximated form. We investigate the dipole radiation in [7] as motivation to study the propagation of radiation in curved spacetimes, I.INTRODUCTION giving an analysis of the geodesics that light and therefore The general relativistic treatment of electromagnetism is electromagnetic waves travels along in a general, spherically typically developed by promoting Minkowski spacetime to a symmetric spacetime. We further the study of null ray solutions general, possibly curved, spacetime, which is more physically to the calculation for non-radial null rays. This investigation relevant since matter curves spacetime in the universe. allows us to delve into gravitational redshift, re-deriving a well- The current pedagogical works on relativistic electrodynam- known approximation [3], the redshift of light sent between two ics come from two points of view: the first is for those with a stationary observers. We then derive results for the redshift of background in electromagnetism and deals with introductory light experiencing both gravitational redshift and the Doppler special relativity; the second is for those with a background in shift, sent between two observers which may be in 1) radial gravity, which studies curved spacetimes and general relativity, motion, or 2) a combination of both radial and azimuthal with examples involving electromagnetism. The aim of this motion. article is to review these two areas as a joint theory for the applied electromagnetism audience. II.ELECTROMAGNETISMIN FLAT SPACETIME We begin by briefly recalling Maxwell’s equations with We begin by introducing index notation; the language √ of relativity relies on this to condense equations. Without c = 1/ ε0µ0 = 1 and with metric convention (−, +, +, +) this notation, the Einstein field equations would instead be in flat spacetime, meaning that all line elements/metrics will 2 2 2 written as 10 separate partial differential equations, as opposed be of the form ds = −F (x)dt + dsspace, where F (x) is to a single equation. We then present the reformulation of some function of potentially all coordinates. We use c = 1, electrodynamics in terms of the Faraday tensor, an object key natural units, to drop prefactors of c, which can be restored to all relativistic electrodynamics. After a small exposition via dimensional analysis. In particular for electromagnetism in of electromagnetism in flat spacetime, we begin with general flat spacetime, we consider Minkowski space in coordinates with line element ds2 = −dt2 + dx2 + dy2 + dz2. This leads arXiv:1912.11851v3 [gr-qc] 25 Apr 2021 relativity. Generalisations of the derivative to the covariant derivative and the process of minimal coupling are introduced, to ρ bringing Maxwell’s equations in terms of the Faraday tensor to ∇ · E = , ∇ · B = 0, general relativistic form. Curvature is introduced and motivated ε0 by comparison of Schwarz’s theorem for second partial ∂B ∂E ∇ × E = − , ∇ × B = µ J + . derivatives to its covariant form, which doesn’t hold in general. ∂t 0 ∂t We present the energy-momentum tensor, an object at the heart of general relativity, coupling matter to electrodynamics. This A. Index Notation tensor, when specialised to electromagnetism, includes the By defining a so-called four-current, J µ = (ρ, J)µ, introduc- Poynting vector and the Maxwell stress tensor, among other ing the derivative notation ∂µ = (∂t, ∇)µ and the Levi-Civita physical quantities that are relevant in the Einstein equations. symbol, one reformulates Maxwell’s equations into a compact Due to properties of the electromagnetic energy-momentum form. The Levi-Civita symbol, also known as the alternating tensor is defined to be Cameron R. D. Bunney is with the School of Mathematical Sciences,  University of Nottingham (e-mail: [email protected]). 1 if ijk an even permutation of 123, Gabriele Gradoni is with the School of Mathematical Sciences, University ijk  of Nottingham and George Green Institute for Electromagnetics Research, ε = −1 if ijk an odd permutation of 123, University of Nottingham (e-mail: [email protected]). 0 otherwise. 2

ρ An even permutation is where cyclic permutations of (123) where δµ is the Kronecker delta, serving as the index equivalent are performed, whereas odd permutations involve flipping two of the identity matrix. One can intuitively think of this in elements, then followed by cyclic permutation, (132). Index, or terms of matrices. Consider matrices A and B, then (AB)ij = P suffix notation, is a way of writing equations without writing k AikBkj. The use of Einstein summation follows exactly out sums explicitly. The general rule is that if an index is in what we would expect. both the upper and lower positions, one goes ahead ans sums In relativity, keeping track of indices is important such as over it. Usually, Roman indices run from 1 to 3, and Greek whether an index is up or down. The act of moving an index run from 0 to 3. The cross product between two vectors can from down to up is called raising and vice versa is called µ be written as lowering. Consider A and Bµ, we raise and lower as follows: i ijk (a × b) = ε ajbk. µ ν µν Aν := gµν A ,B := g Bµ. This can be verified using the definition of the Levi-Civita symbol above. One may then write the differential operations The power of writing the potentials in one four-potential is in of curl and divergence using this convention. These can be an object defined as applied to rewrite out the Maxwell equations in index form µν µ ν ν µ 0 F = ∂ A − ∂ A . (3) i J i ∂iE = , ∂iB = 0, ε0 ijk i ijk i i This (2, 0)-antisymmetric tensor is called the Faraday tensor ε ∂jEk = −∂0B , ε ∂jBk = µ0J + ∂0E . and has components i Remark. Partial derivatives are denoted ∂ = ∂/∂x . µν i  0 E1 E2 E3  1 3 2 µν −E 0 B −B  νµ F =   = −F , B. Faraday Tensor −E2 −B3 0 B1  3 2 1 Since the divergence of the magnetic field is always zero, −E B −B 0 one may write this in terms of a vector potential, and Maxwell’s equations may be reduced down to two using B = ∇ × A. (1) the Faraday tensor,

νµ ν Similarly, one may write ∂µF = µ0J , (4) ∂ ∂ F = 0. (5) E = −∇φ − A. (2) [σ µν] ∂t Both fields are fully described by using only a vector and a Here, the square brackets represent the anti-symmetrisation of scalar potential. One may consider combining both potentials the tensor. On a (0, n)-tensor, this is defined as to make a four-potential to be a possible tool to describing 1 electromagnetism in a relativistic notation: A[α α ...α ] = (Aα1α2...αn + Aα2...αnα1 1 2 n n! µ µ + ... All even permutations A = (φ, A) ,Aµ = (−φ, A)µ. − Aα α ...α − Aα ...α α Minkowski space, as before, has the line element (in coordinates 2 1 n 1 n 2 (t, x, y, z)) ds2 = −dt2 + dx2 + dy2 + dz2, which can be − ... All odd permutations). compactly written as ds2 = η dxµdxν with η a diagonal µν µν The symmetrisation can be similarly defined by writing matrix, ηµν = diag(−1, 1, 1, 1)µν . One often calls ηµν the 2 A(α1α2...αn) and changing the signs on the odd permutations metric as ds and ηµν can be derived from each other. 2 µ ν to "+" signs. In particular, in the equation above, we have For general spacetimes, one writes ds = gµν dx dx as a convention. Here, gµν when considered as a matrix is a 1 ∂ F = (∂ F + ∂ F + ∂ F symmetric matrix. Objects with two indices can be written [σ µν] 3! σ µν µ νσ ν σµ as matrices and objects with one index can be written as −∂ν Fµσ − ∂µFσν − ∂σFνµ). row or column vectors. An object with indices is a tensor if it obeys a certain transformation law which can be found Remark. It is very common to use ηµν = in any introductory general relativity textbook, including [5]. diag(1, −1, −1, −1)µν instead as a metric. In this case Whether a tensor has indices up, down, or both, determines this would have its own Faraday tensor. The relationship µ the nomenclature surround it. For example an example C between the Faraday tensor used here and the other convention µ is a (1, 0)-tensor, whereas Cµ is a (0, 1)-tensor and Cν is a is (1, 1)-tensor, if these three objects obey the right transformation µν µν µν F = −F , rules. By using the metric gµν or the inverse metric g , there 1 3 is a way of moving indices up and down. The inverse metric is defined as the (2, 0)-tensor such that where the subscript 1 represents the single minus sign in our convention and the 3 represents the three minus signs in the νρ ρ gµν g = δµ, other convention. 3

III.ELECTROMAGNETISMIN CURVED SPACETIMES A. Curvature and Einstein’s Field Equations The theory of electromagnetism in Minkowski space gen- A curved spacetime is characterized by its curvature. To eralises to the case where the spacetime has curvature. The start, consider Schwarz’s theorem (partial derivatives commute) method by which the equations are coupled to gravity in the and use our substitution rule curved spacetime is (often) by replacing partial derivatives by ∂µ∂ν 7→ ∇µ∇ν , covariant derivatives ∇µ and the Minkowski metric ηµν by a general metric gµν in a process called minimal coupling. When with ∂µ∂ν = ∂ν ∂µ. If we act on some covector field Bα, it acting on vectors, the partial derivatives become covariant, can be shown that defined as α α β α ∇µ∇ν Bα − ∇ν ∇µBα = (∂ν Γµρ − ∂µΓνρ + ΓµρΓνβ ν ν ν ν ρ ∂µA 7→ ∇µA = ∂µA + ΓµρA , β α −ΓνρΓµα)Bα. and when acting on covectors, one has which in general is non zero. In fact, we define a new tensor,

ρ the Riemann curvature tensor by ∂µBν 7→ ∇µBν = ∂µBν − Γµν Bρ. α ∇µ∇ν Bρ − ∇ν ∇µBρ = R Bα, (6) ν ρνµ Here, Γµρ are the Christoffel symbols defined by which is also equivalent to ν 1 να Γµρ = g (∂µgρα + ∂ρgµα − ∂αgµρ) . ρ ρ ρ α 2 ∇µ∇ν B − ∇ν ∇µB = Rαµν B . (7) The covariant derivative arises since the standard partial This equivalence can be shown using the following symmetries derivative does not have the good transformation rules of a of the Riemann tensor, tensor when applied to anything other than functions in a β curved spacetime. However, for scalars, which have no free gβαRµνρ = Rαµνρ = −Rµανρ = −Rαµρν = Rνραµ. indices, we define the covariant derivative to be The Riemann curvature tensor is the description of the lack of commutativity of covariant derivatives. One may contract ∇µχ = ∂µχ. the Riemann curvature tensor with itself to form the Ricci Furthermore, one demands that the covariant derivative satisfies: curvature tensor, α Rραµ = Rρµ. i) Metric Compatibility: ∇µgνρ = 0, By raising an index and contracting again, one forms a scalar ii) Symmetry: Γν = Γν . µρ ρµ quantity, the (Ricci) scalar curvature: This is called the Levi-Civita connection and it is uniquely R gρµ = Rρ = R. defined by the Christoffel symbols. One may apply minimal ρµ ρ coupling first to the Faraday tensor as such This process of index raising and contraction is also known as taking the trace of the tensor — think of taking the trace of µν µρ νσ F = g g Fρσ, Pn some matrix M, the trace is simply T r(M) = i=1 Mii. = gµρgνσ (∇ A − ∇ A ) , ρ σ σ ρ Definition III.1. A spacetime is flat if and only if the compo- µρ νσ α α
= g g ∂ρAσ − ∂σAρ − Γρσ + Γσρ , nents of the Riemann curvature tensor vanish everywhere. = gµρgνσ (∂ A − ∂ A ) , ρ σ σ ρ Since the Riemann tensor is a tensor, if all components µ ν ν µ = ∂ A − ∂ A . vanish everywhere in one coordinate system, they vanish in all coordinate systems. Here we once again retrieve the flat spacetime definition of In order to state the Einstein field equations, we just need the Faraday tensor as in (3). one last object, the energy-momentum tensor. This is an object Remark. One may apply minimal coupling to the same gauge T µν which we require satisfy: symmetry on the vector potential as before 1) Symmetry: T µν = T νµ, µ µ µ 2 0µν 0 µ ν ρσ A 7→ A + ∇ χ, χ ∈ C (M) 2) Tensor Transformation: T (x ) = Jρ (x)Jσ (x)T (x), 3) Conservation Law: ∇ T µν = 0. where C2 (M) is the space of twice differentiable functions on µ our spacetime M, and due to the equivalence of partial and µ ν Where Jρ (x) and Jσ (x) are the Jacobi matrices introduced covariant derivatives on scalar functions, this is exactly the when transforming a (2, 0)-tensor [5]. The energy-momentum same transformation as before. tensor describes the matter in spacetime. We will work with By the process of minimal coupling, one also receives the the electromagnetic energy-momentum tensor defined as [9] curved spacetime generalisations of the Maxwell equations in   µν 1 µσ ν 1 µν σρ (4) and (5). T = F F σ − g FσρF . µ0 4 4

This has some interesting properties, and in matrix form this is Now, in order to use this gauge, the contraction must happen  1 2 2
µν between the vector potential and the first covariant derivative. E + B S1 S2 S3 2µ0 However, looking at the first term in µν  S1 −σ11 −σ12 −σ13 T =   . α β β α α  S2 −σ21 −σ22 −σ32 ∇β∇ A − ∇β∇ A = µ0J , (9) S3 −σ31 −σ32 −σ33 we can commute the two covariant derivatives. Recalling (7), One may recognise some familiar quantities, T 00 describes the this shows how we are able to commute covariant derivatives: 0i energy, T is the Poynting vector, describing the energy flux,   µν ασ β ασ β β ρ and σ is the Maxwell stress tensor. Also, T has the property g ∇β∇σA = g ∇σ∇βA + RρβσA , 2 2
P3 that it is trace-free since 1/(2µ0) E + B − σii = 0. α β ασ ρ i=1 = ∇ ∇βA + g RρσA , The Maxwell stress tensor is given explicitly by = ∇α∇ Aβ + R αAρ.   β ρ 1 1 2 1 2 σij = EiEj − δijE + BiBj − δijB . µ0 2 2 We define the Einstein tensor as Substituting this equation into (9) yields 1 α β α ρ β α α Gµν = Rµν − gµν R. ∇ ∇βA + Rρ A − ∇β∇ A = µ0J , 2 Finally, the Einstein field equations (without cosmological which under the Lorenz gauge reduces down to constant) are β α α ρ α − ∇β∇ A + Rρ A = µ0J . (10) Gµν = 8πG T µν . In flat spacetimes with Minkowski coordinates, the covariant The left hand side describes the geometry of spacetime and the derivatives reduce to partial derivatives with vanishing Ricci right hand side describes the matter in the spacetime. G is the tensor so this yields the flat wave equation, Newton gravitation constant. This equation can be written in a β α α different way by taking the trace of the equation (contracting ∂β∂ A = −µ0J . with gµν ). One then finds R = −8πGT , thus The left hand side of Eq. (10) has a name of its own, the de 1 Rµν = 8πG (T µν − gµν T ), Rham ’Laplacian’ [1] or the de Rham wave operator and one 2 writes this wave equation as where T = T α is the trace of the energy-momentum tensor. α  α α α ∆(dR)A = µ0J . For our energy-momentum tensor, Tα = 0 so we arrive at the Maxwell-Einstein equations, A crucial point to notice is that this wave equation has an Rµν = 8πG T µν . (8) explicit curvature term, it does not simply follow the minimal coupling method since that would require covariant derivatives Remark. If the electromagnetic field is so weak that it does to commute – the Ricci tensor term explicitly comes from this affect the curvature of spacetime significantly, it suffices to use not occurring. The equation we derived requires this term to Eqs. (4) and (5). However, if the electromagnetic field does imply conservation of charge. We now show a wave solution meaningfully affect the curvature of spacetime, one must use using the geometrical optics approximation in [18, p. 71]. If the full Einstein-Maxwell equations (8). one considers a vacuum spacetime, i.e. where J α = 0, one considers a solution to

B. Electromagnetic Wave Equation in Curved Spacetime β α α ρ − ∇β∇ A + Rρ A = 0, (11) We can rewrite (4) in terms of the four-potential (despite the Faraday tensor being the same in flat and curved spacetimes, in situations where the spacetime scale of variation of the we shall now write it explicitly with the covariant derivative electromagnetic field is much smaller than that of the curvature, for reasons which will become clear later on) i.e. where spacetime does not curve significantly to affect a αβ α β β α α wavelike solution. One would then expect to find solutions ∇βF = ∇β∇ A − ∇β∇ A = µ0J . with nearly constant amplitude, One may again employ the Lorenz gauge in its covariant form, µ Aα = Cαeikµx , (12) ∇ Aµ = 0. µ where derivatives of Cα are "small". One then substitutes this This, as in the flat case, does not fully specify the gauge Ansatz into (11) whilst neglecting the "small" derivatives of condition because one may still transform the vector potential Cα and the Ricci curvature term, leaving

0µ µ µ ! β α ∇µA = ∇µ(A + ∇ χ) = 0. − ∇β∇ A = 0. (13) So, we require the harmonic condition on the scalar function This calculation, reported in Appendix A, yields the condition χ that that for a non-trivial wave solution, the wave four-vector is µ µ ∇µ∇ χ = 0. lightlike, i.e. k kµ = 0. 5

C. Electrodynamics in Curved Spacetimes Using our convention for metric signature (−, +, +, +), and To achieve the dynamical theory, one must first familiarise writing this parametrisation explicitly oneself with a few concepts. The worldline is the trajectory of µ ν 2 2 2 gµν dx dx = ds = −c dτ , a particle. In our spacetimes, this is an assignment leads to the expression [5, p. 9] 1,3 ∼ 4 x : R → R = R , r   Z dxµ dxν t (λ) ∆τ = −gµν dλ. x1 (λ) dλ dλ λ 7→ x (λ) =   , x2 (λ) This is the time interval experienced by an observed moving x3 (λ) along this timelike worldline. The four-velocity is the nor- malised tangent vector (living in the future light cone to the where λ is a parametrisation and is not necessarily the worldline of a particle. That is, the particle’s four-velocity is (coordinate) time. Given a worldline, we can use the metric to classify these in terms of their tangent vectors. We define dxµ uµ = , tµ = dxµ(λ)/dλ to be the tangent vector to the worldline. We dλ say that and if the trajectory is parametrised by its proper time, λ = τ, µ µ ν then the four-velocity is normalised, viz., t is timelike if gµν t t < 0, µ µ ν µ ν t is lightlike if gµν t t = 0, gµν u u = −1. µ µ ν t is spacelike if gµν t t > 0. Therefore, the four-velocity is normalised if and only if the trajectory is parametrised by its proper time. Particles whose trajectories have timelike tangent vectors travel Before writing slower than the speed of light (they have mass); have lightlike down the dynamical equation for electromagnetism, we must tangent vectors travel at the speed of light (they are massless); cast ourselves back to Newtonian mechanics as a motivation. have spacelike tangent vector travel faster than the speed of Its dynamics is governed by derivatives with respect to its light. This can be demonstrated using the light cone. The light parameter. In Newtonian mechanics, time is the parameter, cone allows the decomposition of spacetime at a point in our however in relativistic mechanics we have a choice of parameter, τ spacetime into five sections: those in the future which can be in particular, we can choose proper time . Consider now by reached by a lightlike worldline; the points in the future which the chain rule can be reached by a timelike worldline; the points in the past dAµ(x(τ)) dxν (τ) := ∂ Aµ(x(τ)). which can be reached by a lightlike worldline; the points in dτ dτ ν the past which can be reached by a timelike worldline; and Now, we can apply minimal coupling and define a covariant the points which can be reached by a spacelike worldline. One parameter derivative, may typically suppress one or two spatial dimensions to ease µ ν the drawing of the light cone however the light cone is to ∇A (x(τ)) dx (τ) µ := ∇ν A (x(τ)), represent the four dimensional spacetime. dτ dτ ν dx (τ) µ µ ρ
= ∂ν A (x(τ)) + Γνρ(x(τ))A (x(τ)) , t dτ µ ν dA (x(τ)) µ dx (τ) ρ F uture A = + Γ (x(τ)) A (x(τ)). B dτ νρ dτ y

x This applies to all (1, 0)-tensors and in particular, to the four- C velocity uµ, P ast ∇uµ(τ) duµ(τ) = + Γµ (x(τ))uν (τ)uρ(τ). dτ dτ νρ Fig. 1. Representation of the light cone (suppressing one spatial dimension). In general, the dynamical equation (the generalisation of F = ma) is ∇uµ(τ) The boundary of the light cone is formed by the lightlike m = F µ(x(τ)), (14) worldlines from the centre point so anything on the light cone’s dτ boundary is on a lightlike worldline such as point B. Anything where F µ is a so-called four-force. Using the concept of a inside the light cone (and not on the boundary) is represented four-velocity parametrised by proper time, one can now write by a timelike worldline, point A, and anything outside the down the dynamical equation for electromagnetism: C light cone is on a spacelike worldline, point ! In Newtonian µ ∇u (τ) µν mechanics, the parametrisation is the coordinate time, however m = qF uν . in relativistic mechanics because this coordinate time is not dτ invariant under transformations, one must use a so called proper The right hand side of the equation describes the Lorentz force, time. This is the time experienced by a timelike particle in its m is the mass, and q is the charge of the particle, F µν is again µ frame of reference and this parametrisation is often called τ. the Faraday tensor, and uν = gµν u is the four-velocity. 6

D. Geodesics can be approximated by the electric dipole approximation in cylindrical coordinates, which reads An important example of (14) is when F µ = 0. In Newtonian µ p ω mechanics, we know that ma = 0 yields solutions of the form φ = − 0 0 cos θ sin (ω (t − r)), x = a + bt, straight lines. These are the shortest distances 4πr µ0p0ω between two points in Euclidean geometry; the name for the A = − sin (ω (t − r))zˆ, 4πr shortest distance is a geodesic. Now, when F µ = 0, we recover the geodesic equation which is usually written in terms of the or, in compact tensorial notation worldlines cos θ 2 µ ν ρ µ µ0p0ω  0  d x (τ) µ dx (τ) dx (τ) A = − sin (ω (t − r))   . + Γ = 0. (15) 4πr 0 dτ 2 νρ dτ dτ   1 This plays an important role in differential geometry as it tells One should note that in the argument of the trigonometric you the shortest path between two points. An example would functions, t − r = t − |x − x0| = t , the retarded time. Then be if one used the metric for a sphere of radius 1, r using zˆ = cos θ rˆ − sin θ θˆ, one can calculate the fields that ds2 = dθ2 + sin2 θdϕ2, these potentials create, which reads µ p ω2 E = − 0 0 sin θ cos (ω (t − r))θˆ, solving this equation would yield solutions of great circles — 4πr circles that pass through the centre of the sphere, diving it µ p ω2 B = − 0 0 sin θ cos (ω (t − r))ϕˆ. into two equal parts. Solving this equation is in general not 4πr easy. The geodesic equation is more familiar to the reader than The Faraday tensor here is at first glance, it is actually the Euler-Lagrange equation for two different action integrals. Classically, the dynamics of a  0 0 1 0 2 system described by a Lagrangian L[q(t), q˙(t), t] can be found µν µ0p0ω  0 0 1 0 F = − sin θ cos (ω (t − r))   . by looking at the minimum of its action. The action S[q(t)] is 4πr −1 −1 0 0 defined as an integral 0 0 0 0 Z These fields correspond to waves of frequency ω, travelling S[q(t)] = L[q(t), q˙(t), t]dt. radially at the speed of light. The case of dipole radiation is −iωt considered in [10, p. 271] with q(t) = q0e , leading to the In Newtonian mechanics, the parameter used is time. In complex fields relativity, time is on an equal footing to space so we cannot µ p ω2 µ p ω2 single it out and integrate over it, therefore we integrate over E = − 0 0 sin θ eiω(r−t)θˆ, B = − 0 0 sin θ eiω(r−t)ϕˆ. our parameter λ. The two actions relevant to the geodesic 4πr 4πr equation are The strength of this notation is that it allows one to recognise that both E and B are spherical waves of the form u(r, t) = r Z dxµ dxν 1 Z dxµ dxν 1/r F (r − t) + 1/r G(r + t). S = −g dλ, and S = g dλ. 1 µν dλ dλ 2 2 µν dλ dλ √ The first has a geometric meaning since the integrand, ds2 = ds, is the line element. Therefore, by computing the Euler- Lagrange equations minimizes the path. The second is however computationally easier to work with.

E. Flat Dipole Radiation For fullness of this review, we present the flat spacetime tensorial form of dipole radiation. Dipole approximations are Fig. 2. Spherical Wave with θ = π/2, from above and from the side. known in curved spacetimes, however these are more involved and for compactness, we do not cover this. For the case of As shown in Fig. 2, as r → ∞ increases, the waves become Schwarzschild, see [11] which utilises the Newman-Penrose weaker in amplitude and tend to plane waves as r → ∞ formalism. Consider the case of electric dipole radiation by imagining two small metal spheres, connected by a thin wire, separated F. Propagation of Radiation in Curved Spacetimes by a distance d, with a charge oscillating between the two We have seen that an oscillating dipole emanates spherical balls, q(t) = q0 cos(ωt), where ω is the angular frequency electromagnetic waves and it is known that electromagnetic and the maximum value of the dipole moment is p0 = q0d. waves are propagated by the photon which travels on lightlike Following [7, p. 469-471], the scalar and vector potentials worldlines. Another term for lightlike is null. Light and, more 7

generally, photons travel on null geodesics. Geodesics were pt is the momentum conjugate to the time t in the context of described in section III-D; the null part refers to the lightlike Lagrangian mechanics. Taking the derivative of L with respect 2 µ ν condition. Consider the general metric ds = gµν dx dx . to t˙, we have 2 This would be null when ds = 0. Therefore, one considers pt = −F (r)t˙ = −E. the ray paths under this condition. A useful geometry is the Schwarzschild geometry, described by the metric Using the notation, dt = E/F (r) dλ, then substituting into our metric, we have   2 2 2GM 2 dr 2 2 ds = − 1 − dt + + r dΩ , E2 dr2 r 1 − 2GM 0 = −F (r) dλ2 + + r2dΩ2, (17) r F 2(r) F (r) which describes how a spherically symmetric object (such as E2 1   a perfectly spherically symmetric star) curves the spacetime 2 2 ˙ 2 =⇒ 0 = −F (r) 2 + r˙ + r θ + sin θ ϕ˙ , about it. M is a mass parameter, related to the Newtonian F (r) F (r) (18) mass by M = GMNewt, where G is the gravitational constant 1 and MNewt is the Newtonian measured mass. One notes that = −E2 +r ˙2
. (19) the case where M = 0 gives the Minkowski metric of special F (r) relativity so this encompasses the case of no mass. The spherical This yields our second condition on our radial null rays,r˙ = 2 symmetry is described by dΩ2 = dθ2 + sin θdϕ2, the metric ±E. This clearly has the solution r = ±Eλ + const. One can for the sphere1. Though, this is not the most general form of solve for t in terms of r by realising that this metric t˙ dt dλ dt 1 dr2 = = = ∓ . ds2 = −F (r)dt2 + + r2dΩ2, (16) r˙ dλ dr dr F (r) F (r) which encompasses metrics such as the Reissner-Nordström Then one solves for this ODE to see  r  metric describing a non-rotating charged mass and also the anti ±t = r + 2GM ln − 1 + const. de-Sitter spacetime metric, a spacetime with constant negative 2GM curvature - its Ricci scalar is constant and negative. Such a These were derived from the geodesic Lagrangian using the spacetime is of interest for string theorists thanks to its role in Euler-Lagrange equations so we can clearly see that these are the AdS/CFT correspondence [13]. geodesics but it is not always so easy to see that these are still When studying the geodesics of null rays, one considers null (lightlike). both radial, and non-radial null rays. The radial case is when 2) Non-Radial Null Rays: By the spherical symmetry of there is no motion within the angles, and the non-radial case the geometry, one choose coordinates such that, for example, is when one does allow movement within the angles. θ = π/2 since a geodesic is described in an invariant plane [16, 1) Radial Null Rays: One has ds2 = 0 for the null condition p. 11]. Thus dθ = 0, θ˙ = 0, and sin θ = 1. This simplifies and the radial condition comes into play in the Euler-Lagrange the metric greatly, equations for the geodesic Lagrangian. From section III-D, dr2 one has ds2 = −F (r)dt2 + + r2dϕ2. Z 1  1  F (r) S = −F (r)t˙2 + r˙2 + r2θ˙2 + r2 sin2 θϕ˙ 2 dλ. 2 F (r) In similar fashion to before, the geodesic Lagrangian in this case is Since these rays are purely radial one has θ˙ =ϕ ˙ = 0 where 1  1  a dot refers to differentiation with respect to λ. The geodesic L = −F (r)t˙2 + r˙2 + r2ϕ˙ 2 . Lagrangian is then 2 F (r) 1  1  Upon inspection, one sees that the Lagrangian is independent L = −F (r)t˙2 + r˙2 . 2 F (r) of both t and ϕ, leading to one of the same conservation laws as before F (r) t˙ = E, as well as To refresh the reader, the Euler-Lagrange equations are d ∂L ∂L ∂L 2 = . pϕ = = r ϕ˙ = `, (20) dλ ∂x˙ µ ∂xµ ∂ϕ˙ One can see that the geodesic Lagrangian is independent of t, which describes conservation of angular momentum. Returning therefore ∂L/∂t = const, and we call this constant −E. More to the metric, explicitly, by setting µ = 0, we have, 1 0 = ds2 = −F (r)dt2 + dr2 + r2dϕ2, d ∂L d ∂L ∂L F (r) = pt = = = 0. dλ ∂x˙ 0 dλ ∂x0 ∂t ˙2 1 2 2 2 | {z } | {z } =⇒ 0 = −F (r)t + r˙ + r ϕ˙ , LHS of Euler-Lagrange Equations RHS of Euler-Lagrange Equations F (r) E2 r˙2 `2 = − + + , 1This is the definition of dΩ2 in this work, however it is important to F (r) F (r) r2 know that this is also used for spherical symmetry in higher dimensions. More `2 explicitly, the dimension of our spacetime is 4: we use 1 for time, 1 for radial 2 2 =⇒ r˙ = E − 2 F (r), distance and 2 for angles. One also considers a spacetime of dimension n, then r dΩ2 would then be the metric of Sn−2, the sphere of dimension (n − 2). 8 leaving the radial equation ODE, our gravitational redshift formula. s In our calculations, we use only spacetimes whose metrics `2  2GM  r˙ = ± E2 − 1 − . are Schwarzschild-like, i.e. of the form (16), describing well r2 r the exterior regions of spherical objects such as planets, stars, This same equation can be found [17, p. 1] with differing or non-rotating black holes. However, not all spacetime metrics constants, we have α = `, and β = E/`. However, if one is are of this form. A ubiquitous cosmological model of the simply interested in the orbit that the ray takes, one looks at universe which encompasses its expansion is the Friedmann- dϕ/dr =ϕ/ ˙ r˙. Lemaître-Robertson-Walker (FLRW) model, with metric ds2 = −dt2 + a2(t) dx2 + dy2 + dz2
, G. Gravitational Redshift where a(t) is referred to as the scale factor of the universe. The effect of a mass curving spacetime causes a shift in The specific details of the universe it describes are of a the frequency of light/electromagnetic radiation emitted by homogeneous, isotropic, expanding universe. This model is one particle observed by another and this has applications explored in [2], where the resulting electric field to the within communication, such as two satellites communicating. sourced Maxwell equations is found. By introducing the We present here a geometric derivation based on Schrödinger’s t variable t˜ = R 1/a(t0)dt0 and noting dt˜ = 1/a(t)dt via the derivation [15, p. 47-53]. t0 fundamental theorem of calculus, the metric is transformed into 2 2 2 2 2 2 2 2 γ0 ˜ q0 ds = a (t)(−dt + dx + dy + dz ) = a (t)ds . p0 Minkowski f µ Since the Minkowski metric and FLRW metric differ by only q p γ an overall positive scale factor, they are said to be conformally equivalent. The redshift of monochromatic light in FLRW µ µ u (τ) v (˜τ) spacetime is calculated in [2] for the case a(t) = a0 + bt as Source Observer ν 0 observed = e−b|x−x |, Fig. 3. Two observers and their respective worldlines, interacting via emission νemitted of photons. where x and x0 are the positions of the observer and source respectively. Consider Fig. 3 above. This shows two particles travelling 1) Stationary Observers: Consider the case where both the along timelike worldlines, parametrised by their own proper source and observer are separated only radially and neither times, τ and τ˜, each having their own four-velocities, uµ(τ) are moving spatially. The light being emitted from the source and vµ(˜τ). Naming the left hand side particle the source and arrives at the observer and travels along a radial, null geodesic the right hand side the observer, consider when the source is so that we have explicit formulae for the elements of the tangent at p along its worldline and it emits a photon that the observer vector of this geodesic. Choosing an energy scale such that observes at point q. Then, (very) shortly after, the source emits E = 1, one has a second photon at point p0, being observed at point q0. The  1  world lines for each photon are wµ = (t,˙ r,˙ 0, 0) = , ±1, 0, 0 .23 F γ = xµ(λ), The four-velocities of the source and observer also take simple γ0 = yµ(λ) = xµ(λ) + f µ(λ), forms, simply being uµ = (u0, 0, 0, 0) and vµ = (v0, 0, 0, 0). where both γ and γ0 are null geodesics and f µ is "small". We have one final piece of information to use — since the four- 0 µ µ From p to p , the source moves a distance of dτ(u )p = (f )p velocities are parametrised by proper time, each is normalised, µ µ and the observer moves a distance of dτ˜(v )q = (f )q where thus one has that dτ γ γ0 is the proper time elapsed for the source between and µ 0 2 0 1 −1 = u uµ = −F (u ) =⇒ u = √ . and dτ˜ is similarly defined. It is shown in Appendix B that F the quantity f x˙ µ is constant over γ. As a corollary of this, µ The positive sign is chosen since the four-velocity points and writing x˙ = dx/dτ, we have towards future regarding time; in terms of the lightcone, the µ µ (x ˙ fµ)p = (x ˙ fµ)q, four-velocity lives in the future null cone rather than the past. µ µ (x ˙ uµ)pdτ = (x ˙ vµ)qdτ,˜ We can now use the redshift formula (21) to calculate what the µ redshift due to gravity is. Using the notation Fi = 1−2GM/ri dτ (x ˙ vµ)q =⇒ = ν , where the source is located at r1 and the observer at r2, dτ˜ (x ˙ uν )p µ 1 √1 r µ ν (w vµ) w0v (−F2 ) F where x˙ is the tangent of the null geodesic that the photon emitted = q = 0 = F2 F2 = 1 . µ σ 0 1 1 takes and one might call this instead w . Calling ν the νobserved (w uσ) w u0 (−F1 √ ) F2 p F1 F1 frequency of the observed light, we now have 2 µ In r˙, we have two possibilities for the cases whether the light’s r-coordinate ν (w vµ) emitted = q , (21) is increasing or decreasing. σ 3Note also since ϕ˙ = 0, that ` = 0, rendering r˙ = ±1. νobserved (w uσ)p 9

Thus the final formula is given by 3) Radial and Azimuthal Motion: We can also go one step further with our calculations and allow azimuthal motion of v u 2GM our source and observer. We will use units such that c = 1, νemitted u1 − r 1 E = 1, and ` = 1.4 To do this, we require that our light travels = t 2GM , (22) νobserved 1 − r2 on a null, non-radial geodesic given by ! 1 r 1 1 showing the gravitational redshift. Then by rearranging and wµ = (t,˙ r,˙ 0, ϕ˙) = , ± 1 − F, 0, , observing at infinity, one has F r2 r2

1 and that our four-velocities are of the forms µ 0 r ϕ µ 0 r ϕ νobserved = q νemitted, (23) u = (u , u , 0, u ) and v = (v , v , 0, v ). We write out 1 − 2GM 0 1 r 2 ϕ r1 explicitly that uµ = (−F1u , u , 0, r u ) to be clear on F1 1  GM  the lowered index form of the four-velocity. Again, the four- ≈ 1 + νemitted. (24) velocity must be normalised such that r1 µ 0 2 1 r 2 2 ϕ 2 u uµ = −F1(u ) + (u ) + r1(u ) = −1, This approximation is often found in textbooks to describe F1 gravitational redshift. This formula can be found in, for example giving [3, p. 620] by noting νobserved = v∞ and νemitted = v0. q 0 1 r 2 2 ϕ 2 u = F1 + (u ) + r1F1(u ) . 2) Radial Motion: Now we have the groundwork set, one F1 can introduce the case where both the source and observer are The redshift formula gives us moving radially. The change to the setup is in the four-velocity µ 0 r s of the source and observer. We now have u = (u , u , 0, 0) µ 0 1 r 1 ϕ µ 0 r µ µ w uµ = −u ± u 1 − F1 + u , and v = (v , v , 0, 0) such that u uµ = v vµ = −1, i.e. that 2 F1 r1 u and v are timelike. This translates to and 1 s 0 2 r 2 µ 0 1 r 1 ϕ −F1(u ) + (u ) = −1. w v = −v ± v 1 − F + v . F µ 2 2 1 F2 r2 where we want to solve for the zeroth component. Solving and Taking the quotient and substituting in for u0, v0, one sees simplifying, we therefore have ν F observed = 2 × νemitted F1 0 1 p r 2 u = F1 + (u ) , p r 2 2 ϕ 2 rq 1 ϕ  F1 F1 + (u ) + r1F1(u ) ∓ u 1 − 2 F1 − F1u r1  q  , 0 p r 2 2 ϕ 2 r 1 ϕ v F2 + (v ) + r2F2(v ) ∓ v 1 − 2 F2 − F2v with the same calculation being made for . The redshift r2 formula now tells us that (26)

0 1 1 0 1 r giving the equation for the redshift due to gravity and both ν w u + w u (−F1u ) ± u observed = 0 1 = F1 F1 , 0 1 1 0 1 r radial and azimuthal motion. νemitted w v0 + w v1 (−F2v ) ± v F2 F2 Taking the non-azimuthal limit requires care as the units  0 r  p r 2 r ! F2 F1u ∓ u F2 F1 + (u ) ∓ u with ` = 1 can no longer hold. Therefore, to do this, one must = = . ϕ 0 2 p r 2 r first place ` back into the equation and then substitute u = 0, F1 F2v ∓ v F1 F2 + (v ) ∓ v vϕ = 0, ` = 0, which recovers Eqn. (25). Filling out this equation, one retrieves Together, Eqs. (22), (25), and (26) give the observed redshifts for different set ups between observer and source, q and are relevant wherever the Schwarzschild geometry, or a 2GM  2GM r 2 r  1 − 1 − r + (u ) ∓ u νobserved r2 1 Schwarzschild-like geometry, can describe a spacetime such = 2GM  q  , (25) νemitted 1 − 2GM r 2 r as outside the Schwarzschild radius of a Schwarzschild black r1 1 − + (v ) ∓ v r2 hole, or outside the surface of a spherically symmetric star. The equations show the combined effects of both gravitational for the redshift due to both radial motion and a gravitational redshift and the Doppler shift due to the motion of both source field described by the Schwarzschild metric. Again, under and observer. the stationary assumption with r at infinity, we recover (23). 2 We can see the well-known gravitational redshift by plotting Taking the Newtonian limit of this formula, i.e. such that (22). We assume that r > r > 2GM and that the observer gravitational effects are negligible F ≈ 1, we use r , r >> 2 1 i 2 1 at r is not moving. Gravitational effects will be strongest 2GM, ur, vr << 1, and keep to leading order in ur and vr, 2 close to the mass. With appropriate units, we can normalise we find ν /ν = (1∓ur)/(1∓vr). This is exactly the observed emitted 2GM = 1. We then look at the region 1 < r < 3/2, that is classical Doppler effect. This calculation tells us that the ∓-sign 1 terms in (25) give the leading Doppler contributions whereas 4We can do this as long as ϕ˙ > 0, since we can form dimensions of time, the Fi terms tell us about the gravitational contributions. mass, and length using E, `, and c. 10

r2 = 3/2. Since in most physical situations, both source and Simplifying by using the Kronecker deltas to change indices observer will be far away from the Schwarzschild radius, the and the metrics to raise or lower indices, we find effect of the gravitational redshift is negligible, however it is  ρ  ρ α βµ ikρx σ α βµ σ ikρx still of interest for purposes of reviewing phenomena arising −iC kµδ ikσe ∂βx = C δ δβ kµkσ e , ρ from the general relativistic formalism. α ikρx βν µ σ = C e g δν δβ kµkσ, α βν σ µ = A g δν kσ δν kµ, α βν = A g kν kβ, α ν = A k kν .

APPENDIX B µ PROOF THAT fµx˙ IS CONSTANT OVER γ 0 µ ν Given that γ is null, we have 0 = gµν (y)y ˙ y˙ . By Taylor expansion to first order in f,     0 = (g (x) + f σ∂ g (x)) x˙ µ + f˙µ x˙ ν + f˙ν , Fig. 4. Gravitational Redshift. µν σ µν µ ν µ ν σ µ ˙ν = gµν x˙ x˙ +x ˙ x˙ (∂σgµν ) f + 2gµν x˙ f . Figure (4) asymptotes at r = 1, physically representing the µ ν 1 We also have gµν x˙ x˙ = 0 because γ is null, and also that Schwarzschild radius. This radius represents the point at which µ ˙ν d µ ν σ d µ all future-directed non-spacelike curves must go inwards to the 2gµν x˙ f = 2 (gµν x˙ f ) − 2f (gµσx˙ ) . mass at the centre. At r1 = r2 = 3/2, we find νobs/νemit = 1 dλ dλ as expected. Thus returning to the above equality and continuing, we have d  d 1  0 = 2 (f x˙ µ) − 2f σ (g x˙ µ) − (∂ g )x ˙ µx˙ ν . dλ µ dλ σµ 2 σ µν IV. CONCLUSION But the term within the bracket is zero since γ is a geodesic 1 µ ν We reviewed the main relevant concepts of flat spacetime as this is the Euler-Lagrange equation for L = 2 gµν x˙ x˙ . electromagnetism in order to introduce the general relativistic d 1 d (g x˙ µ) − (∂ g )x ˙ µx˙ ν = 0 =⇒ 0 = 2 (f x˙ µ) . electromagnetism. We apply this formalism to derive two dλ σµ 2 σ µν dλ µ equations determining the redshift of light detected by a moving Which is exactly the condition telling us that f x˙ µ is constant. observer, as emitted by a moving observer. These equations µ couple together the effects of the Doppler shift and gravitational REFERENCES redshift to give a single equation, depending on the geometry of the two moving observers. Further research in this area would [1] Donato Bini, Christian Cherubini, Robert Jantzen, and Remo Ruffini. De rham wave equation for tensor valued p-forms. In- be in the direction of how interference in curved spacetime ternational Journal of Modern Physics D.8., 12:1363–1384, 2003. might affect communication between two antennas. DOI:10.1142/S0218271803003785. [2] Neil V. Budko. Exponential cosmological redshift in a linearly expanding universe. 2009. [3] Bradley W. Carroll and Dale A. Ostlie. An Introduction to Modern APPENDIX A Astrophysics. Pearson Education Inc., 2013. GEOMETRIC OPTIC WAVE EQUATION SOLUTION [4] Sean Carroll. Lecture notes on general relativity. https://ned.ipac. caltech.edu/level5/March01/Carroll3/Carroll_contents.html, December 1997. Accessed September 17th 2019. Substituting Eq. (12) into (13), we find [5] Sean Carroll. Spacetime and Geometry: An Introduction to General Relativity. Pearson, first edition, 2013. β α [6] Steven Errede. A brief history of the development of classical 0 = −∇β∇ A , electrodynamics. http://web.hep.uiuc.edu/home/serrede/P435/Lecture_  µ  β α ikµx Notes/A_Brief_History_of_Electromagnetism.pdf. Accessed November = −∇β∇ C e . 8th 2019. [7] David J. Griffiths. Introduction to Electrodynamics. Pearson Education Inc., fourth edition, 2013. The covariant derivative on a scalar is just a partial derivative [8] F. Gronwald and J. Nitsch. The structure of the electromagnetic field as ν ν ν µ and note that ∂µx = ∂x /∂x = δν . Applying this twice derived from first principles. IEEE Antennas and Propagation Magazine, yields 43:64–79, 2001. 10.1109/MAP.2001.951560. [9] S.W. Hawking and G.F.R. Ellis. The large scale structure of space-time. Cambridge University Press, 1973.  σ  α ikσ x β µ = −C ∇β ikµe ∂ x , [10] John David Jackson. Classical Electrodynamics. John Wiley & Sons, Inc., third edition, 1962. σ α βµ ikσ x [11] B Linet. Dipole field solution of maxwell's equations in the schwarzschild = −iC kµδ ∇βe , metric. Journal of Physics A: Mathematical and General, 8(3):328–333,  ρ  α βµ ikρx σ mar 1975. = −iC kµδ ikσe ∂βx . [12] Jorma Louko. Black holes, 2019. 11

[13] Juan Maldacena. International Journal of Theoretical Physics, 38(4):1113–1133, 1999. [14] Charles W. Misner, Kip S. Thorne, and John Archibald Wheeler. Gravitation. Princeton University Press, 2017. [15] Erwin Schrödinger. Expanding Universes. Cambridge University Press, Cambridge, second edition, 1956. [16] G. Slezakova. Geodesic Geometry of Black Holes. PhD dissertation, The University of Waikato, Hamilton, New Zealand, 2006. Accessed September 17th 2019. [17] J. K. Synge. The escape of photons from gravitationally intense stars. Monthly Notices of the Royal Astronomical Society, 131:463, 1965. [18] Robert M. Wald. General Relativity. The University of Chicago Press, 1984. [19] Steven Weinberg. Gravitation and Cosmology: Principles and Appli- cations of the General Theory of Relativity. John Wiley & Sons, Inc., 1972.