WASHINGTON UNIVERSITY Department of Physics ...

3 downloads 0 Views 805KB Size Report
Thanks to Clifford Will for his help and support. Thanks also to my mother Betsey ...... [38] C. R. Johnson, Phys. Rev. D 31 (1985) 1252. [39] V. V. Narlikar, and ...
WASHINGTON UNIVERSITY Department of Physics

Dissertation Examination Committee: Clifford M. Will, Chair Mark Alford Ramanath Cowsick Renato Feres Barry Spielman Wai-Mo Suen

¨ EXTENSIONS OF THE EINSTEIN-SCHRODINGER NON-SYMMETRIC THEORY OF GRAVITY by James A. Shifflett

A dissertation presented to the Graduate School of Arts and Sciences of Washington University in partial fulfillment of the requirements for the degree of Doctor of Philosophy May 2008 Saint Louis, Missouri

Acknowledgements Thanks to Clifford Will for his help and support. Thanks also to my mother Betsey Shifflett for her encouragement during my graduate studies, and to my late father John Shifflett for his encouragement long ago. This work was funded in part by the National Science Foundation under grants PHY 03-53180 and PHY 06-52448.

ii

Contents Acknowledgements

ii

Abstract

v

1 Introduction

1

2 Extension of the Einstein-Schr¨ odinger theory 2.1 The Lagrangian density . . . . . . . . . . . . 2.2 The Einstein equations . . . . . . . . . . . . . 2.3 Maxwell’s equations . . . . . . . . . . . . . . . 2.4 The connection equations . . . . . . . . . . .

for Abelian fields . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

. . . .

7 7 12 16 17

3 Exact Solutions 3.1 An exact electric monopole solution . . . . . . . . . . . . . . . . . . . 3.2 An exact electromagnetic plane-wave solution . . . . . . . . . . . . .

28 28 32

4 The 4.1 4.2 4.3

equations of motion The Lorentz force equation . . . . . . . . . . . . . . . . . . . . . . . . Equations of motion of the electric monopole solution . . . . . . . . . The Einstein-Infeld-Hoffmann equations of motion . . . . . . . . . . .

34 34 36 41

5 Observational consequences 5.1 Pericenter advance . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5.2 Deflection and time delay of light . . . . . . . . . . . . . . . . . . . . 5.3 Shift in Hydrogen atom energy levels . . . . . . . . . . . . . . . . . .

51 51 60 69

6 Application of Newman-Penrose methods 6.1 Newman-Penrose methods applied to the exact field equations . . . . 6.2 Newman-Penrose asymptotically flat O(1/r2 ) expansion of the field equations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

71 71

7 Extension of the Einstein-Schr¨ odinger theory for 7.1 The Lagrangian density . . . . . . . . . . . . . . 7.2 Invariance properties of the Lagrangian density . 7.3 The field equations . . . . . . . . . . . . . . . . . 8 Conclusions

non-Abelian . . . . . . . . . . . . . . . . . . . . . . . .

90

fields107 . . . 107 . . . 112 . . . 117 122

iii

A A divergence identity

123

˜σ = 0 B Variational derivatives for fields with the symmetry Γ [µσ]

124

C Approximate solution for Nνµ in terms of gνµ and fνµ

126

˜ α in terms of gνµ and fνµ D Approximate solution for Γ νµ

129

E Derivation of the generalized contracted Bianchi identity

134

F Validation of the EIH method to post-Coulombian order

139

G Application of point-particle post-Newtonian methods

143

H Alternative derivation of the Lorentz force equation

147

I

Alternative derivation of the O(Λ−1 b ) field equations

151

J A weak field Lagrangian density

153

K Proca-waves as Pauli-Villars ghosts?

161

L Lm , Tµν , j µ and kinetic equations for spin-0 and spin-1/2 sources

165

M Alternative ways to derive the Einstein-Schr¨ odinger theory

172

N Derivation of the electric monopole solution

179

O The function Vˆ (r) in the electric monopole solution

185

P The electric monopole solution in alternative coordinates

190

Q The electromagnetic plane-wave solution in alternative coordinates192 R Some properties of the non-symmetric Ricci tensor

197

S Calculation of the non-symmetric Ricci tensor in tetrad form

201

T Proof of a nonsymmetric matrix decomposition theorem

204

U Calculation of the exact Υbca in Newman-Penrose form

206

V Check of the approximate Υbca in Newman-Penrose form

210

W Kursunoglu’s theory with sources and non-Abelian fields

214

X Possible extension of the theory to non-Abelian symmetric fields

220

Bibliography

247 iv

Abstract We modify the Einstein-Schr¨odinger theory to include a cosmological constant Λz which multiplies the symmetric metric. The cosmological constant Λz is assumed to be nearly cancelled by Schr¨odinger’s cosmological constant Λb which multiplies the nonsymmetric fundamental tensor, such that the total Λ = Λz +Λb matches measurement. The resulting theory becomes exactly Einstein-Maxwell theory in the limit as |Λz | → ∞. For |Λz | ∼ 1/(Planck length)2 the field equations match the ordinary Einstein and Maxwell equations except for extra terms which are < 10−16 of the usual terms for worst-case field strengths and rates-of-change accessible to measurement. Additional fields can be included in the Lagrangian, and these fields may couple to the symmetric metric and the electromagnetic vector potential, just as in Einstein-Maxwell theory. The ordinary Lorentz force equation is obtained by taking the divergence of the Einstein equations when sources are included. The Einstein-Infeld-Hoffmann (EIH) equations of motion match the equations of motion for Einstein-Maxwell theory to Newtonian/Coulombian order, which proves the existence of a Lorentz force without requiring sources. An exact charged solution matches the Reissner-Nordstr¨om solution except for additional terms which are ∼ 10−66 of the usual terms for worst-case radii accessible to measurement. An exact electromagnetic v

plane-wave solution is identical to its counterpart in Einstein-Maxwell theory. Pericenter advance, deflection of light and time delay of light have a fractional difference of < 10−56 compared to Einstein-Maxwell theory for worst-case parameters. When a spin-1/2 field is included in the Lagrangian, the theory gives the ordinary Dirac equation, and the charged solution results in fractional shifts of < 10−50 in Hydrogen atom energy levels. Newman-Penrose methods are used to derive an exact solution of the connection equations, and to show that the charged solution is Petrov type-D like the Reissner-Nordstr¨om solution. The Newman-Penrose asymptotically flat O(1/r2 ) expansion of the field equations is shown to match Einstein-Maxwell theory. Finally we generalize the theory to non-Abelian fields, and show that a special case of the resulting theory closely approximates Einstein-Weinberg-Salam theory.

vi

Chapter 1 Introduction Einstein-Maxwell theory is the standard theory which couples general relativity with electrodynamics. In this theory, space-time geometry and gravity are described by a metric gµν which is symmetric (gµν = gνµ ), and the electromagnetic field Fµν is antisymmetric (Fµν = −Fνµ ). The fact that these two fields could be combined together into one second rank tensor was noticed long ago by researchers looking for a more unified description of the physical laws. The Einstein-Schr¨odinger theory is a generalization of vacuum general relativity which allows a nonsymmetric field Nµν in place of the symmetric gµν . The theory without a cosmological constant was first proposed by Einstein and Straus[1, 2, 3, 4, 5]. Schr¨odinger later showed that it could be derived from a very simple Lagrangian density if a cosmological constant Λb was included[6, 7, 8]. Einstein and Schr¨odinger suspected that the theory might include electrodynamics, where the nonsymmetric “fundamental tensor” Nµν contained both the metric and electromagnetic field. However, this hope was dashed when it was found that the theory did not predict a Lorentz force between charged particles[9, 10].

1

In this dissertation we describe a simple modification of the Einstein-Schr¨odinger theory[11, 12, 13, 14] which closely approximates Einstein-Maxwell theory, and where the Lorentz force does occur. The modification involves the addition of a second cosmological term Λz gµν to the field equations, where gµν is the symmetric metric. We assume this term is nearly canceled by Schr¨odinger’s “bare” cosmological term Λb Nµν , where Nµν is the nonsymmetric fundamental tensor. The total “physical” cosmological constant Λ = Λb + Λz can then be made to match cosmological measurements of the accelerating universe. The origin of our Λz is unknown. One possibility is that Λz could arise from vacuum fluctuations, an idea discussed by many authors[15, 16, 17, 18]. Zero-point fluctuations are essential to both quantum electrodynamics and the Standard Model, and are thought to be the cause of the Casimir force[16] and other effects. With this interpretation, the fine tuning of cosmological constants is not so objectionable because it resembles mass/charge/field-strength renormalization in quantum electrodynamics. For example, to cancel electron self-energy in quantum electrodynamics, the “bare” electron mass becomes large for a cutoff frequency ωc ∼ 1/(Planck length), and infinite if ωc → ∞, but the total “physical” mass remains small. In a similar manner, to cancel zero-point energy in our theory, the “bare” cosmological constant Λb ∼ ωc4 × (Planck length)2 becomes large if ωc ∼ 1/(Planck length), and infinite if ωc → ∞, but the total “physical” Λ remains small. There are other possible origins of Λz . For example Λz could arise dynamically, related to the minimum of a potential of some additional field in the theory. Apart from the discussion above, speculation about the origin of Λz is outside the scope of this dissertation. Our main 2

goal is to show that the theory closely approximates Einstein-Maxwell theory, and for non-Abelian fields the Einstein-Weinberg-Salam theory (general relativity coupled to electro-weak theory). Like Einstein-Maxwell theory, our theory can be coupled to additional fields using a symmetric metric gµν and vector potential Aµ , and it is invariant under a U(1) gauge transformation. The theory does not enlarge the invariance group. When coupled to the Standard Model, the combined Lagrangian is invariant under the usual U(1) ⊗ SU(2) ⊗ SU(3) gauge group. The usual U(1) gauge term F µνFµν is incorporated together with the geometry, and is not explicitly in the Lagrangian. The non-Abelian version of the theory can also be coupled to the Standard Model, in which case both the U(1) and SU(2) gauge terms are incorporated together with the geometry. This is done much as it is done in [19, 20] with Bonnor’s theory. Whether the SU(3) gauge term of the Standard model could also be incorporated with a larger gauge group, or by using higher space-time dimensions, is beyond the scope of this dissertation. The Abelian version of our theory is similar to [21, 22] but with the opposite sign of Λb and Λz . Because of this difference our theory involves Hermitian fields instead of real fields, and the spherically symmetric solutions have much different properties near the origin and do not come in an infinite set. The Abelian version of our theory is also roughly the electromagnetic dual of another theory[23, 24, 25, 26]. Compared to all of these other theories, our theory also allows coupling to additional fields (sources), and it allows Λ 6= 0, and it is derived from a Lagrangian density which incorporates a new type of non-symmetric Ricci tensor with different invariance properties. Many other modifications of the Einstein-Schr¨odinger theory have been consid3

ered. For example in Bonnor’s theory[27, 28] the antisymmetric part of the fundamental tensor N[τ ρ] or its dual is taken to be the electromagnetic field, and a Lorentz √ force is derived, but only because a −N N a[ρτ ] N[τ ρ] term is appended onto the usual Lagrangian density. Other theories include an assortment of additional terms in the Lagrangian density[29, 30]. Such theories lack the mathematical simplicity of the original Einstein-Schr¨odinger theory, and for that reason they seem unsatisfying. This criticism seems less applicable to our theory because there are such good motivations √ for including a Λz −g term in the Lagrangian density. Some previous work[31, 32, 33] shows that the original Einstein-Schr¨odinger theory has problems with negative energy “ghosts”. As will be seen in §2.4, this problem is avoided in our theory in an unusual way. In [31, 32, 33] referenced above, the electromagnetic field is assumed to be an independent field added onto the Lagrangian, and it is unrelated to N[νµ] . Because of the coupling of Nνµ to the electromagnetic field in such theories, there would be observable violations[34, 35, 36] of the principle of equivalence for values of N[νµ] which occur in the theory. Such problems do not apply in our theory, mainly because we assume a symmetric metric which is defined in terms of Nνµ , and it is this symmetric metric which appears in Maxwell’s equations, and any coupling to additional fields. Such problems are also avoided in our theory partly because of the small values of N[νµ] which occur. In most previous work on the original Einstein-Schr¨odinger theory, the electromagnetic field is assumed to be the dual of N[τ ρ] . Even though this is the same definition used in [9, 10] to show there is no Lorentz force, several authors claim that a Lorentzlike force can be demonstrated[37, 38, 39]. However, the solutions[40, 41, 42] that 4

must be used for test particles have bad asymptotic behavior, such as a radial electric field which is independent of radius at large distances. Our theory uses a different definition of the electromagnetic field, and it has satisfactory exact solutions for both an electric monopole as in §3.1, and an electromagnetic plane-wave as in §3.2. Many others have contributed to the Einstein-Schr¨odinger theory. Of particular significance to our modified theory are contributions related to the choice of metric[43, 44, 45, 37], the generalized contracted Bianchi identity[43, 44, 45], the inclusion of sources[45, 19, 37], and exact solutions with a cosmological constant[46, 47]. This dissertation is organized as follows. In §2.1 we discuss the Lagrangian density. In §2.2-§2.4 we derive the field equations and quantify how closely they approximate the field equations of Einstein-Maxwell theory. In §3.1 we present an exact charged solution and show that it closely approximates the Reissner-Nordstr¨om solution. In §3.2 we present an exact electromagnetic plane-wave solution which is identical to its counterpart in Einstein-Maxwell theory. In §4.1 we derive the ordinary Lorentz force equation by taking the divergence of the Einstein equations when sources are included. In §4.2 we use the Lorentz force equation to derive the equations of motion for charged and neutral particles around the charged solution. In §4.3 we derive the Lorentz force using the EIH method, which requires no sources in the Lagrangian. In §5.1-§5.2 we calculate pericenter advance, deflection of light, and time delay of light, and compare the results to Einstein-Maxwell theory. In §5.3 we include a spin1/2 field in the Lagrangian and estimate the shift in Hydrogen atom energy levels for this theory as compared with Einstein-Maxwell theory. In §6.1 we represent the exact field equations in Newman-Penrose tetrad form, and use this to derive an exact 5

solution of the connection equations, and to show that the charged solution is Petrov type-D like the Reissner-Nordstr¨om solution. In §6.2 we derive the Newman-Penrose asymptotically flat O(1/r2 ) expansion of the field equations, and compare the results to Einstein-Maxwell theory. In §7.1-§7.3 we consider a generalization of our theory to non-Abelian fields, and show that a special case of the theory closely approximates Einstein-Weinberg-Salam theory.

6

Chapter 2 Extension of the Einstein-Schr¨ odinger theory for Abelian fields

2.1

The Lagrangian density

Einstein-Maxwell theory can be derived from a Palatini Lagrangian density, meaning that it depends on a connection Γλρτ as well as the metric gρτ , 1 √ −g [ g µν Rνµ (Γ) + 2Λb ] 16π 1√ + −gA[α,ρ] g αµ g ρνA[µ,ν] + Lm (uν , ψ, gµν , Aν · · · ). 4π

L(Γλρτ , gρτ , Aν ) = −

(2.1)

Here Λb is a bare cosmological constant. The Lm term couples the metric gµν and electromagnetic potential Aµ to additional fields, such as a hydrodynamic velocity vector uν , spin-1/2 wavefunction ψ, or perhaps the other fields of the Standard Model. The

7

bλ in place of original Einstein-Schr¨odinger theory allows a nonsymmetric Nµν and Γ ρτ the symmetric gµν and Γλρτ , and excludes the



−gA[α,ρ] g αµg ρνA[µ,ν] term (see Ap-

pendix M). Our “Λ-renormalized” Einstein-Schr¨odinger (LRES) theory introduces an additional cosmological term



−gΛz ,

i h 1 √ b + (n−2)Λb −N N aµν Rνµ (Γ) 16π 1 √ − −g (n−2)Λz + Lm (uν , ψ, gµν , Aν . . . ), 16π

bλ , Nρτ ) = − L(Γ ρτ

(2.2)

where Λb ≈ −Λz so that the total Λ matches astronomical measurements[48], Λ = Λb + Λz ≈ 10−56 cm−2 ,

(2.3)

and the physical metric and electromagnetic potential are defined to be √

√ −g g µν = −N N a(µν) ,

p bσ /[(n−1) −2Λb ]. Aν = Γ [νσ]

Equation (2.4) defines g µν unambiguously because



(2.4)

√ −g = [−det( −g g µν )]1/(n−2) .

Here and throughout this paper we use geometrized units with c = G = 1, the symbols

( )

and

[ ]

around indices indicate symmetrization and antisymmetrization, g =

det(gµν ), N = det(Nµν ), and N aσν is the inverse of Nνµ such that N aσνNνµ = δµσ . The dimension is assumed to be n=4, but “n” is retained in the equations to show how easily the theory can be generalized. The Lm term is not to include a



−gA[α,β] g αµ g βνA[µ,ν]

b is a form of part but may contain the rest of the Standard Model. In (2.2), Rνµ (Γ) non-symmetric Ricci tensor with special invariance properties to be discussed later, b =Γ bα − Γ bα bσ bα bσ bα bτ bα Rνµ (Γ) νµ,α (α(ν),µ) + Γνµ Γ(ασ) − Γνα Γσµ − Γ[τ ν] Γ[αµ] /(n−1).

(2.5)

bανµ is the Christoffel connection This tensor reduces to the ordinary Ricci tensor when Γ bα = 0 and Γ bα with Γ [νµ] α[ν,µ] = 0, as occurs in ordinary general relativity. 8

bα into a new connection Γ ˜ α , and Aσ from (2.4), It is helpful to decompose Γ νµ νµ p bανµ = Γ ˜ ανµ + (δµα Aν − δνα Aµ ) −2Λb , Γ σ σ b[σµ] b[σν] ˜ ανµ = Γ bανµ + (δµα Γ )/(n−1). − δνα Γ where Γ

(2.6) (2.7)

˜ α has the symmetry By contracting (2.7) on the right and left we see that Γ νµ ˜α , bα = Γ ˜α = Γ Γ αν (να) να

(2.8)

bανα had n3 . Substituting (2.6) so it has only n3 −n independent components whereas Γ into (2.5) as in R.17 gives p b = Rνµ (Γ) ˜ + 2A[ν,µ] −2Λb . Rνµ (Γ)

(2.9)

˜ α and Aσ , Using (2.9), the Lagrangian density (2.2) can be written in terms of Γ νµ h i p 1 √ ˜ νµ + 2A[ν,µ] −2Λb ) + (n−2)Λb −N N aµν (R 16π 1 √ − −g (n−2)Λz + Lm (uν , ψ, gµν , Aσ . . . ). (2.10) 16π

λ bρτ L(Γ , Nρτ ) = −

˜ νµ = Rνµ (Γ), ˜ and from (2.8,2.5) we have Here R ˜ νµ = Γ ˜ ανµ,α − Γ ˜ αα(ν,µ) + Γ ˜ σνµ Γ ˜ ασα − Γ ˜ σνα Γ ˜ ασµ . R

(2.11)

˜ α and Aν fully parameterize Γ bα and can be treated as independent From (2.6,2.8), Γ νµ νµ ˜ ανµ = 0 and variables. It is simpler to calculate the field equations by setting δL/δ Γ bα = 0, so we will follow this method. δL/δAν = 0 instead of setting δL/δ Γ νµ To do quantitative comparisons of this theory to Einstein-Maxwell theory we will need to use some value for Λz . One possibility is that Λz results from zero-point

9

fluctuations[15, 16, 17, 18], in which case using (2.3) we get Λb ≈ −Λz ∼ Cz ωc4 lP2 ∼ 1066 cm−2 ,

(2.12)

ωc = (cutoff frequency) ∼ 1/lP ,

(2.13)

Cz

³ ´ 60 1 fermion boson = − ∼ spin states spin states 2π 2π

(2.14)

where lP = (Planck length) = 1.6 × 10−33 cm. We will also consider the limit ωc → ∞, |Λz | → ∞, Λb → ∞ as in quantum electrodynamics, and we will prove that ³ ´ ³ ´ lim Λ-renormalized Einstein-Maxwell = . theory |Λz | → ∞ Einstein-Schr¨odinger theory

(2.15)

The non-symmetric Ricci tensor (2.5) has the following invariance properties bα ) = Rµν (Γ bα ), Rνµ (Γ ρτ τρ

(2.16)

bα + δ α ϕ,τ ] ) = Rνµ (Γ bα ) for an arbitrary ϕ(xσ ). Rνµ (Γ ρτ [ρ ρτ

(2.17)

From (2.16,2.17), the Lagrangians (2.2,2.10) are invariant under charge conjugation, ˜ ανµ → Γ ˜ αµν , Γ bανµ → Γ bαµν , Nνµ→ Nµν , N aνµ → N aµν, Q → −Q, Aσ → −Aσ , Γ

(2.18)

and also under an electromagnetic gauge transformation ψ → ψeiφ , Aα → Aα −

p h ¯ ˜α → Γ ˜α , Γ bα → Γ bα + 2¯h δ α φ,τ ] −2Λb , φ,α , Γ ρτ ρτ ρτ ρτ [ρ Q Q

(2.19)

bα , Nνµ ˜α , Γ assuming that Lm is invariant. With Λb > 0, Λz < 0 as in (2.12) then Γ νµ νµ ˜ νµ and Rνµ (Γ) b are Hermitian from (2.16), and gνµ , Aσ and N aνµ are all Hermitian, R and L are real from (2.4,2.2,2.10). In this theory the metric gµν from (2.4) is used for measuring space-time intervals, for calculating geodesics, and for raising and lowering of indices. The covariant 10

derivative “;” is always done using the Christoffel connection formed from gµν , Γανµ =

1 ασ g (gµσ,ν + gσν,µ − gνµ,σ ). 2

(2.20)

We will see that taking the divergence of the Einstein equations using (2.20,2.4) gives the ordinary Lorentz force equation. The electromagnetic field is defined in terms of the vector potential (2.4) Fµν = Aν,µ − Aµ,ν .

(2.21)

However, we will also define another field f µν √

√ √ 1/2 −g f µν = −N N a[νµ] Λb / 2 i.

(2.22)

√ −1/2 Then from (2.4), g µν and f µν 2 iΛb are parts of a total field, √ √ √ −1/2 ( −N / −g )N aνµ = g µν +f µν 2 iΛb .

(2.23)

We will see that the field equations require fµν ≈ Fµν to a very high precision. The bλ may definitions (2.4) of gµν and Aν in terms of the “fundamental” fields Nρτ , Γ ρτ seem unnatural from an empirical viewpoint. On the other hand, our Lagrangian density (2.2) seems simpler than (2.1) of Einstein-Maxwell theory, it contains fewer fields, and these fields have no symmetry restrictions. However, these are all very subjective considerations. It is much more important that our theory closely matches Einstein-Maxwell theory, and hence measurement. Note that there are many nonsymmetric generalizations of the Ricci tensor besides b from (2.5) and the ordinary Ricci tensor Rνµ (Γ). b For example, our version Rνµ (Γ) b Rµν (Γ), b Rνµ (Γ bT ) and Rµν (Γ bT ), and we could form any weighted average of Rνµ (Γ), 11

bα , Γ bα bα bσ bα bσ then add any linear combination of the tensors Γ α[ν,µ] [ν|α,|µ] , Γ[νµ] Γ[σα] , Γ[νσ] Γ[µα] , and bα Γ bσ bα Γ [αν] [σµ] . All of these generalized Ricci tensors would be linear in Γνµ,σ , quadratic in bα , and would reduce to the ordinary Ricci tensor when Γ bα = 0 and Γ bα Γ νµ [νµ] α[ν,µ] = 0, as occurs in ordinary general relativity. Even if we limit the tensor to only four terms, there are still eight possibilities. We assert that invariance properties like (2.16,2.17) are the most sensible way to choose among the different alternatives, not criteria such as the number of terms in the expression. α Finally, let us discuss some notation issues. We use the symbol Γνµ for the Christof-

˜ α and Γ bα fel connection (2.20) whereas Einstein and Schr¨odinger used it for our Γ νµ νµ respectively. We use the symbol gµν for the symmetric metric (2.4) whereas Einstein and Schr¨odinger used it for our Nµν , the nonsymmetric fundamental tensor. Also, to represent the inverse of Nαµ we use N aσα instead of the more conventional N ασ , because this latter notation would be ambiguous when using g µν to raise indices. While our notation differs from previous literature on the Einstein-Schr¨odinger theory, this change is required by our explicit metric definition, and it is necessary to be consistent with the much larger body of literature on Einstein-Maxwell theory.

2.2

The Einstein equations

√ To set δL/δ( −N N aµν ) = 0 we need some initial results. Using (2.4) and the identities det(sM ) = sn det(M ), det(M −1 ) = 1/det(M ) gives √ √ −N = (−det( −N N a.. ))1/(n−2) , √

√ √ −g = (−det( −g g .. ))1/(n−2) = (−det( −N N a(..) ))1/(n−2) . 12

(2.24) (2.25)

−1 Using (2.24,2.25,2.4) and the identity ∂(det(M ·· ))/∂M µν = Mνµ det(M ·· ) gives

√ ∂ −N Nνµ √ = , (n−2) ∂( −N N aµν )

√ ∂ −g gνµ √ = . (n−2) ∂( −N N aµν )

(2.26)

√ Setting δL/δ( −N N aµν ) = 0 using (2.10,2.26) gives the field equations, ·

µ ¶ ¸ ∂L ∂L √ √ 0 = −16π − ,ω ∂( −N N aµν ) ∂( −N N aµν ), ω √ ˜ νµ + 2A[ν,µ] 2 iΛ1/2 + Λb Nνµ + Λz gνµ − 8πSνµ , = R b

(2.27) (2.28)

where Sνµ and the energy-momentum tensor Tνµ are defined by δLm δLm Sνµ ≡ 2 √ =2 √ , aµν δ( −gg µν ) δ( −N N ) 1 1 Tνµ ≡ Sνµ − gνµ Sαα , Sνµ = Tνµ − gνµ Tαα . 2 (n − 2)

(2.29) (2.30)

The second equality in (2.29) results because Lm in (2.2) contains only the metric √

√ √ −g g µν = −N N a(µν) from (2.4), and not −N N a[µν] . Taking the symmetric and

antisymmetric parts of (2.28) and using (2.21) gives µ ˜ R(νµ) + Λb N(νµ) + Λz gνµ = 8π Tνµ − √ −1/2 ˜ [νµ] Λ−1 . N[νµ] = Fνµ 2 iΛb − R b

¶ 1 α gνµ Tα , (n − 2)

(2.31) (2.32)

Also from the curl of (2.32) we get ˜ [νµ,σ] + Λb N[νµ,σ] = 0. R

(2.33)

To put (2.31) into a form which looks more like the ordinary Einstein equations, we need some preliminary results. The definitions (2.4,2.22) of gνµ and fνµ can be inverted exactly to give Nνµ in terms of gνµ and fνµ . An expansion in powers of Λ−1 b

13

will better serve our purposes, and is derived in Appendix C, N(νµ)

µ = gνµ − 2 fν σ fσµ −

¶ 1 ρσ 4 −2 gνµ f fσρ Λ−1 b + (f )Λb . . . 2(n−2)

√ −1/2 −3/2 N[νµ] = fνµ 2 iΛb + (f 3 )Λb . . . .

(2.34) (2.35)

Here the notation (f 3 ) and (f 4 ) refers to terms like fνα f α σ f σ µ and fνα f α σ f σ ρ f ρ µ . Let us consider the size of these higher order terms relative to the leading order term for worst-case fields accessible to measurement. In geometrized units an elementary charge has r Qe = e

G = c4

r

e2 G¯h √ = α lP = 1.38 × 10−34 cm h ¯ c c3

where α = e2 /¯ hc is the fine structure constant and lP =

p

(2.36)

G¯ h/c3 is the Planck

length. If we assume that charged particles retain f 1 0 ∼ Q/r2 down to the smallest radii probed by high energy particle physics experiments (10−17 cm) we have from (2.36,2.12), |f 1 0 |2 /Λb ∼ (Qe /(10−17 )2 )2 /Λb ∼ 10−66 .

(2.37)

Here |f 1 0 | is assumed to be in some standard spherical or cartesian coordinate system. If an equation has a tensor term which can be neglected in one coordinate system, it can be neglected in any coordinate system, so it is only necessary to prove it in one coordinate system. The fields at 10−17 cm from an elementary charge would be larger than near any macroscopic charged object, and would also be larger than the strongest plane-wave fields. Therefore the higher order terms in (2.34-2.35) must be < 10−66 of the leading order terms, so they will be completely negligible for most purposes. 14

˜ α = 0. In §2.4 we will calculate the connection equations resulting from δL/δ Γ νµ Solving these equations gives (2.62,2.63,2.67,2.69), which can be abbreviated as −1 α ˜α Γ (νµ) = Γνµ + O(Λb ),

˜ α = O(Λ−1/2 ), Γ [νµ] b

(2.38)

˜ νµ = Gνµ + O(Λ−1 ), G b

˜ [νµ] = O(Λ−1/2 ), R b

(2.39)

˜ νµ = Rνµ (Γ), ˜ Rνµ = Rνµ (Γ) and where Γανµ is the Christoffel connection (2.20), R ˜ νµ = R ˜ (νµ) − 1 gνµ R ˜ ρ, G ρ 2 −1/2

In (2.39) the notation O(Λ−1 b ) and O(Λb −1/2

f[νµ,α]; α Λb

Gνµ = Rνµ −

1 gνµ R. 2

(2.40)

) indicates terms like f σ ν;α f α µ;σ Λ−1 b and

.

From the antisymmetric part of the field equations (2.32) and (2.35,2.39) we get fνµ = Fνµ + O(Λ−1 b ).

(2.41)

So fνµ and Fνµ only differ by terms with Λb in the denominator, and the two become identical in the limit as Λb → ∞. Combining the symmetric part of the field equations (2.31) with its contraction, and substituting (2.40,2.34,2.3) µ ¶ 1 1 σ ρ ρσ N(νµ) − gνµ Nρ = gνµ − 2 fν fσµ − gνµ f fσρ Λ−1 b 2 2(n−2) µ ¶ 1 1 ρσ ρσ 4 −2 − gνµ n + gνµ f fσρ − nf fσρ Λ−1 b + (f )Λb . . . 2 2(n−2) ³ n´ − 2fν σ fσµ Λ−1 = gνµ 1 − b 2 µ ¶ 1 n 4 −2 +gνµ +1− f ρσfσρ Λ−1 b + (f )Λb . . . (n−2) 2(n−2) µ ¶ ³n ´ 1 σ −1 ρσ = −2 fν fσµ − gνµ f fσρ Λb − − 1 gνµ + (f 4 )Λ−2 b ... 4 2

15

gives the Einstein equations µ ¶ ³n ´ 1 ρ ˜ Gνµ = 8πTνµ − Λb N(νµ) − gνµ Nρ +Λz − 1 gνµ , (2.42) 2 2 µ ¶ ³n ´ 1 σ ρσ = 8πTνµ + 2 fν fσµ − gνµ f fσρ + Λ − 1 gνµ + (f 4 )Λ−1 b . . . . (2.43) 4 2 From (2.29,2.30) we see that Tνµ will be the same as in ordinary general relativity, for example when we include classical hydrodynamics or spin-1/2 fields as in [49] or Appendix L. Therefore from (2.41,2.39), equation (2.43) differs from the ordinary Einstein equations only by terms with Λb in the denominator, and it becomes identical to the ordinary Einstein equations in the limit as Λb → ∞ (with an observationally valid total Λ). In §2.4 we will examine how close the approximation is for Λb from (2.12).

2.3

Maxwell’s equations

Setting δL/δAτ = 0 and using (2.10,2.22) gives · µ ¶ ¸ ∂L ∂L 4π 0 = √ − ,ω −g ∂Aτ ∂Aτ,ω √ √ 1/2 2 iΛb √ ( −gf ωτ ), ω a[ωτ ] τ √ √ = ( −N N ), ω − 4πj = − 4πj τ , 2 −g −g

(2.44) (2.45)

where j

τ

µ ¶ ¸ · −1 ∂Lm ∂Lm = √ − ,ω . −g ∂Aτ ∂Aτ,ω

(2.46)

From (2.45,2.21) we get Maxwell’s equations, f ωτ ; ω = 4πj τ ,

(2.47)

F[νµ,α] = 0.

(2.48)

16

µ where fνµ = Fνµ + O(Λ−1 will be b ) from (2.41). From (2.2,2.46) we see that j

the same as in ordinary general relativity, for example when we include classical hydrodynamics or spin-1/2 fields as in [49] or Appendix L. From (2.41), we see that equations (2.47,2.48) differ from the ordinary Maxwell equations only by terms with Λb in the denominator, and these equations become identical to the ordinary Maxwell equations in the limit as Λb → ∞. In §2.4 we will examine how close the approximation is for Λb from (2.12). Because Lm couples to additional fields only through gµν and Aµ , any equations associated with additional fields will be the same as in ordinary general relativity. For example in the spin-1/2 case, setting δL/δ ψ¯ = 0 will give the ordinary Dirac equation in curved space as in [49] or Appendix L. It would be interesting to investigate what ˜ α in Lm , although there does not appear to be any results if one includes fµν , Nµν or Γ µν empirical reason for doing so. A continuity equation follows from (2.47) regardless of the type of source, j ρ;ρ =

1 τρ f ;[τ ;ρ] = 0. 4π

(2.49)

Note that the covariant derivative in (2.47,2.49) is done using the Christoffel connection (2.20) formed from the symmetric metric (2.4).

2.4

The connection equations

˜ ανµ = 0 requires some preliminary calculations. With the definition Setting δL/δ Γ ∆L ∂L = − β ˜τρ ˜ βτρ ∆Γ ∂Γ

µ

∂L ˜ βτρ,ω ∂Γ

17

¶ ,ω

...

(2.50)

and (2.10,2.11) we can calculate, −16π

√ ∆L ρ ˜α ˜ σν[µ| δβα δστ δ ρ ) −N N aµν (δβσ δντ δ[µ| Γσ|α] + Γ = 2 |α] β ˜ ∆Γτ ρ √ √ ρ ω ρ ω −2( −N N aµν δβα δντ δ[µ δα] ), ω − ( −N N aµν δβα δατ δ[ν δµ] ),ω √ √ √ √ ˜ ρ −N N aµτ − Γ ˜ τ −N N aρν + Γ ˜ α −N N aρτ = −( −N N aρτ ), β − Γ νβ βα βµ √ √ √ ˜ τ −N N aµν ) + δ τ ( −N N a[ρω] ),ω , +δβρ (( −N N a ωτ ), ω + Γ νµ β

√ ∆L = (n−2)( −N N a[ρω] ), ω , ˜α ∆Γ αρ √ √ √ ∆L ˜ τ −N N aµν ) + ( −N N a[τ ω] ), ω . −16π = (n−1)(( −N N a ωτ ), ω + Γ νµ ˜α ∆Γ −16π

(2.51) (2.52) (2.53)

τα

In these last two equations, the index contractions occur after the derivatives. At ˜ α has the symmetry (2.8), it has only n3− n this point we must be careful. Because Γ νµ independent components, so there can only be n3 − n independent field equations associated with it. It is shown in Appendix B that instead of just setting (2.51) to zero, the field equations associated with such a field are given by the expression, # " δβρ δβτ ∆L ∆L ∆L − (2.54) 0 = 16π − ˜ ααρ (n−1) ∆Γ ˜ ατα ˜ βτρ (n−1) ∆Γ ∆Γ √ √ √ √ ˜ ρ −N N aµτ + Γ ˜ τ −N N aρν − Γ ˜ α −N N aρτ = ( −N N aρτ ), β + Γ νβ βα βµ √ −δβτ ( −N N a[ρω] ),ω +

√ √ 1 ((n−2)δβτ ( −N N a[ρω] ),ω +δβρ ( −N N a[τ ω] ),ω ) (n−1) √ √ √ √ ˜ τνβ −N N aρν + Γ ˜ ρ −N N aµτ − Γ ˜ αβα −N N aρτ = ( −N N aρτ ), β + Γ βµ √ √ 1 (δβτ ( −N N a[ρω] ),ω − δβρ ( −N N a[τ ω] ),ω ) (n−1) √ √ √ √ ˜ τσβ −N N aρσ + Γ ˜ ρ −N N aστ − Γ ˜ αβα −N N aρτ = ( −N N aρτ ), β + Γ βσ √ 8π 2 i √ τ] − −gj [ρ δβ . (2.55) 1/2 (n−1)Λb √ These are the connection equations, analogous to ( −gg ρτ );β = 0 in the symmet−

ric case. Note that we can also derive Ampere’s law (2.45) by antisymmetrizing 18

and contracting these equations.

From the definition of matrix inverse N aρτ =

(1/N )∂N/∂Nτ ρ , N aρτ Nτ µ = δµρ we get the identity √ √ √ √ ∂ −N −N aρτ −N aρτ ( −N ),σ = N Nτ ρ,σ = − N ,σ Nτ ρ . Nτ ρ,σ = ∂Nτ ρ 2 2

(2.56)

Contracting (2.55) with Nτ ρ using (2.8,2.56), and dividing this by (n−2) gives, √ √ ˜ α −N = − ( −N ), β − Γ αβ

√ 8π 2 i 1/2 (n−1)(n−2)Λb



−gj ρ N[ρβ] .

(2.57)

From (2.57) we get ˜α Γ α[ν,µ] −

√ 8π 2 i 1/2

(n−1)(n−2)Λb

µ√ ¶ √ −g ρ √ j N[ρ[ν] ,µ] = (ln −N ),[ν,µ] = 0. −N

(2.58)

From (2.55,2.57) we get the contravariant connection equations, N

aρτ

˜ τ aρσ + Γ ˜ ρ N aστ ,β + Γσβ N βσ

√ 8π 2 i

¶ √ µ −g [ρ τ ] 1 α aρτ √ = j δβ + j N[αβ] N .(2.59) 1/2 (n−2) −N (n−1)Λb

Multiplying this by −Nνρ Nτ µ gives the covariant connection equations, √ ¶ √ µ −8π 2 i −g 1 α α ˜ ˜ √ Nνµ,β − Γνβ Nαµ − Γβµ Nνα = N[αβ] Nνµ j α . (2.60) Nν[α Nβ]µ + 1/2 (n−2) −N (n−1)Λb Equation (2.60) together with (2.31,2.33,2.8) are often used to define the EinsteinSchr¨odinger theory, particularly when Tνµ = 0, j α = 0. Equations (2.55) or (2.60) can be solved exactly as in [50, 51] or §6.1, similar to the way gρτ ;β = 0 can be solved to get the Christoffel connection. An expansion in powers of Λ−1 b will better serve our purposes, and such an expansion is derived in

19

Appendix D, and is also stated without derivation in [45], ˜ ανµ = Γανµ + Υανµ , Γ · 1 α α Υ(νµ) = −2 f τ(ν fµ)α;τ + f ατfτ (ν;µ) + ((f ρσfσρ ), α gνµ − 2(f ρσfσρ ),(ν δµ) ) 4(n−2) µ ¶¸ 4π ρ 2 α 40 −2 α + j f ρ gνµ + fρ(ν δµ) Λ−1 b + (f )Λb . . . , (n−2) (n−1) · ¸ √ 8π 1 −1/2 −3/2 α α α α α Υ[νµ] = (fνµ; + f µ;ν − f ν;µ ) + j[ν δµ] 2 iΛb + (f 30 )Λb . . . 2 (n−1) · ¸ 1 8π α ρσ α 40 −2 Υαν = 2 (f fσρ ),ν − j fαν Λ−1 b + (f )Λb . . . . 2(n−2) (n−1)(n−2)

(2.61)

(2.62) ,(2.63) (2.64)

In (2.61), Γανµ is the Christoffel connection (2.20). The notation (f 30 ) and (f 40 ) refers to terms like f α τ f τ σ f σ [ν;µ] and f α τ f τ σ f σ ρ f ρ (ν;µ) . As in (2.34,2.35), we see from (2.37) that the higher order terms in (2.62-2.64) must be < 10−66 of the leading order terms, so they will be completely negligible for most purposes. Extracting Υτσβ of (2.61) from (2.11) gives (R.6,R.7), α σ α σ α σ α ˜ (νµ) = Rνµ + Υα R (νµ);α −Υα(ν;µ) −Υ(να) Υ(σµ) −Υ[να] Υ[σµ] +Υ(νµ) Υσα ,

(2.65)

˜ [νµ] = Υα −Υσ Υα −Υσ Υα +Υσ Υα . R [νµ];α (να) [σµ] [να] (σµ) [νµ] σα

(2.66)

20

Substituting (2.61-2.64,2.47) into (2.65) using ` = f ρσfσρ gives α α σ α ˜ (νµ) = Rνµ + Υ(νµ); R α − Υα(ν;µ) − Υ[να] Υ[σµ] . . . ·µ ¶ 1 τ α ατ α α = Rνµ − 2 f (ν fµ) ; τ + f fτ (ν; µ) + ( `, gνµ − 2 `,(ν δµ) ) ; α 4(n−2) µ ¶ 4π ρ 2 α α + j ;α f ρ gνµ + fρ(ν δµ) (n−2) (n−1) ¶ µ 2 4π ρ α α + j f ρ gνµ + fρ(ν δµ) ;α (n−2) (n−1) 1 8π + `,(ν; µ) − (j α fα(ν );µ) 2(n−2) (n−1)(n−2) µ ¶ 1 16π σ σ σ σ fνα; +f α;ν −f ν; α + − j[ν δα] 4 (n−1) ¶¸ µ 16π α α α α × fσµ; +f µ; σ −f σ; µ + j[σ δµ] Λ−1 b ... (n−1) · 1 = Rνµ − 2f τ (ν fµ) α ;τ ;α + 2f ατ fτ (ν; µ); α + `, α ;α gνµ 2(n−2) 1 −f σ ν;α f α µ;σ + f σ ν; α fσµ; α + f σ α;ν f α σ; µ 2 ¸ 2 32π 32π 2 ρ 8π ρ α τ +8πj fτ (ν;µ) − jν jµ + j jρ gνµ + j ;α f ρ gνµ Λ−1 b ..., (n−1) (n−2) (n−2) · 1 n ˜ ρ = R − 2f τ β fβ α ;τ ;α + `, α ;α − f σβ ;α f α β;σ + f σβ ; α fσβ; α R ρ 2(n−2) 2 µ ¶ ¸ 1 n 8πn ρ α ατ 2 ρ −8πf jτ ;α −32π 1+ − j jρ + j ;α f ρ Λ−1 b ... (n−1) (n−2) (n−2) · 3 n `, α ;α − f[σβ;α] f [σβ ; α] = R + −2f τ β fβ α ;τ ;α − 2(n−2) 2 ¸ 2 32π n 16π ατ ρ − j jρ − f jτ ;α Λ−1 b ... . (n−1)(n−2) (n−2)

and using (2.40) gives µ ˜ (Gνµ − Gνµ ) = − 2f τ (ν fµ) α ;τ ;α + 2f ατ fτ (ν; µ);α 1 − f σ ν;α f α µ;σ + f σ ν;α fσµ; α + f σ α;ν f α σ; µ 2 1 3 − gνµ f τ β fβ α ;τ ;α − gνµ (f ρσfσρ ), α ;α − gνµ f[σβ;α] f [σβ ; α] 4 4 ¶ 2 2 32π 16π τ ρ +8πj fτ (ν;µ) − jν jµ + gνµ j jρ Λ−1 b . . . . (2.67) (n−1) (n−1) 21

From (2.43) we can define an “effective” energy momentum tensor T˜νµ which applies when Gνµ is used in the Einstein equations and Lm = 0, 8π T˜νµ

¶ µ 1 σ ρσ ˜ νµ − Gνµ ). = 2 fν fσµ − gνµ f fσρ − (G 4

(2.68)

Substituting (2.63,2.47) into (2.66) gives −3/2 ˜ [νµ] = Υα )... R [νµ];α + O(Λb µ ¶ √ 1 8π −1/2 α α α = (fνµ; +f µ;ν −f ν;µ );α + j[ν,µ] 2 iΛb . . . 2 (n−1) µ ¶ √ 3 8π −1/2 α α α f[νµ,α]; + f µ;ν;α − f ν;µ;α + j[ν,µ] 2 iΛb . . . = 2 (n−1) µ ¶ √ 3 8π(n−2) −1/2 α α α = 2 iΛb . . . . (2.69) f[νµ,α]; +2f µ;[ν;α] −2f ν;[µ;α] − j[ν,µ] 2 (n−1)

As we have already noted in §2.2 and §2.3, the Λb in the denominator of (2.67,2.69) causes our Einstein and Maxwell equations (2.43,2.47,2.48) to become the ordinary Einstein and Maxwell equations in the limit as ωc → ∞, |Λz | → ∞, Λb → ∞, and it also causes the relation fνµ ≈ Fνµ from (2.41) to become exact in this limit. Let us examine how close these approximations are when Λb ∼ 1066 cm−2 as in (2.12). We will start with the Einstein equations (2.43). Let us consider worst-case values ˜ νµ − Gνµ accessible to measurement, and compare these to the ordinary electroof G magnetic term in the Einstein equations (2.43). If we assume that charged particles retain f 1 0 ∼ Q/r2 down to the smallest radii probed by high energy particle physics experiments (10−17 cm) we have, |f 1 0;1 /f 1 0 |2 /Λb ∼ 4/Λb (10−17 )2 ∼ 10−32 ,

(2.70)

|f 1 0;1;1 /f 1 0 |/Λb ∼ 6/Λb (10−17 )2 ∼ 10−32 .

(2.71)

22

−1 ατ So for electric monopole fields, terms like f σ ν;α f α µ;σ Λ−1 b and f fτ (ν; µ);α Λb in (2.67)

must be < 10−32 of the ordinary electromagnetic term in (2.43). And regarding j τ as a substitute for (1/4π)f ωτ ; ω from (2.47), the same is true for the jν terms. For an electromagnetic plane-wave in a flat background space we have Aµ = A²µ sin(kα xα ) , ²α ²α = −1 , k α kα = k α ²α = 0,

(2.72)

fνµ = 2A[µ,ν] = 2A²[µ kν] cos(kα xα ),

(2.73)

j σ = 0.

Here A is the magnitude, k α is the wavenumber, and ²α is the polarization. Substituting (2.72,2.73) into (2.67), all of the terms vanish for a flat background space. Also, for the highest energy gamma rays known in nature (1020 eV, 1034 Hz) we have from (2.12), |f 1 0;1 /f 1 0 |2 /Λb ∼ (E/¯ hc)2 /Λb ∼ 10−16 ,

(2.74)

|f 1 0;1;1 /f 1 0 |/Λb ∼ (E/¯ hc)2 /Λb ∼ 10−16 .

(2.75)

So for electromagnetic plane-wave fields, even if some of the extra terms in (2.67) were non-zero because of spatial curvatures, they must still be < 10−16 of the ordinary electromagnetic term in (2.68). Therefore even for the most extreme worst-case fields accessible to measurement, the extra terms in the Einstein equations (2.43) must all be < 10−16 of the ordinary electromagnetic term. Now let us look at the approximation fνµ ≈ Fνµ from (2.41), and Maxwell’s equations (2.47,2.48). From the covariant derivative commutation rule, the cyclic identity 2Rν[τ α]µ = Rνµατ , the definition of the Weyl tensor Cνµατ , and the Einstein equations

23

Rνµ = −Λgνµ + (f 2 ) . . . from (2.43) we get 1 2f α ν;[µ;α] = Rτ νµα f α τ + Rτ α µα f τ ν = Rνµατ f ατ + Rτ µ fτ ν 2 µ ¶ 1 4 2 [α τ ] [α τ ] ατ = Cνµ + δ R − δ δ R fατ − Rτ µ fντ 2 (n−2) [ν µ] (n−1)(n−2) [ν µ] (n−2)Λ 1 ατ = f Cατ νµ + fνµ + (f 3 ) . . . . (2.76) 2 (n−1) Substituting (2.35) into the antisymmetric field equations (2.32) gives √ ˜ [νµ] 2 iΛ−1/2 /2 + (f 3 )Λ−1 . . . , fνµ = Fνµ + R b b

(2.77)

and using (2.69,2.76) we get fνµ

µ ¶ 8π(n−2) 2(n−2)Λ 3 τα ατ = Fνµ + θ[τ,α] ενµ +f Cατ νµ + fνµ + j[ν,µ] +(f ) Λ−1 .. b .(2.78) (n−1) (n−1)

where ετ νµα = (Levi−Civita tensor), Cατ νµ = (Weyl tensor), and θτ =

1 2 f[νµ,α] ετ νµα , f[νµ,α] = − θτ ετ νµα . 4 3

(2.79)

The θ[τ,α] ενµ τ α Λ−1 b term in (2.78) is divergenceless so that it has no effect on Ampere’s law (2.47). The fνµ Λ/Λb term is ∼ 10−122 of fνµ from (2.3,2.12). The (f 3 )Λ−1 b term is < 10−66 of fνµ from (2.37). The largest observable values of the Weyl tensor might be expected to occur near the Schwarzschild radius, rs = 2Gm/c2 , of black holes, where it takes on values around rs /r3 . The largest value of rs /r3 would occur near the lightest black holes, which would be of about one solar mass, where from (2.12), 1 1 C0101 ∼ = 2 Λb Λb rs Λb

µ

c2 2Gm¯

¶2 ∼ 10−77 .

(2.80)

And regarding j τ as a substitute for (1/4π)f ωτ ; ω from (2.47), the j[ν,µ] Λ−1 b term is < 10−32 of fνµ from (2.71). Therefore, the last four terms in (2.78) must all be 24

< 10−32 of fνµ . Consequently, even for the most extreme worst-case fields accessible to measurement, the extra terms in Maxwell’s equations (2.47,2.48) must be < 10−32 of the ordinary terms. The divergenceless term θ[τ,α] ενµ τ α Λ−1 b of (2.78) should also be expected to be < 10−32 of fνµ from (2.70,2.71,2.79). However, we need to consider the possibility where θτ changes extremely rapidly. Taking the curl of (2.78), the Fνµ and j[ν,µ] terms drop out, µ f[νµ,σ] = θτ ;α;[σ ενµ]

τα

¶ 2(n−2)Λ 30 f[νµ,σ] + (f ) Λ−1 + (f Cατ [νµ ),σ] + b .... (n−1) ατ

Contracting this with Λb ερσνµ /2 and using (2.79) gives, 1 4(n−2)Λ ρ 2Λb θρ = −2 θ[ρ ; σ]; σ + ερσνµ (f ατ Cατ [νµ ),σ] + θ + (f 30 ) . . . 2 (n−1) Using θσ ;σ = 0 from the definition (2.79), the covariant derivative commutation rule, and the Einstein equations Rνµ = −Λgνµ + (f 2 ) . . . from (2.43), gives θσ ;ρ;σ = Rσρ θσ = −θρ Λ + (f 30 ) . . . , and we get something similar to the Proca equation[52, 53], µ θρ =

− θρ; σ;

σ

¶ 1 σνµ ατ 1 (3n−7)Λ 30 + ερ (f Cατ [νµ ),σ] + θρ + (f ) .... 2 (n−1) 2Λb

(2.81)

Here we are using a (1, −1, −1, −1) metric signature. Equation (2.81) suggests that θρ Proca-wave solutions might exist in this theory. Assuming that the magnitude of Cατ νµ is roughly proportional to θρ for such waves, and assuming that fµν goes according to (2.78) with Fµν = 0, the extra terms in (2.81) could perhaps be neglected in the weak field approximation. Using (2.81) and Λb ≈ −Λz = Cz ωc4 lP2 from (2.12), such Proca-wave solutions would have an extremely high minimum frequency ωP roca =

p p 2Λb ≈ 2Cz ωc2 lP ∼ 1043 rad/s, 25

(2.82)

where the cutoff frequency ωc and Cz come from (2.13,2.14). There are several points to make about (2.81,2.82). 1) A particle associated with a θρ field would have mass h ¯ ωP roca , which is much greater than could be produced by particle accelerators, and so it would presumably not conflict with high energy physics experiments. 2) We have recently shown that sin[ kr−ωt] Proca-wave solutions do not exist in the theory, using an asyptotically flat Newman-Penrose 1/r expansion similar to [54, 55]. However, it is still possible that wave-packet solutions could exist.

3) Substituting the k = 0 flat space Proca-wave solution

θρ = (0, 1, 0, 0)sin[ωP roca t] and Fµν = 0 into (2.78,2.68,2.67), and assuming a flat background space gives T˜00 = −2/Λb < 0. This suggests that Proca-wave solutions might have negative energy, but because sin[ kr − ωt] solutions do not exist, and because of the other approximations used, this calculation is extremely uncertain. 4) With a cutoff frequency ωc ∼ 1/lP from (2.13) we have ωP roca > ωc from (2.82,2.13,2.14), so Proca-waves would presumably be cut off. More precisely, (2.82) says that Proca√ waves would be cut off if ωc > 1/(lP 2Cz ). Whether ωc is caused by a discreteness, uncertainty or foaminess of spacetime near the Planck length[56, 57, 58, 59, 60], or by some other effect, the same ωc which cuts off Λz in (2.12) should also cut off very high frequency electromagnetic and gravitational waves, and Proca-waves. 5) If wavepacket Proca-wave solutions do exist, and they have negative energy, it is possible that θρ could function as a kind of built-in Pauli-Villars field. Pauli-Villars regularization in quantum electrodynamics requires a negative energy Proca field with a mass h ¯ ωP roca that goes to infinity as ωc → ∞, as we have from (2.82). This idea is supported by the effective weak field Lagrangian derived in Appendix J, and is discussed more 26

fully in Appendix K. 6) As mentioned initially, it might be more correct to take the limit of this theory as ωc → ∞, |Λz | → ∞, Λb → ∞, as in quantum electrodynamics. In this limit (2.81,2.82) require that θρ → 0 or ωP roca → ∞, and the theory becomes exactly Einstein-Maxwell theory as in (2.15). 7) Finally, we should emphasize that Proca-wave solutions are only a possibility suggested by equation (2.81). Their existence and their possible interpretation are just speculation at this point. We are continuing to pursue these questions.

27

Chapter 3 Exact Solutions

3.1

An exact electric monopole solution

Here we present an exact charged solution for this theory which closely approximates the Reissner-Nordstr¨om solution[61, 62] of Einstein-Maxwell theory. The solution is derived in Appendix N, and a MAPLE program[63] which checks the solution is available. The solution is 1 2 dr − cˇr2 dθ2 − cˇr2 sin2 θdφ2 , cˇa √ √ Q 2 f 10 = , −N = r sin θ, −g = cˇr2 sin θ, cˇr2 · µ ¶µ ¶¸ Q 4M 4Λ Q2 Vˆ Λ 0 F01 = −A0 = 2 1 + − + 2 cˇ − 1 − 1− , r Λb r3 3Λb Λb r4 Λb µ ¶ 2M Λr2 Q2 Vˆ Λ a = 1− − + 2 1− , r 3 r Λb ds2 = cˇadt2 −

(3.1) (3.2) (3.3) (3.4)

where (0 ) means ∂/∂r, and cˇ and Vˆ are very close to one for ordinary radii, s µ 2 ¶i 2Q2 Q2 (2i)! 2Q cˇ = 1− =1− ··· − 2 i , (3.5) 4 4 Λb r Λb r [i!] 4 (2i−1) Λb r4 µZ ¶ µ 2 ¶i rΛb r3 (2i)! Q2 2Q 2 ˆ V = r cˇ dr − ··· + , (3.6) =1+ 2 4 i Q 3 10Λb r i!(i+1)! 4 (4i+1) Λb r4 28

and the nonzero connections are aa0 cˇ2 4a2 Q2 a0 −a0 1 0 0 1 ˜ ˜ ˜ ˜ Γ00 = , Γ10 = Γ01 = − , Γ11 = , 2 Λb r5 2a 2a ˜2 = Γ ˜2 = Γ ˜3 = Γ ˜3 = 1 , Γ 12 21 13 31 r

(3.7)

˜ 1 = −ar , Γ ˜ 1 = −ar sin2 θ , Γ ˜3 = Γ ˜ 3 = cot θ , Γ ˜ 2 = −sin θcos θ, Γ 22 33 23 32 33 √ √ a 2 iQ 2a 2 2 3 3 1 1 ˜ = −Γ ˜ =Γ ˜ = −Γ ˜ = −√ ˜ = −Γ ˜ = − √ 2 iQ . Γ , Γ 02 20 03 30 10 01 Λb r 3 Λb r3 With Λz = 0, Λb = Λ we get the Papapetrou solution[46, 47] of the unmodified Einstein-Schr¨odinger theory. In this case the M/Λb r3 term in (3.3) would be huge from (2.3), and the Q2 /r2 term in (3.4) disappears, which is why the Papapetrou solution was found to be unsatisfactory[46]. However, we are instead assuming Λb ≈ −Λz from (2.12). In this case the solution matches the Reissner-Nordstr¨om solution except for terms which are negligible for ordinary radii. To see this, first recall that Λ/Λb ∼ 10−122 from (2.3,2.12), so the Λ terms are all extremely tiny. Ignoring the Λ terms and keeping only the O(Λ−1 b ) terms in (3.3,3.4,3.5,3.6) gives F01 = A0 = a = cˇ =

· ¸ Q 4M 4Q2 1+ − + O(Λ−2 b ), r2 Λb r 3 Λb r 4 · ¸ Q M 4Q2 1+ − + O(Λ−2 b ), r Λb r3 5Λb r4 · ¸ 2M Q2 Q2 1− + 2 1+ + O(Λ−2 b ), r r 10Λb r4 Q2 + O(Λ−2 1− b ). 4 Λb r

(3.8) (3.9) (3.10) (3.11)

For the smallest radii probed by high-energy particle physics we get from (2.37), Q2 ∼ 10−66 . Λb r4

(3.12)

The worst-case value of M/Λb r3 might be near the Schwarzschild radius rs of black holes where r = rs = 2M and M/Λb r3 = 1/2Λb rs2 . This value will be largest for the 29

lightest black holes, and the lightest black hole that we can expect to observe would be of about one solar mass, where we have M 1 1 ∼ = 3 2 Λb r 2Λb rs 2Λb

µ

c2 2Gm¯

¶2 ∼ 10−77 .

(3.13)

Also, an electron has M = Gme /c2 = 7 × 10−56 cm, and using (2.12) and the smallest radii probed by high-energy particle physics (10−17 cm) we have M 7 × 10−56 ∼ ∼ 10−70 . Λb r3 1066 (10−17 )3

(3.14)

From (3.12,3.13,3.14,2.3,2.12) we see that our electric monopole solution (3.1-3.4) has a fractional difference from the Reissner-Nordstr¨om solution of at most 10−66 for worst-case radii accessible to measurement. Clearly our solution does not have the deficiencies of the Papapetrou solution[46, 47] in the original theory, and it is almost certainly indistinguishable from the Reissner-Nordstr¨om solution experimentally. Also, from (6.142-6.147) the solution is of Petrov Type-D. And the solution reduces to the Schwarzschild solution for Q = 0. And from (3.8-3.11) we see that the solution goes to the Reissner-Nordstr¨om solution exactly in the limit as Λb → ∞. The only significant difference between our electric monopole solution and the Reissner-Nordstr¨om solution occurs on the Planck scale. From (3.1,3.5), the surface area of the solution is[64], µ

surface area



Z

Z

π

=



dθ 0

s √ dφ gθθ gφφ = 4πr2 cˇ = 4πr2

0

1−

2Q2 . Λb r4

(3.15)

The origin of the solution is where the surface area vanishes, so in our coordinates the origin is not at r = 0 but rather at r0 =

p Q(2/Λb )1/4 . 30

(3.16)

An alternative coordinate system is investigated in Appendix P where the origin is at ρ = 0, but it is less satisfactory in most respects than the one we are using. From (2.36,2.12) we have r0 ∼ lP ∼ 10−33 cm for an elementary charge, and r0 ¿ 2M for any realistic astrophysical black hole. For Q/M < 1 the behavior at the origin is hidden behind an event horizon nearly identical to that of the Reissner-Nordstr¨om solution. For Q/M > 1 where there is no event horizon, the behavior at the origin differs markedly from the simple naked singularity of the Reissner-Nordstr¨om solution. For the Reissner-Nordstr¨om solution all of the relevant fields have singularities at the origin, with g00 ∼ Q2 /r2 , A0 = Q/r, F01 = Q/r2 , R00 ∼ 2Q4 /r6 and R11 ∼ 2/r2 . For our √ √ solution the metric has a less severe singularity at the origin, with g11 ∼ − r/ r − r0 . Also, the fields Nµν , N aνµ ,

√ √ √ √ √ −N , Aν , −gf νµ , −gfνµ , −gg νµ , −ggνµ , and the

functions “a” and Vˆ all have finite nonzero values and derivatives at the origin, because it can be shown (see Appendix O) that Vˆ (r0 ) =



√ 2 [Γ(1/4)]2 /6 π −2/3 =

˜ αµν and √−g R ˜ νµ are also finite and nonzero at the origin, 1.08137. The fields Fνµ , Γ so if we use the tensor density form of the field equations (2.28,2.47), there is no ambiguity as to whether the field equations are satisfied at this location. −1/2

Finally let us consider the result from (2.37) that |f µ σ Λb

| < 10−33 for worst-case

electromagnetic fields accessible to measurement. The “smallness” of this value may √ −1/2 are part of the total seem unappealing at first, considering that g µν and f µν 2 iΛb √ √ √ −1/2 field ( −N / −g)N aνµ = g µν +f µν 2 iΛb as in (2.23). However, for an elementary −1/2

charge, |f µν Λb

| is not really small if one compares it to g µν − η µν instead of g µν .

Our charged solution (3.1,3.2,3.4) has g 00 ≈ 1+ 2M/r+ Q2/r2 and f 01 ≈ Q/r2 . So for −1/2

an elementary charge, we see from (2.36,2.12) that |f 01 Λb 31

| ∼ Q2/r2 for any radius.

3.2

An exact electromagnetic plane-wave solution

Here we present an exact electromagnetic plane-wave solution for this theory which is identical to the electromagnetic plane-wave solution in Einstein-Maxwell theory, sometimes called the Baldwin-Jeffery solution[65, 66, 67]. We will not do a full derivation, but a MAPLE program[63] which checks the solution is available. We present the solution in the form of a pp-wave solution[68], and a gravitational wave component is included for generality. The solution is expressed in terms of null coordinates √ √ x, y, u = (t − z)/ 2, v = (t + z)/ 2,  

gµν

−1   0  =   0    0

fµν = 2A[ν,µ]





0 0 fˇx   0      0  ˇ 0 0 f √  y √ µν  , −gf = 2  (3.17)   0  0 0 0      −fˇx −fˇy 0 0    0 0 −fˇx 0       √  0 0 −fˇy 0 √ √  , −g = −N = 1, (3.18) = 2A,[ν kµ] = 2   fˇ fˇ 0 0  x y      0 0 0 0 0

0 0   −1 0 0  ,  0 H 1    0 1 0

where √ kµ = (0, 0, −1, 0), Aµ = (0, 0, A, 0), A = − 2(xfˇx + y fˇy ),

(3.19)

ˆ + A2 H = 2H

(3.20)

= 2(h+ x2 + h× xy − h+ y 2 ) + 2(fˇx2 + fˇy2 )(x2 + y 2 ), ˆ = h+ x2 + h× xy − h+ y 2 + (y fˇx − xfˇy )2 , H 32

(3.21) (3.22)

and the nonzero connections are ˜ 1 = 1 ∂H , Γ 33 2 ∂x ˜ 4 = 1 ∂H − Γ 13 2 ∂x ˜ 4 = 1 ∂H − Γ 23 2 ∂y

1 ∂H 1 ∂H 2 ∂(fˇx2 + fˇy2 ) 2 4 ˜ ˜ Γ33 = , Γ33 = − , 2 ∂y 2 ∂u Λb ∂u ˇ 2i ∂ fˇx ˜ 4 = 1 ∂H + √2i ∂ fx , √ , Γ 31 2 ∂x Λb ∂u Λb ∂u ˇ 2i ∂ fˇy ˜ 4 = 1 ∂H + √2i ∂ fy . √ , Γ 32 2 ∂y Λb ∂u Λb ∂u

(3.23)

Here h+ (u), h× (u) characterize the gravitational wave component, fˇx (u), fˇy (u) characterize the electromagnetic wave component, and all of these are arbitrary functions √ of the coordinate u = (t − z)/ 2. ˜ µν = Rµν , and the elecFor the parameterization (3.17-3.20), it happens that R tromagnetic field is a null field[68, 66] with f σ µ f µ σ = det(f µ ν ) = 0. For this case, as shown in §6.1, all of the higher order terms in (2.34,2.35,2.43) vanish so that 1/2 √ Fµν = fµν = N[µν] Λb / 2 i and our Einstein and Maxwell equations are identical to

those of Einstein-Maxwell theory. Maxwell’s equations (2.47,2.48) are satisfied automatically from (3.18,3.17), and the Einstein equations reduce to, 2ˆ 2ˆ ˜ 33 + Λb (N33 − g33 ) = ∂ H + ∂ H − 2(fˇ2 + fˇ2 ). 0 = R x y ∂x2 ∂y 2

(3.24)

This has the solution (3.22,3.21). In Appendix Q the solution above is transformed to ordinary x, y, z, t coordinates, and also to the alternative x, y, u, v coordinates of [66]. The solution has been discussed extensively in the literature on Einstein-Maxwell theory[65, 68, 66, 67] so we will not interpret it further. It is the same solution which forms the incoming waves for the Bell-Szekeres colliding plane-wave solution[67], although the full Bell-Szekeres solution does not satisfy our theory because the electromagnetic field is not null after the collision. 33

Chapter 4 The equations of motion

4.1

The Lorentz force equation

A generalized contracted Bianchi identity for this theory can be derived using only ˜α , the connection equations (2.55) and the symmetry (2.8) of Γ νµ √ √ √ ˜ σλ + −N N aσν R ˜ λσ ),ν − −N N aνσ R ˜ σν,λ = 0. ( −N N aνσ R

(4.1)

This identity can also be written in a manifestly covariant form √ √ √ ˜ σλ + −N N aσν R ˜ λσ );ν − −N N aνσ R ˜ σν;λ = 0, ( −N N aνσ R

(4.2)

˜ νµ from (2.4,2.22,2.40), or in terms of g ρτ , f ρτ and G ³ ´√ ˜ σ = 3 f σρ R ˜ [σρ,ν] + f ασ;α R ˜ [σν] 2 iΛ−1/2 . G ν; σ b 2

(4.3)

The identity was originally derived[3, 7] assuming j ν = 0 in (2.55), and later expressed in terms of the metric (2.4) by [23, 43, 44, 37]. The derivation for j ν 6= 0 was first done[45] by applying an infinitesimal coordinate transformation to an invariant 34

integral, and it is also done in Appendix E using a much different direct computation method. Clearly (4.1,4.3) are generalizations of the ordinary contracted Bianchi √ √ identity 2( −g Rν λ ),ν− −g g νσRσν,λ = 0 or Gσν;σ = 0, which is also valid in this theory. Another useful identity[23] is derived in Appendix A using only the definitions (2.4,2.22) of gµν and fµν , ³

1 N ν) − δνµ Nρρ 2 (µ

´

³



´√ −1/2 3 σρ σρ = f N[σρ,ν] + f ;σ N[ρν] 2 iΛb . 2

(4.4)

The ordinary Lorentz force equation results from taking the divergence of the Einstein equations (2.42) using (4.3,2.47,2.32,4.4,2.21) ´ ´√ ³ 3 σρ ˜ ˜ [σν] 2 iΛ−1/2 + Λb N (µ ν) − 1 δ µ N ρ ; µ f R[σρ,ν] + 4πj σ R b 2 2 ν ρ ³ ´ ´√ ³ ˜ [σν] − Λb 3 f σρ N[σρ,ν] 2 iΛ−1/2 + Λb N (µ ν) − 1 δνµ Nρρ ; µ = 4πj σ R b 2 2 √ ˜ [σν] + Λb f ρσ ;ρ N[σν] ) 2 iΛ−1/2 = (4πj σ R b

8πTν;σ σ =

³

(4.5) (4.6) (4.7)

√ ˜ [σν] + Λb N[σν] ) 2 iΛ−1/2 = 4πj σ (R b

(4.8)

= 16πj σ A[σ,ν] ,

(4.9)

Tν;σ σ = Fνσ j σ .

(4.10)

See Appendix H for an alternative derivation of the Lorentz force equation. In Appendix L we also show that the Lorentz force equation and the continuity equation can be derived from the Klein-Gordon equation for spin-0 fields. Note that the covariant derivatives in (4.2,4.3,4.4,4.10) are all done using the Christoffel connection (2.20) formed from the symmetric metric (2.4).

35

4.2

Equations of motion of the electric monopole solution

Here we calculate the equations of motion when one body is much heavier than the other so this body remains approximately stationary and is represented by the charged solution (3.1-3.7). We ignore radiation reaction effects. The Lorentz-force equation (4.10) for the classical hydrodynamics case is duα Q2 F µ u = + Γαµν uµ uν , dλ α

µ

dxν u = . dλ ν

(4.11)

The stationary and moving bodies have masses M , M2 and charges Q, Q2 . We are using dλ = ds/M2 instead of ds because the unitless parameter λ is still meaningful for the null geodesics of photons where ds → 0 and M2 → 0. Using the metric (3.1) and the relation (r2 cˇ)0 = 2r/ˇ c from (3.5), the non-zero Christoffel connections (2.20) are Γ100 =

aˇ c (aˇ c)0 (aˇ c)0 1 (aˇ c)0 , Γ010 = , Γ111 = − , Γ212 = Γ313 = 2 , 2 2aˇ c 2aˇ c cˇ r

(4.12)

Γ122 = −ar , Γ133 = −arsin2 θ , Γ323 = cot θ , Γ233 = −sin θcos θ. The equations of motion (4.11) are then c)0 r2 aˇ c(aˇ c)0 t2 dur (aˇ − u − aruθ2 − arsin2 θuφ2 + u , dλ 2aˇ c 2 duθ 2 0 = + 2 ur uθ − sin θcosθ uφ2 , dλ rˇ c φ du (r2 cˇ)0 0 = + 2 ur uφ + 2cotθ uθ uφ , dλ r cˇ t du (aˇ c)0 r t Q2 F01 r u = + uu. aˇ c dλ aˇ c

aˇ cQ2 F01 ut =

36

(4.13) (4.14) (4.15) (4.16)

For motion in the equatorial plane we may put uθ = 0, θ = π/2, and (4.14) is identically satisfied. Then from (4.15) we get 0 =

1 d(uφ r2 cˇ) , r2 cˇ dλ

(4.17)

uφ r2 cˇ = (constant) = L = (angular momentum).

(4.18)

From (4.16,3.3) we get µ



¢ 1 d¡ t u aˇ c + Q2 A0 , aˇ c dλ

(4.19)

ut aˇ c + Q2 A0 = (constant) = E = (total energy).

(4.20)

1 0= aˇ c

d(ut aˇ c) + Q2 A00 ur dλ

=

Recalling that dλ = ds/M2 and uθ = 0, θ = π/2 we also have M22 = uα uα = aˇ cut2 −

1 r2 u −ˇ cr2 uφ2 . aˇ c

(4.21)

Eliminating t and λ from (4.21) using (4.18,4.20) gives µ

M22 aˇ c

dr φ = (aˇ c) u − u dφ 2 t2

¶2

µ

2 2 φ2

−aˇ cr u

dr L = (E − Q2 A0 ) − dφ r2 cˇ 2

¶2

aL2 − 2 . (4.22) r

This can be rewritten as an integral Z φ =

Ldr/r2 p . cˇ (E − Q2 A0 )2 −aL2 /r2 − M22 aˇ c

(4.23)

For Λb → ∞, a = 1 we have flat-space electrodynamics, and the integral can be done analytically. For Λb → ∞ we have Einstein-Maxwell theory, and the integral becomes an elliptic integral. For a finite Λb the integral is more complicated, but using (3.9,3.10,3.11) for A0 , a, cˇ and neglecting powers higher than 1/r4 also leads to an elliptic integral. The time dependence can be obtained using (4.18,4.20) to get dt/dφ = ut /uφ = (r2 cˇ/L)(E − Q2 A0 )/aˇ c = (E − Q2 A0 )r2 /aL 37

(4.24)

so that from (4.23), Z

(E − Q2 A0 )dr p aˇ c (E − Q2 A0 )2 −aL2 /r2 − M22 aˇ c

t =

(4.25)

We can also obtain the results (4.23,4.25) using the Hamilton-Jacobi approach as in [75], p. 94-95 and 306-308. From (4.21,3.1), the Hamilton-Jacobi equation is µ

M22

¶µ ¶ ∂S ∂S = g + Q2 A µ + Q2 Aν ∂xµ ∂xν µ ¶2 µ ¶2 µ ¶2 1 1 ∂S ∂S ∂S − 2 . = + Q2 A0 − aˇ c aˇ c ∂t ∂r cˇr ∂φ µν

(4.26) (4.27)

The solution is Z S = −Et + Lφ + Sr (r),

Sr =

dr aˇ c

q (E − Q2 A0 )2 −aL2 /r2 − M22 aˇ c . (4.28)

Then (4.23) is obtained from the equation ∂S/∂L = (constant), and (4.25) is obtained from the equation ∂S/∂E = (constant). Let us analyze the special case L = Q2 = 0 using the effective potential method of [69, 66]. From (4.22,4.18) and the definitions E˜ = E/M2 , ds = M2 dλ we have µ ¶2 dr 2 ˜ aˇ c=E − . ds

(4.29)

This equation can be expressed as a non-relativistic potential problem, 1 2

µ

dr ds

¶2 =

E˜ 2 − 1 ˜ −V, 2

(4.30)

where (dr/ds)2 /2 corresponds to the kinetic energy per mass, and V˜ is the so-called “effective potential”, aˇ c−1 V˜ = . 2 38

(4.31)

Since Vˆ (r0 ) =



√ 2 [Γ(1/4)]2 /6 π − 2/3 = 1.08137, we can assume that Vˆ ≈ 1 for

present purposes, and the effective potential becomes, µ ¶s 2 2M Q 2Q2 1 1 1− + 2 1− − . V˜ ≈ 4 2 r r Λb r 2

(4.32)

Let us consider the case for elementary particles where Q À M . This case is more interesting than astronomical objects because there is no event horizon to hide the behavior close to the origin at r0 = (2Q2 /Λb )1/4 = 3.16 × 10−34 cm where cˇ = 0. Assuming an electron charge and mass we have Q = Qe = e

p

G/c4 =



α lP = 1.38 ×

10−34 cm and M = Gme /c2 = 7×10−56 cm. In this case the mass term in “a” is smaller that the charge term for r < Q2 /2M = 1.36 × 10−13 cm, which is close to the classical electron radius. The following table shows the rough behavior of V˜ , V˜ vs. r for three Q/Qe values Our charged solution r/Qe \Q/Qe 1.68 2.66 3.76 4.60 5.32 5.94 1021 1025 ∞

Reissner-Nordstr¨om solution

.265 1.06 1.86 −.05602 − − −.00519 −.15032 − −.00002 −.00478 −.05315 .00055 .00769 .01525 .00062 .00953 .02611 .00060 .00937 .02722 −.00001 .00000 .00001 = − ≈ = Using (3.16) and Qe =



Q2e

α lP from (2.36), the Me term is insignificant and we get

√ √ µ 2¶ √ 4Q2e 2Q2e Λb 4Q2e 2Qe Qe 48 2α lP √ √ . =− ∆E0 ≈ − 3 − 3 √ a0 Λb 5 a0 2a0 Λb 5a20 Λb 2Qe

(5.150)

The term in the parenthesis is the ground state energy of a Hydrogen atom. With E0 = e2 /2a0 ∼ 13.6eV , lP = 1.6 × 10−33 cm, Λb ∼ 1066 cm−2 , h ∼ 4 × 10−15 eV · s, and a0 = h ¯ 2 /me e2 ∼ 5 × 10−9 cm we get √ 48 2α lP ∆E0 e2 −50 −50 √ 10 ∼ 10−49 eV, ∆f0 ∼ ∼ 10−34 Hz. (5.151) ∼ 10 , ∆E ∼ 0 2 2a0 h 5a0 Λb This is clearly unmeasurable.

70

Chapter 6 Application of Newman-Penrose methods

6.1

Newman-Penrose methods applied to the exact field equations

Here we use Newman-Penrose tetrad formalism to derive several results. In particular, we derive an exact solution for Nσµ in terms of gσµ and fσµ , and an exact solution of the connection equations (2.59), and we confirm the approximate solutions (2.34,2.35) and (2.62,2.63). We also derive the spin coefficients and Weyl tensor components for our charged solution (3.1-3.7), and show that it has Petrov type-D classification. Throughout this section, Latin letters indicate tetrad indices and Greek letters indicate tensor indices, and we assume n = 4 and the definitions √ −1/2 fˆνµ = f νµ 2 iΛb ,

√ ˆj ν = j ν 2 iΛ−1/2 , b

71

√ ˆ = Q 2 iΛ−1/2 . Q b

(6.1)

Using the definitions (2.4,2.22) we have W

σµ

√ −N = √ N aµσ = g σµ + fˆσµ . −g

(6.2)

Let us consider the following theorem which is similar to one in [51]: Theorem:Assume W σµ is a real tensor, fˆσµ = W [σµ] , and g σµ = W (σµ) is an invertible metric which can be put into Newman-Penrose tetrad form 

0   0  ,  −1    0

(6.3)

lσ = e1 σ , nσ = e2 σ , mσ = e3 σ , m∗ σ = e4 σ ,

(6.4)

lσ = e2 σ , nσ = e1 σ , mσ = −e4 σ , m∗σ = −e3 σ ,

(6.5)

δµσ = ea σ ea µ

, δba = eb σ ea σ ,

(6.6)

e = det(ea ν ) = ²αβσµ lα nβ mσ m∗µ

(6.7)

gab =

g ab = g αβ ea α eb β

e∗ = −e.

0 1 0   1 0 0  =  0 0 0    0 0 −1



(6.8)

where lσ and nσ are real, mσ and m∗σ are complex conjugates. Then tetrads ea ν may

72

be chosen such that 

W ab = W αβ ea α eb β



fˆab

0   −ˇ  u =   0    0



(1+ uˇ) 0 0   0     (1− uˇ) 0 0 0    , =    0  0 0 −(1+i` u )       0 0 −(1−i` u) 0    

0 uˇ 0      u 0 0 0  a  0 −ˇ , fˆ b =    0 0 0 0 −i` u       0 0 0 i` u 0

uˇ 0

0  0     uˇ 0 0    , fˆab =    0 i` u 0        0 0 −i` u 0

(6.9)



−ˇ u

0

0   0 0 0  , (6.10)  0 0 i` u    0 −i` u 0

where u` and uˇ are real, except for null fields with fˆσ µ fˆµ σ = det(fˆµ ν ) = 0, in which case tetrads may be chosen such that 

W ab = W αβ ea α eb β



fˆab

0   0  =   0    0

0 0

0   1  =  0    0 

0   0 −´ u −´ u  ,  u´ 0 0     u´ 0 0

1



0

0   0 −´ u −´ u  ,  u´ 0 −1     u´ −1 0  0   0  fˆa b =   u´    u´

(6.11)



0 0 0   0 u´ u´  ,  0 0 0    0 0 0

  0    0  fˆab =   −´  u   −´ u

 0 u´ u´   0 0 0  , (6.12)  0 0 0    0 0 0

where u´ is real. If W σµ is instead Hermitian, things are unchanged except that the scalars “` u” (u grave) and “ˇ u” (u check) are imaginary instead of real. The above theorem is proven in Appendix T. 73

One difference from the usual Newman-Penrose formalism is that gauge freedom is restricted so that only type III tetrad transformations can be used. Covariant derivative is done in the usual fashion, T a b|c = ea σ eb µ T σ µ;τ ec τ = T a b,c + γ a dc T d b − γ d bc T a d .

(6.13)

For the spin coefficients we will follow the conventions of Chandrasekhar[64], 1 γabc = (λabc + λcab − λbca ) = ea µ ebµ;σ ec σ , 2

(6.14)

γabc = −γbac , γ a ac = 0,

(6.15)

λabc = (ebσ,µ − ebµ,σ )ea σ ec µ = ebσ,µ (ea σ ec µ − ea µ ec σ ) = γabc − γcba ,

(6.16)

λabc = −λcba ,

(6.17)

ρ = γ314 , µ = γ243 , τ = γ312 , π = γ241 ,

(6.18)

κ = γ311 , σ = γ313 , λ = γ244 , ν = γ242 ,

(6.19)

² = (γ211 + γ341 )/2

, γ = (γ212 + γ342 )/2,

(6.20)

α = (γ214 + γ344 )/2

, β = (γ213 + γ343 )/2.

(6.21)

With these coefficients and with other tetrad quantities, complex conjugation causes the exchange 3 → 4, 4 → 3. As usual we may also define directional derivative operators, D = e1 α

∂ ∂xα

,

∆ = e2 α

∂ ∂xα

,

δ = e3 α

∂ ∂xα

,

δ ∗ = e4 α

∂ . (6.22) ∂xα

Substituting (2.65,2.66) into (2.28) gives the Einstein equations and antisymmetric

74

field equations in tetrad form µ

Rbd

1 = 8πG Tbd − gbd Taa 2

¶ − Λb N(bd) − Λe gbd

a c a c a c −Υ(bd)|a + Υaa(b|d) + Υ(ba) Υ(cd) + Υ[ba] Υ[cd] − Υ(bd) Υaca ,

(6.23)

√ 1/2 a c a c a c Λb N[bd] = 2A[d|b] 2 iΛb − Υ[bd]|a + Υ(ba) Υ[cd] + Υ[ba] Υ(cd) − Υ[bd] Υaca . (6.24) The usual Ricci identities will be valid if we define Φab values in terms of the righthand side of (6.23), Φ00 = −R11 /2, Φ22 = −R22 /2, Φ02 = −R33 /2, Φ20 = −R44 /2,

(6.25)

Φ01 = −R13 /2, Φ10 = −R14 /2, Φ12 = −R23 /2, Φ21 = −R24 /2,

(6.26)

ˆ = R/24 = (R12 − R34 )/12. Φ11 = −(R12 + R34 )/4, Λ

(6.27)

First let us consider the case where we do not have fˆσ µ fˆµ σ = det(fˆµ ν ) = 0. It is easily verified from (6.10) that the scalars are given by[51] u` = uˇ =

q√ q√

$ − `/4 ,

(6.28)

$ + `/4 ,

(6.29)

where $ = ( `/4)2 − fˆ/g,

(6.30)

fˆ = det(fˆµν ) , g = det(gµν ),

(6.31)

fˆ/g = −ˇ u2 u`2 , u2 − u`2 ). ` = fˆσ µ fˆµ σ = 2(ˇ

75

(6.32) (6.33)

From (6.2), the fundamental tensor of the Einstein-Schr¨odinger theory is,  

N aab

Nbc

0 0   0 (1− uˇ)     (1+ uˇ)  √ 0 0 0   −g¦  , = √  −N¦  0  0 0 −(1−i` u )       0 0 −(1+i` u) 0 

(6.34)



0 (1− uˇ)ˇ c/` c 0 0        (1+ uˇ)ˇ c/` c 0 0 0   , =     0 0 0 −(1−i` u)` c/ˇ c       0 0 −(1+i` u)` c/ˇ c 0

(6.35)

where c` = √

√ 1 = 1−` s2 , 1+ u`2

u` = s`/` c, cˇ = √

√ 1 s2 , = 1+ˇ 1− uˇ2

uˇ = sˇ/ˇ c, p

p i −det(Nab ) = , cˇc` p = −det(gab ) = i.

−N¦ =



−g¦

(6.36) (6.37) (6.38) (6.39) (6.40) (6.41)

Note the correspondence of s`, c`, u` and sˇ, cˇ, uˇ to circular and hyperbolic trigonometry functions.

76

From (6.10,6.13), Ampere’s law (2.47) becomes, 4π ˆc b ˆac c ˆba j = fˆbc ,b + γab f +γab f , c 4π ˆ2 3 ˆ12 4 ˆ12 2 ˆ43 2 ˆ34 f + γ14 f + γ34 f + γ43 f j = fˆ12 ,1 + γ13 c

(6.42) (6.43)

= Dˇ u − ρ∗ uˇ − ρˇ u − ρi` u + ρ∗ i` u

(6.44)

= Dˇ u − ρw − ρ∗ w∗ ,

(6.45)

4π ˆ1 4 ˆ21 1 ˆ43 1 ˆ34 3 ˆ21 f + γ24 f + γ34 f + γ43 f j = fˆ21 ,2 + γ23 c

(6.46)

= −∆ˇ u − µˇ u − µ∗ uˇ + µ∗ i` u − µi` u

(6.47)

= −∆ˇ u − µw − µ∗ w∗ ,

(6.48)

4π ˆ4 4 ˆ12 4 ˆ21 2 ˆ34 1 ˆ34 f f + γ21 f + γ12 f + γ32 j = fˆ34 ,3 + γ31 c

(6.49)

= −iδ` u − π ∗ i` u + τ i` u + τ uˇ + π ∗ uˇ

(6.50)

= −iδ` u + τ w + π ∗ w∗ ,

(6.51)

where w = uˇ + i` u

(6.52)

The connection equations are easier to work with in contravariant form (2.59) than √ √ in covariant form (2.55). Multiplying (2.59) by −N / −g and using (6.34,6.13,2.61)

77

and (6.40,6.41) gives 0=

Obcd

√ ¢ −N¦ ¡ acd d c d c = √ N ,b + γab N aca + γab N aad + Υab N aca + Υba N aad −g¦ ¶ µ 8π ˆ[d c] 1 ˆa acd + j δb − j N[ab] N , 3c 2

(6.53)

0 = Ob11 =

1 1 Υ2b (1− uˇ) + Υb2 (1+ uˇ),

(6.54)

0 = Ob22 =

2 2 (1− uˇ), (1+ uˇ) + Υb1 Υ1b

(6.55)

3 3 0 = Ob33 = −Υ4b (1−i` u) − Υb4 (1+i` u), √ µ ¶ ( −N¦ ),b ± 2 ± 1 ± 12 0 = Ob = ∓ˇ u,b + (1∓ uˇ) − √ + Υ2b + Υb1 −N¦ µ ¶ 1 ˆa 8π [2 1] ˆ + ±j δb − j N[ab] c`cˇ(1∓ uˇ) 3c 2 ¡ ¢ 2 1 = (1∓ uˇ) ∓ˇ u,b cˇ2 − u`u`,b c`2 + ±Υ2b + ±Υb1 µ ¶ ±1 ˆ[2 1] 8π (1∓ uˇ) ˆ[4 3] + j δb + i` u j δb , 3c (1± uˇ) (1+ u`2 ) µ √ ¶ ( −N¦ ),b ± 4 ± 3 ± 34 √ 0 = Ob = ±i` u,b + (1∓i` u) − Υ4b − Υb3 −N¦ µ ¶ 1 ˆa 8π [4 3] ˆ ±j δb + j N[ab] c`cˇ(1∓i` u) + 3c 2 ¡ ¢ 4 3 = (1∓i` u) ±i` u,b c`2 − uˇuˇ,b cˇ2 − ±Υ4b − ±Υb3 µ ¶ ±1 ˆ[4 3] 8π (1∓i` u) ˆ[2 1] + j δb + uˇ j δb , 3c (1±i` u) (1− uˇ2 )

(6.56)

(6.57)

(6.58)

(6.59)

(6.60)

4 2 0 = ±Ob24 = γ31b (−(1± uˇ) + (1∓i` u)) + ±Υ1b (1± uˇ) − ±Υb3 (1∓i` u) ± 4 2 = ∓γ31b w + ±Υ1b (1± uˇ) − ±Υb3 (1∓i` u) ±

8π ˆ[4 2] j δb , 3c

3 1 0 = ±Ob13 = γ24b ((1∓ uˇ) − (1±i` u)) + ±Υ2b (1∓ uˇ) − ±Υb4 (1±i` u) ±

u) ± = ∓γ24b w + ±Υ32b (1∓ uˇ) − ±Υ1b4 (1±i`

8π ˆ[4 2] j δb (6.61) 3c (6.62)

8π ˆ[3 1] j δb 3c

8π ˆ[3 1] j δb , 3c

(6.63) (6.64)

To save space in the equations above we are using the notation, −

Obdc =

+

Obcd = Obcd , 78



Υdcb = + Υdbc = Υdbc .

(6.65)

The connection equations (6.54-6.64) can be solved by forming linear combinations of them where all of the Υbca terms cancel except for the desired one. The calculations are done in Appendix U. Splitting the result into symmetric and antisymmetric components gives 2 = Υ(12) 1 = Υ(12) 4 Υ(34) = 1 Υ(11) = 2 Υ(22) = 3 Υ(33) =

4πˇ c2 uˇ ˆ2 j , 3c 4πˇ c2 uˇ ˆ1 cˇ2 uˇ∆ˇ u+ j , 3c 4π` c2 i` u ˆ4 −` c2 u`δ` u+ j , 3c 4πˇ ucˇ2 ˆ2 u`D` uc`2 − uˇDˇ ucˇ2 + j , 3c 4πˇ ucˇ2 ˆ1 j , u`∆` uc`2 − uˇ∆ˇ ucˇ2 − 3c 4πi` uc`2 ˆ4 u`δ` uc`2 − uˇδˇ ucˇ2 − j , 3c cˇ2 uˇDˇ u−

2 1 3 Υ(11) = Υ(22) = Υ(44) = 0, 2 Υ(23) = 1 Υ(13) = 3 Υ(13) = 3 Υ(23) = 4 Υ(12) = 2 = Υ(34) 1 Υ(43) = 2 Υ(13) = 4 = Υ(13) 4 = Υ(11)

2πi` uc`2 ˆ4 i` u (δˇ ucˇ2 − iδ` uc`2 ) − j , 2 3c i` u 2πi` uc`2 ˆ4 − (δˇ ucˇ2 + iδ` uc`2 ) − j , 2 3c uˇ 2πˇ ucˇ2 ˆ2 − (Dˇ ucˇ2 + iD` uc`2 ) + j , 2 3c uˇ 2πˇ ucˇ2 ˆ1 − (∆ˇ ucˇ2 − i∆` uc`2 ) − j , 2 3c µ ¶ uˇcˇ2 cˇ2 − δˇ u 2 + τ w − π∗ w∗ , 2 c` µ ¶ 2 i` uc` c`2 ∗ ∗ − iD` u 2 + ρw − ρ w , 2 cˇ µ ¶ 2 i` uc` c`2 ∗ ∗ − i∆` u 2 − µw + µ w , 2 cˇ νwˇ u κwˇ u 1 , Υ(24) =− , zˇ zˇ σwi` u λwi` u 3 =− , Υ(24) , z` z` κw2 νw2 3 =− , Υ(22) , zˇ zˇ 79

(6.66) (6.67) (6.68) (6.69) (6.70) (6.71) (6.72) (6.73) (6.74) (6.75) (6.76) (6.77) (6.78) (6.79) (6.80) (6.81) (6.82)

2 Υ(33) =

σw2 z`

1 , Υ(44) =−

λw2 , z`

4πˇ c2 ˆ2 j , 3c 4πˇ c2 ˆ1 −ˇ c2 ∆ˇ u− j , 3c 4π` c2 ˆ4 −i` c2 δ` u− j , 3c 1 2π` c2 ˆ4 j , (δˇ ucˇ2 − iδ` uc`2 ) − 2 3c 1 2π` c2 ˆ4 − (δˇ ucˇ2 + iδ` uc`2 ) − j , 2 3c 1 2πˇ c2 ˆ2 (Dˇ ucˇ2 + iD` uc`2 ) − j , 2 3c 1 2πˇ c2 ˆ1 − (∆ˇ ucˇ2 − i∆` uc`2 ) − j , 2 3c µ ¶ cˇ2 cˇ2 δˇ u 2 + τ w − π ∗ w∗ , 2 c` µ ¶ 2 c` c`2 ∗ ∗ iD` u 2 + ρw − ρ w , 2 cˇ µ ¶ 2 c` c`2 ∗ ∗ − i∆` u 2 − µw + µ w , 2 cˇ κw νw 1 − , Υ[24] =− , zˇ zˇ σw λw 3 , Υ[24] = , z` z`

(6.83)

2 Υ[12] = −ˇ c2 Dˇ u+

(6.84)

1 Υ[12] =

(6.85)

4 Υ[34] = 2 = Υ[23] 1 Υ[13] = 3 Υ[13] = 3 Υ[23] = 4 Υ[12] = 2 Υ[34] = 1 Υ[43] = 2 Υ[13] = 4 Υ[13] =

(6.86) (6.87) (6.88) (6.89) (6.90) (6.91) (6.92) (6.93) (6.94) (6.95)

where z` = [(1±i` u)2 (1± uˇ) + (1∓i` u)2 (1∓ uˇ)]/2 = 1 + 2iˇ uu` − u`2 ,

(6.96)

zˇ = [(1± uˇ)2 (1±i` u) + (1∓ uˇ)2 (1∓i` u)]/2 = 1 + 2iˇ uu` + uˇ2 .

(6.97)

As an error check, it is easy to verify that these results agree with (2.57) and (2.8), a Υ(ba)

´ 8π ³ 2 [1ˆ2] 2 [3 ˆ4] = u`u`,b c` − uˇuˇ,b cˇ + uˇcˇ δb j − i` uc` δb j 3c √ √ √ ( −g¦ ),b ( −N¦ ),b 4π −g¦ ˆa + √ + √ j N[ab] , = − √ −g¦ 3c −N¦ −N¦ 2

2

a Υ[ba] = 0.

(6.98) (6.99) (6.100)

80

The tetrad formalism allows the approximation |fˆν µ | ¿ 1 to be stated somewhat more rigorously as |` u| ¿ 1, |ˇ u| ¿ 1. From (6.28-6.33), a charged particle will have ˆ 2 , u` = 0. From (2.37) we have |ˇ uˇ ≈ Q/r u|2 ∼ 10−66 for worst-case fields accessible to measurement, so the approximation |fˆν µ | ¿ 1 is quite valid for almost all cases of interest. Let us consider the tetrad version of the approximate solution of the connection equations (2.62-2.63), which is calculated in Appendix V. This solution differs from the exact solution (6.66-6.95) only by the factors c`,ˇ c,` z ,ˇ z . From (6.96,6.97) and cˇ ≈ 1 + uˇ2 /2,

c` ≈ 1 − u`2 /2,

(6.101)

these factors will induce terms which are two orders higher in u` and uˇ than the leading order terms. This confirms that the next higher order terms in (2.62-2.63) will be two orders higher in f µ ν than the leading order terms, and from (2.37) these terms must be < 10−66 of the leading order terms. Now consider the tetrad version of the approximation (2.34,2.35) for Nνµ in terms of gνµ and fνµ . From (6.36,6.38,6.33,6.35,6.3,6.10) we have, to second order in u` and uˇ, 1 + uˇ2 /2 + u`2 /2 = 1 + uˇ2 − `/4,

(6.102)

−` c/ˇ c ≈ −1 + uˇ2 /2 + u`2 /2 = −1 + u`2 + `/4,  

(6.103)

cˇ/` c ≈

N(ab)

c 0 0   0 cˇ/`     cˇ/` 0 0    c 0  ≈ gab + fˆa c fˆcb − 1 gab `, =    4  0 0 0 −` c/ˇ c       0 0 −` c/ˇ c 0

81

(6.104)



N[ab]



−ˇ ucˇ/` c 0 0   0     uˇcˇ/` 0 0 0   c   ≈ fˆab . =     0 0 0 i` uc`/ˇ c       0 0 −i` uc`/ˇ c 0

(6.105)

These results match the order fˆ2 approximations (2.34,2.35). The next higher order terms of (6.102,6.103) will be two orders higher in u` and uˇ than the leading order terms. This confirms that the next higher order terms in (2.34,2.35) will be two orders higher in fˆµ ν than the leading order terms, and from (2.37) these terms must be < 10−66 of the leading order terms. Now let us consider the tetrad version of the charged solution (3.1-3.7). The tetrads are similar to those of the Reissner-Nordstr¨om solution[64], except for the cˇ factors, , e1 α = lα = (1/aˇ c, 1, 0, 0),

e1α = lα = (1, −1/aˇ c, 0, 0)

(6.106)

1 1 (aˇ c, 1, 0, 0) , e2 α = nα = (1, −aˇ c, 0, 0), (6.107) 2 2 p 1 = −r cˇ/2 (0, 0, 1, i sin θ), e3 α = mα = √ (0, 0, 1, i csc θ), (6.108) r 2ˇ c

e2α = nα = e3α = mα

where “a” is defined with (3.4) and from (6.1,6.36-6.39,3.5) we have u` = 0 , s` = 0 , c` = 1, ˆ Q sˇ = 2, cˇ cˇ r ˆ Q sˇ = 2 , r

(6.109) (6.110)

uˇ =

cˇ = √

(6.111) s

√ 1 = 1 + sˇ2 = 1 − uˇ2 82

1+

ˆ2 Q r4

(6.112)

From (6.2,6.106-6.108,6.1,6.40,6.41,6.7), the tetrad solution matches the solution (3.1,3.2) derived previously, W σµ = eσ a W ab eb µ 



1 1  0 1+ uˇ 0 0  aˇc     1− uˇ 0 0 0  1 − aˇc  2  2  = eσ a      0 0 0 −1 0  0      0 0 −1 0 0 0   1

1

 aˇc 2    1 − aˇc  2 =   0 0    0 0 

0 0 √1 r 2ˇ c i csc √ θ r 2ˇ c

0

(1+ˇ u)  2

  (1−ˇu) 0   aˇc    0 √1  r 2ˇ c    √ θ − ircsc 0 2ˇ c

(6.113)

 0 0 √1 r 2ˇ c √1 r 2ˇ c

0

   0      i csc √ θ r 2ˇ c 

(6.114)

√ θ − ircsc 2ˇ c



− (1+ˇu2)aˇc

0

(1−ˇ u)

0

0

− r√12ˇc

0

√ θ − r√12ˇc − ircsc 2ˇ c



s 0 0 1/a −ˇ       sˇ −aˇ  2 c 0 0   1 , =  cˇ   0  2 0 −1/r 0       0 0 0 −1/r2 sin2 θ

0

   0    (6.115)   i csc √ θ r 2ˇ c 

(6.116)

1/e = det(ea ν ) = i csc θ/ˇ cr 2 ,

(6.117)

e = det(ea ν ) = −iˇ cr2 sin θ,

(6.118)

p √ −N = −N¦ e = ie/ˇ cc` = r2 sin θ,

(6.119)



−g =



−g¦ e = ie = cˇr2 sin θ.

(6.120)

Let us calculate the spin coefficients and Wely tensor components for our charged

83

solution, so that it may classified. The nonzero tetrad derivatives are, µ

e32,1

¶0

p (aˇ c)0 , e33,2 = −r cˇ/2 i cos θ , 2 0 2 p rˇ c −ˇ c + sˇ2 −1 = − cˇ/2 − √ = √ =√ , e33,1 = e32,1 i sin θ . 2 2ˇ c 2ˇ c cˇ 2ˇ c cˇ

e11,1 = −

1 aˇ c

,

e20,1 =

(6.121) (6.122)

From these and (6.16), the λabc coefficients are λa1b = e11,1 (ea 1 eb 1 − ea 1 eb 1 ) = 0, λ221 = e20,1 (e2 0 e1 1 − e2 1 e1 0 ) =

(6.123)

(aˇ c)0 , 2

(6.124)

λ123 = e20,1 (e1 0 e3 1 − e1 1 e3 0 ) = 0,

(6.125)

λ223 = e20,1 (e2 0 e3 1 − e2 1 e3 0 ) = 0,

(6.126)

λ324 = e20,1 (e3 0 e4 1 − e3 1 e4 0 ) = 0,

(6.127)

λ132 = e30,1 (e1 0 e2 1 − e1 1 e2 0 ) = 0,

(6.128)

λ233 = −e32,1 e2 1 e3 2 − e33,1 e2 1 e3 3 = 0,

(6.129) µ

1

2

1

3

λ243 = −e42,1 e2 e3 − e43,1 e2 e3 = −2 λ441 = e42,1 e4 2 e1 1 + e43,1 e4 3 e1 1 = 0,

λ334

¶µ

−aˇ c 2



1 a √ =− , 2rˇ c r 2ˇ c

(6.130) (6.131)

µ

λ431

−1 √ 2ˇ c cˇ

¶ −1 1 1 √ = − 2, = e32,1 e4 e1 + e33,1 e4 e1 = 2 √ (6.132) rˇ c 2ˇ c cˇ r 2ˇ cµ ¶ p i csc θ 1 cot θ 3 2 2 3 √ √ = √ . (6.133) = e33,2 (e3 e4 − e3 e4 ) = 2(−r cˇ/2 i cos θ) r 2ˇ c r 2ˇ c r 2ˇ c 2

1

3

1

84

From (6.14), the spin coefficients are similar to those of the Reissner-Nordstr¨om solution[64], except for the cˇ factors, 1 , rˇ c2 a γ243 = λ243 = − , 2rˇ c 1 1 cot θ (γ213 +γ343 ) = λ334 = √ , 2 2 2r 2ˇ c 1 1 −cot θ (γ214 +γ344 ) = λ344 = √ , 2 2 2r 2ˇ c 0 1 1 (aˇ c) (γ212 +γ342 ) = λ221 = , 2 2 4 1 (γ211 +γ341 ) = 0, 2

ρ = γ314 = λ431 = −

(6.134)

µ =

(6.135)

β = α = γ = ² =

(6.136) (6.137) (6.138) (6.139)

τ = γ312 = 0,

(6.140)

π = γ241 = 0,

(6.141)

κ = γ311 = 0,

(6.142)

σ = γ313 = 0,

(6.143)

λ = γ244 = 0,

(6.144)

ν = γ242 = 0.

(6.145)

The Weyl tensor components calculated with MAPLE are Ψ2

1 = − cˇ

Ã

ˆ2 2Q 1+ 4 r



ˆ 2 Vˆ m Λe Q Λe Λe cˇ − + − 3 4 r 4r 6 6

Ψ0 = Ψ1 = Ψ3 = Ψ4 = 0.

! +

ˆ2 ˆ2 Λe Q Q + , 6r4 2ˇ cr6

(6.146) (6.147)

The results κ = σ = λ = ν = ² = 0 and Ψ0 = Ψ1 = Ψ3 = Ψ4 = 0 prove that the charged solution (3.1,3.7) has the classification of Petrov type-D, the same as the Reissner-Nordstr¨om solution. 85

Next let us find the exact solution for Nσµ in terms of gσµ and fσµ . Using (6.10,6.35) we have 

 2

fˆac fˆc b

Nab

 0 uˇ   uˇ2 0  =   0 0    0 0 





0 0 u3 0 0   0 −ˇ       uˇ3 0 0 0 0 0     , fˆac fˆc d fˆd b =  ,     0 0 0 −i` 3 0 u`2  u          2 3 u` 0 0 0 i` u 0 

0 (1− uˇ)/` c2 0 0        (1+ uˇ)/` c2 0 0 0   ,  = cˇc`    0 0 0 −(1−i` u)/ˇ c2        0 0 −(1+i` u)/ˇ c2 0 

(6.148)

(6.149)



0 (1− uˇ)(1 + u`2 ) 0 0       (1+ uˇ)(1 + u`2 )  0 0 0   . (6.150) = cˇc`     0 0 0 −(1−i` u)(1 − uˇ2 )       0 0 −(1+i` u)(1 − uˇ2 ) 0 Using (6.36,6.38,6.31,6.33) gives 1 1 1 cˇc` = p =p =q . (6.151) (1 − uˇ2 )(1 + u`2 ) 1 − (ˇ u2 − u`2 ) − uˇ2 u`2 ˆ 1 − `/2 + f /g The exact solution for Nµν is then N(µν) = cˇc`((1 − `/2)gµν + fˆµρ fˆρ ν ) ,

(6.152)

N[µν] = cˇc`((1 − `/2)fˆµν + fˆµρ fˆρ α fˆα ν ) ,

(6.153)

Nµν = cˇc`((1 − `/2)δµα + fˆµρ fˆρα )(gαν + fˆαν ) .

86

(6.154)

These can be written so that they are approximately correct for any dimension, N(µν) =

(1 − `/(n−2))gµν + fˆµρ fˆρ ν q , ˆ 1 − `/(n−2) + f /g

(6.155)

N[µν] =

(1 − `/(n−2))fˆµν + fˆµρ fˆρ α fˆα ν q , ˆ 1 − `/(n−2) + f /g

(6.156)

(1 − `/(n−2))δµα + fˆµρ fˆρα q (gαν + fˆαν ) , ˆ 1 − `/(n−2) + f /g

(6.157)

Nµν =

Nαα = q

n − 2`/(n−2)

,

(6.158)

gµν (1−n/2) + fˆµρ fˆρ ν 1 N(µν) − g(µν) Nαα = q . 2 ˆ 1 − `/(n−2) + f /g

(6.159)

1 − `/(n−2) + fˆ/g

Now let us consider the null field case where fˆσ µ fˆµ σ = det(fˆµ ν ) = 0. Using (6.11) we have



N aab

p

−N¦





0 0 0 1     1 0 u´ u´    , =    0 −´ u 0 −1       0 −´ u −1 0 p = −det(Nab ) = i,

−g¦ =

p

−det(gab ) = i.

 u2 2´    1  Nba =    −´  u   −´ u

 1



u´    0 0 0  ,  0 0 −1    0 −1 0

(6.160)

(6.161) (6.162)

In terms of ordinary Newman-Penrose formalism, we are representing null fields with the three complex Maxwell scalars set to φ0 = fˆ13 = u´, φ1 = 0, φ2 = 0. With a type III tetrad transformation we have φ0 → φ0 eiθ/A, φ1 → φ1 , φ2 → φ2 e−iθA, for arbitrary real functions θ and A. Therefore by performing a type III transformation we may always choose u´ to be a real constant representing the magnitude of the field. This ˜ α , and when ˆj ν = 0 is sometimes helpful because it reduces the number of terms in Γ σµ 87

it makes Ampere’s law just a relationship between spin coefficients. From (6.12,6.13), Ampere’s law (2.47) for null fields is 4π ˆc b ˆac c ˆba j = fˆbc ,b + γab f +γab f , c 4π ˆ2 1 ˆ42 2 ˆ42 3 ˆ42 4 ˆ42 2 ˆ42 2 ˆ24 j = fˆ42 ,4 + γ41 f + γ42 f + γ43 f + γ44 f + γ24 f + γ42 f c

(6.163)

1 ˆ32 2 ˆ32 3 ˆ32 4 ˆ32 2 ˆ32 2 ˆ23 + fˆ32 ,3 + γ31 f + γ32 f + γ33 f + γ34 f + γ23 f + γ32 f

(6.164)

= u´,4 + u´,3 + (γ241 − γ344 + γ124 + γ231 − γ433 + γ123 )´ u,

(6.165)

= δ ∗ u´ + δ´ u + 2Re(π − 2α)´ u,

(6.166)

4π ˆ1 1 ˆ42 1 ˆ24 1 ˆ32 1 ˆ23 j = γ24 f + γ42 f + γ23 f + γ32 f c

(6.167)

= (−γ242 − γ232 )´ u

(6.168)

= −2Re(ν)´ u,

(6.169)

4π ˆ4 4 ˆ24 4 ˆ42 4 ˆ24 3 ˆ24 2 ˆ24 1 ˆ24 f f + γ42 f + γ24 f + γ24 f + γ23 f + γ22 j = fˆ24 ,2 + γ21 c 4 ˆ23 4 ˆ32 f f + γ32 + γ23

(6.170)

= −´ u,2 + (−γ122 + γ423 + γ342 − γ323 )´ u

(6.171)

= −∆´ u + (2γ − µ + λ∗ )´ u,

(6.172)

The connection equations (2.55) can be solved exactly for null fields. Using ˆ4 ˜ ˇ (2) (D.9,D.7,D.11) and letting Uαστ = Υ αστ from (D.12), the order f solution for Γανµ is given by, ˜ ανµ = Γανµ + Υ ˇ [αµ]τ fˆτ ν + Υ ˇ ανµ ,(6.173) ˇ [αν]τ fˆτ µ + Υ ˇ (νµ)τ fˆτ α − 1 Υσ gµν +Υσ gν)α + Υ Γ σ(µ 2 σα where Υσσα =

2 ˇ Υστ α fˆτ σ , (n−2)

(6.174) 88

ˇ ανµ = Uανµ + Uαστ fˆσ µ fˆτ ν + U(µσ)τ fˆσ α fˆτ ν − U(νσ)τ fˆσ α fˆτ µ + U[νµ]σ fˆσ τ fˆτ α Υ 1 2 Uστ α fˆτ σ fˆµν + Uσρτ fˆρσ fˆτ [µ gν]α , (n−2) (n−2) 1 ˆ 8π ˆ = (fνµ;α + fˆαµ;ν − fˆαν;µ ) + j[ν gµ]α . 2 (n−1) +

Uανµ

(6.175) (6.176)

It happens that for the special case of null fields (6.173-6.176) is exact instead of approximate. This can be proven by substituting (6.175,6.176) into (D.11) and using (6.3,6.10) with constant u´. The only properties of Uamn needed to prove this are Uamn = −Uanm and U223 = U224 = 0, and these are easy to see from its definition, 1 ˆ 8π ˆ j[n gm]a . (fnm;a + fˆam;n − fˆan;m ) + 2 (n−1) 1³ ˆ = fnm,a − γ b na fˆbm − γ b ma fˆnb 2

(6.177)

Uanm =

+ fˆam,n − γ b an fˆbm − γ b mn fˆab ´ −fˆan,m + γ b am fˆbn + γ b nm fˆab +

8π ˆ j[n gm]a . (n−1)

(6.178)

It is unclear whether we can ever have ˆj σ = 6 0 for null fields, but the solution of the connection equations works even for this case. Finally, from (6.160,6.3,6.12) we see that the solution (2.34,2.35) for Nσµ in terms of gσµ and fσµ is exact instead of approximate for null fields,   

fˆac fˆc b

u 2´    0  =    0    0

2

0 0 0   0 0 0  ,  0 0 0    0 0 0

 0   2´  u2 c a fˆ c fˆ b =    0    0

1 N(ba) = gba + fˆbc fˆc a − gab fˆa c fˆc a , 4

89

 0 0 0   0 0 0  ,  0 0 0    0 0 0

N[ba] = fˆba .

(6.179)

(6.180)

6.2

Newman-Penrose asymptotically flat O(1/r2) expansion of the field equations

Here we solve the LRES field equations to O(1/r2 ) in a Newman-Penrose tetrad frame, assuming an asymptotically flat 1/r expansion of the unknowns, and assuming a retarded time coordinate which remains constant on a surface moving along with any radial propagation of radiation. We consider two main cases. For propagation at the speed-of-light with k = ω we show that LRES theory and Einstein-Maxwell theory are the same. We demonstrate radiation in the form of electromagnetic and gravitational waves, and peeling behavior of the Weyl scalars, and we show that the Proca equation (2.81) has the trivial solution θν = 0 corresponding to Faraday’s law. For propagation different than the speed-of-light with k < ω and 2Λb = ω 2 −k 2 , the Proca equation could potentially have Proca-wave solutions, and this analysis could determine whether such solutions have positive or negative energy. In fact what we find is that no Proca-wave solutions exist. This work emulates the analysis of EinsteinMaxwell theory in [54, 55] and to a lesser extent in [84]. It is all implemented in a REDUCE symbolic algebra program[63] called LRES 1OR RETARDED.TXT. In the following, Latin letters a, b...h indicate tetrad indices, and Greek letters indicate tensor indices. Let us ignore the θ, φ coordinates for the moment. Following

90

[54] we assume that in t, r coordinates the flat-space tetrads are     0 ea

ν

1/2 −1/2 , =   1 1  

 1 −1  , =   1/2 1/2   

b 0e ν

1/2 −1/2 1 0  1/2 1/2   = . e = 0 aν      1 1 0 −1 1 −1

(6.181)

(6.182)

We can check that these satisfy the requirements for Newman-Penrose tetrads      b ν 0e ν 0e a

 1 −1   1/2 1 1 0 = ,  =      −1/2 1 0 1 1/2 1/2     

 1 1/2 1/2 1/2 1 0  b   = , g = e e = 0 µν 0 µ 0 bν      −1 1/2 1 −1 0 −1      0 gab

1/2 1/2  1/2 1 0 1  = . = 0 eaν 0 eν b =       1 −1 −1/2 1 1 0

(6.183)

(6.184)

(6.185)

The calculations are done using a retarded time coordinate u = t − kr/ω

(6.186)

where k = (wavenumber), ω = (f requency), r = (radius). The transformation from t, r coordinates to u, r coordinates has the transformation matrix     ∂u

 ∂t T =   

∂r ∂t ∂t

T

−1

 ∂u =   ∂r ∂u

∂u ∂r 

1 −k/ω  = ,    ∂r 0 1 ∂r   

∂t ∂r 

1 k/ω  . =    ∂r 0 1 ∂r

91

(6.187)

(6.188)

Transforming the flat-space tetrads and metric to u, r coordinates gives      k/ω − 1   1 −1  1 k/ω   1 b     ,  e = = 0 ν      1/2 1/2 0 1 1/2 k/2ω + 1/2     0 ea

ν

1/2 −1/2  1  =    1 1 −k/ω   

(6.189) 

0 k/2ω + 1/2 −1/2 = , (6.190)    1 1 − k/ω 1    

 1 0 1 0  1 k/ω   1 0 1 k/ω    =   (6.191)  g = 0 νµ        k/ω 1 0 −1 0 1 k/ω 1 0 −1   k/ω   1 , =    k/ω k 2 /ω 2 − 1   0g

νµ

(6.192)

1 − k 2 /ω 2 k/ω  . =    k/ω −1

(6.193)

We assume that in cartesian coordinates the 1st approximation beyond flat-space has a 1/r falloff, the 2nd approximation has a 1/r2 falloff, etcetera. Considering that 0 gµν = diag(1, −1, −r2 , −r2 sin2 θ) and 0 g µν = diag(1, −1, −1/r2 , −1/r2 sin2 θ) in spherical coordinates, we can conclude that in spherical coordinates a 1/r falloff should look like (1/r, 1/r, 1, 1) for a covariant vector and (1/r, 1/r, 1/r2 , 1/r2 ) for a contravariant vector. Following [54] for the θ, φ part of the flat-space tetrads, and using the results above, the covariant tetrads are assumed to be of the form eb ν = 0 eb ν + 1 eb ν + 2 eb ν

92

(6.194)

where 

b 0e ν

b 1e ν



k/ω − 1 0 0  1      1/2 k/2ω + 1/2  0 0    , =   √ √  0  2 −ir sinθ/ 2 0 −r/      √ √  0 0 −r/ 2 ir sinθ/ 2   a0 /r    b /r  0 = ²   c /r  0   d0 /r 

a1 /r a2 a3    b1 /r b2 b3   ,  c1 /r c2 c3     d1 /r d2 d3 2

b 2e ν

 A0 /r    B /r2  0 = ²2    C /r2  0   D0 /r2

(6.195)

(6.196)



A1 /r2 A2 /r A3 /r    B1 /r2 B2 /r B3 /r   .  C1 /r2 C2 /r C3 /r     D1 /r2 D2 /r D3 /r

(6.197)

and the contravariant tetrads are assumed to be of the form ea ν = 0 ea ν + 1 ea ν + 2 ea ν

93

(6.198)

where 

0 ea

1 ea

ν

ν



0 0 k/2ω + 1/2 −1/2       1 − k/ω  1 0 0   , =    √ √  0 0 −1/(r 2) i/( 2r sinθ)        √ √ 0 0 −1/(r 2) −i/( 2r sinθ)   0 a /r    b0 /r  = ²   c0 /r    d0 /r  0

2 ea

ν

a1 /r a2 /r2 a3 /r2    1 2 2 3 2 b /r b /r b /r  ,  c1 /r c2 /r2 c3 /r2     d1 /r d2 /r2 d3 /r2

(6.199)

(6.200)

 2

 A /r   B 0 /r2  2 = ²   C 0 /r2    D0 /r2

1

A /r

2

2

A /r

3

B 1 /r2 B 2 /r3 C 1 /r2 C 2 /r3 D1 /r2 D2 /r3

3

3

A /r    3 3 B /r  .  C 3 /r3     D3 /r3

(6.201)

These tetrads are a generalization of those used in [54], reducing to the same form for speed-of-light propagation with k = ω. The functions di , di , Di , Di are complex conjugates of the functions ci , ci , C i , Ci . The ² parameter is included in the program to keep track of the order of terms, but we set ² = 1 in any final result. Using a symbolic algebra linear equation solver, the program calculates the 16 coefficients ai , bi , ci , di in terms of the coefficients ai , bi , ci , di by solving to O(²) the set of 16 equations (0 ea ν + 1 ea ν )(0 eb ν + 1 eb ν ) = δab .

(6.202)

Then it calculates the 16 coefficients Ai , B i , C i , Di in terms of the coefficients ai , bi , ci , di , 94

Ai , Bi , Ci , Di by solving to O(²2 ) the set of 16 equations ea ν eb ν = δab .

(6.203)

Since the contravariant antisymmetric field f µν satisfies Ampere’s law (2.47) exactly in both LRES theory and Einstein-Maxwell theory, we require the dual field ∗ = εµναβ f αβ /2 fµν

(6.204)

to be the curl of a dual potential ∗ = A∗ν,µ − A∗µ,ν . fµν

(6.205)

This ensures that Ampere’s law is satisfied automatically. Using the same considerations as with the tetrads regarding a 1/r falloff in spherical coordinates, the dual potential is assumed to be of the form A∗ν = 1 A∗ν + 2 A∗ν

(6.206)

where µ



∗ 1 Aν

µ ∗ 2 Aν

= ²2

= ² h0 /r, h1 /r, h2 , h3 , ¶ H0 /r2 , H1 /r2 , H2 /r, H3 /r .

(6.207) (6.208)

This dual potential is also used to make the system of equations well defined, because A∗ν contains only 4 unknowns, and our Proca equation contains only 4 equations. Note that [54] does not use either an ordinary potential or a dual potential, but instead uses the 6 components of the electromagnetic field for unknowns. It is unclear why he does this, since he obtains the same result but with more calculations. 95

Our unknowns are then the dual potential components hi , Hi and the tetrad components ai , bi , ci , Ai , Bi , Ci . Note that the unknowns di , Di are complex conjugates of ci , Ci , so c0i , c00i , Ci0 , Ci00 determine di , Di . All of the unknowns are assumed to depend only on the three coordinates u, θ, φ, and not on r. The goal is to calculate the field equations and then solve them for these unknowns. The first step in calculating the field equations is to calculate the λabc coefficients, spin coefficients and Riemann tensor, which are found from the equations λabc = ebσ,µ (ea σ ec µ − ec σ ea µ ), γabc =

1 (λabc + λcab − λbca ) = ebµ;ν ea µ ec ν , 2

Rmnpq = −γmnp,q + γmnq,p − γmnr λp r q + γmrp γ r nq − γmrq γ r np .

(6.209) (6.210) (6.211)

The calculation of λabc to O(1/r2 ) would ordinarily be very time consuming because of the O(r) components in ea ν . To speed things up we calculate λabc to O(²2 ) and then truncate the result to O(1/r2 ). We checked that this gives the same result as doing the calculation the long way. ˜ mn we use the method described in To calculate the nonsymmetric Ricci tensor R Appendix S. To do this calculation we use Υαµν from the solution (2.61-2.64) to the connection equations, and ∗ f ab = −εabcd fcd /2

(6.212)

from above. The fundamental tensor is calculated in tetrad form using the following

96

relations from (2.23,C.12,C.14), p √ −N¦ = (1 − fˆab fˆba /4) −g¦ ,

(6.213)

p √ 1/ −N¦ = (1 + fˆab fˆba /4)/ −g¦ ,

(6.214)

p √ N aab = (g ab − fˆab ) −g¦ / −N¦ , p √ Nbc = (gbc + fˆbc + fˆbd fˆd c ) −N¦ / −g¦ , where



gab = g ab

fˆab √

p −det(gab ) = i,

p p −N¦ = −det(Nab ).

(6.216)



0 0 1 0     1 0 0 0   , =    0 0 0 −1       0 0 −1 0 √ −1/2 = f ab 2 i Λb ,

−g¦ =

(6.215)

(6.217)

(6.218) (6.219) (6.220)

Then we calculate the tetrad equivalent of the source-free field equations, which consist of the symmetric part of the Einstein equations (2.31), and the Proca equation derived from (2.33), ˜ (ab) + Λb N(ab) + Λz gab = 0, R

(6.221)

˜ [ab|c] + Λb N[ab|c] ) = 0. εabcd (R

(6.222)

Great care is taken to ensure that everything is calculated to O(1/r2 ). To get practical computation time and memory usage when multiplying two expressions, it was essential to determine the min and max powers of 1/r in each expression, and then 97

to truncate them to the lowest power of 1/r required for their product to be accurate to O(1/r2 ). The calculation of the symmetric field equations (6.221) was checked by p doing the calculation another way, using the ordinary Ricci tensor Rnq = Rnpq from

˜ ab − Rab similar to G ˜ ab − Gab from (2.67). (6.211), together with an expression for R The calculation of the Proca equation (6.222) was also checked in a similar manner, using the approximate Proca equation (2.81). Now let us consider the solution of the field equations for the speed-of-light propagation case where k = ω and we do not require 2Λb = ω 2 − k 2 . The solution is implemented in the subroutine solvekeqw(). Let us call the O(1/r) and O(1/r2 ) Einstein equations 1 Eab and 2 Eab and the Proca equations 1 Pa and 2 Pa . Looking first at the O(1/r) field equations we find that 1 E12

⇒ ∂ 2 a1 /∂u2 = 0,

(6.223)

1 E13

√ − 1 E14 ⇒ ∂ 2 a3 /∂u2 = − 2sin(θ)∂ 2 c001 /∂u2 ,

(6.224)

1 E13

√ + 1 E14 ⇒ ∂ 2 a2 /∂u2 = − 2∂ 2 c01 /∂u2 ,

(6.225)

1 E11 1 P1

⇒ ∂ 2 c003 /∂u2 = −sin(θ)∂ 2 c02 /∂u2 ,

(6.226)

⇒ ∂h1 /∂u = 0.

(6.227)

This solves the field equations to O(1/r). The O(1/r2 ) equations impose similar requirements as the O(1/r) equations, but

98

are more restrictive, 2 E22

⇒ ∂a1 /∂u = 0,

(6.228)

2 E23

√ − 2 E24 ⇒ ∂a3 /∂u = − 2sin(θ)∂c001 /∂u,

(6.229)

2 E23

√ + 2 E24 ⇒ ∂a2 /∂u = − 2∂c01 /∂u.

(6.230)

Applying these requirements, E33 then requires that either a1 = 0 or ∂ 2 c02 /∂u2 = 0 and ∂ 2 c03 /∂u2 = −sin(θ)∂ 2 c002 /∂u2 .

(6.231)

Following [54] we will concentrate on the case a1 = 0. Then we find that 2 E34

⇒ ∂c003 /∂u = −sin(θ)∂c02 /∂u.

(6.232)

The remaining field equations do not put any constraints on O(²) parameters (ai , bi , ci , hi ), but they can instead be solved to get complicated expressions for O(²2 ) parameters (Ai , Bi , Ci , Hi ) in terms of O(²) parameters 2 P2

⇒ ∂ 2 H10 /∂u2 ,

(6.233)

2 E12

⇒ ∂ 2 A1 /∂u2 ,

(6.234)

2 E11

⇒ ∂ 2 C300 /∂u2 ,

(6.235)

2 E13

− 2 E14 ⇒ ∂ 2 C100 /∂u2 ,

(6.236)

2 E13

+ 2 E14 ⇒ ∂ 2 C10 /∂u2 .

(6.237)

Substituting these expressions solves all of the field equations to O(1/r2 ). Note that the only requirement on the dual potential is that h1 is a constant, and we are free to set h1 = 0 since hν is a potential. The remaining components h0 , h2 , h3 99

are all arbitrary functions of u, θ, φ, which is to be expected since plane waves in flat space can have any shape and angular pattern. Also, we find that Fab = fab to O(1/r2 ). The Proca field and the electric and magnetic fields are, θρ = 0,

(6.238)

1 ∂h3 1 ∂h2 , Eφ = − , rsin θ ∂u r ∂u 1 ∂h2 1 ∂h3 = 0, Bθ = , Bφ = . r ∂u rsin θ ∂u

Er = 0, Eθ =

(6.239)

Br

(6.240)

We define the “effective” energy-momentum tensor as 8πTab = Gab where Gab is the Einstein tensor formed from the symmetric metric. With T µν = ea µ T ab eb ν we set P0 = T 00 , Pr = T 01 , Pθ = T 02 r, Pφ = T 03 rsin θ, where the factors r and rsin θ account for basis vector scaling. The resulting energy and power densities are P0

1 = Pr = 2 4πr sin2 θ



∂h2 ∂u

¶2

µ sin2 θ +

∂h3 ∂u

¶2 # , Pθ = Pφ = 0. (6.241)

The Ψ0 Weyl tensor component is O(1/r), indicating the presence of gravitational radiation, 1 1 Ψ0 = √ 2 rsin θ

µ

∂ 2 c00 ∂ 2 c0 ∂ 2 c0 2sin θ 22 − i sin θ 22 − i 23 ∂u ∂u ∂u

¶ .

(6.242)

The functions c02 , c002 , c03 are arbitrary functions of u, θ, φ. The Weyl tensor component Ψ1 is O(1/r2 ), and Ψ2 , Ψ3 , Ψ4 are of higher order in 1/r, which indicates the start of peeling behavior. Our calculations were only done to O(1/r2 ), so we could not verify the peeling behavior beyond this order. Note that this peeling behavior is opposite to the usual behavior because we have made our tetrads consistent with [54]. As shown in [64], the tetrad transformation e1 ν ↔ e2 ν has no effect on the metric but causes 100

the exchanges Ψ4 ↔ Ψ0 , Ψ3 ↔ Ψ1 , and this transformation would make our results conform to the usual convention. Finally, it happens that our O(1/r2 ) solution for the tetrads and electromagnetic field solves the O(1/r2 ) Einstein-Maxwell field equations. Therefore, from the standpoint of a Newman-Penrose 1/r expansion of the field equations, LRES theory is identical to Einstein-Maxwell theory to O(1/r2 ) for speed-of-light propagation. Now let us consider the solution of the field equations for the case k < ω where we require 2Λb = ω 2 −k 2 , and 2Λb = mass2 of possible Proca radiation. The solution is implemented in the subroutine solvekltw(). Again we call the O(1/r) and O(1/r2 ) Einstein equations 1 Eab and 2 Eab and the Proca equations 1 Pa and 2 Pa . We will start with the O(1/r) equations. In the following, the requirements actually involve 2nd derivatives with respect to u but we are integrating them once with zero constant of integration. This has the effect of excluding possible tetrad solutions involving linear functions of u, which is justified by the bad behavior of such solutions as t → ∞.

1 E33

1 E33

+ 1 E44 +21 E34 ⇒ ∂c02 /∂u = 0,

(6.243)

1 E33

+ 1 E44 −21 E34 ⇒ ∂c003 /∂u = 0,

(6.244)

− 1 E44 ⇒ ∂c03 /∂u = −sin θ∂c002 /∂u, 1 E12

⇒ a1 = (2b1 − 4¯ a)(ω − k)/(ω + k) for some a ¯(θ, φ),

(6.245) (6.246)

1 E13

√ + 1 E14 ⇒ ∂c01 /∂u = [2(w−k)∂b2 /∂u−(w+k)∂a2 /∂u]/(2 2 ω),

1 E13

√ − 1 E14 ⇒ ∂c001 /∂u = [2(w−k)∂b3 /∂u−(w+k)∂a3 /∂u]/(2 2 ωsinθ). (6.248)

101

(6.247)

In the following I am not integrating with respect to u because hν is a potential, 21 P2 − 1 P1 or 21 P2 + 1 P1 ⇒ ∂ 4 h1 /∂u4 + ω 2 ∂ 2 h1 /∂u2 = 0, ¯ 1 sin(ωu + h ˇ 1) + h ˜1 + h ˆ 1 u, ⇒ h1 = h 1 P3

+ 1 P4 ⇒ ∂ 4 h2 /∂u4 + ω 2 ∂ 2 h2 /∂u2 = 0, ¯ 2 sin(ωu + h ˇ 2) + h ˜2 + h ˆ 2 u, ⇒ h2 = h

1 P3

− 1 P4 ⇒ ∂ 4 h3 /∂u4 + ω 2 ∂ 2 h3 /∂u2 = 0, ¯ 3 sin(ωu + h ˇ 3) + h ˜3 + h ˆ 3 u. ⇒ h3 = h

(6.249) (6.250) (6.251) (6.252) (6.253) (6.254)

ˇ 1, h ˇ 2, h ˇ 3 are constants, and h ¯ 1, h ˜ 1, h ˆ 1 ,h ¯ 2, h ˜ 2, h ˆ 2 ,h ¯ 3, h ˜ 3, h ˆ 3 are dependent only on Here h ∗ θ, φ and not on “u”. But with the terms linear in u, some components of fµν become à ! µ ¶ ¯1 ˜1 ˆ1 ∂ h ∂ h ∂ h 1 ∂H1 1 ∗ ˇ sin(ωu + h1 ) + +u − 2 + H2 , (6.255) f23 = − r ∂θ ∂θ ∂θ r ∂θ à ! µ ¶ ¯1 ˜1 ˆ1 1 ∂ h ∂ h ∂ h 1 ∂H1 ∗ ˇ f24 = − sin(ωu + h1 ) + +u − 2 + H3 , (6.256) r ∂φ ∂φ ∂φ r ∂φ ∗ f34 =

˜ ˆ ¯3 ∂h ˇ 3 ) + ∂ h3 + u ∂ h3 sin(ωu + h ∂θ ∂θ ∂θ ¯ ˜ ˆ 2 1 µ ∂H3 ∂H2 ¶ ∂ h2 ∂ h ∂ h 2 ˇ 2) − − sin(ωu + h −u + − . ∂φ ∂φ ∂φ r ∂θ ∂φ

(6.257)

To get good behavior as t → ∞ we have ∗ ∗ ˆ1 = h ` 1 = constant, f23 and f24 ⇒ h

(6.258)

∗ ˆ 2 = ω∂s(θ, φ)/∂θ + h ´ 2 (θ), f34 ⇒ h

ˆ 3 = ω∂s(θ, φ)/∂φ + h ´ 3 (φ), (6.259) h

´ 2 (θ), h ´ 3 (φ), h ` 1 . We may also let h ˇ3 = 0 where we are using the new variables s(θ, φ), h without any loss of generality. This solves the field equations to O(1/r). For the O(1/r2 ) equations if we let 1 s(θ, φ) = ω

µ

Z f (θ, φ) −

Z ´ 2 (θ)dθ − h 102

¶ ´h3 (φ)dφ ,

(6.260)

then the following combination of the Proca equations gives µ

¶ ωr2 sin2 θ 22 P2 0 = + ω − k ω + k c(ω 2 − k 2 ) µ ¶ ∂ 2 f (θ, φ) ∂f (θ, φ) ∂ ˘ 1 ω sin2 θ, = + sinθ sinθ +h ∂φ2 ∂θ ∂θ 2 P1

(6.261) (6.262)

˘1 = h ` 1 ω/(ω 2 − k 2 ). This is solved by assuming where h f (θ, φ) = v(θ, φ) + v¯(θ),

(6.263)

˘ 1 cos θ + h ´ 1 )/sin θ, ∂¯ v (θ)/∂θ = (ω h

(6.264)

´ 1 is another constant. Then v(θ, φ) must satisfy the generalized Legendre where h equation for l = m = 0 ∂ 2 v(θ, φ) ∂ + sinθ 2 ∂φ ∂θ

µ

∂v(θ, φ) sinθ ∂θ

¶ = 0,

(6.265)

which has the unique solution v(θ, φ) = (Y00 spherical harmonic) = constant.

(6.266)

Therefore we have ∂f (θ, φ)/∂φ = 0,

(6.267)

˘ 1 cos θ + h ´ 1 )/sin θ. ∂f (θ, φ)/∂θ = (ω h

(6.268)

∗ Now let us look at another component of fµν ,

∗ f13

˘ ´ ¯ 2 ωcos(ωu + h ˇ 2 ) + (h1 ω cos θ + h1 ) + 1 = h sin θ r

µ

∂H2 ∂h0 − ∂u ∂θ

¶ −

1 ∂H0 . (6.269) r2 ∂θ

∗ ˘ 1 = 0 and h ´ 1 = 0. finite for θ = 0 or π (the z-axis) we must have h To make f13

Applying these results and forming a combination of the Einstein equations gives, µ 0=−

22 E11 2 E34 + ω+k ω−k



¯ 2 ω 2 sin2 θ h 4r2 sin2 θ 1 ¯ 2 sin2 θ + h ¯ 2. = +h 2 3 c2 (ω + k)(ω 2 − k 2 ) ω2 − k2 103

(6.270)

¯1 = h ¯2 = h ¯ 3 = 0. The sum of positive numbers must be positive so this requires h Therefore h1 , h2 , h3 have no wavelike component but are just functions of θ and φ. From the above result we can tell that there are not going to be any Proca wave solutions. However let us continue to solve the equations and see what we get. Next we look at some additional combinations of the Proca equation 2 P1 2 P3

+ 2 P4

2 P3

− 2 P4

2 ∂ 4 H1 ∂ 3 h0 ∂h0 2 ∂ H1 + ω + + ω2 = 0, 4 2 3 ∂u ∂u ∂u ∂u 2 2 ∂ 4 H2 ∂ 4 h0 2 ∂ H2 2 ∂ h0 ⇒ + ω − − ω = 0, ∂u4 ∂u2 ∂θ∂u3 ∂θ∂u 2 2 ∂ 4 h0 ∂ 4 H3 2 ∂ H3 2 ∂ h0 ⇒ + ω − − ω = 0. ∂u4 ∂u2 ∂φ∂u3 ∂φ∂u



(6.271) (6.272) (6.273)

These have the general solution, ¯ 0 sin(ωu + h ˇ 0) + h ˜ 0, h0 = −∂H1 /∂u + h

(6.274)

¯ 2 sin(ωu + H ˇ 2) + H ˜2 + H ˆ 2 u, H2 = −∂H1 /∂θ + H

(6.275)

¯ 3 sin(ωu + H ˇ 3) + H ˜3 + H ˆ 3 u. H3 = −∂H1 /∂φ + H

(6.276)

ˇ 0, H ¯ 0, h ˜ 0, h ˆ 0, H ˇ 2, H ˇ 3 are constants, and h ¯ 2, H ˜ 2, H ˆ 2, H ¯ 3, H ˜ 3, H ˆ 3 are functions of Here h ∗ only θ, φ, and not “u”. Looking again at fµν gives

˜1 1 ∂h 1 ¯ ˇ ˜ ˆ − 2 (H 2 sin(ωu + H2 ) + H2 + H2 u), r ∂θ r ˜1 1 ∂h 1 ¯ ˇ ˜ ˆ =− − 2 (H 3 sin(ωu + H3 ) + H3 + H3 u). r ∂φ r

∗ f23 =−

(6.277)

∗ f24

(6.278)

ˆ2 = H ˆ 3 = 0. To get good asymptotic behavior as t → ∞ requires H The remaining field equations do not put any constraints on O(²) parameters (ai , bi , ci , hi ), but they can instead be solved to get complicated expressions for O(²2 ) parameters (Ai , Bi , Ci , Hi ) in terms of O(²) parameters. Substituting these expressions solves all of the field equations to O(1/r2 ). 104

Again we define the “effective” energy-momentum tensor as 8πTab = Gab where Gab is the Einstein tensor formed from the symmetric metric. With T µν = ea µ T ab eb ν we set P0 = T 00 , Pr = T 01 , Pθ = T 02 r, Pφ = T 03 rsin θ, where the factors r and rsin θ account for basis vector scaling. The Proca field, electric and magnetic fields, power densities, and Weyl scalars all indicate no radiation, θν = 0

to O(1/r),

(6.279)

Er = Eθ = Eφ = Br = Bθ = Bφ = 0

to O(1/r),

(6.280)

P0 = Pr = Pθ = Pφ = 0

to O(1/r2 ),

(6.281)

Ψ0 = Ψ1 = Ψ2 = Ψ3 = Ψ4 = 0

to O(1/r2 ).

(6.282)

The fact that all of the Weyl scalars vanish indicates that there is no gravitational radiation, which is to be expected because we are requiring propagation at a speed different than the speed-of-light. The lack of any 1/r component of θρ or 1/r2 component of Pr indicates that there is no propagating Proca radiation. The Proca field ¯ 0 cos(ωu + h ˇ 0 )/r2 components, but this does not correspond to propaθa does have h gating radiation and seems of little interest. Also, many of the functions contained in these terms would probably be determined if we were to solve the field equations to a higher order, and it is likely that this would cause these higher order terms to vanish. So the final result is that LRES theory does not have Proca-wave solutions, given the assumed form of the solution. However, there is some uncertainty as to whether this analysis really rules out Proca-wave solutions, because we may have put too strong of a constraint on the form of the solution. For a wavepacket type of solu105

tion, the wavepacket should be expected to spread out as a function of radius, and it seems unlikely that this behavior could be represented by a simple 1/r expansion. Also, for a continuous-wave type of solution, the analysis assumes a constant speed of propagation, whereas one might expect the speed of propagation to slow down due to the energy of the wave at smaller radii. There is another issue regarding the propagation speed of Proca waves. The Proca equation is 2Λb θρ = − θρα; α using a (1,−1,−1,−1) signature. The continuous-wave solution in flat space goes as sin(ωt − kx), where k < ω and 2Λb = ω 2 − k 2 . From a quantum mechanical viewpoint we have mass2 = h ¯ 2 2Λb = (¯hω)2 − (¯hk)2 = energy 2 − momentum2 , so the √ √ “particle” velocity would be v ≈ momentum/mass = h ¯ k/(¯h 2Λb ) = k/ ω 2 −k 2 < 1, which is below the speed of light. This is consistent with the group velocity which √ is vgroup = dω/dk = k/ 2Λb +k 2 < 1. However with k < ω the phase velocity is vphase = ω/k =



2Λb +k 2 /k > 1, so the wavefront (and our retarded coordinate) is

travelling at greater than the speed of light! This does not seem right. It makes one wonder if we are not finding a continuous-wave Proca wave solution simply because they are somehow inconsistent with even ordinary general relativity.

106

Chapter 7 Extension of the Einstein-Schr¨ odinger theory for non-Abelian fields

7.1

The Lagrangian density

Here we generalize LRES theory to non-Abelian fields. The resulting theory incorporates the U(1) and SU(2) gauge terms of the Weinberg-Salam Lagrangian, and when the rest of the Weinberg-Salam Lagrangian is included in a matter term, we get a close approximation to ordinary Einstein-Weinberg-Salam theory. Einstein-WeinbergSalam theory can be derived from a Palatini Lagrangian density, 1 √ −g [ g µν Rνµ (Γ) + 2Λb ] 16π 1 √ + −g tr(Fρα g αµ g ρνFνµ ) + Lm (gµν , Aν , ψ, φ · · · ), 32π

L(Γλρτ , gρτ , Aν ) = −

107

(7.1)

where the electro-weak field tensor is defined as Fνµ = 2A[µ,ν] +

ie [Aν , Aµ ]. 2¯hsinθw

(7.2)

The Hermitian vector potential Aσ can be decomposed into a real U(1) gauge vector Aσ , and the three real SU(2) gauge vectors biν , Aν = IAν + σi biν , where the σi are the Pauli spin matrices,    



(7.3)







1 0 0 1 0 −i 1 0   , σ1 =   , σ2 =   , σ3 =  , I=         0 1 1 0 i 0 0 −1

(7.4)

[σi , σj ] = 2i²ijk σk , σi† = σi , tr(σi ) = 0, tr(σi σj ) = 2δji .

(7.5)

The Lm term couples the metric gµν and vector potential Aµ to a spin-1/2 wavefunction ψ, scalar function φ, and perhaps the additional fields of the Standard Model. Here and throughout this paper we use geometrized units with c = G = 1, the symbols

( )

and

[ ]

around indices indicate symmetrization and antisymmetrization, and

[A, B] = AB−BA. The constant θw is the weak mixing angle and Λb is a bare cosmological constant. The factor of 1/2 in (7.2) results because we are including Aν and biν in one gauge term tr(Fρα g αµ g ρνFνµ ), and because we are using σi instead of the usual τi = σi /2. bλ in The original Einstein-Schr¨odinger theory allows a nonsymmetric Nµν and Γ ρτ place of the symmetric gµν and Γλρτ , and excludes the tr(Fρα g αµ g ρνFνµ ) term. Our “non-Abelian Λ-renormalized Einstein-Schr¨odinger theory” introduces an additional bρ and Nνµ to have d×d matrix cosmological term g1/2d Λz as in (2.2), and also allows Γ νµ 108

components, 1 ˆ νµ ) + d(n−2)Λb ] N 1/2d [ tr(N aµν R 16π 1 1/2d − g d(n−2)Λz + Lm (gµν , Aν , ψ, φ · · · ), 16π

bα , Nρτ ) = − L(Γ ρτ

(7.6)

where Λb ≈ −Λz so that the total Λ matches astronomical measurements[48] Λ = Λb +Λz ≈ 10−56 cm−2 ,

(7.7)

and the vector potential is defined to be p σ b[νσ] Aν = Γ /[(n −1) −2Λb ].

(7.8)

The Lm term is not to include a tr(Fρα g αµ g ρνFνµ ) term but may contain the rest of the Weinberg-Salam theory. Matrix indices are assumed to have size d = 2, and tensor indices are assumed to have dimension n=4, but we will retain “d” and “n” in the equations to show how easily the theory can be generalized. The non-Abelian Ricci tensor is ˆ νµ = R

bανµ,α Γ



α b(α(ν),µ) Γ

bτ Γ bρ Γ [τ ν] [ρµ] 1 bα bσ 1 bσ bα α σ bσµ − bνα Γ + Γνµ Γ(σα) + Γ(σα) Γνµ − Γ . 2 2 (n−1)

(7.9)

For Abelian fields the third and fourth terms are the same, and this tensor reduces to the Abelian version (2.5). This tensor reduces to the ordinary Ricci tensor for bα = 0 and Γ bα Γ [νµ] α[ν,µ] = 0, as occurs in ordinary general relativity. Let us define the symmetric tensor gµν by g1/2d gµν = N 1/2d N a(µν) .

(7.10)

Note that (7.10) defines gµν unambiguously because g = [det(g1/2d gµν )]2/(n−2) . The “physical” metric is denoted with a different symbol gµν , and in this paper we will 109

just be assuming the special case gµν = Igµν . The symmetric metric is used for measuring space-time intervals, covariant derivatives, and for raising and lowering indices. If we did not assume gµν = Igµν , we would need to choose between several metric definitions which all reduce to the definition (2.4) for Abelian fields, √

−ggµν = tr(g1/2d gµν )/d

or gµν = tr(gµν )/d

or gµν = tr(gµν )/d,

(7.11)

and we would also need to choose between g1/2d Λz

or



−g Λz

(7.12)

in the Lagrangian density (7.6). These definitions are all the same with the assumption gµν = Igµν , so we will not choose between them here. The determinants g = det(gνµ ) and N = det(Nνµ ) are defined as usual but where Nνµ and gνµ are taken to be nd × nd matrices. The inverse of Nνµ is defined to be N aµkνi = (1/N )∂N/∂Nνiµk where i,k are matrix indices, or N aµν = (1/N )∂N/∂Nνµ using matrix notation. The field N aµν satisfies the relation N aµkνiNνiσj = δσµ δjk , or N aµνNνσ = δσµ I using matrix notation. Likewise gνσ is the inverse of gµν such that ¯ατ = T ν Nνµ T µ for some coordinate transformation T ν = gµνgνσ = δσµ I. Assuming N α τ α ¯ = det(N ¯ατ ) will contain d times as many T ν ∂xν /∂ x¯α , the transformed determinant N α factors as it would if Nατ had no matrix components, so N and g are scalar densities √ √ of weight 2d. The factors N 1/2d and g1/2d are used in (7.6) instead of −N and −g to make the Lagrangian density a scalar density of weight 1 as required. Note that with an even d, we do not want the factor of −1. For our theory the electro-weak field tensor f νµ is defined by √ 1/2 g1/2d f νµ = iN 1/2d N a[νµ] Λb / 2. 110

(7.13)

√ −1/2 Then from (7.10), gµν and f µν 2 iΛb are parts of a total field, √ −1/2 (N/g)1/2d N aνµ = gµν +f µν 2 iΛb . We will see that the field equations require fνµ ≈ 2A[µ,ν] +

(7.14) √

−2Λb [Aν , Aµ ] to a very

high precision. From (7.2,2.3) we see that this agrees with Einstein-Weinberg-Salam theory when µ ¶2 1 e α = 2 = 1.457 × 1063 cm−2 , −Λz ≈ Λb = 2 2 2¯hsinθw 8lP sin θw where lP =

(7.15)

p

G¯ h/c3 = 1.616×10−33 cm, α = e2 /¯ hc = 1/137 and sin2 θw = .2397.

ρ α bνµ ˜νµ It is helpful to decompose Γ into a new connection Γ , and Aν from (7.8),

bα = Γ ˜ α + (δ α Aν − δ α Aµ ) Γ νµ νµ µ ν

p

−2Λb ,

α bσ bσ ˜α = Γ bα + (δ α Γ where Γ νµ νµ µ [σν] − δν Γ[σµ] )/(n−1).

(7.16) (7.17)

˜ ανµ has the symmetry By contracting (7.17) on the right and left we see that Γ ˜α = Γ bα = Γ ˜α , Γ να (να) αν

(7.18)

so it has only n3 −n independent components. Substituting the decomposition (7.16) into (7.9) gives from (R.16), b = Rνµ (Γ) ˜ + 2A[ν,µ] Rνµ (Γ)

p −2Λb + 2Λb [Aν , Aµ ]

p ˜ α ] − [A(ν , Γ ˜ α ]) −2Λb . + ([Aα , Γ νµ µ)α

(7.19)

˜ ανµ and Aσ Using (7.19), the Lagrangian density (7.6) can be rewritten in terms of Γ

111

from (7.17,7.8), L = −

h p 1 ˜ νµ + 2A[ν,µ] −2Λb + 2Λb [Aν , Aµ ] N 1/2d tr(N aµν (R 16π i p ˜ α ] − [A(ν , Γ ˜ α ]) −2Λb )) + d(n−2)Λb + ([Aα , Γ νµ µ)α



1 1/2d g d(n−2)Λz + Lm (gµν , Aσ , ψe , φ . . . ). 16π

(7.20)

˜ νµ = Rνµ (Γ), ˜ and from (7.18) our non-Abelian Ricci tensor (7.9) reduces to Here R 1 ˜σ ˜α 1 ˜α ˜σ ˜α ˜σ ˜α ˜ νµ = Γ ˜α − Γ R α(ν,µ) + Γνµ Γσα + Γσα Γνµ − Γνα Γσµ . νµ,α 2 2

(7.21)

˜ ανµ and Aν fully parameterize Γ bανµ and can be treated as indepenFrom (7.16,7.18), Γ dent variables. The fields N 1/2d N a(νµ) and N 1/2d N a[νµ] (or gνµ and f νµ ) fully parameterize Nνµ and can also be treated as independent variables. It is simpler to calculate ˜ ανµ = 0, δL/δAν = 0, δL/δ(N 1/2d N a(µν) ) = 0 and the field equations by setting δL/δ Γ bα = 0 and δL/δNνµ = 0, so we will δL/δ(N 1/2d N a[µν] ) = 0 instead of setting δL/δ Γ νµ follow this method.

7.2

Invariance properties of the Lagrangian density

Here we show that the Lagrangian density is real (invariant under complex conjugation), and is also invariant under U(1) and SU(2) gauge transformations. The Abelian Lambda-renormalized Einstein-Schr¨odinger theory comes in two versions, one where bρ and Nνµ are real, and one where they are Hermitian. The non-Abelian theory Γ νµ bρνµ and Nνµ are real, and one where they have also comes in two versions, one where Γ 112

bα ∗ = Γ bα and N ∗ = Nµkνi , where i, k are matrix nd×nd Hermitian symmetry, Γ νiµk µkνi νiµk indices. Using matrix notation these symmetries become bα ∗ = Γ bα T , Γ νµ µν

˜α ∗ = Γ ˜α T , Γ νµ µν

∗ T Nνµ = Nµν ,

N aµν ∗ = N aνµ T ,

(7.22)

where “T” indicates matrix transpose (not transpose over tensor indices). We will assume this Hermitian case because it results from Λz < 0, Λb > 0 as in (2.12). From (7.22,7.10,7.13,7.8) the physical fields are all composed of d×d Hermitian matrices, ∗ T bα ∗ = Γ bα T , A∗ = AT . = fνµ , Γ gνµ ∗ = gνµ T, g∗νµ = gTνµ , f νµ ∗ = f νµ T, fνµ ν (νµ) (νµ) ν

(7.23)

Hermitian fνµ and Aν are just what we need to approximate Einstein-WeinbergSalam theory. And of course gνµ and gνµ will be Hermitian if we assume the special case where they are multiples of the identity matrix. Writing the symmetries as ∗ Nνiµk = Nµkνi , g∗νiµk = gνkµi = gµkνi , and using the result that the determinant of a

Hermitian matrix is real, we see that the nd× nd matrix determinants are real N ∗ = N,

g∗ = g,

g ∗ = g.

(7.24)

Also, using (7.22) and the identity M1T M2T = (M2 M1 )T we can deduce a remarkable property of our non-Abelian Ricci tensor (7.9), which is that it has the same nd×nd bανµ and Nνµ , Hermitian symmetry as Γ ˆ∗ = R ˆT . R νµ µν

(7.25)

From the properties (7.25,7.22,7.24) and the identities tr(M1 M2 ) = tr(M2 M1 ), tr(M T ) = tr(M ) we see that our Lagrangian density (7.6) or (7.20) is real.

113

With an SU(2) gauge transformation we assume a transformation matrix U that is special (det(U ) = 1) and unitary (U † U = I). Taking into account (7.3,7.8,7.16), we assume that under an SU(2) gauge transformation the fields transform as follows, 1 Bν → U Bν U −1 − √ U,ν U −1 , −2Λb 1 Aν → U Aν U −1 − √ U,ν U −1 , −2Λb

(7.26) (7.27)

Aν → Aν ,

(7.28)

bα U −1 + 2δ α U,µ] U −1 , bα → U Γ Γ [ν νµ νµ

(7.29)

−1 bα bα Γ (νµ) → U Γ(νµ) U ,

(7.30)

−1 bα bα Γ + (n−1) U,µ U −1 , [αµ] → U Γ[αµ] U

(7.31)

˜ ανµ → U Γ ˜ ανµ U −1 , Γ

(7.32)

Nνµ → U Nνµ U −1 ,

gνµ → U gνµ U −1 ,

fνµ → U fνµ U −1 ,

(7.33)

N aµν → U N aµν U −1 ,

gµν → U gµν U −1 ,

f µν → U f µν U −1 .

(7.34)

Under a U(1) gauge transformation all of the fields are unchanged except 1 Aν → Aν + √ ϕ,ν , 2Λb I A ν → Aν + √ ϕ,ν , 2Λb bα → Γ bα − 2iI δ α ϕ,µ] , Γ νµ νµ [ν bα bα Γ [αµ] → Γ[αµ] − (n−1)iI ϕ,µ .

114

(7.35) (7.36) (7.37) (7.38)

Writing the SU(2) gauge transformation (7.33) as   

0 Nνµ

U 0 0 0 N00      0 U 0 0 N   10  =    0 0 U 0 N   20     0 0 0 U N30



N01 N02 N03 U −1 0 0 0       0 U −1 0 0  N11 N12 N13          −1 N21 N22 N23  0 0 U 0      N31 N32 N33 0 0 0 U −1

(7.39)

and using the identity det(M1 M2 ) = det(M1 )det(M2 ), we see that the nd× nd matrix determinants are invariant under an SU(2) gauge transformation, N → N,

g → g,

g → g.

(7.40)

Another remarkable property of our non-Abelian Ricci tensor (7.9) is that it transforms the same as Nνµ under an SU(2) gauge transformation (7.29), as in (R.11), α bαρτ U −1 +2δ[ρ bαρτ )U −1 Rνµ (U Γ U,τ ] U −1 ) = U Rνµ (Γ

for any matrix U (xσ ).

(7.41)

The results (7.40,7.41) actually apply for a general matrix U , and do not require that det(U ) = 1 or U † U = I. Using the special case U = Ie−iϕ in (7.41) we see that our non-Abelian Ricci tensor (7.9) is also invariant under a U(1) gauge transformation, bα − 2iI δ α ϕ,τ ] ) = Rνµ (Γ bα ) Rνµ (Γ ρτ [ρ ρτ

for any ϕ(xσ ).

(7.42)

From (7.41,7.33,7.40,7.42) and the identity tr(M1 M2 ) = tr(M2 M1 ) we see that our Lagrangian density (7.6) or (7.20) is invariant under both U(1) and SU(2) gauge transformations, thus satisfying an important requirement to approximate EinsteinWeinberg-Salam theory. One of the motivations for this theory is that the Λz = 0, Lm = 0 version can be derived from a purely affine Lagrangian density as well as a Palatini Lagrangian 115

density, the same as with the Abelian theory in Appendix M. The purely affine Lagrangian density is bα ) = [ det(Nνµ )]1/2d , L(Γ ρτ

(7.43)

where Nνµ is simply defined to be ˆ νµ /Λb . Nνµ = −R

(7.44)

bαρτ = 0 gives the Considering that N aµν = (1/N )∂N/∂Nνµ , we see that setting δL/δ Γ same result obtained from (7.6) with Λz = 0, Lm = 0, ˆ νµ /δ Γ bα ] = 0. tr[N aµν δ R ρτ

(7.45)

bαρτ , there are no δL/δ(N 1/2d N aµν ) = 0 field equations. Since (7.43) depends only on Γ However, the definition (7.44) exactly matches the δL/δ(N 1/2d N aµν ) = 0 field equations obtained from (7.6) with Λz = 0, Lm = 0. Note that there are other definitions of N and g which would make the Lagrangian density (7.6) real and gauge invariant, for example we could have defined N = tr(det(Nνµ )) or N = Det(det(Nνµ )), where det() is done only over the tensor indices. However, with these definitions the field N aµν = (1/N )∂N/∂Nνµ would not be a matrix inverse such that N aσνNνµ = δµσ I. Calculations would be very unwieldy in a theory where N aµν = (1/N )∂N/∂Nνµ appeared in the field equations but was not a genuine inverse of Nνµ . In addition, it would be impossible to derive the Λz = 0, Lm = 0 version of the theory from a purely affine Lagrangian density, thus removing a motivation for the theory. Note that we also cannot use the definition N = det(tr(Nνµ )) as in [19] because (−det(tr(Nνµ )))1/2 and ˆ νµ )))1/2 would not depend on the traceless part of the fields. (−det(tr(R 116

7.3

The field equations

Let us calculate the field equations for the following special case, ˜ ανµ = tr(Γ ˜ ανµ )I/d, Γ

gνµ = tr(gνµ )I/d.

(7.46)

In this case Aν and N 1/2d N a[νµ] are the only independent variables in (7.20) which are not just multiples of the identity matrix I. This assumption is both coordinate independent and gauge independent, considering (7.32,7.34). We assume this special case because it gives us Einstein-Weinberg-Salam theory, and because it greatly simplifies ˜ νµ = tr(R ˜ νµ )I/d, and the the theory. With the assumption (7.46) we also have R √ ˜ α ] − [A(ν , Γ ˜ α ]) −2Λb vanishes in the Lagrangian density (7.20). And term ([Aα , Γ νµ µ)α as mentioned initially, with the assumption (7.46) several metric definitions such as (7.11) are the same, so we need not choose one or the other. It would be interesting to investigate the more general theory described by the Lagrangian density (7.6,7.20) without the restriction (7.46). However, it is important to emphasize that any solution of the restricted theory will also be a solution of any of the more general theories which use one of the metric definitions (7.11). Setting δL/δAτ = 0 and using the definition (7.13) of f νµ gives the ordinary Weinberg-Salam equivalent of Ampere’s law, (g1/2d f ωτ ), ω −

p

−2Λb g1/2d [f ωτ, Aω ] = 4πg1/2d j τ ,

(7.47)

where the source current j τ is defined by jτ =

−1 δLm . g1/2d δAτ 117

(7.48)

˜ β = 0 using a Lagrange multiplier term tr[Ωρ Γ ˜ α ] to enforce the Setting δL/δ Γ τρ [αρ] symmetry (7.18), and using the result tr[(g1/2d f ωτ ), ω ] = 4πg1/2d tr[j τ ] derived from (7.47,7.3,7.5) gives the connection equations, ˜ τσβ N 1/2d N aρσ + Γ ˜ ρ N 1/2d N aστ − Γ ˜ αβα N 1/2d N aρτ ] tr[(N 1/2d N aρτ ), β + Γ βσ √ 8π 2 i τ] = g1/2d tr[j [ρ ]δβ . 1/2 (n−1)Λb

(7.49)

Setting δL/δ(N 1/2d N a(µν) ) = 0 using the identities N = [det(N 1/2d N aµν )]2/(n−2) and g = [det(N 1/2d N a(µν) )]2/(n−2) gives our equivalent of the Einstein equations, ˜ (νµ) + Λb N(νµ) + Λz gνµ ] = 8πtr[Sνµ ], tr[R

(7.50)

where Sνµ is defined by Sνµ ≡ 2

δLm 1/2d δ(N N (µν) )

=2

δLm . δ(g1/2d gµν )

(7.51)

Setting δL/δ(N 1/2d N a[µν] ) = 0 using the identities N = [det(N 1/2d N aµν )]2/(n−2) and g = [det(N 1/2d N a(µν) )]2/(n−2) gives, ˜ [νµ] + 2A[ν,µ] R

p −2Λb + 2Λb [Aν , Aµ ] + Λb N[νµ] = 0.

(7.52)

Note that the antisymmetric field equations (7.52) lack a source term because Lm in (7.20) contains only g1/2d gµν = N 1/2d N a(νµ) from (7.10), and not N 1/2d N a[νµ] . The trace operations in (7.49,7.50) occur because we are assuming the special case (7.46). The ˜ β and δL/δ(N 1/2d N a(µν) ) vanish because off-diagonal matrix components of δL/δ Γ τρ with (7.46), the Lagrangian density contains no off-diagonal matrix components of ˜ β and N 1/2d N a(µν) . The trace operation sums up the contributions from the diagonal Γ τρ ˜ βτρ and N 1/2d N a(µν) because (7.46) means that for a given set matrix components of Γ of tensor indices, all of the diagonal matrix components are really the same variable. 118

To put (7.47-7.52) into a form which looks more like the ordinary EinsteinWeinberg-Salam field equations we need to do some preliminary calculations. The definitions (7.10,7.13) of gνµ and fνµ can be inverted to give Nνµ in terms of gνµ and fνµ . An expansion in powers of Λ−1 b is derived in Appendix C, N(νµ)

µ = gνµ − 2 f σ (ν fµ)σ −

¶ 1 tr(f ρ σf σ ρ ) −1 −3/2 gνµ Λb + (f 3 )Λb . . . 2(n−2) d

√ −1/2 N[νµ] = fνµ 2 iΛb + (f 2 )Λ−1 b .... −3/2

Here (f 3 )Λb

(7.53) (7.54)

ˆρ ˆσ ˆ −3/2 and f σ [ν fµ]σ Λ−1 . and (f 2 )Λ−1 b b refer to terms like f σ f (µ fν)ρ Λb

Because of the assumption (7.46) and the trace operation in (7.49), the connection equations (7.49) are the same as with the Abelian theory (2.55) but with the substitution of tr[fνµ ]/d and tr[j ν ]/d instead of fνµ and j ν . Therefore the solution of the connection equations from (2.38) can again be abbreviated as ˜ α(νµ) = IΓανµ + (f 0 f )Λ−1 . . . Γ b

˜ α[νµ] = (f 0 )Λ−1 . . . , Γ b

(7.55)

where Γανµ is the Christoffel connection, Γανµ =

1 ασ g (gµσ,ν + gσν,µ − gνµ,σ ). 2

(7.56)

Substituting (7.55) using (R.4) shows that as in (2.39), the Non-symmetric Ricci tensor (7.21) can again be abbreviated as ˜ (νµ) = IRνµ + (f 0 f 0 )Λ−1 + (f f 00 )Λ−1 . . . , R b b

˜ [νµ] = (f 00 )Λ−1/2 . . . , R b

(7.57)

00 −1 and where Rνµ = Rνµ (Γ) is the ordinary Ricci tensor. Here (f 0 f 0 )Λ−1 b , (f f )Λb −1/2

(f 00 )Λb

−1/2

−1 ατ α indicate terms like tr(f σ ν;α )tr(f α µ;σ )Λ−1 b , tr(f )tr(fτ (ν; µ);α )Λb and tr(f[νµ,α]; )Λb

119

.

Combining (7.53,7.57,2.3) with the symmetric field equations (7.50) and their contraction gives µ ¶ tr(f σ(ν fµ)σ ) 1 tr(Tνµ ) tr(f ρσfσρ ) = 8π +2 − gνµ d d 4 d ³n ´ −1/2 00 −1 +Λ − 1 gνµ + (f 3 )Λb + (f 0 f 0 )Λ−1 b + (f f )Λb . . . , 2

Gνµ

(7.58)

where the Einstein tensor and energy-momentum tensor are 1 Gνµ = Rνµ − gνµ Rαα , 2 −1/2

Here (f 3 )Λb

Tνµ = Sνµ −

1 gνµ Sαα . 2

(7.59) −1/2

00 −1 ρ σ , (f 0 f 0 )Λ−1 b and (f f )Λb indicate terms like f σ f (µ fν)ρ Λb

, tr(f σ ν;α )tr(f α µ;σ )Λ−1 b

and tr(f ατ )tr(fτ (ν; µ);α )Λ−1 b . This shows that the Einstein equations (7.58) match those of Einstein-Weinberg-Salam theory except for extra terms which will be very small relative to the leading order terms because of the large value Λb ∼ 1063 cm−2 from (2.12). Combining (7.54,7.57) with the antisymmetric field equations (7.52) gives fνµ = 2A[µ,ν] + −1/2

Here (f 2 )Λb

p

−1/2

−2Λb [Aν , Aµ ] + (f 2 )Λb

+ (f 00 )Λ−1 b ....

(7.60) −1/2

and (f 00 )Λ−1 indicate terms like f σ [ν fµ]σ Λ−1 and tr(f[νµ,α]; α )Λb b b

.

From (2.12) we see that the fνµ in Ampere’s law (7.47) matches the electro-weak tensor (7.2) except for extra terms which will be very small relative to the leading order terms because of the large value Λb ∼ 1063 cm−2 from (2.12). Finally, let us do a quantitative comparison of our non-Abelian LRES theory to Einstein-Weinberg-Salam theory. If Λz is due to zero-point fluctuations we would usually expect Λb ∼ ωc4 lP2 ∼ 1066 cm−2 with cutoff frequency ωc ∼ 1/lP as in (2.12,2.13). Our Λb from (7.15) is consistent with this interpretation with a cutoff frequency 120

ωc ∼ α1/4/lP , which is just as reasonable as ωc ∼ 1/lP as far as anyone knows. For the Abelian LRES theory with Λb ∼ 1066 cm−2 , we showed in §2.4 that the higher order terms in the Einstein-Maxwell field equations were < 10−16 of the ordinary terms for worst-case field strengths and rates of change accessible to measurement. Therefore, for non-Abelian LRES theory with Λb = 1.457×1063 cm−2 from (7.15), the higher order terms in the field equations will be < 10−13 of the ordinary terms for worst-case field strengths and rates of change accessible to measurement. This is far below the level that could be detected by experiment. One aspect of this theory which might differ from Einstein-Weinberg-Salam theory is the possible existence of Proca waves, as discussed at the end of §2.4 for the purely electromagnetic case. The only change for the non-Abelian case is that Λb is fixed, so we cannot use the argument that the potential ghost goes away in the limit as ωc → ∞, Λb → ∞. If Proca-waves really do exist in the theory, it is possible that they could be interpreted as a built-in Pauli-Villars field as discussed in §2.4 and Appendix K. Finally, we should mention again that this theory would differ from Einstein-Weinberg-Salam theory if we do not assume the special case (7.46) where ˜ α are restricted to be multiples of the identity matrix. Further work is gνµ and Γ νµ necessary to compare this more general theory to experiment for reasonable choices of the metric definition (7.11). Some preliminary work on this topic can be found in Appendix X.

121

Chapter 8 Conclusions The Einstein-Schr¨odinger theory is modified to include a cosmological constant Λz which multiplies the symmetric metric. This cosmological constant is assumed to be nearly cancelled by Schr¨odinger’s “bare” cosmological constant Λb which multiplies the nonsymmetric fundamental tensor, such that the total “physical” cosmological constant Λ = Λb + Λz matches measurement. The resulting Λ-renormalized Einstein-Schr¨odinger theory closely approximates ordinary Einstein-Maxwell theory when |Λz | ∼ 1/(Planck length)2 , and it becomes exactly Einstein-Maxwell theory in the limit as |Λz | → ∞. In a similar manner, when the theory is generalized to nonAbelian fields, a special case closely approximates Einstein-Weinberg-Salam theory.

122

Appendix A A divergence identity Here we derive (4.4) using only the definitions (2.4,2.22) of gνµ and fνµ , and the identity (2.56), µ ¶ √ 3 1 µ ρ −1/2 (µ N ν) − δν Nρ ; µ − f σρ N[σρ,ν] 2 iΛb 2 2 √ 1 3 −1/2 = g σρ (N(ρν);σ + N(νσ);ρ − N(ρσ);ν ) − f σρ N[σρ;ν] 2 iΛb 2√ 2 ¤ 1 −N £ a(σρ) = √ N (N(ρν);σ + N(νσ);ρ − N(ρσ);ν ) − 3N a[ρσ] N[σρ;ν] 2 −g √ ¤ 1 −N £ aσρ = √ N (N(ρν);σ + N(νσ);ρ − N(ρσ);ν ) + 3N aσρ N[ρν;σ] 2 −g √ 1 −N aσρ = √ N (Nρν;σ + Nνσ;ρ − Nρσ;ν ) 2 −g √ ¤ 1 −N £ aσρ = √ N (Nρν;σ + Nνσ;ρ ) − N aσρ (Nρσ,ν − Γαρν Nασ − Γασν Nρα ) 2 −g √ 1 −N aσρ 1 √ = − √ (N ;σ Nρν + N aσρ ;ρ Nνσ ) − √ ( −N );ν 2 −g −g ¶ µ√ ¶ ¸ ·µ √ 1 −N aσρ −N aσρ √ N √ N =− ;σ Nρν + ;ρ Nνσ 2 −g −g i √ √ 1 h ρσ −1/2 −1/2 ρσ ρσ ρσ = − (g + f 2 iΛb );σ Nρν + (g + f 2 iΛb );ρ Nνσ 2 √ −1/2 = f σρ ;σ N[ρν] 2 iΛb .

123

(A.1) (A.2) (A.3) (A.4) (A.5) (A.6) (A.7) (A.8) (A.9) (A.10)

Appendix B Variational derivatives for fields ˜σ = 0 with the symmetry Γ [µσ] The field equations associated with a field with symmetry properties must have the same number of independent components as the field. For a field with the symmetry ˜ σ = 0, the field equations can be found by introducing a Lagrange multiplier Ωµ , Γ [µσ] Z 0=δ

˜ σ )dn x. (L + Ωµ Γ [µσ]

(B.1)

Minimizing the integral with respect to Ωµ shows that the symmetry is enforced. Using the definition, ∆L ∂L = − β ˜τρ ˜ βτρ ∆Γ ∂Γ

µ

∂L ˜ βτρ,ω ∂Γ

¶ ,ω

... ,

(B.2)

˜ β gives and minimizing the integral with respect to Γ τρ 0=

∆L ∆L 1 τ ρ δσ] = + Ωµ δβσ δ[µ + (Ωτ δβρ − δβτ Ωρ ). β β ˜ ˜ ∆Γτ ρ ∆Γτ ρ 2

124

(B.3)

Contracting this on the left and right gives Ωρ =

2 ∆L 2 ∆L =− . α ˜ ˜α (n−1) ∆Γαρ (n−1) ∆Γ ρα

(B.4)

Substituting (B.4) back into (B.3) gives δβρ δβτ ∆L ∆L ∆L 0 = − − . β ˜α ˜α ˜ τ ρ (n−1) ∆Γ (n−1) ∆Γ ∆Γ αρ τα

(B.5)

In (B.4,B.5) the index contractions occur after the derivatives. Contracting (B.5) on the right and left gives the same result, so it has the same number of independent ˜ αµν . This is a general expression for the field equations associated components as Γ ˜ σ = 0. with a field having the symmetry Γ [µσ]

125

Appendix C Approximate solution for Nνµ in terms of gνµ and fνµ Here we invert the definitions (7.11,7.13) of gνµ and fνµ to obtain (7.53,7.54), the approximation of Nνµ in terms of gνµ and fνµ , and we also do the same for Abelian fields as in (2.4,2.22) and (2.34,2.35). First let us define the notation √ −1/2 fˆνµ = f νµ 2 i Λb .

(C.1)

We assume that |fˆν µ | ¿ 1 for all components of the unitless field fˆν µ , and find a solution in the form of a power series expansion in fˆν µ . We will first consider the problem for non-Abelian fields. For the following calculations we will treat the fields as nd×nd matrices but we will only show the tensor indices explicitly. Lowering an index on the right side of the equation (N/g)1/2d N aνµ = gµν+fˆµν from (7.14) we get (N/g)1/2d N aµ α = δαµ I − fˆµ α . 126

(C.2)

Using fˆα α = 0, the well known formula det(eM ) = exp (tr(M )), and the power series ln(1−x) = −x − x2 /2 − x3 /3 . . . we get[85], 1 ln(det(I − fˆ)) = tr(ln(I − fˆ)) = − tr(fˆρ σ fˆσ ρ ) + (fˆ3 ) . . . 2

(C.3)

Here the notation (fˆ3 ) refers to terms like tr(fˆτ α fˆα σ fˆσ τ ). Taking ln(det()) on both sides of (C.2) using the result (C.3) and the identities det(sM ) = snd det(M ) and det(M −1 ) = 1/det(M ) gives 1 ln(det[(N/g)1/2d N aµ α ]) = ln((N/g)n/2−1 ) = − tr(fˆρ σ fˆσ ρ ) + (fˆ3 ) . . . , 2 1 ln[(N/g)1/2d ] = − tr(fˆρ σ fˆσ ρ ) + (fˆ3 ) . . . . 2d(n−2)

(C.4) (C.5)

Taking ex on both sides of this and using ex = 1 + x + x2 /2 . . . gives (N/g)1/2d = 1−

1 tr(fˆρ σfˆσ ρ ) + (fˆ3 ) . . . . 2d(n−2)

(C.6)

Using the power series (1−x)−1 = 1 + x + x2 + x3 . . . , or multiplying by (C.2) on the right we can calculate the inverse of (C.2) to get[85] (g/N )1/2d N ν µ = δµν I + fˆν µ + fˆν σ fˆσ µ + (fˆ3 ) . . . .

(C.7)

Lowering this on the left gives, Nνµ = (N/g)1/2d (gνµ + fˆνµ + fˆνσ fˆσ µ + (fˆ3 ) . . . ).

(C.8)

Here (fˆ3 ) refers to terms like fˆνα fˆα σ fˆσ µ . Using (7.46,C.8,C.6,C.1) we get the result (7.53,7.54). Now let us redo the calculation for Abelian fields. Lowering an index on the √ √ equation ( −N / −g )N aµν = g νµ + fˆνµ from (2.4,2.22) gives √ −N √ N aµ α = δαµ − fˆµ α . −g 127

(C.9)

Let us consider the tensor fˆµ α = fˆµν gνα . Because gνα is symmetric and fˆµν is antisymmetric, it is clear that fˆα α = 0. Also because fˆνσ fˆσ µ is symmetric it is clear that fˆν σ fˆσ µ fˆµ ν = 0. In matrix language therefore tr(fˆ) = 0, tr(fˆ3 ) = 0, and in fact tr(fˆp ) = 0 for any odd p. Using the well known formula det(eM ) = exp (tr(M )) and the power series ln(1−x) = −x − x2 /2 − x3 /3 − x4 /4 . . . we then get[85], 1 ln(det(I − fˆ)) = tr(ln(I − fˆ)) = − fˆρ σ fˆσ ρ + (fˆ4 ) . . . 2

(C.10)

Here the notation (fˆ4 ) refers to terms like fˆτ α fˆα σ fˆσ ρ fˆρ τ . Taking ln(det()) on both sides of (C.9) using the result (C.10) and the identities det(sM ) = sn det(M ) and det(M −1 ) = 1/det(M ) gives µ√ ¶ µ (n/2−1) ¶ 1 1 N −N fˆρ σ fˆσ ρ + (fˆ4 ) . . . (C.11) ln √ = ln (n/2−1) = − −g (n−2) g 2(n−2) Taking ex on both sides of (C.11) and using ex = 1 + x + x2 /2 . . . gives √ 1 −N √ fˆρσfˆσρ + (fˆ4 ) . . . = 1− −g 2(n−2)

(C.12)

Using the power series (1−x)−1 = 1 + x + x2 + x3 . . . , or multiplying (C.9) term by term, we can calculate the inverse of (C.9) to get[85] √ −g ν N µ = δµν + fˆν µ + fˆν σ fˆσ µ + fˆν ρ fˆρ σ fˆσ µ + (fˆ4 ) . . . √ −N √ −N Nνµ = √ (gνµ + fˆνµ + fˆνσ fˆσ µ + fˆνρ fˆρ σ fˆσ µ + (fˆ4 ) . . . ). −g

(C.13) (C.14)

Here the notation (fˆ4 ) refers to terms like fˆνα fˆα σ fˆσ ρ fˆρ µ . Since fˆνσ fˆσ µ is symmetric and fˆνρ fˆρ σ fˆσ µ is antisymmetric, we obtain from (C.14,C.12,C.1) the final result (2.34,2.35).

128

Appendix D ˜α Approximate solution for Γ νµ in terms of gνµ and fνµ Here we derive the approximate solution (2.62,2.63) to the connection equations (2.55). First let us define the notation √ −1/2 fˆνµ = f νµ 2 i Λb ,

√ ˆj σ = j σ 2 i Λ−1/2 , b

¯ ανµ = Υα(νµ) , Υ

ˇ ανµ = Υα[νµ] . Υ

(D.1)

We assume that |fˆν µ | ¿ 1 for all components of the unitless field fˆν µ , and find a solution in the form of a power series expansion in fˆν µ . Using (2.59) and √ √ ( −N ),α 8π −g ˆσ σ ˜ √ √ Γσα = + j N[σα] (n−2)(n−1) −N −N

129

(D.2)

√ ˜ α = Γα +Υα , ( −N /√−g )N aµν = g νµ + fˆνµ from (2.61,2.4,2.22) from (2.57) and Γ νµ νµ νµ we get √ −N ˜ ν N aµτ + Γ ˜ µ N aτ ν ) 0 = √ (N aµν ,α + Γ τα ατ −g µ ¶ 8π 1 ˆτ [µ ν] aµν ˆ − j δα + j N[τ α] N (D.3) (n−1) (n−2) √ µ√ ¶ −N N aµν −N ˜ ν aµτ ˜ µ aτ ν ˜ σσα −Γσσα )N aµν ) √ = (Γτ α N + Γατ N − (Γ ,α + √ −g −g 8π ˆ[µ ν] − j δα (D.4) (n−1) = (g νµ + fˆνµ );α + Υντ α (g τ µ + fˆτ µ ) + Υµατ (g ντ + fˆντ ) − Υσσα (g νµ + fˆνµ ) −

8π ˆ[µ ν] j δα (n−1)

(D.5)

= fˆνµ ;α + Υντ α g τ µ + Υντ α fˆτ µ + Υµατ g ντ + Υµατ fˆντ − Υσσα g νµ − Υσσα fˆνµ +

4π ˆν µ ˆµ ν (j δα − j δα ). (n−1)

(D.6)

Contracting this with gνµ gives ˇ σατ fˆτ σ ⇒ Υσσα = 0 = (2 − n)Υσσα − 2Υ

2 ˇ Υστ α fˆτ σ . (n−2)

(D.7)

Lowering the indices of (D.6) and making linear combinations of its permutations gives Υανµ = Υανµ

µ fˆνµ;α +Υνµα +Υντ α fˆτ µ +Υµαν +Υµατ fˆν τ −Υσσα gνµ −Υσσα fˆνµ ¶ 4π ˆ + (jν gαµ − ˆjµ gνα ) (n−1) µ 1 ˆ − fµα;ν +Υµαν +Υµτ ν fˆτ α +Υανµ +Υαντ fˆµ τ −Υσσν gµα −Υσσν fˆµα 2 ¶ 4π ˆ + (jµ gνα − ˆjα gµν ) (n−1) µ 1 ˆ − fαν;µ +Υανµ +Υατ µ fˆτ ν +Υνµα +Υνµτ fˆα τ −Υσσµ gαν −Υσσµ fˆαν 2 ¶ 4π ˆ ˆ + (jα gµν − jν gαµ ) . (D.8) (n−1) 1 + 2

130

Cancelling out the Υανµ terms on the right-hand side, collecting terms, and separating out the symmetric and antisymmetric parts gives, ¯ ανµ = Υ ˇ [αµ]τ fˆτ ν + Υ ˇ [αν]τ fˆτ µ + Υ ˇ (νµ)τ fˆτ α − 1 Υσσα gνµ + Υσσ(ν gµ)α Υ (D.9) 2 ¯ [νµ]τ fˆτ α − 1 Υσ fˆνµ + Υσ fˆµ]α ¯ (αν)τ fˆτ µ + Υ ˇ ανµ = −Υ ¯ (αµ)τ fˆτ ν + Υ Υ σ[ν 2 σα 1 8π ˆ + (fˆνµ;α + fˆαµ;ν − fˆαν;µ ) + j[ν gµ]α . (D.10) 2 (n−1) Substituting (D.9) into (D.10) ˇ ανµ Υ

µ ¶ 1 ˇ 1 σ σ σ σ σ ˇ [αµ]σ fˆ τ + Υ ˇ (µτ )σ fˆ α − Υ gµτ + Υ gτ )α fˆτ ν = − Υ[ατ ]σ fˆ µ + Υ σ(µ 2 2 σα ¶ µ 1 ˇ 1 σ σ σ σ σ ˇ [µα]σ fˆ τ + Υ ˇ (ατ )σ fˆ µ − Υ gατ + Υ gτ )µ fˆτ ν − Υ[µτ ]σ fˆ α + Υ σ(α 2 2 σµ µ ¶ 1 ˇ 1 σ σ σ σ σ ˇ [αν]σ fˆ τ + Υ ˇ (ντ )σ fˆ α − Υσα gντ + Υσ(ν gτ )α fˆτ µ + Υ[ατ ]σ fˆ ν + Υ 2 2 µ ¶ 1 ˇ 1 σ σ σ σ σ ˆ ˆ ˆ ˇ [να]σ f τ + Υ ˇ (ατ )σ f ν − Υσν gατ + Υσ(α gτ )ν fˆτ µ Υ[ντ ]σ f α + Υ + 2 2 µ ¶ 1 σ 1 ˇ σ σ σ σ ˆ ˆ ˆ ˇ ˇ Υ[ντ ]σ f µ + Υ[νµ]σ f τ + Υ(µτ )σ f ν − Υσν gµτ + Υσ(µ gτ )ν fˆτ α + 2 2 µ ¶ 1 σ 1 ˇ σ σ σ σ ˆ ˆ ˆ ˇ ˇ Υ[µτ ]σ f ν + Υ[µν]σ f τ + Υ(ντ )σ f µ − Υσµ gντ + Υσ(ν gτ )µ fˆτ α − 2 2 1 − Υσσα fˆνµ + Υσσ[ν fˆµ]α 2 1 8π ˆ + (fˆνµ;α + fˆαµ;ν − fˆαν;µ ) + j[ν gµ]α 2 (n−1) ´ 1 ³ˇ σ σ ˆ ˆ ˆτ ˇ = − Υ f + Υ f ατ σ µ µτ σ α f ν 2 ´ 1 ³ˇ σ σ ˆ ˆ ˆτ ˇ f + Υ f + Υ ατ σ ν ντ σ α f µ 2 ´ 1 ³ˇ σ σ σ ˆ ˆτ ˆ ˆ ˇ ˇ f f + 2 Υ f − Υ + Υ µ f α ν τ τ νσ τ µσ [νµ]σ 2 1 + Υσσα fˆµν + Υσστ fˆτ [µ gν]α 2 1 8π ˆ + (fˆνµ;α + fˆαµ;ν − fˆαν;µ ) + j[ν gµ]α , 2 (n−1)

131

and using (D.7) gives, ˇ ανµ = Υ ˇ αστ fˆσ µ fˆτ ν + Υ ˇ (µσ)τ fˆσ α fˆτ ν − Υ ˇ (νσ)τ fˆσ α fˆτ µ + Υ ˇ [νµ]σ fˆσ τ fˆτ α Υ 1 ˇ 2 ˇ Υστ α fˆτ σ fˆµν + Υσρτ fˆρσ fˆτ [µ gν]α (n−2) (n−2) 1 8π ˆ + (fˆνµ;α + fˆαµ;ν − fˆαν;µ ) + j[ν gµ]α . 2 (n−1) +

(D.11)

Equation (D.11) is useful for finding exact solutions to the connection equations ˇ ανµ . because it consists of only n2 (n − 1)/2 equations in the n2 (n − 1)/2 unknowns Υ Also, from (D.11) we can immediately see that ˇ ανµ = 1 (fˆνµ;α + fˆαµ;ν − fˆαν;µ ) + 4π (ˆjν gµα − ˆjµ gνα ) + (fˆ30 ) . . . . Υ 2 (n−1)

(D.12)

Here the notation (fˆ30 ) refers to terms like fˆατ fˆτ σ fˆσ [ν;µ] . With (D.12) as a starting ˇ ανµ by recursively substituting the current point, one can calculate more accurate Υ ˇ ανµ into (D.11). Then this Υ ˇ ανµ can be substituted into (D.7,D.9) to get Υ ¯ ανµ . For Υ our purposes (D.12) will be accurate enough. Substituting (D.12) into (D.7) we get 4π ˆ (j(ν gα)µ − ˆjµ gνα ) + (fˆ30 ) . . . , (n−1) 1ˆ 4π ˆ = fνα;µ + j[ν gα]µ + (fˆ30 ) . . . , 2 (n−1) µ ¶ 2 1ˆ 4π ˆ = fτ σ;α + j[τ gσ]α fˆτ σ + (fˆ40 ) . . . (n−2) 2 (n−1) 8π −1 ˆj τ fˆτ α + (fˆ40 ) . . . . (fˆρσfˆσρ ),α + = 2(n−2) (n−1)(n−2)

ˇ (αν)µ = −fˆµ(ν;α) + Υ

(D.13)

ˇ [αν]µ Υ

(D.14)

Υσσα

132

(D.15)

Substituting these equations into (D.9) gives µ

¶ 1ˆ 2π ˆ ˆ fαµ;τ + (jα gµτ − jµ gατ ) fˆτ ν 2 (n−1) µ ¶ 1ˆ 2π ˆ ˆ + fνα;τ + (jν gατ − jα gντ ) fˆτ µ 2 (n−1) µ ¶ 2π ˆ ˆ ˆ ˆ + −fτ (µ;ν) + (jµ gντ + jν gµτ − 2jτ gµν ) fˆτ α (n−1) µ ¶ 1 −1 8π τ ˆ ρσ ˆ ˆ ˆ j fτ α gνµ − (f fσρ ),α + 2 2(n−2) (n−1)(n−2) µ ¶ −1 8π ρσ ˆ τ ˆ ˆ ˆ + (f fσρ ),(ν + j fτ (ν gµ)α + (fˆ40 ) . . . 2(n−2) (n−1)(n−2) ³ ´ 1 τ τ ˆ ρσ ˆ ρσ ˆ ˆ ˆ ˆ ˆ ˆ = f (ν fµ)α;τ + fα fτ (ν;µ) + (f fσρ ),α gνµ − 2(f fσρ ),(ν gµ)α 4(n−2) µ ¶ 2 ˆ 4π ˆτ ˆ j fατ gνµ + fτ (ν gµ)α + (fˆ40 ) . . . . (D.16) + (n−2) (n−1)

¯ ανµ = − Υ

Here the notation (fˆ40 ) refers to terms like fˆατ fˆτ σ fˆσ ρ fˆρ (ν;µ) . Raising the indices on (D.16,D.12,D.15) and using (D.1) gives the final result (2.62,2.63,2.64). 1 α ¯ α = fˆτ (ν fˆµ) α ;τ + fˆατ fˆτ (ν;µ) + Υ ( `, α gνµ − 2 `,(ν δµ) ) νµ 4(n−2) µ ¶ 4π ˆρ ˆα 2 ˆ α + j f ρ gνµ + fρ(ν δµ) + (fˆ40 ) . . . , (n−2) (n−1) α ˇ ανµ = 1 (fˆνµ; α + fˆα µ;ν − fˆα ν;µ ) + 8π ˆj[ν δµ] + (fˆ30 ) . . . . Υ 2 (n−1)

133

(D.17) (D.18)

Appendix E Derivation of the generalized contracted Bianchi identity Here we derive the generalized contracted Bianchi identity (4.3) from the connection ˜ α . Whereas [45] derived the equations (2.55), and from the symmetry (2.8) of Γ νµ identity by performing an infinitesimal coordinate transformation on an invariant integral, we will instead use a direct method similar to [3], but generalized to include charge currents. First we make the following definitions, Wτ ρ =



√ √ −g W τ ρ = −N N aρτ = −g (g τ ρ + fˆτ ρ ),

√ −1/2 fˆνµ = f νµ 2 iΛb ,

√ ˆjα = √−gj α 2 iΛ−1/2 , b

τ ˜σ ˜ τ ναµ = Γ ˜τ − Γ ˜τ + Γ ˜σ Γ ˜τ ˜σ ˜τ R νµ,α να,µ νµ σα − Γνα Γσµ + δν Γσ[α,µ] ,

˜σ ˜σ ˜α ˜α ˜σ Γ ˜α + Γ ˜ νµ = R ˜ α ναµ = Γ ˜α − Γ R νµ σα − Γνα Γσµ + Γσ[ν,µ] . να,µ νµ,α

134

(E.1) (E.2) (E.3) (E.4)

˜ νµ is our non-symmetric Ricci tensor (2.11), which has the property from Here R (2.16), ˜T ) = R ˜ µν . Rνµ (Γ

(E.5)

˜ νµ and R ˜ τ ναµ reduce to the ordinary Ricci and Riemann tensors for The tensors R symmetric fields where Γσσ[ν,µ] = Rσσµν /2 = 0. Rewriting the connection equations (2.55) in terms of the definitions above gives, ˜ ρ Wτ σ − Γ ˜ σ Wτ ρ − ˜ τ Wσρ + Γ 0 = Wτ ρ ,λ + Γ σλ σλ λσ

4π ˆρ τ ˆτ ρ (j δλ − j δλ ). (n−1)

(E.6)

Differentiating (E.6), antisymmetrizing, and substituting (E.6) for Wτ ρ ,λ gives, µ 0 =

˜ τ Wσρ + Γ ˜ ρ Wτ σ − Γ ˜ σ Wτ ρ − W ,[λ + Γ σ[λ σ[λ [λ|σ τρ

σρ τσ ˜τ ˜ρ ˜ σ Wτ ρ − = Γ +Γ −Γ σ[λ,ν] W σ[λ,ν] [λ|σ,|ν] W

¶ 4π ˆρ τ ˆτ ρ (j δ[λ − j δ[λ ) , ν] (E.7) (n−1)

4π ˆρ τ ˆτ ρ (j ,[ν δλ] − j ,[ν δλ] ) (n−1)

˜ τ Wσρ ,ν] + Γ ˜ ρ Wτ σ ,ν] − Γ ˜ σ Wτ ρ ,ν] +Γ σ[λ σ[λ [λ|σ ρ τσ τ ˜ τσ[λ,ν] Wσρ + Γ ˜ρ ˜ σσ[λ,ν] Wτ ρ − 4π (ˆjρ ,[ν δλ] = Γ −Γ − ˆjτ ,[ν δλ] ) [λ|σ,|ν] W (n−1) ¶ µ 4π ˆρ σ ˆσ ρ ρ σ αρ σα α σρ τ ˜ ˜ ˜ ˜ (j δν] − j δν] ) −Γσ[λ Γα|ν] W + Γν]α W − Γν]α W − (n−1) µ ¶ 4π ˆσ τ ˆτ σ ρ τ ασ σ τα α τσ ˜ ˜ ˜ ˜ −Γ[λ|σ Γα|ν] W + Γν]α W − Γν]α W − (j δν] − j δν] ) (n−1) ¶ µ 4π ˆρ τ ˆτ ρ ρ τ αρ τα α τρ σ ˜ ˜ ˜ ˜ (j δν] − j δν] ) . +Γσ[λ Γα|ν] W + Γν]α W − Γν]α W − (n−1)

(E.8)

(E.9)

Cancelling the terms 2B-3A, 2C-4A, 3C-4B and using (E.3) gives, 0 =

h i i 1 h σρ ˜ τ ˜ τ ˆjρ − Γ ˜ ρ ˆjτ ˜ T ) + 4π W R σνλ + Wτ σ Rρ σνλ (Γ Γ [λν] 2 (n−1) [νλ] i 4π h ˆτ ˜ τ ˆjσ − Γ ˜ σ ˆjτ )δ ρ − (ˆjρ ,[ν + Γ ˜ ρ ˆjσ − Γ ˜ σ ˆjρ )δ τ . (E.10) + (j ,[ν + Γ σ[ν σ[ν σ[ν λ] λ] [ν|σ (n−1)

135

Multiplying by 2, contracting over ρν , and using (E.5) and ˆjν,ν = 0 from (2.49) gives, ˜ τ σνλ + Wτ σ Rν σνλ (Γ ˜T ) + 0 = W R σν

8π h ˜ τ ˆν ˜ ν ˆτ i Γ j − Γ[λν] j (n−1) [νλ]

i 8π h ˆτ τ ˆσ σ ˆτ ν ν ν ˆσ σ ˆν τ ˆ ˜ ˜ ˜ ˜ + (j ,[ν + Γσ[ν j − Γσ[ν j )δλ] − (j ,[ν + Γ[ν|σ j − Γσ[ν j )δλ] (E.11) (n−1) ˜ τ σνλ + Wτ σ R ˜ λσ − 4π(n−2) (ˆjτ ,λ + Γ ˜ τ ˆjσ − Γ ˜ σ ˆjτ ). = Wσν R (E.12) σλ σλ (n−1) This is a generalization of the symmetry Rτ λ = Rλ τ of the ordinary Ricci tensor. Next we will use the generalized uncontracted Bianchi identity[3], which can be verified by direct computation, +

+

+

˜ τ σ ν α ;λ + R ˜ τ σ α λ ;ν + R ˜ τ σ λ ν ;α = 0. R +−+

+++

+−−

(E.13)

The +/− notation is from [3] and indicates that covariant derivative is being done ˜ ανµ instead of the usual Γανµ . A plus by an index means that the associated with Γ derivative index is to be placed on the right side of the connection, and a minus means that it is to be placed on the left side. Note that the identity (E.13) is true for either the ordinary Riemann tensor or for our definition (E.3). This is because the ˜σ two tensors differ by the term δντ Γ σ[α,µ] , so that the expression (E.13) would differ by ˜ ρρ [ ν ,α ];λ + Γ ˜ ρρ [α , λ ];ν + Γ ˜ ρρ [ λ , ν ];α ). But this difference vanishes because for the term δστ (Γ −+

++

−−

an arbitrary curl Y[α,λ] we have ˜ σαλ Y[ν,σ] ˜ σλν Y[σ,α] − Γ Y [ ν ,α ];λ + Y [α , λ ];ν + Y [ λ , ν ];α = Y[ν,α],λ − Γ −+

++

−−

˜ σλν Y[α,σ] ˜ σαν Y[σ,λ] − Γ + Y[α,λ],ν − Γ ˜ σ Y[σ,ν] − Γ ˜ σ Y[λ,σ] = 0. + Y[λ,ν],α − Γ αλ αν

(E.14)

A simple form of the generalized contracted Bianchi identity results if we contract 136

+−

˜ τ σνλ and (E.6) for W σ ν ;τ , (E.13) over Wσν and τα , then substitute (E.12) for Wσν R +

+

+

˜ τ σ ν τ ;λ + R ˜ τ σ τ λ ;ν + R ˜ τ σ λ ν ;τ ) 0 = Wσν (R +−+

+++

(E.15)

+−−

+

˜ σ λ ;ν − Wσν R ˜ τ σ ν λ ;τ ˜ σ ν ;λ + Wσν R = −Wσν R ++

+−

+−−



(E.16)

+

˜ σ ν ;λ + (Wσ ν R ˜ σ λ );ν − (Wσν R ˜ τ σν λ );τ = −Wσν R +−

+

−W

+− σν



˜

;ν Rσλ

+W

+− σν

˜τ ;τ R σνλ

(E.17)



˜ σ ν ;λ + (Wσ ν R ˜ σ λ );ν = −Wσν R +− + ¶ µ + + + 4π(n−2) ˆ+τ τ σ σ τσ ˜ τ ˜ ˆj − Γ ˜ ˆj ) ; τ + W Rλσ − (j , λ + Γ σλ σλ − − (n−1) − − 4π ˆν σ ˆσ ν ˜ 4π ˆν σ ˆσ ν ˜ τ − (j δν − j δν )Rσλ + (j δτ − j δτ )R σνλ (n−1) (n−1) −

(E.18)

+

˜ σ λ );ν + (W ν σ R ˜ σ ν ;λ + (Wσ ν R ˜ λ σ );ν = −Wσν R +−

+



+ + 4π(n−2) ˆ+τ ˜ στ ˆjσ − Γ ˜ σσ ˆj τ );τ (j , λ + Γ λ λ − (n−1) − − 4π ˆν ˜ σ 4π(n−2)ˆσ ˜ j Rσλ + j R σνλ + (n−1) (n−1)



(E.19)

˜ ˜α ˜ ˜ σν,λ − Γ ˜α R = −Wσν (R σλ αν − Γλν Rσα ) ˜ σλ ),ν + Γ ˜ ν Wσα R ˜ σλ − Γ ˜ α Wσν R ˜ σα − Γ ˜ α Wσν R ˜ σλ +(Wσν R να λν αν ˜ λσ ),ν + Γ ˜ ν Wασ R ˜ λσ − Γ ˜ α Wνσ R ˜ ασ − Γ ˜ α Wνσ R ˜ λσ +(Wνσ R αν νλ αν −

4π(n−2) ˆτ ˜ τσλ,τ ˆjσ + Γ ˜ τσλˆjσ,τ − Γ ˜ σσλ,τ ˆjτ − Γ ˜ σσλˆjτ,τ [j ,λ,τ + Γ (n−1) ˜ τατ (ˆjα ,λ + Γ ˜ ασλˆjσ − Γ ˜ σσλˆjα ) +Γ ˜ σσαˆjτ ) ˜ τσαˆjσ − Γ ˜ ατλ (ˆjτ ,α + Γ −Γ ˜ α (ˆjτ ,λ + Γ ˜ τ ˆjσ − Γ ˜ σ ˆjτ ) −Γ ατ σλ σλ ˜ σλ − Γ ˜ α ) − ˆjσ (Γ ˜α −Γ ˜ α )] − ˆjσ (R α[σ,λ] ασ,λ αλ,σ

(E.20)

With the ˆjσ terms of (E.20), 4C-6A,4D-8D,5A-7A,5B-7B,5C-7C all cancel, 4A and 137

4E are zero because ˆjν,ν = 0 from (2.49), and 4B,6B,6C,8C cancel the Ricci tensor term ˜ ανµ factor cancel, which are 8A,8B. With the Wτ σ terms of (E.20), all those with a Γ the terms 1C-2C,1B-3C,2B-2D,3B-3D. Doing the cancellations and using (E.1) we get √ √ √ ˜ σλ + −N N aσν R ˜ λσ ),ν − −N N aνσ R ˜ σν,λ . 0 = ( −N N aνσ R

(E.21)

Equation (E.21) is a simple generalization of the ordinary contracted Bianchi identity √ √ ˜α 2( −g Rν λ ),ν − −g g νσRσν,λ = 0, and it applies even when j τ 6= 0. Because Γ νµ has cancelled out of (E.21), the Christoffel connection Γανµ would also cancel, so a manifestly tensor relation can be obtained by replacing the ordinary derivatives with covariant derivatives done with Γανµ , √ √ √ ˜ σλ + −N N aσν R ˜ λσ );ν − −N N aνσ R ˜ σν;λ . 0 = ( −N N aνσ R

(E.22)

Rewriting the identity in terms of g ρτ and fˆρτ as defined by (E.1,E.2) gives, √ ˜ σλ + √−g (g νσ + fˆνσ )R ˜ λσ );ν − √−g (g σν + fˆσν )R ˜ σν;λ (E.23) 0 = ( −g (g σν + fˆσν )R = =



˜ (ν λ);ν − R ˜ σ ] + √−g [2(fˆνσ R ˜ [λσ] );ν + fˆνσ R ˜ [σν];λ ] −g [2R σ;λ

(E.24)



˜ (ν λ);ν − R ˜ σσ;λ ] + √−g [3fˆνσ R ˜ [σν,λ] + 2fˆνσ ;ν R ˜ [λσ] ]. −g [2R

(E.25)

√ Dividing by 2 −g gives another form of the generalized contracted Bianchi identity µ

˜ σσ ˜ (ν λ) − 1 δλν R R 2

¶ ;ν

3 ˜ [νσ,λ] + fˆνσ;ν R ˜ [σλ] . = fˆνσ R 2

From (E.2,2.40) we get the final result (4.3).

138

(E.26)

Appendix F Validation of the EIH method to post-Coulombian order Here we state the post-Coulombian equations of motion of Einstein-Maxwell theory obtained by two authors[72, 74] using the EIH method, and show that they match the equations of motion obtained from the Darwin Lagrangian[53]. For two particles the Darwin Lagrangian takes the form La =

ma va2 1 ma va4 eb ea eb + − ea + 2 [va · vb + (va · nab )(vb · nab )] . 2 2 8 c Rab 2c Rab

(F.1)

Here we are using the notation i i 2 i i i rab . = rab /Rab , Rab = vai − vbi , niab = rab = rai − rbi , vab r˙ai = vai , r˙bi = vbi , rab

139

(F.2)

From this we get the equations of motion ¶ µ ∂La ∂ ∂La 0 = − (F.3) ∂rai ∂t ∂vai µ i ¶ i i rab ea eb rab s s 3rab vai s s vbi s s s s u u = eb eb 3 + 2 − 3 va vb − 5 va rab vb rab + 3 vb rab + 3 va rab Rab 2c Rab Rab Rab Rab µ ¶ i s ma ea eb v r −ma v˙ ai − 2 (v˙ ai va2 + 2vai vas v˙ as ) − 2 v˙ bi − vbi ab 2 ab 2c 2c Rab Rab µ ¶ s s ea eb i s s i s s i s s i u u vab rab − 2 3 vab vb rab + rab v˙ b rab + rab vb vab − 3rab vb rab 2 (F.4) 2c Rab Rab · 2 ¸ i i ea eb va rab vb2 rab i s s = −mv˙ a + ea eb 3 + 2 − − va vb + 3 Rab c 2 2 Rab u s i i ¤ rs ea eb £ 3ea eb u s rab rab rab e2a e2b rab + 2 −vas vai + vas vbi ab − + . (F.5) v v 3 5 4 c Rab 2c2 b b Rab mb c2 Rab Let us first compare the notation used in the various references, i Landau/Lif shitz rai rbi rab Rab i i W allace η ζ βi r i i Gorbatenko ξ η −Ri R Bazanski ξ i η i −Ri r Anderson xiA xiB xiAB xAB i Jackson r1i r2i r12 R

ea e1 Q e1 qA q1

eb e2 q e2 qB q2

ma m1 M m1 mA m1

mb m2 m m2 mB m2

(F.6)

The Wallace[72] equations of motion (including radiation reaction term) are ¶ µ ¶ µ ¶ ·µ ∂ ∂ 1 1 1 s s s ˙s m1 η¨ + e1 e2 m = e1 e2 η˙ η˙ + η˙ ζ ∂η r 2 ∂η m r m

¸ µ ¶ ∂ 1 1 ∂ 3r ˙r ˙s + (η˙ η˙ − η˙ ζ + ζ ζ ) s − ζ ζ ∂η r 2 ∂η m η r η s µ ¶ e21 e22 1 ∂ 1 2 − + e1 (e1 ˙¨ η m + e2 ˙¨ ζ m ). (F.7) m m2 r ∂η r 3 s m

s ˙m

˙s ˙m

Using µ ¶ 1 ∂r 1 ∂ 2r βr βs 1 βm ∂ , = β , = − 3 + δsr = − s m 3 s r s ∂η r r ∂η r ∂η η r r 3 ∂ r βs βr 3βr βs βm βm = −δ − δ + − 3 δsr rm sm ∂η m η r η s r3 r3 r5 r 1 ∂ 3 r ˙ r ˙ s ζ˙ m ζ˙ s βs 3ζ˙ r ζ˙ s βr βs βm βm ζ˙ s ζ˙ s − ζ ζ = − + , 2 ∂η m η r η s r3 2r5 2r3 140

(F.8) (F.9) (F.10)

we get · µ ¶ βm 1 s s βm s ˙s m1 η¨ − e1 e2 3 = e1 e2 − η˙ η˙ + η˙ ζ r 2 r3 m

¸ 1 ∂ 3r ˙r ˙s s ˙ m βs ˙ − (η˙ η˙ − η˙ ζ + ζ ζ ) 3 − ζ ζ r 2 ∂η m η r η s e2 e2 1 βm 2 + 1 2 η m + e2 ˙¨ ζ m) (F.11) + e1 (e1 ˙¨ m2 r r 3 3 " Ã ! 1 s s ζ˙ s ζ˙ s βm s ˙s = e1 e2 − η˙ η˙ + η˙ ζ − 2 2 r3 # 3ζ˙ r ζ˙ s βr βs βm s m s ˙ m βs − (η˙ η˙ − η˙ ζ ) 3 − r 2r5 s m

+

s ˙m

e21 e22 βm 2 + e1 (e1 ˙¨ η m + e2 ˙¨ ζ m ). m2 r 4 3

(F.12)

Translating this into the Landau/Lifshitz notation we see that it agrees with (F.5), ¶ m · µ 2 m rab va vb2 rab s s ma v˙ − ea eb 3 = ea eb − + va vb − 3 Rab 2 2 Rab ¸ s r s m 3vbr vbs rab rab rab s m s m rab − (va va − va vb ) 3 − 5 Rab 2Rab e2 e2 r m 2 + a b ab + ea (ea v¨am + eb v¨bm ). 4 mb Rab 3 m

(F.13)

The Gorbatenko[74] equations of motion (including radiation reaction term) are M ξ¨k

" qQ (ξ˙l η˙ l ) (Rl ξ˙l ) (Rl ξ˙l ) ˙ (ξ˙l ξ˙l ) = − 3 Rk + qQ R − η ˙ + ξ + Rk k k k R R3 R3 R3 2R3 ¸ η¨k (Rl η¨l ) 3 (Rl η˙ l )2 (η˙ l η˙ l ) 2 − − Rk − Rk − Rk + (Q ˙¨ ξk + q ˙¨ ηk )Q. (F.14) 3 5 3 2R 2R 2 R 2R 3

The Coulombian order equations for the η k particle are the first two terms but with ξ k → η k ,M → m, Q → q, q → Q, Rk → −Rk . Using these equations we have m¨ ηk ≈

qQ Rk R3

⇒ mRl η¨l ≈

qQ R

(Rl η¨l ) qQ Rk Rk ≈ − 3 2R 2m R4 η¨k (Rl η¨l ) qQRk ⇒− − Rk ≈ − . 3 2R 2R mR4 141

⇒−

(F.15) (F.16)

Substituting this last equation into (F.14) and assuming 1/(mR) is O(λ1 ) gives M ξ¨k

" qQ (ξ˙l η˙ l ) (Rl ξ˙l ) (Rl ξ˙l ) ˙ (ξ˙l ξ˙l ) = − 3 Rk + qQ R − η ˙ + ξ + Rk k k k R R3 R3 R3 2R3 ¸ qQRk 3 (Rl η˙ l )2 (η˙ l η˙ l ) 2 − ξk + q ˙¨ ηk )Q. (F.17) − Rk − Rk + (Q ˙¨ 4 5 3 mR 2 R 2R 3

Translating this into the Landau/Lifshifz notation we see that it agrees with (F.5), ma v˙ ak

· l l l l l l va vb k va k rab va k eb ea k rab va2 k − = r + e e r + v − v − r b a ab ab b a 3 3 3 3 3 ab Rab Rab Rab Rab 2Rab ¸ k l l 2 2 eb ea rab 3 (rab vb ) k vb2 k + − rab + r + (ea v¨ak + eb v¨bk )ea . (F.18) 4 5 3 ab mb Rab 2 Rab 2Rab 3

142

Appendix G Application of point-particle post-Newtonian methods Here we apply point-particle post-Newtonian methods to LRES theory in order to calculate what the theory predicts for the Kreuzer experiment[86]. The Kreuzer experiment is an experiment which can distinguish between active gravitational mass and inertial mass. Active gravitational mass is the mass which is the source of the Newtonian gravitational potential mA /r. Inertial mass is the mass which relates the acceleration of a body to an applied force. In particular inertial mass is the mass in the Lorentz force equation muν uµ;ν = −Q(Aν,µ −Aµ,ν )uν , which is exactly the same in our theory (4.10,4.11) as in Einstein-Maxwell theory. In [86] the computations only require the lowest-order post-Newtonian version of the Lorentz force equation, mp dvpi = mp Ui∗ + Qp Ei dt

143

(G.1)

where U∗ =

X mp

(G.2)

rp

p

Ei = A0,i .

(G.3)

In [86] the computations also only require the lowest-order post-Newtonian approximation of the electromagnetic field, fi0 = Ei = −∇i ψ,

fik = 0,

(G.4)

where ψ=

X Qp p

rp

.

(G.5)

Here we will be using the notation x = (position of observer), rp = x−xp ,

rp = |rp |,

xp = (position of particle p)

rpq = xp −xq ,

rpq = |rpq |.

(G.6) (G.7)

From our electric monopole solution (3.1,3.2,3.8) we have f01 F01

µ Q = 2 1− r µ Q ≈ 2 1+ r

¶−1/2 µ ¶ 2Q2 Q Q2 ≈ 2 1+ , Λb r 4 r Λb r4 ¶ 4m 4Q2 − . Λb r3 Λb r4

(G.8) (G.9)

The extra terms in (G.8,G.9) fall off as 1/r5 or 1/r6 , and they all include a factor of 1/Λb and are < 10−66 of the Q/r2 term for worst-case radii accessible to measurement. Based upon this result and the close approximation of equations (2.47,2.48) to the ordinary Maxwell equations, we will also assume the approximation (G.4). Therefore

144

we have f[νβ;α] = 0, fνµ; ν = 0, and the 00 component of our effective electromagnetic energy momentum tensor (2.68,2.67) is µ ¶ 1 ν EM ρν ˜ 8π T00 ≈ 2 f0 fν0 − g00 f fνρ 4 µ 1 + 2f τ (0 f0) α ;τ ;α + 2f ατ fτ (0; 0);α − f ν 0;α f α 0;ν + f ν 0;α fν0; α + f ν α;0 f α ν; 0 2 ¶ 1 3 − g00 f τ β fβ α ;τ ;α − (f ρν fνρ ), α ;α g00 − g00 f[νβ;α] f [νβ ; α] Λ−1 (G.10) b 4 4 µ ¶ 1 0i i ≈ 2 f0 fi0 − f fi0 2 ¶ µ 1 0i α 0i k i k α 0 i + 2f fi(0; 0);0 −f 0;i f 0;k +f 0;α fk0; +f i;0 f 0; 0 − (f fi0 ), ;α Λ−1 b . (G.11) 2 Time derivatives result in higher order post-Newtonian terms, raising and lowering with ηµν differs from gµν only by higher order post-Newtonian terms, and covariant derivative differs from ordinary derivative only by higher order post-Newtonian terms. Therefore we have µ ¶ µ ¶ 1 2 1 2 EM 2 ˜ 8π T00 ≈ +2 E − E + (E ),k,k − Ek,i Ei,k + Ek,i Ek,i Λ−1 b . (G.12) 2 2 From (G.4), the last two terms cancel and we have EM 8π T˜00 ≈ E2 +

Using ∇2 ψ = −4π

P p

1 2 2 1 2 ∇ (E ) = |∇ψ|2 + ∇ |∇ψ|2 . 2Λb 2Λb

(G.13)

Qp δ(rp ) from (G.5) we get the identity

X 1 2 2 1 ∇ ψ = ∇ · ∇(ψ 2 ) = ∇ · (ψ∇ψ) = |∇ψ|2 − 4π Qp δ(rp )ψ. 2 2 p

145

(G.14)

Then we have X X Qq 1 2 2 1 2 + ∇ |∇ψ|2 ∇ ψ + 4π Qp δ(rp ) 2 rq 2Λb p q XX Qq 1 2 2 1 2 = ∇ ψ + 4π Qp δ(rp ) + ∇ |∇ψ|2 2 r 2Λ q b p q X X Qq 1 2 2 1 2 = ∇ ψ + 4π Qp δ(rp ) + ∇ |∇ψ|2 2 r 2Λ pq b p q6=p ! Ã 1 2 1 2 X Qp X Qq + ∇ |∇ψ|2 . = ∇2 ψ − 2 rp q6=p rpq 2Λb p

EM 8π T˜00 =

(G.16) (G.17) (G.18)

P

mp δ(rp ) from (G.2) and including the mass part of T˜00 we get à ! X X Qp X Qq 1 2 2 1 2 = 4π mp δ(rp ) + ∇ ψ − + ∇ |∇ψ|2 (G.19) 2 rp q6=p rpq 2Λb p p à ! 1 2 X Qp X Qq 1 2 = −∇2 U ∗ + ∇2 ψ − + ∇ |∇ψ|2 . (G.20) 2 r r 2Λ p pq b p q6=p

Using ∇2 U ∗ = −4π 8π T˜00

(G.15)

p

From [86], G00 = ∇2 g00 /2 to lowest post-Newtonian order, so the Einstein equations are 1 G00 = ∇2 g00 = −∇2 U ∗ + ∇2 2

Ã

1 2 X Qp X Qq ψ − 2 rp q6=p rpq p

! +

1 2 ∇ |∇ψ|2 . (G.21) 2Λb

This has the solution g00 = 1 − 2U ∗ + ψ 2 − 2

X Qp X Qq 1 + |∇ψ|2 . rp q6=p rpq Λb p

(G.22)

Using (G.2,G.5) we get g00 = 1 − 2

X mp p

rp

+

à !2 X Qp p

rp

¯ ¯2 X Qp X Qq 1 ¯¯X Qp rpi ¯¯ −2 + ¯ ¯ . (G.23) rp q6=p rpq Λb ¯ p rp3 ¯ p

The only difference between this expression and that of ordinary Einstein-Maxwell theory is the last term, and this term falls off as 1/r4 . Since the difference between gravitational mass and inertial mass in [86] depends only on terms which fall off as 1/r, these two masses are the same for our theory as for Einstein-Maxwell theory. 146

Appendix H Alternative derivation of the Lorentz force equation Here we check the results in §4.1 by deriving the Lorentz force equation (4.10) in a different way, using the field equations (2.28) and a simple form of the generalized contracted Bianchi identity (4.1). Let us make the definitions √ √ √ √ −1/2 −1/2 Wνσ = −N N aσν = −g (g νσ + fˆνσ ), fˆνµ = f νµ 2 iΛb , ˆj ν = j ν 2 iΛb . (H.1) The generalized contracted Bianchi identity (4.1) then becomes ˜ νλ + Wσν R ˜ λν );σ − Wνσ R ˜ νσ;λ 0 = (Wνσ R

(H.2)

˜ νσ,λ . ˜ λν ),σ − Wνσ R ˜ νλ + Wσν R = (Wνσ R

(H.3)

147

From (2.56,2.22,2.4) we have √ √ −N N aρτ ;ν Nτ ρ = −N (N aρτ ,ν + Γραν N aατ + Γταν N aρα )Nτ ρ =

√ √ √ −N (N aρτ ,ν Nτ ρ + 2Γααν ) = −2( −N ),ν + 2 −N Γααν (H.5)

√ = −2( −N );ν , √

−g fˆτ ρ ;ν N[τ ρ] = =

(H.4)



(H.6)

√ −g fˆτ ρ ;ν Nτ ρ = ( −N N aρτ );ν Nτ ρ

√ √ −N N aρτ ;ν Nτ ρ + n( −N );ν

√ = (n−2)( −N );ν .

(H.7) (H.8) (H.9)

Making linear combinations of the field equations (2.28) gives, 0 =

³ ´ √ √ ˜ νλ + 2A[ν,λ] 2 iΛ1/2 + Λb Nνλ + Λe gνλ − 8πSνλ −N N aσν R b

+

³ ´ √ √ ˜ λν + 2A[λ,ν] 2 iΛ1/2 + Λb Nλν + Λe gλν − 8πSλν −N N aνσ R b



³ ´ √ √ ˜ να + 2A[ν,α] 2 iΛ1/2 + Λb Nνα + Λe gνα − 8πSνα δ σ −N N aαν R λ b

(H.10)

˜ νλ + Wσν R ˜ λν − Wνα R ˜ να δ σ − 16π √−g T σ = Wνσ R λ λ √ √ √ − 4 −g(2f νσ A[ν,λ] + f αν A[ν,α] δλσ )+ Λb −N (2−n)δλσ + Λe −g (2−n)δλσ . (H.11) Using (H.2,H.9) the divergence of this equation gives the Lorentz force equation, 0 =



˜ να − 16π √−g T σ − √−gΛb fˆτ ρ ;λ N[τ ρ] −g fˆαν ;λ R λ;σ

¢ √ ¡ − 4 −g −8πj ν A[ν,λ] + 2f νσ A[ν,λ];σ + f αν ;λ A[ν,α] + f αν A[ν,α];λ =



(H.12)

√ σ ˜ [να] + 2A[ν,α] 2 iΛ1/2 + Λb N[να] )− 16π √−g Tλ;σ −g fˆαν ;λ (R b

¢ √ ¡ − 4 −g −8πj ν A[ν,λ] + 3f αν A[ν,α,λ]

(H.13)

¢ √ ¡ σ = −16π −g Tλ;σ + 2j ν A[λ,ν] .

(H.14)

148

Contracting H.11 and dividing by (2−n) gives ˜ να − 16π √−g 0 = Wνα R

1 Tσ (2−n) σ √ √ √ + 4 −g f αν A[ν,α] + Λb −N n+ Λe −g n.

(H.15)

Adding this to H.11 gives 0 = W

νσ

˜ νλ + W R

σν



˜ λν − 16π −g R

µ Tλσ

1 − T σ δσ (n−2) σ λ

√ √ √ − 8 −g f νσ A[ν,λ] + 2Λb −N δλσ + 2Λe −g δλσ .



(H.16)

Taking the non-covariant divergence of this also gives the Lorentz force equation, µ √ ˜ Rσν,λ − 16π ( −g Tλσ ),σ −

¶ √ 1 σ 0 = W ( −g Tσ ),λ (n−2) √ ¢ √ √ ¡ −8 −g −4πj ν A[ν,λ] + f νσ A[ν,λ],σ + 2Λb ( −N ),λ + 2Λe ( −g),λ (H.17) µ ¶ √ 1√ 1 √ σν ˜ σ σν σ = W Rσν,λ − 16π ( −g Tλ );σ − −g g ,λ Tσν + ( −g Sσ ),λ 2 2 µ ¶ √ 1 −8 −g 4πj ν A[λ,ν] + f νσ (3A[ν,λ,σ] − A[σ,ν],λ ) 2 √ √ +Λb −N N aνσ Nσν,λ + Λe −gg νσ gσν,λ (H.18) σν

√ ˜ σν,λ +2A[σ,ν],λ 2 iΛ1/2 + Λb Nσν,λ + Λe gσν,λ ) = Wσν (R b µ ¶ µ ¶ √ √ 1 σν α σ −8π − −g g ,λ Sσν − gσν Sα + ( −g Sσ ),λ 2 √ √ σ −8 −g 4πj ν A[λ,ν] −16π −g Tλ;σ

(H.19)

√ ˜ σν +2A[σ,ν] 2 iΛ1/2 + Λb Nσν + Λe gσν ),λ = Wσν (R b µ ¶ √ √ 1√ σν σν α σ −8π − −g g ,λ Sσν + −g g ,λ gσν Sα + ( −g Sσ ),λ 2 ¢ √ ¡ σ −16π −g Tλ;σ +2j ν A[λ,ν]

(H.20)

149

= Wσν

³

´ √ ˜ σν +2A[σ,ν] 2 iΛ1/2 + Λb Nσν + Λe gσν − 8πSσν ,λ R b

¡ √ ¢ √ √ σ −8π − −g Sσ,λ − ( −g),λ Sαα + ( −g Sσσ ),λ ¢ √ ¡ σ −16π −g Tλ;σ +2j ν A[λ,ν] ¢ √ ¡ σ = −16π −g Tλ;σ +2j ν A[λ,ν] .

(H.21) (H.22)

It is interesting to consider the antisymmetric part of (H.11), ˜ [λ]ν] = 4f ν [σ A[λ],ν] . fˆν [σ R

(H.23)

Comparing this to the antisymmetric part of the normal field equations (2.28), we see that it implies f ν [σ N[λ]ν] = 0. This identity can be proven as follows, √ 8 −g ˆν[σ [λ] √ f N ν] = (N aνσ −N aσν )(N λ ν −Nν λ )−(N aνλ −N aλν )(N σ ν −Nν σ ) (H.24) −N = −N aνσ Nν λ −N aσν N λ ν +N aνλ Nν σ +N aλν N σ ν (H.25) √ 1 −N √ [−N aνσ Nνρ (N aρλ +N aλρ )−N aσν Nρν (N aρλ +N aλρ ) = 2 −g +N aνλ Nνρ (N aρσ +N aσρ )+N aλν Nρν (N aρσ +N aσρ )] (H.26) √ 1 −N √ [−N aνσ Nνρ N aλρ −N aσν Nρν N aρλ = 2 −g +N aνλ Nνρ N aσρ +N aλν Nρν N aρσ ] = 0.

150

(H.27)

Appendix I Alternative derivation of the O(Λ−1 b ) field equations Here we check the results in §2.2-§2.3 by deriving the approximate Einstein and Maxwell equations in different way, using an O(Λ−1 b ) approximation to the Lagrangian density. Inserting the order O(Λ−1 b ) result

√ √ −N ≈ −g (1 + f ρνfνρ /Λb (n−2)) from

(C.12) into the Lagrangian density (2.10) and using (2.3) gives, √ 1 √ ˜ σµ + 2A[σ,µ] 2 iΛ1/2 ) −N N aµσ (R b 16π µ ¶ √ √ 1 1 1 ρν − (n−2)Λb −g 1+ f fνρ − (n−2)Λz −g + Lm 16π Λb (n−2) 16π h √ √ √ i 1 1/2 aµσ ˜ ≈ − −N N (Rσµ + 2A[σ,µ] 2 iΛb ) + (n−2)Λ −g 16π √ 1 ρν f fνρ + Lm . − −g 16π

L ≈ −

(I.1)

(I.2)

The f ρνfνρ term looks superficially like the ordinary electromagnetic term, except that fσµ is defined by (2.22) rather than as 2A[µ,σ] , and there is also a sign difference. ˜ α is identical to (2.10), so Ampere’s law and The dependence of (I.2) on Aµ and Γ σµ 151

the connection equations will be the same as usual. To take the variational derivative √ of (I.2) with respect to −N N aµσ it is convenient to rewrite it as, i √ 1 h√ ˜ σµ + 2A[σ,µ] 2 iΛ1/2 ) + (n−2)Λ√−g −N N aµσ (R b 16π √ 1 √ √ gνα gρτ − −g −gf ρν −gf ατ √ √ + Lm . 16π −g −g

L ≈ −

(I.3)

Using √ gτ (ν gµ)σ ∂(gτ σ / −g) √ √ √ = − −g −g ∂( −N N aµν )

µ

¶ √ √ ∂ ( −gg ρτ gτ σ / −g) √ because =0 , ∂( −N N aµν )

(I.4)

and (2.26,2.22) we have δL 0 = −16π √ δ( −N N aµσ ) √ ˜ σµ + 2A[σ,µ] 2 iΛ1/2 + Λgσµ = R b ¢ gνα gρτ √ √ √ ¡ − −g δµ[ρ δσν] −gf ατ + −gf ρν δµ[α δστ ] √ √ −g −g µ ¶ √ √ √ gν(σ gµ)α gρτ gνα gρ(σ gµ)τ ρν ατ √ √ √ − −g −gf −gf +√ √ √ −g −g −g −g −g −g √ √ gνα gρτ 1 −gf ρν −gf ατ √ √ +gσµ (n−2) −g −g δLm −16π √ δ( −N N aµσ ) √ ˜ σµ + 2A[σ,µ] 2 iΛ1/2 + Λgσµ + Λb fσµ = R b µ ¶ ¶ µ 1 1 ν ρν α −2 fσ fνµ − gσµ f fνρ − 8π Tσµ − gσµ Tα . 2(n−2) (n − 2)

(I.5)

(I.6)

(I.7)

Symmetrizing (I.7) and combining it with its contraction gives the approximate Einstein equations (2.43). Antisymmetrizing (I.7) and doing the same analysis as before gives (2.78) and Maxwell equations (2.47,2.48).

152

Appendix J A weak field Lagrangian density Here we derive an O(Λ−1 b ) Lagrangian density for this theory which depends only on the fields gµν , Aα and θα . The Lagrangian density is valid in the sense that it gives the same field equations to O(Λ−1 b ) as calculated in §2.2-§2.4. However, the weak field Lagrangian density is derived with a somewhat ad-hoc procedure, and since it does not describe LRES theory exactly, and exact solutions would be different, one should exercise caution when using it to make definite conclusions about the theory. Inserting the O(Λ−1 b ) result (C.12) into our Palatini Lagrangian density and using

153

the O(Λ−1 b ) results (2.32,2.3,2.67,2.78,2.76) gives, √ √ √ 1 √ −1/2 ˜ σµ + 2A[σ,µ] 2 iΛ1/2 ) ( −g g σµ + −gf σµ 2 iΛb )(R b 16π µ ¶ √ √ 1 1 1 − (n−2)Λb −g 1+ f ρνfνρ − (n−2)Λz −g 16π (n−2)Λb 16π 1 √ ˜ µ + (n−2)Λ − 2f σµfµσ )− √−g 1 f ρνfνρ −g(R − µ 16π 16π √ −g − (R + (n−2)Λ−f σµfµσ 16π 2 τβ α n 3 − f fβ ;τ ;α − `, α ;α − f[νβ;α] f [νβ ; α] Λb 2(n−2)Λb 2Λb ¶ 32π 2 n 16π ρ ατ − j jρ − f jτ ;α (n−1)(n−2)Λb (n−2)Λb √ −g − (R + (n−2)Λ−f σµfµσ 16π ³ ´ n 4 1 `, α ;α + θρ θρ − f βτ C´βτ αρ f αρ +8πj[β,τ ] − Λb 2(n−2)Λb Λb ¶ 2 32π n 16π − j ρ jρ − f ατ jτ ;α (n−1)(n−2)Λb (n−2)Λb µ √ −g n 4 R + (n−2)Λ−f σµ–Fµσ − `, α ;α + θρ θρ − 16π 2(n−2)Λb Λb ¶ 2 32π n 8π(n−4) ατ ρ − j jρ + f jτ ;α (n−1)(n−2)Λb (n−2)Λb µ √ −g n 4 R + (n−2)Λ−f σµ–Fµσ − `, α ;α + θρ θρ − 16π 2(n−2)Λb Λb ¶ 2 32π n 8π(n−4) ρ ατ τ − j jρ + ((f jτ );α − 4πj jτ ) (n−1)(n−2)Λb (n−2)Λb µ √ −g 32π 2 (n−2) ρ 4 − R + (n−2)Λ−f σµ–Fµσ − j jρ + θρ θρ 16π (n−1)Λb Λ ¶ b n 8π(n−4) ατ − `, α ;α + (f jτ );α . 2(n−2)Λb (n−2)Λb

L2 ≈ −

≈ ≈









(J.1) (J.2)

(J.3)

(J.4)

(J.5)

(J.6)

(J.7)

where we write (2.78) as 1 ρτ ´ f Cρτ σµ , Λb 1 8π(n−2) = 2A[µ,σ] + θ[τ,α] εσµ τ α + j[σ,µ] , Λb Λb (n−1)

fσµ = –Fσµ +

(J.8)

–Fσµ

(J.9)

C´σµαρ = Rσµαρ − gσ[α Rρ]µ + gµ[α Rρ]σ . 154

(J.10)

Removing the total divergences gives √

L2

−g ≈ − 16π

µ

¶ √ 2π(n−2) ρ 4 ρ R + (n−2)Λ−f Fµσ + θ θρ + −g j jρ . (J.11) Λb (n−1)Λb σµ–

Now let us redefine the electromagnetic potential 4π(n−2) A˘µ = Aµ − jµ . Λb (n−1)

(J.12)

In terms of this shifted potential, fµν and –Fµν from (2.78) lose their j[σ,µ] terms and Maxwell’s equations (2.47,2.48) become more exact. This redefinition brings the same j σ jσ term out of Lm for all of the Lm cases. For the classical hydrodynamics case (L.1), µQ σ µ u Aσ − uα gαν uν , m 2 µQ α √ u , = m −g √ 4π(n−2) σ = − −g j jσ . Λb (n−1)

Lm = − jα ∆Lm

(J.13) (J.14) (J.15)

For the spin-0 case (L.2), Lm = Dµ = jα = ∆Lm =

¶ − µ h ¯ 2 ¯← ¯ ψ Dµ D ψ − mψψ , m ← − ← − ∂ iQ ∂ iQ + Aµ , Dµ = − Aµ , µ µ ∂x h ¯ ∂x h ¯ ← − i¯hQ ¯ α (ψD ψ − ψ¯ D α ψ), 2m µ ¶ √ 1h √ 4π(n−2) σ ¯2 iQ 2m 4π(n−2) σ −g j jσ = − −g j jσ . − 2m h ¯ i¯hQ Λb (n−1) Λb (n−1) √

1 −g 2

µ

155

(J.16) (J.17) (J.18) (J.19)

For the spin-1/2 case (L.4), Lm Dµ

¶ ← − ν i¯h ¯ ν 2¯ ¯ = −g (ψγ Dν ψ − ψ D ν γ ψ) − mc ψψ , 2 ← − ← − ∂ ∂ iQ ˜ ˜ †µ − iQ Aµ , = + Γµ + Aµ , Dµ = +Γ µ µ ∂x h ¯ ∂x h ¯ √

µ

¯ α ψ, j α = Qψγ ∆Lm = Here lP =



(J.21) (J.22)

µ ¯ σψ −g i¯ hψγ

(J.20)

iQ h ¯



√ 4π(n−2) σ 4π(n−2) jσ = − −g j jσ . Λb (n−1) Λb (n−1)

(J.23)

p h ¯ G/c3 is the Planck length. The ∆Lm contribution is halved by the j σ jσ

term in the original Lagrangian density so that the total shifted Lagrangian density is ¶ √ µ √ 2π(n−2) ρ −g 4 ρ σµ– ˘ ˘ ˘ L2 ≈ − R + (n−2)Λ− f Fµσ + θ θρ − −g j jρ 16π Λb Λb (n−1) +L˘hydrodynamics + L˘spin−0 + L˘spin−1/2 . . . .

(J.24)

˘ µσ are (J.8,J.9) but without the jα terms, where f˘σµ and –F ˘ σµ + 1 f˘αρ C´σµαρ , f˘σµ = –F Λb ˘ σµ = 2A˘[µ,σ] + 1 θ[τ,α] εσµ τ α . –F Λb

(J.25) (J.26)

Expanding things out and ignoring higher order powers of Cσµαρ /Λb and total divergences, our effective weak-field Lagrangian density becomes L˘2

¶ µ √ gravitational −g ≈ − (R + (n−2)Λ) terms 16π µ ¶ √ √ 2π(n−2) ρ electromagnetic −g ˘[µ σ] ˘ + A ; A[σ,µ] − −g j jρ terms 4π Λb (n−1) µ ¶ √ ¢ −g 1 ¡ [µ σ] Proca ρ − θ ; θ[σ,µ] + Λb θ θρ terms 4π Λ2b ¶ µσρν µ ¶ µ ¶ √ µ −g ˘ 1 coupling 1 τα C βλ − θ[τ,α] εσµ A˘[ρ,ν] + θ[β,λ] ερν A[µ,σ] + terms 4π 2Λb Λb 2Λb +L˘hydrodynamics + L˘spin−0 + L˘spin−1/2 . . . . 156

(J.27)

Even though we call this a weak-field Lagrangian density, it only neglects terms which are < 10−64 of the leading order terms for worst-case field strengths. We can probably assume θα ≈ 0 because Proca plane waves would have a minimum frequency ωP roca =



2Λb ∼



2ωc2 lP ∼



2/lP which exceeds the zero-point cut-off frequency

ωc ∼ 1/lP . Alternatively this field could function as a built-in Pauli-Villars field as discussed in Appendix K. Finally, note that if we take the variational derivative of our shifted Lagrangian density (J.27) with respect to ψ¯ we get the unshifted Dirac equation, · µ ¶ ¸ 1 ∂L ∂L 0 = √ − ,λ c −g ∂ ψ¯ ∂ ψ¯,λ 2πc2 (n−2)Q ρ ˘ ν ψ − mcψ ψγ jρ + i¯hγ ν D = G(n−1) = i¯ hγ ν Dν ψ − mcψ.

(J.28) (J.29) (J.30)

Presumably the same thing occurs with the Klein-Gordon equation. This should be expected because our shift in the electromagnetic potential is only a redefinition and should not result in different field equations. Now let us derive the field equations from the weak field Lagrangian density for the source-free case. The O(Λ−1 b ) Maxwell equations (2.47,2.48) can be derived by setting δL2 /δAν = 0. Using (J.10,J.9,J.11) and f σµ = I–aσµαρ –Fαρ ,

(J.31)

–Iαρσµ = gα[σ gµ]ρ − C´αρσµ /Λb

(J.32)

δ(f σµ –Fµσ ) = 2f σµ δ –Fµσ + f σµ f αρ δ –Iσµαρ ,

(J.33)

157

gives √ δL2 δ –Fµσ 0 = −8π = − −gf σµ = δAν δAν

µ

√ ∂A[σµ] 2 −gf σµ ∂Aν,ω

¶ ,ω

√ = 2( −gf νω ),ω . (J.34)

From this and (2.21) we get Maxwell’s equations (2.47,2.48) as before. The O(Λ−1 b ) Proca equation (2.81) can be derived by setting δL2 /δθν = 0. Using (J.33,J.32,J.10,J.9,J.8,J.11) gives 0 = −8πΛ2b

δL2 δθν

(J.35)

√ δ –Fµσ √ = −Λ2b −gf σµ + −g 4Λb θν δθν µ ¶ ρτ √ √ σµ ∂(εµσ θρ,τ ) = Λb −gf −g 4Λb θν ,ω + ∂θν,ω √ √ = Λb ( −gf σµ εµσ νω ),ω + −g 4Λb θν

(J.36) (J.37) (J.38)

√ √ √ = 2Λb −g(Aµ,σ εµσνω );ω +( −g (4θ[ν ; ω] + f ρτ C´ρτ σµ εµσνω )),ω + −g 4Λb θν (J.39) √ √ √ = 4( −gθ[ν ; ω] ),ω − −g ενωσµ (f ρτ C´ρτ [σµ ),ω] + −g 4Λb θν .

(J.40)

From this we get the Proca equation (2.81) as before. The O(Λ−1 b ) Einstein equations (2.42) can be derived by setting δL2 /δgσµ = 0. First we will deal with the C´λραν term in (J.33,J.32). From (J.10,2.76) we have f λρ f αν δ C´λραν = f λρ f αν δ(Rλραν − gλ[α Rν]ρ + gρ[α Rν]λ )

(J.41)

= fτ ρ f αν δ(Rτ ραν −2δατ Rνρ ) + f λρ f αν δgλτ (Rτ ραν − 2δατ Rνρ ) (J.42) = fτ ρ f αν δ(Rτ ραν −2δατ Rνρ ) + 2δgλτ f λρ f ατ ;ρ;α .

(J.43)

To calculate fτ ρ f αν δ(Rτ ραν−2δατ Rνρ ) we will assume locally geodesic coordinates where Γρσµ = 0. With this method, terms with a Γρσµ factor can be ignored, and covariant derivatives are equivalent to ordinary derivatives, as long as they are not inside a 158

derivative. Then from the definition of the Ricci tensor we have fτ ρ f αν δ(Rτ ραν −2δατ Rνρ ) = 2fτ ρ f αν δΓτρ[ν,α] + 4 `ρν δΓαν[ρ,α] = 2fτ ρ f αν δΓτρν,α + 2 `ρν δΓανρ,α − 2 `ρν δΓανα,ρ

(J.44) (J.45)

= 2fτ ρ f αν (δΓτρν );α + 2 `ρν (δΓανρ );α − 2 `ρν (δΓανα );ρ (J.46) = 2(fτ ρ f αν δΓτρν );α + 2( `ρν δΓανρ );α − 2( `ρν δΓανα );ρ −2(fτ ρ f αν );α δΓτρν − 2 `ρν ;α δΓανρ + 2 `ρν ;ρ δΓανα ,

(J.47)

where `ρν = f ρ τ f τ ν .

(J.48)

The first line of (J.47) is the divergence of a vector, so assuming that δΓανρ = 0 on the boundary of integration, these terms can be dropped. Substituting the Christoffel connection (2.20) into the remaining terms gives fτ ρ f αν δ(Rτ ραν −2δατ Rνρ ) = −(fτ ρ f αν );α g τ β (δgνβ,ρ + δgβρ,ν − δgρν,β ) − `ρν ;α g αβ (δgρβ,ν + δgβν,ρ − δgνρ,β ) + `ρν ;ρ g αβ (δgαβ,ν + δgβν,α − δgνα,β ) = −(f βρ f αν );α 2δgνβ,ρ − `ρν ; β (2δgρβ,ν − δgνρ,β ) + `ρν ;ρ g αβ δgαβ,ν

(J.49) (J.50)

= −(f βρ f αν );α 2(δgνβ );ρ − `ρν ; β (2(δgρβ );ν − (δgνρ );β ) + `ρν ;ρ g αβ (δgαβ );ν (J.51) = −2((f βρ f αν );α δgνβ );ρ − 2( `ρν ; β δgρβ );ν + ( `ρν ; β δgνρ );β + ( `ρν ;ρ g αβ δgαβ );ν + 2(f βρ f αν );α;ρ δgνβ + 2 `ρν ; β ;ν δgρβ − `ρν ; β ;β δgνρ − `ρν ;ρ;ν g αβ δgαβ .

(J.52)

The first line of (J.52) is the divergence of a vector, so assuming that δgσµ = 0 on the boundary of integration, these terms can be dropped. Using (J.52,J.43,J.33,J.32,J.9, 159

2.79,2.76,J.11) and assuming a covariant θρ , the terms of δL2 /δgσµ are then −1 λρ αν δ C´λραν −1 f f = ((f µρ f ασ );α;ρ + (f σρ f αµ );α;ρ + `σν ; µ ;ν + `µν ; σ ;ν 2Λb δgσµ 2Λb ¢ − `µσ ; β ;β − `ρν ;ρ;ν g σµ + f σρ f αµ ;ρ;α + f µρ f ασ ;ρ;α , 1 2

µ 2f

– δgσµ

αρ δ Fρα



(J.53)

1 σρ µτ ν 1 (f ερ θτ, ν + f µρ ερ στ ν θτ, ν − g σµ f αρ ερα τ ν θτ, ν ) (J.54) Λb 2 ¶ µ 3 µ [ρσ ν] 3 σµ 1 3 σ [ρµ ν] [ρα ν] , = − f ρ f ; ; ν − f ρ f ; ; ν + g fαρ f ; ; ν (J.55) Λb 2 2 4 =

µ ¶ 2 δ(θρ θρ ) 2 σ µ 1 9 σρ [νβ µ] σµ ρ − = θ θ = g f[νβ;ρ] f ; + 2g θ θρ (J.56) Λb δgσµ Λb Λb 4 µ ¶ 1 1 σ σ νβ µ µν β βµ ν σµ ρ = (fνβ; +2f ν;β )(f ; +f ; +f ; ) + 2g θ θρ (J.57) Λb 4 µ 1 1 1 1 1 = fνβ; σ f νβ ; µ + fνβ; σ f βµ ; ν + f σ ν;β f νβ ; µ + f σ ν;β f µν ; β Λb 4 2 2 2 ¶ 1 σµ 1 σµ 1 σ βµ ν νβ α αν β + f ν;β f ; − g fνβ; α f ; − g fνβ; α f ; , (J.58) 2 4 2 √ 2 1 δ −g ρ 1 3 σµ − √ θ θρ = − g σµ θρ θρ = g f[νβ; α] f [νβ ; α] , Λb −g δgσµ Λb 8Λb µ ¶ √ 1 δ −g αρ 1 1 √ f –Fρα = g σµ f αρ –Fρα 2 −g δgσµ 4 µ ¶ 1 σµ αρ 1 τν ´ = g f fρα − f Cτ νρα 4 Λb 1 σµ αρ β 1 = g σµ f αρ fρα − g f f [ρ;α];β , 4 2Λb 1 λρ αν δ(gλ[α gν]ρ ) f f = f σρ f µ ρ , 2 δgσµ √ 1 δ( −g((n−2)Λ+R)) (n − 2) σµ √ = Λg − Gσµ . −g δgσµ 2

(J.59) (J.60) (J.61) (J.62) (J.63) (J.64)

Setting the sum of these terms to zero and applying Ampere’s law (2.47) yields the order O(Λ−1 b ) Einstein equations (2.42).

160

Appendix K Proca-waves as Pauli-Villars ghosts? Here we investigate the possibility that the θν field in (2.81,J.27) could function as a built-in Pauli-Villars field. Recall that in quantum electrodynamics, a cutoff wavenumber is often not implemented by simply substituting kc for ∞ in the upper limit of integrals. Instead, Pauli-Villars regularization is often used because it is Lorentz invariant. With the Pauli-Villars method, a ficticious particle is introduced into the Lagrangian which has a huge mass, M say, and which has the opposite sign in the Lagrangian, meaning that it is a ghost, with negative norm or negative energy. This has the same effect as a cutoff wavenumber where kc = M (in quantum electrodynamics natural units). To calculate the electron self-energy for example, this ficticious particle is a ghost Proca particle. When calculating the amplitude of a process involving a photon, the integral associated with every Feynman graph will

161

contain a Feynman propagator of the form Dij =

gij . (p − k)2 + i²

(K.1)

Introducing the ghost Proca particle means that for every Feynman graph that one had originally, there will be a new one where the associated integral contains instead Dij = −

(p −

k)2

gij − M 2 + i²

(note the minus sign).

(K.2)

Because of the way that amplitudes are combined, this has the effect of replacing the original photon propagator in all Feynman graphs with the sum of the two above Dij =

gij gij − . 2 2 (p − k) + i² (p − k) − M 2 + i²

(K.3)

When one integrates over k, the two parts cancel each other for k >> M , effectively cutting off the integral. However, because the additional term of the propagator has the huge mass M in its denominator, this term has virtually no effect for ordinary momenta that can be produced in accelerators and astrophysical phenomena. The point of this is that the ghost Proca particle that seems to come out of the theory is just what is needed for Pauli-Villars regularization. This can be seen from the effective weak-field Lagrangian density for the theory in Appendix J. All that would necessary to do would be to include a coupling of this particle with the electron when spin-1/2 particles are added to the theory. So this ghost Proca particle could actually be a blessing in disguise, because it potentially frees the theory from divergences which must be removed artificially in ordinary quantum electrodynamics. This also illustrates that a ghost particle with mass near the inverse Planck length is a whole different animal than one with an ordinary mass. 162

With further work, it might even be possible to free the theory of its reliance on an externally imposed cutoff frequency. The Pauli-Villars cutoff caused by the ghost Proca particle should also cutoff the calculation of Λb (for photons anyway). Then combing (2.12,2.81) we get ωc = kc =

p 2Λb ,

(K.4)

Λb = Cz lP2 ωc4 = Cz lP2 (2Λb )2 , ⇒ Λb =

1 , 4Cz lP2

(K.5) (K.6)

where lP is the Planck length, ωc is the cutoff frequency (2.13), and Cz comes from (2.14). But we have calculated in (7.15) what Λb must be in the non-Abelian theory, and presumably this should apply for the Abelian theory also. Equating the values from (7.15) and (K.6) gives 1 α Λb = = 2 2 4Cz lP 8lP sin2 θw



2sin2 θw Cz = . α

(K.7)

Using α = e2 /¯ hc = 1/137.036, sin2 θw = .2397 ±.0013 and the definition (2.14) gives ³

´ 4πsin2 θ w fermion boson = 412.8 ± 2. − = spin states spin states α

(K.8)

The result (K.8) is interesting, partly because the theory predicts that the difference should be an integer, and this potentially allows the theory to be proven or disproven. At present the weak mixing angle θw cannot be measured accurately enough to determine whether we are seeing an integer or not. The issue is also complicated because the value of sin2 θw /α “runs”, meaning that its value depends logarithmically on the energy at which it is measured. To really do an accurate calculation we would need to use its value at the same cutoff frequency ωc used to calculate Λz , but this could be 163

done. It is likely that measurement accuracy will improve enough in the near future so that we can determine whether (K.8) is consistent with an integer. The result (K.8) is also interesting because it might select among the different possibilities of matrix size for the non-Abelian version of our theory. For the non-Abelian theory we used 2×2 matrices in order to get Einstein-Weinberg-Salam theory, but we had the choice of using any matrix size “d”, corresponding to U(1)⊗SU(d) instead of U(1)⊗SU(2). Each choice of “d” will result in different numbers of fermion and boson spin states. It would be very nice if some choice of “d” agreed with (K.8). The choice d = 5 with U(1)⊗SU(5) is particularly interesting because SU(5) has long been considered as a way of unifying the strong and weak forces in the U(1)⊗SU(2)⊗SU(3) gauge structure of the Standard Model. However, the calculation of the left-hand side of (K.8) is complicated, and it is unclear whether to include scalar particles and gravitons, and it is even more unclear how to account for possible additional particles associated with a non-Abelian gνµ . For the Standard Model the left-hand side of (K.8) works out to about 60. To do this rigorously, we would also need to include Pauli-Villars ghosts corresponding to electrons. Of course the theory only approximates electro-vac EinsteinMaxwell theory, so we must add in spin-1/2 particles (one can think of this as 1st quantization of an electric monopole solution). In any case, for every spin-1/2 particle that are added to the theory, it would be easy to also add in a corresponding Pauli-Villars ghost. Surely having Pauli-Villars ghosts as an inherent part of quantum electrodynamics can’t be any worse than introducing ficticious particles just to make divergent integrals come out finite, and then forgetting about them. 164

Appendix L Lm, Tµν , j µ and kinetic equations for spin-0 and spin-1/2 sources Here we display the matter Lagrangian Lm for the classical hydrodynamics, spin-0 and spin-1/2 cases, and we derive the energy-momentum tensors and charge currents for each. Then we derive the Klein-Gordon equation and Dirac equation, and we derive the continuity equation and Lorentz force equation from the Klein-Gordon equation. All of these results are shown to be identical to ordinary Einstein-Maxwell theory and one-particle quantum mechanics. For the classical hydrodynamics case we can form a rather artificial Lm which depends on a mass scalar density µ and a velocity vector uν , neither of which is constrained (that is we will not require δL/δµ = 0 or δL/δuσ = 0), Lm = −

µ µQ ν u Aν − uα gασ uσ . m 2

165

(L.1)

For the spin-0 case as in [49], matter is represented with a scalar wave-function ψ, ¶ − µ h ¯ 2 ¯← ¯ = ψ Dµ D ψ − mψψ , m ← − ← − ∂ ∂ iQ iQ = + Aµ , Dµ = − Aµ . µ µ ∂x h ¯ ∂x h ¯ √

Lm Dµ

1 −g 2

µ

(L.2) (L.3)

For the spin-1/2 case as in [49], matter is represented by a four-component wavefunction ψ, and things are defined using tetrads e(a) σ , Lm =



µ −g

¶ ← − i¯ h ¯ σ σ ¯ (ψγ Dσ ψ − ψ¯ D σ γ ψ) − mψψ , 2   

γ σ = γ (a) e(a) σ

g

(a)(b)

= e

(a)

τ

e

(b)

,

σ

g

τσ

(L.4) 

I 0   0 σi  ,  , γ (i) =  γ (0) =      −σi 0 0 −I    1 = (γ (a) γ (b) + γ (b) γ (a) ) =   2

1

(L.5)

  , (L.6) 

−1 −1 −1

1 g τ σ = e(a) τ e(b) σ g (a)(b) = (γ τ γ σ + γ σ γ τ ), 2 (a)

e(a) τ e(c) τ = δ(c) Dµ ˇµ Γ

,

e(a) τ e(a) σ = δτσ ,

(L.7) (L.8)

← − ← − ∂ iQ ∂ ˇ ˇ † − iQ Aµ , = + Γµ + Aµ , Dµ = +Γ (L.9) µ µ µ ∂x h ¯ ∂x h ¯ 1 (a)(b) σ i = Σ e(a) e(b)σ;µ , Σ(a)(b) = (γ (a) γ (b) − γ (b) γ (a) ). (L.10) 2 2

In the equations above, m is mass, Q is charge and the σi are the Pauli spin matrices. ← − In (L.3,L.9) the conjugate derivative operator Dµ is made to operate from right to left to simplify subsequent calculations. The spin-0 and spin-1/2 Lm 0 s are the ordinary expressions for quantum fields in curved space[49]. To calculate Tµν for the spin-1/2 case we will need to first calculate the derivative √ √ √ ∂( −g e(a) τ )/∂( −N N aµν ). Multiplying (L.7) by −g e(b) σ and taking its derivative

166

with respect to



−g g µν gives √

−g g τ σ e(b) σ =

τ σ (b) δ(µ δν) e σ

+





−g e(a) τ g (a)(b) ,

√ ∂( −g e(a) τ ) (a)(b) ∂e(b) σ √ √ −g g = g . ∂( −g g µν ) ∂( −g g µν ) τσ

(L.11) (L.12)

√ Taking the derivative of (L.8) with respect to −g g µν and using (2.26) gives √ τ gνµ (a) √ ∂e(a) τ (a) ∂( −g e(c) ) √ 0=e τ √ − δ + −g e(c) τ . (L.13) ∂( −g g µν ) (n−2) (c) ∂( −g g µν ) Substituting (L.13) into (L.12) we finally get √ √ ∂( −g e(a) τ ) (a)(b) ∂e(b) σ (b) τ τσ √ √ δ = −g g g − e (L.14) (ν µ) ∂( −g g µν ) ∂( −g g µν ) √ ¶ µ τ gνµ (b) τσ (b) ∂( −g e(a) ) (a) = g e σ−e τ √ e σ (L.15) (n−2) ∂( −g g µν ) √ τ gνµ (b)τ (a)(b) ∂( −g e(a) ) √ = e −g , (L.16) (n−2) ∂( −g g µν ) √ √ ∂( −g e(a) τ ) ∂( −g e(a) τ ) 1 gνµ τ √ √ = = e(a)(ν δµ) + e(a) τ . (L.17) µν aµν ∂( −gg ) 2 2(n−2) ∂( −N N ) From (2.30) we see that Sνµ and Tνµ are different for each Lm case. For the classical hydrodynamics case (L.1), Sνµ Tνµ

µ ¶ µ 1 α = √ uν uµ − gνµ u uα , −g (n − 2) µ = √ uν uµ . −g

For the spin-0 case (L.2) as in [49], µ ¶ − 1 1 2 ¯← 2¯ Sνµ = h ¯ ψ D(ν Dµ) ψ − gνµ m ψψ , m (n−2) ´ − − σ 1 ³ 2 ¯← 1 2 ¯← 2¯ Tνµ = h ¯ ψ D(ν Dµ) ψ − gνµ (¯ h ψ Dσ D ψ−m ψψ) . m 2

(L.18) (L.19)

(L.20) (L.21)

For the spin-1/2 case (L.4) as in [49], ¶ µ ← − σ ← − i¯h ¯ 1 σ ¯ ¯ ¯ gνµ (ψγ Dσ ψ − ψ Dσ γ ψ) ,(L.22) Sνµ = ψγ(ν Dµ) ψ − ψ D(µ γν) ψ − 2 (n−2) ´ ← − i¯h ³ ¯ Tνµ = ψγ(ν Dµ) ψ − ψ¯ D(µ γν) ψ . (L.23) 2 167

← − ¯ σ , the energyNote that in the purely classical limit as i¯hDσ ψ → pσ ψ , −i¯hψ¯ Dσ → ψp momentum tensors (L.21) for spin-0 and (L.23) for spin-1/2 both go to the classical hydrodynamics case (L.19). From (2.46) we see that j τ is different for each Lm case. For the classical hydrodynamics case (L.1), µQ α √ u . m −g

(L.24)

← − i¯hQ ¯ α (ψD ψ − ψ¯ D α ψ). 2m

(L.25)

jα = For the spin-0 case (L.2) as in [49], jα =

For the spin-1/2 case (L.4) as in [49], ¯ α ψ. j α = Qψγ

(L.26)

For the spin-0 case, the Klein-Gordon equation is obtained by setting δL/δ ψ¯ = 0, 0 = = = =

· µ ¶ ¸ ∂L −2 ∂L √ − ,λ −g ∂ ψ¯ ∂ ψ¯,λ · 2µ ¶ ¸ √ h ¯ iQ h ¯2 µ λ − − Aµ D ψ − mψ − √ ( −gD ψ),λ m h ¯ m −g ¸ · ¶ 2 µ √ h ∂ −1 −¯ iQ µ 2 √ −gD ψ − m ψ + Aµ m −g ∂xµ h ¯ · 2 ¸ √ 1 h ¯ µ 2 √ Dµ −gD + m ψ. m −g

(L.27) (L.28) (L.29) (L.30)

The conjugate Klein-Gordon equation is found by setting δL/δψ = 0, µ ¶ ¸ · −2 ∂L ∂L 0 = √ − ,λ −g ∂ψ ∂ψ,λ · ¸ −µ √ ← − h 1 ¯ ← ¯2 2 = ψ D −g Dµ √ +m . m −g 168

(L.31) (L.32)

This is just the complex conjugate of the Klein-Gordon equation (L.30) if ψ¯ = ψ ∗ . For the spin-1/2 case, the Dirac equation is found in a similar manner, · µ ¶ ¸ ∂L 1 ∂L 0 = √ − ,λ −g ∂ ψ¯ ∂ ψ¯,λ = i¯ hγ σ Dσ ψ − mψ.

(L.33) (L.34)

The conjugate Dirac equation is, · µ ¶ ¸ 1 ∂L ∂L 0 = √ − ,λ −g ∂ψ ∂ψ,λ ← − ¯ = −i¯hψ¯ Dσ γ σ − mψ.

(L.35) (L.36)

Both the Klein-Gordon and Dirac equations match those of ordinary one-particle quantum mechanics in curved space[49]. Note that for the spin-0 case, instead of deriving the continuity equation (2.49,L.25) from the divergence of Ampere’s law, it can also be derived from the Klein-Gordon equation. Using (L.30,L.32,L.3,L.25) we get, 0 = = = = = =

¶ ¸ ¶ µ · µ one side of conjugate one side of iQ ¯ ψ (L.37) − ψ Klein−Gordon equation Klein−Gordon equation 2¯h µ 2 ¶ √ ← −µ √ ← − h iQ ¯ h ¯ ¯2 µ 2 2 ψ √ Dµ −gD + m − D −g Dµ √ −m ψ (L.38) 2m¯ h −g −g µµ ¶ µ ¶¶ √ − ← − √ ← − 1 1 i¯hQ ¯ ← ψ Dµ + √ Dµ −g Dµ − D µ −g Dµ √ +Dµ ψ (L.39) 2m −g −g Ãà ← ! !! à − ← − √ √ ← − i¯hQ ¯ ∂ 1 ∂ 1 ∂ ∂ ψ +√ −g Dµ − D µ −g µ √ + µ ψ (L.40) µ µ 2m ∂x −g ∂x ∂x −g ∂x à ! ← − √ ← − ∂ 1 i¯hQ 1 ∂ ¯√ √ (ψ −gDµ ψ) − (ψ¯ D µ −gψ) µ √ (L.41) µ 2m −g ∂x ∂x −g µ ¶ √ i¯ ← − 1 ∂ hQ ¯ µ µ √ (ψD ψ − ψ¯ D ψ) (L.42) −g −g ∂xµ 2m

= j µ ;µ .

(L.43)

169

Similarly, instead of deriving the Lorentz-force equation (4.10,L.21) from the divergence of the Einstein equations, it can also be derived from the Klein-Gordon equation. Using (L.30,L.32,L.3,L.21) we get, µ 0 = = =

=

=

¶ ¶ ← − µ one side of conjugate Dρ ψ ψ¯ Dρ one side of (L.44) + Klein−Gordon equation Klein−Gordon equation 2 2 µ ¶ ¶ − µ 2 2 ¯← √ √ ← − ← − D ψ ψ D h ¯ h ¯ ρ ρ λ 2 λ 2 √ Dλ −gD + m ψ (L.45) ψ¯ D −g Dλ √ +m + −g 2m 2m −g " √ √ ← − ← − 1 ∂( −gDλ ψ) h ¯ 2 ∂(ψ¯ D λ −g) 1 ¯ √ D ψ + ψ D ρ ρ√ 2m ∂xλ −g −g ∂xλ µ ¶ µ ¶ ¸ ← − ← − iQ iQ λ λ +ψ¯ D − Aλ Dρ ψ + ψ¯ Dρ Aλ D ψ h ¯ h ¯ − m ¯ ¯← + (ψD (L.46) ρ ψ + ψ Dρ ψ) 2 " √ ← − √ ← − 1 ∂( −gDλ ψ) h ¯ 2 ∂(ψ¯ D λ −g) 1 ¯ √ Dρ ψ + ψ Dρ √ 2m ∂xλ −g −g ∂xλ µ µ ¶ ¶µ ¶ ← −λ ← −λ iQ ∂ψ iQ iQ ¯ ¯ + ψ D − Aλ + ψ D − Aλ Aρ ψ h ¯ ∂xρ h ¯ h ¯ µ ¶ µ ¶ µ ¶ ¸ iQ ∂ ψ¯ iQ iQ λ λ ¯ + ρ Aλ D ψ + − Aρ ψ Aλ D ψ ∂x h ¯ h ¯ h ¯ µ ¶ ∂ ψ¯ m ¯ ∂ψ ψ ρ + ρψ (L.47) + 2 ∂x ∂x " √ ← − √ ← − h ¯ 2 ∂(ψ¯ D λ −g) 1 1 ∂( −gDλ ψ) ¯ √ Dρ ψ + ψ Dρ √ 2m ∂xλ −g −g ∂xλ µ µ µ ¶ ¶ ¶ ← − ← − ← − iQ ∂ψ iQ ∂ψ iQ λ λ λ − ψ¯ D Aλ − ψ¯ D − Aρ + ψ¯ D − Aρ Dλ ψ ρ λ h ¯ ∂x h ¯ ∂x h ¯ µ µ ¶ µ ¶ ¯ ¶ ¸ ¯ ← −λ iQ iQ ∂ψ λ iQ ∂ψ λ λ ¯ D ψ − ψ D − Aρ D ψ − ρ − Aλ D ψ + − Aρ ∂x h ¯ h ¯ ∂xλ h ¯ ¯ m ∂(ψψ) (L.48) + 2 ∂xρ

170

" √ ← − √ ← − 1 ∂( −gDλ ψ) h ¯ 2 ∂(ψ¯ D λ −g) 1 ¯ √ Dρ ψ + ψ Dρ √ = 2m ∂xλ −g −g ∂xλ µ ¶ µ ¶ ¶ µ ← −λ ∂ iQ ∂ iQ 2iQ ¯ Aρ ψ − ρ Aλ ψ − A[ρ,λ] ψ + ψD ∂xλ h ¯ ∂x h ¯ h ¯ µ µ µ ¶ ¸ ¶ ¶ ∂ iQ ¯ iQ ¯ ∂ 2iQ λ ¯ + − Aρ ψ − ρ − Aλ ψ + A[ρ,λ] ψ D ψ ∂xλ h ¯ ∂x h ¯ h ¯ ¯ m ∂(ψψ) (L.49) + ∂xρ " 2← √ − √ ← − h ¯ 2 ∂(ψ¯ D λ −g) 1 1 ∂( −gDλ ψ) ¯ √ √ = D ψ + ψ D ρ ρ 2m ∂xλ −g −g ∂xλ ← − ∂g νλ ← − ∂g µλ ¯ ¯ + ψ Dν Dλ ψ − ψ Dµ ρ Dλ ψ ∂xρ ∂x µ ¶ Ã ¯← − − ! ¯← ← − ∂(D ψ) ∂(D ψ) ∂( ψ D ) ∂( ψ D ρ λ ρ λ) + ψ¯ D λ − + − Dλ ψ λ ρ λ ρ ∂x ∂x ∂x ∂x ¸ ¯ ← −λ 2iQ ¯ λ m ∂(ψψ) ¯ + (ψD ψ − ψ D ψ)A[ρ,λ] + (L.50) h ¯ 2 ∂xρ · ¸ ´ λν ← − ← − ← − ← − 1 h ¯2 ∂ ³√ ∂g ∂ λ λ λ √ = −g(ψ¯ D Dρ ψ+ ψ¯ Dρ D ψ) + ψ¯ Dλ ρ Dν ψ − ρ ψ¯ Dλ D ψ 2m −g ∂xλ ∂x ∂x ¯ ← − m ∂(ψψ) iQ¯h ¯ λ + (ψD ψ − ψ¯ D λ ψ)A[ρ,λ] (L.51) + ρ 2 ∂x m i − ← − ← − 1 h 2 ¯← ¯ ,ρ = h ¯ (ψ D λ Dρ ψ + ψ¯ Dρ Dλ ψ);λ − (¯h2 ψ¯ Dλ Dλ ψ − m2 ψψ) 2m ← − iQ¯ h ¯ λ + (ψD ψ − ψ¯ D λ ψ)A[ρ,λ] (L.52) m = Tρ;λ λ + 2j λ A[ρ,λ] .

(L.53)

Presumably, similar results occur for the spin-1/2 case, but this was not verified.

171

Appendix M Alternative ways to derive the Einstein-Schr¨ odinger theory The original Einstein-Schr¨odinger theory can be derived from many different Lagrangian densities. In fact it results from any Lagrangian density of the form, bλ , Nρτ ) = − L(Γ ρτ

h √ 1 √ ˜ c1 Γ ˜ α + 2A[ν,µ] 2 iΛ1/2 ) −N N aµν (Rνµ (Γ)+ α[ν,µ] b 16π i + (n−2)Λb , (M.1)

where c1 , c2 , c3 are arbitrary constants and σ ˜α ˜ ασα − Γ ˜ να ˜ = Γ ˜ ανµ,α − Γ ˜ ανα,µ + Γ ˜ σνµ Γ Γσµ , Rνµ (Γ)

˜α = Γ bα + Γ νµ νµ

h i 2 bσ +(c2 − 1) δ α Γ bσ , c2 δµα Γ [σν] ν [σµ] (n−1)

bσ /c3 . Aν = Γ [σν]

(M.2) (M.3) (M.4)

Contracting (M.3) on the right and left gives ˜ αβα = Γ

h i 1 bααβ − (c2 n + c2 − n)Γ bαβα = Γ ˜ ααβ , (c2 n + c2 − 1)Γ (n−1) 172

(M.5)

˜ α has only n3 − n independent components. Also, from (M.3,M.4) we have so Γ νµ ¤ 2c3 £ α c2 δµ Aν + (c2 − 1)δνα Aµ , (n−1)

bα = Γ ˜α − Γ νµ νµ

(M.6)

˜ ανµ and Aν fully parameterize Γ bανµ and can be treated as independent variables. so Γ ˜ α = 0 and δL/δAν = 0 must give the same field equations as Therefore setting δL/δ Γ νµ bανµ = 0. Because the field equations can be derived in this way, the constants δL/δ Γ c2 and c3 are clearly arbitrary, and because of (2.58) with j σ = 0, the constant c1 is also arbitrary. √ 1/2 For c1 = 1, c2 = 1/2, c3 = −(n−1) 2 iΛb , (M.1) reduces to bλ , Nρτ ) = − L(Γ ρτ

h i 1 √ b + (n−2)Λb , −N N aµν Rνµ (Γ) 16π

(M.7)

√ which is our original Lagrangian density (2.2) without the Λz −g and Lm terms, where we have the invariance properties from (2.18,2.19), ˜α → Γ ˜α , Γ bα → Γ bα , Nνµ → Nµν , N aµν→ N aνµ ⇒ L → L, Aν → −Aν , Γ νµ µν νµ µν

(M.8)

√ h ¯ 1/2 α ˜ αρτ → Γ ˜ αρτ , Γ bαρτ → Γ bαρτ + 2¯h δ[ρ φ,α , Γ φ,τ ] 2 iΛb ⇒ L → L. Q Q

(M.9)

Aα → Aα −

˜α = Γ ˜α = Γ bα , and from (2.57,2.28,M.1) the field equations For this case we have Γ σα ασ (ασ) require a generalization of the result L,σ −Γαασ L = 0 that occurs with the Lagrangian density of ordinary vacuum general relativity, that is bα L = 0 or L,σ −Re(Γ bα )L = 0. L,σ − Γ (ασ) ασ

(M.10)

√ 1/2 For the alternative choice, c1 = 0, c2 = n/(n + 1), c3 = −(n − 1) 2 iΛb /2, we have bα and from (2.57,2.28,M.1) the field equations require ˜α = Γ ˜α = Γ Γ ασ ασ σα bα L = 0. L,σ − Γ ασ 173

(M.11)

√ 1/2 For the alternative choice c1 = 1, c2 = 0, c3 = −(n−1) 2 iΛb , (M.1) reduces to h i 1 √ λ aµν b b L(Γρτ , Nρτ ) = − −N N