Tuesday, July 23, 2024

Variation of an action with the speed of light limit: Einstein-Hilbert action does not allow gravitational waves?

UPDATE August 3, 2024: We cannot use an infinitesimal variation to slow down time by dt as explained below. The problem is that if we try to "glue" the time-translated part of the history to the part which is stretched by dt, the gluing point will have a discontinuous time derivative of the metric g.

That is, if one wants to displace a segment of a sine wave a little bit, then gluing it to the old wave makes the derivative discontinuous. A discontinuous derivative of a metric may mean, for example, matter of infinite density at the point of discontinuity.

We have to check what happens if we do the gluing at the birth place of a gravitational wave.

See our blog post August 4, 2024. Our claim below is correct: there is no history where a gravitational wave is generated. But the proof is a bit different.

----

Let us have a physical system M which is initially static or "stationary", but then some local event happens, and the state starts to change. We cannot allow signals that propagate faster than the speed of light. The variation of the action S of the system must respect this principle.

What implications does the speed limit have on the variation?


An example: some of the mass M suddenly disappears


This is our favorite example of fooling the Einstein field equations. The system M is a static spherically symmetric mass. Then we use a magic trick to make some of the mass disappear suddenly. Let us denote the new mass by

       M'  <  M.

We assume weak fields.

If we were allowed to alter the metric g infinitely fast, then we could adjust the metric to that Schwarzschild metric which satisfies the field equations for the reduced mass M'. But we cannot do that.















Let us try to find a stationary point of the Einstein-Hilbert action S for this example.


                 ρ
               ###      r      ●       r      ###
             shell             M'             shell
   ^ y
   |
    -----> x


Let us imagine that we right now have a metric g which matches the mass M' at the center, and g also matches an imagined spherical shell of a mass M - M' at some radius r from the center. The shell does not have any matter, because the matter disappeared through magic, but the metric g is adjusted as if there were matter there.

The shell is not thin. The Ricci tensor R for the shell is obtained through trace reversing the mass stress-energy tensor T, which only contains the mass density ρ of the shell as its 00 component:

       R  =
              ρ / 2      0           0               0
              0           ρ / 2      0               0
              0           0            ρ / 2         0
              0           0            0         ρ / 2,

and the Ricci scalar is:

       R  =  ρ.

The metric g has the diagonal almost (-1, 1, 1, 1) because the field is weak. The metric has g₀₀ slightly larger than -1, and the x coordinate metric g₁₁ slightly larger than 1.

The Wikipedia formula gives:

       δR  =  R₀₀  *  -δg₀₀,

       δR  =  R₁₁  *  -δg₁₁,

and so on. 

If we make g₀₀ larger, that is, slow down time, then δR is negative. Also,

       sqrt(-det(g))

in the action S formula becomes smaller. Thus, the metric g is not a stationary point of the action S. We can increase g₀₀, and S declines linearly with g₀₀.

Can we make R = 0 simply by slowing down time? Hardly. Various components of R depend in a complex way on g₀₀.


Replacing the shell with gravitational waves which carry the mass-energy M - M' away


If it is possible to find a metric which describes a realistic gravitational wave, and where R = 0, then we do not need a magic trick in the example above. A spherical gravitational wave front carries the mass-energy M - M' away. The Ricci tensor is everywhere zero outside M'.

However, we can ban gravitational wave solutions by demanding that the metric has to be perfectly symmetric. There are no longitudinal gravitational waves.


Ivor Robinson and Andrzej Trautman (1960) found spherical wave solutions of the Einstein field equations. Let us look at their metric:






The wave equation is:






The authors say that in the linear approximation they have waves whose wave fronts are at

       σ  =  constant.

Apparently,

       σ  =  r  -  c t,

using more familiar coordinates, r for the radius and t for time.

The authors write that if they choose K > 0, then there is a "singularity" on each wave front, which "might represent a flow of energy" to the wave, so that energy does not run out when the wave carries it away.

The description is very far from a practical gravitational wave whose source would be a binary star. The authors neither mention anything about the metric change when the central mass M loses mass-energy in gravitational waves.


A spherically symmetric metric of time plus a shell of gravitational waves makes R ≠ 0 inside the gravitational waves


Proof sketch. We assume that the metric of time, g₀₀, is spherically symmetric. We assume that the mass-energy of a gravitational wave pulls test masses just like the mass-energy of matter does.


          M - M'
            ###            ●            ###
          shell           M'           shell


The central mass M' radiated away the mass M - M' in a shell of gravitational waves.

Let us have initially static (static in the coordinates) test masses m. Outside the shell, they feel a gravitational pull which is relatively larger than inside the shell. Therefore, test masses will become denser at the shell. Thus,

       R₀₀ > 1,

inside the shell, unless the spatial metric can cancel the focusing effect of the metric of time.


                         ^
                         |
                         •
            <---- •   ×   • ---->                               ● M'
                         •    m
                         |
                         v

             ##############
                      shell
                     M - M'


Let us consider test masses m sent from inside the shell from one static point, at low velocities to every direction. The cloud of test masses will have a smaller volume than in euclidean space, if we just consider the metric of time. The mass-energy density of the gravitational wave pulls the test masses together.

Let us reset R₀₀ to zero by making the spatial metric saddle-shaped, that is, the average spatial curvature is negative. Then the volume of the cloud of test masses is surprisingly large: that can cancel the focusing effect of the metric of time.

Let us then send a tight cone of geodesics to a spatial direction where the spatial metric has a negative Ricci curvature. Since the spatial metric overall has a negative curvature, such a direction must exist. Can we use the metric of time to cancel the negative Ricci curvature in that direction?

We can probably ignore the field of M' altogether, because its field has the Ricci tensor R = 0 outside M'.



   ^  g₀₀
   |
   |                   __--------------------
   |        ___----               
-1
                        P
                      ###
                     shell
                        r <------ M'


Let us then look at the metric of time g₀₀ created by the shell. At some r inside the shell, we have a point P where

        d²g₀₀ / dr²  <  0.


spatial
direction
  ^
  |          \         /   geodesics
  |            \     /  s
  |              \ /
  |               ×  P start point
   --------------------------> t


Let us send a cone of geodesics from P to a spatial direction where the spatial Ricci curvature is negative.

As the cone of geodesics to a spatial direction opens up, its "diameter" in terms of time changes according to g₀₀. Let s be the spatial distance from the start point P. We have:

       d²g₀₀(s) / ds²  <  0.

The metric of time adds to the negative Ricci curvature into the spatial direction. We cannot use the metric of time to cancel a negative Ricci curvature to the spatial direction.

It could still be that if the metric of time is not spherically symmetric, we might be able to cancel the negative Ricci curvature to the spatial direction.

Conjecture 1. If a central mass M magically loses a mass M - M', or radiates it away as gravitational waves, then the Ricci tensor R cannot be zero outside the remaining central mass M'.


The mass-energy density of a sine gravitational wave packet in nature is its Ricci scalar -R / κ?



Chris Hirata (2019) calculated the Ricci tensor R of a sine gravitational wave

       g  =  η  +  h₊,

and noticed that the Ricci scalar R for it is zero. He derives the mass-energy density by determining what kind of a stress-energy tensor T produces R in the Einstein field equations. The component 

       R₀₀  =  κ T₀₀

has a negative mass-energy density T₀₀. Hirata assumes that -T₀₀ is the energy density carried by the gravitational wave. Here,

       κ  =  8 π G / c⁴

is the Einstein constant.

Let us have a gravitational wave packet g which inside the packet is roughly a sine wave.

Let us "correct" the wave packet metric g by superposing the temporal (time) metric perturbation of mass-energy of the density -T₀₀ moving at a speed of light along with the wave packet. In this way, we increase R₀₀ to zero for the new metric g'.

Now we have a new metric g', for which the the Ricci scalar R is

       R  =  κ T₀₀.

The energy density is

       -R / κ.

For mass-energy of matter with no pressure, we have the energy density

       R / κ.

Is this reasonable? Could this be a stationary point of the Einstein-Hilbert action S if we take into account that faster-than-light signals are not allowed?


Superposing a matter wave (electromagnetic wave) on a gravitational wave packet


Chris Hirata calculated that the R₀₀ Ricci tensor component in a sine gravitational wave is:






Let us for a while forget about the oscillation of R₀₀ between positive and negative values. We will analyze that later.

We assume weak fields. The metric is

       g  =  η  +  h  =

                -1                0                 0                0

                 0                1 + 2 h₊      0                0

                 0                0                 1 - 2 h₊      0

                 0                0                 0                0,

where |h₊| is very small and h₊ describes a sine wave.

Let us denote 

       R₀₀  =  -κ ρ.

The Ricci tensor R of a sine gravitational wave is approximately

       R  =

            -κ ρ        0            0         κ ρ

             0            0            0            0

             0            0            0            0

             κ ρ        0            0        -κ ρ.

The Einstein tensor G is the same as R because the Ricci scalar R is approximately zero.

The stress-energy tensor T which corresponds to G is approximately

       T  =

            -ρ          0           0           ρ

             0          0           0           0

             0          0           0           0
 
             ρ          0           0          -ρ.

The stress-energy tensor T which corresponds to a particle moving at the speed of light is







where







Chris Hirata in his paper has the speed of light c set to 1, as well as Newton's gravity constant.

Let us superpose matter moving at the speed of light onto the sine gravitational wave.

If we have a mass-energy density ρ moving to the z direction at the speed of light, then its stress-energy tensor is:

       T'  =

             ρ          0           0           ρ

             0          0           0           0

             0          0           0           0
 
             ρ          0           0           ρ.

Let us trace reverse T'. We obtain

       R'  =

             κ ρ        0            0        κ ρ

             0           0            0            0

             0           0            0            0

             κ ρ        0            0        κ ρ.

Let

       η + k 

be the metric whose Ricci tensor is R' (does such a perturbation exist?).

Can we assume that the Ricci tensor R is approximately linear in perturbations of the metric? If yes, then the Ricci tensor for

       g''  =  η  +  h  +  k

is approximately

       R''  =

             0             0             0      2 κ ρ

             0             0             0             0

             0             0             0             0

             2 κ ρ      0             0             0,

and the Ricci scalar R'' is approximately zero, assuming that g'' is approximately orthogonal.

We found a metric g'', for which the variation of the Einstein-Hilbert action S is approximately zero! Since the off-diagonal corner elements of R'' are different from zero, this is not a solution of the full Einstein field equations, but it is an almost stationary solution.

Can we now use Noether's approximate theorem from our July 9, 2024 blog post to prove that energy is approximately conserved?


Noether's theorem DOES prove that energy is conserved for stationary points of the Einstein-Hilbert action S; but the next section proves that we do not have a stationary point



The claim in the headline of this section follows almost trivially if we assume that the start state and the end state are such that we can determine the energy of the system from the matter lagrangian LM alone and ignore gravity!


                        history,
                        gravitational waves
                        which get absorbed back
                ●  ---------------------------------------------  ●
            start state                          end state
            where energy                   where energy
            well defined                      well defined


For example, the end and the start state might be static matter whose density is very low.

We assume that there exists metrics g, g' for the start and the end states.

We assume that we can obtain the Ricci tensor R for the metrics g and g' by trace reversing the stress-energy tensor. Then the Ricci scalar at a location is

      R  =  κ ρ,

where κ is the Einstein constant and ρ is the mass density at the location.








In the Einstein-Hilbert action S for the history start and end states, the volume element sqrt(-det(g)) is almost exactly 1, and

       1 / (2 κ)  *  R  +  LM  =  1/2 ρ  -  ρ  =  -1/2 ρ.

At the start and end states, varying the the history of the system, as explained in the Wikipedia article about Noether's theorem, by dt at the start and -dt at the end must keep the action S constant. This implies that the integral of ρ over the entire space must be the same at the start and end states. That is, energy is conserved.

The same argument should work if the start and the end states are spherically symmetric spheres of incompressible fluid, just as in the Schwarzschild interior and exterior solutions.

Conservation of energy between the start and end states does not require us to assume anything more about the history between the states, except that there is a continuos time time symmetry which allows us to "shift" the middle part of the history by dt.

The middle part of the history can include gravitational waves which get absorbed back into the matter content M.

NOTE: No, it cannot! See the next section.

We, naturally, have to assume that there exists a stationary point of the action S. Our results with rogue variations may show that there is no stationary point?


Noether's time variation shows that the Einstein-Hilbert action does NOT allow gravitational waves?


            oscillating masses
                       <--->
                   ●\/\/\/●    variation: before emitting
                                   the wave, speed up
                                   everything by dt

 
                   ----------    quadrupole wave
                   ----------    packet, R = 0
                   ----------
                   ----------
                   ----------    variation: slow down
                                   everything by dt


The experiment:

1.   Initially, the oscillating masses are static. Let their mass-energy be M.

2.   The masses then oscillate, and radiate a large amount of the energy of the system away in a gravitational wave pulse. Let the metric be g.

3.   The masses are again static at the end of the process. The remaining mass-energy is M' < M.


The Noether time variation:

A.   Before emitting the waves, we "speed up time" by dt.

B.   When the wave packet is receding from the system, we "slow down time" by dt.


Since M > M', the action S of the mass system has a non-zero change when A and B are applied. This is because in LM, the "potential energy" is larger at A than at B.


But the change to S from the gravitational wave packet at B is zero because the Ricci tensor R is zero and:








The "slowing down of time" by dt can be done by changing the metric g of the gravitational wave infinitesimally. The formula above from the Wikipedia article about the Einstein-Hilbert action says that the Ricci scalar R stays zero in the variation.  The integral of R over the wave packet stays zero.

Assumption. In the formula for δR above, we have to assume that the error term, which is of the form ~ (δg^μν)², "behaves well" and the integral of the error over the entire wave packet is ~ dt², when we slow down time by dt. Also, the variation must be well approximated as the sum of individual variations of each gμν.


We showed that the action S changes in an infinitesimal variation of the metric g. Thus, g is not a stationary point of the action. The Einstein-Hilbert action has no solution for an experiment where a gravitational wave packet is generated.

Let us still check that the Einstein field equations really have to hold in this experiment. Could it be that the light-speed limit of signals spoils the derivation of the Einstein field equations in the Wikipedia article?

The contradiction which we derived shows that a lagrangian density is not allowed to "hide" energy. The energy visible at the start of the experiment must be visible in the lagrangian at all times. One is not allowed to sweep it under the rug. The contradiction also shows that our result above, about energy conservation in the Einstein-Hilbert action is void if there are gravitational waves: there is no stationary point of the action S at all.


The lagrangian density of electromagnetic waves: if we treat B as a "dependent" field from E, then it can handle electromagnetic waves


In the preceding section, we slowed down a packet of gravitational waves by a time dt. We cannot speed them up by dt, because then the signal would propagate faster than light. Or could we allow such an infinitesimal breach of the light speed limit?

Let us check what people say about the electromagnetic lagrangian and electromagnetic waves. We cannot find anything in the literature.

If the system is, as in the preceding section, such that it does not contain any light-speed waves in a certain period of time, then we can "speed up time" by dt during that period, and "slow down time" by dt when a wave is propagating.


(Jim Branson, 2013)

In electromagnetism, the lagrangian density for the electromagnetic wave does not "hide" its energy?








If we slow down time, the "kinetic" part B² of the wave should be reduced. Does that happen? If B is not treated as an independent field, but is born from the electric field through the formula

        B ~ dE / dt,

then B² really is reduced as time is slowed. It is like a vibrating string: if we slow down time, the kinetic energy of the string is reduced.

For a metric g, we do not have a similar division of which part of the field is "dependent". That is why the Einstein-Hilbert action cannot handle a gravitational wave.

Linearized Einstein equations would work better, because they probably have properties which are analogous to electromagnetism.


Conclusions


Let us close this long blog post.

We derived a very important result: the Einstein Hilbert action does not have a stationary point for a gravitational wave packet. That is, general relativity does not allow gravitational waves to exist – even though we have observed in nature such waves.

The reason for the failure is that the Ricci tensor R is zero in a wave, and "hides" the energy of a wave.

Fixing the problem would require making general relativity more analogous to electromagnetism, so that we can separate the "kinetic part" (B²) of the lagrangian density. 

A sine gravitational wave derived from the linearized Einstein field equations has certain components of R non-zero. It may behave more reasonably than a R = 0 wave if we "slow down time" by dt.






Above we have the energy density of a sine gravitational wave, calculated by Chris Hirata (2019). The density is ~ square of the time derivative of the metric perturbation h₊. This is reasonable, just like 1/2 m v² is the kinetic energy of an oscillating string.

It is not clear if the light-speed limit on signals affects the variational calculus at all. We might even allow infinitesimal breaches of the limit in variations, though not in histories themselves.

We did not yet get more insight into why the average stress-tensor component -T₀₀ gives the correct energy density of a gravitational wave occurring in nature.

Wednesday, July 17, 2024

Average Ricci tensor inside a gravitational wave

In our previous blog post we were stumped by the question why R₀₀ < 0 "on the average" in a gravitational wave. The Gauss-Bonnet theorem would suggest that the "average" of R₀₀ has to be zero.














Intuitively, Γkij measures how much a vector pointing to a coordinate direction i changes to the k direction when the vector is parallel transported "one unit distance" to the j direction.


















R₀₀ in a gravitational wave



Chris Hirata (2019) calculated the average R₀₀ for a gravitational wave of the h₊ polarization in the TT gauge.

We want to analyze what contributes to the 00 component of the Ricci tensor R in the case of a gravitational wave. The coordinates t, x, y, z are numbered 0, 1, 2, 3. The gravitational wave perturbation of the metric is a sine wave h₊:

       g  =

              -1              0               0              0

               0              1 + 2 h₊    0              0

               0              0               1 - 2 h₊    0

               0              0               0              1,

where

       h(t, x, y, z)  =  ε  *  sin(ω t  -  x / λ),

and ε > 0 is very small.

The 00 component of the Ricci tensor R is

       R₀₀  =  R⁰₀₀₀  +  R¹₀₁₀  +  R¹₀₂₀  +  R³₀₃₀

in terms of the Riemann tensor R. Let us study R¹₀₁₀.


               1               B
                 -------------- 
                |               |
                |               |
                 --------------
               A               2
     ^ t
     | 
      ------> x
              

An observer carries a vector, originally pointing to the time direction t, along two alternative routes from A to B, either going through 1 or 2, and measures the difference of the carried vectors at B.

The derivative of the metric component g₁₁ differs along the two routes.








Chris Hirata (2019) calculated above the values of the Christoffel symbols. The Ricci tensor component:

       R¹₀₁₀  =  dΓ¹₁₀ / dt 

                     + ...

The value Hirata calculated for R₀₀ is:







where the two last terms cause the average of R₀₀ to be negative. Let us explain where they come from.

Imagine a second square and a second observer placed on the right side of the diagram above. The first observer carries a vector of the form (a, b, 0, 0) from 2 to B, where a ≈ 1 and 0 < b is very small. Since the x metric changes between 2 and B, the value of b slightly changes.

But the second observer carries a vector (1, 0, 0, 0) from 2 to B. The contributions of the first and the second observers do not "cancel" each other for the line 2 to B! This is the reason why a Gauss-Bonnet-like theorem is not true for R₀₀. The sum of contributions for many adjacent loops is not the same as the contribution of the large loop which encloses the small loops.

Now we have the explanation for the counterintuitive fact that a wave metric g has R₀₀ which, on the average, differs from 0.

Note that static test masses placed initially at fixed coordinates (x, y, z) stay at those coordinates. There is no "focusing" of the test masses even though the average R₀₀ differs from zero!


Conjectures


Definition. Focusing of initially static test masses in a spatial volume V: the test masses start at fixed spatial coordinate positions in V at the same coordinate time, and their coordinate velocity is initially zero. Let D be the density of test masses per cubic meter, where the volume is measured in proper distances. There is focusing if

       d²D / dt² < 0,

where t is the coordinate time.


Conjecture 1. If the Ricci tensor R = 0 everywhere, then there is no focusing or defocusing of initially static test masses.


We used Conjecture 1 in our June 17, 2024 blog post where we showed that the Einstein field equations do not allow the existence of gravitational wave packets. The gravity of a wave packet focuses initially static test masses, but R = 0 inside the wave and outside it. This is a contradiction.

Conjecture 2. If R₀₀ ≥ 0 everywhere, then there is focusing of initially static test masses. If R₀₀ ≤ 0 everywhere, then there is defocusing.


Conjecture 3. Let us start from the Schwarzschild interior and exterior solution for a spherical mass M. However we change the metric g inside M, the integral of the Ricci scalar R, weighed by the "volume element" of the Einstein-Hilbert action, over M and its immediate neighborhood, does not change.


We are not sure if Conjecture 3 is true. Its intention is to say that a "total curvature" does not change inside M when the metric g does not change outside M.

The total curvature comes from the mass density, the pressure, and possible gravitational waves inside M.

Conjecture 3 implies a conservation law for "energy", which includes also gravitational waves.

Conjecture 3 implies that the Ricci scalar R cannot be zero for gravitational waves. The waves are never solutions of the full Einstein field equations.


Proof sketch of Conjecture 1: failure


UPDATE July 18, 2024: the proof sketch does not work. If the metric of x, y, z has a negative curvature, then the test masses will disperse faster than in euclidean space. If we put a focusing effect to the metric of time by slowing down time in the middle of the wave packet, we can make R₀₀ zero, but still focus initially static test masses, because the metric of time focuses them.

However, does this necessarily make R₁₁ or R₂₂ or R₃₃ non-zero?

----

                             ^
                             |
                         •      •
                  •      •      •      •
         <--    •      •  × P      •    -->
                  •      •      •      •
                         •      •       
                             |      C = cloud of dispersing
                             v            test masses m


Let us try to prove Conjecture 1.


We assume that the geometric definition in Wikipedia holds: R₀₀ at a point P = (t, x, y, z) tells us how fast test masses m sent at low velocities from P disperse. A zero curvature means that the proper volume of the test mass cloud C grows just like it would grow in euclidean space.

Let the metric be g.

Let us assume that the entire volume of space at some time coordinate t₀ is filled with distinct (i.e., intersections are empty) clouds C of dispersing test masses sent from various points P.

For each test mass m at (t₀, x, y, z), we assume another test mass m' at the same point, but whose coordinate velocity is zero. For each cloud C, we obtain a cloud C' of initially static test masses.

Assumption. The metric g does not make the clouds C' to intersect for some time after t₀.


The Assumption is true for the sine gravitational wave, because test masses m' stay at fixed coordinate positions. The Assumption probably is true for sufficiently "well behaved" metrics g.

Since the clouds C' fill the entire volume of space, and there is focusing, there has to be a cloud C' whose volume V' has a negative second time derivative d²V' / dt² at t₀.

Conjecture 4. If C' has d²V' / dt² < 0, then the volume of the corresponding cloud C has to grow slower than C would grow in euclidean space. That is, R₀₀ > 0 for the point P of the cloud C.


Conjecture 4 should be true at least for sufficiently "well behaved" metrics g. No, it is not true! See the UPDATE we wrote above.


A method of making R₀₀ zero in a metric, but still focusing initially static test masses


Let us make the spatial metric such that it has negative curvature. It is "saddle-shaped". If we send test masses from a spacetime point to all spatial directions, the volume of the test mass cloud grows faster than in euclidean metric. That is, R₀₀ < 0.

We can then reset R₀₀ to zero by making time to flow slower near the center of the system. This focuses test masses and can cancel tge negative R₀₀.

If we have initially static test masses, they are pulled toward the center by the slower flow of time there. Thus, the system does focus initially static test masses, even though R₀₀ = 0.

A gravitational wave in the TT gauge has an undulating spatial metric. Apparently, the saddle shape (negative curvature) dominates there over the spherical shape (positive curvature). This explains how, on the average, R₀₀ < 0.

One could try to reset R₀₀ in a wave packet to zero by making time to flow slower at the center of the wave packet. We are not sure if this can succeed.

Also, one could try to make R = 0 in the packet, just like the Einstein field equations require.


The paper in the link states that exact wavelike solutions for the Einstein field equations are known only for some unnatural, symmetric configurations.

Our own work in June 2024 with rogue variations revealed that general relativity has no solutions at all for "dynamic systems". Thus, it is hopeless to try to find an exact solution which also includes the source that produces the wave.


A gravitational wave packet must have R₀₀ zero but focus initially static test masses => then the energy density can be calculated with the method above


Let us have a spherically symmetric mass M. If M would be static, then the metric g around M would have the Ricci tensor R zero. In particular, R₀₀ is zero.

Let us then imagine that we make M to oscillate back and forth. Then M produces a gravitational wave, which is kind of a "retarded" field which cannot fully keep up with the accelerated movement of M.

It might be that the retarded field has the Ricci tensor R zero, just like the static field of M.

If a gravitational wave packet carries a mass-energy m, we expect it to pull test masses toward itself just like any mass of m. We conclude that the wave must focus intially static test masses.

Let us assume that the linearized Einstein equations give an approximate form for a gravitational wave packet, except that R₀₀ fluctuates while it should not. 

Chris Hirata (2019) calculated R₀₀ for a sine wave.

We can make R₀₀ approximately zero by modifying the metric of time in such a way that we put a "negative mass density" at those locations of the wave where R₀₀ > 0, and a "positive mass density" at locations where R₀₀ < 0. The added masses will travel at the speed of light.

If the new metric with added masses approximates a real gravitational wave packet, then the Chris Hirata (2019) method of calculation really gives the energy density of a gravitational wave.

There are a lot of assumptions in the reasoning above. If they happen to be correct, then we have an explanation for the energy density of a gravitational wave.

Our reasoning did not use the full Einstein field equations at all. The following are enough:

1.   the field around a spherical mass M is approximately Schwarzschild; in particular, the stretching of the spatial metric is like in Schwarzschild;

2.   when M is oscillated, the retarded field forms a gravitational wave, and the retarded field inherits the property R₀₀ = 0 from the static field;

3.   the spatial metric of the retarded field is approximately what we get from the linearized Einstein equations; we probably can derive the metric also without using any Einstein equations, by just considering retardation of the field;

4.  the metric of time in a wave packet is like for an electromagnetic wave packet of the same mass-energy E.


Note that above, the metric of time reveals the energy content of a gravitational wave to observers. The spatial metric is what actually carries the energy: we can extract the energy by using the fact that the wave stretches proper distances.


A gravitational wave must not shorten proper distances if it does slow down time as much: otherwise, we get superluminal communication


In this blog we have remarked that if a gravitational wave would shorten proper distances, and not slow down time enough, then the coordinate speed of light would exceed the speed of light in the underlying Minkowski space.

We cannot accept superluminal communication. Therefore, the metric η + h₊ produced by the linearized Einstein equations cannot be the correct metric.

Could it be that the energy density of a wave packet slows down time enough, so that the coordinate speed of light is not too fast? It does not look like that. The energy density is not large enough.

Note that in the static field of a spherical mass M, the spatial metric is always stretched, never contracted.

In our earlier blog posts, we speculated that the stretched spatial metric in a gravitational wave comes from energy shipping in an underlying field which is analogous to an electromagnetic wave.


Conclusions


We finally found a possible explanation for the fact that R₀₀ calculated from a sine gravitational wave gives the energy density of the wave.

A gravitational sine wave derived from the linearized Einstein equations has R₀₀ < 0, but its Ricci scalar R = 0. By adding the wave metric perturbation to the interior Schwarzschild solution we can reset R₀₀ to zero locally. But we cannot reset R. This suggests that gravitational wave packets always have R > 0, which reflects the fact that they carry energy.

Conjecture. A gravitational wave packet has R₀₀ = 0, but R > 0.


The sine wave has both R₀₀ and R₃₃ negative because the metric to the x and y directions is negatively curved, "saddle-shaped". A narrow conical beam of geodesics to the t or z direction spreads out faster than in euclidean space.

We have to check what the Einstein-Hilbert action says about this conjecture.

We have to think about waves which only stretch, not shorten, proper distances. How to model them?

Tuesday, July 9, 2024

Energy density of a gravitational wave

The energy density of a gravitational wave is in literature calculated by determining how much the Einstein tensor G of a wave metric g differs from zero, and interpreting the 00 component of the difference as the mass-energy density of the wave.

https://en.wikipedia.org/wiki/Binary_pulsar

Empirically, we know that this is the correct energy density. Binary pulsars proved this in the 1980s and 1990s.

Why does this method give the right energy density?


A gravitational wave must act as a source of gravity


Let us introduce an example configuration which we call a "convertible mass" M:

Let us have a mass M which is almost spherically symmetric. We can send energy radially between different layers of M using gravitational waves. The waves cannot be exactly spherically symmetric, because then they would be longitudinal waves. Presumably, we can make the waves almost spherically symmetric.

It is natural to assume that the almost Schwarzschild gravity field of M far away does not change appreciably in the process.

We conclude that the waves must act as sources of gravity, exactly according to their energy content.

One may imagine the following conversion: the mass M is a spherical vessel filled with photon gas. We use some method to convert the photons into gravitons. It is plausible that the large-scale metric inside M should remain unchanged, except that now there are mini gravitational waves bouncing around.

On June 17, 2024 we heuristically proved that the Ricci tensor component R₀₀ of a gravitational wave must differ from zero, because the wave acts as a source of gravity.


The gravitational wave "derived" from the Einstein equations


On May 21, 2024 we showed, using rogue variations, that the Einstein field equations do not seem to have a solution for any dynamic system at all. Thus, one cannot really "solve" the Einstein equations for a binary system where two masses M₁ and M₂ orbit each other.

The metric g of the quadrupole wave has to be derived in some way, but g is not a solution of the Einstein field equations.

The question is why does this "some way" produce:

1.   the right metric g for a wave, and

2.   why does the Einstein tensor G then give the correct energy density for g?


The gravitational wave metric g is derived from linearized Einstein field equations.


In the link someone (probably Chris Hirata, 2019) has derived the energy density of a gravitational wave.






We want to solve the metric g, where h is a perturbation and η is the flat metric.







The Einstein tensor G is decomposed into a linear part (1) and nonlinear parts (2), etc.







Then it is guessed that the energy density can be obtained from a nonlinear part (2) of the equation. We want to find out why this works.


The electromagnetic lagrangian














Let us analyze what is "potential energy" and what is "kinetic energy" in the lagrangian density.

If we have a static charge Q, it corresponds to a current j into the "direction of time". The magnetic field B in a static configuration is zero.

We can make a potential pit for Q by introducing an electric field E. We can then put Q into a lower electric potential, but the integral of E² grows. 

The lagrangian L(x) seems to be of the type:

       potential energy  -  kinetic energy.

Let us analyze with this lagrangian the process:

1.   we make a charge Q to move back and forth, attached to a spring;

2.   Q radiates a dipole wave which travels over a large distance;

3.   the wave is absorbed by many dipole antennae far away.


The start configuration of the action is the charge Q oscillating back and forth. The end configuration has Q almost stopped and the charges in the receiving antennae oscillating back and forth.

Why does the electromagnetic wave behave in the well known way?


The Einstein-Hilbert action in "canonical coordinates"











The May 21, 2024 rogue variations can be eliminated if we introduce canonical coordinates, against which we measure the kinetic energy of particles. That is, the kinetic energy in the Einstein-Hilbert action matter lagrangian LM is calculated against these coordinates.

Good canonical coordinates might be standard coordinates of the underlying Minkowski space, where the particles move.

The Einstein-Hilbert action should then work much better?


The energy of gravitational waves


In the formula above,

      LM  =  kinetic energy  -  potential energy.

The formula treats R as "kinetic energy". The wave g does act as a source of gravity, just like an electromagnetic wave does in LM.

What exactly is the reason why we believe that the Einstein-Hilbert action does not allow R ≠ 0 in space empty of ordinary matter? Obviously, it comes from the variational calculus, which is based on the work of Riemann, Christoffel, and Ricci-Curbastro. Could there be some error there, such that the error only is manifested when a wave moves at the speed of light?

The Einstein-Hilbert action gives a reasonable solution for a spherical static mass M. It is the Schwarzschild metric. But for the "kinetic energy" of a gravity field (= wave) we get a strange result.

It might be that the correctness of the action for a static, or an almost static field, makes the wave metric g in literature approximately correct, when it is derived in literature using whatever methods. But the action is not correct for a field moving at the speed of light?

Could it be that the Einstein-Hilbert action actually says that R ≠ 0 inside a gravitational wave?


The carpet analogy: a gravitational wave necessarily breaks the Einstein field equations


Our November 9, 2023 carpet model may cast light on the variation problem.

In a section above we introduced the "convertible mass" M, of which a part can temporarily be converted to an almost exactly spherically symmetric gravitational wave, and then M absorbs the wave again.

If the Ricci curvature R₀₀ inside M would change in the "conversion", then outside M, the Ricci curvature R₀₀ has to differ from zero, too. That is, the Einstein field equations are broken in every attempted solution of the new metric g.

Now it looks natural that the Ricci scalar R of a gravitational wave can be > R. The Einstein equations will be broken anyway.

This is just like the carpet analogy. One can press the carpet to the floor at any chosen location, but cannot press it to the floor everywhere.

If we try to make the Ricci scalar R zero in a zone which contains just gravitational waves and no matter, then, apparently, R will inevitably grow elsewhere. The variation which is used to derive the Einstein field equations, may ignore the fact that R just moves around. Just like if we have a wrinkle in a carpet which is too large: the wrinkle just moves around. Once again, we have to check the details in the variation used at:


It may be that the variation calculation is otherwise correct, except that it does not take into account that faster-than-light signals are not allowed. If they were allowed, then we could instantaneously adjust the metric throughout the entire space, to reflect the mass-energy "lost" in M into gravitational waves.

The metric g in a gravitational wave really is close to the metric which we get from the linearized Einstein equations. Furthermore, the Ricci tensor R inside the wave is not zero, thus breaking the full Einstein field equations.


C. Denson Hill and Pawel Nurowski (2017) seem to claim that gravitational waves must satisfy the full Einstein equations. Above we obtained the opposite result: a real-world gravitational wave must necessarily break the full Einstein field equations.

The history in a gravitational wave would be a stationary point of the Einstein-Hilbert action, in the allowed histories where superluminal signals are banned.


The ADM formalism may have got this right, since it starts from the Einstein-Hilbert action. Has anyone noticed that the ADM formalism breaks the Einstein field equations?


Why does an approximate solution from the linearized Einstein equations give the right Ricci tensor component R₀₀?


Above we reasoned that gravitational waves are not solutions of the full Einstein field equations. But why are they solutions of the linearized equations, and on top of that, give even the right value for R₀₀, such that it tells the energy density right?

A possible solution: we must add to the Einstein-Hilbert action the energy of gravitational waves as mass-energy?


Noether's approximate theorem


Let us have a particle moving in a potential V. Let the path x(t) of the particle be a stationary point of a newtonian action S. We then slightly modify the potential. Let the new potential be V'.

The new lagrangian density is

       kinetic energy 1/2 m x'(t)²  -  potential energy in V'.

The old path x(t) is an almost stationary point of the action, in the sense that varying x(t) by some fixed amount only slightly changes the value of the action S.

The total energy

       kinetic energy  +  potential energy

is almost conserved in the path x(t).

Noether's approximate theorem. If a history is an almost stationary point of the action S, in the sense that changing the history by some fixed amount only changes the value of S very little, then energy is almost conserved in that history.


A model with a vibrating string


                 waves     waves
         ------ ~~~~ ------ ~~~ -------  string
                    |               |
                    v               v
          gravity pulls on energy


Let us consider a vibrating string in newtonian mechanics, where newtonian gravity also pulls on the kinetic and elastic energy of the vibration.

If we ignore newtonian gravity, then the string is governed by a linear wave equation. This might be analogous to the linearized Einstein equations.

If we solve the vibration using the linear wave equation, we get a solution which approximately solves the behavior of the string. This is analogous to solving a gravitational wave g from linearized Einstein equations.

Let us then apply the full equation of the system to the approximate solution. The full equation is derived through variational calculus, and the full equation calculates how much the action S of a history changes when the history is varied. In this case, the variation we are interested in is if the string is pulled slightly lower at the locations where there is wave energy.

The full equation calculates how much the potential energy of the system is reduced if we pull the string 1 meter lower. We see that the value of the full equation measures the energy of a wave!

Similarly, the component T₀₀ from the full Einstein equation might calculate the energy of a gravitational wave g.

We still have to check the details. This may explain why the full Einstein equations produce the energy density of a gravitational wave g.


The Chris Hirata (2019) calculation of R₀₀









The metric g in the TT gauge makes the spatial metric of x and y coordinates to oscillate according to the perturbation h₊:







The average R₀₀ over many wavelengths is:








The calculation claims that R₀₀ then has a value which, on the average, is negative. In this blog we have remarked that R₀₀ measures how initially static test masses become denser as time flows. But why would a wave which only manipulates spatial distances, in an oscillating fashion, cause any change in the density of test masses in the long run?

Could it be that a wave packet manipulates the density of the test masses? But why would it decrease the density of the test masses, rather than increase?








Let us use the geodesic equation to determine the orbit of a test mass which is initially at some fixed coordinates x, y, z, and only "moves" through time, t.

Above, we have μ any of 1, 2, or 3. Let μ be 1. The rightmost term is zero because dx¹ / dt is zero. The first term on the right side of the equation is zero because Γ¹₀₀ is zero, as calculated by Hirata.

Thus, test masses stay at their initial coordinate positions x, y, z. In the metric, geodesics going to the time direction are not "focused" or "defocused".

We probably found a solution to the mystery:

R₀₀ really measures the focusing of test masses which start from the same point and move to every direction at very slow speeds. It is not enough to study test masses which start static.

Let the test masses start to every direction at a very slow coordinate speed v at a moment when the spatial metric is flat. The spatial metric inside a gravitational wave undulates. Then the test masses travel about half of the time at a coordinate speed which is

        v / (1 + ε),

and a half of the time at

        v / (1 - ε).

The average coordinate speed is larger than v.

This might mean that the average R₀₀ < 0 inside the wave.

Then we can try to "compensate" the negative average value of R₀₀ by adding another metric perturbation h', which has R₀₀ > 0, on the average.

The faraway metric in h' will act like mass pulling test masses toward the wave packet. The wave appears to have mass-energy.

Maybe we should make the compensation h' to move along the gravitational wave? Then we can probably make R₀₀ zero or almost zero inside the gravitational wave. Outside the wave packet, h' pulls test masses toward the packet.

However, can this really work? If the wave packet pulls test masses toward itself, is it possible that R₀₀ = 0 everywhere?

In two dimensions, if initially parallel geodesic lines going to the t direction approach each other in the proper distance in the x direction, then certainly R₀₀ > 0?

What about t and two spatial dimensions x, y?


Does a gravitational wave have R₀₀ = 0 in some coordinate system?


Initially static test masses never move away from their initial spatial coordinates x, y, z. Thus the coordinate system which we use is "freely falling". Above we argued, and the literature says that R₀₀ > 0 there. The wave is not a solution of the full Einstein field equations.

Somewhere someone claimed that we can choose coordinates where R = 0. We do not see how that could be possible.


Equivalence of two notions of Ricci curvature


If initially static test masses approach each other, in this blog we have claimed that then R₀₀ > 0 somewhere.

This is equivalent to saying that initially parallel lines approach each other.

But could it be that using the algebraic definition of R this is not always true?


Let us determine curvatures using vector parallel transports. Let us assume that we just have the coordinates t and x.

We divide a square 0 < t < 1 and 0 < x < 1 into small subsquares. Is it so that the rotation of a transported vector is obtained by summing the rotations for each subsquare? That seems to be the case.

Let us have a gravitational wave g which squeezes and stretches alternately the x and y directions. We put initially static test masses at various coordinates x, y, z. The masses remain at the same coordinates.

Let us parallel transport a vector pointing to the z direction around a subsquare whose sides are in the t and x directions.

If the wave packet is localized and the surrounding space has a flat metric, then no vector is rotated at all in a parallel transport. This means that all components of the Riemann tensor are zero, on the average, inside the wave packet. Then all components of R are zero, on the average.

Then R₀₀ = 0, on the average, in a wave packet, contrary to what the calculation of Chris Hirata (2019) claims.

Does the difference come from different weights that we give to volumes inside a wave packet?

The averaging by dividing a cubic volume into subcubes is the logical way to do the averaging, of course.


The Gauss-Bonnet theorem says that we obtain the geodesic curvature of the edge of a 2-dimensional region by integrating the gaussian curvature K inside the region. 

Does the analogous theorem hold for Riemann curvature?


Conclusions


Let us close this long blog post. We were not yet able to determine why R₀₀ averaged over a gravitational wave happens to give the correct energy density of the wave.

We will next look at Riemann curvature. Are the formulae with Christoffel symbols correct? Chris Hirata (2019) gets an average R₀₀ < 0 over a gravitational wave, while the Gauss-Bonnet theorem suggests that the average of R₀₀ over a wave packet should be exactly zero. Real-world gravitational waves are wave packets.