Using the Zwicky Transient Facility alert stream, we are conducting a large spectroscopic campaign to construct a complete, volume-limited sample of transients brighter than 20 mag, and coincident ...within 100″ of galaxies in the Census of the Local Universe catalog. We describe the experiment design and spectroscopic completeness from the first 16 months of operations, which have classified 754 supernovae. We present results from a systematic search for calcium-rich gap transients in the sample of 22 low-luminosity (peak absolute magnitude M > −17), hydrogen-poor events found in the experiment. We report the detection of eight new events, and constrain their volumetric rate to 15% 5% of the SN Ia rate. Combining this sample with 10 previously known events, we find a likely continuum of spectroscopic properties ranging from events with SN Ia-like features (Ca-Ia objects) to those with SN Ib/c-like features (Ca-Ib/c objects) at peak light. Within the Ca-Ib/c events, we find two populations distinguished by their red (g − r 1.5 mag) or green ( mag) colors at the r-band peak, wherein redder events show strong line blanketing features and slower light curves (similar to Ca-Ia objects), weaker He lines, and lower Ca ii/O i in the nebular phase. We find that all together the spectroscopic continuum, volumetric rates, and striking old environments are consistent with the explosive burning of He shells on low-mass white dwarfs. We suggest that Ca-Ia and red Ca-Ib/c objects arise from the double detonation of He shells, while green Ca-Ib/c objects are consistent with low-efficiency burning scenarios like detonations in low-density shells or deflagrations.
The discovery of a transient kilonova following the gravitational-wave (GW) event GW170817 highlighted the critical need for coordinated rapid and wide-field observations, inference, and follow-up ...across the electromagnetic spectrum. In the southern hemisphere, the Dark Energy Camera (DECam) on the Blanco 4 m telescope is well suited to this task, as it is able to cover wide fields quickly while still achieving the depths required to find kilonovae like the one accompanying GW170817 to ∼500 Mpc, the binary neutron star (NS) horizon distance for current generation of LIGO/Virgo collaboration (LVC) interferometers. Here, as part of the multi-facility follow-up by the Global Relay of Observatories Watching Transients Happen collaboration, we describe the observations and automated data movement, data reduction, candidate discovery, and vetting pipeline of our target-of-opportunity DECam observations of S190426c, the first possible NS-black hole merger detected in GWs. Starting 7.5 hr after S190426c, over 11.28 hr of observations, we imaged an area of 525 deg2 (r band) and 437 deg2 (z band); this was 16.3% of the total original localization probability, and nearly all of the probability visible from the southern hemisphere. The machine-learning-based pipeline was optimized for fast turnaround, delivering transients for human vetting within 17 minutes, on average, of shutter closure. We reported nine promising counterpart candidates 2.5 hr before the end of our observations. One hour after our data-taking ended (roughly 20 hr after the announcement of S190426c), LVC released a refined skymap that reduced the probability coverage of our observations to 8.0%, demonstrating a critical need for localization updates on shorter (∼hour) timescales. Our observations yielded no detection of a bona fide counterpart to mz = 21.7 and mr = 22.2 at the 5 level of significance, consistent with the refined LVC positioning. We view these observations and rapid inferencing as an important real-world test for this novel end-to-end wide-field pipeline.
Abstract
AT2019wey (SRGA J043520.9+552226, SRGE J043523.3+552234) is a transient first reported by the ATLAS optical survey in 2019 December. It rose to prominence upon detection, three months later, ...by the Spektrum-Roentgen-Gamma (SRG) mission in its first all-sky survey. X-ray observations reported in Yao et al. suggest that AT2019wey is a Galactic low-mass X-ray binary (LMXB) with a black hole (BH) or neutron star (NS) accretor. Here we present ultraviolet, optical, near-infrared, and radio observations of this object. We show that the companion is a short-period (
P
≲ 16 hr) low-mass (<1
M
⊙
) star. We consider AT2019wey to be a candidate BH system since its locations on the
L
radio
–
L
X
and
L
opt
–
L
X
diagrams are closer to BH binaries than NS binaries. We demonstrate that from 2020 June to August, despite the more than 10 times brightening at radio and X-ray wavelengths, the optical luminosity of AT2019wey only increased by 1.3–1.4 times. We interpret the UV/optical emission before the brightening as thermal emission from a truncated disk in a hot accretion flow and the UV/optical emission after the brightening as reprocessing of the X-ray emission in the outer accretion disk. AT2019wey demonstrates that combining current wide-field optical surveys and SRG provides a way to discover the emerging population of short-period BH LMXB systems with faint X-ray outbursts.
The first two months of the third Advanced LIGO and Virgo observing run (2019 April-May) showed that distant gravitational-wave (GW) events can now be readily detected. Three candidate mergers ...containing neutron stars (NS) were reported in a span of 15 days, all likely located more than 100 Mpc away. However, distant events such as the three new NS mergers are likely to be coarsely localized, which highlights the importance of facilities and scheduling systems that enable deep observations over hundreds to thousands of square degrees to detect the electromagnetic counterparts. On 2019 May 10 02:59:39.292 UT the GW candidate S190510g was discovered and initially classified as a binary neutron star (BNS) merger with 98% probability. The GW event was localized within an area of 3462 deg2, later refined to 1166 deg2 (90%) at a distance of 227 92 Mpc. We triggered Target-of-Opportunity observations with the Dark Energy Camera (DECam), a wide-field optical imager mounted at the prime focus of the 4 m Blanco Telescope at Cerro Tololo Inter-American Observatory in Chile. This Letter describes our DECam observations and our real-time analysis results, focusing in particular on the design and implementation of the observing strategy. Within 24 hr of the merger time, we observed 65% of the total enclosed probability of the final skymap with an observing efficiency of 94%. We identified and publicly announced 13 candidate counterparts. S190510g was reclassified 1.7 days after the merger, after our observations were completed, with a "BNS merger" probability reduced from 98% to 42% in favor of a "terrestrial classification.
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 deg2 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 Mej < 0.04 M at polar viewing angles, or Mej < 0.03 M if the opacity is κ < 2 cm2g−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 χ < 0.7 for mass ratios Q < 6, with weaker constraints for more compact NSs.
Abstract
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 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
M
ej
< 0.03
M
⊙
,
X
lan
> 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%.
Palomar Gattini-IR is a new wide-field, near-infrared (NIR) robotic time domain survey operating at Palomar Observatory. Using a 30 cm telescope mounted with a H2RG detector, Gattini-IR achieves a ...field of view (FOV) of 25 sq. deg. with a pixel scale of 8 7 in J-band. Here, we describe the system design, survey operations, data processing system and on-sky performance of Palomar Gattini-IR. As a part of the nominal survey, Gattini-IR scans 7500 square degrees of the sky every night to a median 5 depth of 15.7 AB mag outside the Galactic plane. The survey covers 15,000 square degrees of the sky visible from Palomar with a median cadence of 2 days. A real-time data processing system produces stacked science images from dithered raw images taken on sky, together with point-spread function (PSF)-fit source catalogs and transient candidates identified from subtractions within a median delay of 4 hr from the time of observation. The calibrated data products achieve an astrometric accuracy (rms) of 0 7 with respect to Gaia DR2 for sources with signal-to-noise ratio > 10, and better than 0 35 for sources brighter than 12 Vega mag. The photometric accuracy (rms) achieved in the PSF-fit source catalogs is better than 3% for sources brighter than 12 Vega mag and fainter than the saturation magnitude of 8.5 Vega mag, as calibrated against the Two Micron All Sky Survey catalog. The detection efficiency of transient candidates injected into the images is better than 90% for sources brighter than the 5 limiting magnitude. The photometric recovery precision of injected sources is 3% for sources brighter than 13 mag, and the astrometric recovery rms is 0 9. Reference images generated by stacking several field visits achieve depths of 16.5 AB mag over 60% of the sky, while it is limited by confusion in the Galactic plane. With a FOV 40× larger than any other existing NIR imaging instrument, Gattini-IR is probing the reddest and dustiest transients in the local universe such as dust obscured supernovae in nearby galaxies, novae behind large columns of extinction within the galaxy, reddened microlensing events in the Galactic plane and variability from cool and dust obscured stars. We present results from transients and variables identified since the start of the commissioning period.
Abstract
The nova rate in the Milky Way remains largely uncertain, despite its vital importance in constraining models of Galactic chemical evolution as well as understanding progenitor channels for ...Type Ia supernovae. The rate has been previously estimated to be in the range of ≈10–300 yr
−1
, either based on extrapolations from a handful of very bright optical novae or the nova rates in nearby galaxies; both methods are subject to debatable assumptions. The total discovery rate of optical novae remains much smaller (≈5–10 yr
−1
) than these estimates, even with the advent of all-sky optical time-domain surveys. Here, we present a systematic sample of 12 spectroscopically confirmed Galactic novae detected in the first 17 months of Palomar Gattini-IR (PGIR), a wide-field near-infrared time-domain survey. Operating in the
J
band (≈1.2
μ
m), which is significantly less affected by dust extinction compared to optical bands, the extinction distribution of the PGIR sample is highly skewed to a large extinction values (>50% of events obscured by
A
V
≳ 5 mag). Using recent estimates for the distribution of Galactic mass and dust, we show that the extinction distribution of the PGIR sample is commensurate with dust models. The PGIR extinction distribution is inconsistent with that reported in previous optical searches (null-hypothesis probability <0.01%), suggesting that a large population of highly obscured novae have been systematically missed in previous optical searches. We perform the first quantitative simulation of a 3
π
time-domain survey to estimate the Galactic nova rate using PGIR, and derive a rate of
≈
43.7
−
8.7
+
19.5
yr
−1
. Our results suggest that all-sky near-infrared time-domain surveys are well poised to uncover the Galactic nova population.
Interaction-powered supernovae (SNe) explode within an optically thick circumstellar medium (CSM) that could be ejected during eruptive events. To identify and characterize such pre-explosion ...outbursts, we produce forced-photometry light curves for 196 interacting SNe, mostly of Type IIn, detected by the Zwicky Transient Facility between early 2018 and 2020 June. Extensive tests demonstrate that we only expect a few false detections among the 70,000 analyzed pre-explosion images after applying quality cuts and bias corrections. We detect precursor eruptions prior to 18 Type IIn SNe and prior to the Type Ibn SN 2019uo. Precursors become brighter and more frequent in the last months before the SN and month-long outbursts brighter than magnitude −13 occur prior to 25% (5-69%, 95% confidence range) of all Type IIn SNe within the final three months before the explosion. With radiative energies of up to 1049 erg, precursors could eject ∼1 M of material. Nevertheless, SNe with detected precursors are not significantly more luminous than other SNe IIn, and the characteristic narrow hydrogen lines in their spectra typically originate from earlier, undetected mass-loss events. The long precursor durations require ongoing energy injection, and they could, for example, be powered by interaction or by a continuum-driven wind. Instabilities during the neon- and oxygen-burning phases are predicted to launch precursors in the final years to months before the explosion; however, the brightest precursor is 100 times more energetic than anticipated.
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.