I
DYNAMICS OF CONTINUA AND PARTICLES FROM GENERAL CQVARIANCE OF NEWTONIAN GRAVITATION THEORY
Ch. DUVAL * and H.P. KUNZLE
Centre de Physique Théorique, C.N.P.,S. Marseille
ABSTRACT : The principle of general «variance, which states that the total action functional in General Relativity is independent of coordinate transformations, is shown to be also applicable to the four-di:nensional geometric theory of Newtonian gravitation. It leads to the cor rect conservation (or balance) equations of continuum mecnanics as well as the equations of motion of test particles in c gravitational field. The degeneracy of the "metric" of Newtonian space-time forces us to intro duce a "gauge field" which fixes the connection and leads to a conserved current, the mass flow. The particle equati'oreare also derived from an invariant Hamiltonian structure on the extended Galilei group and a minimal interaction principle. One not only finds the same equations of notion but even the same gauge fields.
JULY 1976 76/P.849 x U.E.R. Scientifique de Marseille-Luminy
** Supported in part by the National Research Council of Canada
+ On leave from Department of Mathematics, • University of Alberta, EDMONTON (Canada T6G 2G1) POSTAL ADD5ESS : Centre de Physique Théorique, CNRS 31, chemin Joseph Aiguier F-13274 MARSEILLE CEDEX 2 France -2-
1. INTRODUCTION
The covan'ant space-time formulation of Newtonian gravi tation theory goes back to Frank [f| and Cartan P2j and has since been developed by various authors £3-8 J • The Newtonian gravitational field equations and the balance equations of continuum mechanics are" here formulated in terms of a connection and its curvature on a 4- dimensional space-time exactly as in General Relativity theory, except t.iat the Lorentz metric is replaced by a Galilei structure [7l, i.e. essentially a degenerate "metric" that singles out a class of local reference frames related to each other by homogeneous Galilei trans formations. Neither the equations for particles, nor for fields were obtained from a covariant variational principle, but rather they were found by analogy with the relativistic theory.
In fact, it seems that such a Lagrangian formulation can not exist, at least not in terms of the given variables. For example, on the tangent bundle of flat space-time taken as evolution space E of an isolated 1-particle system, there does not exist a Galilei invariant Lagrangian, while a Lorentz invariant one does [91. In par ticle mechanics the problem can be solved in two ways. Either one keeps the original evolution space but gives up the Lagrangian and requires only an invariant Hamiltonlan structure to exist (i.e. a presyrcplectic 2-form C on TlR. which defines- both equations of motion and Poisson brackets), or one enlarges the evolution space E to tome E such that a projection mapTTrE —•* E exists and attempts to find an invariant
• 1-form 9 such that el & B TV*V . Such a Cartan 1-form 6 is normally obtained from a Lagrangian (e.g. &= |^=^ dc\ L9>10J)>
but it may exist in cases where no Lagrangian is^defined (like on certai evolution spaces for particles with spin [_11 » 12J ) and replaces tne la ter for almost all purposes. In particular, the minimal interaction principle can be directly applied to Q (though not to G~ ), In
6/P.849 the last section we will illustrate this procedure by constructing a model of particle with spin in a Newtonian gravitational field starting from group theoretical considerations and the minimal interaction prin ciple. The same procedure starting with the Poincaré group would lead tc the corresponding relativistic model in a curved space-time.
For a classical field theory, it is not yet too well under stood what local differential geometric object corresponds to the Hamil- tonian structure in mechanics (for attempts to identify this-object see [l3| ). On the other hand it is always possible to enlarge the number of variables on which the Lagrar.gian density may depend (i.e. to intro duce some sort of "g^uge fields"). It might give new insights into the structure of field theories to be able to derive the covariant form of tr.ô Newtonian equations of gravitation this way. But this appears to be a non trivial open problem.
What will be solved in this paper is the related problem of how to derive the correct conservation (balance) equations of non relativistic continuum theory from the principle of general covariance. That means we require some action functional to be invariant under arbi trary space-time coordinate transformations as well as gauge transforma tions for the introduced gauge fields. This approach is well known in General Relativity (see, e.g. [}{]). We wish to show that the same principle applies just as well to the classical Newtonian theories.
The procedure is, however, somewhat less straightforward. It cannot be based directly on the Galilei structure and the connection as primary variables. Vie find it necessary to introduce two "potentials"
A rt and LA that together with the Galilei structure determine the connection but not in a one-to-one way. Thus there appears a "gauge group", invariance with respect to which implies that the corresponding current is conserved. In this way the concept of mass current appears in the Newtonian theory separately from the stress-energy tensor, in
76/?.349
f -<;-
contradistinction to the raiativistic case, similarly as the mass appears in a different role in Galilei and Poincaré invariant mechanics, res pectively.
The principle of general covariance has been cast into a particularly elegant and effective form by Souriau [l5j which permits it to be directly applied also in the limiting case of matter distribu tion with support on a single worldline. Then it yields the equations of motion for test particles in a given external field. He obtain this way the equations of a particle with mass, spin and quadrupole moment in a Newtonian gravitational field, and find that they agree with those deri ved by means of the minimal interaction principle.
Section 2 reviews Souriau's approach to the principle of general covariance and test particle motion in General Relativity; section 3 the basic definitions for Newtonian structure. In the next two sections we derive the balance equations for continua and the equa tions of motion for particles, respectively, in the Newtonian theory. The last section deals with the Hamiltonian model of spinning particles.
2. PRINCIPLE OF GENERAL COVARIANCE
The Palatini method of considering the stress-energy tensor I ' of rn.it.tcran d non gravitational fields as' the variational derivative of a Lagrangian density with respect to the space-time metric ^.v has been restated by Souriau |_15j in the following form : Consider the set S of all Lorentz metrics â«/i on a given 4-dimonsional manifold V and a 1-form J" en the manifold
' ' S is an open subset of the set of sections of a vector bundle and as such has the structure of an infinite dimensional manifold. Ke will r.ot make this more precise than is usually done in the context of variational calculus used in physical applications (see, e.g., ricH).
76/P.849 -5-
*/% ^ (2.1) 'V ' where the variation °*J is interpreted as a tangent vectorto 5 at g and the integral required to exist.
The principle of general cqvariance now demands that J be a 1-form not only on S but rather on the quotient "manifold"
H =$/$*& (v ) (where "*>'$0(V ) denotes the group of all C °° diffeomorphisms of V with compact support) since any two'metrics related by a diffeomorphism are physically equivalent. This means thatf J~' must vanish on all "vectors" H. tangent to the oroit of g under $;# [v) in S, in other words that
for all vector fields fc on V with compact support (where «fç denotes the Lie derivative v/i th respect to J" ). This is well known to be equivalent to the conservation equation of the stress-energy ten sor T71 ,
v^-r-y» = „ (2.2)
Here the connection has been taken to be the unique sym metric Lorentz connection. One could replace S by the set ^ftV'^/O lV*J/*r = ° } without requiring that F^ be symmetric and (2.1) by '
TcrCv^ffjrrjj^ e'C IV^ }^3 (2.3) then ?
M;i^v©T^w, V
75/P.849 -6-
(21 leads again to equation (2.2) if one defines* '
T*.f*+ w^©•CM -r? -.*{•<-© -\ &' : y <2-4> n 1 r\~A and now requires that 1 «A and Vyi are symmetric. If, however, torsion is admitted as, for example, in the Einstein-Cartan theory flT], \^„ can be decomposed into the Levi-Cività connection
ant a cont W/s î * orsion tensor. Part of ©y is again ab sorbed in T i and the rest interpreted as a spin density which will appear as a source term in (2.2). (The ©-' has been called hyper- momentum by Hehl et al. [l8]). Since in this principle only the (global) functional is considered as a physical quantity it is clear tnat in (2.3) ~T"fl and ©»' could only be individually well defined local fields if there were no constraint at all between %*e/± anc* ^«/i • ^ l~v/i is to be a symmetric lorentz connection, then either of T'y1 or G^J/1 can be dropped and the physical stress-energy tensor defined by (2-4) (with either "T"*/1 or ©*/" equal to zero). Then equa tion (2.2) will still be a consequence of the general principle- of co- variance.
The advantage of Souriau's approach is that the same prin ciple also leads very easily to the equations of motion of test particles in the gravitational field (and other external fields). The method is summarized as follows : Suppose the functional J has its support on a worldline A in V and is a distribution of, say, no higher than second order. Then it can be assumed to be of the fora
7Lhl,\i.T\.©firmer'» }* "••' A '
(2) ( ) and f "1 denote symmetrization and antisymraetri2ation respectively.
76/P.849 -7-
where t is any curve parameter of -A and \ ' s o < ^-^z - » and Q} /** has the same symmetries as the curvature of a symme tric Lorentz connection. By the principle of general covariance we have H«^V^r?-0 VW*
for all vector fields v on V with compact support.
^ In particular, for J =• /\^i, £ where A and B are t* -functions vanishing on j\ and £ any C vector field with compact support equation (2.6) leads to o = f ©7120* £* d* A £ wnich is equivalent to \ÙJ beinç, of the form for some M . -^ , where V is the tangent vector to A with respect to t. Using this and now putting Ç^r A Ç * with /{|j\ = O we have from (2.6)
A where • :- V Vj, and integration by parts has been used. Equa tion (2.7) is equivalent to
76/P.849 for some In defined along -A. . Substituting again into (2.6) gives
X Vr Q \?0 (2.9) - M + i< - < T"f Define pi* 6*: r .2 then (2.9) becomes
x î --jR% v 5 Z. * '^ Equations (2.10) and (2.11) are now fairly accepted as equations of motion of a test particle if O ' is interpreted as its spin, x as its 4-momentum and Q / as its (mass) quadrupole moment (e.g. F19» 15, 201 arid references given therein). It is well known that they do not form a complete set of equations for all dynamical quantities necessary to describe the state of motion of a test particle, unless additional conditions are imposed, the most natural being (when çy* r 0 ) [l9. 12, 20] . Note that again not all of the coefficients T^1 , £jV • Q7i'r' ln equation (2.5) can 76/P.049 -9- be interpreted physically in an unambiguous way. To obtain the equa tion of motion of a monopole particle ©_/^ and Q r can be assu- med zero. For the dipole or multipole case the term T" "/ is not necessary. uanoxwi STRUCTURES A Galilei structure on a four-dimensional C -nr-nfold V can be given by a pair of tensor fields (fA & ) "here Y^ is symmetric and % spans the kernel of "/"A at every point of V. (See L?J for details and proofs for all the material of this section.) A vector U is called unit timelike if U H\, = ' > it is spacelike if U ^ = o . The three-dimensional distribution 1 of hyperplanes Sx fom.ed by a ! spacelike tangent vectors at each K. G. V* is completely integrate, in particular, if the 1-form !/< c p^ cix.'1 is closed : dy - o • The leaves of this folia tion then carry a natural Riemannian metric induced by J7 . Linear symmetric Galilei connections (i.e. satisfying V* ^ =o and 7* ^ a 0 ) (3.1) exist if and only if Tp is closed and are determined up to an arbi trary 2-form F = '/_; r«„ lx "A eU P • In fact, the general solution of (3.3) for syxnetric p"f can be written in the form c;*<. £yr.j,^,, ,r, r *v^ (3.2) where U is an arbitrary given unit tinelike vector field and 76/?.849 -10- n ' •* I *n is the umq-.c symmetric Galilei connection such that U . is geodesic and «rotational, i.e. such that ([3, p. 197]) U V U = o and V U1 = o (3.3) (where indices are as always raised by means =f yr,. The curvature tensor of any symmetric Galilei connection has the properties We will only consider ? subclass of sytrmetric Galilei connections, called Newtonian connections which are characterized by satisfying in addition or. equivalently, by dF : 0 (3.5) -/» any Galilei connection that is obtained as a limit of symmetric Lorentz ' connections [21^. It means that we confine ourselves to conservative forces. We will call the triple f"Y /, ip , V\,„ ) on V a Newtonian structure if (3.4) holds. A unit timelike vector field U (which can be con sidered as an observer field) is also necessary to define a "generalised inverse of If 7* by means of ? y u r o «d \J^c S,f- c^uf* (3.6) "'J -J 76/P.049 -11- Any orientable Galilei manifold carries a canonical volume element vol which can be written as It has been shown in [7] that an asymptotically Euclidean Newtonian space-time satisfying the field equations is completely equivalent to the classical Newtonian gravitation theory {with a potential *u satisfying the Poisson equation ^^ ?ttnkp in Euclidean S-sp^ce). He end this section by deriving some formulae that will later oe needed. Covariant derivatives will always refer to P^» given by {3.2) except when specified otherwise. From {3.2) w. have ' f F ^"H--VVf1 P X > -i î thus, in view of (3.3), Wr V l/ -- - F* U? (3.7) VC- « ^ 1 F*/1 s 2. Using (3.6) we can derive from thi>se equation.. X F, = .A Vxr„ V, U (3-9) The general variation of "^ is obtained from (3.6) directly as / 76/P.349 -12- (3.10 M Y(< p.) A a particular case of which is v i (3.11) ^v^ V^ 4. BALA.-;C£ EQUATIONS FOR CONTINUA IN THE NEWTONIAN THEORY Taking S now to be the set of all Galilei structures (-%T tyg } on V , one might think of replacing (2.1) by But as shown in [211 the balance equations following from the general covariance of this functional lack the acceleration term of the classical equations. This is not surprising because (if ^ • PK ) do not fully determine the connection. Even if we require the connection to be sym metric, we must therefore take this functional to be f. ^u*, m-J^y^ *si* ®?K} } "* (4-i) where ^ r In the relati- vistic theory the o Pj« can be expressed in terms of o'2~s, cor a Newtonian connection p ^ compatible with the Galilei struc ture WÛ have by (3.2,5) ' v1 - rT ip V 76/P.E19 -13- where 1 ^A is characterized by (3.3) in terms ^f derivatives of [V r, 4"a~ 1 . and a unit vector field u . A ^ is an arbit. ary 1-form. The variation ^Hx/i can thus be expressed in terns of ile W could be considered fixed. However, there is no canonical choice for U. and it seems therefore best to regard the connection fixed (in a some what redundant way) by ( *$ J / if/^ ) and the pair (A^, U )• Thus (^.1) can be replaced by %{M^)y-[\-i^u'^^st* J"f/I«+ K.^J.W (4.2) We now require that be a 1-fonn on the manifold of equivalence classes of Newtonian structures. Thus J must in prrticular be left invariant by any "gauge -transformation", i.e. a change ("f^ £ A'U"J"* C 7"^£/^*, ù*J that does not change the Newtonian structure (TA^, Pj%) • To find the infinitesimal gauge transformation suppose therefore that »?/* $"<&. *• S^JL * o . Then IU Su. : o and varying equation (3.9) gives {in view of (3.10, u)) V U - 27 P AI A[* pip a H f - 2 Tf VSue *** v V ft r Pi - ur *\ (w )• But 2.1 è a r therefore we find "y F- \i - Ï* Vi In other words variât ons such that S "Y ' T o , atfj ° ° anc! 4 **.. •r-/ 76/ = .849 -14- for any W and Therefore «J must vanish for all such variations and we have from (4.2) for arbitrary spacelike gU and ^ with compact support. This implies that VJ*--. (4.1) where K. is an arbitrary scalar. Using (4.5) we now state that sj must vanish on varia tions that are infinitesimal diffeomorphisms with compact support, that is by Stokes' theorem and (4.4), where 76/P.819 -15- Since (4.6) must hold for all vector fields Ç with compact support, we have 7p^--^pVft V Jfi- (4-8! Equation (4.4) states that O is a conserved current which is most naturally interpreted as the mass flow. Assume therefore * that -X*r « V* } V* tfa = j (4.9) and confine considerations to a region where the mass density y° is positive. Now, with the help of (3.9) expression (4.3) can be written as Vp^ .. . j,£A ^ «4.10» with ' : a V V . But we wish to replace all references to * it *i the arbitrary field U by one to the physical matter flow V . It is easily seen that V " V * * ^ v/»> ~v*f- ^ where v^y/" n 5.-^V . Then fro.-n (4.10) V 76/P.849 -16- where (3.11) has been used as well as But by (4.4) and (4.9) this can be written as v-- - x V^ (412) provided one defines Equation (4.12) is the closest analogue to the relativistic conservation law (2.2) one can find in the Nev/tonian theory, since it is not possible to absorb the right hand side in the divergence. The relativistic stress energy tensor therefore decouples in the Newtonian theory in a tensor describing stresses and energy flow, and the matter flow vector J" - a V . where TVK has the property which follows at once from (4.7) and (4.13) since Xryy^-i - c> This description seems quite satisfactory since it is rather analogous to the situation in particle mechanics. Contravariant syrrr.etric tensors "T y have also been used in Newtonian theory (e.g. [8])» but they do not describe stresses and energy (flow). Just as in the relativistic theory the tensor can be decomposed in different reference frames and then interpreted physi cally. The most canonical one is the "matter frame". Therein view of (4.14) 76/P.849 / -17- where °| Yu = o and ^^3 = ° = J^ V , ana one can interpret £ as the specific (non mechanical) internal energy den se . sity, <* US the energy lor heat) flow and b- . as the stress tensor. This decomposition corresponds to the well known one in Relativity where a nutter flow vector p V is introduced ad hoc, namely to Kith CJ Vj r o , f- I V» « O It was shown in [2lT that equations (4.12) for the decora- position (4.15) are equivalent in a flat Galilei space-time to the clas sical, balance equations of continuum mechanics, namely to A where V = ( 1, V ) , /I = I,?, 3 and k -- - PDC, - 2 ?0^ V ' is a body force density expressed in terms of the only non-vanishing components of I~I