Showing 1–50 of 158 results for author: Bellm, C

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

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)

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

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

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

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

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

Supernova Siblings and their Parent Galaxies in the Zwicky Transient Facility Bright Transient Surve

Authors: M. L. Graham, C. Fremling, D. A. Perley, R. Biswas, C. A. Phillips, J. Sollerman, P. E. Nugent, S. Nance, S. Dhawan, J. Nordin, A. Goobar, A. Miller, J. D. Neill, X. J. Hall, M. J. Hankins, D. A. Duev, M. M. Kasliwal, M. Rigault, E. C. Bellm, D. Hale, P. Mróz, S. R. Kulkarni

Abstract: Supernova (SN) siblings -- two or more SNe in the same parent galaxy -- are useful tools for exploring progenitor stellar populations as well as properties of the host galaxies such as distance, star formation rate, dust extinction, and metallicity. Since the average SN rate for a Milky Way-type galaxy is just one per century, a large imaging survey is required to discover an appreciable sample of… ▽ More Supernova (SN) siblings -- two or more SNe in the same parent galaxy -- are useful tools for exploring progenitor stellar populations as well as properties of the host galaxies such as distance, star formation rate, dust extinction, and metallicity. Since the average SN rate for a Milky Way-type galaxy is just one per century, a large imaging survey is required to discover an appreciable sample of SN siblings. From the wide-field Zwicky Transient Facility (ZTF) Bright Transient Survey (BTS; which aims for spectroscopic completeness for all transients which peak brighter than $r{<}$18.5 mag) we present 10 SN siblings in 5 parent galaxies. For each of these families we analyze the SN's location within the host and its underlying stellar population, finding agreement with expectations that SNe from more massive progenitors are found nearer to their host core and in regions of more active star formation. We also present an analysis of the relative rates of core collapse and thermonuclear SN siblings, finding a significantly lower ratio than past SN sibling samples due to the unbiased nature of the ZTF. △ Less

Submitted 29 December, 2021; originally announced December 2021.

Comments: 14 pages, 4 tables, 7 figures, accepted to MNRAS

Microlensing Events in the Galactic Plane Using the Zwicky Transient Facility

Authors: Antonio C. Rodriguez, Przemek Mróz, Shrinivas R. Kulkarni, Igor Andreoni, Eric C. Bellm, Richard Dekany, Andrew J. Drake, Dmitry A. Duev, Frank J. Masci, Thomas A. Prince, Reed Riddle, David L. Shupe

Abstract: Microlensing is a powerful technique to study the Galactic population of "dark" objects such as exoplanets both bound and unbound, brown dwarfs, low-luminosity stars, old white dwarfs, neutron stars, and almost the only way to study isolated stellar-mass black holes. The majority of previous efforts to search for gravitational microlensing events have concentrated towards high-density fields such… ▽ More Microlensing is a powerful technique to study the Galactic population of "dark" objects such as exoplanets both bound and unbound, brown dwarfs, low-luminosity stars, old white dwarfs, neutron stars, and almost the only way to study isolated stellar-mass black holes. The majority of previous efforts to search for gravitational microlensing events have concentrated towards high-density fields such as the Galactic bulge. Microlensing events in the Galactic plane have the advantage of closer proximity and better constrained relative proper motions, leading to better constrained lens mass estimates at the expense of a lower optical depth compared to events towards the Galactic bulge. We use the Zwicky Transient Facility (ZTF) Data Release 5 (DR5) compiled from 2018--2021 to survey the Galactic plane in the region of $|b| < 20^\circ$. We find a total of 60 candidate microlensing events including three that show a strong microlensing parallax effect. The rate of events traces Galactic structure, decreasing exponentially as a function Galactic longitude with scale length $\ell_0 \sim 37^\circ$. On average, we find Einstein timescales of our microlensing events to be about three times as long ($\sim60$ days) compared to those towards the Galactic bulge ($\sim20$ days). This pilot project demonstrates that microlensing towards the Galactic plane shows strong promise for characterization of dark objects within the Galatic disk. △ Less

Submitted 14 December, 2021; originally announced December 2021.

Comments: submitted to ApJ

Characterizing Sparse Asteroid Light Curves with Gaussian Processes

Authors: Christina Willecke Lindberg, Daniela Huppenkothen, R. Lynne Jones, Bryce T. Bolin, Mario Juric, V. Zach Golkhou, Eric C. Bellm, Andrew J. Drake, Matthew J. Graham, Russ R. Laher, Ashish A. Mahabal, Frank J. Masci, Reed Riddle, Kyung Min Shin

Abstract: In the era of wide-field surveys like the Zwicky Transient Facility and the Rubin Observatory's Legacy Survey of Space and Time, sparse photometric measurements constitute an increasing percentage of asteroid observations, particularly for asteroids newly discovered in these large surveys. Follow-up observations to supplement these sparse data may be prohibitively expensive in many cases, so to ov… ▽ More In the era of wide-field surveys like the Zwicky Transient Facility and the Rubin Observatory's Legacy Survey of Space and Time, sparse photometric measurements constitute an increasing percentage of asteroid observations, particularly for asteroids newly discovered in these large surveys. Follow-up observations to supplement these sparse data may be prohibitively expensive in many cases, so to overcome these sampling limitations, we introduce a flexible model based on Gaussian Processes to enable Bayesian parameter inference of asteroid time series data. This model is designed to be flexible and extensible, and can model multiple asteroid properties such as the rotation period, light curve amplitude, changing pulse profile, and magnitude changes due to the phase angle evolution at the same time. Here, we focus on the inference of rotation periods. Based on both simulated light curves and real observations from the Zwicky Transient Facility, we show that the new model reliably infers rotational periods from sparsely sampled light curves, and generally provides well-constrained posterior probability densities for the model parameters. We propose this framework as an intermediate method between fast, but very limited period detection algorithms and much more comprehensive, but computationally expensive shape modeling based on ray-tracing codes. △ Less

Submitted 24 November, 2021; originally announced November 2021.

Comments: 27 pages, 18 figures, accepted for publication in AJ, associated software available at https://github.com/dirac-institute/asterogap/tree/v0.1

A WC/WO star exploding within an expanding carbon-oxygen-neon nebula

Authors: A. Gal-Yam, R. Bruch, S. Schulze, Y. Yang, D. A. Perley, I. Irani, J. Sollerman, E. C. Kool, M. T. Soumagnac, O. Yaron, N. L. Strotjohann, E. Zimmerman, C. Barbarino, S. R. Kulkarni, M. M. Kasliwal, K. De, Y. Yao, C. Fremling, L. Yan, E. O. Ofek, C. Fransson, A. V. Filippenko, W. Zheng, T. G. Brink, C. M. Copperwheat , et al. (24 additional authors not shown)

Abstract: The final explosive fate of massive stars, and the nature of the compact remnants they leave behind (black holes and neutron stars), are major open questions in astrophysics. Many massive stars are stripped of their outer hydrogen envelopes as they evolve. Such Wolf-Rayet (W-R) stars emit strong and rapidly expanding (v_wind>1000 km/s) winds indicating a high escape velocity from the stellar surfa… ▽ More The final explosive fate of massive stars, and the nature of the compact remnants they leave behind (black holes and neutron stars), are major open questions in astrophysics. Many massive stars are stripped of their outer hydrogen envelopes as they evolve. Such Wolf-Rayet (W-R) stars emit strong and rapidly expanding (v_wind>1000 km/s) winds indicating a high escape velocity from the stellar surface. A fraction of this population is also helium depleted, with spectra dominated by highly-ionized emission lines of carbon and oxygen (Types WC/WO). Evidence indicates that the most commonly-observed supernova (SN) explosions that lack hydrogen and helium (Types Ib/Ic) cannot result from massive WC/WO stars, leading some to suggest that most such stars collapse directly into black holes without a visible supernova explosions. Here, we present observations of supernova SN 2019hgp, discovered about a day after explosion. The short rise time and rapid decline place it among an emerging population of rapidly-evolving transients (RETs). Spectroscopy reveals a rich set of emission lines indicating that the explosion occurred within a nebula composed of carbon, oxygen, and neon. Narrow absorption features show that this material is expanding at relatively high velocities (>1500 km/s) requiring a compact progenitor. Our observations are consistent with an explosion of a massive WC/WO star, and suggest that massive W-R stars may be the progenitors of some rapidly evolving transients. △ Less

Submitted 24 November, 2021; originally announced November 2021.

Comments: Unedited author version, Nature in press

The Type Icn SN 2021csp: Implications for the Origins of the Fastest Supernovae and the Fates of Wolf-Rayet Stars

Authors: Daniel A. Perley, Jesper Sollerman, Steve Schulze, Yuhan Yao, Christoffer Fremling, Avishay Gal-Yam, Anna Y. Q. Ho, Yi Yang, Erik C. Kool, Ido Irani, Lin Yan, Igor Andreoni, Dietrich Baade, Eric C. Bellm, Thomas G. Brink, Ting-Wan Chen, Aleksandar Cikota, Michael W. Coughlin, Richard Dekany, Dmitry A. Duev, Alexei V. Filippenko, Peter Hoeflich, Mansi M. Kasliwal, S. R. Kulkarni, Ragnhild Lunnan , et al. (9 additional authors not shown)

Abstract: We present observations of SN 2021csp, the second example of a newly-identified type of supernova (Type Icn) hallmarked by strong, narrow, P Cygni carbon features at early times. The SN appears as a fast and luminous blue transient at early times, reaching a peak absolute magnitude of -20 within 3 days due to strong interaction between fast SN ejecta (v ~ 30000 km/s) and a massive, dense, fast-mov… ▽ More We present observations of SN 2021csp, the second example of a newly-identified type of supernova (Type Icn) hallmarked by strong, narrow, P Cygni carbon features at early times. The SN appears as a fast and luminous blue transient at early times, reaching a peak absolute magnitude of -20 within 3 days due to strong interaction between fast SN ejecta (v ~ 30000 km/s) and a massive, dense, fast-moving C/O wind shed by the WC-like progenitor months before explosion. The narrow line features disappear from the spectrum 10-20 days after explosion and are replaced by a blue continuum dominated by broad Fe features, reminiscent of Type Ibn and IIn supernovae and indicative of weaker interaction with more extended H/He-poor material. The transient then abruptly fades ~60 days post-explosion when interaction ceases. Deep limits at later phases suggest minimal heavy-element nucleosynthesis, a low ejecta mass, or both, and imply an origin distinct from that of classical Type Ic supernovae. We place SN 2021csp in context with other fast-evolving interacting transients, and discuss various progenitor scenarios: an ultrastripped progenitor star, a pulsational pair-instability eruption, or a jet-driven fallback supernova from a Wolf-Rayet star. The fallback scenario would naturally explain the similarity between these events and radio-loud fast transients, and suggests a picture in which most stars massive enough to undergo a WR phase collapse directly to black holes at the end of their lives. △ Less

Submitted 11 January, 2022; v1 submitted 23 November, 2021; originally announced November 2021.

Comments: Accepted to ApJ 29th December 2021. Originally submitted to ApJ on 6th August 2021

Blazar variability with the Vera C. Rubin Legacy Survey of Space and Time (LSST)

Authors: Claudia M. Raiteri, Maria I. Carnerero, Barbara Balmaverde, Eric C. Bellm, William Clarkson, Filippo D'Ammando, Maurizio Paolillo, Gordon T. Richards, Massimo Villata, Peter Yoachim, Ilsang Yoon

Abstract: With their emission mainly coming from a relativistic jet pointing towards us, blazars are fundamental sources to study extragalactic jets and their central engines, consisting of supermassive black holes (SMBHs) fed by accretion discs. They are also candidate sources of high-energy neutrinos and cosmic rays. Because of the jet orientation, the non-thermal blazar emission is Doppler beamed; its va… ▽ More With their emission mainly coming from a relativistic jet pointing towards us, blazars are fundamental sources to study extragalactic jets and their central engines, consisting of supermassive black holes (SMBHs) fed by accretion discs. They are also candidate sources of high-energy neutrinos and cosmic rays. Because of the jet orientation, the non-thermal blazar emission is Doppler beamed; its variability is unpredictable and occurs on time-scales from less than one hour to years. The comprehension of the diverse mechanisms producing the flux and spectral changes requires well-sampled multiband light curves on long time periods. In particular, outbursts are the best test bench to shed light on the underlying physics, especially when studied in a multiwavelength context. The Vera C. Rubin Legacy Survey of Space and Time (Rubin-LSST) will monitor the southern sky for ten years in six photometric bands, offering a formidable tool to study blazar variability features in a statistical way. The alert system will allow us to trigger follow-up observations of outstanding events, especially at high (keV-to-GeV) and very high (TeV) energies. We here examine the simulated Rubin-LSST survey strategies with the aim of understanding which cadences are more suitable for the blazar variability science. Our metrics include light curve and colour sampling. We also investigate the problem of saturation, which will affect the brightest and many flaring sources, and will have a detrimental impact on follow-up observations. △ Less

Submitted 19 November, 2021; originally announced November 2021.

Comments: 18 pages, 12 figures, accepted for publication in the ApJS Rubin Cadence Focus Issue

Establishing accretion flares from massive black holes as a source of high-energy neutrinos

Authors: Sjoert van Velzen, Robert Stein, Marat Gilfanov, Marek Kowalski, Kimitake Hayasaki, Simeon Reusch, Yuhan Yao, Simone Garrappa, Anna Franckowiak, Suvi Gezari, Jakob Nordin, Christoffer Fremling, Yashvi Sharma, Lin Yan, Erik C. Kool, Daniel Stern, Patrik M. Veres, Jesper Sollerman, Pavel Medvedev, Rashid Sunyaev, Eric C. Bellm, Richard G. Dekany, Dimitri A. Duev, Matthew J. Graham, Mansi M. Kasliwal , et al. (4 additional authors not shown)

Abstract: The origin of cosmic high-energy neutrinos remains largely unexplained. For high-energy neutrino alerts from IceCube, a coincidence with time-variable emission has been seen for three different types of accreting black holes: (1) a gamma-ray flare from a blazar (TXS 0506+056), (2) an optical transient following a stellar tidal disruption event (TDE; AT2019dsg), and (3) an optical outburst from an… ▽ More The origin of cosmic high-energy neutrinos remains largely unexplained. For high-energy neutrino alerts from IceCube, a coincidence with time-variable emission has been seen for three different types of accreting black holes: (1) a gamma-ray flare from a blazar (TXS 0506+056), (2) an optical transient following a stellar tidal disruption event (TDE; AT2019dsg), and (3) an optical outburst from an active galactic nucleus (AGN; AT2019fdr). For the latter two sources, infrared follow-up observations revealed a powerful reverberation signal due to dust heated by the flare. This discovery motivates a systematic study of neutrino emission from all supermassive black hole with similar dust echoes. Because dust reprocessing is agnostic to the origin of the outburst, our work unifies TDEs and high-amplitude flares from AGN into a population that we dub accretion flares. Besides the two known events, we uncover a third flare that is coincident with a PeV-scale neutrino (AT2019aalc). Based solely on the optical and infrared properties, we estimate a significance of 3.6$σ$ for this association of high-energy neutrinos with three accretion flares. Our results imply that at least ~10% of the IceCube high-energy neutrino alerts could be due to accretion flares. This is surprising because the sum of the fluence of these flares is at least three orders of magnitude lower compared to the total fluence of normal AGN. It thus appears that the efficiency of high-energy neutrino production in accretion flares is increased compared to non-flaring AGN. We speculate that this can be explained by the high Eddington ratio of the flares. △ Less

Submitted 3 April, 2024; v1 submitted 17 November, 2021; originally announced November 2021.

Comments: Accepted for publication in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society (2024), Volume 529, Issue 3, 2559-2576

Candidate Tidal Disruption Event AT2019fdr Coincident with a High-Energy Neutrino

Authors: Simeon Reusch, Robert Stein, Marek Kowalski, Sjoert van Velzen, Anna Franckowiak, Cecilia Lunardini, Kohta Murase, Walter Winter, James C. A. Miller-Jones, Mansi M. Kasliwal, Marat Gilfanov, Simone Garrappa, Vaidehi S. Paliya, Tomas Ahumada, Shreya Anand, Cristina Barbarino, Eric C. Bellm, Valery Brinnel, Sara Buson, S. Bradley Cenko, Michael W. Coughlin, Kishalay De, Richard Dekany, Sara Frederick, Avishay Gal-Yam , et al. (24 additional authors not shown)

Abstract: The origins of the high-energy cosmic neutrino flux remain largely unknown. Recently, one high-energy neutrino was associated with a tidal disruption event (TDE). Here we present AT2019fdr, an exceptionally luminous TDE candidate, coincident with another high-energy neutrino. Our observations, including a bright dust echo and soft late-time X-ray emission, further support a TDE origin of this flar… ▽ More The origins of the high-energy cosmic neutrino flux remain largely unknown. Recently, one high-energy neutrino was associated with a tidal disruption event (TDE). Here we present AT2019fdr, an exceptionally luminous TDE candidate, coincident with another high-energy neutrino. Our observations, including a bright dust echo and soft late-time X-ray emission, further support a TDE origin of this flare. The probability of finding two such bright events by chance is just 0.034%. We evaluate several models for neutrino production and show that AT2019fdr is capable of producing the observed high-energy neutrino, reinforcing the case for TDEs as neutrino sources. △ Less

Submitted 10 June, 2022; v1 submitted 17 November, 2021; originally announced November 2021.

Comments: 20 pages, 6 figures, 6 tables

Journal ref: Phys. Rev. Lett. 128, 221101 (2022)

Variability-selected intermediate mass black hole candidates in dwarf galaxies from ZTF and WISE

Authors: Charlotte Ward, Suvi Gezari, Peter Nugent, Eric C. Bellm, Richard Dekany, Andrew Drake, Dmitry A. Duev, Matthew J. Graham, Mansi M. Kasliwal, Erik C. Kool, Frank J. Masci, Reed L. Riddle

Abstract: While it is difficult to observe the first black hole seeds in the early Universe, we can study intermediate mass black holes (IMBHs) in local dwarf galaxies for clues about their origins. In this paper we present a sample of variability--selected AGN in dwarf galaxies using optical photometry from the Zwicky Transient Facility (ZTF) and forward--modeled mid--IR photometry of time--resolved Wide--… ▽ More While it is difficult to observe the first black hole seeds in the early Universe, we can study intermediate mass black holes (IMBHs) in local dwarf galaxies for clues about their origins. In this paper we present a sample of variability--selected AGN in dwarf galaxies using optical photometry from the Zwicky Transient Facility (ZTF) and forward--modeled mid--IR photometry of time--resolved Wide--field Infrared Survey Explorer ({\it WISE}) coadded images. We found that 44 out of 25,714 dwarf galaxies had optically variable AGN candidates, and 148 out of 79,879 dwarf galaxies had mid--IR variable AGN candidates, corresponding to active fractions of $0.17\pm0.03$\% and $0.19\pm0.02$\% respectively. We found that spectroscopic approaches to AGN identification would have missed 81\% of our ZTF IMBH candidates and 69\% of our {\it WISE} IMBH candidates. Only $9$ candidates have been detected previously in radio, X-ray, and variability searches for dwarf galaxy AGN. The ZTF and {\it WISE} dwarf galaxy AGN with broad Balmer lines have virial masses down to $10^{5.5}M_\odot$ and for the rest of the sample, BH masses predicted from host galaxy mass range between $10^{5.2}M_\odot<M_{\text{BH}}<10^{7.3}M_\odot$. We found that only 5 of 152 previously reported variability--selected AGN candidates from the Palomar Transient Factory in common with our parent sample were variable in ZTF. We also determined a nuclear supernova fraction of $0.05\pm0.01$\% year$^{-1}$ for dwarf galaxies in ZTF. Our ZTF and {\it WISE} IMBH candidates show the promise of variability searches for the discovery of otherwise hidden low mass AGN. △ Less

Submitted 15 July, 2022; v1 submitted 25 October, 2021; originally announced October 2021.

Comments: Submitted to ApJ. 25 pages, 10 figures, 4 tables

Discovery of a double detonation thermonuclear supernova progenitor

Authors: Thomas Kupfer, Evan B. Bauer, Jan van Roestel, Eric C. Bellm, Lars Bildsten, Jim Fuller, Thomas A. Prince, Ulrich Heber, Stephan Geier, Matthew J. Green, Shrinivas R. Kulkarni, Steven Bloemen, Russ R. Laher, Ben Rusholme, David Schneider

Abstract: We present the discovery of a new double detonation progenitor system consisting of a hot subdwarf B (sdB) binary with a white dwarf companion with an P=76.34179(2) min orbital period. Spectroscopic observations are consistent with an sdB star during helium core burning residing on the extreme horizontal branch. Chimera light curves are dominated by ellipsoidal deformation of the sdB star and a we… ▽ More We present the discovery of a new double detonation progenitor system consisting of a hot subdwarf B (sdB) binary with a white dwarf companion with an P=76.34179(2) min orbital period. Spectroscopic observations are consistent with an sdB star during helium core burning residing on the extreme horizontal branch. Chimera light curves are dominated by ellipsoidal deformation of the sdB star and a weak eclipse of the companion white dwarf. Combining spectroscopic and light curve fits we find a low mass sdB star, $M_{\rm sdB}=0.383\pm0.028$ M$_\odot$ with a massive white dwarf companion, $M_{\rm WD}=0.725\pm0.026$ M$_\odot$. From the eclipses we find a blackbody temperature for the white dwarf of 26,800 K resulting in a cooling age of $\approx$25 Myrs whereas our MESA model predicts an sdB age of $\approx$170 Myrs. We conclude that the sdB formed first through stable mass transfer followed by a common envelope which led to the formation of the white dwarf companion $\approx$25 Myrs ago. Using the MESA stellar evolutionary code we find that the sdB star will start mass transfer in $\approx$6 Myrs and in $\approx$60 Myrs the white dwarf will reach a total mass of $0.92$ M$_\odot$ with a thick helium layer of $0.17$ M$_\odot$. This will lead to a detonation that will likely destroy the white dwarf in a peculiar thermonuclear supernova. PTF1 2238+7430 is only the second confirmed candidate for a double detonation thermonuclear supernova. Using both systems we estimate that at least $\approx$1% of white dwarf thermonuclear supernovae originate from sdB+WD binaries with thick helium layers, consistent with the small number of observed peculiar thermonuclear explosions. △ Less

Submitted 6 January, 2022; v1 submitted 22 October, 2021; originally announced October 2021.

Comments: accepted for publication in ApJL, 14 pages, 5 figures, 1 table

Give Me a Few Hours: Exploring Short Timescales in Rubin Observatory Cadence Simulations

Authors: Eric C. Bellm, Colin J. Burke, Michael W. Coughlin, Igor Andreoni, Claudia M. Raiteri, Rosaria Bonito

Abstract: The limiting temporal resolution of a time-domain survey in detecting transient behavior is set by the time between observations of the same sky area. We analyze the distribution of visit separations for a range of Vera C. Rubin Observatory cadence simulations. Current simulations are strongly peaked at the 22 minute visit pair separation and provide effectively no constraint on temporal evolution… ▽ More The limiting temporal resolution of a time-domain survey in detecting transient behavior is set by the time between observations of the same sky area. We analyze the distribution of visit separations for a range of Vera C. Rubin Observatory cadence simulations. Current simulations are strongly peaked at the 22 minute visit pair separation and provide effectively no constraint on temporal evolution within the night. This choice will necessarily prevent Rubin from discovering a wide range of astrophysical phenomena in time to trigger rapid followup. We present a science-agnostic metric to supplement detailed simulations of fast-evolving transients and variables and suggest potential approaches for improving the range of timescales explored. △ Less

Submitted 5 October, 2021; originally announced October 2021.

Comments: 7 pages, 2 figures. Submitted to the ApJS Rubin Cadence Focus Issue