Related papers: Elliptica: a new pseudo-spectral code for the cons…
In this work, we introduce a spectral-infinite element method for solving Einstein's constraint equations in hyperbolic form. As an application of this, we use this method for computing asymptotically flat perturbations of a Kerr black hole…
We develop a novel technique through spectral decompositions to study the gravitational perturbations of a black hole, without needing to decouple the linearized field equations into master equations and separate their radial and angular…
We construct new models of black hole-neutron star binaries in quasiequilibrium circular orbits by solving Einstein's constraint equations in the conformal thin-sandwich decomposition together with the relativistic equations of…
A numerical solution scheme for the Einstein field equations based on generalized harmonic coordinates is described, focusing on details not provided before in the literature and that are of particular relevance to the binary black hole…
Numerical relativity is an essential tool in studying the coalescence of binary black holes (BBHs). It is still computationally prohibitive to cover the BBH parameter space exhaustively, making phenomenological fitting formulas for BBH…
We investigate an analytical framework for reconstructing bulk geometries from pole-skipping data. Previously, this method enabled the recursive recovery of near-horizon metric derivatives in static, planar-symmetric black holes. Building…
For highly compact neutron stars, constructing numerical initial data for black hole-neutron star binary evolutions is very difficult. We describe improvements to an earlier method that enable it to handle these more challenging cases. We…
Traditional black-hole binary puncture initial data is conformally flat. This unphysical assumption is coupled with a lack of radiation signature from the binary's past life. As a result, waveforms extracted from evolutions of this data…
It is well known that multigrid methods are optimally efficient for solution of elliptic equations (O(N)), which means that effort is proportional to the number of points at which the solution is evaluated). Thus this is an ideal method to…
In this work we have used for the first time pseudo-spectral methods to perform numerical simulations of spherically symmetric black hole formations on a Friedman-Robertson-Walker universe. With these methods, the differential equations…
Binary black hole simulations starting from quasi-circular (i.e., zero radial velocity) initial data have orbits with small but non-zero orbital eccentricities. In this paper the quasi-equilibrium initial-data method is extended to allow…
We present a new pseudo-spectral code for the simulation of evolution systems that are second order in space. We test this code by evolving a non-linear scalar wave equation. These non-linear waves can be stably evolved using very simple…
In a recent paper [{\em F. Bernal, J. Mor\'on-Vidal and J.A. Acebr\'on, Comp.$\&$ Math. App. 146:294-308 (2023)}] an hybrid supercomputing algorithm for elliptic equations has been put forward. The idea is that the interfacial nodal…
Using a post-Newtonian diagnostic tool developed by Mora and Will, we examine numerically generated quasiequilibrium initial data sets that have been used in recently successful numerical evolutions of binary black holes through plunge,…
We describe a generalised method for ellipsoid fitting against a minimum set of data points. The proposed method is numerically stable and applies to a wide range of ellipsoidal shapes, including highly elongated and arbitrarily oriented…
We present the first numerical simulations in null coordinates of the collapse of nonspherical regular initial data to a black hole. We restrict to twist-free axisymmetry, and re-investigate the critical collapse of a non-spherical massless…
We consider the initial data problem for several black holes in vacuum with arbitrary momenta and spins on a three space with punctures. We compactify the internal asymptotically flat regions to obtain a computational domain without inner…
We obtain an explicit solution of the momentum constraint for conformally flat, maximal slicing, initial data which gives an alternative to the purely longitudinal extrinsic curvature of Bowen and York. The new solution is related, in a…
Motivated by the goal for high accuracy modeling of gravitational radiation emitted by isolated systems, recently, there has been renewed interest in the numerical solution of the hyperboloidal initial value problem for Einstein's field…
By applying a parabolic-hyperbolic formulation of the constraints and superposing Kerr-Schild black holes, a simple method is introduced to initialize time evolution of binary systems. As the input parameters are essentially the same as…