Showing 1–50 of 261 results for author: Margutti, R

Constraints on Relativistic Jets from the Fast X-ray Transient 210423 using Prompt Radio Follow-Up Observations

Authors: Dina Ibrahimzade, R. Margutti, J. S. Bright, P. Blanchard, K. Paterson, D. Lin, H. Sears, A. Polzin, I. Andreoni, G. Schroeder, K. D. Alexander, E. Berger, D. L. Coppejans, A. Hajela, J. Irwin, T. Laskar, B. D. Metzger, J. C. Rastinejad, L. Rhodes

Abstract: Fast X-ray Transients (FXTs) are a new observational class of phenomena with no clear physical origin. This is at least partially a consequence of limited multi-wavelength follow up of this class of transients in real time. Here we present deep optical ($g-$ and $i-$ band) photometry with Keck, and prompt radio observations with the VLA of FXT 210423 obtained at ${δt \approx 14-36}$ days since the… ▽ More Fast X-ray Transients (FXTs) are a new observational class of phenomena with no clear physical origin. This is at least partially a consequence of limited multi-wavelength follow up of this class of transients in real time. Here we present deep optical ($g-$ and $i-$ band) photometry with Keck, and prompt radio observations with the VLA of FXT 210423 obtained at ${δt \approx 14-36}$ days since the X-ray trigger. We use these multi-band observations, combined with publicly available data sets, to constrain the presence and physical properties of on-axis and off-axis relativistic jets such as those that can be launched by neutron-star mergers and tidal disruption events, which are among the proposed theoretical scenarios of FXTs. Considering a wide range of possible redshifts $z\le3.5$, circumstellar medium (CSM) density $n={10^{-6}-10^{-1}\,\rm{cm^{-3}}}$, isotropic-equivalent jet kinetic energy $E_{k,iso}={10^{48}-10^{55}\,\rm{erg}}$, we find that we can rule out wide jets with opening angle ${θ_{j}=15^{\circ}}$ viewed within ${10^{\circ}}$ off-axis. For more collimated jets (${θ_{j}=3^{\circ}}$) we can only rule out on-axis (${θ_{obs}=0^{\circ}}$) orientations. This study highlights the constraining power of prompt multi-wavelength observations of FXTs discovered in real time by current (e.g., Einstein Probe) and future facilities. △ Less

Submitted 11 July, 2024; v1 submitted 9 July, 2024; originally announced July 2024.

Comments: 14 pages, 6 figures

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

SN 2024ggi in NGC 3621: Rising Ionization in a Nearby, CSM-Interacting Type II Supernova

Authors: W. V. Jacobson-Galán, K. W. Davis, C. D. Kilpatrick, L. Dessart, R. Margutti, R. Chornock, R. J. Foley, P. Arunachalam, K. Auchettl, C. R. Bom, R. Cartier, D. A. Coulter, G. Dimitriadis, D. Dickinson, M. R. Drout, A. T. Gagliano, C. Gall, B. Garretson, L. Izzo, D. O. Jones, N. LeBaron, H. -Y. Miao, D. Milisavljevic, Y. -C. Pan, A. Rest , et al. (6 additional authors not shown)

Abstract: We present UV/optical/NIR observations and modeling of supernova (SN) 2024ggi, a type II supernova (SN II) located in NGC 3621 at 7.2 Mpc. Early-time ("flash") spectroscopy of SN 2024ggi within +0.8 days of discovery shows emission lines of H I, He I, C III, and N III with a narrow core and broad, symmetric wings (i.e., IIn-like) arising from the photoionized, optically-thick, unshocked circumstel… ▽ More We present UV/optical/NIR observations and modeling of supernova (SN) 2024ggi, a type II supernova (SN II) located in NGC 3621 at 7.2 Mpc. Early-time ("flash") spectroscopy of SN 2024ggi within +0.8 days of discovery shows emission lines of H I, He I, C III, and N III with a narrow core and broad, symmetric wings (i.e., IIn-like) arising from the photoionized, optically-thick, unshocked circumstellar material (CSM) that surrounded the progenitor star at shock breakout. By the next spectral epoch at +1.5 days, SN 2024ggi showed a rise in ionization as emission lines of He II, C IV, N IV/V and O V became visible. This phenomenon is temporally consistent with a blueward shift in the UV/optical colors, both likely the result of shock breakout in an extended, dense CSM. The IIn-like features in SN 2024ggi persist on a timescale of $t_{\rm IIn} = 3.8 \pm 1.6$ days at which time a reduction in CSM density allows the detection of Doppler broadened features from the fastest SN material. SN 2024ggi has peak UV/optical absolute magnitudes of $M_{\rm w2} = -18.7$ mag and $M_{\rm g} = -18.1$ mag that are consistent with the known population of CSM-interacting SNe II. Comparison of SN 2024ggi with a grid of radiation hydrodynamics and non-local thermodynamic equilibrium (nLTE) radiative-transfer simulations suggests a progenitor mass-loss rate of $\dot{M} = 10^{-2}$M$_{\odot}$ yr$^{-1}$ ($v_w$ = 50 km/s), confined to a distance of $r < 5\times 10^{14}$ cm. Assuming a wind velocity of $v_w$ = 50 km/s, the progenitor star underwent an enhanced mass-loss episode in the last ~3 years before explosion. △ Less

Submitted 25 June, 2024; v1 submitted 29 April, 2024; originally announced April 2024.

Comments: Accepted for publication in ApJ. arXiv admin note: text overlap with arXiv:2306.04721, arXiv:2403.02382

The Peculiar Radio Evolution of the Tidal Disruption Event ASASSN-19bt

Authors: Collin T. Christy, Kate D. Alexander, Yvette Cendes, Ryan Chornock, Tanmoy Laskar, Raffaella Margutti, Edo Berger, Michael Bietenholz, Deanne Coppejans, Fabio De Colle, Tarraneh Eftekhari, Thomas W. -S. Holoien, Tatsuya Matsumoto, James C. A. Miller-Jones, Enrico Ramirez-Ruiz, Richard Saxton, Sjoert van Velzen, Mark Wieringa

Abstract: We present detailed radio observations of the tidal disruption event (TDE) ASASSN-19bt/AT2019ahk, obtained with the Australia Telescope Compact Array (ATCA), the Atacama Large Millimeter/submillimeter Array (ALMA), and the MeerKAT radio telescopes, spanning 40 to 1464 days after the onset of the optical flare. We find that ASASSN-19bt displays unusual radio evolution compared to other TDEs, as the… ▽ More We present detailed radio observations of the tidal disruption event (TDE) ASASSN-19bt/AT2019ahk, obtained with the Australia Telescope Compact Array (ATCA), the Atacama Large Millimeter/submillimeter Array (ALMA), and the MeerKAT radio telescopes, spanning 40 to 1464 days after the onset of the optical flare. We find that ASASSN-19bt displays unusual radio evolution compared to other TDEs, as the peak brightness of its radio emission increases rapidly until 457 days post-optical discovery and then plateaus. Using a generalized approach to standard equipartition techniques, we estimate the energy and corresponding physical parameters for two possible emission geometries: a non-relativistic spherical outflow and a relativistic outflow observed from an arbitrary viewing angle. We find that the non-relativistic solution implies a continuous energy rise in the outflow from $E\sim10^{46}$ erg to $E\sim10^{49}$ erg with $β\approx 0.05$, while the off-axis relativistic jet solution instead suggests $E\approx10^{52}$ erg with $Γ\sim10$ erg at late times in the maximally off-axis case. We find that neither model provides a holistic explanation for the origin and evolution of the radio emission, emphasizing the need for more complex models. ASASSN-19bt joins the population of TDEs that display unusual radio emission at late times. Conducting long-term radio observations of these TDEs, especially during the later phases, will be crucial for understanding how these types of radio emission in TDEs are produced. △ Less

Submitted 18 April, 2024; originally announced April 2024.

Comments: 25 pages. Submitted to ApJ

Late-time X-ray Observations of the Jetted Tidal Disruption Event AT2022cmc: The Relativistic Jet Shuts Off

Authors: T. Eftekhari, A. Tchekhovskoy, K. D. Alexander, E. Berger, R. Chornock, T. Laskar, R. Margutti, Y. Yao, Y. Cendes, S. Gomez, A. Hajela, D. R. Pasham

Abstract: The tidal disruption event (TDE) AT2022cmc represents the fourth known example of a relativistic jet produced by the tidal disruption of a stray star providing a unique probe of the formation and evolution of relativistic jets in otherwise dormant supermassive black holes (SMBHs). Here we present deep, late-time Chandra observations of AT2022cmc extending to $t_{\rm obs} \approx 400$ days after di… ▽ More The tidal disruption event (TDE) AT2022cmc represents the fourth known example of a relativistic jet produced by the tidal disruption of a stray star providing a unique probe of the formation and evolution of relativistic jets in otherwise dormant supermassive black holes (SMBHs). Here we present deep, late-time Chandra observations of AT2022cmc extending to $t_{\rm obs} \approx 400$ days after disruption. Our observations reveal a sudden decrease in the X-ray brightness by a factor of $\gtrsim 14$ over a factor of $\approx 2.3$ in time, and a deviation from the earlier power-law decline with a steepening $α\gtrsim 3.2$ ($F_X \propto t^{-α}$), steeper than expected for a jet break, and pointing to the cessation of jet activity at $t_{\rm obs} \approx 215$ days. Such a transition has been observed in two previous TDEs (Swift J1644+57 and Swift J2058+05). From the X-ray luminosity and the timescale of jet shutoff, we parameterize the mass of the SMBH in terms of unknown jet efficiency and accreted mass fraction parameters. Motivated by the disk-jet connection in AGN, we favor black hole masses $\lesssim 10^5 \ \rm M_{\odot}$ (where the jet and disk luminosities are comparable), and disfavor larger black holes (in which extremely powerful jets are required to outshine their accretion disks). We additionally estimate a total accreted mass of $\approx 0.1 \rm \ M_{\odot}$. Applying the same formalism to Swift J1644+57 and Swift J2058+05, we favor comparable black hole masses for these TDEs of $\lesssim$ a few $\times 10^5 \ \rm M_{\odot}$, and suggest that jetted TDEs may preferentially form from lower mass black holes when compared to non-relativistic events, owing to generally lower jet and higher disk efficiencies at higher black hole masses. △ Less

Submitted 15 April, 2024; originally announced April 2024.

Comments: Submitted to ApJ; 12 pages, 4 figures

Final Moments II: Observational Properties and Physical Modeling of CSM-Interacting Type II Supernovae

Authors: W. V. Jacobson-Galán, L. Dessart, K. W. Davis, C. D. Kilpatrick, R. Margutti, R. J. Foley, R. Chornock, G. Terreran, D. Hiramatsu, M. Newsome, E. Padilla Gonzalez, C. Pellegrino, D. A. Howell, A. V. Filippenko, J. P. Anderson, C. R. Angus, K. Auchettl, K. A. Bostroem, T. G. Brink, R. Cartier, D. A. Coulter, T. de Boer, M. R. Drout, N. Earl, K. Ertini , et al. (30 additional authors not shown)

Abstract: We present ultraviolet/optical/near-infrared observations and modeling of Type II supernovae (SNe II) whose early-time ($δt < 2$ days) spectra show transient, narrow emission lines from shock ionization of confined ($r < 10^{15}$ cm) circumstellar material (CSM). The observed electron-scattering broadened line profiles (i.e., IIn-like) of HI, He I/II, C III/IV, and N III/IV/V from the CSM persist… ▽ More We present ultraviolet/optical/near-infrared observations and modeling of Type II supernovae (SNe II) whose early-time ($δt < 2$ days) spectra show transient, narrow emission lines from shock ionization of confined ($r < 10^{15}$ cm) circumstellar material (CSM). The observed electron-scattering broadened line profiles (i.e., IIn-like) of HI, He I/II, C III/IV, and N III/IV/V from the CSM persist on a characteristic timescale ($t_{\rm IIn}$) that marks a transition to a lower-density CSM and the emergence of Doppler-broadened features from the fast-moving SN ejecta. Our sample, the largest to date, consists of 39 SNe with early-time IIn-like features in addition to 35 "comparison" SNe with no evidence of early-time IIn-like features, all with ultraviolet observations. The total sample consists of 50 unpublished objects with 474 previously unpublished spectra and 50 multiband light curves, collected primarily through the Young Supernova Experiment and Global Supernova Project collaborations. For all sample objects, we find a significant correlation between peak ultraviolet brightness and both $t_{\rm IIn}$ and the rise time, as well as evidence for enhanced peak luminosities in SNe II with IIn-like features. We quantify mass-loss rates and CSM density for the sample through matching of peak multiband absolute magnitudes, rise times, $t_{\rm IIn}$ and optical SN spectra with a grid of radiation hydrodynamics and non-local thermodynamic equilibrium (nLTE) radiative-transfer simulations. For our grid of models, all with the same underlying explosion, there is a trend between the duration of the electron-scattering broadened line profiles and inferred mass-loss rate: $t_{\rm IIn} \approx 3.8[\dot{M}/(0.01 \textrm{M}_{\odot} \textrm{yr}^{-1})]$ days. △ Less

Submitted 4 March, 2024; originally announced March 2024.

Comments: 58 pages, 24 figures, submitted to ApJ. Supplementary figures available on Github (https://github.com/wynnjacobson-galan/Flash_Spectra_Sample). Data release following publication

SN2023ixf in Messier 101: the twilight years of the progenitor as seen by Pan-STARRS

Authors: Conor L. Ransome, V. Ashley Villar, Anna Tartaglia, Sebastian Javier Gonzalez, Wynn V. Jacobson-Galán, Charles D. Kilpatrick, Raffaella Margutti, Ryan J. Foley, Matthew Grayling, Yuan Qi Ni, Ricardo Yarza, Christine Ye, Katie Auchettl, Thomas de Boer, Kenneth C. Chambers, David A. Coulter, Maria R. Drout, Diego Farias, Christa Gall, Hua Gao, Mark E. Huber, Adaeze L. Ibik, David O. Jones, Nandita Khetan, Chien-Cheng Lin , et al. (6 additional authors not shown)

Abstract: The nearby type II supernova, SN2023ixf in M101 exhibits signatures of early-time interaction with circumstellar material in the first week post-explosion. This material may be the consequence of prior mass loss suffered by the progenitor which possibly manifested in the form of a detectable pre-supernova outburst. We present an analysis of the long-baseline pre-explosion photometric data in $g$,… ▽ More The nearby type II supernova, SN2023ixf in M101 exhibits signatures of early-time interaction with circumstellar material in the first week post-explosion. This material may be the consequence of prior mass loss suffered by the progenitor which possibly manifested in the form of a detectable pre-supernova outburst. We present an analysis of the long-baseline pre-explosion photometric data in $g$, $w$, $r$, $i$, $z$ and $y$ filters from Pan-STARRS as part of the Young Supernova Experiment, spanning $\sim$5,000 days. We find no significant detections in the Pan-STARRS pre-explosion light curve. We train a multilayer perceptron neural network to classify pre-supernova outbursts. We find no evidence of eruptive pre-supernova activity to a limiting absolute magnitude of $-7$. The limiting magnitudes from the full set of $gwrizy$ (average absolute magnitude $\approx$-8) data are consistent with previous pre-explosion studies. We use deep photometry from the literature to constrain the progenitor of SN2023ixf, finding that these data are consistent with a dusty red supergiant (RSG) progenitor with luminosity $\log\left(L/L_\odot\right)$$\approx$5.12 and temperature $\approx$3950K, corresponding to a mass of 14-20 M$_\odot$ △ Less

Submitted 7 December, 2023; originally announced December 2023.

Comments: 19 pages, 8 figures, 1 table

Astronomy as a Field: A Guide for Aspiring Astrophysicists

Authors: Ava Polzin, Yasmeen Asali, Sanah Bhimani, Madison Brady, Mandy C. Chen, Lindsay DeMarchi, Michelle Gurevich, Emily Lichko, Emma Louden, Julie Malewicz, Samantha Pagan, Malena Rice, Zili Shen, Emily Simon, Candice Stauffer, J. Luna Zagorac, Katie Auchettl, Katelyn Breivik, Hsiao-Wen Chen, Deanne Coppejans, Sthabile Kolwa, Raffaella Margutti, Priyamvada Natarajan, Erica Nelson, Kim L. Page , et al. (3 additional authors not shown)

Abstract: This book was created as part of the SIRIUS B VERGE program to orient students to astrophysics as a broad field. The 2023-2024 VERGE program and the printing of this book is funded by the Women and Girls in Astronomy Program via the International Astronomical Union's North American Regional Office of Astronomy for Development and the Heising-Simons Foundation; as a result, this document is written… ▽ More This book was created as part of the SIRIUS B VERGE program to orient students to astrophysics as a broad field. The 2023-2024 VERGE program and the printing of this book is funded by the Women and Girls in Astronomy Program via the International Astronomical Union's North American Regional Office of Astronomy for Development and the Heising-Simons Foundation; as a result, this document is written by women in astronomy for girls who are looking to pursue the field. However, given its universal nature, the material covered in this guide is useful for anyone interested in pursuing astrophysics professionally. △ Less

Submitted 26 December, 2023; v1 submitted 7 December, 2023; originally announced December 2023.

Comments: Introductory guide for students interested in pursuing astrophysics; to be submitted to BAAS

The High Energy X-ray Probe (HEX-P): Sensitive broadband X-ray observations of transient phenomena in the 2030s

Authors: Murray Brightman, Raffaella Margutti, Ava Polzin, Amruta Jaodand, Kenta Hotokezaka, Jason A. J. Alford, Gregg Hallinan, Elias Kammoun, Kunal Mooley, Megan Masterson, Lea Marcotulli, Arne Rau, George A. Younes, Daniel Stern, Javier A. García, Kristin Madsen

Abstract: HEX-P will launch at a time when the sky is being routinely scanned for transient gravitational wave, electromagnetic and neutrino phenomena that will require the capabilities of a sensitive, broadband X-ray telescope for follow up studies. These include the merger of compact objects such as neutron stars and black holes, stellar explosions, and the birth of new compact objects. \hexp\ will probe… ▽ More HEX-P will launch at a time when the sky is being routinely scanned for transient gravitational wave, electromagnetic and neutrino phenomena that will require the capabilities of a sensitive, broadband X-ray telescope for follow up studies. These include the merger of compact objects such as neutron stars and black holes, stellar explosions, and the birth of new compact objects. \hexp\ will probe the accretion and ejecta from these transient phenomena through the study of relativistic outflows and reprocessed emission, provide unique capabilities for understanding jet physics, and potentially revealing the nature of the central engine. △ Less

Submitted 8 November, 2023; originally announced November 2023.

Comments: Submitted to Frontiers in Astronomy and Space Sciences

Roaring to softly whispering: Persistent X-ray emission at the location of the Fast Blue Optical Transient AT2018cow $\sim$3.7 yrs after discovery and implications on accretion-powered scenarios

Authors: G. Migliori, R. Margutti, B. D. Metzger, R. Chornock, C. Vignali, D. Brethauer, D. L. Coppejans, T. Maccarone, L. Rivera Sandoval, J. S. Bright, T. Laskar, D. Milisavljevic, E. Berger, J. Nayana

Abstract: We present the first deep X-ray observations of a luminous FBOT AT2018cow, at $\sim3.7\,\rm{yr}$ since discovery, together with the re-analysis of the observation at $δt\sim 220$ d. X-ray emission is significantly detected at a location consistent with AT2018cow. The very soft X-ray spectrum and sustained luminosity are distinct from the spectral and temporal behavior of the LFBOT in the first… ▽ More We present the first deep X-ray observations of a luminous FBOT AT2018cow, at $\sim3.7\,\rm{yr}$ since discovery, together with the re-analysis of the observation at $δt\sim 220$ d. X-ray emission is significantly detected at a location consistent with AT2018cow. The very soft X-ray spectrum and sustained luminosity are distinct from the spectral and temporal behavior of the LFBOT in the first $\sim100$ d, and would possibly signal the emergence of a new emission component, although a robust association with AT2018cow can only be claimed at $δt \sim220$ d, while at $δt \sim1350$ d contamination of the host galaxy cannot be excluded. We interpret these findings in the context of the late-time panchromatic emission from AT2018cow, which includes the detection of persistent, slowly-fading UV emission with $νL_ν\approx 10^{39}\,\rm{erg\,s^{-1}}$. Similar to previous works, (and in analogy with arguments for Ultra-Luminous X-ray sources --ULXs), these late-time observations are consistent with thin-disks around Intermediate Mass Black Holes (IMBHs, with $M_{\bullet}\approx 10^3-10^4\, \rm{M_{\odot}}$) accreting at sub-Eddington rates. However, differently from previous studies, we find that smaller-mass BHs with $M_{\bullet}\approx 10-100\,\rm{M_{\odot}}$ accreting at $\gtrsim$ the Eddington rate cannot be ruled out, and provide a natural explanation for the inferred compact size ($R_{\rm out} \approx 40\,R_{\odot}$) of the accretion disk years after the optical flare. Most importantly, irrespective of the accretor mass, our study lends support to the hypothesis that LFBOTs are accretion-powered phenomena and that, specifically, LFBOTs constitute electromagnetic manifestations of super-Eddington accreting systems that evolve to $\lesssim$ Eddington over a $\approx 100$ days time scale. △ Less

Submitted 6 February, 2024; v1 submitted 27 September, 2023; originally announced September 2023.

Comments: 17 pages, 3 figures, accepted for publication in ApJL

Constraints on the $z\sim5$ Star-Forming Galaxy Luminosity Function From $\textit{Hubble Space Telescope}$ Imaging of an Unbiased and Complete Sample of Long Gamma-ray Burst Host Galaxies

Authors: Huei Sears, Ryan Chornock, Jay Strader, Daniel A. Perley, Peter K. Blanchard, Raffaella Margutti, Nial R. Tanvir

Abstract: We present rest-frame UV \textit{Hubble Space Telescope} imaging of the largest and most complete sample of 23 long duration gamma-ray burst (GRB) host galaxies between redshifts 4 and 6. Of these 23, we present new WFC3/F110W imaging for 19 of the hosts, which we combine with archival WFC3/F110W and WFC3/F140W imaging for the remaining four. We use the photometry of the host galaxies from this sa… ▽ More We present rest-frame UV \textit{Hubble Space Telescope} imaging of the largest and most complete sample of 23 long duration gamma-ray burst (GRB) host galaxies between redshifts 4 and 6. Of these 23, we present new WFC3/F110W imaging for 19 of the hosts, which we combine with archival WFC3/F110W and WFC3/F140W imaging for the remaining four. We use the photometry of the host galaxies from this sample to characterize both the rest-frame UV luminosity function (LF) and the size-luminosity relation of the sample. We find that when assuming the standard Schechter-function parameterization for the UV LF, the GRB host sample is best fit with $α= -1.30^{+0.30}_{-0.25}$ and $M_* = -20.33^{+0.44}_{-0.54}$ mag, which is consistent with results based on $z\sim5$ Lyman-break galaxies. We find that $\sim68\%$ of our size-luminosity measurements fall within or below the same relation for Lyman-break galaxies at $z\sim4$. This study observationally confirms expectations that at $z\sim5$ Lyman-break and GRB host galaxies should trace the same population and demonstrates the utility of GRBs as probes of hidden star-formation in the high-redshift universe. Under the assumption that GRBs unbiasedly trace star formation at this redshift, our non-detection fraction of 7/23 is consistent at the $95\%$-confidence level with $13 - 53\%$ of star formation at redshift $z\sim5$ occurring in galaxies fainter than our detection limit of $M_{1600 A} \sim -18.3$ mag. △ Less

Submitted 27 August, 2023; originally announced August 2023.

Comments: 26 pages, 5 figures

JWST Observations of the Extraordinary GRB 221009A Reveal an Ordinary Supernova Without Signs of $r$-Process Enrichment in a Low-Metallicity Galaxy

Authors: Peter K. Blanchard, V. Ashley Villar, Ryan Chornock, Tanmoy Laskar, Yijia Li, Joel Leja, Justin Pierel, Edo Berger, Raffaella Margutti, Kate D. Alexander, Jennifer Barnes, Yvette Cendes, Tarraneh Eftekhari, Daniel Kasen, Natalie LeBaron, Brian D. Metzger, James Muzerolle Page, Armin Rest, Huei Sears, Daniel M. Siegel, S. Karthik Yadavalli

Abstract: Identifying the astrophysical sites of the $r$-process, one of the primary mechanisms by which heavy elements are formed, is a key goal of modern astrophysics. The discovery of the brightest gamma-ray burst of all time, GRB 221009A, at a relatively nearby redshift, presented the first opportunity to spectroscopically test the idea that $r$-process elements are produced following the collapse of ra… ▽ More Identifying the astrophysical sites of the $r$-process, one of the primary mechanisms by which heavy elements are formed, is a key goal of modern astrophysics. The discovery of the brightest gamma-ray burst of all time, GRB 221009A, at a relatively nearby redshift, presented the first opportunity to spectroscopically test the idea that $r$-process elements are produced following the collapse of rapidly rotating massive stars. Here we present spectroscopic and photometric $\textit{James Webb Space Telescope}$ (JWST) observations of GRB 221009A obtained $+168$ and $+170$ rest-frame days after the initial gamma-ray trigger, and demonstrate they are well-described by a supernova (SN) and power-law afterglow, with no evidence for an additional component from $r$-process emission, and that the SN component strongly resembles the near-infrared spectra of previous SNe, including SN 1998bw. We further find that the SN associated with GRB 221009A is slightly fainter than the expected brightness of SN 1998bw at this phase, concluding that the SN is therefore not an unusual GRB-SN. We infer a nickel mass of $\approx0.09$ M$_{\odot}$, consistent with the lack of an obvious SN detection in the early-time data. We find that the host galaxy of GRB 221009A has a very low metallicity of $\approx0.12$ Z$_{\odot}$ and our resolved host spectrum shows that GRB 221009A occurred in a unique environment in its host characterized by strong H$_2$ emission lines consistent with recent star formation, which may hint at environmental factors being responsible for its extreme energetics. △ Less

Submitted 27 August, 2023; originally announced August 2023.

Comments: 32 pages, 14 Figures, Submitted to Nature Astronomy

Ubiquitous Late Radio Emission from Tidal Disruption Events

Authors: Yvette Cendes, Edo Berger, Kate D. Alexander, Ryan Chornock, Raffaella Margutti, Brian Metzger, Mark H. Wieringa, Michael F. Bietenholz, Aprajita Hajela, Tanmoy Laskar, Michael C. Stroh, Giacomo Terreran

Abstract: We present radio observations of 23 optically-discovered tidal disruption events (TDEs) on timescales of about 500-3200 days post-discovery. We detect 9 new TDEs that did not have detectable radio emission at earlier times, indicating a late-time brightening after several hundred (and up to 2300 days); an additional 6 TDEs exhibit radio emission whose origin is ambiguous or may be attributed to th… ▽ More We present radio observations of 23 optically-discovered tidal disruption events (TDEs) on timescales of about 500-3200 days post-discovery. We detect 9 new TDEs that did not have detectable radio emission at earlier times, indicating a late-time brightening after several hundred (and up to 2300 days); an additional 6 TDEs exhibit radio emission whose origin is ambiguous or may be attributed to the host galaxy or an AGN. We also report new rising components in two TDEs previously detected in the radio (iPTF16fnl and AT2019dsg) at ~1000 days. While the radio emission in some of the detected TDEs peaked on a timescale of ~2-4 years, more than half of the sample still shows rising emission. The range of luminosities for the sample is 10^37-10^39 erg/s, about two orders of magnitude below the radio luminosity of the relativistic TDE Sw1644+57. Our data set indicates that about 40% of all optical TDEs are detected in the radio hundreds to thousands of days after discovery, and that this is probably more common than early radio emission peaking at ~100 days. Using an equipartition analysis, we find evidence for a delayed launch of the radio-emitting outflows, with delay timescales of ~500-2000 days, inferred velocities of ~0.02-0.15c, and kinetic energies of ~10^47-10^49 erg. We rule out off-axis relativistic jets as a viable explanation for this population, and conclude delayed outflows are a more likely explanation, such as from delayed disk formation. Finally, we find comparable densities in the circumnuclear environments of these TDEs as for those with early radio emission, and find the TDEs still rising in luminosity are consistent with free expansion. We conclude that late radio emission marks a fairly ubiquitous but heretofore overlooked phase of TDE evolution. △ Less

Submitted 25 August, 2023; originally announced August 2023.

Comments: 30 pages. Submitted to ApJ

SN 2022oqm: A Bright and Multi-peaked Calcium-rich Transient

Authors: S. Karthik Yadavalli, V. Ashley Villar, Luca Izzo, Yossef Zenati, Ryan J. Foley, J. Craig Wheeler, Charlotte R. Angus, Dominik Bánhidi, Katie Auchettl, Barna Imre Bíró, Attila Bódi, Zsófia Bodola, Thomas de Boer, Kenneth C. Chambers, Ryan Chornock, David A. Coulter, István Csányi, Borbála Cseh, Srujan Dandu, Kyle W. Davis, Connor Braden Dickinson, Diego Farias, Joseph Farah, Christa Gall, Hua Gao , et al. (38 additional authors not shown)

Abstract: We present the photometric and spectroscopic evolution of SN 2022oqm, a nearby multi-peaked hydrogen- and helium-weak calcium-rich transient (CaRT). SN 2022oqm was detected 13.1 kpc from its host galaxy, the face-on spiral galaxy NGC 5875. Extensive spectroscopic coverage reveals an early hot (T >= 40,000 K) continuum and carbon features observed $\sim$1~day after discovery, SN Ic-like photospheri… ▽ More We present the photometric and spectroscopic evolution of SN 2022oqm, a nearby multi-peaked hydrogen- and helium-weak calcium-rich transient (CaRT). SN 2022oqm was detected 13.1 kpc from its host galaxy, the face-on spiral galaxy NGC 5875. Extensive spectroscopic coverage reveals an early hot (T >= 40,000 K) continuum and carbon features observed $\sim$1~day after discovery, SN Ic-like photospheric-phase spectra, and strong forbidden calcium emission starting 38 days after discovery. SN 2022oqm has a relatively high peak luminosity (MB = -17 mag) for (CaRTs), making it an outlier in the population. We determine that three power sources are necessary to explain the light curve (LC), with each corresponding to a distinct peak. The first peak is powered by an expanding blackbody with a power law luminosity, suggesting shock cooling by circumstellar material (CSM). Subsequent LC evolution is powered by a double radioactive decay model, consistent with two sources of photons diffusing through optically thick ejecta. From the LC, we derive an ejecta mass and 56Ni mass of ~0.6 solar masses and ~0.09 solar masses. Spectroscopic modeling suggests 0.6 solar masses of ejecta, and with well-mixed Fe-peak elements throughout. We discuss several physical origins for SN 2022oqm and find either a surprisingly massive white dwarf progenitor or a peculiar stripped envelope model could explain SN 2022oqm. A stripped envelope explosion inside a dense, hydrogen- and helium-poor CSM, akin to SNe Icn, but with a large 56Ni mass and small CSM mass could explain SN 2022oqm. Alternatively, helium detonation on an unexpectedly massive white dwarf could also explain SN 2022oqm. △ Less

Submitted 4 April, 2024; v1 submitted 24 August, 2023; originally announced August 2023.

Comments: 35 pages, 17 figures, 7 tables, Accepted for Publication in ApJ

SN 2022joj: A Potential Double Detonation with a Thin Helium shell

Authors: E. Padilla Gonzalez, D. A. Howell, G. Terreran, C. McCully, M. Newsome, J. Burke, J. Farah, C. Pellegrino, K. A. Bostroem, G. Hosseinzadeh, J. Pearson, D. J. Sand, M. Shrestha, N. Smith, Y. Dong, N. Meza Retamal, S. Valenti, S. Boos, K. J. Shen, D. Townsley, L. Galbany, L. Piscarreta, R. J. Foley, M. J. Bustamante-Rosell, D. A. Coulter , et al. (12 additional authors not shown)

Abstract: We present photometric and spectroscopic data for SN 2022joj, a nearby peculiar Type Ia supernova (SN Ia) with a fast decline rate ($\rm{Δm_{15,B}=1.4}$ mag). SN 2022joj shows exceedingly red colors, with a value of approximately ${B-V \approx 1.1}$ mag during its initial stages, beginning from $11$ days before maximum brightness. As it evolves the flux shifts towards the blue end of the spectrum,… ▽ More We present photometric and spectroscopic data for SN 2022joj, a nearby peculiar Type Ia supernova (SN Ia) with a fast decline rate ($\rm{Δm_{15,B}=1.4}$ mag). SN 2022joj shows exceedingly red colors, with a value of approximately ${B-V \approx 1.1}$ mag during its initial stages, beginning from $11$ days before maximum brightness. As it evolves the flux shifts towards the blue end of the spectrum, approaching ${B-V \approx 0}$ mag around maximum light. Furthermore, at maximum light and beyond, the photometry is consistent with that of typical SNe Ia. This unusual behavior extends to its spectral characteristics, which initially displayed a red spectrum and later evolved to exhibit greater consistency with typical SNe Ia. We consider two potential explanations for this behavior: double detonation from a helium shell on a sub-Chandrasekhar-mass white dwarf and Chandrasekhar-mass models with a shallow distribution of $\rm{^{56}Ni}$. The shallow nickel models could not reproduce the red colors in the early light curves. Spectroscopically, we find strong agreement between SN 2022joj and double-detonation models with white dwarf masses around 1 $\rm{M_{\odot}}$ and thin He-shell between 0.01 and 0.02 $\rm{M_{\odot}}$. Moreover, the early red colors are explained by line-blanketing absorption from iron-peak elements created by the double detonation scenario in similar mass ranges. However, the nebular spectra composition in SN 2022joj deviates from expectations for double detonation, as we observe strong [Fe III] emission instead of [Ca II] lines as anticipated from double detonation models. More detailed modeling, e.g., including viewing angle effects, is required to test if double detonation models can explain the nebular spectra. △ Less

Submitted 11 August, 2023; originally announced August 2023.

Transients stemming from collapsing massive stars: The missing pieces to advance joint observations of photons and high-energy neutrinos

Authors: Ersilia Guarini, Irene Tamborra, Raffaella Margutti, Enrico Ramirez-Ruiz

Abstract: Collapsing massive stars lead to a broad range of astrophysical transients, whose multi-wavelength emission is powered by a variety of processes including radioactive decay, activity of the central engine, and interaction of the outflows with a dense circumstellar medium. These transients are also candidate factories of neutrinos with energy up to hundreds of PeV. We review the energy released by… ▽ More Collapsing massive stars lead to a broad range of astrophysical transients, whose multi-wavelength emission is powered by a variety of processes including radioactive decay, activity of the central engine, and interaction of the outflows with a dense circumstellar medium. These transients are also candidate factories of neutrinos with energy up to hundreds of PeV. We review the energy released by such astrophysical objects across the electromagnetic wavebands as well as neutrinos, in order to outline a strategy to optimize multi-messenger follow-up programs. We find that, while a significant fraction of the explosion energy can be emitted in the infrared-optical-ultraviolet (UVOIR) band, the optical signal alone is not optimal for neutrino searches. Rather, the neutrino emission is strongly correlated with the one in the radio band, if a dense circumstellar medium surrounds the transient, and with X-rays tracking the activity of the central engine. Joint observations of transients in radio, X-rays, and neutrinos will crucially complement those in the UVOIR band, breaking degeneracies in the transient parameter space. Our findings call for heightened surveys in the radio and X-ray bands to warrant multi-messenger detections. △ Less

Submitted 19 September, 2023; v1 submitted 7 August, 2023; originally announced August 2023.

Comments: 26 pages, including 7 figures and 2 appendices. Minor revisions, conclusions unchanged. Matches version accepted for publication in PRD

Early-time Spectropolarimetry of the Aspherical Type II Supernova SN 2023ixf

Authors: Sergiy S. Vasylyev, Yi Yang, Alexei V. Filippenko, Kishore Patra, Thomas G. Brink, Lifan Wang, Ryan Chornock, Rafaella Margutti, Elinor L. Gates, Adam J. Burgasser, Preethi R. Karpoor, Natalie LeBaron, Emma Softich, Christopher A. Theissen, Eli Wiston, WeiKang Zheng

Abstract: We present six epochs of optical spectropolarimetry of the Type II supernova (SN) 2023ixf ranging from $\sim$ 2 to 15 days after the explosion. Polarimetry was obtained with the Kast double spectrograph on the Shane 3 m telescope at Lick Observatory, representing the earliest such observations ever captured for an SN. We observe a high continuum polarization $p_{\text{cont}} \approx 1$ % on days +… ▽ More We present six epochs of optical spectropolarimetry of the Type II supernova (SN) 2023ixf ranging from $\sim$ 2 to 15 days after the explosion. Polarimetry was obtained with the Kast double spectrograph on the Shane 3 m telescope at Lick Observatory, representing the earliest such observations ever captured for an SN. We observe a high continuum polarization $p_{\text{cont}} \approx 1$ % on days +1.4 and +2.5 before dropping to 0.5 % on day +3.5, persisting at that level up to day +14.5. Remarkably, this change coincides temporally with the disappearance of highly ionized "flash" features. The decrease of the continuum polarization is accompanied by a $\sim 70^\circ$ rotation of the polarization position angle ($PA$) as seen across the continuum. The early evolution of the polarization may indicate different geometric configurations of the electron-scattering atmosphere as seen before and after the disappearance of the emission lines associated with highly-ionized species (e.g., He II, C IV, N III), which are likely produced by elevated mass loss shortly prior to the SN explosion. We interpret the rapid change of polarization and $PA$ from days +2.5 to +4.5 as the time when the SN ejecta emerge from the dense asymmetric circumstellar material (CSM). The temporal evolution of the continuum polarization and the $PA$ is consistent with an aspherical SN explosion that exhibits a distinct geometry compared to the CSM. The rapid follow-up spectropolarimetry of SN 2023ixf during the shock ionization phase reveals an exceptionally asymmetric mass-loss process leading up to the explosion. △ Less

Submitted 12 October, 2023; v1 submitted 3 July, 2023; originally announced July 2023.

Comments: Accepted to ApJL

Luminous Radio Emission from the Superluminous Supernova 2017ens at 3.3 years after explosion

Authors: Raffaella Margutti, J. S. Bright, D. J. Matthews, D. L. Coppejans, K. D. Alexander, E. Berger, M. Bietenholz, R. Chornock, L. DeMarchi, M. R. Drout, T. Eftekhari, W. V. Jacobson-Galan, T. Laskar, D. Milisavljevic, K. Murase, M. Nicholl, C. M. B. Omand, M. Stroh, G. Terreran, A. Z. VanderLey

Abstract: We present the results from a multi-year radio campaign of the superluminous supernova (SLSN) 2017ens, which yielded the earliest radio detection of a SLSN to date at the age of $\sim$3.3 years after explosion. SN2017ens was not detected at radio frequencies in the first $\sim$300\,d of evolution but reached $L_ν\approx 10^{28}\,\rm{erg\,s^{-1}\,cm^{-2}}$ at $ν\sim 6$ GHz, $\sim1250$ days post-exp… ▽ More We present the results from a multi-year radio campaign of the superluminous supernova (SLSN) 2017ens, which yielded the earliest radio detection of a SLSN to date at the age of $\sim$3.3 years after explosion. SN2017ens was not detected at radio frequencies in the first $\sim$300\,d of evolution but reached $L_ν\approx 10^{28}\,\rm{erg\,s^{-1}\,cm^{-2}}$ at $ν\sim 6$ GHz, $\sim1250$ days post-explosion. Interpreting the radio observations in the context of synchrotron radiation from the supernova shock interaction with the circumstellar medium (CSM), we infer an effective mass-loss rate of $\approx 10^{-4}\,\rm{M_{\odot}yr^{-1}}$ at $r\sim 10^{17}$ cm from the explosion's site, for a wind speed of $v_w=50-60\,\rm{km\,s^{-1}}$ measured from optical spectra. These findings are consistent with the spectroscopic metamorphosis of SN2017ens from hydrogen-poor to hydrogen-rich $\sim190$ d after explosion reported by Chen et al., 2018. SN2017ens is thus an addition to the sample of hydrogen-poor massive progenitors that explode shortly after having lost their hydrogen envelope. The inferred circumstellar densities, implying a CSM mass up to $\sim0.5\,\rm{M_{\odot}}$, and low velocity of the ejection point at binary interactions (in the form of common envelope evolution and subsequent envelope ejection) playing a role in shaping the evolution of the stellar progenitors of SLSNe in the $\lesssim 500$ yr preceding core collapse. △ Less

Submitted 23 June, 2023; originally announced June 2023.

Comments: 9 pages, 4 figures, 1 table

Millimeter Observations of the Type II SN2023ixf: Constraints on the Proximate Circumstellar Medium

Authors: Edo Berger, Garrett K. Keating, Raffaella Margutti, Keiichi Maeda, Kate D. Alexander, Yvette Cendes, Tarraneh Eftekhari, Mark Gurwell, Daichi Hiramatsu, Anna Y. Q. Ho, Tanmoy Laskar, Ramprasad Rao, Peter K. G. Williams

Abstract: We present 1.3 mm (230 GHz) observations of the recent and nearby Type II supernova, SN2023ixf, obtained with the Submillimeter Array (SMA) at 2.6-18.6 days after explosion. The observations were obtained as part the SMA Large Program POETS (Pursuit of Extragalactic Transients with the SMA). We do not detect any emission at the location of SN2023ixf, with the deepest limits of… ▽ More We present 1.3 mm (230 GHz) observations of the recent and nearby Type II supernova, SN2023ixf, obtained with the Submillimeter Array (SMA) at 2.6-18.6 days after explosion. The observations were obtained as part the SMA Large Program POETS (Pursuit of Extragalactic Transients with the SMA). We do not detect any emission at the location of SN2023ixf, with the deepest limits of $L_ν(230\,{\rm GHz})\lesssim 8.6\times 10^{25}$ erg s$^{-1}$ Hz$^{-1}$ at 2.7 and 7.7 days, and $L_ν(230\,{\rm GHz})\lesssim 3.4\times 10^{25}$ erg s$^{-1}$ Hz$^{-1}$ at 18.6 days. These limits are about a factor of 2 times dimmer than the mm emission from SN2011dh (IIb), about an order of magnitude dimmer compared to SN1993J (IIb) and SN2018ivc (IIL), and about 30 times dimmer than the most luminous non-relativistic SNe in the mm-band (Type IIb/Ib/Ic). Using these limits in the context of analytical models that include synchrotron self-absorption and free-free absorption we place constraints on the proximate circumstellar medium around the progenitor star, to a scale of $\sim 2\times 10^{15}$ cm, excluding the range $\dot{M}\sim {\rm few}\times 10^{-6}-10^{-2}$ M$_\odot$ yr$^{-1}$ (for a wind velocity, $v_w=115$ km s$^{-1}$, and ejecta velocity, $v_{\rm eje}\sim (1-2)\times 10^4$ km s$^{-1}$). These results are consistent with an inference of the mass loss rate based on optical spectroscopy ($\sim 2\times 10^{-2}$ M$_\odot$ yr$^{-1}$ for $v_w=115$ km s$^{-1}$), but are in tension with the inference from hard X-rays ($\sim 7\times 10^{-4}$ M$_\odot$ yr$^{-1}$ for $v_w=115$ km s$^{-1}$). This tension may be alleviated by a non-homogeneous and confined CSM, consistent with results from high-resolution optical spectroscopy. △ Less

Submitted 15 June, 2023; originally announced June 2023.

Comments: Submitted

Early hard X-rays from the nearby core-collapse supernova SN2023ixf

Authors: Brian W. Grefenstette, Murray Brightman, Hannah P. Earnshaw, Fiona A. Harrison, Raffaella Margutti

Abstract: We present NuSTAR observations of the nearby SN 2023ixf in M101 (d=6.9 Mpc) which provide the earliest hard X-ray detection of a non-relativistic stellar explosion to date at $δt\approx$4-d and $δt\approx$11-d. The spectra are well described by a hot thermal bremsstrahlung continuum with $T>25 \rm{keV}$ shining through a thick neutral medium with a neutral hydrogen column that decreases with time… ▽ More We present NuSTAR observations of the nearby SN 2023ixf in M101 (d=6.9 Mpc) which provide the earliest hard X-ray detection of a non-relativistic stellar explosion to date at $δt\approx$4-d and $δt\approx$11-d. The spectra are well described by a hot thermal bremsstrahlung continuum with $T>25 \rm{keV}$ shining through a thick neutral medium with a neutral hydrogen column that decreases with time (initial $N_{\rm{Hint}}=2.6\times 10^{23} \rm{cm^{-2}}$). A prominent neutral Fe K$α$ emission line is clearly detected, similar to other strongly interacting SNe such as SN2020jl. The rapidly decreasing intrinsic absorption with time suggests the presence of a dense but confined circumstellar medium (CSM). The absorbed broadband X-ray luminosity (0.3--79 keV) is $L_{X} \approx 2.5 \times 10^{40}$ erg s$^{-1}$ during both epochs, with the increase in overall X-ray flux related to the decrease in the absorbing column. Interpreting these observations in the context of thermal bremsstrahlung radiation originating from the interaction of the SN shock with a dense medium we infer large particle densities in excess of $n_{\rm{CSM}}\approx 4\times 10^{8} \rm{cm^{-3}}$ at $r<10^{15} \rm{cm}$, corresponding to an enhanced progenitor mass-loss rate of $\dot M \approx 3\times 10^{-4}$ M$_{\odot}$ yr$^{-1}$ for an assumed wind velocity of $v_w=50$ km s$^{-1}$. △ Less

Submitted 7 June, 2023; originally announced June 2023.

Comments: 8 pages, 2 figures. Submitted

SN2023ixf in Messier 101: A Variable Red Supergiant as the Progenitor Candidate to a Type II Supernova

Authors: Charles D. Kilpatrick, Ryan J. Foley, Wynn V. Jacobson-Galán, Anthony L. Piro, Stephen J. Smartt, Maria R. Drout, Alexander Gagliano, Christa Gall, Jens Hjorth, David O. Jones, Kaisey S. Mandel, Raffaella Margutti, Conor L. Ransome, V. Ashley Villar, David A. Coulter, Hua Gao, David Jacob Matthews, Yossef Zenati

Abstract: We present pre-explosion optical and infrared (IR) imaging at the site of the type II supernova (SN II) 2023ixf in Messier 101 at 6.9 Mpc. We astrometrically registered a ground-based image of SN 2023ixf to archival Hubble Space Telescope (HST), Spitzer Space Telescope (Spitzer), and ground-based near-IR images. A single point source is detected at a position consistent with the SN at wavelengths… ▽ More We present pre-explosion optical and infrared (IR) imaging at the site of the type II supernova (SN II) 2023ixf in Messier 101 at 6.9 Mpc. We astrometrically registered a ground-based image of SN 2023ixf to archival Hubble Space Telescope (HST), Spitzer Space Telescope (Spitzer), and ground-based near-IR images. A single point source is detected at a position consistent with the SN at wavelengths ranging from HST $R$-band to Spitzer 4.5 $μ$m. Fitting to blackbody and red supergiant (RSG) spectral-energy distributions (SEDs), we find that the source is anomalously cool with a significant mid-IR excess. We interpret this SED as reprocessed emission in a 8600 $R_{\odot}$ circumstellar shell of dusty material with a mass $\sim$5$\times10^{-5} M_{\odot}$ surrounding a $\log(L/L_{\odot})=4.74\pm0.07$ and $T_{\rm eff}=3920\substack{+200\\-160}$ K RSG. This luminosity is consistent with RSG models of initial mass 11 $M_{\odot}$, depending on assumptions of rotation and overshooting. In addition, the counterpart was significantly variable in pre-explosion Spitzer 3.6 $μ$m and 4.5 $μ$m imaging, exhibiting $\sim$70% variability in both bands correlated across 9 yr and 29 epochs of imaging. The variations appear to have a timescale of 2.8 yr, which is consistent with $κ$-mechanism pulsations observed in RSGs, albeit with a much larger amplitude than RSGs such as $α$ Orionis (Betelgeuse). △ Less

Submitted 7 June, 2023; originally announced June 2023.

Comments: 14 pages, 5 figures, submitted to ApJL, comments welcome

SN 2023ixf in Messier 101: Photo-ionization of Dense, Close-in Circumstellar Material in a Nearby Type II Supernova

Authors: W. V. Jacobson-Galan, L. Dessart, R. Margutti, R. Chornock, R. J. Foley, C. D. Kilpatrick, D. O. Jones, K. Taggart, C. R. Angus, S. Bhattacharjee, L. A. Braff, D. Brethauer, A. J. Burgasser, F. Cao, C. M. Carlile, K. C. Chambers, D. A. Coulter, E. Dominguez-Ruiz, C. B. Dickinson, T. de Boer, A. Gagliano, C. Gall, H. Gao, E. L. Gates, S. Gomez , et al. (43 additional authors not shown)

Abstract: We present UV/optical observations and models of supernova (SN) 2023ixf, a type II SN located in Messier 101 at 6.9 Mpc. Early-time ("flash") spectroscopy of SN 2023ixf, obtained primarily at Lick Observatory, reveals emission lines of H I, He I/II, C IV, and N III/IV/V with a narrow core and broad, symmetric wings arising from the photo-ionization of dense, close-in circumstellar material (CSM) l… ▽ More We present UV/optical observations and models of supernova (SN) 2023ixf, a type II SN located in Messier 101 at 6.9 Mpc. Early-time ("flash") spectroscopy of SN 2023ixf, obtained primarily at Lick Observatory, reveals emission lines of H I, He I/II, C IV, and N III/IV/V with a narrow core and broad, symmetric wings arising from the photo-ionization of dense, close-in circumstellar material (CSM) located around the progenitor star prior to shock breakout. These electron-scattering broadened line profiles persist for $\sim$8 days with respect to first light, at which time Doppler broadened features from the fastest SN ejecta form, suggesting a reduction in CSM density at $r \gtrsim 10^{15}$ cm. The early-time light curve of SN2023ixf shows peak absolute magnitudes (e.g., $M_{u} = -18.6$ mag, $M_{g} = -18.4$ mag) that are $\gtrsim 2$ mag brighter than typical type II supernovae, this photometric boost also being consistent with the shock power supplied from CSM interaction. Comparison of SN 2023ixf to a grid of light curve and multi-epoch spectral models from the non-LTE radiative transfer code CMFGEN and the radiation-hydrodynamics code HERACLES suggests dense, solar-metallicity, CSM confined to $r = (0.5-1) \times 10^{15}$ cm and a progenitor mass-loss rate of $\dot{M} = 10^{-2}$ M$_{\odot}$yr$^{-1}$. For the assumed progenitor wind velocity of $v_w = 50$ km s$^{-1}$, this corresponds to enhanced mass-loss (i.e., ``super-wind'' phase) during the last $\sim$3-6 years before explosion. △ Less

Submitted 21 August, 2023; v1 submitted 7 June, 2023; originally announced June 2023.

Comments: 19 pages, 8 figures. Accepted for publication in ApJL

Unprecedented X-ray Emission from the Fast Blue Optical Transient AT2022tsd

Authors: D. J. Matthews, R. Margutti, B. D. Metzger, D. Milisavljevic, G. Migliori, T. Laskar, D. Brethauer, E. Berger, R. Chornock, M. Drout, E. Ramirez-Ruiz

Abstract: We present the X-ray monitoring campaign of AT2022tsd in the time range $δt_{rest} = 23 - 116$ d rest-frame since discovery. With an initial 0.3-10 keV X-ray luminosity of $L_x \approx 10^{44}$ erg s$^{-1}$ at $δt_{rest}\approx$ 23 d, AT2022tsd is the most luminous FBOT to date and rivals even the most luminous GRBs. We find no statistical evidence for spectral evolution. The average X-ray spectru… ▽ More We present the X-ray monitoring campaign of AT2022tsd in the time range $δt_{rest} = 23 - 116$ d rest-frame since discovery. With an initial 0.3-10 keV X-ray luminosity of $L_x \approx 10^{44}$ erg s$^{-1}$ at $δt_{rest}\approx$ 23 d, AT2022tsd is the most luminous FBOT to date and rivals even the most luminous GRBs. We find no statistical evidence for spectral evolution. The average X-ray spectrum is well described by an absorbed simple power-law spectral model with best-fitting photon index $Γ= 1.89 ^{+0.09}_{-0.08}$ and marginal evidence at the 3$σ$ confidence level for intrinsic absorption $NH_{int}\approx 4\times10^{19}$ cm$^{-2}$. The X-ray light-curve behavior can be either interpreted as a power-law decay $L_x\propto t^α$ with $α\approx -2$ and superimposed X-ray variability, or as a broken power-law with a steeper post-break decay as observed in other FBOTs such as AT2018cow. We briefly compare these results to accretion models of TDEs and GRB afterglow models. △ Less

Submitted 1 June, 2023; originally announced June 2023.

Comments: 5 pages, 1 figure, 1 table. Submitted to RNAAS

A standard siren measurement of the Hubble constant using GW170817 and the latest observations of the electromagnetic counterpart afterglow

Authors: A. Palmese, R. Kaur, A. Hajela, R. Margutti, A. McDowell, A. MacFadyen

Abstract: We present a new constraint on the Hubble constant $H_0$ using the latest measurements of the electromagnetic counterpart to the gravitational wave (GW) event GW170817. We use the latest optical, X-ray and radio observations of the afterglow up to $\sim 3.5$ years after the GW detection, and properly take into account the impact of the host galaxy peculiar velocity. We find… ▽ More We present a new constraint on the Hubble constant $H_0$ using the latest measurements of the electromagnetic counterpart to the gravitational wave (GW) event GW170817. We use the latest optical, X-ray and radio observations of the afterglow up to $\sim 3.5$ years after the GW detection, and properly take into account the impact of the host galaxy peculiar velocity. We find $75.46^{+5.34}_{-5.39}$ km s$^{-1}$ Mpc$^{-1}$ (68\% Credible Interval), a $\sim7\%$ precision measurement, which is a significant improvement compared to the $14\%$ precision of the first standard siren measurement. Our result is consistent within $1σ$ with the Cepheid-anchored Supernova and within $1.5 σ$ with the Cosmic Microwave Background measurements of the Hubble constant. We also explore the impact of the various assumptions made when fitting for the afterglow on the Hubble constant estimate. △ Less

Submitted 31 May, 2023; originally announced May 2023.

Comments: 9 pages, 5 figures

Scary Barbie: An Extremely Energetic, Long-Duration Tidal Disruption Event Candidate Without a Detected Host Galaxy at z = 0.995

Authors: Bhagya M. Subrayan, Dan Milisavljevic, Ryan Chornock, Raffaella Margutti, Kate D. Alexander, Vandana Ramakrishnan, Paul C. Duffell, Danielle A. Dickinson, Kyoung-Soo Lee, Dimitrios Giannios, Geoffery Lentner, Mark Linvill, Braden Garretson, Matthew J. Graham, Daniel Stern, Daniel Brethauer, Tien Duong, Wynn Jacobson-Galán, Natalie LeBaron, David Matthews, Huei Sears, Padma Venkatraman

Abstract: We report multi-wavelength observations and characterization of the ultraluminous transient AT 2021lwx (ZTF20abrbeie; aka ``Barbie'') identified in the alert stream of the Zwicky Transient Facility (ZTF) using a Recommender Engine For Intelligent Transient Tracking (REFITT) filter on the ANTARES alert broker. From a spectroscopically measured redshift of 0.995, we estimate a peak observed pseudo-b… ▽ More We report multi-wavelength observations and characterization of the ultraluminous transient AT 2021lwx (ZTF20abrbeie; aka ``Barbie'') identified in the alert stream of the Zwicky Transient Facility (ZTF) using a Recommender Engine For Intelligent Transient Tracking (REFITT) filter on the ANTARES alert broker. From a spectroscopically measured redshift of 0.995, we estimate a peak observed pseudo-bolometric luminosity of log (L$_{\text{max}} / [\text{erg}/\text{s}]$) = 45.7 from slowly fading ztf-$\it{g}$ and ztf-$r$ light curves spanning over 1000 observer-frame days. The host galaxy is not detected in archival Pan-STARRS observations ($g > 23.3$ mag), implying a lower limit to the outburst amplitude of more than 5 mag relative to the quiescent host galaxy. Optical spectra from Lick and Keck Observatories exhibit strong emission lines with narrow cores from the H Balmer series and ultraviolet semi-forbidden lines of Si III] $λ$1892, C III] $λ$1909, and C II] $λ$2325. Typical nebular lines in AGN spectra from ions such as [O II] and [O III] are not detected. These spectral features, along with the smooth light curve that is unlike most AGN flaring activity, and the luminosity that exceeds any observed or theorized supernova, lead us to conclude that AT 2021lwx is most likely an extreme tidal disruption event (TDE). Modeling of ZTF photometry with MOSFiT suggests that the TDE was between a $\approx 14 M_{\odot}$ star and a supermassive black hole of mass $M_{\text{BH}} \sim$ $10^{8} M_{\odot}$. Continued monitoring of the still-evolving light curve along with deep imaging of the field once AT 2021lwx has faded can test this hypothesis and potentially detect the host galaxy. △ Less

Submitted 8 June, 2023; v1 submitted 21 February, 2023; originally announced February 2023.

Comments: 15 pages, 4 figures, 1 Table; Version as published in The Astrophysical Journal Letters. Observations of AT 2021lwx published in the paper can be found at https://bsubraya.github.io/research/

The Radio to GeV Afterglow of GRB 221009A

Authors: Tanmoy Laskar, Kate D. Alexander, Raffaella Margutti, Tarraneh Eftekhari, Ryan Chornock, Edo Berger, Yvette Cendes, Anne Duerr, Daniel A. Perley, Maria Edvige Ravasio, Ryo Yamazaki, Eliot H. Ayache, Thomas Barclay, Rodolfo Barniol Duran, Shivani Bhandari, Daniel Brethauer, Collin T. Christy, Deanne L. Coppejans, Paul Duffell, Wen-fai Fong, Andreja Gomboc, Cristiano Guidorzi, Jamie A. Kennea, Shiho Kobayashi, Andrew Levan , et al. (5 additional authors not shown)

Abstract: GRB 221009A ($z=0.151$) is one of the closest known long $γ$-ray bursts (GRBs). Its extreme brightness across all electromagnetic wavelengths provides an unprecedented opportunity to study a member of this still-mysterious class of transients in exquisite detail. We present multi-wavelength observations of this extraordinary event, spanning 15 orders of magnitude in photon energy from radio to… ▽ More GRB 221009A ($z=0.151$) is one of the closest known long $γ$-ray bursts (GRBs). Its extreme brightness across all electromagnetic wavelengths provides an unprecedented opportunity to study a member of this still-mysterious class of transients in exquisite detail. We present multi-wavelength observations of this extraordinary event, spanning 15 orders of magnitude in photon energy from radio to $γ$-rays. We find that the data can be partially explained by a forward shock (FS) from a highly-collimated relativistic jet interacting with a low-density wind-like medium. Under this model, the jet's beaming-corrected kinetic energy ($E_K \sim 4\times10^{50}$ erg) is typical for the GRB population. The radio and mm data provide strong limiting constraints on the FS model, but require the presence of an additional emission component. From equipartition arguments, we find that the radio emission is likely produced by a small amount of mass ($\lesssim6\times10^{-7} M_\odot$) moving relativistically ($Γ\gtrsim9$) with a large kinetic energy ($\gtrsim10^{49}$ erg). However, the temporal evolution of this component does not follow prescriptions for synchrotron radiation from a single power-law distribution of electrons (e.g. in a reverse shock or two-component jet), or a thermal electron population, perhaps suggesting that one of the standard assumptions of afterglow theory is violated. GRB 221009A will likely remain detectable with radio telescopes for years to come, providing a valuable opportunity to track the full lifecycle of a powerful relativistic jet. △ Less

Submitted 22 February, 2023; v1 submitted 8 February, 2023; originally announced February 2023.

Comments: Accepted for publication in the Astrophysical Journal Letters

GRB minimum variability timescale with Insight-HXMT and Swift: implications for progenitor models, dissipation physics and GRB classifications

Authors: A. E. Camisasca, C. Guidorzi, L. Amati, F. Frontera, X. Y. Song, S. Xiao, S. L. Xiong, S. N. Zhang, R. Margutti, S. Kobayashi, C. G. Mundell, M. Y. Ge, A. Gomboc, S. M. Jia, N. Jordana-Mitjans, C. K. Li, X. B. Li, R. Maccary, M. Shrestha, W. C. Xue, S. Zhang

Abstract: The dissipation process of GRB prompt emission is still unknown. Study of temporal variability may provide a unique way to discriminate the imprint of the inner engine activity from geometry and propagation related effects. We define the minimum variability timescale (MVT) as the shortest duration of individual pulses that shape a light curve for a sample of GRBs and test correlations with peak lu… ▽ More The dissipation process of GRB prompt emission is still unknown. Study of temporal variability may provide a unique way to discriminate the imprint of the inner engine activity from geometry and propagation related effects. We define the minimum variability timescale (MVT) as the shortest duration of individual pulses that shape a light curve for a sample of GRBs and test correlations with peak luminosity, Lorentz factor, and jet opening angle. We compare these correlations with predictions from recent numerical simulations for a relativistic structured -- possibly wobbling -- jet and assess the value of MTV as probe of prompt-emission physics. We used the peak detection algorithm mepsa to identify the shortest pulse within a GRB time history and estimate its full width half maximum (FWHM). We applied this framework to two sets of GRBs: Swift (from 2005 to July 2022) and Insight-HXMT (from June 2017 to July 2021, including 221009A). We then selected 401 GRBs with measured z to test for correlations. On average short GRBs have significantly shorter MVT than long GRBs. The MVT distribution of short GRBs with extended emission such as 060614 and 211211A is compatible only with that of short GRBs. This provides a new clue on the progenitor's nature. The MVT for long GRBs anticorrelates with peak luminosity. We confirm the anticorrelation with the Lorentz factor and find a correlation with the jet opening angle as estimated from the afterglow, along with an inverse correlation with the number of pulses. The MVT can identify the emerging putative new class of long GRBs that are suggested to be produced by compact binary mergers. For otherwise typical long GRBs, the different correlations between MVT and peak luminosity, Lorentz factor, jet opening angle, and number of pulses can be explained within the context of structured, possibly wobbling, weakly magnetised relativistic jets. (summarised) △ Less

Submitted 3 January, 2023; originally announced January 2023.

Comments: 18 pages, 15 figures, accepted by A&A

Journal ref: A&A 671, A112 (2023)

The Birth of a Relativistic Jet Following the Disruption of a Star by a Cosmological Black Hole

Authors: Dheeraj R. Pasham, Matteo Lucchini, Tanmoy Laskar, Benjamin P. Gompertz, Shubham Srivastav, Matt Nicholl, Stephen J. Smartt, James C. A. Miller-Jones, Kate D. Alexander, Rob Fender, Graham P. Smith, Michael D. Fulton, Gulab Dewangan, Keith Gendreau, Eric R. Coughlin, Lauren Rhodes, Assaf Horesh, Sjoert van Velzen, Itai Sfaradi, Muryel Guolo, N. Castro Segura, Aysha Aamer, Joseph P. Anderson, Iair Arcavi, Sean J. Brennan , et al. (41 additional authors not shown)

Abstract: A black hole can launch a powerful relativistic jet after it tidally disrupts a star. If this jet fortuitously aligns with our line of sight, the overall brightness is Doppler boosted by several orders of magnitude. Consequently, such on-axis relativistic tidal disruption events (TDEs) have the potential to unveil cosmological (redshift $z>$1) quiescent black holes and are ideal test beds to under… ▽ More A black hole can launch a powerful relativistic jet after it tidally disrupts a star. If this jet fortuitously aligns with our line of sight, the overall brightness is Doppler boosted by several orders of magnitude. Consequently, such on-axis relativistic tidal disruption events (TDEs) have the potential to unveil cosmological (redshift $z>$1) quiescent black holes and are ideal test beds to understand the radiative mechanisms operating in super-Eddington jets. Here, we present multi-wavelength (X-ray, UV, optical, and radio) observations of the optically discovered transient \target at $z=1.193$. Its unusual X-ray properties, including a peak observed luminosity of $\gtrsim$10$^{48}$ erg s$^{-1}$, systematic variability on timescales as short as 1000 seconds, and overall duration lasting more than 30 days in the rest-frame are traits associated with relativistic TDEs. The X-ray to radio spectral energy distributions spanning 5-50 days after discovery can be explained as synchrotron emission from a relativistic jet (radio), synchrotron self-Compton (X-rays), and thermal emission similar to that seen in low-redshift TDEs (UV/optical). Our modeling implies a beamed, highly relativistic jet akin to blazars but requires extreme matter-domination, i.e, high ratio of electron-to-magnetic field energy densities in the jet, and challenges our theoretical understanding of jets. △ Less

Submitted 29 November, 2022; originally announced November 2022.

Comments: To appear in Nature Astronomy on 30th November 2022. Also see here for an animation explaining the result: https://youtu.be/MQHdSbxuznY

Inferencing Progenitor and Explosion Properties of Evolving Core-collapse Supernovae from Zwicky Transient Facility Light Curves

Authors: Bhagya M. Subrayan, Danny Milisavljevic, Takashi J. Moriya, Kathryn E. Weil, Geoffrey Lentner, Mark Linvill, John Banovetz, Braden Garretson, Jack Reynolds, Niharika Sravan, Ryan Chornock, Rafaella Margutti

Abstract: We analyze a sample of 45 Type II supernovae from the Zwicky Transient Facility (ZTF) public survey using a grid of hydrodynamical models in order to assess whether theoretically-driven forecasts can intelligently guide follow up observations supporting all-sky survey alert streams. We estimate several progenitor properties and explosion physics parameters including zero-age-main-sequence (ZAMS) m… ▽ More We analyze a sample of 45 Type II supernovae from the Zwicky Transient Facility (ZTF) public survey using a grid of hydrodynamical models in order to assess whether theoretically-driven forecasts can intelligently guide follow up observations supporting all-sky survey alert streams. We estimate several progenitor properties and explosion physics parameters including zero-age-main-sequence (ZAMS) mass, mass-loss rate, kinetic energy, 56Ni mass synthesized, host extinction, and the time of explosion. Using complete light curves we obtain confident characterizations for 34 events in our sample, with the inferences of the remaining 11 events limited either by poorly constraining data or the boundaries of our model grid. We also simulate real-time characterization of alert stream data by comparing our model grid to various stages of incomplete light curves (t less than 25 days, t less than 50 days, all data), and find that some parameters are more reliable indicators of true values at early epochs than others. Specifically, ZAMS mass, time of explosion, steepness parameter beta, and host extinction are reasonably constrained with incomplete light curve data, whereas mass-loss rate, kinetic energy and 56Ni mass estimates generally require complete light curves spanning greater than 100 days. We conclude that real-time modeling of transients, supported by multi-band synthetic light curves tailored to survey passbands, can be used as a powerful tool to identify critical epochs of follow up observations. Our findings are relevant to identify, prioritize, and coordinate efficient follow up of transients discovered by Vera C. Rubin Observatory. △ Less

Submitted 28 November, 2022; originally announced November 2022.

Comments: 27 pages, 14 figures, Accepted to The Astrophysical Journal

The Young Supernova Experiment Data Release 1 (YSE DR1): Light Curves and Photometric Classification of 1975 Supernovae

Authors: P. D. Aleo, K. Malanchev, S. Sharief, D. O. Jones, G. Narayan, R. J. Foley, V. A. Villar, C. R. Angus, V. F. Baldassare, M. J. Bustamante-Rosell, D. Chatterjee, C. Cold, D. A. Coulter, K. W. Davis, S. Dhawan, M. R. Drout, A. Engel, K. D. French, A. Gagliano, C. Gall, J. Hjorth, M. E. Huber, W. V. Jacobson-Galán, C. D. Kilpatrick, D. Langeroodi , et al. (58 additional authors not shown)

Abstract: We present the Young Supernova Experiment Data Release 1 (YSE DR1), comprised of processed multi-color Pan-STARRS1 (PS1) griz and Zwicky Transient Facility (ZTF) gr photometry of 1975 transients with host-galaxy associations, redshifts, spectroscopic/photometric classifications, and additional data products from 2019 November 24 to 2021 December 20. YSE DR1 spans discoveries and observations from… ▽ More We present the Young Supernova Experiment Data Release 1 (YSE DR1), comprised of processed multi-color Pan-STARRS1 (PS1) griz and Zwicky Transient Facility (ZTF) gr photometry of 1975 transients with host-galaxy associations, redshifts, spectroscopic/photometric classifications, and additional data products from 2019 November 24 to 2021 December 20. YSE DR1 spans discoveries and observations from young and fast-rising supernovae (SNe) to transients that persist for over a year, with a redshift distribution reaching z~0.5. We present relative SN rates from YSE's magnitude- and volume-limited surveys, which are consistent with previously published values within estimated uncertainties for untargeted surveys. We combine YSE and ZTF data, and create multi-survey SN simulations to train the ParSNIP and SuperRAENN photometric classification algorithms; when validating our ParSNIP classifier on 472 spectroscopically classified YSE DR1 SNe, we achieve 82% accuracy across three SN classes (SNe Ia, II, Ib/Ic) and 90% accuracy across two SN classes (SNe Ia, core-collapse SNe). Our classifier performs particularly well on SNe Ia, with high (>90%) individual completeness and purity, which will help build an anchor photometric SNe Ia sample for cosmology. We then use our photometric classifier to characterize our photometric sample of 1483 SNe, labeling 1048 (~71%) SNe Ia, 339 (~23%) SNe II, and 96 (~6%) SNe Ib/Ic. YSE DR1 provides a training ground for building discovery, anomaly detection, and classification algorithms, performing cosmological analyses, understanding the nature of red and rare transients, exploring tidal disruption events and nuclear variability, and preparing for the forthcoming Vera C. Rubin Observatory Legacy Survey of Space and Time. △ Less

Submitted 21 February, 2023; v1 submitted 14 November, 2022; originally announced November 2022.

Comments: Accepted to ApJS; 64 pages; 35 figures; 10 tables

The Luminosity Phase Space of Galactic and Extragalactic X-ray Transients Out to Intermediate Redshifts

Authors: Ava Polzin, Raffaella Margutti, Deanne Coppejans, Katie Auchettl, Kim L. Page, Georgios Vasilopoulos, Joe S. Bright, Paolo Esposito, Peter K. G. Williams, Koji Mukai, Edo Berger

Abstract: We present a detailed compilation and analysis of the X-ray phase space of low- to intermediate-redshift ($ 0\le z \le 1$) transients that consolidates observed light curves (and theory where necessary) for a large variety of classes of transient/variable phenomena in the 0.3--10 keV energy band. We include gamma-ray burst afterglows, supernovae, supernova shock breakouts and shocks interacting wi… ▽ More We present a detailed compilation and analysis of the X-ray phase space of low- to intermediate-redshift ($ 0\le z \le 1$) transients that consolidates observed light curves (and theory where necessary) for a large variety of classes of transient/variable phenomena in the 0.3--10 keV energy band. We include gamma-ray burst afterglows, supernovae, supernova shock breakouts and shocks interacting with the environment, tidal disruption events and active galactic nuclei, fast blue optical transients, cataclysmic variables, magnetar flares/outbursts and fast radio bursts, cool stellar flares, X-ray binary outbursts, and ultraluminous X-ray sources. Our overarching goal is to offer a comprehensive resource for the examination of these ephemeral events, extending the X-ray duration-luminosity phase space (DLPS) to show luminosity evolution. We use existing observations (both targeted and serendipitous) to characterize the behavior of various transient/variable populations. Contextualizing transient signals in the larger DLPS serves two primary purposes: to identify areas of interest (i.e., regions in the parameter space where one would expect detections, but in which observations have historically been lacking) and to provide initial qualitative guidance in classifying newly discovered transient signals. We find that while the most luminous (largely extragalactic) and least luminous (largely Galactic) part of the phase space is well-populated at $t > 0.1$ days, intermediate luminosity phenomena (L$_x = 10^{34} - 10^{42}$ erg s$^{-1}$) represent a gap in the phase space. We thus identify L$_x = 10^{34} - 10^{42}$ erg s$^{-1}$ and $t = 10^{-4} - 0.1$ days as a key discovery phase space in transient X-ray astronomy. △ Less

Submitted 5 September, 2023; v1 submitted 2 November, 2022; originally announced November 2022.

Comments: 12 figures, 13 tables; version accepted to ApJ

The Jet Opening Angle and Event Rate Distributions of Short Gamma-ray Bursts from Late-time X-ray Afterglows

Authors: Alicia Rouco Escorial, Wen-fai Fong, Edo Berger, Tanmoy Laskar, Raffaella Margutti, Genevieve Schroeder, Jillian C. Rastinejad, Dylaan Cornish, Sarah Popp, Maura Lally, Anya E. Nugent, Kerry Paterson, Brian D. Metzger, Ryan Chornock, Kate Alexander, Yvette Cendes, Tarraneh Eftekhari

Abstract: We present a comprehensive study of 29 short gamma-ray bursts (SGRBs) observed $\approx 0.8-60$ days post-burst using $Chandra$ and $XMM-Newton$. We provide the inferred distributions of SGRB jet opening angles and true event rates to compare against neutron star merger rates. We perform uniform analysis and modeling of their afterglows, obtaining 10 opening angle measurements and 19 lower limits.… ▽ More We present a comprehensive study of 29 short gamma-ray bursts (SGRBs) observed $\approx 0.8-60$ days post-burst using $Chandra$ and $XMM-Newton$. We provide the inferred distributions of SGRB jet opening angles and true event rates to compare against neutron star merger rates. We perform uniform analysis and modeling of their afterglows, obtaining 10 opening angle measurements and 19 lower limits. We report on two new opening angle measurements (SGRBs 050724A and 200411A) and eight updated values, obtaining a median value of $\langle θ_{\rm j} \rangle \approx 6.1^{\circ}$ [-3.2$^{\circ}$,+9.3$^{\circ}$] (68\% confidence on the full distribution) from jet measurements alone. For the remaining events, we infer $θ_{\rm j}\gtrsim 0.5-26^{\circ}$. We uncover a population of SGRBs with wider jets of $θ_{\rm j} \gtrsim 10^{\circ}$ (including two measurements of $θ_{\rm j} \gtrsim 15^{\circ}$), representing $\sim 28\%$ of our sample. Coupled with multi-wavelength afterglow information, we derive a total true energy of $\langle E_{\rm true, tot} \rangle \approx 10^{49}-10^{50}$\,erg which is consistent with MHD jet launching mechanisms. Furthermore, we determine a range for the beaming-corrected event rate of $\mathfrak{R}_{\rm true} \approx360-1800$ Gpc$^{-3}$ yr$^{-1}$, set by the inclusion of a population of wide jets on the low end, and the jet measurements alone on the high end. From a comparison with the latest merger rates, our results are consistent with the majority of SGRBs originating from binary neutron star mergers. However, our inferred rates are well above the latest neutron star-black hole merger rates, consistent with at most a small fraction of SGRBs originating from such mergers. △ Less

Submitted 11 October, 2022; originally announced October 2022.

Comments: 32 pages, 8 figures, 4 tables (Submitted to ApJ)

Rubin Observatory LSST Transients and Variable Stars Roadmap

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

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

Submitted 8 August, 2022; originally announced August 2022.

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

Science Cases for the Keck Wide-Field Imager

Authors: J. Cooke, C. Angus, K. Auchettl, J. Bally, B. Bolin, S. Brough, J. N. Burchett, R. Foley, G. Foran, D. Forbes, J. Gannon, R. Hirai, G. G. Kacprzak, R. Margutti, C. Martinez-Lombilla, U. Mestric, A. Moller, A. Rest, J. Rhodes, R. M. Rich, F. Schussler, R. Wainscoat, J. Walawender, I. Wold, J. Zhang

Abstract: The Keck Wide-Field Imager (KWFI) is a proposed 1-degree diameter field of view UV-sensitive optical camera for Keck prime focus. KWFI will be the most powerful optical wide-field camera in the world and the only such 8m-class camera sensitive down to ~3000 A for the foreseeable future. Twenty science cases are described for KWFI compiled largely during 2019-2021, preceded by a brief discussion of… ▽ More The Keck Wide-Field Imager (KWFI) is a proposed 1-degree diameter field of view UV-sensitive optical camera for Keck prime focus. KWFI will be the most powerful optical wide-field camera in the world and the only such 8m-class camera sensitive down to ~3000 A for the foreseeable future. Twenty science cases are described for KWFI compiled largely during 2019-2021, preceded by a brief discussion of the instrument, components, and capabilities for context. △ Less

Submitted 24 July, 2022; originally announced July 2022.

Comments: 56 pages, 31 figures

Evidence for Extended Hydrogen-Poor CSM in the Three-Peaked Light Curve of Stripped Envelope Ib Supernova

Authors: Yossef Zenati, Qinan Wang, Alexey Bobrick, Lindsay DeMarchi, Hila Glanz, Mor Rozner, Armin Rest, Brian D. Metzger, Raffaella Margutti, Sebastian Gomez, Nathan Smith, Silvia Toonen, Joe S. Bright, Colin Norman, Ryan J. Foley, Alexander Gagliano, Julian H. Krolik, Stephen J. Smartt, Ashley V. Villar, Gautham Narayan, Ori Fox, Katie Auchettl, Daniel Brethauer, Alejandro Clocchiatti, Sophie V. Coelln , et al. (18 additional authors not shown)

Abstract: We present multi-band ATLAS photometry for SN 2019tsf, a stripped-envelope Type Ib supernova (SESN). The SN shows a triple-peaked light curve and a late (re-)brightening, making it unique among stripped-envelope systems. The re-brightening observations represent the latest photometric measurements of a multi-peaked Type Ib SN to date. As late-time photometry and spectroscopy suggest no hydrogen, t… ▽ More We present multi-band ATLAS photometry for SN 2019tsf, a stripped-envelope Type Ib supernova (SESN). The SN shows a triple-peaked light curve and a late (re-)brightening, making it unique among stripped-envelope systems. The re-brightening observations represent the latest photometric measurements of a multi-peaked Type Ib SN to date. As late-time photometry and spectroscopy suggest no hydrogen, the potential circumstellar material (CSM) must be H-poor. Moreover, late (>150 days) spectra show no signs of narrow emission lines, further disfavouring CSM interaction. On the contrary, an extended CSM structure is seen through a follow-up radio campaign with Karl G. Jansky Very Large Array (VLA), indicating a source of bright optically thick radio emission at late times, which is highly unusual among H-poor SESNe. We attribute this phenomenology to an interaction of the supernova ejecta with spherically-asymmetric CSM, potentially disk-like, and we present several models that can potentially explain the origin of this rare Type Ib supernova. The warped disc model paints a novel picture, where the tertiary companion perturbs the progenitors CSM, that can explain the multi-peaked light curves of SNe, and here we apply it to SN 2019tsf. This SN 2019tsf is likely a member of a new sub-class of Type Ib SNe and among the recently discovered class of SNe that undergo mass transfer at the moment of explosion △ Less

Submitted 23 July, 2022; v1 submitted 14 July, 2022; originally announced July 2022.

Comments: 23 pages, Comments are welcome, Submitted to ApJ

A Mildly Relativistic Outflow Launched Two Years after Disruption in the Tidal Disruption Event AT2018hyz

Authors: Yvette Cendes, Edo Berger, Kate Alexander, Sebastian Gomez, Aprajita Hajela, Ryan Chornock, Tanmoy Laskar, Raffaella Margutti, Brian Metzger, Michael Bietenholz, Daniel Brethauer, Mark Wieringa

Abstract: We present late-time radio/millimeter (as well as optical/UV and X-ray) detections of the tidal disruption event (TDE) AT2018hyz, spanning $970 - 1300$ d after optical discovery. In conjunction with earlier deeper limits, including at $\approx 700$ d, our observations reveal rapidly rising emission at $0.8-240$ GHz, steeper than $F_ν\propto t^5$ relative to the time of optical discovery. Such a st… ▽ More We present late-time radio/millimeter (as well as optical/UV and X-ray) detections of the tidal disruption event (TDE) AT2018hyz, spanning $970 - 1300$ d after optical discovery. In conjunction with earlier deeper limits, including at $\approx 700$ d, our observations reveal rapidly rising emission at $0.8-240$ GHz, steeper than $F_ν\propto t^5$ relative to the time of optical discovery. Such a steep rise cannot be explained in any reasonable scenario of an outflow launched at the time of disruption (e.g., off-axis jet, sudden increase in the ambient density), and instead points to a delayed launch. Our multi-frequency data allow us to directly determine the radius and energy of the radio-emitting outflow, showing that it was launched $\approx 750$ d after optical discovery. The outflow velocity is mildly relativistic, with $β\approx 0.25$ and $\approx 0.6$ for a spherical and a $10^\circ$ jet geometry, respectively, and the minimum kinetic energy is $E_K\approx 5.8\times 10^{49}$ and $\approx 6.3\times 10^{49}$ erg, respectively. This is the first definitive evidence for the production of a delayed mildly-relativistic outflow in a TDE; a comparison to the recently-published radio light curve of ASASSN-15oi suggests that the final re-brightening observed in that event (at a single frequency and time) may be due to a similar outflow with a comparable velocity and energy. Finally, we note that the energy and velocity of the delayed outflow in AT2018hyz are intermediate between those of past non-relativistic TDEs (e.g., ASASSN-14li, AT2019dsg) and the relativistic TDE Sw\,J1644+57. We suggest that such delayed outflows may be common in TDEs. △ Less

Submitted 28 June, 2022; originally announced June 2022.

Comments: Submitted to ApJ

Short GRB Host Galaxies. II. A Legacy Sample of Redshifts, Stellar Population Properties, and Implications for their Neutron Star Merger Origins

Authors: Anya E. Nugent, Wen-fai Fong, Yuxin Dong, Joel Leja, Edo Berger, Michael Zevin, Ryan Chornock, Bethany E. Cobb, Luke Zoltan Kelley, Charles D. Kilpatrick, Andrew Levan, Raffaella Margutti, Kerry Paterson, Daniel Perley, Alicia Rouco Escorial, Nathan Smith, Nial Tanvir

Abstract: We present the stellar population properties of 69 short gamma-ray burst (GRB) host galaxies, representing the largest uniformly-modeled sample to-date. Using the Prospector stellar population inference code, we jointly fit photometry and/or spectroscopy of each host galaxy. We find a population median redshift of $z=0.64^{+0.83}_{-0.32}$ ($68\%$ confidence), including 10 new or revised photometri… ▽ More We present the stellar population properties of 69 short gamma-ray burst (GRB) host galaxies, representing the largest uniformly-modeled sample to-date. Using the Prospector stellar population inference code, we jointly fit photometry and/or spectroscopy of each host galaxy. We find a population median redshift of $z=0.64^{+0.83}_{-0.32}$ ($68\%$ confidence), including 10 new or revised photometric redshifts at $z\gtrsim1$. We further find a median mass-weighted age of $t_m=0.8^{+2.71}_{-0.53}$Gyr, stellar mass of $\log(M_*/M_\odot)=9.69^{+0.75}_{-0.65}$, star formation rate of SFR=$1.44^{+9.37}_{-1.35}M_\odot$yr$^{-1}$, stellar metallicity of $\log(Z_*/Z_\odot)=-0.38^{+0.44}_{-0.42}$, and dust attenuation of $A_V=0.43^{+0.85}_{-0.36}$~mag (68\% confidence). Overall, the majority of short GRB hosts are star-forming ($\approx84\%$), with small fractions that are either transitioning ($\approx6\%$) or quiescent ($\approx10\%$); however, we observe a much larger fraction ($\approx40\%$) of quiescent and transitioning hosts at $z\lesssim0.25$, commensurate with galaxy evolution. We find that short GRB hosts populate the star-forming main sequence of normal field galaxies, but do not include as many high-mass galaxies, implying that their binary neutron star (BNS) merger progenitors are dependent on a combination of host star formation and stellar mass. The distribution of ages and redshifts implies a broad delay-time distribution, with a fast-merging channel at $z>1$ and a decreased BNS formation efficiency at lower redshifts. If short GRB hosts are representative of BNS merger hosts within the horizon of current gravitational wave detectors, these results can inform future searches for electromagnetic counterparts. All of the data and modeling products are available on the BRIGHT website. △ Less

Submitted 22 November, 2022; v1 submitted 3 June, 2022; originally announced June 2022.

Comments: 32 pages, 15 figures, 3 tables, published with ApJ

Journal ref: ApJ 940 57 (2022)

Short GRB Host Galaxies I: Photometric and Spectroscopic Catalogs, Host Associations, and Galactocentric Offsets

Authors: Wen-fai Fong, Anya E. Nugent, Yuxin Dong, Edo Berger, Kerry Paterson, Ryan Chornock, Andrew Levan, Peter Blanchard, Kate D. Alexander, Jennifer Andrews, Bethany E. Cobb, Antonino Cucchiara, Derek Fox, Chris L. Fryer, Alexa C. Gordon, Charles D. Kilpatrick, Ragnhild Lunnan, Raffaella Margutti, Adam Miller, Peter Milne, Matt Nicholl, Daniel Perley, Jillian Rastinejad, Alicia Rouco Escorial, Genevieve Schroeder , et al. (3 additional authors not shown)

Abstract: We present a comprehensive optical and near-infrared census of the fields of 90 short gamma-ray bursts (GRBs) discovered in 2005-2021, constituting all short GRBs for which host galaxy associations are feasible ($\approx$ 60% of the total Swift short GRB population). We contribute 245 new multi-band imaging observations across 49 distinct GRBs and 25 spectra of their host galaxies. Supplemented by… ▽ More We present a comprehensive optical and near-infrared census of the fields of 90 short gamma-ray bursts (GRBs) discovered in 2005-2021, constituting all short GRBs for which host galaxy associations are feasible ($\approx$ 60% of the total Swift short GRB population). We contribute 245 new multi-band imaging observations across 49 distinct GRBs and 25 spectra of their host galaxies. Supplemented by literature and archival survey data, the catalog contains 335 photometric and 40 spectroscopic data sets. The photometric catalog reaches $3σ$ depths of $\gtrsim 24-27$ mag and $\gtrsim 23-26$ mag for the optical and near-infrared bands, respectively. We identify host galaxies for 84 bursts, in which the most robust associations make up 54% (49/90) of events, while only a small fraction, 6.7%, have inconclusive host associations. Based on new spectroscopy, we determine 17 host spectroscopic redshifts with a range of $z\approx 0.15-1.6$ and find that $\approx$ 25-44% of Swift short GRBs originate from $z>1$. We also present the galactocentric offset catalog for 83 short GRBs. Taking into account the large range of individual measurement uncertainties, we find a median of projected offset of $\approx 7.9$ kpc, for which the bursts with the most robust associations have a smaller median of $\approx 4.9$ kpc. Our catalog captures more high-redshift and low-luminosity hosts, and more highly-offset bursts than previously found, thereby diversifying the population of known short GRB hosts and properties. In terms of locations and host luminosities, the populations of short GRBs with and without detectable extended emission are statistically indistinguishable. This suggests that they arise from the same progenitors, or from multiple progenitors which form and evolve in similar environments. All of the data products are available on the BRIGHT website. △ Less

Submitted 3 June, 2022; originally announced June 2022.

Comments: 53 pages, 9 figures, 6 tables, submitted

Seven years of coordinated Chandra-NuSTAR observations of SN2014C unfold the extreme mass-loss history of its stellar progenitor

Authors: Daniel Brethauer, Raffaella Margutti, Danny Milisavljevic, Michael F. Bietenholz, Ryan Chornock, Deanne L. Coppejans, Fabio De Colle, Aprajita Hajela, Giacomo Terreran, Felipe Vargas, Lindsay DeMarchi, Chelsea Harris, Wynn V. Jacobson-Galán, Atish Kamble, Daniel Patnaude, Michael C. Stroh

Abstract: We present the results from our seven-year long broad-band X-ray observing campaign of SN\,2014C with \emph{Chandra} and \emph{NuSTAR}. These coordinated observations represent the first look at the evolution of a young extragalactic SN in the 0.3-80 keV energy range in the years after core collapse. We find that the spectroscopic metamorphosis of SN\,2014C from an ordinary type Ib SN into an inte… ▽ More We present the results from our seven-year long broad-band X-ray observing campaign of SN\,2014C with \emph{Chandra} and \emph{NuSTAR}. These coordinated observations represent the first look at the evolution of a young extragalactic SN in the 0.3-80 keV energy range in the years after core collapse. We find that the spectroscopic metamorphosis of SN\,2014C from an ordinary type Ib SN into an interacting SN with copious hydrogen emission is accompanied by luminous X-rays reaching $L_x\approx 5.6\times10^{40}\, \rm{erg\,s^{-1}}$ (0.3--100 keV) at $\sim 1000$ days post explosion and declining as $L_x\propto t^{-1}$ afterwards. The broad-band X-ray spectrum is of thermal origin and shows clear evidence for cooling after peak, with $T(t)\approx 20 \,{\rm keV}(t/t_{\rm pk})^{-0.5}$. Soft X-rays of sub-keV energy suffer from large photoelectric absorption originating from the local SN environment with $NH_{\rm int}(t)\approx3\times 10^{22}(t/400 \,\rm{days})^{-1.4}\,\rm{cm^{-2}}$. We interpret these findings as the result of the interaction of the SN shock with a dense ($n\approx 10^{5}-10^{6}\,\rm{cm^{-3}}$), H-rich disk-like circumstellar medium (CSM) with inner radius $\sim2\times 10^{16}$ cm and extending to $\sim 10^{17}$ cm. Based on the declining $NH_{\rm int}(t)$ and X-ray luminosity evolution, we infer a CSM mass of $\sim(1.2\,f$--2.0$\sqrt{f}) \rm{M_{\odot}}$, where $f$ is the volume filling factor. Finally, we place SN\,2014C in the context of 119 core-collapse SNe with evidence for strong shock interaction with a thick circumstellar medium and we highlight the challenges that the current mass-loss theories (including wave-driven mass loss, binary interaction and line-driven winds) face when interpreting the wide dynamic ranges of CSM parameters inferred from observations. △ Less

Submitted 1 June, 2022; originally announced June 2022.

Comments: We recognize that the supernova within our sample is non-exhaustive, please email daniel_brethauer@berkeley.edu for suggestions of additional core-collapse interacting supernovae

Neutrino Emission from Luminous Fast Blue Optical Transients

Authors: Ersilia Guarini, Irene Tamborra, Raffaella Margutti

Abstract: Mounting evidence suggests that Luminous Fast Blue Optical Transients (LFBOTs) are powered by a compact object, launching an asymmetric and fast outflow responsible for the radiation observed in the ultraviolet, optical, infrared, radio, and X-ray bands. Proposed scenarios aiming to explain the electromagnetic emission include an inflated cocoon, surrounding a jet choked in the extended stellar en… ▽ More Mounting evidence suggests that Luminous Fast Blue Optical Transients (LFBOTs) are powered by a compact object, launching an asymmetric and fast outflow responsible for the radiation observed in the ultraviolet, optical, infrared, radio, and X-ray bands. Proposed scenarios aiming to explain the electromagnetic emission include an inflated cocoon, surrounding a jet choked in the extended stellar envelope. In alternative, the observed radiation may arise from the disk formed by the delayed merger of a black hole with a Wolf-Rayet star. We explore the neutrino production in these scenarios, i.e. internal shocks in a choked jet and interaction between the outflow and the circumstellar medium (CSM). If observed on-axis, the choked jet provides the dominant contribution to the neutrino fluence. Intriguingly, the IceCube upper limit on the neutrino emission inferred from the closest LFBOT, AT2018cow, excludes a region of the parameter space otherwise allowed by electromagnetic observations. After correcting for the Eddington bias on the observation of cosmic neutrinos, we conclude that the emission from an on-axis choked jet and CSM interaction is compatible with the detection of two track-like neutrino events observed by the IceCube Neutrino Observatory in coincidence with AT2018cow, and otherwise considered to be of atmospheric origin. While the neutrino emission from LFBOTs does not constitute the bulk of the diffuse background of neutrinos observed by IceCube, detection prospects of nearby LFBOTs with IceCube and the upcoming IceCube-Gen2 are encouraging. Follow-up neutrino searches will be crucial for unravelling the mechanism powering this emergent transient class. △ Less

Submitted 7 July, 2022; v1 submitted 24 May, 2022; originally announced May 2022.

Comments: 23 pages, including 9 figures and 1 appendix. Minor changes, matches version accepted for publication in ApJ

Journal ref: Astrophys. J. 935 (2022) 157

The First Short GRB Millimeter Afterglow: The Wide-Angled Jet of the Extremely Energetic SGRB 211106A

Authors: Tanmoy Laskar, Alicia Rouco Escorial, Genevieve Schroeder, Wen-fai Fong, Edo Berger, Péter Veres, Shivani Bhandari, Jillian Rastinejad, Charles D. Kilpatrick, Aaron Tohuvavohu, Raffaella Margutti, Kate D. Alexander, James DeLaunay, Jamie A. Kennea, Anya Nugent, K. Paterson, Peter K. G. Williams

Abstract: We present the discovery of the first millimeter afterglow of a short-duration $γ$-ray burst (SGRB) and the first confirmed afterglow of an SGRB localized by the GUANO system on Swift. Our Atacama Large Millimeter/Sub-millimeter Array (ALMA) detection of SGRB 211106A establishes an origin in a faint host galaxy detected in Hubble Space Telescope (HST) imaging at $0.7\lesssim z\lesssim1.4$. From th… ▽ More We present the discovery of the first millimeter afterglow of a short-duration $γ$-ray burst (SGRB) and the first confirmed afterglow of an SGRB localized by the GUANO system on Swift. Our Atacama Large Millimeter/Sub-millimeter Array (ALMA) detection of SGRB 211106A establishes an origin in a faint host galaxy detected in Hubble Space Telescope (HST) imaging at $0.7\lesssim z\lesssim1.4$. From the lack of a detectable optical afterglow, coupled with the bright millimeter counterpart, we infer a high extinction, $A_{\rm V}\gtrsim2.6$ mag along the line of sight, making this the one of the most highly dust-extincted SGRBs known to date. The millimeter-band light curve captures the passage of the synchrotron peak from the afterglow forward shock and reveals a jet break at $t_{\rm jet}=29.2^{+4.5}_{-4.0}$~days. For a presumed redshift of $z=1$, we infer an opening angle, $θ_{\rm jet}=(15.5\pm1.4)$~degrees, and beaming-corrected kinetic energy of $\log(E_{\rm K}/{\rm erg})=51.8\pm0.3$, making this one of the widest and most energetic SGRB jets known to date. Combining all published millimeter-band upper limits in conjunction with the energetics for a large sample of SGRBs, we find that energetic outflows in high density environments are more likely to have detectable millimeter counterparts. Concerted afterglow searches with ALMA should yield detection fractions of 24-40% on timescales of $\gtrsim2$~days at rates $\approx0.8$-1.6 per year, outpacing the historical discovery rate of SGRB centimeter-band afterglows. △ Less

Submitted 1 August, 2022; v1 submitted 6 May, 2022; originally announced May 2022.

Comments: 20 pages, 9 figures, 4 tables. Version accepted for publication in ApJ Letters

Radio Analysis of SN 2004C Reveals an Unusual CSM Density Profile as a Harbinger of Core Collapse

Authors: Lindsay DeMarchi, R. Margutti, J. Dittman, A. Brunthaler, D. Milisavljevic, Michael F. Bietenholz, C. Stauffer, D. Brethauer, D. Coppejans, K. Auchettl, K. D. Alexander, C. D. Kilpatrick, Joe S. Bright, L. Z. Kelley, Michael C. Stroh, W. V. Jacobson-Galan

Abstract: We present extensive multi-frequency VLA and VLBA observations of the radio-bright supernova (SN) IIb SN 2004C that span $\sim(40-2793)$ days post-explosion. We interpret the temporal evolution of the radio spectral energy distribution (SED) in the context of synchrotron self-absorbed (SSA) emission from the explosion's forward shock as it expands in the circumstellar medium (CSM) previously sculp… ▽ More We present extensive multi-frequency VLA and VLBA observations of the radio-bright supernova (SN) IIb SN 2004C that span $\sim(40-2793)$ days post-explosion. We interpret the temporal evolution of the radio spectral energy distribution (SED) in the context of synchrotron self-absorbed (SSA) emission from the explosion's forward shock as it expands in the circumstellar medium (CSM) previously sculpted by the mass-loss history of the stellar progenitor. VLBA observations and modeling of the VLA data point to a blastwave with average velocity $\sim0.06c$ that carries an energy of $\sim 10^{49}$ erg. Our modeling further reveals a flat CSM density profile $ρ_{\rm{CSM}} \propto R^{-0.03 \pm0.22}$ up to a break radius $R_{br} \approx (1.96 \pm 0.10) \times 10^{16}$ cm, with a steep density gradient following $ρ_{\rm{CSM}} \propto R^{-2.3 \pm 0.5}$ at larger radii. We infer that the flat part of the density profile corresponds to a CSM shell with mass $\sim0.021 M_{\odot}$, and that the progenitor's effective mass-loss rate varied with time over the range $(50-500) \times 10^{-5} M_{\odot} \rm{yr}^{-1}$ for an adopted wind velocity $v_w =1000$ km $s^{-1}$ and shock microphysical parameters $ε_e = 0.1, ε_B = 0.01$. These results add to the mounting observational evidence for departures from the traditional single-wind mass-loss scenarios in evolved, massive stars in the centuries leading up to core collapse. Potentially viable scenarios include mass loss powered by gravity waves and/or interaction with a binary companion. △ Less

Submitted 14 March, 2022; originally announced March 2022.

Comments: 32 pages, 10 figures, submitted to ApJ

The Circumstellar Environments of Double-Peaked, Calcium-strong Supernovae 2021gno and 2021inl

Authors: Wynn Jacobson-Galán, Padma Venkatraman, Raffaella Margutti, David Khatami, Giacomo Terreran, Ryan J. Foley, Rodrigo Angulo, Charlotte R. Angus, Katie Auchettl, Peter K. Blanchard, Alexey Bobrick, Joe S. Bright, Cirilla D. Couch, David A. Coulter, Karoli Clever, Kyle W. Davis, Thomas de Boer, Lindsay DeMarchi, Sierra A. Dodd, David O. Jones, Jessica Johnson, Charles D. Kilpatrick, Nandita Khetan, Zhisen Lai, Danial Langeroodi , et al. (20 additional authors not shown)

Abstract: We present panchromatic observations and modeling of calcium-strong supernovae (SNe) 2021gno in the star-forming host galaxy NGC 4165 (D = 30.5 Mpc) and 2021inl in the outskirts of elliptical galaxy NGC 4923 (D = 80 Mpc), both monitored through the Young Supernova Experiment (YSE) transient survey. The multi-color light curves of both SNe show two peaks, the former peak being derived from shock co… ▽ More We present panchromatic observations and modeling of calcium-strong supernovae (SNe) 2021gno in the star-forming host galaxy NGC 4165 (D = 30.5 Mpc) and 2021inl in the outskirts of elliptical galaxy NGC 4923 (D = 80 Mpc), both monitored through the Young Supernova Experiment (YSE) transient survey. The multi-color light curves of both SNe show two peaks, the former peak being derived from shock cooling emission (SCE) and/or shock interaction with circumstellar material (CSM). The primary peak in SN 2021gno is coincident with luminous, rapidly decaying X-ray emission ($L_x = 5 \times 10^{41}$ erg s$^{-1}$) detected by Swift-XRT at $δt = 1$ day after explosion, this observation being the second ever detection of X-rays from a calcium-strong transient. We interpret the X-ray emission from SN 2021gno in the context of shock interaction with dense CSM that extends to $r < 3 \times 10^{14}$ cm. Based on modeling of the SN 2021gno X-ray spectrum, we calculate a CSM mass range of $M_{\rm CSM} = (0.3 - 1.6) \times 10^{-3}$ M$_{\odot}$ and particle densities of $n = (1-4) \times 10^{10}$ cm$^{-3}$. Radio non-detections of SN 2021gno indicate a low-density environment at larger radii ($r > 10^{16}$ cm) and a progenitor mass loss rate of $\dot{M} < 10^{-4}$ M$_{\odot}$ yr$^{-1}$, for $v_w = 500$ km s$^{-1}$. For radiation derived from SCE, modeling of the primary light curve peak in both SNe indicates an extended progenitor envelope mass and radius of $M_e = 0.02 - 0.05$ M$_{\odot}$ and $R_e = 30 - 230$ R$_{\odot}$. The explosion properties of SNe 2021gno and 2021inl suggest progenitor systems containing either a low-mass massive star or a white dwarf (WD), the former being unlikely for either object given the lack of star formation at both explosion sites. Furthermore, the progenitor environments of both SNe are consistent with explosion models for low-mass hybrid He/C/O WD + C/O WD binaries. △ Less

Submitted 7 March, 2022; originally announced March 2022.

Comments: 33 pages, 19 figures. Submitted to ApJ. Comments welcome!

Physical Properties of the Host Galaxies of Ca-rich Transients

Authors: Y. Dong, D. Milisavljevic, J. Leja, S. K. Sarbadhicary, A. E. Nugent, R. Margutti, W. V. Jacobson-Galan, A. Polin, J. Banovetz, J. M. Reynolds, B. Subrayan

Abstract: Calcium-rich (Ca-rich) transients are a new class of supernovae (SNe) that are known for their comparatively rapid evolution, modest peak luminosities, and strong nebular calcium emission lines. Currently, the progenitor systems of Ca-rich transients remain unknown. Although they exhibit spectroscopic properties not unlike core-collapse Type Ib/c SNe, nearly half are found in the outskirts of thei… ▽ More Calcium-rich (Ca-rich) transients are a new class of supernovae (SNe) that are known for their comparatively rapid evolution, modest peak luminosities, and strong nebular calcium emission lines. Currently, the progenitor systems of Ca-rich transients remain unknown. Although they exhibit spectroscopic properties not unlike core-collapse Type Ib/c SNe, nearly half are found in the outskirts of their host galaxies that are predominantly elliptical, suggesting a closer connection to the older stellar populations of SNe Ia. In this paper, we present a compilation of publicly available multiwavelength observations of all known and/or suspected host galaxies of Ca-rich transients ranging from FUV to IR, and use these data to characterize their stellar populations with prospector. We estimate several galaxy parameters including integrated star formation rate, stellar mass, metallicity, and age. For nine host galaxies, the observations are sensitive enough to obtain nonparametric star formation histories, from which we recover SN rates and estimate probabilities that the Ca-rich transients in each of these host galaxies originated from a core-collapse vs. Type Ia-like explosion. Our work supports the notion that the population of Ca-rich transients do not come exclusively from core-collapse explosions, and must either be only from white dwarf stars or a mixed population of white dwarf stars with other channels, potentially including massive star explosions. Additional photometry and explosion site spectroscopy of larger samples of Ca-rich host galaxies will improve these estimates and better constrain the ratio of white dwarf vs. massive star progenitors of Ca-rich transients. △ Less

Submitted 7 February, 2022; originally announced February 2022.

Comments: 22 pages, 7 figures, 5 tables, accepted to ApJ

Shocked jets in CCSNe can power the zoo of fast blue optical transients

Authors: Ore Gottlieb, Alexander Tchekhovskoy, Raffaella Margutti

Abstract: Evidence is mounting that recent multiwavelength detections of fast blue optical transients (FBOTs) in star-forming galaxies comprise a new class of transients, whose origin is yet to be understood. We show that hydrogen-rich collapsing stars that launch relativistic jets near the central engine can naturally explain the entire set of FBOT observables. The jet-star interaction forms a mildly-relat… ▽ More Evidence is mounting that recent multiwavelength detections of fast blue optical transients (FBOTs) in star-forming galaxies comprise a new class of transients, whose origin is yet to be understood. We show that hydrogen-rich collapsing stars that launch relativistic jets near the central engine can naturally explain the entire set of FBOT observables. The jet-star interaction forms a mildly-relativistic shocked jet (inner cocoon) component, which powers cooling emission that dominates the high velocity optical signal during the first few weeks, with a typical energy of $ \sim 10^{50}-10^{51} $ erg. During this time, the cocoon radial energy distribution implies that the optical lightcurve exhibits a fast decay of $ L \propto t^{-2.4} $. After a few weeks, when the velocity of the emitting shell is $ \sim 0.01 $ c, the cocoon becomes transparent, and the cooling envelope governs the emission. The interaction between the cocoon and the dense circumstellar winds generates synchrotron self-absorbed emission in the radio bands, featuring a steady rise on a month timescale. After a few months the relativistic outflow decelerates, enters the observer's line of sight, and powers the peak of the radio lightcurve, which rapidly decays thereafter. The jet (and the inner cocoon) become optically thin to X-rays $ \sim $ day after the collapse, allowing X-ray photons to diffuse from the central engine that launched the jet to the observer. Cocoon cooling emission is expected at higher volumetric rates than gamma-ray bursts (GRBs) by a factor of a few, similar to FBOTs. We rule out uncollimated outflows, however both GRB jets and failed collimated jets are compatible with all observables. △ Less

Submitted 30 March, 2022; v1 submitted 12 January, 2022; originally announced January 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

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

Progenitor and Close-In Circumstellar Medium of Type II Supernova 2020fqv from High-Cadence Photometry and Ultra-Rapid UV Spectroscopy

Authors: Samaporn Tinyanont, Ryan Ridden-Harper, Ryan Foley, Viktoriya Morozova, Charles Kilpatrick, Georgios Dimitriadis, Lindsay DeMarchi, Alexander Gagliano, Wynn V. Jacobson-Galán, Alexander Messick, Justin Pierel, Anthony Piro, Enrico Ramirez-Ruiz, Matthew Siebert, Kenneth Chambers, Karoli Clever, David Coulter, Kishalay De, Matthew Hankins, Tiara Hung, Saurabh Jha, Camilo Jiménez-Ángel, David Jones, Mansi Kasliwal, Chien-Cheng Lin , et al. (13 additional authors not shown)

Abstract: We present observations of SN 2020fqv, a Virgo-cluster Type II core-collapse supernova (CCSN) with a high temporal resolution light curve from the Transiting Exoplanet Survey Satellite (TESS) covering the time of explosion; ultraviolet (UV) spectroscopy from the Hubble Space Telescope (HST) starting 3.3 days post-explosion; ground-based spectroscopic observations starting 1.1~days post-explosion;… ▽ More We present observations of SN 2020fqv, a Virgo-cluster Type II core-collapse supernova (CCSN) with a high temporal resolution light curve from the Transiting Exoplanet Survey Satellite (TESS) covering the time of explosion; ultraviolet (UV) spectroscopy from the Hubble Space Telescope (HST) starting 3.3 days post-explosion; ground-based spectroscopic observations starting 1.1~days post-explosion; along with extensive photometric observations. Massive stars have complicated mass-loss histories leading up to their death as CCSNe, creating circumstellar medium (CSM) with which the SNe interact. Observations during the first few days post-explosion can provide important information about the mass-loss rate during the late stages of stellar evolution. Model fits to the quasi-bolometric light curve of SN 2020fqv reveal ~0.23 $M_{\odot}$ of CSM confined within ~1450 $R_{\odot}$ ($10^{14}$ cm) from its progenitor star. Early spectra (<4 days post-explosion), both from HST and ground-based observatories, show emission features from high-ionization metal species from the outer, optically thin part of this CSM. We find that the CSM is consistent with an eruption caused by the injection of $\sim$$5\times 10^{46}$ erg into the stellar envelope $\sim$300 days pre-explosion, potentially from a nuclear burning instability at the onset of oxygen burning. Light-curve fitting, nebular spectroscopy, and pre-explosion \textit{HST} imaging consistently point to a red supergiant (RSG) progenitor with $M_{\rm ZAMS}$$\approx$$13.5$--$15 \, M_{\odot}$, typical for SN~II progenitor stars. This finding demonstrates that a typical RSG, like the progenitor of SN 2020fqv, has a complicated mass-loss history immediately before core collapse. △ Less

Submitted 20 October, 2021; originally announced October 2021.

Comments: Published in MNRAS

Radio and X-ray observations of the luminous Fast Blue Optical Transient AT2020xnd

Authors: Joe S. Bright, Raffaella Margutti, David Matthews, Daniel Brethauer, Deanne Coppejans, Mark H. Wieringa, Brian D. Metzger, Lindsay DeMarchi, Tanmoy Laskar, Charles Romero, Kate D. Alexander, Assaf Horesh, Giulia Migliori, Ryan Chornock, E. Berger, Michael Bietenholz, Mark J. Devlin, Simon R. Dicker, W. V. Jacobson-Galán, Brian S. Mason, Dan Milisavljevic, Sara E. Motta, Tony Mroczkowski, Enrico Ramirez-Ruiz, Lauren Rhodes , et al. (3 additional authors not shown)

Abstract: We present deep X-ray and radio observations of the Fast Blue Optical Transient (FBOT) AT2020xnd/ZTF20acigmel at $z=0.2433$ from $13$d to $269$d after explosion. AT2020xnd belongs to the category of optically luminous FBOTs with similarities to the archetypal event AT2018cow. AT2020xnd shows luminous radio emission reaching $L_ν\approx8\times10^{29}$ergs$^{-1}$Hz$^{-1}$ at 20GHz and $75$d post exp… ▽ More We present deep X-ray and radio observations of the Fast Blue Optical Transient (FBOT) AT2020xnd/ZTF20acigmel at $z=0.2433$ from $13$d to $269$d after explosion. AT2020xnd belongs to the category of optically luminous FBOTs with similarities to the archetypal event AT2018cow. AT2020xnd shows luminous radio emission reaching $L_ν\approx8\times10^{29}$ergs$^{-1}$Hz$^{-1}$ at 20GHz and $75$d post explosion, accompanied by luminous and rapidly fading soft X-ray emission peaking at $L_{X}\approx6\times10^{42}$ergs$^{-1}$. Interpreting the radio emission in the context of synchrotron radiation from the explosion's shock interaction with the environment we find that AT2020xnd launched a high-velocity outflow ($v\sim$0.1-0.2$c$) propagating into a dense circumstellar medium (effective $\dot M\approx10^{-3}M_{\rm{sol}}$yr$^{-1}$ for an assumed wind velocity of $v_w=1000$kms$^{-1}$). Similar to AT2018cow, the detected X-ray emission is in excess compared to the extrapolated synchrotron spectrum and constitutes a different emission component, possibly powered by accretion onto a newly formed black hole or neutron star. These properties make AT2020xnd a high-redshift analog to AT2018cow, and establish AT2020xnd as the fourth member of the class of optically-luminous FBOTs with luminous multi-wavelength counterparts. △ Less

Submitted 11 October, 2021; originally announced October 2021.

Comments: 21 pages, 6 figures, 6 tables. Submitted to ApJ

Final Moments I: Precursor Emission, Envelope Inflation, and Enhanced Mass loss Preceding the Luminous Type II Supernova 2020tlf

Authors: Wynn Jacobson-Galán, Luc Dessart, David Jones, Raffaella Margutti, Deanne Coppejans, Georgios Dimitriadis, Ryan J. Foley, Charles D. Kilpatrick, David J. Matthews, Sofia Rest, Giacomo Terreran, Patrick D. Aleo, Katie Auchettl, Peter K. Blanchard, David A. Coulter, Kyle W. Davis, Thomas de Boer, Lindsay DeMarchi, Maria R. Drout, Nicholas Earl, Alexander Gagliano, Christa Gall, Jens Hjorth, Mark E. Huber, Adaeze L. Ibik , et al. (11 additional authors not shown)

Abstract: We present panchromatic observations and modeling of supernova (SN) 2020tlf, the first normal type II-P/L SN with confirmed precursor emission, as detected by the Young Supernova Experiment transient survey with the Pan-STARRS1 telescope. Pre-explosion emission was detected in $riz-$bands at 130 days prior to SN 2020tlf and persisted at relatively constant flux until first light. Soon after discov… ▽ More We present panchromatic observations and modeling of supernova (SN) 2020tlf, the first normal type II-P/L SN with confirmed precursor emission, as detected by the Young Supernova Experiment transient survey with the Pan-STARRS1 telescope. Pre-explosion emission was detected in $riz-$bands at 130 days prior to SN 2020tlf and persisted at relatively constant flux until first light. Soon after discovery, "flash" spectroscopy of SN 2020tlf revealed prominent narrow symmetric emission lines ($v_w < 300$ km s$^{-1}$) that resulted from the photo-ionization of unshocked circumstellar material (CSM) shedded in progenitor mass loss episodes in the final weeks to months before explosion. Surprisingly, this novel display of pre-SN emission and associated mass loss occurred in a RSG progenitor with ZAMS mass of only 10-12 M$_{\odot}$, as inferred from nebular spectra. Modeling of the light curve and multi-epoch spectra with the non-LTE radiative transfer code CMFGEN and radiation-hydrodynamical (RHD) code HERACLES suggests a dense CSM limited to $r \approx 10^{15}$ cm, and mass loss rate of $10^{-2}$ M$_{\odot}$ yr$^{-1}$. The subsequent luminous light-curve plateau and persistent blue excess indicates an extended progenitor, compatible with a RSG model with $R_{\star} = 1100$ R$_{\odot}$. Inferences from the limits on the shock-powered X-ray and radio luminosity are consistent with these conclusions and suggest a CSM density of $ρ< 2 \times 10^{-16}$ g cm$^{-3}$ for distances of $r \approx 5 \times 10^{15}$ cm, as well as a mass loss rate of $\dot M<1.3 \times 10^{-5}\,\rm{M_{\odot}\,yr^{-1}}$ at larger distances. A promising power source for the observed precursor emission is the ejection of stellar material following energy disposition into the stellar envelope as a result of gravity waves emitted during either neon/oxygen burning or a nuclear flash from silicon combustion. △ Less

Submitted 13 December, 2021; v1 submitted 24 September, 2021; originally announced September 2021.

Comments: 36 pages, 20 figures. Accepted for publication in ApJ