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

Divergence of curvature scalars * metric

  1. Apr 17, 2010 #1
    How can one work out what terms like:
    are in terms of the divergence of the Ricci curvature or Ricci scalar?

    One student noted that since:
    [tex]G^{ab} = R^{ab} - \frac12 g^{ab}R[/tex]
    [tex]{G^{ab}}_{;b} = 0[/tex]
    that we could maybe use the fact that
    [tex]G^{ab}G_{ab} = R^{ab}R_{ab} - \frac12 R^{ab}g_{ab}R - \frac12 R R_{ab}g^{ab} + \frac14 RRg^{ab}g_{ab} = R^{ab}R_{ab} [/tex]
    to help? We weren't sure where to go next.

    Can someone better with tensor manipulations show how we could work this out?
    Last edited: Apr 17, 2010
  2. jcsd
  3. Apr 17, 2010 #2


    User Avatar
    Science Advisor
    Homework Helper

    You have divergence (Ricci tensor) = 1/2 divergence (Ricci scalar). For your original expression just apply the product rule using that rule. Remember that the metric is divergence free.
  4. Apr 17, 2010 #3
    Basically, you are saying I should be able to use
    [tex] {R^{ab}}_{;b} = \frac12 g^{ab}R_{,b}[/tex]
    to help here.

    I feel like an idiot, but I really don't see how to use that.

    If I expand out, I get:
    [tex](g^{cd}R^{ab}R_{ab})_{;d} = {g^{cd}}_{;d}R^{ab}R_{ab} + g^{cd}{R^{ab}}_{;d}R_{ab} + g^{cd}R^{ab}R_{ab;d}[/tex]
    As you mention, the first of those is just zero, since the metric covariant derivatives are zero.

    But how do I relate terms like
    to terms like
  5. Apr 17, 2010 #4


    User Avatar
    Gold Member

    Applying Liebnitz


    might help. I have to dash ...

    [I just noticed this is what nicksauce said ]
  6. Apr 17, 2010 #5
    Yes, that is what I wrote above. I just don't know how to relate terms containing
    to the divergence
    which after that, I can of course relate to the divergence of the Ricci scalar.

    What I'm actually trying to get at, is a friend asked what is the next highest order term which we can add to the EFE such that:
    1) still shows T_uv is divergenceless
    2) has the same vacuum equations (R_uv = 0)

    I thought it would be a neat thing to work out for myself, but I'd need to know the divergence of more complicated terms. So I could put in something simple and hopefully a couple counter terms to make the divergence zero. But I quickly end up with terms I can't figure out how to simplify further so I can't tell if stuff cancels or not.

    I'll admit I tried to sneak a peak at the answer, but searching papers I only found examples of people manipulating the Hilbert action. There's probably someone who already worked it out for the field equations directly, but if so I didn't find it. The reason I want to focus on the field equations is this sometimes allows simplier changes, for example imagine adding magnetic monopoles from the Maxwell E and B field equations point of view -- which is easy, vs trying to add it in an action formalism which no one has found a consistent local way to do it yet. So the simplest change in the action formalism, isn't necessarily the simplest change to the field equations, to meet the desired conditions.
  7. Apr 17, 2010 #6


    User Avatar
    Science Advisor

    Don't look if you want to have the fun :bugeye: of working it out for yourself (I believe these still give local energy conservation, but the vacuum equations are not the same):

    http://www.gravityresearchfoundation.org/pdf/awarded/2008/Padmanabhan_2008.pdf [Broken]
    Last edited by a moderator: May 4, 2017
  8. Apr 19, 2010 #7
    Thanks atyy.
    That first paper makes me a bit uncomfortable. For instance "Horizons are always observer dependent" and other claims; it all seems kind of sloppy. Especially declaring at the end that there have to be "atoms of spacetime".
    As for the others, I've seen f(R) and Lovelock gravity before while searching, but those papers are from the point of view of modifying the action. Which I commented on my thoughts about above. It does give me an idea though...

    if the action doesn't have explicit time or spatial dependence, I guess I can use it as a tool to "generate" more terms I know (from physical principles) have zero divergence. Right?

    So going along that idea (and using some hints from wikipedia), does this look right?
    (interspersed with blatant copying from wikipedia, to speed up the math)

    Consider the action
    [tex] S = \int (\frac{1}{2\kappa}R_{ab}R^{ab} + \mathcal{L}_m) \sqrt{-g} d^4x[/tex]

    The variation of the determinant is:
    [tex]\delta \sqrt{-g}= -\frac{1}{2} \sqrt{-g} g^{ab} \delta g_{ab}[/tex]

    The variation of the curvature scalar is:
    \delta (R_{ab}R^{ab}) &= R^{ab} \delta R_{ab} + R_{ab} \delta R^{ab} \\
    &= 2 R^{ab} \delta R_{ab} \\
    &= 2 R^{ab}(\nabla_c \delta \Gamma^c_{ab} - \nabla_b \delta \Gamma^c_{c a})

    Since [tex]\delta \Gamma^\lambda_{\mu\nu}\,[/tex] is actually the difference of two connections, it should transform as a tensor. Therefore, it can be written as
    [tex]\delta \Gamma^\lambda_{\mu\nu}=\frac{1}{2}g^{\lambda d}\left(\nabla_\mu\delta g_{d\nu}+\nabla_\nu\delta g_{d\mu}-\nabla_d\delta g_{\mu\nu} \right)[/tex]

    and substituting in the equation above one finds:
    \delta R_{ab} &= [\nabla_c \delta \Gamma^c_{ab} - \nabla_b \delta \Gamma^c_{c a}] \\
    &= [\nabla_c \frac{1}{2}g^{c d}\left(\nabla_a\delta g_{db}+\nabla_b\delta g_{da}-\nabla_d\delta g_{ab} \right)
    -\nabla_b \frac{1}{2}g^{c d}\left(\nabla_c\delta g_{da}+\nabla_a\delta g_{dc}-\nabla_d\delta g_{ca} \right)] \\
    &= \frac{1}{2}g^{c d}[\left(\nabla_c\nabla_a\delta g_{db}+\nabla_c\nabla_b\delta g_{da}-\nabla_c\nabla_d\delta g_{ab} \right)
    -\left(\nabla_b\nabla_c\delta g_{da}+\nabla_b\nabla_a\delta g_{dc}-\nabla_b\nabla_d\delta g_{ca} \right)] \\
    &= \frac{1}{2}g^{c d}[\nabla_c\nabla_a\delta g_{db}-\nabla_c\nabla_d\delta g_{ab} - \nabla_b\nabla_a\delta g_{dc} + \nabla_b\nabla_d\delta g_{ca}]
    (hmm... somehow the overall sign is off compared to the f(R) gravity wiki article. But magnitudes match.)

    In the action, this term in integrated, and via integration by parts, we can move the derivatives from on the metric piece to on the R_ab piece.

    After the integration by parts, we have:

    [tex] \begin{align*}
    \delta [\sqrt{-g} R_{ab} R^{ab}] &= \sqrt{-g}[ - \delta g_{ab}\frac12 g^{ab} R^{cd}R_{cd} +
    [ \delta g_{db}\nabla_c\nabla_a g^{c d} R^{ab}-\delta g_{ab}\nabla_c\nabla_d g^{c d} R^{ab} - \delta g_{dc}\nabla_b\nabla_a g^{c d} R^{ab} + \delta g_{ca}\nabla_b\nabla_d g^{c d} R^{ab}
    ]] \\
    &= \sqrt{-g}[ - \delta g_{ab}\frac12 g^{ab} R^{cd}R_{cd} +
    [ \delta g_{ab}\nabla_c\nabla_d g^{c a} R^{db}-\delta g_{ab}\nabla_c\nabla_d g^{c d} R^{ab} - \delta g_{ab}\nabla_c\nabla_d g^{b a} R^{dc} + \delta g_{ba}\nabla_c\nabla_d g^{b d} R^{ac}
    ]] \\
    &= \sqrt{-g}\delta g_{ab}[ - \frac12 g^{ab} R^{cd}R_{cd} +
    \nabla_c\nabla_d g^{c a} R^{db} - \nabla_c\nabla_d g^{c d} R^{ab} - \nabla_c\nabla_d g^{b a} R^{dc} + \nabla_c\nabla_d g^{b d} R^{ac}
    ] \\
    &= \sqrt{-g}\delta g_{ab}[ - \frac12 g^{ab} R^{cd}R_{cd} +
    \nabla^a\nabla_d R^{db} - \nabla^d\nabla_d R^{ab} - g^{b a} \nabla_c\nabla_d R^{dc} + \nabla_c\nabla^b R^{ac}
    ] \\
    &= \sqrt{-g}\delta g_{ab}[ - \frac12 g^{ab} R^{cd}R_{cd} +
    \nabla^a\frac12 g^{db}R_{;d} - \nabla^d\nabla_d R^{ab} - g^{b a} \nabla_c\frac12 g^{dc}R_{;d} + \nabla^b \frac12 g^{ac} R_{;c}
    ] \\
    &= \sqrt{-g}\delta g_{ab}[ - \frac12 g^{ab} R^{cd}R_{cd} +
    \frac12\nabla^a\nabla^b R - \nabla^d\nabla_d R^{ab} - g^{b a} \frac12 \nabla^d\nabla_d R + \frac12 \nabla^a\nabla^b R
    ] \\
    &= \sqrt{-g}\delta g_{ab}[ - \frac12 g^{ab} R^{cd}R_{cd} +
    \nabla^a\nabla^b R - \nabla^d\nabla_d ( R^{ab} + \frac12 g^{ab} R)

    So, the field equations for this "theory" are:
    [tex]- \frac12 g^{ab} R^{cd}R_{cd} +
    \nabla^a\nabla^b R - \nabla^d\nabla_d ( R^{ab} + \frac12 g^{ab} R) = \kappa T^{ab}[/tex]

    The whole point of this was to find the divergence of metric * scalar invariant. So this method seems to have generated the result that:

    [tex]- (\frac12 g^{ab} R^{cd}R_{cd})_{;b} +
    \nabla^a\nabla^b R_{;b} - \nabla^d\nabla_d ( {R^{ab}}_{;b} + \frac12 g^{ab} R_{;b}) = 0[/tex]
    [tex]- (\frac12 g^{ab} R^{cd}R_{cd})_{;b} +
    \nabla^a\nabla^b R_{;b} - \nabla^d\nabla_d ( \frac12 g^{ab} R_{;b} + \frac12 g^{ab} R_{;b}) = 0[/tex]
    which simplifying is
    [tex] \begin{align*} (g^{ab} R^{cd}R_{cd})_{;b}
    &= 2 (\nabla^a\nabla^b R_{;b} - \nabla^d\nabla_d g^{ab}R_{;b}) \\
    &= 2 (\nabla^a\nabla^b\nabla_b R - \nabla^d\nabla_d \nabla^a R) \\
    &= 0

    Is this really correct?
    That's a pretty neat way of generating divergence equations. But it seems pretty involved.

    As an aside, does:
    [tex]- \frac12 g^{ab} R^{cd}R_{cd} + \nabla^a\nabla^b R - \nabla^d\nabla_d ( R^{ab} + \frac12 g^{ab} R) = \kappa T^{ab}[/tex]
    have compatible vacuum solutions with:
    [tex]R^{ab} - \frac12 g^{ab} R = \kappa T^{ab}[/tex]
    If so, I guess I answered the other goal of playing with these as well. But I'm not sure how to say difinitively either way. I mean
    [tex]R^{ab} = 0[/tex]
    clearly solves
    [tex]- \frac12 g^{ab} R^{cd}R_{cd} + \nabla^a\nabla^b R - \nabla^d\nabla_d ( R^{ab} + \frac12 g^{ab} R) = 0[/tex]
    [tex]- 2 R^{cd}R_{cd} + \nabla^a\nabla_a R - \nabla^d\nabla_d ( R + 2 R) = 0[/tex]
    [tex]R^{cd}R_{cd} = -\nabla^a\nabla_a R[/tex]
    but I don't know if boundary conditions, or something else I'm ignoring, would screw things up.

    EDIT: Needed to fix quite a few mistakes and typos.
    Last edited: Apr 20, 2010
  9. Apr 20, 2010 #8


    User Avatar
    Science Advisor

    Are you looking for Gauss-Bonnet terms in the action, perhaps?
  10. Apr 20, 2010 #9
    I was using that action merely as a way to calculate the divergence of some terms. Just an idea that popped in my head after thinking about the f(R) theories some more after atyy mentioned them. I'm not even sure if that is a valid method to calculate divergences ... but if so, it is an interesting method for generating divergence relations.

    Can someone skim through my reasoning in the last post and tell me if that is even a valid method?

    And even better, while it may be tedious, would someone mind checking my math there? It seems really strange to me that:
    [tex](g^{cd}R^{ab}R_{ab})_{;d} = 0[/tex]
    and therefore [tex]\nabla^c(R^{ab}R_{ab}) = 0[/tex]
    The curvature scalar is covariantly constant!? How can that be true? And if it is, what does that mean? Should this have been obvious from a different line of attack (maybe something along the lines of the openning post, which I couldn't figure out what to do with)?
  11. Apr 20, 2010 #10
    If we calculate the variation of the scalar curvature [tex]R=R_{\mu\nu}g^{\mu\nu}[/tex], one gets

    \delta (R_{\mu\nu}g^{\mu\nu}) &= R_{\mu\nu} \delta g^{\mu\nu}+g^{\mu\nu}\delta R_{\mu\nu} \\
    &= R_{\mu\nu} \delta g^{\mu\nu}+\nabla_{\kappa}(g^{\mu\nu}\delta\Gamma_{\mu\nu}^{\kappa}-g^{\mu\kappa}\delta\Gamma_{\mu\rho}^{\rho}). \\
    So the further calculations all have to be incorrect!

  12. Apr 21, 2010 #11
    R is not the only curvature scalar. It should be very clear from what I wrote that I was talking about the curvature scalar [tex]R_{ab}R^{ab}[/tex] not R. So I think you misunderstood what I trying to do.

    Anyway, I did find some errors. There are a couple places where I was sloppy with order of operations, which is important since covariant derivatives don't commute. I also dropped a term in the variation of R_ab R^ab. I also tracked down where that sign error was coming from, which is a minor fix. So I need to rewrite those calculations.

    I realize tracking though my math is probably tedious, but the general line of argument is fairly simple:
    1) I take a Lagrangian that doesn't have explicit time or spatial dependence.
    2) I therefore expect the stress energy should be divergence-less, and hence can get a new relation on the divergence of some curvature terms.

    So even without trudging through the math, I hope some people here can comment.
    Is this method valid? If not, I don't really want to take the time to add in those extra terms to fix the calculation.

    Either way, the point of this thread was to learn how to calculate the divergence of terms like (as I said in the openning post)
    if someone knows another way to generate divergence relations of curvature terms besides starting with the second Bianchi identity, please do let me know.
  13. Apr 21, 2010 #12
    Yet again it is not true, to the best of my knowledge! The variations of [tex]R^{ab}[/tex] and its covariant form are not the same at all! You are putting both equal and this gives rise to an incorrect expression for the so-called "curvature scalar" you've defined here!

  14. Apr 22, 2010 #13
    I thank you for responding to my threads, but I wish you would take a bit more time reading what I wrote. As I said in the last post: "I also dropped a term in the variation of R_ab R^ab." Also, I am not sure why you feel the need to refer to it as a 'so-called "curvature scalar"' ... it IS a curvature scalar. There are many curvature scalars

    Regarding my comment about dropping a term, if we look at:
    [tex] \delta (R^{ab}R_{ab}) = \delta (R_{cd}g^{ca}g^{db}R_{ab}) = 2 R_{cd}g^{ca}g^{db} \delta R_{ab} + 2 R_{cd}g^{ca} R_{ab} \delta g^{db} = 2 R^{ab}\delta R_{ab} + 2 g^{ca}R_{cd}R_{ab} \delta g^{db}[/tex]
    you can see that I dropped the second term, which being already proportional to the metric variation adds a term directly to the field equations.

    Anyway, I found an article which states the result of considering the curvature scalars:
    [tex]R, R^{ab}R_{ab}, R^{abcd}R_{abcd}[/tex]

    So let's skip the math and return to my question:
    The general line of argument is fairly simple,
    1) I take a Lagrangian that doesn't have explicit time or spatial dependence.
    2) I therefore expect the stress energy should be divergence-less, and hence can get a new relation on the divergence of some curvature terms.

    So even without trudging through the math, I hope some people here can comment.
    Is this method valid?
    Last edited: Apr 22, 2010
  15. Apr 22, 2010 #14
    Yes, but look at the the notation! Later on in your calculations, you make use of [tex]R[/tex] as the curvature scalar within the field equations and I can't understand how that is going to be possible to introduce the fundamental laws of GR with this [tex]R.[/tex] Either you've written the post in a bad way or I'm getting it all wrong, this confusion in the text must be somehow cleared out so I can tell you if the method works or not!

    Your problem is that you are just on a roll when doing all calculations and seems like you don't even trouble yourself to check what you've written clearly right after everything is done! I don't see in your post #7 anything impying the fact that you dropped the second term in the equations!

    Anyways, if the math is seamless, why not? We can admot a new field equation with a new curvature invariant as once Einstein published his own and we read it and started to believe in it! But there are so many other problems with the math: the following is completely erroneous and misleading:

    [tex]- (\frac12 g^{ab} R^{cd}R_{cd})_{;b} + \nabla^a\nabla^b R_{;b} - \nabla^d\nabla_d ( {R^{ab}}_{;b} + \frac12 g^{ab} R_{;b}) = 0.[/tex]

    This is the first time I see something like this.

  16. Apr 23, 2010 #15
    Ugh. Altabeh, again, I appreciate that you take the time to write a response... but it is incredibly frustrating when you reply without apparently taking the time to read what I have even written or am asking. It is very frustrating. Every post of yours in this thread, and a couple in other threads have been like that.

    Let me give examples in this thread:
    From the opening post, and from the very start of my calculations in post #7, it should have been clear that I was trying to work out a divergence relation involving the divergence of a term metric * the scalar R_ab R^ab.

    You write a post saying everything is wrong because you feel the scalar I want is R. And therefore dismiss everything. Clearly showing you didn't even read enough to get the basics of what I am even asking.

    I then respond explaining again that the curvature scalar I'm focussing on is R_ab R^ab. I also point out that I've redone some of the calculations and found multiple errors. I list some of those errors. I then point out that one can ignore the math for now, as the more helpful information is whether this method for obtaining divergence relations is valid at all.

    You again seemingly ignore most of that and focus on the old math instead ... and again ignoring what I wrote, you even complain about an error I already pointed out. And yet again complain about the scalar I'm focussing on, even casting doubt on it by referring to it as a 'so-called "curvature scalar"'.

    I respond again. Stressing yet again the curvature scalar I want to focus on is R_ab R^ab. Since you complained about errors, I stated more explicitly this time one of the errors. I again stress that focussing on the math is not as helpful as telling me whether the general method is viable. I summarize this method yet again. I even post a paper I found that gives the result I should find if I go back and fix my errors, so that we can stop focussing on the details of working out the field equations.

    What do you do? You yet AGAIN complain about the scalar. You AGAIN refer to errors in my math all the way back in post 7 (which I've already pointed out myself I found errors in), and again ignoring the important parts of the discussion.

    Do you see the pattern here?
    You seem very knowledgeable. And again, I thank you for taking the time to respond. But it is not useful, and in fact incredibly frustrating, when you seemingly ignore what I write. So please, read a bit more carefully otherwise we'll just be talking past each other every time.

    I felt the need to respond to this though. Yes, I've made mistakes. But please don't take this as "you don't even trouble yourself to check what you've written". I put a lot of time into writing up those calculations. And I have checked them and I've caught many mistakes myself, and I can see how to fix them. Since I found a paper that gave the final answer (albeit without showing the work), and I found where my mistakes were, I didn't feel the need to retype up all the calculations.

    So, back to the discussion:
    I've checked by hand for f(R) gravity that the field equations are divergenceless in general (and not just in solutions to the equations), as I expected. But it would be nice to know that this method works in general. To put in explicitly...

    If I do the following steps:
    1) put a function of any scalars I want in replace of f(R) in the f(R) lagrangian
    2) solve for the field equations
    3) look at the vacuum case (T_uv = 0)
    4) take the divergence of the vacuum field equation, thus obtaining a relationship between the divergence of many curvature parts

    Is the relation I obtain from this method true in general or only true for solutions of those field equations. I've already made heuristic arguments for why I expect them to be true in general, and therefore this is a valid method to generate relationships between the divergence of curvature terms ... but it would be nice to know for sure (and explicitly why) this is the case.
  17. Apr 23, 2010 #16
    You seem to have forgotten to take a quick peep at your post #7 in this thread: All you said there was that you introduced an scalar but called it inattentively "scalar curvature" and assigned to it the same symbol we use for the "scalar curvature" [tex]g^{\mu\nu}R_{\mu\nu}[/tex] in GR so when I asked for the reason behind defining R in this way, you suggested me to take a look at this http://en.wikipedia.org/wiki/Curvature_invariant_(general_relativity)" [Broken] where the curvature invariants are discussed briefly without any reference to the fact that all of them use distinctive symbols other than [tex]R.[/tex] So am I given this right to get confused when seeing your posts and misleading notations!?

    The reason is that you looked like you got us led to the EFE starting from this [tex]R[/tex] ending up with the scalar curvature included in the equations. I don't know what makes it all silly to you to have questions like this formed in mind after reading your convoluted post!

    Right after you were born, what was your language to be understandable to your parents!? Nothing and it is the reason the baby language is just an inspiration of the inward feelings under the circumstances the baby sees at that moment! Speaking of which, for physics to be understood well, at this level of prestige, the standard language is mathematics but your calculations are all leaky and all we can do about the procedure is to lean on our inward feelings about it unless you provide us with a great indefectible mathematical framework because you're not talking about something really elementary but fully scientific and professional in the real sense of the words! I think the answer I gave in an early post gets everything straight:

    Anyways, if the math is seamless, why not? We can admit a new field equation with a new curvature invariant as once Einstein published his own and we read it and started to believe in it! But there are so many other problems with the math:

    As the title of this thread says all we are going to talk about is mathematical so if I got the time to spend on checking the math, I'd be glad to announce the results soon! Remember that all you're doing here is strongly backed up by mathematics and if this part isn't seamless, I can't even predict what'll happen to the theory in the end!

    Your theory has to be divergence-less and this can be obtained by taking the covariant derivative of the left-hand side of the field equations without even looking at the other side. This step compels the theory 1) go astray from the main road, i.e. the generality of equations for any material distribution given, and 2) nullifies the effect of "divergence" because it's already zero! It is mandatory for this step to be modified to give the correct equations! Also

    will only be valid then in a vacuum spacetime.

    If your theory makes the following general demands met, then it works fine:

    1- The modification I'm taking about is a necessity for your theory to be "general";
    2- The relation you obtained, i.e.

    [tex] \begin{align*} (g^{ab} R^{cd}R_{cd})_{;b} &= 2 (\nabla^a\nabla^b R_{;b} - \nabla^d\nabla_d g^{ab}R_{;b}) \\&= 2 (\nabla^a\nabla^b\nabla_b R - \nabla^d\nabla_d \nabla^a R) \\&= 0\end{align*}[/tex]

    has many wrong things with it! If this was flawless mathematically, your theory at least in a limited sense would rock!

    Are these helpful to you or yet you think I'm stuck on nagging at your [tex]R[/tex]?

    Last edited by a moderator: May 4, 2017
  18. Apr 23, 2010 #17
    I am using standard notation.
    [tex]R_{ab}[/tex] is used to refer to the Ricci curvature. The curvature scalar I am talking about is, once again,

    which is of course not
    the Ricci scalar.

    It is clear from the very first line in this thread what curvature scalar I was referring to. It baffles me that you are complaining that I said in post #7 "variation of the curvature scalar" and then I wrote the curvature scalar I was referring to once again ... how can you complain that is unclear even now? Fine, you misunderstood. But I've tried to correct your misunderstanding several times now.

    The fact that you are now starting to use "R" to refer to the curvature scalar I was referring to, is only making me worry you still don't understand. Especially since I am using standard tensor notation here. There should be no ambiguity.

    There are many curvature scalars. If you need to give this one a name, wikipedia refers to it as "The principal quadratic invariant of the Ricci tensor" in that link I gave you. R^{ab}R_{ab} is unambiguous and shorter.

    I hope we have finally settled any notation issues.

    That is completely missing the point of what I'm trying to get here. I'm generating field equations from a Lagrangian, and using the method I've described several times now, in hopes to generate general relations between curvature terms. I'm not using the Lagrangian to suggest new physics.

    And in that post, as in even your latest post, you again continue to complain about the specific math in post 7 which is frustrating, considering I already pointed out myself that there were errors and therefore we should focus on whether the method is viable. Heck I even posted a paper containing the correct field equations that I would obtain if I went back and fixed the errors.

    If you insist on seeing the correct field equations, let me link the paper again:
    equation 15 are the field equations

    Going through all the math to get the field equations is unimportant if the method I'm proposing to use these to get general divergence relations is wrong. I have done no math to prove the method I suggest is valid, only given heuristic arguments. I don't know how to prove one way or the other (besides getting the equations and checking by hand that they have the property I expected). If you know how to do either of these mathematically:
    1) How to calculate the covariant gradient of a non-tensor object (like a Christoffel symbol), in which case I could calculate any gradient of a curvature scalar built from the Reimann curvature directly.
    2) How to prove that the method I'm proposing to generate divergence relations, which hold in general (and not just in some specific theory), will work.
    Then that would be of great help, and I'd love to hear those math details.
    Complaining about mistakes in post #7, that I've already pointed out myself and furthermore linked to the correct equations later, is not helpful.

    I'm sorry. I've reread this multiple times. I cannot figure out what you are trying to say here.

    Are you saying the field equations obtained this way "have to be divergence-less" due to the mathematical way in which they were derived? If so, can you please show me how you know that? That is what I've been asking about this method.

    But then you go on to say that looking at the vacuum case "compels the theory 1) go astray from the main road". And that "It is mandatory for this step to be modified to give the correct equations!" Modified to what? Astray from what? I think you are still misunderstanding what this thread is about. I am trying to obtain relations between divergences of curvature terms. And I hit about this some-what bizarre method to generate them, and I'm pursuing it because I don't know how to find them another way (and no one has suggested an alternate way).

    The final comment you make there is that what I obtain will NOT be relations which hold in general. So you seem to strongly disagree with my method, but this seems to contradict your first statement that the equations must be divergence-less.

    Did you only mean they "have to be divergence-less" in the sense that for the field equations to make physical sense the solutions to the equations must be divergence-less?
    Last edited by a moderator: May 4, 2017
  19. Apr 23, 2010 #18


    User Avatar
    Gold Member

    I believe this is the case but I can't get corroboration. I'm pretty sure the Einstein tensor obeys [itex]{G^{\mu\nu}}_{;\nu}=0[/itex] identically. Any canonical EMT obeys this, so I think G must also.

    [Edit]Got it. It's confirmed here


    It follows from the Bianchi identity of course !
    Last edited: Apr 23, 2010
  20. Apr 23, 2010 #19
    Yeah, it really seems like it should, and every case I've been able to work out by hand has worked so far, but I don't know how to prove it.

    Yeah, that's probably the easiest case. Unfortunately that series of contractions on the Bianchi identity is the only useful one. So to get anymore equations relating the divergence of curvature terms, I need to resort to something more drastic.

    The [itex]{G^{\mu\nu}}_{;\nu}=0[/itex] is the relation we get when using the gravitational action:
    [tex]S= - {1 \over 2\kappa}\int R \sqrt{-g} d^4x[/tex]

    I'm hoping that replacing R with any function of curvature scalars will yield "field equations" that are in general divergence-free. In vacuum the solutions must be, so I'm worried that all this is showing is that solutions are divergence free ... and thus I can't use this to get a new set of divergence relations for curvature terms.

    However, since the action would have no explicit time or spatial dependence regardless of the function of curvature scalars chosen, I'd expect that the field equations from such a modified lagrangian would be divergence free. As I mentioned previously, I've check it a couple cases and it has been true so far. Maybe it is a simple step to show that this is always the case, so any guidance along that direction would be much appreciated.
  21. Apr 24, 2010 #20
    Let's face it: You as someone who does not deal with physics on a regular basis (at least I can smell it from the extraordinarily flawed mathematics of yours), just defined something in your own sense and kept talking about it in some post! Okay and there is no problem with this, but what do you expect from the reader when you bring up the Einstein's field equations (EFE) with a completely different [tex]R[/tex] in the same post just a few lines right below an streak of equations involved with your own [tex]R[/tex]!? What does it mean and what should we do there to find a comparison between them while there is such huge ambigiuty and a nerve-wracking violation of energy\momentum conservation?

    So your notation is NOT satandard and to keep notations in agreement with the standards established in that paper you just cited in your latest post, you would have used [tex]P[/tex] instead of [tex]R[/tex] which, as in the same paper, refers basically to the standard scalar curvature [tex]g^{ab}R_{ab}.[/tex]

    No quite contrary! You seem to first have gone way far from the topic, and secondly I don't assume you're trying to insult me by making a duffer out of me! People need to tease out things that sound unclear to them and if this bothers you, I'll be out of here because I don't want to be gone off on this way anymore! I'm here to help so just stay cool and keep to the subject!

    Considering all the mathematics in your post #7 has serious problems in it, the relation you're giving, i.e.

    [tex]- (\frac12 g^{ab} R^{cd}R_{cd})_{;b} + \nabla^a\nabla^b R_{;b} - \nabla^d\nabla_d ( {R^{ab}}_{;b} + \frac12 g^{ab} R_{;b}) = 0[/tex]

    has only one way of being interpreted, say, logically and physically when compared to the field equations obtained in the paper "The Cosmology of Generalized Modified Gravity Models" as follows: The mathematics, as I examined this discreetly myself, does not prove the relation right and so you have to go with some other alternative which is the same method the authors of the article above take advantage of. The method is strongly based on a new parameter [tex]\mu[/tex] appearing in a sort of famous Lagrangian used earlier by many other authors (see, for instance, the equation (1) in http://arxiv.org/PS_cache/astro-ph/pdf/0306/0306438v2.pdf" [Broken] This value since is present as a multiple in the correction terms, can be adjusted the way that

    1) the amount of non-conserved quantities must get really small (this is shown to be of order 10^{-33}eV for a slow-expanding universe so the violation is really negligible);

    2) the model keeps to be on a par with both observational data and the property of self-consistency of EFE (this means that invoking the conservation laws to hold in the new theory has nothing to do with the energy tensor being zero or not. They must be satisfied in any form).

    Besides these, the EFE must be extractable from the new Lagrangian by setting the new parameter equal to some constant, basically, zero! But in your Lagrangian none of these are obeyed. Specially that your Lagrangian does not reduce to the Einstein-Hilbert Lagrangian.Also, we don't have a degree of freedom to adjust it according to the model under study. I feel like you changed your mind about using the primitive idea of getting some free values involved in the field equations.

    For what reason do we need to calculate all of them when we know that the theory is vacuous? The basic Lagrangian must be modified in such a way that the EFE be extracted from it when all the correction terms vanish. Yours is not even right in the sign of the second term on the left hand side of the EFE.

    If you have no knowledge of math, why did you write all those void equations in post #7? Was that a practice essay or a play-around with equations? If you don't know that we spend time to read them, you better never think of writing nonsense again!

    They have to be, at the very very least, minimally divergence-less as they are in Carroll's papers. But in your Lagrangian this is not possible because first a mathematical derivation cannot be obtained and secondly there is nothing to be adjusted and the Lagrangian itself does not reduce to the fundamental Lagrangian of GR. This property is indeed essential to constructing a self-consistent field equation and if it is not cared about, you have to force the material distribution to consist of dark energy to compensate for the probable lost mass\momentum in interactions with ordinary matter. Looking for a relation between scalar curvature terms in the way you are into, is not feasible unless a necessary modification of the proposed Lagrangian is carried out!

    Going astray from keeping the theory at minimal level of the violation of the conservations laws.

    Talked about it above.

    What do you mean? Explain elaborately.

    As I said, being divergence-free must be at least MINIMALLY satisfied. If not, the theory, the relations to be obtained and anything else you want to get from this method will be incorrect. How does your theory guarantee the occurring of such thing within it!?

    Last edited by a moderator: May 4, 2017
Share this great discussion with others via Reddit, Google+, Twitter, or Facebook