Go to previous page Go up Go to next page

9.2 Lagrangian and Hamiltonian formulations

The conservative part of the equations of motion in harmonic coordinates (131View Equation) is derivable from a generalized Lagrangian, depending not only on the positions and velocities of the bodies, but also on their accelerations: ai = dvi/dt 1 1 and ai= dvi/dt 2 2. As shown by Damour and Deruelle [55], the accelerations in the harmonic-coordinates Lagrangian occur already from the 2PN order. This fact is in accordance with a general result of Martin and Sanz [100] that N-body equations of motion cannot be derived from an ordinary Lagrangian beyond the 1PN level, provided that the gauge conditions preserve the Lorentz invariance. Note that we can always arrange for the dependence of the Lagrangian upon the accelerations to be linear, at the price of adding some so-called “multi-zero” terms to the Lagrangian, which do not modify the equations of motion (see, e.g., Ref. [63]). At the 3PN level, we find that the Lagrangian also depends on accelerations. It is notable that these accelerations are sufficient - there is no need to include derivatives of accelerations. Note also that the Lagrangian is not unique because we can always add to it a total time derivative dF/dt, where F depends on the positions and velocities, without changing the dynamics. We find [66Jump To The Next Citation Point]
harm L = ... (137)


Witness the accelerations occuring at the 2PN and 3PN orders; see also the gauge-dependent logarithms of ' r12/r1 and ' r12/r2, and the single term containing the regularization ambiguity c. We refer to [66Jump To The Next Citation Point] for the explicit expressions of the ten conserved quantities corresponding to the integrals of energy (also given in Eq. (133View Equation)), linear and angular momenta, and center-of-mass position. Notice that while it is strictly forbidden to replace the accelerations by the equations of motion in the Lagrangian, this can and should be done in the final expressions of the conserved integrals derived from that Lagrangian.

Now we want to exhibit a transformation of the particles dynamical variables - or contact transformation, as it is called in the jargon - which transforms the 3PN harmonic-coordinates Lagrangian (137View Equation) into a new Lagrangian, valid in some ADM or ADM-like coordinate system, and such that the associated Hamiltonian coincides with the 3PN Hamiltonian that has been obtained by Damour, Jaranowski and Schäfer [60Jump To The Next Citation Point]. In ADM coordinates the Lagrangian will be “ordinary”, depending only on the positions and velocities of the bodies. Let this contact transformation be Y i(t) = yi(t) + dyi(t) 1 1 1 and 1 <--> 2, where i Y1 and i y 1 denote the trajectories in ADM and harmonic coordinates, respectively. For this transformation to be able to remove all the accelerations in the initial Lagrangian Lharm up to the 3PN order, we determine [66Jump To The Next Citation Point] it to be necessarily of the form

1 [@Lharm @F 1 ] ( 1) dyi1 = --- ----i--+ --i-+ -6Xi1 + O -8 (138) m1 @a 1 @v1 c c
(and idem 1 <--> 2), where F is a freely adjustable function of the positions and velocities, made of 2PN and 3PN terms, and where Xi 1 represents a special correction term, that is purely of order 3PN. The point is that once the function F is specified there is a unique determination of the correction term i X 1 for the contact transformation to work (see Ref. [66Jump To The Next Citation Point] for the details). Thus, the freedom we have is entirely coded into the function F, and the work then consists in showing that there exists a unique choice of F for which our Lagrangian Lharm is physically equivalent, via the contact transformation (138View Equation), to the ADM Hamiltonian of Ref. [60Jump To The Next Citation Point]. An interesting point is that not only the transformation must remove all the accelerations in harm L, but it should also cancel out all the logarithms ln(r12/r'1) and ln(r12/r'2), because there are no logarithms in ADM coordinates. The result we find, which can be checked to be in full agreement with the expression of the gauge vector in Eq. (132View Equation), is that F involves the logarithmic terms
22 G3m m [ (r ) (r )] F = -----6-12--2- m21(n12v1)ln -1'2 - m22(n12v2)ln -12'- + ..., (139) 3 c r12 r1 r2
together with many other non-logarithmic terms (indicated by dots) that are entirely specified by the isometry of the harmonic and ADM descriptions of the motion. For this particular choice of F the ADM Lagrangian reads as
harm harm ( ) LADM = Lharm + dL-----dyi1 + dL-----dyi2 + dF-+ O 1- . (140) dyi1 dyi2 dt c8
Inserting into this equation all our explicit expressions we find
LADM = ... (141)


The notation is the same as in Eq. (137View Equation), except that we use upper-case letters to denote the ADM-coordinates positions and velocities; thus, for instance N12 = (Y1 - Y2)/R12 and (N V ) = N .V 12 1 12 1. The Hamiltonian is simply deduced from the latter Lagrangian by applying the usual Legendre transformation. Posing i ADM i P 1 = @L /@V 1 and 1 <--> 2, we get [87888960Jump To The Next Citation Point66]21
ADM H = ... (142)


Arguably, the results given by the ADM-Hamiltonian formalism (for the problem at hand) look simpler than their harmonic-coordinate counterparts. Indeed, the ADM Lagrangian is ordinary - no accelerations - and there are no logarithms nor associated gauge constants ' r1 and ' r2. But of course, one is free to describe the binary motion in whatever coordinates one likes, and the two formalisms, harmonic (137View Equation) and ADM (141View Equation, 142View Equation), describe rigorously the same physics. On the other hand, the higher complexity of the harmonic-coordinates Lagrangian (137View Equation) enables one to perform more tests of the computations, notably by inquiring about the future of the constants r' 1 and r' 2, that we know must disappear from physical quantities such as the center-of-mass energy and the total gravitational-wave flux.
  Go to previous page Go up Go to next page