$\mathcal{H}$-EFTCAMB: A Cobaya-Integrated, Python-Wrapped Extension of EFTCAMB for Covariant Horndeski Gravity

The paper introduces $\mathcal{H}$-EFTCAMB, a publicly available, Cobaya-integrated Python extension of EFTCAMB that enables robust cosmological testing of arbitrary covariant Horndeski gravity theories through both EFT mapping and direct scalar field equation solving, the latter of which remains valid even in regimes where the EFT approach fails.

The Gist

Imagine the universe is a giant, complex video game. For decades, scientists have been playing this game using a specific set of rules called General Relativity (Einstein's theory of gravity). In this game, the universe is expanding faster and faster, and the standard explanation is that there's a mysterious "Dark Energy" fueling this speed.

But what if the rules of the game are actually different? What if gravity itself changes over time or space, rather than just being fueled by a mysterious energy? This is the world of Modified Gravity.

This paper introduces a new, super-charged tool called H-EFTCAMB. Think of it as a massive software upgrade for the "Cosmic Simulator" that scientists use to test these new theories against real-world data.

Here is a breakdown of what this paper is about, using simple analogies:

1. The Old Tool vs. The New Tool

• The Old Tool (EFTCAMB): Imagine a very smart, reliable calculator that can solve equations for a specific list of gravity theories. It's great, but it's like a calculator that only knows how to do math for "Standard Gravity" and a few specific "Modified Gravity" flavors. If you wanted to test a brand-new, weird theory, you'd have to rewrite the calculator's code from scratch.
• The New Tool (H-EFTCAMB): This is the "Pro Version." The "H" stands for Horndeski, which is the "Master Key" or the "Ultimate Library" of all possible gravity theories that keep the math stable (no explosions or nonsense).
  • The Analogy: If the old tool was a locked box with a few specific keys, the new tool is a universal keyring. It can now open any door in the "Horndeski Library" without needing to rebuild the lock. You can feed it a raw recipe (a mathematical formula), and it will instantly cook up the predictions for how the universe should look.

2. Two Ways to Cook the Meal

The paper highlights two ways this new tool handles these theories:

• Method A (The Map): It translates the complex, raw recipe (the "Covariant Lagrangian") into a simplified map (the "EFT framework") that the old tool already knew how to read. This is like translating a foreign language into English so you can read it.
• Method B (The Direct Route): Sometimes, the "map" breaks down. Imagine a road that loops back on itself or a field that oscillates (wiggles back and forth). The old tool would crash here. The new tool has a Direct Route mode. It solves the equations directly, even when the math gets messy or the field "turns over." It's like having a drone that can fly through a storm where a car (the old tool) would get stuck.

3. The "Plug-and-Play" Feature

One of the biggest upgrades is how easy it is to use.

• The Analogy: In the past, connecting this scientific calculator to the data analysis software (called Cobaya) was like trying to wire a vintage radio to a modern smart home system. It required a lot of custom wiring and technical know-how.
• The Upgrade: H-EFTCAMB now comes with a Python Wrapper. Think of this as a universal adapter plug. You can now plug this powerful gravity simulator directly into the modern data analysis systems scientists use every day. It's "seamless," meaning you don't need to be a coding wizard to run complex simulations.

4. Why Do We Need This?

Scientists are currently facing a crisis in cosmology. Different ways of measuring the universe's expansion rate give slightly different answers (the "Hubble Tension").

• The Detective Work: To solve this, they need to test thousands of "what-if" scenarios. "What if gravity gets stronger in the future?" "What if it oscillates?"
• The Validation: The authors didn't just build the tool; they tested it rigorously. They compared their new "Universal Keyring" against the old calculator and against other top-tier simulators (like hi class). They found that H-EFTCAMB gives the exact same answers as the trusted tools but can also handle the weird, oscillating theories that the others can't touch.

5. A Real-World Example: "Thawing Gravity"

The paper ends with a demo. They used their new tool to test a theory called "Thawing Gravity with Screening."

• The Scenario: Imagine a theory where gravity is modified, but there's a "shield" (screening) that hides this modification in our solar system (so we don't notice it here) but lets it show up in the vast emptiness of the universe.
• The Result: Using H-EFTCAMB, they proved this theory is still a viable candidate to explain our data. They showed that the new "Direct Route" mode was essential because the field in this theory was wiggling (oscillating), which would have caused the old tools to crash.

Summary

H-EFTCAMB is the new, Swiss Army Knife for cosmologists.

1. It's Universal: It handles the most general class of gravity theories (Horndeski).
2. It's Robust: It can solve equations even when the math gets weird or oscillating.
3. It's User-Friendly: It plugs easily into modern data analysis software.
4. It's Verified: It works perfectly and is ready for the next generation of telescopes (like Euclid and DESI) to help us figure out if Einstein's gravity needs a software update.

In short, it's the tool that allows scientists to stop guessing and start rigorously testing the most exotic ideas about how our universe works.

Technical Summary

Here is a detailed technical summary of the paper "H−EFTCAMB: A Cobaya-Integrated, Python-Wrapped Extension of EFTCAMB for Covariant Horndeski Gravity."

1. Problem Statement

The field of observational cosmology has entered an era of precision, driven by data from the Cosmic Microwave Background (CMB), Large Scale Structure (LSS), and distance ladders. While the standard $\Lambda$CDM model explains observations well, tensions in data (e.g., the Hubble tension) suggest potential deviations from General Relativity (GR) or the need for dynamical Dark Energy (DE).

Existing numerical tools like EFTCAMB and hi class have been instrumental in testing modified gravity (MG) theories. However, they face specific limitations:

• EFTCAMB (Original): Primarily designed around the Effective Field Theory (EFT) of Dark Energy framework. It maps covariant theories to EFT operators. This approach breaks down when the scalar field undergoes "turn-overs" (e.g., oscillations where $\dot{\phi} \to 0$), as the unitary gauge becomes singular. It also lacked a general module to handle arbitrary covariant Lagrangians directly without pre-defined mappings.
• Integration: Older versions relied heavily on the native CAMB interface and CosmoMC, limiting seamless integration with modern Bayesian inference pipelines like Cobaya.
• Flexibility: There was a need for non-parametric reconstruction methods and the ability to define arbitrary functional forms for Lagrangian components directly from covariant theories, specifically within the Horndeski class (the most general scalar-tensor theory with second-order equations of motion).

2. Methodology

The authors present H−EFTCAMB, the official successor to EFTCAMB, which introduces a hybrid architecture combining the stability of the EFT framework with the directness of covariant Lagrangian evolution.

A. The Horndeski Module (Direct Covariant Evolution)

The core innovation is a new module (Flag 5) that accepts an arbitrary covariant Horndeski Lagrangian as input:
$$S_H = \int d^4x \sqrt{-g} (L_2 + L_3 + L_4 + L_5)$$
where $L_i$ are functions of the scalar field $\phi$ and its kinetic term $X$.

• Dual Solving Modes:
  1. EFT Mapping Mode: Maps the covariant Lagrangian to EFT operators and solves using the standard EFT framework (efficient for non-oscillating fields).
  2. Direct $\delta\phi$ Mode: Directly integrates the equations of motion for the scalar field perturbation $\delta\phi$. This bypasses the unitary gauge singularity, allowing the code to handle oscillating fields and turn-overs where $\dot{\phi}$ crosses zero, a scenario where the original EFT approach fails.
• Implementation: The Lagrangian functions ($G_i(\phi, X)$) are defined in a single subroutine (EFTCAMBHdskLagrangianCoefficients), allowing users to plug in arbitrary analytic functions, basis expansions, or spline interpolations.

B. Flexible Function Specification

H−EFTCAMB introduces advanced methods for defining free functions:

• Non-parametric Reconstruction: Users can define functions via cubic or quintic spline interpolations on a dense grid. This allows for data-driven reconstruction of EFT functions or Lagrangian components without assuming a specific parametric form.
• Multi-variable Dependence: EFT functions can now be defined as combinations of the scale factor $a$ and the dark energy density fraction $\Omega_{DE}$ (e.g., $\alpha_i = c_i \Omega_{DE}$), a feature previously unsupported in the standard EFT framework.

C. Python Wrapper and Cobaya Integration

• A robust Python wrapper exposes the Fortran backend, enabling seamless integration with Cobaya, the standard cosmological sampler.
• This allows for efficient Bayesian parameter estimation and nested sampling using modern datasets (Planck, DESI, Pantheon+).
• The wrapper supports outputting intermediate variables (background, perturbations, source functions) for debugging and theoretical analysis.

D. Validation Strategy

The code was validated through:

1. Internal Consistency: Comparing the new Horndeski module against the original EFTCAMB implementation of the Scaling Cubic Galileon (SCG) model.
2. External Cross-Checks: Comparing results with hi class for:
   • Jordan-Brans-Dicke (JBD) theory.
   • The $\alpha_i = c_i \Omega_{DE}$ parametrization.

3. Key Contributions

1. Direct Horndeski Solver: The first public Einstein-Boltzmann solver capable of directly evolving general Horndeski theories from their covariant Lagrangians, specifically handling oscillating scalar fields via the $\delta\phi$ formalism.
2. Non-Parametric Capabilities: Introduction of spline interpolation for arbitrary function definition, enabling model-independent reconstruction of DE/MG dynamics.
3. Modern Workflow Integration: Full integration with Cobaya via a Python wrapper, facilitating large-scale data analysis with next-generation surveys.
4. Expanded Model Library: Inclusion of built-in mappings for various DE models (e.g., $f(R)$, K-mouflage, Acoustic DE) and new initial conditions consistent with non-zero EFT functions in the early universe, alongside positivity bounds checks.
5. Unified Flag Structure: A reorganized flag tree (Fig. 1 & 2) that clarifies model selection, distinguishing between "Designer" (inputting background) and "Solving Background" (deriving from Lagrangian) approaches.

4. Results

• Validation:
  • SCG Comparison: The new Horndeski module produced CMB temperature/polarization spectra and matter power spectra that agreed perfectly (sub-percent level) with the original EFTCAMB SCG implementation.
  • JBD Comparison: Results for Jordan-Brans-Dicke theory matched hi class across a range of coupling parameters ($\omega_{BD}$), confirming the accuracy of the direct $\delta\phi$ solver.
  • $\alpha_i$ Parametrization: Successfully reproduced results for the $\alpha_i = c_i \Omega_{DE}$ parametrization, which was previously impossible in EFTCAMB due to its restriction to scale-factor-dependent functions.
• Application Example (Thawing Gravity with Screening - TGs):
  • The authors applied H−EFTCAMB to a "Thawing Gravity" model augmented with a Galileon screening operator ($X\Box\phi$).
  • Using nested sampling with Planck, DESI, and Pantheon+ data, they demonstrated that the screening mechanism allows the model to satisfy solar system constraints while maintaining the ability to explain phantom crossing ($w_{DE} < -1$).
  • Crucially, the model involved oscillating fields near recombination, a regime where hi class fails, demonstrating the unique capability of H−EFTCAMB.

5. Significance

H−EFTCAMB represents a major leap forward in the computational toolkit for testing gravity. By bridging the gap between high-energy theoretical Lagrangians and observational data, it allows cosmologists to:

• Test general Horndeski theories without relying on specific parametrizations or approximations that break down in oscillating regimes.
• Perform model-independent reconstructions of cosmic acceleration using non-parametric methods.
• Leverage the full power of Cobaya for Bayesian inference with upcoming precision data (e.g., Euclid, DESI, CMB-S4).

The code is publicly available and serves as the new standard for confronting modified gravity and dark energy models with current and future cosmological surveys.