An efficient higher-order WKB code for quasinormal modes and greybody factors

This paper presents an optimized and significantly faster version of a higher-order WKB Mathematica code for computing quasinormal modes and greybody factors, which achieves substantial speed gains by expanding the effective potential in a Taylor series around its maximum rather than evaluating the full analytic WKB formula for each specific potential.

The Gist

Imagine a black hole not as a silent, dark vacuum, but as a giant, cosmic bell. When you ring this bell (by dropping a star into it or smashing two black holes together), it doesn't just sit there; it "rings" with specific tones. In physics, these tones are called Quasinormal Modes (QNMs). They tell us the black hole's size, shape, and how stable it is.

There's also a related concept called Grey-body Factors. Think of the black hole's gravity as a steep hill. If you roll a ball (a wave of light or gravity) toward the hill, some of it rolls back, and some rolls over the top. The "Grey-body Factor" is simply the percentage of the ball that makes it over the hill.

The Problem: The Old Calculator Was Too Slow

For years, physicists have used a mathematical tool called the WKB approximation to predict these tones and transmission rates. It's like a recipe for calculating how a bell rings without having to build the bell first.

However, the old version of this "recipe" had a major flaw: it was incredibly slow and clumsy.

• The Old Way: Imagine trying to bake a complex cake. The old method required you to write out the entire mathematical formula for every single ingredient, every single step, and every possible variation of the recipe before you could even turn on the oven. If the recipe was complicated (like a black hole with weird physics), the computer would spend hours just writing down the instructions, often getting lost in the details.
• The Result: For many interesting black holes, the calculation would take hours, or the computer would simply give up because the math was too messy.

The Solution: A Smarter, Faster Code

The authors of this paper (Konoplya, Matyjasek, and Zhidenko) have released a brand-new, super-charged version of this calculator. Here is how they made it faster, using a simple analogy:

The "Taylor Series" Trick
Instead of trying to write out the entire, massive, complicated formula for the whole universe, the new code does something much smarter. It zooms in on the very top of the "hill" (the peak of the potential barrier) and asks: "What does the hill look like right here, right now?"

• The Old Way: Trying to map the entire mountain range from space before you can climb it.
• The New Way: Standing at the very peak, taking a few quick measurements of the slope right under your feet, and using those local measurements to predict the rest.

By switching from "writing the whole book" to "taking a quick snapshot of the peak," the new code skips the boring, time-consuming parts.

The Results: From Hours to a Blink

The paper presents a table comparing the old code to the new one. The difference is staggering:

• Old Code: Could take hours to calculate a single tone for a complex black hole. Sometimes, it would lose precision (like a blurry photo) because it had to do so many calculations.
• New Code: Does the same job in a fraction of a second. It's like switching from walking across the country to taking a teleportation beam.

Even for the most complicated black holes (the ones with "non-rational functions," which are just fancy math words for "very messy shapes"), the new code is thousands of times faster.

Why Does This Matter?

1. Speed: Physicists can now test hundreds of different theories about black holes in the time it used to take to test one.
2. Accuracy: Because the code is faster, it can use more advanced math (up to the 16th order of precision) without crashing the computer. This means the predictions are sharper and more reliable.
3. Connection: The code also helps link the "ringing" of the black hole (QNMs) to the "transmission" of waves (Grey-body factors). It's like realizing that if you know the exact pitch of a bell, you can instantly know how much sound escapes through a window.

In a Nutshell

This paper is about giving physicists a Formula 1 race car to replace their horse-drawn carriage. They didn't invent a new way to race (the WKB method is still the same), but they built a much better engine. Now, they can explore the mysterious physics of black holes at the speed of light, rather than waiting for the sun to rise and set while the computer thinks.

Technical Summary

1. Problem Statement

Quasinormal modes (QNMs) and grey-body factors (GBFs) are fundamental characteristics of black hole spacetimes, governing the ringdown phase of perturbations and the transmission probability of waves through effective potential barriers, respectively. The Wentzel–Kramers–Brillouin (WKB) approximation is the standard tool for calculating these quantities due to its efficiency.

However, existing implementations of the WKB method (specifically the version up to the 13th order with Padé resummation) face significant computational bottlenecks:

• Symbolic Differentiation: Previous codes relied on constructing full analytic expressions for high-order derivatives of the effective potential. For complex potentials involving non-rational functions or complicated metric structures, symbolic differentiation becomes extremely time-consuming or computationally infeasible.
• Scalability: As the order of the WKB expansion increases, the complexity of the analytic expressions grows exponentially, rendering high-precision calculations impractical for many modern black hole models (e.g., regular black holes, modified gravity theories).

2. Methodology

The authors present an updated and optimized Mathematica® package that fundamentally changes how the WKB expansion is computed. The core methodological shift involves moving from symbolic to numerical evaluation of derivatives.

• Taylor Series Expansion at the Maximum: Instead of deriving analytic formulas for the $k$-th derivative of the potential $U(x)$, the new code:
  1. Numerically determines the location of the potential maximum, $x_m$.
  2. Expands the effective potential $U(x)$ in a Taylor series around $x_m$.
  3. Evaluates the necessary derivatives ($U_0, U_2, U_3, \dots$) numerically at this peak.
• Higher-Order Extension: The code extends the WKB expansion from the previous 13th order up to the 16th order, incorporating recent theoretical developments.
• Padé Resummation: The package retains the use of balanced Padé approximants to sum the divergent WKB series, ensuring stability and accuracy, particularly for the dominant modes where $|Re(\omega)| \geq |Im(\omega)|$.
• Analytic Expansion Beyond Eikonal: The package includes routines to derive analytic expressions for QNMs and GBFs as series expansions in inverse powers of the multipole number $\kappa = \ell + 1/2$. This allows for systematic corrections beyond the leading eikonal limit ($\ell \to \infty$).
• Correspondence Implementation: It implements the analytic correspondence between QNM frequencies and S-matrix coefficients ($\mathcal{K}$), allowing GBFs to be estimated directly from QNM frequencies.

3. Key Contributions

• Algorithmic Optimization: The primary contribution is the replacement of symbolic differentiation with numerical evaluation of derivatives at the potential peak. This decouples the computation time from the algebraic complexity of the potential function.
• Performance Acceleration: The new code achieves a dramatic reduction in computation time. For complex potentials, the speed gain reaches several orders of magnitude. Calculations that previously took hours (or were impossible) are now completed in fractions of a second.
• Precision Improvement: By avoiding the generation of massive symbolic expressions, the new method significantly reduces precision loss associated with floating-point operations.
• Extended Order: The implementation of the 16th-order WKB formula allows for higher precision in spectral calculations.
• Unified Framework: The package provides a unified environment for calculating QNMs, GBFs, and their analytic correspondence, including tools for rotating black holes and various scalar/vector/tensor perturbations.

4. Results

The authors benchmarked the new package against the previous version using two test cases:

1. Schwarzschild-de Sitter (Dirac field): A standard test case with a relatively simple potential.
2. Regular Black Hole (Scalar field): A complex case arising from proper-time flow in quantum gravity, featuring a cumbersome analytic structure.

Key Findings from Benchmarks (Table 1):

• Time Efficiency: For the 16th-order calculation of the dominant mode in the Regular Black Hole case:
  • Old Code: ~17 seconds (and significantly longer for higher orders).
  • New Code: ~0.03 seconds.
  • The speedup is even more pronounced for higher orders (e.g., at order 13, the old code took ~32 seconds vs. ~0.28 seconds for the new code).
• Precision: The new code exhibits reduced precision loss (fewer binary places lost) for higher orders compared to the old symbolic approach.
• Scalability: In the new package, computation time is nearly independent of the complexity of the effective potential, whereas the old package's time increased drastically with potential complexity.

5. Significance

• Enabling Complex Models: This tool makes it feasible to study QNMs and GBFs for a wide class of "exotic" black hole solutions (regular black holes, wormholes, modified gravity models) that were previously too computationally expensive to analyze at high WKB orders.
• Gravitational Wave Phenomenology: By providing rapid and accurate calculations of ringdown frequencies and transmission coefficients, the code supports the interpretation of gravitational wave signals from black hole mergers and the testing of General Relativity against alternative theories.
• Theoretical Insight: The ability to perform high-order expansions and analytic series in $\ell$ allows for a deeper understanding of the convergence properties of the WKB method and the relationship between the eikonal limit and lower multipole numbers.
• Open Source Utility: The authors provide an ancillary Mathematica notebook with examples, facilitating immediate adoption by the research community for black hole spectroscopy and scattering problems.

In summary, this paper delivers a critical software update that transforms high-order WKB calculations from a computationally prohibitive task into a routine, efficient procedure, thereby expanding the scope of black hole perturbation theory research.