Adiabatic evolution of asymmetric binaries on generic orbits with new fundamental fields I: characterization of gravitational wave fluxes

This paper investigates the dynamics of asymmetric binaries in scalar-tensor extensions of General Relativity by employing a new C++ code, STORM, to compute scalar-field fluxes for generic bound orbits, thereby advancing the construction of accurate waveform models for next-generation gravitational-wave detectors.

The Gist

Imagine the universe as a giant, dark ocean. For over a decade, we've been listening to the ripples on this ocean caused by massive objects crashing together. These ripples are gravitational waves, and they have been our primary way of "seeing" the universe since the first detection by LIGO.

But what if the ocean isn't just water? What if there's an invisible, ghostly mist swirling around the black holes that we can't see, but that changes how the waves ripple?

This paper is about learning how to listen for that invisible mist.

The Big Picture: The "Ghost" in the Machine

In our current understanding of physics (General Relativity), black holes are simple: they are defined only by their mass, spin, and electric charge. They are like perfect, featureless spheres.

However, many scientists suspect that the universe might be more complex. There might be extra "fields" (like the invisible mist mentioned above) that interact with gravity. If these fields exist, they would change how black holes dance and how they scream as they spiral into each other.

The authors of this paper are building a new set of ears to listen for these changes. They are specifically looking at "asymmetric binaries"—systems where a tiny, heavy object (like a small black hole or a neutron star) is spiraling into a giant supermassive black hole.

The Analogy: The Skater and the Ice Rink

To understand what they did, imagine a figure skater (the small object) spinning on a giant, slightly tilted ice rink (the massive black hole).

1. The Old Way (General Relativity): In the standard model, the ice rink is perfectly smooth. The skater's path is predictable, and the sound of their blades on the ice (the gravitational waves) follows a strict, known pattern.
2. The New Way (This Paper): The authors imagine the ice rink is covered in a layer of scented mist (the new scalar field). As the skater spins, they interact with this mist. The mist pushes back, changes the skater's speed, and alters the sound of the blades.
   • If the skater is spinning in a perfect circle, the mist might just make a low hum.
   • If the skater is wobbling (eccentric orbit) or spinning on a tilted axis (inclined orbit), the mist creates a complex, chaotic symphony of sounds.

What They Actually Did

The authors didn't just guess; they built a super-computer simulation to calculate exactly how this "mist" changes the sound.

• The "STORM" Code: They wrote a new, incredibly precise computer program called STORM (Scalar Tensor Orbital Radiation from EMRIs). Think of this as a high-fidelity audio synthesizer. Instead of just playing a simple beep, it can generate the complex, multi-layered music of a black hole spiraling through this invisible mist.
• The "Generic" Orbit: Previous studies mostly looked at simple, flat, circular orbits. But in reality, these cosmic dances are messy. The skaters wobble, tilt, and move in weird loops. This paper is the first to map out the "sound" of these messy, 3D, wobbly orbits interacting with the mist.
• The "Flux": They calculated how much energy is lost to the mist as the skater spirals inward. This energy loss is what eventually causes the two objects to crash together.

The Key Findings (The "Tune")

After running thousands of simulations, they found some fascinating patterns:

1. The "Dipole" Dominance: The strongest signal from the mist comes from a specific type of vibration (called the dipole). It's like the bass note in a song; it's the loudest and most important part.
2. Wobble Matters: The more the orbit is tilted or the more it wobbles (eccentricity), the more complex the "song" becomes. New notes (harmonics) appear that weren't there in the simple, flat orbits.
3. Spin Doesn't Change Much: Surprisingly, how fast the giant black hole is spinning (its "spin") doesn't change the mist's song as much as the wobble of the orbit does. The shape of the dance matters more than the speed of the dancer.
4. The "Golden" Systems: They found that certain combinations of tilt and wobble make the "mist signal" much louder. Future telescopes (like LISA, a space-based gravitational wave detector) should look for these specific systems because they are the best places to catch a glimpse of this new physics.

Why This Matters

We are currently in the "listening" phase of gravitational wave astronomy. We hear the crash, but we don't fully understand the music yet.

This paper provides the sheet music for a new kind of song. If the next generation of detectors (like LISA) hears a black hole spiraling in a way that matches this "misty" sheet music, it will be a revolutionary discovery. It would prove that Einstein's theory of gravity is incomplete and that there are new, fundamental forces at play in the universe.

In short: They built a better microphone and a new set of musical scores so that when the universe sings a new song, we'll finally know what it's saying.

Technical Summary

Here is a detailed technical summary of the paper "Adiabatic evolution of asymmetric binaries on generic orbits with new fundamental fields I: characterization of gravitational wave fluxes" by Gliorio et al.

1. Problem Statement

The paper addresses the challenge of modeling asymmetric binary systems (specifically Extreme and Intermediate Mass Ratio Inspirals, or EMRIs/IMRIs) within extensions of General Relativity (GR) that include a massless scalar field non-minimally coupled to gravity.

• Context: While GR has been tested extensively, next-generation gravitational wave (GW) detectors (like LISA and the Einstein Telescope) require high-precision waveform models to test fundamental physics. Asymmetric binaries are ideal laboratories for this due to their long signal duration and rich harmonic structure.
• Gap: Previous studies in modified gravity (beyond GR) were limited to specific orbital configurations, such as equatorial eccentric or spherical inclined orbits. There was no comprehensive framework for fully generic bound orbits (simultaneously eccentric and inclined) that could compute the full spectrum of scalar radiation fluxes.
• Goal: To characterize the scalar energy and angular momentum fluxes emitted by a scalar-charged compact object inspiraling into a Kerr black hole on a generic orbit, providing the necessary data to build accurate adiabatic waveform models.

2. Methodology

Theoretical Framework

The authors utilize an Effective Field Theory (EFT) approach tailored for the small mass ratio $\epsilon = m_p/M \ll 1$.

• Background: The primary black hole is treated as a Kerr metric (justified because deviations from Kerr are negligible at leading order in the mass ratio for the considered theories). The background scalar field is constant due to no-hair theorems.
• Perturbations: The secondary object is modeled as a point particle with a scalar charge $d$, which acts as the sole new parameter distinguishing the theory from GR. The scalar field perturbation $\psi^{(1)}$ is sourced by this charge.
• Equations: The scalar field equation is a wave equation on the Kerr background (Eq. 11), decoupled from metric perturbations at leading order. The source term is derived from the particle's worldline.

Numerical Implementation: STORM

The authors developed a new C++ code named STORM (Scalar Tensor Orbital Radiation from EMRIs) to solve the perturbation equations.

• Arbitrary Precision: The code uses the Boost multiprecision library to handle numerical stability in the extreme parameter space of EMRIs.
• Geodesic Integration: It integrates the geodesic motion of the secondary on generic Kerr orbits using elliptic functions (Jacobi elliptic functions) and Mino time, parametrized by semi-latus rectum ($p$), eccentricity ($e$), and inclination ($\theta_{\text{inc}}$).
• Mode Decomposition: The scalar field is decomposed into spin-weighted spheroidal harmonics. The code solves:
  • Angular Teukolsky Equation: Using a continued-fraction method and spectral expansion to find eigenvalues and eigenfunctions.
  • Radial Teukolsky Equation: Using a hybrid numerical-analytical scheme based on the Mano-Suzuki-Takasagi (MST) formalism to set boundary conditions at the horizon and infinity, followed by numerical integration in the intermediate region.
• Flux Calculation: The code computes the asymptotic amplitudes of the scalar field modes to derive the energy ($\dot{E}$) and angular momentum ($\dot{L}$) fluxes at both the event horizon and spatial infinity.

3. Key Contributions

1. Generalization to Generic Orbits: This is the first study within this EFT framework to compute scalar fluxes for fully generic bound orbits, simultaneously incorporating eccentricity and inclination.
2. Development of STORM: The release of a robust, arbitrary-precision numerical code capable of evolving perturbations along generic Kerr geodesics, applicable to both scalar and tensor modes.
3. Comprehensive Flux Characterization: The authors mapped the complete harmonic structure of scalar fluxes across a broad parameter space, analyzing the dependence on orbital geometry ($p, e, \theta_{\text{inc}}$) and black hole spin ($a$).
4. Harmonic Hierarchy Analysis: They identified the specific hierarchy of modes $(\ell, m, n, k)$ that dominate the radiation, providing practical guidelines for constructing flux grids for waveform generation.

4. Key Results

Mode Hierarchy and Structure

• Dominance of Dipole Radiation: The radiation is predominantly dipolar ($\ell=1$). The $\ell=m=1$ mode provides the leading contribution to the total flux.
• Multipole Suppression: Higher multipoles ($\ell > 1$) are suppressed but remain non-negligible, especially for highly eccentric and inclined configurations where the spectral content broadens.
• Harmonic Indices ($n, k$):
  • Eccentricity ($e$): Increasing eccentricity significantly broadens the distribution of the radial harmonic index $n$. For low $e$, the peak is at $n=0$; for high $e$, the peak shifts to $|n| \neq 0$.
  • Inclination ($\theta_{\text{inc}}$): Increasing inclination broadens the distribution of the polar harmonic index $k$. While low-$k$ modes remain dominant, the number of contributing harmonics increases with inclination.
  • Peak Locations: For dominant modes ($\ell=m$), the peak in the $k$-spectrum is at $k=0$, while the peak in the $n$-spectrum shifts to higher values as $\ell$ increases. For subdominant modes ($\ell \neq m$), the peak in $k$ occurs at $k = \ell - m$.

Parameter Dependence

• Orbital Radius ($p$): As expected, the total flux decreases as the orbital separation increases.
• Black Hole Spin ($a$): The scalar flux harmonic structure is remarkably insensitive to the black hole spin. While spin induces modest shifts in the location of spectral peaks at high eccentricity, the overall distribution is far more strongly influenced by eccentricity and inclination than by spin.
• Flux Magnitude: The total scalar energy flux varies only moderately across the sampled parameter space (spanning orders of magnitude in $p$ and $e$), suggesting that the cumulative effect of scalar radiation over the full inspiral is the critical factor for detectability.

Tables and Data

The paper provides extensive tables (Tables I–V) listing:

• Total fluxes summed over all harmonics for various $(\ell, m)$ modes.
• Locations of spectral peaks ($n_{\text{peak}}, k_{\text{peak}}$) for different orbital configurations.
• Total flux values for a grid of $a, e, \theta_{\text{inc}}, p$ values, serving as a reference for future waveform modelers.

5. Significance and Outlook

• Foundation for Waveform Models: This work provides the essential "building blocks" (flux grids and amplitude data) required to construct accurate adiabatic inspiral waveforms for scalar-charged EMRIs/IMRIs. These waveforms are critical for the LISA mission and future ground-based detectors to test deviations from GR.
• Precision Tests of Fundamental Physics: By quantifying how scalar charge imprints on the GW signal, the study enables future parameter estimation studies to constrain the scalar charge $d$ and the coupling scale of modified gravity theories (e.g., Gauss-Bonnet gravity).
• Future Work: The authors plan to integrate these flux grids into inspiral evolution frameworks (like FEW) to simulate the full orbital evolution and perform parameter estimation studies. They also aim to extend the analysis to the merger phase and larger mass ratios.

In summary, this paper bridges a critical gap in theoretical gravitational physics by providing the first high-precision, generic-orbit characterization of scalar radiation in modified gravity, laying the groundwork for the next era of precision gravitational wave astronomy.