---
title: 'Fastrot-spec: Neutron Star Oscillation Code'
url: https://www.emergentmind.com/topics/fastrot-spec
type: topic
---

# Fastrot-spec: Neutron Star Oscillation Code

“fastrot-spec” is, as an *Editor’s term*, a concise designation for the ROXAS framework in its intended regime of **fast, rigidly rotating isolated neutron stars**: a **new pseudospectral, nonlinear, general-relativistic code** for computing stellar oscillations and extracting **axisymmetric and non-axisymmetric mode frequencies** relevant to the **kHz gravitational-wave band** expected from binary-neutron-star post-merger remnants [2410.16764]. The code is formulated for **cold barotropic equations of state**, uses **primitive-variable hydrodynamics** coupled to the **extended conformal-flatness condition (xCFC)**, and has the notable capability to **follow the stellar surface at every step**, which is central for non-spherical rotating equilibria [2410.16764].

## 1. Scientific setting and target regime

ROXAS, short for **Relativistic Oscillations of non-aXisymmetric neutron stArS**, is aimed at the oscillation problem for **isolated rotating neutron stars**, with explicit focus on the regime most relevant to **hypermassive neutron stars** formed after binary neutron star mergers [2410.16764]. The underlying motivation is observational: next-generation gravitational-wave detectors are expected to improve sensitivity in the **kHz band**, where such remnants are expected to emit, so robust mode-frequency predictions are needed for gravitational-wave interpretation and equation-of-state constraints [2410.16764].

Within its current implementation, the framework supports **fast, rigidly rotating isolated neutron stars**, **cold barotropic equations of state**, and both **axisymmetric** and **non-axisymmetric oscillations** [2410.16764]. The paper also emphasizes that the code is **very lightweight** and can run on **office computers** [2410.16764]. This combination of target physics and modest computational footprint places the method between perturbative oscillation calculations and full dynamical spacetime evolutions: it is constructed for precision frequency extraction in a constrained but astrophysically important sector of parameter space.

A common misunderstanding is to read the “fastrot” label as implying differential rotation or merger-scale microphysics. The present implementation does not do so. The supported rotation law is **rigid rotation**, and the matter model is **cold**, **barotropic**, and **\(\beta\)-equilibrated** [2410.16764]. Thermal effects and out-of-\(\beta\)-equilibrium physics are identified instead as future extensions.

## 2. Relativistic formulation and equilibrium-preserving evolution

The hydrodynamics is written in terms of **primitive variables** rather than conserved variables, following the formalism cited in the paper, with the perfect-fluid stress-energy tensor
$$
T^{\mu\nu} = (e+p)u^\mu u^\nu + pg^{\mu\nu}.
$$
The matter model is **zero-temperature, \(\beta\)-equilibrated baryonic matter** with a **barotropic EoS**. The evolved variables are the **log-enthalpy**,
$$
H = \ln\left(\frac{e+p}{m_B n_B}\right),
$$
and the **Eulerian velocity** \(U_i\), with sound speed
$$
c_s^2 = \frac{dp}{de}.
$$
The evolution equations for \(H\) and \(U_i\) are given explicitly in the paper as Eqs. \((\ref{eq:evolH})\) and \((\ref{eq:evolUi})\) [2410.16764].

For the spacetime sector, ROXAS adopts the **xCFC** reduction of Einstein’s equations in the fully constrained formalism. The gauge and geometric conditions are
$$
K = 0,\qquad \bar{D}_i\tilde{\gamma}^{ij} = 0,\qquad \tilde{\gamma}^{ij}=f^{ij},
$$
corresponding respectively to **maximal slicing**, **Dirac gauge**, and **conformal flatness** [2410.16764]. Under these assumptions, the metric problem reduces to a hierarchical system of **two linear vector** and **two nonlinear scalar Poisson-like PDEs** [2410.16764]. In practice, this is the key approximation that makes the method substantially cheaper than full GR evolution while remaining accurate enough for the frequency studies reported.

A central numerical feature is the **well-balanced formulation**. Every evolved quantity is decomposed into equilibrium plus perturbation, for example
$$
U_i = U_{i,\mathrm{eq}} + \bar{U}_i,\qquad
H = H_\mathrm{eq} + \bar{H},\qquad
\beta^i = \beta^i_\mathrm{eq} + \bar{\beta}^i,\qquad
N = N_\mathrm{eq} + \bar{N}.
$$
For stationary axisymmetric equilibrium, the first integral is written as
$$
D_i\left(H_\mathrm{eq} + \ln N_\mathrm{eq} - \ln \Gamma_\mathrm{eq}\right)=0.
$$
The rewritten perturbative evolution equations are designed to **remove cancellation errors**, **preserve equilibrium exactly**, and **exploit double precision more effectively** [2410.16764]. This is significant because oscillation calculations are particularly sensitive to spurious drift away from equilibrium; the formalism is therefore tuned not merely for stability, but for long-time fidelity of small-amplitude mode content.

## 3. Pseudospectral architecture and moving-surface hydrodynamics

ROXAS is built on **LORENE** and uses **Chebyshev pseudospectral methods** in spherical coordinates for both metric and hydrodynamic sectors [2410.16764]. The code employs separate grids for the metric solve and the fluid evolution because the two subsystems impose different geometric requirements.

| Subsystem | Domain structure | Spectral representation |
|---|---|---|
| Metric grid | spherical nucleus, multiple spherical shells, compactified external domain | Chebyshev in radius; spherical harmonics / Fourier in angle |
| Hydro grid | spherical nucleus, one or more shells, deformed outer boundary following the stellar surface | Chebyshev in radius; spherical harmonics / Fourier in angle |

The metric grid is used to solve the xCFC equations for \(\beta^i\), \(N\), \(\Psi_c\), and \(\Psi_\mathrm{ev}\), and its use of spherical domains is aligned with the requirements of the elliptic spacetime solver [2410.16764]. The hydrodynamic grid, by contrast, must conform to the actual stellar surface, which is non-spherical for rotating stars. The split-grid architecture is therefore not incidental; it is a direct response to the incompatibility between a metric solver that prefers spherical domains and a fluid domain that must adapt to stellar oblateness and oscillatory deformation.

The code further includes **origin regularization**, **axis regularization**, and **exponential filters** to suppress aliasing [2410.16764]. For axis regularization, the paper notes that scalar fields are transformed from Fourier representation to spherical-harmonic representation and back. These are standard but essential ingredients in spectral computations of rotating configurations, where coordinate singularities and modal contamination can otherwise dominate the error budget.

The most distinctive hydrodynamic feature is the **moving outer boundary**. The surface radius \(R(t,\theta,\phi)\) is evolved from the **impermeable-boundary condition** written as a time-evolution equation for the boundary itself [2410.16764]. The paper stresses why this matters: a fixed spherical outer boundary would poorly represent the matter-vacuum interface, whereas a surface-following boundary maintains a clean hydrodynamic domain and improves the representation of rotating equilibria and their perturbations. The appendix validates this boundary-evolution law with a rotating ellipsoid test [2410.16764].

## 4. Oscillation excitation, spectral diagnostics, and gravitational-wave extraction

Mode identification in ROXAS is performed from the **stellar surface deformation**, specifically through the **spherical-harmonic decomposition of the radius** [2410.16764]. The paper uses the following surface coefficients as diagnostics: **\(\ell=m=0\)** for quasi-radial modes, **\(\ell=2,m=0\)** for axisymmetric quadrupolar modes, and **\(\ell=m=2\)** for non-axisymmetric quadrupolar modes [2410.16764].

The perturbation used to excite oscillations is
$$
\delta H = (x^2 \pm y^2)\left(1-\frac{r^2}{R_S(\theta)^2}\right).
$$
The **plus sign** excites a mixture of **\(\ell=m=0\)** and **\(\ell=2,m=0\)** content, whereas the **minus sign** excites a mixture of **\(\ell=m=0\)** and **\(\ell=m=2\)** content [2410.16764]. Frequency extraction is then performed by Fourier transforming the time series of the selected surface-mode coefficient, with Fourier-space resolution
$$
\delta f = \frac{1}{t_\mathrm{max}}.
$$
The mode labels used in the paper are \(F\) for the fundamental quasi-radial mode, \(H_1\) for its first overtone, \({}^2f_0\) for the fundamental \(\ell=2,m=0\) mode, and \({}^2f_2\) and \({}^2f_{-2}\) for the co- and counterrotating \(\ell=m=2\) modes [2410.16764].

Although xCFC suppresses true radiative degrees of freedom, the code also computes **approximate gravitational-wave signals** via the quadrupole formula and its relativistic corrections. The paper writes
$$
h^{TT}_{ij}(\mathbf{x}, t) =
\frac{2G}{c^4 r}
P_{ij}^{kl}(\mathbf{n})
\ddot{Q}_{kl}\left(t-\frac{r}{c}\right),
$$
together with the stress formula
$$
\ddot{Q}_{ij} =
\int \rho(\mathbf{x},t)
\left(2\dot{x}_i\dot{x}_j
- x_i\partial_j\Phi
- x_j\partial_i\Phi\right)\, d^3x.
$$
The practical amplitudes \(A_+^e\) and \(A_\times^e\) are built using the weighted potential \(\Phi_w\) and the relativistic density
$$
\rho^* = \Psi^6 \Gamma \rho.
$$
These diagnostics do not convert ROXAS into a full waveform generator; rather, they provide a consistent approximate radiative readout of the fluid oscillations represented within the xCFC approximation [2410.16764].

## 5. Validation, recovered spectra, and accuracy claims

The principal validation strategy is comparison against **semi-analytic perturbative methods**, including **axisymmetric perturbative/Cowling studies**, **CFC results**, and **full GR results** for non-axisymmetric modes [2410.16764]. The general outcome reported is that most extracted frequencies agree to **about 1% or better**, while some low-frequency modes show larger relative discrepancies, up to **3–4%**, and in a few cases **more than 10%** [2410.16764]. The paper associates the larger errors with the **lowest-frequency modes** and with **short simulation lengths (~25 ms)**, which limit Fourier resolution.

For the realistic **SLy4-like** rotating star used as an example, the paper quotes
- \({}^2f_2\): ROXAS **\(1.675\) kHz** versus literature **\(1.655\) kHz**, a **1.2%** difference;
- \({}^2f_{-2}\): ROXAS **\(2.185\) kHz** versus literature **\(2.105\) kHz**, a **3.7%** difference [2410.16764].

The code also recovers **spherically symmetric modes** in the non-rotating **BU0** model. The extracted \(F\) and \(H_1\) frequencies are stated to be consistent with the earlier **spherically symmetric version** of the code [2410.16764]. This is methodologically important: it shows that the 2D/3D extension preserves continuity with the simpler radial implementation rather than introducing an incompatible discretization artifact.

The broader interpretive point is that these comparisons are used to support the adequacy of **CFC/xCFC for oscillation-frequency extraction** in this setting. The data do not claim exact equivalence to full GR, but they do indicate that the approximation is sufficiently accurate for the targeted mode studies while being computationally much cheaper [2410.16764].

## 6. Limitations, planned extensions, and terminological ambiguity

The current formulation is explicitly limited to **isolated stars**, **fast rigid rotation**, **cold barotropic EoS**, **\(\beta\)-equilibrium**, **no shocks**, and the **xCFC approximation rather than full GR** [2410.16764]. The paper also notes the use of **spatial symmetries** that reduce angular ranges, identifies **interpolation** as the main computational cost, and remarks that **longer runs** are needed to improve low-frequency precision [2410.16764]. These restrictions are not peripheral: they define the exact domain in which the reported accuracy statements should be interpreted.

Planned future extensions include **out-of-\(\beta\)-equilibrium matter**, **finite-temperature effects**, **relaxed spatial symmetries**, **beyond-CFC spacetime evolution**, and **differential rotation profiles relevant to HMNSs** [2410.16764]. This suggests a staged development trajectory in which the present code serves as a controlled high-accuracy baseline for fast-rotator oscillation calculations before moving toward more merger-realistic physics.

A plausible source of confusion is the suffix **“spec.”** In plasma physics, **SPEC** ordinarily denotes the **Stepped-Pressure Equilibrium Code**, extended to **field-aligned flow and rigid rotation** within **MRxMHD** [2001.06984]. In programming-language research, a **spec** is a behavior-descriptive entity combining **interface contracts**, **test rules**, and **test cases** [1107.3193]. Neither of those usages denotes ROXAS. This suggests that “fastrot-spec” is best treated not as a standardized acronym, but as an informal shorthand for the ROXAS fast-rotator oscillation setting [2410.16764].

Source: https://www.emergentmind.com/topics/fastrot-spec