Motivated by experimental probes of general relativity, we adopt methods from
perturbative (quantum) field theory to compute, up to certain integrals, the
effective lagrangian for its n-body problem. Perturbation theory is performed
about a background Minkowski spacetime to O[(v/c)^4] beyond Newtonian gravity,
where v is the typical speed of these n particles in their center of energy
frame. For the specific case of the 2 body problem, the major efforts underway
to measure gravitational waves produced by in-spiraling compact astrophysical
binaries require their gravitational interactions to be computed beyond the
currently known O[(v/c)^7]. We argue that such higher order post-Newtonian
calculations must be automated for these field theoretic methods to be applied
successfully to achieve this goal. In view of this, we outline an algorithm
that would in principle generate the relevant Feynman diagrams to an arbitrary
order in v/c and take steps to develop the necessary software. The Feynman
diagrams contributing to the n-body effective action at O[(v/c)^6] beyond
Newton are derived.Comment: 39 pages. The Mathematica code used in this paper can be found at
http://www.stargazing.net/yizen/PN.html Version 2: Slight re-wording of
section on removal of accelerations in 2 PN lagrangian; comments added in
conclusion; and typographical errors fixed. Article is similar to that
published in PR