Showing 1–50 of 434 results for author: Bloom, J

Interpreting Attention Layer Outputs with Sparse Autoencoders

Authors: Connor Kissane, Robert Krzyzanowski, Joseph Isaac Bloom, Arthur Conmy, Neel Nanda

Abstract: Decomposing model activations into interpretable components is a key open problem in mechanistic interpretability. Sparse autoencoders (SAEs) are a popular method for decomposing the internal activations of trained transformers into sparse, interpretable features, and have been applied to MLP layers and the residual stream. In this work we train SAEs on attention layer outputs and show that also h… ▽ More Decomposing model activations into interpretable components is a key open problem in mechanistic interpretability. Sparse autoencoders (SAEs) are a popular method for decomposing the internal activations of trained transformers into sparse, interpretable features, and have been applied to MLP layers and the residual stream. In this work we train SAEs on attention layer outputs and show that also here SAEs find a sparse, interpretable decomposition. We demonstrate this on transformers from several model families and up to 2B parameters. We perform a qualitative study of the features computed by attention layers, and find multiple families: long-range context, short-range context and induction features. We qualitatively study the role of every head in GPT-2 Small, and estimate that at least 90% of the heads are polysemantic, i.e. have multiple unrelated roles. Further, we show that Sparse Autoencoders are a useful tool that enable researchers to explain model behavior in greater detail than prior work. For example, we explore the mystery of why models have so many seemingly redundant induction heads, use SAEs to motivate the hypothesis that some are long-prefix whereas others are short-prefix, and confirm this with more rigorous analysis. We use our SAEs to analyze the computation performed by the Indirect Object Identification circuit (Wang et al.), validating that the SAEs find causally meaningful intermediate variables, and deepening our understanding of the semantics of the circuit. We open-source the trained SAEs and a tool for exploring arbitrary prompts through the lens of Attention Output SAEs. △ Less

Submitted 25 June, 2024; originally announced June 2024.

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 Gravity Collective: A Comprehensive Analysis of the Electromagnetic Search for the Binary Neutron Star Merger GW190425

Authors: D. A. Coulter, C. D. Kilpatrick, D. O. Jones, R. J. Foley, A. V. Filippenko, W. Zheng, J. J. Swift, G. S. Rahman, H. E. Stacey, A. L. Piro, C. Rojas-Bravo, J. Anais Vilchez, N. Muñoz-Elgueta, I. Arcavi, G. Dimitriadis, M. R. Siebert, J. S. Bloom, M. J. Bustamante-Rosell, K. E. Clever, K. W. Davis, J. Kutcka, P. Macias, P. McGill, P. J. Quiñonez, E. Ramirez-Ruiz , et al. (12 additional authors not shown)

Abstract: We present an ultraviolet-to-infrared search for the electromagnetic (EM) counterpart to GW190425, the second-ever binary neutron star (BNS) merger discovered by the LIGO-Virgo-KAGRA Collaboration (LVK). GW190425 was more distant and had a larger localization area than GW170817, therefore we use a new tool teglon to redistribute the GW190425 localization probability in the context of galaxy catalo… ▽ More We present an ultraviolet-to-infrared search for the electromagnetic (EM) counterpart to GW190425, the second-ever binary neutron star (BNS) merger discovered by the LIGO-Virgo-KAGRA Collaboration (LVK). GW190425 was more distant and had a larger localization area than GW170817, therefore we use a new tool teglon to redistribute the GW190425 localization probability in the context of galaxy catalogs within the final localization volume. We derive a 90th percentile area of 6,688 deg$^{2}$, a $\sim$1.5$\times$ improvement relative to the LIGO/Virgo map, and show how teglon provides an order of magnitude boost to the search efficiency of small ($\leq$1 deg$^{2}$) field-of-view instruments. We combine our data with all publicly reported imaging data, covering 9,078.59 deg$^2$ of unique area and 48.13% of the LIGO/Virgo-assigned localization probability, to calculate the most comprehensive kilonova, short gamma-ray burst (sGRB) afterglow, and model-independent constraints on the EM emission from a hypothetical counterpart to GW190425 to date under the assumption that no counterpart was found in these data. If the counterpart were similar to AT 2017gfo, there was a 28.4% chance that it would have been detected in the combined dataset. We are relatively insensitive to an on-axis sGRB, and rule out a generic transient with a similar peak luminosity and decline rate as AT 2017gfo to 30% confidence. Finally, across our new imaging and all publicly-reported data, we find 28 candidate optical counterparts that we cannot rule out as being associated with GW190425, finding that 4 such counterparts discovered within the localization volume and within 5 days of merger exhibit luminosities consistent with a kilonova. △ Less

Submitted 23 April, 2024; originally announced April 2024.

Comments: 41 pages, 11 figures, Submitted to ApJ

Testing a Galactic Lensing Hypothesis with the Prompt Emission of GRB 221009A

Authors: Joshua S. Bloom

Abstract: Even at modest amplification, the optical depth to gravitational lensing through the Galaxy is $<10^{-5}$. However, the large apparent isotropic-equivalent energy of GRB 221009A coupled with a path through low Galactic latitude suggests that the conditional probability that this particular GRB was lensed is greater than the very low a priori expectation. With the extreme brightness of the prompt e… ▽ More Even at modest amplification, the optical depth to gravitational lensing through the Galaxy is $<10^{-5}$. However, the large apparent isotropic-equivalent energy of GRB 221009A coupled with a path through low Galactic latitude suggests that the conditional probability that this particular GRB was lensed is greater than the very low a priori expectation. With the extreme brightness of the prompt emission, this Galactic lensing hypothesis can be constrained by autocorrelation analysis of Fermi photons on 0.1-1000 ms timescales. In relating lensing mass, magnification, and autocorrelation timescale, I show that a lensed-induced autocorrelation signature by stellar lenses falls below the minimal variability timescale (MVT) expected from a black hole central engine. However, lensing by Galactic dark matter MACHOs ($M_l > 10-1000\,M_\odot$) could be confirmed with this approach. Regardless, at a peak $γ$-ray photon rate of $>30$ ms$^{-1}$, GRB 221009A represents a prime opportunity to measure the smallest MVTs of GRBs. △ Less

Submitted 15 March, 2024; originally announced March 2024.

Comments: 4 pages, 1 figure. Published October 2022 in Research Notes of the AAS, Volume 6, Number 10

Journal ref: Res. Notes AAS 6 220 (2022)

Iterative assembly of $^{171}$Yb atom arrays with cavity-enhanced optical lattices

Authors: M. A. Norcia, H. Kim, W. B. Cairncross, M. Stone, A. Ryou, M. Jaffe, M. O. Brown, K. Barnes, P. Battaglino, T. C. Bohdanowicz, A. Brown, K. Cassella, C. -A. Chen, R. Coxe, D. Crow, J. Epstein, C. Griger, E. Halperin, F. Hummel, A. M. W. Jones, J. M. Kindem, J. King, K. Kotru, J. Lauigan, M. Li , et al. (25 additional authors not shown)

Abstract: Assembling and maintaining large arrays of individually addressable atoms is a key requirement for continued scaling of neutral-atom-based quantum computers and simulators. In this work, we demonstrate a new paradigm for assembly of atomic arrays, based on a synergistic combination of optical tweezers and cavity-enhanced optical lattices, and the incremental filling of a target array from a repeti… ▽ More Assembling and maintaining large arrays of individually addressable atoms is a key requirement for continued scaling of neutral-atom-based quantum computers and simulators. In this work, we demonstrate a new paradigm for assembly of atomic arrays, based on a synergistic combination of optical tweezers and cavity-enhanced optical lattices, and the incremental filling of a target array from a repetitively filled reservoir. In this protocol, the tweezers provide microscopic rearrangement of atoms, while the cavity-enhanced lattices enable the creation of large numbers of optical traps with sufficient depth for rapid low-loss imaging of atoms. We apply this protocol to demonstrate near-deterministic filling (99% per-site occupancy) of 1225-site arrays of optical traps. Because the reservoir is repeatedly filled with fresh atoms, the array can be maintained in a filled state indefinitely. We anticipate that this protocol will be compatible with mid-circuit reloading of atoms into a quantum processor, which will be a key capability for running large-scale error-corrected quantum computations whose durations exceed the lifetime of a single atom in the system. △ Less

Submitted 18 June, 2024; v1 submitted 29 January, 2024; originally announced January 2024.

Comments: 8 pages, 6 figures

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)

Stellar Spectra Fitting with Amortized Neural Posterior Estimation and nbi

Authors: Keming Zhang, Tharindu Jayasinghe, Joshua S. Bloom

Abstract: Modern surveys often deliver hundreds of thousands of stellar spectra at once, which are fit to spectral models to derive stellar parameters/labels. Therefore, the technique of Amortized Neural Posterior Estimation (ANPE) stands out as a suitable approach, which enables the inference of large number of targets as sub-linear/constant computational costs. Leveraging our new nbi software package, we… ▽ More Modern surveys often deliver hundreds of thousands of stellar spectra at once, which are fit to spectral models to derive stellar parameters/labels. Therefore, the technique of Amortized Neural Posterior Estimation (ANPE) stands out as a suitable approach, which enables the inference of large number of targets as sub-linear/constant computational costs. Leveraging our new nbi software package, we train an ANPE model for the APOGEE survey and demonstrate its efficacy on both mock and real APOGEE stellar spectra. Unique to the nbi package is its out-of-the-box functionality on astronomical inverse problems with sequential data. As such, we have been able to acquire the trained model with minimal effort. We introduce an effective approach to handling the measurement noise properties inherent in spectral data, which utilizes the actual uncertainties in the observed data. This allows training data to resemble observed data, an aspect that is crucial for ANPE applications. Given the association of spectral data properties with the observing instrument, we discuss the utility of an ANPE "model zoo," where models are trained for specific instruments and distributed under the nbi framework to facilitate real-time stellar parameter inference. △ Less

Submitted 9 December, 2023; originally announced December 2023.

Comments: Accepted to ICML 2023 Workshop on Machine Learning for Astrophysics. 7 pages 3 figures

nbi: the Astronomer's Package for Neural Posterior Estimation

Authors: Keming Zhang, Joshua S. Bloom, Stéfan van der Walt, Nina Hernitschek

Abstract: Despite the promise of Neural Posterior Estimation (NPE) methods in astronomy, the adaptation of NPE into the routine inference workflow has been slow. We identify three critical issues: the need for custom featurizer networks tailored to the observed data, the inference inexactness, and the under-specification of physical forward models. To address the first two issues, we introduce a new framewo… ▽ More Despite the promise of Neural Posterior Estimation (NPE) methods in astronomy, the adaptation of NPE into the routine inference workflow has been slow. We identify three critical issues: the need for custom featurizer networks tailored to the observed data, the inference inexactness, and the under-specification of physical forward models. To address the first two issues, we introduce a new framework and open-source software nbi (Neural Bayesian Inference), which supports both amortized and sequential NPE. First, nbi provides built-in "featurizer" networks with demonstrated efficacy on sequential data, such as light curve and spectra, thus obviating the need for this customization on the user end. Second, we introduce a modified algorithm SNPE-IS, which facilities asymptotically exact inference by using the surrogate posterior under NPE only as a proposal distribution for importance sampling. These features allow nbi to be applied off-the-shelf to astronomical inference problems involving light curves and spectra. We discuss how nbi may serve as an effective alternative to existing methods such as Nested Sampling. Our package is at https://github.com/kmzzhang/nbi. △ Less

Submitted 21 December, 2023; v1 submitted 6 December, 2023; originally announced December 2023.

Comments: Update references. Accepted to NeurIPS 2023 Workshop on Deep Learning and Inverse Problems. Initially appeared at ICML 2023 Workshop on Machine Learning for Astrophysics. Code at https://github.com/kmzzhang/nbi

The ZTF Source Classification Project: III. A Catalog of Variable Sources

Authors: Brian F. Healy, Michael W. Coughlin, Ashish A. Mahabal, Theophile Jegou du Laz, Andrew Drake, Matthew J. Graham, Lynne A. Hillenbrand, Jan van Roestel, Paula Szkody, LeighAnna Zielske, Mohammed Guiga, Muhammad Yusuf Hassan, Jill L. Hughes, Guy Nir, Saagar Parikh, Sungmin Park, Palak Purohit, Umaa Rebbapragada, Draco Reed, Daniel Warshofsky, Avery Wold, Joshua S. Bloom, Frank J. Masci, Reed Riddle, Roger Smith

Abstract: The classification of variable objects provides insight into a wide variety of astrophysics ranging from stellar interiors to galactic nuclei. The Zwicky Transient Facility (ZTF) provides time series observations that record the variability of more than a billion sources. The scale of these data necessitates automated approaches to make a thorough analysis. Building on previous work, this paper re… ▽ More The classification of variable objects provides insight into a wide variety of astrophysics ranging from stellar interiors to galactic nuclei. The Zwicky Transient Facility (ZTF) provides time series observations that record the variability of more than a billion sources. The scale of these data necessitates automated approaches to make a thorough analysis. Building on previous work, this paper reports the results of the ZTF Source Classification Project (SCoPe), which trains neural network and XGBoost machine learning (ML) algorithms to perform dichotomous classification of variable ZTF sources using a manually constructed training set containing 170,632 light curves. We find that several classifiers achieve high precision and recall scores, suggesting the reliability of their predictions for 209,991,147 light curves across 77 ZTF fields. We also identify the most important features for XGB classification and compare the performance of the two ML algorithms, finding a pattern of higher precision among XGB classifiers. The resulting classification catalog is available to the public, and the software developed for SCoPe is open-source and adaptable to future time-domain surveys. △ Less

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

Comments: 23 pages, 12 figures, accepted to ApJS. Catalog available at https://zenodo.org/doi/10.5281/zenodo.8410825

Detectability of Self-Lensing Flares of White Dwarfs with Compact Companions

Authors: Guy Nir, Joshua S. Bloom

Abstract: Binaries containing compact objects, if viewed close to edge on, can produce periodic brightening events under certain conditions on the masses, radii, and binary separation. Such flares are caused by one object gravitational lensing the other, in what is known as self-lensing flares. We present a simulation tool that efficiently reproduces the main features of self-lensing flares and facilitates… ▽ More Binaries containing compact objects, if viewed close to edge on, can produce periodic brightening events under certain conditions on the masses, radii, and binary separation. Such flares are caused by one object gravitational lensing the other, in what is known as self-lensing flares. We present a simulation tool that efficiently reproduces the main features of self-lensing flares and facilitates a detection sensitivity analysis for various sky surveys. We estimate the detection prospects for a handful of representative surveys when searching for systems of either two white dwarfs, or a white dwarf with other compact objects, i.e., neutron stars and black holes. We find only a marginal ability to detect such systems in existing surveys. However, we estimate many such systems could be detectable by surveys in the near future, including the Vera Rubin observatory. We provide a quantitative analysis of the detectability of double-compact object self-lensing flares across the landscape of system parameters, and a qualitative discussion of survey and followup approaches to distinguish such flares from confounding events, such as stellar flares, satellite glints, and cosmic rays. We estimate 0.3, 3 and 247 double white dwarf systems could be detected by TESS, ZTF, and LSST, respectively. A similar number of systems with a neutron star or black hole companion could be detected, but we caution that the number densities of such binaries is model dependent and so our detection estimates. Such binaries can be used to constrain models of the end states of binary evolution. △ Less

Submitted 24 November, 2023; originally announced November 2023.

Directional Gaussian spatial processes for South African wind data

Authors: Jacobus S. Blom, Priyanka Nagar, Andriette Bekker

Abstract: Accurate wind pattern modelling is crucial for various applications, including renewable energy, agriculture, and climate adaptation. In this paper, we introduce the wrapped Gaussian spatial process (WGSP), as well as the projected Gaussian spatial process (PGSP) custom-tailored for South Africa's intricate wind behaviour. Unlike conventional models struggling with the circular nature of wind dire… ▽ More Accurate wind pattern modelling is crucial for various applications, including renewable energy, agriculture, and climate adaptation. In this paper, we introduce the wrapped Gaussian spatial process (WGSP), as well as the projected Gaussian spatial process (PGSP) custom-tailored for South Africa's intricate wind behaviour. Unlike conventional models struggling with the circular nature of wind direction, the WGSP and PGSP adeptly incorporate circular statistics to address this challenge. Leveraging historical data sourced from meteorological stations throughout South Africa, the WGSP and PGSP significantly increase predictive accuracy while capturing the nuanced spatial dependencies inherent to wind patterns. The superiority of the PGSP model in capturing the structural characteristics of the South African wind data is evident. As opposed to the PGSP, the WGSP model is computationally less demanding, allows for the use of less informative priors, and its parameters are more easily interpretable. The implications of this study are far-reaching, offering potential benefits ranging from the optimisation of renewable energy systems to the informed decision-making in agriculture and climate adaptation strategies. The WGSP and PGSP emerge as robust and invaluable tools, facilitating precise modelling of wind patterns within the dynamic context of South Africa. △ Less

Submitted 10 November, 2023; originally announced November 2023.

Comments: 9 pages, 9 Figures, conference proceedings paper

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

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

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

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

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

Transitive generalized toggle groups containing a cycle

Authors: Jonathan S. Bloom, Dan Saracino

Abstract: In \cite{striker2018rowmotion} Striker generalized Cameron and Fon-Der-Flaass's notion of a toggle group. In this paper we begin the study of transitive generalized toggle groups that contain a cycle. We first show that if such a group has degree $n$ and contains a transposition or a 3-cycle then the group contains $A_n$. Using the result about transpositions, we then prove that a transitive gener… ▽ More In \cite{striker2018rowmotion} Striker generalized Cameron and Fon-Der-Flaass's notion of a toggle group. In this paper we begin the study of transitive generalized toggle groups that contain a cycle. We first show that if such a group has degree $n$ and contains a transposition or a 3-cycle then the group contains $A_n$. Using the result about transpositions, we then prove that a transitive generalized toggle group that contains a short cycle must be primitive. Employing a result of Jones \cite{jones2014primitive}, which relies on the classification of the finite simple groups, we conclude that any transitive generalized toggle group of degree $n$ that contains a cycle with at least 3 fixed points must also contain $A_n$. Finally, we look at imprimitive generalized toggle groups containing a long cycle and show that they decompose into a direct product of primitive generalized toggle groups each containing a long cycle. △ Less

Submitted 17 October, 2023; originally announced October 2023.

Comments: 11 pages

MSC Class: 05E18

Resolving the explosion of supernova 2023ixf in Messier 101 within its complex circumstellar environment

Authors: E. A. Zimmerman, I. Irani, P. Chen, A. Gal-Yam, S. Schulze, D. A. Perley, J. Sollerman, A. V. Filippenko, T. Shenar, O. Yaron, S. Shahaf, R. J. Bruch, E. O. Ofek, A. De Cia, T. G. Brink, Y. Yang, S. S. Vasylyev, S. Ben Ami, M. Aubert, A. Badash, J. S. Bloom, P. J. Brown, K. De, G. Dimitriadis, C. Fransson , et al. (32 additional authors not shown)

Abstract: Observing a supernova explosion shortly after it occurs can reveal important information about the physics of stellar explosions and the nature of the progenitor stars of supernovae (SNe). When a star with a well-defined edge explodes in vacuum, the first photons to escape from its surface appear as a brief shock-breakout flare. The duration of this flare can extend to at most a few hours even for… ▽ More Observing a supernova explosion shortly after it occurs can reveal important information about the physics of stellar explosions and the nature of the progenitor stars of supernovae (SNe). When a star with a well-defined edge explodes in vacuum, the first photons to escape from its surface appear as a brief shock-breakout flare. The duration of this flare can extend to at most a few hours even for nonspherical breakouts from supergiant stars, after which the explosion ejecta should expand and cool. Alternatively, for stars exploding within a distribution of sufficiently dense optically thick circumstellar material, the first photons escape from the material beyond the stellar edge, and the duration of the initial flare can extend to several days, during which the escaping emission indicates photospheric heating. The difficulty in detecting SN explosions promptly after the event has so far limited data regarding supergiant stellar explosions mostly to serendipitous observations that, owing to the lack of ultraviolet (UV) data, were unable to determine whether the early emission is heating or cooling, and hence the nature of the early explosion event. Here, we report observations of SN 2023ixf in the nearby galaxy M101, covering the early days of the event. Using UV spectroscopy from the Hubble Space Telescope (HST) as well as a comprehensive set of additional multiwavelength observations, we trace the photometric and spectroscopic evolution of the event and are able to temporally resolve the emergence and evolution of the SN emission. △ Less

Submitted 27 March, 2024; v1 submitted 16 October, 2023; originally announced October 2023.

Journal ref: Nature 627, 759 (2024)

Machine-Learning Enhanced Photometric Analysis of the Extremely Bright GRB 210822A

Authors: Camila Angulo-Valdez, Rosa L. Becerra, Margarita Pereyra, Keneth Garcia-Cifuentes, Felipe Vargas, Alan M. Watson, Fabio De Colle, Nissim Fraija, Nathaniel R. Butler, Maria G. Dainotti, Simone Dichiara, William H. Lee, Eleonora Troja, Joshua S. Bloom, J. Jesús González, Alexander S. Kutyrev, J. Xavier Prochaska, Enrico Ramirez-Ruiz, Michael G. Richer

Abstract: We present analytical and numerical models of the bright long GRB 210822A at $z=1.736$. The intrinsic extreme brightness exhibited in the optical, which is very similar to other bright GRBs (e.g., GRBs 080319B, 130427A, 160625A 190114C, and 221009A), makes GRB 210822A an ideal case for studying the evolution of this particular kind of GRB. We use optical data from the RATIR instrument starting at… ▽ More We present analytical and numerical models of the bright long GRB 210822A at $z=1.736$. The intrinsic extreme brightness exhibited in the optical, which is very similar to other bright GRBs (e.g., GRBs 080319B, 130427A, 160625A 190114C, and 221009A), makes GRB 210822A an ideal case for studying the evolution of this particular kind of GRB. We use optical data from the RATIR instrument starting at $T+315.9$ s, with publicly available optical data from other ground-based observatories, as well as Swift/UVOT, and X-ray data from the Swift/XRT instrument. The temporal profiles and spectral properties during the late stages align consistently with the conventional forward shock model, complemented by a reverse shock element that dominates optical emissions during the initial phases ($T<300$ s). Furthermore, we observe a break at $T=80000$s that we interpreted as evidence of a jet break, which constrains the opening angle to be about $θ_\mathrm{j}=(3-5)$ degrees. Finally, we apply a machine-learning technique to model the multi-wavelength light curve of GRB 210822A using the AFTERGLOWPY library. We estimate the angle of sight $θ_{obs}=(6.4 \pm 0.1) \times 10^{-1}$ degrees, the energy $E_0=(7.9 \pm 1.6)\times 10^{53}$ ergs, the electron index $p=2.54 \pm 0.10$, the thermal energy fraction in electrons $ε_\mathrm{e}=(4.63 \pm 0.91) \times 10^{-5}$ and in the magnetic field $ε_\mathrm{B}= (8.66 \pm 1.01) \times 10^{-6}$, the efficiency $χ= 0.89 \pm 0.01$, and the density of the surrounding medium $n_\mathrm{0} = 0.85 \pm 0.01 cm^{-3}$. △ Less

Submitted 17 November, 2023; v1 submitted 18 September, 2023; originally announced September 2023.

Comments: Resubmitted to MNRAS after moderate revision, 12 pages, 6 figures

The Progenitor Star of SN 2023ixf: A Massive Red Supergiant with Enhanced, Episodic Pre-Supernova Mass Loss

Authors: Yu-Jing Qin, Keming Zhang, Joshua Bloom, Jesper Sollerman, Erez A. Zimmerman, Ido Irani, Steve Schulze, Avishay Gal-Yam, Mansi Kasliwal, Michael W. Coughlin, Daniel A. Perley, Christoffer Fremling, Shrinivas Kulkarni

Abstract: We identify the progenitor star of SN 2023ixf in the nearby galaxy Messier 101 using Keck/NIRC2 adaptive optics imaging and pre-explosion HST/ACS images. The supernova position, localized with diffraction-spike pattern and high precision relative astrometry, unambiguously coincides with a single progenitor candidate of m_F814W=24.96(-0.04)(+0.05). Forced photometry further recovers 2-sigma detecti… ▽ More We identify the progenitor star of SN 2023ixf in the nearby galaxy Messier 101 using Keck/NIRC2 adaptive optics imaging and pre-explosion HST/ACS images. The supernova position, localized with diffraction-spike pattern and high precision relative astrometry, unambiguously coincides with a single progenitor candidate of m_F814W=24.96(-0.04)(+0.05). Forced photometry further recovers 2-sigma detections in the F673N and F675W bands and imposes robust flux limits on the bluer bands. Given the reported infrared excess and semi-regular variability of the progenitor, we fit a time-dependent spectral energy distribution (SED) model of a dusty red supergiant (RSG) to a combined dataset of HST photometry, as well as ground-based near-infrared and Spitzer/IRAC [3.6], [4.5] photometry from the literature. The progenitor closely resembles a RSG of T_eff=3343+/-27 K and logL=5.10+/-0.02, with a 0.11+/-0.01 dex (25.2+/-1.7 per cent) variation over the mean luminosity at a period of P=1128.3+/-6.5 days, heavily obscured by a dust envelope with an optical depth of tau=2.83+/-0.03 at 1 micron (or A_V=10.28+/-0.11 mag). Such observed signatures match a post-main sequence star of 18.1(-1.2)(+0.7) Msun, close to the most massive SN II progenitor, with a pulsation-enhanced mass-loss rate of M_dot=(3.58+/-0.15) x 10^(-4) Msun/yr. The dense and confined circumstellar material is likely ejected during the last episode of radial pulsation before the explosion. Notably, we find strong evidence for periodic variation of tau (or both T_eff and tau) along with luminosity, a necessary assumption to reproduce the wavelength dependence of the variability, which implies dust sublimation and condensation during radial pulsations. Given the observed SED, partial dust obscuration remains a possible scenario, but any unobstructed binary companion over 7.1 Msun can be ruled out. △ Less

Submitted 18 September, 2023; originally announced September 2023.

Comments: submitted

Keck Infrared Transient Survey I: Survey Description and Data Release 1

Authors: S. Tinyanont, R. J. Foley, K. Taggart, K. W. Davis, N. LeBaron, J. E. Andrews, M. J. Bustamante-Rosell, Y. Camacho-Neves, R. Chornock, D. A. Coulter, L. Galbany, S. W. Jha, C. D. Kilpatrick, L. A. Kwok, C. Larison, J. R. Pierel, M. R. Siebert, G. Aldering, K. Auchettl, J. S. Bloom, S. Dhawan, A. V. Filippenko, K. D. French, A. Gagliano, M. Grayling , et al. (13 additional authors not shown)

Abstract: We present the Keck Infrared Transient Survey (KITS), a NASA Key Strategic Mission Support program to obtain near-infrared (NIR) spectra of astrophysical transients of all types, and its first data release, consisting of 105 NIR spectra of 50 transients. Such a data set is essential as we enter a new era of IR astronomy with the James Webb Space Telescope (JWST) and the upcoming Nancy Grace Roman… ▽ More We present the Keck Infrared Transient Survey (KITS), a NASA Key Strategic Mission Support program to obtain near-infrared (NIR) spectra of astrophysical transients of all types, and its first data release, consisting of 105 NIR spectra of 50 transients. Such a data set is essential as we enter a new era of IR astronomy with the James Webb Space Telescope (JWST) and the upcoming Nancy Grace Roman Space Telescope (Roman). NIR spectral templates will be essential to search JWST images for stellar explosions of the first stars and to plan an effective Roma} SN Ia cosmology survey, both key science objectives for mission success. Between 2022 February and 2023 July, we systematically obtained 274 NIR spectra of 146 astronomical transients, representing a significant increase in the number of available NIR spectra in the literature. The first data release includes data from the 2022A semester. We systematically observed three samples: a flux-limited sample that includes all transients $<$17 mag in a red optical band (usually ZTF r or ATLAS o bands); a volume-limited sample including all transients within redshift $z < 0.01$ ($D \approx 50$ Mpc); and an SN Ia sample targeting objects at phases and light-curve parameters that had scant existing NIR data in the literature. The flux-limited sample is 39% complete (60% excluding SNe Ia), while the volume-limited sample is 54% complete and is 79% complete to $z = 0.005$. All completeness numbers will rise with the inclusion of data from other telescopes in future data releases. Transient classes observed include common Type Ia and core-collapse supernovae, tidal disruption events (TDEs), luminous red novae, and the newly categorized hydrogen-free/helium-poor interacting Type Icn supernovae. We describe our observing procedures and data reduction using Pypeit, which requires minimal human interaction to ensure reproducibility. △ Less

Submitted 13 September, 2023; originally announced September 2023.

Decentralized Multi-Agent Reinforcement Learning with Global State Prediction

Authors: Joshua Bloom, Pranjal Paliwal, Apratim Mukherjee, Carlo Pinciroli

Abstract: Deep reinforcement learning (DRL) has seen remarkable success in the control of single robots. However, applying DRL to robot swarms presents significant challenges. A critical challenge is non-stationarity, which occurs when two or more robots update individual or shared policies concurrently, thereby engaging in an interdependent training process with no guarantees of convergence. Circumventing… ▽ More Deep reinforcement learning (DRL) has seen remarkable success in the control of single robots. However, applying DRL to robot swarms presents significant challenges. A critical challenge is non-stationarity, which occurs when two or more robots update individual or shared policies concurrently, thereby engaging in an interdependent training process with no guarantees of convergence. Circumventing non-stationarity typically involves training the robots with global information about other agents' states and/or actions. In contrast, in this paper we explore how to remove the need for global information. We pose our problem as a Partially Observable Markov Decision Process, due to the absence of global knowledge on other agents. Using collective transport as a testbed scenario, we study two approaches to multi-agent training. In the first, the robots exchange no messages, and are trained to rely on implicit communication through push-and-pull on the object to transport. In the second approach, we introduce Global State Prediction (GSP), a network trained to forma a belief over the swarm as a whole and predict its future states. We provide a comprehensive study over four well-known deep reinforcement learning algorithms in environments with obstacles, measuring performance as the successful transport of the object to the goal within a desired time-frame. Through an ablation study, we show that including GSP boosts performance and increases robustness when compared with methods that use global knowledge. △ Less

Submitted 28 August, 2023; v1 submitted 22 June, 2023; originally announced June 2023.

Mid-circuit qubit measurement and rearrangement in a $^{171}$Yb atomic array

Authors: M. A. Norcia, W. B. Cairncross, K. Barnes, P. Battaglino, A. Brown, M. O. Brown, K. Cassella, C. -A. Chen, R. Coxe, D. Crow, J. Epstein, C. Griger, A. M. W. Jones, H. Kim, J. M. Kindem, J. King, S. S. Kondov, K. Kotru, J. Lauigan, M. Li, M. Lu, E. Megidish, J. Marjanovic, M. McDonald, T. Mittiga , et al. (20 additional authors not shown)

Abstract: Measurement-based quantum error correction relies on the ability to determine the state of a subset of qubits (ancillae) within a processor without revealing or disturbing the state of the remaining qubits. Among neutral-atom based platforms, a scalable, high-fidelity approach to mid-circuit measurement that retains the ancilla qubits in a state suitable for future operations has not yet been demo… ▽ More Measurement-based quantum error correction relies on the ability to determine the state of a subset of qubits (ancillae) within a processor without revealing or disturbing the state of the remaining qubits. Among neutral-atom based platforms, a scalable, high-fidelity approach to mid-circuit measurement that retains the ancilla qubits in a state suitable for future operations has not yet been demonstrated. In this work, we perform imaging using a narrow-linewidth transition in an array of tweezer-confined $^{171}$Yb atoms to demonstrate nondestructive state-selective and site-selective detection. By applying site-specific light shifts, selected atoms within the array can be hidden from imaging light, which allows a subset of qubits to be measured while causing only percent-level errors on the remaining qubits. As a proof-of-principle demonstration of conditional operations based on the results of the mid-circuit measurements, and of our ability to reuse ancilla qubits, we perform conditional refilling of ancilla sites to correct for occasional atom loss, while maintaining the coherence of data qubits. Looking towards true continuous operation, we demonstrate loading of a magneto-optical trap with a minimal degree of qubit decoherence. △ Less

Submitted 2 October, 2023; v1 submitted 30 May, 2023; originally announced May 2023.

Comments: 12 pages, 7 figures

A data science platform to enable time-domain astronomy

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

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

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

Comments: Accepted to ApJS

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

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

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

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

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

Collapsars as Sites of r-process Nucleosynthesis: Systematic Near-Infrared Follow-up of Type Ic-BL Supernovae

Authors: Shreya Anand, Jennifer Barnes, Sheng Yang, Mansi M. Kasliwal, Michael W. Coughlin, Jesper Sollerman, Kishalay De, Christoffer Fremling, Alessandra Corsi, Anna Y. Q. Ho, Arvind Balasubramanian, Conor Omand, Gokul P. Srinivasaragavan, S. Bradley Cenko, Tomas Ahumada, Igor Andreoni, Aishwarya Dahiwale, Kaustav Kashyap Das, Jacob Jencson, Viraj Karambelkar, Harsh Kumar, Brian D. Metzger, Daniel Perley, Nikhil Sarin, Tassilo Schweyer , et al. (19 additional authors not shown)

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

Submitted 12 February, 2024; v1 submitted 17 February, 2023; originally announced February 2023.

Comments: 37 pages, 8 figures, 4 tables

Journal ref: ApJ, 962, 68 (2024)

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

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

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

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

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

Comments: Matches published version in Nature Astronomy

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

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

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

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

Comments: 62 pages, 10 figures, 4 tables

Journal ref: Nature 617 (2023), 477

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

The Tidal Disruption Event AT2021ehb: Evidence of Relativistic Disk Reflection, and Rapid Evolution of the Disk-Corona System

Authors: Yuhan Yao, Wenbin Lu, Muryel Guolo, Dheeraj R. Pasham, Suvi Gezari, Marat Gilfanov, Keith C. Gendreau, Fiona Harrison, S. Bradley Cenko, S. R. Kulkarni, Jon M. Miller, Dominic J. Walton, Javier A. García, Sjoert van Velzen, Kate D. Alexander, James C. A. Miller-Jones, Matt Nicholl, Erica Hammerstein, Pavel Medvedev, Daniel Stern, Vikram Ravi, R. Sunyaev, Joshua S. Bloom, Matthew J. Graham, Erik C. Kool , et al. (7 additional authors not shown)

Abstract: We present X-ray, UV, optical, and radio observations of the nearby ($\approx78$ Mpc) tidal disruption event (TDE) AT2021ehb/ZTF21aanxhjv during its first 430 days of evolution. AT2021ehb occurs in the nucleus of a galaxy hosting a $\approx 10^{7}\,M_\odot$ black hole ($M_{\rm BH}$ inferred from host galaxy scaling relations). High-cadence Swift and NICER monitoring reveals a delayed X-ray brighte… ▽ More We present X-ray, UV, optical, and radio observations of the nearby ($\approx78$ Mpc) tidal disruption event (TDE) AT2021ehb/ZTF21aanxhjv during its first 430 days of evolution. AT2021ehb occurs in the nucleus of a galaxy hosting a $\approx 10^{7}\,M_\odot$ black hole ($M_{\rm BH}$ inferred from host galaxy scaling relations). High-cadence Swift and NICER monitoring reveals a delayed X-ray brightening. The spectrum first undergoes a gradual ${\rm soft }\rightarrow{\rm hard}$ transition and then suddenly turns soft again within 3 days at $δt\approx 272$ days during which the X-ray flux drops by a factor of ten. In the joint NICER+NuSTAR observation ($δt =264$ days, harder state), we observe a prominent non-thermal component up to 30 keV and an extremely broad emission line in the iron K band. The bolometric luminosity of AT2021ehb reaches a maximum of $6.0^{+10.4}_{-3.8}\% L_{\rm Edd}$ when the X-ray spectrum is the hardest. During the dramatic X-ray evolution, no radio emission is detected, the UV/optical luminosity stays relatively constant, and the optical spectra are featureless. We propose the following interpretations: (i) the ${\rm soft }\rightarrow{\rm hard}$ transition may be caused by the gradual formation of a magnetically dominated corona; (ii) hard X-ray photons escape from the system along solid angles with low scattering optical depth ($\sim\,$a few) whereas the UV/optical emission is likely generated by reprocessing materials with much larger column density -- the system is highly aspherical; (iii) the abrupt X-ray flux drop may be triggered by the thermal-viscous instability in the inner accretion flow leading to a much thinner disk. △ Less

Submitted 24 August, 2022; v1 submitted 25 June, 2022; originally announced June 2022.

Comments: 35 pages, 21 figures, accepted by ApJ

SN2020qlb: A hydrogen-poor superluminous supernova with well-characterized light curve undulations

Authors: S. L. West, R. Lunnan, C. M. B. Omand, T. Kangas, S. Schulze, N. Strotjohann, S. Yang, C. Fransson, J. Sollerman, D. Perley, L. Yan, T. -W. Chen, Z. H. Chen, K. Taggart, C. Fremling, J. S. Bloom, A. Drake, M. J. Graham, M. M. Kasliwal, R. Laher, M. S. Medford, J. D. Neill, R. Riddle, D. Shupe

Abstract: SN\,2020qlb (ZTF20abobpcb) is a hydrogen-poor superluminous supernova (SLSN-I) that is among the most luminous (maximum M$_{g} = -22.25$ mag) and that has one of the longest rise times (77 days from explosion to maximum). We estimate the total radiated energy to be $>2.1\times10^{51}$ erg. SN\,2020qlb has a well-sampled light curve that exhibits clear near and post peak undulations, a phenomenon s… ▽ More SN\,2020qlb (ZTF20abobpcb) is a hydrogen-poor superluminous supernova (SLSN-I) that is among the most luminous (maximum M$_{g} = -22.25$ mag) and that has one of the longest rise times (77 days from explosion to maximum). We estimate the total radiated energy to be $>2.1\times10^{51}$ erg. SN\,2020qlb has a well-sampled light curve that exhibits clear near and post peak undulations, a phenomenon seen in other SLSNe, whose physical origin is still unknown. We discuss the potential power source of this immense explosion as well as the mechanisms behind its observed light curve undulations. We analyze photospheric spectra and compare them to other SLSNe-I. We constructed the bolometric light curve using photometry from a large data set of observations from the Zwicky Transient Facility (ZTF), Liverpool Telescope (LT), and Neil Gehrels Swift Observatory and compare it with radioactive, circumstellar interaction and magnetar models. Model residuals and light curve polynomial fit residuals are analyzed to estimate the undulation timescale and amplitude. We also determine host galaxy properties based on imaging and spectroscopy data, including a detection of the [O III]$λ$4363, auroral line, allowing for a direct metallicity measurement. We rule out the Arnett $^{56}$Ni decay model for SN\,2020qlb's light curve due to unphysical parameter results. Our most favored power source is the magnetic dipole spin-down energy deposition of a magnetar. Two to three near peak oscillations, intriguingly similar to those of SN\,2015bn, were found in the magnetar model residuals with a timescale of $32\pm6$ days and an amplitude of 6$\%$ of peak luminosity. We rule out centrally located undulation sources due to timescale considerations; and we favor the result of ejecta interactions with circumstellar material (CSM) density fluctuations as the source of the undulations. △ Less

Submitted 7 December, 2022; v1 submitted 23 May, 2022; originally announced May 2022.

Comments: 22 pages, 25 figures, submitted to A&A

Journal ref: A&A 670, A7 (2023)

Strictly increasing and decreasing sequences in subintervals of words and a conjecture of Guo and Poznanović

Authors: Jonathan S. Bloom, Dan Saracino

Abstract: We prove a conjecture of Guo and Poznanović concerning chains in certain 01-fillings of moon polyominoes. A key ingredient of our proof is a correspondence between words $w$ and pairs $(\mathcal{W}(w), \mathcal{M}(w))$ of increasing tableaux such that $\mathcal{M}(w)$ determines the lengths of the longest strictly increasing and strictly decreasing sequences in every subinterval of $w$. We define… ▽ More We prove a conjecture of Guo and Poznanović concerning chains in certain 01-fillings of moon polyominoes. A key ingredient of our proof is a correspondence between words $w$ and pairs $(\mathcal{W}(w), \mathcal{M}(w))$ of increasing tableaux such that $\mathcal{M}(w)$ determines the lengths of the longest strictly increasing and strictly decreasing sequences in every subinterval of $w$. We define this correspondence by using Thomas and Yong's K-infusion operator and then use it to obtain the bijections that prove the conjecture of Guo and Poznanović. In constructing our bijections we introduce new variants of the RSK correspondence and Knuth equivalence. △ Less

Submitted 11 November, 2023; v1 submitted 11 April, 2022; originally announced April 2022.

Comments: 36 pages

MSC Class: 05A05; 05A19

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

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

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

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

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

A Study of Reinforcement Learning Algorithms for Aggregates of Minimalistic Robots

Authors: Joshua Bloom, Apratim Mukherjee, Carlo Pinciroli

Abstract: The aim of this paper is to study how to apply deep reinforcement learning for the control of aggregates of minimalistic robots. We define aggregates as groups of robots with a physical connection that compels them to form a specified shape. In our case, the robots are pre-attached to an object that must be collectively transported to a known location. Minimalism, in our setting, stems from the ba… ▽ More The aim of this paper is to study how to apply deep reinforcement learning for the control of aggregates of minimalistic robots. We define aggregates as groups of robots with a physical connection that compels them to form a specified shape. In our case, the robots are pre-attached to an object that must be collectively transported to a known location. Minimalism, in our setting, stems from the barebone capabilities we assume: The robots can sense the target location and the immediate obstacles, but lack the means to communicate explicitly through, e.g., message-passing. In our setting, communication is implicit, i.e., mediated by aggregated push-and-pull on the object exerted by each robot. We analyze the ability to reach coordinated behavior of four well-known algorithms for deep reinforcement learning (DQN, DDQN, DDPG, and TD3). Our experiments include robot failures and different types of environmental obstacles. We compare the performance of the best control strategies found, highlighting strengths and weaknesses of each of the considered training algorithms. △ Less

Submitted 28 March, 2022; originally announced March 2022.

Comments: 7 pages, submitted to IROS 2022

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.

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

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

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

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

Comments: 30 pages, 20 figures, Accepted by APJ

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

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

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

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

Comments: 38 pages, 25 figures, Accepted by APJ

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

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

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

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

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

HEALPix Alchemy: Fast All-Sky Geometry and Image Arithmetic in a Relational Database for Multimessenger Astronomy Brokers

Authors: Leo P. Singer, B. Parazin, Michael W. Coughlin, Joshua S. Bloom, Arien Crellin-Quick, Daniel A. Goldstein, Stéfan van der Walt

Abstract: Efficient searches for electromagnetic counterparts to gravitational wave, high-energy neutrino, and gamma-ray burst events demand rapid processing of image arithmetic and geometry set operations in a database to cross-match galaxy catalogs, observation footprints, and all-sky images. Here we introduce HEALPix Alchemy, an open-source, pure Python implementation of a set of methods that enables rap… ▽ More Efficient searches for electromagnetic counterparts to gravitational wave, high-energy neutrino, and gamma-ray burst events demand rapid processing of image arithmetic and geometry set operations in a database to cross-match galaxy catalogs, observation footprints, and all-sky images. Here we introduce HEALPix Alchemy, an open-source, pure Python implementation of a set of methods that enables rapid all-sky geometry calculations. HEALPix Alchemy is built upon HEALPix, a spatial indexing strategy that is widely used in astronomical databases as well as the native format of LIGO-Virgo-KAGRA gravitational-wave sky localization maps. Our approach leverages new multirange types built into the PostgreSQL 14 database engine. This enables fast all-sky queries against probabilistic multimessenger event localizations and telescope survey footprints. Questions such as "What are the galaxies contained within the 90% credible region of an event?" and "What is the rank-ordered list of the fields within an observing footprint with the highest probability of containing the event?" can be performed in less than a few seconds on commodity hardware using off-the-shelf cloud-managed database implementations without server-side database extensions. Common queries scale roughly linearly with the number of telescope pointings. As the number of fields grows into the hundreds or thousands, HEALPix Alchemy is orders of magnitude faster than other implementations. HEALPix Alchemy is now used as the spatial geometry engine within SkyPortal, which forms the basis of the Zwicky Transient Facility transient marshal, called Fritz. △ Less

Submitted 27 April, 2022; v1 submitted 13 December, 2021; originally announced December 2021.

Comments: For source code, see https://github.com/skyportal/healpix-alchemy

Report number: LIGO-P2100433-v6

Journal ref: AJ 163 209 (2022)

Science with the Ultraviolet Explorer (UVEX)

Authors: S. R. Kulkarni, Fiona A. Harrison, Brian W. Grefenstette, Hannah P. Earnshaw, Igor Andreoni, Danielle A. Berg, Joshua S. Bloom, S. Bradley Cenko, Ryan Chornock, Jessie L. Christiansen, Michael W. Coughlin, Alexander Wuollet Criswell, Behnam Darvish, Kaustav K. Das, Kishalay De, Luc Dessart, Don Dixon, Bas Dorsman, Kareem El-Badry, Christopher Evans, K. E. Saavik Ford, Christoffer Fremling, Boris T. Gansicke, Suvi Gezari, Y. Goetberg , et al. (31 additional authors not shown)

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

Submitted 17 January, 2023; v1 submitted 30 November, 2021; originally announced November 2021.

Comments: 69 pages, 43 figures

GRB 191016A: The onset of the forward shock and evidence of late energy injection

Authors: M. Pereyra, N. Fraija, A. M. Watson, R. L. Becerra, N. R. Butler, F. De Colle, E. Troja, S. Dichiara, E. Fraire-Bonilla, W. H. Lee, A. S. Kutyrev, J. X. Prochaska, J. S. Bloom, J. J. González, E. Ramirez-Ruiz, M. G. Richer

Abstract: We present optical and near-infrared photometric observations of GRB 191016 with the COATLI, DDOTI and RATIR ground-based telescopes over the first three nights. We present the temporal evolution of the optical afterglow and describe 5 different stages that were not completely characterized in previous works, mainly due to scarcity of data points to accurately fit the different components of the o… ▽ More We present optical and near-infrared photometric observations of GRB 191016 with the COATLI, DDOTI and RATIR ground-based telescopes over the first three nights. We present the temporal evolution of the optical afterglow and describe 5 different stages that were not completely characterized in previous works, mainly due to scarcity of data points to accurately fit the different components of the optical emission. After the end of the prompt gamma-ray emission, we observed the afterglow rise slowly in the optical and near-infrared (NIR) wavelengths and peak at around T+1450s in all filters. This was followed by an early decay, a clear plateau from T+5000s to T+11000s, and then a regular late decay. We also present evidence of the jet break at later times, with a temporal index in good agreement with the temporal slope obtained from X-ray observations. Although many of the features observed in the optical light curves of GRBs are usually well explained by a reverse shock (RS) or forward shock(FS), the shallowness of the optical rise and enhanced peak emission in the GRB191016A afterglow is not well-fitted by only a FS or a RS. We propose a theoretical model which considers both of these components and combines an evolving FS with a later embedded RS and a subsequent late energy injection from the central engine activity. We use this model to successfully explain the temporal evolution of the light curves and discuss its implications on the fireball properties. △ Less

Submitted 30 November, 2021; originally announced November 2021.

Comments: 15 pages, 5 figures, 7 Tables, Submitted to MNRAS

A Ubiquitous Unifying Degeneracy in Two-Body Microlensing Systems

Authors: Keming Zhang, B. Scott Gaudi, Joshua S. Bloom

Abstract: While gravitational microlensing by planetary systems provides unique vistas on the properties of exoplanets, observations of a given 2-body microlensing event can often be interpreted with multiple distinct physical configurations. Such ambiguities are typically attributed to the close-wide and inner-outer types of degeneracies that arise from transformation invariances and symmetries of microlen… ▽ More While gravitational microlensing by planetary systems provides unique vistas on the properties of exoplanets, observations of a given 2-body microlensing event can often be interpreted with multiple distinct physical configurations. Such ambiguities are typically attributed to the close-wide and inner-outer types of degeneracies that arise from transformation invariances and symmetries of microlensing caustics. However, there remain unexplained inconsistencies between aforementioned theories and observations. Here, leveraging a fast machine learning inference framework, we present the discovery of the offset degeneracy, which concerns a magnification-matching behaviour on the lens-axis and is formulated independent of caustics. This offset degeneracy unifies the close-wide and inner-outer degeneracies, generalises to resonant topologies, and upon reanalysis, not only appears ubiquitous in previously published planetary events with 2-fold degenerate solutions, but also resolves prior inconsistencies. Our analysis demonstrates that degenerate caustics do not strictly result in degenerate magnifications and that the commonly invoked close-wide degeneracy essentially never arises in actual events. Moreover, it is shown that parameters in offset degenerate configurations are related by a simple expression. This suggests the existence of a deeper symmetry in the equations governing 2-body lenses than previously recognised. △ Less

Submitted 10 May, 2022; v1 submitted 26 November, 2021; originally announced November 2021.

Comments: 21 pages, 8 figures, submitted

Journal ref: Nat Astron 6, 782-787 (2022)

Target of Opportunity Observations of Gravitational Wave Events with Vera C. Rubin Observatory

Authors: Igor Andreoni, Raffaella Margutti, Om Sharan Salafia, B. Parazin, V. Ashley Villar, Michael W. Coughlin, Peter Yoachim, Kris Mortensen, Daniel Brethauer, S. J. Smartt, Mansi M. Kasliwal, Kate D. Alexander, Shreya Anand, E. Berger, Maria Grazia Bernardini, Federica B. Bianco, Peter K. Blanchard, Joshua S. Bloom, Enzo Brocato, Mattia Bulla, Regis Cartier, S. Bradley Cenko, Ryan Chornock, Christopher M. Copperwheat, Alessandra Corsi , et al. (30 additional authors not shown)

Abstract: The discovery of the electromagnetic counterpart to the binary neutron star merger GW170817 has opened the era of gravitational-wave multi-messenger astronomy. Rapid identification of the optical/infrared kilonova enabled a precise localization of the source, which paved the way to deep multi-wavelength follow-up and its myriad of related science results. Fully exploiting this new territory of exp… ▽ More The discovery of the electromagnetic counterpart to the binary neutron star merger GW170817 has opened the era of gravitational-wave multi-messenger astronomy. Rapid identification of the optical/infrared kilonova enabled a precise localization of the source, which paved the way to deep multi-wavelength follow-up and its myriad of related science results. Fully exploiting this new territory of exploration requires the acquisition of electromagnetic data from samples of neutron star mergers and other gravitational wave sources. After GW170817, the frontier is now to map the diversity of kilonova properties and provide more stringent constraints on the Hubble constant, and enable new tests of fundamental physics. The Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST) can play a key role in this field in the 2020s, when an improved network of gravitational-wave detectors is expected to reach a sensitivity that will enable the discovery of a high rate of merger events involving neutron stars (about tens per year) out to distances of several hundred Mpc. We design comprehensive target-of-opportunity observing strategies for follow-up of gravitational-wave triggers that will make the Rubin Observatory the premier instrument for discovery and early characterization of neutron star and other compact object mergers, and yet unknown classes of gravitational wave events. △ Less

Submitted 20 April, 2022; v1 submitted 2 November, 2021; originally announced November 2021.

Comments: Accepted for publication in ApJS. arXiv admin note: text overlap with arXiv:1812.04051

Assembly and coherent control of a register of nuclear spin qubits

Authors: Katrina Barnes, Peter Battaglino, Benjamin J. Bloom, Kayleigh Cassella, Robin Coxe, Nicole Crisosto, Jonathan P. King, Stanimir S. Kondov, Krish Kotru, Stuart C. Larsen, Joseph Lauigan, Brian J. Lester, Mickey McDonald, Eli Megidish, Sandeep Narayanaswami, Ciro Nishiguchi, Remy Notermans, Lucas S. Peng, Albert Ryou, Tsung-Yao Wu, Michael Yarwood

Abstract: We introduce an optical tweezer platform for assembling and individually manipulating a two-dimensional register of nuclear spin qubits. Each nuclear spin qubit is encoded in the ground $^{1}S_{0}$ manifold of $^{87}$Sr and is individually manipulated by site-selective addressing beams. We observe that spin relaxation is negligible after 5 seconds, indicating that $T_1\gg5$ s. Furthermore, utilizi… ▽ More We introduce an optical tweezer platform for assembling and individually manipulating a two-dimensional register of nuclear spin qubits. Each nuclear spin qubit is encoded in the ground $^{1}S_{0}$ manifold of $^{87}$Sr and is individually manipulated by site-selective addressing beams. We observe that spin relaxation is negligible after 5 seconds, indicating that $T_1\gg5$ s. Furthermore, utilizing simultaneous manipulation of subsets of qubits, we demonstrate significant phase coherence over the entire register, estimating $T_2^\star = \left(21\pm7\right)$ s and measuring $T_2^\text{echo}=\left(42\pm6\right)$ s. △ Less

Submitted 10 August, 2021; originally announced August 2021.

Comments: 10 pages, 4 figures

The Gravity Collective: A Search for the Electromagnetic Counterpart to the Neutron Star-Black Hole Merger GW190814

Authors: Charles D. Kilpatrick, David A. Coulter, Iair Arcavi, Thomas G. Brink, Georgios Dimitriadis, Alexei V. Filippenko, Ryan J. Foley, D. Andrew Howell, David O. Jones, Martin Makler, Anthony L. Piro, César Rojas-Bravo, David J. Sand, Jonathan J. Swift, Douglas Tucker, WeiKang Zheng, Sahar S. Allam, James T. Annis, Juanita Antilen, Tristan G. Bachmann, Joshua S. Bloom, Clecio R. Bom, K. Azalee Bostroem, Dillon Brout, Jamison Burke , et al. (57 additional authors not shown)

Abstract: We present optical follow-up imaging obtained with the Katzman Automatic Imaging Telescope, Las Cumbres Observatory Global Telescope Network, Nickel Telescope, Swope Telescope, and Thacher Telescope of the LIGO/Virgo gravitational wave (GW) signal from the neutron star-black hole (NSBH) merger GW190814. We searched the GW190814 localization region (19 deg$^{2}$ for the 90th percentile best localiz… ▽ More We present optical follow-up imaging obtained with the Katzman Automatic Imaging Telescope, Las Cumbres Observatory Global Telescope Network, Nickel Telescope, Swope Telescope, and Thacher Telescope of the LIGO/Virgo gravitational wave (GW) signal from the neutron star-black hole (NSBH) merger GW190814. We searched the GW190814 localization region (19 deg$^{2}$ for the 90th percentile best localization), covering a total of 51 deg$^{2}$ and 94.6% of the two-dimensional localization region. Analyzing the properties of 189 transients that we consider as candidate counterparts to the NSBH merger, including their localizations, discovery times from merger, optical spectra, likely host-galaxy redshifts, and photometric evolution, we conclude that none of these objects are likely to be associated with GW190814. Based on this finding, we consider the likely optical properties of an electromagnetic counterpart to GW190814, including possible kilonovae and short gamma-ray burst afterglows. Using the joint limits from our follow-up imaging, we conclude that a counterpart with an $r$-band decline rate of 0.68 mag day$^{-1}$, similar to the kilonova AT 2017gfo, could peak at an absolute magnitude of at most $-17.8$ mag (50% confidence). Our data are not constraining for ''red'' kilonovae and rule out ''blue'' kilonovae with $M>0.5 M_{\odot}$ (30% confidence). We strongly rule out all known types of short gamma-ray burst afterglows with viewing angles $<$17$^{\circ}$ assuming an initial jet opening angle of $\sim$$5.2^{\circ}$ and explosion energies and circumburst densities similar to afterglows explored in the literature. Finally, we explore the possibility that GW190814 merged in the disk of an active galactic nucleus, of which we find four in the localization region, but we do not find any candidate counterparts among these sources. △ Less

Submitted 12 June, 2021; originally announced June 2021.

Comments: 86 pages, 9 figures

Real-Time Likelihood-Free Inference of Roman Binary Microlensing Events with Amortized Neural Posterior Estimation

Authors: Keming Zhang, Joshua S. Bloom, B. Scott Gaudi, Francois Lanusse, Casey Lam, Jessica R. Lu

Abstract: Fast and automated inference of binary-lens, single-source (2L1S) microlensing events with sampling-based Bayesian algorithms (e.g., Markov Chain Monte Carlo; MCMC) is challenged on two fronts: high computational cost of likelihood evaluations with microlensing simulation codes, and a pathological parameter space where the negative-log-likelihood surface can contain a multitude of local minima tha… ▽ More Fast and automated inference of binary-lens, single-source (2L1S) microlensing events with sampling-based Bayesian algorithms (e.g., Markov Chain Monte Carlo; MCMC) is challenged on two fronts: high computational cost of likelihood evaluations with microlensing simulation codes, and a pathological parameter space where the negative-log-likelihood surface can contain a multitude of local minima that are narrow and deep. Analysis of 2L1S events usually involves grid searches over some parameters to locate approximate solutions as a prerequisite to posterior sampling, an expensive process that often requires human-in-the-loop domain expertise. As the next-generation, space-based microlensing survey with the Roman Space Telescope is expected to yield thousands of binary microlensing events, a new fast and automated method is desirable. Here, we present a likelihood-free inference (LFI) approach named amortized neural posterior estimation, where a neural density estimator (NDE) learns a surrogate posterior $\hat{p}(θ|x)$ as an observation-parametrized conditional probability distribution, from pre-computed simulations over the full prior space. Trained on 291,012 simulated Roman-like 2L1S simulations, the NDE produces accurate and precise posteriors within seconds for any observation within the prior support without requiring a domain expert in the loop, thus allowing for real-time and automated inference. We show that the NDE also captures expected posterior degeneracies. The NDE posterior could then be refined into the exact posterior with a downstream MCMC sampler with minimal burn-in steps. △ Less

Submitted 30 March, 2021; v1 submitted 10 February, 2021; originally announced February 2021.

Comments: 15 pages, 8 figures, 3 tables. Accepted for publication in the Astronomical Journal. This article supersedes arXiv:2010.04156

Estimating Effectiveness of Identifying Human Trafficking via Data Envelopment Analysis

Authors: Geri L. Dimas, Malak El Khalkhali, Alex Bender, Kayse Lee Maass, Renata Konrad, Jeffrey S. Blom, Joe Zhu, Andrew C. Trapp

Abstract: Transit monitoring is a preventative approach used to identify possible cases of human trafficking prior to exploitation while an individual is in transit or before one crosses a border. Transit monitoring is often conducted by non-governmental organizations (NGOs) who train staff to identify and intercept suspicious activity. Love Justice International (LJI) is a well-established NGO that has bee… ▽ More Transit monitoring is a preventative approach used to identify possible cases of human trafficking prior to exploitation while an individual is in transit or before one crosses a border. Transit monitoring is often conducted by non-governmental organizations (NGOs) who train staff to identify and intercept suspicious activity. Love Justice International (LJI) is a well-established NGO that has been conducting transit monitoring for years along the Nepal-India border at multiple monitoring stations. In partnership with LJI, we developed a system that uses data envelopment analysis (DEA) to help LJI decision-makers evaluate the performance of these stations at intercepting potential human-trafficking victims given the amount of resources (e.g. staff, etc.) available and make specific operational improvement recommendations. Our model consists of 91 decision-making units (DMUs) from 7 stations over 13 quarters and considers three inputs, four outputs, and 3 homogeneity criteria. Using this model we identified efficient stations, compared rankings of station performance, and recommended strategies to improve efficiency. To the best of our knowledge, this is the first application of DEA in the anti-human trafficking domain. △ Less

Submitted 22 June, 2022; v1 submitted 9 December, 2020; originally announced December 2020.

Informing the Cataclysmic Variable Donor Sequence from Gaia DR2 Color-Magnitude and Inferred Variability Metrics

Authors: Ellianna S. Abrahams, Joshua S. Bloom, Nami Mowlavi, Paula Szkody, Hans-Walter Rix, Jean-Paul Ventura, Thomas G. Brink, Alexei V. Filippenko

Abstract: Short-period cataclysmic variables (spCVs), with orbital periods below the period gap ($P_{orb}$ < 2 hr), offer insight into the evolutionary models of CVs and can serve as strong emitters of gravitational waves (GWs). To identify new spCV candidates, we crossmatch a catalog of known CVs to sources with robust parallaxes in the Gaia second data release (DR2). We uncover and fit an apparently monot… ▽ More Short-period cataclysmic variables (spCVs), with orbital periods below the period gap ($P_{orb}$ < 2 hr), offer insight into the evolutionary models of CVs and can serve as strong emitters of gravitational waves (GWs). To identify new spCV candidates, we crossmatch a catalog of known CVs to sources with robust parallaxes in the Gaia second data release (DR2). We uncover and fit an apparently monotonic relationship between the color--absolute-magnitude diagram (CMD) position and $P_{orb}$ of these CVs, revealed in DR2. To supplement this relation, we develop a method for identifying sources with large photometric variability, a characteristic trait of spCVs. Using all available Gaia light curves, we construct a machine-learned regression model to predict variability metrics for sources in the CMD locus of known spCVs based solely on time-averaged covariates present in DR2. Using this approach we identify 3,253 candidate spCVs, of which $\sim$95% are previously unknown. Inspection of archival SDSS spectra of these candidates suggests that $>$82% are likely to be spCVs: a noticeably higher recovery rate than previous light curve searches, which bias toward active systems. We obtain optical spectra of 9 new systems at Lick Observatory and confirm that all objects are CV systems. We measure $P_{orb}$ for 7 systems using archival Gaia and Palomar Transient Factory light curves, 3 of which do not have previous $P_{orb}$ measurements. We use the CMD-$P_{orb}$ relation to infer the detectability of these systems to the upcoming LISA mission, and find that six sources may be coherent LISA verification binaries, with an estimated SNR > 5 in the 4 yr mission. This paper demonstrates that the time-averaged Gaia catalog is a powerful tool in the methodical discovery and characterization of time-varying objects, making it complementary to missions like ZTF, TESS, and the Vera Rubin LSST. △ Less

Submitted 24 November, 2020; originally announced November 2020.

Comments: 35 pages, 24 figures. Submitted to AAS Journals

Classification of Periodic Variable Stars with Novel Cyclic-Permutation Invariant Neural Networks

Authors: Keming Zhang, Joshua S. Bloom

Abstract: Neural networks (NNs) have been shown to be competitive against state-of-the-art feature engineering and random forest (RF) classification of periodic variable stars. Although previous work utilising NNs has made use of periodicity by period folding multiple-cycle time-series into a single cycle -- from time-space to phase-space -- no approach to date has taken advantage of the fact that network p… ▽ More Neural networks (NNs) have been shown to be competitive against state-of-the-art feature engineering and random forest (RF) classification of periodic variable stars. Although previous work utilising NNs has made use of periodicity by period folding multiple-cycle time-series into a single cycle -- from time-space to phase-space -- no approach to date has taken advantage of the fact that network predictions should be invariant to the initial phase of the period-folded sequence. Initial phase is exogenous to the physical origin of the variability and should thus be factored out. Here, we present cyclic-permutation invariant networks, a novel class of NNs for which invariance to phase shifts is guaranteed through polar coordinate convolutions, which we implement by means of "Symmetry Padding." Across three different datasets of variable star light curves, we show that two implementations of the cyclic-permutation invariant network: the iTCN and the iResNet, consistently outperform non-invariant baselines and reduce overall error rates by between 4% to 22%. Over a 10-class OGLE-III sample, the iTCN/iResNet achieves an average per-class accuracy of 93.4%/93.3%, compared to RNN/RF accuracies of 70.5%/89.5% in a recent study using the same data. Finding improvement on a non-astronomy benchmark, we suggest that the methodology introduced here should also be applicable to a wide range of science domains where periodic data abounds due to physical symmetries. △ Less

Submitted 28 April, 2021; v1 submitted 2 November, 2020; originally announced November 2020.

Comments: Accepted to MNRAS. Extended abstract accepted to NeurIPS 2020: Machine Learning for the Physical Sciences. Initially appeared in ICLR 2020: Fundamental Science in the Era of AI. Open source implementation at https://github.com/kmzzhang/periodicnetwork. 8 pages, 4 figures

Automating Inference of Binary Microlensing Events with Neural Density Estimation

Authors: Keming Zhang, Joshua S. Bloom, B. Scott Gaudi, Francois Lanusse, Casey Lam, Jessica Lu

Abstract: Automated inference of binary microlensing events with traditional sampling-based algorithms such as MCMC has been hampered by the slowness of the physical forward model and the pathological likelihood surface. Current analysis of such events requires both expert knowledge and large-scale grid searches to locate the approximate solution as a prerequisite to MCMC posterior sampling. As the next gen… ▽ More Automated inference of binary microlensing events with traditional sampling-based algorithms such as MCMC has been hampered by the slowness of the physical forward model and the pathological likelihood surface. Current analysis of such events requires both expert knowledge and large-scale grid searches to locate the approximate solution as a prerequisite to MCMC posterior sampling. As the next generation, space-based microlensing survey with the Roman Space Observatory is expected to yield thousands of binary microlensing events, a new scalable and automated approach is desired. Here, we present an automated inference method based on neural density estimation (NDE). We show that the NDE trained on simulated Roman data not only produces fast, accurate, and precise posteriors but also captures expected posterior degeneracies. A hybrid NDE-MCMC framework can further be applied to produce the exact posterior. △ Less

Submitted 15 February, 2021; v1 submitted 8 October, 2020; originally announced October 2020.

Comments: 7 pages, 1 figure. This article is superseded by arXiv:2102.05673. Accepted to the ML4PS workshop at NeurIPS 2020

Modelling the prompt optical emission of GRB 180325A: the evolution of a spike from the optical to gamma-rays

Authors: Rosa L. Becerra, Fabio De Colle, Jorge Cantó, Susana Lizano, Ricardo F. González, Jonathan Granot, Alain Klotz, Alan M. Watson, Nissim Fraija, Anabella T. Araudo, Eleonora Troja, Jean Luc Atteia, William H. Lee, Damien Turpin, Joshua S. Bloom, Michael Boer, Nathaniel R. Butler, José J. González, Alexander S. Kutyrev, J. Xavier Prochaska, Enrico Ramírez-Ruíz, Michael G. Richer, Carlos G. Román Zúñiga

Abstract: The transition from prompt to the afterglow emission is one of the most exciting and least understood phases in gamma-ray bursts (GRBs). Correlations among optical, X-ray and gamma-ray emission in GRBs have been explored, to attempt to answer whether the earliest optical emission comes from internal and/or external shocks. We present optical photometric observations of GRB 180325A collected with t… ▽ More The transition from prompt to the afterglow emission is one of the most exciting and least understood phases in gamma-ray bursts (GRBs). Correlations among optical, X-ray and gamma-ray emission in GRBs have been explored, to attempt to answer whether the earliest optical emission comes from internal and/or external shocks. We present optical photometric observations of GRB 180325A collected with the TAROT and RATIR ground-based telescopes. These observations show two strong optical flashes with separate peaks at $\sim50\;$s and $\sim120\;$s, followed by a temporally extended optical emission. We also present X-rays and gamma-ray observations of GRB 180325A, detected by the Burst Alert Telescope (BAT) and X-ray Telescope (XRT), on the Neil Gehrels Swift observatory, which both observed a narrow flash at $\sim80\;$s. We show that the prompt gamma- and X-ray early emission shares similar temporal and spectral features consistent with internal dissipation within the relativistic outflow (e.g. by internal shocks or magnetic reconnection), while the early optical flashes are likely generated by the reverse shock that decelerates the ejecta as it sweeps up the external medium. △ Less

Submitted 3 December, 2020; v1 submitted 28 September, 2020; originally announced September 2020.

Comments: Accepted for publication in ApJ

The Palomar Transient Factory Core-Collapse Supernova Host-Galaxy Sample. I. Host-Galaxy Distribution Functions and Environment-Dependence of CCSNe

Authors: Steve Schulze, Ofer Yaron, Jesper Sollerman, Giorgos Leloudas, Amit Gal, Angus H. Wright, Ragnhild Lunnan, Avishay Gal-Yam, Eran O. Ofek, Daniel A. Perley, Alexei V. Filippenko, Mansi M. Kasliwal, Shri R. Kulkarni, Peter E. Nugent, Robert M. Quimby, Mark Sullivan, Nora Linn Strothjohann, Iair Arcavi, Sagi Ben-Ami, Federica Bianco, Joshua S. Bloom, Kishalay De, Morgan Fraser, Christoffer U. Fremling, Assaf Horesh , et al. (29 additional authors not shown)

Abstract: Several thousand core-collapse supernovae (CCSNe) of different flavors have been discovered so far. However, identifying their progenitors has remained an outstanding open question in astrophysics. Studies of SN host galaxies have proven to be powerful in providing constraints on the progenitor populations. In this paper, we present all CCSNe detected between 2009 and 2017 by the Palomar Transient… ▽ More Several thousand core-collapse supernovae (CCSNe) of different flavors have been discovered so far. However, identifying their progenitors has remained an outstanding open question in astrophysics. Studies of SN host galaxies have proven to be powerful in providing constraints on the progenitor populations. In this paper, we present all CCSNe detected between 2009 and 2017 by the Palomar Transient Factory. This sample includes 888 SNe of 12 distinct classes out to redshift $z\approx1$. We present the photometric properties of their host galaxies from the far-ultraviolet to the mid-infrared and model the host-galaxy spectral energy distributions to derive physical properties. The galaxy mass functions of Type Ic, Ib, IIb, II, and IIn SNe ranges from $10^{5}$ to $10^{11.5}~M_\odot$, probing the entire mass range of star-forming galaxies down to the least-massive star-forming galaxies known. Moreover, the galaxy mass distributions are consistent with models of star-formation-weighted mass functions. Regular CCSNe are hence direct tracers of star formation. Small but notable differences exist between some of the SN classes. Type Ib/c SNe prefer galaxies with slightly higher masses (i.e., higher metallicities) and star-formation rates than Type IIb and II SNe. These differences are less pronounced than previously thought. H-poor SLSNe and SNe~Ic-BL are scarce in galaxies above $10^{10}~M_\odot$. Their progenitors require environments with metallicities of $<0.4$ and $<1$ solar, respectively. In addition, the hosts of H-poor SLSNe are dominated by a younger stellar population than all other classes of CCSNe. Our findings corroborate the notion that low-metallicity \textit{and} young age play an important role in the formation of SLSN progenitors. △ Less

Submitted 13 August, 2020; originally announced August 2020.

Comments: 24 pages main text, 14 figures, 9 Tables, catalogue available at http://www.github.com/steveschulze/PTF