A note on Gravitational radiation in generalized Brans-Dicke theory: compact binary systems

This paper corrects a previous error regarding the lower bound on the Brans-Dicke coupling parameter $\omega_0$ in generalized Brans-Dicke-$f(R)$ theories by reanalyzing gravitational radiation data from the binary system PSR J1012+5307, thereby establishing revised constraints that highlight the distinct role of $\omega_0$ in the presence of both massless and massive scalar fields.

The Gist

Imagine the universe as a giant, invisible trampoline. In our standard understanding of physics (Einstein's General Relativity), heavy objects like stars and black holes warp this trampoline, creating ripples we call gravitational waves.

But there's a popular "what if" theory called Brans-Dicke theory. It suggests that the trampoline isn't just made of fabric; it's also influenced by an invisible "ghost" field (a scalar field) that can stretch or shrink the fabric on its own.

Recently, a group of scientists (the authors of the paper you referenced) tried to mix this "ghost field" idea with another complex theory called $f(R)$ gravity. They claimed to have found a new, super-strict rule about how strong this ghost field can be. They said, "We've calculated that the ghost field must be incredibly weak, much weaker than we thought!"

Here is the twist: The authors of this note (Jesus, Velten, and Toniato) are saying, "Hold on a minute. We think you made a tiny mistake in your math code, and it completely flipped your results upside down."

Here is the breakdown of what happened, using some everyday analogies:

1. The Two Ghosts

In this specific theory, there are actually two invisible forces at play:

• The Classic Ghost ($\phi$): This is the original Brans-Dicke field. Think of it as a gentle breeze that can change the weight of the trampoline.
• The Heavy Ghost ($\Phi$): This comes from the $f(R)$ part. Think of this as a heavy backpack the trampoline is wearing. If the backpack is light, the trampoline wobbles a lot. If the backpack is heavy, the trampoline acts more like normal.

The scientists wanted to know: How strong is the breeze (the first ghost) allowed to be, given how heavy the backpack (the second ghost) is?

2. The "Typo" That Flipped the World

The previous study (Ref [1]) looked at data from a cosmic dance between two dead stars (a binary pulsar system called PSR J1012+5307). They measured how fast the stars were spiraling toward each other.

They claimed that for the theory to match reality, the "breeze" parameter ($\omega_0$) had to be huge (over 6 million). In their graph, they showed a "Forbidden Zone" where the theory would break the laws of physics.

The Correction:
The new authors looked at the computer code the previous team used. They found a tiny typo.

• The Mistake: The code had a math exponent of $+3/2$ where it should have been $-3/2$.
• The Analogy: Imagine you are baking a cake and the recipe says "Add 3/2 cups of sugar" but you accidentally typed "Subtract 3/2 cups." The result is a very different cake.
• The Result: Because of this sign error, the previous team's graph was inverted.
  • Their Graph: Showed that low values of the "breeze" were allowed, but high values were forbidden.
  • The Real Graph: Shows that high values are allowed, but low values are forbidden.

3. What the New Graph Actually Means

The new authors redrew the map (Figure 1b and 1c in the paper). Here is what they found:

• If the "Backpack" is very light: The "breeze" (the Brans-Dicke field) must be extremely weak (the number $\omega_0$ must be huge, over 6 million). This is actually stricter than the old Solar System tests, but for a different reason.
• If the "Backpack" gets heavier: The rules relax. The "breeze" can be stronger.
• The Limit: There is a point where the "Backpack" gets so heavy that the theory stops behaving like this special "Brans-Dicke" version and just becomes the standard Einstein gravity we already know.

4. Why This Matters

You might ask, "Why does a typo in a math code matter?"

• Scientific Integrity: It shows that even in complex, high-level physics, a single minus sign can change a conclusion from "This theory is dead" to "This theory is alive but constrained."
• The Real Limit: The previous authors claimed their result was "two orders of magnitude stricter" than anything we knew from the Solar System. The new authors say, "Actually, your result was based on a glitch. The real limit is consistent with what we expect, but it depends heavily on how heavy that second 'ghost' field is."

The Bottom Line

The paper is a correction note. It's like a referee blowing the whistle in a soccer game and saying, "Wait, the goal was offside because the linesman made a mistake."

The authors of this note fixed the math, flipped the graph, and gave us the correct boundaries for how this specific theory of gravity can exist in our universe. They confirmed that while the theory is still interesting, the previous claim of a "super-strict" new limit was an illusion caused by a coding error.

Technical Summary

1. Problem Statement

The paper addresses a discrepancy in the constraints derived for the coupling parameter $\omega_0$ in Brans-Dicke-$f(R)$ (BDf) theory, a "first generation" scalar-tensor theory of gravity.

• Context: BDf theory extends General Relativity (GR) by introducing two scalar fields: a massless Brans-Dicke field ($\phi$) and an effectively massive geometric field ($\Phi$, derived from $f(R)$).
• The Conflict: A recent study (Ref. [1]) claimed to derive a lower bound on $\omega_0$ of $\omega_0 > 6.09723 \times 10^6$ for massless BD fields with a geometric field mass $m_f < 10^{-29}$ GeV. The authors of Ref. [1] asserted this bound was "two orders of magnitude stricter" than constraints derived from Solar System tests (specifically the Cassini mission, which sets $\omega_0 \gtrsim 40,000$).
• The Issue: The authors of this note argue that the result in Ref. [1] is physically inconsistent and mathematically flawed, particularly regarding how the massive field $m_f$ influences the constraints and the behavior of the theory in the limit where it should recover standard Brans-Dicke results.

2. Methodology

The authors re-analyze the gravitational radiation emitted by compact binary systems (specifically the pulsar system PSR J1012+5307) within the BDf framework.

• Theoretical Framework: They utilize the action for BDf theory involving the scalar field $\phi$ and the function $f(R)$. They examine the trace of the field equations and the evolution equation for the scalar field.
• Key Equation: The analysis focuses on the fractional period decay ($\dot{T}$) of a binary system due to gravitational wave emission. The authors utilize the expression provided in Ref. [1] (Eq. 120 in Ref. [1], labeled as Eq. 4 in this note), which includes terms dependent on the Heaviside function $\Theta$, the geometric mass $m_f$, and the coupling $\omega_0$.
• Reproduction and Verification:
  1. The authors reproduced the plot from Ref. [1] (Fig. 1a) using the same data and parameters.
  2. They performed an independent calculation of the bounds on $\omega_0$ as a function of the geometric mass $m_f$.
  3. Error Identification: By comparing their results with the original, they identified a specific numerical error in the code used by Ref. [1]. The term $(1 - T^2 m_f^2 / 4\pi^2)^{3/2}$ in the decay rate equation was incorrectly implemented with an exponent of $-3/2$ in the original code, rather than the correct $3/2$.

3. Key Contributions

• Correction of a Numerical Error: The primary contribution is the identification of a typo in the numerical implementation of the period decay formula in Ref. [1]. This error inverted the relationship between the geometric mass $m_f$ and the coupling parameter $\omega_0$.
• Re-evaluation of Constraints: The authors provide corrected bounds on $\omega_0$ that restore physical consistency.
• Clarification of Physical Limits: They demonstrate that the corrected model correctly recovers the standard Brans-Dicke limit (and the Cassini bound) as the mass $m_f$ increases, whereas the original result in Ref. [1] failed to do so logically.

4. Results

The corrected analysis yields a different relationship between the coupling parameter $\omega_0$ and the geometric scalar field mass $m_f$:

• Inversion of Excluded Regions: The original plot (Fig. 1a) showed a region of exclusion that suggested extremely high $\omega_0$ values were required for low masses. The corrected plot (Fig. 1b) inverts this, showing that the constraints behave differently.
• Specific Bounds:
  • For very low geometric masses ($m_f \gtrsim 1 \times 10^{-29}$ GeV), the bound is $\omega_0 \gtrsim 6 \times 10^6$.
  • However, as the mass increases to $m_f \gtrsim 7.6 \times 10^{-29}$ GeV, the bound relaxes to $\omega_0 \gtrsim 40,000$.
• Consistency with Solar System Data: The corrected result shows that for larger mass values (approaching the traditional Brans-Dicke limit), the constraint on $\omega_0$ weakens, eventually recovering the well-established Cassini limit of $\omega_0 \gtrsim 40,000$.
• Upper Limit on Mass: The analysis notes an upper bound on $m_f$ ($< 8 \times 10^{-29}$ GeV) imposed by the orbital period of PSR J1012+5307, beyond which the BD limit cannot be investigated using this specific binary system.

5. Significance

• Restoration of Physical Consistency: The paper corrects a significant overestimation of the constraints on BDf theory. The claim that binary pulsar data provides a bound "two orders of magnitude stricter" than Solar System tests is invalidated by this correction.
• Methodological Rigor: It highlights the importance of verifying numerical codes in theoretical astrophysics, as a simple sign or exponent error can lead to drastically different physical conclusions.
• Theoretical Implications: The corrected results clarify the role of the massive geometric field in scalar-tensor theories. It confirms that the influence of the massive field is distinct from the massless BD field and that the theory behaves consistently with known limits (Cassini) when the massive field is sufficiently heavy.
• Impact on Future Research: Researchers relying on the previous bounds for BDf theory constraints must update their models using the corrected relationship between $\omega_0$ and $m_f$ presented in this note.