Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

Efe Pde

  1. Jan 29, 2007 #1
    Has anyone evaluated the Einstein Field Equation purely in partial derivatives wrt x,y,z,t? What does it look like?
  2. jcsd
  3. Jan 29, 2007 #2
    It looks like this:

    [tex]R_{ab} - \frac{1}{2}g_{ab}R = 8\pi T_{ab}.[/tex]

    You're free to use Cartesian coordinates if you wish, but you may run into some difficulty eventually. :rolleyes:

    On a more serious note, plenty of work has been done to cast Einstein's equations in various PDE forms. Obviously, most of this has been centered on exploiting conformal equivalence classes of three-metrics so that one reaches something along the lines of the Lichnerowicz-York equation. More specifically, if we perform the usual 3+1 decomposition of a spacetime so that [itex]\mathcal{M}\simeq\Sigma\times I[/itex] (where [itex]\Sigma[/itex] is a space-like three-manifold and I is some interval in [itex]\mathbb{R}[/itex]) then the question of solving Einstein's equations essentially reduces to that of trying to find solutions to the following quasilinear elliptic equation for [itex]\phi[/itex]:

    [tex]8\Delta\phi - R\phi + K_{ij}^{TT}K^{ij}_{TT}\phi^{-7} - \frac{2}{3}K^2\phi^5 = 0[/tex]

    Here [itex]R[/itex] is the Ricci scalar on [itex]\Sigma[/itex] and [itex]K_{ij}^{TT}[/itex] are the transverse-traceless components of the extrinsic curvature. There are vast amounts of literature on this; search prola.aps.org for more information. There is a lot known about equations like this, but it's also a really active area of research given the current interest in generating realistic initial data for simulations of gravitational wave production in things like binary black hole systems.

    There are also some rather nice attempts at putting the system in first-order hyperbolic form, but I don't know a great deal about those.
    Last edited: Jan 29, 2007
  4. Jan 29, 2007 #3
    But how many degree and how many order is this EFE PDE?
    I know is has just 4 independent variables (x,y,z,t) and 10 (16-6 due to symmetric matrix) dependent variables right?
  5. Jan 29, 2007 #4
    Well, let's take a look at it:

    [tex]8\Delta\phi - R\phi + K_{ij}^{TT}K^{ij}_{TT}\phi^{-7} - \frac{2}{3}K^2\phi^5 = 0[/tex]

    This is a quasilinear elliptic PDE. What we want to solve it for is the conformal factor [itex]\phi[/itex]. Thus, it's a PDE which, if we're very clever, we can solve for a scalar function [itex]\phi[/itex]. So, in your language, it's a second-order PDE (it's a second-order PDE because of the Laplacian [itex]\Delta=\nabla_k\nabla^k[/itex] which is a second-order differential operator that acts on [itex]\phi[/itex]).

    It might be helpful to state explicitly that the Einstein field equations, when written in the form

    [tex]G_{ab} = R_{ab} - \frac{1}{2}g_{ab}R = 0[/tex]

    are a system of PDEs which one wishes to solve for the metric tensor [itex]g_{ab}[/itex]. The Einstein field equations involve the Ricci tensor, which is itself a contraction of the Riemann tensor (the usual convention is that [itex]R_{ab} = R^{c}_{\phantom{c}acb}[/itex]). However, the Riemann tensor involves second-derivatives of the metric, so one concludes that the EFEs are actually a second-order system of nonlinear, coupled PDEs. You can indeed, as you say, use the symmetries of the Riemann tensor to recognise that there are not in fact ten separate PDEs contained within the EFEs. But if you actually want to find solutions to the field equations, attempting to solve [itex]G_{ab}=0[/itex] is crazy.
    Last edited: Jan 29, 2007
  6. Jan 30, 2007 #5


    User Avatar
    Staff Emeritus
    Science Advisor

    I think we need to step back some, here.

    I think what the topic of disucsson is what's known in the literature as "the intial value formulation of General relativity".

    One starts by specifying a 3-d spatial metric (this is the metric induced on some particular spatial hypersurface by the full 4-d metric for the Lorentz interval) on some particular 3-d hypersurface which is a "surface of constant time". One also needs to specify something equivalent to the "time derivative" of the metric on this hypersurface. Defining this is a bit tricky, but it's discussed in the textbooks.

    So the dynamical variables in this theory are the metric coefficients - more particularly, the spatial metric coefficients.

    Given these initial conditions, one can calculate how the metric coefficients change as a function of time (which is a global coordinate).

    I do see from my textbook that one winds up with results similar to coalquay404's, but I'm not quite sure where [itex]\Phi[/itex] enters the picture or what it represents.
  7. Jan 30, 2007 #6

    Chris Hillman

    User Avatar
    Science Advisor

    EFE written out explicity in terms of partials?

    It depends upon what coordinate chart and what spacetime model you are using, and what kind of solution you are looking for (what you plug in on the right hand side for the stress-energy tensor).

    To keep things simple, I'll assume you are asking about the vacuum EFE [itex]G^{\hat{a}\hat{b}} = 0[/tex], where the hats suggest that I plan to write out the components, not in a coordinate basis, but with respect to a suitable frame field.

    Next, you can write down some general coordinate chart for an arbitrary Lorentzian spacetime, and then, writing down a suitable coframe and computing the Einstein tensor, you will obtain an explicit version of the vacuum EFE such that, in principle, the solution space consists of all local vacuum solutions in gtr. This will look fairly awful, although you can organize the terms so that it takes only two pages or so to write out in all its glorious detail (allowing for algebraic symmetries); see for example Frankel, Gravitational Curvature.

    A much more common procedure would be to make additional assumptions, such as symmetry assumptions, which make greatly reduce the complexity of the EFE. For example, if you assume that your solution is static spherically symmetric, and write down a simple chart consistent with those assumptions, typically the Schwarzschild chart, the EFE reduces to a simple system of ODEs, which is easily solved, yielding the familiar Schwarzschild vacuum.

    Here is a somewhat more typical example: as it happens, I solved the vacuum EFE a few hours ago (not for the first time!), while reading a new eprint (see http://www.arxiv.org/abs/gr-qc/0701161) in which Gowdy and Edmonds introduce a new type of Gowdy vacuum. Following Gowdy's classic papers from 1971-4, they start with the chart (actually I'll make some minor changes in notation)
    [tex]ds^2 = -2 \, e^{2 \, a(u,v)} \, du \, dv + b(u,v) \, \left( e^{2 \, h(u,v)} \, dx^2 + e^{-2 \, h(u,v)} \, dy^2 \right) [/tex]
    [tex]-1 < u,v < 1, \; -\infty < x,y < \infty [/tex]
    For those who like NP tetrads, one can take as an alternative definition of the Gowdy model the tetrad
    [tex]\vec{\ell} = \partial_u, \; \; \vec{n} = e^{-2a} \, \partial_v, \; \; \vec{m} = \frac{1}{\sqrt{2 \, b}} \, \left( e^{-h} \, \partial_x + i \, e^{h} \, \partial_y \right)[/tex]
    (plus mbar).

    Then, computing either the Einstein tensor or (in the Newman-Penrose formalism) the Ricci spinor components, the vacuum EFE reduces to
    [tex]b_{uv} = 0[/tex]
    [tex]h_{uv} = \frac{-1}{2 b} \, \left( b_u \, h_v + b_v \, h_u \right)[/tex]
    [tex]a_u = \frac{b_{uu}}{2 b_u} - \frac{b_u}{4 b} + \frac{b \, h_u^2}{b_u}[/tex]
    [tex]a_v = \frac{b_{vv}}{2 b_v} - \frac{b_v}{4 b} + \frac{b \, h_v^2}{b_v}[/tex]
    where the subscripts indicate partial derivatives, a very common notion in the literature when PDEs are mentioned. (Yes, for those studying the original eprint, there is a missing + in Equation (4) of the current version as of this morning.)

    Here, the first equation is just the usual one-dimensional wave equation in double null coordinates [itex]u,v[/itex]. Once one chooses a solution of the wave equation and plugs into the second equation, one can solve for h. Finally, a is then completely determined by quadrature.

    It turns out that the pattern of a pair of master equations (second order coupled nonlinear PDEs), plus two first order equations giving a third function in terms of the first two, is fairly typical for spacetimes with two independent Killing vectors, except that here the master equations also partially uncouple.

    (Computing the Weyl spinor components, the Gowdy vacuums turn out to be Petrov I in general.)

    In the new eprint, Gowdy and Edmonds suggest starting with the following solution of the wave equation:
    [tex]b(u,v) = u^2-v^2[/tex]
    (In Gowdy's previous papers, he started with some other simple solutions of the wave equation.) Then a "three parameter solution" (some of this alleged freedom actually involves choice of coordinate chart) of the remaining equations is given by separation of variables:
    [tex]h(u,v) = \frac{1}{\sqrt{ (A \, u^2 + B) \, (A \, v^2 + B) }} [/tex]
    [tex]a(u,v) = \log \left( \frac{\sqrt{uv}}{ (u+v)^{1/4} \, (u-v)^{1/4} } \right) + \frac{A^2 \, (u+v)^2 \, (u-v)^2}{8 \, (A \, u^2 + B)^2 \, (A \, v^2 + B)^2} + C[/tex]
    So this is a simple explicit example of a Gowdy-Edmonds vacuum.

    This only takes a few seconds (casesplit followed by pdsolve in GRTensorII), so we've come a long way since Einstein gloomily predicted in 1915 that his new gravitational field equations were so complicated that they would probably never be solved in closed form. This prognostication was almost immediately disproven by Karl Schwarzschild; to be fair, Einstein no doubt meant that he doubted that the general vacuum solution of his field equation could be written down in closed form. In recent decades, however, Penrose and Connes have invested considerable effort in attempts to write down the general vacuum solution, in some sense, in closed form, so far with limited success.

    Gowdy and Edmonds point out that their new family of solutions model a linear polarized gravitational wave propagating on a Kasner vacuum background. Some Gowdy vacuums also turn out to admit interpretations as the interaction zone of a colliding plane wave (CPW) model, although I haven't yet had a chance to check whether the neccessary condition for this to happen is satisfied here.

    One last point: the approach I have taken here illustrates the most elementary thing one might mean by "EFE written out in terms of partial derivatives". The initial value formulation provides a very important and somewhat more sophisticated approach. One reason why casting the EFE into "hyperbolic" form is so important is that this allows us to quote some general results in order to conclude that the initial value formulation of general relativity is (in a technical sense) well-posed.
    Last edited: Jan 30, 2007
  8. Jan 30, 2007 #7
    This is reasonably correct apart from the claim that "one starts by specifying a 3-d spatial metric." The entire point of the initial value formulation of general relativity is that GR can be completely determined from a singular Lagrangian. Since the Lagrangian is singular, this implies that general relativity is, in some not very clear sense, a gauge theory. The gauge freedom in general relativity arises from one's ability to employ four-dimensional diffeomorphisms to map solutions of the EFEs into equivalent solutions. Since the theory is based on a singular Lagrangian, it is also constrained. This is the crucial point: since there are constraints in the theory, one cannot just specify a 3-d metric and hope that everything works out. The three-metric needs to satisfy the constraint equations before it can be deemed to be satisfactory initial data.

    Look, in the 3+1 approach to GR, one has four constraint equations. Three of these are encapsulated in what is called the momentum constriaint:

    [tex] \mathcal{H}^i = -2\nabla_j\pi^{ij} = 0[/tex]

    where [itex]\pi^{ij}[/itex] is the momentum conjugate to the three-metric, and [itex]\nabla[/itex] is a [itex]g_{ij}[/itex]-compatible covariant derivative. The momentum constraint is well understood and really presents no difficulty when trying to, for example, study numerical evolutions of the field equations. The remaining constraint, known as the Hamiltonian constraint, is given by

    [tex]\mathcal{H} = \frac{1}{\sqrt{g}}(\pi_{ij}\pi^{ij} - \frac{1}{2}\pi^2) - gR = 0[/tex].

    It is the Hamiltonian constraint which is the very essence of general relativity: it encapsulates the statement that one is free to slice a four-dimensional solution to Einstein's equations in (almost) whatever way one likes -- this is general covariance by another name. (The Hamiltonian constraint is also responsible for almost every single difficulty involved in attempts to quantize general relativity, but that's a different story. :smile: )

    The point is that, in practice, you do not begin by specifying a three-metric as initial data. You can do that if you wish, but you need to make sure that the initial data satisfies the constraint equations. If it does, you can then plug it into the evolution equations and be guaranteed that what comes out is equivalent to a well-posed IV problem.

    There are various ways of trying to attack the IV problem. You can approach it from the (thick) sandwich point of view where you specify two three-metrics, so that they bound a four-dimensional region of spacetime. When Wheeler proposed this idea he thought it would be only a matter of time before somebody figured out how to make it work. Unfortunately, forty years later, we still know nothing about the (thick) sandwich approach. We do know that it probably doesn't work for general relativity since a similar approach to solving the IV problem for electromagnetism gives rise to divergences which can't be gotten rid of.

    An alternative approach is called the thin sandwich formulation. Here one specifies a three-metric [itex]g_{ij}[/itex] and another tensor [itex]\dot{g}_{ij}[/itex] on a spatial surface and evolves them forward. (It's crucial to note that despite our suggestive notation, [itex]\dot{g}_{ij}[/itex] is not the "velocity" of the three-metric.) The thin sandwich approach is very useful since we know plenty of nice existence and uniqueness proofs for solutions. But there are two problems with it: (1) it's not at all clear what sort of solution you'll get for given initial data; this makes it difficult to identify suitable initial data which will give you, for example, a binary black hole solution; (2) You're still dealing with tensorial quantities - this is, in general, a difficult way of going about things.

    The solution to these problems was recognised years ago by Jim York (I think he may have been at UNC at the time, with later refinement coming after he moved to Cornell), and is called the conformal thin sandwich formalism. Instead of working directly with a three-metric as initial data, you conformally transform the three-metric and work instead with conformal equivalence classes of three-metrics. Ploughing through some quite impressive algebra leads one to the conclusion that the entire problem can be reduced to that of solving the so-called Lichnerowicz-York equation, the quasilinear elliptic PDE that I mentioned earlier. This is a vast simplification of the problem and, indeed, is the only known robust method of generating solutions to the EFEs.
  9. Jan 30, 2007 #8


    User Avatar
    Staff Emeritus
    Science Advisor

    OK, I see that by totally neglecting the constraint equations I glossed over some points which turn out to be crucial.

    On some further reading, it looks like [tex]\phi^2 = \Omega[/tex], where the conformal transformation is [tex]g'_{ab} = \Omega^2 g_{ab}[/tex], so your equation for [itex]\phi[/itex] (which also appears in Wald) is the equation for the square root of the conformal factor.

    Given that the the conformal factor satisfies the above diffeq, I guess it's supposed to be "easy" to solve for the time evolution of [itex]\dot{g}_{ij}[/itex], though I'm not quite sure I understand all the details. If I'm interpreting the equations correctly, the requirement is that [itex]K_{ab} = \dot{g}_{ab}[/itex] be traceles and also satisfy [itex]\nabla^a K_{ab} = 0[/itex]. Wald writes a D rather than an [itex]\nabla[/itex], so I might be misinterpreting something.
  10. Jan 30, 2007 #9
    That's essentially correct. Unfortunately, Wald defines a conformally transformed metric as

    [tex]\overline{g}_{ij} = \Omega^2g_{ij}.[/tex]

    However, when you plough through the calculations (and they're very long!) you find that using the fourth power on the conformal factor is far more convenient. Ultimately, of course, this is just a clever way of minimizing pain when having to deal with calculations involving a metric on a three-dimensional manifold.

    Essentially this is correct. Note, however, that the extrinsic curvature is given more correctly in terms of [itex]\dot{g}_{ij}[/itex] as

    [tex]K_{ij} = -\frac{1}{2N}(\dot{g}_{ij} - 2\nabla_{(i}N_{j)})[/tex]

    where [itex]N[/itex] is the lapse and [itex]N^i[/itex] is the shift. (I think Wald actually uses a different sign to this but the expression here is the one which most frequently occurs in the literature.) The expression you gave is equivalent to imposing the gauge choice [itex](N,N^i)=(1,0)[/itex], which some people call a stable "geodesic slicing".

    The other point to note is that since we're dealing with elliptic PDEs, none of this is "easy". :smile:

    BTW, a nice, easy introduction to this sort of thing is http://relativity.livingreviews.org/Articles/lrr-2000-5/index.html [Broken].
    Last edited by a moderator: May 2, 2017
  11. Jan 31, 2007 #10
    D is the covariant derivative with respect to the spatial metric [itex]h_{ij}[/itex], whereas [itex]\nabla[/itex] is the covariant derivative with respect to the full spacetime metric. Up to a choice of sign convention, [itex]K_{ij} = (1/2N)(\dot{h}_{ij}-D_i N_j - D_j N_i)[/itex]
    Last edited: Jan 31, 2007
  12. Jan 31, 2007 #11


    User Avatar
    Staff Emeritus
    Science Advisor

    Aha! - that clears up some of the lose ends that were puzzling me, Wald mentioned that the shift & lapse functions weren't considered to be part of the dynamics, now I understand better how they are eliminated by the right gauge choice.

    And that clears up the D vs [itex]\nabla[/itex] issue as well.

    Thanks, Wandering & coalquay.
    Last edited: Jan 31, 2007
  13. Jan 31, 2007 #12
    Actually, I have a basic question of my own. Within this ADM formalism, what sort of gauge choices are possible? With and without coupling to matter? For GR, how does one decide whether a particular choice of gauge is legitimate, ADM or not?

    Also, I've thought that the lapse function and shift vector are not part of the dynamics simply because there are no time derivatives of them in the lagrangian. I think it's not necessary to gauge them away to see this?
    Last edited: Jan 31, 2007
  14. Feb 1, 2007 #13
    There are two separate answers to the allowed sorts of gauge choice.

    First, since foliating a spacetime is a purely geometrical exercise, it's pretty easy to convince oneself that essentially any choice of gauge is allowed. In the ADM formalism the most obvious gauge freedom is contained within the pair [itex](N,N^i)[/itex]. Since the lapse is just a function and the shift is a three-vector, this means that you have (at least) four gauge degrees of freedom per space point. The only geometrical restriction on the choice of gauge is that the lapse function must always be positive; apart from this one is free to choose any value of lapse and shift one likes.

    Secondly, while you can make any choice of lapse and shift that you want, the real test of whether a gauge choice is "legitimate" is whether or not it makes the equations of motion for the gravitational field easier to understand. For example, recall that the extrinsic curvature is defined by

    [tex] K_{ij} = -\frac{1}{2N}(\dot{g}_{ij} - 2\nabla_{(i}N_{j)})[/tex]

    Rearranging this immediately gives

    [tex]\dot{g}_{ij} = -2NK_{ij} + 2\nabla_{(i}N_{j)}[/tex].

    So, on the face of it, you would think that simply choosing [itex]N=1[/itex] and [itex]N^i=(0,0,0)[/itex] would make your life a lot simpler when analysing the dynamics since then the equation of motion for the three-metric reduces to

    [tex]\dot{g}_{ij} = -2K_{ij}[/tex].

    For a very, very long time, people thought that this choice of gauge was a smart move; after all, getting rid of the lapse and shift from your equations of motion by making this choice of gauge appears to simplify things a great deal, right?

    Wrong. It turns out that for many practical applications, the gauge choice [itex](N,N^i)=(1,\vec{0})[/itex] is actually a bad idea since it makes it difficult to interpret the physical evolution. Recently, in other ADM-like formulations (like the BSSN formulation), people have realised that a smart non-zero choice of shift actually makes numerical simulations run for much longer and become much more stable. A good outline of some smart choices of gauge is given in

    "Kinematics and Dynamics of General Relativity", James W. York, Jr., in Sources of Gravitational Radiation, (Larry Smarr, ed.)

    A guy at the NumRel group at Penn State has a copy of this paper available for download http://gravity.psu.edu/~ksmith/research/research_links.php [Broken]. If you're interested in these sorts of things, it's pretty much required reading.

    EDIT: There is one important caveat here: above, I said that essentially any choice of lapse and shift are allowed. This is true if you add to it the phrase "on the initial Cauchy surface". A common choice for the initial hypersurface is that it be constant mean curvature (CMC), which means that the trace of the extrinsic curvature is constant everywhere on it. If the initial data satisfies the constraints you can then evolve the data in any way you want by choosing any values of the lapse and shift at later time (this is in fact the statement that general relativity has four-dimensional general covariance).

    In general, however, subsequent hypersurfaces will not be CMC with such arbitrary choices of gauge. If you want these later hypersurfaces to be CMC also, then the lapse function (note: not the shift) generally has to satisfy a reasonably complicated PDE. Apart from this one point, everything I've said above holds true.

    This is correct. Apart from the purely geometrical process of foliating a spacetime (which really implies that the lapse and shift must be gauge), there is another way to see this based on the action. If you write down the 3+1 action (I'm assuming that [itex]\partial\Sigma=\emptyset[/itex] here, otherwise we'd need to include boundary terms):

    [tex]S = \int dt\int_\Sigma d^3x \sqrt{g}N(R + K_{ij}K^{ij} - K^2)[/tex]

    then no "time" derivatives of either the lapse of shift appear. This, coupled with the Euler-Lagrange equations, immediately tells you that they are gauge variables.

    However, there's a third way to see it too. The initially amazing thing about the ADM formalism was that they showed that the Hamiltonian for general relativity vanishes on-shell (by on-shell I mean on that subspace of the phase space for which the constraint equations are satisfied). Suppose that you have a phase space [itex]\Gamma[/itex] for a field theory with phase space coordinates [itex](\phi,\pi)[/itex]. Suppose further that you have some constraints in this theory (for argument's sake, let there be [itex]m[/itex] constraints). The constraints can generally be written as some functions

    [tex]\Psi^A(\phi,\pi) = 0,\qquad A=1,\ldots,m.[/tex]

    A moment's thought about this will convince you that the constraint equations therefore actually define a subspace, [itex]\Gamma_1[/itex], of the phase space called the primary constraint surface. A little more thought then convinces you that any physical evolution of the system, since it has to satisfy the constraint equations, must actually trace out a path in phase space which is confined completely to [itex]\Gamma_1\subset\Gamma[/itex].

    Now, there's a very nice little theorem (see, for example, Dirac's Lectures on Quantum Mechanics or Henneaux & Teitelboim's Quantization of Gauge Systems) which says that if you have some function(al) [itex]F[\phi,\pi;x)[/itex] which vanishes on the primary constraint surface then you can in fact write it as a linear combination of the constraints together with some Lagrange multipliers:

    [tex]F[\phi,\pi;x)|_{\Gamma_1} = 0
    \Rightarrow F[\phi,\pi;x) = \int_\Sigma \mu_A\Psi^A[\phi,\pi;x)[/tex].

    Now apply this to general relativity: GR is a constrained system so any function that vanishes on-shell in GR (i.e., one which satisfies the Hamiltonian and momentum constraints) can be written as a linear combination of these constraints with Lagrange multipliers. However, crucially, GR is an example of something called a reparametrization-invariant theory, so its Hamiltonian vanishes. This means that its Hamiltonian density can be written

    [tex]\mathcal{H} = \mu\mathcal{H} + \mu_i\mathcal{H}^i[/tex]

    where [itex]\mathcal{H}[/itex] and [itex]\mathcal{H}^i[/itex] are the Hamiltonain and momentum constraints, respectively, and the [itex]\mu[/itex]s are Lagrange multipliers. If you go through things to find the equations of motion, you eventually reach the surprising result that these Lagrange multipliers must be the lapse and shift in order for the equations of motion to agree with those obtained from the Lagrangian picture. Thus, working purely from a symplectic geometry perspective, we've shown that the lapse and shift are pure gauge variables.
    Last edited by a moderator: May 2, 2017
  15. Aug 23, 2007 #14
    I would appreciate if anyone can give me a tip of advice.I have to find the metric tensor g from the EFE.Particularly I want to implement the problem which is presented in the Eric Paquet's paper "Representation of 3D and 4D Objects Based on an Associated Curved Space and a General Coordinate Transformation Invariant Description",
    [www.hindawi.com/GetPDF.aspx?doi=10.1155/2007/42505 ] . I wanted to know if there are available solutions for that problem or how I can find them using a c++ library.I think I don't need time in this approach,just only space.Finally, this type of EFE has a specific name?Any ideas please?
Share this great discussion with others via Reddit, Google+, Twitter, or Facebook