Pith. sign in

REVIEW 5 major objections 7 minor 97 references

Celephais: efficient spectral initial data code for precessing compact binaries

T0 review · 5 major / 7 minor · reviewed 2026-08-11 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read Celephais constructs spectrally accurate compact-binary initial data with arbitrary spin orientations, using a sparse, adaptively refined spectral solve that cuts unknowns by about three in high-mass-ratio benchmarks.

desk verdict A genuinely useful spectral-initial-data methods paper whose main caveat is an unquantified spin-curl approximation in the tilted-neutron-star hydrostatic equilibrium. read the letter →

arxiv 2608.07653 v1 pith:4KMRQLR2 submitted 2026-08-07 gr-qc astro-ph.HEastro-ph.IM

Hao-Jui Kuan This is my paper · ORCID
classification gr-qcastro-ph.HEastro-ph.IM MSC 83-0865N3565F5083C05 PACS 04.25.Dm04.30.-w04.40.Dg
keywords spectralmethodsinitialdatacompactbinariesprecessingspinsneutronstarsblackhole-neutronstaradaptivehp-refinementsparseJacobian
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

The paper introduces Celephais, a spectral initial-data code that solves the Einstein constraint equations for binary neutron stars and black-hole–neutron-star binaries without assuming the spins are aligned with the orbital axis. The aim is to make spectrally accurate, precessing compact-binary initial data cheap enough to produce in bulk for numerical-relativity surveys. To do so it assembles only the nonzero entries of the global multi-domain Jacobian, factors the sparse matrix once, and reuses that factorisation as a preconditioner in a matrix-free Newton–Krylov solve, then adaptively refines only the spectral directions that are underresolved. For a mass-ratio-20 black-hole–neutron-star benchmark, the adaptive scheme reaches the uniform-grid constraint accuracy with about three times fewer unknowns. The paper validates the data by binding-energy comparisons with post-Newtonian theory, an eccentricity-reduction test, and a full evolution whose waveform-reconstructed precession axis follows a post-Newtonian simple-precession model.

What carries the argument

The load-bearing object is the sparse, globally coupled multi-domain Jacobian of the extended conformal thin-sandwich residual system. The code builds it column-wise by scanning each residual's symbolic expression through a dependency graph, evaluating only residual rows whose supports contain the varied field–domain pair, and propagates several compatible Jacobian columns per forward-mode automatic-differentiation traversal. The resulting matrix is factored once by a sparse direct solver and reused as a right preconditioner in a matrix-free Newton–Krylov solve, with the factors refreshed periodically or when residual growth demands it; when the binary admits a reflection symmetry about the orbital plane with spins in that plane, the Jacobian block-diagonalises into even and odd parity sectors that are factored separately. Adaptive hp-refinement then inspects spectral-tail ratios in each domain and coordinate direction, marking pairs whose tail ratios exceed tolerances and choosing between radial h moves and anisotropic p moves. The final piece is the spin-aware post-Newtonian force balance: the orbital frequency is solved from the instantaneous radial balance of the nonspinning circular baseline plus the leading spin–orbit and spin–spin accelerations treated as vectors, with the nonradial part left to precessional dynamics.

What would settle it

Recompute one of the precessing binary-neutron-star configurations keeping the spin-curl term in the fluid Euler equation instead of dropping it, and compare the enthalpy or rest-mass-density distribution plus the measured eccentricity and waveform phase against the current data; a visible drift or a systematic tilt-dependent offset in the stellar profile would falsify the claim that arbitrary spin orientations are handled to the stated accuracy.

Watch

Extended reading notes

Core claim

The central claim is that the cost barrier to spectral compact-binary initial data with generic spin orientations can be broken by explicitly exploiting the sparsity of the Newton Jacobian rather than treating it as dense. The paper shows that a structurally filtered, batched automatic-differentiation assembly yields a sparse matrix with only a few percent nonzero entries for a representative precessing binary neutron star; factoring that matrix once with a sparse direct solver and reusing the factors as a refreshed right preconditioner keeps the memory bounded while the matrix-free Newton–Krylov iteration retains the full Jacobian action at every step. On top of this, an hp-refinement scheme driven by spectral-tail ratios concentrates resolution in the domains and directions that limit convergence, recovering the uniform-grid constraint accuracy with about three times fewer unknowns for a mass-ratio-20 black-hole–neutron-star system and about 1.6 times fewer for mass ratio 10. The author further extends the post-Newtonian force-balance condition to arbitrary spin orientations and shows that, for a precessing binary neutron star with one star tilted at 81 degrees, the post-Newtonian-initialised data have eccentricity below $10^{-2}$ and one evolution-based correction removes most of the orbital-frequency oscillation; the evolved waveform's reconstructed radiation axis tracks the 3.5PN simple-precession model through most of the inspiral.

Load-bearing premise

The neutron-star equilibrium in tilted-spin configurations is computed with an approximate Bernoulli relation that neglects the spin-curl term in the Euler equation; if spin-driven fluid currents significantly change the enthalpy balance, the stellar models are not truly quasiequilibrium.

Editorial extensions

If this is right

  • Precessing binary-neutron-star and black-hole–neutron-star initial data can be produced spectrally at laptop scale: for the production precessing binary-neutron-star case the first Jacobian assembly takes under 10 minutes and the full job uses under 17 GB of memory.
  • High-mass-ratio black-hole–neutron-star data can be brought to uniform-grid constraint accuracy with about three times fewer unknowns, and mass-ratio-10 data with about 1.6 times fewer, easing the cost of large parameter surveys.
  • Post-Newtonian force-balance initialisation gives precessing binary-neutron-star eccentricity below 10^-2 before correction and a usable starting point after one evolution-based iteration.
  • A full evolution of the precessing binary-neutron-star system is stable through merger, and the waveform-reconstructed precession axis follows the 3.5PN simple-precession model through most of the inspiral, linking prescribed spin tilts to observed dynamics.
  • The code is planned for public release under the GNU General Public License, which would let other groups reproduce the configurations and build extensions such as puncture black-hole data or alternative theories of gravity.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • If the approximate Bernoulli relation is the binding limitation, then retaining the spin-curl term in the fluid Euler equation would let the same machinery claim true quasiequilibrium for arbitrary spin tilts; the paper's own diagnostics (virial error, eccentricity, waveform phase) provide the template for testing that upgrade.
  • The sparse-Jacobian-plus-refreshed-preconditioner strategy is not tied to binary-neutron-star physics; the same dependency-filtered batched differentiation and parity-sector splitting could lower memory for other global spectral elliptic solvers in numerical relativity, such as puncture black-hole initial data or scalar-tensor field equations.
  • The aligned-spin approximation in the radial-infall estimate is the weakest part of the eccentricity-control pipeline for strongly precessing systems; a testable extension is a fully vectorial 3.5PN infall-rate estimate with the spin projection evolved along the precession cone.
  • The spectral-tail ratios used for adaptive refinement could double as an automated resolution-selection criterion in production waveform campaigns, since the paper's results suggest the final accuracy is largely independent of the refinement path once the limiting direction reaches sufficient order.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

5 major / 7 minor

Summary. The paper presents Celephais, a Kadath-based spectral initial-data code for binary neutron star (BNS) and black hole-neutron star (BHNS) systems with generic spin orientations. The key numerical innovations are a dependency-filtered, batched sparse Jacobian assembler; a MUMPS factorization reused as a refreshed right preconditioner in a Jacobian-free Newton-Krylov solve; and an anisotropic hp-refinement scheme driven by spectral tails. The validation suite includes compatibility with FUKA data, binding-energy comparisons with post-Newtonian sequences, a precessing BNS eccentricity-reduction test, and a full evolution of the same system from which the precession axis is reconstructed and compared with a 3.5PN simple-precession model. The main claims are that the code solves a representative precessing BNS on a laptop in minutes, that adaptive refinement recovers uniform-grid constraint accuracy with about three times fewer unknowns for a mass-ratio-20 BHNS, and that the code delivers spectrally resolved BNS/BHNS initial data with arbitrary spin orientations.

Significance. If the claims hold, Celephais would be a valuable public tool: it extends the Kadath/FUKA ecosystem toward generic spin orientations, addresses the memory bottleneck of dense spectral Jacobians, and provides a concrete AMR strategy for high-mass-ratio binaries. The paper is commendably concrete about its numerical machinery and about several of its own limitations, and it makes good use of external checks rather than relying on internal self-consistency alone. The efficiency benchmarks and the sparse-preconditioning strategy are plausible and well described. However, the headline claim of arbitrary spin orientations rests on an explicitly approximate treatment of spinning neutron-star hydrostatics whose error is never quantified, and the tilted-spin validation is limited to a single BNS configuration. These issues are load-bearing for the paper's central claim and require additional quantitative evidence.

major comments (5)
  1. [Sec. II A, Eqs. (17)-(21)] The approximate Bernoulli relation (21) drops the spin-curl term V^j(D_j \hat u_i - D_i \hat u_j) from Eq. (17). For the tilted-spin BNS of Table I (81 degrees spin-orbit tilt, chi=0.3), this term does not vanish by construction, yet the manuscript never estimates its magnitude or reports a stellar enthalpy-residual diagnostic. Because hydrostatic equilibrium of the spinning star is the physical basis for the claim of arbitrary spin orientations, please provide a quantitative bound on the neglected term for the configurations presented, e.g. by evaluating the full enthalpy balance (17) on the solved data or by comparing against a solve that retains the spin-curl term. The dynamical checks in Secs. V B-V D are not sensitive enough to exclude a small equilibrium bias, since a mildly non-equilibrium star can still inspiral and precess approximately as expected.
  2. [Sec. V A, Fig. 7] The binding-energy comparison with post-Newtonian sequences is presented as a validation of quasiequilibrium quality, but the vertical bars in Fig. 7 are the virial discrepancy (36), not statistical or systematic uncertainties in E_b. The numerical sequences are said to follow the PN trends, yet no resolution study of E_b or an error budget is given. To make this check quantitative, please report E_b as a function of resolution or provide an independent uncertainty estimate; with the present presentation, the agreement cannot be distinguished from a curve-fit effect.
  3. [Sec. IV A, Figs. 5 and 6] The adaptive-refinement gain stated as 'about three times fewer unknowns' is measured exclusively by the Hamiltonian-constraint norm ||H||_2. The XCTS system is a coupled set of elliptic equations, and the momentum constraint and the stellar-surface conditions are equally part of the solved problem. Reporting only the Hamiltonian norm may overstate the efficiency gain. Please report at least the momentum-constraint norm or the full nonlinear residual for the same refinement sequences.
  4. [Sec. V] No tilted-spin BHNS test is presented. The BHNS examples in Secs. IV and V A are either nonspinning or have aligned black-hole spin, while the only precessing configuration is the BNS of Table I and Secs. V B-V D. Since the abstract and conclusion claim arbitrary spin orientations for both BNS and BHNS initial data, the BHNS part of that claim is currently unvalidated. Please add at least one misaligned-spin BHNS case, or explicitly restrict the claim of arbitrary spin orientations to BNS systems.
  5. [Secs. V C and V D] The full evolution is performed at a single resolution, and the precession-axis comparison is conditional on the measured (2,2) frequency input to the PN model (Eqs. (59)-(65)). The manuscript explicitly acknowledges that waveform-phase convergence is not established, and that the simple-precession model is only valid in the slow-motion regime. This is acceptable as a usability check, but it means the comparison cannot isolate a possible bias in the initial stellar model. A second resolution and an early-time diagnostic of the fluid state (e.g. the enthalpy residual or the spin-curl term of Eq. (17)) would materially strengthen the claim that the generic-spin initial data are true quasiequilibrium data.
minor comments (7)
  1. [Sec. III A] The first paragraph contains typos: 'estalishing' should be 'establishing' and 'machienry' should be 'machinery'.
  2. [Sec. III C] The sentence 'Therefore, we intented to use one fourth of np' contains a typo: 'intented' should be 'intended'.
  3. [Sec. VI] The sentence 'the black-hole is treated with an excision horizon with' is grammatically incomplete; it appears to be missing a word or clause.
  4. [Fig. 5 caption] The caption contains a garbled phrase 'coe/cientenergy'; it should likely read 'coefficient energy'.
  5. [Fig. 10 legend] The legend text 'orbital-axisdirection(O'ShaughnessyhLaLbi)' is garbled and should refer cleanly to the O'Shaughnessy et al. radiation-axis construction.
  6. [Sec. II D / Table I] The text uses 'Macbook Pro' in one place and 'MacBook Pro' in Table I; please standardize the spelling.
  7. [Sec. V B, Eq. (51)] The precessing eccentricity estimate uses an aligned-spin approximation for \dot{a}, which the text acknowledges as a 'principal limitation'. Consider reporting how much this approximation affects the initial eccentricity for a moderately misaligned case, or at least note that the measured eccentricity below 10^-2 already includes this approximation.

Circularity Check

0 steps flagged · score 1.0 of 10

No significant circularity: the central solver and validation rest on independent equations, external codes, and transparently labeled calibration steps.

full rationale

No circular step can be exhibited from the paper's own text. The derivation chain is self-contained: the XCTS system (Eqs. 6-8), matter sources (Eqs. 14-16), the approximate Bernoulli relation (Eq. 21), and the boundary conditions define a fixed elliptic problem, and the code's performance is measured by constraint residuals and unknown counts, not by the quantities those residuals are used to produce. The validation is external at the points that matter: FUKA-imported solutions (Section V A) are a compatibility check against an independent public code, even though references [70,71] include the present author; PN binding-energy sequences (Fig. 7) are independent analytic estimates; and the precession-axis comparison (Section V D) explicitly uses the measured waveform frequency to drive the PN simple-precession model, stating that it does not independently predict frequency evolution. The eccentricity-reduction test (Section V B) is presented honestly as a calibration: PN estimates initialize, then 'dynamical measurements determine the final correction,' so the final low eccentricity is a measured outcome, not a prediction derived from the fitted parameters. The paper also discloses its limitations: the spin-curl term neglected in Eq. (21), the aligned-spin approximation for adot in precessing configurations, the absence of waveform-phase convergence, and the ill-conditioned spin-tilt inversion. These are modeling or validation caveats, not circular reductions. The only self-citations supply test data or prior applications and are not load-bearing for the central claim.

Assumptions & free parameters 4 free parameters · 5 assumptions · 0 invented entities

The central physics rests on standard XCTS assumptions and two modeling approximations (helical symmetry and the approximate Bernoulli relation); the efficiency claim depends on hand-chosen AMR thresholds. No new physical entities, fields, or conserved quantities are introduced.

free parameters (4)
  • AMR tail tolerances = tau2=1e-8, tau_inf=1e-7
    Default tolerances in Eq. (49); chosen by hand. The reported efficiency gains depend on these values but they are not fitted to external data.
  • AMR policy constants = dA=1, Gamma_p=2, Gamma_p,h=1.1, n_max=15
    Marking and execution-policy thresholds in Section IV B; hand-chosen to balance refinement cost. The 3x unknown-count claim depends on these choices.
  • Tail width w = w=2
    Tail definition in Eq. (47); arbitrary choice of the highest two modes for the truncation estimate.
  • Krylov and forcing bounds = eta in [1e-8,1e-3], m<=48
    Eisenstat-Walker forcing caps in Eq. (46); algorithmic parameters not fitted to data.
assumptions (5)
  • domain assumption Quasiequilibrium and approximate helical Killing vector
    Section II, before Eq. (11); binary assumed stationary in a corotating frame on orbital timescale. Standard in initial-data codes but an approximation for precessing binaries.
  • domain assumption Conformal flatness, maximal slicing, and instantaneous stationarity of freely specifiable fields
    Section II: tilde gamma_ij = f_ij, K=0, tilde u_ij=0, partial_t K=0. These reduce XCTS to the Isenberg-Wilson-Mathews form.
  • domain assumption Approximate Bernoulli relation with spin-curl terms neglected
    Section II A, Eq. (21); hydrostatic equilibrium for spinning neutron stars relies on this approximation.
  • domain assumption Parity-sector separation under y-reflection
    Section III B; assumes spin axes lie in the x-z plane and \ dot a = 0 for exact sectoring. For nonzero \ dot a, only the preconditioner uses the mask.
  • domain assumption p/rho rescaling gives smooth stellar surfaces
    Section II A; standard practice following FUKA/SGRID to avoid Gibbs phenomena at the stellar surface.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Celephais: efficient spectral initial data code for precessing compact binaries." pith.science (2026). https://pith.science/paper/4KMRQLR2

@misc{pith2026260807653,
  author       = {Pith},
  title        = {Pith review of: Celephais: efficient spectral initial data code for precessing compact binaries},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/4KMRQLR2}},
  note         = {Machine review of arXiv:2608.07653}
}
abstract

Large numerical-relativity surveys require compact-binary initial data that are both spectrally accurate and inexpensive to construct, including for systems with unequal masses and misaligned spins. We present \celephais, a compact-object initial-data code built on the \texttt{Kadath} spectral library, that constructs binary-neutron-star and black-hole--neutron-star initial data without imposing equatorial symmetry. The method exploits the sparse structure of the globally coupled multi-domain Jacobian and the approximate parity separation of fields. The assembled matrix is factored with \texttt{MUMPS} and reused as a refreshed right preconditioner in a Jacobian-free Newton--Krylov iteration, thereby avoiding dense storage. An adaptive $hp$--refinement scheme then concentrates resolution where the spectral tails are not yet resolved. For a mass-ratio-$20$ black-hole--neutron-star benchmark, the adaptive schemes recover the uniform-grid constraint accuracy with about three times fewer unknowns. We also extend the post-Newtonian force-balance estimate to arbitrary spin orientations and use it to initialise eccentricity reduction. Validation comprises binding-energy comparisons with post-Newtonian sequences, a precessing binary-neutron-star eccentricity-reduction test, and a full evolution whose waveform-reconstructed precession axis follows a post-Newtonian simple-precession model. These results establish an efficient route to spectrally resolved binary-neutron-star and black-hole--neutron-star initial data with arbitrary spin orientations.

Figures

Figures reproduced from arXiv: 2608.07653 by the authors.

Figure 1
Figure 1. FIG. 1. Binary multi-domain decomposition and admissible refinement moves (schematic, not to scale). Two surface-fitted stars are connected [PITH_FULL_IMAGE:figures/full_fig_p005_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2. Dependency filter for a residual operator [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. FIG. 3. Sparse Jacobian for a low-resolution equal-mass BNS with [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figures from the paper (7 more)
Figure 4
Figure 4. Figure 4: FIG. 4. The two admissible radial refinement moves on one domain [PITH_FULL_IMAGE:figures/full_fig_p007_4.png]
Figure 5
Figure 5. Figure 5: FIG. 5. Cost–accuracy relation for a mass-ratio- [PITH_FULL_IMAGE:figures/full_fig_p008_5.png]
Figure 6
Figure 6. Figure 6: FIG. 6. Local and global e [PITH_FULL_IMAGE:figures/full_fig_p008_6.png]
Figure 7
Figure 7. Figure 7: then compares quasicircular Celephaïs binding￾energy sequences with their PN estimates before any evolution￾based eccentricity correction. For the most demanding configu￾rations tested, agreement requires either a uniform N = 13 grid or an N = 11 baseline followed by o…
Figure 8
Figure 8. Figure 8: FIG. 8. Eccentricity reduction for the same precessing BNS as in [PITH_FULL_IMAGE:figures/full_fig_p010_8.png]
Figure 9
Figure 9. Figure 9: FIG. 9. Angle conventions for the single-spin precessing frame. At [PITH_FULL_IMAGE:figures/full_fig_p011_9.png]
Figure 10
Figure 10. Figure 10: FIG. 10. Radiation-axis precession reconstructed from the waveform [PITH_FULL_IMAGE:figures/full_fig_p012_10.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

97 extracted references · 27 canonical work pages

  1. [1]

    Abbottet al.(KAGRA, VIRGO, LIGO Scientific), GWTC-3: Compact Binary Coalescences Observed by LIGO and Virgo during the Second Part of the Third Observing Run, Phys

    R. Abbottet al.(KAGRA, VIRGO, LIGO Scientific), GWTC-3: Compact Binary Coalescences Observed by LIGO and Virgo during the Second Part of the Third Observing Run, Phys. Rev. X13, 041039 (2023), arXiv:2111.03606 [gr-qc]

    arXiv 2023

  2. [2]

    R. W. Kiendrebeogoet al., Updated Observing Scenarios and Multimessenger Implications for the International Gravitational- wave Networks O4 and O5, Astrophys. J.958, 158 (2023), arXiv:2306.09234 [astro-ph.HE]

    arXiv 2023

  3. [3]

    Spectral methods in general relativistic astrophysics

    S. Bonazzola, E. Gourgoulhon, and J. A. Marck, Spectral meth- ods in general relativistic astrophysics, J. Comput. Appl. Math. 109, 433 (1999), arXiv:gr-qc/9811089

    work page Pith review arXiv 1999

  4. [4]

    H. P. Pfeiffer, L. E. Kidder, M. A. Scheel, and S. A. Teukolsky, A Multidomain spectral method for solving elliptic equations, Comput. Phys. Commun.152, 253 (2003), arXiv:gr-qc/0202096

    arXiv 2003

  5. [5]

    Grandclement and J

    P. Grandclement and J. Novak, Spectral methods for numerical relativity, Living Rev. Rel.12, 1 (2009), arXiv:0706.2286 [gr- qc]

    arXiv 2009

  6. [6]

    LORENE website: Langage Objet pour la RElativité Numériqué, https://gitlab.in2p3.fr/lorene/Lorene

  7. [7]

    Gourgoulhon, P

    E. Gourgoulhon, P. Grandclement, K. Taniguchi, J.-A. Marck, and S. Bonazzola, Quasiequilibrium sequences of synchronized and irrotational binary neutron stars in general relativity: 1. Method and tests, Phys. Rev. D63, 064029 (2001), arXiv:gr- qc/0007028

    arXiv 2001

  8. [8]

    Grandclement, E

    P. Grandclement, E. Gourgoulhon, and S. Bonazzola, Binary black holes in circular orbits. 2. Numerical methods and first results, Phys. Rev. D65, 044021 (2002), arXiv:gr-qc/0106016

    arXiv 2002

Show all 97 references
  1. [9]

    Grandclement, Kadath: A Spectral solver for theoretical physics, J

    P. Grandclement, Kadath: A Spectral solver for theoretical physics, J. Comput. Phys.229, 3334 (2010), arXiv:0909.1228 [gr-qc]

    2010 arXiv

  2. [10]

    L. J. Papenfort, S. D. Tootle, P. Grandclément, E. R. Most, and L. Rezzolla, New public code for initial data of unequal-mass, spinning compact-object binaries, Phys. Rev. D104, 024057 (2021), arXiv:2103.09911 [gr-qc]

    2021 arXiv

  3. [11]

    Dietrich, N

    T. Dietrich, N. Moldenhauer, N. K. Johnson-McDaniel, S. Bernuzzi, C. M. Markakis, B. Brügmann, and W. Tichy, Binary Neutron Stars with Generic Spin, Eccentricity, Mass ratio, and Compactness - Quasi-equilibrium Sequences and First Evolutions, Phys. Rev. D92, 124007 (2015), arX...

    2015 arXiv

  4. [12]

    Tichy, A

    W. Tichy, A. Rashti, T. Dietrich, R. Dudi, and B. Brügmann, Constructing binary neutron star initial data with high spins, high compactnesses, and high mass ratios, Phys. Rev. D100, 124046 (2019), arXiv:1910.09690 [gr-qc]

    2019 arXiv

  5. [13]

    Rashti, F

    A. Rashti, F. M. Fabbri, B. Brügmann, S. V . Chaurasia, T. Diet- rich, M. Ujevic, and W. Tichy, New pseudospectral code for the construction of initial data, Phys. Rev. D105, 104027 (2022), arXiv:2109.14511 [gr-qc]

    2022 arXiv

  6. [14]

    Foucart, L

    F. Foucart, L. E. Kidder, H. P. Pfeiffer, and S. A. Teukolsky, Initial data for black hole-neutron star binaries: A Flexible, high-accuracy spectral method, Phys. Rev. D77, 124051 (2008), arXiv:0804.3787 [gr-qc]

    2008 arXiv

  7. [15]

    Taciket al., Binary Neutron Stars with Arbitrary Spins in Numerical Relativity, Phys

    N. Taciket al., Binary Neutron Stars with Arbitrary Spins in Numerical Relativity, Phys. Rev. D92, 124012 (2015), [Erratum: Phys.Rev.D 94, 049903 (2016)], arXiv:1508.06986 [gr-qc]. 14

    2015 arXiv

  8. [16]

    Assumpcao, L

    T. Assumpcao, L. R. Werneck, T. P. Jacques, and Z. B. Etienne, Fast hyperbolic relaxation elliptic solver for numerical relativity: Conformally flat, binary puncture initial data, Phys. Rev. D105, 104037 (2022), arXiv:2111.02424 [gr-qc]

    2022 arXiv

  9. [17]

    N. L. Vuet al., A scalable elliptic solver with task-based paral- lelism for the SpECTRE numerical relativity code, Phys. Rev. D 105, 084027 (2022), arXiv:2111.06767 [gr-qc]

    2022 arXiv

  10. [18]

    Vincent, H

    T. Vincent, H. P. Pfeiffer, and N. L. Fischer, hp-adaptive dis- continuous Galerkin solver for elliptic equations in numerical relativity, Phys. Rev. D100, 084052 (2019), arXiv:1907.01572 [physics.comp-ph]

    2019 arXiv

  11. [19]

    Gourgoulhon, 3+1 formalism and bases of numerical relativity (2007), arXiv:gr-qc/0703035

    E. Gourgoulhon, 3+1 formalism and bases of numerical relativity (2007), arXiv:gr-qc/0703035

    2007 arXiv

  12. [20]

    J. W. York, Jr., Conformal ’thin sandwich’ data for the initial- value problem, Phys. Rev. Lett.82, 1350 (1999), arXiv:gr- qc/9810051

    1999

  13. [21]

    H. P. Pfeiffer and J. W. York, Jr., Extrinsic curvature and the Einstein constraints, Phys. Rev. D67, 044022 (2003), arXiv:gr- qc/0207095

    2003

  14. [22]

    G. B. Cook, Corotating and irrotational binary black holes in quasicircular orbits, Phys. Rev. D65, 084003 (2002), arXiv:gr- qc/0108076

    2002

  15. [23]

    Bonazzola, E

    S. Bonazzola, E. Gourgoulhon, P. Grandclement, and J. Novak, A Constrained scheme for Einstein equations based on Dirac gauge and spherical coordinates, Phys. Rev. D70, 104007 (2004), arXiv:gr-qc/0307082

    2004 arXiv

  16. [24]

    J. A. Isenberg, Waveless approximation theories of gravity, Int. J. Mod. Phys. D17, 265 (2008), arXiv:gr-qc/0702113

    2008 arXiv

  17. [25]

    J. R. Wilson and G. J. Mathews, Instabilities in Close Neutron Star Binaries, Phys. Rev. Lett.75, 4161 (1995)

    1995

  18. [26]

    J. R. Wilson, G. J. Mathews, and P. Marronetti, Relativistic numerical model for close neutron star binaries, Phys. Rev. D 54, 1317 (1996), arXiv:gr-qc/9601017

    1996 arXiv

  19. [27]

    Shibata, K

    M. Shibata, K. Uryu, and J. L. Friedman, Deriving formula- tions for numerical computation of binary neutron stars in qua- sicircular orbits, Phys. Rev. D70, 044044 (2004), [Erratum: Phys.Rev.D 70, 129901 (2004)], arXiv:gr-qc/0407036

    2004 arXiv

  20. [28]

    K. Uryu, F. Limousin, J. L. Friedman, E. Gourgoulhon, and M. Shibata, Binary neutron stars in a waveless approximation, Phys. Rev. Lett.97, 171101 (2006), arXiv:gr-qc/0511136

    2006 arXiv

  21. [29]

    Bonazzola, E

    S. Bonazzola, E. Gourgoulhon, and J.-A. Marck, A Relativis- tic formalism to compute quasiequilibrium configurations of nonsynchronized neutron star binaries, Phys. Rev. D56, 7740 (1997), arXiv:gr-qc/9710031

    1997 arXiv

  22. [30]

    Tichy, B

    W. Tichy, B. Bruegmann, and P. Laguna, Gauge conditions for binary black hole puncture data based on an approximate helical Killing vector, Phys. Rev. D68, 064008 (2003), arXiv:gr- qc/0306020

    2003

  23. [31]

    Ossokine, F

    S. Ossokine, F. Foucart, H. P. Pfeiffer, M. Boyle, and B. Szilágyi, Improvements to the construction of binary black hole initial data, Class. Quant. Grav.32, 245010 (2015), arXiv:1506.01689 [gr-qc]

    2015 arXiv

  24. [32]

    Moldenhauer, C

    N. Moldenhauer, C. M. Markakis, N. K. Johnson-McDaniel, W. Tichy, and B. Brügmann, Initial data for binary neutron stars with adjustable eccentricity, Phys. Rev. D90, 084043 (2014), arXiv:1408.4136 [gr-qc]

    2014 arXiv

  25. [33]

    Buonanno, L

    A. Buonanno, L. E. Kidder, A. H. Mroue, H. P. Pfeiffer, and A. Taracchini, Reducing orbital eccentricity of precessing black- hole binaries, Phys. Rev. D83, 104034 (2011), arXiv:1012.1549 [gr-qc]

    2011 arXiv

  26. [34]

    Tichy, A New numerical method to construct binary neu- tron star initial data, Class

    W. Tichy, A New numerical method to construct binary neu- tron star initial data, Class. Quant. Grav.26, 175018 (2009), arXiv:0908.0620 [gr-qc]

    2009 arXiv

  27. [35]

    Tichy, Initial data for binary neutron stars with arbitrary spins, Phys

    W. Tichy, Initial data for binary neutron stars with arbitrary spins, Phys. Rev. D84, 024041 (2011), arXiv:1107.1440 [gr-qc]

    2011 arXiv

  28. [36]

    Tichy, Constructing quasi-equilibrium initial data for binary neutron stars with arbitrary spins, Phys

    W. Tichy, Constructing quasi-equilibrium initial data for binary neutron stars with arbitrary spins, Phys. Rev. D86, 064024 (2012), arXiv:1209.5336 [gr-qc]

    2012 arXiv

  29. [37]

    Shibata, A Relativistic formalism for computation of irrota- tional binary stars in quasiequilibrium states, Phys

    M. Shibata, A Relativistic formalism for computation of irrota- tional binary stars in quasiequilibrium states, Phys. Rev. D58, 024012 (1998), arXiv:gr-qc/9803085

    1998 arXiv

  30. [38]

    Gourgoulhon, Relations between three formalisms for irrota- tional binary neutron stars in general relativity (1998), arXiv:gr- qc/9804054

    E. Gourgoulhon, Relations between three formalisms for irrota- tional binary neutron stars in general relativity (1998), arXiv:gr- qc/9804054

    1998

  31. [39]

    Bonazzola, E

    S. Bonazzola, E. Gourgoulhon, and J.-A. Marck, Numerical models of irrotational binary neutron stars in general relativity, Phys. Rev. Lett.82, 892 (1999), arXiv:gr-qc/9810072

    1999 arXiv

  32. [40]

    Foucart, M

    F. Foucart, M. D. Duez, L. E. Kidder, and S. A. Teukolsky, Black hole-neutron star mergers: effects of the orientation of the black hole spin, Phys. Rev. D83, 024005 (2011), arXiv:1007.4203 [astro-ph.HE]

    2011 arXiv

  33. [41]

    Foucart, M

    F. Foucart, M. B. Deaton, M. D. Duez, L. E. Kidder, I. Mac- Donald, C. D. Ott, H. P. Pfeiffer, M. A. Scheel, B. Szilagyi, and S. A. Teukolsky, Black hole-neutron star mergers at realistic mass ratios: Equation of state and spin orientation effects, Phys. Rev. D87, 084006 (2013...

    2013 arXiv

  34. [42]

    Tsatsin and P

    P. Tsatsin and P. Marronetti, Initial data for neutron star bi- naries with arbitrary spins, Phys. Rev. D88, 064060 (2013), arXiv:1303.6692 [gr-qc]

    2013 arXiv

  35. [43]

    Kawaguchi, K

    K. Kawaguchi, K. Kyutoku, H. Nakano, H. Okawa, M. Shibata, and K. Taniguchi, Black hole-neutron star binary merger: Depen- dence on black hole spin orientation and equation of state, Phys. Rev. D92, 024014 (2015), arXiv:1506.05473 [astro-ph.HE]

    2015 arXiv

  36. [44]

    Dietrich, S

    T. Dietrich, S. Bernuzzi, B. Brügmann, M. Ujevic, and W. Tichy, Numerical Relativity Simulations of Precessing Bi- nary Neutron Star Mergers, Phys. Rev. D97, 064002 (2018), arXiv:1712.02992 [gr-qc]

    2018 arXiv

  37. [45]

    Kyutoku, S

    K. Kyutoku, S. Fujibayashi, K. Hayashi, K. Kawaguchi, K. Kiuchi, M. Shibata, and M. Tanaka, On the Possibility of GW190425 Being a Black Hole–Neutron Star Binary Merger, Astrophys. J. Lett.890, L4 (2020), arXiv:2001.04474 [astro- ph.HE]

    2020 arXiv

  38. [46]

    Bonazzola, E

    S. Bonazzola, E. Gourgoulhon, M. Salgado, and J. A. Marck, Axisymmetric rotating relativistic bodies: A new numerical ap- proach for ’exact’ solutions, Astron. Astrophys.278, 421 (1993)

    1993

  39. [47]

    Shibata, Rotating black hole surrounded by self-gravitating torus in the puncture framework, Phys

    M. Shibata, Rotating black hole surrounded by self-gravitating torus in the puncture framework, Phys. Rev. D76, 064035 (2007)

    2007

  40. [48]

    E. Gourgoulhon, An Introduction to the theory of rotating rel- ativistic stars, inCompStar 2010: School and Workshop on Computational Tools for Compact Star Astrophysics(2010) arXiv:1003.5015 [gr-qc]

    2010 arXiv

  41. [49]

    Gourgoulhon, P

    E. Gourgoulhon, P. Grandclement, and S. Bonazzola, Binary black holes in circular orbits. 1. A Global space-time approach, Phys. Rev. D65, 044020 (2002), arXiv:gr-qc/0106015

    2002 arXiv

  42. [50]

    G. B. Cook and H. P. Pfeiffer, Excision boundary conditions for black hole initial data, Phys. Rev. D70, 104016 (2004), arXiv:gr-qc/0407078

    2004 arXiv

  43. [51]

    J. L. Jaramillo, E. Gourgoulhon, and G. A. Mena Marugan, Inner boundary conditions for black hole initial data derived from isolated horizons, Phys. Rev. D70, 124036 (2004), arXiv:gr- qc/0407063

    2004

  44. [52]

    Grandclément and J

    P. Grandclément and J. Nicoules, Boundary conditions for sta- tionary black holes: Application to Kerr, Martínez-Troncoso- Zanelli, and hairy black holes, Phys. Rev. D105, 104011 (2022), arXiv:2203.09341 [gr-qc]. 15

    2022 arXiv

  45. [53]

    Arnowitt, S

    R. Arnowitt, S. Deser, and C. W. Misner, Energy and the Criteria for Radiation in General Relativity, Phys. Rev.118, 1100 (1960)

    1960

  46. [54]

    B. S. DeWitt, Quantum Theory of Gravity. 1. The Canonical Theory, Phys. Rev.160, 1113 (1967)

    1967

  47. [55]

    O Murchadha and J

    N. O Murchadha and J. W. York, Gravitational energy, Phys. Rev. D10, 2345 (1974)

    1974

  48. [56]

    J. W. York, Jr., Kinematics and Dynamics of General Relativity, inWorkshop on Sources of Gravitational Radiation(1978) pp. 83–126

    1978

  49. [57]

    R. L. Arnowitt, S. Deser, and C. W. Misner, The Dynamics of general relativity, Gen. Rel. Grav.40, 1997 (2008), arXiv:gr- qc/0405109

    2008

  50. [58]

    Regge and C

    T. Regge and C. Teitelboim, Role of Surface Integrals in the Hamiltonian Formulation of General Relativity, Annals Phys. 88, 286 (1974)

    1974

  51. [59]

    Gourgoulhon,3+1 Formalism in General Relativity, Lecture Notes in Physics (Springer, 2012)

    E. Gourgoulhon,3+1 Formalism in General Relativity, Lecture Notes in Physics (Springer, 2012)

    2012

  52. [60]

    Uryu and Y

    K. Uryu and Y . Eriguchi, A New numerical method for con- structing quasiequilibrium sequences of irrotational binary neu- tron stars in general relativity, Phys. Rev. D61, 124023 (2000), arXiv:gr-qc/9908059

    2000 arXiv

  53. [61]

    Ansorg, A

    M. Ansorg, A. Kleinwachter, and R. Meinel, Highly accurate calculation of rotating neutron stars: detailed description of the numerical methods, Astron. Astrophys.405, 711 (2003), arXiv:astro-ph/0301173

    2003 arXiv

  54. [62]

    Griewank and A

    A. Griewank and A. Walther,Evaluating Derivatives: Principles and Techniques of Algorithmic Differentiation, 2nd ed. (Society for Industrial and Applied Mathematics, Philadelphia, 2008)

    2008

  55. [63]

    Karypis and V

    G. Karypis and V . Kumar, A Fast and High Quality Multilevel Scheme for Partitioning Irregular Graphs, SIAM Journal on Scientific Computing20, 359 (1998)

    1998

  56. [64]

    R. S. Dembo, S. C. Eisenstat, and T. Steihaug, Inexact newton methods, SIAM Journal on Numerical Analysis19, 400 (1982), https://doi.org/10.1137/0719025

    1982 doi

  57. [65]

    S. C. Eisenstat and H. F. Walker, Choosing the forcing terms in an inexact newton method, SIAM Journal on Scientific Comput- ing17, 16 (1996), https://doi.org/10.1137/0917003

    1996 doi

  58. [66]

    Renkhoff, D

    S. Renkhoff, D. Cors, D. Hilditch, and B. Brügmann, Adaptive hp refinement for spectral elements in numerical relativity, Phys. Rev. D107, 104043 (2023), arXiv:2302.00575 [gr-qc]

    2023 arXiv

  59. [67]

    Ansorg, A Double-domain spectral method for black hole excision data, Phys

    M. Ansorg, A Double-domain spectral method for black hole excision data, Phys. Rev. D72, 024018 (2005), arXiv:gr- qc/0505059

    2005

  60. [68]

    Persson and J

    P.-O. Persson and J. Peraire, Sub-cell shock capturing for dis- continuous galerkin methods, in44th AIAA Aerospace Sciences Meeting and Exhibit(American Institute of Aeronautics and Astronautics, 2006) p. 112

    2006

  61. [69]

    Dörfler, A convergent adaptive algorithm for poisson’s equa- tion, SIAM Journal on Numerical Analysis33, 1106 (1996), https://doi.org/10.1137/0733054

    W. Dörfler, A convergent adaptive algorithm for poisson’s equa- tion, SIAM Journal on Numerical Analysis33, 1106 (1996), https://doi.org/10.1137/0733054

    1996 doi

  62. [70]

    H.-J. Kuan, K. Kiuchi, and M. Shibata, Tidal Resonance in Binary Neutron Star Inspirals: A High-Precision Study in Numerical Relativity, Phys. Rev. Lett.135, 141403 (2025), arXiv:2411.16850 [hep-ph]

    2025

  63. [71]

    H.-J. Kuan, I. Markin, M. Ujevic, T. Dietrich, K. Kiuchi, M. Shi- bata, and W. Tichy, The error budget of binary neutron star merger simulations for configurations with high spin (2025), arXiv:2506.02115 [gr-qc]

    2025 arXiv

  64. [72]

    Blanchet, Post-Newtonian Theory for Gravitational Waves, Living Rev

    L. Blanchet, Post-Newtonian Theory for Gravitational Waves, Living Rev. Rel.17, 2 (2014), arXiv:1310.1528 [gr-qc]

    2014 arXiv

  65. [73]

    L. E. Kidder, Coalescing binary systems of compact objects to postNewtonian 5/2 order. 5. Spin effects, Phys. Rev. D52, 821 (1995), arXiv:gr-qc/9506022

    1995 arXiv

  66. [74]

    A. Bohe, S. Marsat, G. Faye, and L. Blanchet, Next-to-next- to-leading order spin-orbit effects in the near-zone metric and precession equations of compact binaries, Class. Quant. Grav. 30, 075017 (2013), arXiv:1212.5520 [gr-qc]

    2013 arXiv

  67. [75]

    H. P. Pfeiffer, D. A. Brown, L. E. Kidder, L. Lindblom, G. Lovelace, and M. A. Scheel, Reducing orbital eccentricity in binary black hole simulations, Class. Quant. Grav.24, S59 (2007), arXiv:gr-qc/0702106

    2007 arXiv

  68. [76]

    Purrer, S

    M. Purrer, S. Husa, and M. Hannam, An Efficient iterative method to reduce eccentricity in numerical-relativity simula- tions of compact binary inspiral, Phys. Rev. D85, 124051 (2012), arXiv:1203.4258 [gr-qc]

    2012 arXiv

  69. [77]

    Ramos-Buades, S

    A. Ramos-Buades, S. Husa, and G. Pratten, Simple procedures to reduce eccentricity of binary black hole simulations, Phys. Rev. D99, 023003 (2019), arXiv:1810.00036 [gr-qc]

    2019 arXiv

  70. [78]

    Habib, M

    S. Habib, M. A. Scheel, and S. A. Teukolsky, Eccentricity re- duction for quasicircular binary evolutions, Phys. Rev. D111, 084059 (2025), arXiv:2410.05531 [gr-qc]

    2025 arXiv

  71. [79]

    Kyutoku, M

    K. Kyutoku, M. Shibata, and K. Taniguchi, Reducing orbital eccentricity in initial data of binary neutron stars, Phys. Rev. D 90, 064006 (2014), arXiv:1405.6207 [gr-qc]

    2014 arXiv

  72. [80]

    Kyutoku, K

    K. Kyutoku, K. Kawaguchi, K. Kiuchi, M. Shibata, and K. Taniguchi, Reducing orbital eccentricity in initial data of black hole–neutron star binaries in the puncture framework, Phys. Rev. D103, 023002 (2021), arXiv:2009.03896 [gr-qc]

    2021 arXiv

  73. [81]

    Blanchet, Innermost circular orbit of binary black holes at the third postNewtonian approximation, Phys

    L. Blanchet, Innermost circular orbit of binary black holes at the third postNewtonian approximation, Phys. Rev. D65, 124009 (2002), arXiv:gr-qc/0112056

    2002 arXiv

  74. [82]

    M.-Z. Han, K. Kiuchi, and M. Shibata, SACRA-K: A Performance-Portable Numerical Relativity Code with Kokkos (2026), arXiv:2607.08743 [astro-ph.HE]

    2026 arXiv

  75. [83]

    Yamamoto, M

    T. Yamamoto, M. Shibata, and K. Taniguchi, Simulating coa- lescing compact binaries by a new code SACRA, Phys. Rev. D 78, 064054 (2008), arXiv:0806.4007 [gr-qc]

    2008 arXiv

  76. [84]

    Kiuchi, K

    K. Kiuchi, K. Kawaguchi, K. Kyutoku, Y . Sekiguchi, M. Shibata, and K. Taniguchi, Sub-radian-accuracy gravitational waveforms of coalescing binary neutron stars in numerical relativity, Phys. Rev. D96, 084060 (2017), arXiv:1708.08926 [astro-ph.HE]

    2017 arXiv

  77. [85]

    Kiuchi, K

    K. Kiuchi, K. Kawaguchi, K. Kyutoku, Y . Sekiguchi, and M. Shi- bata, Sub-radian-accuracy gravitational waves from coalescing binary neutron stars in numerical relativity. II. Systematic study on the equation of state, binary mass, and mass ratio, Phys. Rev. D101, 084006 (2020...

    2020 arXiv

  78. [86]

    Kiuchi, L

    K. Kiuchi, L. E. Held, Y . Sekiguchi, and M. Shibata, Imple- mentation of advanced Riemann solvers in a neutrino-radiation magnetohydrodynamics code in numerical relativity and its ap- plication to a binary neutron star merger, Phys. Rev. D106, 124041 (2022), arXiv:2205.04487 ...

    2022 arXiv

  79. [87]

    Kiuchi, Toward high-precision inspiral gravitational wave- forms from binary neutron star mergers in numerical relativity, Phys

    K. Kiuchi, Toward high-precision inspiral gravitational wave- forms from binary neutron star mergers in numerical relativity, Phys. Rev. D112, 084008 (2025), arXiv:2508.10981 [astro- ph.HE]

    2025 arXiv

  80. [88]

    Nakano, J

    H. Nakano, J. Healy, C. O. Lousto, and Y . Zlochower, Perturbative extraction of gravitational waveforms generated with Numerical Relativity, Phys. Rev. D91, 104022 (2015), arXiv:1503.00718 [gr-qc]

    2015 arXiv

  81. [89]

    Reisswig and D

    C. Reisswig and D. Pollney, Notes on the integration of numeri- cal relativity waveforms, Class. Quant. Grav.28, 195015 (2011), arXiv:1006.1632 [gr-qc]

    2011 arXiv

  82. [90]

    O’Shaughnessy, B

    R. O’Shaughnessy, B. Vaishnav, J. Healy, Z. Meeks, and D. Shoe- maker, Efficient asymptotic frame selection for binary black hole spacetimes using asymptotic radiation, Phys. Rev. D84, 124002 (2011), arXiv:1109.5224 [gr-qc]. 16

    2011 arXiv

  83. [91]

    Kawaguchi, K

    K. Kawaguchi, K. Kyutoku, H. Nakano, and M. Shibata, Extract- ing the orbital axis from gravitational waves of precessing binary systems, Phys. Rev. D97, 024017 (2018), arXiv:1705.07459 [gr-qc]

    2018 arXiv

  84. [92]

    T. A. Apostolatos, C. Cutler, G. J. Sussman, and K. S. Thorne, Spin induced orbital precession and its modulation of the grav- itational wave forms from merging binaries, Phys. Rev. D49, 6274 (1994)

    1994

  85. [93]

    Eça and M

    L. Eça and M. Hoekstra, A procedure for the estimation of the numerical uncertainty of CFD calculations based on grid refinement studies, Journal of Computational Physics262, 104 (2014)

    2014

  86. [94]

    T. A. Oliver, N. Malaya, R. Ulerich, and R. D. Moser, Estimating uncertainties in statistics computed from direct numerical simu- lation, Physics of Fluids26, 035101 (2014), arXiv:1311.0828 [physics.flu-dyn]

    2014 arXiv

  87. [95]

    Dietrich, A

    T. Dietrich, A. Samajdar, S. Khan, N. K. Johnson-McDaniel, R. Dudi, and W. Tichy, Improving the NRTidal model for bi- nary neutron star systems, Phys. Rev. D100, 044003 (2019), arXiv:1905.06011 [gr-qc]

    2019 arXiv

  88. [96]

    Hotokezaka, K

    K. Hotokezaka, K. Kyutoku, H. Okawa, and M. Shibata, Explor- ing tidal effects of coalescing binary neutron stars in numerical relativity. II. Long-term simulations, Phys. Rev. D91, 064060 (2015), arXiv:1502.03457 [gr-qc]

    2015 arXiv

  89. [97]

    H.-J. Kuan, K. Van Aelst, A. T.-L. Lam, and M. Shibata, Binary neutron star mergers in massive scalar-tensor theory: Quasiequi- librium states and dynamical enhancement of the scalarization, Phys. Rev. D108, 064057 (2023), arXiv:2309.01709 [gr-qc]

    2023 arXiv

Pith tools

Reviewed August 11, 2026 · model on record in the stance chip above.