Source author record

B. D. Metzger

B. 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

32works
8topics
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

32 published item(s)

preprint2022arXiv

Evidence for X-ray Emission in Excess to the Jet Afterglow Decay 3.5 yrs After the Binary Neutron Star Merger GW 170817: A New Emission Component

For the first $\sim3$ years after the binary neutron star merger event GW 170817 the radio and X-ray radiation has been dominated by emission from a structured relativistic off-axis jet propagating into a low-density medium with n $< 0.01\,\rm{cm^{-3}}$. We report on observational evidence for an excess of X-ray emission at $δt>900$ days after the merger. With $L_x\approx5\times 10^{38}\,\rm{erg\,s^{-1}}$ at 1234 days, the recently detected X-ray emission represents a $\ge 3.2\,σ$ (Gaussian equivalent) deviation from the universal post jet-break model that best fits the multi-wavelength afterglow at earlier times. In the context of JetFit afterglow models, current data represent a departure with statistical significance $\ge 3.1\,σ$, depending on the fireball collimation, with the most realistic models showing excesses at the level of $\ge 3.7\,σ$. A lack of detectable 3 GHz radio emission suggests a harder broad-band spectrum than the jet afterglow. These properties are consistent with the emergence of a new emission component such as synchrotron radiation from a mildly relativistic shock generated by the expanding merger ejecta, i.e. a kilonova afterglow. In this context, we present a set of ab-initio numerical-relativity BNS merger simulations that show that an X-ray excess supports the presence of a high-velocity tail in the merger ejecta, and argues against the prompt collapse of the merger remnant into a black hole. Radiation from accretion processes on the compact-object remnant represents a viable alternative. Neither a kilonova afterglow nor accretion-powered emission have been observed before, as detections of BNS mergers at this phase of evolution are unprecedented.

preprint2022arXiv

Extragalactic Millimeter Transients in the Era of Next Generation CMB Surveys

The next generation of wide-field cosmic microwave background (CMB) surveys are uniquely poised to open a new window for time-domain astronomy in the millimeter band. Here we explore the discovery phase space for extragalactic transients with near-term and future CMB experiments to characterize the expected population. We use existing millimeter-band light curves of known transients (gamma-ray bursts, tidal disruption events, fast blue optical transients, neutron star mergers) and theoretical models, in conjunction with known and estimated volumetric rates. Using Monte Carlo simulations of various CMB survey designs (area, cadence, depth, duration) we estimate the detection rates and the resulting light curve characteristics. We find that existing and near-term surveys will find tens to hundreds of long-duration gamma-ray bursts (LGRBs), driven primarily by detections of the reverse shock emission, and including off-axis LGRBs. Next-generation experiments (CMB-S4, CMB-HD) will find tens of fast blue optical transients (FBOTs) in the nearby universe and will detect a few tidal disruption events. CMB-HD will additionally detect a small number of short gamma-ray bursts, where these will be discovered within the detection volume of next generation gravitational wave experiments like Cosmic Explorer.

preprint2022arXiv

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

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

preprint2021arXiv

Late-Time Evolution and Modeling of the Off-Axis Gamma-ray Burst Candidate FIRST J141918.9+394036

We present new radio and optical data, including very long baseline interferometry, as well as archival data analysis, for the luminous decades-long radio transient FIRST J141918.9+394036. The radio data reveal a synchrotron self-absorption peak around 0.3 GHz and a radius of around 1.3 mas (0.5 pc) 26 years post-discovery, indicating a blastwave energy $\sim5 \times 10^{50}$ erg. The optical spectrum shows a broad [OIII]$λ$4959,5007 emission-line that may indicate collisional-excitation in the host galaxy, but its association with the transient cannot be ruled out. The properties of the host galaxy are suggestive of a massive stellar progenitor that formed at low metallicity. Based on the radio light curve, blastwave velocity, energetics, nature of the host galaxy and transient rates we find that the properties of FIRST J1419+39 are most consistent with long gamma-ray burst (LGRB) afterglows. Other classes of (optically-discovered) stellar explosions as well as neutron star mergers are disfavored, and invoking any exotic scenario may not be necessary. It is therefore likely that FIRST J1419+39 is an off-axis LGRB afterglow (as suggested by Law et al. and Marcote et al.), and under this premise the inverse beaming fraction is found to be $f_b^{-1}\simeq280^{+700}_{-200}$, corresponding to an average jet half-opening angle $<θ_j>\simeq5^{+4}_{-2}$ degrees (68% confidence), consistent with previous estimates. From the volumetric rate we predict that surveys with the VLA, ASKAP and MeerKAT will find a handful of FIRST J1419+39-like events over the coming years.

preprint2020arXiv

The Broad-band Counterpart of the Short GRB 200522A at $z=0.5536$: A Luminous Kilonova or a Collimated Outflow with a Reverse Shock?

We present the discovery of the radio afterglow and near-infrared (NIR) counterpart of the Swift short GRB 200522A, located at a small projected offset of $\approx 1$ kpc from the center of a young, star-forming host galaxy at $z=0.5536$. The radio and X-ray luminosities of the afterglow are consistent with those of on-axis cosmological short GRBs. The NIR counterpart, revealed by our HST observations at a rest-frame time of $\approx2.3$ days, has a luminosity of $\approx (1.3-1.7) \times 10^{42}$ erg s$^{-1}$. This is substantially lower than on-axis short GRB afterglow detections, but is a factor of $\approx 8$-$17$ more luminous than the kilonova of GW170817, and significantly more luminous than any kilonova candidate for which comparable observations exist. The combination of the counterpart's color ($i-y = -0.08\pm 0.21$; rest-frame) and luminosity cannot be explained by standard radioactive heating alone. We present two scenarios to interpret the broad-band behavior of GRB 200522A: a synchrotron forward shock with a luminous kilonova (potentially boosted by magnetar energy deposition), or forward and reverse shocks from a $\approx14^{\circ}$, relativistic ($Γ_0 \gtrsim 80$) jet. Models which include a combination of enhanced radioactive heating rates, low-lanthanide mass fractions, or additional sources of heating from late-time central engine activity may provide viable alternate explanations. If a stable magnetar was indeed produced in GRB 200522A, we predict that late-time radio emission will be detectable starting $\approx 0.3$-$6$ years after the burst for a deposited energy of $\approx 10^{53}$ erg. Counterparts of similar luminosity to GRB 200522A associated with gravitational wave events will be detectable with current optical searches to $\approx\!250$ Mpc.

preprint2020arXiv

Wandering Massive Black Holes or Analogs of the First Repeating Fast Radio Burst?

The discovery of a persistent radio source coincident with the first repeating fast radio burst, FRB 121102, and offset from the center of its dwarf host galaxy has been used as evidence for a link with young millisecond magnetars born in superluminous supernovae (SLSNe) or long-duration gamma-ray bursts (LGRBs). A prediction of this scenario is that compact radio sources offset from the centers of dwarf galaxies may serve as signposts for at least some FRBs. Recently, Reines et al. 2019 presented the discovery of 20 such radio sources in nearby ($z\lesssim 0.055$) dwarf galaxies, and argued that these cannot be explained by emission from HII regions, normal supernova remnants, or normal radio supernovae. Instead, they attribute the emission to accreting wandering massive black holes. Here, we explore the alternative possibility that these sources are analogs of FRB 121102. We compare their properties -- radio luminosities, spectral energy distributions, light curves, ratios of radio-to-optical flux, and spatial offsets -- to FRB 121102, a few other well-localized FRBs, and potentially related systems, and find that these are all consistent as arising from the same population. We further compare their properties to the magnetar nebula model used to explain FRB 121102, as well as to theoretical off-axis LGRB light curves, and find overall consistency. Finally, we find a consistent occurrence rate relative to repeating FRBs and LGRBs. We outline key follow-up observations to further test these possible connections.

preprint2019arXiv

Two years of non-thermal emission from the binary neutron star merger GW170817: rapid fading of the jet afterglow and first constraints on the kilonova fastest ejecta

We present Chandra and VLA observations of GW170817 at ~521-743 days post merger, and a homogeneous analysis of the entire Chandra data set. We find that the late-time non-thermal emission follows the expected evolution from an off-axis relativistic jet, with a steep temporal decay $F_ν\propto t^{-1.95\pm0.15}$ and a simple power-law spectrum $F_ν\propto ν^{-0.575\pm0.007}$. We present a new method to constrain the merger environment density based on diffuse X-ray emission from hot plasma in the host galaxy and we find $n\le 9.6 \times 10^{-3}\,\rm{cm^{-3}}$. This measurement is independent from inferences based on the jet afterglow modeling and allows us to partially solve for model degeneracies. The updated best-fitting model parameters with this density constraint are a fireball kinetic energy $E_0 = 1.5_{-1.1}^{+3.6}\times 10^{49}\,\rm{erg}$ ($E_{iso}= 2.1_{-1.5}^{+6.4}\times10^{52}\, \rm{erg}$), jet opening angle $θ_{0}= 5.9^{+1.0}_{-0.7}\,\rm{deg}$ with characteristic Lorentz factor $Γ_j = 163_{-43}^{+23}$, expanding in a low-density medium with $n_0 = 2.5_{-1.9}^{+4.1} \times 10^{-3}\, \rm{cm^{-3}}$ and viewed $θ_{obs} = 30.4^{+4.0}_{-3.4}\, \rm{deg}$ off-axis. The synchrotron emission originates from a power-law distribution of electrons with $p=2.15^{+0.01}_{-0.02}$. The shock microphysics parameters are constrained to $ε_{\rm{e}} = 0.18_{-0.13}^{+0.30}$ and $ε_{\rm{B}}=2.3_{-2.2}^{+16.0} \times 10^{-3}$. We investigate the presence of X-ray flares and find no statistically significant evidence of $\ge2.5σ$ of temporal variability at any time. Finally, we use our observations to constrain the properties of synchrotron emission from the deceleration of the fastest kilonova ejecta with energy $E_k^{KN}\propto (Γβ)^{-α}$ into the environment, finding that shallow stratification indexes $α\le6$ are disfavored.

preprint2016arXiv

A Dark Energy Camera Search for an Optical Counterpart to the First Advanced LIGO Gravitational Wave Event GW150914

We report initial results of a deep search for an optical counterpart to the gravitational wave event GW150914, the first trigger from the Advanced LIGO gravitational wave detectors. We used the Dark Energy Camera (DECam) to image a 102 deg$^2$ area, corresponding to 38% of the initial trigger high-probability sky region and to 11% of the revised high-probability region. We observed in i and z bands at 4-5, 7, and 24 days after the trigger. The median $5σ$ point-source limiting magnitudes of our search images are i=22.5 and z=21.8 mag. We processed the images through a difference-imaging pipeline using templates from pre-existing Dark Energy Survey data and publicly available DECam data. Due to missing template observations and other losses, our effective search area subtends 40 deg$^{2}$, corresponding to 12% total probability in the initial map and 3% of the final map. In this area, we search for objects that decline significantly between days 4-5 and day 7, and are undetectable by day 24, finding none to typical magnitude limits of i= 21.5,21.1,20.1 for object colors (i-z)=1,0,-1, respectively. Our search demonstrates the feasibility of a dedicated search program with DECam and bodes well for future research in this emerging field.

preprint2016arXiv

A Dark Energy Camera Search for Missing Supergiants in the LMC After the Advanced LIGO Gravitational Wave Event GW150914

The collapse of the core of a star is expected to produce gravitational radiation. While this process will usually produce a luminous supernova, the optical signatue could be subluminous and a direct collapse to a black hole, with the star just disappearing, is possible. The gravitational wave event GW150914 reported by the LIGO Virgo Collaboration (LVC) on 2015 September 16, was detected by a burst analysis and whose high probability spatial localization included the Large Magellanic Cloud. Shortly after the announcement of the event, we used the Dark Energy Camera to observe 102 deg$^2$ of the localization area, including a 38 deg$^2$ area centered on the LMC. Using a catalog of 152 LMC luminous red supergiants, candidates to undergo a core collapse without a visible supernova, we find that the positions of 144 of these are inside our images, and that all are detected - none have disappeared. There are other classes of candidates: we searched existing catalogs of red supergiants, yellow supergiants, Wolf-Rayet stars, and luminous blue variable stars, recovering all that were inside the imaging area. Based on our observations, we conclude that it is unlikely that GW150914 was caused by the core collapse of a supergiant in the LMC, consistent with the LIGO Collaboration analyses of the gravitational wave form as best described by a binary black hole merger. We discuss how to generalize this search for future very nearby core collapse candidates.

preprint2016arXiv

A DECam Search for an Optical Counterpart to the LIGO Gravitational Wave Event GW151226

We report the results of a Dark Energy Camera (DECam) optical follow-up of the gravitational wave (GW) event GW151226, discovered by the Advanced LIGO detectors. Our observations cover 28.8 deg$^2$ of the localization region in the $i$ and $z$ bands (containing 3% of the BAYESTAR localization probability), starting 10 hours after the event was announced and spanning four epochs at $2-24$ days after the GW detection. We achieve $5σ$ point-source limiting magnitudes of $i\approx21.7$ and $z\approx21.5$, with a scatter of $0.4$ mag, in our difference images. Given the two day delay, we search this area for a rapidly declining optical counterpart with $\gtrsim 3σ$ significance steady decline between the first and final observations. We recover four sources that pass our selection criteria, of which three are cataloged AGN. The fourth source is offset by $5.8$ arcsec from the center of a galaxy at a distance of 187 Mpc, exhibits a rapid decline by $0.5$ mag over $4$ days, and has a red color of $i-z\approx 0.3$ mag. These properties roughly match the expectations for a kilonova. However, this source was detected several times, starting $94$ days prior to GW151226, in the Pan-STARRS Survey for Transients (dubbed as PS15cdi) and is therefore unrelated to the GW event. Given its long-term behavior, PS15cdi is likely a Type IIP supernova that transitioned out of its plateau phase during our observations, mimicking a kilonova-like behavior. We comment on the implications of this detection for contamination in future optical follow-up observations.

preprint2016arXiv

Localization and broadband follow-up of the gravitational-wave transient GW150914

A gravitational-wave (GW) transient was identified in data recorded by the Advanced Laser Interferometer Gravitational-wave Observatory (LIGO) detectors on 2015 September 14. The event, initially designated G184098 and later given the name GW150914, is described in detail elsewhere. By prior arrangement, preliminary estimates of the time, significance, and sky location of the event were shared with 63 teams of observers covering radio, optical, near-infrared, X-ray, and gamma-ray wavelengths with ground- and space-based facilities. In this Letter we describe the low-latency analysis of the GW data and present the sky localization of the first observed compact binary merger. We summarize the follow-up observations reported by 25 teams via private Gamma-ray Coordinates Network circulars, giving an overview of the participating facilities, the GW sky localization coverage, the timeline and depth of the observations. As this event turned out to be a binary black hole merger, there is little expectation of a detectable electromagnetic (EM) signature. Nevertheless, this first broadband campaign to search for a counterpart of an Advanced LIGO source represents a milestone and highlights the broad capabilities of the transient astronomy community and the observing strategies that have been developed to pursue neutron star binary merger events. Detailed investigations of the EM data and results of the EM follow-up campaign are being disseminated in papers by the individual teams.

preprint2016arXiv

Numerical simulations of the jetted tidal disruption event Swift J1644+57

In this work we focus on the technical details of the numerical simulations of the non- thermal transient Swift J1644+57, whose emission is probably produced by a two-component jet powered by a tidal disruption event. In this context we provide details of the coupling between the relativistic hydrodynamic simulations and the radiative transfer code. First, we consider the technical demands of one-dimensional simulations of a fast relativistic jet, and show to what extent (for the same physical parameters of the model) do the computed light curves depend on the numerical parameters of the different codes employed. In the second part we explain the difficulties of computing light curves from axisymmetric two dimensonal simulations and discuss a procedure that yields an acceptable tradeoff between the computational cost and the quality of the results.

preprint2016arXiv

SN 2015bn: a detailed multi-wavelength view of a nearby superluminous supernova

We present observations of SN 2015bn (= PS15ae = CSS141223-113342+004332 = MLS150211-113342+004333), a Type I superluminous supernova (SLSN) at redshift $z=0.1136$. As well as being one of the closest SLSNe I yet discovered, it is intrinsically brighter ($M_U\approx-23.1$) and in a fainter galaxy ($M_B\approx-16.0$) than other SLSNe at $z\sim0.1$. We used this opportunity to collect the most extensive dataset for any SLSN I to date, including densely-sampled spectroscopy and photometry, from the UV to the NIR, spanning $-$50 to +250 days from optical maximum. SN 2015bn fades slowly, but exhibits surprising undulations in the light curve on a timescale of 30-50 days, especially in the UV. The spectrum shows extraordinarily slow evolution except for a rapid transformation between +7 and +20-30 days. No narrow emission lines from slow-moving material are observed at any phase. We derive physical properties including the bolometric luminosity, and find slow velocity evolution and non-monotonic temperature and radial evolution. A deep radio limit rules out a healthy off-axis gamma-ray burst, and places constraints on the pre-explosion mass loss. The data can be consistently explained by a $\gtrsim10\,{\rm M}_\odot$ stripped progenitor exploding with $\sim 10^{51}\,$erg kinetic energy, forming a magnetar with a spin-down timescale of $\sim20$ days (thus avoiding a gamma-ray burst) that reheats the ejecta and drives ionization fronts. The most likely alternative scenario -- interaction with $\sim20\,{\rm M}_\odot$ of dense, inhomogeneous circumstellar material -- can be tested with continuing radio follow-up.

preprint2016arXiv

Supplement: Localization and broadband follow-up of the gravitational-wave transient GW150914

This Supplement provides supporting material for arXiv:1602.08492 . We briefly summarize past electromagnetic (EM) follow-up efforts as well as the organization and policy of the current EM follow-up program. We compare the four probability sky maps produced for the gravitational-wave transient GW150914, and provide additional details of the EM follow-up observations that were performed in the different bands.

preprint2016arXiv

The influence of circumnuclear environment on the radio emission from TDE jets

Dozens of stellar tidal disruption events (TDEs) have been identified at optical, UV and X-ray wavelengths. A small fraction of these, most notably Swift J1644+57, produce radio synchrotron emission, consistent with a powerful, relativistic jet shocking the surrounding circumnuclear gas. The dearth of similar non-thermal radio emission in the majority of TDEs may imply that powerful jet formation is intrinsically rare, or that the conditions in galactic nuclei are typically unfavorable for producing a detectable signal. Here we explore the latter possibility by constraining the radial profile of the gas density encountered by a TDE jet using a one-dimensional model for the circumnuclear medium which includes mass and energy input from a stellar population. Near the jet Sedov radius of 10$^{18}$ cm, we find gas densities in the range of $n_{18} \sim$ 0.1$-$1000 cm$^{-3}$ across a wide range of plausible star formation histories. Using one- and two-dimensional relativistic hydrodynamical simulations, we calculate the synchrotron radio light curves of TDE jets (as viewed both on and off-axis) across the allowed range of density profiles. We find that bright radio emission would be produced across the plausible range of nuclear gas densities by jets as powerful as Swift J1644+57, and we quantify the relationship between the radio luminosity and jet energy. We use existing radio detections and upper limits to constrain the energy distribution of TDE jets. Radio follow up observations several months to several years after the TDE candidate will strongly constrain the energetics of any relativistic flow.

preprint2015arXiv

Late time multi wavelength observations of Swift J1644+5734: A luminous optical/IR bump and quiescent X-ray emission

We present late-time multi-wavelength observations of Swift J1644+57, suggested to be a relativistic tidal disruption flare (TDF). Our observations extend to >4 years from discovery, and show that 1.4 years after outburst the relativistic jet switched-off on a timescale less than tens of days, corresponding to a power-law decay faster than $t^{-70}$. Beyond this point weak X-rays continue to be detected at an approximately constant luminosity of $L_X \sim 5 \times 10^{42}$ erg s$^{-1}$, and are marginally inconsistent with a continuing decay of $t^{-5/3}$, similar to that seen prior to the switch-off. Host photometry enables us to infer a black hole mass of $M_{BH}=3 \times 10^6$ M$_{\odot}$, consistent with the late time X-ray luminosity arising from sub-Eddington accretion onto the black hole in the form of either an unusually optically faint AGN or a slowly varying phase of the transient. Optical/IR observations show a clear bump in the light curve at timescales of 30-50 days, with a peak magnitude (corrected for host galaxy extinction) of $M_R \sim -22-23$. The luminosity of the bump is significantly higher than seen in other, non-relativistic TDFs and does not match any re-brightening seen at X-ray or radio wavelengths. Its luminosity, light curve shape and spectrum are broadly similar to those seen in superluminous SNe, although subject to large uncertainties in the correction of the significant host extinction. We discuss these observations in the context of both TDF and massive star origins for Swift J1644+5734 and other candidate relativistic tidal flares.

preprint2015arXiv

The radio afterglow of Swift J1644+57 reveals a powerful jet with fast core and slow sheath

We model the non-thermal transient Swift J1644+57 as resulting from a relativistic jet powered by the accretion of a tidally-disrupted star onto a super-massive black hole. Accompanying synchrotron radio emission is produced by the shock interaction between the jet and the dense circumnuclear medium, similar to a gamma-ray burst afterglow. An open mystery, however, is the origin of the late-time radio rebrightening, which occurred well after the peak of the jetted X-ray emission. Here, we systematically explore several proposed explanations for this behavior by means of multi-dimensional hydrodynamic simulations coupled to a self-consistent radiative transfer calculation of the synchrotron emission. Our main conclusion is that the radio afterglow of Swift J1644+57 is not naturally explained by a jet with a one-dimensional top-hat angular structure. However, a more complex angular structure comprised of an ultra-relativistic core (Lorentz factor $Γ\sim 10$) surrounded by a slower ($Γ\sim $ 2) sheath provides a reasonable fit to the data. Such a geometry could result from the radial structure of the super-Eddington accretion flow or as the result of jet precession. The total kinetic energy of the ejecta that we infer of $\sim$ few $10^{53}\,$erg requires a highly efficient jet launching mechanism. Our jet model providing the best fit to the light curve of the on-axis event Swift J1644+57 is used to predict the radio light curves for off-axis viewing angles. Implications for the presence of relativistic jets from TDEs detected via their thermal disk emission, as well as the prospects for detecting orphan TDE afterglows with upcoming wide-field radio surveys and resolving the jet structure with long baseline interferometry, are discussed.

preprint2013arXiv

Signatures of magnetar central engines in short GRB lightcurves

A significant fraction of the Long Gamma-ray Bursts (LGRBs) in the Swift sample have a plateau phase showing evidence of ongoing energy injection. We suggest that many Short Gamma-ray Bursts (SGRBs) detected by the Swift satellite also show evidence of energy injection. Explaining this observation within the typical SGRB progenitor model is challenging as late time accretion, often used to explain plateaus in LGRBs, is likely to be absent from the SGRB population. Alternatively, it is predicted that the remnant of NS-NS mergers may not collapse immediately to a BH (or even collapse at all), forming instead an unstable millisecond pulsar (magnetar) which powers a plateau phase in the X-ray lightcurve. By fitting the magnetar model to all of the Swift SGRBs observed until May 2012, we find that about half can be clearly fitted with a magnetar plateau phase while the rest are consistent with forming a magnetar but the data are insufficient to prove a plateau phase. More data, both at early times and a larger sample, are required to confirm this. This model can be tested by detecting the gravitational wave emission from events using the next generation gravitational wave observatories.

preprint2013arXiv

Superluminous X-rays from a superluminous supernova

The discovery of a population of superluminous supernovae (SLSNe), with peak luminosities a factor of ~100 brighter than normal SNe (typically SLSNe have M_V <-21), has shown an unexpected diversity in core-collapse supernova properties. Numerous models have been postulated for the nature of these events, including a strong interaction of the shockwave with a dense circumstellar environment, a re-energizing of the outflow via a central engine, or an origin in the catastrophic destruction of the star following a loss of pressure due to pair production in an extremely massive stellar core (so-called pair instability supernovae). Here we consider constraints that can be placed on the explosion mechanism of Hydrogen-poor SLSNe (SLSNe-I) via X-ray observations, with XMM-Newton, Chandra and Swift, and show that at least one SLSNe-I is likely the brightest X-ray supernovae ever observed, with L_X ~ 10^45 ergs/s, ~150 days after its initial discovery. This is a luminosity 3 orders of magnitude higher than seen in other X-ray supernovae powered via circumstellar interactions. Such high X-ray luminosities are sufficient to ionize the ejecta and markedly reduce the optical depth, making it possible to see deep into the ejecta and any source of emission that resides there. Alternatively, an engine could have powered a moderately relativistic jet external to the ejecta, similar to those seen in gamma-ray bursts. If the detection of X-rays does require an engine it implies that these SNe do create compact objects, and that the stars are not completely destroyed in a pair instability event. Future observations will determine which, if any, of these mechanisms are at play in superluminous supernovae.

preprint2011arXiv

Nickel-Rich Outflows Produced by the Accretion-Induced Collapse of White Dwarfs: Lightcurves and Spectra

The accretion-induced collapse (AIC) of a white dwarf to form a neutron star can leave behind a rotationally supported disk with mass of up to ~ 0.1 M_sun. The disk is initially composed of free nucleons but as it accretes and spreads to larger radii, the free nucleons recombine to form helium, releasing sufficient energy to unbind the remaining disk. Most of the ejected mass fuses to form Ni56 and other iron group elements. We present spherically symmetric radiative transfer calculations of the transient powered by the radioactive heating of this ejecta. For an ejecta mass of 1e-2 M_sun (3e-3 M_sun), the lightcurve peaks after <~ 1 day with a peak bolometric luminosity ~ 2e41 erg/s (~ 5e40 erg/s), i.e., a "kilonova"; the decay time is ~ 4 (2) days. Overall, the spectra redden with time reaching U-V ~ 4 after ~ 1 day; the optical colors (B-V) are, however, somewhat blue. Near the peak in the lightcurve, the spectra are dominated by Doppler broadened Nickel features, with no distinct spectral lines present. At ~ 3-5 days, strong Calcium lines are present in the infrared, although the Calcium mass fraction is only ~ 1e-4.5. If rotationally supported disks are a common byproduct of AIC, current and upcoming transient surveys such as the Palomar Transient Factory should detect a few AIC per year for an AIC rate of ~ 1e-2 of the Type Ia rate. We discuss ways of distinguishing AIC from other rapid, faint transients, including .Ia's and the ejecta from binary neutron star mergers.

preprint2011arXiv

Polarimetry of the transient relativistic jet of GRB 110328 / Swift J164449.3+573451

We present deep infrared (Ks band) imaging polarimetry and radio (1.4 and 4.8 GHz) polarimetry of the enigmatic transient Swift J164449.3+573451. This source appears to be a short lived jet phenomenon in a galaxy at redshift z = 0.354, activated by a sudden mass accretion onto the central massive black hole, possibly caused by the tidal disruption of a star. We aim to find evidence for this scenario through linear polarimetry, as linear polarisation is a sensitive probe of jet physics, source geometry and the various mechanisms giving rise to the observed radiation. We find a formal Ks band polarisation measurement of P_lin = 7.4 +/- 3.5 % (including systematic errors). Our radio observations show continuing brightening of the source, which allows sensitive searches for linear polarisation as a function of time. We find no evidence of linear polarisation at radio wavelengths of 1.4 GHz and 4.8 GHz at any epoch, with the most sensitive 3 sigma limits as deep as 2.1%. These upper limits are in agreement with expectations from scenarios in which the radio emission is produced by the interaction of a relativistic jet with a dense circumsource medium. We further demonstrate how the polarisation properties can be used to derive properties of the jet in Swift J164449.3+573451, exploiting the similarities between this source and the afterglows of gamma-ray bursts.

preprint2011arXiv

Short GRBs with Extended Emission from Magnetar Birth: Jet Formation and Collimation

Approximately 1/4-1/2 of short duration Gamma-Ray Bursts are followed by variable X-ray emission lasting ~ 100 s with a fluence comparable or exceeding that of the initial burst itself. The long duration and significant energy of this `extended emission'(EE) poses a major challenge to the standard binary neutron star merger model. metzger08 recently proposed that the EE is powered by the spin-down of a strongly magnetized neutron star. However, the effects of surrounding material on the magnetar outflow have not yet been considered. Here we present time-dependent axisymmetric relativistic magnetohydrodynamic simulations of the interaction of the relativistic proto-magnetar wind with a surrounding 10^-1 10^-3 M_\odot envelope, which represents material ejected during the merger; the supernova following AIC; or via outflows from the accretion disk. The collision between the relativistic magnetar wind and the expanding ejecta produces a magnetized nebula inside the ejecta. A strong toroidal magnetic field builds up in the nebula, which drives a bipolar jet out through the ejecta, similar to the magnetar model developed in the case of long duration GRBs. We quantify the `break-out' time and opening angle of the jet theta_j as a function of the wind energy flux dot{E} and ejecta mass M_ej. We show that dot{E} and theta_j are inversely correlated, such that the beaming-corrected (isotropic) luminosity of the jet is primarily a function of M_ej. Both variability arguments, and the lower limit on the power of magnetar outflows capable of producing bright emission, suggest that the true opening angle of the magnetar jet must be relatively large. The model thus predicts a class of events for which the EE is observable with no associated short GRB. These may appear as long-duration GRBs or X-Ray Flashes, which may be detected by future all-sky X-ray survey missions.

preprint2011arXiv

The formation of high-field magnetic white dwarfs from common envelopes

The origin of highly-magnetized white dwarfs has remained a mystery since their initial discovery. Recent observations indicate that the formation of high-field magnetic white dwarfs is intimately related to strong binary interactions during post-main-sequence phases of stellar evolution. If a low-mass companion, such as a planet, brown dwarf, or low-mass star is engulfed by a post-main-sequence giant, the hydrodynamic drag in the envelope of the giant leads to a reduction of the companion's orbit. Sufficiently low-mass companions in-spiral until they are shredded by the strong gravitational tides near the white dwarf core. Subsequent formation of a super-Eddington accretion disk from the disrupted companion inside a common envelope can dramatically amplify magnetic fields via a dynamo. Here, we show that these disk-generated fields are sufficiently strong to explain the observed range of magnetic field strengths for isolated, high-field magnetic white dwarfs. A higher-mass binary analogue may also contribute to the origin of magnetar fields.

preprint2011arXiv

The Proto-Magnetar Model for Gamma-Ray Bursts

Long duration Gamma-Ray Bursts (GRBs) originate from the core collapse of massive stars, but the identity of the central engine remains elusive. Previous work has shown that rapidly spinning, strongly magnetized proto-neutron stars (`millisecond proto-magnetars') produce outflows with energies, timescales, and magnetizations sigma_0 (maximum Lorentz factor) that are consistent with those required to produce long GRBs. Here we extend this work in order to construct a self-consistent model that directly connects the properties of the central engine to the observed prompt emission. Just after the launch of the supernova shock, a wind heated by neutrinos is driven from the proto-magnetar. The outflow is collimated into a bipolar jet by its interaction with the star. As the magnetar cools, the wind becomes ultra-relativistic and Poynting-flux dominated (sigma_0 >> 1) on a timescale comparable to that required for the jet to clear a cavity through the star. Although the site and mechanism of the prompt emission are debated, we calculate the emission predicted by two models: magnetic dissipation and internal shocks. Our results favor the magnetic dissipation model in part because it predicts a relatively constant `Band' spectral peak energy E_peak with time during the GRB. The jet baryon loading decreases abruptly when the neutron star becomes transparent to neutrinos at t ~ 10-100 seconds. Jets with ultra-high magnetization cannot effectively accelerate and dissipate their energy, suggesting this transition ends the prompt emission and may explain the steep decay phase that follows. We assess several phenomena potentially related to magnetar birth, including low luminosity GRBs, thermal-rich GRBs/X-ray Flashes, very luminous supernovae, and short duration GRBs with extended emission.

preprint2010arXiv

A faint type of supernova from a white dwarf with a helium-rich companion

Supernovae (SNe) are thought to arise from two different physical processes. The cores of massive, short-lived stars undergo gravitational core collapse and typically eject a few solar masses during their explosion. These are thought to appear as as type Ib/c and II SNe, and are associated with young stellar populations. A type Ia SN is thought to arise from the thermonuclear detonation of a white dwarf star composed mainly of carbon and oxygen, whose mass approaches the Chandrasekhar limit. Such SNe are observed in both young and old stellar environments. Here we report our discovery of the faint type Ib SN 2005E in the halo of the nearby isolated galaxy, NGC 1032. The lack of any trace of recent star formation near the SN location (Fig. 1), and the very low derived ejected mass (~0.3 M_sun), argue strongly against a core-collapse origin for this event. Spectroscopic observations and the derived nucleosynthetic output show that the SN ejecta have high velocities and are dominated by helium-burning products, indicating that SN 2005E was neither a subluminous nor a regular SN Ia (Fig. 2). We have therefore found a new type of stellar explosion, arising from a low-mass, old stellar system, likely involving a binary with a primary white dwarf and a helium-rich secondary. The SN ejecta contain more calcium than observed in any known type of SN and likely additional large amounts of radioactive 44Ti. Such SNe may thus help resolve fundamental physical puzzles, extending from the composition of the primitive solar system and that of the oldest stars, to the Galactic production of positrons.

preprint2010arXiv

Electromagnetic Counterparts of Compact Object Mergers Powered by the Radioactive Decay of R-process Nuclei

The most promising astrophysical sources of kHz gravitational waves (GWs) are the inspiral and merger of binary neutron star(NS)/black hole systems. Maximizing the scientific return of a GW detection will require identifying a coincident electro-magnetic (EM) counterpart. One of the most likely sources of isotropic EM emission from compact object mergers is a supernova-like transient powered by the radioactive decay of heavy elements synthesized in ejecta from the merger. We present the first calculations of the optical transients from compact object mergers that self-consistently determine the radioactive heating by means of a nuclear reaction network; using this heating rate, we model the light curve with a one dimensional Monte Carlo radiation transfer calculation. For an ejecta mass ~1e-2 M_sun[1e-3 M_sun] the resulting light curve peaks on a timescale ~ 1 day at a V-band luminosity nu L_nu ~ 3e41[1e41] ergs/s (M_V = -15[-14]); this corresponds to an effective "f" parameter ~3e-6 in the Li-Paczynski toy model. We argue that these results are relatively insensitive to uncertainties in the relevant nuclear physics and to the precise early-time dynamics and ejecta composition. Due to the rapid evolution and low luminosity of NS merger transients, EM counterpart searches triggered by GW detections will require close collaboration between the GW and astronomical communities. NS merger transients may also be detectable following a short-duration Gamma-Ray Burst or "blindly" with present or upcoming optical transient surveys. Because the emission produced by NS merger ejecta is powered by the formation of rare r-process elements, current optical transient surveys can directly constrain the unknown origin of the heaviest elements in the Universe.

preprint2010arXiv

Gamma-Ray Burst Central Engines: Black Hole Vs. Magnetar

Discovered over forty years ago, Gamma-Ray Bursts (GRBs) remain a forefront topic in modern astrophysics. Perhaps the most fundamental question associated with GRBs is the nature of the astrophysical agent (or agents) that ultimately powers them: the central engine. In this review, I focus on the possible central engines of long-duration GRBs, and the constraints that present observations place on these models. Long GRBs are definitively associated with the deaths of massive stars, but whether the central engine is an accreting black hole or a rapidly-spinning, highly-magnetized neutron star (a "proto-magnetar") remains unsettled. This distinction has been brought into particular focus by recent MHD simulations of the core-collapse of massive, rotating "collapsar progenitors," which suggest that powerful magneto-centrifugal outflows from the proto-neutron star may stave off black hole formation entirely. Although both black hole and magnetar GRB models remain viable, I argue that the magnetar model is more mature in the sense that it provides quantitative explanations for the durations, energies, Lorentz factors, and collimation of long GRB outflows. Given these virtues, one promising strategy to break the present stalemate is to further develop the magnetar model until inescapable (and falsifiable) predictions emerge. This course of action signals a renewed challenge to translate time-dependent jet properties (power, magnetization, and Lorentz factor) into observables (gamma-ray light curves and spectra).

preprint2009arXiv

GRB 080503: Implications of a Naked Short Gamma-Ray Burst Dominated by Extended Emission

We report on observations of GRB 080503, a short gamma-ray burst with very bright extended emission (about 30 times the gamma-ray fluence of the initial spike) in conjunction with a thorough comparison to other short Swift events. In spite of the prompt-emission brightness, however, the optical counterpart is extraordinarily faint, never exceeding 25 mag in deep observations starting at ~1 hr after the BAT trigger. The optical brightness peaks at ~1 day and then falls sharply in a manner similar to the predictions of Li & Paczynski (1998) for supernova-like emission following compact-binary mergers. However, a shallow spectral index and similar evolution in X-rays inferred from Chandra observations are more consistent with an afterglow interpretation. The extreme faintness of this probable afterglow relative to the bright gamma-ray emission argues for a very low-density medium surrounding the burst (a "naked" GRB), consistent with the lack of a coincident host galaxy down to 28.5 mag in deep HST imaging. Our observations reinforce the notion that short GRBs generally occur outside regions of active star formation, but demonstrate that in some cases the luminosity of the extended prompt emission can greatly exceed that of the short spike, which may constrain theoretical interpretation of this class of events. Because most previous BAT short bursts without observed extended emission are too faint for this signature to have been detectable even if it were present at typical level, conclusions based solely on the observed presence or absence of extended emission in the existing Swift sample are premature. (abridged)

preprint2009arXiv

Magnetized Relativistic Jets and Long-Duration GRBs from Magnetar Spindown during Core-Collapse Supernovae

We use ideal axisymmetric relativistic magnetohydrodynamic simulations to calculate the spindown of a newly formed millisecond, B ~ 10^{15} G, magnetar and its interaction with the surrounding stellar envelope during a core-collapse supernova (SN) explosion. The mass, angular momentum, and rotational energy lost by the neutron star are determined self-consistently given the thermal properties of the cooling neutron star's atmosphere and the wind's interaction with the surrounding star. The magnetar drives a relativistic magnetized wind into a cavity created by the outgoing SN shock. For high spindown powers (~ 10^{51}-10^{52} ergs/s), the magnetar wind is super-fast at almost all latitudes, while for lower spindown powers (~ 10^{50} erg/s), the wind is sub-fast but still super-Alfvenic. In all cases, the rates at which the neutron star loses mass, angular momentum, and energy are very similar to the corresponding free wind values (<~ 30% differences), in spite of the causal contact between the neutron star and the stellar envelope. In addition, in all cases that we consider, the magnetar drives a collimated (~5-10 deg.) relativistic jet out along the rotation axis of the star. Nearly all of the spindown power of the neutron star escapes via this polar jet, rather than being transferred to the more spherical SN explosion. The properties of this relativistic jet and its expected late-time evolution in the magnetar model are broadly consistent with observations of long duration gamma-ray bursts (GRBs) and their associated broad-lined Type Ic SN.

preprint2009arXiv

Nickel-Rich Outflows from Accretion Disks Formed by the Accretion-Induced Collapse of White Dwarfs

A white dwarf (WD) approaching the Chandrasekhar mass may in several cases undergo accretion-induced collapse (AIC) to a neutron star (NS) before a thermonuclear explosion ensues. It has generally been assumed that AIC does not produce a detectable supernova (SN). If, however, the progenitor WD is rapidly rotating (as may be expected due to its prior accretion), a centrifugally supported disk forms around the NS upon collapse. We calculate the subsequent evolution of this accretion disk using time-dependent height-integrated simulations with initial conditions taken from the AIC calculations of Dessart et al. (2006). Initially, the disk is cooled by neutrinos and its composition is driven neutron-rich (electron fraction Ye ~ 0.1) by electron captures. However, as the disk viscously spreads, it is irradiated by neutrinos from the central proto-NS, which dramatically alters its neutron-to-proton ratio. We find that electron neutrino captures increase Ye to ~ 0.5 by the time that weak interactions in the disk freeze out. Because the disk becomes radiatively inefficient and begins forming alpha-particles soon after freeze out, powerful winds blow away most of the disk's remaining mass. These Ye ~ 0.5 outflows synthesize up to a few times 1e-2 Msun in 56Ni. As a result, AIC may be accompanied by a radioactively powered SN-like transient that peaks on a timescale of ~ 1 day. Since few intermediate mass elements are likely synthesized, these Ni-rich explosions should be spectroscopically distinct from other SNe. PanSTARRs and the Palomar Transient Factory should detect a few AIC transients per year if their true rate is ~1/100 of the Type Ia rate, and LSST should detect hundreds per year. High cadence observations (< 1 day) are optimal for the detection and follow-up of AIC (abridged).

preprint2009arXiv

The effects of r-process heating on fall-back accretion in compact object mergers

We explore the effects of r-process nucleosynthesis on fall-back accretion in neutron star(NS)-NS and black hole-NS mergers, and the resulting implications for short-duration gamma-ray bursts (GRBs). Though dynamically important, the energy released during the r-process is not yet taken into account in merger simulations. We use a nuclear reaction network to calculate the heating (due to beta-decays and nuclear fission) experienced by material on the marginally-bound orbits nominally responsible for late-time fall-back. Since matter with longer orbital periods t_orb experiences lower densities, for longer periods of time, the total r-process heating rises rapidly with t_orb, such that material with t_orb > 1 seconds can become completely unbound. Thus, r-process heating fundamentally changes the canonical prediction of an uninterrupted power-law decline in the fall-back rate dM/dt at late times. When the timescale for r-process to complete is > 1 second, the heating produces a complete cut-off in fall-back accretion after ~ 1 second; if robust, this would imply that fall-back accretion cannot explain the late-time X-ray flaring observed following some short GRBs. However, for a narrow, but physically plausible, range of parameters, fall-back accretion can resume after ~ 10 s, despite having been strongly suppressed for ~ 1-10 s after the merger. This suggests the intriguing possibility that the gap observed between the prompt and extended emission in short GRBs is a manifestation of r-process heating.

preprint2008arXiv

Neutron-Rich Freeze-Out in Viscously Spreading Accretion Disks Formed from Compact Object Mergers

Accretion disks with masses ~0.001-0.1 Msun form during the merger of neutron star (NS)-NS and black hole-NS binaries. Initially, such hyper-accreting disks cool efficiently by neutrino emission and their composition is driven neutron-rich by pair captures under degenerate conditions. However, as the disk viscously spreads and its temperature drops, cooling becomes inefficient and the disk becomes advective. Analytic arguments and numerical simulations suggest that once this occurs, powerful winds likely drive away most of the disk's remaining mass. We calculate the thermal evolution and nuclear composition of viscously spreading accretion disks formed from compact object mergers using one-dimensional height-integrated simulations. We show that freeze-out from weak equilibrium necessarily accompanies the disk's late-time transition to an advective state. As a result, hyper-accreting disks generically freeze out neutron-rich (with electron fraction Ye ~ 0.2-0.4), and their late-time outflows robustly synthesize rare neutron-rich isotopes. Using the measured abundances of these isotopes in our solar system, we constrain the compact object merger rate in the Milky Way to be < 1e-5 (M_d,0/0.1 Msun)^(-1) per year, where M_d,0 is the average initial mass of the accretion disk. Thus, either the NS-NS merger rate is at the low end of current estimates or the average disk mass produced during a typical merger is << 0.1 Msun. We also show that if most short duration gamma-ray bursts (GRBs) are produced by compact object mergers, their beaming fraction must exceed f_b ~ 0.13(M_d,0/0.1 Msun), corresponding to a jet half-opening angle > 30(M_d,0/0.1 Msun)^(1/2) degrees. This is consistent with other evidence that short duration GRB outflows are less collimated than those produced in long duration GRBs.