Showing 1–50 of 137 results for author: Laher, R R

ZTF SN Ia DR2: The secondary maximum in Type Ia supernovae

Authors: M. Deckers, K. Maguire, L. Shingles, G. Dimitriadis, M. Rigault, M. Smith, A. Goobar, J. Nordin, J. Johansson, M. Amenouche, U. Burgaz, S. Dhawan, M. Ginolin, L. Harvey, W. D. Kenworthy, Y. -L. Kim, R. R. Laher, N. Luo, S. R. Kulkarni, F. J. Masci, T. E. Müller-Bravo, P. E. Nugent, N. Pletskova, J. Purdum, B. Racine , et al. (2 additional authors not shown)

Abstract: Type Ia supernova (SN Ia) light curves have a secondary maximum that exists in the $r$, $i$, and near-infrared filters. The secondary maximum is relatively weak in the $r$ band, but holds the advantage that it is accessible, even at high redshift. We used Gaussian Process fitting to parameterise the light curves of 893 SNe Ia from the Zwicky Transient Facility's (ZTF) second data release (DR2), an… ▽ More Type Ia supernova (SN Ia) light curves have a secondary maximum that exists in the $r$, $i$, and near-infrared filters. The secondary maximum is relatively weak in the $r$ band, but holds the advantage that it is accessible, even at high redshift. We used Gaussian Process fitting to parameterise the light curves of 893 SNe Ia from the Zwicky Transient Facility's (ZTF) second data release (DR2), and we were able to extract information about the timing and strength of the secondary maximum. We found $>5σ$ correlations between the light curve decline rate ($Δm_{15}(g)$) and the timing and strength of the secondary maximum in the $r$ band. Whilst the timing of the secondary maximum in the $i$ band also correlates with $Δm_{15}(g)$, the strength of the secondary maximum in the $i$ band shows significant scatter as a function of $Δm_{15}(g)$. We found that the transparency timescales of 97 per cent of our sample are consistent with double detonation models, and that SNe Ia with small transparency timescales ($<$ 32 d) reside predominantly in locally red environments. We measured the total ejected mass for the normal SNe Ia in our sample using two methods, and both were consistent with medians of $1.3\ \pm \ 0.3$ and $1.2\ \pm\ 0.2$ solar masses. We find that the strength of the secondary maximum is a better standardisation parameter than the SALT light curve stretch ($x_1$). Finally, we identified a spectral feature in the $r$ band as Fe II, which strengthens during the onset of the secondary maximum. The same feature begins to strengthen at $<$ 3 d post maximum light in 91bg-like SNe. Finally, the correlation between $x_1$ and the strength of the secondary maximum was best fit with a broken line, with a split at $x_1^0\ =\ -0.5\ \pm\ 0.2$, suggestive of the existence of two populations of SNe Ia. △ Less

Submitted 27 June, 2024; originally announced June 2024.

Comments: 18 pages, 14 figures

ZTF SN Ia DR2: Impact of the galaxy cluster environment on the stretch distribution of Type Ia supernovae

Authors: F. Ruppin, M. Rigault, M. Ginolin, G. Dimitriadis, A. Goobar, J. Johansson, K. Maguire, J. Nordin, M. Smith, M. Aubert, J. Biedermann, Y. Copin, U. Burgaz, B. Carreres, F. Feinstein, D. Fouchez, T. E. Muller-Bravo, L. Galbany, S. L. Groom, W. D. Kenworthy, Y. -L. Kim, R. R. Laher, P. Nugent, B. Popovic, J. Purdum , et al. (5 additional authors not shown)

Abstract: Understanding the impact of the astrophysical environment on Type Ia supernova (SN Ia) properties is crucial to minimize systematic uncertainties in cosmological analyses based on this probe. We investigate the dependence of the SN Ia SALT2.4 light-curve stretch on the distance from their nearest galaxy cluster to study a potential effect of the intracluster medium (ICM) environment on SN Ia intri… ▽ More Understanding the impact of the astrophysical environment on Type Ia supernova (SN Ia) properties is crucial to minimize systematic uncertainties in cosmological analyses based on this probe. We investigate the dependence of the SN Ia SALT2.4 light-curve stretch on the distance from their nearest galaxy cluster to study a potential effect of the intracluster medium (ICM) environment on SN Ia intrinsic properties. We use the largest SN Ia sample to date and cross-match it with existing X-ray, Sunyaev-Zel'dovich, and optical cluster catalogs in order to study the dependence between stretch and distance to the nearest detected cluster from each SN Ia. We model the underlying stretch distribution with a Gaussian mixture with relative amplitudes that depend on redshift and cluster-centric distance. We find a significant improvement of the fit quality of the stretch distribution if we include the distance-dependant term in the model with a variation of the Akaike information criterion $\rm{ΔAIC} = -10.2$. Because of the known correlation between galaxy age and distance from cluster center, this supports previous evidence that the age of the stellar population is the underlying driver of the bimodial shape of the SN Ia stretch distribution. We further compute the evolution of the fraction of quenched galaxies as a function of distance with respect to cluster center from our best-fit model of the SNe Ia stretch distribution and compare it to previous results obtained from $Hα$ line measurements, optical broadband photometry, and simulations. We find our estimate to be compatible with these results. The results of this work indicate that SNe Ia searches at high redshift targeted towards clusters to maximize detection probability should be considered with caution as the stretch distribution of the detected sample would be strongly biased towards the old sub-population of SNe Ia. △ Less

Submitted 3 June, 2024; originally announced June 2024.

Comments: Submitted to Astronomy & 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

Sub-relativistic Outflow and Hours-Timescale Large-amplitude X-ray Dips during Super-Eddington Accretion onto a Low-mass Massive Black Hole in the Tidal Disruption Event AT2022lri

Authors: Yuhan Yao, Muryel Guolo, Francesco Tombesi, Ruancun Li, Suvi Gezari, Javier A. García, Lixin Dai, Ryan Chornock, Wenbin Lu, S. R. Kulkarni, Keith C. Gendreau, Dheeraj R. Pasham, S. Bradley Cenko, Erin Kara, Raffaella Margutti, Yukta Ajay, Thomas Wevers, Tom M. Kwan, Igor Andreoni, Joshua S. Bloom, Andrew J. Drake, Matthew J. Graham, Erica Hammerstein, Russ R. Laher, Natalie LeBaron , et al. (10 additional authors not shown)

Abstract: We present the tidal disruption event (TDE) AT2022lri, hosted in a nearby ($\approx\!144$ Mpc) quiescent galaxy with a low-mass massive black hole ($10^4\,M_\odot < M_{\rm BH} < 10^6\,M_\odot$). AT2022lri belongs to the TDE-H+He subtype. More than 1 Ms of X-ray data were collected with NICER, Swift, and XMM-Newton from 187 d to 672 d after peak. The X-ray luminosity gradually declined from… ▽ More We present the tidal disruption event (TDE) AT2022lri, hosted in a nearby ($\approx\!144$ Mpc) quiescent galaxy with a low-mass massive black hole ($10^4\,M_\odot < M_{\rm BH} < 10^6\,M_\odot$). AT2022lri belongs to the TDE-H+He subtype. More than 1 Ms of X-ray data were collected with NICER, Swift, and XMM-Newton from 187 d to 672 d after peak. The X-ray luminosity gradually declined from $1.5\times 10^{44}\,{\rm erg\,s^{-1}}$ to $1.5\times 10^{43}\,{\rm erg\,s^{-1}}$ and remains much above the UV and optical luminosity, consistent with a super-Eddington accretion flow viewed face-on. Sporadic strong X-ray dips atop a long-term decline are observed, with variability timescale of $\approx\!0.5$ hr--1 d and amplitude of $\approx\!2$--8. When fitted with simple continuum models, the X-ray spectrum is dominated by a thermal disk component with inner temperature going from $\sim\! 146$ eV to $\sim\! 86$ eV. However, there are residual features that peak around 1 keV, which, in some cases, cannot be reproduced by a single broad emission line. We analyzed a subset of time-resolved spectra with two physically motivated models describing either a scenario where ionized absorbers contribute extra absorption and emission lines or where disk reflection plays an important role. Both models provide good and statistically comparable fits, show that the X-ray dips are correlated with drops in the inner disk temperature, and require the existence of sub-relativistic (0.1--0.3$c$) ionized outflows. We propose that the disk temperature fluctuation stems from episodic drops of the mass accretion rate triggered by magnetic instabilities or/and wobbling of the inner accretion disk along the black hole's spin axis. △ Less

Submitted 18 May, 2024; originally announced May 2024.

Comments: 35 pages, 20 figures, submitted

The Zwicky Transient Facility Bright Transient Survey. III. $\texttt{BTSbot}$: Automated Identification and Follow-up of Bright Transients with Deep Learning

Authors: Nabeel Rehemtulla, Adam A. Miller, Theophile Jegou Du Laz, Michael W. Coughlin, Christoffer Fremling, Daniel A. Perley, Yu-Jing Qin, Jesper Sollerman, Ashish A. Mahabal, Russ R. Laher, Reed Riddle, Ben Rusholme, Shrinivas R. Kulkarni

Abstract: The Bright Transient Survey (BTS) aims to obtain a classification spectrum for all bright ($m_\mathrm{peak}\,\leq\,18.5\,$mag) extragalactic transients found in the Zwicky Transient Facility (ZTF) public survey. BTS critically relies on visual inspection ("scanning") to select targets for spectroscopic follow-up, which, while effective, has required a significant time investment over the past… ▽ More The Bright Transient Survey (BTS) aims to obtain a classification spectrum for all bright ($m_\mathrm{peak}\,\leq\,18.5\,$mag) extragalactic transients found in the Zwicky Transient Facility (ZTF) public survey. BTS critically relies on visual inspection ("scanning") to select targets for spectroscopic follow-up, which, while effective, has required a significant time investment over the past $\sim5$ yr of ZTF operations. We present $\texttt{BTSbot}$, a multi-modal convolutional neural network, which provides a bright transient score to individual ZTF detections using their image data and 25 extracted features. $\texttt{BTSbot}$ is able to eliminate the need for daily human scanning by automatically identifying and requesting spectroscopic follow-up observations of new bright transient candidates. $\texttt{BTSbot}$ recovers all bright transients in our test split and performs on par with scanners in terms of identification speed (on average, $\sim$1 hour quicker than scanners). We also find that $\texttt{BTSbot}$ is not significantly impacted by any data shift by comparing performance across a concealed test split and a sample of very recent BTS candidates. $\texttt{BTSbot}$ has been integrated into Fritz and $\texttt{Kowalski}$, ZTF's first-party marshal and alert broker, and now sends automatic spectroscopic follow-up requests for the new transients it identifies. During the month of October 2023, $\texttt{BTSbot}$ selected 296 sources in real-time, 93% of which were real extragalactic transients. With $\texttt{BTSbot}$ and other automation tools, the BTS workflow has produced the first fully automatic end-to-end discovery and classification of a transient, representing a significant reduction in the human-time needed to scan. Future development has tremendous potential for creating similar models to identify and request follow-up observations for specific types of transients. △ Less

Submitted 26 January, 2024; originally announced January 2024.

Comments: 26 pages, 12 figures; to be submitted to ApJ; comments welcome

A Joint SRG/eROSITA + ZTF Search: Discovery of a 97-min Period Eclipsing Cataclysmic Variable with Evidence of a Brown Dwarf Secondary

Authors: Ilkham Galiullin, Antonio C. Rodriguez, Shrinivas R. Kulkarni, Rashid Sunyaev, Marat Gilfanov, Ilfan Bikmaev, Lev Yungelson, Jan van Roestel, Boris T. Gänsicke, Irek Khamitov, Paula Szkody, Kareem El-Badry, Mikhail Suslikov, Thomas A. Prince, Mikhail Buntov, Ilaria Caiazzo, Mark Gorbachev, Matthew J. Graham, Rustam Gumerov, Eldar Irtuganov, Russ R. Laher, Pavel Medvedev, Reed Riddle, Ben Rusholme, Nail Sakhibullin , et al. (2 additional authors not shown)

Abstract: Cataclysmic variables (CVs) that have evolved past the period minimum during their lifetimes are predicted to be systems with a brown dwarf donor. While population synthesis models predict that around $\approx 40-70\%$ of the Galactic CVs are post-period minimum systems referred to as "period bouncers", only a few dozen confirmed systems are known. We report the study and characterisation of a new… ▽ More Cataclysmic variables (CVs) that have evolved past the period minimum during their lifetimes are predicted to be systems with a brown dwarf donor. While population synthesis models predict that around $\approx 40-70\%$ of the Galactic CVs are post-period minimum systems referred to as "period bouncers", only a few dozen confirmed systems are known. We report the study and characterisation of a new eclipsing CV, SRGeJ041130.3+685350 (SRGeJ0411), discovered from a joint SRG/eROSITA and ZTF program. The optical spectrum of SRGeJ0411 shows prominent hydrogen and helium emission lines, typical for CVs. We obtained optical high-speed photometry to confirm the eclipse of SRGeJ0411 and determine the orbital period to be $P_\textrm{orb} \approx 97.530$ minutes. The spectral energy distribution suggests that the donor has an effective temperature of $\lesssim 1,800$ K. We constrain the donor mass with the period--density relationship for Roche-lobe-filling stars and find that $M_\textrm{donor} \lesssim 0.04\ M_\odot$. The binary parameters are consistent with evolutionary models for post-period minimum CVs, suggesting that SRGeJ0411 is a new period bouncer. The optical emission lines of SRGeJ0411 are single-peaked despite the system being eclipsing, which is typically only seen due to stream-fed accretion in polars. X-ray spectroscopy hints that the white dwarf in SRGeJ0411 could be magnetic, but verifying the magnetic nature of SRGeJ0411 requires further investigation. The lack of optical outbursts has made SRGeJ0411 elusive in previous surveys, and joint X-ray and optical surveys highlight the potential for discovering similar systems in the near future. △ Less

Submitted 8 January, 2024; originally announced January 2024.

Comments: 17 pages, 16 figures, accepted to MNRAS

SN 2021adxl: A luminous nearby interacting supernova in an extremely low metallicity environment

Authors: S. J. Brennan, S. Schulze, R. Lunnan, J. Sollerman, L. Yan, C. Fransson, I. Irani, J. Melinder, T. -W. Chen, K. De, C. Fremling, Y. -L. Kim, D. Perley, P. J. Pessi, A. J. Drake, M. J. Graham, R. R. Laher, F. J. Masci, J. Purdum, H. Rodriguez

Abstract: SN 2021adxl is a slowly evolving, luminous, Type IIn supernova with asymmetric emission line profiles, similar to the well-studied SN 2010jl. We present extensive optical, near-ultraviolet, and near-infrared photometry and spectroscopy covering ~1.5 years post discovery. SN 2021adxl occurred in an unusual environment, atop a vigorously star-forming region that is offset from its host galaxy core.… ▽ More SN 2021adxl is a slowly evolving, luminous, Type IIn supernova with asymmetric emission line profiles, similar to the well-studied SN 2010jl. We present extensive optical, near-ultraviolet, and near-infrared photometry and spectroscopy covering ~1.5 years post discovery. SN 2021adxl occurred in an unusual environment, atop a vigorously star-forming region that is offset from its host galaxy core. The appearance of Lyα, O II, as well as the compact core, would classify the host of SN 2021adxl as a Blueberry galaxy, analogous to the higher redshift Green Pea galaxies. Using several abundance indicators, we find a metallicity of the explosion environment of only 10% solar, the lowest reported metallicity for a Type IIn SN environment. SN 2021adxl reaches a peak magnitude of r ~ -20.2 mag and since discovery, SN 2021adxl has faded by only ~4 magnitudes in the r band with a cumulative radiated energy of ~1.5e50 erg over 18 months. SN 2021adxl shows strong signs of interaction with a complex circumstellar medium, seen by the detection of X-rays, revealed by the detection of coronal emission lines, and through multi-component hydrogen and helium profiles. In order to further understand this interaction, we model the Hα profile using a Monte-Carlo electron scattering code. The blueshifted high-velocity component is consistent with emission from a radially thin, spherical shell resulting in the broad emission components due to electron scattering. Using the velocity evolution of this emitting shell, we find that the SN ejecta collide with circumstellar material of at least 5 Msun, assuming a steady-state mass-loss rate of 4-6e-3 Msun per year for the first ~200 days of evolution. Continuing the observations of SN 2021adxl may reveal signatures of dust formation or an infrared excess, similar to that seen for SN 2010jl. △ Less

Submitted 20 December, 2023; originally announced December 2023.

Comments: 21 pages with 27 Figures. Submitted to A&A. MC line scattering code can be found at https://github.com/Astro-Sean/escatter

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)

$\texttt{tdescore}$: An Accurate Photometric Classifier for Tidal Disruption Events

Authors: Robert Stein, Ashish Mahabal, Simeon Reusch, Matthew Graham, Mansi M. Kasliwal, Marek Kowalski, Suvi Gezari, Erica Hammerstein, Szymon J. Nakoneczny, Matt Nicholl, Jesper Sollerman, Sjoert van Velzen, Yuhan Yao, Russ R. Laher, Ben Rusholme

Abstract: Optical surveys have become increasingly adept at identifying candidate Tidal Disruption Events (TDEs) in large numbers, but classifying these generally requires extensive spectroscopic resources. Here we present $\texttt{tdescore}$, a simple binary photometric classifier that is trained using a systematic census of $\sim$3000 nuclear transients from the Zwicky Transient Facility (ZTF). The sample… ▽ More Optical surveys have become increasingly adept at identifying candidate Tidal Disruption Events (TDEs) in large numbers, but classifying these generally requires extensive spectroscopic resources. Here we present $\texttt{tdescore}$, a simple binary photometric classifier that is trained using a systematic census of $\sim$3000 nuclear transients from the Zwicky Transient Facility (ZTF). The sample is highly imbalanced, with TDEs representing $\sim$2% of the total. $\texttt{tdescore}$ is nonetheless able to reject non-TDEs with 99.6% accuracy, yielding a sample of probable TDEs with recall of 77.5% for a precision of 80.2%. $\texttt{tdescore}$ is thus substantially better than any available TDE photometric classifier scheme in the literature, with performance not far from spectroscopy as a method for classifying ZTF nuclear transients, despite relying solely on ZTF data and multi-wavelength catalogue cross-matching. In a novel extension, we use `SHapley Additive exPlanations' ($\texttt{SHAP}$) to provide a human-readable justification for each individual $\texttt{tdescore}$ classification, enabling users to understand and form opinions about the underlying classifier reasoning. $\texttt{tdescore}$ can serve as a model for photometric identification of TDEs with time-domain surveys, such as the upcoming Rubin observatory. △ Less

Submitted 4 April, 2024; v1 submitted 30 November, 2023; originally announced December 2023.

Comments: Accepted for publication, 16 pages, 7 figures, 2 tables

Journal ref: ApJL 965 L14 2024

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

Panic at the ISCO: time-varying double-peaked broad lines from evolving accretion disks are common amongst optically variable AGN

Authors: Charlotte Ward, Suvi Gezari, Peter Nugent, Matthew Kerr, Michael Eracleous, Sara Frederick, Erica Hammerstein, Matthew J. Graham, Sjoert van Velzen, Mansi M. Kasliwal, Russ R. Laher, Frank J. Masci, Josiah Purdum, Benjamin Racine, Roger Smith

Abstract: About 3-10\% of Type I active galactic nuclei (AGN) have double-peaked broad Balmer lines in their optical spectra originating from the motion of gas in their accretion disk. Double-peaked profiles arise not only in AGN, but occasionally appear during optical flares from tidal disruption events and changing-state AGN. In this paper we identify 250 double-peaked emitters (DPEs) amongst a parent sam… ▽ More About 3-10\% of Type I active galactic nuclei (AGN) have double-peaked broad Balmer lines in their optical spectra originating from the motion of gas in their accretion disk. Double-peaked profiles arise not only in AGN, but occasionally appear during optical flares from tidal disruption events and changing-state AGN. In this paper we identify 250 double-peaked emitters (DPEs) amongst a parent sample of optically variable broad-line AGN in the Zwicky Transient Facility (ZTF) survey, corresponding to a DPE fraction of 19\%. We model spectra of the broad H$α$ emission line regions and provide a catalog of the fitted accretion disk properties for the 250 DPEs. Analysis of power spectra derived from the 5 year ZTF light curves finds that DPE light curves have similar amplitudes and power law indices to other broad-line AGN. Follow-up spectroscopy of 12 DPEs reveals that $\sim$50\% display significant changes in the relative strengths of their red and blue peaks over long $10-20$ year timescales, indicating that broad-line profile changes arising from spiral arm or hotspot rotation are common amongst optically variable DPEs. Analysis of the accretion disk parameters derived from spectroscopic modeling provides evidence that DPEs are not in a special accretion state, but are simply normal broad-line AGN viewed under the right conditions for the accretion disk to be easily visible. We include inspiraling SMBH binary candidate SDSSJ1430+2303 in our analysis, and discuss how its photometric and spectroscopic variability is consistent with the disk-emitting AGN population in ZTF. △ Less

Submitted 9 January, 2024; v1 submitted 5 September, 2023; originally announced September 2023.

Comments: Accepted to ApJ. 34 pages, 10 figures, 4 tables

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)

SRGeJ045359.9+622444: A 55-min Period Eclipsing AM CVn Discovered from a Joint SRG/eROSITA + ZTF Search

Authors: Antonio C. Rodriguez, Ilkham Galiullin, Marat Gilfanov, Shrinivas R. Kulkarni, Irek Khamitov, Ilfan Bikmaev, Jan van Roestel, Lev Yungelson, Kareem El-Badry, Rashid Sunayev, Thomas A. Prince, Mikhail Buntov, Ilaria Caiazzo, Andrew Drake, Mark Gorbachev, Matthew J. Graham, Rustam Gumerov, Eldar Irtuganov, Russ R. Laher, Frank J. Masci, Pavel Medvedev, Josiah Purdum, Nail Sakhibullin, Alexander Sklyanov, Roger Smith , et al. (2 additional authors not shown)

Abstract: AM CVn systems are ultra-compact binaries where a white dwarf accretes from a helium-rich degenerate or semi-degenerate donor. Some AM CVn systems will be among the loudest sources of gravitational waves for the upcoming Laser Interferometer Space Antenna (LISA), yet the formation channel of AM CVns remains uncertain. We report the study and characterisation of a new eclipsing AM CVn, SRGeJ045359.… ▽ More AM CVn systems are ultra-compact binaries where a white dwarf accretes from a helium-rich degenerate or semi-degenerate donor. Some AM CVn systems will be among the loudest sources of gravitational waves for the upcoming Laser Interferometer Space Antenna (LISA), yet the formation channel of AM CVns remains uncertain. We report the study and characterisation of a new eclipsing AM CVn, SRGeJ045359.9+622444 (hereafter SRGeJ0453), discovered from a joint SRG/eROSITA and ZTF program to identify cataclysmic variables (CVs). We obtained optical photometry to confirm the eclipse of SRGeJ0453 and determine the orbital period to be $P_\textrm{orb} = 55.0802 \pm 0.0003$ min. We constrain the binary parameters by modeling the high-speed photometry and radial velocity curves and find $M_\textrm{donor} = 0.044 \pm0.024 M_{\odot}$ and $R_\textrm{donor}=0.078 \pm 0.012 R_{\odot}$. The X-ray spectrum is approximated by a power-law model with an unusually flat photon index of $Γ\sim 1$ previously seen in magnetic CVs with SRG/eROSITA, but verifying the magnetic nature of SRGeJ0453 requires further investigation. Optical spectroscopy suggests that the donor star of SRGeJ0453 could have initially been a He star or a He white dwarf. SRGeJ0453 is the ninth eclipsing AM CVn system published to date, and its lack of optical outbursts have made it elusive in previous surveys. The discovery of SRGeJ0453 using joint X-ray and optical surveys highlights the potential for discovering similar systems in the near future. △ Less

Submitted 22 June, 2023; originally announced June 2023.

Comments: Submitted to ApJ

A New Forced Photometry Service for the Zwicky Transient Facility

Authors: Frank J. Masci, Russ R. Laher, Benjamin Rusholme, David Shupe, Roberta Paladini, Steve Groom, Avery Wold, Adam A. Miller, Andrew Drake

Abstract: We describe the Zwicky Transient Facility (ZTF) Forced Photometry Service (ZFPS) as developed and maintained by the ZTF Science Data System Team at IPAC/Caltech. The service is open for public use following a subscription. The ZFPS has been operational since early 2020 and has been used to generate publication quality lightcurves for a myriad of science programs. The ZFPS has been recently upgrade… ▽ More We describe the Zwicky Transient Facility (ZTF) Forced Photometry Service (ZFPS) as developed and maintained by the ZTF Science Data System Team at IPAC/Caltech. The service is open for public use following a subscription. The ZFPS has been operational since early 2020 and has been used to generate publication quality lightcurves for a myriad of science programs. The ZFPS has been recently upgraded to allow users to request forced-photometry lightcurves for up to 1500 sky positions per request in a single web-application submission. The underlying software has been recoded to take advantage of a parallel processing architecture with the most compute-intensive component rewritten in C and optimized for the available hardware. The ZTF processing cluster consists of 66 compute nodes, each hosting at least 16 physical cores. The compute nodes are generally idle following nightly real-time processing of the ZTF survey data and when other ad hoc processing tasks have been completed. The ZFPS and associated infrastructure at IPAC/Caltech therefore enable thousands of forced-photometry lightcurves to be generated along with a wealth of quality metrics to facilitate analyses and filtering of bad quality data prior to scientific use. △ Less

Submitted 11 July, 2023; v1 submitted 25 May, 2023; originally announced May 2023.

Comments: 21 pages, 11 figures. Living document. Submitted nowhere but here

SN 2019odp: A Massive Oxygen-Rich Type Ib Supernova

Authors: T. Schweyer, J. Sollerman, A. Jerkstrand, M. Ergon, T. -W. Chen, C. M. B. Omand, S. Schulze, M. W. Coughlin, I. Andreoni, C. Fremling, A. Rau, Y. Sharma, N. L. Strotjohann, L. Yan, M. J. Graham, M. M. Kasliwal, R. R. Laher, J. Purdum, P. Rosnet, B. Rusholme, R. Smith

Abstract: We present and analyze observations of the Type Ib supernova (SN) 2019odp (a.k.a ZTF19abqwtfu) covering epochs within days of the explosion to the nebular phase at 360 d post-explosion. We discuss them in the context of recombination cooling emission for the early excess emission and consider progenitor models based on the nebular phase spectra. Our observations include photometric observations ma… ▽ More We present and analyze observations of the Type Ib supernova (SN) 2019odp (a.k.a ZTF19abqwtfu) covering epochs within days of the explosion to the nebular phase at 360 d post-explosion. We discuss them in the context of recombination cooling emission for the early excess emission and consider progenitor models based on the nebular phase spectra. Our observations include photometric observations mainly in the optical and low to medium-resolution spectroscopic observations covering the complete observable time-range. We expand on existing methods to derive oxygen mass estimates from nebular phase spectroscopy. Our spectroscopic observations confirm the presence of He in the SN ejecta and we thus (re)classify it as a Type Ib supernova. From the pseudo-bolometric lightcurve we estimate a high ejecta mass $M_\text{ej} \sim 4 - 7~M_\odot$. The high ejecta mass, large nebular [O I]/[Ca II] line flux ratio ($1.2 - 1.9$) and an oxygen mass above $\gtrapprox 0.5\, M_\odot$ point towards a progenitor with pre-explosion mass higher than $18\,M_\odot$. The compact nature of the progenitor ($\lesssim 10\,R_\odot$) suggests a Wolf-Rayet (WR) star as progenitor. △ Less

Submitted 24 March, 2023; originally announced March 2023.

Comments: Submitted to A&A

Zwicky Transient Facility and Globular Clusters: The Period-Luminosity and Period-Wesenheit Relations for SX Phoenicis Variables in the gri-Band

Authors: Chow-Choong Ngeow, Anupam Bhardwaj, Matthew J. Graham, Brian F. Healy, Russ R. Laher, Reed Riddle, Avery Wold

Abstract: SX Phoenicis (SXP) variables are short period pulsating stars that exhibit a period-luminosity (PL) relation. We derived the gri-band PL and extinction-free period-Wesenheit (PW) relations, as well as the period-color (PC) and reddening-free period-Q-index (PQ) relations for 47 SXP stars in located in 21 globular clusters using the optical light curves taken from Zwicky Transient Facility (ZTF). T… ▽ More SX Phoenicis (SXP) variables are short period pulsating stars that exhibit a period-luminosity (PL) relation. We derived the gri-band PL and extinction-free period-Wesenheit (PW) relations, as well as the period-color (PC) and reddening-free period-Q-index (PQ) relations for 47 SXP stars in located in 21 globular clusters using the optical light curves taken from Zwicky Transient Facility (ZTF). These empirically relations were derived for the first time in the gri filters except for the g-band PL relation. We used our gi band PL and PW relations to derive a distance modulus to Crater II dwarf spheroidal which hosts one SXP variable. Assuming that the fundamental and first-overtone pulsation mode for the SXP variable in Crater II, we found distance moduli of $20.03 \pm 0.23$ mag and $20.37 \pm 0.24$ mag, respectively, using the PW relation, where the latter is in excellent agreement with independent RR Lyrae based distance to Crater II dwarf galaxy. △ Less

Submitted 13 March, 2023; originally announced March 2023.

Comments: 12 pages, 1 Table and 10 Figures; AJ accepted

A Distance Measurement to M33 Using Optical Photometry of Mira Variables

Authors: Jia-Yu Ou, Chow-Choong Ngeow, Anupam Bhardwaj, Matthew J. Graham, Russ R. Laher, Frank J. Masci, Reed Riddle

Abstract: We present a systematic analysis to determine and improve the pulsation periods of 1637 known long-period Mira variables in M33 using $gri$-band light curves spanning $\sim18$~years from several surveys, including M33 variability survey, Panoramic Survey Telescope and Rapid Response System, Palomar Transient Factory (PTF), intermediate PTF, and Zwicky Transient Facility. Based on these collections… ▽ More We present a systematic analysis to determine and improve the pulsation periods of 1637 known long-period Mira variables in M33 using $gri$-band light curves spanning $\sim18$~years from several surveys, including M33 variability survey, Panoramic Survey Telescope and Rapid Response System, Palomar Transient Factory (PTF), intermediate PTF, and Zwicky Transient Facility. Based on these collections of light curves, we found that optical band light curves that are as complete as possible are crucial to determine the periods of distant Miras. We demonstrated that the machine learning techniques can be used to classify Miras into O-rich and C-rich based on the $(J-K_s)$ period--color plane. Finally, We derived the distance modulus to M33 using O-rich Miras at maximum light together with our improved periods as $24.67 \pm 0.06$~mag, which is in good agreement with the recommended value given in the literature. △ Less

Submitted 6 February, 2023; originally announced February 2023.

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

Volumetric rates of Luminous Red Novae and Intermediate Luminosity Red Transients with the Zwicky Transient Facility

Authors: Viraj R. Karambelkar, Mansi M. Kasliwal, Nadejda Blagorodnova, Jesper Sollerman, Robert Aloisi, Shreya G. Anand, Igor Andreoni, Thomas G. Brink, Rachel Bruch, David Cook, Kaustav Kashyap Das, Kishalay De, Andrew Drake, Alexei V. Filippenko, Christoffer Fremling, George Helou, Anna Ho, Jacob Jencson, David Jones, Russ R. Laher, Frank J. Masci, Kishore C. Patra, Josiah Purdum, Alexander Reedy, Tawny Sit , et al. (5 additional authors not shown)

Abstract: Luminous red novae (LRNe) are transients characterized by low luminosities and expansion velocities, and are associated with mergers or common envelope ejections in stellar binaries. Intermediate-luminosity red transients (ILRTs) are an observationally similar class with unknown origins, but generally believed to either be electron capture supernovae (ECSN) in super-AGB stars, or outbursts in dust… ▽ More Luminous red novae (LRNe) are transients characterized by low luminosities and expansion velocities, and are associated with mergers or common envelope ejections in stellar binaries. Intermediate-luminosity red transients (ILRTs) are an observationally similar class with unknown origins, but generally believed to either be electron capture supernovae (ECSN) in super-AGB stars, or outbursts in dusty luminous blue variables (LBVs). In this paper, we present a systematic sample of 8 LRNe and 8 ILRTs detected as part of the Census of the Local Universe (CLU) experiment on the Zwicky Transient Facility (ZTF). The CLU experiment spectroscopically classifies ZTF transients associated with nearby ($<150$ Mpc) galaxies, achieving 80% completeness for m$_{r}<20$\,mag. Using the ZTF-CLU sample, we derive the first systematic LRNe volumetric-rate of 7.8$^{+6.5}_{-3.7}\times10^{-5}$ Mpc$^{-3}$ yr$^{-1}$ in the luminosity range $-16\leq$M$_{\rm{r}}$$\leq -11$ mag. We find that in this luminosity range, the LRN rate scales as dN/dL $\propto L^{-2.5\pm0.3}$ - significantly steeper than the previously derived scaling of $L^{-1.4\pm0.3}$ for lower luminosity LRNe (M$_{V}\geq-10$). The steeper power law for LRNe at high luminosities is consistent with the massive merger rates predicted by binary population synthesis models. We find that the rates of the brightest LRNe (M$_{r}\leq-13$ mag) are consistent with a significant fraction of them being progenitors of double compact objects (DCOs) that merge within a Hubble time. For ILRTs, we derive a volumetric rate of $2.6^{+1.8}_{-1.4}\times10^{-6}$ Mpc$^{-3}$yr$^{-1}$ for M$_{\rm{r}}\leq-13.5$, that scales as dN/dL $\propto L^{-2.5\pm0.5}$. This rate is $\approx1-5\%$ of the local core-collapse supernova rate, and is consistent with theoretical ECSN rate estimates. △ Less

Submitted 9 November, 2022; originally announced November 2022.

Comments: 32 pages, 16 figures. Submitted to ApJ, comments welcome

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

Probing the low-mass end of core-collapse supernovae using a sample of strongly-stripped Calcium-rich Type IIb Supernovae from the Zwicky Transient Facility

Authors: Kaustav K. Das, Mansi M. Kasliwal, Christoffer Fremling, Sheng Yang, Steve Schulze, Jesper Sollerman, Tawny Sit, Kishalay De, Anastasios Tzanidakis, Daniel A. Perley, Shreya Anand, Igor Andreoni, C. Barbarino, K. Brudge, Avishay Gal-Yam, Viraj Karambelkar, S. R. Kulkarni, Yashvi Sharma, Yi Yang, Yuhan Yao, Andrew Drake, Russ R. Laher, Frank J. Masci, Michael S. Medford, Harrison Reedy , et al. (1 additional authors not shown)

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

Submitted 2 October, 2023; v1 submitted 11 October, 2022; originally announced October 2022.

Comments: Accepted for publication in ApJ

SN 2022oqm -- a Ca-rich explosion of a compact progenitor embedded in C/O circumstellar material

Authors: I. Irani, Ping Chen, Jonathan Morag, S. Schulze, A. Gal-Yam, Nora L. Strotjohann, Ofer Yaron, E. A. Zimmerman, Amir Sharon, Daniel A. Perley, J. Sollerman, Aaron Tohuvavohu, Kaustav K. Das, Mansi M. Kasliwal, Rachel Bruch, Thomas G. Brink, WeiKang Zheng, Kishore C. Patra, Sergiy S. Vasylyev, Alexei V. Filippenko, Yi Yang, Matthew J. Graham, Joshua S. Bloom, Paolo Mazzali, Josiah Purdum , et al. (5 additional authors not shown)

Abstract: We present the discovery and analysis of SN\,2022oqm, a Type Ic supernova (SN) detected $<1$\,day after explosion. The SN rises to a blue and short-lived (2\,days) initial peak. Early-time spectral observations of SN\,2022oqm show a hot (40,000\,K) continuum with high-ionization C and O absorption features at velocities of 4000\,km\,s$^{-1}$, while its photospheric radius expands at 20,000\,\kms,… ▽ More We present the discovery and analysis of SN\,2022oqm, a Type Ic supernova (SN) detected $<1$\,day after explosion. The SN rises to a blue and short-lived (2\,days) initial peak. Early-time spectral observations of SN\,2022oqm show a hot (40,000\,K) continuum with high-ionization C and O absorption features at velocities of 4000\,km\,s$^{-1}$, while its photospheric radius expands at 20,000\,\kms, indicating a pre-existing distribution of expanding C/O material. After $\sim2.5$\,days, both the spectrum and light curves evolve into those of a typical SN Ic, with line velocities of $\sim10,000$\,km\,s$^{-1}$, in agreement with the photospheric radius evolution. The optical light curves reach a second peak at $t\approx15$\,days. By $t=60$\,days, the spectrum of \oqm\ becomes nearly nebular, displaying strong \ion{Ca}{2} and [\ion{Ca}{2}] emission with no detectable [\ion{O}{1}], marking this event as Ca-rich. The early behavior can be explained by $10^{-3}$\,\msun\ of optically thin circumstellar material (CSM) surrounding either (1) a massive compact progenitor such as a Wolf-Rayet star, (2) a massive stripped progenitor with an extended envelope, or (3) a binary system with a white dwarf. We propose that the early-time light curve is powered by both interaction of the ejecta with the optically thin CSM and shock cooling (in the massive-star scenario). The observations can be explained by CSM that is optically thick to X-ray photons, is optically thick in the lines as seen in the spectra, and is optically thin to visible-light continuum photons that come either from downscattered X-rays or from the shock-heated ejecta. Calculations show that this scenario is self-consistent. △ Less

Submitted 21 September, 2023; v1 submitted 5 October, 2022; originally announced October 2022.

Comments: 36 pages, 23 figures. Comments are welcome ido.irani@weizmann.ac.il or idoirani@gmail.com. Accepted to ApJ

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)

Zwicky Transient Facility and Globular Clusters: Calibration of the gr-Band Absolute Magnitudes for the Yellow Post-Asymptotic-Giant-Branch Stars

Authors: Chow-Choong Ngeow, Anupam Bhardwaj, Daniel Reiley, Russ R. Laher, Josiah Purdum, Ben Rusholme

Abstract: We present the first absolute calibration for the yellow post-asymptotic-giant-branch (PAGB) stars in the g- and r-band based on time-series observations from the Zwicky Transient Facility. These absolute magnitudes were calibrated using four yellow PAGB stars (one non-varying star and three Type II Cepheids) located in the globular clusters. We provide two calibrations of the gr-band absolute mag… ▽ More We present the first absolute calibration for the yellow post-asymptotic-giant-branch (PAGB) stars in the g- and r-band based on time-series observations from the Zwicky Transient Facility. These absolute magnitudes were calibrated using four yellow PAGB stars (one non-varying star and three Type II Cepheids) located in the globular clusters. We provide two calibrations of the gr-band absolute magnitudes for the yellow PAGB stars, by using an arithmetic mean and a linear regression. We demonstrate that the linear regression provides a better fit to the g-band absolute magnitudes for the yellow PAGB stars. These calibrated gr-band absolute magnitudes have a potential to be used as population II distance indicators in the era of time-domain synoptic sky surveys. △ Less

Submitted 15 August, 2022; originally announced August 2022.

Comments: 6 pages, 3 Figures and 1 Table; AJ accepted

Zwicky Transient Facility and Globular Clusters: The Period-Luminosity and Period-Wesenheit Relations for Type II Cepheids

Authors: Chow-Choong Ngeow, Anupam Bhardwaj, Jing-Yi Henderson, Matthew J. Graham, Russ R. Laher, Michael S. Medford, Josiah Purdum, Ben Rusholme

Abstract: We present the first gri-band period-luminosity (PL) and period-Wesenheit (PW) relations for 37 Type II Cepheids (hereafter TIIC) located in 18 globular clusters based on photometric data from the Zwicky Transient Facility. We also updated BV IJHK-band absolute magnitudes for 58 TIIC in 24 globular clusters using the latest homogeneous distances to the globular clusters. The slopes of g/r/i and B/… ▽ More We present the first gri-band period-luminosity (PL) and period-Wesenheit (PW) relations for 37 Type II Cepheids (hereafter TIIC) located in 18 globular clusters based on photometric data from the Zwicky Transient Facility. We also updated BV IJHK-band absolute magnitudes for 58 TIIC in 24 globular clusters using the latest homogeneous distances to the globular clusters. The slopes of g/r/i and B/V/I band PL relations are found to be statistically consistent when using the same sample of distance and reddening. We employed the calibration of ri-band PL/PW relations in globular clusters to estimate a distance to M31 based on a sample of ~270 TIIC from the PAndromeda project. The distance modulus to M31, obtained using calibrated ri-band PW relation, agrees well with the recent determination based on classical Cepheids. However, distance moduli derived using the calibrated r- and i-band PL relations are systematically smaller by ~0.2 mag, suggesting there are possible additional systematic error on the PL relations. Finally, we also derive the period-color (PC) relations and for the first time the period-Q-index (PQ) relations, where the Q-index is reddening-free, for our sample of TIIC. The PC relations based on (r-i) and near-infrared colors and the PQ relations are found to be relatively independent of the pulsation periods. △ Less

Submitted 5 August, 2022; originally announced August 2022.

Comments: 14 pages, 8 Figures and 5 Tables, AJ accepted

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

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

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

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

Submitted 25 February, 2022; originally announced February 2022.

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

Observing the Sun as a star: Design and early results from the NEID solar feed

Authors: Andrea S. J. Lin, Andrew Monson, Suvrath Mahadevan, Joe P. Ninan, Samuel Halverson, Colin Nitroy, Chad F. Bender, Sarah E. Logsdon, Shubham Kanodia, Ryan C. Terrien, Arpita Roy, Jacob K. Luhn, Arvind F. Gupta, Eric B. Ford, Fred Hearty, Russ R. Laher, Emily Hunting, William R. McBride, Noah Isaac Salazar Rivera, Jayadev Rajagopal, Marsha J. Wolf, Paul Robertson, Jason T. Wright, Cullen H. Blake, Caleb I. Canas , et al. (5 additional authors not shown)

Abstract: Efforts with extreme-precision radial velocity (EPRV) instruments to detect small-amplitude planets are largely limited, on many timescales, by the effects of stellar variability and instrumental systematics. One avenue for investigating these effects is the use of small solar telescopes which direct disk-integrated sunlight to these EPRV instruments, observing the Sun at high cadence over months… ▽ More Efforts with extreme-precision radial velocity (EPRV) instruments to detect small-amplitude planets are largely limited, on many timescales, by the effects of stellar variability and instrumental systematics. One avenue for investigating these effects is the use of small solar telescopes which direct disk-integrated sunlight to these EPRV instruments, observing the Sun at high cadence over months or years. We have designed and built a solar feed system to carry out "Sun-as-a-star" observations with NEID, a very high precision Doppler spectrometer recently commissioned at the WIYN 3.5m Telescope at Kitt Peak National Observatory. The NEID solar feed has been taking observations nearly every day since December 2020; data is publicly available at the NASA Exoplanet Science Institute (NExScI) NEID Solar Archive: \url{https://neid.ipac.caltech.edu/search_solar.php}. In this paper, we present the design of the NEID solar feed and explanations behind our design intent. We also present early radial velocity (RV) results which demonstrate NEID's RV stability on the Sun over 4 months of commissioning: 0.66~m/s RMS under good sky conditions and improving to 0.41~m/s RMS under best conditions. △ Less

Submitted 15 February, 2022; v1 submitted 10 December, 2021; originally announced December 2021.

Comments: Accepted to AJ

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

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)

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

Faintest of them all : ZTF 21aaoryiz/SN 2021fcg -- Discovery of an extremely low luminosity type Iax supernova

Authors: Viraj R. Karambelkar, Mansi M. Kasliwal, Kate Maguire, Shreya G. Anand, Igor Andreoni, Kishalay De, Andrew Drake, Dmitry A. Duev, Matthew J. Graham, Erik C. Kool, Russ R. Laher, Mark R. Magee, Ashish A. Mahabal, Michael S. Medford, Daniel Perley, Mickael Rigault, Ben Rusholme, Steve Schulze, Yashvi Sharma, Jesper Sollerman, Anastasios Tzanidakis, Richard Walters, Yuhan Yao

Abstract: We present the discovery of ZTF 21aaoryiz/SN 2021fcg -- an extremely low-luminosity Type Iax supernova. SN 2021fcg was discovered by the Zwicky Transient Facility in the star-forming galaxy IC0512 at a distance of $\approx$ 27 Mpc. It reached a peak absolute magnitude of $M_{r} =$ $-12.66\pm0.20$ mag, making it the least luminous thermonuclear supernova discovered to date. The E(B-V) contribution… ▽ More We present the discovery of ZTF 21aaoryiz/SN 2021fcg -- an extremely low-luminosity Type Iax supernova. SN 2021fcg was discovered by the Zwicky Transient Facility in the star-forming galaxy IC0512 at a distance of $\approx$ 27 Mpc. It reached a peak absolute magnitude of $M_{r} =$ $-12.66\pm0.20$ mag, making it the least luminous thermonuclear supernova discovered to date. The E(B-V) contribution from the underlying host galaxy is unconstrained. However, even if it were as large as 0.5 mag, the peak absolute magnitude would be $M_{r} = -13.78\pm0.20$ mag -- still consistent with being the lowest luminosity SN. Optical spectra of SN 2021fcg taken at 37 and 65 days post maximum show strong [Ca II], Ca II and Na I D emission and several weak [Fe II] emission lines. The [Ca II] emission in the two spectra has extremely low velocities of $\approx 1300$ and $1000$ km s$^{-1}$ respectively. The spectra very closely resemble those of the very low luminosity Type Iax supernovae SN 2008ha, SN 2010ae and SN 2019gsc taken at similar phases. The peak bolometric luminosity of SN 2021fcg is $\approx$ $2.5^{+1.5}_{-0.3}\times10^{40}$ erg s$^{-1}$ which is a factor of three lower than that for SN 2008ha. The bolometric lightcurve of SN 2021fcg is consistent with a very low ejected nickel mass (M$_{\rm{Ni}} \approx 0.8^{+0.4}_{-0.5}\times10^{-3}$ M$_{\odot}$). The low luminosity and nickel mass of SN 2021fcg pose a challenge to the picture that low luminosity SNe Iax originate from deflagrations of near M$_{\rm{ch}}$ hybrid carbon-oxygen-neon white dwarfs. Instead, the merger of a carbon-oxygen and oxygen-neon white dwarf is a promising model to explain SN 2021fcg. △ Less

Submitted 8 October, 2021; originally announced October 2021.

Comments: accepted for publication in ApJL. 10 pages, 5 figures

Discovery and characterization of five new eclipsing AM CVn systems

Authors: J. van Roestel, T. Kupfer, M. J. Green, S. Wong, L. Bildsten, K. Burdge, T. Prince, T. R. Marsh, P. Szkody, C. Fremling, M. J. Graham, V. S. Dhillon, S. P. Littlefair, E. C. Bellm, M. Coughlin, D. A. Duev, D. A. Goldstein, R. R. Laher, B. Rusholme, R. Riddle, R. Dekany, S. R. Kulkarni

Abstract: AM CVn systems are ultra-compact, helium-rich, accreting binaries with degenerate or semi-degenerate donors. We report the discovery of five new eclipsing AM CVn systems with orbital periods of 61.5, 55.5, 53.3, 37.4, and 35.4 minutes. These systems were discovered by searching for deep eclipses in the Zwicky Transient Facility (ZTF) lightcurves of white dwarfs selected using Gaia parallaxes. We o… ▽ More AM CVn systems are ultra-compact, helium-rich, accreting binaries with degenerate or semi-degenerate donors. We report the discovery of five new eclipsing AM CVn systems with orbital periods of 61.5, 55.5, 53.3, 37.4, and 35.4 minutes. These systems were discovered by searching for deep eclipses in the Zwicky Transient Facility (ZTF) lightcurves of white dwarfs selected using Gaia parallaxes. We obtained phase-resolved spectroscopy to confirm that all systems are AM CVn binaries, and we obtained high-speed photometry to confirm the eclipse and characterize the systems. The spectra of two long-period systems (61.5 and 53.3 minutes) show many emission and absorption lines, indicating the presence of N, O, Na, Mg, Si, and Ca, and also the K and Zn, elements which have never been detected in AM CVn systems before. By modelling the high-speed photometry, we measured the mass and radius of the donor star, potentially constraining the evolutionary channel that formed these AM CVn systems. We determined that the average mass of the accreting white dwarf is $\approx0.8$$\mathrm{M_{\odot}}$, and that the white dwarfs in long-period systems are hotter than predicted by recently updated theoretical models. The donors have a high entropy and are a factor of $\approx$ 2 more massive compared to zero-entropy donors at the same orbital period. The large donor radius is most consistent with He-star progenitors, although the observed spectral features seem to contradict this. The discovery of 5 new eclipsing AM~CVn systems is consistent with the known observed AM CVn space density and estimated ZTF recovery efficiency. Based on this estimate, we expect to find another 1--4 eclipsing AM CVn systems as ZTF continues to obtain data. This will further increase our understanding of the population, but will require high precision data to better characterize these 5 systems and any new discoveries. △ Less

Submitted 30 August, 2021; v1 submitted 15 July, 2021; originally announced July 2021.

Comments: submitted to MNRAS

Inferring kilonova population properties with a hierarchical Bayesian framework I : Non-detection methodology and single-event analyses

Authors: Siddharth R. Mohite, Priyadarshini Rajkumar, Shreya Anand, David L. Kaplan, Michael W. Coughlin, Ana Sagués-Carracedo, Muhammed Saleem, Jolien Creighton, Patrick R. Brady, Tomás Ahumada, Mouza Almualla, Igor Andreoni, Mattia Bulla, Matthew J. Graham, Mansi M. Kasliwal, Stephen Kaye, Russ R. Laher, Kyung Min Shin, David L. Shupe, Leo P. Singer

Abstract: We present ${\tt nimbus}$ : a hierarchical Bayesian framework to infer the intrinsic luminosity parameters of kilonovae (KNe) associated with gravitational-wave (GW) events, based purely on non-detections. This framework makes use of GW 3-D distance information and electromagnetic upper limits from multiple surveys for multiple events, and self-consistently accounts for finite sky-coverage and pro… ▽ More We present ${\tt nimbus}$ : a hierarchical Bayesian framework to infer the intrinsic luminosity parameters of kilonovae (KNe) associated with gravitational-wave (GW) events, based purely on non-detections. This framework makes use of GW 3-D distance information and electromagnetic upper limits from multiple surveys for multiple events, and self-consistently accounts for finite sky-coverage and probability of astrophysical origin. The framework is agnostic to the brightness evolution assumed and can account for multiple electromagnetic passbands simultaneously. Our analyses highlight the importance of accounting for model selection effects, especially in the context of non-detections. We show our methodology using a simple, two-parameter linear brightness model, taking the follow-up of GW190425 with the Zwicky Transient Facility (ZTF) as a single-event test case for two different prior choices of model parameters -- (i) uniform/uninformative priors and (ii) astrophysical priors based on surrogate models of Monte Carlo radiative transfer simulations of KNe. We present results under the assumption that the KN is within the searched region to demonstrate functionality and the importance of prior choice. Our results show consistency with ${\tt simsurvey}$ -- an astronomical survey simulation tool used previously in the literature to constrain the population of KNe. While our results based on uniform priors strongly constrain the parameter space, those based on astrophysical priors are largely uninformative, highlighting the need for deeper constraints. Future studies with multiple events having electromagnetic follow-up from multiple surveys should make it possible to constrain the KN population further. △ Less

Submitted 15 November, 2021; v1 submitted 15 July, 2021; originally announced July 2021.

Comments: 18 pages, 5 figures, accepted for publication in ApJ

Year 1 of the ZTF high-cadence Galactic Plane Survey: Strategy, goals, and early results on new single-mode hot subdwarf B-star pulsators

Authors: Thomas Kupfer, Thomas A. Prince, Jan van Roestel, Eric C. Bellm, Lars Bildsten, Michael W. Coughlin, Andrew J. Drake, Matthew J. Graham, Courtney Klein, Shrinivas R. Kulkarni, Frank J. Masci, Richard Walters, Igor Andreoni, Rahul Biswas, Corey Bradshaw, Dmitry A. Duev, Richard Dekany, Joseph A. Guidry, JJ Hermes, Russ R. Laher, Reed Riddle

Abstract: We present the goals, strategy and first results of the high-cadence Galactic plane survey using the Zwicky Transient Facility (ZTF). The goal of the survey is to unveil the Galactic population of short-period variable stars, including short period binaries and stellar pulsators with periods less than a few hours. Between June 2018 and January 2019, we observed 64 ZTF fields resulting in 2990 deg… ▽ More We present the goals, strategy and first results of the high-cadence Galactic plane survey using the Zwicky Transient Facility (ZTF). The goal of the survey is to unveil the Galactic population of short-period variable stars, including short period binaries and stellar pulsators with periods less than a few hours. Between June 2018 and January 2019, we observed 64 ZTF fields resulting in 2990 deg$^2$ of high stellar density in ZTF-$r$ band along the Galactic Plane. Each field was observed continuously for 1.5 to 6 hrs with a cadence of 40 sec. Most fields have between 200 and 400 observations obtained over 2-3 continuous nights. As part of this survey we extract a total of $\approx$230 million individual objects with at least 80 epochs obtained during the high-cadence Galactic Plane survey reaching an average depth of ZTF-$r$ $\approx$20.5 mag. For four selected fields with 2 million to 10 million individual objects per field we calculate different variability statistics and find that $\approx$1-2% of the objects are astrophysically variable over the observed period. We present a progress report on recent discoveries, including a new class of compact pulsators, the first members of a new class of Roche Lobe filling hot subdwarf binaries as well as new ultracompact double white dwarfs and flaring stars. Finally we present a sample of 12 new single-mode hot subdwarf B-star pulsators with pulsation amplitudes between ZTF-$r$ = 20-76 mmag and pulsation periods between $P$ = 5.8-16 min with a strong cluster of systems with periods $\approx$ 6 min. All of the data have now been released in either ZTF Data Release 3 or data release 4. △ Less

Submitted 6 May, 2021; originally announced May 2021.

Comments: accepted for publication in MNRAS, 16 pages, 10 figures and 5 tables

A systematic search for outbursting AM CVn systems with the Zwicky Transient Facility

Authors: Jan van Roestel, Leah Creter, Thomas Kupfer, Paula Szkody, Jim Fuller, Matthew J. Green, R. Michael Rich, John Sepikas, Kevin Burdge, Ilaria Caiazzo, Przemek Mroz, Thomas A. Prince, Dmitry A. Duev, Matthew J. Graham, David L. Shupe, Russ R. Laher, Ashish A. Mahabal, Frank J. Masci

Abstract: AM CVn systems are a rare type of accreting binary that consists of a white dwarf and a helium-rich, degenerate donor star. Using the Zwicky Transient Facility (ZTF), we searched for new AM CVn systems by focusing on blue, outbursting stars. We first selected outbursting stars using the ZTF alerts. We cross-matched the candidates with $Gaia$ and Pan-STARRS catalogs. The initial selection of candid… ▽ More AM CVn systems are a rare type of accreting binary that consists of a white dwarf and a helium-rich, degenerate donor star. Using the Zwicky Transient Facility (ZTF), we searched for new AM CVn systems by focusing on blue, outbursting stars. We first selected outbursting stars using the ZTF alerts. We cross-matched the candidates with $Gaia$ and Pan-STARRS catalogs. The initial selection of candidates based on the $Gaia$ $BP$-$RP$ contains 1751 unknown objects. We used the Pan-STARRS $g$-$r$ and $r$-$i$ color in combination with the $Gaia$ color to identify 59 high-priority candidates. We obtained identification spectra of 35 sources, of which 18 are high priority candidates, and discovered 9 new AM CVn systems and one magnetic CV which shows only He-II lines. Using the outburst recurrence time, we estimate the orbital periods which are in the range of 29 to 50 minutes. We conclude that targeted followup of blue, outbursting sources is an efficient method to find new AM CVn systems, and we plan to followup all candidates we identified to systematically study the population of outbursting AM CVn systems. △ Less

Submitted 7 May, 2021; v1 submitted 5 May, 2021; originally announced May 2021.

Comments: Submitted, comments and feedback welcome

Fast-transient Searches in Real Time with ZTFReST: Identification of Three Optically-discovered Gamma-ray Burst Afterglows and New Constraints on the Kilonova Rate

Authors: Igor Andreoni, Michael W. Coughlin, Erik C. Kool, Mansi M. Kasliwal, Harsh Kumar, Varun Bhalerao, Ana Sagués Carracedo, Anna Y. Q. Ho, Peter T. H. Pang, Divita Saraogi, Kritti Sharma, Vedant Shenoy, Eric Burns, Tomás Ahumada, Shreya Anand, Leo P. Singer, Daniel A. Perley, Kishalay De, U. C. Fremling, Eric C. Bellm, Mattia Bulla, Arien Crellin-Quick, Tim Dietrich, Andrew Drake, Dmitry A. Duev , et al. (10 additional authors not shown)

Abstract: While optical surveys regularly discover slow transients like supernovae on their own, the most common way to discover extragalactic fast transients, fading away in a few nights, is via follow-up observations of gamma-ray burst and gravitational-wave triggers. However, wide-field surveys have the potential to also identify rapidly fading transients independently of such external triggers. The volu… ▽ More While optical surveys regularly discover slow transients like supernovae on their own, the most common way to discover extragalactic fast transients, fading away in a few nights, is via follow-up observations of gamma-ray burst and gravitational-wave triggers. However, wide-field surveys have the potential to also identify rapidly fading transients independently of such external triggers. The volumetric survey speed of the Zwicky Transient Facility (ZTF) makes it sensitive to faint and fast-fading objects as kilonovae, the optical counterparts to binary neutron stars and neutron star-black hole mergers, out to almost 200Mpc. We introduce an open-source software infrastructure, the ZTF REaltime Search and Triggering, ZTFReST, designed to identify kilonovae and fast optical transients in ZTF data. Using the ZTF alert stream combined with forced photometry, we have implemented automated candidate ranking based on their photometric evolution and fitting to kilonova models. Automated triggering of follow-up systems, such as Las Cumbres Observatory, has also been implemented. In 13 months of science validation, we found several extragalactic fast transients independent of any external trigger (though some counterparts were identified later), including at least one supernova with post-shock cooling emission, two known afterglows with an associated gamma-ray burst, two known afterglows without any known gamma-ray counterpart, and three new fast-declining sources (ZTF20abtxwfx, ZTF20acozryr, and ZTF21aagwbjr) that are likely associated with GRB200817A, GRB201103B, and GRB210204A. However, we have not found any objects which appear to be kilonovae; therefore, we constrain the rate of GW170817-like kilonovae to $R < 900$Gpc$^{-3}$yr$^{-1}$. A framework such as ZTFReST could become a prime tool for kilonova and fast transient discovery with the Vera C. Rubin Observatory. △ Less

Submitted 13 April, 2021; originally announced April 2021.

Comments: Submitted for publication

Real-time Discovery of AT2020xnd: A Fast, Luminous Ultraviolet Transient with Minimal Radioactive Ejecta

Authors: Daniel A. Perley, Anna Y. Q. Ho, Yuhan Yao, Christoffer Fremling, Joseph P. Anderson, Steve Schulze, Harsh Kumar, G. C. Anupama, Sudhanshu Barway, Eric C. Bellm, Varun Bhalerao, Ting-Wan Chen, Dmitry A. Duev, Lluís Galbany, Matthew J. Graham, Mariusz Gromadzki, Claudia P. Gutiérrez, Nada Ihanec, Cosimo Inserram, Mansi M. Kasliwal, Erik C. Kool, S. R. Kulkarni, Russ R. Laher, Frank J. Masci, James D. Neill , et al. (7 additional authors not shown)

Abstract: The many unusual properties of the enigmatic AT2018cow suggested that at least some subset of the empirical class of fast blue optical transients (FBOTs) represents a genuinely new astrophysical phenomenon. Unfortunately, the intrinsic rarity and fleeting nature of these events have made it difficult to identify additional examples early enough to acquire the observations necessary to constrain th… ▽ More The many unusual properties of the enigmatic AT2018cow suggested that at least some subset of the empirical class of fast blue optical transients (FBOTs) represents a genuinely new astrophysical phenomenon. Unfortunately, the intrinsic rarity and fleeting nature of these events have made it difficult to identify additional examples early enough to acquire the observations necessary to constrain theoretical models. We present here the Zwicky Transient Facility discovery of AT2020xnd (ZTF20acigmel, the "Camel") at z=0.243, the first unambiguous AT2018cow analog to be found and confirmed in real time. AT2018cow and AT2020xnd share all key observational properties: a fast optical rise, sustained high photospheric temperature, absence of a second peak attributable to ejection of a radioactively-heated stellar envelope, extremely luminous radio, millimetre, and X-ray emission, and a dwarf-galaxy host. This supports the argument that AT2018cow-like events represent a distinct phenomenon from slower-evolving radio-quiet supernovae, likely requiring a different progenitor or a different central engine. The sample properties of the four known members of this class to date disfavour tidal disruption models but are consistent with the alternative model of an accretion powered jet following the direct collapse of a massive star to a black hole. Contextual filtering of alert streams combined with rapid photometric verification using multi-band imaging provides an efficient way to identify future members of this class, even at high redshift. △ Less

Submitted 14 October, 2021; v1 submitted 2 March, 2021; originally announced March 2021.

Comments: Accepted to MNRAS

Time-series and Phasecurve Photometry of Episodically-Active Asteroid (6478) Gault in a Quiescent State Using APO, GROWTH, P200 and ZTF

Authors: Josiah N. Purdum, Zhong-Yi Lin, Bryce T. Bolin, Kritti Sharma, Philip I. Choi, Varun Bhalerao, Harsh Kumar, Robert Quimby, Joannes C. Van Roestel, Chengxing Zhai, Yanga R. Fernandez, Josef Hanuš, Carey M. Lisse, Dennis Bodewits, Christoffer Fremling, Nathan Ryan Golovich, Chen-Yen Hsu, Wing-Huen Ip, Chow-Choong Ngeow, Navtej S. Saini, Michael Shao, Yuhan Yao, Tomás Ahumada, Shreya Anand, Igor Andreoni , et al. (27 additional authors not shown)

Abstract: We observed Episodically Active Asteroid (6478) Gault in 2020 with multiple telescopes in Asia and North America and have found that it is no longer active after its recent outbursts at the end of 2018 and start of 2019. The inactivity during this apparation allowed us to measure the absolute magnitude of Gault of H_r = 14.63 +/- 0.02, G_r = 0.21 +/- 0.02 from our secular phasecurve observations.… ▽ More We observed Episodically Active Asteroid (6478) Gault in 2020 with multiple telescopes in Asia and North America and have found that it is no longer active after its recent outbursts at the end of 2018 and start of 2019. The inactivity during this apparation allowed us to measure the absolute magnitude of Gault of H_r = 14.63 +/- 0.02, G_r = 0.21 +/- 0.02 from our secular phasecurve observations. In addition, we were able to constrain Gault's rotation period using time-series photometric lightcurves taken over 17 hours on multiple days in 2020 August, September and October. The photometric lightcurves have a repeating $\lesssim$0.05 magnitude feature suggesting that (6478) Gault has a rotation period of ~2.5 hours and may have a semi-spherical or top-like shape, much like Near-Earth Asteroids Ryugu and Bennu. The rotation period of ~2.5 hours is near to the expected critical rotation period for an asteroid with the physical properties of (6478) Gault suggesting that its activity observed over multiple epochs is due to surface mass shedding from its fast rotation spun up by the Yarkovsky-O'Keefe-Radzievskii-Paddack effect. △ Less

Submitted 27 March, 2021; v1 submitted 25 February, 2021; originally announced February 2021.

Comments: 23 pages, 18 figures; Accepted by ApJ Letters

The ZTF Source Classification Project: I. Methods and Infrastructure

Authors: Jan van Roestel, Dmitry A. Duev, Ashish A. Mahabal, Michael W. Coughlin, Przemek Mróz, Kevin Burdge, Andrew Drake, Matthew J. Graham, Lynne Hillenbrand, C. Fremling, David Hale, Russ R. Laher, Frank J. Masci, Reed Riddle, Philippe Rosnet, Ben Rusholme, Roger Smith, Maayane T. Soumagnac, Richard Walters, Thomas A. Prince, S. R. Kulkarni

Abstract: The Zwicky Transient Facility (ZTF) has been observing the entire northern sky since the start of 2018 down to a magnitude of 20.5 ($5 σ$ for 30s exposure) in $g$, $r$, and $i$ filters. Over the course of two years, ZTF has obtained light curves of more than a billion sources, each with 50-1000 epochs per light curve in $g$ and $r$, and fewer in $i$. To be able to use the information contained in… ▽ More The Zwicky Transient Facility (ZTF) has been observing the entire northern sky since the start of 2018 down to a magnitude of 20.5 ($5 σ$ for 30s exposure) in $g$, $r$, and $i$ filters. Over the course of two years, ZTF has obtained light curves of more than a billion sources, each with 50-1000 epochs per light curve in $g$ and $r$, and fewer in $i$. To be able to use the information contained in the light curves of variable sources for new scientific discoveries, an efficient and flexible framework is needed to classify them. In this paper, we introduce the methods and infrastructure which will be used to classify all ZTF light curves. Our approach aims to be flexible and modular and allows the use of a dynamical classification scheme and labels, continuously evolving training sets, and the use of different machine learning classifier types and architectures. With this setup, we are able to continuously update and improve the classification of ZTF light curves as new data becomes available, training samples are updated, and new classes need to be incorporated. △ Less

Submitted 22 February, 2021; originally announced February 2021.

HO Puppis: Not a Be Star but a Newly Confirmed IW And-Type Star

Authors: Chien-De Lee, Jia-Yu Ou, Po-Chieh Yu, Chow-Choong Ngeow, Po-Chieh Huang, Wing-Huen Ip, Franz-Josef Hambsch, Hyun-il Sung, Jan van Roestel, Richard Dekany, Andrew J. Drake, Matthew J. Graham, Dmitry A. Duev, Stephen Kaye, Thomas Kupfer, Russ R. Laher, Frank J. Masci, Przemek Mroz, James D. Neill, Reed Riddle, Ben Rusholme, Richard Walters

Abstract: HO Puppis (HO Pup) was considered as a Be-star candidate based on its gamma-Cassiopeiae-type light curve, but lacked spectroscopic confirmation. Using distance measured from Gaia Data Release 2 and the spectral-energy-distribution (SED) fit on broadband photometry, the Be-star nature of HO Pup is ruled out. Furthermore, based on the 28,700 photometric data points collected from various time-domain… ▽ More HO Puppis (HO Pup) was considered as a Be-star candidate based on its gamma-Cassiopeiae-type light curve, but lacked spectroscopic confirmation. Using distance measured from Gaia Data Release 2 and the spectral-energy-distribution (SED) fit on broadband photometry, the Be-star nature of HO Pup is ruled out. Furthermore, based on the 28,700 photometric data points collected from various time-domain surveys and dedicated intensive-monitoring observations, the light curves of HO Pup closely resemble IW And-type stars (as pointed out in Kimura et al. 2020a), exhibiting characteristics such as quasi-standstill phase, brightening, and dips. The light curve of HO Pup displays various variability timescales, including brightening cycles ranging from 23 to 61 days, variations with periods between 3.9 days and 50 minutes during the quasi-standstill phase, and a semi-regular ~14-day period for the dip events. We have also collected time-series spectra (with various spectral resolutions), in which Balmer emission lines and other expected spectral lines for an IW And-type star were detected (even though some of these lines were also expected to be present for Be stars). We detect Bowen fluorescence near the brightening phase, and that can be used to discriminate between IW And-type stars and Be stars. Finally, despite only observing for four nights, the polarization variation was detected, indicating that HO Pup has significant intrinsic polarization. △ Less

Submitted 19 February, 2021; originally announced February 2021.

Comments: 24 pages, 7 tables and 12 figures; ApJ accepted

Initial Characterization of Active Transitioning Centaur, P/2019 LD2 (ATLAS), using Hubble, Spitzer, ZTF, Keck, APO and GROWTH Visible & Infrared Imaging and Spectroscopy

Authors: Bryce T. Bolin, Yanga R. Fernandez, Carey M. Lisse, Timothy R. Holt, Zhong-Yi Lin, Josiah N. Purdum, Kunal P. Deshmukh, James M. Bauer, Eric C. Bellm, Dennis Bodewits, Kevin B. Burdge, Sean J. Carey, Chris M. Copperwheat, George Helou, Anna Y. Q. Ho, Jonathan Horner, Jan van Roestel, Varun Bhalerao, Chan-Kao Chang, Christine Chen, Chen-Yen Hsu, Wing-Huen Ip, Mansi M. Kasliwal, Frank J. Masci, Chow-Choong Ngeow , et al. (21 additional authors not shown)

Abstract: We present visible and mid-infrared imagery and photometry of temporary Jovian co-orbital comet P/2019 LD$_2$ taken with HST/WFC3, Spitzer/IRAC, the GROWTH telescope network, visible spectroscopy from Keck/LRIS and archival ZTF observations taken between 2019 April and 2020 August. Our observations indicate that the nucleus of LD$_2$ has a radius between 0.2-1.8 km assuming a 0.08 albedo and a com… ▽ More We present visible and mid-infrared imagery and photometry of temporary Jovian co-orbital comet P/2019 LD$_2$ taken with HST/WFC3, Spitzer/IRAC, the GROWTH telescope network, visible spectroscopy from Keck/LRIS and archival ZTF observations taken between 2019 April and 2020 August. Our observations indicate that the nucleus of LD$_2$ has a radius between 0.2-1.8 km assuming a 0.08 albedo and a coma dominated by $\sim$100$μ$ m-scale dust ejected at $\sim$1 m/s speeds with a $\sim$1'' jet pointing in the SW direction. LD$_2$ experienced a total dust mass loss of $\sim$10$^8$ kg at a loss rate of $\sim$6 kg/s with Af$ρ$/cross-section varying between $\sim$85 cm/125 km$^2$ and $\sim$200 cm/310 km$^2$ from 2019 April 9 to 2019 Nov 8. If the increase in Af$ρ$/cross-section remained constant, it implies LD$_2$'s activity began $\sim$2018 November when within 4.8 au of the Sun, implying the onset of H$_2$O sublimation. We measure CO/CO$_2$ gas production of $\lesssim$10$^{27}$ mol/s /$\lesssim$10$^{26}$ mol/s from our 4.5 $μ$m Spitzer observations, $g$-$r$ = 0.59$\pm$0.03, $r$-$i$ = 0.18$\pm$0.05, $i$-$z$ = 0.01$\pm$0.07 from GROWTH observations, H$_2$O gas production of $\lesssim$80 kg/s scaling from our estimated $C_2$ production of $Q_{C_2}\lesssim$7.5$\times10^{24}$ mol/s from Keck/LRIS spectroscopy. We determine that the long-term orbit of LD$_2$ is similar to Jupiter family comets having close encounters with Jupiter within $\sim$0.5 Hill radius in the last $\sim$3 y, within 0.8 Hill radius in $\sim$9 y. Additionally, 78.8$\%$ of our orbital clones are ejected from the Solar System within $1 \times 10^{6}$ years having a dynamical half-life of 3.4 $\times 10^5$ years. △ Less

Submitted 5 January, 2021; v1 submitted 7 November, 2020; originally announced November 2020.

Comments: 28 pages, 8 figures, 2 Tables, accepted for publication in AJ

TDE Hosts are Green and Centrally Concentrated: Signatures of a Post-Merger System

Authors: Erica Hammerstein, Suvi Gezari, Sjoert van Velzen, S. Bradley Cenko, Nathaniel Roth, Charlotte Ward, Sara Frederick, Tiara Hung, Matthew Graham, Ryan J. Foley, Eric C. Bellm, Christopher Cannella, Andrew J. Drake, Thomas Kupfer, Russ R. Laher, Ashish A. Mahabal, Frank J. Masci, Reed Riddle, César Rojas-Bravo, Roger Smith

Abstract: We study the properties of the galaxies hosting the first 19 tidal disruption events (TDEs) detected with the Zwicky Transient Facility (ZTF) within the context of a carefully constructed, representative host galaxy sample. We find that the ZTF sample of TDE hosts is dominated by compact "green valley" galaxies. After we restrict the comparison sample to galaxies with a similar concentration, as m… ▽ More We study the properties of the galaxies hosting the first 19 tidal disruption events (TDEs) detected with the Zwicky Transient Facility (ZTF) within the context of a carefully constructed, representative host galaxy sample. We find that the ZTF sample of TDE hosts is dominated by compact "green valley" galaxies. After we restrict the comparison sample to galaxies with a similar concentration, as measured by Sersic index, we find this green valley over representation is even larger. That is, concentrated red sequence galaxies are not producing TDEs at elevated levels. We present host galaxy spectra which show that E+A galaxies are overrepresented in the ZTF sample by a factor of $\approx$22, which is lower than previous TDE host galaxy studies have found. We find that this overrepresentation can be fully accounted for when taking into account the masses, colors, and Sérsic indices of the ZTF TDE hosts. The combination of both green colors and high Sérsic index of the typical TDE host galaxy could be explained if the TDE rate is temporarily enhanced following a merger that leads to a higher central concentration of stars. △ Less

Submitted 21 October, 2020; v1 submitted 20 October, 2020; originally announced October 2020.

Comments: 12 pages, 5 figures, submitted to ApJL

An 8.8 minute orbital period eclipsing detached double white dwarf binary

Authors: Kevin B. Burdge, Michael W. Coughlin, Jim Fuller, David L. Kaplan, S. R. Kulkarni, Thomas R. Marsh, Thomas A. Prince, Eric C. Bellm, Richard G. Dekany, Dmitry A. Duev, Matthew J. Graham, Ashish A. Mahabal, Frank J. Masci, Russ R. Laher, Reed Riddle, Maayane T. Soumagnac

Abstract: We report the discovery of ZTF J2243+5242, an eclipsing double white dwarf binary with an orbital period of just $8.8$ minutes, the second known eclipsing binary with an orbital period less than ten minutes. The system likely consists of two low-mass white dwarfs, and will merge in approximately 400,000 years to form either an isolated hot subdwarf or an R Coronae Borealis star. Like its… ▽ More We report the discovery of ZTF J2243+5242, an eclipsing double white dwarf binary with an orbital period of just $8.8$ minutes, the second known eclipsing binary with an orbital period less than ten minutes. The system likely consists of two low-mass white dwarfs, and will merge in approximately 400,000 years to form either an isolated hot subdwarf or an R Coronae Borealis star. Like its $6.91\, \rm min$ counterpart, ZTF J1539+5027, ZTF J2243+5242 will be among the strongest gravitational wave sources detectable by the space-based gravitational-wave detector The Laser Space Interferometer Antenna (LISA) because its gravitational-wave frequency falls near the peak of LISA's sensitivity. Based on its estimated distance of $d=2120^{+131}_{-115}\,\rm pc$, LISA should detect the source within its first few months of operation, and should achieve a signal-to-noise ratio of $87\pm5$ after four years. We find component masses of $M_A= 0.349^{+0.093}_{-0.074}\,M_\odot$ and $M_B=0.384^{+0.114}_{-0.074}\,M_\odot$, radii of $R_A=0.0308^{+0.0026}_{-0.0025}\,R_\odot$ and $R_B = 0.0291^{+0.0032}_{-0.0024}\,R_\odot$, and effective temperatures of $T_A=22200^{+1800}_{-1600}\,\rm K$ and $T_B=16200^{+1200}_{-1000}\,\rm K$. We determined all of these properties, and the distance to this system, using only photometric measurements, demonstrating a feasible way to estimate parameters for the large population of optically faint ($r>21 \, m_{\rm AB}$) gravitational-wave sources which the Vera Rubin Observatory (VRO) and LISA should identify. △ Less

Submitted 9 October, 2020; v1 submitted 7 October, 2020; originally announced October 2020.

Comments: 12 pages, 7 figures, 3 tables, submitted