Bjerrum-Bohr, Donoghue, and Holstein (2002) calculate some kind of vertex correction diagrams for gravity.
What is the energy density of the gravity field?
^ E
|
----------- + capacitor
----------- - plates
If we have a static electric field E, we can extract energy from it locally by using capacitor plates which have opposite charges.
• m
------------- scaffolding
| |
●
M
To extract energy from a static gravity field, using a small mass m, we must use a scaffolding which encloses the mass M. Then we can lower a small mass m and extract energy. Energy cannot be extracted locally.
This is the old problem about where is the energy of a gravity field located. Gravitational waves certainly carry energy. But is there energy in a static gravity field and what is the energy density, if not zero?
Electromagnetic/gravitational waves: a quantitative model
F M/Q M/Q -F
------> ● ● <------
a d
----------------- ----------------
accelerate decelerate
time t time t'
1. First we accelerate to the right a particle of a charge Q or a mass M, with a constant force F for a time t. The particle moves a distance a.
2. Then we decelerate it with a constant force -F. It stops after a distance d, after some time t'.
3. The work lost,
W = F (a - d)
escaped in radiation. The momentum
p = F (t - t')
escaped in radiation.
Since gravitational waves carry 16X the energy of analogous electromagnetic waves, the difference a - d has to be 16X for gravity.
\ rubber plate
\
● ----> F
/
/
Let us analyze the "effective inertia" of the particle in the process. Let us first look at a "rubber plate" model of the field of the particle. In phase 1, the far field of the particle does not have time to react. The effective inertia of the particle may be reduced because of that.
But the effective inertia also increases because the rubber plate in the diagram pulls the particle to the left.
\ rubber plate
\
-F <----- ●
/
/
In phase 2, the bending of the rubber plate actually helps -F to pull the particle to the left. The effective inertia is reduced. This explains why d < a.
How do we obtain a 16X radiated energy?
A. If the rubber plate has a zero mass and resists stretching very much, then it carries away a lot of energy.
B. If the rubber plate has mass, and and is easily stretched, then it carries away a lot of energy.
In the case of gravity, A is a more beautiful option. Assigning mass-energy to a static gravity field does not look nice.
The effects of A and B on the effective inertia of an electric charge
Let us estimate the effects of A and B of the previous section when an electric charge Q is accelerated. The effective inertia of the particle is reduced because the far field does not follow it, but the electric field lines resist stretching, and add effective inertia.
Let a mildly relativistic electron pass a proton at a distance R. In the previous blog post we calculated from the Larmor formula that the radiated energy is
W = 4/3 * 1 / (4 π ε₀)³ * e⁶ / c⁴ * 1 / R³
* 1 / me².
The transit past the proton lasts a time
t ~ 2 R / c.
The pulling force feels the inertia of the electron, and it also must do the work to create the radiation.
Let us assume that the electron is static and that the proton flies past it in the time t. A typical acceleration of the electron is
ae ~ 1 / (4 π ε₀) * e² / R² * 1 / me.
The distance r which the electron moves is very crudely
r ~ 1/2 ae t²
= 2 * 1 / (4 π ε₀) * e² / R² * 1 / me
* R² / c²
= 2 * 1 / (4 π ε₀) * e² / c² * 1 / me.
The force required to create the radiation is
FW ~ W / r
= 1/2 * 4/3 * 1 / (4 π ε₀)² * e⁴ / c²
* 1 / R³ * 1 / me.
Let us compare this to the typical pulling force by the proton on the electron:
Fe = 1 / (4 π ε₀) * e² / R²,
FW / Fe ~ 2/3 * 1 / (4 π ε₀) * e² / c² * 1 / R
* 1 / me
≈ 2/3 re / R.
We see that when R is roughly re, the creation of radiation involves a force roughly the same as the pulling force of proton.
The far field of the electron does not have time to follow the electron in the process. This reduces the effective inertia of the electron by
~ re / R.
We see that the the inertia effect of the radiation is roughly of the same magnitude as of the electric field lagging behind. What is the role and the relation of these two effects?
mildly relativistic
proton+
● ---->
○------------- • -------------○ far field
e-
If the "spring" connecting the far field to the electron is very loose, then the the effective inertia reduction from the far field lagging behind dominates.
Question. Is the inertia of the far field relevant at all? The electron communicates with it through the electromagnetic field. Could it be that the inertia of the far field can be ignored? Is the inertia of the electron always me, even if the far field does not have time to react?
In a hydrogen atom, the effective mass of the electron is not reduced by its far field, which, presumably, cannot keep up with the rapid orbit of the electron around the proton.
On October 22, 2025 we presented a hypothesis that the vertex correction in QED comes from almost-bremsstahlung which almost can escape from the electron, but cannot become a real photon because it does not have h f of energy. Then we can say that the vertex correction in QED does not have anything to do with the far field, but is solely about bremsstrahlung. Also, there would be no purely classical vertex correction.
The vertex correction in gravity
k
~~~~~~~ graviton
p / \ p + q
e- ----------------------------------------
|
| q
| photon or graviton
|
proton+ ----------------------------------------
The virtual graviton k probably couples to the mass me of the electron, and to the kinetic energy of the electron. Then the integral looks much like the QED vertex correction. Can we renormalize it?
The propagator for the graviton should be the impulse response of the gravity field to a Dirac delta disturbance. What is it like?
Since gravity is nonlinear, we do not even know if there exists a solution which corresponds to a Dirac delta disturbance of the field.
The scattering process above is dynamic. In May 2024 we tentatively proved that the Einstein field equations have no dynamic solutions at all.
To make progress, we must try to estimate what the "correct" propagator for the graviton might be like.
The coupling constant for gravity, and the graviton propagator
In a pure momentum exchange, like a planet passing by a star, the graviton propagator, and the associated coupling constants, look like the photon propagator in QED. That is, Newton's force is like Coulomb's force:
F ~ m M / r².
But in an emission of gravitational waves, the power is 16-fold compared to electromagnetism. The coupling constant should be tuned to reflect this?
In this blog we have had the idea of "teleportation" to describe a gravitational wave. If we have a binary black hole, there are large dynamic changes in the metric around it. One can use the changes in the metric to drain energy from the gravity field at a 16-fold pace, compared to the analogous electric field.
A gravitational wave "teleports" the metric to a distant place, attenuating it by 1 / R², where R is the distance, and only preserving components which are normal to the radius between the observer and the binary black hole, since there are no longitudinal waves. In the teleportation interpretation, the propagator is similar to electromagnetism, but the coupling is complex, since the metric is complicated close to a binary black hole.
Scattering processes are short-distance phenomena. We do not need teleportation in them.
Hypothesis. The Green's function for the gravity field is a "hammer hit" to each of the 16 components of the metric tensor separately. Each component is capable of carrying away the energy of a hammer hit to the analogous electromagnetic field of the analogous electric charge. This explains why gravitational waves carry 16 times the energy of equivalent electromagnetic waves. The propagator for each component is like the one for the photon.
The following could be a counterargument to the hypothesis: off-diagonal components of the metric g are not independent. E.g., g₀₁ = g₁₀. But it could be that such a component receives a double share of the hammer energy.
What about a hammer hitting just the component g₀₀, i.e., creating a transient mass? The problem may be that the homogeneous gravity differential equation, that is, the vacuum equation, does not have solutions for such a configuration. That is, there is no Green's function.
A hammer strike in QED is allowed to break the rules: the 4-momentum can be arbitrary. Could it be that the hammer strike in gravity can do the same?
Corollary. Why does gravity imitate the Coulomb force in a system of orbiting masses? Because the only component relevant to the orbital motion is the metric of time. Its propagator is just like the photon propagator.
What is a Green's function for a field with many interconnected components?
The Einstein field equations tie together the 16 components of the metric tensor. The equation in vacuum can be called the homogeneous equation. How do we disturb the homogeneous equation with a Dirac delta function, and what is the Green's function for it?
For a scalar field this is clear. But what about a matrix equation?
Worldline quantum field theory of Jakobsen et al. (2021)
Jakobsen, Mogull, Plefka, Steinhoff (2021) calculate the classical bremsstrahlung in gravity.
The authors of the paper have developed a worldline quantum field theory, to be able to calculate the bremsstrahlung waveform perturbatively.
The authors introduce retarded propagators. This looks promising!
Above is another paper by Mogull, Plefka, and Steinhoff (2021).
The setting is that a gravitational wave h is emitted by a particle, or absorbed by a particle.
h particle
~~~~~~
~~~~~~ • z(τ)
~~~~~~
---> τ proper time
z location
We may assume that the particle is initially static at the origin (0, 0, 0) of the spatial coordinates. Its position vector is denoted by z(τ), where τ is the proper time of the particle.
The authors calculate a Fourier decomposition of the metric perturbation (gravitational wave) hμν.
Let us assume that the particle only moves on the y axis. Then the second Fourier decomposition above is simply the Fourier decomposition of the real-valued function y(τ).
Above is the position space propagator for the graviton.
The authors give a strange "propagator" for the position of the particle. The propagator should be used to calculate the probability amplitude of a particle moving from one state to another. Why would the amplitude be ~ τ₁ - τ₂ ? Let us take that at a face value and check how the authors will use this "propagator".
#################
An aside: a loop with two photons breaks the classical limit in QED?
e- ----------------------------------
| | p - k virtual
| q + k | photons
e+ ----------------------------------
k = arbitrary 4-momentum
in the loop
Let us assume that we have been able to calculate the Feynman integral for the above loopy diagram – possibly renormalizing the result.
Let us upscale the system, making the electron mass me N²-fold and the elementary charge e N-fold, where N is a large positive number. Then 4-momenta scale by the factor N² and the coupling constant e² / (4 π) scales by the factor N².
We consider a small part of the Feynman integral:
d⁴k [photon and electron propagators]
* [coupling constant for 4 vertices].
1. The 4-momentum volume element d⁴k scales by (N²)⁴ = N⁸.
2. The coupling constants e² / (4 π) scale by (N²)⁴ = N⁸.
3. A photon propagator scales by 1 / (N²)² = N⁴.
4. An electron propagator scales by 1 / N².
5. The combined scaling of the propagators is 1 / N⁴ * 1 / N² * 1 / N⁴ * 1 / N² = 1 / N¹².
6. The total scaling of the integral is
N⁸ * N⁸ / N¹² = N⁴.
e- ----------------------------
| q
e+ ----------------------------
The scaling of the simplest tree level scattering diagram is from the coupling constants (N²)² and from the photon propagator 1 / (N²)². That is, the scaling is 1.
The scaling of the loopy diagram is N⁴. The loopy diagram would dominate for large N. This only makes sense if the integral of the loopy diagram is zero.
M. Consoli (1979) calculates box diagram integrals (Figure 7. (a) and (b) in the paper). Consoli does not mention that the integrals would be zero.
We conclude that the classical limit of the Feynman box integral is nonsensical. This suggests that the integral is nonsensical also if we use the normal electron mass me and charge e. What is the problem? Apparently, the Feynman integral is not the right way to analyze the "fine details" of a particle orbit in a Coulomb field. The simplest tree level diagram works very well, but a loop with two virtual photon exchanges is not sensible.
*** WORK IN PROGRESS ***








