SXS Catalog: Extreme Spacetimes Simulations

• The SXS Catalog is a comprehensive, curated archive offering high-accuracy numerical and analytic solutions to Einstein’s equations for simulating extreme spacetimes.
• It leverages advanced spectral methods, adaptive mesh refinement, and systematic COM and memory corrections to generate robust gravitational-waveform data.
• The catalog spans extensive binary black hole, neutron star, and exotic spacetime simulations, underpinning parameter estimation and strong-field general relativity tests.

The Simulating eXtreme Spacetimes Catalog (often abbreviated as the SXS or eXtreme Spacetimes Catalog) is a comprehensive, systematically curated, publicly accessible archive of numerical and analytic solutions to Einstein’s equations in dynamical and exotic spacetimes. It provides high-accuracy gravitational waveforms, remnant properties, and geometric data from binary black hole (BBH), neutron star, and other extreme-orbit astrophysical systems, as well as computational frameworks and explicit algorithms for generating, parameterizing, and manipulating highly symmetric, singular, or quantum-spacetime scenarios. The SXS Catalog is widely used as the principal dataset underpinning waveform modeling, gravitational-wave (GW) parameter estimation, and strong-field tests of general relativity with current and next-generation detector networks.

1. Genesis, Scope, and Catalog Evolution

The initial SXS Catalog was released in 2013, comprising 174 distinct BBH simulations generated by the Spectral Einstein Code (SpEC) using the generalized harmonic formulation, with up to 35 inspiral orbits and mass ratios $q=m_1/m_2\in[1,8]$ and spin magnitudes up to $|\chi_{1,2}|=0.98$ (Mroue et al., 2013). An order-of-magnitude increase occurred with the 2019 update to 2,018 configurations, expanding mass ratios to $q\in[1,10]$ and including 1,426 precessing systems (Boyle et al., 2019). The latest catalog (version 3.0.0, 2025) nearly doubles the entries to 3,756, reaching $q=8$ for dense precessing coverage and up to $q\approx 20$ in nonprecessing cases with $|\vec\chi|\le0.8$ (Scheel et al., 19 May 2025). Eccentric, aligned, anti-aligned, and precessing configurations are represented, with more than 250 eccentric and approximately 1,000 fully precessing systems in the 2025 release.

Besides core BBH simulations, the catalog includes waveform corrections such as consistent center-of-mass (COM) alignment and nonlinear gravitational-wave memory, as well as explicit constructions of exotic and singular spacetimes (e.g., Taub-NUT, NHEK, swirling universes (Colléaux et al., 26 Sep 2025); Cauchy-compact Minkowski/BTZ manifolds (Brunswic, 2021)), high-efficiency code-based parameter explorations for extreme-mass-ratio inspirals (EMRIs) (East et al., 2013), and laboratory quantum-spacetime analog models (Chen et al., 2021).

2. Methodologies and Data Structures

Binary Black Hole and Compact Object Simulations

All waveforms are produced by evolving the full nonlinear Einstein field equations using SpEC’s multidomain spectral approach, adaptive mesh refinement (AMR), and advanced initial data techniques. The domain comprises overlapping curved hexahedra and spherical shells, dynamically refined according to truncation error. Inner boundaries enforce outflow at apparent horizons; outer boundaries use constraint-preserving radiative conditions (Mroue et al., 2013, Boyle et al., 2019, Scheel et al., 19 May 2025). Waveform extraction is performed on $\gtrsim 24$ coordinate spheres at $R\ge100M$ and extrapolated in powers of $1/R$, typically to order $N=2$ or $|\chi_{1,2}|=0.98$0, to obtain the asymptotic strain $|\chi_{1,2}|=0.98$1 and Newman–Penrose scalar $|\chi_{1,2}|=0.98$2 (Iozzo et al., 2020).

Key improvements across versions:

• COM Correction: Removal of residual coordinate drifts via Newtonian COM fitting and supertranslation, suppressing spurious mode mixing in precessing systems and yielding physically smooth amplitudes.
• Memory Correction: Addition of displacement memory via BMS balance-law-derived corrections, implemented in the sxs Python package and satisfying consistency with CCE and balance laws to $|\chi_{1,2}|=0.98$3 fractional error (Mitman et al., 2020).
• Spectral Efficiency: Through multidomain spectral methods, SXS simulations are $|\chi_{1,2}|=0.98$4 more efficient in CPU-hours per waveform cycle compared to finite-difference codes, enabling extended runs (up to 148 orbits in the latest release) at high accuracy (Scheel et al., 19 May 2025).
• Metadata and Organization: Each simulation includes metadata (mass ratio, spin vectors, eccentricity, waveform length, remnant properties), multiple HDF5 and ASCII data formats, and explicit file versioning for reproducibility.

Mathematical Frameworks for Exotic Spacetimes

Extensions to maximally symmetric manifolds (de Sitter, anti-de Sitter), near-horizon extreme Kerr (NHEK), swirling universes, and Cauchy-compact flat 3-manifolds with BTZ singularities are handled by parametric constructions, group-theoretic identifications, and explicit coordinate-analytic continuation maps (Kopczynski, 2023, Colléaux et al., 26 Sep 2025, Brunswic, 2021). The catalog documents explicit embedding formulas, isometry group generators, extension (“cusp-compactification”) procedures, and algorithms for constructing fundamental domains and singularity insertions (see Section 5).

Extreme Mass-Ratio and Quantum Simulation Extensions

For extreme mass-ratio (e.g., star–supermassive black hole orbits), methods decompose the metric as $|\chi_{1,2}|=0.98$5, with $|\chi_{1,2}|=0.98$6 the known background and $|\chi_{1,2}|=0.98$7 the dynamical perturbation, employing error subtraction at each timestep to suppress background truncation error (East et al., 2013). For quantum-spacetime analogs, Dirac equations in exotic metrics (Alcubierre, FLRW, etc.) are mapped to cold ion Hamiltonians, enabling simulation of superluminal propagation and relativity effects via laser and parametric trap controls (Chen et al., 2021).

3. Accuracy, Validation, and Error Characterization

Catalog waveforms undergo stringent convergence tests across resolutions and extraction radii:

• Mismatch: Formally, $|\chi_{1,2}|=0.98$8. Median mismatch between highest SpEC resolutions is $|\chi_{1,2}|=0.98$9; waveform-extraction uncertainty is $q\in[1,10]$0. Extrapolation inaccuracies between $q\in[1,10]$1 and $q\in[1,10]$2 are $q\in[1,10]$3 for most cases (Boyle et al., 2019, Scheel et al., 19 May 2025).
• Remnant Error: 90th percentile uncertainties on remnant mass and spin are $q\in[1,10]$4 and $q\in[1,10]$5, an order of magnitude below fit-based discrepancies (Boyle et al., 2019).
• Gauge and Balance Law Compliance: Use of the full set of Weyl scalars ($q\in[1,10]$6) and BMS law-based corrections ensures waveform consistency to O($q\in[1,10]$7) in the Bondi frame at null infinity (Mitman et al., 2020, Iozzo et al., 2020).

Validation against post-Newtonian and perturbative approximations demonstrates phase agreement within $q\in[1,10]$8–$q\in[1,10]$9 rad for inspiral precession, radiated GW energy matching point-particle predictions to few-percent accuracy, and full waveform compatibility for parameter estimation (Mroue et al., 2013, East et al., 2013).

4. Parameter-Space Coverage and Data Access

The SXS Catalog densely samples the following parameters (see Table):

Parameter	2013 Catalog	2019 Catalog	2025 Catalog
Mass ratio ($q=8$0)	1–8	1–10	1–8 (dense), $q=8$120 (some)
Spin $q=8$2	up to 0.98	up to 0.998	$q=8$3 (main grid)
Precessing runs	91	1,426	$q=8$41,000
Eccentric runs	43	—	$q=8$5250
Quasi-circular, aligned	40+	592+	$q=8$62,700
Orbits/waveform	up to 35	7–351 cycles	up to 148 orbits

All data, including time series for $q=8$7, multipolar $q=8$8, horizon diagnostics, and metadata, are accessible via the sxs Python package and at https://data.black-holes.org (Scheel et al., 19 May 2025). Versioned DOIs are provided for reproducibility.

5. Exotic and Singular Spacetimes Subcatalogs

The catalog includes algorithmic “recipes” and parameterization theorems for a wide class of extreme and exotic spacetimes:

• Double Wick Rotations: Complex coordinate continuations relate families such as Taub-NUT, NHEK, swirling, and Melvin spacetimes (with variants for $q=8$9 curvature), via closed-form metric transformations and Lie algebra mappings among Killing fields, independent of the underlying theory (Colléaux et al., 26 Sep 2025). The procedure is systematically codified, including analytic expressions for resulting metrics and their isometry generators.
• Cauchy-Compact Flat Spacetimes with Extreme BTZ: The global homeomorphism between moduli spaces of Cauchy-maximal locally Minkowski manifolds with extreme BTZ singularities and the tangent bundle of Teichmüller space is established; explicit extension (cusp-compactification) procedures are provided (Brunswic, 2021).
• Simulations in Maximally Symmetric and Embedded Spaces: Real-time, numerically stable frameworks for simulating observer kinematics, isometry actions, and causality in (anti-)de Sitter universes are implemented via tessellations, isometry composition, and mesh projection (Kopczynski, 2023).
• Quantum Simulation Analogs: Trapped-ion systems and Hamiltonian engineering enable exploration of Dirac propagation, light-cone tilting, and quantum relativistic effects in curved backgrounds (Chen et al., 2021).

6. Applications in Gravitational Wave Astronomy and Beyond

The SXS Catalog is the reference standard for the construction and validation of GW template families, effective-one-body (EOB), phenomenological and surrogate models, and hybrid (PN-NR) waveforms for LIGO, Virgo, and future detectors (Mroue et al., 2013, Boyle et al., 2019, Scheel et al., 19 May 2025). It underpins direct parameter estimation, detector response studies, systematic bias quantification, and strong-field GR tests. The inclusion of memory, consistent COM frames, and extended parameter coverage enables high-fidelity injections, tests of BMS law satisfaction, and the modeling of signals from BNS and mixed binary systems (Mitman et al., 2020, Shankar et al., 2022). The catalog facilitates research in quantum spacetime simulation, relativistic visualization, and analytic geometry construction, supporting theoretical explorations at the interface of GR, quantum field theory, and geometric topology.

7. Future Directions and Expansions

Ongoing catalog expansion targets systematic coverage of the full BBH parameter space to $q\approx 20$0, higher spin magnitude, diverse eccentricity, and more detailed remnant and matter sector diagnostics. Forthcoming work emphasizes:

• Automated inclusion of higher-order memory effects (spin, center-of-mass), full Weyl scalar datasets, and multiple gauge/CCE extrapolations (Mitman et al., 2020, Iozzo et al., 2020).
• Cross-catalog compatibility, e.g., workflow pipelines for RIT, Georgia Tech, ET AL, and waveform modelers.
• More extensive neutron star and magnetohydrodynamic (MHD) entries via codes such as GRaM-X (Shankar et al., 2022).
• Extension to cosmological and quantum-gravitational parameter spaces through analog experiments and explicit coordinate constructions (Chen et al., 2021, Kopczynski, 2023).

The SXS Catalog thus forms the backbone of computational and analytic research on extreme, dynamical, and singular spacetimes, driving methodological development and applications across gravitational-wave astrophysics, geometric topology, quantum analog simulation, and numerical relativity.