Welcome
to my page on using the techniques usually associated with perturbative
(quantum) field theory -- see here and here -- to calculate general
relativistic corrections to the Newtonian gravitational lagrangian for
point masses.
- Yi-Zen
Update Dec 2011: I have written the Mathematica package TensoriaCalc to compute geometric tensors - Riemann, Ricci, Christoffel symbols, etc. - within a Riemannian geometry framework. This is closely related to a parallel effort of mine to write a package that can manipulate abstract tensors in physicists' index notation. (I already have a beta version of TensoriaLite, which reproduces some of the tensor contraction capabilities of FeynCalc, including the computation of the N-graviton Feynman rules for the Einstein-Hilbert lagrangian.) Stay tuned! Motivation To measure gravitational waves produced by binary black holes and/or neutron stars as they spiral towards each other, gravitational wave observatories such as GEO, LIGO, TAMA, and VIRGO need the effective lagrangian of 2 point masses up to 3 PN and beyond: O[(v/c) ^{6}],
and higher, relative to Newtonian gravity, where v is the typical speed of the
binaries in its center of energy frame and c
is the speed of light. The effective lagrangian is used to construct
theoretical templates so that the raw data can be compared against them
to see if there is significant enough correlation to claim detection;
see Blanchet
for a review. For now, I am focusing on the integer PN order,
conservative portion
of the dynamics.On the other hand, experiments such as APOLLO, GTDM, LATOR, BEACON, etc. are beginning to probe the non-Euclidean nature of the solar system's geometry beyond the currently well tested O[(v/c) ^{2}],
by measuring the timing and deflection of light propagation within the
solar system to higher accuracy than before. Knowing the geometry of
the
solar system at a theoretical level amounts to knowing the effective
lagrangian for n point masses, and vice versa.In view of the gravitational wave detection efforts, because the calculation at such a high order is very arduous, it is imperative that the necessary software be developed to automate the computation as far as possible. Code for post
Newtonian diagrams, related tools
This section contains the code I used to compute the n-body effective
lagrangian up to O[(v/c)^{4}]. Please do let me know if
you find errors / bugs.Update 28 Jan 2009: I have added in the 0 thru 2 PN Mathematica notebooks some description for FeynD, the notation used to represent Feynman diagrams. - FeynCalc I assume you have Mathematica. Install the high energy physics package
- Generator of fully connected diagrams To generate the fully connected Feynman diagrams whose evaluation would yield the effective lagrangian at a given PN order, first load the following and do Kernel > Evaluation > Evaluate Initialization:
- Sample files for Newtonian thru 2 PN Next, look at the following files for Newtonian thru 2 PN to understand the format / notation of the "Feynman diagrams" generated. A portion of the diagrams can be calculated automatically.
- N-graviton Feynman rules The diagrams that need to be evaluated by a human being can be done so manually within each cell. To calculate diagrams involving graviton vertices, we need the N-graviton Feynman rules for N > 1. Here I have used the linearized de Donder gauge. Simply open the following file, do Kernel > Evaluation > Evaluate Initialization and the Feynman rule command GravitonNVert is ready for use.
- Differentiation and integration For manipulation and differentiation of Euclidean distances, vectors, velocities, accelerations, etc., as well as the computation of the Feynman integrals, it helps to have some tools handy. Again, do Kernel > Evaluation > Evaluate Initialization after opening the file.
Future developments
To push our knowledge of binary dynamics beyond the currently
known 3.5 PN, O[(v/c)^{7}], using field theory, the first steps
to take would be to improve the efficiency of the codes in IntegerPN_FullyConnectedDiagrams.nb
and EinsteinHilbertLagrangianFeynmanRules_LinearizeddeDonderGauge.nb. |