GRTresna: An open-source code to solve the initial data constraints in numerical relativity

GRTresna is an open-source, multigrid solver designed to generate initial data for numerical relativity simulations involving fundamental fields around black holes and inhomogeneous cosmological spacetimes, based on the formalism presented in arXiv:2207.03125.

The Gist

Imagine you are trying to build a massive, realistic simulation of the universe on a computer. Maybe you want to watch two black holes collide, or see how the very first moments of the Big Bang played out. To do this, you need to write a "script" for the universe. But before the movie can start rolling, you have to set the scene perfectly.

This paper introduces a new tool called GRTresna (which sounds like "tool" in the Basque language, fitting for a Swiss-army-knife of physics code). Think of GRTresna as a specialized "Scene Setter" for these universe simulations.

Here is how it works, broken down into simple concepts:

1. The Problem: The "Perfect Start" is Hard to Find

In the real world, gravity is described by Einstein's equations. To simulate this on a computer, scientists have to slice time into tiny steps. Before they can take the first step (simulating the next second), they must provide the computer with a snapshot of the universe at "Time Zero."

However, this snapshot isn't just a random picture. It has to obey strict rules, like a puzzle where every piece must fit perfectly. If the pieces don't fit, the simulation crashes or produces nonsense. These rules are called constraints.

• The Analogy: Imagine trying to build a house of cards. You can't just throw cards down randomly; the bottom layer must be perfectly balanced, or the whole thing collapses before you even start. GRTresna is the tool that helps you build that perfect, stable bottom layer.

2. What Makes GRTresna Special?

There are other tools out there that build these "Time Zero" snapshots, but they are often like specialized tools: great for building houses of cards with just two specific types of cards (like black holes), but bad at handling other materials.

GRTresna is different because it is flexible:

• It handles "Weird" Stuff: Most tools struggle when you add "fundamental fields" (like invisible energy fields or scalar fields) into the mix. GRTresna is designed specifically to handle these complex ingredients alongside black holes.
• It works in the "Big Picture": While other tools focus on small, dense objects (like merging stars), GRTresna can also set up the stage for the entire universe, including uneven, lumpy cosmic landscapes (inhomogeneous cosmology).
• It's a "Lego" Kit: The code is built so that scientists can easily swap out the "cards" they are using. If a researcher wants to test a new type of matter or a new theory of gravity, they can plug it into GRTresna without rebuilding the whole engine.

3. How It Works (The Magic Trick)

The paper explains that GRTresna uses a method called Multigrid.

• The Analogy: Imagine you are trying to smooth out a crumpled piece of paper. You could try to fix every tiny wrinkle one by one (which takes forever). Or, you could look at the big folds, smooth those out, then look at the medium folds, and finally fix the tiny ones.
• GRTresna does this mathematically. It solves the big problems first, then the medium ones, then the small ones. This makes it incredibly fast and efficient, even when dealing with complex, jagged shapes in the simulation.

4. What Has It Already Done?

The paper lists several "test drives" where GRTresna has successfully set the scene for other simulations:

• Dark Matter & Black Holes: It helped simulate how dark matter behaves around pairs of colliding black holes.
• The Early Universe: It helped study how the universe expanded and heated up right after the Big Bang (a phase called "preheating").
• Spinning Black Holes: It helped create the initial conditions for black holes that are spinning.
• New Gravity Theories: It helped test theories of gravity that are slightly different from Einstein's.

5. Who Is It For?

GRTresna is open-source, meaning anyone can download it, look at the code, and use it for free. It is designed to work seamlessly with a family of other physics codes (like GRChombo), acting as a reliable partner that prepares the data so the main simulation can run smoothly.

In summary: GRTresna is a versatile, fast, and open-source "Scene Setter" that allows physicists to build the perfect starting point for their universe simulations, especially when those simulations involve complex energy fields or the vast, uneven structure of the cosmos. It doesn't run the whole movie, but it ensures the first frame is perfect so the rest of the story can play out correctly.

Technical Summary

Technical Summary of GRTresna

Problem Statement
Numerical relativity (NR) is essential for solving the Einstein Equations in strong-field regimes, primarily focusing on the time evolution of spacetime metrics. However, before evolution can begin, the initial data set must satisfy four coupled, non-linear elliptic partial differential equations (PDEs) known as the Hamiltonian and momentum constraints. While these constraints can be solved under restrictive assumptions, such limitations significantly narrow the range of physical scenarios accessible to NR, particularly those involving fundamental fields (e.g., scalar fields) and inhomogeneous cosmological spacetimes. Existing solvers are often specialized for compact object mergers (neutron stars and black holes) or rely on spectral methods that, while highly accurate, may lack the flexibility required for complex, non-standard matter configurations or cosmological applications.

Methodology
The paper introduces GRTresna, an open-source multigrid solver designed to address these constraint equations. The code is built upon the formalism established in Aurrekoetxea, Clough & Lim [1] and is implemented in C++ using object-oriented programming and templating. Key methodological components include:

• Solving Framework: GRTresna implements two variations of the Conformal Transverse-Traceless (CTT) method: CTTK and CTTK-Hybrid. These methods are specifically noted for their suitability in handling fundamental fields within the matter content.
• Numerical Approach: The solver utilizes a multigrid method to efficiently reduce errors across a hierarchy of discretizations. This approach is integrated with the Chombo library, inheriting its Adaptive Mesh Refinement (AMR) capabilities, which use block-structured Berger-Rigoutsos grid generation.
• Parallelism: The code employs hybrid OpenMP/MPI parallelism to handle computational demands.
• Flexibility in Matter Sources: While the current version supports scalar fields, the templated architecture allows users to easily substitute other matter types (e.g., vector fields) by modifying the methods that compute energy and momentum densities.
• Boundary and Data Handling: The code supports extrapolating, reflective, and periodic boundary conditions, compatible with the NR evolution code GRChombo. It can generate analytical initial data or read grids from existing HDF5 files, facilitating the upgrade of matter configurations evolved on fixed backgrounds (e.g., via GRDzhadzha) to full NR simulations with backreaction.

Key Contributions

• General-Purpose Solver: GRTresna provides a flexible, open-source tool capable of handling both black hole spacetimes (including superpositions of boosted and/or spinning black holes via Bowen-York data) and inhomogeneous cosmological spacetimes with periodic boundary conditions.
• Integration with the GRTL Ecosystem: The code is designed to be fully compatible with the GRTL Collaboration's family of codes, particularly GRChombo, allowing for seamless data exchange via checkpoint files for $t=0$ restarts.
• Extensibility: The codebase is structured to allow the implementation of new solver methods, additional matter types, and extensions to theories beyond General Relativity (GR).
• Diagnostics: The solver computes Hamiltonian and momentum constraint errors at each iteration, outputting the norm of these values to text files for validation.

Results and Applications
The paper does not present new simulation results but rather demonstrates the code's utility through a summary of research projects where GRTresna has already been successfully employed. These applications include:

• Investigating the robustness of inflation against scalar field inhomogeneities.
• Studying the formation of oscillons during inflationary preheating.
• Analyzing the formation of spinning primordial black holes.
• Examining the effects of scalar dark matter environments around binary black holes.
• Modeling the general relativistic evolution of polarized Proca stars.
• Solving initial condition problems for modified gravity theories.

Significance
The paper positions GRTresna as a necessary complement to existing high-precision spectral solvers (such as TwoPunctures or SpECTRE), which are often optimized for waveform generation in compact object mergers. While GRTresna is limited to second-order accuracy due to its multigrid approach, its primary significance lies in its flexibility and generality. It addresses a gap in the field by providing a publicly available solver capable of handling inhomogeneous cosmological spacetimes and fundamental field matter types, scenarios for which fully general solvers were previously lacking. The code aims to expand the physics that NR evolutions can probe, particularly in scenarios involving theories beyond GR and the Standard Model, while remaining accessible for integration into other NR codes via standard HDF5 output.