Search

We present nimbus: a hierarchical Bayesian framework to infer the intrinsic luminosity parameters of kilonovae (KNe) associated with gravitational-wave (GW) events, based purely on nondetections. This framework makes use of GW 3D distance information and electromagnetic upper limits from multiple surveys for multiple events and self-consistently accounts for the finite sky coverage and probability of astrophysical origin. The framework is agnostic to the brightness evolution assumed and can account for multiple electromagnetic passbands simultaneously. Our analyses highlight the importance of accounting for model selection effects, especially in the context of nondetections. We show our methodology using a simple, two-parameter linear brightness model, taking the follow-up of GW190425 with the Zwicky Transient Facility as a single-event test case for two different prior choices of model parameters: (i) uniform/uninformative priors and (ii) astrophysical priors based on surrogate models of Monte Carlo radiative-transfer simulations of KNe. We present results under the assumption that the KN is within the searched region to demonstrate functionality and the importance of prior choice. Our results show consistency with simsurveyâ€"an astronomical survey simulation tool used previously in the literature to constrain the population of KNe. While our results based on uniform priors strongly constrain the parameter space, those based on astrophysical priors are largely uninformative, highlighting the need for deeper constraints. Future studies with multiple events having electromagnetic follow-up from multiple surveys should make it possible to constrain the KN population further. [ABSTRACT FROM AUTHOR]

Searches for electromagnetic counterparts of gravitational-wave signals have redoubled since the first detection in 2017 of a binary neutron star merger with a gamma-ray burst, optical/infrared kilonova, and panchromatic afterglow. Yet, one LIGO/Virgo observing run later, there has not yet been a second, secure identification of an electromagnetic counterpart. This is not surprising given that the localization uncertainties of events in LIGO and Virgo's third observing run, O3, were much larger than predicted. We explain this by showing that improvements in data analysis that now allow LIGO/Virgo to detect weaker and hence more poorly localized events have increased the overall number of detections, of which well-localized, gold-plated events make up a smaller proportion overall. We present simulations of the next two LIGO/Virgo/KAGRA observing runs, O4 and O5, that are grounded in the statistics of O3 public alerts. To illustrate the significant impact that the updated predictions can have, we study the follow-up strategy for the Zwicky Transient Facility. Realistic and timely forecasting of gravitational-wave localization accuracy is paramount given the large commitments of telescope time and the need to prioritize which events are followed up. We include a data release of our simulated localizations as a public proposal planning resource for astronomers. [ABSTRACT FROM AUTHOR]

In the past few years, new observations of neutron stars (NSs) and NS mergers have provided a wealth of data that allow one to constrain the equation of state (EOS) of nuclear matter at densities above nuclear saturation density. However, most observations were based on NSs with masses of about 1.4 M ⊙, probing densities up to ∼three to four times the nuclear saturation density. Even higher densities are probed inside massive NSs such as PSR J0740+6620. Very recently, new radio observations provided an update to the mass estimate for PSR J0740+6620, and X-ray observations by the NICER and XMM telescopes constrained its radius. Based on these new measurements, we revisit our previous nuclear physics multimessenger astrophysics constraints and derive updated constraints on the EOS describing the NS interior. By combining astrophysical observations of two radio pulsars, two NICER measurements, the two gravitational-wave detections GW170817 and GW190425, detailed modeling of the kilonova AT 2017gfo, and the gamma-ray burst GRB 170817A, we are able to estimate the radius of a typical 1.4 M ⊙ NS to be 11.94 âˆ' 0.87 + 0.76 km at 90% confidence. Our analysis allows us to revisit the upper bound on the maximum mass of NSs and disfavors the presence of a strong first-order phase transition from nuclear matter to exotic forms of matter, such as quark matter, inside NSs. [ABSTRACT FROM AUTHOR]

We study a sample of 16 Type Ia supernovae (SNe Ia) having both spectroscopic and photometric observations within 2–3 days after the first light. The early B − V colors of such a sample tend to show a continuous distribution. For objects with normal ejecta velocity (NV), the C ii λ6580 feature is always visible in the early spectra, while it is absent or very weak in the high-velocity (HV) counterpart. Moreover, the velocities of the detached high-velocity features (HVFs) of the Ca II near-IR triplet (CaIR3) above the photosphere are found to be much higher in HV objects than in NV objects, with typical values exceeding 30,000 km s−1 at 2–3 days. We further analyze the relation between the velocity shift of late-time [Fe II ] lines (v[Fe II ]) and host galaxy mass. We find that all HV objects have redshifted v[Fe II ], while NV objects have both blue- and redshifted v[Fe II ]. It is interesting to point out that the objects with redshifted v[Fe II ] are all located in massive galaxies, implying that HV and a portion of NV objects may have similar progenitor metallicities and explosion mechanisms. We propose that, with a geometric/projected effect, the He-detonation model may account for the similarity in birthplace environment and the differences seen in some SNe Ia, including B − V colors, C II features, CaIR3 HVFs at early times, and v[Fe II ] in the nebular phase. Nevertheless, some features predicted by He-detonation simulation, such as the rapidly decreasing light curve, deviate from the observations, and some NV objects with blueshifted nebular v[Fe II ] may involve other explosion mechanisms. [ABSTRACT FROM AUTHOR]

We present a systematic search for optical counterparts to 13 gravitational wave (GW) triggers involving at least one neutron star during LIGO/Virgo's third observing run (O3). We searched binary neutron star (BNS) and neutron star black hole (NSBH) merger localizations with the Zwicky Transient Facility (ZTF) and undertook follow-up with the Global Relay of Observatories Watching Transients Happen (GROWTH) collaboration. The GW triggers had a median localization area of 4480 deg2, a median distance of 267 Mpc, and false-alarm rates ranging from 1.5 to 10−25 yr−1. The ZTF coverage in the g and r bands had a median enclosed probability of 39%, median depth of 20.8 mag, and median time lag between merger and the start of observations of 1.5 hr. The O3 follow-up by the GROWTH team comprised 340 UltraViolet/Optical/InfraRed (UVOIR) photometric points, 64 OIR spectra, and three radio images using 17 different telescopes. We find no promising kilonovae (radioactivity-powered counterparts), and we show how to convert the upper limits to constrain the underlying kilonova luminosity function. Initially, we assume that all GW triggers are bona fide astrophysical events regardless of false-alarm rate and that kilonovae accompanying BNS and NSBH mergers are drawn from a common population; later, we relax these assumptions. Assuming that all kilonovae are at least as luminous as the discovery magnitude of GW170817 (−16.1 mag), we calculate that our joint probability of detecting zero kilonovae is only 4.2%. If we assume that all kilonovae are brighter than −16.6 mag (the extrapolated peak magnitude of GW170817) and fade at a rate of 1 mag day−1 (similar to GW170817), the joint probability of zero detections is 7%. If we separate the NSBH and BNS populations based on the online classifications, the joint probability of zero detections, assuming all kilonovae are brighter than −16.6 mag, is 9.7% for NSBH and 7.9% for BNS mergers. Moreover, no more than <57% (<89%) of putative kilonovae could be brighter than −16.6 mag assuming flat evolution (fading by 1 mag day−1) at the 90% confidence level. If we further take into account the online terrestrial probability for each GW trigger, we find that no more than <68% of putative kilonovae could be brighter than −16.6 mag. Comparing to model grids, we find that some kilonovae must have Mej < 0.03 M⊙, Xlan > 10−4, or ϕ > 30° to be consistent with our limits. We look forward to searches in the fourth GW observing run; even 17 neutron star mergers with only 50% coverage to a depth of −16 mag would constrain the maximum fraction of bright kilonovae to <25%. [ABSTRACT FROM AUTHOR]

The first binary neutron star merger, GW170817, was accompanied by a radioactivity-powered optical/infrared transient called a kilonova. To date, no compelling kilonova has been found in all-sky optical surveys, independently of short gamma-ray burst and gravitational-wave triggers. In this work, we searched the first 23 months of the Zwicky Transient Facility (ZTF) data stream for candidate kilonovae in the form of rapidly evolving transients. We combined ZTF alert queries with forced point-spread-function photometry and nightly flux stacking to increase our sensitivity to faint and fast transients. Automatic queries yielded >11,200 candidates, 24 of which passed quality checks and selection criteria based on a grid of kilonova models tailored for both binary neutron star and neutron star–black hole mergers. None of the candidates in our sample was deemed a possible kilonova after thorough vetting. The sources that passed our selection criteria are dominated by Galactic cataclysmic variables. We identified two fast transients at high Galactic latitude, one of which is the confirmed afterglow of long-duration GRB 190106A, the other is a possible cosmological afterglow. Using a survey simulation code, we constrained the kilonova rate for a range of models including top-hat, linearly decaying light curves, and synthetic light curves obtained with radiative transfer simulations. For prototypical GW170817-like kilonovae, we constrain the rate to be Gpc−3 yr−1 (95% confidence). By assuming a population of kilonovae with the same geometry and composition of GW170817 observed under a uniform viewing angle distribution, we obtained a constraint on the rate of R < 4029 Gpc−3 yr−1. [ABSTRACT FROM AUTHOR]

Colors of Type Ia supernovae (SNe Ia) in the first few days after explosion provide a potential discriminant between different models. In this paper, we present g − r colors of 65 SNe Ia discovered within 5 days from first light by the Zwicky Transient Facility in 2018, a sample that is about three times larger than that in the literature. We find that g − r colors are intrinsically rather homogeneous at early phases, with about half of the dispersion attributable to photometric uncertainties (0.18 mag). Colors are nearly constant starting from 6 days after first light (g − r ∼ −0.15 mag), while the time evolution at earlier epochs is characterized by a continuous range of slopes, from events rapidly transitioning from redder to bluer colors (slope of ∼−0.25 mag day−1) to events with a flatter evolution. The continuum in the slope distribution is in good agreement both with models requiring some amount of 56Ni mixed in the outermost regions of the ejecta and with "double-detonation" models having thin helium layers () and varying carbon–oxygen core masses. At the same time, six events show evidence for a distinctive "red bump" signature predicted by double-detonation models with larger helium masses. We finally identify a significant correlation between the early-time g − r slopes and supernova brightness, with brighter events associated to flatter color evolution (p-value = 0.006). The distribution of slopes, however, is consistent with being drawn from a single population, with no evidence for two components as claimed in the literature based on B − V colors. [ABSTRACT FROM AUTHOR]

We present our spectropolarimetric observations of SN 2017egm, a Type I superluminous supernova (SLSN-I) in a nearby galaxy NGC 3191, with the Subaru telescope at +185.0 days after the g-band maximum light. This is the first spectropolarimetric observation for SLSNe at late phases. We find that the degree of the polarization in the late phase significantly changes from that measured at the earlier phase. The spectrum at the late phase shows a strong Ca emission line and therefore we reliably estimate the interstellar polarization (ISP) component assuming that the emission line is intrinsically unpolarized. By subtracting the estimated ISP, we find that the intrinsic polarization at the early phase is only ∼0.2%, which indicates an almost spherical photosphere, with an axial ratio ∼1.05. The intrinsic polarization at the late phase increases to ∼0.8%, which corresponds to the photosphere with an axial ratio ∼1.2. A nearly constant position angle of the polarization suggests the inner ejecta are almost axisymmetric. By these observations, we conclude that the inner ejecta are more aspherical than the outer ejecta. This may suggest the presence of a central energy source producing aspherical inner ejecta. [ABSTRACT FROM AUTHOR]