Showing 1–50 of 220 results for author: Bellm, E

ZTF SN Ia DR2: Study of Type Ia Supernova lightcurve fits

Authors: M. Rigault, M. Smith, N. Regnault, D. W. Kenworthy, K. Maguire, A. Goobar, G. Dimitriadis, M. Amenouche, M. Aubert, C. Barjou-Delayre, C. E. Bellm, U. Burgaz, B. Carreres, Y. Copin, M. Deckers, T. de Jaeger, S. Dhawan, F. Feinstein, D. Fouchez, L. Galbany, M. Ginolin, J. M. Graham, Y. -L. Kim, M. Kowalski, D. Kuhn , et al. (12 additional authors not shown)

Abstract: Type Ia supernova (SN Ia) cosmology relies on the estimation of lightcurve parameters to derive precision distances that leads to the estimation of cosmological parameters. The empirical SALT2 lightcurve modeling that relies on only two parameters, a stretch x1, and a color c, has been used by the community for almost two decades. In this paper we study the ability of the SALT2 model to fit the ne… ▽ More Type Ia supernova (SN Ia) cosmology relies on the estimation of lightcurve parameters to derive precision distances that leads to the estimation of cosmological parameters. The empirical SALT2 lightcurve modeling that relies on only two parameters, a stretch x1, and a color c, has been used by the community for almost two decades. In this paper we study the ability of the SALT2 model to fit the nearly 3000 cosmology-grade SN Ia lightcurves from the second release of the Zwicky Transient Facility (ZTF) cosmology science working group. While the ZTF data was not used to train SALT2, the algorithm is modeling the ZTF SN Ia optical lightcurves remarkably well, except for lightcurve points prior to -10 d from maximum, where the training critically lacks statistics. We find that the lightcurve fitting is robust against the considered choice of phase-range, but we show the [-10; +40] d range to be optimal in terms of statistics and accuracy. We do not detect any significant features in the lightcurve fit residuals that could be connected to the host environment. Potential systematic population differences related to the SN Ia host properties might thus not be accountable for by the addition of extra lightcurve parameters. However, a small but significant inconsistency between residuals of blue- and red-SN Ia strongly suggests the existence of a phase-dependent color term, with potential implications for the use of SNe Ia in precision cosmology. We thus encourage modellers to explore this avenue and we emphasize the importance that SN Ia cosmology must include a SALT2 retraining to accurately model the lightcurves and avoid biasing the derivation of cosmological parameters. △ Less

Submitted 4 June, 2024; originally announced June 2024.

Comments: 10 pages, 9 figures. Submitted to Astronomy and Astrophysics

ZTF SN Ia DR2: Environmental dependencies of stretch and luminosity of a volume limited sample of 1,000 Type Ia Supernovae

Authors: M. Ginolin, M. Rigault, M. Smith, Y. Copin, F. Ruppin, G. Dimitriadis, A. Goobar, J. Johansson, K. Maguire, J. Nordin, M. Amenouche, M. Aubert, C. Barjou-Delayre, M. Betoule, U. Burgaz, B. Carreres, M. Deckers, S. Dhawan, F. Feinstein, D. Fouchez, L. Galbany, C. Ganot, L. Harvey, T. de Jaeger, W. D. Kenworthy , et al. (21 additional authors not shown)

Abstract: To get distances, Type Ia Supernovae magnitudes are corrected for their correlation with lightcurve width and colour. Here we investigate how this standardisation is affected by the SN environment, with the aim to reduce scatter and improve standardisation. We first study the SN Ia stretch distribution, as well as its dependence on environment, as characterised by local and global (g-z) colour and… ▽ More To get distances, Type Ia Supernovae magnitudes are corrected for their correlation with lightcurve width and colour. Here we investigate how this standardisation is affected by the SN environment, with the aim to reduce scatter and improve standardisation. We first study the SN Ia stretch distribution, as well as its dependence on environment, as characterised by local and global (g-z) colour and stellar mass. We then look at the standardisation parameter $α$, which accounts for the correlation between residuals and stretch, along with its environment dependence and linearity. We finally compute magnitude offsets between SNe in different astrophysical environments after colour and stretch standardisation, aka steps. This analysis is made possible due to the unprecedented statistics of the ZTF SN Ia DR2 volume-limited sample. The stretch distribution exhibits a bimodal behaviour, as previously found in literature. However, we find the distribution means to decrease with host stellar mass at a 9.0$σ$ significance. We demonstrate, at the 14.3$σ$ level, that the stretch-magnitude relation is non-linear, challenging the usual linear stretch-residuals relation. Fitting for a broken-$α$ model, we indeed find two different slopes between stretch regimes ($x_1<-0.49\pm0.06$): $α_{low}=0.28\pm0.01$ and $α_{high}=0.09\pm0.01$, a $Δ_α=-0.19\pm0.01$ difference. As the relative proportion of SNe Ia in the high-/low-stretch modes evolves with redshift and environment, this implies that a linear $α$ also evolves with redshift and environment. Concerning the environmental magnitude offset $γ$, we find it to be greater than 0.14 mag regardless of the considered environmental tracer used (local or global colour and stellar mass), all measured at the $\geq 6σ$ level, increased to $\sim0.18\pm0.01$ mag when accounting for the stretch-non linearity. △ Less

Submitted 31 May, 2024; originally announced May 2024.

Comments: 15 pages, 9 figures, submitted to Astronomy and Astrophysics

Candidate strongly-lensed Type Ia supernovae in the Zwicky Transient Facility archive

Authors: A. Townsend, J. Nordin, A. Sagués Carracedo, M. Kowalski, N. Arendse, S. Dhawan, A. Goobar, J. Johansson, E. Mörtsell, S. Schulze, I. Andreoni, E. Fernández, A. G. Kim, P. E. Nugent, F. Prada, M. Rigault, N. Sarin, D. Sharma, E. C. Bellm, M. W. Coughlin, R. Dekany, S. L. Groom, L. Lacroix, R. R. Laher, R. Riddle , et al. (39 additional authors not shown)

Abstract: Gravitationally lensed Type Ia supernovae (glSNe Ia) are unique astronomical tools for studying cosmological parameters, distributions of dark matter, the astrophysics of the supernovae and the intervening lensing galaxies themselves. Only a few highly magnified glSNe Ia have been discovered by ground-based telescopes, such as the Zwicky Transient Facility (ZTF), but simulations predict the existe… ▽ More Gravitationally lensed Type Ia supernovae (glSNe Ia) are unique astronomical tools for studying cosmological parameters, distributions of dark matter, the astrophysics of the supernovae and the intervening lensing galaxies themselves. Only a few highly magnified glSNe Ia have been discovered by ground-based telescopes, such as the Zwicky Transient Facility (ZTF), but simulations predict the existence of a fainter, undetected population. We present a systematic search in the ZTF archive of alerts from 1 June 2019 to 1 September 2022. Using the AMPEL platform, we developed a pipeline that distinguishes candidate glSNe Ia from other variable sources. Initial cuts were applied to the ZTF alert photometry before forced photometry was obtained for the remaining candidates. Additional cuts were applied to refine the candidates based on their light curve colours, lens galaxy colours, and the resulting parameters from fits to the SALT2 SN Ia template. Candidates were also cross-matched with the DESI spectroscopic catalogue. Seven transients passed all the cuts and had an associated galaxy DESI redshift, which we present as glSN Ia candidates. While superluminous supernovae (SLSNe) cannot be fully rejected, two events, ZTF19abpjicm and ZTF22aahmovu, are significantly different from typical SLSNe and their light curves can be modelled as two-image glSN Ia systems. From this two-image modelling, we estimate time delays of 22 $\pm$ 3 and 34 $\pm$ 1 days for the two events, respectively, which suggests that we have uncovered a population with longer time delays. The pipeline is efficient and sensitive enough to parse full alert streams. It is currently being applied to the live ZTF alert stream to identify and follow-up future candidates while active. This pipeline could be the foundation for glSNe Ia searches in future surveys, like the Vera C. Rubin Observatory's Legacy Survey of Space and Time. △ Less

Submitted 28 May, 2024; originally announced May 2024.

Comments: 21 pages, 15 figures

Searching for gravitational wave optical counterparts with the Zwicky Transient Facility: summary of O4a

Authors: Tomás Ahumada, Shreya Anand, Michael W. Coughlin, Vaidehi Gupta, Mansi M. Kasliwal, Viraj R. Karambelkar, Robert D. Stein, Gaurav Waratkar, Vishwajeet Swain, Theophile Jegou du Laz, Akash Anumarlapudi, Igor Andreoni, Mattia Bulla, Gokul P. Srinivasaragavan, Andrew Toivonen, Avery Wold, Eric C. Bellm, S. Bradley Cenko, David L. Kaplan, Jesper Sollerman, Varun Bhalerao, Daniel Perley, Anirudh Salgundi, Aswin Suresh, K-Ryan Hinds , et al. (27 additional authors not shown)

Abstract: During the first half of the fourth observing run (O4a) of the International Gravitational Wave Network (IGWN), the Zwicky Transient Facility (ZTF) conducted a systematic search for kilonova (KN) counterparts to binary neutron star (BNS) and neutron star-black hole (NSBH) merger candidates. Here, we present a comprehensive study of the five high-significance (FAR < 1 per year) BNS and NSBH candida… ▽ More During the first half of the fourth observing run (O4a) of the International Gravitational Wave Network (IGWN), the Zwicky Transient Facility (ZTF) conducted a systematic search for kilonova (KN) counterparts to binary neutron star (BNS) and neutron star-black hole (NSBH) merger candidates. Here, we present a comprehensive study of the five high-significance (FAR < 1 per year) BNS and NSBH candidates in O4a. Our follow-up campaigns relied on both target-of-opportunity observations (ToO) and re-weighting of the nominal survey schedule to maximize coverage. We describe the toolkit we have been developing, Fritz, an instance of SkyPortal, instrumental in coordinating and managing our telescope scheduling, candidate vetting, and follow-up observations through a user-friendly interface. ZTF covered a total of 2841 deg$^2$ within the skymaps of the high-significance GW events, reaching a median depth of g~20.2 mag. We circulated 15 candidates, but found no viable KN counterpart to any of the GW events. Based on the ZTF non-detections of the high-significance events in O4a, we used a Bayesian approach, nimbus, to quantify the posterior probability of KN model parameters that are consistent with our non-detections. Our analysis favors KNe with initial absolute magnitude fainter than -16 mag. The joint posterior probability of a GW170817-like KN associated with all our O4a follow-ups was 64%. Additionally, we use a survey simulation software, simsurvey, to determine that our combined filtered efficiency to detect a GW170817-like KN is 36%, when considering the 5 confirmed astrophysical events in O3 (1 BNS and 4 NSBH), along with our O4a follow-ups. Following Kasliwal et al. (2020), we derived joint constraints on the underlying KN luminosity function based on our O3 and O4a follow-ups, determining that no more than 76% of KNe fading at 1 mag/day can peak at a magnitude brighter than -17.5 mag. △ Less

Submitted 20 May, 2024; originally announced May 2024.

Comments: submitted

Dramatic rebrightening of the type-changing stripped-envelope supernova SN 2023aew

Authors: Yashvi Sharma, Jesper Sollerman, Shrinivas R. Kulkarni, Takashi J. Moriya, Steve Schulze, Stan Barmentloo, Michael Fausnaugh, Avishay Gal-Yam, Anders Jerkstrand, Tomás Ahumada, Eric C. Bellm, Kaustav K. Das, Andrew Drake, Christoffer Fremling, Saarah Hall, K. R. Hinds, Theophile Jegou du Laz, Viraj Karambelkar, Mansi M. Kasliwal, Frank J. Masci, Adam A. Miller, Guy Nir, Daniel A. Perley, Josiah N. Purdum, Yu-Jing Qin , et al. (10 additional authors not shown)

Abstract: Multi-peaked supernovae with precursors, dramatic light-curve rebrightenings, and spectral transformation are rare, but are being discovered in increasing numbers by modern night-sky transient surveys like the Zwicky Transient Facility (ZTF). Here, we present the observations and analysis of SN 2023aew, which showed a dramatic increase in brightness following an initial luminous (-17.4 mag) and lo… ▽ More Multi-peaked supernovae with precursors, dramatic light-curve rebrightenings, and spectral transformation are rare, but are being discovered in increasing numbers by modern night-sky transient surveys like the Zwicky Transient Facility (ZTF). Here, we present the observations and analysis of SN 2023aew, which showed a dramatic increase in brightness following an initial luminous (-17.4 mag) and long (~100 days) unusual first peak (possibly precursor). SN 2023aew was classified as a Type IIb supernova during the first peak but changed its type to resemble a stripped-envelope supernova (SESN) after the marked rebrightening. We present comparisons of SN 2023aew's spectral evolution with SESN subtypes and argue that it is similar to SNe Ibc during its main peak. P-Cygni Balmer lines are present during the first peak, but vanish during the second peak's photospheric phase, before H$α$ resurfaces again during the nebular phase. The nebular lines ([O I], [Ca II], Mg I], H$α$) exhibit a double-peaked structure which hints towards a clumpy or non-spherical ejecta. We analyze the second peak in the light curve of SN 2023aew and find it to be broader than normal SESNe as well as requiring a very high $^{56}$Ni mass to power the peak luminosity. We discuss the possible origins of SN 2023aew including an eruption scenario where a part of the envelope is ejected during the first peak which also powers the second peak of the light curve through SN-CSM interaction. △ Less

Submitted 5 February, 2024; originally announced February 2024.

Comments: 22 pages, 11 figures, 5 tables

AT2019pim: A Luminous Orphan Afterglow from a Moderately Relativistic Outflow

Authors: Daniel A. Perley, Anna Y. Q. Ho, Michael Fausnaugh, Gavin P. Lamb, Mansi M. Kasliwal, Tomas Ahumada, Shreya Anand, Igor Andreoni, Eric Bellm, Varun Bhalerao, Bryce Bolin, Thomas G. Brink, Eric Burns, S. Bradley Cenko, Alessandra Corsi, Alexei V. Filippenko, Dmitry Frederiks, Adam Goldstein, Rachel Hamburg, Rahul Jayaraman, Peter G. Jonker, Erik C. Kool, Shrinivas Kulkarni, Harsh Kumar, Russ Laher , et al. (12 additional authors not shown)

Abstract: Classical gamma-ray bursts (GRBs) have two distinct emission episodes: prompt emission from ultra-relativistic ejecta and afterglow from shocked circumstellar material. While both components are extremely luminous in known GRBs, a variety of scenarios predict the existence of luminous afterglow emission with little or no associated high-energy prompt emission. We present AT 2019pim, the first secu… ▽ More Classical gamma-ray bursts (GRBs) have two distinct emission episodes: prompt emission from ultra-relativistic ejecta and afterglow from shocked circumstellar material. While both components are extremely luminous in known GRBs, a variety of scenarios predict the existence of luminous afterglow emission with little or no associated high-energy prompt emission. We present AT 2019pim, the first secure example of this phenomenon to be identified. Serendipitously discovered during follow-up observations of a gravitational-wave trigger and located in a contemporaneous TESS sector, it is hallmarked by a fast-rising (t ~ 2 hr), luminous (M_UV,peak ~ -24.4 mag) optical transient with accompanying luminous X-ray and radio emission. No gamma-ray emission consistent with the time and location of the transient was detected by Fermi-GBM or by Konus, placing strong limits on an accompanying GRB. We investigate several independent observational aspects of the afterglow in the context of constraints on relativistic motion and find all of them are consistent with an initial Lorentz factor of Gamma_0 ~ 30-50, significantly lower than in any well-observed GRB and consistent with the theoretically-predicted "dirty fireball" scenario in which the high-energy prompt emission is stifled by pair production. However, we cannot rule out a structured jet model in which only the line-of-sight material was ejected at low-Gamma, off-axis from a classical high-Gamma jet core. This event represents a milestone in orphan afterglow searches, demonstrating that luminous afterglows with weak or no detectable gamma-ray radiation exist in nature and can be discovered by high-cadence optical surveys. △ Less

Submitted 29 January, 2024; originally announced January 2024.

Comments: Submitted to MNRAS

Four new eclipsing accreting ultracompact white dwarf binaries found with the Zwicky Transient Facility

Authors: J. M. Khalil, J. van Roestel, E. C. Bellm, J. S. Bloom, R. Dekany, A. J. Drake, M. J. Graham, S. L. Groom, S. R. Kulkarni, R. R. Laher, A. A. Mahabal, T. Prince, R. Riddle

Abstract: Context. Accreting ultracompact binaries contain a white dwarf that is accreting from a degenerate object and have orbital periods shorter than 65 minutes. Aims. The aims of this letter are to report the discovery and the orbital period of four new eclipsing accreting ultracompact binaries found using the Zwicky Transient Facility, and to discuss their photometric properties. Methods. We searc… ▽ More Context. Accreting ultracompact binaries contain a white dwarf that is accreting from a degenerate object and have orbital periods shorter than 65 minutes. Aims. The aims of this letter are to report the discovery and the orbital period of four new eclipsing accreting ultracompact binaries found using the Zwicky Transient Facility, and to discuss their photometric properties. Methods. We searched through a list of 4171 dwarf novae compiled using the Zwicky Transient Facility and used the Box Least Square method to search for periodic signals in the data. Results. We found four new eclipsing accreting ultracompact binaries with orbital periods between 25.9-56 minutes, one of which is previously published as an AM CVn, while the other three systems are new discoveries. The other two shorter period systems are likely also AM CVn systems, while the longest period system with a period of 56 minutes shows multiple super-outbursts observed in two years which is more consistent with it being a Helium-CV. △ Less

Submitted 15 December, 2023; originally announced December 2023.

Comments: Submitted to A&A

Journal ref: A&A 683, L10 (2024)

Microlensing Events in Five Years of Photometry from the Zwicky Transient Facility

Authors: Ruocheng Zhai, Antonio C. Rodriguez, Casey Y. Lam, Eric C. Bellm, Josiah Purdum, Frank J. Masci, Avery Wold

Abstract: Microlensing has a unique advantage for detecting dark objects in the Milky Way, such as free floating planets, neutron stars, and stellar-mass black holes. Most microlensing surveys focus towards the Galactic bulge, where higher stellar density leads to a higher event rate. However, microlensing events in the Galactic plane are closer, and take place over longer timescales. This enables a better… ▽ More Microlensing has a unique advantage for detecting dark objects in the Milky Way, such as free floating planets, neutron stars, and stellar-mass black holes. Most microlensing surveys focus towards the Galactic bulge, where higher stellar density leads to a higher event rate. However, microlensing events in the Galactic plane are closer, and take place over longer timescales. This enables a better measurement of the microlensing parallax, which serves as an independent constraint on the mass of the dark lens. In this work, we systematically searched for microlensing events in Zwicky Transient Facility (ZTF) Data Release 17 from 2018--2023 in the Galactic plane region $|b| < 20^\circ$. We find 124 high-confidence microlensing events and 54 possible events. In the event selection, we use the efficient \texttt{EventFinder} algorithm to detect microlensing signals, which could be used for large datasets such as future ZTF data releases or data from the Rubin Observatory Legacy Survey of Space and Time (LSST). With detection efficiencies of ZTF fields from catalog-level simulations, we calculate the mean Einstein timescale to be $\langle t_\mathrm{E}\rangle = 51.7 \pm 3.3$ days, smaller than previous results of the Galactic plane to within 1.5-$σ$. We calculate optical depths and event rates, which we interpret with caution due to the use of visual inspection in creating our final sample. With two years of additional ZTF data in DR17, we have more than doubled the amount of microlensing events (60) found in the three-year DR5 search and found events with longer Einstein timescales than before. △ Less

Submitted 30 November, 2023; originally announced November 2023.

Comments: 9 figures, 3 tables. Submitted to ApJ

An Optical Search for New Outbursting Low Mass X-Ray Binaries

Authors: Yuankun Wang, Eric C. Bellm, Allison Crossland, William I. Clarkson, Alessandro Mazzi, Reed Riddle, Russ R. Laher, Ben Rusholme

Abstract: Transient Low-Mass X-ray binaries (LMXBs) are discovered largely by X-ray and gamma-ray all-sky monitors. The X-ray outburst is also accompanied by an optical brightening, which empirically can precede detection of X-rays. Newly sensitive optical synoptic surveys may offer a complementary pathway for discovery, and potential for insight into the initial onset and propagation of the thermal instabi… ▽ More Transient Low-Mass X-ray binaries (LMXBs) are discovered largely by X-ray and gamma-ray all-sky monitors. The X-ray outburst is also accompanied by an optical brightening, which empirically can precede detection of X-rays. Newly sensitive optical synoptic surveys may offer a complementary pathway for discovery, and potential for insight into the initial onset and propagation of the thermal instability that leads to the ionization of the accretion disk. We use the Zwicky Transient Facility (ZTF) alert stream to perform a comprehensive search at optical wavelengths for previously undiscovered outbursting LMXBs. Our pipeline first crossmatches the positions of the alerts to cataloged X-ray sources, and then analyzes the 30-day lightcurve of matched alerts by thresholding on differences with an 8-day exponentially weighted moving average. In addition to an nineteen month-long live search, we ran our pipeline over four years of ZTF archival data, recovering 4 known LMXBs. We also independently detected an outburst of MAXI J1957+032 in the live search and found the first outburst of Swift J1943.4+0228, an unclassified X-ray transient, in 10 years. Using Monte Carlo simulations of the Galactic LMXB population, we estimate that 29% of outbursting LMXBs are detectable by ZTF and that 4.4% of LMXBs would be present in the crossmatched X-ray catalogs, giving an estimated Galactic population of $3390^{+3980}_{-1930}$. We estimate that our current pipeline can detect 1.3% of all outbursting LMXBs, including those previously unknown, but that Rubin Observatory's Legacy Survey of Space and Time (LSST) will be able to detect 43% of outbursting LMXBs. △ Less

Submitted 29 November, 2023; originally announced November 2023.

Comments: 17 pages, 10 figures, Accepted for publication in ApJ

A Pilot Search for Gravitational Self-Lensing Binaries with the Zwicky Transient Facility

Authors: Allison Crossland, Eric C. Bellm, Courtney Klein, James R. A. Davenport, Thomas Kupfer, Steven L. Groom, Russ R. Laher, Reed Riddle

Abstract: Binary systems containing a compact object may exhibit periodic brightening episodes due to gravitational lensing as the compact object transits the companion star. Such "self-lensing" signatures have been detected before for white dwarf binaries. We attempt to use these signatures to identify detached stellar-mass neutron star and black hole binaries using data from the Zwicky Transient Facility… ▽ More Binary systems containing a compact object may exhibit periodic brightening episodes due to gravitational lensing as the compact object transits the companion star. Such "self-lensing" signatures have been detected before for white dwarf binaries. We attempt to use these signatures to identify detached stellar-mass neutron star and black hole binaries using data from the Zwicky Transient Facility (ZTF). We present a systematic search for self-lensing signals in Galactic binaries from a subset of high-cadence ZTF data taken in 2018. We identify 19 plausible candidates from the search, although because each candidate is observed to only brighten once, other origins such as stellar flares are more likely. We discuss prospects for more comprehensive future searches of the ZTF data. △ Less

Submitted 29 November, 2023; originally announced November 2023.

Comments: 12 pages. Submitted to the Open Journal of Astrophysics

Minutes-duration Optical Flares with Supernova Luminosities

Authors: Anna Y. Q. Ho, Daniel A. Perley, Ping Chen, Steve Schulze, Vik Dhillon, Harsh Kumar, Aswin Suresh, Vishwajeet Swain, Michael Bremer, Stephen J. Smartt, Joseph P. Anderson, G. C. Anupama, Supachai Awiphan, Sudhanshu Barway, Eric C. Bellm, Sagi Ben-Ami, Varun Bhalerao, Thomas de Boer, Thomas G. Brink, Rick Burruss, Poonam Chandra, Ting-Wan Chen, Wen-Ping Chen, Jeff Cooke, Michael W. Coughlin , et al. (52 additional authors not shown)

Abstract: In recent years, certain luminous extragalactic optical transients have been observed to last only a few days. Their short observed duration implies a different powering mechanism from the most common luminous extragalactic transients (supernovae) whose timescale is weeks. Some short-duration transients, most notably AT2018cow, display blue optical colours and bright radio and X-ray emission. Seve… ▽ More In recent years, certain luminous extragalactic optical transients have been observed to last only a few days. Their short observed duration implies a different powering mechanism from the most common luminous extragalactic transients (supernovae) whose timescale is weeks. Some short-duration transients, most notably AT2018cow, display blue optical colours and bright radio and X-ray emission. Several AT2018cow-like transients have shown hints of a long-lived embedded energy source, such as X-ray variability, prolonged ultraviolet emission, a tentative X-ray quasiperiodic oscillation, and large energies coupled to fast (but subrelativistic) radio-emitting ejecta. Here we report observations of minutes-duration optical flares in the aftermath of an AT2018cow-like transient, AT2022tsd (the "Tasmanian Devil"). The flares occur over a period of months, are highly energetic, and are likely nonthermal, implying that they arise from a near-relativistic outflow or jet. Our observations confirm that in some AT2018cow-like transients the embedded energy source is a compact object, either a magnetar or an accreting black hole. △ Less

Submitted 16 November, 2023; originally announced November 2023.

Comments: 79 pages, 3 figures (main text) + 7 figures (extended data) + 2 figures (supplementary information). Published online in Nature on 15 November 2023

The Early Ultraviolet Light-Curves of Type II Supernovae and the Radii of Their Progenitor Stars

Authors: Ido Irani, Jonathan Morag, Avishay Gal-Yam, Eli Waxman, Steve Schulze, Jesper Sollerman, K-Ryan Hinds, Daniel A. Perley, Ping Chen, Nora L. Strotjohann, Ofer Yaron, Erez A. Zimmerman, Rachel Bruch, Eran O. Ofek, Maayane T. Soumagnac, Yi Yang, Steven L. Groom, Frank J. Masci, Reed Riddle, Eric C. Bellm, David Hale

Abstract: We present a sample of 34 normal SNe II detected with the Zwicky Transient Facility, with multi-band UV light-curves starting at $t \leq 4$ days after explosion, as well as X-ray detections and upper limits. We characterize the early UV-optical colors and provide prescriptions for empirical host-extinction corrections. We show that the $t > 2\,$days UV-optical colors and the blackbody evolution of… ▽ More We present a sample of 34 normal SNe II detected with the Zwicky Transient Facility, with multi-band UV light-curves starting at $t \leq 4$ days after explosion, as well as X-ray detections and upper limits. We characterize the early UV-optical colors and provide prescriptions for empirical host-extinction corrections. We show that the $t > 2\,$days UV-optical colors and the blackbody evolution of the sample are consistent with the predictions of spherical phase shock-cooling (SC), independently of the presence of `flash ionization" features. We present a framework for fitting SC models which can reproduce the parameters of a set of multi-group simulations without a significant bias up to 20% in radius and velocity. Observations of about half of the SNe II in the sample are well-fit by models with breakout radii $<10^{14}\,$cm. The other half are typically more luminous, with observations from day 1 onward that are better fit by a model with a large $>10^{14}\,$cm breakout radius. However, these fits predict an early rise during the first day that is too slow. We suggest these large-breakout events are explosions of stars with an inflated envelope or a confined CSM with a steep density profile, at which breakout occurs. Using the X-ray data, we derive constraints on the extended ($\sim10^{15}$ cm) CSM density independent of spectral modeling, and find most SNe II progenitors lose $<10^{-4} M_{\odot}\, \rm yr^{-1}$ a few years before explosion. This provides independent evidence the CSM around many SNe II progenitors is confined. We show that the overall observed breakout radius distribution is skewed to higher radii due to a luminosity bias. We argue that the $66^{+11}_{-22}\%$ of red supergiants (RSG) explode as SNe II with breakout radii consistent with the observed distribution of field RSG, with a tail extending to large radii, likely due to the presence of CSM. △ Less

Submitted 14 April, 2024; v1 submitted 25 October, 2023; originally announced October 2023.

Comments: Submitted to ApJ. Comments are welcome at ido.irani@weizmann.ac.il or idoirani@gmail.com. Accepted version

Characterizing the Ordinary Broad-lined Type Ic SN 2023pel from the Energetic GRB 230812B

Authors: Gokul P. Srinivasaragavan, Vishwajeet Swain, Brendan M. O'Connor, Shreya Anand, Tomás Ahumada, Daniel A. Perley, Robert Stein, Jesper Sollerman, Christoffer Fremling, S. Bradley Cenko, Sarah Antier, Nidhal Guessoum, Thomas Hussenot-Desenonges, Patrice Hello, Stephen Lesage, Erica Hammerstein, M. Coleman Miller, Igor Andreoni, Varun Bhalerao, Joshua S. Bloom, Anirban Dutta, Avishay Gal-Yam, K-Ryan Hinds, Amruta D. Jaodand, Mansi M. Kasliwal , et al. (17 additional authors not shown)

Abstract: We report observations of the optical counterpart of the long gamma-ray burst (LGRB) GRB 230812B, and its associated supernova (SN) SN 2023pel. The proximity ($z = 0.36$) and high energy ($E_{γ, \rm{iso}} \sim 10^{53}$ erg) make it an important event to study as a probe of the connection between massive star core-collapse and relativistic jet formation. With a phenomenological power-law model for… ▽ More We report observations of the optical counterpart of the long gamma-ray burst (LGRB) GRB 230812B, and its associated supernova (SN) SN 2023pel. The proximity ($z = 0.36$) and high energy ($E_{γ, \rm{iso}} \sim 10^{53}$ erg) make it an important event to study as a probe of the connection between massive star core-collapse and relativistic jet formation. With a phenomenological power-law model for the optical afterglow, we find a late-time flattening consistent with the presence of an associated SN. SN 2023pel has an absolute peak $r$-band magnitude of $M_r = -19.46 \pm 0.18$ mag (about as bright as SN 1998bw) and evolves on quicker timescales. Using a radioactive heating model, we derive a nickel mass powering the SN of $M_{\rm{Ni}} = 0.38 \pm 0.01$ $\rm{M_\odot}$, and a peak bolometric luminosity of $L_{\rm{bol}} \sim 1.3 \times 10^{43}$ $\rm{erg}$ $\rm{s^{-1}}$. We confirm SN 2023pel's classification as a broad-lined Type Ic SN with a spectrum taken 15.5 days after its peak in $r$ band, and derive a photospheric expansion velocity of $v_{\rm{ph}} = 11,300 \pm 1,600$ $\rm{km}$ $\rm{s^{-1}}$ at that phase. Extrapolating this velocity to the time of maximum light, we derive the ejecta mass $M_{\rm{ej}} = 1.0 \pm 0.6$ $\rm{M_\odot}$ and kinetic energy $E_{\rm{KE}} = 1.3^{+3.3}_{-1.2} \times10^{51}$ $\rm{erg}$. We find that GRB 230812B/SN 2023pel has SN properties that are mostly consistent with the overall GRB-SN population. The lack of correlations found in the GRB-SN population between SN brightness and $E_{γ, \rm{iso}}$ for their associated GRBs, across a broad range of 7 orders of magnitude, provides further evidence that the central engine powering the relativistic ejecta is not coupled to the SN powering mechanism in GRB-SN systems. △ Less

Submitted 9 December, 2023; v1 submitted 22 October, 2023; originally announced October 2023.

Comments: 21 pages, 7 Figures, 1 Table, Accepted to ApJ Letters

A 12.4 day periodicity in a close binary system after a supernova

Authors: Ping Chen, Avishay Gal-Yam, Jesper Sollerman, Steve Schulze, Richard S. Post, Chang Liu, Eran O. Ofek, Kaustav K. Das, Christoffer Fremling, Assaf Horesh, Boaz Katz, Doron Kushnir, Mansi M. Kasliwal, Shri R. Kulkarni, Dezi Liu, Xiangkun Liu, Adam A. Miller, Kovi Rose, Eli Waxman, Sheng Yang, Yuhan Yao, Barak Zackay, Eric C. Bellm, Richard Dekany, Andrew J. Drake , et al. (15 additional authors not shown)

Abstract: Neutron stars and stellar-mass black holes are the remnants of massive star explosions. Most massive stars reside in close binary systems, and the interplay between the companion star and the newly formed compact object has been theoretically explored, but signatures for binarity or evidence for the formation of a compact object during a supernova explosion are still lacking. Here we report a stri… ▽ More Neutron stars and stellar-mass black holes are the remnants of massive star explosions. Most massive stars reside in close binary systems, and the interplay between the companion star and the newly formed compact object has been theoretically explored, but signatures for binarity or evidence for the formation of a compact object during a supernova explosion are still lacking. Here we report a stripped-envelope supernova, SN 2022jli, which shows 12.4-day periodic undulations during the declining light curve. Narrow H$α$ emission is detected in late-time spectra with concordant periodic velocity shifts, likely arising from hydrogen gas stripped from a companion and accreted onto the compact remnant. A new Fermi/LAT $γ$-ray source is temporally and positionally consistent with SN 2022jli. The observed properties of SN 2022jli, including periodic undulations in the optical light curve, coherent H$α$ emission shifting, and evidence for association with a $γ$-ray source, point to the explosion of a massive star in a binary system leaving behind a bound compact remnant. Mass accretion from the companion star onto the compact object powers the light curve of the supernova and generates the $γ$-ray emission. △ Less

Submitted 11 October, 2023; originally announced October 2023.

Comments: Accepted for publication in Nature

The first systematically identified repeating partial tidal disruption event

Authors: Jean J. Somalwar, Vikram Ravi, Yuhan Yao, Muryel Guolo, Matthew Graham, Erica Hammerstein, Wenbin Lu, Matt Nicholl, Yashvi Sharma, Robert Stein, Sjoert van Velzen, Eric C. Bellm, Michael W. Coughlin, Steven L. Groom, Frank J. Masci, Reed Riddle

Abstract: Tidal disruption events (TDEs) occur when a star enters the tidal radius of a supermassive black hole (SMBH). If the star only grazes the tidal radius, a fraction of the stellar mass will be accreted in a partial TDE (pTDE). The remainder can continue orbiting and may re-disrupted at pericenter, causing a repeating pTDE. pTDEs may be as or more common than full TDEs (fTDEs), yet few are known. In… ▽ More Tidal disruption events (TDEs) occur when a star enters the tidal radius of a supermassive black hole (SMBH). If the star only grazes the tidal radius, a fraction of the stellar mass will be accreted in a partial TDE (pTDE). The remainder can continue orbiting and may re-disrupted at pericenter, causing a repeating pTDE. pTDEs may be as or more common than full TDEs (fTDEs), yet few are known. In this work, we present the discovery of the first repeating pTDE from a systematically-selected sample, AT\,2020vdq. AT\,2020vdq was originally identified as an optically- and radio-flaring TDE. Around $3$ years after its discovery, it rebrightened dramatically and rapidly in the optical. The optical flare was remarkably fast and luminous compared to previous TDEs. It was accompanied by extremely broad (${\sim}0.1c$) optical/UV spectral features and faint X-ray emission ($L_X \sim 3\times10^{41}$\,erg\,s$^{-1}$), but no new radio-emitting component. Based on the transient optical/UV spectral features and the broadband light curve, we show that AT\,2020vdq is a repeating pTDE. We then use it to constrain TDE models; in particular, we favor a star originally in a very tight binary system that is tidally broken apart by the Hills mechanism. We also constrain the repeating pTDE rate to be $10^{-6}$ to $10^{-5}$ yr$^{-1}$ galaxy$^{-1}$, with uncertainties dominated by the unknown distribution of pTDE repeat timescales. In the Hills framework, this means the binary fraction in the galactic nucleus is of the order few percent. △ Less

Submitted 5 October, 2023; originally announced October 2023.

Comments: 24 pages, 13 figures, submitted to ApJ

An Optically-Discovered Outburst from XTE J1859+226

Authors: Eric C. Bellm, Yuankun Wang, Jan van Roestel, Rebecca A. Phillipson, Michael W. Coughlin, John A. Tomsick, Steven L. Groom, Brian Healy, Josiah Purdum, Ben Rusholme, Jesper Sollerman, Peter Bealo, Stefano Lora, Eddy Muyllaert, Ivo Peretto, Erik J. Schwendeman

Abstract: Using the Zwicky Transient Facility, in 2021 February we identified the first known outburst of the Black Hole X-ray Transient XTE J1859+226 since its discovery in 1999. The outburst was visible at X-ray, UV, and optical wavelengths for less than 20 days, substantially shorter than its 320-day full outburst in 1999, and the observed peak luminosity was two orders of magnitude lower. Its peak bolom… ▽ More Using the Zwicky Transient Facility, in 2021 February we identified the first known outburst of the Black Hole X-ray Transient XTE J1859+226 since its discovery in 1999. The outburst was visible at X-ray, UV, and optical wavelengths for less than 20 days, substantially shorter than its 320-day full outburst in 1999, and the observed peak luminosity was two orders of magnitude lower. Its peak bolometric luminosity was only $2\times 10^{35}$ erg s$^{-1}$, implying an Eddington fraction of about $3\times10^{-4}$. The source remained in the hard spectral state throughout the outburst. From optical spectroscopy measurements we estimate an outer disk radius of 10$^{11}$ cm. The low observed X-ray luminosity is not sufficient to irradiate the entire disk, but we observe a surprising exponential decline in the X-ray lightcurve. These observations highlight the potential of optical and infrared (O/IR) synoptic surveys to discover low-luminosity activity from X-ray transients. △ Less

Submitted 19 September, 2023; originally announced September 2023.

Comments: 12 pages, 6 figures, accepted for publication in ApJ

A rotating white dwarf shows different compositions on its opposite faces

Authors: Ilaria Caiazzo, Kevin B. Burdge, Pier-Emmanuel Tremblay, James Fuller, Lilia Ferrario, Boris T. Gaensicke, J. J. Hermes, Jeremy Heyl, Adela Kawka, S. R. Kulkarni, Thomas R. Marsh, Przemek Mroz, Thomas A. Prince, Harvey B. Richer, Antonio C. Rodriguez, Jan van Roestel, Zachary P. Vanderbosch, Stephane Vennes, Dayal Wickramasinghe, Vikram S. Dhillon, Stuart P. Littlefair, James Munday, Ingrid Pelisoli, Daniel Perley, Eric C. Bellm , et al. (13 additional authors not shown)

Abstract: White dwarfs, the extremely dense remnants left behind by most stars after their death, are characterised by a mass comparable to that of the Sun compressed into the size of an Earth-like planet. In the resulting strong gravity, heavy elements sink toward the centre and the upper layer of the atmosphere contains only the lightest element present, usually hydrogen or helium. Several mechanisms comp… ▽ More White dwarfs, the extremely dense remnants left behind by most stars after their death, are characterised by a mass comparable to that of the Sun compressed into the size of an Earth-like planet. In the resulting strong gravity, heavy elements sink toward the centre and the upper layer of the atmosphere contains only the lightest element present, usually hydrogen or helium. Several mechanisms compete with gravitational settling to change a white dwarf's surface composition as it cools, and the fraction of white dwarfs with helium atmospheres is known to increase by a factor ~2.5 below a temperature of about 30,000 K; therefore, some white dwarfs that appear to have hydrogen-dominated atmospheres above 30,000 K are bound to transition to be helium-dominated as they cool below it. Here we report observations of ZTF J203349.8+322901.1, a transitioning white dwarf with two faces: one side of its atmosphere is dominated by hydrogen and the other one by helium. This peculiar nature is likely caused by the presence of a small magnetic field, which creates an inhomogeneity in temperature, pressure or mixing strength over the surface. ZTF J203349.8+322901.1 might be the most extreme member of a class of magnetic, transitioning white dwarfs -- together with GD 323, a white dwarf that shows similar but much more subtle variations. This new class could help shed light on the physical mechanisms behind white dwarf spectral evolution. △ Less

Submitted 14 August, 2023; originally announced August 2023.

Comments: 45 pages, 11 figures

Journal ref: Nature 620, 61-66 (2023)

Swift Deep Galactic Plane Survey Classification of Swift J170800$-$402551.8 as a Candidate Intermediate Polar Cataclysmic Variable

Authors: B. O'Connor, E. Gogus, J. Hare, K. Mukai, D. Huppenkothen, J. Brink, D. A. H. Buckley, A. Levan, M. G. Baring, R. Stewart, C. Kouveliotou, P. Woudt, E. Bellm, S. B. Cenko, P. A. Evans, J. Granot, C. Hailey, F. Harrison, D. Hartmann, A. J. van der Horst, L. Kaper, J. A. Kennea, S. B. Potter, P. O. Slane, D. Stern , et al. (2 additional authors not shown)

Abstract: Here, we present the results of our multi-wavelength campaign aimed at classifying \textit{Swift} J170800$-$402551.8 as part of the \textit{Swift} Deep Galactic Plane Survey (DGPS). We utilized Target of Opportunity (ToO) observations with \textit{Swift}, \textit{NICER}, \textit{XMM-Newton}, \textit{NuSTAR}, and the Southern African Large Telescope (SALT), as well as multi-wavelength archival obse… ▽ More Here, we present the results of our multi-wavelength campaign aimed at classifying \textit{Swift} J170800$-$402551.8 as part of the \textit{Swift} Deep Galactic Plane Survey (DGPS). We utilized Target of Opportunity (ToO) observations with \textit{Swift}, \textit{NICER}, \textit{XMM-Newton}, \textit{NuSTAR}, and the Southern African Large Telescope (SALT), as well as multi-wavelength archival observations from \textit{Gaia}, VPHAS, and VVV. The source displays a periodicity of 784 s in our \textit{XMM-Newton} observation. The X-ray spectrum (\textit{XMM-Newton} and \textit{NuSTAR}) can be described by thermal bremsstrahlung radiation with a temperature of $kT$\,$\approx$\,$30$ keV. The phase-folded X-ray lightcurve displays a double-peaked, energy-dependent pulse-profile. We used \textit{Chandra} to precisely localize the source, allowing us to identify and study the multi-wavelength counterpart. Spectroscopy with SALT identified a Balmer H$α$ line, and potential HeI lines, from the optical counterpart. The faintness of the counterpart ($r$\,$\approx$\,$21$ AB mag) favors a low-mass donor star. Based on these criteria, we classify \textit{Swift} J170800$-$402551.8 as a candidate intermediate polar cataclysmic variable, where the spin period of the white dwarf is 784 s. △ Less

Submitted 28 August, 2023; v1 submitted 26 July, 2023; originally announced July 2023.

Comments: Accepted in MNRAS

The Swift Deep Galactic Plane Survey (DGPS) Phase-I Catalog

Authors: B. O'Connor, C. Kouveliotou, P. A. Evans, N. Gorgone, A. J. van Kooten, S. Gagnon, H. Yang, M. G. Baring, E. Bellm, P. Beniamini, J. Brink, D. A. H. Buckley, S. B. Cenko, O. D. Egbo, E. Gogus, J. Granot, C. Hailey, J. Hare, F. Harrison, D. Hartmann, A. J. van der Horst, D. Huppenkothen, L. Kaper, O. Kargaltsev, J. A. Kennea , et al. (8 additional authors not shown)

Abstract: The \textit{Swift} Deep Galactic Plane Survey is a \textit{Swift} Key Project consisting of 380 tiled pointings covering 40 deg$^{2}$ of the Galactic Plane between longitude $10$\,$<$\,$|l|$\,$<$\,$30$ deg and latitude $|b|$\,$<$\,$0.5$ deg. Each pointing has a $5$ ks exposure, yielding a total of 1.9 Ms spread across the entire survey footprint. Phase-I observations were carried out between March… ▽ More The \textit{Swift} Deep Galactic Plane Survey is a \textit{Swift} Key Project consisting of 380 tiled pointings covering 40 deg$^{2}$ of the Galactic Plane between longitude $10$\,$<$\,$|l|$\,$<$\,$30$ deg and latitude $|b|$\,$<$\,$0.5$ deg. Each pointing has a $5$ ks exposure, yielding a total of 1.9 Ms spread across the entire survey footprint. Phase-I observations were carried out between March 2017 and May 2021. The Survey is complete to depth $L_X$\,$>$\,$10^{34}$ erg s$^{-1}$ to the edge of the Galaxy. The main Survey goal is to produce a rich sample of new X-ray sources and transients, while also covering a broad discovery space. Here, we introduce the Survey strategy and present a catalog of sources detected during Phase-I observations. In total, we identify 928 X-ray sources, of which 348 are unique to our X-ray catalog. We report on the characteristics of sources in our catalog and highlight sources newly classified and published by the DGPS team. △ Less

Submitted 21 September, 2023; v1 submitted 25 June, 2023; originally announced June 2023.

Comments: Accepted in ApJS. This is the final version

Discovery of Gaia17bpp, a Giant Star with the Deepest and Longest Known Dimming Event

Authors: Anastasios Tzanidakis, James R. A. Davenport, Eric C. Bellm, Yuankun Wang

Abstract: We report the serendipitous discovery of Gaia17bpp/2MASS J19372316+1759029, a star with a deep single large-amplitude dimming event of $\sim$4.5 magnitudes that lasted over 6.5 years. Using the optical to IR spectral energy distribution (SED), we constrain the primary star to be a cool giant M0-III star with effective temperature $T_{\text{eff}}$=3,850 K and radius R=58 R$_{\odot}$. Based on the S… ▽ More We report the serendipitous discovery of Gaia17bpp/2MASS J19372316+1759029, a star with a deep single large-amplitude dimming event of $\sim$4.5 magnitudes that lasted over 6.5 years. Using the optical to IR spectral energy distribution (SED), we constrain the primary star to be a cool giant M0-III star with effective temperature $T_{\text{eff}}$=3,850 K and radius R=58 R$_{\odot}$. Based on the SED fitting, we obtained a bimodal posterior distribution of primary stellar masses at 1.5 M${\odot}$ and 3.7 M${\odot}$. Within the last 66 years of photometric coverage, no other significant dimming events of this depth and duration were identified in the optical light curves. Using a Gaussian Process, we fit a high-order Gaussian model to the optical and IR light curves and conclude the dimming event exhibits moderate asymmetries from optical to IR. At the minimum of the dimming event, the (W$_{1}$-W$_{2}$) color was bluer by $\sim$0.2 mag relative to the primary star outside the dimming event. The ingress and egress colors show a shallow reddening profile. We suggest that the main culprit of the dimming event is likely due to the presence of a large, optically thick disk transiting the primary giant star. By fitting a monochromatic transit model of an oblate disk transiting a star, we found good agreement with a slow-moving, 0.005 km sec$^{-1}$, disk with a $\sim$1.4 AU radius. We propose that Gaia17bpp belongs to a rare binary star population similar to the Epsilon Aurigae system, which consists of a secondary star enshrouded by an optically thick debris disk. △ Less

Submitted 21 June, 2023; originally announced June 2023.

Comments: 16 pages, 12 figures, submitted to ApJ, comments welcome!

NANCY: Next-generation All-sky Near-infrared Community surveY

Authors: Jiwon Jesse Han, Arjun Dey, Adrian M. Price-Whelan, Joan Najita, Edward F. Schlafly, Andrew Saydjari, Risa H. Wechsler, Ana Bonaca, David J Schlegel, Charlie Conroy, Anand Raichoor, Alex Drlica-Wagner, Juna A. Kollmeier, Sergey E. Koposov, Gurtina Besla, Hans-Walter Rix, Alyssa Goodman, Douglas Finkbeiner, Abhijeet Anand, Matthew Ashby, Benedict Bahr-Kalus, Rachel Beaton, Jayashree Behera, Eric F. Bell, Eric C Bellm , et al. (184 additional authors not shown)

Abstract: The Nancy Grace Roman Space Telescope is capable of delivering an unprecedented all-sky, high-spatial resolution, multi-epoch infrared map to the astronomical community. This opportunity arises in the midst of numerous ground- and space-based surveys that will provide extensive spectroscopy and imaging together covering the entire sky (such as Rubin/LSST, Euclid, UNIONS, SPHEREx, DESI, SDSS-V, GAL… ▽ More The Nancy Grace Roman Space Telescope is capable of delivering an unprecedented all-sky, high-spatial resolution, multi-epoch infrared map to the astronomical community. This opportunity arises in the midst of numerous ground- and space-based surveys that will provide extensive spectroscopy and imaging together covering the entire sky (such as Rubin/LSST, Euclid, UNIONS, SPHEREx, DESI, SDSS-V, GALAH, 4MOST, WEAVE, MOONS, PFS, UVEX, NEO Surveyor, etc.). Roman can uniquely provide uniform high-spatial-resolution (~0.1 arcsec) imaging over the entire sky, vastly expanding the science reach and precision of all of these near-term and future surveys. This imaging will not only enhance other surveys, but also facilitate completely new science. By imaging the full sky over two epochs, Roman can measure the proper motions for stars across the entire Milky Way, probing 100 times fainter than Gaia out to the very edge of the Galaxy. Here, we propose NANCY: a completely public, all-sky survey that will create a high-value legacy dataset benefiting innumerable ongoing and forthcoming studies of the universe. NANCY is a pure expression of Roman's potential: it images the entire sky, at high spatial resolution, in a broad infrared bandpass that collects as many photons as possible. The majority of all ongoing astronomical surveys would benefit from incorporating observations of NANCY into their analyses, whether these surveys focus on nearby stars, the Milky Way, near-field cosmology, or the broader universe. △ Less

Submitted 20 June, 2023; originally announced June 2023.

Comments: Submitted to the call for white papers for the Roman Core Community Survey (June 16th, 2023), and to the Bulletin of the AAS

Identification of 1RXS J165424.6-433758 as a polar cataclysmic variable

Authors: B. O'Connor, J. Brink, D. A. H. Buckley, K. Mukai, C. Kouveliotou, E. Gogus, S. B. Potter, P. Woudt, A. Lien, A. Levan, O. Kargaltsev, M. G. Baring, E. Bellm, S. B. Cenko, P. A. Evans, J. Granot, C. Hailey, F. Harrison, D. Hartmann, A. J. van der Horst, D. Huppenkothen, L. Kaper, J. A. Kennea, P. O. Slane, D. Stern , et al. (3 additional authors not shown)

Abstract: We present the results of our X-ray, ultraviolet, and optical follow-up campaigns of 1RXS J165424.6-433758, an X-ray source detected with the \textit{Swift} Deep Galactic Plane Survey (DGPS). The source X-ray spectrum (\textit{Swift} and \textit{NuSTAR}) is described by thermal bremsstrahlung radiation with a temperature of $kT=10.1\pm1.2$ keV, yielding an X-ray ($0.3-10$ keV) luminosity… ▽ More We present the results of our X-ray, ultraviolet, and optical follow-up campaigns of 1RXS J165424.6-433758, an X-ray source detected with the \textit{Swift} Deep Galactic Plane Survey (DGPS). The source X-ray spectrum (\textit{Swift} and \textit{NuSTAR}) is described by thermal bremsstrahlung radiation with a temperature of $kT=10.1\pm1.2$ keV, yielding an X-ray ($0.3-10$ keV) luminosity $L_X=(6.5\pm0.8)\times10^{31}$ erg s$^{-1}$ at a \textit{Gaia} distance of 460 pc. Spectroscopy with the Southern African Large Telescope (SALT) revealed a flat continuum dominated by emission features, demonstrating an inverse Balmer decrement, the $\lambda4640$ Bowen blend, almost a dozen HeI lines, and HeII $\lambda4541$, $\lambda4686$ and $λ5411$. Our high-speed photometry demonstrates a preponderance of flickering and flaring episodes, and revealed the orbital period of the system, $P_\textrm{orb}=2.87$ hr, which fell well within the cataclysmic variable (CV) period gap between $2-3$ hr. These features classify 1RXS J165424.6-433758 as a nearby polar magnetic CV. △ Less

Submitted 10 September, 2023; v1 submitted 8 June, 2023; originally announced June 2023.

Comments: Accepted for publication in ApJ

LSST Survey Strategy in the Galactic Plane and Magellanic Clouds

Authors: R. A. Street, X. Li, S. Khakpash, E. Bellm, L. Girardi, L. Jones, N. S. Abrams, Y. Tsapras, M. P. G. Hundertmark, E. Bachelet, P. Gandhi, P. Szkody, W. I. Clarkson, R. Szabo, L. Prisinzano, R. Bonito, D. A. H. Buckley, J. P. Marais, R. Di Stefano

Abstract: Galactic science encompasses a wide range of subjects in the study of the Milky Way and Magellanic Clouds, from Young Stellar Objects to X-ray Binaries. Mapping these populations, and exploring transient phenomena within them, are among the primary science goals of the Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST). While early versions of the survey strategy dedicated relative… ▽ More Galactic science encompasses a wide range of subjects in the study of the Milky Way and Magellanic Clouds, from Young Stellar Objects to X-ray Binaries. Mapping these populations, and exploring transient phenomena within them, are among the primary science goals of the Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST). While early versions of the survey strategy dedicated relatively few visits to the Galactic Plane region, more recent strategies under consideration envision higher cadence within selected regions of high scientific interest. The range of galactic science presents a challenge in evaluating which strategies deliver the highest scientific returns. Here we present metrics designed to evaluate Rubin survey strategy simulations based on the cadence of observations they deliver within regions of interest to different topics in galactic science, using variability categories defined by timescale. We also compare the fractions of exposures obtained in each filter with those recommended for the different science goals. We find that the baseline_v2.x simulations deliver observations of the high-priority regions at sufficiently high cadence to reliably detect variability on timescales >10 d or more. Follow-up observations may be necessary to properly characterize variability, especially transients, on shorter timescales. Combining the regions of interest for all the science cases considered, we identify those areas of the Galactic Plane and Magellanic Clouds of highest priority. We recommend that these refined survey footprints be used in future simulations to explore rolling cadence scenarios, and to optimize the sequence of observations in different bandpasses. △ Less

Submitted 19 May, 2023; originally announced May 2023.

Comments: Accepted for publication in ApJS. 26 pages and 36 figures

A data science platform to enable time-domain astronomy

Authors: Michael W. Coughlin, Joshua S. Bloom, Guy Nir, Sarah Antier, Theophile Jegou du Laz, Stéfan van der Walt, Arien Crellin-Quick, Thomas Culino, Dmitry A. Duev, Daniel A. Goldstein, Brian F. Healy, Viraj Karambelkar, Jada Lilleboe, Kyung Min Shin, Leo P. Singer, Tomas Ahumada, Shreya Anand, Eric C. Bellm, Richard Dekany, Matthew J. Graham, Mansi M. Kasliwal, Ivona Kostadinova, R. Weizmann Kiendrebeogo, Shrinivas R. Kulkarni, Sydney Jenkins , et al. (28 additional authors not shown)

Abstract: SkyPortal is an open-source software package designed to efficiently discover interesting transients, manage follow-up, perform characterization, and visualize the results. By enabling fast access to archival and catalog data, cross-matching heterogeneous data streams, and the triggering and monitoring of on-demand observations for further characterization, a SkyPortal-based platform has been oper… ▽ More SkyPortal is an open-source software package designed to efficiently discover interesting transients, manage follow-up, perform characterization, and visualize the results. By enabling fast access to archival and catalog data, cross-matching heterogeneous data streams, and the triggering and monitoring of on-demand observations for further characterization, a SkyPortal-based platform has been operating at scale for 2 yr for the Zwicky Transient Facility Phase II community, with hundreds of users, containing tens of millions of time-domain sources, interacting with dozens of telescopes, and enabling community reporting. While SkyPortal emphasizes rich user experiences (UX) across common frontend workflows, recognizing that scientific inquiry is increasingly performed programmatically, SkyPortal also surfaces an extensive and well-documented API system. From backend and frontend software to data science analysis tools and visualization frameworks, the SkyPortal design emphasizes the re-use and leveraging of best-in-class approaches, with a strong extensibility ethos. For instance, SkyPortal now leverages ChatGPT large-language models (LLMs) to automatically generate and surface source-level human-readable summaries. With the imminent re-start of the next-generation of gravitational wave detectors, SkyPortal now also includes dedicated multi-messenger features addressing the requirements of rapid multi-messenger follow-up: multi-telescope management, team/group organizing interfaces, and cross-matching of multi-messenger data streams with time-domain optical surveys, with interfaces sufficiently intuitive for the newcomers to the field. (abridged) △ Less

Submitted 14 June, 2023; v1 submitted 28 April, 2023; originally announced May 2023.

Comments: Accepted to ApJS

SN 2020udy: a SN Iax with strict limits on interaction consistent with a helium-star companion

Authors: Kate Maguire, Mark R. Magee, Giorgos Leloudas, Adam A. Miller, Georgios Dimitriadis, Miika Pursiainen, Mattia Bulla, Kishalay De, Avishay Gal-Yam, Daniel A. Perley, Christoffer Fremling, Viraj R. Karambelkar, Jakob Nordin, Simeon Reusch, Steve Schulze, Jesper Sollerman, Giacomo Terreran, Yi Yang, Eric C. Bellm, Steven L. Groom, Mansi M. Kasliwal, Shrinivas R. Kulkarni, Leander Lacroix, Frank J. Masci, Josiah N. Purdum , et al. (2 additional authors not shown)

Abstract: Early observations of transient explosions can provide vital clues to their progenitor origins. In this paper we present the nearby Type Iax (02cx-like) supernova (SN), SN 2020udy that was discovered within hours ($\sim$7 hr) of estimated first light. An extensive dataset of ultra-violet, optical, and near-infrared observations was obtained, covering out to $\sim$150 d after explosion. SN 2020udy… ▽ More Early observations of transient explosions can provide vital clues to their progenitor origins. In this paper we present the nearby Type Iax (02cx-like) supernova (SN), SN 2020udy that was discovered within hours ($\sim$7 hr) of estimated first light. An extensive dataset of ultra-violet, optical, and near-infrared observations was obtained, covering out to $\sim$150 d after explosion. SN 2020udy peaked at -17.86$\pm$0.43 mag in the r band and evolved similarly to other 'luminous' SNe Iax, such as SNe 2005hk and 2012Z. Its well-sampled early light curve allows strict limits on companion interaction to be placed. Main-sequence companion stars with masses of 2 and 6 M$_\odot$ are ruled out at all viewing angles, while a helium-star companion is allowed from a narrow range of angles (140-180$^\circ$ away from the companion). The spectra and light curves of SN2020udy are in good agreement with those of the 'N5def' deflagration model of a near Chandrasekhar-mass carbon-oxygen white dwarf. However, as has been seen in previous studies of similar luminosity events, SN 2020udy evolves slower than the model. Broad-band linear polarisation measurements taken at and after peak are consistent with no polarisation, in agreement with the predictions of the companion-star configuration from the early light curve measurements. The host galaxy environment is low metallicity and is consistent with a young stellar population. Overall, we find the most plausible explosion scenario to be the incomplete disruption of a CO white dwarf near the Chandrasekhar-mass limit, with a helium-star companion. △ Less

Submitted 24 April, 2023; originally announced April 2023.

Comments: 18 pages, 14 figures, submitted to MNRAS

The broad-lined Type-Ic supernova SN 2022xxf with extraordinary two-humped light curves

Authors: H. Kuncarayakti, J. Sollerman, L. Izzo, K. Maeda, S. Yang, S. Schulze, C. R. Angus, M. Aubert, K. Auchettl, M. Della Valle, L. Dessart, K. Hinds, E. Kankare, M. Kawabata, P. Lundqvist, T. Nakaoka, D. Perley, S. I. Raimundo, N. L. Strotjohann, K. Taguchi, Y. -Z. Cai, P. Charalampopoulos, Q. Fang, M. Fraser, C. P. Gutierrez , et al. (38 additional authors not shown)

Abstract: We report on our study of supernova (SN) 2022xxf based on observations obtained during the first four months of its evolution. The light curves (LCs) display two humps of similar maximum brightness separated by 75 days, unprecedented for a broad-lined (BL) Type Ic supernova (SN IcBL). SN 2022xxf is the most nearby SN IcBL to date (in NGC 3705, $z = 0.0037$, at a distance of about 20 Mpc). Optical… ▽ More We report on our study of supernova (SN) 2022xxf based on observations obtained during the first four months of its evolution. The light curves (LCs) display two humps of similar maximum brightness separated by 75 days, unprecedented for a broad-lined (BL) Type Ic supernova (SN IcBL). SN 2022xxf is the most nearby SN IcBL to date (in NGC 3705, $z = 0.0037$, at a distance of about 20 Mpc). Optical and near-infrared photometry and spectroscopy are used to identify the energy source powering the LC. Nearly 50 epochs of high signal-to-noise-ratio spectroscopy were obtained within 130 days, comprising an unparalleled dataset for a SN IcBL, and one of the best-sampled SN datasets to date. The global spectral appearance and evolution of SN 2022xxf points to typical SN Ic/IcBL, with broad features (up to $\sim14000$ km s$^{-1}$) and a gradual transition from the photospheric to the nebular phase. However, narrow emission lines (corresponding to $\sim1000-2500$ km s$^{-1}$) are present in the spectra from the time of the second rise, suggesting slower-moving circumstellar material (CSM). These lines are subtle, in comparison to the typical strong narrow lines of CSM-interacting SNe, for example, Type IIn, Ibn, and Icn, but some are readily noticeable at late times such as in Mg I $λ$5170 and [O I] $λ$5577. Unusually, the near-infrared spectra show narrow line peaks in a number of features formed by ions of O and Mg. We infer the presence of CSM that is free of H and He. We propose that the radiative energy from the ejecta-CSM interaction is a plausible explanation for the second LC hump. This interaction scenario is supported by the color evolution, which progresses to the blue as the light curve evolves along the second hump, and the slow second rise and subsequent rapid LC drop. (Abstract abridged) △ Less

Submitted 14 August, 2023; v1 submitted 29 March, 2023; originally announced March 2023.

Comments: Accepted version

Journal ref: A&A 678, A209 (2023)

Orbital decay in an accreting and eclipsing 13.7 minute orbital period binary with a luminous donor

Authors: Kevin B. Burdge, Kareem El-Badry, Saul Rappaport, Tin Long Sunny Wong, Evan B. Bauer, Lars Bildsten, Ilaria Caiazzo, Deepto Chakrabarty, Emma Chickles, Matthew J. Graham, Erin Kara, S. R. Kulkarni, Thomas R. Marsh, Melania Nynka, Thomas A. Prince, Robert A. Simcoe, Jan van Roestel, Zach Vanderbosch, Eric C. Bellm, Richard G. Dekany, Andrew J. Drake, George Helou, Frank J. Masci, Jennifer Milburn, Reed Riddle , et al. (2 additional authors not shown)

Abstract: We report the discovery of ZTF J0127+5258, a compact mass-transferring binary with an orbital period of 13.7 minutes. The system contains a white dwarf accretor, which likely originated as a post-common envelope carbon-oxygen (CO) white dwarf, and a warm donor ($T_{\rm eff,\,donor}= 16,400\pm1000\,\rm K$). The donor probably formed during a common envelope phase between the CO white dwarf and an e… ▽ More We report the discovery of ZTF J0127+5258, a compact mass-transferring binary with an orbital period of 13.7 minutes. The system contains a white dwarf accretor, which likely originated as a post-common envelope carbon-oxygen (CO) white dwarf, and a warm donor ($T_{\rm eff,\,donor}= 16,400\pm1000\,\rm K$). The donor probably formed during a common envelope phase between the CO white dwarf and an evolving giant which left behind a helium star or helium white dwarf in a close orbit with the CO white dwarf. We measure gravitational wave-driven orbital inspiral with $\sim 35σ$ significance, which yields a joint constraint on the component masses and mass transfer rate. While the accretion disk in the system is dominated by ionized helium emission, the donor exhibits a mixture of hydrogen and helium absorption lines. Phase-resolved spectroscopy yields a donor radial-velocity semi-amplitude of $771\pm27\,\rm km\, s^{-1}$, and high-speed photometry reveals that the system is eclipsing. We detect a {\it Chandra} X-ray counterpart with $L_{X}\sim 3\times 10^{31}\,\rm erg\,s^{-1}$. Depending on the mass-transfer rate, the system will likely evolve into either a stably mass-transferring helium CV, merge to become an R Crb star, or explode as a Type Ia supernova in the next million years. We predict that the Laser Space Interferometer Antenna (LISA) will detect the source with a signal-to-noise ratio of $24\pm6$ after 4 years of observations. The system is the first \emph{LISA}-loud mass-transferring binary with an intrinsically luminous donor, a class of sources that provide the opportunity to leverage the synergy between optical and infrared time domain surveys, X-ray facilities, and gravitational-wave observatories to probe general relativity, accretion physics, and binary evolution. △ Less

Submitted 23 March, 2023; originally announced March 2023.

Comments: 13 pages, 7 figures, 2 tables, submitted to ApJL

Tidal Disruption Event Demographics with the Zwicky Transient Facility: Volumetric Rates, Luminosity Function, and Implications for the Local Black Hole Mass Function

Authors: Yuhan Yao, Vikram Ravi, Suvi Gezari, Sjoert van Velzen, Wenbin Lu, Steve Schulze, Jean J. Somalwar, S. R. Kulkarni, Erica Hammerstein, Matt Nicholl, Matthew J. Graham, Daniel A. Perley, S. Bradley Cenko, Robert Stein, Angelo Ricarte, Urmila Chadayammuri, Eliot Quataert, Eric C. Bellm, Joshua S. Bloom, Richard Dekany, Andrew J. Drake, Steven L. Groom, Ashish A. Mahabal, Thomas A. Prince, Reed Riddle , et al. (4 additional authors not shown)

Abstract: 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 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 sa… ▽ More 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 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 $φ(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.3}\lesssim (M_{\rm BH}/M_\odot) \lesssim 10^{7.3}$, the TDE mass function follows $φ(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. △ Less

Submitted 7 September, 2023; v1 submitted 11 March, 2023; originally announced March 2023.

Comments: Replaced following peer-review process. 38 pages, 23 figures. Accepted for publication in ApJL

Young stellar objects, accretion disks, and their variability with Rubin Observatory LSST

Authors: R. Bonito, L. Venuti, S. Ustamujic, P. Yoachim, R. A. Street, L. Prisinzano, P. Hartigan, M. G. Guarcello, K. G. Stassun, T. Giannini, E. D. Feigelson, A. Caratti o Garatti, S. Orlando, W. I. Clarkson, P. McGehee, E. C. Bellm, J. E. Gizis

Abstract: Vera C. Rubin Observatory, through the Legacy Survey of Space and Time (LSST), will allow us to derive a panchromatic view of variability in young stellar objects (YSOs) across all relevant timescales. Indeed, both short-term variability (on timescales of hours to days) and long-term variability (months to years), predominantly driven by the dynamics of accretion processes in disk-hosting YSOs, ca… ▽ More Vera C. Rubin Observatory, through the Legacy Survey of Space and Time (LSST), will allow us to derive a panchromatic view of variability in young stellar objects (YSOs) across all relevant timescales. Indeed, both short-term variability (on timescales of hours to days) and long-term variability (months to years), predominantly driven by the dynamics of accretion processes in disk-hosting YSOs, can be explored by taking advantage of the multi-band filters option available in Rubin LSST, in particular the $u,g,r,i$ filters that enable us to discriminate between photospheric stellar properties and accretion signatures. The homogeneity and depth of sky coverage that will be achieved with LSST will provide us with a unique opportunity to characterize the time evolution of disk accretion as a function of age and varying environmental conditions (e.g. field crowdedness, massive neighbors, metallicity), by targeting different star-forming regions. In this contribution to the Rubin LSST Survey Strategy Focus Issue, we discuss how implementing a dense observing cadence to explore short-term variability in YSOs represents a key complementary effort to the Wide-Fast-Deep observing mode that will be used to survey the sky over the full duration of the main survey ($\approx$10 years). The combination of these two modes will be vital to investigate the connection between the inner disk dynamics and longer-term eruptive variability behaviors, such as those observed on EXor objects. △ Less

Submitted 2 February, 2023; originally announced February 2023.

Comments: 11 pages, 4 figures, 1 table; accepted for publication in The Astrophysical Journal Supplement Series

A Systematic Study of Ia-CSM Supernovae from the ZTF Bright Transient Survey

Authors: Yashvi Sharma, Jesper Sollerman, Christoffer Fremling, Shrinivas R. Kulkarni, Kishalay De, Ido Irani, Steve Schulze, Nora Linn Strotjohann, Avishay Gal-Yam, Kate Maguire, Daniel A. Perley, Eric C. Bellm, Erik C. Kool, Thomas Brink, Rachel Bruch, Maxime Deckers, Richard Dekany, Alison Dugas, Samantha Goldwasser, Matthew J. Graham, Melissa L. Graham, Steven L. Groom, Matt Hankins, Jacob Jencson, Joel P. Johansson , et al. (13 additional authors not shown)

Abstract: Among the supernovae (SNe) that show strong interaction with the circumstellar medium, there is a rare subclass of Type Ia supernovae, SNe Ia-CSM, that show strong narrow hydrogen emission lines much like SNe IIn but on top of a diluted over-luminous Type Ia spectrum. In the only previous systematic study of this class (Silverman et al. 2013), 16 objects were identified, 8 historic and 8 from the… ▽ More Among the supernovae (SNe) that show strong interaction with the circumstellar medium, there is a rare subclass of Type Ia supernovae, SNe Ia-CSM, that show strong narrow hydrogen emission lines much like SNe IIn but on top of a diluted over-luminous Type Ia spectrum. In the only previous systematic study of this class (Silverman et al. 2013), 16 objects were identified, 8 historic and 8 from the Palomar Transient Factory (PTF). Now using the successor survey to PTF, the Zwicky Transient Facility (ZTF), we have classified 12 additional objects of this type through the systematic Bright Transient Survey (BTS). In this study, we present and analyze the optical and mid-IR light curves, optical spectra, and host galaxy properties of this sample. Consistent with previous studies, we find the objects to have slowly evolving light curves compared to normal SNe Ia with peak absolute magnitudes between -19.1 and -21, spectra having weak H$β$, large Balmer decrements of ~7 and strong Ca NIR emission. Out of 10 SNe from our sample observed by NEOWISE, 9 have $3σ$ detections, along with some showing a clear reduction in red-wing of H$α$, indicative of newly formed dust. We do not find our SN Ia-CSM sample to have a significantly different distribution of equivalent width of He I $\lambda5876$ than SNe IIn as observed in Silverman et al. 2013. The hosts tend to be late-type galaxies with recent star formation. We also derive a rate estimate of 29$^{+27}_{-21}$ Gpc$^{-3}$ yr$^{-1}$ for SNe Ia-CSM which is ~0.02--0.2 % of the SN Ia rate. This work nearly doubles the sample of well-studied Ia-CSM objects in Silverman et al. 2013, increasing the total number to 28. △ Less

Submitted 11 January, 2023; originally announced January 2023.

Comments: 26 pages, 14 figures, 5 tables. Submitted to ApJ

The prevalence and influence of circumstellar material around hydrogen-rich supernova progenitors

Authors: Rachel J. Bruch, Avishay Gal-Yam, Ofer Yaron, Ping Chen, Nora L. Strotjohann, Ido Irani, Erez Zimmerman, Steve Schulze, Yi Yang, Young-Lo Kim, Mattia Bulla, Jesper Sollerman, Mickael Rigault, Eran Ofek, Maayane Soumagnac, Frank J. Masci, Christoffer Fremling, Daniel Perley, Jakob Nordin, S. Bradley Cenko, Anna Y. Q. Ho, S. Adams, Igor Adreoni, Eric C. Bellm, Nadia Blagorodnova , et al. (22 additional authors not shown)

Abstract: Narrow transient emission lines (flash-ionization features) in early supernova (SN) spectra trace the presence of circumstellar material (CSM) around the massive progenitor stars of core-collapse SNe. The lines disappear within days after the SN explosion, suggesting that this material is spatially confined, and originates from enhanced mass loss shortly (months to a few years) prior to explosion.… ▽ More Narrow transient emission lines (flash-ionization features) in early supernova (SN) spectra trace the presence of circumstellar material (CSM) around the massive progenitor stars of core-collapse SNe. The lines disappear within days after the SN explosion, suggesting that this material is spatially confined, and originates from enhanced mass loss shortly (months to a few years) prior to explosion. We performed a systematic survey of H-rich (Type II) SNe discovered within less than two days from explosion during the first phase of the Zwicky Transient Facility (ZTF) survey (2018-2020), finding thirty events for which a first spectrum was obtained within $< 2$ days from explosion. The measured fraction of events showing flash ionisation features ($>36\%$ at $95\%$ confidence level) confirms that elevated mass loss in massive stars prior to SN explosion is common. We find that SNe II showing flash ionisation features are not significantly brighter, nor bluer, nor more slowly rising than those without. This implies that CSM interaction does not contribute significantly to their early continuum emission, and that the CSM is likely optically thin. We measured the persistence duration of flash ionisation emission and find that most SNe show flash features for $\approx 5 $ days. Rarer events, with persistence timescales $>10$ days, are brighter and rise longer, suggesting these may be intermediate between regular SNe II and strongly-interacting SNe IIn. △ Less

Submitted 13 December, 2022; v1 submitted 6 December, 2022; originally announced December 2022.

A very luminous jet from the disruption of a star by a massive black hole

Authors: Igor Andreoni, Michael W. Coughlin, Daniel A. Perley, Yuhan Yao, Wenbin Lu, S. Bradley Cenko, Harsh Kumar, Shreya Anand, Anna Y. Q. Ho, Mansi M. Kasliwal, Antonio de Ugarte Postigo, Ana Sagues-Carracedo, Steve Schulze, D. Alexander Kann, S. R. Kulkarni, Jesper Sollerman, Nial Tanvir, Armin Rest, Luca Izzo, Jean J. Somalwar, David L. Kaplan, Tomas Ahumada, G. C. Anupama, Katie Auchettl, Sudhanshu Barway , et al. (56 additional authors not shown)

Abstract: Tidal disruption events (TDEs) are bursts of electromagnetic energy released when supermassive black holes (SMBHs) at the centers of galaxies violently disrupt a star that passes too close. TDEs provide a new window to study accretion onto SMBHs; in some rare cases, this accretion leads to launching of a relativistic jet, but the necessary conditions are not fully understood. The best studied jett… ▽ More Tidal disruption events (TDEs) are bursts of electromagnetic energy released when supermassive black holes (SMBHs) at the centers of galaxies violently disrupt a star that passes too close. TDEs provide a new window to study accretion onto SMBHs; in some rare cases, this accretion leads to launching of a relativistic jet, but the necessary conditions are not fully understood. The best studied jetted TDE to date is Swift J1644+57, which was discovered in gamma-rays, but was too obscured by dust to be seen at optical wavelengths. Here we report the optical discovery of AT2022cmc, a rapidly fading source at cosmological distance (redshift z=1.19325) whose unique lightcurve transitioned into a luminous plateau within days. Observations of a bright counterpart at other wavelengths, including X-rays, sub-millimeter, and radio, supports the interpretation of AT2022cmc as a jetted TDE containing a synchrotron "afterglow", likely launched by a SMBH with spin $a \gtrsim 0.3$. Using 4 years of Zwicky Transient Facility (ZTF) survey data, we calculate a rate of $0.02 ^{+ 0.04 }_{- 0.01 }$ Gpc$^{-3}$ yr$^{-1}$ for on-axis jetted TDEs based on the luminous, fast-fading red component, thus providing a measurement complementary to the rates derived from X-ray and radio observations. Correcting for the beaming angle effects, this rate confirms that about 1% of TDEs have relativistic jets. Optical surveys can use AT2022cmc as a prototype to unveil a population of jetted TDEs. △ Less

Submitted 29 November, 2022; originally announced November 2022.

Comments: Accepted for publication in Nature

Deep Drilling in the Time Domain with DECam: Survey Characterization

Authors: Melissa L. Graham, Robert A. Knop, Thomas Kennedy, Peter E. Nugent, Eric Bellm, Márcio Catelan, Avi Patel, Hayden Smotherman, Monika Soraisam, Steven Stetzler, Lauren N. Aldoroty, Autumn Awbrey, Karina Baeza-Villagra, Pedro H. Bernardinelli, Federica Bianco, Dillon Brout, Riley Clarke, William I. Clarkson, Thomas Collett, James R. A. Davenport, Shenming Fu, John E. Gizis, Ari Heinze, Lei Hu, Saurabh W. Jha , et al. (19 additional authors not shown)

Abstract: This paper presents a new optical imaging survey of four deep drilling fields (DDFs), two Galactic and two extragalactic, with the Dark Energy Camera (DECam) on the 4 meter Blanco telescope at the Cerro Tololo Inter-American Observatory (CTIO). During the first year of observations in 2021, $>$4000 images covering 21 square degrees (7 DECam pointings), with $\sim$40 epochs (nights) per field and 5… ▽ More This paper presents a new optical imaging survey of four deep drilling fields (DDFs), two Galactic and two extragalactic, with the Dark Energy Camera (DECam) on the 4 meter Blanco telescope at the Cerro Tololo Inter-American Observatory (CTIO). During the first year of observations in 2021, $>$4000 images covering 21 square degrees (7 DECam pointings), with $\sim$40 epochs (nights) per field and 5 to 6 images per night per filter in $g$, $r$, $i$, and/or $z$, have become publicly available (the proprietary period for this program is waived). We describe the real-time difference-image pipeline and how alerts are distributed to brokers via the same distribution system as the Zwicky Transient Facility (ZTF). In this paper, we focus on the two extragalactic deep fields (COSMOS and ELAIS-S1), characterizing the detected sources and demonstrating that the survey design is effective for probing the discovery space of faint and fast variable and transient sources. We describe and make publicly available 4413 calibrated light curves based on difference-image detection photometry of transients and variables in the extragalactic fields. We also present preliminary scientific analysis regarding Solar System small bodies, stellar flares and variables, Galactic anomaly detection, fast-rising transients and variables, supernovae, and active galactic nuclei. △ Less

Submitted 16 November, 2022; originally announced November 2022.

Comments: 22 pages, 17 figures, 2 tables. Accepted to MNRAS

Uncovering a population of gravitational lens galaxies with magnified standard candle SN Zwicky

Authors: Ariel Goobar, Joel Johansson, Steve Schulze, Nikki Arendse, Ana Sagués Carracedo, Suhail Dhawan, Edvard Mörtsell, Christoffer Fremling, Lin Yan, Daniel Perley, Jesper Sollerman, Rémy Joseph, K-Ryan Hinds, William Meynardie, Igor Andreoni, Eric Bellm, Josh Bloom, Thomas E. Collett, Andrew Drake, Matthew Graham, Mansi Kasliwal, Shri Kulkarni, Cameron Lemon, Adam A. Miller, James D. Neill , et al. (13 additional authors not shown)

Abstract: Detecting gravitationally lensed supernovae is among the biggest challenges in astronomy. It involves a combination of two very rare phenomena: catching the transient signal of a stellar explosion in a distant galaxy and observing it through a nearly perfectly aligned foreground galaxy that deflects light towards the observer. High-cadence optical observations with the Zwicky Transient Facility, w… ▽ More Detecting gravitationally lensed supernovae is among the biggest challenges in astronomy. It involves a combination of two very rare phenomena: catching the transient signal of a stellar explosion in a distant galaxy and observing it through a nearly perfectly aligned foreground galaxy that deflects light towards the observer. High-cadence optical observations with the Zwicky Transient Facility, with an unparalleled large field of view, led to the detection of a multiply-imaged Type Ia supernova (SN Ia), ``SN Zwicky", a.k.a. SN 2022qmx. Magnified nearly twenty-five times, the system was found thanks to the ``standard candle" nature of SNe Ia. High-spatial resolution imaging with the Keck telescope resolved four images of the supernova with very small angular separation, corresponding to an Einstein radius of only $θ_E =0.167"$ and almost identical arrival times. The small $θ_E$ and faintness of the lensing galaxy is very unusual, highlighting the importance of supernovae to fully characterise the properties of galaxy-scale gravitational lenses, including the impact of galaxy substructures. △ Less

Submitted 14 June, 2023; v1 submitted 1 November, 2022; originally announced November 2022.

Comments: Matches published version in Nature Astronomy

A search for relativistic ejecta in a sample of ZTF broad-lined Type Ic supernovae

Authors: Alessandra Corsi, Anna Y. Q. Ho, S. Bradley Cenko, Shrinivas R. Kulkarni, Shreya Anand, Sheng Yang, Jesper Sollerman, Gokul P. Srinivasaragavan, Conor M. B. Omand, Arvind Balasubramanian, Dale A. Frail, Christoffer Fremling, Daniel A. Perley, Yuhan Yao, Aishwarya S. Dahiwale, Kishalay De, Alison Dugas, Matthew Hankins, Jacob Jencson, Mansi M. Kasliwal, Anastasios Tzanidakis, Eric C. Bellm, Russ R. Laher, Frank J. Masci, Josiah N. Purdum , et al. (1 additional authors not shown)

Abstract: 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 r… ▽ More 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$σ$ 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. △ Less

Submitted 17 October, 2022; originally announced October 2022.

Comments: 28 pages, 11 figures, submitted to AAS journal

A radio-detected Type Ia supernova with helium-rich circumstellar material

Authors: Erik C. Kool, Joel Johansson, Jesper Sollerman, Javier Moldón, Takashi J. Moriya, Steve Schulze, Laura Chomiuk, Chelsea Harris, Miguel Pérez-Torres, Seppo Mattila, Peter Lundqvist, Matthew Graham, Sheng Yang, Daniel A. Perley, Nora Linn Strotjohann, Christoffer Fremling, Avishay Gal-Yam, Jeremy Lezmy, Kate Maguire, Conor Omand, Mathew Smith, Igor Andreoni, Eric C. Bellm, Kishalay De, Joshua S. Bloom , et al. (12 additional authors not shown)

Abstract: Type Ia supernovae (SNe Ia) are thermonuclear explosions of degenerate white dwarf (WD) stars destabilized by mass accretion from a companion star, but the nature of their progenitors remains poorly understood. A way to discriminate between progenitor systems is through radio observations; a non-degenerate companion star is expected to lose material through winds or binary interaction prior to exp… ▽ More Type Ia supernovae (SNe Ia) are thermonuclear explosions of degenerate white dwarf (WD) stars destabilized by mass accretion from a companion star, but the nature of their progenitors remains poorly understood. A way to discriminate between progenitor systems is through radio observations; a non-degenerate companion star is expected to lose material through winds or binary interaction prior to explosion, and the SN ejecta crashing into this nearby circumstellar material (CSM) should result in radio synchrotron emission. However, despite extensive efforts, no SN Ia has ever been detected at radio wavelengths, which suggests a clean environment and a companion star that is itself a degenerate WD star. Here we report on the study of SN 2020eyj, a SN Ia showing helium-rich CSM, as revealed by its spectral features, infrared emission and, for the first time in a SN Ia, a radio counterpart. Based on our modeling, we conclude the CSM likely originates from a single-degenerate (SD) binary system where a WD accretes material from a helium donor star, an often hypothesized formation channel for SNe Ia. We describe how comprehensive radio follow-up of SN 2020eyj-like SNe Ia can improve the constraints on their progenitor systems. △ Less

Submitted 17 May, 2023; v1 submitted 14 October, 2022; originally announced October 2022.

Comments: 62 pages, 10 figures, 4 tables

Journal ref: Nature 617 (2023), 477

A dense $\mathbf{0.1 M_{\rm \odot}}$ star in a 51-minute orbital period eclipsing binary

Authors: Kevin B. Burdge, Kareem El-Badry, Thomas R. Marsh, Saul Rappaport, Warren R. Brown, Ilaria Caiazzo, Deepto Chakrabarty, V. S. Dhillon, Jim Fuller, Boris T. Gänsicke, Matthew J. Graham, Erin Kara, S. R. Kulkarni, S. P. Littlefair, Przemek Mróz, Pablo Rodríguez-Gil, Jan van Roestel, Robert A. Simcoe, Eric C. Bellm, Andrew J. Drake, Richard G. Dekany, Steven L. Groom, Russ R. Laher, Frank J. Masci, Reed Riddle , et al. (2 additional authors not shown)

Abstract: In over a thousand known cataclysmic variables (CVs), where a white dwarf is accreting from a hydrogen-rich star, only a dozen have orbital periods below 75 minutes. One way to achieve these short periods requires the donor star to have undergone substantial nuclear evolution prior to interacting with the white dwarf, and it is expected that these objects will transition to helium accretion. These… ▽ More In over a thousand known cataclysmic variables (CVs), where a white dwarf is accreting from a hydrogen-rich star, only a dozen have orbital periods below 75 minutes. One way to achieve these short periods requires the donor star to have undergone substantial nuclear evolution prior to interacting with the white dwarf, and it is expected that these objects will transition to helium accretion. These transitional CVs have been proposed as progenitors of helium CVs. However, no known transitional CV is expected to reach an orbital period short enough to account for most of the helium CV population, leaving the role of this evolutionary pathway unclear. Here we report observations of ZTF J1813+4251, a 51-minute orbital period, fully eclipsing binary system consisting of a star with a temperature comparable to that of the Sun but a density 100 times greater due to its helium-rich composition, accreting onto a white dwarf. Phase-resolved spectra, multi-band light curves and the broadband spectral energy distribution allow us to obtain precise and robust constraints on the masses, radii and temperatures of both components. Evolutionary modeling shows that ZTF J1813+4251 is destined to become a helium CV binary, reaching an orbital period under 20 minutes, rendering ZTF J1813+4251 a previously missing link between helium CV binaries and hydrogen-rich CVs. △ Less

Submitted 4 October, 2022; originally announced October 2022.

Comments: 48 Pages, 12 figures, 2 tables, Published online by Nature on Oct 5, 2022

Journal ref: Nature 610 467-471 (2022)

A light in the dark: searching for electromagnetic counterparts to black hole-black hole mergers in LIGO/Virgo O3 with the Zwicky Transient Facility

Authors: Matthew J. Graham, Barry McKernan, K. E. Saavik Ford, Daniel Stern, S. G. Djorgovski, Michael Coughlin, Kevin B. Burdge, Eric C. Bellm, George Helou, Ashish A. Mahabal, Frank J. Masci, Josiah Purdum, Philippe Rosnet, Ben Rusholme

Abstract: The accretion disks of active galactic nuclei (AGN) are promising locations for the merger of compact objects detected by gravitational wave (GW) observatories. Embedded within a baryon-rich, high density environment, mergers within AGN are the only GW channel where an electromagnetic (EM) counterpart must occur (whether detectable or not). Considering AGN with unusual flaring activity observed by… ▽ More The accretion disks of active galactic nuclei (AGN) are promising locations for the merger of compact objects detected by gravitational wave (GW) observatories. Embedded within a baryon-rich, high density environment, mergers within AGN are the only GW channel where an electromagnetic (EM) counterpart must occur (whether detectable or not). Considering AGN with unusual flaring activity observed by the Zwicky Transient Facility (ZTF), we describe a search for candidate EM counterparts to binary black hole (BBH) mergers detected by LIGO/Virgo in O3. After removing probable false positives, we find nine candidate counterparts to BBH mergers mergers during O3 (seven in O3a, two in O3b) with a $p$-value of 0.019. Based on ZTF sky coverage, AGN geometry, and merger geometry, we expect $\approx 3(N_{\rm BBH}/83)(f_{\rm AGN}/0.5)$ potentially detectable EM counterparts from O3, where $N_{\rm BBH}$ is the total number of observed BBH mergers and $f_{\rm AGN}$ is the fraction originating in AGN. Further modeling of breakout and flaring phenomena in AGN disks is required to reduce our false positive rate. Two of the events are also associated with mergers with total masses $> 100M_\odot$, which is the expected rate for O3 if hierarchical (large mass) mergers occur in the AGN channel. Candidate EM counterparts in future GW observing runs can be better constrained by coverage of the Southern sky as well as spectral monitoring of unusual AGN flaring events in LIGO/Virgo alert volumes. A future set of reliable AGN EM counterparts to BBH mergers will yield an independent means of measuring cosmic expansion ($H_0$) as a function of redshift. △ Less

Submitted 26 September, 2022; originally announced September 2022.

Comments: 24 pages, 6 figures, submitted to ApJ

Rubin Observatory LSST Transients and Variable Stars Roadmap

Authors: Kelly M. Hambleton, Federica B. Bianco, Rachel Street, Keaton Bell, David Buckley, Melissa Graham, Nina Hernitschek, Michael B. Lund, Elena Mason, Joshua Pepper, Andrej Prsa, Markus Rabus, Claudia M. Raiteri, Robert Szabo, Paula Szkody, Igor Andreoni, Simone Antoniucci, Barbara Balmaverde, Eric Bellm, Rosaria Bonito, Giuseppe Bono, Maria Teresa Botticella, Enzo Brocato, Katja Bucar Bricman, Enrico Cappellaro , et al. (57 additional authors not shown)

Abstract: The Vera C. Rubin Legacy Survey of Space and Time holds the potential to revolutionize time domain astrophysics, reaching completely unexplored areas of the Universe and mapping variability time scales from minutes to a decade. To prepare to maximize the potential of the Rubin LSST data for the exploration of the transient and variable Universe, one of the four pillars of Rubin LSST science, the T… ▽ More The Vera C. Rubin Legacy Survey of Space and Time holds the potential to revolutionize time domain astrophysics, reaching completely unexplored areas of the Universe and mapping variability time scales from minutes to a decade. To prepare to maximize the potential of the Rubin LSST data for the exploration of the transient and variable Universe, one of the four pillars of Rubin LSST science, the Transient and Variable Stars Science Collaboration, one of the eight Rubin LSST Science Collaborations, has identified research areas of interest and requirements, and paths to enable them. While our roadmap is ever-evolving, this document represents a snapshot of our plans and preparatory work in the final years and months leading up to the survey's first light. △ Less

Submitted 8 August, 2022; originally announced August 2022.

Comments: 202 pages (in book format) 34 figures plus chapter heading figures (13)

Prospects of Gravitational Wave Follow-up Through a Wide-field Ultra-violet Satellite: a Dorado Case Study

Authors: Bas Dorsman, Geert Raaijmakers, S. Bradley Cenko, Samaya Nissanke, Leo P. Singer, Mansi M. Kasliwal, Anthony L. Piro, Eric C. Bellm, Dieter H. Hartmann, Kenta Hotokezaka, Kamilė Lukošiūtė

Abstract: The detection of gravitational waves from binary neuron star merger GW170817 and electromagnetic counterparts GRB170817 and AT2017gfo kick-started the field of gravitational wave multimessenger astronomy. The optically red to near infra-red emission (`red' component) of AT2017gfo was readily explained as produced by the decay of newly created nuclei produced by rapid neutron capture (a kilonova).… ▽ More The detection of gravitational waves from binary neuron star merger GW170817 and electromagnetic counterparts GRB170817 and AT2017gfo kick-started the field of gravitational wave multimessenger astronomy. The optically red to near infra-red emission (`red' component) of AT2017gfo was readily explained as produced by the decay of newly created nuclei produced by rapid neutron capture (a kilonova). However, the ultra-violet to optically blue emission (`blue' component) that was dominant at early times (up to 1.5 days) received no consensus regarding its driving physics. Among many explanations, two leading contenders are kilonova radiation from a lanthanide-poor ejecta component or shock interaction (cocoon emission). In this work, we simulate AT2017gfo-like light curves and perform a Bayesian analysis to study whether an ultra-violet satellite capable of rapid gravitational wave follow-up, could distinguish between physical processes driving the early `blue' component. We find that a Dorado-like ultra-violet satellite, with a 50 sq. deg. field of view and a limiting magnitude (AB) of 20.5 for a 10 minute exposure is able to distinguish radiation components up to at least 160 Mpc if data collection starts within 3.2 or 5.2 hours for two possible AT2017gfo-like light curve scenarios. We also study the degree to which parameters can be constrained with the obtained photometry. We find that, while ultra-violet data alone constrains parameters governing the outer ejecta properties, the combination of both ground-based optical and space-based ultra-violet data allows for tight constraints for all but one parameter of the kilonova model up to 160 Mpc. These results imply that an ultra-violet mission like Dorado would provide unique insights into the early evolution of the post-merger system and its driving emission physics. △ Less

Submitted 20 June, 2022; originally announced June 2022.

Comments: 13 pages, 2 tables, 7 figures. Comments welcome. To be submitted to ApJ

Discovery of Two Polars from a Crossmatch of ZTF and the SRG/eFEDS X-ray Catalog

Authors: Antonio C. Rodriguez, Shrinivas R. Kulkarni, Thomas A. Prince, Paula Szkody, Kevin B. Burdge, Ilaria Caiazzo, Jan van Roestel, Zachary P. Vanderbosch, Kareem El-Badry, Eric C. Bellm, Boris T. Gänsicke, Matthew J. Graham, Ashish A. Mahabal, Frank J. Masci, Przemek Mróz, Reed Riddle, Ben Rusholme

Abstract: Magnetic CVs are luminous Galactic X-ray sources but have been difficult to find in purely optical surveys due to their lack of outburst behavior. The eROSITA telescope aboard the Spektr-RG (SRG) mission is conducting an all-sky X-ray survey and recently released the public eROSITA Final Equatorial Depth Survey (eFEDS) catalog. We crossmatched the eFEDS catalog with photometry from the Zwicky Tran… ▽ More Magnetic CVs are luminous Galactic X-ray sources but have been difficult to find in purely optical surveys due to their lack of outburst behavior. The eROSITA telescope aboard the Spektr-RG (SRG) mission is conducting an all-sky X-ray survey and recently released the public eROSITA Final Equatorial Depth Survey (eFEDS) catalog. We crossmatched the eFEDS catalog with photometry from the Zwicky Transient Facility (ZTF) and discovered two new magnetic cataclysmic variables (CVs). We obtained high-cadence optical photometry and phase-resolved spectroscopy for each magnetic CV candidate and found them both to be polars. Among the newly discovered magnetic CVs is ZTFJ0850+0443, an eclipsing polar with orbital period $P_\textrm{orb} = 1.72$ hr, white dwarf mass $M_\textrm{WD} = 0.81 \pm 0.08 M_\odot$ and accretion rate $\dot{M} \sim 10^{-11} M_\odot$/yr. We suggest that ZTFJ0850+0443 is a low magnetic field strength polar, with $B_\textrm{WD} \lesssim 10$ MG. We also discovered a non-eclipsing polar, ZTFJ0926+0105, with orbital period $P_\textrm{orb} = 1.48$ hr, magnetic field strength $B_\textrm{WD} \gtrsim 26$ MG, and accretion rate $\dot{M} \sim 10^{-12} M_\odot$/yr. △ Less

Submitted 9 June, 2022; originally announced June 2022.

Comments: Submitted to ApJ

A 62-minute orbital period black widow binary in a wide hierarchical triple

Authors: Kevin B. Burdge, Thomas R. Marsh, Jim Fuller, Eric C. Bellm, Ilaria Caiazzo, Deepto Chakrabarty, Michael W. Coughlin, Kishalay De, V. S. Dhillon, Matthew J. Graham, Pablo Rodrí guez-Gil, Amruta D. Jaodand, David L. Kaplan, Erin Kara, Albert K. H. Kong, S. R. Kulkarni, Kwan-Lok Li, S. P. Littlefair, Walid A. Majid, Przemek Mróz, Aaron B. Pearlman, E. S. Phinney, Jan van Roestel, Robert A. Simcoe, Igor Andreoni , et al. (8 additional authors not shown)

Abstract: Over a dozen millisecond pulsars are ablating low-mass companions in close binary systems. In the original "black widow", the 8-hour orbital period eclipsing pulsar PSR J1959+2048 (PSR B1957+20), high energy emission originating from the pulsar is irradiating and may eventually destroy a low-mass companion. These systems are not only physical laboratories that reveal the dramatic result of exposin… ▽ More Over a dozen millisecond pulsars are ablating low-mass companions in close binary systems. In the original "black widow", the 8-hour orbital period eclipsing pulsar PSR J1959+2048 (PSR B1957+20), high energy emission originating from the pulsar is irradiating and may eventually destroy a low-mass companion. These systems are not only physical laboratories that reveal the dramatic result of exposing a close companion star to the relativistic energy output of a pulsar, but are also believed to harbour some of the most massive neutron stars, allowing for robust tests of the neutron star equation of state. Here, we report observations of ZTF J1406+1222, a wide hierarchical triple hosting a 62-minute orbital period black widow candidate whose optical flux varies by a factor of more than 10. ZTF J1406+1222 pushes the boundaries of evolutionary models, falling below the 80 minute minimum orbital period of hydrogen-rich systems. The wide tertiary companion is a rare low metallicity cool subdwarf star, and the system has a Galactic halo orbit consistent with passing near the Galactic center, making it a probe of formation channels, neutron star kick physics, and binary evolution. △ Less

Submitted 4 May, 2022; originally announced May 2022.

Comments: 63 pages, 15 figures, 2 tables, published in Nature on May 4, 2022

Journal ref: Nature 605, 41-45 (2022)

Neutrino follow-up with the Zwicky Transient Facility: Results from the first 24 campaigns

Authors: Robert Stein, Simeon Reusch, Anna Franckowiak, Marek Kowalski, Jannis Necker, Sven Weimann, Mansi M. Kasliwal, Jesper Sollerman, Tomas Ahumada, Pau Amaro-Seoane, Shreya Anand, Igor Andreoni, Eric C. Bellm, Joshua S. Bloom, Michael Coughlin, Kishalay De, Christoffer Fremling, Suvi Gezari, Matthew Graham, Steven L. Groom, George Helou, David L. Kaplan, Viraj Karambelkar, Albert K. H. Kong, Erik C. Kool , et al. (11 additional authors not shown)

Abstract: The Zwicky Transient Facility (ZTF) performs a systematic neutrino follow-up program, searching for optical counterparts to high-energy neutrinos with dedicated Target-of-Opportunity (ToO) observations. Since first light in March 2018, ZTF has taken prompt observations for 24 high-quality neutrino alerts from the IceCube Neutrino Observatory, with a median latency of 12.2 hours from initial neutri… ▽ More The Zwicky Transient Facility (ZTF) performs a systematic neutrino follow-up program, searching for optical counterparts to high-energy neutrinos with dedicated Target-of-Opportunity (ToO) observations. Since first light in March 2018, ZTF has taken prompt observations for 24 high-quality neutrino alerts from the IceCube Neutrino Observatory, with a median latency of 12.2 hours from initial neutrino detection. From two of these campaigns, we have already reported tidal disruption event (TDE) AT 2019dsg and likely TDE AT 2019fdr as probable counterparts, suggesting that TDEs contribute >7.8% of the astrophysical neutrino flux. We here present the full results of our program through to December 2021. No additional candidate neutrino sources were identified by our program, allowing us to place the first constraints on the underlying optical luminosity function of astrophysical neutrino sources. Transients with optical absolutes magnitudes brighter that $-21$ can contribute no more than 87% of the total, while transients brighter than $-22$ can contribute no more than 58% of the total, neglecting the effect of extinction and assuming they follow the star formation rate. These are the first observational constraints on the neutrino emission of bright populations such as superluminous supernovae. None of the neutrinos were coincident with bright optical AGN flares comparable to that observed for TXS 0506+056/IC170922A, with such optical blazar flares producing no more than 26% of the total neutrino flux. We highlight the outlook for electromagnetic neutrino follow-up programs, including the expected potential for the Rubin Observatory. △ Less

Submitted 4 April, 2024; v1 submitted 31 March, 2022; originally announced March 2022.

Journal ref: MNRAS, 521, 4, 2023, 5046-5063

In search of short gamma-ray burst optical counterpart with the Zwicky Transient Facility

Authors: Tomás Ahumada, Shreya Anand, Michael W. Coughlin, Igor Andreoni, Erik C. Kool, Harsh Kumar, Simeon Reusch, Ana Sagués-Carracedo, Robert Stein, S. Bradley Cenko, Mansi M. Kasliwal, Leo P. Singer, Rachel Dunwoody, Joseph Mangan, Varun Bhalerao, Mattia Bulla, Eric Burns, Matthew J. Graham, David L. Kaplan, Daniel Perley, Mouza Almualla, Joshua S. Bloom, Virginia Cunningham, Kishalay De, Pradip Gatkine , et al. (24 additional authors not shown)

Abstract: The Fermi Gamma-ray Burst Monitor (GBM) triggers on-board in response to $\sim$ 40 short gamma-ray bursts (SGRBs) per year; however, their large localization regions have made the search for optical counterparts a challenging endeavour. We have developed and executed an extensive program with the wide field of view of the Zwicky Transient Facility (ZTF) camera, mounted on the Palomar 48 inch Oschi… ▽ More The Fermi Gamma-ray Burst Monitor (GBM) triggers on-board in response to $\sim$ 40 short gamma-ray bursts (SGRBs) per year; however, their large localization regions have made the search for optical counterparts a challenging endeavour. We have developed and executed an extensive program with the wide field of view of the Zwicky Transient Facility (ZTF) camera, mounted on the Palomar 48 inch Oschin telescope (P48), to perform target-of-opportunity (ToO) observations on 10 Fermi-GBM SGRBs during 2018 and 2020-2021. Bridging the large sky areas with small field of view optical telescopes in order to track the evolution of potential candidates, we look for the elusive SGRB afterglows and kilonovae (KNe) associated with these high-energy events. No counterpart has yet been found, even though more than 10 ground based telescopes, part of the Global Relay of Observatories Watching Transients Happen (GROWTH) network, have taken part in these efforts. The candidate selection procedure and the follow-up strategy have shown that ZTF is an efficient instrument for searching for poorly localized SGRBs, retrieving a reasonable number of candidates to follow-up and showing promising capabilities as the community approaches the multi-messenger era. Based on the median limiting magnitude of ZTF, our searches would have been able to retrieve a GW170817-like event up to $\sim$ 200 Mpc and SGRB afterglows to z = 0.16 or 0.4, depending on the assumed underlying energy model. Future ToOs will expand the horizon to z = 0.2 and 0.7 respectively. △ Less

Submitted 22 March, 2022; originally announced March 2022.

A Uniform Type Ia Supernova Distance Ladder with the Zwicky Transient Facility: Absolute Calibration Based on the Tip of the Red Giant Branch (TRGB) Method

Authors: Suhail Dhawan, Ariel Goobar, Joel Johansson, In Sung Jang, Mickael Rigault, Luke Harvey, Kate Maguire, Wendy L. Freedman, Barry F. Madore, Mathew Smith, Jesper Sollerman, Young-Lo Kim, Igor Andreoni, Eric C. Bellm, Michael W. Coughlin, R. Dekany, Matthew J. Graham, Shrinivas R. Kulkarni, Russ R. Laher, Michael S. Medford, James D. Neill, Guy Nir, Reed Riddle, Ben Rusholme

Abstract: The current Cepheid-calibrated distance ladder measurement of $H_0$ is reported to be in tension with the values inferred from the cosmic microwave background (CMB), assuming standard cosmology. However, some tip of the red giant branch (TRGB) estimates report $H_0$ in better agreement with the CMB. Hence, it is critical to reduce systematic uncertainties in local measurements to understand the Hu… ▽ More The current Cepheid-calibrated distance ladder measurement of $H_0$ is reported to be in tension with the values inferred from the cosmic microwave background (CMB), assuming standard cosmology. However, some tip of the red giant branch (TRGB) estimates report $H_0$ in better agreement with the CMB. Hence, it is critical to reduce systematic uncertainties in local measurements to understand the Hubble tension. In this paper, we propose a uniform distance ladder between the second and third rungs, combining SNe~Ia observed by the Zwicky Transient Facility (ZTF) with a TRGB calibration of their absolute luminosity. A large, volume-limited sample of both calibrator and Hubble flow SNe~Ia from the \emph{same} survey minimizes two of the largest sources of systematics: host-galaxy bias and non-uniform photometric calibration. We present results from a pilot study using existing TRGB distance to the host galaxy of ZTF SN~Ia SN 2021rhu (aka ZTF21abiuvdk) in NGC7814. Combining the ZTF calibrator with a volume-limited sample from the first data release of ZTF Hubble flow SNe~Ia, we infer $H_0 = 76.94 \pm 6.4\, {\rm km}\,{\rm s^{-1}}\,{\rm Mpc^{-1}}$, an $8.3 \%$ measurement. The error budget is dominated by the single object calibrating the SN~Ia luminosity in this pilot study. However, the ZTF sample includes already five other SNe~Ia within $\sim$ 20 Mpc for which TRGB distances can be obtained with HST. Finally, we present the prospects of building this distance ladder out to 80 Mpc with JWST observations of more than one hundred ZTF SNe~Ia. △ Less

Submitted 11 July, 2022; v1 submitted 8 March, 2022; originally announced March 2022.

Comments: 10 pages (incl. references), 4 Figures, 2 Tables, accepted to ApJ

The Final Season Reimagined: 30 Tidal Disruption Events from the ZTF-I Survey

Authors: Erica Hammerstein, Sjoert van Velzen, Suvi Gezari, S. Bradley Cenko, Yuhan Yao, Charlotte Ward, Sara Frederick, Natalia Villanueva, Jean J. Somalwar, Matthew J. Graham, Shrinivas R. Kulkarni, Daniel Stern, Igor Andreoni, Eric C. Bellm, Richard Dekany, Suhail Dhawan, Andrew J. Drake, Christoffer Fremling, Pradip Gatkine, Steven L. Groom, Anna Y. Q. Ho, Mansi M. Kasliwal, Viraj Karambelkar, Erik C. Kool, Frank J. Masci , et al. (8 additional authors not shown)

Abstract: Tidal disruption events (TDEs) offer a unique way to study dormant black holes. While the number of observed TDEs has grown thanks to the emergence of wide-field surveys in the past few decades, questions regarding the nature of the observed optical, UV, and X-ray emission remain. We present a uniformly selected sample of 30 spectroscopically classified TDEs from the Zwicky Transient Facility Phas… ▽ More Tidal disruption events (TDEs) offer a unique way to study dormant black holes. While the number of observed TDEs has grown thanks to the emergence of wide-field surveys in the past few decades, questions regarding the nature of the observed optical, UV, and X-ray emission remain. We present a uniformly selected sample of 30 spectroscopically classified TDEs from the Zwicky Transient Facility Phase I survey operations with follow-up \textit{Swift} UV and X-ray observations. Through our investigation into correlations between light curve properties, we recover a shallow positive correlation between the peak bolometric luminosity and decay timescales. We introduce a new spectroscopic class of TDE, TDE-featureless, which are characterized by featureless optical spectra. The new TDE-featureless class shows larger peak bolometric luminosities, peak blackbody temperatures, and peak blackbody radii. We examine the differences between the X-ray bright and X-ray faint populations of TDEs in this sample, finding that X-ray bright TDEs show higher peak blackbody luminosities than the X-ray faint sub-sample. This sample of optically selected TDEs is the largest sample of TDEs from a single survey yet, and the systematic discovery, classification, and follow-up of this sample allows for robust characterization of TDE properties, an important stepping stone looking forward toward the Rubin era. △ Less

Submitted 28 November, 2022; v1 submitted 2 March, 2022; originally announced March 2022.

Comments: 44 pages, 22 figures, 9 tables, accepted to ApJ

Constraining Type Ia supernova explosions and early flux excesses with the Zwicky Transient Factory

Authors: M. Deckers, K. Maguire, M. R. Magee, G. Dimitriadis, M. Smith, A. Sainz de Murieta, A. A. Miller, A. Goobar, J. Nordin, M. Rigault, E. Bellm, M. W. Coughlin, R. R. Laher, D. Shupe, M. J. Graham, M. M. Kasliwal, R. Walters

Abstract: In the new era of time-domain surveys Type Ia supernovae are being caught sooner after explosion, which has exposed significant variation in their early light curves. Two driving factors for early time evolution are the distribution of nickel in the ejecta and the presence of flux excesses of various causes. We perform an analysis of the largest young SN Ia sample to date. We compare 115 SN Ia lig… ▽ More In the new era of time-domain surveys Type Ia supernovae are being caught sooner after explosion, which has exposed significant variation in their early light curves. Two driving factors for early time evolution are the distribution of nickel in the ejecta and the presence of flux excesses of various causes. We perform an analysis of the largest young SN Ia sample to date. We compare 115 SN Ia light curves from the Zwicky Transient Facility to the turtls model grid containing light curves of Chandrasekhar-mass explosions with a range of nickel masses, nickel distributions and explosion energies. We find that the majority of our observed light curves are well reproduced by Chandrasekhar-mass explosion models with a preference for highly extended nickel distributions. We identify six SNe Ia with an early-time flux excess in our g- and r-band data (four `blue' and two `red' flux excesses). We find an intrinsic rate of 18+/-11 per cent of early flux excesses in SNe Ia at z < 0.07, based on three detected flux excesses out of 30 (10 per cent) observed SNe Ia with a simulated efficiency of 57 per cent. This is comparable to rates of flux excesses in the literature but also accounts for detection efficiencies. Two of these events are mostly consistent with CSM interaction, while the other four have longer lifetimes in agreement with companion interaction and nickel-clump models. We find a higher frequency of flux excesses in 91T/99aa-like events (44+/-13 per cent). △ Less

Submitted 25 February, 2022; originally announced February 2022.

The Hydrogen-Poor Superluminous Supernovae from the Zwicky Transient Facility Phase-I Survey: II. Light Curve Modeling and Characterization of Undulations

Authors: Z. H. Chen, Lin Yan, T. Kangas, R. Lunnan, J. Sollerman, S. Schulze, D. A. Perley, T. -W. Chen, K. Taggart, K. R. Hinds, A. Gal-Yam, X. F. Wang, K. De, E. Bellm, J. S. Bloom, R. Dekany, M. Graham, M. Kasliwal, S. Kulkarni, R. Laher, D. Neill, B. Rusholme

Abstract: We present analysis of the light curves (LCs) of 77 hydrogen-poor superluminous supernovae (SLSNe-I) discovered during the Zwicky Transient Facility Phase-I operation. We find that the majority (67\%) of the sample can be fit equally well by both magnetar and ejecta-circumstellar medium (CSM) interaction plus $^{56}$Ni decay models. This implies that LCs alone can not unambiguously constrain the p… ▽ More We present analysis of the light curves (LCs) of 77 hydrogen-poor superluminous supernovae (SLSNe-I) discovered during the Zwicky Transient Facility Phase-I operation. We find that the majority (67\%) of the sample can be fit equally well by both magnetar and ejecta-circumstellar medium (CSM) interaction plus $^{56}$Ni decay models. This implies that LCs alone can not unambiguously constrain the physical power sources for a SLSN-I. However, 23\% of the sample show inverted V-shape, steep declining LCs or features of long rise and fast post-peak decay, which are better described by the CSM+Ni model. The remaining 10\% of the sample favor the magnetar model. Moreover, our analysis shows that the LC undulations are quite common, with a fraction of $18-44\% $ in our gold sample. Among those strongly undulating events, about 62\% of them are found to be CSM-favored, implying that the undulations tend to occur in the CSM-favored events. Undulations show a wide range in energy and duration, with median values (and 1$σ$ errors) being as $1.7\%^{+1.5\%}_{-0.7\%}\,\rm E_{\rm rad,total}$ and $28.8^{+14.4}_{-9.1}$\,days, respectively. Our analysis of the undulation time scales suggests that intrinsic temporal variations of the central engine can explain half of the undulating events, while CSM interaction can account for the majority of the sample. Finally, all of the well-observed He-rich SLSNe-Ib have either strongly undulating LCs or the LCs are much better fit by the CSM+Ni model. These observations imply that their progenitor stars have not had enough time to lose all of the He-envelopes before supernova explosions, and H-poor CSM are likely to present in these events. △ Less

Submitted 3 November, 2022; v1 submitted 4 February, 2022; originally announced February 2022.

Comments: 30 pages, 20 figures, Accepted by APJ

The Hydrogen-Poor Superluminous Supernovae from the Zwicky Transient Facility Phase-I Survey: I. Light Curves and Measurements

Authors: Z. H. Chen, Lin Yan, T. Kangas, R. Lunnan, S. Schulze, J. Sollerman, D. A. Perley, T. -W. Chen, K. Taggart, K. R. Hinds, A. Gal-Yam, X. F. Wang, I. Andreoni, E. Bellm, J. S. Bloom, K. Burdge, A. Burgos, D. Cook, A. Dahiwale, K. De, R. Dekany, A. Dugas, S. Frederik, C. Fremling, M. Graham , et al. (18 additional authors not shown)

Abstract: During the Zwicky Transient Facility (ZTF) Phase-I operation, 78 hydrogen-poor superluminous supernovae (SLSNe-I) were discovered in less than three years, making up the largest sample from a single survey. This paper (Paper I) presents the data, including the optical/ultraviolet light curves and classification spectra, while Paper II in this series will focus on the detailed analysis of the light… ▽ More During the Zwicky Transient Facility (ZTF) Phase-I operation, 78 hydrogen-poor superluminous supernovae (SLSNe-I) were discovered in less than three years, making up the largest sample from a single survey. This paper (Paper I) presents the data, including the optical/ultraviolet light curves and classification spectra, while Paper II in this series will focus on the detailed analysis of the light curves and modeling. Our photometry is primarily taken by the ZTF in the $g,r,i$ bands, and with additional data from other ground-based facilities and Swift. The events of our sample cover a redshift range of $z = 0.06 - 0.67$, with a median and $1σ$ error (16\% and 84\% percentiles) $z_{\rm med} = 0.265^{+0.143}_{-0.135}$. The peak luminosity covers $-22.8\,{\rm mag} \leq M_{g,\rm peak} \leq -19.8$\,mag, with a median value of $-21.48^{+1.13}_{-0.61}$\,mag. Their light curves evolve slowly with the mean rest-frame rise time of $t_{\rm rise} = 41.9\pm17.8$\,days. The luminosity and time scale distributions suggest that low luminosity SLSNe-I with peak luminosity $\sim -20$\,mag or extremely fast rising events ($<10$\,days) exist but are rare. We confirm previous findings that slowly rising SLSNe-I also tend to fade slowly. The rest-frame color and temperature evolution show large scatters, suggesting that the SLSN-I population may have diverse spectral energy distributions. The peak rest-frame color shows a moderate correlation with the peak absolute magnitude, i.e. brighter SLSNe-I tend to have bluer colors. With optical and ultraviolet photometry, we construct bolometric luminosity and derive a bolometric correction relation generally applicable for converting $g,r$-band photometry to bolometric luminosity for SLSNe-I. △ Less

Submitted 3 November, 2022; v1 submitted 4 February, 2022; originally announced February 2022.

Comments: 38 pages, 25 figures, Accepted by APJ

Cosmological Fast Optical Transients with the Zwicky Transient Facility: A Search for Dirty Fireballs

Authors: Anna Y. Q. Ho, Daniel A. Perley, Yuhan Yao, Dmitry Svinkin, A. de Ugarte Postigo, R. A. Perley, D. Alexander Kann, Eric Burns, Igor Andreoni, Eric C. Bellm, Elisabetta Bissaldi, Joshua S. Bloom, Richard Dekany, Andrew J. Drake, José Feliciano Agüí Fernández, Dmitry Frederiks, Matthew J. Graham, Boyan A. Hristov, Mansi M. Kasliwal, S. R. Kulkarni, Harsh Kumar, Russ R. Laher, Alexandra L. Lysenko, Bagrat Mailyan, Christian Malacaria , et al. (11 additional authors not shown)

Abstract: Dirty fireballs are a hypothesized class of relativistic massive-star explosions with an initial Lorentz factor $Γ_\mathrm{init}$ below the $Γ_\mathrm{init}\sim100$ required to produce a long-duration gamma-ray burst (LGRB), but which could still produce optical emission resembling LGRB afterglows. Here we present the results of a search for on-axis optical afterglows using the Zwicky Transient Fa… ▽ More Dirty fireballs are a hypothesized class of relativistic massive-star explosions with an initial Lorentz factor $Γ_\mathrm{init}$ below the $Γ_\mathrm{init}\sim100$ required to produce a long-duration gamma-ray burst (LGRB), but which could still produce optical emission resembling LGRB afterglows. Here we present the results of a search for on-axis optical afterglows using the Zwicky Transient Facility (ZTF). Our search yielded seven optical transients that resemble on-axis LGRB afterglows in terms of their red colors ($g-r>0$ mag), faint host galaxy ($r>23$ mag), and rapid fading ($dr/dt>1$ mag/day). Spectroscopy of the transient emission within a few days of discovery established cosmological distances ($z=0.876$ to $z=2.9$) for six events, tripling the number of afterglows with redshift measurements discovered by optical surveys without a $γ$-ray trigger. Upon a retrospective search, four events (ZTF20abbiixp/AT2020kym, ZTF21aagwbjr/AT2021buv, ZTF21aakruew/AT2021cwd, ZTF21abfmpwn/AT2021qbd) turned out to have a likely associated LGRB (GRB200524A, GRB210204A, GRB210212B, GRB210610B), while three did not (ZTF20aajnksq/AT2020blt, ZTF21aaeyldq/AT2021any, ZTF21aayokph/AT2021lfa). Our search revealed no definitive new class of events: the simplest explanation for the apparently "orphan" events is that they were regular LGRBs missed by high-energy satellites due to detector sensitivity and duty cycle, although it is possible that they were intrinsically faint in $γ$-rays or viewed slightly off-axis. We rule out a scenario in which dirty fireballs have a similar energy per solid angle to LGRBs and are an order of magnitude more common. In addition, we set the first direct constraint on the ratio of the opening angles of the material producing $γ$-rays and the material producing early optical afterglow emission, finding that they must be comparable. △ Less

Submitted 22 August, 2022; v1 submitted 28 January, 2022; originally announced January 2022.

Comments: 37 pages, 12 figures. Accepted for publication in ApJ