skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Analysis of the Conformally Flat Approximation for Binary Neutron Star Initial Conditions

Abstract

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 in this paper, we analyze the viability 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 quasicircular orbits. We illustrate the use of this approximation for orbiting neutron stars in the quasicircular 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.

Authors:
 [1]; ORCiD logo [2];  [2];  [3]
  1. Univ. of Notre Dame, IN (United States). Center for Astrophysics, Dept. of Physics and Center for Research Computing
  2. Univ. of Notre Dame, IN (United States). Center for Astrophysics, Dept. of Physics
  3. Univ. of Notre Dame, IN (United States). Center for Astrophysics, Dept. of Physics ; Hanoi National Univ. of Education (Vietnam); Univ. of Notre Dame, IN (United States). Joint Inst. for Nuclear Astrophysics (JINA)
Publication Date:
Research Org.:
Univ. of Notre Dame, IN (United States)
Sponsoring Org.:
USDOE
OSTI Identifier:
1366506
Grant/Contract Number:
FG02-95ER40934
Resource Type:
Journal Article: Accepted Manuscript
Journal Name:
Advances in Astronomy
Additional Journal Information:
Journal Volume: 2017; Journal ID: ISSN 1687-7969
Country of Publication:
United States
Language:
English
Subject:
79 ASTRONOMY AND ASTROPHYSICS

Citation Formats

Suh, In-Saeng, Mathews, Grant J., Haywood, J. Reese, and Lan, N. Q.. Analysis of the Conformally Flat Approximation for Binary Neutron Star Initial Conditions. United States: N. p., 2017. Web. doi:10.1155/2017/6127031.
Suh, In-Saeng, Mathews, Grant J., Haywood, J. Reese, & Lan, N. Q.. Analysis of the Conformally Flat Approximation for Binary Neutron Star Initial Conditions. United States. doi:10.1155/2017/6127031.
Suh, In-Saeng, Mathews, Grant J., Haywood, J. Reese, and Lan, N. Q.. Mon . "Analysis of the Conformally Flat Approximation for Binary Neutron Star Initial Conditions". United States. doi:10.1155/2017/6127031. https://www.osti.gov/servlets/purl/1366506.
@article{osti_1366506,
title = {Analysis of the Conformally Flat Approximation for Binary Neutron Star Initial Conditions},
author = {Suh, In-Saeng and Mathews, Grant J. and Haywood, J. Reese and Lan, N. Q.},
abstractNote = {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 in this paper, we analyze the viability 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 quasicircular orbits. We illustrate the use of this approximation for orbiting neutron stars in the quasicircular 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.},
doi = {10.1155/2017/6127031},
journal = {Advances in Astronomy},
number = ,
volume = 2017,
place = {United States},
year = {Mon Jan 09 00:00:00 EST 2017},
month = {Mon Jan 09 00:00:00 EST 2017}
}

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record

Save / Share:
  • Using the technique of the conformal deformation of initial data, we study the problem of mapping flat data onto solutions of Ashtekar's constraints. The procedure followed is that of York and we show that it is useful even for complex general relativity.
  • A new numerical method to construct binary black hole/neutron star initial data is presented. The method uses three spherical coordinate patches; two of these are centered at the binary compact objects and cover a neighborhood of each object; the third patch extends to the asymptotic region. As in the Komatsu-Eriguchi-Hachisu method, nonlinear elliptic field equations are decomposed into a flat space Laplacian and a remaining nonlinear expression that serves in each iteration as an effective source. The equations are solved iteratively, integrating a Green's function against the effective source at each iteration. Detailed convergence tests for the essential part ofmore » the code are performed for a few types of selected Green's functions to treat different boundary conditions. Numerical computation of the gravitational potential of a fluid source, and a toy model for a binary black hole field, are carefully calibrated with the analytic solutions to examine accuracy and convergence of the new code. As an example of the application of the code, an initial data set for binary black holes in the Isenberg-Wilson-Mathews formulation is presented, in which the apparent horizons are located using a method described in Appendix 1.« less
  • By asymptotically matching a post-Newtonian (PN) metric to two perturbed Schwarzschild metrics, we generate approximate initial data (in the form of an approximate 4-metric) for a nonspinning black hole binary in a circular orbit. We carry out this matching through O(v{sup 4}) in the binary's orbital velocity v, and thus the resulting data, like the O(v{sup 4}) PN metric, are conformally curved. The matching procedure also fixes the quadrupole and octupole tidal deformations of the holes, including the 1PN corrections to the quadrupole fields. Far from the holes, we use the appropriate PN metric that accounts for retardation, which wemore » construct using the highest-order PN expressions available to compute the binary's past history. The data set's uncontrolled remainders are thus O(v{sup 5}) throughout the time slice; we also generate an extension to the data set that has uncontrolled remainders of O(v{sup 6}) in the purely PN portion of the time slice (i.e., not too close to the holes). This extension also includes various other readily available higher-order terms. The addition of these terms decreases the constraint violations in certain regions, even though it does not increase the data's formal accuracy. The resulting data are smooth, since we join all the metrics together by smoothly interpolating between them. We perform this interpolation using transition functions constructed to avoid introducing excessive additional constraint violations. Because of their inclusion of tidal deformations and outgoing radiation, these data should substantially reduce both the high- and low-frequency components of the initial spurious ('junk') radiation observed in current simulations that use conformally flat initial data. Such reductions in the nonphysical components of the initial data will be necessary for simulations to achieve the accuracy required to supply Advanced LIGO and LISA with the templates necessary for parameter estimation.« less
  • The central engine of short gamma-ray bursts (sGRBs) is hidden from direct view, operating at a scale much smaller than that probed by the emitted radiation. Thus we must infer its origin not only with respect to the formation of the trigger—the actual astrophysical configuration that is capable of powering an sGRB—but also from the consequences that follow from the various evolutionary pathways that may be involved in producing it. Considering binary neutron star mergers we critically evaluate, analytically and through numerical simulations, whether the neutrino-driven wind produced by the newly formed hyper-massive neutron star can allow the collimated relativisticmore » outflow that follows its collapse to actually produce an sGRB or not. Upon comparison with the observed sGRB duration distribution, we find that collapse cannot be significantly delayed (≤100 ms) before the outflow is choked, thus limiting the possibility that long-lived hyper-massive remnants can account for these events. In the case of successful breakthrough of the jet through the neutrino-driven wind, the energy stored in the cocoon could contribute to the precursor and extended emission observed in sGRBs.« less
  • The initial mass function (IMF), binary fraction, and distributions of binary parameters (mass ratios, separations, and eccentricities) are indispensable inputs for simulations of stellar populations. It is often claimed that these are poorly constrained, significantly affecting evolutionary predictions. Recently, dedicated observing campaigns have provided new constraints on the initial conditions for massive stars. Findings include a larger close binary fraction and a stronger preference for very tight systems. We investigate the impact on the predicted merger rates of neutron stars and black holes. Despite the changes with previous assumptions, we only find an increase of less than a factor ofmore » 2 (insignificant compared with evolutionary uncertainties of typically a factor of 10–100). We further show that the uncertainties in the new initial binary properties do not significantly affect (within a factor of 2) our predictions of double compact object merger rates. An exception is the uncertainty in IMF (variations by a factor of 6 up and down). No significant changes in the distributions of final component masses, mass ratios, chirp masses, and delay times are found. We conclude that the predictions are, for practical purposes, robust against uncertainties in the initial conditions concerning binary parameters, with the exception of the IMF. This eliminates an important layer of the many uncertain assumptions affecting the predictions of merger detection rates with the gravitational wave detectors aLIGO/aVirgo.« less