Constraints on a fifth force from the stellar orbits around the central supermassive black hole of the Milky Way

This study utilizes Markov chain Monte Carlo simulations of the S2 star's orbit around the Milky Way's central supermassive black hole to constrain the strength and range of a potential fifth force within a Yukawa gravity model, finding that the resulting parameters are consistent with previous research and compatible with recent GRAVITY Collaboration measurements of Schwarzschild precession.

The Gist

Imagine the center of our galaxy, the Milky Way, as a cosmic dance floor. In the middle of this floor sits a massive, invisible partner: a Supermassive Black Hole (Sgr A*). Dancing around it are several stars, but the most famous one is S2. S2 is a speedy dancer, whipping around the black hole in a tight, elliptical loop.

For decades, physicists have watched S2 to see if it follows the rules of General Relativity (Einstein's theory of gravity). So far, it mostly does. But scientists are always looking for a "ghost" in the machine—a tiny deviation that might hint at a Fifth Force.

What is this "Fifth Force"?

Think of the four forces we know (gravity, electromagnetism, strong nuclear, weak nuclear) as the standard rules of the universe. A "Fifth Force" would be a new, invisible rulebook.

• The Analogy: Imagine gravity is a magnet pulling the star in. A Fifth Force would be like a gentle, invisible fan blowing the star slightly away.
• Why care? If this force exists, it might explain why galaxies spin the way they do without needing "Dark Matter" (the invisible glue we assume holds them together). It might also explain "Dark Energy" (the force pushing the universe apart).

The Experiment: A Cosmic Detective Story

The authors of this paper acted like cosmic detectives. They asked: "If this invisible fan (Fifth Force) exists, how would it change S2's dance steps?"

They didn't just guess; they built a simulation.

1. The Model: They used a mathematical model called "Yukawa gravity." Think of this as a recipe that adds a specific ingredient (the Fifth Force) to the standard gravity soup. This ingredient has two settings:
   • Strength ($\delta$): How hard the invisible fan blows.
   • Range ($\lambda$): How far the fan's breeze reaches.
2. The Test: They ran the simulation three times, changing the "Range" of the fan:
   • Case 1: The breeze is short (a few hundred AU). It only affects the star when it's very close to the black hole.
   • Case 2: The breeze is medium (about the size of the whole orbit). It affects the star throughout its dance.
   • Case 3: The breeze is huge (thousands of AU). It covers the whole dance floor and beyond.

The Method: The "Guess and Check" Machine

To find the right settings for the fan, they used a computer method called MCMC (Markov Chain Monte Carlo).

• The Analogy: Imagine trying to find the perfect temperature for a shower. You turn the knob a little, check the water, turn it a bit more, check again. You do this thousands of times until you find the exact spot that feels just right.
• The computer did this with the Fifth Force settings, comparing its simulated star orbits against the actual photos and measurements astronomers have taken of S2 over the last 30 years.

What Did They Find?

The results were a mix of "nothing new found" and "interesting possibilities."

1. The Fan Gets Stronger the Further it Reaches:
   The computer found that if the Fifth Force has a long range (Case 3), it needs to be stronger to fit the data. If the range is short (Case 1), it can be very weak.

   • Simple takeaway: The longer the invisible breeze reaches, the harder it has to blow to explain the star's path.

2. The "Ghost" is Still Elusive:
   In all three cases, the computer found a "best fit" for the Fifth Force, but the uncertainty was huge.

   • The Analogy: It's like trying to hear a whisper in a hurricane. The data suggests the whisper might be there, but the noise (measurement errors) is so loud that we can't be sure. The results are consistent with "No Fifth Force at all" (General Relativity is still winning).

3. Consistency with Other Studies:
   Even though they used a slightly different mathematical recipe than other teams (like the GRAVITY Collaboration), they got similar results. This suggests that if a Fifth Force exists, its strength is likely very small (around 0.5% to 15% of normal gravity, depending on the range).

4. The "Relativity Check":
   The team also checked if their results messed up a specific prediction of Einstein's theory (called Schwarzschild precession). They found that their "Fifth Force" models didn't break Einstein's rules; they fit right within the margin of error.

The Bottom Line

This paper is a sophisticated way of saying: "We looked very hard for a new force of nature using the star S2 as our test subject. We didn't find definitive proof, but we narrowed down the possibilities."

• Did they find the Fifth Force? Not yet. The data is still too fuzzy.
• Did they learn anything? Yes. They confirmed that if this force exists, it behaves in a specific way (stronger if it reaches further), and their results match other scientists' findings.
• What's next? We need better telescopes and more precise measurements of S2. Just like you need a quieter room to hear a whisper, we need clearer data to hear the "whisper" of the Fifth Force.

In short, the universe is still keeping its secrets, but the astronomers are getting better at listening.

Technical Summary

1. Problem Statement

The paper investigates the potential existence of a fifth force (an additional repulsive or attractive interaction beyond Newtonian gravity and General Relativity) within the extreme gravitational environment of the Galactic Center (GC). Specifically, the authors aim to constrain the parameters of this force using the observed orbital dynamics of the S2 star orbiting the supermassive black hole (SMBH) Sgr A*.

Theoretical frameworks, such as massive gravity and Extended Theories of Gravity (ETG), predict a Yukawa-like correction to the Newtonian gravitational potential. This correction introduces two key parameters:

• $\delta$ (Strength): The coupling constant representing the intensity of the fifth force.
• $\lambda$ (Range): The characteristic length scale over which the force acts.

The study seeks to determine if deviations in the S2 star's orbit from General Relativity (GR) predictions can be attributed to such a force, and to establish constraints on $\delta$ and $\lambda$ across different scales relative to the S2 orbit.

2. Methodology

2.1. Theoretical Framework: Extended PPN Formalism

The authors utilize an extended Parametrized Post-Newtonian (PPN) formalism tailored for Yukawa gravity. The equation of motion for the star includes standard GR terms plus a specific Yukawa correction term:
$$\vec{\ddot{r}} = -\frac{GM\vec{r}}{r^3} + \text{GR terms} + \delta \cdot \frac{GM}{1+\delta} \left[ 1 - \left(1 + \frac{r}{\lambda}\right)e^{-r/\lambda} \right] \frac{\vec{r}}{r^3}$$
This formulation differs from standard PPN used in GR and from potentials used in recent studies by the GRAVITY Collaboration, allowing for a comparison of model dependencies. The potential used is derived from $f(R)$ theories:
$$\Phi(r) = -\frac{GM}{(1+\delta)r} \left( 1 + \delta e^{-r/\lambda} \right)$$

2.2. Data and Simulation

• Observational Data: The study fits simulated orbits to astrometric observations of the S2 star from reference [56]. Radial velocity data was excluded from the primary fit due to sparse cadence, though velocities were calculated post-fit for validation.
• Numerical Integration: The equations of motion were solved using a Runge-Kutta 5(4) algorithm (SciPy solve_ivp).
• Parameter Fitting: A Markov Chain Monte Carlo (MCMC) analysis (using the emcee ensemble sampler) was employed to estimate posterior probability distributions for the SMBH mass ($M$), six orbital elements of S2, and the fifth force parameters ($\delta, \lambda$).

2.3. Experimental Cases and Priors

The authors analyzed three distinct cases based on the range $\lambda$ relative to the S2 semi-major axis ($a$):

1. Case 1: $\lambda \sim$ a few hundred AU (deep inside the S2 orbit).
2. Case 2: $\lambda \sim$ a thousand AU (comparable to the S2 orbit size).
3. Case 3: $\lambda \sim$ several thousand AU (much larger than the S2 orbit).

Two types of priors were tested to assess sensitivity:

• Uniform Priors: Broad distributions to explore a wide parameter space.
• Gaussian Priors: Narrower distributions centered on known values to test robustness and reduce uncertainty.

2.4. Consistency Check with $f_{SP}$

The study evaluated whether the derived orbits were consistent with the Schwarzschild precession parameter ($f_{SP}$), recently measured by the GRAVITY Collaboration as $1.10 \pm 0.19$. Theoretical expressions relating $\lambda, \delta, f_{SP}$, and orbital elements were used to estimate the compatibility of the MCMC results with this measurement.

3. Key Results

3.1. Constraints on Fifth Force Parameters

The MCMC analysis yielded the following best-fit values for the strength $\delta$ across the three cases (using uniform priors):

• Case 1 ($\lambda \approx 357$ AU): $\delta \approx 0.005$
• Case 2 ($\lambda \approx 1766$ AU): $\delta \approx 0.02$
• Case 3 ($\lambda \approx 6021$ AU): $\delta \approx 0.15$

Key Trend: As the range $\lambda$ increases, the required strength $\delta$ to fit the data also increases. Conversely, the relative error ($\Delta\delta/\delta$) decreases as $\lambda$ increases.

3.2. Sensitivity to Priors

• Uniform Priors: Produced broader posterior distributions with larger uncertainties. The SMBH mass ($M$) estimates were slightly lower than standard values.
• Gaussian Priors: Resulted in significantly narrower distributions and smaller uncertainties. The SMBH mass estimates converged closer to standard values ($\sim 4.35 \times 10^6 M_\odot$).
• Sign Reversal: Under Gaussian priors, the best-fit $\delta$ values became negative ($\delta \approx -0.007, -0.03, -0.11$), suggesting an attractive rather than repulsive force in those specific fits, though the magnitude remained consistent with the uniform prior results.
• Conclusion: The results show a dependence on prior choice, particularly for $M$ and the sign of $\delta$, indicating that current data precision is insufficient to robustly constrain these parameters without strong prior assumptions.

3.3. Correlations

• $M$ vs. $\delta$: A significant correlation was found between the SMBH mass and the fifth force strength, particularly with Gaussian priors (up to ~57%). This is attributed to the rescaling factor $\delta \cdot M / (1+\delta)$ in the equation of motion, which can cause parameter degeneracy.
• $\delta$ vs. $\lambda$: No strong linear correlation was found between strength and range. However, the authors suggest a complex non-linear correlation exists, supported by theoretical relationships derived from the precession formula.

3.4. Compatibility with $f_{SP}$

The derived orbital fits in all three cases are compatible within error intervals with the GRAVITY Collaboration's measured value of $f_{SP} = 1.10 \pm 0.19$. This implies that the observed deviations from GR (if any) could potentially be explained by a fifth force, but the current uncertainties are too large to confirm this definitively.

4. Key Contributions

1. Novel Potential Application: The study applies an extended PPN formalism based on $f(R)$ gravity (Potential Eq. 4) to S2 star orbits, distinct from the potentials used in previous GRAVITY Collaboration studies.
2. Multi-Scale Analysis: It systematically investigates the fifth force across three distinct range scales (deep inside, comparable to, and larger than the S2 orbit), providing a more comprehensive view of constraints.
3. Prior Sensitivity Analysis: The paper rigorously tests the dependence of results on Uniform vs. Gaussian priors, highlighting the limitations of current data in distinguishing between model parameters without strong external constraints.
4. Model Independence: The results suggest that the constraints on the fifth force strength are largely model-independent, as similar $\delta$ values are obtained despite using different Yukawa-like potentials compared to other recent studies.

5. Significance and Conclusion

The paper demonstrates that while current astrometric data of the S2 star allows for the existence of a fifth force, the uncertainties remain too large to rule out General Relativity or definitively detect the force.

• Current Status: The constraints obtained ($\delta \sim 0.005$ to $0.15$) are consistent with previous studies but lack the precision to confirm a non-zero fifth force.
• Future Outlook: The authors conclude that future, more precise observations (specifically reducing the uncertainty in $f_{SP}$ to $\pm 0.1$ or $\pm 0.05$) are essential to tighten these constraints.
• Physical Insight: The study reinforces the idea that if a fifth force exists at the Galactic Center, its strength likely scales with its range, and its detection requires disentangling it from the SMBH mass via high-precision orbital tracking.

In summary, this work provides a robust framework for testing modified gravity theories using stellar dynamics but underscores the necessity of next-generation observational data to move from "compatible with" to "constrained by" the fifth force hypothesis.