Showing 1–50 of 267 results for author: Metzger, B D

Fragmention in Gravitationally-Unstable Collapsar Disks and Sub-Solar Neutron Star Mergers

Authors: Brian D. Metzger, Lam Hui, Matteo Cantiello

Abstract: Although stable neutron stars (NS) can in principle exist down to masses Mns ~ 0.1Msun, standard models of stellar core-collapse predict a robust lower limit Mns >~ 1.2Msun, roughly commensurate with the Chandrasekhar mass Mch of the progenitor's iron core (electron fraction Ye ~ 0.5). However, this limit may be circumvented in sufficiently dense neutron-rich environments (Ye << 0.5) for which Mch… ▽ More Although stable neutron stars (NS) can in principle exist down to masses Mns ~ 0.1Msun, standard models of stellar core-collapse predict a robust lower limit Mns >~ 1.2Msun, roughly commensurate with the Chandrasekhar mass Mch of the progenitor's iron core (electron fraction Ye ~ 0.5). However, this limit may be circumvented in sufficiently dense neutron-rich environments (Ye << 0.5) for which Mch ~ Ye^2 is reduced to < Msun. Such physical conditions could arise in the black hole accretion disks formed from the collapse of rapidly-rotating stars (``collapsars''), as a result of gravitational instabilities and cooling-induced fragmentation, similar to models for planet formation in protostellar disks. We confirm that the conditions to form sub-solar mass NS (ssNS) may be marginally satisfied in the outer regions of massive neutrino-cooled collapsar disks. If the disk fragments into multiple ssNS, their subsequent coalescence offers a channel for precipitating sub-solar mass LIGO/Virgo gravitational-wave mergers that does not implicate primordial black holes. The model makes several additional predictions: (1) ~Hz frequency Doppler modulation of the ssNS-merger gravitational wave signals due to the binary's orbital motion in the disk; (2) at least one additional gravitational wave event (coincident within <~ hours), from the coalescence of the ssNS-merger remnant(s) with the central black hole; (3) an associated gamma-ray burst and supernova counterpart, the latter boosted in energy and enriched with r-process elements from the NS merger(s) embedded within the exploding stellar envelope (``kilonovae inside a supernova''). △ Less

Submitted 10 July, 2024; originally announced July 2024.

Comments: 9 pages, 2 figures

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

High-energy neutrino signatures from pulsar remnants of binary neutron-star mergers: coincident detection prospects with gravitational waves

Authors: Mainak Mukhopadhyay, Shigeo S. Kimura, Brian D. Metzger

Abstract: Binary neutron-star (BNS) mergers are accompanied by multi-messenger emissions, including gravitational wave (GW), neutrino, and electromagnetic signals. Some fraction of BNS mergers may result in a rapidly spinning magnetar as a remnant, which can enhance both the EM and neutrino emissions. In this study, we model the neutrino emissions from such systems and discuss the prospects for detecting th… ▽ More Binary neutron-star (BNS) mergers are accompanied by multi-messenger emissions, including gravitational wave (GW), neutrino, and electromagnetic signals. Some fraction of BNS mergers may result in a rapidly spinning magnetar as a remnant, which can enhance both the EM and neutrino emissions. In this study, we model the neutrino emissions from such systems and discuss the prospects for detecting the neutrinos coincident with GW signatures. We consider a scenario where a magnetar remnant drives a pulsar wind using its spin energy. The wind interacts with the surrounding kilonova ejecta, forming a nebula filled with non-thermal photons. Ions and nuclei extracted from the magnetar's surface can be accelerated in the polar-cap and the termination-shock regions. We investigate the neutrino fluences resulting from photomeson interactions, where accelerated CR protons interact with the photons in the nebula. Our findings indicate that the peak neutrino fluence is $\sim 10^{-2}\rm GeV~cm^{-2}$ for a source at $40$ Mpc, which is reached approximately $\mathcal{O}( 1-10\ {\rm days})$ post merger. Finally, we examine the potential for GW-triggered stacking searches with IceCube-Gen2 using next-generation GW detectors such as the Cosmic Explorer (CE) and the Einstein Telescope (ET). We conclude that, assuming an optimistic neutrino emission model, a combination of CE+ET would offer a high probability of neutrino detection from these sources within an operational timescale of $\sim 20$ years. In case of non-detection, $2 σ$ level constraints on model parameters can be established within similar joint operation timescales. △ Less

Submitted 5 July, 2024; originally announced July 2024.

Comments: 21 pages, 8 figures

A Radio Study of Persistent Radio Sources in Nearby Dwarf Galaxies: Implications for Fast Radio Bursts

Authors: Y. Dong, T. Eftekhari, W. Fong, S. Bhandari, E. Berger, O. S. Ould-Boukattine, J. W. T. Hessels, N. Sridhar, A. Reines, B. Margalit, J. Darling, A. C. Gordon, J. E. Greene, C. D. Kilpatrick, B. Marcote, B. D. Metzger, K. Nimmo, A. E. Nugent, Z. Paragi, P. K. G. Williams

Abstract: We present 1 - 12 GHz Karl G. Jansky Very Large Array observations of 9 off-nuclear persistent radio sources (PRSs) in nearby (z < 0.055) dwarf galaxies, along with high-resolution European very-long baseline interferometry (VLBI) Network (EVN) observations for one of them at 1.7GHz. We explore the plausibility that these PRSs are associated with fast radio burst (FRB) sources by examining their p… ▽ More We present 1 - 12 GHz Karl G. Jansky Very Large Array observations of 9 off-nuclear persistent radio sources (PRSs) in nearby (z < 0.055) dwarf galaxies, along with high-resolution European very-long baseline interferometry (VLBI) Network (EVN) observations for one of them at 1.7GHz. We explore the plausibility that these PRSs are associated with fast radio burst (FRB) sources by examining their properties, physical sizes, host-normalized offsets, spectral energy distributions (SEDs), radio luminosities, and light curves, and compare them to those of the PRSs associated with FRBs 20121102A and 20190520B, two known active galactic nuclei (AGN), and one likely AGN in our sample with comparable data, as well as other radio transients exhibiting characteristics analogous to FRB-PRSs. We identify a single source in our sample, J1136+2643, as the most promising FRB- PRS, based on its compact physical size and host-normalized offset. We further identify two sources, J0019+1507 and J0909+5955, with physical sizes comparable to FRB-PRSs, but which exhibit large offsets and flat spectral indices potentially indicative of a background AGN origin. We test the viability of neutron star wind nebulae and hypernebulae models for J1136+2643, and find that the physical size, luminosity, and SED of J1136+2643 are broadly consistent with these models. Finally, we discuss the alternative interpretation that the radio sources are instead powered by accreting massive black holes and outline future prospects and follow-up observations for differentiating between these scenarios. △ Less

Submitted 1 May, 2024; originally announced May 2024.

Comments: 24 pages, 7 figures, 3 tables

Coupled Disk-Star Evolution in Galactic Nuclei and the Lifetimes of QPE Sources

Authors: Itai Linial, Brian D. Metzger

Abstract: A modest fraction of the stars in galactic nuclei fed towards the central supermassive black hole (SMBH) approach on low-eccentricity orbits driven by gravitational-wave radiation (extreme mass ratio inspiral, EMRI). In the likely event that a gaseous accretion disk is created in the nucleus during this slow inspiral (e.g., via an independent tidal-disruption event; TDE), star-disk collisions gene… ▽ More A modest fraction of the stars in galactic nuclei fed towards the central supermassive black hole (SMBH) approach on low-eccentricity orbits driven by gravitational-wave radiation (extreme mass ratio inspiral, EMRI). In the likely event that a gaseous accretion disk is created in the nucleus during this slow inspiral (e.g., via an independent tidal-disruption event; TDE), star-disk collisions generate regular short-lived flares consistent with the observed quasi-periodic eruption (QPE) sources. We present a model for the coupled star-disk evolution which self-consistently accounts for mass and thermal energy injected into the disk from stellar collisions and associated mass ablation. For weak collision/ablation heating, the disk is thermally-unstable and undergoes limit-cycle oscillations which modulate its properties and lead to accretion-powered outbursts on timescales of years to decades, with a time-averaged accretion rate $\sim 0.1 \dot{M}_{\rm Edd}$. Stronger collision/ablation heating acts to stabilize the disk, enabling roughly steady accretion at the EMRI-stripping rate. In either case, the stellar destruction time through ablation, and hence the maximum QPE lifetime, is $\sim 10^{2}-10^{3}$ yr, far longer than fall-back accretion after a TDE. The quiescent accretion disks in QPE sources may at the present epoch be self-sustaining and fed primarily by EMRI ablation. Indeed, the observed range of secular variability broadly match those predicted for collision-fed disks. Changes in the QPE recurrence pattern following such outbursts, similar to that observed in GSN 069, could arise from temporary misalignment between the EMRI-fed disk and the SMBH equatorial plane as the former regrows its mass after a state transition. △ Less

Submitted 18 April, 2024; originally announced April 2024.

Comments: 24 pages, 6 figures. Submitted to ApJ. Comments welcome!

The origin of the coherent radio flash potentially associated with GRB 201006A

Authors: Nikhil Sarin, Teagan A. Clarke, Spencer J. Magnall, Paul D. Lasky, Brian D. Metzger, Edo Berger, Navin Sridhar

Abstract: Rowlinson et al. 2023 recently claimed the detection of a coherent radio flash 76.6 minutes after a short gamma-ray burst. They proposed that the radio emission may be associated with a long-lived neutron star engine. We show through theoretical and observational arguments that the coherent radio emission, if real and indeed associated with GRB 201006A and at the estimated redshift, is unlikely to… ▽ More Rowlinson et al. 2023 recently claimed the detection of a coherent radio flash 76.6 minutes after a short gamma-ray burst. They proposed that the radio emission may be associated with a long-lived neutron star engine. We show through theoretical and observational arguments that the coherent radio emission, if real and indeed associated with GRB 201006A and at the estimated redshift, is unlikely to be due to the collapse of the neutron star, ruling out a blitzar-like mechanism. Instead, we show if a long-lived engine was created, it must have been stable with the radio emission likely linked to the intrinsic magnetar activity. However, we find that the optical upper limits require fine-tuning to be consistent with a magnetar-driven kilonova: we show that neutron-star engines that do satisfy the optical constraints would have produced a bright kilonova afterglow that should already be observable by the VLA or MeerKAT (for ambient densities typical for short GRBs). Given the optical limits and the current lack of a kilonova afterglow, we instead posit that no neutron star survived the merger, and the coherent radio emission was produced far from a black hole central engine via mechanisms such as synchrotron maser or magnetic reconnection in the jet -- a scenario consistent with all observations. We encourage future radio follow-up to probe the engine of this exciting event and continued prompt radio follow-up of short GRBs. △ Less

Submitted 11 April, 2024; originally announced April 2024.

Comments: Submitted. 10 pages, 4 figures

Birefringence tests of gravity with multi-messenger binaries

Authors: Macarena Lagos, Leah Jenks, Maximiliano Isi, Kenta Hotokezaka, Brian D. Metzger, Eric Burns, Will M. Farr, Scott Perkins, Kaze W. K. Wong, Nicolas Yunes

Abstract: Extensions to General Relativity (GR) allow the polarization of gravitational waves (GW) from astrophysical sources to suffer from amplitude and velocity birefringence, which respectively induce changes in the ellipticity and orientation of the polarization tensor. We introduce a multi-messenger approach to test this polarization behavior of GWs during their cosmological propagation using binary s… ▽ More Extensions to General Relativity (GR) allow the polarization of gravitational waves (GW) from astrophysical sources to suffer from amplitude and velocity birefringence, which respectively induce changes in the ellipticity and orientation of the polarization tensor. We introduce a multi-messenger approach to test this polarization behavior of GWs during their cosmological propagation using binary sources, for which the initial polarization is determined by the inclination and orientation angles of the orbital angular momentum vector with respect to the line of sight. In particular, we use spatially-resolved radio imaging of the jet from a binary neutron star (BNS) merger to constrain the orientation angle and hence the emitted polarization orientation of the GW signal at the site of the merger, and compare to that observed on Earth by GW detectors. For GW170817 we constrain the deviation from GR due to amplitude birefringence to $κ_A = -0.12^{+0.60}_{-0.61}$, while the velocity birefringence parameter $κ_V$ remains unconstrained. The inability to constrain $κ_V$ is due to the fact that Virgo did not detect GW170817, and measurements of the polarization orientation require information from a combination of multiple detectors with different alignments. For this reason, we also mock future BNS mergers with resolved afterglow proper motion and project that $κ_V$ could be constrained to a precision of $5\,$rad (corresponding to an angular shift of the GW polarization of $δφ_V\approx 0.2\,$rad for a BNS at $100\,$Mpc) by a future network of third-generation ground-based GW detectors such as Cosmic Explorer and the radio High Sensitivity Array. Crucially, this velocity birefringence effect cannot be constrained with dark binary mergers as it requires polarization information at the emission time, which can be provided only by electromagnetic emission. △ Less

Submitted 7 February, 2024; originally announced February 2024.

Long Plateau Doth So: How Internal Heating Sources Affect Hydrogen-Rich Supernova Light Curves

Authors: Tatsuya Matsumoto, Brian D. Metzger, Jared A. Goldberg

Abstract: Some hydrogen-rich core-collapse supernovae (type IIP SNe) exhibit evidence for a sustained energy source powering their light curves, resulting in a brighter and/or longer-lasting hydrogen-recombination plateau phase. We present a semi-analytic SNIIP light curve model that accounts for the effects of an arbitrary internal heating source, considering as special cases $^{56}$Ni/$^{56}$Co decay, a c… ▽ More Some hydrogen-rich core-collapse supernovae (type IIP SNe) exhibit evidence for a sustained energy source powering their light curves, resulting in a brighter and/or longer-lasting hydrogen-recombination plateau phase. We present a semi-analytic SNIIP light curve model that accounts for the effects of an arbitrary internal heating source, considering as special cases $^{56}$Ni/$^{56}$Co decay, a central engine (millisecond magnetar or accreting compact object), and shock interaction with a dense circumstellar disk. While a sustained internal power source can boost the plateau luminosity commensurate with the magnitude of the power, the duration of the recombination plateau can typically be increased by at most a factor $\sim 2-3$ compared to the zero-heating case. For a given ejecta mass and initial kinetic energy, the longest plateau duration is achieved for a constant heating rate at the highest magnitude that does not appreciably accelerate the ejecta. This finding has implications for the minimum ejecta mass required to explain particularly long-lasting supernovae such as iPTF14hls, and for confidently identifying rare explosions of the most-massive hydrogen-rich (e.g. population III) stars. We present a number of analytic estimates which elucidate the key features of the detailed model. △ Less

Submitted 24 January, 2024; originally announced January 2024.

Comments: 27 pages, 15 figures, 1 table, submitted to ApJ, comments welcome!!!

The Effects of $r$-Process Enrichment in Hydrogen-Rich Supernovae

Authors: Anirudh Patel, Jared A. Goldberg, Mathieu Renzo, Brian D. Metzger

Abstract: Core-collapse supernovae are candidate sites for the rapid neutron capture process ($r$-process). We explore the effects of enrichment from $r$-process nuclei on the light-curves of hydrogen-rich supernovae (SNe IIP) and assess the detectability of these signatures. We modify the radiation transport code $\texttt{SNEC}$ to include the approximate effects of opacity and radioactive heating from… ▽ More Core-collapse supernovae are candidate sites for the rapid neutron capture process ($r$-process). We explore the effects of enrichment from $r$-process nuclei on the light-curves of hydrogen-rich supernovae (SNe IIP) and assess the detectability of these signatures. We modify the radiation transport code $\texttt{SNEC}$ to include the approximate effects of opacity and radioactive heating from $r$-process elements in the SN ejecta. We present models spanning a range of total $r$-process masses $M_{\rm r}$ and their assumed radial distribution within the ejecta, finding that $M_{\rm r} \gtrsim 10^{-2} M_\odot$ is sufficient to induce appreciable differences in their light-curves as compared to ordinary SNe IIP (without any $r$-process elements). The primary photometric signatures of $r$-process enrichment include a shortening of the plateau phase, coinciding with the hydrogen-recombination photosphere retreating to the $r$-process-enriched layers, and a steeper post-plateau decline associated with a reddening of the SN colors. We compare our $r$-process-enriched models to ordinary IIP models and observational data, showing that yields of $M_{\rm r} \gtrsim 10^{-2} M_\odot$ are potentially detectable across several of the metrics used by transient observers, provided that the $r$-process rich layers are mixed $\gtrsim$ halfway to the ejecta surface. This detectability threshold can roughly be reproduced analytically using a two-zone ("kilonova within a supernova") picture. Assuming that a small fraction of SNe produce a detectable $r$-process yield $M_{\rm r} \gtrsim 10^{-2}M_\odot$, and respecting constraints on the total Galactic production rate, we estimate that $\gtrsim 10^{3}-10^4$ SNe need be observed to find one $r$-enriched event, a feat that may become possible with the Vera Rubin Observatory. △ Less

Submitted 28 April, 2024; v1 submitted 23 January, 2024; originally announced January 2024.

Comments: 15 pages, 7 figures, 3 appendices

A Hubble Space Telescope Search for r-Process Nucleosynthesis in Gamma-ray Burst Supernovae

Authors: J. C. Rastinejad, W. Fong, A. J. Levan, N. R. Tanvir, C. D. Kilpatrick, A. S. Fruchter, S. Anand, K. Bhirombhakdi, S. Covino, J. P. U. Fynbo, G. Halevi, D. H. Hartmann, K. E. Heintz, L. Izzo, P. Jakobsson, G. P. Lamb, D. B. Malesani, A. Melandri, B. D. Metzger, B. Milvang-Jensen, E. Pian, G. Pugliese, A. Rossi, D. M. Siegel, P. Singh , et al. (1 additional authors not shown)

Abstract: The existence of a secondary (in addition to compact object mergers) source of heavy element ($r$-process) nucleosynthesis, the core-collapse of rapidly-rotating and highly-magnetized massive stars, has been suggested by both simulations and indirect observational evidence. Here, we probe a predicted signature of $r$-process enrichment, a late-time ($\gtrsim 40$ days post-burst) distinct red color… ▽ More The existence of a secondary (in addition to compact object mergers) source of heavy element ($r$-process) nucleosynthesis, the core-collapse of rapidly-rotating and highly-magnetized massive stars, has been suggested by both simulations and indirect observational evidence. Here, we probe a predicted signature of $r$-process enrichment, a late-time ($\gtrsim 40$ days post-burst) distinct red color, in observations of GRB-supernovae (GRB-SNe) which are linked to these massive star progenitors. We present optical to near-IR color measurements of four GRB-SNe at $z \lesssim 0.4$, extending out to $> 500$ days post-burst, obtained with the Hubble Space Telescope and large-aperture ground-based telescopes. Comparison of our observations to models indicates that GRBs 030329, 100316D and 130427A are consistent with both no enrichment and producing $0.01 - 0.15 M_{\odot}$ of $r$-process material if there is a low amount of mixing between the inner $r$-process ejecta and outer SN layers. GRB 190829A is not consistent with any models with $r$-process enrichment $\geq 0.01 M_{\odot}$. Taken together the sample of GRB-SNe indicates color diversity at late times. Our derived yields from GRB-SNe may be underestimated due to $r$-process material hidden in the SN ejecta (potentially due to low mixing fractions) or the limits of current models in measuring $r$-process mass. We conclude with recommendations for future search strategies to observe and probe the full distribution of $r$-process produced by GRB-SNe. △ Less

Submitted 9 April, 2024; v1 submitted 7 December, 2023; originally announced December 2023.

Comments: Resubmission after comments. Accepted to ApJ. 36 pages, 7 figures, 2 tables

Ultraviolet Quasi-periodic Eruptions from Star-Disk Collisions in Galactic Nuclei

Authors: Itai Linial, Brian D. Metzger

Abstract: ``Quasi-periodic eruptions'' (QPE) are recurrent nuclear transients with periods of several hours to almost a day, which thus far have been detected exclusively in the X-ray band. We have shown that many of the key properties of QPE flares (period, luminosity, duration, emission temperature, alternating long-short recurrence time behavior, source rates) are naturally reproduced by a scenario invol… ▽ More ``Quasi-periodic eruptions'' (QPE) are recurrent nuclear transients with periods of several hours to almost a day, which thus far have been detected exclusively in the X-ray band. We have shown that many of the key properties of QPE flares (period, luminosity, duration, emission temperature, alternating long-short recurrence time behavior, source rates) are naturally reproduced by a scenario involving twice-per-orbit collisions between a solar-type star on a mildly eccentric orbit, likely brought into the nucleus as an extreme mass-ratio inspiral (EMRI), and the gaseous accretion disk of a supermassive black hole (SMBH). The flare is generated by the hot shocked debris expanding outwards from either side of the disk midplane, akin to dual miniature supernovae. Here, we consider the conditions necessary for disk-star collisions to generate lower-temperature flares which peak in the ultraviolet (UV) instead of the X-ray band. We identify a region of parameter space at low SMBH mass $M_{\bullet} \sim 10^{5.5}M_{\odot}$ and QPE periods $P \gtrsim 10$ hr for which the predicted flares are sufficiently luminous $L_{\rm UV} \sim 10^{41}$ erg s$^{-1}$ to outshine the quiescent disk emission at these wavelengths. The prospects to discover such ``UV QPEs'' with future satellite missions such as ULTRASAT and UVEX depends on the prevalence of very low-mass SMBH and the occurrence rate of stellar EMRIs onto them. For gaseous disks produced by the tidal disruption of stars, we predict that X-ray QPEs will eventually shut off, only to later reappear as UV-QPEs as the accretion rate continues to drop. △ Less

Submitted 27 November, 2023; originally announced November 2023.

Comments: 9 pages, 3 figures. Comments welcome

Dynamics of baryon ejection in magnetar giant flares: implications for radio afterglows, r-process nucleosynthesis, and fast radio bursts

Authors: Jakub Cehula, Todd A. Thompson, Brian D. Metzger

Abstract: We explore the impact of a magnetar giant flare (GF) on the neutron star (NS) crust, and the associated baryon mass ejection. We consider that sudden magnetic energy dissipation creates a thin high-pressure shell above a portion of the NS surface, which drives a relativistic shockwave into the crust, heating a fraction of these layers sufficiently to become unbound along directions unconfined by t… ▽ More We explore the impact of a magnetar giant flare (GF) on the neutron star (NS) crust, and the associated baryon mass ejection. We consider that sudden magnetic energy dissipation creates a thin high-pressure shell above a portion of the NS surface, which drives a relativistic shockwave into the crust, heating a fraction of these layers sufficiently to become unbound along directions unconfined by the magnetic field. We explore this process using spherically-symmetric relativistic hydrodynamical simulations. For an initial shell pressure $P_{\rm GF}$ we find the total unbound ejecta mass roughly obeys the relation $M_{\rm{ej}}\sim4-9\times10^{24}\:\rm{g}\:(P_{\rm GF}/10^{30}\:\rm{ergs}\:\rm{cm}^{-3})^{1.43}$. For $P_{\rm{GF}}\sim10^{30}-10^{31}\:\rm{ergs}\:\rm{cm}^{-3}$ corresponding to the dissipation of a magnetic field of strength $\sim10^{15.5}-10^{16}\:\rm{G}$, we find $M_{\rm{ej}}\sim10^{25}-10^{26}\:\rm{g}$ with asymptotic velocities $v_{\rm{ej}}/c\sim0.3-0.6$ compatible with the ejecta properties inferred from the afterglow of the December 2004 GF from SGR 1806-20. Because the flare excavates crustal material to a depth characterized by an electron fraction $Y_e\approx0.40-0.46$, and is ejected with high entropy and rapid expansion timescale, the conditions are met for heavy element $r$-process nucleosynthesis via the alpha-rich freeze-out mechanism. Given an energetic GF rate of roughly once per century in the Milky Way, we find that magnetar GFs could be an appreciable heavy $r$-process source that tracks star formation. We predict that GFs are accompanied by short $\sim$minutes long, luminous $\sim10^{39}\:\rm{ergs}\:\rm{s}^{-1}$ optical transients powered by $r$-process decay ("nova brevis"), akin to scaled-down kilonovae. Our findings also have implications for the synchrotron nebulae surrounding some repeating fast radio burst sources. △ Less

Submitted 1 February, 2024; v1 submitted 9 November, 2023; originally announced November 2023.

Comments: 24 pages, 12 figures, 2 tables. v2: matches the accepted version in MNRAS (only minor edits were made). Videos showing our Fiducial simulation are available at https://dx.doi.org/10.5281/zenodo.10593900 and at https://www.youtube.com/playlist?list=PLtfS1xXqUmtC2A48QhRvGXxbeMK1vkNnC

TESS photometry of the nova eruption in V606 Vul: asymmetric photosphere and multiple ejections?

Authors: Kirill V. Sokolovsky, Elias Aydi, Konstantin Malanchev, Colin J. Burke, Koji Mukai, J. L. Sokoloski, Brian D. Metzger, Kirill E. Atapin, Aleksandre A. Belinski, Yu-Ching Chen, Laura Chomiuk, Pavol A. Dubovsky, Claude-Andre Faucher-Giguere, Rebekah A. Hounsell, Natalia P. Ikonnikova, Vsevolod Yu. Lander, Junyao Li, Justin D. Linford, Amy J. Mioduszewski, Isabella Molina, Ulisse Munari, Sergey A. Potanin, Robert M. Quimby, Michael P. Rupen, Simone Scaringi , et al. (46 additional authors not shown)

Abstract: Lightcurves of many classical novae deviate from the canonical "fast rise - smooth decline" pattern and display complex variability behavior. We present the first TESS-space-photometry-based investigation of this phenomenon. We use TESS Sector 41 full-frame images to extract a lightcurve of the slow Galactic nova V606 Vul that erupted nine days prior to the start of the TESS observations. The ligh… ▽ More Lightcurves of many classical novae deviate from the canonical "fast rise - smooth decline" pattern and display complex variability behavior. We present the first TESS-space-photometry-based investigation of this phenomenon. We use TESS Sector 41 full-frame images to extract a lightcurve of the slow Galactic nova V606 Vul that erupted nine days prior to the start of the TESS observations. The lightcurve covers the first of two major peaks of V606 Vul that was reached 19 days after the start of the eruption. The nova reached its brightest visual magnitude V=9.9 in its second peak 64 days after the eruption onset, following the completion of Sector 41 observations. To increase the confidence level of the extracted lightcurve, we performed the analysis using four different codes implementing the aperture photometry (Lightkurve, VaST) and image subtraction (TESSreduce, tequila_shots) and find good agreement between them. We performed ground-based photometric and spectroscopic monitoring to complement the TESS data. The TESS lightcurve reveals two features: periodic variations (0.12771d, 0.01mag average peak-to-peak amplitude) that disappeared when the source was within 1mag of peak optical brightness and a series of isolated mini-flares (with peak-to-peak amplitudes of up to 0.5mag) appearing at seemingly random times. We interpret the periodic variations as the result of azimuthal asymmetry of the photosphere engulfing the nova-hosting binary that was distorted by and rotating with the binary. Whereas we use spectra to associate the two major peaks in the nova lightcurve with distinct episodes of mass ejection, the origin of mini-flares remains elusive. △ Less

Submitted 8 November, 2023; originally announced November 2023.

Comments: 31 pages, 10 figures, submitted to ApJ; comments welcome

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

The effect of relativistic precession on light curves of tidal disruption events

Authors: Diego Calderón, Ondřej Pejcha, Brian D. Metzger, Paul C. Duffell

Abstract: The disruption of a star by the tidal forces of a spinning black hole causes the stellar stream to precess affecting the conditions for triggering the tidal disruption event (TDE). In this work, we study the effect that precession imprints on TDE light curves due to the interaction of the TDE wind and luminosity with the stream wrapped around the black hole. We perform two-dimensional radiation-hy… ▽ More The disruption of a star by the tidal forces of a spinning black hole causes the stellar stream to precess affecting the conditions for triggering the tidal disruption event (TDE). In this work, we study the effect that precession imprints on TDE light curves due to the interaction of the TDE wind and luminosity with the stream wrapped around the black hole. We perform two-dimensional radiation-hydrodynamic simulations using the moving-mesh hydrodynamic code JET with its radiation treatment module. We study the impact of black hole mass, accretion efficiency, and inclination between the orbital and spin planes. From our results, we identified two behaviours: $i)$ models with low-mass black holes ($M_\text{h}\sim10^6~\text{M}_{\odot}$), low inclination ($i\sim0$), and low accretion efficiency ($η\sim0.01$) show light curves with a short early peak caused by the interaction of the wind with the inner edge of the stream. The line of sight has little effect on the light curve, since the stream covers a small fraction of the solid angle due to the precession occurring in the orbital plane; $ii)$ models with high-mass black holes ($M_\text{h}\gtrsim10^7~\text{M}_{\odot}$), high inclination ($i\sim90^{\circ}$), and high accretion efficiency ($η\sim0.1$) produce light curves with luminosity peaks that can be delayed by up to 50-100 d depending on the line of sight due to presence of the precessed stream blocking the radiation in the early phase of the event. Our results show that black hole spin and misalignment do not imprint recognisable features on the light curves but rather can add complications to their analysis. △ Less

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

Comments: 19 pages (+1 in Appendix), 13 figures (+1 in Appendix), 1 table. Accepted for publication in MNRAS

Revisiting the classics: On the evolutionary origin of the "Fe II" and "He/N" spectral classes of novae

Authors: E. Aydi, L. Chomiuk, J. Strader, K. V. Sokolovsky, R. E. Williams, D. A. H. Buckley, A. Ederoclite, L. Izzo, R. Kyer, J. D. Linford, A. Kniazev, B. D. Metzger, J. Mikolajewska, P. Molaro, I. Mollina, K. Mukai, U. Munari, M. Orio, T. Panurach, B. J. Shappee, K. J. Shen, J. L. Sokoloski, R. Urquhart, F. M. Walter

Abstract: The optical spectra of novae are characterized by emission lines from the hydrogen Balmer series and either Fe II or He/N, leading to their traditional classification into two spectral classes: "Fe II" and "He/N". For decades, the origins of these spectral features were discussed in the literature in the contexts of different bodies of gas or changes in the opacity of the ejecta, particularly asso… ▽ More The optical spectra of novae are characterized by emission lines from the hydrogen Balmer series and either Fe II or He/N, leading to their traditional classification into two spectral classes: "Fe II" and "He/N". For decades, the origins of these spectral features were discussed in the literature in the contexts of different bodies of gas or changes in the opacity of the ejecta, particularly associated with studies by R. E. Williams and S. N. Shore. Here, we revisit these major studies with dedicated, modern data sets, covering the evolution of several novae from early rise to peak all the way to the nebular phase. Our data confirm previous suggestions in the literature that the "Fe II" and "He/N" spectral classes are phases in the spectroscopic evolution of novae driven primarily by changes in the opacity, ionization, and density of the ejecta, and most if not all novae go through at least three spectroscopic phases as their eruptions evolve: an early He/N (phase 1; observed during the early rise to visible peak and characterized by P Cygni lines of He I, N II, and N III), then an Fe II (phase 2; observed near visible peak and characterized by P Cygni lines of Fe II and O I), and then a later He/N (phase 3; observed during the decline and characterized by emission lines of He I. He II, N II, and N III), before entering the nebular phase. This spectral evolution seems to be ubiquitous across novae, regardless of their speed class; however the duration of each of these phase differs based on the speed class of the nova. △ Less

Submitted 27 October, 2023; v1 submitted 13 September, 2023; originally announced September 2023.

Comments: 21 pages, 14 figures, 11 tables, Accepted in MNRAS

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

A Unified Theory of Jetted Tidal Disruption Events: From Promptly Escaping Relativistic to Delayed Transrelativistic Jets

Authors: Odelia Teboul, Brian D. Metzger

Abstract: Only a tiny fraction ~ 1% of stellar tidal disruption events (TDE) generate powerful relativistic jets evidenced by luminous hard X-ray and radio emissions. We propose that a key property responsible for both this surprisingly low rate and a variety of other observations is the typically large misalignment ψ between the orbital plane of the star and the spin axis of the supermassive black hole (SM… ▽ More Only a tiny fraction ~ 1% of stellar tidal disruption events (TDE) generate powerful relativistic jets evidenced by luminous hard X-ray and radio emissions. We propose that a key property responsible for both this surprisingly low rate and a variety of other observations is the typically large misalignment ψ between the orbital plane of the star and the spin axis of the supermassive black hole (SMBH). Such misaligned disk/jet systems undergo Lense-Thirring precession together about the SMBH spin axis. We find that TDE disks precess sufficiently rapidly that winds from the accretion disk will encase the system on large scales in a quasi-spherical outflow. We derive the critical jet efficiency η > ηcrit for both aligned and misaligned precessing jets to successfully escape from the disk-wind ejecta. As ηcrit is higher for precessing jets, less powerful jets only escape after alignment with the SMBH spin. Alignment can occur through magneto-spin or hydrodynamic mechanisms, which we estimate occur on typical timescales of weeks and years, respectively. The dominant mechanism depends on η and the orbital penetration factor \b{eta}. Hence depending only on intrinsic parameters of the event {ψ,η,\b{eta}}, we propose that each TDE jet can either escape prior to alignment, thus exhibiting erratic X-ray light curve and two-component radio afterglow (e.g., Swift J1644+57) or escape after alignment. Relatively rapid magneto-spin alignments produce relativistic jets exhibiting X-ray power-law decay and bright afterglows (e.g., AT2022cmc), while long hydrodynamic alignments give rise to late jet escape and delayed radio flares (e.g., AT2018hyz). △ Less

Submitted 9 August, 2023; originally announced August 2023.

Comments: 17 pages, 7 figures, 1 table

Tidal Disruption Events through the Lens of the Cooling Envelope Model

Authors: Nikhil Sarin, Brian D. Metzger

Abstract: The cooling envelope model for tidal disruption events (TDE) postulates that while the stellar debris streams rapidly dissipate their bulk kinetic energy (``circularize"), this does not necessarily imply rapid feeding of the supermassive black hole (SMBH). The bound material instead forms a large pressure-supported envelope which powers optical/UV emission as it undergoes gradual Kelvin-Helmholtz… ▽ More The cooling envelope model for tidal disruption events (TDE) postulates that while the stellar debris streams rapidly dissipate their bulk kinetic energy (``circularize"), this does not necessarily imply rapid feeding of the supermassive black hole (SMBH). The bound material instead forms a large pressure-supported envelope which powers optical/UV emission as it undergoes gradual Kelvin-Helmholtz contraction. We present results interpreting a sample of 15 optical TDE within the cooling envelope model in order to constrain the SMBH mass $M_{\rm BH}$, stellar mass $M_{\star}$, and orbital penetration factor $β$. The distributions of inferred properties from our sample broadly follow the theoretical expectations of loss-cone analysis assuming a standard stellar initial mass function. However, we find a deficit of events with $M_{\rm BH} \lesssim 5\times 10^{5}M_{\odot}$ and $M_{\star} \lesssim 0.5M_{\odot}$, which could result in part from the reduced detectability of TDEs with these properties. Our model fits also illustrate the predicted long delay between the optical light curve peak and when the SMBH accretion rate reaches its maximum. The latter occurs only once the envelope contracts to the circularization radius on a timescale of months to years, consistent with delayed-rising X-ray and non-thermal radio flares seen in a growing number of TDE. △ Less

Submitted 19 January, 2024; v1 submitted 27 July, 2023; originally announced July 2023.

Comments: Published in ApJL. 19 pages, 9 figures. Model and inference workflow implemented in Redback at https://github.com/nikhil-sarin/redback

An X-ray Census of Fast Radio Burst Host Galaxies: Constraints on AGN and X-ray Counterparts

Authors: T. Eftekhari, W. Fong, A. C. Gordon, N. Sridhar, C. D. Kilpatrick, S. Bhandari, A. T. Deller, Y. Dong, A. Rouco Escorial, K. E. Heintz, J. Leja, B. Margalit, B. D. Metzger, A. B. Pearlman, J. X. Prochaska, S. D. Ryder, P. Scholz, R. M. Shannon, N. Tejos

Abstract: We present the first X-ray census of fast radio burst (FRB) host galaxies to conduct the deepest search for AGN and X-ray counterparts to date. Our sample includes seven well-localized FRBs with unambiguous host associations and existing deep Chandra observations, including two events for which we present new observations. We find evidence for AGN in two FRB host galaxies based on the presence of… ▽ More We present the first X-ray census of fast radio burst (FRB) host galaxies to conduct the deepest search for AGN and X-ray counterparts to date. Our sample includes seven well-localized FRBs with unambiguous host associations and existing deep Chandra observations, including two events for which we present new observations. We find evidence for AGN in two FRB host galaxies based on the presence of X-ray emission coincident with their centers, including the detection of a luminous ($L_X\approx\,5\times\,10^{42}\,\rm\,erg\,s^{-1}$) X-ray source at the nucleus of FRB20190608B's host, for which we infer an SMBH mass of $\rm{M_{BH}\sim\,10^{8}\,M_{\odot}}$ and an Eddington ratio $\rm{L_{bol}/L_{Edd}\approx\,0.02}$, characteristic of geometrically thin disks in Seyfert galaxies. We also report nebular emission line fluxes for 24 highly secure FRB hosts (including 10 hosts for the first time), and assess their placement on a BPT diagram, finding that FRB hosts trace the underlying galaxy population. We further find that the hosts of repeating FRBs are not confined to the star-forming locus, contrary to previous findings. Finally, we place constraints on associated X-ray counterparts to FRBs in the context of ultraluminous X-ray sources (ULXs), and find that existing X-ray limits for FRBs rule out ULXs brighter than $L_X\gtrsim\,10^{40}\,\rm\,erg\,s^{-1}$. Leveraging the CHIME/FRB catalog and existing ULX catalogs, we search for spatially coincident ULX-FRB pairs. We identify a total of 28 ULXs spatially coincident with the localization regions for 17 FRBs, but find that the DM-inferred redshifts for the FRBs are inconsistent with the ULX redshifts, disfavoring an association between these specific ULX-FRB pairs. △ Less

Submitted 27 November, 2023; v1 submitted 7 July, 2023; originally announced July 2023.

Comments: Accepted to ApJ

JWST detection of heavy neutron capture elements in a compact object merger

Authors: A. Levan, B. P. Gompertz, O. S. Salafia, M. Bulla, E. Burns, K. Hotokezaka, L. Izzo, G. P. Lamb, D. B. Malesani, S. R. Oates, M. E. Ravasio, A. Rouco Escorial, B. Schneider, N. Sarin, S. Schulze, N. R. Tanvir, K. Ackley, G. Anderson, G. B. Brammer, L. Christensen, V. S. Dhillon, P. A. Evans, M. Fausnaugh, W. -F. Fong, A. S. Fruchter , et al. (58 additional authors not shown)

Abstract: The mergers of binary compact objects such as neutron stars and black holes are of central interest to several areas of astrophysics, including as the progenitors of gamma-ray bursts (GRBs), sources of high-frequency gravitational waves and likely production sites for heavy element nucleosynthesis via rapid neutron capture (the r-process). These heavy elements include some of great geophysical, bi… ▽ More The mergers of binary compact objects such as neutron stars and black holes are of central interest to several areas of astrophysics, including as the progenitors of gamma-ray bursts (GRBs), sources of high-frequency gravitational waves and likely production sites for heavy element nucleosynthesis via rapid neutron capture (the r-process). These heavy elements include some of great geophysical, biological and cultural importance, such as thorium, iodine and gold. Here we present observations of the exceptionally bright gamma-ray burst GRB 230307A. We show that GRB 230307A belongs to the class of long-duration gamma-ray bursts associated with compact object mergers, and contains a kilonova similar to AT2017gfo, associated with the gravitational-wave merger GW170817. We obtained James Webb Space Telescope mid-infrared (mid-IR) imaging and spectroscopy 29 and 61 days after the burst. The spectroscopy shows an emission line at 2.15 microns which we interpret as tellurium (atomic mass A=130), and a very red source, emitting most of its light in the mid-IR due to the production of lanthanides. These observations demonstrate that nucleosynthesis in GRBs can create r-process elements across a broad atomic mass range and play a central role in heavy element nucleosynthesis across the Universe. △ Less

Submitted 5 July, 2023; originally announced July 2023.

Comments: Submitted. Comments welcome! Nature (2023)

Large-scale Evolution of Seconds-long Relativistic Jets from Black Hole-Neutron Star Mergers

Authors: Ore Gottlieb, Danat Issa, Jonatan Jacquemin-Ide, Matthew Liska, Francois Foucart, Alexander Tchekhovskoy, Brian D. Metzger, Eliot Quataert, Rosalba Perna, Daniel Kasen, Matthew D. Duez, Lawrence E. Kidder, Harald P. Pfeiffer, Mark A. Scheel

Abstract: We present the first numerical simulations that track the evolution of a black hole-neutron star (BH-NS) merger from pre-merger to $r\gtrsim10^{11}\,{\rm cm}$. The disk that forms after a merger of mass ratio $q=2$ ejects massive disk winds ($3-5\times10^{-2}\,M_{\odot}$). We introduce various post-merger magnetic configurations and find that initial poloidal fields lead to jet launching shortly a… ▽ More We present the first numerical simulations that track the evolution of a black hole-neutron star (BH-NS) merger from pre-merger to $r\gtrsim10^{11}\,{\rm cm}$. The disk that forms after a merger of mass ratio $q=2$ ejects massive disk winds ($3-5\times10^{-2}\,M_{\odot}$). We introduce various post-merger magnetic configurations and find that initial poloidal fields lead to jet launching shortly after the merger. The jet maintains a constant power due to the constancy of the large-scale BH magnetic flux until the disk becomes magnetically arrested (MAD), where the jet power falls off as $L_j\sim t^{-2}$. All jets inevitably exhibit either excessive luminosity due to rapid MAD activation when the accretion rate is high or excessive duration due to delayed MAD activation compared to typical short gamma-ray bursts (sGRBs). This provides a natural explanation for long sGRBs such as GRB 211211A but also raises a fundamental challenge to our understanding of jet formation in binary mergers. One possible implication is the necessity of higher binary mass ratios or moderate BH spins to launch typical sGRB jets. For post-merger disks with a toroidal magnetic field, dynamo processes delay jet launching such that the jets break out of the disk winds after several seconds. We show for the first time that sGRB jets with initial magnetization $σ_0>100$ retain significant magnetization ($σ\gg1$) at $r>10^{10}\,{\rm cm}$, emphasizing the importance of magnetic processes in the prompt emission. The jet-wind interaction leads to a power-law angular energy distribution by inflating an energetic cocoon whose emission is studied in a companion paper. △ Less

Submitted 18 August, 2023; v1 submitted 26 June, 2023; originally announced June 2023.

Comments: For movies of the simulations, see https://oregottlieb.com/bhns.html

Hours-long Near-UV/Optical Emission from Mildly Relativistic Outflows in Black Hole-Neutron Star Mergers

Authors: Ore Gottlieb, Danat Issa, Jonatan Jacquemin-Ide, Matthew Liska, Alexander Tchekhovskoy, Francois Foucart, Daniel Kasen, Rosalba Perna, Eliot Quataert, Brian D. Metzger

Abstract: The ongoing LIGO-Virgo-KAGRA observing run O4 provides an opportunity to discover new multi-messenger events, including binary neutron star (BNS) mergers such as GW170817, and the highly anticipated first detection of a multi-messenger black hole-neutron star (BH-NS) merger. While BNS mergers were predicted to exhibit early optical emission from mildly relativistic outflows, it has remained uncert… ▽ More The ongoing LIGO-Virgo-KAGRA observing run O4 provides an opportunity to discover new multi-messenger events, including binary neutron star (BNS) mergers such as GW170817, and the highly anticipated first detection of a multi-messenger black hole-neutron star (BH-NS) merger. While BNS mergers were predicted to exhibit early optical emission from mildly relativistic outflows, it has remained uncertain whether the BH-NS merger ejecta provides the conditions for similar signals to emerge. We present the first modeling of early near-ultraviolet/optical emission from mildly relativistic outflows in BH-NS mergers. Adopting optimal binary properties: a mass ratio of $q=2$ and a rapidly rotating BH, we utilize numerical relativity and general relativistic magnetohydrodynamic (GRMHD) simulations to follow the binary's evolution from pre-merger to homologous expansion. We use an M1 neutrino transport GRMHD simulation to self-consistently estimate the opacity distribution in the outflows and find a bright near-ultraviolet/optical signal that emerges due to jet-powered cocoon cooling emission, outshining the kilonova emission at early time. The signal peaks at an absolute magnitude of $\sim -15$ a few hours after the merger, longer than previous estimates, which did not consider the first principles-based jet launching. By late 2024, the Rubin Observatory will have the capability to track the entire signal evolution or detect its peak up to distances of $\gtrsim1$ Gpc. In 2026, ULTRASAT will conduct all-sky surveys within minutes, detecting some of these events within $\sim 200$ Mpc. The BH-NS mergers with higher mass ratios or lower BH spins would produce shorter and fainter signals. △ Less

Submitted 8 August, 2023; v1 submitted 26 June, 2023; originally announced June 2023.

Comments: For movies of the simulations, see https://oregottlieb.com/bhns.html

Long-Term Evolution of Massive-Star Post-Common Envelope Circumbinary Disks and the Environments of Fast Luminous Transients

Authors: Semih Tuna, Brian D. Metzger

Abstract: If the envelope of a massive star is not entirely removed during common envelope (CE) interaction with an orbiting compact (e.g., black hole [BH] or neutron star [NS]) companion, the residual bound material eventually cools, forming a centrifugally-supported disk around the binary containing the stripped He core. We present a time-dependent height-integrated model for the long-term evolution of po… ▽ More If the envelope of a massive star is not entirely removed during common envelope (CE) interaction with an orbiting compact (e.g., black hole [BH] or neutron star [NS]) companion, the residual bound material eventually cools, forming a centrifugally-supported disk around the binary containing the stripped He core. We present a time-dependent height-integrated model for the long-term evolution of post-CE circumbinary disks (CBD), accounting for mass and angular momentum exchange with the binary and irradiation heating by the He core and photoevaporation wind mass-loss. A large fraction of the CBD's mass is accreted prior to its outwards viscous spreading and wind-dispersal on a timescale ~10^{4}-10^{5} yr, driving significant changes in the binary separation, even for disks containing ~ 10% of the original envelope mass. Insofar that the CBD lifetime is comparable to the thermal (and, potentially, nuclear) timescale of the He core, over which a second mass-transfer episode onto the companion can occur, the presence of the CBD could impact the stability of this key phase. Disruption of the He core by the BH/NS would result in a jetted energetic explosion into the dense gaseous CBD (<~10^{15} cm) and its wind (>~ 10^{16} cm), consistent with the environments of luminous fast blue optical transients like AT2018cow. Evolved He cores which undergo core-collapse still embedded in their CBD could generate Type Ibn/Icn supernovae. Thousands of dusty wind-shrouded massive-star CBD may be detectable as extragalactic luminous infrared sources with the Roman Space Telescope; synchrotron radio nebulae powered by the CBD-fed BH/NS may accompany these systems. △ Less

Submitted 16 June, 2023; originally announced June 2023.

Comments: 26 pages, 11 figures

Three-Dimensional General-Relativistic Simulations of Neutrino-Driven Winds from Magnetized Proto-Neutron Stars

Authors: Dhruv K. Desai, Daniel M. Siegel, Brian D. Metzger

Abstract: Formed in the aftermath of a core-collapse supernova or neutron star merger, a hot proto-neutron star (PNS) launches an outflow driven by neutrino heating lasting for up to tens of seconds. Though such winds are considered potential sites for the nucleosynthesis of heavy elements via the rapid neutron capture process ($r$-process), previous work has shown that unmagnetized PNS winds fail to achiev… ▽ More Formed in the aftermath of a core-collapse supernova or neutron star merger, a hot proto-neutron star (PNS) launches an outflow driven by neutrino heating lasting for up to tens of seconds. Though such winds are considered potential sites for the nucleosynthesis of heavy elements via the rapid neutron capture process ($r$-process), previous work has shown that unmagnetized PNS winds fail to achieve the necessary combination of high entropy and/or short dynamical timescale in the seed nucleus formation region. We present three-dimensional general-relativistic magnetohydrodynamical (GRMHD) simulations of PNS winds which include the effects of a dynamically strong ($B \gtrsim 10^{15}$ G) dipole magnetic field. After initializing the magnetic field, the wind quickly develops a helmet-streamer configuration, characterized by outflows along open polar magnetic field lines and a ``closed'' zone of trapped plasma at lower latitudes. Neutrino heating within the closed zone causes the thermal pressure of the trapped material to rise in time compared to the polar outflow regions, ultimately leading to the expulsion of this matter from the closed zone on a timescale of $\sim$60 ms, consistent with the predictions of \citet{Thompson03}. The high entropies of these transient ejecta are still growing at the end of our simulations and are sufficient to enable a successful 2nd-peak $r$-process in at least a modest $\gtrsim 1\%$ of the equatorial wind ejecta. △ Less

Submitted 6 June, 2023; originally announced June 2023.

Comments: 22 pages, 16 figures

Many Roads Lead to Lithium: Formation Pathways For Lithium-Rich Red Giants

Authors: Maryum Sayeed, Melissa K. Ness, Benjamin T. Montet, Matteo Cantiello, Andrew R. Casey, Sven Buder, Megan Bedell, Katelyn Breivik, Brian D. Metzger, Sarah L. Martell, Leah McGee-Gold

Abstract: Stellar models predict that lithium (Li) inside a star is destroyed during the first dredge-up phase, yet 1.2% of red giant stars are Li-rich. We aim to uncover possible origins of this population, by analysing 1155 Li-rich giants (A(Li) $\geq$ 1.5) in GALAH DR3. To expose peculiar traits of Li-rich stars, we construct a reference sample of Li-normal (doppelgänger) stars with matched evolutionary… ▽ More Stellar models predict that lithium (Li) inside a star is destroyed during the first dredge-up phase, yet 1.2% of red giant stars are Li-rich. We aim to uncover possible origins of this population, by analysing 1155 Li-rich giants (A(Li) $\geq$ 1.5) in GALAH DR3. To expose peculiar traits of Li-rich stars, we construct a reference sample of Li-normal (doppelgänger) stars with matched evolutionary state and fiducial supernova abundances. Comparing Li-rich and doppelgänger spectra reveals systematic differences in the H-$α$ and Ca-triplet line profiles associated with the velocity broadening measurement. We also find twice as many Li-rich stars appear to be fast rotators (2% with $v_\textrm{broad} \gtrsim 20$ km s$^{-1}$) compared to doppelgängers. On average, Li-rich stars have higher abundances than their doppelgängers, for a subset of elements, and Li-rich stars at the base of RGB have higher mean $s-$process abundances ($\geq 0.05$ dex for Ba, Y, Zr), relative to their doppelgängers. External mass-transfer from intermediate-mass AGB companions could explain this signature. Additional companion analysis excludes binaries with mass ratios $\gtrsim$ 0.5 at $\gtrsim$ 7 AU. We also discover that highly Ba-enriched stars are missing from the Li-rich population, possibly due to low-mass AGB companions which preclude Li-enrichment. Finally, we confirm a prevalence of Li-rich stars on the red clump that increases with lithium, which supports an evolutionary state mechanism for Li-enhancement. Multiple culprits, including binary spin-up and mass-transfer, are therefore likely mechanisms of Li-enrichment. △ Less

Submitted 5 June, 2023; originally announced June 2023.

Comments: 29 pages, 19 figures, 6 tables. Submitted to ApJ

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

Catching a nova X-ray/UV flash in the visible? Early spectroscopy of the extremely slow Nova Velorum 2022 (Gaia22alz)

Authors: E. Aydi, L. Chomiuk, J. Mikołajewska, J. Brink, B. D. Metzger, J. Strader, D. A. H. Buckley, E. J. Harvey, T. W. -S. Holoien, L. Izzo, A. Kawash, J. D. Linford, P. Molaro, B. Mollina, P. Mróz, K. Mukai, M. Orio, T. Panurach, P. Senchyna, B. J. Shappee, K. J. Shen, J. L. Sokoloski, K. V. Sokolovsky, R. Urquhart, R. E. Williams

Abstract: We present early spectral observations of the very slow Galactic nova Gaia22alz, over its gradual rise to peak brightness that lasted 180 days. During the first 50 days, when the nova was only 3--4 magnitudes above its normal brightness, the spectra showed narrow (FWHM $\approx$ 400 km s$^{-1}$) emission lines of H Balmer, He I, He II, and C IV, but no P Cygni absorption. A few weeks later, the hi… ▽ More We present early spectral observations of the very slow Galactic nova Gaia22alz, over its gradual rise to peak brightness that lasted 180 days. During the first 50 days, when the nova was only 3--4 magnitudes above its normal brightness, the spectra showed narrow (FWHM $\approx$ 400 km s$^{-1}$) emission lines of H Balmer, He I, He II, and C IV, but no P Cygni absorption. A few weeks later, the high-excitation He II and C IV lines disappeared, and P Cygni profiles of Balmer, He I, and eventually Fe II lines emerged, yielding a spectrum typical of classical novae before peak. We propose that the early spectra of Gaia22alz are produced in the white dwarf's envelope or accretion disk, reprocessing X-ray and ultraviolet emission from the white dwarf after a dramatic increase in the rate of thermonuclear reactions, during a phase known as the ``early X-ray/UV flash''. If true, this would be one of the rare times that the optical signature of the early X-ray/UV flash has been detected. While this phase might last only a few hours in other novae and thus be easily missed, it was possible to detect in Gaia22alz due to its very slow and gradual rise and thanks to the efficiency of new all-sky surveys in detecting transients on their rise. We also consider alternative scenarios that could explain the early spectral features of Gaia22alz and its unusually slow rise. △ Less

Submitted 9 April, 2023; originally announced April 2023.

Comments: 20 pages, 12 figures, 2 tables. Submitted to MNRAS

EMRI + TDE = QPE: Periodic X-ray Flares from Star-Disk Collisions in Galactic Nuclei

Authors: Itai Linial, Brian D. Metzger

Abstract: Roughly half of the quasi-periodic eruption (QPE) sources in galactic nuclei exhibit a remarkably regular alternating "long-short'' pattern of recurrence times between consecutive flares. We show that a main-sequence star (brought into the nucleus as an extreme mass-ratio inspiral; EMRI) which passes twice per orbit through the accretion disk of the supermassive black-hole (SMBH) on a mildly eccen… ▽ More Roughly half of the quasi-periodic eruption (QPE) sources in galactic nuclei exhibit a remarkably regular alternating "long-short'' pattern of recurrence times between consecutive flares. We show that a main-sequence star (brought into the nucleus as an extreme mass-ratio inspiral; EMRI) which passes twice per orbit through the accretion disk of the supermassive black-hole (SMBH) on a mildly eccentric inclined orbit, each time shocking and ejecting optically-thick gas clouds above and below the midplane, naturally reproduces observed properties of QPE flares. Inefficient photon production in the ejecta renders the QPE emission much harder than the blackbody temperature, enabling the flares to stick out from the softer quiescent disk spectrum. Destruction of the star via mass ablation limits the QPE lifetime to decades, precluding a long-lived AGN as the gaseous disk. By contrast, a tidal disruption event (TDE) naturally provides a transient gaseous disk on the requisite radial scale, with a rate exceeding the EMRI inward migration rate, suggesting that many TDEs should host a QPE. This picture is consistent with the X-ray TDE observed several years prior to the QPE appearance from GSN 069. Remarkably, a second TDE-like flare was observed from this event, starting immediately after detectable QPE activity ceased; this event could plausibly result from the (partial or complete) destruction of the QPE-generating star triggered by runaway mass-loss, though other explanations cannot be excluded. Our model can also be applied to black hole-disk collisions, such as those invoked in the context of the candidate SMBH binary OJ 287. △ Less

Submitted 1 September, 2023; v1 submitted 28 March, 2023; originally announced March 2023.

Comments: 21 pages, 4 figures, 2 tables, Accepted to ApJ

A long-duration gamma-ray burst of dynamical origin from the nucleus of an ancient galaxy

Authors: Andrew J. Levan, Daniele B. Malesani, Benjamin P. Gompertz, Anya E. Nugent, Matt Nicholl, Samantha Oates, Daniel A. Perley, Jillian Rastinejad, Brian D. Metzger, Steve Schulze, Elizabeth R. Stanway, Anne Inkenhaag, Tayyaba Zafar, J. Feliciano Agui Fernandez, Ashley Chrimes, Kornpob Bhirombhakdi, Antonio de Ugarte Postigo, Wen-fai Fong, Andrew S. Fruchter, Giacomo Fragione, Johan P. U. Fynbo, Nicola Gaspari, Kasper E. Heintz, Jens Hjorth, Pall Jakobsson , et al. (7 additional authors not shown)

Abstract: The majority of long duration ($>2$ s) gamma-ray bursts (GRBs) are believed to arise from the collapse of massive stars \cite{Hjorth+03}, with a small proportion created from the merger of compact objects. Most of these systems are likely formed via standard stellar evolution pathways. However, it has long been thought that a fraction of GRBs may instead be an outcome of dynamical interactions in… ▽ More The majority of long duration ($>2$ s) gamma-ray bursts (GRBs) are believed to arise from the collapse of massive stars \cite{Hjorth+03}, with a small proportion created from the merger of compact objects. Most of these systems are likely formed via standard stellar evolution pathways. However, it has long been thought that a fraction of GRBs may instead be an outcome of dynamical interactions in dense environments, channels which could also contribute significantly to the samples of compact object mergers detected as gravitational wave sources. Here we report the case of GRB 191019A, a long GRB (T_90 = 64.4 +/- 4.5 s) which we pinpoint close (<100 pc projected) to the nucleus of an ancient (>1~Gyr old) host galaxy at z=0.248. The lack of evidence for star formation and deep limits on any supernova emission make a massive star origin difficult to reconcile with observations, while the timescales of the emission rule out a direct interaction with the supermassive black hole in the nucleus of the galaxy, We suggest that the most likely route for progenitor formation is via dynamical interactions in the dense nucleus of the host, consistent with the centres of such galaxies exhibiting interaction rates up to two orders of magnitude larger than typical field galaxies. The burst properties could naturally be explained via compact object mergers involving white dwarfs (WD), neutron stars (NS) or black holes (BH). These may form dynamically in dense stellar clusters, or originate in a gaseous disc around the supermassive black hole. Future electromagnetic and gravitational-wave observations in tandem thus offer a route to probe the dynamical fraction and the details of dynamical interactions in galactic nuclei and other high density stellar systems. △ Less

Submitted 22 March, 2023; originally announced March 2023.

Comments: Accepted to Nature Astronomy. This is the submitted version and will differ from the published version due to modifications in the refereeing process

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

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

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

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

Comments: 37 pages, 8 figures, 4 tables

Journal ref: ApJ, 962, 68 (2024)

The first JWST spectrum of a GRB afterglow: No bright supernova in observations of the brightest GRB of all time, GRB 221009A

Authors: A. J. Levan, G. P. Lamb, B. Schneider, J. Hjorth, T. Zafar, A. de Ugarte Postigo, B. Sargent, S. E. Mullally, L. Izzo, P. D'Avanzo, E. Burns, J. F. Agüí Fernández, T. Barclay, M. G. Bernardini, K. Bhirombhakdi, M. Bremer, R. Brivio, S. Campana, A. A. Chrimes, V. D'Elia, M. Della Valle, M. De Pasquale, M. Ferro, W. Fong, A. S. Fruchter , et al. (35 additional authors not shown)

Abstract: We present JWST and Hubble Space Telescope (HST) observations of the afterglow of GRB 221009A, the brightest gamma-ray burst (GRB) ever observed. This includes the first mid-IR spectra of any GRB, obtained with JWST/NIRSPEC (0.6-5.5 micron) and MIRI (5-12 micron), 12 days after the burst. Assuming that the intrinsic spectral slope is a single power-law, with $F_ν \propto ν^{-β}$, we obtain… ▽ More We present JWST and Hubble Space Telescope (HST) observations of the afterglow of GRB 221009A, the brightest gamma-ray burst (GRB) ever observed. This includes the first mid-IR spectra of any GRB, obtained with JWST/NIRSPEC (0.6-5.5 micron) and MIRI (5-12 micron), 12 days after the burst. Assuming that the intrinsic spectral slope is a single power-law, with $F_ν \propto ν^{-β}$, we obtain $β\approx 0.35$, modified by substantial dust extinction with $A_V = 4.9$. This suggests extinction above the notional Galactic value, possibly due to patchy extinction within the Milky Way or dust in the GRB host galaxy. It further implies that the X-ray and optical/IR regimes are not on the same segment of the synchrotron spectrum of the afterglow. If the cooling break lies between the X-ray and optical/IR, then the temporal decay rates would only match a post jet-break model, with electron index $p<2$, and with the jet expanding into a uniform ISM medium. The shape of the JWST spectrum is near-identical in the optical/nIR to X-shooter spectroscopy obtained at 0.5 days and to later time observations with HST. The lack of spectral evolution suggests that any accompanying supernova (SN) is either substantially fainter or bluer than SN 1998bw, the proto-type GRB-SN. Our HST observations also reveal a disc-like host galaxy, viewed close to edge-on, that further complicates the isolation of any supernova component. The host galaxy appears rather typical amongst long-GRB hosts and suggests that the extreme properties of GRB 221009A are not directly tied to its galaxy-scale environment. △ Less

Submitted 22 March, 2023; v1 submitted 15 February, 2023; originally announced February 2023.

Comments: Accepted for publication to the Astrophysical Journal Letters for the GRB 221009A Special Issue. The results of this paper are under press embargo until March 28, 18 UT. 19 pages, 8 figures, 2 tables

VERITAS and Fermi-LAT constraints on the Gamma-ray Emission from Superluminous Supernovae SN2015bn and SN2017egm

Authors: A. Acharyya, C. B. Adams, P. Bangale, W. Benbow, J. H. Buckley, M. Capasso, V. V. Dwarkadas, M. Errando, A. Falcone, Q. Feng, J. P. Finley, G. M. Foote, L. Fortson, A. Furniss, G. Gallagher, A. Gent, W. F Hanlon, O. Hervet, J. Holder, T. B. Humensky, W. Jin, P. Kaaret, M. Kertzman, M. Kherlakian, D. Kieda , et al. (34 additional authors not shown)

Abstract: Superluminous supernovae (SLSNe) are a rare class of stellar explosions with luminosities ~10-100 times greater than ordinary core-collapse supernovae. One popular model to explain the enhanced optical output of hydrogen-poor (Type I) SLSNe invokes energy injection from a rapidly spinning magnetar. A prediction in this case is that high-energy gamma rays, generated in the wind nebula of the magnet… ▽ More Superluminous supernovae (SLSNe) are a rare class of stellar explosions with luminosities ~10-100 times greater than ordinary core-collapse supernovae. One popular model to explain the enhanced optical output of hydrogen-poor (Type I) SLSNe invokes energy injection from a rapidly spinning magnetar. A prediction in this case is that high-energy gamma rays, generated in the wind nebula of the magnetar, could escape through the expanding supernova ejecta at late times (months or more after optical peak). This paper presents a search for gamma-ray emission in the broad energy band from 100 MeV to 30 TeV from two Type I SLSNe, SN2015bn, and SN2017egm, using observations from Fermi-LAT and VERITAS. Although no gamma-ray emission was detected from either source, the derived upper limits approach the putative magnetar's spin-down luminosity. Prospects are explored for detecting very-high-energy (VHE; 100 GeV - 100 TeV) emission from SLSNe-I with existing and planned facilities such as VERITAS and CTA. △ Less

Submitted 13 February, 2023; originally announced February 2023.

Comments: 20 pages, 7 figures, 2 tables

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

The multi-wavelength view of shocks in the fastest nova V1674 Her

Authors: K. V. Sokolovsky, T. J. Johnson, S. Buson, P. Jean, C. C. Cheung, K. Mukai, L. Chomiuk, E. Aydi, B. Molina, A. Kawash, J. D. Linford, A. J. Mioduszewski, M. P. Rupen, J. L. Sokoloski, M. N. Williams, E. Steinberg, I. Vurm, B. D. Metzger, K. L. Page, M. Orio, R. M. Quimby, A. W. Shafter, H. Corbett, S. Bolzoni, J. DeYoung , et al. (19 additional authors not shown)

Abstract: Classical novae are shock-powered multi-wavelength transients triggered by a thermonuclear runaway on an accreting white dwarf. V1674 Her is the fastest nova ever recorded (time to declined by two magnitudes is t_2=1.1 d) that challenges our understanding of shock formation in novae. We investigate the physical mechanisms behind nova emission from GeV gamma-rays to cm-band radio using coordinated… ▽ More Classical novae are shock-powered multi-wavelength transients triggered by a thermonuclear runaway on an accreting white dwarf. V1674 Her is the fastest nova ever recorded (time to declined by two magnitudes is t_2=1.1 d) that challenges our understanding of shock formation in novae. We investigate the physical mechanisms behind nova emission from GeV gamma-rays to cm-band radio using coordinated Fermi-LAT, NuSTAR, Swift and VLA observations supported by optical photometry. Fermi-LAT detected short-lived (18 h) 0.1-100 GeV emission from V1674 Her that appeared 6 h after the eruption began; this was at a level of (1.6 +/- 0.4)x10^-6 photons cm^-2 s^-1. Eleven days later, simultaneous NuSTAR and Swift X-ray observations revealed optically thin thermal plasma shock-heated to kT_shock = 4 keV. The lack of a detectable 6.7 keV Fe K_alpha emission suggests super-solar CNO abundances. The radio emission from V1674 Her was consistent with thermal emission at early times and synchrotron at late times. The radio spectrum steeply rising with frequency may be a result of either free-free absorption of synchrotron and thermal emission by unshocked outer regions of the nova shell or the Razin-Tsytovich effect attenuating synchrotron emission in dense plasma. The development of the shock inside the ejecta is unaffected by the extraordinarily rapid evolution and the intermediate polar host of this nova. △ Less

Submitted 21 March, 2023; v1 submitted 6 February, 2023; originally announced February 2023.

Comments: 20 pages, 9 figures, 3 tables. Accepted to MNRAS

Synchrotron Afterglow Model for AT 2022cmc: Jetted Tidal Disruption Event or Engine-Powered Supernova?

Authors: Tatsuya Matsumoto, Brian D. Metzger

Abstract: AT 2022cmc is a luminous optical transient ($νL_ν \gtrsim 10^{45}$ erg s$^{-1}$) accompanied by decaying non-thermal X-rays (peak duration $t_{\rm X} \lesssim$ days and isotropic energy $E_{\rm X,iso} \gtrsim 10^{53}$ erg) and a long-lived radio/mm synchrotron afterglow, which has been interpreted as a jetted tidal disruption event (TDE). Both an equipartition analysis and a detailed afterglow mod… ▽ More AT 2022cmc is a luminous optical transient ($νL_ν \gtrsim 10^{45}$ erg s$^{-1}$) accompanied by decaying non-thermal X-rays (peak duration $t_{\rm X} \lesssim$ days and isotropic energy $E_{\rm X,iso} \gtrsim 10^{53}$ erg) and a long-lived radio/mm synchrotron afterglow, which has been interpreted as a jetted tidal disruption event (TDE). Both an equipartition analysis and a detailed afterglow model reveals the radio/mm emitting plasma to be expanding mildly relativistically (Lorentz factor $Γ\gtrsim\,few$) with an opening angle $θ_{\rm j}\simeq0.1$ and roughly fixed energy $E_{\rm j,iso} \gtrsim few \times 10^{53}$ erg into an external medium of density profile $n \propto R^{-k}$ with $k \simeq 1.5-2$, broadly similar to that of the first jetted TDE candidate Swift J1644+57 and consistent with Bondi accretion at a rate $\sim 10^{-3}\dot{M}_{\rm Edd}$ onto a $10^{6}M_{\odot}$ black hole before the outburst. The rapidly decaying optical emission over the first days is consistent with fast-cooling synchrotron radiation from the same forward shock as the radio/mm emission, while the bluer slowly decaying phase to follow likely represents a separate thermal emission component. Emission from the reverse shock may have peaked during the first days, but whose non-detection in the optical band places an upper bound $Γ_{\rm j} \lesssim 100$ on the Lorentz factor of the unshocked jet. Although a TDE origin for AT 2022cmc is indeed supported by some observations, the vast difference between the short-lived jet activity phase $t_{\rm X} \lesssim$ days relative to the months-long thermal optical emission, also challenges this scenario. A stellar core-collapse event giving birth to a magnetar or black hole engine of peak duration $\sim 1$ day offers an alternative model also consistent with the circumburst environment, if interpreted as a massive-star wind. △ Less

Submitted 20 April, 2023; v1 submitted 27 January, 2023; originally announced January 2023.

Comments: 11 pages, 9 figures, 2 tables, accepted for publication in MNRAS

A collapsar origin for GRB 211211A is (just barely) possible

Authors: Jennifer Barnes, Brian D. Metzger

Abstract: Gamma-ray bursts (GRBs) have historically been divided into two classes. Short-duration GRBs are associated with binary neutron-star mergers (NSMs), while long-duration bursts are connected to a subset of core-collapse supernovae (SNe). GRB 211211A recently made headlines as the first long-duration burst purportedly generated by an NSM. The evidence for an NSM origin was excess optical and near-in… ▽ More Gamma-ray bursts (GRBs) have historically been divided into two classes. Short-duration GRBs are associated with binary neutron-star mergers (NSMs), while long-duration bursts are connected to a subset of core-collapse supernovae (SNe). GRB 211211A recently made headlines as the first long-duration burst purportedly generated by an NSM. The evidence for an NSM origin was excess optical and near-infrared emission consistent with the kilonova observed after the gravitational wave-detected NSM GW170817. Kilonovae derive their unique electromagnetic signatures from the properties of the heavy elements synthesized by rapid neutron capture (the r-process) following the merger. Recent simulations suggest that the "collapsar" SNe that trigger long GRBs may also produce r-process elements. While observations of GRB 211211A and its afterglow ruled out an SN typical of those that follow long GRBs, an unusual collapsar could explain both the duration of GRB 211211A and the r-process-powered excess in its afterglow. We use semianalytic radiation transport modeling to evaluate low-mass collapsars as the progenitors of GRB 211211A-like events. We compare a suite of collapsar models to the afterglow-subtracted emission that followed GRB 211211A, and find the best agreement for models with high kinetic energies and an unexpected pattern of Ni-56 enrichment. We discuss how core-collapse explosions could produce such ejecta, and how distinct our predictions are from those generated by more straightforward kilonova models. We also show that radio observations can distinguish between kilonovae and the more massive collapsar ejecta we consider here. △ Less

Submitted 3 January, 2023; originally announced January 2023.

Comments: Submitted to ApJ. Comments welcome

High-Energy Neutrinos from Gamma-Ray-Faint Accretion-Powered Hypernebulae

Authors: Navin Sridhar, Brian D. Metzger, Ke Fang

Abstract: Hypernebulae are inflated by accretion-powered winds accompanying hyper-Eddington mass transfer from an evolved post-main sequence star onto a black hole or neutron star companion. The ions accelerated at the termination shock -- where the collimated fast disk winds/jet collide with the slower, wide-angled wind-fed shell -- can generate high-energy neutrinos via hadronic ($pp$) reactions, and phot… ▽ More Hypernebulae are inflated by accretion-powered winds accompanying hyper-Eddington mass transfer from an evolved post-main sequence star onto a black hole or neutron star companion. The ions accelerated at the termination shock -- where the collimated fast disk winds/jet collide with the slower, wide-angled wind-fed shell -- can generate high-energy neutrinos via hadronic ($pp$) reactions, and photohadronic ($pγ$) interactions with the disk thermal and Comptonized nonthermal background photons. It has been suggested that some fast radio bursts (FRBs) may be powered by such short-lived jetted hyper-accreting engines. Although neutrino emission associated with the ms-duration bursts themselves is challenging to detect, the persistent radio counterparts of some FRB sources -- if associated with hypernebulae -- could contribute to the high energy neutrino diffuse background flux. If the hypernebula birth rate follows that of steller-merger transients and common envelope events, we find that their volume-integrated neutrino emission -- depending on the population-averaged mass-transfer rates -- could explain up to $\sim25\%$ of the high-energy diffuse neutrino flux observed by the IceCube Observatory and the Baikal-Gigaton Volume Detector (GVD) Telescope. The time-averaged neutrino spectrum from hypernebula -- depending on the population parameters -- can also reproduce the observed diffuse neutrino spectrum. The neutrino emission could in some cases furthermore extend to >100 PeV, detectable by future ultra-high-energy neutrino observatories. The large optical depth through the nebula to Breit-Wheeler ($γγ$) interaction attenuates the escape of GeV-PeV gamma-rays co-produced with the neutrinos, rendering these gamma-ray-faint neutrino sources, consistent with the \textit{Fermi} observations of the isotropic gamma-ray background. △ Less

Submitted 1 January, 2024; v1 submitted 21 December, 2022; originally announced December 2022.

Comments: 20 pages, 10 figures, 1 table; accepted for publication in The ApJ

Journal ref: The Astrophysical Journal 960, 74 (Published on 2023-December-27)

A short gamma-ray burst from a proto-magnetar remnant

Authors: N. Jordana-Mitjans, C. G. Mundell, C. Guidorzi, R. J. Smith, E. Ramirez-Ruiz, B. D. Metzger, S. Kobayashi, A. Gomboc, I. A. Steele, M. Shrestha, M. Marongiu, A. Rossi, B. Rothberg

Abstract: The contemporaneous detection of gravitational waves and gamma rays from the GW170817/GRB 170817A, followed by kilonova emission a day after, confirmed compact binary neutron-star mergers as progenitors of short-duration gamma-ray bursts (GRBs), and cosmic sources of heavy r-process nuclei. However, the nature (and lifespan) of the merger remnant and the energy reservoir powering these bright gamm… ▽ More The contemporaneous detection of gravitational waves and gamma rays from the GW170817/GRB 170817A, followed by kilonova emission a day after, confirmed compact binary neutron-star mergers as progenitors of short-duration gamma-ray bursts (GRBs), and cosmic sources of heavy r-process nuclei. However, the nature (and lifespan) of the merger remnant and the energy reservoir powering these bright gamma-ray flashes remains debated, while the first minutes after the merger are unexplored at optical wavelengths. Here, we report the earliest discovery of bright thermal optical emission associated with the short GRB 180618A with extended gamma-ray emission, with ultraviolet and optical multicolour observations starting as soon as 1.4 minutes post-burst. The spectrum is consistent with a fast-fading afterglow and emerging thermal optical emission at 15 minutes post-burst, which fades abruptly and chromatically (flux density $F_ν \propto t^{-α}$, $α=4.6 \pm 0.3$) just 35 minutes after the GRB. Our observations from gamma rays to optical wavelengths are consistent with a hot nebula expanding at relativistic speeds, powered by the plasma winds from a newborn, rapidly-spinning and highly magnetized neutron star (i.e. a millisecond magnetar), whose rotational energy is released at a rate $L_{\rm th} \propto t^{-(2.22\pm 0.14)}$ to reheat the unbound merger-remnant material. These results suggest such neutron stars can survive the collapse to a black hole on timescales much larger than a few hundred milliseconds after the merger, and power the GRB itself through accretion. Bright thermal optical counterparts to binary merger gravitational wave sources may be common in future wide-field fast-cadence sky surveys. △ Less

Submitted 10 November, 2022; originally announced November 2022.

Journal ref: 2022 ApJ 939 106

Limits on Simultaneous and Delayed Optical Emission from Well-localized Fast Radio Bursts

Authors: Daichi Hiramatsu, Edo Berger, Brian D. Metzger, Sebastian Gomez, Allyson Bieryla, Iair Arcavi, D. Andrew Howell, Ryan Mckinven, Nozomu Tominaga

Abstract: We present the largest compilation to date of optical observations during and following fast radio bursts (FRBs). The data set includes our dedicated simultaneous and follow-up observations, as well as serendipitous archival survey observations, for a sample of 15 well-localized FRBs: eight repeating and seven one-off sources. Our simultaneous (and nearly simultaneous with a $0.4$ s delay) optical… ▽ More We present the largest compilation to date of optical observations during and following fast radio bursts (FRBs). The data set includes our dedicated simultaneous and follow-up observations, as well as serendipitous archival survey observations, for a sample of 15 well-localized FRBs: eight repeating and seven one-off sources. Our simultaneous (and nearly simultaneous with a $0.4$ s delay) optical observations of 13 (1) bursts from the repeating FRB 20220912A provide the deepest such limits to date for any extragalactic FRB, reaching a luminosity limit of $νL_ν\lesssim 10^{42}$ erg s$^{-1}$ ($\lesssim 2\times10^{41}$ erg s$^{-1}$) with $15-400$ s exposures; an optical-flux-to-radio-fluence ratio of $f_{\rm opt}/F_{\rm radio}\lesssim 10^{-7}$ ms$^{-1}$ ($\lesssim 10^{-8}$ ms$^{-1}$); and flux ratio of $f_{\rm opt}/f_{\rm radio}\lesssim 0.02-\lesssim 2\times 10^{-5}$ ($\lesssim 10^{-6}$) on millisecond to second timescales. These simultaneous limits provide useful constraints in the context of FRB emission models, such as the pulsar magnetosphere and pulsar nebula models. Interpreting all available optical limits in the context of the synchrotron maser model, we find that they constrain the flare energies to $\lesssim 10^{43}-10^{49}$ erg (depending on the distances of the various repeating FRBs, with $\lesssim 10^{39}$ erg for the Galactic SGR 1935+2154). These limits are generally at least an order of magnitude larger than those inferred from the FRBs themselves, although in the case of FRB 20220912A our simultaneous and rapid follow-up observations severely restrict the model parameter space. We conclude by exploring the potential of future simultaneous and rapid-response observations with large optical telescopes. △ Less

Submitted 6 May, 2023; v1 submitted 7 November, 2022; originally announced November 2022.

Comments: Updated to match the published letter in ApJL, 2023 April 26

Journal ref: The Astrophysical Journal Letters, 947:L28 (17pp), 2023 April 20

Evidence for multiple shocks from the $γ$-ray emission of RS Ophiuchi

Authors: Rebecca Diesing, Brian D. Metzger, Elias Aydi, Laura Chomiuk, Indrek Vurm, Siddhartha Gupta, Damiano Caprioli

Abstract: In August of 2021, Fermi-LAT, H.E.S.S., and MAGIC detected GeV and TeV $γ$-ray emission from an outburst of recurrent nova RS Ophiuchi. This detection represents the first very high energy $γ$-rays observed from a nova, and opens a new window to study particle acceleration. Both H.E.S.S. and MAGIC described the observed $γ$-rays as arising from a single, external shock. In this paper, we perform d… ▽ More In August of 2021, Fermi-LAT, H.E.S.S., and MAGIC detected GeV and TeV $γ$-ray emission from an outburst of recurrent nova RS Ophiuchi. This detection represents the first very high energy $γ$-rays observed from a nova, and opens a new window to study particle acceleration. Both H.E.S.S. and MAGIC described the observed $γ$-rays as arising from a single, external shock. In this paper, we perform detailed, multi-zone modeling of RS Ophiuchi's 2021 outburst including a self-consistent prescription for particle acceleration and magnetic field amplification. We demonstrate that, contrary to previous work, a single shock cannot simultaneously explain RS Ophiuchi's GeV and TeV emission, particularly the spectral shape and distinct light curve peaks. Instead, we put forward a model involving multiple shocks that reproduces the observed $γ$-ray spectrum and temporal evolution. The simultaneous appearance of multiple distinct velocity components in the nova optical spectrum over the first several days of the outburst supports the presence of distinct shocks, which may arise either from the strong latitudinal dependence of the density of the external circumbinary medium (e.g., in the binary equatorial plane versus the poles) or due to internal collisions within the white dwarf ejecta (as powers the $γ$-ray emission in classical novae). △ Less

Submitted 9 May, 2023; v1 submitted 3 November, 2022; originally announced November 2022.

Comments: 19 pages, 10 figures, published in 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)

Dust Echoes from Luminous Fast Blue Optical Transients

Authors: Brian D. Metzger, Daniel Perley

Abstract: Luminous fast blue optical transients (LFBOTs) such as AT2018cow form a rare class of engine-powered explosions of uncertain origin. A hallmark feature of these events is radio/millimeter synchrotron emission powered by the interaction of fast (v > 0.1 c) ejecta and dense circumstellar material (CSM) extending to large radii > 1e16 cm surrounding the progenitor. Assuming this CSM to be an outflow… ▽ More Luminous fast blue optical transients (LFBOTs) such as AT2018cow form a rare class of engine-powered explosions of uncertain origin. A hallmark feature of these events is radio/millimeter synchrotron emission powered by the interaction of fast (v > 0.1 c) ejecta and dense circumstellar material (CSM) extending to large radii > 1e16 cm surrounding the progenitor. Assuming this CSM to be an outflow from the progenitor, we show that dust grains up to ~1 micron in size can form in the outflow in the years before the explosion. This dusty CSM would attenuate the transient's ultraviolet (UV) emission prior to peak light, before being destroyed by the rising luminosity, reddening the pre-maximum colors (consistent with the pre-maximum red-to-blue color evolution of the LFBOT candidate MUSSES2020J). Re-radiation by the dust before being destroyed generates an near-infrared (NIR) "echo" of luminosity ~1e41-1e42 erg/s lasting weeks, which is detectable over the transient's rapidly fading blue continuum. We show that this dust echo is compatible with the previously unexplained NIR excess observed in AT2018cow. The gradual decay of the early NIR light curve can result from CSM which is concentrated in wide-angle outflow or torus, consistent with the highly aspherical geometry of AT2018cow's ejecta. Pre-maximum optical/UV and NIR follow-up of LFBOTs provide an new probe of their CSM environments and place additional constraints on their progenitors. △ Less

Submitted 4 October, 2022; originally announced October 2022.

Comments: 8 pages, 4 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

Cooling Envelope Model for Tidal Disruption Events

Authors: Brian D. Metzger

Abstract: We present a toy model for the thermal optical/UV/X-ray emission from tidal disruption events (TDE). Motivated by recent hydrodynamical simulations, we assume the debris streams promptly and rapidly circularize (on the orbital period of the most tightly bound debris), generating a hot quasi-spherical pressure-supported envelope of radius R_v ~ 1e14 cm (photosphere radius ~1e15 cm) surrounding the… ▽ More We present a toy model for the thermal optical/UV/X-ray emission from tidal disruption events (TDE). Motivated by recent hydrodynamical simulations, we assume the debris streams promptly and rapidly circularize (on the orbital period of the most tightly bound debris), generating a hot quasi-spherical pressure-supported envelope of radius R_v ~ 1e14 cm (photosphere radius ~1e15 cm) surrounding the supermassive black hole (SMBH). As the envelope cools radiatively, it undergoes Kelvin-Helmholtz contraction R_v ~ t^(-1), its temperature rising T_eff ~ t^(1/2) while its total luminosity remains roughly constant; the optical luminosity decays as nu L_nu ~ R_v^2 T_eff ~ t^(-3/2). Despite this similarity to the mass fall-back rate Mdot_fb ~ t^(-5/3), envelope heating from fall-back accretion is sub-dominant compared to the envelope cooling luminosity except near optical peak (where they are comparable). Envelope contraction can be delayed by energy injection from accretion from the inner envelope onto the SMBH in a regulated manner, leading to a late-time flattening of the optical/X-ray light curves, similar to those observed in some TDEs. Eventually, as the envelope contracts to near the circularization radius, the SMBH accretion rate rises to its maximum, in tandem with the decreasing optical luminosity. This cooling-induced (rather than circularization-induced) delay of up to several hundred days, may account for the delayed onset of thermal X-rays, late-time radio flares, and high-energy neutrino generation, observed in some TDEs. We compare the model predictions to recent TDE light curve correlation studies, finding agreement as well as points of tension. △ Less

Submitted 14 July, 2022; originally announced July 2022.

Comments: 12 pages, 4 figures

Rejuvenated accretors have less bound envelopes: Impact of Roche lobe overflow on subsequent common envelope events

Authors: M. Renzo, E. Zapartas, S. Justham, K. Breivik, M. Lau, R. Farmer, M. Cantiello, B. D. Metzger

Abstract: Common-envelope (CE) evolution is an outstanding open problem in stellar evolution, critical to the formation of compact binaries including gravitational-wave sources. In the ``classical'' isolated binary evolution scenario for double compact objects, the CE is usually the second mass transfer phase. Thus, the donor star of the CE is the product of a previous binary interaction, often stable Roche… ▽ More Common-envelope (CE) evolution is an outstanding open problem in stellar evolution, critical to the formation of compact binaries including gravitational-wave sources. In the ``classical'' isolated binary evolution scenario for double compact objects, the CE is usually the second mass transfer phase. Thus, the donor star of the CE is the product of a previous binary interaction, often stable Roche-lobe overflow (RLOF). Because of the accretion of mass during the first RLOF, the main-sequence core of the accretor star grows and is ``rejuvenated''. This modifies the core-envelope boundary region and decreases significantly the envelope binding energy for the remaining evolution. Comparing accretor stars from self-consistent binary models to stars evolved as single, we demonstrate that the rejuvenation can lower the energy required to eject a CE by $\sim 42-96\%$ for both black hole and neutron star progenitors, depending on the evolutionary stage and final orbital separation. Therefore, binaries experiencing first stable mass transfer may more easily survive subsequent CE events and result in possibly wider final separations compared to current predictions. Despite their high mass, our accretors also experience extended ``blue loops'', which may have observational consequences for low-metallicity stellar populations and asteroseismology. △ Less

Submitted 19 November, 2022; v1 submitted 30 June, 2022; originally announced June 2022.

Comments: accepted for publication in ApJL, reproducible at https://github.com/mathren/CE_accretors

Dissecting the microphysics behind the metallicity-dependence of massive stars radii

Authors: Chengcheng Xin, Mathieu Renzo, Brian D. Metzger

Abstract: Understanding the radii of massive stars throughout their evolution is important to answering numerous questions about stellar physics, from binary interactions on the main sequence to the pre-supernova radii. One important factor determining a star's radius is the fraction of its mass in elements heavier than Helium (metallicity, $Z$). However, the metallicity enters stellar evolution through sev… ▽ More Understanding the radii of massive stars throughout their evolution is important to answering numerous questions about stellar physics, from binary interactions on the main sequence to the pre-supernova radii. One important factor determining a star's radius is the fraction of its mass in elements heavier than Helium (metallicity, $Z$). However, the metallicity enters stellar evolution through several distinct microphysical processes, and which dominates can change throughout stellar evolution and with the overall magnitude of $Z$. We perform a series of numerical experiments with 15M$_{\odot}$ MESA models computed doubling separately the metallicity entering the radiative opacity, the equation of state, and the nuclear reaction network to isolate the impact of each on stellar radii. We explore separately models centered around two metallicity values: one near solar $Z=0.02$ and another sub-solar $Z\sim10^{-3}$, and consider several key epochs from the end of the main sequence to core carbon depletion. We find that the metallicity entering the opacity dominates at most epochs for the solar metallicity models, contributing to on average $\sim$60 - 90% of the total change in stellar radius. Nuclear reactions have a larger impact ($\sim$50 - 70%) during most epochs in the subsolar $Z$ models. The methodology introduced here can be employed more generally to propagate known microphysics errors into uncertainties on macrophysical observables including stellar radii. △ Less

Submitted 24 June, 2022; v1 submitted 22 June, 2022; originally announced June 2022.

Angular Momentum Transport in Proto-Neutron Stars and the Fate of Neutron Star Merger Remnants

Authors: Ben Margalit, Adam S. Jermyn, Brian D. Metzger, Luke F. Roberts, Eliot Quataert

Abstract: Both the core collapse of rotating massive stars, and the coalescence of neutron star (NS) binaries, result in the formation of a hot, differentially rotating NS remnant. The timescales over which differential rotation is removed by internal angular-momentum transport processes (`viscosity') has key implications for the remnant's long-term stability and the NS equation-of-state (EOS). Guided by a… ▽ More Both the core collapse of rotating massive stars, and the coalescence of neutron star (NS) binaries, result in the formation of a hot, differentially rotating NS remnant. The timescales over which differential rotation is removed by internal angular-momentum transport processes (`viscosity') has key implications for the remnant's long-term stability and the NS equation-of-state (EOS). Guided by a non-rotating model of a cooling proto-NS, we estimate the dominant sources of viscosity using an externally imposed angular velocity profile $Ω(r)$. Although the magnetorotational instability provides the dominant source of effective viscosity at large radii, convection and/or the Spruit-Tayler dynamo dominate in the core of merger remnants where $dΩ/dr \geq 0$. Furthermore, the viscous timescale in the remnant core is sufficiently short that solid body rotation will be enforced faster than matter is accreted from rotationally-supported outer layers. Guided by these results, we develop a toy model for how the merger remnant core grows in mass and angular momentum due to accretion. We find that merger remnants with sufficiently massive and slowly rotating initial cores may collapse to black holes via envelope accretion, even when the total remnant mass is less than the usually considered threshold $\approx 1.2 M_{\rm TOV}$ for forming a stable solid-body rotating NS remnant (where $M_{\rm TOV}$ is the maximum non-rotating NS mass supported by the EOS). This qualitatively new picture of the post-merger remnant evolution and stability criterion has important implications for the expected electromagnetic counterparts from binary NS mergers and for multi-messenger constraints on the NS EOS. △ Less

Submitted 22 August, 2022; v1 submitted 21 June, 2022; originally announced June 2022.

Comments: ApJ accepted version

Radio Nebulæ from Hyper-Accreting X-ray Binaries as Common Envelope Precursors and Persistent Counterparts of Fast Radio Bursts

Authors: Navin Sridhar, Brian D. Metzger

Abstract: Roche lobe overflow from a donor star onto a black hole or neutron star binary companion can evolve to a phase of unstable runaway mass-transfer, lasting as short as hundreds of orbits ($\lesssim 10^{2}$ yr for a giant donor), and eventually culminating in a common envelope event. The highly super-Eddington accretion rates achieved during this brief phase (… ▽ More Roche lobe overflow from a donor star onto a black hole or neutron star binary companion can evolve to a phase of unstable runaway mass-transfer, lasting as short as hundreds of orbits ($\lesssim 10^{2}$ yr for a giant donor), and eventually culminating in a common envelope event. The highly super-Eddington accretion rates achieved during this brief phase ($\dot{M} \gtrsim 10^{5}\dot{M}_{\rm Edd})$ are accompanied by intense mass-loss in disk winds, analogous but even more extreme than ultra-luminous X-ray (ULX) sources in the nearby universe. Also in analogy with observed ULX, this expanding outflow will inflate an energetic `bubble' of plasma into the circumbinary medium. Embedded within this bubble is a nebula of relativistic electrons heated at the termination shock of the faster $v \gtrsim 0.1 c$ wind/jet from the inner accretion flow. We present a time-dependent, one-zone model for the synchrotron radio emission and other observable properties of such ULX `hyper-nebulae'. If ULX jets are sources of repeating fast radio bursts (FRB), as recently proposed, such hyper-nebulae could generate persistent radio emission and contribute large and time-variable rotation measure to the bursts, consistent with those seen from FRB 20121102 and FRB 190520B. ULX hyper-nebulae can be discovered independent of an FRB association in radio surveys such as VLASS, as off-nuclear point-sources whose fluxes can evolve significantly on timescales as short as years, possibly presaging energetic transients from common envelope mergers. △ Less

Submitted 18 September, 2022; v1 submitted 21 June, 2022; originally announced June 2022.

Comments: 27 pages, 9 figures, 1 table; accepted for publication in ApJ

Journal ref: ApJ 937 5 (2022)

Supernova Precursor Emission and the Origin of Pre-Explosion Stellar Mass-Loss

Authors: Tatsuya Matsumoto, Brian D. Metzger

Abstract: A growing number of core collapse supernovae (SNe) which show evidence for interaction with dense circumstellar material (CSM) are accompanied by "precursor" optical emission rising weeks to months prior to the explosion. The precursor luminosities greatly exceed the Eddington limit of the progenitor star, implying they are accompanied by substantial mass-loss. Here, we present a semi-analytic mod… ▽ More A growing number of core collapse supernovae (SNe) which show evidence for interaction with dense circumstellar material (CSM) are accompanied by "precursor" optical emission rising weeks to months prior to the explosion. The precursor luminosities greatly exceed the Eddington limit of the progenitor star, implying they are accompanied by substantial mass-loss. Here, we present a semi-analytic model for SN precursor light curves which we apply to constrain the properties and mechanisms of the pre-explosion mass-loss. We explore two limiting mass-loss scenarios: (1) an "eruption" arising from shock break-out following impulsive energy deposition below the stellar surface; (2) a steady "wind" due to sustained heating of the progenitor envelope. The eruption model, which resembles a scaled-down version of Type IIP SNe, can explain the luminosities and timescales of well-sampled precursors, for ejecta masses $\sim 0.1-1\,M_{\odot}$ and velocities $\sim 100-1000\,\rm km\,s^{-1}$. By contrast, the steady-wind scenario cannot explain the highest precursor luminosities $\gtrsim10^{41}\,\rm erg\,s^{-1}$, under the constraint that the total ejecta mass not exceed the entire progenitor mass (though the less-luminous SN 2020tlf precursor can be explained by a mass-loss rate $\sim1\,M_{\odot}\,\rm yr^{-1}$). However, shock interaction between the wind and pre-existing (earlier ejected) CSM may boost its radiative efficiency and mitigate this constraint. In both eruption and wind scenarios the precursor ejecta forms compact ($\lesssim10^{15}$ cm) optically-thick CSM at the time of core collapse; though only directly observable via rapid post-explosion spectroscopy ($\lesssim$ few days before being overtaken by the SN ejecta), this material can boost the SN luminosity via shock interaction. △ Less

Submitted 23 August, 2022; v1 submitted 16 June, 2022; originally announced June 2022.

Comments: 21 pages, 8 figures, 2 tables, accepted for publication in ApJ