skip to main content

Title: Comparing inclination-dependent analyses of kilonova transients
ABSTRACT The detection of the optical transient AT2017gfo proved that binary neutron star mergers are progenitors of kilonovae (KNe). Using a combination of numerical-relativity and radiative-transfer simulations, the community has developed sophisticated models for these transients for a wide portion of the expected parameter space. Using these simulations and surrogate models made from them, it has been possible to perform Bayesian inference of the observed signals to infer properties of the ejected matter. It has been pointed out that combining inclination constraints derived from the KN with gravitational-wave measurements increases the accuracy with which binary parameters can be estimated, in particular breaking the distance-inclination degeneracy from gravitational wave inference. To avoid bias from the unknown ejecta geometry, constraints on the inclination angle for AT2017gfo should be insensitive to the employed models. In this work, we compare different assumptions about the ejecta and radiative reprocesses used by the community and we investigate their impact on the parameter inference. While most inferred parameters agree, we find disagreement between posteriors for the inclination angle for different geometries that have been used in the current literature. According to our study, the inclusion of reprocessing of the photons between different ejecta types improves the modeling more » fits to AT2017gfo and, in some cases, affects the inferred constraints. Our study motivates the inclusion of large ∼ 1-mag uncertainties in the KN models employed for Bayesian analysis to capture yet unknown systematics, especially when inferring inclination angles, although smaller uncertainties seem appropriate to capture model systematics for other intrinsic parameters. We can use this method to impose soft constraints on the ejecta geometry of the KN AT2017gfo. « less
Authors:
; ; ; ; ; ; ; ; ;
Award ID(s):
1806990 1911206
Publication Date:
NSF-PAR ID:
10279145
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
502
Issue:
2
Page Range or eLocation-ID:
3057 to 3065
ISSN:
0035-8711
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT The joint detection of the gravitational wave GW170817, of the short γ-ray burst GRB170817A and of the kilonova AT2017gfo, generated by the the binary neutron star (NS) merger observed on 2017 August 17, is a milestone in multimessenger astronomy and provides new constraints on the NS equation of state. We perform Bayesian inference and model selection on AT2017gfo using semi-analytical, multicomponents models that also account for non-spherical ejecta. Observational data favour anisotropic geometries to spherically symmetric profiles, with a log-Bayes’ factor of ∼104, and favour multicomponent models against single-component ones. The best-fitting model is an anisotropic three-component composed ofmore »dynamical ejecta plus neutrino and viscous winds. Using the dynamical ejecta parameters inferred from the best-fitting model and numerical–relativity relations connecting the ejecta properties to the binary properties, we constrain the binary mass ratio to q < 1.54 and the reduced tidal parameter to $120\lt \tilde{\Lambda }\lt 1110$. Finally, we combine the predictions from AT2017gfo with those from GW170817, constraining the radius of a NS of 1.4 M⊙ to 12.2 ± 0.5 km (1σ level). This prediction could be further strengthened by improving kilonova models with numerical-relativity information.« less
  2. Context. The PSR J2222−0137 binary system has a set of features that make it a unique laboratory for tests of gravity theories. Aims. To fully exploit the system’s potential for these tests, we aim to improve the measurements of its physical parameters, spin and orbital orientation, and post-Keplerian parameters, which quantify the observed relativistic effects. Methods. We describe an improved analysis of archival very long baseline interferometry (VLBI) data, which uses a coordinate convention in full agreement with that used in timing. We have also obtained much improved polarimetry of the pulsar with the Five hundred meter Aperture Spherical Telescopemore »(FAST). We provide an improved analysis of significantly extended timing datasets taken with the Effelsberg, Nançay, and Lovell radio telescopes; this also includes previous timing data from the Green Bank Telescope. Results. From the VLBI analysis, we have obtained a new estimate of the position angle of the ascending node, Ω = 189 −18 +19 deg (all uncertainties are 68% confidence limits), and a new reference position for the pulsar with an improved and more conservative uncertainty estimate. The FAST polarimetric results, and in particular the detection of an interpulse, yield much improved estimates for the spin geometry of the pulsar, in particular an inclination of the spin axis of the pulsar of ∼84 deg. From the timing, we obtain a new ∼1% test of general relativity (GR) from the agreement of the Shapiro delay parameters and the rate of advance of periastron. Assuming GR in a self-consistent analysis of all effects, we obtain much improved masses: 1.831(10)  M ⊙ for the pulsar and 1.319(4)  M ⊙ for the white dwarf companion; the total mass, 3.150(14)  M ⊙ , confirms this as the most massive double degenerate binary known in the Galaxy. This analysis also yields the orbital orientation; in particular, the orbital inclination is 85.27(4) deg – indicating a close alignment between the spin of the pulsar and the orbital angular momentum – and Ω = 187.7(5.7) deg, which matches our new VLBI estimate. Finally, the timing also yields a precise measurement of the variation in the orbital period, Ṗ b = 0.251(8) × 10 −12 ss −1 ; this is consistent with the expected variation in the Doppler factor plus the orbital decay caused by the emission of gravitational waves predicted by GR. This agreement introduces stringent constraints on the emission of dipolar gravitational waves.« less
  3. Abstract Binary neutron star mergers (NSMs) have been confirmed as one source of the heaviest observable elements made by the rapid neutron-capture ( r -) process. However, modeling NSM outflows—from the total ejecta masses to their elemental yields—depends on the unknown nuclear equation of state (EOS) that governs neutron star structure. In this work, we derive a phenomenological EOS by assuming that NSMs are the dominant sources of the heavy element material in metal-poor stars with r -process abundance patterns. We start with a population synthesis model to obtain a population of merging neutron star binaries and calculate their EOS-dependentmore »elemental yields. Under the assumption that these mergers were responsible for the majority of r -process elements in the metal-poor stars, we find parameters representing the EOS for which the theoretical NSM yields reproduce the derived abundances from observations of metal-poor stars. For our proof-of-concept assumptions, we find an EOS that is slightly softer than, but still in agreement with, current constraints, e.g., by the Neutron Star Interior Composition Explorer, with R 1.4 = 12.25 ± 0.03 km and M TOV = 2.17 ± 0.03 M ⊙ (statistical uncertainties, neglecting modeling systematics).« less
  4. ABSTRACT Searches for gravitational-wave counterparts have been going in earnest since GW170817 and the discovery of AT2017gfo. Since then, the lack of detection of other optical counterparts connected to binary neutron star or black hole–neutron star candidates has highlighted the need for a better discrimination criterion to support this effort. At the moment, low-latency gravitational-wave alerts contain preliminary information about binary properties and hence whether a detected binary might have an electromagnetic counterpart. The current alert method is a classifier that estimates the probability that there is a debris disc outside the black hole created during the merger as wellmore »as the probability of a signal being a binary neutron star, a black hole–neutron star, a binary black hole, or of terrestrial origin. In this work, we expand upon this approach to both predict the ejecta properties and provide contours of potential light curves for these events, in order to improve the follow-up observation strategy. The various sources of uncertainty are discussed, and we conclude that our ignorance about the ejecta composition and the insufficient constraint of the binary parameters by low-latency pipelines represent the main limitations. To validate the method, we test our approach on real events from the second and third Advanced Laser Interferometer Gravitational-Wave Observatory (LIGO)–Virgo observing runs.« less
  5. ABSTRACT The structure of the broad-line region (BLR) is an essential ingredient in the determination of active galactic nucleus (AGN) virial black hole masses, which in turn are important to study the role of black holes in galaxy evolution. Constraints on the BLR geometry and dynamics can be obtained from velocity-resolved studies using reverberation mapping data (i.e. monitoring data). However, monitoring data are observationally expensive and only available for a limited sample of AGNs, mostly confined to the local Universe. Here, we explore a new version of a Bayesian inference, physical model of the BLR that uses an individual spectrummore »and prior information on the BLR size from the radius–luminosity relation, to model the AGN BLR geometry and dynamics. We apply our model to a sample of 11 AGNs, which have been previously modelled using monitoring data. Our single-epoch BLR model is able to constrain some of the BLR parameters with inferred parameter values that agree within the uncertainties with those determined from the modelling of monitoring data. We find that our model is able to derive stronger constraints on the BLR for AGNs with broad emission lines that qualitatively have more substructure and more asymmetry, presumably as they contain more information to constrain the physical model. The performance of this model makes it a practical and cost-effective tool to determine some of the BLR properties of a large sample of low- and high-redshift AGNs, for which monitoring data are not available.« less