Next Article in Journal
The Approximation of the Nonlinear Singular System with Impulses and Sliding Mode Control via a Singular Polynomial Fuzzy Model Approach
Next Article in Special Issue
The σ Cohomology Analysis for Symmetric Higher-Spin Fields
Previous Article in Journal
Comparative Study of Peridynamics and Finite Element Method for Practical Modeling of Cracks in Topology Optimization
Previous Article in Special Issue
Circuit Complexity from Cosmological Islands
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Noether Currents and Maxwell-Type Equations of Motion in Higher Derivative Gravity Theories

Yukawa Institute for Theoretical Physics, Kyoto University, Kyoto 606-8502, Japan
Symmetry 2021, 13(8), 1408; https://doi.org/10.3390/sym13081408
Submission received: 29 June 2021 / Revised: 23 July 2021 / Accepted: 27 July 2021 / Published: 2 August 2021
(This article belongs to the Special Issue Manifest and Hidden Symmetries in Field and String Theories)

Abstract

:
In general coordinate invariant gravity theories whose Lagrangians contain arbitrarily high order derivative fields, the Noether currents for the global translation and for the Nakanishi’s I O S p ( 8 | 8 ) choral symmetry containing the BRS symmetry as its member are constructed. We generally show that for each of these Noether currents, a suitable linear combination of equations of motion can be brought into the form of a Maxwell-type field equation possessing the Noether current as its source term.

1. Introduction

The equation of motion for the Yang-Mills field A μ a in the covariant gauge is given in the form
D ν F μ ν a + μ B a i g f a b c μ c ¯ b · c c = g j μ a ,
where D ν F μ ν a is the covariant divergence of the field strength F μ ν a ; j μ a is the color current from the matter field; and B a , c a and c ¯ a are Nakanishi-Lautrup (NL), Faddeev-Popov (FP) ghost, and anti-ghost fields, respectively. This equation was first noted by Ojima [1] to be rewritten into the form of the Maxwell-type equation of motion:
ν F a μ ν + { Q B , D μ c ¯ a } = g J μ a .
Here, Q B is the BRS charge, and J μ a in the RHS is the Noether current for the global gauge transformation (=color rotation) under which all the gauge fields A μ a , NL, and FP ghost fields, B a , c a , c ¯ a , transform as adjoint representations, given by
J μ a = ( A ν × F ν μ ) a + j μ a + ( A μ × B ) a i ( c ¯ × D μ c ) a + i ( μ c ¯ × c ) a ,
with ( A × B ) a f a b c A b B c . This form of YM field Equation (2) is particular, firstly, in the simple divergence form for the field strength, ν F a μ ν , and secondly, in the BRS exact form for the NL and FP ghost contribution terms.
This form of YM field equation, which we call the “Maxwell-type equation of motion”, played very important roles in discussing [2,3]
  • the existence of an elementary BRS quartet of asymptotic fields;
  • the spontaneous breaking of color symmetry and the Higgs phenomenon;
  • the unbroken color symmetry and color confinement.
From the technical viewpoint, it was also useful to simplify the computations of equal time commutation (ETC) relations for some field variables, as well as to derive Ward–Takahashi identities.
In addition, in gravity theory, there is a beautiful canonical formulation given by Nakanishi in a series of papers [4,5,6,7,8,9,10] based on the Einstein–Hilbert action with BRS gauge fixing in de Donder gauge. It is summarized in his textbook [11], co-authored with Ojima. He remarked there that the Einstein gravity field equation can also be rewritten in the form of the Maxwell-type [12]. In this formulation, he also found a beautiful theorem [11] together with Ojima in which the graviton can be identified with a Nambu–Goldstone (NG) massless tensor particle accompanying the spontaneous breaking of G L ( 4 ) symmetry down to S O ( 1 , 3 ) Lorentz symmetry, thus proving the exact masslessness of the graviton in Einstein gravity theory. (Ogievetsky, independently, identified the graviton with the Nambu–Goldstone tensor in his non-linear realization theory for G L ( 4 ) / S O ( 1 , 3 ) [13]). This is a gravitational extension of the Ferrari–Picasso theorem [14], which proves that the photon is an NG vector boson accompanying the spontaneous breaking of a vector-charge Q μ symmetry, corresponding to the gauge symmetry with a transformation parameter linear in x μ . Nakanishi also found in his G L ( 4 ) -invariant de Donder gauge that there exists an 8 + 8 dimensional Poincaré-like I O S p ( 8 | 8 ) supersymmetry, which he called choral symmetry, containing (as its member) BRS and FP ghost scale symmetries as well as the G L ( 4 ) and rigid translation corresponding to the GC transformation with transformation parameter ε μ linear in x μ , ε μ = a ν μ x ν + b μ .
However, this work is a formal theory based on the Einstein–Hilbert action. It is perturbative non-renormalizable and may not give a well-defined theory, although there is a possibility that it may satisfy the so-called asymptotic safety [15] and give a UV complete theory. (There recently appeared an interesting paper [16] which proposes a novel perturbative approach to the Einstein–Hilbert gravity using the quadratic gravity terms as regulators which, the authors claim, can eventually be removed without harm.)
On the other hand, however, there are many investigations of higher derivative gravity theories. In particular, quadratic gravity [17,18,19] attracted much attention in connection with the perturbative renormalizability [20], Weyl invariant theory [21,22,23], and asymptotic freedom [24,25,26].
These higher derivative theories suffer from the massive (negative metric) ghost problem in the perturbative regime, although there have been many proposals for possible ways out (see, e.g., [19] for a review). This ghost problem is, however, outside the scope of this paper.
Even if we are much less ambitious than making gravity theory UV-complete, we still have several motivations to consider higher derivative gravity theories.
From the low energy effective field theory viewpoint, it is quite natural to consider the actions containing higher and higher order derivative fields, successively, from low to high energies. The Einstein–Hilbert action is the lowest derivative order, the quadratic gravity actions are the next derivative order, and so on.
Or, alternatively, one may simply want a gravity theory with a UV cut-off M valid only in the low energy region E < M . A simple momentum cut-off does not work here since it breaks the GC-invariance. Pauli–Villars regulators respecting the GC-invariance can be supplied by considering the covariant higher derivative terms. As noted by Stelle, the gravity field propagator behaves as ∼ 1 / p 4 in the quadratic gravity and sufficiently cuts off the UV contribution to make the theory renormalizable in 4D.
For regulators to work sufficiently enough to make all the quantities finite in 4D, however, the propagator must drop as fast as ∼ 1 / p 6 . Such behavior would be supplied, for instance, by the quadratic term of covariant quantities which contain third order derivatives of the gravity field.
In this paper, we will consider a general gravity theory which is invariant under the general coordinate (GC) transformation and contains arbitrarily high order derivatives of gravity and matter fields, and we
  • derive a concrete form of the Noether current for the rigid translation, i.e., energy momentum tensor;
  • derive the Maxwell-type gravity equation of motion in a gauge-unfixed, i.e., classical system;
  • derive the Maxwell-type equation analogous to Equation (2) in a gauge-fixed quantum system in the de Donder–Nakanishi gauge;
  • and derive the Noether currents of the I O S p ( 8 | 8 ) symmetry present in the de Donder-Nakanishi gauge.
The original motivation for the present author to consider this problem is to give a sound proof for the existence theorem [27] of the massless graviton, claiming that there should exist a spin 2 massless graviton in any GC invariant theory, insofar as it realizes a translational invariant vacuum with a flat Minkowski metric. This is a generalization of the Ferrari–Picasso theorem for the massless photon and the Nakanishi–Ojima theorem for the massless graviton. Those theorems were proved explicitly assuming the renormalizable QED and Einstein gravity, respectively. To prove the existence theorem generally, however, it is necessary to have the Maxwell-type gravity equation of motion in any GC invariant system assuming no particular form of action.
This paper is organized as follows. In Section 2, we present a totally general classical system containing arbitrarily high order derivative fields, which is only assumed to be GC transformation invariant. To treat such a system, we introduce a series of generalized both-side derivatives and prove some formulas they satisfy. Based on these, we derive an expression for the energy-momentum tensor for such a general system as the Noether current for the translation invariance and show that the gravity field equation of motion can be cast into the form of the Maxwell-type equation. In Section 3, these results are generalized in the gauge-fixed system by adopting the G L ( 4 ) -invariant de Donder gauge à la Nakanishi. In Section 4, using the same technique, we show that each of the Noether currents of the I O S p ( 8 | 8 ) symmetry can be written in a form of the source current of a Maxwell-type equation. Section 5 is devoted to the conclusion. Some technical points on O S p transformations are discussed in Appendix A and Appendix B. In Appendix A, the O S p ( 8 | 8 ) transformation of the gauge-fixing plus Faddeev–Popov term is computed for the x μ -dependent transformation parameter. In Appendix B, to obtain some familiarity with the O S p -symmetry, we briefly study the simplest model, a O S p ( 2 | 2 ) -invariant scalar field system on flat Minkowski background; the O S p ( 2 | 2 ) Noether current is derived, and the I O S p ( 2 | 2 ) algebra is confirmed from the canonical (anti-)commutation relations.

2. Gravity Equation of Motion in a Generic Higher Derivative System

We consider a generic system whose action contains higher order derivative fields up to the N-th order μ 1 μ N :
S [ ϕ ] = d 4 x L ϕ , μ ϕ , μ ν ϕ , μ ν ρ ϕ , , μ 1 μ N ϕ ,
where ϕ stands for a collection of fields { ϕ j } (whose index i may be suppressed when unimportant), and we use abbreviations such as
μ 1 μ 2 μ n ϕ j μ 1 μ 2 μ n ϕ j ,
L j ; μ 1 μ 2 μ n L ( μ 1 μ 2 μ n ϕ j ) | weight 1 .
The suffix ‘weight 1’ in the latter means that we always keep the weight as one, irrespective of whether the n indices μ 1 , μ 2 , , μ n take the same values or not; namely, for the case L = a μ ν μ ν ϕ , for instance, L / ( 11 ϕ ) = a 11 and L / ( 12 ϕ ) = a 12 + a 21 , but we define  L / ( μ ν ϕ ) | weight 1 = ( a μ ν + a ν μ ) / 2 ! always. The functional derivative of the action S with respect to ϕ j is given by
δ S δ ϕ j = L ϕ j μ L j ; μ + μ ν L j ; μ ν μ ν ρ L j ; μ ν ρ + = n = 0 N ( ) n μ 1 μ n L j ; μ 1 μ n ,
where N is the highest order n of the derivative fields μ 1 μ n ϕ contained in L (so that L j ; μ 1 μ n = 0 for n N + 1 ), and, for the n = 0 case of empty set { μ 1 μ n } , μ 1 μ n = 1 and L j ; μ 1 μ n = L / ϕ j are understood. The Euler–Lagrange equations are given by δ S / δ ϕ j = 0 .
The Lagrangian generally changes under an infinitesimal transformation ϕ ϕ + δ ϕ , as
δ L = L ϕ j δ ϕ j + L j ; μ μ δ ϕ j + L j ; μ ν μ ν δ ϕ j + = n = 0 N L j ; μ 1 μ n μ 1 μ n δ ϕ j ,
where summation over the repeated j is also implied. We consider the system which is invariant under the gauge transformation taking the form
δ ϕ j ( x ) = G ρ j ε ρ ( x ) + T ρ j μ μ ε ρ ( x ) .
For the GC transformation x ρ x ρ = x ρ ε ρ ( x ) , this field transformation reads more explicitly
δ ϕ j ( x ) = ε ρ ( x ) ρ ϕ j + ϕ j ρ μ μ ε ρ ( x ) i . e . , G ρ j = ρ ϕ j , T ρ j μ = ϕ j ρ μ .
(Here, we are taking ε ρ with an opposite sign to Nakanishi’s so that the definitions of G ρ j , T ρ j μ and ϕ j ρ μ all have an opposite sign to Nakanishi’s.) For a general tensor field ϕ j = T ν 1 ν p σ i σ q , the symbol ϕ j ρ μ is defined by [11]:
T ν 1 ν p σ 1 σ q ρ μ = i = 1 p δ ν i μ T ν 1 ν i 1 ρ ν i + 1 ν p σ 1 σ q j = 1 q δ ρ σ j T ν 1 ν p σ 1 σ i 1 μ σ i + 1 σ q .
The GC invariance of the system implies that the Lagrangian is a scalar density so that the change of L is given by a total divergence:
δ L = μ ( L ε μ ) .
For the GC transformation δ ϕ in Equation (8), we can equate expression (11) for δ L with Equation (7) and obtain an identity:
n = 0 N L j ; μ 1 μ n μ 1 μ n G ρ j ε ρ ( x ) + T ρ j μ μ ε ρ ( x ) μ ( L ε μ ) = 0 .
This equation, if expanded in a power of derivatives μ on the gauge transformation parameter (function) ε ρ ( x ) , yields
n = 0 N L j ; μ 1 μ n μ 1 μ n G ρ j ρ L ε ρ + n = 0 N 1 ( n + 1 ) L j ; μ α 1 α n α 1 α n G ρ j + n = 0 N L j ; α 1 α n α 1 α n T ρ j μ L δ ρ μ μ ε ρ + k = 1 N K k ν 1 ν k μ ρ ν 1 ν k μ ε ρ ( x ) = 0 ,
where K k ’s ( k = 0 , 1 , 2 , ) are defined by
K k ν 1 ν k μ ρ n = 0 N k 1 n + k + 1 C k + 1 L j ; ν 1 ν k μ α 1 α n α 1 α n G ρ j + n = 0 N k n + k C k L j ; ν 1 ν k α 1 α n α 1 α n T ρ j μ ,
with n C k denoting the binomial coefficient n k = n ! / k ! ( n k ) ! . Note that the first summation term n = 0 N 1 ( ) + n = 0 N ( ) in the coefficient of μ ε ρ in the second line in Equation (13) is only identical with the quantity K k ν 1 ν k μ ρ for the case k = 0 .
Since the functions ε ρ , μ ε ρ , , μ 1 μ k ε ρ are mutually independent, the coefficients should vanish separately, implying the following N + 2 identities, which we shall refer to as δ L -identities below:
n = 0 N L j ; μ 1 μ n μ 1 μ n G ρ j ρ L = 0
K 0 ρ μ L δ ρ μ = 0
Sym { ν 1 ν k μ } K k ν 1 ν k μ ρ = 0 for k = 1 , 2 , , N .
Here, Sym { ν 1 ν k μ } implies the totally symmetric part with respect to the k + 1 indices ν 1 , , ν k and μ . Note that only the totally symmetric part of K k ν 1 ν k μ ρ should vanish since it vanishes when multiplied by the totally symmetric function ν 1 ν k μ ε ρ ( x ) . Note also that K k ν 1 ν k μ ρ is manifestly symmetric with respect to the first k indices ν 1 , , ν k , as is clear from the defining Equation (14).
Now, we can derive useful identities for rewriting the suitable linear combination of equations of motion (6)
δ S δ ϕ j = n = 0 N ( ) n μ 1 μ n L j ; μ 1 μ n .
(We call the quantities δ S / δ ϕ j ‘equations of motion’, although it is an abuse of terminology, since the equation of motion itself is the equation δ S / δ ϕ j = 0 .)
First, a linear combination ( δ S / δ ϕ j ) G ρ j of the equations of motion is rewritten into the following form by adding the first δ L -identity (15):
( δ S / δ ϕ j ) G ρ j = n = 0 N L j ; μ 1 μ n μ 1 μ n G ρ j ( ) n μ 1 μ n L j ; μ 1 μ n · G ρ j ρ L .
To rewrite this more concisely, we introduce a generalized ‘both-side’ derivative defined for n 0 by [28]:
F μ 1 μ 2 μ n G F μ 1 μ 2 μ n G μ 1 F · μ 2 μ n G + μ 1 μ 2 F · μ 3 μ n G + ( ) n μ 1 μ 2 μ n F · G ,
for two arbitrary functions F and G, with understanding μ 1 μ 2 μ n = 1 when n = 0 . This derivative is no longer symmetric under permutation of the indices but satisfies a useful formula [28]:
μ F μ α 1 α n α 1 α n G = F μ α 1 α n μ α 1 α n G + ( ) n μ α 1 α n F μ α 1 α n · G ,
for any totally symmetric function F μ α 1 α n with respect to the n + 1 indices { μ , α 1 , , α n } . Applying this formula, we can rewrite the identity (19) as
δ S δ ϕ j G ρ j = μ J ρ μ J ρ μ = n = 0 N 1 L j ; μ α 1 α n α 1 α n G ρ j δ ρ μ L .
This J ρ μ is the Noether current for the global GC transformation with x-independent ε ρ (=translation), i.e., energy-momentum tensor, for the higher derivative system. This identity shows that it is indeed conserved when the equations of motion δ S / δ ϕ j = 0 are satisfied.
Now, in order to derive various identities from the rest of the δ L -identities, (16) and (17), we need to introduce generalized both-side derivatives and some formulas for them.
We define k-th both-side derivative k by induction both in the number k and the differential order n:
α 1 α n k α 1 α n 1 k · α n + α 1 α n k 1 with initial condition
k = 1 : α 1 α n 1 = ( ) n α 1 α n for n 0 n = 0 : α 1 α n k | n = 0 = 1 for k 1 .
It is easy to see that the k = 0 is just the same as the original ‘both-side’ derivative ⟷ introduced above in Equation (20); indeed, it satisfies the above recursive defining relation (23) for k = 0 as follows:
α 1 α n = = 0 n ( ) α 1 α · α + 1 α n = = 0 n 1 ( ) α 1 α · α + 1 α n 1 α n + ( ) n α 1 α n = α 1 α n 1 · α n + α 1 α n k = 1 .
Then, as a generalization of the k = 0 formula (21), we have the following formula, which holds for all k 0 , n 0 and for any totally symmetric function F μ α 1 α n with respect to the n + 1 indices { μ , α 1 , , α n } :
μ F μ α 1 α n α 1 α n k G = F μ α 1 α n n + k + 1 C k μ α 1 α n μ α 1 α n k 1 G .
The proof easily goes by induction in the number N k + n in the region k 0 and n 0 . First, note that this formula holds at k = 0 boundary as shown above for n 0 and also clearly holds at n = 0 boundary with k 0 since the relevant k-th both-side derivatives there are α 1 α n k | n = 0 = 1 and that of a single derivative μ k 1 which is simply, by Equation (23),
μ k 1 = μ + μ k 2 = k μ + μ 1 = k μ μ .
Thus, it is sufficient to prove the formula only for k 1 and n 1 . Now, assume that formula (26) holds for all k 0 and n 0 values in the region k + n N , and let us evaluate the LHS of the formula for any k 1 and n 1 with k + n = N + 1 . If we use the defining Equation (23),
F μ α 1 α n α 1 α n k G = F μ α 1 α n α 1 α n 1 k ( α n G ) + F μ α 1 α n α 1 α n k 1 G ,
the two terms on the RHS have lower values k + n = N by one than the LHS, to which we can apply the formula by the induction assumption, so that
μ F μ α 1 α n α 1 α n k G = F μ α 1 α n n + k C k μ α 1 α n 1 μ α 1 α n 1 k 1 ( α n G ) + F μ α 1 α n n + k C k 1 μ α 1 α n μ α 1 α n k 2 G = F μ α 1 α n [ n + k C k + n + k C k 1 μ α 1 α n G μ α 1 α n 1 k 1 α n + μ α 1 α n k 2 G ] .
If we note an identity (of Pascal’s triangle) n + k C k + n + k C k 1 = n + k + 1 C k and again apply the defining Equation (23) with k k 1 , then we see that the last expression is simply reproducing the RHS of formula (26), finishing the proof.
Now, we are ready to derive the Maxwell-type form of gravity equation of motion. For this purpose, let us introduce the following quantity J k for k = 0 , 1 , 2 , :
J k ν 1 ν k μ ρ n = 0 N k 1 L j ; ν 1 ν k μ α 1 α n α 1 α n k G ρ j + n = 0 N k L j ; ν 1 ν k α 1 α n α 1 α n k 1 T ρ j μ .
The first of this quantity J 0 ρ μ with k = 0 is a combination of the equation of motion, ( δ S / δ ϕ j ) T ρ j μ , Lagrangian L , and the energy-momentum tensor J ρ μ :
J 0 ρ μ = J ρ μ + δ ρ μ L + δ S δ ϕ j T ρ j μ .
This can be seen from Equations (22) and (18) which are rewritten by using the definition of the k-th both-side derivative with k = 0 and 1 , respectively, into
n = 0 N 1 L j ; μ α 1 α n α 1 α n 0 G ρ j = J ρ μ + δ ρ μ L n = 0 N L j ; α 1 α n α 1 α n 1 T ρ j μ = δ S δ ϕ j T ρ j μ .
Owing to the general formula (26), the two quantities, J k introduced here (30) and K k defined previously in Equation (14), satisfy the following recurrence relation:
K k ν 1 ν k μ ρ J k ν 1 ν k μ ρ = τ J k + 1 τ ν 1 ν k μ ρ .
When applying the formula (26) to derive this equality, we should note that the summation over the set of n + 1 dummy indices { τ , α 1 α n } contained in the RHS quantity τ J k + 1 τ ν 1 ν k μ ρ is identified with the summation over the set { α 1 α n + 1 } contained in the LHS quantities K k ν 1 ν k μ ρ J k ν 1 ν k μ ρ by identifying α n + 1 as τ . This implies that the n = 0 terms existing in the summations n = 0 in K k and J k on the LHS do not appear on the RHS. However, the n = 0 terms in K k and J k are seen to be the same, thus canceling themselves on the LHS.
From this relation (33), we find, suppressing the tensor indices,
J k = K k J k + 1 = K k K k + 1 + 2 J k + 2 = = = 0 K ( ) K k + + ( ) K + 1 J k + K + 1 .
Since L ; ν 1 ν k = 0 for k > N , J k and K k vanish for k N + 1 . Thus, we find the following expression for J 0 ρ μ , reviving the tensor indices:
J 0 μ ρ = k = 0 N ( ) k ν 1 ν k K k ν 1 ν k μ ρ .
We now insert Equation (31) into J 0 on the LHS; then, noting that the term δ ρ μ L there cancels the K k = 0 term on the RHS due to the second δ L -identity (16), K 0 ρ μ = δ ρ μ L , we find the gravity field equation in the form
δ S δ ϕ j T ρ j μ = J ρ μ + k = 1 N ( ) k ν 1 ν k K k ν 1 ν k μ ρ .
This is still not the final form. The last summation term can be written as a divergence form of a ‘field-strength’ tensor F ρ ν μ , but it is not yet ν μ antisymmetric:
k = 1 N ( ) k ν 1 ν k K k ν 1 ν k μ ρ ν F ρ ν μ F ρ ν μ = k = 0 N 1 ( ) k ν 1 ν k K k + 1 ν 1 ν k ν μ ρ .
However, thanks to the remaining δ L -identities (17), we can modify it into an ν μ antisymmetric field strength F ˜ ρ ν μ , satisfying
ν F ρ ν μ = ν F ˜ ρ ν μ .
As noted before, the tensor K k ν 1 ν k μ ρ defined in Equation (14) is manifestly totally symmetric with respect to the first k-indices { ν 1 , , ν k } . The δ L -identities (17) say that it vanishes if further symmetrized, including the last index μ ; namely, taking the cyclic permutation of the k + 1 indices { ν 1 , , ν k , μ } ,
K k ν 1 ν k μ ρ + K k ν 2 ν k μ ν 1 ρ + K k ν 3 ν k μ ν 1 ν 2 ρ + + K k μ ν 1 ν k ρ = 0 .
If we apply k-ple divergence ν 1 ν k to this, the ν j indices become dummy, and, due to the manifest total symmetry among the first k indices of K k , the k terms from the second to the last yield the same quantity, and we obtain
ν 1 ν k K k ν 1 ν k μ ρ + k ν 1 ν k K k ν 2 ν k μ ν 1 ρ = 0 .
Or, taking k k + 1 and renaming ν k + 1 ν in the first term and ν 1 ν in the second term, we have
ν ν 1 ν k K k + 1 ν 1 ν k ν μ ρ + ( k + 1 ) ν 1 ν k K k + 1 ν 1 ν k μ ν ρ = 0 . .
This means that K k + 1 ν 1 ν k ν μ ρ can be replaced by a ν μ antisymmetric tensor, which we can define as
K ˜ k + 1 ν 1 ν k ν μ ρ = k + 1 k + 2 K k + 1 ν 1 ν k ν μ ρ K k + 1 ν 1 ν k μ ν ρ .
Indeed, the difference between K k + 1 and K ˜ k + 1 is given by
K k + 1 ν 1 ν k ν μ ρ K ˜ k + 1 ν 1 ν k ν μ ρ = 1 k + 2 K k + 1 ν 1 ν k ν μ ρ + ( k + 1 ) K k + 1 ν 1 ν k μ ν ρ ,
whose ( k + 1 ) -ple divergence ν ν 1 ν k is guaranteed to vanish by Equation (41). Thus, we find that the ‘field-strength’ F ρ ν μ in Equation (37) can be replaced by the ν μ antisymmetric one:
F ˜ ρ ν μ = k = 0 N 1 ( ) k ν 1 ν k K ˜ k + 1 ν 1 ν k ν μ ρ .
With this antisymmetric field strength, the gravity equation of motion is finally written in the desired form of the Maxwell-type equation:
δ S δ ϕ j T ρ j μ = J ρ μ ν F ˜ ρ μ ν .
This is an equation for the gauge-unfixed classical system.
Here, we note a more explicit expression for K ˜ k + 1 in terms of the Lagrangian. Substituting expression (14) for K k into Definition (42), we note that the G ρ j -proportional part contained in K k + 1 ν 1 ν k ν μ is ν μ symmetric so that only the T ρ j μ -proportional part contributes to K ˜ k + 1 , and obtain
K ˜ k + 1 ν 1 ν k ν μ ρ k + 1 k + 2 n = 0 N k 1 n + k + 1 C k + 1 L j ; ν 1 ν k ν α 1 α n α 1 α n T ρ j μ ( ν μ ) for k = 0 , 1 , 2 , .

3. Quantum Theory with de Donder Gauge

Let us now consider the quantum system. We add the gauge-fixing and corresponding Faddeev–Popov (FP) term to the classical GC invariant Lagrangian L cl . (We call the Lagrangian L in the previous section L cl hereafter.) We actually adopt Nakanishi’s simpler form of L GF + L FP = L GF + FP [4,5]:
L GF + FP = h δ N ( i κ 1 g μ ν μ c ¯ ν ) = κ 1 g ˜ μ ν μ b ν i g ˜ μ ν · μ c ¯ ρ · ν c ρ = δ B ( i κ 1 g ˜ μ ν μ c ¯ ν ) μ i g ˜ λ ν λ c ¯ ν · c μ ,
with h g and g ˜ μ ν h g μ ν . Here, the usual BRS transformation δ B (obtained by replacing ε μ ( x ) κ c μ ( x ) for the usual gravity/matter fields) is given by a sum of Nakanishi’s BRS δ N and the translation κ c λ λ :
δ B Φ = δ N Φ δ N ( x λ ) λ Φ , δ N ( x λ ) = κ c λ , δ N Φ = κ μ c α · Φ α μ , δ B g ˜ μ ν = κ λ c μ · g ˜ λ ν + λ c ν · g ˜ μ λ λ ( c λ g ˜ μ ν ) δ B c ¯ μ = i B μ , δ N c ¯ μ = i b μ , B μ = b μ + i κ c λ λ c ¯ μ , δ N c μ = 0 , δ B c μ = κ c λ λ c μ .
We call this gauge specified by the gauge-fixing and FP term L GF + FP in (47) “de Donder-Nakanishi gauge”. It corresponds to the de Donder–Landau gauge possessing no α η μ ν B μ B ν term, violating G L ( 4 ) invariance by the use of η μ ν . Since the present L GF + FP for the de Donder–Nakanishi gauge is given in a usual BRS exact form for the de Donder–Landau gauge up to a total derivative term as shown in the last expression in Equation (47), it is also invariant under the usual BRS transformation δ B . The use of Nakanishi’s BRS δ N , which represents the tensorial transformation part of the usual BRS transformation δ B , and the use of the b μ = i 1 δ N c ¯ μ field, in particular, make manifest the existence of much larger I O S p ( 8 | 8 ) symmetry, called choral symmetry by Nakanishi, which contains symmetries of energy-momentum, G L ( 4 ) , BRS, FP-ghost scale transformation, etc., as will be discussed explicitly in the next section.
We still consider the GC transformation x ρ x ρ = x ρ ε ρ ( x ) in this quantum theory to derive identities. The gravity/matter fields ϕ j are transformed in the same way as before:
δ ϕ j ( x ) = G ρ j ε ρ ( x ) + T ρ j μ μ ε ρ ( x ) , with G ρ j = ρ ϕ j , T ρ j μ = ϕ j ρ μ .
We call the newly added fields b μ , c ¯ μ and c μ ghost fields collectively, and treat them all as scalar fields under GC transformation, namely denoting ghost fields by ϕ M = ( b μ , c ¯ μ , c μ ) collectively:
δ ϕ M ( x ) = G ρ M ε ρ ( x ) + T ρ M μ μ ε ρ ( x ) , with G ρ M = ρ ϕ M , T ρ M μ = 0 .
Of course, L GF + FP is not invariant under the GC transformation, but we can easily calculate the change by noting the structure of the L GF + FP , which is written formally as a scalar density:
L GF + FP = κ 1 g ˜ μ ν E μ ν , E μ ν = μ b ν + i κ · μ c ¯ ρ · ν c ρ .
If the ghost part tensor E μ ν truly behaved as a μ ν covariant tensor, L GF + FP would be a scalar density transforming only into the total divergence μ ( L GF + FP ε μ ) . This is actually true for the FP ghost part i κ · μ c ¯ ρ · ν c ρ in E μ ν since c ρ and c ¯ ρ are regarded as scalars, so that their simple derivatives μ c ¯ ρ and ν c ρ behave as μ and ν vectors, giving the desired μ ν tensor as a product. But the NL field part μ b ν transforms only as a μ vector since b ν is regarded as a scalar, so that the ν leg rotation part of the transformation of g ˜ μ ν , i.e., δ g ˜ μ ν g ˜ μ ρ ρ ε ν , is not canceled. We thus see
δ L GF + FP = κ 1 g ˜ μ ρ μ b ν · ρ ε ν + μ ( L GF + FP ε μ ) .
Thus, the total Lagrangian in our quantum gravity theory
L = L cl + L GF + FP
changes under the GC transformation as
δ L = μ ( L ε μ ) + κ 1 g ˜ μ ρ μ b ν · ρ ε ν = ρ L · ε ρ + ( L δ ρ μ + κ 1 g ˜ μ ν ν b ρ ) μ ε ρ .
This namely differs from Equation (11) in the classical system case, only in the point that the κ 1 g ˜ μ ν ν b ρ term is added in the first order derivative term μ ε ρ . Therefore, the δ L -identities in the previous section almost all remain the same and only the first order μ ε ρ identity (16) is slightly changed into
K 0 ρ μ ( L δ ρ μ + κ 1 g ˜ μ ν ν b ρ ) = 0 .
Note that we should now understand that L is the total Lagrangian containing the ghost part L GF + FP , and the fields ϕ j cover not only the gravity/matter fields ϕ j but also the ghost fields ϕ M = ( b μ , c ¯ μ , c μ ) . The equation of motion δ S / δ ϕ , of course, takes the same form (18) as before. The zeroth order δ L -identity (15), in particular, remains the same, and the global translation current (energy-momentum tensor) is given by the same form of equation as Equation (22):
δ S δ ϕ j G ρ j = μ J ρ μ J ρ μ = n = 0 N 1 L j ; μ α 1 α n α 1 α n G ρ j δ ρ μ L .
Thus, Equation (31) for J 0 holds unchanged. The identity (35) also holds as it stands. In going from Equation (35) to the gravity Equation (36), however, the L δ ρ μ term from J 0 now does not totally cancel the first term K k = 0 ρ μ = L δ ρ μ + κ 1 g ˜ μ ν ν b ρ but leaves the κ 1 g ˜ μ ν ν b ρ term. Thus, Equation (36) is now replaced by
δ S δ ϕ j T ρ j μ = J ρ μ + κ 1 g ˜ μ ν ν b ρ + k = 1 N ( ) k ν 1 ν k K k ν 1 ν k μ ρ .
Note, here, that the implicit summation over ϕ j also contains the ghost fields ϕ M which contribute only to the n = 0 terms since the ghost fields appear only in the first order derivatives in the de Donder–Nakanishi gauge Lagrangian (47).
The final form of the Maxwell-type gravity field equation is therefore given by
δ S δ ϕ j T ρ j μ = J ρ μ κ 1 g ˜ μ ν ν b ρ ν F ˜ ρ μ ν .
in place of the previous classical one (45). The expressions in Equation (44) for the field-strength F ˜ μ ν ρ and Equation (46) for the quantities K ˜ k + 1 ν 1 ν k ν μ ρ remain the same as before. Here, L is understood to be the total Lagrangian, but only the classical Lagrangian part actually contributes in this case, because all the ghost fields ϕ M = { b μ , c ¯ μ , c μ } have vanishing contributions since T M ρ μ = 0 for them. That is, the field-strength is in fact the same as that in the classical theory with Lagrangian L cl .
One may wonder why the final Maxwell-type gravity equation of motion (58) is slightly different from the Yang–Mills case, since the present ghost field term κ 1 g ˜ μ ν ν b ρ is not written in a BRS exact form such as { Q B , D μ c ¯ } in the latter. It is actually possible to rewrite Equation (58) into such a form. Indeed, the term κ 1 g ˜ μ ν ν b ρ is in fact BRS exact up to a divergence of an antisymmetric tensor:
κ 1 g ˜ μ ν ν b ρ = δ B i κ 1 g ˜ μ ν ν c ¯ ρ ν i ( c μ g ˜ ν σ c ν g ˜ μ σ ) σ c ¯ ρ .
The gravity field Equation (58), therefore, can be rewritten into quite a similar form as the Maxwell-type YM equation:
ν F ˜ ˜ ρ μ ν + Q B , κ 1 g ˜ μ ν ν c ¯ ρ = J ρ μ ,
where we have written δ B ( ) = { i Q B , } in terms of the BRS charge Q B and defined a modified field strength F ˜ ˜ ρ μ ν :
F ˜ ˜ ρ μ ν = F ˜ ρ μ ν + i ( c μ g ˜ ν σ c ν g ˜ μ σ ) σ c ¯ ρ .
This form of the Maxwell-type equation (60) with the BRS exact term was also derived for the Einstein theory case by Nakanishi [12].

4. Noether Current for the Choral Symmetries in a Generic Higher Derivative System

BRS symmetry or, more generally, choral symmetries I O S p ( 8 | 8 ) exist for any GC transformation invariant systems if one adopts the gauge-fixing Lagrangian (47) of the de Donder–Nakanishi gauge [11]. This is because the currents of the choral symmetries are conserved as far as the equations of motion
μ ( g ˜ μ ν ν X M ) = 0
hold for the 16 component ‘fields’ (=4d coordinate x μ and three fields) [9,10].
X M = ( x ^ μ , b μ , c μ , c ¯ μ ) , x ^ μ x μ / κ .
Indeed, this equation of motion for the coordinate X M x ρ actually implies the de Donder condition on the gravity field:
μ ( g ˜ μ ν ν x ρ ) = μ g ˜ μ ρ = 0 .
In addition, the FP ghost equations of motion
L FP c ¯ ρ = i μ ( g ˜ μ ν ν c ρ ) = 0 , L FP c ρ = + i μ ( g ˜ μ ν ν c ¯ ρ ) = 0
directly follow from the gauge-fixing Lagrangian (47), implying the equations for X M = c ρ and c ¯ ρ . The equation for X M = b ρ may be a bit non-trivial, but we now already know the Maxwell-type gravity equation of motion (58), the divergence μ of which immediately leads to
μ ( g ˜ μ ν ν b ρ ) = 0 .
These 16 components’ d’Alembert’s equations of motion hold if and only if the gauge-fixing Lagrangian L GF + FP is given by that of the de Donder–Nakanishi (47), which can be written in a manifestly O S p ( 8 | 8 ) invariant form. (The following discussion on the O S p ( 8 | 8 ) invariance may be viewed as a mere recapitulation of Nakanishi’s paper [10,11], but we have simplified and made in particular the signs and i factors more tractable by introducing a hermitian O S p ( 8 | 8 ) metric (68). The derivation of the O S p ( 8 | 8 ) Noether current in the higher derivative system is, of course, new.)
L GF + FP = κ 1 h E = κ 1 h g μ ν E μ ν , E μ ν = 1 2 μ b ν + i κ μ c ¯ ρ · ν c ρ + ( μ ν ) = 1 2 μ b ρ · ν x ρ + i κ μ c ¯ ρ · ν c ρ + ( μ ν ) = κ 2 η N M μ X M ν X N .
where η M N is the O S p ( 8 | 8 ) metric given by
η N M = δ μ ν δ ν μ i δ μ ν i δ ν μ = η N M ( inverse ) .
Note the symmetry property of this (c-number) metric:
η M N = ( ) | M | · | N | η N M = ( ) | M | η N M = ( ) | N | η N M η ˜ N M ,
where the statistics index | M | is 0 or 1 when X M is bosonic or fermionic, respectively. This property (69) is because η M N is ‘diagonal’ in the sense that its off-diagonal, bose-fermi, and fermi-bose matrix elements vanish, i.e., η M N = 0 when | M | | N | , so that | M | = | N | = | M | · | N | in front of η M N . Note also that the η ˜ N M introduced here is simply the transposed metric η N M T = η M N . Thus, we have
κ 2 η N M X M X N = κ 2 X M η ˜ M N X N = κ 2 x ^ μ b μ c μ c ¯ μ δ μ ν δ ν μ i δ μ ν + i δ ν μ x ^ ν b ν c ν c ¯ ν = κ 2 x ^ μ b μ + b μ x ^ μ i c μ c ¯ μ + i c ¯ μ c μ = b μ x μ + i κ c ¯ μ c μ .
Noting the d’Alembert’s equations of motion for X M , Nakanishi constructed the conserved currents.
M M N μ g ˜ μ ν X M ν X N
P M μ g ˜ μ ν ν X M = g ˜ μ ν 1 ν X M .
(Our current M M N μ presented here is not exactly equal to Nakanishi’s original one [11], M μ ( X M , X N ) , but the precise relation reads M μ ( X M , X N ) = i | M | · | N | M M N μ . )
He showed from the equal-time commutation relations (ETCR) derived in the Einstein gravity theory that their charge operators
M M N d 3 x M M N 0 = ( ) 1 + | M | · | N | M N M , P M d 3 x P M 0
generate the following transformations on all the fields Φ , gravity and matter fields ϕ j , and the O S p ( 8 | 8 ) ghost-‘fields’ X M = ( x ^ μ , b μ , c μ , c ¯ μ ) :
[ i M M N , Φ } δ M N Φ = δ N M N Φ κ ( δ N M N x ^ ρ ) ρ Φ ,
where the Nakanishi transformation δ N M N is an O S p rotation for the O S p ( 8 | 8 ) ghost-fields X L , given by
δ N M N X L = η ˜ N L X M + ( ) | M | · | N | η ˜ M L X N .
This transformation, in particular, gives for the coordinate x ^ ρ
δ N M N x ^ ρ = η ˜ N x ^ ρ X M + η ˜ M x ^ ρ X N ( κ ) 1 E M N ρ ,
which is nonvanishing only when X M or/and X N is b μ . Additionally, the Nakanishi transformation of the gravity/matter fields ϕ j is given by
δ N M N ϕ j = μ E M N ν · ϕ j ν μ .
Therefore, if either X M or X N equals b μ , the δ M N transformation is just the GC transformation with transformation parameter ε ρ ( x ) E M N ρ = κ ( δ N M N x ^ ρ ) for the gravity/matter fields ϕ j ,
δ M N ϕ j = E M N ρ G ρ j + μ E M N ρ · T ρ j μ ( G ρ j = ρ ϕ j , T ρ j μ = ϕ j ρ μ ) , with field dependent parameter E M N ρ = κ ( δ N M N x ^ ρ ) ,
and, for the O S p ( 8 | 8 ) ghost-fields X L , the GC transformation as scalar fields plus an O S p rotation:
δ M N X L = E M N ρ G ρ L ( η ˜ N L X M ( ) | M | · | N | η ˜ M L X N ) ( G ρ L = ρ X L ) .
Note, here, that the ‘transformation parameter’ E M N ρ may now be fermionic when | M | + | N | = 1 . We have therefore put the factors G ρ j and T ρ j μ linear in ϕ j behind the parameter E M N ρ in Equations (77) and (78) to avoid the sign factor ( 1 ) | ϕ j | ( | M | + | N | ) .
The choral invariance of our total Lagrangian L = L cl + L GF + FP is now clear; the gravity/matter fields only receive a special GC transformation with parameter ε ρ ( x ) = E M N ρ in Equation (77), and so the L cl part is invariant. The gauge-fixing Lagrangian L GF + FP is also clearly invariant since in this form of the Lagrangian (67), the O S p -vector field components X M , including the coordinate x ^ ρ , are treated as scalar fields, and hence E μ ν is clearly μ ν tensor, and h g μ ν E μ ν is manifestly a GC scalar density. Moreover, E μ ν written in the form (67) is manifestly invariant under (global) O S p ( 8 | 8 ) rotation. Note that this invariance is made manifest by making the mere parameter coordinate x μ transform as if being a field both under the GC transformation and the O S p rotation; in fact, those two transformations on x ^ μ cancel each other, and the coordinate x μ remains intact under δ M N as any non-field parameters should be: indeed, Equation (73) indicates for Φ = x ^ μ ,
δ M N x ^ μ = δ N M N x ^ μ κ ( δ N M N x ^ ρ ) ρ x ^ μ = 0 ,
where the first term is the O S p rotation, and the second term is the GC transformation of x ^ μ regarded as a ‘scalar field’.
Let us now compute the Noether currents corresponding to these choral symmetries in our general higher derivative GC invariant system. We shall show that the Noether currents coincide with the Nakanishi’s simple form (70), aside from the divergence of an antisymmetric tensor.
To do this systematically, we devise a local version of the choral symmetry transformation (73), or (77) and (78). We multiply them by a local graded transformation parameter ε N M ( x ) from the left, so that it reduces to the original O S p ( 8 | 8 ) transformation in the global limit ε N M ( x ) ε N M : const . ; namely, we define the transformation,
δ ϕ j = ( ε N M E M N ρ ) G ρ j + μ ( ε N M E M N ρ ) · T ρ j μ for gravity / matter fields ϕ j , δ X L = ( ε N M E M N ρ ) G ρ L ε N M ( η ˜ N L X M ( ) | M | · | N | η ˜ M L X N ) for O S p coordinate fields X L .
We take our parameter ε N M Grassmann even or odd according to | N | + | M | = 0 or 1, respectively, so that the product ( ε N M E M N ρ ) always becomes an ordinary bosonic ‘parameter’ and hence can be moved to anywhere without worrying about sign changes. Note, however, that, in order to correctly obtain the Noether current corresponding to the O S p transformation δ M N in Equation (73), or Equations (77) and (78), we have to factor out the parameter ε N M from the left since it is multiplied from the left here. However, the general procedure explained in the previous sections to derive the Noether current in the higher derivative theories, which we follow now, has placed the transformation parameter at the most right end, and the troublesome point is that the transformation parameter for the O S p transformation is the graded one, ε N M , but not the bosonic product ( ε N M E M N ρ ) . It is necessary to move these graded quantities separately and freely to apply the general procedure to this case, although the transformation parameter ε N M eventually has to be factored out from the left. The best way to forget about the bothering sign factors appearing in changing the order of the graded quantities is to adopt a convention similar to the so-called ‘implicit grading’ [29]. We take as a natural order of those graded quantities, ε N M first, E M N ρ second, and other graded quantities such as X M , G ρ J , and T j m u ρ third. Initially, these quantities appear in this natural order, since the product factor ( ε N M E M N ρ ) appearing in Equation (80) is bosonic and can be placed on the most left in any case. Then, from this natural order, we freely move those factors separately anywhere without writing any sign factors. The implicit grading scheme means that the correct sign factors should be recovered when necessary; that is, in any terms containing those graded quantities, the necessary sign factor can be found by counting how many instances of changing order are necessary to bring those factors into the natural order. We hereafter adopt this implicit grading scheme.
We should note that the GC transformation part of this transformation (80) now takes exactly the same form as the GC transformation (8) with the (bosonic) transformation parameter ε ρ ( x ) ε N M E M N ρ :
δ Φ j = G ρ j ε ρ + T ρ j μ μ ε ρ
for all the fields Φ j = ( ϕ j , X M ) , and so the total action is still invariant, meaning that the total Lagrangian transforms as a scalar density: δ L = μ ( L ε μ ) = μ ( L ε N M E M N ρ ) .
As for the remaining O S p ( 8 | 8 ) rotation part of the O S p coordinate ‘fields’ X L ,
δ ^ O S p X L ε N M η ˜ N L X M ( ) | M | · | N | η ˜ M L X N ,
However, the L GF + FP in Equation (67) is no longer invariant under the rotation with the x-dependent parameter ε M N . As shown explicitly in Appendix A, we can immediately find the change of L GF + FP as Equation (A6):
δ L GF + FP = g ˜ μ ν μ ε N M · X M ν X N = g ˜ μ ν X M ν X N μ ε N M ,
where it should be noted that we have already used ‘implicit grading’ in the last equality. Following this implicit grading scheme, we can write the change of our total Lagrangian L = L cl + L GF + FP under our transformation (80) in the form
δ L = μ ( L E M N μ ε N M ) + g ˜ μ ν X M ν X N · μ ε N M = μ ( L E M N μ ) · ε N M + L E M N μ + g ˜ μ ν X M ν X N μ ε N M .
Now, we can rewrite our transformation in the same form as the general gauge transformation (8), which contains the zero-th and first order differentiation of the transformation parameter ε N M ; that is, it is unifiedly given for gravity/matter and ghost fields Φ I = ( ϕ j , X L ) in the form
δ Φ I = G I M N ε N M + T I M N μ μ ε N M ,
where the coefficients G I M N and T I M N μ are given as
G I M N = G ρ j E M N ρ + T ρ j μ μ E M N ρ for Φ I = ϕ j : gravity / matter G ρ L E M N ρ η ˜ N L X M ( ) | M | · | N | η ˜ M L X N for Φ I = X L : ghost T I M N μ = T ρ j μ E M N ρ for Φ I = ϕ j : gravity / matter 0 for Φ I = X L : ghost .
Thus, we can now follow the general discussions presented in the previous two sections to derive the Noether currents in this system. We should also note that the Lagrangian change δ L is now given by Equation (84) in place of Equation (11). Then, we see that the previous δ L -identities following from the coefficients of n-th order derivatives μ 1 μ n ε ρ of the transformation parameter ε ρ , Equations (15) and (17), now also hold with the understanding that the coefficients G ρ and T ρ μ of the gauge transformation (8) are now replaced by those of the present transformation: that is, by making the replacement
G ρ j G I M N , T ρ j μ T I M N μ ,
and taking account of the form of δ L given in Equation (84), we find the identities as coefficients of n-th order derivatives μ 1 μ n ε N M with n = 0 , 1 and n 2 , respectively,
n = 0 N L I ; μ 1 μ n μ 1 μ n G I M N ρ ( L E M N ρ ) = 0
K 0 M N μ L E M N μ g ˜ μ ν X M ν X N = 0
Sym { ν 1 ν k μ } K k M N ν 1 ν k μ = 0 for k = 1 , 2 , ,
where K k M N ν 1 ν k μ now reads
K k M N ν 1 ν k μ n = 0 N k 1 n + k + 1 C k + 1 L I ; ν 1 ν k μ α 1 α n α 1 α n G I M N + n = 0 N k n + k C k L I ; ν 1 ν k α 1 α n α 1 α n T I M N μ .
Then, by combining the equation of motion (6), we first obtain a conservation equation for the O S p ( 8 | 8 ) Noether current from the 0-th order δ L -identity (88) as an analogue of Equation (22) or (56):
δ S δ Φ I G I M N = μ J M N μ J M N μ = n = 0 N 1 L I ; μ α 1 α n α 1 α n G I M N L E M N μ .
Second, as an analogue of Equation (45) or Equation (58), we find the desired equation from the first order δ L -identity (89) and second or higher order δ L -identity (90):
δ S δ Φ I T I M N μ = J M N μ g ˜ μ ν X M ν X N ν F ˜ M N μ ν .
where F ˜ M N ν μ is the antisymmetric tensor (as an ambiguity term of O S p Noether current) given by
F ˜ M N ν μ = k = 0 N 1 ( ) k ν 1 ν k K ˜ k + 1 M N ν 1 ν k ν μ
with
K ˜ k + 1 M N ν 1 ν k ν μ k + 1 k + 2 n = 0 N k 1 n + k + 1 C k + 1 L I ; ν 1 ν k ν α 1 α n α 1 α n T I M N μ ( ν μ ) for k = 0 , 1 , 2 , .
Equation (93) shows that the Noether current of O S p ( 8 | 8 ) symmetry takes the form on-shell:
J M N μ = g ˜ μ ν X M ν X N + ν F ˜ M N μ ν ,
so that the charge can be given by Nakanishi’s simple form:
M M N = d 3 x h g 0 ν X M ν X N
as symmetry generators for any local operators. (It is, however, quite another problem whether the RHS volume integral of (97) gives a well-defined charge even when the symmetry is spontaneously unbroken. Generally, the volume integral converges only when the current contains the ambiguity term ν F ˜ M N μ ν with a suitable coefficient, which is, in fact, a key point when discussing the spontaneous breaking or non-breaking of the symmetry of the charge.)
Finally, recall that we are adopting implicit grading. Then, the factor E M N μ or its derivatives contained in G I M N or T I M N μ in the definitions of the current J N M μ (92), the field strength F ˜ M N ν μ (94), and K ˜ k + 1 (95) should be placed at the furthest left, and if they are kept in the place as written, then the grading sign factor should be set as necessary for bringing them there from the furthest left. Fortunately, however, the sign factors actually turn out to be unnecessary here. This is because the other graded quantities to jump over when bringing the factor E M N μ to the furthest left are essentially only the bosonic Lagrangian L as a whole. For instance, the Noether current (92) contains the terms
L I ; μ α 1 α n α 1 α n G I M N L E M N μ .
In the second term, it is the Lagrangian L itself for E M N μ to jump over. In the first term, on the other hand, L I ; μ α 1 α n = L / μ α 1 α n ϕ I is fermionic when the field ϕ I is a fermion. However, it is immediately followed by the factor G ρ I = ρ ϕ I in the G I M N = G ρ I E M N ρ , so that the net factor L I ; μ α 1 α n G ρ I in front of E M N ρ is always bosonic, carrying the same statistics as the original L .

5. Conclusions

In this paper, we considered the general GC invariant theory which contains arbitrarily high order derivative fields. We identified the explicit expression for the energy-momentum as the Noether current corresponding to the rigid case of GC transformation and showed that a linear combination of the equations of motion can be rewritten into the form of a Maxwell-type field equation, which has the total energy-momentum as its source. This was done both for the gauge-unfixed classical system and for the gauge-fixed quantum system in the de Donder–Nakanishi gauge. The Maxwell-type field equation in the latter formally has an additional term coming from the gauge-fixing compared with that in the former, which was shown to take a BRS-exact form, just as in the well-known Yang–Mills case.
By using the same technique, we derived similar expressions of Maxwell-type equations for the Noether currents for the I O S p ( 8 | 8 ) choral symmetries. It confirmed that Nakanishi’s original result persists to any GC invariant system.
These results will be useful for proving the existence theorem [27] of the graviton (photon) in any GC (local U ( 1 ) ) gauge invariant system as far as the translation (global U ( 1 ) ) symmetry is not spontaneously broken. We hope we can report on this matter in the near future.
The techniques used for the higher derivative systems in this paper will also be useful for studying more general problems. For instance, basic problems such as the (non-)equivalence between canonical quantization and path-integral quantization may be discussed in a general fashion for the higher derivative systems. This is because the present technique can easily be combined with Ostrogradsky’s canonical formalism for higher derivative systems [30].
[Note added in proof] We are informed by the author of Ref. [31] in which similar results to ours were obtained in arbitrary higher derivative gravity systems. We see that he essentially proved the gravity equations of motion taking the form of Maxwell-type equations in the case of gauge-unfixed classical system, although the detailed comparison of equations is not easy since rather different notations are used. His assumption on the gauge transformation is, however, more restrictive than ours to prove the anti-symmetric property of the field strength (which he calls quasilocality).

Funding

This work is supported in part by the JSPS KAKENHI Grant Number JP18K03659.

Acknowledgments

The author would like to thank Ichiro Oda for valuable discussions on Weyl invariant gravity theories and reading the early version of this manuscript.

Conflicts of Interest

The author declares no conflict of interest.

Appendix A. OSp Transformation of GF + FP

Note that the parameter ε N M and the metric η N M = η N M have the following statistics and transposition properties:
| ε N M | = | M | + | N | , i . e . , ε N M X L = ( ) | L | ( | M | + | N | ) X L ε N M ε N M = ( ε T ) M N = ( ) 1 + | M | · | N | ε M N η N M = ( η T ) M N = ( ) | M | η M N = ( ) | N | η M N η ˜ M N .
Using these properties and also noting that the metric η M N is a usual c-number (i.e., bosonic) quantity, we can rewrite the two terms in the O S p transformation (82) of the O S p coordinates X L as
ε N M η ˜ N L X M = η L N ε N M X M = ( η 1 ε X ) L + ε N M ( ) | M | · | N | η ˜ M L X N = ( ) | M | · | N | η L M ( ) 1 + | M | · | N | ε M N X N = η L M ε M N X N . = ( η 1 ε X ) L ,
where X, η 1 , and ε are regarded as a column vector and matrices in the final expressions. Thus, the two terms reduce to the same expression, and so the O S p transformation (82) is rewritten concisely into
δ O S p X M = 2 ( η 1 ε X ) M = 2 η M N ε N L X L .
Now, it is easy to see how the gauge-fixing Lagrangian L GF + FP = κ 1 g ˜ μ ν E μ ν changes under the O S p transformation (82). Noting the expression of E μ ν in Equation (67), we find
δ O S p L GF + FP = δ O S p ( κ 1 g ˜ μ ν E μ ν ) = 1 2 g ˜ μ ν η N M δ O S p ( μ X M ν X N ) = g ˜ μ ν η N M ( δ O S p μ X M ) ν X N = g ˜ μ ν η N M μ ( 2 η 1 ε X ) M ν X N = 2 g ˜ μ ν μ ( ε N M X M ) ν X N = 2 g ˜ μ ν μ ε N M · X M ν X N + ε N M μ X M ν X N .
Now, we note that the O S p transformation parameter ε N M must be graded antisymmetric
ε N M = ( ) | N | · | M | ε M N
in order for the Lagrangian to be global O S p invariant. Indeed, only the second term remains in the global transformation, and it vanishes if and only if ε N M is graded antisymmetric since μ X M ν X N , multiplied by g ˜ μ ν , is graded symmetric. The second term thus also vanishes here, and the remaining first term is rewritten into
δ O S p L GF + FP = 2 g ˜ μ ν μ ε N M · X M ν X N = g ˜ μ ν μ ε N M · X M ν X N ( ) | M | · | N | X N ν X M = g ˜ μ ν μ ε N M · X M ν X N ν X M · X N = g ˜ μ ν μ ε N M · X M ν X N ,
where, in going to the second line, we have exchanged the dummy indices N M and used the graded antisymmetry property of ε N M , (A5).

Appendix B. OSp(2|2)-Invariant Scalar Field System

To see the property of the O S p symmetry, let us consider here the O S p -invariant system on a flat Minkowski background in which scalar fields belong to an O S p ( 2 | 2 ) vector representation.
L = 1 2 η N M μ ϕ M μ ϕ N ,
where ϕ M is a 2 + 2 -component O S p -vector whose first 2-components are bosons and the remaining 2-components are fermions, and the O S p ( 2 | 2 ) metric is given by
η N M = σ 1 σ 2 = η N M .
Note that we have added an overall minus sign to our Lagrangian (A7) in order to make it coincide with the convention of the gauge-fixing Lagrangian L GF + FP in Equation (67), although it is not physically important, in any case, since the O S p metric is neither positive- nor negative-definite.
The infinitesimal O S p rotation is parametrized as
δ ϕ M = ε N M ϕ N .
Under this rotation, the quadratic kinetic Lagrangian (A7) is transformed as
δ L = 1 2 η N M ε L M μ ϕ L μ ϕ N + μ ϕ M ε L N μ ϕ L = 1 2 η N M ε L M μ ϕ L μ ϕ N + μ ϕ M η N M ε L N μ ϕ L = 1 2 ε N L μ ϕ L μ ϕ N + μ ϕ M ( ) | M | ε M L μ ϕ L = 1 2 ε N L μ ϕ L μ ϕ N + ε M L μ ϕ L μ ϕ M = ε N M μ ϕ M μ ϕ N
with ε N M η N L ε M L . This can further be rewritten as
ε N M μ ϕ M μ ϕ N = ε M N μ ϕ N μ ϕ M = ( ) | M | · | N | ε M N μ ϕ M μ ϕ N = 1 2 ( ) | M | · | N | ε M N + ε N M μ ϕ M μ ϕ N .
Thus, if the transformation parameter ε N M is graded antisymmetric, i.e.,
ε N M = ( ) | M | · | N | ε M N ,
then the Lagrangian (A7) is O S p -invariant.
We define the canonical conjugate variable π M by the right-derivative
π M L / ϕ ˙ M = η M N ϕ ˙ N .
Then, the ETCR is given by
[ ϕ M , π N } = i δ N M [ ϕ M , ϕ ˙ N } = η N L [ ϕ M , π L } = η N L i δ L M = i η N M ,
so that
[ ϕ M , ϕ ˙ N } = i ( ) N η M N i η ˜ M N , [ ϕ ˙ M , ϕ N } = i ( ) N η M N = i η ˜ M N .
(If the conjugate momentum π M were defined by the left derivative, then the ETCR should be given by [ π M , ϕ N } = i δ M N .)
The Noether current J M N μ for the O S p ( 2 | 2 ) transformation is defined by
1 2 ε N M J M N μ L / ( μ ϕ M ) δ ϕ M = η M N μ ϕ N ε L M ϕ L = ( ) N μ ϕ N ε N L ϕ L = ε N L ϕ L μ ϕ N = ε N M 1 2 ( ϕ M μ ϕ N ( ) | N | · | M | ϕ N μ ϕ M ) = ε N M 1 2 ( ϕ M μ ϕ N μ ϕ M · ϕ N ) = 1 2 ε N M ( ϕ M μ ϕ N ) ,
where we have used the graded-antisymmetry property of ε N M in going from the second to third lines. Thus, we have
J M N μ = ϕ M μ ϕ N ,
so that the charge is given by
M M N = d 3 x J M N 0 = d 3 x ( ϕ M ϕ ˙ N ϕ ˙ M ϕ N ) .
This actually generates the original O S p rotation by the ETCRs:
[ i M M N , ϕ L } = d 3 x i ϕ M [ ϕ ˙ N , ϕ L } i ( ) | N | · | L | [ ϕ ˙ M , ϕ L } ϕ N = ϕ M η ˜ N L + ( ) | N | · | L | η ˜ M L ϕ N = ϕ M η ˜ N L + ( ) | N | · | M | ϕ N η ˜ M L ,
[ i M M N , ϕ ˙ L } = d 3 x i ( ) | N | · | L | [ ϕ M , ϕ ˙ L } ϕ ˙ N i ϕ ˙ M [ ϕ N , ϕ ˙ L } = ( ) | N | · | L | η ˜ M L ϕ ˙ N ϕ ˙ M η ˜ N L = ϕ ˙ M η ˜ N L + ( ) | N | · | M | ϕ ˙ N η ˜ M L .
The I O S p ( 2 | 2 ) algebra is confirmed as
[ i M M N , M R S } = [ i M M N , d 3 x ( ϕ R ϕ ˙ S ϕ ˙ R ϕ S ) } = d 3 x [ ϕ M η ˜ N R ϕ ˙ S + ϕ ˙ M η ˜ N R ϕ S + ( ) ( | M | + | N | ) | R | ϕ R ϕ ˙ M η ˜ N S + ϕ ˙ R ϕ M η ˜ N S ( ) | M | · | N | M N ] = M M S η ˜ N R + ( ) | R | · | S | M M R η ˜ N S ( ) | M | · | N | M N S η ˜ M R + ( ) | R | · | S | M N R η ˜ M S = M M S η ˜ N R + ( 3 graded antisymmetrization terms under ( M N ) and ( R S ) ) ,
[ i M M N , P R } = [ i M M N , d 3 x ( ϕ ˙ R ) } = d 3 x ϕ ˙ M η ˜ N R + ( ) | M | · | N | M N = P M η ˜ N R + ( ) | M | · | N | P N η ˜ M R .

References

  1. Ojima, I. Observables and Quark Confinement in the Covariant Canonical Formalism of the Yang-Mills Theory. Nucl. Phys. B 1978, 143, 340–352. [Google Scholar] [CrossRef] [Green Version]
  2. Kugo, T. The Gauge Invariance and a Possible Mechanism of Color Confinement. Phys. Lett. B 1979, 83, 93–97. [Google Scholar] [CrossRef]
  3. Kugo, T.; Ojima, I. Local Covariant Operator Formalism of Nonabelian Gauge Theories and Quark Confinement Problem. Prog. Theor. Phys. Suppl. 1979, 66, 1–130. [Google Scholar] [CrossRef] [Green Version]
  4. Nakanishi, N. Indefinite-Metric Quantum Field Theory of General Relativity. Prog. Theor. Phys. 1978, 59, 972. [Google Scholar] [CrossRef] [Green Version]
  5. Nakanishi, N. Remarks on the Indefinite Metric Quantum Field Theory of General Relativity. Prog. Theor. Phys. 1978, 59, 2175. [Google Scholar] [CrossRef] [Green Version]
  6. Nakanishi, N. Indefinite Metric Quantum Field Theory of General Relativity. 2. Commutation Relations. Prog. Theor. Phys. 1978, 60, 1190. [Google Scholar] [CrossRef] [Green Version]
  7. Nakanishi, N. Indefinite Metric Quantum Field Theory of General Relativity. 3. Poincare Generators. Prog. Theor. Phys. 1978, 60, 1890. [Google Scholar] [CrossRef] [Green Version]
  8. Nakanishi, N. Indefinite Metric Quantum Field Theory of General Relativity. 5. Vierbein Formalism. Prog. Theor. Phys. 1979, 62, 779. [Google Scholar] [CrossRef] [Green Version]
  9. Nakanishi, N. Indefinite Metric Quantum Field Theory of General Relativity. 9. ‘Choral’ of Symmetries. Prog. Theor. Phys. 1980, 63, 2078. [Google Scholar] [CrossRef] [Green Version]
  10. Nakanishi, N. Indefinite Metric Quantum Field Theory of General Relativity. 10. Sixteen-dimensional Superspace. Prog. Theor. Phys. 1980, 64, 639. [Google Scholar] [CrossRef] [Green Version]
  11. Nakanishi, N.; Ojima, I. Covariant operator formalism of gauge theories and quantum gravity. World Sci. Lect. Notes Phys. 1990, 27, 1–434. [Google Scholar]
  12. Nakanishi, N. Indefinite Metric Quantum Field Theory of General Relativity. 14. Sixteen-dimensional Noether Supercurrents and General Linear Invariance. Prog. Theor. Phys. 1981, 66, 1843. [Google Scholar] [CrossRef] [Green Version]
  13. Ogievetsky, V.I. Nonlinear Realizations of Internal and Space-Time Symmetries. In Proceedings of the X-th Winter School of Theoretical Physics, Karpacz, Wroclow, 19 February–4 March 1973; Volume 1, p. 117. [Google Scholar]
  14. Ferrari, R.; Picasso, L.E. Spontaneous breakdown in quantum electrodynamics. Nucl. Phys. B 1971, 31, 316–330. [Google Scholar] [CrossRef]
  15. Niedermaier, M.; Reuter, M. The Asymptotic Safety Scenario in Quantum Gravity. Living Rev. Rel. 2006, 9, 5–173. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Pottel, S.; Sibold, K. On the Perturbative Quantization of Einstein-Hilbert Gravity. arXiv 2020, arXiv:2012.11450. [Google Scholar]
  17. Alvarez-Gaume, L.; Kehagias, A.; Kounnas, C.; Lüst, D.; Riotto, A. Aspects of Quadratic Gravity. Fortsch. Phys. 2016, 64, 176–189. [Google Scholar] [CrossRef] [Green Version]
  18. Salvio, A.; Strumia, A. Agravity. J. High Energy Phys. 2014, 06, 080. [Google Scholar] [CrossRef]
  19. Salvio, A. Quadratic Gravity. Front. Phys. 2018, 6, 77. [Google Scholar] [CrossRef] [Green Version]
  20. Stelle, K.S. Renormalization of Higher Derivative Quantum Gravity. Phys. Rev. D 1977, 16, 953–969. [Google Scholar] [CrossRef]
  21. Kaku, M.; Townsend, P.K.; van Nieuwenhuizen, P. Gauge Theory of the Conformal and Superconformal Group. Phys. Lett. B 1977, 69, 304–308. [Google Scholar] [CrossRef]
  22. Fradkin, E.S.; Tseytlin, A.A. Conformal Supergravity. Phys. Rep. 1985, 119, 233–362. [Google Scholar] [CrossRef]
  23. Salvio, A.; Strumia, A. Agravity up to infinite energy. Eur. Phys. J. C 2018, 78, 124. [Google Scholar] [CrossRef]
  24. Julve, J.; Tonin, M. Quantum Gravity with Higher Derivative Terms. Nuovo Cim. B 1978, 46, 137–152. [Google Scholar] [CrossRef]
  25. Fradkin, E.S.; Tseytlin, A.A. Renormalizable asymptotically free quantum theory of gravity. Nucl. Phys. B 1982, 201, 469–491. [Google Scholar] [CrossRef]
  26. Holdom, B.; Ren, J. QCD analogy for quantum gravity. Phys. Rev. D 2016, 93, 124030. [Google Scholar] [CrossRef] [Green Version]
  27. Kugo, T.; Terao, H.; Uehara, S. Dynamical Gauge Bosons and Hidden Local Symmetries. Prog. Theor. Phys. Suppl. 1985, 85, 122–135. [Google Scholar] [CrossRef]
  28. Hamazaki, T.; Kugo, T. Defining the Nambu-Jona-Lasinio model by higher derivative kinetic term. Prog. Theor. Phys. 1994, 92, 645–668. [Google Scholar] [CrossRef] [Green Version]
  29. Butter, D. N=1 Conformal Superspace in Four Dimensions. Annals Phys. 2010, 325, 1026–1080. [Google Scholar] [CrossRef] [Green Version]
  30. Ostrogradsky, M. Memoires sur les equations dierentielle relatives au probleme des isoperimetres. Mem. Ac. St. Petersbourg VI 1850, 4, 385. [Google Scholar]
  31. Samokhvalov, S. The problem of motion in gauge theories of gravity. Gen. Rel. Grav. 2021, 53, 44. [Google Scholar] [CrossRef]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kugo, T. Noether Currents and Maxwell-Type Equations of Motion in Higher Derivative Gravity Theories. Symmetry 2021, 13, 1408. https://doi.org/10.3390/sym13081408

AMA Style

Kugo T. Noether Currents and Maxwell-Type Equations of Motion in Higher Derivative Gravity Theories. Symmetry. 2021; 13(8):1408. https://doi.org/10.3390/sym13081408

Chicago/Turabian Style

Kugo, Taichiro. 2021. "Noether Currents and Maxwell-Type Equations of Motion in Higher Derivative Gravity Theories" Symmetry 13, no. 8: 1408. https://doi.org/10.3390/sym13081408

APA Style

Kugo, T. (2021). Noether Currents and Maxwell-Type Equations of Motion in Higher Derivative Gravity Theories. Symmetry, 13(8), 1408. https://doi.org/10.3390/sym13081408

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop