Search

Eruptive mass loss of massive stars prior to supernova (SN) explosion is key to understanding their evolution and end fate. An observational signature of pre-SN mass loss is the detection of an early, short-lived peak prior to the radioactive-powered peak in the lightcurve of the SN. This is usually attributed to the SN shock passing through an extended envelope or circumstellar medium (CSM). Such an early peak is common for double-peaked Type IIb SNe with an extended Hydrogen envelope but is uncommon for normal Type Ibc SNe with very compact progenitors. In this paper, we systematically study a sample of 14 double-peaked Type Ibc SNe out of 475 Type Ibc SNe detected by the Zwicky Transient Facility. The rate of these events is ~ 3-9 % of Type Ibc SNe. A strong correlation is seen between the peak brightness of the first and the second peak. We perform a holistic analysis of this sample's photometric and spectroscopic properties. We find that six SNe have ejecta mass less than 1.5 Msun. Based on the nebular spectra and lightcurve properties, we estimate that the progenitor masses for these are less than ~ 12 Msun. The rest have an ejecta mass > 2.4 Msun and a higher progenitor mass. This sample suggests that the SNe with low progenitor masses undergo late-time binary mass transfer. Meanwhile, the SNe with higher progenitor masses are consistent with wave-driven mass loss or pulsation-pair instability-driven mass loss simulations., Submitted to ApJ. References added. Comments are welcome

Abridged - Stars with ZAMS masses between 140 and $260 M_\odot$ are thought to explode as pair-instability supernovae (PISNe). During their thermonuclear runaway, PISNe can produce up to several tens of solar masses of radioactive nickel, resulting in luminous transients similar to some superluminous supernovae (SLSNe). Yet, no unambiguous PISN has been discovered so far. SN2018ibb is a H-poor SLSN at $z=0.166$ that evolves extremely slowly compared to the hundreds of known SLSNe. Between mid 2018 and early 2022, we monitored its photometric and spectroscopic evolution from the UV to the NIR with 2-10m class telescopes. SN2018ibb radiated $>3\times10^{51} \rm erg$ during its evolution, and its bolometric light curve reached $>2\times10^{44} \rm erg\,s^{-1}$ at peak. The long-lasting rise of $>93$ rest-frame days implies a long diffusion time, which requires a very high total ejected mass. The PISN mechanism naturally provides both the energy source ($^{56}$Ni) and the long diffusion time. Theoretical models of PISNe make clear predictions for their photometric and spectroscopic properties. SN2018ibb complies with most tests on the light curves, nebular spectra and host galaxy, potentially all tests with the interpretation we propose. Both the light curve and the spectra require 25-44 $M_\odot$ of freshly nucleosynthesised $^{56}$Ni, pointing to the explosion of a metal-poor star with a He-core mass of 120-130 $M_\odot$ at the time of death. This interpretation is also supported by the tentative detection of [Co II]$λ$1.025$μ$m, which has never been observed in any other PISN candidate or SLSN before. Powering by a central engine, such as a magnetar or a black hole, can be excluded with high confidence. This makes SN2018ibb by far the best candidate for being a PISN, to date., Submitted to A&A, 44 pages, main text 38 pages, 37 figures, 16 Tables

We conduct a systematic tidal disruption event (TDE) demographics analysis using the largest sample of optically selected TDEs. A flux-limited, spectroscopically complete sample of 33 TDEs is constructed using the Zwicky Transient Facility (ZTF) over three years (from October 2018 to September 2021). We infer the black hole (BH) mass ($M_{\rm BH}$) with host galaxy scaling relations, showing that the sample $M_{\rm BH}$ ranges from $10^{5.1}\,M_\odot$ to $10^{8.2}\,M_\odot$. We developed a survey efficiency corrected maximum volume method to infer the rates. The rest-frame $g$-band luminosity function (LF) can be well described by a broken power-law of $\phi (L_g)\propto [(L_g / L_{\rm bk})^{0.3} + (L_g / L_{\rm bk})^{2.6}]^{-1}$, with $L_{\rm bk}=10^{43.1}\,{\rm erg\,s^{-1}}$. In the BH mass regime of $10^{5.2}\lesssim (M_{\rm BH}/M_\odot) \lesssim 10^{7.2}$, the TDE mass function follows $\phi(M_{\rm BH})\propto M_{\rm BH}^{-0.25}$, which favors a flat local BH mass function ($dn_{\rm BH}/d{\rm log}M_{\rm BH}\approx{\rm constant}$). We confirm the significant rate suppression at the high-mass end ($M_{\rm BH}\gtrsim 10^{7.5}\,M_\odot$), which is consistent with theoretical predictions considering direct capture of hydrogen-burning stars by the event horizon. At a host galaxy mass of $M_{\rm gal}\sim 10^{10}\,M_\odot$, the average optical TDE rate is $\approx 3.2\times 10^{-5}\,{\rm galaxy^{-1}\,yr^{-1}}$. We constrain the optical TDE rate to be [3.7, 7.4, and 1.6$]\times 10^{-5}\,{\rm galaxy^{-1}\,yr^{-1}}$ in galaxies with red, green, and blue colors., Comment: 37 pages, 24 figures, submitted to ApJ

One of the open questions following the discovery of GW170817 is whether neutron star mergers are the only astrophysical sites capable of producing r-process elements. Simulations have shown that 0.01-0.1M$_\odot$ of r-process material could be generated in the outflows originating from the accretion disk surrounding the rapidly rotating black hole that forms as a remnant to both neutron star mergers and collapsing massive stars associated with long-duration gamma-ray bursts (collapsars). The hallmark signature of r-process nucleosynthesis in the binary neutron star merger GW170817 was its long-lasting near-infrared emission, thus motivating a systematic photometric study of the light curves of broadlined stripped-envelope (Ic-BL) supernovae (SNe) associated with collapsars. We present the first systematic study of 25 SNe Ic-BL -- discovered with the Zwicky Transient Facility and from the literature -- in the optical/near-infrared bands to determine what quantity of r-process material, if any, is synthesized in these explosions. Using semi-analytic models designed to account for r-process production in SNe Ic-BL, we perform light curve fitting to derive constraints on the r-process mass for these SNe. We also perform independent light curve fits to models without r-process. We find that the r-process-free models are a better fit to the light curves of the objects in our sample. Thus we conclude that there is no compelling evidence of r-process enrichment in any of our objects. Further high-cadence infrared photometric studies and nebular spectroscopic analysis would be sensitive to smaller quantities of r-process ejecta mass or indicate whether all collapsars are completely devoid of r-process nucleosynthesis., 35 pages, 8 figures, 4 tables. Submitted for publication in The Astrophysical Journal

SN 1987A was an unusual hydrogen-rich core-collapse supernova originating from a blue supergiant star. Similar blue supergiant explosions remain a small family of events, and are broadly characterized by their long rises to peak. The Zwicky Transient Facility (ZTF) Census of the Local Universe (CLU) experiment aims to construct a spectroscopically complete sample of transients occurring in galaxies from the CLU galaxy catalog. We identify 13 long-rising (>40 days) Type II supernovae from the volume-limited CLU experiment during a 3.5 year period from June 2018 to December 2021, approximately doubling the previously known number of these events. We present photometric and spectroscopic data of these 13 events, finding peak r-band absolute magnitudes ranging from -15.6 to -17.5 mag and the tentative detection of Ba II lines in 9 events. Using our CLU sample of events, we derive a long-rising Type II supernova rate of $1.37^{+0.26}_{-0.30}\times10^{-6}$ Mpc$^{-3}$ yr$^{-1}$, $\approx$1.4% of the total core-collapse supernova rate. This is the first volumetric rate of these events estimated from a large, systematic, volume-limited experiment., 32 pages, 17 figures, 5 tables. Submitted to ApJ

Full list of authors: Andreoni, Igor; Coughlin, Michael W.; Perley, Daniel A.; Yao, Yuhan; Lu, Wenbin; Cenko, S. Bradley; Kumar, Harsh; Anand, Shreya; Ho, Anna Y. Q.; Kasliwal, Mansi M.; de Ugarte Postigo, Antonio; Sagués-Carracedo, Ana; Schulze, Steve; Kann, D. Alexander; Kulkarni, S. R.; Sollerman, Jesper; Tanvir, Nial; Rest, Armin; Izzo, Luca; Somalwar, Jean J.; Kaplan, David L.; Ahumada, Tomás; Anupama, G. C.; Auchettl, Katie; Barway, Sudhanshu; Bellm, Eric C.; Bhalerao, Varun; Bloom, Joshua S.; Bremer, Michael; Bulla, Mattia; Burns, Eric; Campana, Sergio; Chandra, Poonam; Charalampopoulos, Panos; Cooke, Jeff; D'Elia, Valerio; Das, Kaustav Kashyap; Dobie, Dougal; Agüí Fernández, José Feliciano; Freeburn, James; Fremling, Cristoffer; Gezari, Suvi; Goode, Simon; Graham, Matthew J. ; Hammerstein, Erica; Karambelkar, Viraj R.; Kilpatrick, Charles D.; Kool, Erik C.; Krips, Melanie; Laher, Russ R.; Leloudas, Giorgos; Levan, Andrew; Lundquist, Michael J.; Mahabal, Ashish A.; Medford, Michael S.; Miller, M. Coleman; Möller, Anais; Mooley, Kunal P.; Nayana, A. J.; Nir, Guy; Pang, Peter T. H.; Paraskeva, Emmy; Perley, Richard A.; Petitpas, Glen; Pursiainen, Miika; Ravi, Vikram; Ridden-Harper, Ryan; Riddle, Reed; Rigault, Mickael; Rodriguez, Antonio C.; Rusholme, Ben; Sharma, Yashvi; Smith, I. A.; Stein, Robert D.; Thöne, Christina; Tohuvavohu, Aaron; Valdes, Frank; van Roestel, Jan; Vergani, Susanna D.; Wang, Qinan; Zhang, Jielai.--Correction to: Nature https://doi.org/10.1038/s41586-022-05465-8 Published online 30 November 2022., Tidal disruption events (TDEs) are bursts of electromagnetic energy that are released when supermassive black holes at the centres of galaxies violently disrupt a star that passes too close1. TDEs provide a window through which to study accretion onto supermassive black holes; in some rare cases, this accretion leads to launching of a relativistic jet2,3,4,5,6,7,8,9, but the necessary conditions are not fully understood. The best-studied jetted TDE so far is Swift J1644+57, which was discovered in γ-rays, but was too obscured by dust to be seen at optical wavelengths. Here we report the optical detection of AT2022cmc, a rapidly fading source at cosmological distance (redshift z = 1.19325) the unique light curve of which transitioned into a luminous plateau within days. Observations of a bright counterpart at other wavelengths, including X-ray, submillimetre and radio, supports the interpretation of AT2022cmc as a jetted TDE containing a synchrotron ‘afterglow’, probably launched by a supermassive black hole with spin greater than approximately 0.3. Using four years of Zwicky Transient Facility10 survey data, we calculate a rate of 0.02+0.04−0.01 Gpc−3 yr−1 for on-axis jetted TDEs on the basis of the luminous, fast-fading red component, thus providing a measurement complementary to the rates derived from X-ray and radio observations11. Correcting for the beaming angle effects, this rate confirms that approximately 1 per cent of TDEs have relativistic jets. Optical surveys can use AT2022cmc as a prototype to unveil a population of jetted TDEs. © 2022. The Author(s), under exclusive licence to Springer Nature Limited., M.W.C. acknowledges support from the National Science Foundation with grant numbers PHY-2010970 and OAC-2117997. E.C.K. acknowledges support from the G.R.E.A.T .research environment and the Wenner-Gren Foundations. M. Bulla acknowledges support from the Swedish Research Council (reg. no. 2020-03330). H.K. and T.A. thank the LSSTC Data Science Fellowship Program, which is funded by LSSTC, NSF Cybertraining Grant no. 1829740, the Brinson Foundation, and the Moore Foundation; their participation in the programme has benefited this work. W.L. was supported by the Lyman Spitzer, Jr. Fellowship at Princeton University. M.R. has received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (grant agreement no. 759194 - USNAC). G.L., P. Charalampopoulos and M.P. were supported by a research grant (19054) from VILLUM FONDEN. P.T.H.P. is supported by the research programme of the Netherlands Organization for Scientific Research (NWO). D.A.K. acknowledges support from Spanish National Research Project RTI2018-098104-J-I00 (GRBPhot). The material is based on work supported by NASA under award no. 80GSFC17M0002. A.J.N. acknowledges DST-INSPIRE Faculty Fellowship (IFA20-PH-259) for supporting this research. I.A. is a Neil Gehrels Fellow. This work has been supported by the research project grant ‘Understanding the Dynamic Universe’ funded by the Knut and Alice Wallenberg Foundation under Dnr KAW 2018.0067. Based on observations obtained with the Samuel Oschin Telescope 48-inch and the 60-inch Telescope at the Palomar Observatory as part of the Zwicky Transient Facility (ZTF) project. ZTF is supported by the National Science Foundation under grant no. AST-1440341 and a collaboration including Caltech, IPAC, the Weizmann Institute for Science, the Oskar Klein Center at Stockholm University, the University of Maryland, the University of Washington (UW), Deutsches Elektronen-Synchrotron and Humboldt University, Los Alamos National Laboratories, the TANGO Consortium of Taiwan, the University of Wisconsin at Milwaukee, and Lawrence Berkeley National Laboratories. Operations are conducted by Caltech Optical Observatories, IPAC, and UW. The work is partly based on the observations made with the Gran Telescopio Canarias (GTC), installed in the Spanish Observatorio del Roque de los Muchachos of the Instituto de Astrofísica de Canarias, on the island of La Palma. We acknowledge all co-investigators of our GTC proposal. SED Machine is based on work supported by the National Science Foundation under grant no. 1106171. The ZTF forced-photometry service was funded under the Heising–Simons Foundation grant no. 12540303 (PI: M.J.G.). The research leading to these results has received funding from the European Union Seventh Framework Programme (FP7/2013-2016) under grant agreement no. 312430 (OPTICON). Based on observations collected at the European Southern Observatory under ESO programme 106.21T6.015. This work made use of data from the GROWTH-India Telescope (GIT) set up by the Indian Institute of Astrophysics (IIA) and the Indian Institute of Technology Bombay (IITB). It is located at the Indian Astronomical Observatory (Hanle), operated by IIA. We acknowledge funding by the IITB alumni batch of 1994, which partially supports operations of the telescope. Telescope technical details are available at https://sites.google.com/view/growthindia. Based on observations made with the Nordic Optical Telescope, owned in collaboration by the University of Turku and Aarhus University, and operated jointly by Aarhus University, the University of Turku and the University of Oslo, representing Denmark, Finland and Norway, the University of Iceland and Stockholm University at the Observatorio del Roque de los Muchachos, La Palma, Spain, of the Instituto de Astrofísica de Canarias. The Liverpool Telescope is operated on the island of La Palma by Liverpool John Moores University in the Spanish Observatorio del Roque de los Muchachos of the Instituto de Astrofísica de Canarias with financial support from the UK Science and Technology Facilities Council. The National Radio Astronomy Observatory is a facility of the National Science Foundation operated under cooperative agreement by Associated Universities, Inc. Partly based on observations collected at Centro Astronómico Hispano en Andalucía (CAHA) at Calar Alto, operated jointly by Instituto de Astrofísica de Andalucía (CSIC) and Junta de Andalucía. The James Clerk Maxwell Telescope is operated by the East Asian Observatory on behalf of The National Astronomical Observatory of Japan, Academia Sinica Institute of Astronomy and Astrophysics, the Korea Astronomy and Space Science Institute, the National Astronomical Research Institute of Thailand and the Center for Astronomical Mega-Science (as well as the National Key R&D Program of China with no. 2017YFA0402700). Additional funding support is provided by the Science and Technology Facilities Council of the UK and participating universities and organizations in the UK and Canada. Additional funds for the construction of SCUBA-2 were provided by the Canada Foundation for Innovation. The JCMT data reported here were obtained under project M22AP030 (principal investigator D.A.P.). The Starlink software is currently supported by the East Asian Observatory. The Submillimeter Array is a joint project between the Smithsonian Astrophysical Observatory and the Academia Sinica Institute of Astronomy and Astrophysics and is funded by the Smithsonian Institution and the Academia Sinica. This work is based on observations carried out under project number W21BK with the IRAM NOEMA Interferometer. IRAM is supported by INSU/CNRS (France), MPG (Germany) and IGN (Spain). Based on observations obtained at the international Gemini Observatory, a programme of NSF’s NOIRLab, which is managed by the Association of Universities for Research in Astronomy (AURA) under a cooperative agreement with the National Science Foundation. On behalf of the Gemini Observatory partnership: the National Science Foundation (USA), National Research Council (Canada), Agencia Nacional de Investigación y Desarrollo (Chile), Ministerio de Ciencia, Tecnología e Innovación (Argentina), Ministério da Ciência, Tecnologia, Inovações e Comunicações (Brazil), and Korea Astronomy and Space Science Institute (Republic of Korea). This work was enabled by observations made from the Gemini North telescope, located within the Maunakea Science Reserve and adjacent to the summit of Maunakea., With funding from the Spanish government through the Severo Ochoa Centre of Excellence accreditation SEV-2017-0709.

The dividing line between gamma-ray bursts (GRBs) and ordinary stripped-envelope core-collapse supernovae (SNe) is yet to be fully understood. Observationally mapping the variety of ejecta outcomes (ultra-relativistic, mildly-relativistic or non-relativistic) in SNe of Type Ic with broad lines (Ic-BL) can provide a key test to stellar explosion models. However, this requires large samples of the rare Ic-BL events with follow-up observations in the radio, where fast ejecta can be probed largely free of geometry and viewing angle effects. Here, we present the results of a radio (and X-ray) follow-up campaign of 16 SNe Ic-BL detected by the Zwicky Transient Facility (ZTF). Our radio campaign resulted in 4 counterpart detections and 12 deep upper limits. None of the events in our sample is as relativistic as SN 1998bw and we constrain the fraction of SN 1998bw-like explosions to $< 19\%$ (3$\sigma$ Gaussian equivalent), a factor of $\approx 2$ smaller than previously established. We exclude relativistic ejecta with radio luminosity densities in between $\approx 5\times10^{27}$ erg s$^{-1}$ Hz$^{-1}$ and $\approx 10^{29}$ erg s$^{-1}$ Hz$^{-1}$ at $t\gtrsim 20$ d since explosion for $\approx 60\%$ of the events in our sample. This shows that SNe Ic-BL similar to the GRB-associated SN 1998bw, SN 2003lw, SN 2010dh, or to the relativistic SN 2009bb and iPTF17cw, are rare. Our results also exclude an association of the SNe Ic-BL in our sample with largely off-axis GRBs with energies $E\gtrsim 10^{50}$ erg. The parameter space of SN2006aj-like events (faint and fast-peaking radio emission) is, on the other hand, left largely unconstrained and systematically exploring it represents a promising line of future research., Comment: 28 pages, 11 figures, submitted to AAS journal

The fate of stars in the zero-age main-sequence (ZAMS) range $\approx 8-12$ Msun is unclear. They could evolve to form white dwarfs or explode as electron-capture supernovae (SNe) or iron core-collapse SNe (CCSNe). Even though the initial mass function indicates that this mass range should account for over 40% of all CCSNe progenitors, few have been observationally confirmed, likely owing to the faintness and rapid evolution of these transients. In this paper, we present a sample of nine Ca-rich/O-poor Type IIb SNe detected by the Zwicky Transient Facility with progenitors likely in this mass range. We perform a holistic analysis of the spectroscopic and photometric properties of the sample. These sources have a flux ratio of [Ca II] $\lambda \lambda$7291, 7324 to [O I] $\lambda \lambda$6300, 6364 of $\gtrsim$ 2 in their nebular spectra. Comparing the measured [O I] luminosity ($\lesssim 10^{39} \mathrm{erg\ s^{-1}}$) and derived oxygen mass ($\lesssim 0.1$ Msun) with theoretical models, we infer that the progenitor ZAMS mass for these explosions is less than 12 Msun. These correspond to He-stars with core masses less than around 3 Msun. We find that the ejecta properties (Mej $\lesssim 1$ Msun) are also consistent with those expected for such low mass He-stars. The low ejecta mass of these sources indicates a class of strongly-stripped SNe that is a transition between the regular stripped-envelope SNe and ultra-stripped SNe. The progenitor could be stripped by a main sequence companion and result in the formation of a neutron star $-$ main sequence binary. Such binaries have been suggested to be progenitors of neutron star $-$ white dwarf systems that could merge within a Hubble time, and be detectable with LISA., Comment: 21 pages, 13 Figures. Submitted to ApJ

We present results from extensive broadband follow-up of GRB 210204A over the period of thirty days. We detect optical flares in the afterglow at 7.6 x 10^5 s and 1.1 x 10^6 s after the burst: the most delayed flaring ever detected in a GRB afterglow. At the source redshift of 0.876, the rest-frame delay is 5.8 x 10^5 s (6.71 d). We investigate possible causes for this flaring and conclude that the most likely cause is a refreshed shock in the jet. The prompt emission of the GRB is within the range of typical long bursts: it shows three disjoint emission episodes, which all follow the typical GRB correlations. This suggests that GRB 210204A might not have any special properties that caused late-time flaring, and the lack of such detections for other afterglows might be resulting from the paucity of late-time observations. Systematic late-time follow-up of a larger sample of GRBs can shed more light on such afterglow behaviour. Further analysis of the GRB 210204A shows that the late time bump in the light curve is highly unlikely due to underlying SNe at redshift (z) = 0.876 and is more likely due to the late time flaring activity. The cause of this variability is not clearly quantifiable due to the lack of multi-band data at late time constraints by the bad weather conditions. The flare of GRB 210204A is the latest flare detected to date., 18 pages, 12 figures

We present a search for extragalactic fast blue optical transients (FBOTs) during Phase I of the Zwicky Transient Facility (ZTF). We identify 38 candidates with durations above half-maximum light 1 d < t1/2 < 12 d, of which 28 have blue (g-r, Replaced following peer-review process. 46 pages, 20 figures. Accepted for publication in ApJ

UVEX is a proposed medium class Explorer mission designed to provide crucial missing capabilities that will address objectives central to a broad range of modern astrophysics. The UVEX design has two co-aligned wide-field imagers operating in the FUV and NUV and a powerful broadband medium resolution spectrometer. In its two-year baseline mission, UVEX will perform a multi-cadence synoptic all-sky survey 50/100 times deeper than GALEX in the NUV/FUV, cadenced surveys of the Large and Small Magellanic Clouds, rapid target of opportunity followup, as well as spectroscopic followup of samples of stars and galaxies. The science program is built around three pillars. First, UVEX will explore the low-mass, low-metallicity galaxy frontier through imaging and spectroscopic surveys that will probe key aspects of the evolution of galaxies by understanding how star formation and stellar evolution at low metallicities affect the growth and evolution of low-metallicity, low-mass galaxies in the local universe. Such galaxies contain half the mass in the local universe, and are analogs for the first galaxies, but observed at distances that make them accessible to detailed study. Second, UVEX will explore the dynamic universe through time-domain surveys and prompt spectroscopic followup capability will probe the environments, energetics, and emission processes in the early aftermaths of gravitational wave-discovered compact object mergers, discover hot, fast UV transients, and diagnose the early stages of stellar explosions. Finally, UVEX will become a key community resource by leaving a large all-sky legacy data set, enabling a wide range of scientific studies and filling a gap in the new generation of wide-field, sensitive optical and infrared surveys provided by the Rubin, Euclid, and Roman observatories. This paper discusses the scientific potential of UVEX, and the broad scientific program., Comment: 69 pages, 43 figures

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. 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 of 4480 deg^2, median distance of 267 Mpc and false alarm rates ranging from 1.5 to 1e-25 per yr. The ZTF coverage had a median enclosed probability of 39%, median depth of 20.8mag, and median response time of 1.5 hr. The O3 follow-up by the GROWTH team comprised 340 UVOIR photometric points, 64 OIR spectra, and 3 radio. We find no promising kilonova (radioactivity-powered counterpart) and we convert the upper limits to constrain the underlying kilonova luminosity function. Assuming that all kilonovae are at least as luminous as GW170817 at discovery (-16.1mag), we calculate our joint probability of detecting zero kilonovae is only 4.2%. If we assume that all kilonovae are brighter than -16.6mag (extrapolated peak magnitude of GW170817) and fade at 1 mag/day (similar to GW170817), the joint probability of zero detections is 7%. If we separate the NSBH and BNS populations, the joint probability of zero detections, assuming all kilonovae are brighter than -16.6mag, is 9.7% for NSBH and 7.9% for BNS mergers. Moreover, 1e-4 or phi>30deg to be consistent with our limits. (Abridged), Submitted to ApJ