Source author record

Brian D. Metzger

Brian D. Metzger appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

76works
13topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

76 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 "luminous" 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 "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.

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'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 ("Blue Ring Nebula") 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'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 "hidden" 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 "quasi-periodic eruptions" (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'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'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 ''red'' kilonovae and rule out ''blue'' kilonovae with $M>0.5 M_{\odot}$ (30% confidence). We strongly rule out all known types of short gamma-ray burst afterglows with viewing angles $<$17$^{\circ}$ assuming an initial jet opening angle of $\sim$$5.2^{\circ}$ and explosion energies and circumburst densities similar to afterglows explored in the literature. Finally, we explore the possibility that GW190814 merged in the disk of an active galactic nucleus, of which we find four in the localization region, but we do not find any candidate counterparts among these sources.

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 "Fast Radio Burst" 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 ("$r$-") 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 "tails" 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's acceleration need not be monotonic, but instead can involve successive transitions between "fast" and "slow" 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', 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.

preprint2016arXiv

A radio jet from the optical and x-ray bright stellar tidal disruption flare ASASSN-14li

The tidal disruption of a star by a supermassive black hole leads to a short-lived thermal flare. Despite extensive searches, radio follow-up observations of known thermal stellar tidal disruption flares (TDFs) have not yet produced a conclusive detection. We present a detection of variable radio emission from a thermal TDF, which we interpret as originating from a newly-launched jet. The multi-wavelength properties of the source present a natural analogy with accretion state changes of stellar mass black holes, suggesting all TDFs could be accompanied by a jet. In the rest frame of the TDF, our radio observations are an order of magnitude more sensitive than nearly all previous upper limits, explaining how these jets, if common, could thus far have escaped detection.

preprint2016arXiv

Assisted Inspirals of Stellar Mass Black Holes Embedded in AGN Disks: Solving the "Final AU Problem"

We explore the evolution of stellar mass black hole binaries (BHBs) which are formed in the self-gravitating disks of active galactic nuclei (AGN). Hardening due to three-body scattering and gaseous drag are effective mechanisms that reduce the semi-major axis of a BHB to radii where gravitational waves take over, on timescales shorter than the typical lifetime of the AGN disk. Taking observationally-motivated assumptions for the rate of star formation in AGN disks, we find a rate of disk-induced BHB mergers ($\mathcal{R} \sim 3~{\rm yr}^{-1}~{\rm Gpc}^{-3}$, but with large uncertainties) that is comparable with existing estimates of the field rate of BHB mergers, and the approximate BHB merger rate implied by the recent Advanced LIGO detection of GW150914. BHBs formed thorough this channel will frequently be associated with luminous AGN, which are relatively rare within the sky error regions of future gravitational wave detector arrays. This channel could also possess a (potentially transient) electromagnetic counterpart due to super-Eddington accretion onto the stellar mass black hole following the merger.

preprint2016arXiv

Binary Stellar Mergers with Marginally-Bound Ejecta: Excretion Disks, Inflated Envelopes, Outflows, and their Luminous Transients

We study mass loss from the outer Lagrange point (L2) in binary stellar mergers and their luminous transients by means of radiative hydrodynamical simulations. Previously, we showed that for binary mass ratios 0.06 < q < 0.8, synchronous L2 mass loss results in a radiatively inefficient, dust-forming unbound equatorial outflow. A similar outflow exists irrespective of q if the ratio of the sound speed to the orbital speed at the injection point is sufficiently large, ε= c_T/v_orb > 0.15. By contrast, for cold L2 mass-loss (ε< 0.15) from binaries with q < 0.06 or q > 0.8, the equatorial outflow instead remains marginally-bound and falls back to the binary over tens to hundreds of binary orbits, where it experiences additional tidal torqueing and shocking. As the bound gas becomes virialized with the binary, the luminosity of the system increases slowly at approximately constant photosphere radius, causing the temperature to rise. Subsequent evolution depends on the efficiency of radiative cooling. If the bound atmosphere is able to cool efficiently, as quantified by radiative diffusion time being shorter than the advection time (t_diff/t_adv < 1), then the virialized gas collapses to an excretion disk, while for t_diff/t_adv > 1 an isotropic wind is formed. Between these two extremes, an inflated envelope transports the heat generated near the binary to the surface by meridional flows. In all cases, the radiated luminosity reaches a fraction ~0.01 to 0.1 of Mdot v_orb^2/2, where Mdot is the mass outflow rate. We discuss the implications of our results for transients in the luminosity gap between classical novae and supernovae, such as V1309 Sco and V838 Mon.

preprint2016arXiv

Pair Fireball Precursors of Neutron Star Mergers

If at least one neutron star (NS) is magnetized in a binary NS merger, then the orbital motion of the conducting companion during the final inspiral induces a strong voltage and current along the magnetic field lines connecting the NSs. If a modest fraction, eta, of the extracted electromagnetic power extracted accelerates relativistic particles, the resulting gamma-ray emission a compact volume will result in the formation of an electron-positron pair fireball. Applying a steady-state pair wind model, we quantify the detectability of the precursor fireball with gamma-ray satellites. For eta ~ 1 the gamma-ray detection horizon of D_max ~ 10(Bd/1e14 G)^{3/4} Mpc is much closer than the Advanced LIGO/Virgo horizon of 200 Mpc, unless the NS surface magnetic field strength is very large, Bd > 1e15 G. Given the quasi-isotropic nature of the emission, mergers with weaker NS fields could contribute a nearby population of short gamma-ray bursts. Power not dissipated close to the binary is carried to infinity along the open field lines by a large scale Poynting flux. Reconnection within this outflow, well outside of the pair photosphere, provides a potential site for non-thermal emission, such as a coherent millisecond radio burst.

preprint2016arXiv

Production of the entire range of r-process nuclides by black hole accretion disc outflows from neutron star mergers

We consider $r$-process nucleosynthesis in outflows from black hole accretion discs formed in double neutron star and neutron star -- black hole mergers. These outflows, powered by angular momentum transport processes and nuclear recombination, represent an important -- and in some cases dominant -- contribution to the total mass ejected by the merger. Here we calculate the nucleosynthesis yields from disc outflows using thermodynamic trajectories from hydrodynamic simulations, coupled to a nuclear reaction network. We find that outflows produce a robust abundance pattern around the second $r$-process peak (mass number $A \sim 130$), independent of model parameters, with significant production of $A < 130$ nuclei. This implies that dynamical ejecta with high electron fraction may not be required to explain the observed abundances of $r$-process elements in metal poor stars. Disc outflows reach the third peak ($ A \sim 195$) in most of our simulations, although the amounts produced depend sensitively on the disc viscosity, initial mass or entropy of the torus, and nuclear physics inputs. Some of our models produce an abundance spike at $A = 132$ that is absent in the Solar system $r$-process distribution. The spike arises from convection in the disc and depends on the treatment of nuclear heating in the simulations. We conclude that disc outflows provide an important -- and perhaps dominant -- contribution to the $r$-process yields of compact binary mergers, and hence must be included when assessing the contribution of these systems to the inventory of $r$-process elements in the Galaxy.

preprint2016arXiv

Radio Constraints on Long-Lived Magnetar Remnants in Short Gamma-Ray Bursts

The merger of a neutron star binary may result in the formation of a rapidly-spinning magnetar. The magnetar can potentially survive for seconds or longer as a supramassive neutron star before collapsing to a black hole if, indeed, it collapses at all. During this process, a fraction of the magnetar's rotational energy of ~10^53 erg is transferred via magnetic spin-down to the surrounding ejecta. The resulting interaction between the ejecta and the surrounding circumburst medium powers a >year-long synchrotron radio transient. We present a search for radio emission with the Very Large Array following nine short-duration gamma-ray bursts (GRBs) at rest-frame times of ~1.3-7.6 years after the bursts, focusing on those events which exhibit early-time excess X-ray emission that may signify the presence of magnetars. We place upper limits of <18-32 microJy on the 6.0 GHz radio emission, corresponding to spectral luminosities of <(0.05-8.3)x10^39 erg s^-1. Comparing these limits to the predicted radio emission from a long-lived remnant and incorporating measurements of the circumburst densities from broad-band modeling of short GRB afterglows, we rule out a stable magnetar with an energy of 10^53 erg for half of the events in our sample. A supramassive remnant that injects a lower rotational energy of 10^52 erg is ruled out for a single event, GRB 050724A. This study represents the deepest and most extensive search for long-term radio emission following short GRBs to date, and thus the most stringent limits placed on the physical properties of magnetars associated with short GRBs from radio observations.

preprint2016arXiv

Rates of Stellar Tidal Disruption as Probes of the Supermassive Black Hole Mass Function

Rates of stellar tidal disruption events (TDEs) by supermassive black holes (SMBHs) due to two-body relaxation are calculated using a large galaxy sample (N=146) in order to explore the sensitivity of the TDE rates to observational uncertainties, such as the parametrization of galaxy light profiles and the stellar mass function. The largest uncertainty arises due to the poorly constrained occupation fraction of SMBHs in low-mass galaxies, which otherwise dominate the total TDE rate. The detection rate of TDE flares by optical surveys is calculated as a function of SMBH mass and other observables for several physically-motivated models of TDE emission. We also quantify the fraction of galaxies that produce deeply penetrating disruption events. If the majority of the detected events are characterized by super-Eddington luminosities (such as disk winds, or synchrotron radiation from an off-axis relativistic jet), then the measured SMBH mass distribution will tightly constrain the low-end SMBH occupation fraction. If Eddington-limited emission channels dominate, however, then the occupation fraction sensitivity is much less pronounced in a flux-limited survey (although still present in a volume-complete event sample). The SMBH mass distribution of the current sample of TDEs, though highly inhomogeneous and encumbered by selection effects, already suggests that Eddington-limited emission channels dominate. Even our most conservative rate estimates appear to be in tension with much lower observationally inferred TDE rates, and we discuss several possible resolutions to this discrepancy.

preprint2016arXiv

Supernovae powered by magnetars that transform into black holes

Rapidly rotating, strongly magnetized neutron stars (magnetars) can release their enormous rotational energy via magnetic spin-down, providing a power source for bright transients such as superluminous supernovae. On the other hand, particularly massive (so-called supramassive) neutron stars require a minimum rotation rate to support their mass against gravitational collapse, below which the neutron star collapses to a black hole. We model the light curves of supernovae powered by magnetars which transform into black holes. Although the peak luminosities can reach high values in the range of superluminous supernovae, their post maximum light curves can decline very rapidly because of the sudden loss of the central energy input. Early black hole transformation also enhances the shock breakout signal from the magnetar-driven bubble relative to the main supernova peak. Our synthetic light curves of supernovae powered by magnetars transforming to black holes are consistent with those of some rapidly evolving bright transients recently reported by Arcavi et al. (2016).

preprint2016arXiv

Time dependent models of accretion disks with nuclear burning following the tidal disruption of a white dwarf by a neutron star

We construct time-dependent one-dimensional (vertically averaged) models of accretion disks produced by the tidal disruption of a white dwarf (WD) by a binary neutron star (NS) companion. Nuclear reactions in the disk midplane burn the WD matter to increasingly heavier elements at sequentially smaller radii, releasing substantial energy which can impact the disk dynamics. A model for disk outflows is employed, by which cooling from the outflow balances other sources of heating (viscous, nuclear) in regulating the Bernoulli parameter of the midplane to a fixed value $\lesssim 0$. We perform a comprehensive parameter study of the compositional yields and velocity distributions of the disk outflows for WDs of different initial compositions. For C/O WDs, the radial composition profile of the disk evolves self-similarly in a quasi-steady-state manner, and is remarkably robust to model parameters. The nucleosynthesis in helium WD disks does not exhibit this behavior, which instead depends sensitively on factors controlling the disk midplane density (e.g. the strength of the viscosity, $α$). By the end of the simulation, a substantial fraction of the WD mass is unbound in outflows at characteristic velocities of $\sim 10^9~{\rm cm~s}^{-1}$. The outflows from WD-NS merger disks contain $10^{-4}-3 \times 10^{-3} M_\odot$ of radioactive $^{56}$Ni, resulting in fast ($\sim$ week long) dim ($\sim10^{40}~{\rm erg~s}^{-1}$) optical transients; shock heating of the ejecta by late time outflows may increase the peak luminosity to $\sim 10^{43}~{\rm erg~s}^{-1}$. The accreted mass onto the neutron star is probably not sufficient to induce gravitational collapse, but may be capable of spinning up the NS to periods of $\sim 10~{\rm ms}$. This is a new possible channel for forming isolated recycled pulsars.

preprint2015arXiv

A Bright Year for Tidal Disruptions?

When a star is tidally disrupted by a supermassive black hole (BH), roughly half of its mass falls back to the BH at super-Eddington rates. Being tenuously gravitationally bound and unable to cool radiatively, only a small fraction f_in << 1 of the returning debris will likely be incorporated into the disk and accrete, with the vast majority instead becoming unbound in an outflow of velocity ~1e4 km/s. This slow outflow spreads laterally, encasing the BH. For months or longer, the outflow remains sufficiently neutral to block hard EUV and X-ray radiation from the hot inner disk, which instead becomes trapped in a radiation-dominated nebula. Ionizing nebular radiation heats the inner edge of the ejecta to temperatures of T > few 1e4 K, converting the emission to optical/near-UV wavelengths where photons more readily escape due to the lower opacity. This can explain the unexpectedly low and temporally constant effective temperatures of optically-discovered TDE flares. For BHs with relatively high masses M_BH > 1e7 M_sun the ejecta can become ionized at an earlier stage, or for a wider range of viewing angles, producing a TDE flare which is instead dominated by thermal X-ray emission. We predict total radiated energies consistent with those of observed TDE flares, and ejecta velocities that agree with the measured emission line widths. The peak optical luminosity for M_BH < 1e6 M_sun is suppressed due to adiabatic losses in the inner disk wind, possibly contributing to the unexpected dearth of optical TDEs in galaxies with low mass BHs. In the classical picture, where f_in ~ 1, TDEs de-spin supermassive BHs and cap their maximum spins well below theoretical accretion physics limits. This cap is greatly relaxed in our model, and existing Fe K-alpha spin measurements provide suggestive preliminary evidence that f_in < 1.

preprint2015arXiv

Circumnuclear Media of Quiescent Supermassive Black Holes

We calculate steady-state, one-dimensional hydrodynamic profiles of hot gas in slowly accreting ("quiescent") galactic nuclei for a range of central black hole masses $M_{\bullet}$, parametrized gas heating rates, and observationally-motivated stellar density profiles. Mass is supplied to the circumnuclear medium by stellar winds, while energy is injected primarily by stellar winds, supernovae, and black hole feedback. Analytic estimates are derived for the stagnation radius (where the radial velocity of the gas passes through zero) and the large scale gas inflow rate, $\dot{M}$, as a function of $M_{\bullet}$ and the gas heating efficiency, the latter being related to the star-formation history. We assess the conditions under which radiative instabilities develop in the hydrostatic region near the stagnation radius, both in the case of a single burst of star formation and for the average star formation history predicted by cosmological simulations. By combining a sample of measured nuclear X-ray luminosities, $L_x$, of nearby quiescent galactic nuclei with our results for $\dot{M}(M_{\bullet})$ we address whether the nuclei are consistent with accreting in a steady-state, thermally-stable manner for radiative efficiencies predicted for radiatively inefficiency accretion flows. We find thermally-stable accretion cannot explain the short average growth times of low mass black holes in the local Universe, which must instead result from gas being fed in from large radii, due either to gas inflows or thermal instabilities acting on larger, galactic scales. Our results have implications for attempts to constrain the occupation fraction of SMBHs in low mass galaxies using the mean $L_x-M_{\bullet}$ correlation, as well as the predicted diversity of the circumnuclear densities encountered by relativistic outflows from tidal disruption events.

preprint2015arXiv

Cool and Luminous Transients from Mass-Losing Binary Stars

We study transients produced by equatorial disk-like outflows from catastrophically mass-losing binary stars with an asymptotic velocity and energy deposition rate near the inner edge which are proportional to the binary escape velocity v_esc. As a test case, we present the first smoothed-particle radiation-hydrodynamics calculations of the mass loss from the outer Lagrange point with realistic equation of state and opacities. The resulting spiral stream becomes unbound for binary mass ratios 0.06 < q < 0.8. For synchronous binaries with non-degenerate components, the spiral-stream arms merge at a radius of ~10a, where a is the binary semi-major axis, and the accompanying shock thermalizes about 10% of the kinetic power of the outflow. The mass-losing binary outflows produce luminosities reaching up to ~10^6 L_Sun and effective temperatures spanning 500 < T_eff < 6000 K, which is compatible with many of the class of recently-discovered red transients such as V838 Mon and V1309 Sco. Dust readily forms in the outflow, potentially in a catastrophic global cooling transition. The appearance of the transient is viewing angle-dependent due to vastly different optical depths parallel and perpendicular to the binary plane. We predict a correlation between the peak luminosity and the outflow velocity, which is roughly obeyed by the known red transients. Outflows from mass-losing binaries can produce luminous (10^5 L_Sun) and cool (T_eff < 1500 K) transients lasting a year or longer, as has potentially been detected by Spitzer surveys of nearby galaxies.

preprint2015arXiv

Does the Collapse of a Supramassive Neutron Star Leave a Debris Disk?

One possible channel for black hole formation is the collapse of a rigidly rotating massive neutron star as it loses its angular momentum or gains excessive mass through accretion. It was proposed that part of the neutron star may form a debris disk around the black hole. Such short-lived massive disks could be the sources of powerful jets emitting cosmological gamma-ray bursts. Whether the collapse creates a disk depends on the equation of state of the neutron star. We survey a wide range of equations of states allowed by observations and find that disk formation is unfeasible. We conclude that this channel of black hole formation is incapable of producing powerful jets, and discuss implications for models of gamma-ray bursts.

preprint2015arXiv

Electromagnetic Signatures of Neutron Star Mergers in the Advanced LIGO Era

The mergers of binaries containing neutron stars and stellar-mass black holes are the most promising sources for direct detection in gravitational waves by the interferometers Advanced LIGO and Virgo over the next few years. The concurrent detection of electromagnetic emission from these events would greatly enhance the scientific return of these discoveries. Here we review the state of the art in modeling the electromagnetic signal of neutron star binary mergers across different phases of the merger and multiple wavelengths. We focus on those observables which provide the most sensitive diagnostics of the merger physics and the contribution to the synthesis of rapid neutron capture ($r$-process) elements in the Galaxy. We also outline expected future developments on the observational and theoretical sides of this rapidly evolving field.

preprint2015arXiv

Extragalactic Transients in the Era of Wide-Field Radio Surveys. I. Detection Rates and Light Curve Characteristics

The impending era of wide-field radio surveys has the potential to revolutionize our understanding of astrophysical transients. Here we evaluate the prospects of a wide range of planned and hypothetical radio surveys using the properties and volumetric rates of known and hypothetical classes of extragalactic synchrotron radio transients (e.g., on- and off-axis gamma-ray bursts [GRB], supernovae, tidal disruption events [TDE], compact object mergers). Utilizing these sources and physically motivated considerations we assess the allowed phase-space of radio luminosity and peak timescale for extragalactic transients. We also include for the first time effects such as redshift evolution of the rates, K-corrections, and non-Euclidean luminosity distance, which affect the detection rates of the most sensitive surveys. The number of detected events is calculated by means of a Monte Carlo method, using the various survey properties (depth, cadence, area) and realistic detection criteria that include a cut on the minimum variability of the transients during the survey and an assessment of host galaxy contamination. Near-term GHz frequency surveys (ASKAP/VAST, Very Large Array Sky Survey) will detect few events: <~30-50 on- and off-axis long GRBs and off-axis tidal disruption events, and ~10-20 neutron star binary mergers if ~1% of the mergers result in a stable millisecond magnetar. Low-frequency surveys (e.g., LOFAR) are unlikely to detect any transients, while a hypothetical large-scale mm survey may detect ~40 on-axis long GRBs. On the other hand, SKA surveys at ~0.1-1 GHz have the potential to uncover thousands of transients, mainly on- and off-axis long GRBs, on-axis short GRBs, off-axis TDEs, and neutron star binary mergers with magnetar remnants.

preprint2015arXiv

Gamma-ray novae as probes of relativistic particle acceleration at non-relativistic shocks

The Fermi LAT discovery that classical novae produce >100 MeV gamma-rays establishes that shocks and relativistic particle acceleration are key features of these events. These shocks are likely to be radiative due to the high densities of the nova ejecta at early times coincident with the gamma-ray emission. Thermal X-rays radiated behind the shock are absorbed by neutral gas and reprocessed into optical emission, similar to Type IIn (interacting) supernovae. Gamma-rays are produced by collisions between relativistic protons with the nova ejecta (hadronic scenario) or Inverse Compton/bremsstrahlung emission from relativistic electrons (leptonic scenario), where in both scenarios the efficiency for converting relativistic particle energy into LAT gamma-rays is at most a few tens of per cent. The ratio of gamma-ray and optical luminosities, L_gam/L_opt, thus sets a lower limit on the fraction of the shock power used to accelerate relativistic particles, e_nth. The measured values of L_gam/L_opt for two classical novae, V1324 Sco and V339 Del, constrains e_nth > 1e-2 and > 1e-3, respectively. Inverse Compton models for the gamma-ray emission are disfavored given the low electron acceleration efficiency, e_nth ~ 1e-4-1e-3, inferred from observations of Galactic cosmic rays and particle-in-cell (PIC) numerical simulations. A fraction > 100(0.01/e_nth) and > 10(0.01/e_nth) per cent of the optical luminosity is powered by shocks in V1324 Sco and V339 Del, respectively. Such high fractions challenge standard models that instead attribute all nova optical emission to the direct outwards transport of thermal energy released near the white dwarf surface.

preprint2015arXiv

Magnetar Driven Shock Breakout and Double Peaked Supernova Light Curves

The light curves of some luminous supernovae are suspected to be powered by the spindown energy of a rapidly rotating magnetar. Here we describe a possible signature of the central engine: a burst of shock breakout emission occurring several days after the supernova explosion. The energy input from the magnetar inflates a high-pressure bubble that drives a shock through the pre-exploded supernova ejecta. If the magnetar is powerful enough, that shock will near the ejecta surface and become radiative. At the time of shock breakout, the ejecta will have expanded to a large radius (~10^{14} cm) so that the radiation released is at optical/ultraviolet wavelengths (T ~ 20,000 K) and lasts for several days. The luminosity and timescale of this magnetar driven shock breakout are similar to the first peak observed recently in the double-peaked light curve of SNLSQ14BDQ. However, for a large region of model parameter space, the breakout emission is predicted to be dimmer than the diffusive luminosity from direct magnetar heating. A distinct double peaked light curve may therefore only be conspicuous if thermal heating from the magnetar is suppressed at early times. We describe how such a delay in heating may naturally result from inefficient dissipation and thermalization of the pulsar wind magnetic energy. Without such suppression, the breakout may only be noticeable as a small bump or kink in the early luminosity or color evolution, or as a small but abrupt rise in the photospheric velocity. A similar breakout signature may accompany other central engines in supernovae, such as a black hole accreting fallback material.

preprint2015arXiv

Non-Thermal Radio Emission from Colliding Flows in Classical Nova V1723 Aql

The importance of shocks in nova explosions has been highlighted by Fermi's discovery of γ-ray producing novae. Over three years of multi-band VLA radio observations of the 2010 nova V1723 Aql show that shocks between fast and slow flows within the ejecta led to the acceleration of particles and the production of synchrotron radiation. Soon after the start of the eruption, shocks in the ejecta produced an unexpected radio flare, resulting in a multi-peaked radio light curve. The emission eventually became consistent with an expanding thermal remnant with mass $2 \times 10^{-4} M_\odot$ and temperature $10^4$ K. However, during the first two months, the $\gtrsim 10^6$ K brightness temperature at low frequencies was too high to be due to thermal emission from the small amount of X-ray producing shock-heated gas. Radio imaging showed structures with velocities of 400 km s$^{-1}$ (d/6 kpc) in the plane of the sky, perpendicular to a more elongated 1500 km s$^{-1}$ (d/6 kpc) flow. The morpho-kinematic structure of the ejecta from V1723 Aql appears similar to nova V959 Mon, where collisions between a slow torus and a faster flow collimated the fast flow and gave rise to γ-ray producing shocks. Optical spectroscopy and X-ray observations of V1723 Aql during the radio flare are consistent with this picture. Our observations support the idea that shocks in novae occur when a fast flow collides with a slow collimating torus. Such shocks could be responsible for hard X-ray emission, γ-ray production, and double-peaked radio light curves from some classical novae.

preprint2015arXiv

Novae as Tevatrons: Prospects for CTA and IceCube

The discovery of novae as sources of ~GeV gamma-rays highlights the key role of shocks and relativistic particle acceleration in these transient systems. Although there is evidence for a spectral cut-off above energies ~1-100 GeV at particular epochs in some novae, the maximum particle energy achieved in these accelerators has remained an open question. The high densities of the nova ejecta (~10 orders of magnitude larger than in supernova remnants) render the gas far upstream of the shock neutral and shielded from ionizing radiation. The amplification of the magnetic field needed for diffusive shock acceleration requires ionized gas, thus confining the acceleration process to a narrow photo-ionized layer immediately ahead of the shock. Based on the growth rate in this layer of the hybrid non-resonant cosmic ray current-driven instability (considering also ion-neutral damping), we quantify the maximum particle energy, Emax, across the range of shock velocities and upstream densities of interest. We find values of Emax ~ 10 GeV - 10 TeV, which are broadly consistent with the inferred spectral cut-offs, but which could also in principle lead to emission extending to higher energies >100 GeV accessible to atmosphere Cherenkov telescopes, such as the planned Cherenkov Telescope Array (CTA). Detecting TeV neutrinos with IceCube in hadronic scenarios appears to be more challenging, although the prospects are improved if the shock power during the earliest, densest phases of the nova outburst is higher than is implied by the observed GeV light curves, due to downscattering of the gamma-rays by electrons within the ejecta. Novae provide ideal nearby laboratories to study magnetic field amplification and the onset of cosmic ray acceleration, because other time-dependent sources (e.g. radio supernovae) typically occur too distant to detect as gamma-ray sources.

preprint2015arXiv

The Diversity of Transients from Magnetar Birth

Strongly-magnetized, rapidly-rotating neutron stars are contenders for the central engines of both long-duration gamma-ray bursts (LGRBs) and hydrogen-poor super-luminous supernovae (SLSNe-I). Models for typical (~minute long) LGRBs invoke magnetars with high dipole magnetic fields (Bd > 1e15 G) and short spin-down times, while models for SLSNe-I invoke neutron stars with weaker fields and longer spin-down times of weeks. Here we identify a transition region in the space of Bd and birth period for which a magnetar can power both a long GRB and a luminous SN. In particular, we show that a 2 ms period magnetar with a spin-down time of ~1e4 s can explain the observations of both the ultra-long GRB 111209 and its associated luminous SN2011kl. For magnetars with longer spin down times, we predict even longer duration (~1e6 s) GRBs and brighter supernovae, a correlation that extends to Swift J2058+05 (commonly interpreted as a tidal disruption event). We further show that previous estimates of the maximum rotational energy of a proto-magnetar were too conservative and energies up to Emax ~1-2e53 erg are possible. The magnetar model can therefore comfortably accommodate the extreme energy requirements recently posed by the most luminous supernova ASASSN-15lh. The high ionization flux from a pulsar wind nebula powering ASASSN-15lh may lead to an "ionization break-out" X-ray burst over the coming months, which would be accompanied by an abrupt change in the optical spectrum. We conclude by briefly contrasting millisecond magnetar and black hole models for SLSNe and ultra-long GRBs.

preprint2014arXiv

Neutrino-heated winds from rotating proto-magnetars

We calculate the steady-state properties of neutrino-driven winds from strongly magnetized, rotating proto-neutron stars (`proto-magnetars') under the assumption that the outflow geometry is set by the force-free magnetic field of an aligned dipole. Our goal is to assess proto-magnetars as sites of r-process nucleosynthesis and gamma-ray burst engines. One dimensional solutions calculated along flux tubes corresponding to different polar field lines are stitched together to determine the global properties of the flow at a given neutrino luminosity and rotation period. Proto-magnetars with rotation periods of P~2-5 ms are shown to produce outflows more favorable for the production of third-peak r-process nuclei due to their much shorter expansion times through the seed nucleus formation region, yet only moderately lower entropies, as compared to normal spherical PNS winds. Proto-magnetars with moderately rapid birth periods P~3-5 ms may thus represent a promising Galactic r-process site which is compatible with a variety of other observations, including the recent discovery of possible magnetar-powered supernovae in metal poor galaxies. We also confirm previous results that the outflows from proto-magnetars with P~1-2 ms can achieve maximum Lorentz factors Gamma ~ 100-1000 in the range necessary to power gamma-ray bursts (GRBs). The implications of GRB jets with a heavy nuclei-dominated composition as sources of ultra-high energy cosmic rays are also addressed.

preprint2014arXiv

Neutron-powered precursors of kilonovae

The merger of binary neutron stars (NSs) ejects a small quantity of neutron rich matter, the radioactive decay of which powers a day to week long thermal transient known as a kilonova. Most of the ejecta remains sufficiently dense during its expansion that all neutrons are captured into nuclei during the r-process. However, recent general relativistic merger simulations by Bauswein and collaborators show that a small fraction of the ejected mass (a few per cent, or ~1e-4 Msun) expands sufficiently rapidly for most neutrons to avoid capture. This matter originates from the shocked-heated interface between the merging NSs. Here we show that the beta-decay of these free neutrons in the outermost ejecta powers a `precursor' to the main kilonova emission, which peaks on a timescale of a few hours following merger at U-band magnitude ~22 (for an assumed distance of 200 Mpc). The high luminosity and blue colors of the neutron precursor render it a potentially important counterpart to the gravitational wave source, that may encode valuable information on the properties of the merging binary (e.g. NS-NS versus NS-black hole) and the NS equation of state. Future work is necessary to assess the robustness of the fast moving ejecta and the survival of free neutrons in the face of neutrino absorptions, although the precursor properties are robust to a moderate amount of leptonization. Our results provide additional motivation for short latency gravitational wave triggers and rapid follow-up searches with sensitive ground based telescopes.

preprint2014arXiv

Outflows from accretion disks formed in neutron star mergers: effect of black hole spin

The accretion disk that forms after a neutron star merger is a source of neutron-rich ejecta. The ejected material contributes to a radioactively-powered electromagnetic transient, with properties that depend sensitively on the composition of the outflow. Here we investigate how the spin of the black hole remnant influences mass ejection on the thermal and viscous timescales. We carry out two-dimensional, time-dependent hydrodynamic simulations of merger remnant accretion disks including viscous angular momentum transport and approximate neutrino self-irradiation. The gravity of the spinning black hole is included via a pseudo-Newtonian potential. We find that a disk around a spinning black hole ejects more mass, up to a factor of several, relative to the non-spinning case. The enhanced mass loss is due to energy release by accretion occurring deeper in the gravitational potential, raising the disk temperature and hence the rate of viscous heating in regions where neutrino cooling is ineffective. The mean electron fraction of the outflow increases moderately with BH spin due to a highly-irradiated (though not neutrino-driven) wind component. While the bulk of the ejecta is still very neutron-rich, thus generating heavy r-process elements, the leading edge of the wind contains a small amount of Lanthanide-free material. This component can give rise to a ~1 day blue optical `bump' in a kilonova light curve, even in the case of prompt BH formation, which may facilitate its detection.

preprint2014arXiv

Red or blue? A potential kilonova imprint of the delay until black hole formation following a neutron star merger

Mergers of binary neutron stars (NSs) usually result in the formation of a hypermassive neutron star (HMNS). Whether- and when this remnant collapses to a black hole (BH) depends primarily on the equation of state and on angular momentum transport processes, both of which are uncertain. Here we show that the lifetime of the merger remnant may be directly imprinted in the radioactively powered kilonova emission following the merger. We employ axisymmetric, time-dependent hydrodynamic simulations of remnant accretion disks orbiting a HMNS of variable lifetime, and characterize the effect of this delay to BH formation on the disk wind ejecta. When BH formation is relatively prompt (~ 100 ms), outflows from the disk are sufficiently neutron rich to form heavy r-process elements, resulting in ~ week-long emission with a spectral peak in the near-infrared (NIR), similar to that produced by the dynamical ejecta. In contrast, delayed BH formation allows neutrinos from the HMNS to raise the electron fraction in the polar direction to values such that potentially Lanthanide-free outflows are generated. The lower opacity would produce a brighter, bluer, and shorter-lived ~ day-long emission (a `blue bump') prior to the late NIR peak from the dynamical ejecta and equatorial wind. This new diagnostic of BH formation should be useful for events with a signal to noise lower than that required for direct detection of gravitational waveform signatures.

preprint2014arXiv

Shocks in nova outflows. I. Thermal emission

Evidence for shocks in nova outflows include (1) multiple velocity components in the optical spectra; (2) keV X-ray emission weeks to months after the outburst; (3) early radio flare on timescales of months, in excess of that predicted from the freely expanding photo-ionized gas; and (4) ~ GeV gamma-rays. We present a 1D model for the shock interaction between the fast nova outflow and a dense external shell (DES) and its associated thermal X-ray, optical, and radio emission. The forward shock is radiative initially when the density of shocked gas is highest, at which times radio emission originates from the dense cooling layer immediately downstream of the shock. The radio light curve is characterized by sharper rises to maximum and later peak times at progressively lower frequencies, with a peak brightness temperature that is approximately independent of frequency. We apply our model to the recent gamma-ray classical nova V1324 Sco, obtaining an adequate fit to the early radio maximum for reasonable assumptions about the fast nova outflow and assuming the DES possesses a velocity ~1e3 km/s and mass ~ 2e-4 M_sun; the former is consistent with the velocities of narrow line absorption systems observed previously in nova spectra, while the total ejecta mass of the DES and fast outflow is consistent with that inferred independently by modeling the late radio peak. Rapid evolution of the early radio light curves require the DES possess a steep outer density profile, which may indicate that the onset of mass loss from the white dwarf was rapid, providing indirect evidence that the DES was expelled by the thermonuclear runaway event. Reprocessed X-rays from the shock absorbed by the DES at early times may contribute significantly to the optical/UV emission, which we speculate is responsible for the previously unexplained `plateaus' and secondary maxima in nova optical light curves.

preprint2013arXiv

Constraints on long-lived remnants of neutron star binary mergers from late-time radio observations of short duration gamma-ray bursts

The coalescence of a binary neutron star (NS) system may in some cases produce a massive NS remnant that is long-lived and, potentially, indefinitely stable to gravitational collapse. Such a remnant has been proposed as an explanation for the late X-ray emission observed following some short duration gamma-ray bursts (GRBs) and as possible electromagnetic counterparts to the gravitational wave chirp. A stable NS merger remnant necessarily possesses a large rotational energy > 1e52 erg, the majority of which is ultimately deposited into the surrounding circumburst medium (CBM) at mildly relativistic velocities. We present Very Large Array radio observations of 7 short GRBs, some of which possessed temporally extended X-ray emission, on timescales of ~1-3 years following the initial burst. No radio sources were detected, with typical upper limits ~0.3 mJy at 1.4 GHz. A basic model for the synchrotron emission from the blast wave is used to constrain the presence of a long-lived NS merger remnant in each system. Depending on the GRB, our non-detections translate into upper limits on the CBM density n < 3e-2 - 3e-3 particles/cm^3 required for consistency with the remnant hypothesis. Our upper limits rule out a long-lived remnant in GRB 050724 and 060505, but cannot rule out such a remnant in other systems due to their lower inferred CBM densities based on afterglow modeling or the lack of such constraints. Additional VLA observations in the near future could place tighter limits on the presence of merger remnants in these system. The population of long-lived NS merger remnants will also be constrained by their (non-)detection with upcoming radio transient surveys.

preprint2013arXiv

Delayed outflows from black hole accretion tori following neutron star binary coalescence

Expulsion of neutron-rich matter following the merger of neutron star (NS) binaries is crucial to the radioactively-powered electromagnetic counterparts of these events and to their relevance as sources of r-process nucleosynthesis. Here we explore the long-term (viscous) evolution of remnant black hole accretion disks formed in such mergers by means of two-dimensional, time-dependent hydrodynamical simulations. The evolution of the electron fraction due to charged-current weak interactions is included, and neutrino self-irradiation is modeled as a lightbulb that accounts for the disk geometry and moderate optical depth effects. Over several viscous times (~1s), a fraction ~10% of the initial disk mass is ejected as a moderately neutron-rich wind (Y_e ~ 0.2) powered by viscous heating and nuclear recombination, with neutrino self-irradiation playing a sub-dominant role. Although the properties of the outflow vary in time and direction, their mean values in the heavy-element production region are relatively robust to variations in the initial conditions of the disk and the magnitude of its viscosity. The outflow is sufficiently neutron-rich that most of the ejecta forms heavy r-process elements with mass number A >130, thus representing a new astrophysical source of r-process nucleosynthesis, distinct from that produced in the dynamical ejecta. Due to its moderately high entropy, disk outflows contain a small residual fraction ~1% of helium, which could produce a unique spectroscopic signature.

preprint2013arXiv

Ionization Break-Out from Millisecond Pulsar Wind Nebulae: an X-ray Probe of the Origin of Superluminous Supernovae

Magnetic spin-down of a millisecond neutron star has been proposed as the power source of hydrogen-poor "superluminous" supernovae (SLSNe-I). However, producing an unambiguous test that can distinguish this model from alternatives, such as circumstellar interaction, has proven challenging. After the supernova explosion, the pulsar wind inflates a hot cavity behind the expanding stellar ejecta: the nascent millisecond pulsar wind nebula. Electron/positron pairs injected by the wind cool through inverse Compton scattering and synchrotron emission, producing a pair cascade and hard X-ray spectrum inside the nebula. These X-rays ionize the inner exposed side of the ejecta, driving an ionization front that propagates outwards with time. Under some conditions this front can breach the ejecta surface within months after the optical supernova peak, allowing ~0.1-1 keV photons to escape the nebula unattenuated with a characteristic luminosity L_X ~ 1e43-1e45 erg/s. This "ionization break-out" may explain the luminous X-ray emission observed from the transient SCP 06F, providing direct evidence that this SLSN was indeed engine-powered. Luminous break-out requires a low ejecta mass and that the spin-down time of the pulsar be comparable to the photon diffusion timescale at optical maximum, the latter condition similar to that required for a supernova with a high optical fluence. These relatively special requirements may explain why most SLSNe-I are not accompanied by detectable X-ray emission. Global asymmetry of the supernova ejecta increases the likelihood of an early break-out along the direction of lowest density. Atomic states with lower threshold energies are more readily ionized at earlier times near optical maximum, allowing UV break-out across a wider range of pulsar and ejecta properties than X-ray break-out, possibly contributing to the blue/UV colors of SLSNe-I.

preprint2013arXiv

Optical and X-ray emission from stable millisecond magnetars formed from the merger of binary neutron stars

The coalescence of binary neutron stars (NSs) may in some cases produce a stable massive NS remnant rather than a black hole. Due to the substantial angular momentum from the binary, such a remnant is born rapidly rotating and likely acquires a strong magnetic field (a `millisecond magnetar'). Magnetic spin-down deposits a large fraction of the rotational energy from the magnetar behind the small quantity of mass ejected during the merger. This has the potential for creating a bright transient that could be useful for determining whether a NS or black hole was formed in the merger. We investigate the expected signature of such an event, including for the first time the important impact of electron/positron pairs injected by the millisecond magnetar into the surrounding nebula. These pairs cool via synchrotron and inverse Compton emission, producing a pair cascade and hard X-ray spectrum. A fraction of these X-rays are absorbed by the ejecta walls and re-emitted as thermal radiation, leading to an optical/UV transient peaking at a luminosity of ~1e43-1e44 erg/s on a timescale of several hours to days. This is dimmer than predicted by simpler analytic models because the large optical depth of electron/positron pairs across the nebula suppresses the efficiency with which the magnetar spin down luminosity is thermalized. Nevertheless, the optical/UV emission is more than two orders of magnitude brighter than a radioactively powered `kilonova.' In some cases nebular X-rays are sufficiently luminous to re-ionize the ejecta, in which case non-thermal X-rays escape the ejecta unattenuated with a similar peak luminosity and timescale as the optical radiation. We discuss the implications of our results for the temporally extended X-ray emission that is observed to follow some short gamma-ray bursts (GRBs), including the kilonova candidates GRB 080503 and GRB 130603B.

preprint2013arXiv

Short GRB 130603B: Discovery of a jet break in the optical and radio afterglows, and a mysterious late-time X-ray excess

We present radio, optical/NIR, and X-ray observations of the afterglow of the short-duration 130603B, and uncover a break in the radio and optical bands at 0.5 d after the burst, best explained as a jet break with an inferred jet opening angle of 4-8 deg. GRB 130603B is only the third short GRB with a radio afterglow detection to date, and the first time that a jet break is evident in the radio band. We model the temporal evolution of the spectral energy distribution to determine the burst explosion properties and find an isotropic-equivalent kinetic energy of (0.6-1.7) x 10^51 erg and a circumburst density of 5 x 10^-3-30 cm^-3. From the inferred opening angle of GRB 130603B, we calculate beaming-corrected energies of Egamma (0.5-2) x 10^49 erg and EK (0.1-1.6) x 10^49 erg. Along with previous measurements and lower limits we find a median short GRB opening angle of 10 deg. Using the all-sky observed rate of 10 Gpc^-3 yr^-1, this implies a true short GRB rate of 20 yr^-1 within 200 Mpc, the Advanced LIGO/VIRGO sensitivity range for neutron star binary mergers. Finally, we uncover evidence for significant excess emission in the X-ray afterglow of GRB 130603B at >1 d and conclude that the additional energy component could be due to fall-back accretion or spin-down energy from a magnetar formed following the merger.

preprint2013arXiv

Swift J1644+57 gone MAD: the case for dynamically-important magnetic flux threading the black hole in a jetted tidal disruption event

The unusual transient Swift J1644+57 likely resulted from a collimated relativistic jet powered by accretion onto a massive black hole (BH) following the tidal disruption (TD) of a star. Several mysteries cloud the interpretation of this event: (1) extreme flaring and `plateau' shape of the X-ray/gamma-ray light curve during the first 10 days after the gamma-ray trigger; (2) unexpected rebrightening of the forward shock radio emission months after trigger; (3) no obvious evidence for jet precession, despite misalignment typically expected between the angular momentum of the accretion disk and BH; (4) recent abrupt shut-off in jet X-ray emission after 1.5 years. Here we show that all of these seemingly disparate mysteries are naturally resolved by one assumption: the presence of strong magnetic flux Phi threading the BH. Initially, Phi is weak relative to high fall-back mass accretion rate, Mdot, and the disk and jets precess about the BH axis = our line of sight. As Mdot drops, Phi becomes dynamically important and leads to a magnetically-arrested disk (MAD). MAD naturally aligns disk and jet axis along the BH spin axis, but only after a violent rearrangement phase (jet wobbling). This explains the erratic light curve at early times and the lack of precession at later times. We use our model for Swift J1644+57 to constrain BH and disrupted star properties, finding that a solar-mass main sequence star disrupted by a relatively low mass, M~10^5-10^6 Msun, BH is consistent with the data, while a WD disruption (though still possible) is disfavored. The magnetic flux required to power Swift J1644+57 is too large to be supplied by the star itself, but it could be collected from a quiescent `fossil' accretion disk present in the galactic nucleus prior to the TD. The presence (lack of) of such a fossil disk could be a deciding factor in what TD events are accompanied by powerful jets.[abridged]

preprint2013arXiv

X-ray Decay Lines from Heavy Nuclei in Supernova Remnants as a Probe of the r-Process Origin and the Birth Periods of Magnetars

The origin of rapid neutron capture (r-process) nuclei remains one of the longest standing mysteries in nuclear astrophysics. Core collapse supernovae (SNe) and neutron star binary mergers are likely r-process sites, but little evidence yet exists for their in situ formation in such environments. Motivated by the advent of sensitive new or planned X-ray telescopes such as the Nuclear Spectroscopic Telescope Array (NuSTAR) and the Large Observatory for X-ray Timing (LOFT), we revisit the prospects for the detection of X-ray decay lines from r-process nuclei in young or nearby supernova remnants. For all remnants planned to be observed by NuSTAR (and several others), we conclude that r-process nuclei are detectable only if the remnant possesses a large overabundance O > 1e3 relative to the average yield per SN. Prospects are better for the next Galactic SN (assumed age of 3 years and distance of 10 kpc), for which an average r-process yield is detectable via the 10.7(9.2) keV line complexes of Os194 by LOFT at 6 sigma (5 sigma) confidence; the 27.3 keV line complex of Sb125 is detectable by NuSTAR at 2 sigma for O > 2. We also consider X-rays lines from the remnants of Galactic magnetars, motivated by the much higher r-process yields of the magneto-rotationally driven SNe predicted to birth magnetars. The ~ 3.6-3.9 keV lines of Sn126 are potentially detectable in the remnants of the magnetars 1E1547.0-5408 and 1E2259+586 by LOFT for an assumed r-process yield predicted by recent simulations. The (non-)detection of these lines can thus probe whether magnetars are indeed born with millisecond periods. Finally, we consider a blind survey of the Galactic plane with LOFT for r-process lines from the most recent binary neutron star merger remnant, concluding that a detection is unlikely without additional information on the merger location.

preprint2012arXiv

Afterglow Model for the Radio Emission from the Jetted Tidal Disruption Candidate Swift J1644+57

The recent transient event Swift J1644+57 has been interpreted as emission from a collimated relativistic jet, powered by the sudden onset of accretion onto a supermassive black hole following the tidal disruption of a star. Here we model the radio-microwave emission as synchrotron radiation produced by the shock interaction between the jet and the gaseous circumnuclear medium (CNM). At early times after the onset of the jet (t < 5-10 days) a reverse shock propagates through and decelerates the ejecta, while at later times the outflow approaches the Blandford-McKee self-similar evolution (possibly modified by additional late energy injection). The achromatic break in the radio light curve of J1644+57 is naturally explained as the transition between these phases. We show that the temporal indices of the pre- and post-break light curve are consistent with those predicted if the CNM has a wind-type radial density profile n ~ 1/r^2. The observed synchrotron frequencies and self-absorbed flux constrain the fraction of the post-shock thermal energy in relativistic electrons epsilon_e ~ 0.03-0.1; the CNM density at 1e18 cm ~ 1-10 1/cm^3; and the initial Lorentz factor Gamma_j ~ 10-20 and opening angle theta_j ~ (0.1-1)Gamma_j^(-1) ~ 0.01-0.1 of the jet. Radio modeling thus provides robust independent evidence for a narrowly collimated outflow. Extending our model to the future evolution of J1644+57, we predict that the radio flux at low frequencies (<~ few GHz) will begin to brighten more rapidly once the characteristic frequency crosses below the self-absorption frequency on a timescale of months (indeed, such a transition may already have begun). Our results demonstrate that relativistic outflows from tidal disruption events provide a unique probe of the conditions in distant, previously inactive galactic nuclei, complementing studies of normal AGN.

preprint2012arXiv

Comparing H-alpha and HI Surveys as Means to a Complete Local Galaxy Catalog in the Advanced LIGO/Virgo Era

Identifying the electromagnetic counterparts of gravitational wave (GW) sources detected by upcoming networks of advanced ground-based interferometers will be challenging due in part to the large number of unrelated astrophysical transients within the ~10-100 square degree sky localizations. A potential way to greatly reduce the number of such false positives is to limit detailed follow-up to only those candidates near galaxies within the GW sensitivity range of ~200 Mpc for binary neutron star mergers. Such a strategy is currently hindered by the fact that galaxy catalogs are grossly incomplete within this volume. Here we compare two methods for completing the local galaxy catalog: (1) a narrow-band H-alpha imaging survey; and (2) an HI emission line radio survey. Using H-alpha fluxes, stellar masses (M_star), and star formation rates (SFR) from galaxies in the Sloan Digital Sky Survey (SDSS), combined with HI data from the GALEX Arecibo SDSS Survey and the Herschel Reference Survey, we estimate that a H-alpha survey with a luminosity sensitivity of L_H-alpha = 1e40 erg/s at 200 Mpc could achieve a completeness of f_SFR ~ 75% with respect to total SFR, but only f_Mstar ~ 33% with respect to stellar mass (due to lack of sensitivity to early-type galaxies). These numbers are significantly lower than those achieved by an idealized spectroscopic survey due to the loss of H-alpha flux resulting from resolving out nearby galaxies and the inability to correct for the underlying stellar continuum. An HI survey with sensitivity similar to the proposed WALLABY survey on ASKAP could achieve f_SFR ~ 80% and f_Mstar ~ 50%, somewhat higher than that of the H-alpha survey. Finally, both H-alpha and HI surveys should achieve > 50% completeness with respect to the host galaxies of short duration gamma-ray bursts, which may trace the population of binary neutron star mergers.

preprint2012arXiv

Global Models of Runaway Accretion in White Dwarf Debris Disks

A growing sample of white dwarfs (WDs) with metal-enriched atmospheres are accompanied by excess infrared emission, indicating that they are encircled by a compact dusty disk of solid debris. Such `WD debris disks' are thought to originate from the tidal disruption of asteroids or other minor bodies, but the precise mechanism(s) responsible for transporting matter to the WD surface remains unclear, especially in those systems with the highest inferred metal accretion rates dM_Z/dt ~ 1e8-1e10 g/s. Here we present global time-dependent calculations of the coupled evolution of the gaseous and solid components of WD debris disks. Solids transported inwards (initially due to PR drag) sublimate at tens of WD radii, producing a source of gas that accretes onto the WD surface and viscously spreads outwards in radius, where it overlaps with the solid disk. If the aerodynamic coupling between the solids and gaseous disks is sufficiently strong (and/or the gas viscosity sufficiently weak), then gas builds up near the sublimation radius faster than it can viscously spread away. Since the rate of drag-induced solid accretion increases with gas density, this results in a runaway accretion process, during which the WD accretion rate reaches values orders of magnitude higher than can be achieved by PR drag alone. We explore the evolution of WD debris disks across a wide range of physical conditions and calculate the predicted distribution of observed accretion rates dM_Z/dt, finding reasonable agreement with the current sample. Although the conditions necessary for runaway accretion are at best marginally satisfied given the minimal level of aerodynamic drag between circular gaseous and solid disks, the presence of other stronger forms of solid-gas coupling---such as would result if the gaseous disk is only mildly eccentric---substantially increase the likelihood of runaway accretion.

preprint2012arXiv

Late Time Radio Emission from X-ray Selected Tidal Disruption Events

We present new observations with the Karl G. Jansky Very Large Array of seven X-ray-selected tidal disruption events (TDEs). The radio observations were carried out between 9 and 22 years after the initial X-ray discovery, and, thus, probe the late-time formation of relativistic jets and jet interactions with the interstellar medium in these systems. We detect a compact radio source in the nucleus of the galaxy IC 3599 and a compact radio source that is a possible counterpart to RX J1420.4+5334. We find no radio counterparts for five other sources with flux density upper limits between 51 and 200 microJy (3 sigma). If the detections truly represent late radio emission associated with a TDE, then our results suggest that a fraction >~ 10% of X-ray-detected TDEs are accompanied by relativistic jets. We explore several models for producing late radio emission, including interaction of the jet with gas in the circumnuclear environment (blast wave model), and emission from the core of the jet itself. Upper limits on the radio flux density from archival observations suggest that the jet formation may have been delayed for years after the TDE, possibly triggered by the accretion rate dropping below a critical threshold of ~ 10^{-2} -- 10^{-3} Mdot_Edd. The non-detections are also consistent with this scenario; deeper radio observations can determine whether relativistic jets are present in these systems. The emission from RX J1420.4+5334 is also consistent with the predictions of the blast wave model, however the radio emission from IC 3599 is substantially underluminous, and its spectral slope is too flat, relative to the blast wave model expectations. Future radio monitoring of IC 3599 and RX J1420.4+5334 will help to better constrain the nature of the jets in these systems.

preprint2012arXiv

Nuclear Dominated Accretion Flows in Two Dimensions. I. Torus Evolution with Parametric Microphysics

We explore the evolution of radiatively inefficient accretion disks in which nuclear reactions are dynamically important (`Nuclear Dominated Accretion Flows', or NuDAFs). Examples of such disks are those generated by the merger of a white dwarf with a neutron star or black hole, or by the collapse of a rotating star. Here we present two-dimensional hydrodynamic simulations that systematically explore the effect of adding a single nuclear reaction to a viscous torus. The equation of state, anomalous shear stress, and nuclear reactions are given parametric forms. Our results point to the existence of two qualitatively different regimes of NuDAF evolution: (1) steady accretion with quiescent burning; or (2) detonation of the disk. These outcomes are controlled primarily by the ratio of the nuclear energy released to the enthalpy at the burning radius. Disks detonate if this ratio exceeds a critical value (~1), and if burning occurs in regions where neutrino cooling is unimportant. Thermonuclear runaways are seeded by the turbulent mixing of hot ash with cold fuel at the burning front. Subcritical disks do not explode, but instead power a persistent collimated outflow of unbound material composed primarily of ash, with a mass-loss rate that increases with the nuclear binding energy released. We discuss the implications of our results for supernova-like counterparts from astrophysical events in the NuDAF regime. In particular, detonations following a white dwarf - neutron star merger could account for some subluminous Type Ia supernovae, such as the class defined by SN 2002cx.

preprint2012arXiv

Nuclear-dominated accretion and subluminous supernovae from the merger of a white dwarf with a neutron star or black hole

We construct one dimensional steady-state models of accretion disks produced by the tidal disruption of a white dwarf (WD) by a neutron star (NS) or stellar mass black hole (BH). At radii r <~ 1e8.5-1e9 cm the midplane density and temperature are sufficiently high to burn the initial white dwarf material into increasingly heavier elements (e.g. Mg, Si, S, Ca, Fe, and Ni) at sequentially smaller radii. When the energy released by nuclear reactions is comparable to that released gravitationally, we term the disk a nuclear-dominated accretion flow (NuDAF). At small radii <~1e7 cm Fe photo-disintegrates into He and then free nuclei, and cooling by neutrinos may be efficient. At the high accretion rates of relevance ~ 0.1-1e-4 Msun/s, most of the disk is radiatively inefficient and prone to outflows powered by viscous dissipation and nuclear burning. Outflow properties are calculated by requiring that material in the midplane be marginally bound (Bernoulli constant <~ 0), due (in part) to cooling by matter escaping the disk. For reasonable assumptions regarding the properties of disk winds, we show that a significant fraction >50-80 per cent of the total WD mass is unbound. The ejecta composition is predominantly O, C, Si, Mg, Ne, Fe, and S [He, C, Si, S, Ar, and Fe], in the case of C-O [He] WDs, respectively, along with a small quantity ~1e-3-1e-2 Msun of radioactive Ni56 and, potentially, a trace amount of H. We use our results to evaluate possible EM counterparts of WD-NS/BH mergers, including optical transients powered by the radioactive decay of Ni56 and radio transients powered by the interaction of the ejecta with the interstellar medium. We address whether recently discovered subluminous Type I supernovae result from WD-NS/BH mergers. Our results also have implications for accretion following the core collapse of massive stars in collapsar models for gamma-ray bursts.

preprint2012arXiv

Optical and X-ray Transients from Planet-Star Mergers

We evaluate the prompt observational signatures of the merger between a massive close-in planet (a `hot Jupiter') and its host star, events with an estimated Galactic rate of ~0.1-1/yr. Depending on the ratio of the mean density of the planet rho_p to that of the star rho_star, a merger results in three possible outcomes. If rho_p/rho_star > 5, then the planet directly plunges below the stellar atmosphere before being disrupted by tidal forces. Dissipation of orbital energy creates a hot wake behind the planet, producing a EUV/soft X-ray transient as the planet sinks below the stellar surface. The peak luminosity L_X ~ 1e36 erg/s is achieved weeks-months prior to merger, after which the stellar surface is enshrouded by an outflow. The final inspiral is accompanied by an optical transient powered by the recombination of hydrogen in the outflow, which peaks at L~1e37-38 erg/s on a timescale ~days. If instead rho_planet/rho_star < 5, then Roche Lobe overflow occurs above the stellar surface. For rho_p/rho_star < 1, mass transfer is stable, resulting the planet being accreted on a relatively slow timescale. However, for 1 < rho_p/rho_star < 5, mass transfer may instead be unstable, resulting in the planet being dynamically disrupted into an accretion disk around the star. Super-Eddington outflows from the disk power an optical transient with L~1e37-38 erg/s and characteristic duration ~week-months. The disk itself becomes visible once the accretion rate become sub-Eddington, resulting in a bolometric brightening and spectral shift to the UV. Optical transients from planet merger events may resemble classical novae, but are distinguished by lower ejecta mass and velocity ~100s km/s, and by hard pre- and post-cursor emission, respectively. Promising search strategies include combined optical, UV, and X-ray surveys of nearby massive galaxies with cadences from days to months.

preprint2011arXiv

A relativistic jetted outburst from a massive black hole fed by a tidally disrupted star

While gas accretion onto some massive black holes (MBHs) at the centers of galaxies actively powers luminous emission, the vast majority of MBHs are considered dormant. Occasionally, a star passing too near a MBH is torn apart by gravitational forces, leading to a bright panchromatic tidal disruption flare (TDF). While the high-energy transient Swift J164449.3+573451 ("Sw 1644+57") initially displayed none of the theoretically anticipated (nor previously observed) TDF characteristics, we show that the observations (Levan et al. 2011) suggest a sudden accretion event onto a central MBH of mass ~10^6-10^7 solar masses. We find evidence for a mildly relativistic outflow, jet collimation, and a spectrum characterized by synchrotron and inverse Compton processes; this leads to a natural analogy of Sw 1644+57 with a smaller-scale blazar. The phenomenologically novel Sw 1644+57 thus connects the study of TDFs and active galaxies, opening a new vista on disk-jet interactions in BHs and magnetic field generation and transport in accretion systems.

preprint2011arXiv

Heavy Nuclei Synthesized in Gamma-Ray Burst Outflows as the Source of UHECRs

Recent measurements by the Pierre Auger Observatory suggest that the composition of ultra-high energy cosmic rays (UHECRs) becomes dominated by heavy nuclei at high energies. However, until now there has been no astrophysical motivation for considering a source highly enriched in heavy elements. Here we demonstrate that the outflows from Gamma-Ray Bursts (GRBs) may indeed be composed primarily of nuclei with masses A ~ 40-200, which are synthesized as hot material expands away from the central engine. In particular, if the jet is magnetically-dominated (rather than a thermally-driven fireball) its low entropy enables heavy elements to form efficiently. Adopting the millisecond proto-magnetar model for the GRB central engine, we show that heavy nuclei are both synthesized in proto-magnetar winds and can in principle be accelerated to energies >1e20 eV in the shocks or regions of magnetic reconnection that are responsible for powering the GRB. Similar results may apply to accretion-powered GRB models if the jet originates from a magnetized disk wind. Depending on the precise distribution of nuclei synthesized, we predict that the average primary mass may continue to increase beyond Fe group elements at the highest energies, possibly reaching the A ~ 90 (Zirconium), A ~ 130 (Tellurium), or even A ~ 195 (Platinum) peaks. Future measurements of the UHECR composition at energies >~ 1e20 eV can thus confirm or constrain our model and, potentially, probe the nature of GRB outflows. The longer attenuation length of ultra-heavy nuclei through the extragalactic background light greatly expands the volume of accesible sources and alleviates the energetic constraints on GRBs as the source of UHECRs.

preprint2011arXiv

Radio transients from stellar tidal disruption by massive black holes

The tidal disruption of a star by a supermassive black hole provides us with a rare glimpse of these otherwise dormant beasts. It has long been predicted that the disruption will be accompanied by a thermal `flare', powered by the accretion of bound stellar debris. Several candidate disruptions have been discovered in this manner at optical, UV and X-ray wavelengths. Here we explore the observational consequences if a modest fraction of the accretion power is channeled into an ultra-relativistic outflow. We show that a relativistic jet decelerates due to its interaction with the interstellar medium at sub-parsec distances from the black hole. Synchrotron radiation from electrons accelerated by the reverse shock powers a bright radio-infrared transient that peaks on a timescale ~1 yr after disruption. Emission from the forward shock may be detectable for several years after the peak. Deep radio follow-up observations of tidal disruption candidates at late times can test for the presence of relativistic ejecta. Upcoming radio transient surveys may independently discover tens to hundreds of tidal disruptions per year, complimenting searches at other wavelengths. Non-thermal emission from tidal disruption probes the physics of jet formation under relatively clean conditions, in which the flow parameters are independently constrained.

preprint2011arXiv

Relic proto-stellar disks and the origin of luminous circumstellar interaction in core collapse supernovae

A small fraction of core collapse supernovae (SNe) show evidence that the outgoing blast wave has encountered a substantial mass ~ 1-10 M_sun of circumstellar matter (CSM) at radii ~100-1000 AU, much more than can nominally be explained by pre-explosion stellar winds. In extreme cases this interaction may power the most luminous, optically-energetic SNe yet discovered. Interpretations for the origin of the CSM have thus far centered on explosive eruptions from the star just ~ years to decades prior to the core collapse. Here we consider an alternative possibility that the inferred CSM is a relic disk left over from stellar birth. We investigate this hypothesis by calculating the evolution of proto-stellar disks around massive stars following their early embedded phase using a self-similar accretion model. We identify an initial gravitationally-unstable ("gravito-turbulent") phase, followed by a much longer period of irradiation-supported accretion during which less effective non-gravitational forms of angular momentum transport dominate. Although external influences, such as the presence of a wide binary companion, may preclude disk survival in many systems, we find that massive (~1-10 M_sun) disks can preferentially survive around the most massive stars. Reasons for this perhaps counterintuitive include (1) the shorter stellar lifetimes and (2) large photo-evaporation radii (~ 1000 AU) of very massive stars; (3) suppression of the magneto-rotational instability due to the shielding from external sources of ionization; and (4) relative invulnerability of massive disks to lower mass stellar collisions and luminous blue variable eruptions. Because very luminous SNe are rare, testing the relic disk model requires constraining the presence of long-lived disks around a small fraction of very massive stars.

preprint2011arXiv

What is the Most Promising Electromagnetic Counterpart of a Neutron Star Binary Merger?

The final inspiral of double neutron star and neutron star-black hole binaries are likely to be detected by advanced networks of ground-based gravitational wave (GW) interferometers. Maximizing the science returns from such a discovery will require the identification and localization of an electromagnetic (EM) counterpart. Here we critically evaluate and compare several possible counterparts, including short-duration gamma-ray bursts (SGRBs), "orphan" optical and radio afterglows, and ~day-long optical transients powered by the radioactive decay of heavy nuclei synthesized in the merger ejecta ("kilonovae"). We assess the promise of each counterpart in terms of four "Cardinal Virtues": detectability, high fraction, identifiability, and positional accuracy. Taking into account the search strategy for typical error regions of ~10s degs sq., we conclude that SGRBs are the most useful to confirm the cosmic origin of a few GW events, and to test the association with NS mergers. However, for the more ambitious goal of localizing and obtaining redshifts for a large sample of GW events, kilonovae are instead preferred. Off-axis optical afterglows will be detectable for at most ~10% of all events, while radio afterglows are promising only for the unique combination of energetic relativistic ejecta in a high density medium, and even then will require hundreds of hours of EVLA time per event. Our main recommendations are:(i) an all-sky gamma-ray satellite is essential for temporal coincidence detections, and for GW searches of gamma-ray triggered events; (ii) LSST should adopt a 1-day cadence follow-up strategy, ideally with ~0.5 hr per pointing to cover GW error regions (the standard 4-day cadence and depth will severely limit the probability of a unique identification); and (iii) radio searches should only focus on the relativistic case, which requires observations for a few months.

preprint2010arXiv

Proto-Magnetars as GRB Central Engines: Uncertainties, Limitations, and Particulars

The millisecond proto-magnetar model for the central engine of long-duration gamma-ray bursts is briefly reviewed. Limitations and uncertainties in the model are highlighted. A short discussion of the maximum energy, maximum duration, radiative efficiency, jet formation mechanism, late-time energy injection, and (non-)association with supernovae of millisecond magnetar-powered GRBs is provided.

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 "proto-magnetars" 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.