No Arabic abstract
Through numerical simulations of boson-star head-on collisions, we explore the quality of binary initial data obtained from the superposition of single-star spacetimes. Our results demonstrate that evolutions starting from a plain superposition of individual boosted boson-star spacetimes are vulnerable to significant unphysical artefacts. These difficulties can be overcome with a simple modification of the initial data suggested in [PRD 99 (2018) 044046] for collisions of oscillatons. While we specifically consider massive complex scalar field boson star models up to a 6th-order-polynomial potential, we argue that this vulnerability is universal and present in other kinds of exotic compact systems and hence needs to be addressed.
We present improvements to construction of binary black hole initial data used in SpEC (the Spectral Einstein Code). We introduce new boundary conditions for the extended conformal thin sandwich elliptic equations that enforce the excision surfaces to be slightly inside rather than on the apparent horizons, thus avoiding extrapolation into the black holes at the last stage of initial data construction. We find that this improves initial data constraint violations near and inside the apparent horizons by about 3 orders of magnitude. We construct several initial data sets that are intended to be astrophysically equivalent but use different free data, boundary conditions, and initial gauge conditions. These include free data chosen as a superposition of two black holes in time-independent horizon-penetrating harmonic and damped harmonic coordinates. We also implement initial data for which the initial gauge satisfies the harmonic and damped harmonic gauge conditions; this can be done independently of the free data, since this amounts to a choice of the time derivatives of the lapse and shift. We compare these initial data sets by evolving them. We show that the gravitational waveforms extracted during the evolution of these different initial data sets agree very well after excluding initial transients. However, we do find small differences between these waveforms, which we attribute to small differences in initial orbital eccentricity, and in initial BH masses and spins, resulting from the different choices of free data. Among the cases considered, we find that superposed harmonic initial data leads to significantly smaller transients, smaller variation in BH spins and masses during these transients, smaller constraint violations, and more computationally efficient evolutions. Finally, we study the impact of initial data choices on the construction of zero-eccentricity initial data.
The construction of accurate and consistent initial data for various binary parameters is a critical ingredient for numerical relativity simulations of the compact binary coalescence. In this article, we present an upgrade of the pseudospectral SGRID code, which enables us to access even larger regions of the binary neutron star parameter space. As a proof of principle, we present a selected set of first simulations based on initial configurations computed with the new code version. In particular, we simulate two millisecond pulsars close to their breakup spin, highly compact neutron stars with masses at about $98%$ of the maximum supported mass of the employed equation of state, and an unequal mass systems with mass ratios even outside the range predicted by population synthesis models ($q = 2.03$). The discussed code extension will help us to simulate previously unexplored binary configurations. This is a necessary step to construct and test new gravitational wave approximants and to interpret upcoming binary neutron star merger observations. When we construct initial data, one has to specify various parameters, such as a rotation parameter for each star. Some of these parameters do not have direct physical meaning, which makes comparisons with other methods or models difficult. To facilitate this, we introduce simple estimates for the initial spin, momentum, mass, and center of mass of each individual star.
This paper reports on our effort in modeling realistic astrophysical neutron star binaries in general relativity. We analyze under what conditions the conformally flat quasiequilibrium (CFQE) approach can generate ``astrophysically relevant initial data, by developing an analysis that determines the violation of the CFQE approximation in the evolution of the binary described by the full Einstein theory. We show that the CFQE assumptions significantly violate the Einstein field equations for corotating neutron stars at orbital separations nearly double that of the innermost stable circular orbit (ISCO) separation, thus calling into question the astrophysical relevance of the ISCO determined in the CFQE approach. With the need to start numerical simulations at large orbital separation in mind, we push for stable and long term integrations of the full Einstein equations for the binary neutron star system. We demonstrate the stability of our numerical treatment and analyze the stringent requirements on resolution and size of the computational domain for an accurate simulation of the system.
Binary neutron stars in circular orbits can be modeled as helically symmetric, i.e., stationary in a rotating frame. This symmetry gives rise to a first integral of the Euler equation, often employed for constructing equilibrium solutions via iteration. For eccentric orbits, however, the lack of helical symmetry has prevented the use of this method, and the numerical relativity community has often resorted to constructing initial data by superimposing boosted spherical stars without solving the Euler equation. The spuriously excited neutron star oscillations seen in evolutions of such data arise because such configurations lack the appropriate tidal deformations and are stationary in a linearly comoving---rather than rotating---frame. We consider eccentric configurations at apoapsis that are instantaneously stationary in a rotating frame. We extend the notion of helical symmetry to eccentric orbits, by approximating the elliptical orbit of each companion as instantaneously circular, using the ellipses inscribed circle. The two inscribed helical symmetry vectors give rise to approximate instantaneous first integrals of the Euler equation throughout each companion. We use these integrals as the basis of a self-consistent iteration of the Einstein constraints to construct conformal thin-sandwich initial data for eccentric binaries. We find that the spurious stellar oscillations are reduced by at least an order of magnitude, compared with those found in evolutions of superposed initial data. The tidally induced oscillations, however, are physical and qualitatively similar to earlier evolutions. Finally, we show how to incorporate radial velocity due to radiation reaction in our inscribed helical symmetry vectors, which would allow one to obtain truly non-eccentric initial data when our eccentricity parameter $e$ is set to zero.
The spatially conformally flat approximation (CFA) is a viable method to deduce initial conditions for the subsequent evolution of binary neutron stars employing the full Einstein equations. Here we review the status of the original formulation of the CFA for the general relativistic hydrodynamic initial conditions of binary neutron stars. We illustrate the stability of the conformally flat condition on the hydrodynamics by numerically evolving ~100 quasi-circular orbits. We illustrate the use of this approximation for orbiting neutron stars in the quasi-circular orbit approximation to demonstrate the equation of state dependence of these initial conditions and how they might affect the emergent gravitational wave frequency as the stars approach the innermost stable circular orbit.