Related papers: Elliptica: a new pseudo-spectral code for the cons…
The focus of this work is on the construction of initial data including a neutron star on a hyperboloidal slice. As simplest scenario for this first step, spherical symmetry is considered together with a polytropic-like equation of state…
We present a multi-domain spectral method to compute initial data of binary systems in General Relativity. By utilizing adapted conformal coordinates, the vacuum region exterior to the gravitational sources is divided up into two subdomains…
Construction of astrophysically realistic initial data remains a central problem when modelling the merger and eventual coalescence of binary black holes in numerical relativity. The objective of this paper is to provide astrophysically…
Spacetime is foliated by spatial hypersurfaces in the 3+1 split of General Relativity. The initial value problem then consists of specifying initial data for all relevant fields on one such a spatial hypersurface. These fields are the…
We solve the elliptic equations associated with the Hamiltonian and momentum constraints, corresponding to a system composed of two black holes with arbitrary linear and angular momentum. These new solutions are based on a Kerr-Schild…
In this paper we investigate the parabolic-hyperbolic formulation of the vacuum constraint equations introduced by R{\'a}cz with a view to constructing multiple black hole initial data sets without spin. In order to respect the natural…
Numerical relativity codes now being developed will evolve initial data representing colliding black holes at a relatively late stage in the collision. The choice of initial data used for code development has been made on the basis of…
The starting point of any general relativistic numerical simulation is a solution of the Hamiltonian and momentum constraints that (ideally) represents an astrophysically realistic scenario. We present a new method to produce initial data…
We present an approximate metric for a binary black hole spacetime to construct initial data for numerical relativity. This metric is obtained by asymptotically matching a post-Newtonian metric for a binary system to a perturbed…
Recent efforts to numerically simulate compact objects in alternative theories of gravity have largely focused on the time-evolution equations. Another critical aspect is the construction of constraint-satisfying initial data with precise…
Next-generation gravitational wave detectors such as Cosmic Explorer, the Einstein Telescope, and LISA, demand highly accurate and extensive gravitational wave (GW) catalogs to faithfully extract physical parameters from observed signals.…
We present approximate analytical solutions to the Hamiltonian and momentum constraint equations, corresponding to systems composed of two black holes with arbitrary linear and angular momentum. The analytical nature of these initial data…
The production of numerical relativity waveforms that describe quasicircular binary black hole mergers requires high-quality initial data, and an algorithm to iteratively reduce residual eccentricity. To date, these tools remain closed…
Standard puncture initial data have been widely used for numerical binary black hole evolutions despite their shortcomings, most notably the inherent lack of gravitational radiation at the initial time that is later followed by a burst of…
We calculate puncture initial data corresponding to both single and binary black hole solutions of the constraint equations by means of a pseudo-spectral method applied in a single spatial domain. Introducing appropriate coordinates, these…
We present a new scheme for constructing initial data for the Einstein field equations using the conformal thin-sandwich formulation that does not assume conformal flatness or approximate Killing vectors. This includes a method for…
Binary boson stars can be used to model the nonlinear dynamics and gravitational wave signals of merging ultracompact, but horizonless, objects. However, doing so requires initial data satisfying the Hamiltonian and momentum constraints of…
We present an algorithm for solving the general relativistic initial value equations for a corotating polytropic star in quasicircular orbit with a nonspinning black hole. The algorithm is used to obtain initial data for cases where the…
We present a new solver for coupled nonlinear elliptic partial differential equations (PDEs). The solver is based on pseudo-spectral collocation with domain decomposition and can handle one- to three-dimensional problems. It has three…
Novel applications of Numerical Relativity demand for more flexible algorithms and tools. In this paper, I develop and test a multigrid solver, based on the infrastructure provided by the Einstein Toolkit, for elliptic partial differential…