Search

An open question in SN Ia research is where the boundary lies between 'normal' Type Ia supernovae (SNe Ia) that are used in cosmological measurements and those that sit off the Phillips relation. We present the spectroscopic modelling of one such '86G-like' transitional SN Ia, SN 2021rhu, that has recently been employed as a local Hubble Constant calibrator using a tip of the red-giant branch measurement. We detail its modelling from -12 d until maximum brightness using the radiative-transfer spectral-synthesis code tardis. Please check and correct this paper accordingly. We base our modelling on literature delayed-detonation and deflagration models of Chandrasekhar mass white dwarfs, as well as the double-detonation models of sub-Chandrasekhar mass white dwarfs. We present a new method for 'projecting' abundance profiles to different density profiles for ease of computation. Due to the small velocity extent and low outer densities of the W7 profile, we find it inadequate to reproduce the evolution of SN 2021rhu as it fails to match the high-velocity calcium components. The host extinction of SN 2021rhu is uncertain but we use modelling with and without an extinction correction to set lower and upper limits on the abundances of individual species. Comparing these limits to literature models we conclude that the spectral evolution of SN 2021rhu is also incompatible with double-detonation scenarios, lying more in line with those resulting from the delayed-detonation mechanism (although there are some discrepancies, in particular a larger titanium abundance in SN 2021rhu compared to the literature). This suggests that SN 2021rhu is likely a lower luminosity, and hence lower temperature, version of a normal SN Ia.

A merger of binary neutron stars creates heavy unstable elements whose radioactive decay produces a thermal emission known as a kilonova. In this paper, we predict the photometric and polarimetric behaviour of this emission by performing 3D Monte Carlo radiative transfer simulations. In particular, we choose three hydrodynamical models for merger ejecta, two including jets with different luminosities and one without a jet structure, to help decipher the impact of jets on the light curve and polarimetric behaviour. In terms of photometry, we find distinct colour evolutions across the three models. Models without a jet show the highest variation in light curves for different viewing angles. In contrast to previous studies, we find models with a jet to produce fainter kilonovae when viewed from orientations close to the jet axis, compared to a model without a jet. In terms of polarimetry, we predict relatively low levels (similar to 0.3-0.4 per cent) at all orientations that, however, remain non-negligible until a few days after the merger and longer than previously found. Despite the low levels, we find that the presence of a jet enhances the degree of polarization at wavelengths ranging from 0.25 to 2.5 mu m, an effect that is found to increase with the jet luminosity. Thus, future photometric and polarimetric campaigns should observe kilonovae in blue and red filters for a few days after the merger to help constrain the properties of the ejecta (e.g. composition) and jet.

LIGO and Virgo's third observing run revealed the first neutron star-black hole (NSBH) merger candidates in gravitational waves. These events are predicted to synthesize r-process elements(1,2)creating optical/near-infrared 'kilonova' emission. The joint gravitational wave and electromagnetic detection of an NSBH merger could be used to constrain the equation of state of dense nuclear matter(3), and independently measure the local expansion rate of the Universe(4). Here, we present the optical follow-up and analysis of two of the only three high-significance NSBH merger candidates detected to date, S200105ae and S200115j, with the Zwicky Transient Facility(5). The Zwicky Transient Facility observed similar to 48% of S200105ae and similar to 22% of S200115j's localization probabilities, with observations sensitive to kilonovae brighter than -17.5 mag fading at 0.5 mag d(-1)in the g- and r-bands; extensive searches and systematic follow-up of candidates did not yield a viable counterpart. We present state-of-the-art kilonova models tailored to NSBH systems that place constraints on the ejecta properties of these NSBH mergers. We show that with observed depths of apparent magnitude similar to 22 mag, attainable in metre-class, wide-field-of-view survey instruments, strong constraints on ejecta mass are possible, with the potential to rule out low mass ratios, high black hole spins and large neutron star radii.

The rise of multimessenger astronomy has brought with it the need to exploit all available data streams and learn more about the astrophysical objects that fall within its breadth. One possible avenue is the search for serendipitous optical/near-infrared counterparts of gamma-ray bursts (GRBs) and gravitational-wave (GW) signals, known as kilonovae. With surveys such as the Zwicky Transient Facility (ZTF), which observes the sky with a cadence of ∼3 d, the existing counterpart locations are likely to be observed; however, due to the significant amount of sky to explore, it is difficult to search for these fast-evolving candidates. Thus, it is beneficial to optimize the survey cadence for realtime kilonova identification and enable further photometric and spectroscopic observations. We explore how the cadence of wide field-of-view surveys like ZTF can be improved to facilitate such identifications. We show that with improved observational choices, e.g. the adoption of three epochs per night on a ∼ nightly basis, and the prioritization of redder photometric bands, detection efficiencies improve by about a factor of two relative to the nominal cadence. We also provide realistic hypothetical constraints on the kilonova rate as a form of comparison between strategies, assuming that no kilonovae are detected throughout the long-term execution of the respective observing plan. These results demonstrate how an optimal use of ZTF increases the likelihood of kilonova discovery independent of GWs or GRBs, thereby allowing for a sensitive search with less interruption of its nominal cadence through Target of Opportunity programs.

Gamma-ray bursts (GRBs) are among the brightest and most energetic events in the Universe. The duration and hardness distribution of GRBs has two clusters(1), now understood to reflect (at least) two different progenitors(2). Short-hard GRBs (SGRBs; T-90 < 2 s) arise from compact binary mergers, and long-soft GRBs (LGRBs; T-90 > 2 s) have been attributed to the collapse of peculiar massive stars (collapsars)(3). The discovery of SN 1998bw/GRB 980425 (ref. (4)) marked the first association of an LGRB with a collapsar, and AT 2017gfo (ref. (5))/GRB 170817A/GW170817 (ref. (6)) marked the first association of an SGRB with a binary neutron star merger, which also produced a gravitational wave. Here, we present the discovery of ZTF20abwysqy (AT2020scz), a fast-fading optical transient in the Fermi satellite and the Interplanetary Network localization regions of GRB 200826A; X-ray and radio emission further confirm that this is the afterglow. Follow-up imaging (at rest-frame 16.5 days) reveals excess emission above the afterglow that cannot be explained as an underlying kilonova, but which is consistent with being the supernova. Although the GRB duration is short (rest-frame T-90 of 0.65 s), our panchromatic follow-up data confirm a collapsar origin. GRB 200826A is the shortest LGRB found with an associated collapsar; it appears to sit on the brink between a successful and a failed collapsar. Our discovery is consistent with the hypothesis that most collapsars fail to produce ultra-relativistic jets.

The detection of the binary neutron star (BNS) merger GW170817 and the associated electromagnetic (EM) counterpart, the ‘kilonova’ (kN) AT2017gfo, opened a new era in multimessenger astronomy. However, despite many efforts, it has been proven very difficult to find additional kNe, even though LIGO/Virgo has reported at least one BNS event during their latest run, O3. The focus of this work is the exploration of the sensitivity of the adopted optical surveys searching for kNe during O3. We propose ways to optimize the choices of filters and survey depth to boost the detection efficiency for these faint and fast-evolving transients in the future. In particular, we use kN models to explore the dependence on ejecta mass, geometry, viewing angle, wavelength coverage, and source distance. We find that the kN detection efficiency has a strong viewing-angle dependence, especially for filters blueward of i-band. This loss of sensitivity can be mitigated by early, deep, observations. Efficient gri counterpart searches for kNe at ∼200 Mpc would require reaching a limiting magnitude mlim = 23 mag, to ensure good sensitivity over a wide range of the model phase-space. We conclude that kN searches during O3 were generally too shallow to detect BNS optical counterparts, even under optimistic assumptions.

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 fits to AT2017gfo and, in some cases, affects the inferred constraints. Our study motivates the inclusion of large similar to 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.

Massive stars die an explosive death as a core-collapse supernova (CCSN). The exact physical processes that cause the collapsing star to rebound into an explosion are not well understood, and the key to resolving this issue may lie in the measurement of the shape of CCSNe ejecta. Spectropolarimetry is the only way to perform this measurement for CCSNe outside the Milky Way and Magellanic Clouds. We present the infrared spectropolarimetric detection of a CCSN enabled by the new highly sensitive WIRC+Pol instrument at Palomar Observatory, which can observe CCSNe (magnitude M = −17 mag) out to 20 Mpc at ~0.1% polarimetric precision. Infrared spectropolarimetry is less affected than optical spectropolarimetry by dust scattering in the circumstellar and interstellar media, thereby providing a less biased probe of the intrinsic geometry of the supernova ejecta. SN 2018hna, a SN 1987A-like explosion, shows 2.0 ± 0.3% continuum polarization in the J band oriented at ~160° on sky 182 days after the explosion. Assuming a prolate geometry as in SN 1987A, we infer an ejecta axis ratio of <0.48 with the axis of symmetry pointing at a 70° position angle. The axis ratio is similar to that of SN 1987A, suggesting that the two CCSNe may share intrinsic geometry and inclination angles. Our data do not rule out oblate ejecta. We also observe one other CCSN and two thermonuclear supernovae in the J band. Supernova 2020oi, a stripped-envelope type Ic SN in Messier 100 has broadband p = 0.37 ± 0.09% at peak light, indicative of either a 10% asymmetry or host interstellar polarization. The type Ia SNe 2019ein and 2020ue have <0.33% and <1.08% polarization near peak light, indicative of asymmetries of less than 10% and 20%, respectively.

The observation of a compact object with a mass of 2.50-2.67M_{99.9%. Even if we weaken previously employed constraints on the maximum mass of neutron stars, the probability of a binary black hole origin is still similar to 81%. Furthermore, we study the impact that this observation has on our understanding of the nuclear equation of state by analyzing the allowed region in the mass-radius diagram of neutron stars for both a binary black hole or neutron star-black hole scenario. We find that the unlikely scenario in which the secondary object was a neutron star requires rather stiff equations of state with a maximum speed of sound c(x) >= root 0.6 times the speed of light, while the binary black hole scenario does not offer any new insight.}

Neutron star mergers eject neutron-rich matter in which heavy elements are synthesized. The decay of these freshly synthesized elements powers electromagnetic transients ('macronovae' or 'kilonovae') whose luminosity and colour strongly depend on their nuclear composition. If the ejecta are very neutron-rich (electron fraction Y-e < 0.25), they contain fair amounts of lanthanides and actinides that have large opacities and therefore efficiently trap the radiation inside the ejecta so that the emission peaks in the red part of the spectrum. Even small amounts of this high-opacity material can obscure emission from lower lying material and therefore act as a 'lanthanide curtain'. Here, we investigate how a relativistic jet that punches through the ejecta can potentially push away a significant fraction of the high opacity material before the macronova begins to shine. We use the results of detailed neutrino-driven wind studies as initial conditions and explore with 3D special relativistic hydrodynamic simulations how jets are propagating through these winds. Subsequently, we perform Monte Carlo radiative transfer calculations to explore the resulting macronova emission. We find that the hole punched by the jet makes the macronova brighter and bluer for on-axis observers during the first few days of emission, and that more powerful jets have larger impacts on the macronova.

We predict linear polarization for a radioactively powered kilonova following the merger of a black hole and a neutron star. Specifically, we perform 3D Monte Carlo radiative transfer simulations for two different models, both featuring a lanthanide-rich dynamical ejecta component from numerical-relativity simulations while only one including an additional lanthanide-free disc-wind component. We calculate polarization spectra for nine different orientations at 1.5, 2.5, and 3.5 d after the merger and in the 0.1-2 mu m wavelength range. We find that both models are polarized at a detectable level 1.5 d after the merger while show negligible levels thereafter. The polarization spectra of the two models are significantly different. The model lacking a disc wind shows no polarization in the optical, while a signal increasing at longer wavelengths and reaching similar to 1-6 per cent at 2 mu m depending on the orientation. The model with a disc-wind component, instead, features a characteristic 'double-peak' polarization spectrum with one peak in the optical and the other in the infrared. Polarimetric observations of future events will shed light on the debated neutron richness of the disc-wind component. The detection of optical polarization would unambiguously reveal the presence of a lanthanide-free disc-wind component, while polarization increasing from zero in the optical to a peak in the infrared would suggest a lanthanide-rich composition for the whole ejecta. Future polarimetric campaigns should prioritize observations in the first similar to 48 h and in the 0.5-2 mu m range, where polarization is strongest, but also explore shorter wavelengths/later times where no signal is expected from the kilonova and the interstellar polarization can be safely estimated.

Time domain astronomy was revolutionized with the discovery of the first kilonova, AT2017gfo, in August 2017, which was associated with the gravitational wave signal GW170817. Since this event, numerous wide-field surveys have been optimizing search strategies to maximize their efficiency of detecting these fast and faint transients. With the Panoramic Survey Telescope and Rapid Response System (Pan-STARRS), we have been conducting a volume-limited survey for intrinsically faint and fast-fading events to a distance of D similar or equal to 200 Mpc. Two promising candidates have been identified from this archival search, with sparse data - PS15cey and PS17cke. Here, we present more detailed analysis and discussion of their nature. We observe that PS15cey was a luminous, fast-declining transient at 320 Mpc. Models of BH-NS mergers with a very stiff equation of state could possibly reproduce the luminosity and decline but the physical parameters are extreme. A more likely scenario is that this was an AT2018kzr-like merger event. PS17cke was a faint and fast-declining event at 15 Mpc. We explore several explosion scenarios of this transient including models of it as a NS-NS and BH-NS merger, the outburst of a massive luminous star, and compare it against other known fast-fading transients. Although there is uncertainty in the explosion scenario due to difficulty in measuring the explosion epoch, we find PS17cke to be a plausible kilonova candidate from the model comparisons.

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 lambda 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.

Context. Gravitational wave (GW) astronomy has rapidly reached maturity, becoming a fundamental observing window for modern astrophysics. The coalescences of a few tens of black hole (BH) binaries have been detected, while the number of events possibly including a neutron star (NS) is still limited to a few. On 2019 August 14, the LIGO and Virgo interferometers detected a high-significance event labelled S190814bv. A preliminary analysis of the GW data suggests that the event was likely due to the merger of a compact binary system formed by a BH and a NS. Aims. In this paper, we present our extensive search campaign aimed at uncovering the potential optical and near infrared electromagnetic counterpart of S190814bv. We found no convincing electromagnetic counterpart in our data. We therefore use our non-detection to place limits on the properties of the putative outflows that could have been produced by the binary during and after the merger. Methods. Thanks to the three-detector observation of S190814bv, and given the characteristics of the signal, the LIGO and Virgo Collaborations delivered a relatively narrow localisation in low latency - a 50% (90%) credible area of 5 deg(2) (23 deg(2)) - despite the relatively large distance of 26752 Mpc. ElectromagNetic counterparts of GRAvitational wave sources at the VEry Large Telescope collaboration members carried out an intensive multi-epoch, multi-instrument observational campaign to identify the possible optical and near infrared counterpart of the event. In addition, the ATLAS, GOTO, GRAWITA-VST, Pan-STARRS, and VINROUGE projects also carried out a search on this event. In this paper, we describe the combined observational campaign of these groups. Results. Our observations allow us to place limits on the presence of any counterpart and discuss the implications for the kilonova (KN), which was possibly generated by this NS-BH merger, and for the strategy of future searches. The typical depth of our wide-field observat

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 R < 1775 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).

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 similar to -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 similar to-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 Ni-56 mixed in the outermost regions of the ejecta and with double-detonation models having thin helium layers (M-He = 0.01 M-circle dot) 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-timeg - rslopes 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.

Joint multimessenger observations with gravitational waves and electromagnetic (EM) data offer new insights into the astrophysical studies of compact objects. The third Advanced LIGO and Advanced Virgo observing run began on 2019 April 1; during the 11 months of observation, there have been 14 compact binary systems candidates for which at least one component is potentially a neutron star. Although intensive follow-up campaigns involving tens of ground and space-based observatories searched for counterparts, no EM counterpart has been detected. Following on a previous study of the first six months of the campaign, we present in this paper the next five months of the campaign from 2019 October to 2020 March. We highlight two neutron star-black hole candidates (S191205ah and S200105ae), two binary neutron star candidates (S191213g and S200213t), and a binary merger with a possible neutron star and a `MassGap' component, S200115j. Assuming that the gravitational-wave (GW) candidates are of astrophysical origin and their location was covered by optical telescopes, we derive possible constraints on the matter ejected during the events based on the non-detection of counterparts. We find that the follow-up observations during the second half of the third observing run did not meet the necessary sensitivity to constrain the source properties of the potential GW candidate. Consequently, we suggest that different strategies have to be used to allow a better usage of the available telescope time. We examine different choices for follow-up surveys to optimize sky localization coverage versus observational depth to understand the likelihood of counterpart detection.

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 similar to 0.2%, which indicates an almost spherical photosphere, with an axial ratio similar to 1.05. The intrinsic polarization at the late phase increases to similar to 0.8%, which corresponds to the photosphere with an axial ratio similar to 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.

Growing evidence suggests that Type Iax supernovae might be the result of thermonuclear deflagrations of Chandrasekhar-mass white dwarfs in binary systems. We carry out Monte Carlo radiative transfer simulations and predict spectropolarimetric features originating from the supernova explosion and subsequent ejecta interaction with the companion star. Specifically, we calculate viewing-angle dependent flux and polarisation spectra for a 3D model simulating the deflagration of a Chandrasekhar-mass white dwarf and, for a second model, simulating the ejecta interaction with a main-sequence star. We find that the intrinsic signal is weakly polarised and only mildly viewing-angle dependent, owing to the overall spherical symmetry of the explosion and the depolarising contribution of iron-group elements dominating the ejecta composition. The interaction with the companion star carves out a cavity in the ejecta and produces a detectable, but modest signal that is significant only at relatively blue wavelengths (less than or similar to 5000 angstrom). In particular, increasingly fainter and redder spectra are predicted for observer orientations further from the cavity, while a modest polarisation signal P similar to 0.2 per cent is found at blue wavelengths for orientations 30 degrees and 45 degrees away from the cavity. We find a reasonable agreement between the interaction model viewed from these orientations and spectropolarimetric data of SN 2005hk and interpret the maximum-light polarisation signal seen at blue wavelengths for this event as a possible signature of the ejecta-companion interaction. We encourage further polarimetric observations of SNe Iax to test whether our results can be extended and generalised to the whole SN Iax class.

On 2019 August 14, the Advanced LIGO and Virgo interferometers detected the high-significance gravitational wave (GW) signal S190814bv. The GW data indicated that the event resulted from a neutron star-black hole (NSBH) merger, or potentially a low-mass binary BH merger. Due to the low false-alarm rate and the precise localization (23 deg(2) at 90%), S190814bv presented the community with the best opportunity yet to directly observe an optical/near-infrared counterpart to an NSBH merger. To search for potential counterparts, the GROWTH Collaboration performed real-time image subtraction on six nights of public Dark Energy Camera images acquired in the 3 weeks following the merger, covering >98% of the localization probability. Using a worldwide network of follow-up facilities, we systematically undertook spectroscopy and imaging of optical counterpart candidates. Combining these data with a photometric redshift catalog, we ruled out each candidate as the counterpart to S190814bv and placed deep, uniform limits on the optical emission associated with S190814bv. For the nearest consistent GW distance, radiative transfer simulations of NSBH mergers constrain the ejecta mass of S190814bv to be M-ej < 0.04 M-circle dot at polar viewing angles, or M-ej < 0.03 M-circle dot if the opacity is kappa < 2 cm(2)g(-1). Assuming a tidal deformability for the NS at the high end of the range compatible with GW170817 results, our limits would constrain the BH spin component aligned with the orbital momentum to be chi < 0.7 for mass ratios Q < 6, with weaker constraints for more compact NSs.

We extend the range of validity of the ARTIS 3D radiative transfer code up to hundreds of days after explosion, when Type Ia supernovae (SNe Ia) are in their nebular phase. To achieve this, we add a non-local thermodynamic equilibrium population and ionization solver, a new multifrequency radiation field model, and a new atomic data set with forbidden transitions. We treat collisions with non-thermal leptons resulting from nuclear decays to account for their contribution to excitation, ionization, and heating. We validate our method with a variety of tests including comparing our synthetic nebular spectra for the well-known one-dimensional W7 model with the results of other studies. As an illustrative application of the code, we present synthetic nebular spectra for the detonation of a sub-Chandrasekhar white dwarf (WD) in which the possible effects of gravitational settling of Ne-22 prior to explosion have been explored. Specifically, we compare synthetic nebular spectra for a 1.06 M-circle dot WD model obtained when 5.5 Gyr of very efficient settling is assumed to a similar model without settling. We find that this degree of Ne-22 settling has only amodest effect on the resulting nebular spectra due to increased Ni-58 abundance. Due to the high ionization in sub-Chandrasekhar models, the nebular [Ni II] emission remains negligible, while the [Ni III] line strengths are increased and the overall ionization balance is slightly lowered in the model with Ne-22 settling. In common with previous studies of sub-Chandrasekhar models at nebular epochs, these models overproduce [Fe III] emission relative to [Fe II] in comparison to observations of normal SNe Ia.

GW170817 showed that neutron star mergers not only emit gravitational waves but also can release electromagnetic signatures in multiple wavelengths. Within the first half of the third observing run of the Advanced LIGO and Virgo detectors, there have been a number of gravitational wave candidates of compact binary systems for which at least one component is potentially a neutron star. In this article, we look at the candidates S190425z, S190426c, S190510g, S190901ap, and S190910h, predicted to have potentially a non-zero remnant mass, in more detail. All these triggers have been followed up with extensive campaigns by the astronomical community doing electromagnetic searches for their optical counterparts; however, according to the released classification, there is a high probability that some of these events might not be of extraterrestrial origin. Assuming that the triggers are caused by a compact binary coalescence and that the individual source locations have been covered during the EM follow-up campaigns, we employ three different kilonova models and apply them to derive possible constraints on the matter ejection consistent with the publicly available gravitational-wave trigger information and the lack of a kilonova detection. These upper bounds on the ejecta mass can be related to limits on the maximum mass of the binary neutron star candidate S190425z and to constraints on the mass-ratio, spin, and NS compactness for the potential black hole-neutron star candidate S190426c. Our results show that deeper electromagnetic observations for future gravitational wave events near the horizon limit of the advanced detectors are essential.

The detection of GW170817 is revolutionizing many areas of astrophysics with the joint observation of gravitational waves and electromagnetic emissions. These multimessenger events provide a new approach to determine the Hubble constant, thus, they are a promising candidate for mitigating the tension between measurements of type-Ia supernovae via the local distance ladder and the cosmic microwave background. In addition to the “standard siren” provided by the gravitational-wave measurement, the kilonova itself has characteristics that allow one to improve existing measurements or to perform yet another, independent measurement of the Hubble constant without gravitational-wave information. Here, we employ standardization techniques borrowed from the type-Ia community and apply them to kilonovae, not using any information from the gravitational-wave signal. We use two versions of this technique, one derived from direct observables measured from the light curve, and the other based on inferred ejecta parameters, e.g., mass, velocity, and composition, for two different models. These lead to Hubble constant measurements of H0=109+49−35 km s−1 Mpc−1 for the measured analysis, and H0=85+22−17 km s−1 Mpc−1 and H0=79+23−15 km s−1 Mpc−1 for the inferred analyses. This measurement has error bars within ∼2 to the gravitational-wave measurements (H0=74+16−8 km s−1 Mpc−1), showing its promise as an independent constraint on H0.

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 deg(2), 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