---
title: 'GRTeclyn: AMReX-Enhanced Numerical Relativity'
url: https://www.emergentmind.com/topics/grteclyn
type: topic
---

# GRTeclyn: AMReX-Enhanced Numerical Relativity

GRTeclyn is a numerical-relativity codebase for fully relativistic simulations of dynamical spacetimes. The published literature presents it most concretely as a **“Port of GRChombo to AMReX”** under active development, and uses it for three-dimensional evolutions of the unstable Ellis–Bronnikov wormhole with gravitational-wave extraction [2604.00071]. A closely related cosmology study describes a tensor-perturbation pipeline implemented in **GRChombo** and states that the implementation is being actively ported to **GRTeclyn**, with the expectation that the GRTeclyn version will follow the GRChombo implementation closely [2412.19731]. Taken together, these works place GRTeclyn within a GRChombo-derived, AMReX-based ecosystem for strong-field general relativity, adaptive mesh refinement, and relativistic perturbation analysis.

## 1. Identity, lineage, and scope

The clearest direct characterization of GRTeclyn appears in the wormhole-collapse study, where it is described as a **codebase** and specifically as a **“Port of GRChombo to AMReX”** [2604.00071]. In that work, GRTeclyn is used to initialize exact isotropic wormhole data on a three-dimensional Cartesian grid, evolve the coupled Einstein–phantom-scalar system, apply **CCZ4** with matter source terms, use **adaptive mesh refinement (AMR)**, track collapse diagnostics, extract gravitational radiation through the Weyl scalar \(\Psi_4\), and monitor constraints.

The inflationary tensor paper uses more qualified language. Its concrete numerical results are obtained with **GRChombo**, but the authors repeatedly state that they are actively porting a dedicated inflationary example to **GRTeclyn** and that the implementation in GRTeclyn is expected to follow closely [2412.19731]. This establishes GRTeclyn not merely as a separate code name, but as the intended recipient of a validated computational pipeline for cosmological tensor initialisation, evolution, and extraction.

A plausible implication is that GRTeclyn is best understood not as a single-purpose application, but as a general infrastructure layer for numerical relativity. The available papers associate it with two distinct regimes: strong-field compact-object dynamics and early-universe tensor perturbations.

## 2. Evolution formalisms and computational infrastructure

The literature associates GRTeclyn with two evolution frameworks. In the wormhole paper, the coupled geometry–matter system is evolved with **CCZ4 with matter**, and the scalar variables are advanced through the **CCZ4RHSWithMatter** module in GRTeclyn [2604.00071]. The evolved state includes conformal metric variables, \(\chi\), \(K\), shift, lapse, and scalar variables \((\phi,\Pi)\). In the inflationary tensor pipeline, the evolution is formulated in **BSSN**, with evolved variables
\[
\{\chi,\tilde{\gamma}_{ij},K,\tilde{A}_{ij},\alpha,\beta^i\},
\]
and matter supplied by the inflaton sector [2412.19731].

The wormhole study gives the more explicit implementation profile. It states that GRTeclyn uses an **AMReX backend** for **block-structured AMR** and **GPU acceleration**, and that production runs use **NVIDIA H100 GPUs** with **one MPI rank per GPU** [2604.00071]. The reported production setup is
\[
L_{\rm full}=64,\qquad N_{\rm full}=256,\qquad dx_{\rm coarse}=0.25,
\]
with up to **5 AMR levels** at **2:1 refinement**, regridded on **\(\chi\)-gradients**, giving
\[
dx_{\rm fine}\approx 7.8\times 10^{-3}.
\]
Time integration uses **4th-order Runge–Kutta** with `dt_multiplier = 0.02`; **octant symmetry** is imposed; inner boundaries use **parity conditions**; outer boundaries use **Sommerfeld boundary conditions**; and the wave-extraction radii satisfy
\[
R_{\rm ext}\in[12,24].
\]

Several numerical controls are also specified in that work: **Kreiss–Oliger dissipation** with
\[
\sigma_{\rm KO}=2.0,
\]
a floor on the conformal variable,
\[
\chi_{\rm min}=10^{-8},
\]
and CCZ4 damping parameters
\[
\kappa_1=3.0,\qquad \kappa_2=0,\qquad \text{covariantZ4}=0.
\]
Constraint monitoring is carried out with volume-weighted \(L_2\) norms on the unrefined level-0 grid [2604.00071].

By contrast, the inflationary tensor study emphasizes the gauge and perturbative-consistency aspects of the GRChombo/GRTeclyn workflow. For the tensor-only runs it uses **geodesic gauge / synchronous gauge** in cosmic time, with initial lapse \(\alpha=1\), initial shift \(\beta^i=0\), and all Bona–Massó coefficients set to zero so that the gauge variables remain fixed [2412.19731]. The same paper also states that the time evolution in GRChombo is **second-order accurate**, that the simulations are performed in a finite **periodic** box implicit in the Fourier-lattice construction, and that **KO dissipation** and **CCZ4-type damping** are disabled for the validation runs because they can distort the physical tensor spectrum.

## 3. Initial data and matter content

One of the main published demonstrations of GRTeclyn is the evolution of an Ellis–Bronnikov wormhole supported by a **massless phantom scalar field** [2604.00071]. The matter model is defined by
\[
T_{\mu\nu} = -\left( \nabla_\mu\phi\,\nabla_\nu\phi -\frac12 g_{\mu\nu}\nabla^\alpha\phi\,\nabla_\alpha\phi -g_{\mu\nu}V(\phi) \right),
\]
with
\[
V(\phi)=0.
\]
The wormhole is initialized in isotropic coordinates so that the conformal variable can be written in algebraically regularized form,
\[
\chi=\left(\frac{4\bar r^2}{4\bar r^2+b_0^2}\right)^2,
\]
while the static initial data satisfy
\[
K_{ij}(t=0)=0,\qquad \Pi(t=0)=0,\qquad \alpha(t=0)=1,\qquad \beta^i(t=0)=0.
\]
The paper stresses that evolving \(\chi\) rather than \(\psi\) makes the compactified origin regular in moving-puncture-type variables.

The same study also details how collapse is forced. To select the compressive branch of the instability, the phantom stress-energy support is globally rescaled,
\[
T_{\mu\nu}\to S_{\rm support}\,T_{\mu\nu},
\]
with
\[
S_{\rm support}=0.5,
\]
and a quadrupolar scalar perturbation is added,
\[
\phi(\bar r,\theta,\varphi)\big|_{t=0} = \phi_{\rm EB}(\bar r) + A_\phi\,Y_{20}(\theta,\varphi)\,\exp\!\left(-\frac{\bar r^2}{\sigma_\phi^2}\right),
\]
using
\[
A_\phi=0.02,\qquad \sigma_\phi=0.5.
\]
The perturbation is applied to \(\phi\) rather than to \(K_{ij}\) or \(\Pi\), so that with \(\Pi(t=0)=0\) and \(K_{ij}(t=0)=0\), the momentum constraint remains exactly satisfied at the initial time [2604.00071].

The cosmological work describes a different initial-data pipeline that is intended for GRChombo/GRTeclyn. It constructs stochastic vacuum tensor perturbations from the Mukhanov–Sasaki mode functions, using a polarization decomposition
\[
h_{ij} = h_+ \epsilon^+_{ij} + h_{\times} \epsilon^{\times}_{ij},
\]
Gaussian random amplitudes in Fourier space, and a crucial **\(k\)-dependent mode-function phase** rather than a random phase alone [2412.19731]. The paper is explicit that the improved method is
\[
h_s(k) = \sqrt{P_{h_s}(k)}e^{i(\theta_{MS}(k)+\theta_s)},
\]
and that this phase is required to reproduce the smooth temporal power evolution predicted by linear theory.

That tensor data are then mapped into the BSSN variables through the tensor-only CPT–BSSN dictionary,
\[
\chi = \frac{1}{a^2},\qquad \tilde{\gamma}_{ij} = \delta_{ij}+h_{ij},\qquad K=-3H,\qquad \tilde A_{ij} = -\frac{1}{2\bar N}\dot h_{ij}.
\]
The authors state that they have developed a new initial-condition class in GRChombo which generates this random tensor perturbation on the lattice, and that this functionality is expected to be made public in a future release of **GRChombo/GRTeclyn** [2412.19731].

## 4. Diagnostics, extraction, and verification

The wormhole paper presents GRTeclyn as a platform for collapse diagnostics and gravitational-wave extraction [2604.00071]. The areal-radius diagnostic is written as
\[
R_{\rm areal}(\bar r)=\frac{\bar r}{\sqrt{\chi(\bar r)}},
\]
and the sign of \(K\) is interpreted through
\[
\partial_t(\sqrt{\gamma}) = -\alpha K \sqrt{\gamma} + \nabla_i\beta^i\sqrt{\gamma},
\]
so that \(K>0\) indicates contraction and \(K<0\) indicates expansion. Because there is **no production elliptic apparent-horizon finder** in GRTeclyn in that work, the authors implement a trapped-surface proxy by evaluating the outgoing null expansion \(\theta_+\) on coordinate spheres and using the criterion
\[
\theta_+\le 0.
\]

Gravitational radiation is extracted through the Newman–Penrose scalar
\[
\Psi_4=\ddot h_+ - i\ddot h_\times.
\]
The extracted signal is decomposed into spin-weighted spherical harmonics, and the dominant mode is found to be \((\ell,m)=(2,0)\), consistent with the imposed \(Y_{20}\) perturbation. A key methodological test is the propagation-speed estimate
\[
v=\frac{R_{i+1}-R_i}{t_{i+1}^{\rm peak}-t_i^{\rm peak}},
\]
used to distinguish physical gravitational waves from superluminal CCZ4 constraint modes. For the perturbed collapse run, the paper reports
\[
v\approx 0.995c
\]
between \(R=12\) and \(R=16\), with mild numerical dispersion further out, while the unperturbed run yields only very weak, superluminal signals interpreted as constraint-related contamination [2604.00071].

The same work also defines a one-sided power spectral density for \(\Psi_4\),
\[
S_{\Psi_4}(f) = \frac{\Delta t^2}{T} \left( \left|\mathcal{F}\{w(t)\,\mathrm{Re}(\Psi_4)\}\right|^2 + \left|\mathcal{F}\{w(t)\,\mathrm{Im}(\Psi_4)\}\right|^2 \right),
\]
uses a Morlet-wavelet continuous wavelet transform for time–frequency analysis, and converts to strain through
\[
S_h(f)=\frac{S_{\Psi_4}(f)}{(2\pi f)^4},\qquad h_{\rm char}=\sqrt{f\,S_h}.
\]
These quantities are then compared to detector sensitivity curves.

The inflationary tensor study develops a different extraction and verification pipeline for the GRChombo/GRTeclyn framework [2412.19731]. In the perturbative tensor-only regime it extracts the tensor perturbation as
\[
h_{ij}(t)=\tilde{\gamma}_{ij}(t)-\delta_{ij}.
\]
For future scalar-plus-tensor runs it also provides a more general transverse-traceless projector,
\[
h_{ij} = \tilde{\gamma}_{ij} - \delta_{ij} - \frac{1}{3} \delta_{ij} g - 2 D^{(i} V^{j)}.
\]
The primary validation observable is the tensor power spectrum, obtained from the polarization fields \(h_+(k)\) and \(h_\times(k)\) and compared against the Mukhanov–Sasaki prediction. Constraint behavior is analyzed through the Hamiltonian and momentum constraints, and the paper states that the tensor-only stochastic initial conditions violate constraints only at **second order or smaller**, with the corresponding numerical checks showing that these conditions are met [2412.19731].

## 5. Demonstrated scientific applications

The most direct GRTeclyn application in the current literature is three-dimensional wormhole dynamics [2604.00071]. With exact initial data, full support
\[
S_{\rm support}=1,\qquad A_\phi=0,
\]
and flat initial lapse, the code keeps the wormhole nearly static up to \(t\lesssim 1.5M\), after which truncation-level effects push the solution onto the **rarefactive expansion branch**. The throat then expands exponentially with measured rate
\[
\lambda\approx 9.012\,M^{-1}.
\]
The paper reports that moving-puncture gauge is much better suited to collapse than to this inflationary expansion, since the lapse “anti-collapses,” coordinate stretching outruns AMR, and constraints grow.

For the perturbed collapse run,
\[
S_{\rm support}=0.5,\qquad A_\phi=0.02,\qquad \sigma_\phi=0.5,
\]
the wormhole undergoes rapid compression. The areal radius drops from \(0.50\) to approximately
\[
R_{\rm areal}\approx 0.14
\]
by \(t\approx 2M\); a trapped surface forms almost immediately; and \(\Pi\), initially exactly zero, is dynamically generated up to \(\sim 0.4\). After horizon formation, the swallowed phantom matter produces a **“phantom bounce”**, the areal radius grows again to about
\[
R_{\rm areal}\approx 0.28,
\]
and an outward curvature shock propagates through the grid. By
\[
t\approx 18.5M,
\]
the trapped-surface proxy disappears and the lapse falls to a floor near \(10^{-10}\). This sequence is significant because it shows GRTeclyn evolving a topologically nontrivial, exotic-matter spacetime through collapse, apparent-horizon-proxy formation, rebound, and wave emission in full \(3\)D.

The same paper also gives an astrophysical scaling argument. For a wormhole with
\[
M\sim 10^3\,M_\odot
\]
at
\[
D=1\,{\rm Mpc},
\]
the dominant signal moves into the \(\sim 100\)–\(300\) Hz band, but for the moderate perturbation amplitude
\[
A_\phi=0.02
\]
the characteristic strain lies **slightly below** Advanced LIGO design sensitivity. The authors state that detectability would require closer sources, larger initial asymmetries, or next-generation detectors [2604.00071].

A second application area is early-universe tensor dynamics, although here the validated numerical results are explicitly presented for **GRChombo**, with GRTeclyn identified as the target of an ongoing port [2412.19731]. That work presents a complete method for the initialisation and extraction of first-order inflationary tensor perturbations with full gravitational backreaction, shows agreement of the background quantities \(a(t)\), \(H(t)\), and \(\epsilon\) with the Friedmann/Klein–Gordon system, and reproduces the expected transition of the tensor spectrum from subhorizon
\[
P_h\propto k^{-1}
\]
to superhorizon
\[
P_h\propto k^{-3}
\]
behavior. It also reports preliminary non-Gaussianity diagnostics based on skewness and kurtosis, with the genuinely nonlinear signal appearing as subtle amplitude-dependent growth for large rescalings \(A\gtrsim 100\).

## 6. Limitations and development trajectory

The current literature also defines the main limitations of GRTeclyn and the near-term development path. In the wormhole study, there is **no production elliptic apparent-horizon finder**, so horizon identification is approximate; **excision is not used**, with moving-puncture gauge relied upon to handle horizon formation; finite-radius extraction is used instead of larger-radius extraction or Cauchy-characteristic extraction; and late-time bounce-driven shock-like gradients degrade convergence and eventually overwhelm the AMR strategy [2604.00071]. The same work emphasizes that the perturbed collapse setup introduces a small Hamiltonian residual at \(t=0\), and explicitly suggests future use of an elliptic solver for fully constraint-satisfying large-amplitude data.

The cosmology pipeline carries a different set of caveats. Its published implementation is still in **GRChombo**, not yet in fully documented GRTeclyn form; the validation runs are **tensor-only**; the simple extraction formula
\[
h_{ij}=\tilde{\gamma}_{ij}-\delta_{ij}
\]
is appropriate only while trace and longitudinal contamination remain negligible; and the work does not yet provide a direct tensor bispectrum estimate [2412.19731]. The planned next steps are to combine scalar and tensor stochastic initialisation methods, evolve scalar and tensor perturbations together in full numerical relativity, study mixed bispectra such as \(\langle\zeta\zeta\gamma\rangle\), and connect the simulations to the **MODAL** bispectrum-estimation pipeline.

These limitations do not define GRTeclyn as incomplete in a generic sense; rather, they identify the specific frontier at which the codebase is being used. The published record shows a framework already capable of GPU-accelerated \(3\)D evolutions with AMR, exotic matter, and waveform extraction, while simultaneously serving as the intended destination for a perturbatively validated cosmological tensor pipeline [2604.00071; 2412.19731].

Source: https://www.emergentmind.com/topics/grteclyn