Researcher profile

Brian D. Metzger

Brian D. Metzger contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
31works
0followers
11topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

31 published item(s)

preprint2026arXiv

No Sign of a Magnetar Remnant Following the Kilonova-Producing Long GRB 211211A $\sim 1.7~$Years Later

In addition to a $γ$-ray burst (GRB), the merger of two neutron stars may produce a temporarily or indefinitely stable neutron star remnant with a strong magnetic field (a "magnetar"). As this magnetar remnant spins down, it can deposit its rotational energy into the surrounding kilonova ejecta, producing synchrotron emission that peaks in the radio bands $\sim$months-years after the merger ("boosted kilonova"). The nearby ($z=0.0763$) long-duration GRB 211211A, which has an apparent kilonova counterpart and likely neutron star merger progenitor, may have produced such a remnant. We observed the location of GRB 211211A at 6 GHz with the NSF's Karl G. Jansky Very Large Array (VLA) spanning $\approx 0.54$-$1.7~$years after the burst. We do not detect any radio emission, placing strong limits on the energy deposited into the ejecta by any remnant to $\lesssim 4.4 \times 10^{52}~{\rm erg}$. Due to the proximity of the event, we are also able to place limits on a kilonova afterglow that did not receive any additional energy deposition, though it is possible such emission will be suppressed until $\sim 4~{\rm years}$ after the burst, when the kilonova is expected to overtake the forward shock of the GRB. Future observations with the VLA and next-generation radio facilities will be able to further constrain the magnetar-boosted kilonova and kilonova afterglow scenarios, as well as directly constrain models in the scenario that GRB 211211A was instead produced by a collapsar.

preprint2026arXiv

Relativistic jets from millisecond proto-magnetars

Rapidly rotating, strongly magnetized neutron stars (``millisecond proto-magnetars'') formed in stellar core-collapse, neutron star mergers, and white dwarf accretion-induced collapse have long been proposed as central engines of gamma-ray bursts (GRB) and accompanying supernovae/kilonovae. However, during the first few seconds after birth, neutrino heating drives baryon-rich winds from the neutron star surface, potentially limiting the magnetization and achievable Lorentz factors of the outflow and casting doubt on whether proto-magnetars can launch ultra-relativistic jets at early times, as needed to power short-duration GRB. We present 3D general-relativistic magnetohydrodynamic simulations of neutrino-heated proto-magnetar winds that incorporate M0 neutrino transport. While the global wind properties broadly agree with previous analytic estimates calibrated to one-dimensional models, our simulations reveal essential multidimensional effects. For rapidly rotating models with spin periods P = 1 ms, centrifugal forces strongly enhance mass loss near the rotational equator, producing a dense, sub-relativistic outflow ( ~0.1c). This equatorial wind naturally confines and collimates less baryon-loaded outflows emerging from higher latitudes, leading to the formation of a structured bipolar jet with a peak magnetization up to ~ 30-100 along the pole, sufficient to reach bulk Lorentz factors ~ 100 on larger scales. The resulting angular stratification of the outflow energy into ultra-relativistic polar and sub-relativistic equatorial components is broadly consistent with the observed partition between beaming-corrected GRB energies and supernova/kilonova ejecta. Our results demonstrate that millisecond proto-magnetars can launch relativistic jets within seconds of formation and highlight their potential role in powering the diverse electromagnetic counterparts of compact-object explosions.

preprint2024arXiv

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

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.

preprint2022arXiv

Bumpy Declining Light Curves Are Common in Hydrogen-poor Superluminous Supernovae

Recent work has revealed that the light curves of hydrogen-poor (Type I) superluminous supernovae (SLSNe), thought to be powered by magnetar central engines, do not always follow the smooth decline predicted by a simple magnetar spin-down model. Here we present the first systematic study of the prevalence and properties of "bumps" in the post-peak light curves of 34 SLSNe. We find that the majority (44-76%) of events cannot be explained by a smooth magnetar model alone. We do not find any difference in supernova properties between events with and without bumps. By fitting a simple Gaussian model to the light-curve residuals, we characterize each bump with an amplitude, temperature, phase, and duration. We find that most bumps correspond with an increase in the photospheric temperature of the ejecta, although we do not see drastic changes in spectroscopic features during the bump. We also find a moderate correlation ($ρ\approx0.5$; $p\approx0.01$) between the phase of the bumps and the rise time, implying that such bumps tend to happen at a certain "evolutionary phase," $(3.7\pm1.4)t_\mathrm{rise}$. Most bumps are consistent with having diffused from a central source of variable luminosity, although sources further out in the ejecta are not excluded. With this evidence, we explore whether the cause of these bumps is intrinsic to the supernova (e.g., a variable central engine) or extrinsic (e.g., circumstellar interaction). Both cases are plausible, requiring low-level variability in the magnetar input luminosity, small decreases in the ejecta opacity, or a thin circumstellar shell or disk.

preprint2022arXiv

Luminous Fast Blue Optical Transients and Type Ibn/Icn SNe from Wolf-Rayet/Black Hole Mergers

Progenitor models for the &#34;luminous&#34; subclass of Fast Blue Optical Transients (LFBOTs; prototype: AT2018cow) are challenged to simultaneously explain all of their observed properties: fast optical rise times < days; peak luminosities >1e44 erg/s; low yields <0.1 Msun of 56Ni; aspherical ejecta with a wide velocity range (<3000 km/s to > 0.1-0.5 c with increasing polar latitude); presence of hydrogen-depleted-but-not-free dense circumstellar material (CSM) on radial scales from ~1e14 cm to ~3e16 cm; embedded variable source of non-thermal X-ray/gamma-rays, suggestive of a compact object. We show that all of these properties are consistent with the tidal disruption and hyper-accretion of a Wolf-Rayet (WR) star by a black hole (BH) or neutron star (NS) binary companion. In contrast with related previous models, the merger occurs with a long delay (>~ 100 years) following the common envelope (CE) event responsible for birthing the binary, as a result of gradual angular momentum loss to a relic circumbinary disk. Disk-wind outflows from the merger-generated accretion flow generate the 56Ni-poor aspherical ejecta with the requisite velocity range. The optical light curve is powered primarily by reprocessing X-rays from the inner accretion flow/jet, though CSM shock interaction also contributes. Primary CSM sources include mass-loss from WR star (e.g., from the L2 point) during the earliest stages of the merger (~<1e14 cm) and the relic CE disk and its photoevaporation-driven wind (>~ 1e16 cm). Longer delayed mergers may instead give rise to supernovae Type Ibn/Icn (depending on the WR evolutionary state), potentially connecting these transient classes with LFBOTs.

preprint2022arXiv

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

A growing number of core collapse supernovae (SNe) which show evidence for interaction with dense circumstellar material (CSM) are accompanied by &#34;precursor&#34; 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 &#34;eruption&#34; arising from shock break-out following impulsive energy deposition below the stellar surface; (2) a steady &#34;wind&#34; 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.

preprint2022arXiv

The first nova eruption in a novalike variable: YZ Ret as seen in X-rays and gamma-rays

Peaking at 3.7 mag on 2020 July 11, YZ Ret was the second-brightest nova of the decade. The nova&#39;s moderate proximity (2.7 kpc from Gaia) provided an opportunity to explore its multi-wavelength properties in great detail. Here we report on YZ Ret as part of a long-term project to identify the physical mechanisms responsible for high-energy emission in classical novae. We use simultaneous Fermi/LAT and NuSTAR observations complemented by XMM-Newton X-ray grating spectroscopy to probe the physical parameters of the shocked ejecta and the nova-hosting white dwarf. The XMM-Newton observations revealed a super-soft X-ray emission which is dominated by emission lines of CV, CVI, NVI, NVII, and OVIII rather than a blackbody-like continuum, suggesting CO-composition of the white dwarf in a high-inclination binary system. Fermi/LAT detected YZ Ret for 15 days with the gamma-ray spectrum best described by a power law with an exponential cut-off at 1.9 +/-0.6 GeV. In stark contrast with theoretical predictions and in keeping with previous NuSTAR observations of Fermi-detected classical novae (V5855 Sgr and V906 Car), the 3.5-78 keV X-ray emission is found to be two orders of magnitude fainter than the GeV emission. The X-ray emission observed by NuSTAR is consistent with a single-temperature thermal plasma. We detect no non-thermal tail of the GeV emission expected to extend down to the NuSTAR band. NuSTAR observations continue to challenge theories of high-energy emission from shocks in novae.

preprint2022arXiv

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

We explore the effects of rapid rotation on the properties of neutrino-heated winds from proto-neutron stars (PNS) formed in core-collapse supernovae or neutron-star mergers by means of three-dimensional general-relativistic hydrodynamical simulations with M0 neutrino transport. We focus on conditions characteristic of a few seconds into the PNS cooling evolution when the neutrino luminosities obey $L_{ν_e} + L_{\barν_e} \approx 7\times 10^{51}$ erg s$^{-1}$, and over which most of the wind mass-loss will occur. After an initial transient phase, all of our models reach approximately steady-state outflow solutions with positive energies and sonic surfaces captured on the computational grid. Our non-rotating and slower-rotating models (angular velocity relative to Keplerian $Ω/Ω_{\rm K} \lesssim 0.4$; spin period $P \gtrsim 2$ ms) generate approximately spherically symmetric outflows with properties in good agreement with previous PNS wind studies. By contrast, our most rapidly spinning PNS solutions ($Ω/Ω_{\rm K} \gtrsim 0.75$; $P \approx 1$ ms) generate outflows focused in the rotational equatorial plane with much higher mass-loss rates (by over an order of magnitude), lower velocities, lower entropy, and lower asymptotic electron fractions, than otherwise similar non-rotating wind solutions. Although such rapidly spinning PNS are likely rare in nature, their atypical nucleosynthetic composition and outsized mass yields could render them important contributors of light neutron-rich nuclei compared to more common slowly rotating PNS birth. Our calculations pave the way to including the combined effects of rotation and a dynamically-important large-scale magnetic field on the wind properties within a 3D GRMHD framework.

preprint2022arXiv

X-ray emission from candidate stellar merger remnant TYC 2597-735-1 and its Blue Ring Nebula

Tight binary or multiple star systems can interact through mass transfer and follow vastly different evolutionary pathways than single stars. The star TYC 2597-735-1 is a candidate for a recent stellar merger remnant resulting from a coalescence of a low-mass companion with a primary star a few thousand years ago. This violent event is evident in a conical outflow (&#34;Blue Ring Nebula&#34;) emitting in UV light and surrounded by leading shock filaments observed in H$α$ and UV emission. From Chandra data, we report the detection of X-ray emission from the location of TYC 2597-735-1 with a luminosity $\log(L_\mathrm{X}/L_\mathrm{bol})=-5.5$. Together with a previously reported period around 14~days, this indicates ongoing stellar activity and the presence of strong magnetic fields on TYC 2597-735-1. Supported by stellar evolution models of merger remnants, we interpret the inferred stellar magnetic field as dynamo action associated with a newly formed convection zone in the atmosphere of TYC 2597-735-1, though internal shocks at the base of an accretion-powered jet cannot be ruled out. We speculate that this object will evolve into an FK Com type source, i.e. a class of rapidly spinning magnetically active stars for which a merger origin has been proposed but for which no relic accretion or large-scale nebula remains visible. We also detect likely X-ray emission from two small regions close to the outer shock fronts in the Blue Ring Nebula, which may arise from either inhomogenities in the circumstellar medium or in the mass and velocity distribution in the merger-driven outflow.

preprint2021arXiv

A Program for Multi-Messenger Standard Siren Cosmology in the Era of LIGO A+, Rubin Observatory, and Beyond

The most promising variation of the standard siren technique combines gravitational-wave (GW) data for binary neutron star (BNS) mergers with redshift measurements enabled by their electromagnetic (EM) counterparts, to constrain cosmological parameters such as $H_0$, $Ω_m$, and $w_0$. Here we evaluate the near- and long-term prospects of multi-messenger cosmology in the era of future GW observatories: Advanced LIGO Plus (A+, 2025), Voyager-like detectors (2030s), and Cosmic Explorer-like detectors (CE, 2035 and beyond). We show that the BNS horizon distance of $\approx 700$ Mpc for A+ is well-matched to the sensitivity of the Vera C. Rubin Observatory (VRO) for kilonova detections. We find that one year of joint A+ and VRO observations will constrain the value of $H_0$ to percent-level precision, given a small investment of VRO time dedicated to target-of-opportunity GW follow-up. In the Voyager era, the BNS-kilonova observations begin to constrain $Ω_m$ with an investment of a few percent of VRO time. With the larger BNS horizon distance in the Cosmic Explorer era, on-axis short gamma-ray bursts (SGRBs) and their afterglows (though accompanying only some of the GW-detected mergers) supplant kilonovae as the most promising counterparts for redshift identification. We show that five years of joint observations with Cosmic Explorer-like facilities and a next-generation gamma-ray satellite with localization capabilities similar to that presently possible with Swift could constrain both $Ω_m$ and $w_0$ to $15-20\%$. We therefore advocate for a robust target-of-opportunity (ToO) program with VRO, and a wide-field gamma-ray satellite with improved sensitivity in the 2030s, to enable standard siren cosmology with next-generation gravitational wave facilities.

preprint2021arXiv

A Toy Model for the Time-Frequency Structure of Fast Radio Bursts: Implications for the CHIME Burst Dichotomy

We introduce a toy model for the time-frequency structure of fast radio bursts (FRB), in which the observed emission is produced as a narrowly-peaked intrinsic spectral energy distribution sweeps down in frequency across the instrumental bandpass as a power-law in time. Though originally motivated by emission models which invoke a relativistic shock, the model could in principle apply to a wider range of emission scenarios. We quantify the burst&#39;s detectability using the frequency bandwidth over which most of its signal-to-noise ratio (SNR) is accumulated. We demonstrate that by varying just a single parameter of the toy model-the power-law index βof the frequency drift rate-one can transform a long (and hence preferentially time-resolved) burst with a narrow time-integrated spectrum into a shorter burst with a broad power-law time-integrated spectrum. We suggest that burst-to-burst diversity in the value of βcould generate the dichotomy between burst duration and frequency-width recently found by CHIME. In shock models, the value of βis related to the radial density profile of external medium, which in light of the preferentially longer duration of bursts from repeating sources may point to diversity in the external environments surrounding repeating versus one-off FRB sources.

preprint2021arXiv

High-Energy Neutrinos and Gamma-Rays from Non-Relativistic Shock-Powered Transients

Shock interaction has been argued to play a role in powering a range of optical transients, including supernovae (particularly the superluminous class), classical novae, stellar mergers, tidal disruption events, and fast blue optical transients. These same shocks can accelerate relativistic ions, generating high-energy neutrino and gamma-ray emission via hadronic pion production. The recent discovery of time-correlated optical and gamma-ray emission in classical novae has revealed the important role of radiative shocks in powering these events, enabling an unprecedented view of the properties of ion acceleration, including its efficiency and energy spectrum, under similar physical conditions to shocks in extragalactic transients. Here we introduce a model for connecting the radiated optical fluence of non-relativistic transients to their maximal neutrino and gamma-ray fluence. We apply this technique to a wide range of extragalactic transient classes in order to place limits on their contributions to the cosmological high-energy gamma-ray and neutrino backgrounds. Based on a simple model for diffusive shock acceleration at radiative shocks, calibrated to novae, we demonstrate that several of the most luminous transients can accelerate protons up to energies $E_{\rm max} \gtrsim 10^{16}$ eV, sufficient to contribute to the IceCube astrophysical background. Furthermore, several of the considered sources$-$particularly hydrogen-poor supernovae$-$may serve as &#34;hidden&#34; gamma-ray sources due to the high gamma-ray opacity of their ejecta, evading constraints imposed by the non-blazar Fermi-LAT background. However, adopting an ion acceleration efficiency $\sim$ 0.3-1$\%$ motivated by nova observations, we find that currently known classes of non-relativistic, potentially shock-powered transients contribute at most a few percent of the total IceCube background.

preprint2021arXiv

Interacting Stellar EMRIs as Sources of Quasi-Periodic Eruptions in Galactic Nuclei

A star that approaches a supermassive black hole (SMBH) on a circular extreme mass ratio inspiral (EMRI) can undergo Roche lobe overflow (RLOF), resulting in a phase of long-lived mass-transfer onto the SMBH. If the interval separating consecutive EMRIs is less than the mass-transfer timescale driven by gravitational wave emission (typically ~1-10 Myr), the semi-major axes of the two stars will approach each another on scales of <~ hundreds to thousands of gravitational radii. Close flybys tidally strip gas from one or both RLOFing stars, briefly enhancing the mass-transfer rate onto the SMBH and giving rise to a flare of transient X-ray emission. If both stars reside in an common orbital plane, these close interactions will repeat on a timescale as short as hours, generating a periodic series of flares with properties (amplitudes, timescales, sources lifetimes) remarkably similar to the &#34;quasi-periodic eruptions&#34; (QPEs) recently observed from galactic nuclei hosting low-mass SMBHs. A cessation of QPE activity is predicted on a timescale of months to years, due to nodal precession of the EMRI orbits out of alignment by the SMBH spin. Channels for generating the requisite coplanar EMRIs include the tidal separation of binaries (Hills mechanism) or Type I inwards migration through a gaseous AGN disk. Alternative scenarios for QPEs, that invoke single stellar EMRIs on an eccentric orbit undergoing a runaway sequence of RLOF events, are strongly disfavored by formation rate constraints.

preprint2021arXiv

Multi-messenger Bayesian parameter inference of a binary neutron-star merger

The combined detection of a binary neutron-star merger in both gravitational waves (GWs) and electromagnetic (EM) radiation spanning the entire spectrum -- GW170817 / AT2017gfo / GRB170817A -- marks a breakthrough in the field of multi-messenger astronomy. Between the plethora of modeling and observations, the rich synergy that exists among the available data sets creates a unique opportunity to constrain the binary parameters, the equation of state of supranuclear density matter, and the physical processes at work during the kilonova and gamma-ray burst. We report, for the first time, Bayesian parameter estimation combining information from GW170817, AT2017gfo, GRB170817 to obtain truly multi-messenger constraints on the tidal deformability $\tildeΛ \in [302,860]$, total binary mass $M \in [2.722,2.751] M_\odot$, the radius of a $1.4$ solar mass neutron star $R \in [11.3,13.5] \rm km$ (with additional $0.2\ \rm km$ systematic uncertainty), and an upper bound on the mass ratio of $q \leq 1.27$, all at 90% confidence. Our joint novel analysis makes use of new phenomenological descriptions of the dynamical ejecta, debris disk mass, and remnant black hole properties, all derived from a large suite of numerical relativity simulations.

preprint2021arXiv

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

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

preprint2021arXiv

Shock-powered radio precursors of neutron star mergers from accelerating relativistic binary winds

During the final stages of a compact object merger, if at least one of the binary components is a magnetized neutron star (NS), then its orbital motion substantially expands the NS&#39;s open magnetic flux -- and hence increases its wind luminosity -- relative to that of an isolated pulsar. As the binary orbit shrinks due to gravitational radiation, the power and speed of this binary-induced inspiral wind may (depending on pair loading) secularly increase, leading to self-interaction and internal shocks in the outflow beyond the binary orbit. The magnetized forward shock can generate coherent radio emission via the synchrotron maser process, resulting in an observable radio precursor to binary NS merger. We perform 1D relativistic hydrodynamical simulations of shock interaction in the accelerating binary NS wind, assuming that the inspiral wind efficiently converts its Poynting flux into bulk kinetic energy prior to the shock radius. This is combined with the shock maser spectrum from particle-in-cell simulations, to generate synthetic radio light curves. The precursor burst with a fluence of $\sim1$ Jy$\cdot$ms at $\sim$GHz frequencies lasts $\sim 1-500$ ms following the merger for a source at $\sim3$ Gpc ($B_{\rm d}/10^{12}$ G)$^{8/9}$, where $B_{\rm d}$ is the dipole field strength of the more strongly-magnetized star. Given an outflow geometry concentrated along the binary equatorial, the signal may be preferentially observable for high-inclination systems, i.e. those least likely to produce a detectable gamma-ray burst.

preprint2021arXiv

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

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

preprint2020arXiv

Collapsars as a major source of r-process elements

The production of elements by rapid neutron capture (r-process) in neutron-star mergers is expected theoretically and is supported by multimessenger observations of gravitational-wave event GW170817: this production route is in principle sufficient to account for most of the r-process elements in the Universe. Analysis of the kilonova that accompanied GW170817 identified delayed outflows from a remnant accretion disk formed around the newly born black hole as the dominant source of heavy r-process material from that event. Similar accretion disks are expected to form in collapsars (the supernova-triggering collapse of rapidly rotating massive stars), which have previously been speculated to produce r-process elements. Recent observations of stars rich in such elements in the dwarf galaxy Reticulum II, as well as the Galactic chemical enrichment of europium relative to iron over longer timescales, are more consistent with rare supernovae acting at low stellar metallicities than with neutron-star mergers. Here we report simulations that show that collapsar accretion disks yield sufficient r-process elements to explain observed abundances in the Universe. Although these supernovae are rarer than neutron-star mergers, the larger amount of material ejected per event compensates for the lower rate of occurrence. We calculate that collapsars may supply more than 80 per cent of the r-process content of the Universe.

preprint2020arXiv

Constraints on the Engines of Fast Radio Bursts

We model the sample of fast radio bursts (FRB), including the newly discovered CHIME repeaters, using the synchrotron blast wave model of Metzger, Margalit & Sironi (2019). This model postulates that FRBs are precursor radiation from ultra-relativistic magnetized shocks generated as flare ejecta from a central engine collide with an effectively stationary external medium. Downward drifting of the burst frequency structure naturally arises from deceleration of the blast-wave. The data are consistent with FRBs being produced by flares of energy $E_{\rm flare} \sim 10^{43}-10^{46}(f_ξ/10^{-3})^{-4/5}$ erg, where $f_ξ$ is the maser efficiency, and minimum bulk Lorentz factors $Γ\approx 10^2-10^3$, which generate the observed FRBs at shock radii $r_{\rm sh} \sim 10^{12}-10^{13}$ cm. We infer upstream densities $n_{\rm ext}(r_{\rm sh}) \sim 10^{2}-10^{4}$ cm$^{-3}$ and radial profiles $n_{\rm ext} \propto r^{-k}$ showing a range of slopes $k \approx [-2,1]$ (which are seen to evolve between bursts), broadly consistent with the upstream medium being the inner edge of an ion-loaded shell released by a recent energetic flare. The burst timescales, energetics, rates, and external medium properties are consistent with repeating FRBs arising from young, hyper-active flaring magnetars, but the methodology presented is generally applicable to any central engine which injects energy impulsively into a dense magnetized medium. Uncertainties and variations of the model are discussed, including the effects of the strong electric field of the FRB wave (strength parameter $a \gg 1$) on the upstream medium. One-dimensional particle-in-cell simulations of magnetized shocks into a pair plasma are presented which demonstrate that high maser efficiency can be preserved, even in the limit $a \gg 1$ in which the FRB wave accelerates the upstream electrons to ultra-relativistic speeds.

preprint2020arXiv

Direct evidence for shock-powered optical emission in a nova

Classical novae are thermonuclear explosions that occur on the surfaces of white dwarf stars in interacting binary systems (Bode & Evans 2008). It has long been thought that the luminosity of classical novae is powered by continued nuclear burning on the surface of the white dwarf after the initial runaway (Gallaher & Starrfield 1978). However, recent observations of GeV $γ$-rays from classical novae have hinted that shocks internal to the nova ejecta may dominate the nova emission. Shocks have also been suggested to power the luminosity of events as diverse as stellar mergers (Metzger & Pejcha 2017), supernovae (Moriya et al. 2018), and tidal disruption events (Roth et al. 2016), but observational confirmation has been lacking. Here we report simultaneous space-based optical and $γ$-ray observations of the 2018 nova V906 Carinae (ASASSN-18fv), revealing a remarkable series of distinct correlated flares in both bands. The optical and $γ$-ray flares occur simultaneously, implying a common origin in shocks. During the flares, the nova luminosity doubles, implying that the bulk of the luminosity is shock-powered. Furthermore, we detect concurrent but weak X-ray emission from deeply embedded shocks, confirming that the shock power does not appear in the X-ray band and supporting its emergence at longer wavelengths. Our data, spanning the spectrum from radio to $γ$-ray, provide direct evidence that shocks can power substantial luminosity in classical novae and other optical transients.

preprint2020arXiv

Electromagnetic transients and gravitational waves from white dwarf disruptions by stellar black holes in triple systems

Mergers of binaries comprised of compact objects can give rise to explosive transient events, heralding the birth of exotic objects which cannot be formed through single star evolution. Using a large number of direct N-body simulations, we explore the possibility that a white dwarf (WD) is dynamically driven to tidal disruption by a stellar-mass black hole (BH) as a consequence of the joint effects of gravitational wave (GW) emission and Lidov-Kozai oscillations imposed by the tidal field of a outer tertiary companion orbiting the inner BH-WD binary. We explore the sensitivity of our results to the distributions of natal kick velocities imparted to the BH and WD upon formation, adiabatic mass loss, semi-major axes and eccentricities of the triples, and stellar mass ratios. We find rates of WD-TDEs in the range $1.2\times 10^{-3}-1.4$ Gpc$^{-3}$ yr$^{-1}$ for $z\leq 0.1$, rarer than stellar TDEs in triples by a factor of $\sim 3$--$30$. The uncertainty in the TDE rates may be greatly reduced in the future using gravitational wave (GW) observations of Galactic binaries and triples with LISA. WD-TDEs may give rise to high energy X-ray or gamma-ray transients of duration similar to long gamma-ray bursts but lacking the signatures of a core-collapse supernova, while being accompanied by a supernova-like optical transient which lasts for only days. WD--BH and WD--NS binaries will also emit GWs in the LISA band before the TDE. The discovery and identification of triple-induced WD-TDE events by future time domain surveys and/or GWs could enable the study of the demographics of BHs in nearby galaxies.

preprint2020arXiv

Fermi-LAT Observations of V549 Vel 2017: a Sub-Luminous Gamma-Ray Nova?

We report on the Fermi-LAT detection (with ~5.7 sigma significance) as well as the multi-wavelength analysis of the 2017 nova eruption V549 Vel. Unlike the recent shock-powered novae ASASSN-16ma and V906 Car, the optical and gamma-ray light curves of V549 Vel show no correlation, likely implying relatively weak shocks in the eruption. Gaia detected a candidate progenitor of V549 Vel and found a parallax measurement of 1.91+/-0.39 mas, equivalent to a mode distance of d=560 pc (90% credible interval of 380-1050 pc). The progenitor was also observed by the 2MASS and WISE surveys. When adopting the Gaia distance, the spectral energy distribution of the progenitor is close to that of a G-type star. Swift/XRT detected the supersoft X-ray emission of the nova (kT=30-40 keV) since day 236, and the inferred blackbody size is comparable to that of other novae assuming d=560 pc (i.e., R_bb~5*10^8 cm). However, there is also an unknown astrometric excess noise of 3.2 mas found in the Gaia data, and the inferred distance becomes controversial. If the Gaia distance is accurate, the gamma-ray luminosity of V549 Vel will be as low as L_gamma~4*10^33 erg/s, making it the least luminous gamma-ray nova known so far. This may imply that the shock properties responsible for the gamma-ray emission in V549 Vel are different from those of the more luminous events. If the nova is located farther away, it is likely a symbiotic system with a giant companion as the observed progenitor.

preprint2020arXiv

Implications of a &#34;Fast Radio Burst&#34; from a Galactic Magnetar

A luminous radio burst was recently detected in temporal coincidence with a hard X-ray flare from the Galactic magnetar SGR 1935+2154 with a time and frequency structure consistent with cosmological fast radio bursts (FRB) and a fluence within a factor of $\lesssim 10$ of the least energetic extragalactic FRB previously detected. Although active magnetars are commonly invoked FRB sources, several distinct mechanisms have been proposed for generating the radio emission which make different predictions for the accompanying higher frequency radiation. We show that the properties of the coincident radio and X-ray flares from SGR 1935+2154, including their approximate simultaneity and relative fluence $E_{\rm radio}/E_{\rm X} \sim 10^{-5}$, as well as the duration and spectrum of the X-ray emission, are consistent with extant predictions for the synchrotron maser shock model. Rather than arising from the inner magnetosphere, the X-rays are generated by (incoherent) synchrotron radiation from thermal electrons heated at the same shocks which produce the coherent maser emission. Although the rate of SGR 1935+2154-like bursts in the local universe is not sufficient to contribute appreciably to the extragalactic FRB rate, the inclusion of an additional population of more active magnetars with stronger magnetic fields than the Galactic population can explain both the FRB rate as well as the repeating fraction, however only if the population of active magnetars are born at a rate that is at least two-orders of magnitude lower than that of SGR 1935+2154-like magnetars. This may imply that the more active magnetar sources are not younger magnetars formed in a similar way to the Milky Way population (e.g. via ordinary supernovae), but instead through more exotic channels such as superluminous supernovae, accretion-induced collapse or neutron star mergers.

preprint2020arXiv

Periodicity in recurrent fast radio bursts and the origin of ultra long period magnetars

The recurrent fast radio burst FRB 180916 was recently shown to exhibit a 16 day period (with possible aliasing) in its bursting activity. Given magnetars as widely considered FRB sources, this period has been attributed to precession of the magnetar spin axis or the orbit of a binary companion. Here, we make the simpler connection to a {\it rotational period}, an idea observationally motivated by the 6.7 hour period of the Galactic magnetar candidate, 1E 161348--5055. We explore three physical mechanisms that could lead to the creation of ultra long period magnetars: (i) enhanced spin-down due to episodic mass-loaded charged particle winds (e.g. as may accompany giant flares), (ii) angular momentum kicks from giant flares and (iii) fallback leading to long lasting accretion disks. We show that particle winds and fallback accretion can potentially lead to a sub-set of the magnetar population with ultra long periods, sufficiently long to accommodate FRB 180916 or 1E 161348--5055. If confirmed, such periods implicate magnetars in relatively mature states (ages $1-10$ kyr) and which possessed large internal magnetic fields at birth $B_{\rm int}\gtrsim 10^{16}$ G. In the low-twist magnetar model for FRBs, such long period magnetars may dominate FRB production for repeaters at lower isotropic-equivalent energies and broaden the energy distribution beyond that expected for a canonical population of magnetars which terminate their magnetic activity at shorter periods $P \lesssim 10$ s.

preprint2020arXiv

Reconstructing Masses of Merging Neutron Stars from Stellar $R$-Process Abundance Signatures

Neutron star mergers (NSMs) are promising astrophysical sites for the rapid neutron-capture (&#34;$r$-&#34;) process, but can their integrated yields explain the majority of heavy-element material in the Galaxy? One method to address this question has utilized a forward approach that propagates NSM rates and yields along with stellar formation rates, in the end comparing those results with observed chemical abundances of $r$-process-rich, metal-poor stars. In this work, we take the inverse approach by utilizing $r$-process-element abundance ratios of metal-poor stars as input to reconstruct the properties---especially the masses---of the neutron star (NS) binary progenitors of the $r$-process stars. This novel analysis provides an independent avenue for studying the population of the original neutron star binary systems that merged and produced the $r$-process material incorporated in Galactic metal-poor halo stars. We use ratios of elements typically associated with the limited-$r$ process and the actinide region to those in the lanthanide region (i.e., Zr/Dy and Th/Dy) to probe the NS masses of the progenitor merger. We find that NSMs can account for all $r$-process material in metal-poor stars that display $r$-process signatures, while simultaneously reproducing the present-day distribution of double-NS (DNS) systems. However, the most $r$-process enhanced stars (the $r$-II stars) on their own would require progenitor NSMs of very asymmetric systems that are distinctly different from present ones in the Galaxy. As this analysis is model-dependent, we also explore variations in line with future expectation regarding potential theoretical and observational updates, and comment on how these variations impact our results.

preprint2020arXiv

Variability in Short Gamma-ray Bursts: Gravitationally Unstable Tidal Tails

Short gamma-ray bursts are thought to result from the mergers of two neutron stars or a neutron star and stellar mass black hole. The final stages of the merger are generally accompanied by the production of one or more tidal &#34;tails&#34; of ejecta, which fall back onto the remnant-disc system at late times. Using the results of a linear stability analysis, we show that if the material comprising these tails is modeled as adiabatic and the effective adiabatic index satisfies $γ\ge 5/3$, then the tails are gravitationally unstable and collapse to form small-scale knots. We analytically estimate the properties of these knots, including their spacing along the tidal tail and the total number produced, and their effect on the mass return rate to the merger remnant. We perform hydrodynamical simulations of the disruption of a polytropic (with the polytropic and adiabatic indices $γ$ equal), $γ=2$ neutron star by a black hole, and find agreement between the predictions of the linear stability analysis and the distribution of knots that collapse out of the instability. The return of these knots to the black hole induces variability in the fallback rate, which can manifest as variability in the lightcurve of the GRB and -- depending on how rapidly the instability operates -- the prompt emission. The late-time variability induced by the return of these knots is also consistent with the extended emission observed in some GRBs.

preprint2020arXiv

X-ray Spectroscopy of the gamma-ray Brightest Nova V906 Car (ASASSN-18fv)

Shocks in gamma-ray emitting classical novae are expected to produce bright thermal and non-thermal X-rays. We test this prediction with simultaneous NuSTAR and Fermi/LAT observations of nova V906 Car, which exhibited the brightest GeV gamma-ray emission to date. The nova is detected in hard X-rays while it is still gamma-ray bright, but contrary to simple theoretical expectations, the detected 3.5-78 keV emission of V906 Car is much weaker than the simultaneously observed >100 MeV emission. No non-thermal X-ray emission is detected, and our deep limits imply that the gamma-rays are likely hadronic. After correcting for substantial absorption (N_H ~ 2 x 10^23 cm^-2), the thermal X-ray luminosity (from a 9 keV optically-thin plasma) is just ~2% of the gamma-ray luminosity. We consider possible explanations for the low thermal X-ray luminosity, including the X-rays being suppressed by corrugated, radiative shock fronts or the X-rays from the gamma-ray producing shock are hidden behind an even larger absorbing column (N_H >10^25 cm^-2). Adding XMM-Newton and Swift/XRT observations to our analysis, we find that the evolution of the intrinsic X-ray absorption requires the nova shell to be expelled 24 days after the outburst onset. The X-ray spectra show that the ejecta are enhanced in nitrogen and oxygen, and the nova occurred on the surface of a CO-type white dwarf. We see no indication of a distinct super-soft phase in the X-ray lightcurve, which, after considering the absorption effects, may point to a low mass of the white dwarf hosting the nova.

preprint2019arXiv

Fast Radio Bursts from Magnetars Born in Binary Neutron Star Mergers and Accretion Induced Collapse

Recently born magnetars are promising candidates for the engines powering fast radio bursts (FRBs). The focus thus far has been placed on millisecond magnetars born in rare core-collapse explosions, motivated by the star forming dwarf host galaxy of the repeating FRB 121102, which is remarkably similar to the hosts of superluminous supernovae (SLSNe) and long gamma-ray bursts (LGRB). However, long-lived magnetars may also be created in binary neutron star (BNS) mergers, in the small subset of cases with a sufficiently low total mass for the remnant to avoid collapse to a black hole, or in the accretion-induced collapse (AIC) of a white dwarf. A BNS FRB channel will be characterized by distinct host galaxy and spatial offset distributions than the SLSNe/LGRB channel; we anticipate a similar host population, although possibly different offset distribution for AIC events. We show that both the BNS and AIC channels are consistent with the recently reported FRB 180924, localized by ASKAP to a massive quiescent host galaxy with an offset of about 1.4 effective radii. Using models calibrated to FRB 121102, we make predictions for the dispersion measure, rotation measure, and persistent radio emission from magnetar FRB sources born in BNS mergers or AIC, and show these are consistent with upper limits from FRB 180924 for reasonable parameters. Depending on the rate of AIC, and the fraction of BNS mergers leaving long-lived stable magnetars, the birth rate of repeating FRB sources associated with older stellar populations could be comparable to that of the core-collapse channel. We also discuss potential differences in the repetition properties of these channels, as a result of differences in the characteristic masses and magnetic fields of the magnetars.

preprint2019arXiv

Internal Shocks from Variable Outflows in Classical Novae

We present one-dimensional hydrodynamical simulations including radiative losses, of internal shocks in the outflows from classical novae, to explore the role of shocks in powering multi-wavelength emission from radio to gamma-ray wavelengths. Observations support a picture in which the initial phases of some novae generate a slow, equatorially-focused outflow (directly from the outer Lagrange point, or from a circumbinary disk), which then transitions to, or is overtaken by, a faster more isotropic outflow from the white dwarf which collides and shocks the slower flow, powering gamma-ray and optical emission through reprocessing by the ejecta. However, the common occurence of multiple peaks in nova light curves suggests that the outflow&#39;s acceleration need not be monotonic, but instead can involve successive transitions between &#34;fast&#34; and &#34;slow&#34; modes. Such a time-fluctuating outflow velocity naturally can reproduce several observed properties of nova, such as correlated gamma-ray and optical flares, expansion of the photosphere coincident with (though lagging slightly) the peak flare luminosity, and complex time-evolution of spectral lines (including accelerating, decelerating, and merging velocity components). While the shocks are still deeply embedded during the gamma-ray emission, the onset of $\sim$ keV X-ray and $\sim 10$ GHz radio synchrotron emission is typically delayed until the forward shock of the outermost monolithic shell (created by merger of multiple internal shock-generated shells) reaches a sufficiently low column through the dense external medium generated by the earliest phase of the outburst.

preprint2019arXiv

Kilonovae

The coalescence of double neutron star (NS-NS) and black hole (BH)-NS binaries are prime sources of gravitational waves (GW) for Advanced LIGO/Virgo and future ground-based detectors. Neutron-rich matter released from such events undergo rapid neutron capture (r-process) nucleosynthesis as it decompresses into space, enriching our universe with rare heavy elements like gold and platinum. Radioactive decay of these unstable nuclei powers a rapidly evolving, approximately isotropic thermal transient known as a ``kilonova&#39;, which probes the physical conditions during the merger and its aftermath. Here I review the history and physics of kilonovae, leading to the current paradigm of day-timescale emission at optical wavelengths from lanthanide-free components of the ejecta, followed by week-long emission with a spectral peak in the near-infrared (NIR). These theoretical predictions, as compiled in the original version of this review, were largely confirmed by the transient optical/NIR counterpart discovered to the first NS-NS merger, GW170817, discovered by LIGO/Virgo. Using a simple light curve model to illustrate the essential physical processes and their application to GW170817, I then introduce important variations about the standard picture which may be observable in future mergers. These include ~hours-long UV precursor emission, powered by the decay of free neutrons in the outermost ejecta layers or shock-heating of the ejecta by a delayed ultra-relativistic outflow; and enhancement of the luminosity from a long-lived central engine, such as an accreting BH or millisecond magnetar. Joint GW and kilonova observations of GW170817 and future events provide a new avenue to constrain the astrophysical origin of the r-process elements and the equation of state of dense nuclear matter.

preprint2007arXiv

Proto-Neutron Star Winds with Magnetic Fields and Rotation

We solve the 1D neutrino-heated non-relativistic MHD wind problem for conditions that range from slowly rotating (spin period P > 10 ms) protoneutron stars (PNSs) with surface field strengths typical of radio pulsars (B < 10^13 G), to &#34;proto-magnetars&#34; with B ~ 10^14-10^15 G in their hypothesized rapidly rotating initial states (P ~ 1 ms). We use the simulations of Bucciantini et al. (2006) to map our monopole results onto a more physical dipole geometry and to estimate the spindown of PNSs when their winds are relativistic. We then quantify the effects of rotation and magnetic fields on the mass loss, energy loss, and r-process nucleosynthesis in PNS winds. We describe the evolution of PNS winds through the Kelvin-Helmholtz cooling epoch, emphasizing the transition between (1) thermal neutrino-driven, (2) non-relativistic magnetically-dominated, and (3) relativistic magnetically-dominated outflows. We find that proto-magnetars with P ~ 1 ms and B > 10^15 G drive relativistic winds with luminosities, energies, and Lorentz factors (magnetization sigma ~ 0.1-1000) consistent with those required to produce long duration gamma-ray bursts and hyper-energetic supernovae (SNe). A significant fraction of the rotational energy may be extracted in only a few seconds, sufficiently rapidly to alter the asymptotic energy of the SN remnant, its morphology, and, potentially, its nucleosynthetic yield. Winds from PNSs with more modest rotation periods (2 - 10 ms) and with magnetar-strength fields produce conditions significantly more favorable for the r-process than winds from slowly rotating PNSs. Lastly, we show that energy and momentum deposition by convectively-excited waves further increase the likelihood of successful r-process in PNS winds.