Researcher profile

Edo Berger

Edo Berger contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
27works
0followers
9topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

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

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

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

Building this graph slice

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

Published work

27 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

The fate of the failed supernova candidate M31-2014-DS1

The fate of massive stars above 20M$_{\odot}$ remains uncertain. Debate persists about whether they die as supernovae (SNe), or if they collapse directly into black holes (BHs) with little or no optical outburst -- so-called ``failed supernovae''. The source M31-2014-DS1 experienced an optical outburst in 2014 and has remained faint at visual wavelengths since then. Due to its persistent faintness, it has been proposed as a failed SN candidate. We present new observations of this candidate obtained using the James Webb Space Telescope (JWST), the Submillimeter Array (SMA), and Chandra. The JWST observations demonstrate that a luminous mid-infrared source persists at the same location a decade after the star faded at visual wavelengths. We model its current spectral energy distribution (SED) as a dust-enshrouded star. No X-ray emission is detected, disfavoring the hypothesis that the late-time luminosity is powered by accretion onto a BH. We find that the remaining source is highly obscured by an asymmetric distribution of circumstellar dust, making it difficult to quantify its physical properties using spherically symmetric radiative transfer codes. The dust geometry requires that the inferred bolometric luminosity is only a lower limit, as a significant fraction of the central source's radiation may escape without being reprocessed by dust. We discuss the implications of these findings in the context of failed SN models and consider the potential overlap with signatures expected from a stellar merger, which also seems to provide a plausible explanation of this source.

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

Constraining the Time of Gravitational Wave Emission from Core-Collapse Supernovae

The advent of sensitive gravitational wave (GW) detectors, coupled with wide-field, high cadence optical time-domain surveys, raises the possibility of the first joint GW-electromagnetic (EM) detections of core-collapse supernovae (CCSNe). For targeted searches of GWs from CCSNe optical observations can be used to increase the sensitivity of the search by restricting the relevant time interval, defined here as the GW search window (GSW). The extent of the GSW is a critical factor in determining the achievable false alarm probability (FAP) for a triggered CCSN search. The ability to constrain the GSW from optical observations depends on how early a CCSN is detected, as well as the ability to model the early optical emission. Here we present several approaches to constrain the GSW, ranging in complexity from model-independent analytical fits of the early light curve, model-dependent fits of the rising or entire light curve, and a new data-driven approach using existing well-sampled CCSN light curves from {\it Kepler} and the Transiting Exoplanet Survey Satellite (TESS). We use these approaches to determine the time of core-collapse and its associated uncertainty (i.e., the GSW). We apply our methods to two Type II SNe that occurred during LIGO/Virgo Observing Run 3: SN\,2019fcn and SN\,2019ejj (both in the same galaxy at $d=15.7$ Mpc). Our approach shortens the duration of the GSW and improves the robustness of the GSW compared to techniques used in past GW CCSN searches.

preprint2022arXiv

Discovery of Three Candidate Magnetar-powered Fast X-ray Transients from Chandra Archival Data

It was proposed that a remnant stable magnetar could be formed in a binary neutron-star merger, leading to a fast X-ray transient (FXT) that can last for thousands of seconds. Recently, Xue et al. suggested that CDF-S XT2 was exactly such a kind of source. If confirmed, such emission can be used to search for electromagnetic counterparts to gravitational wave events from binary neutron-star mergers that have short gamma-ray bursts and the corresponding afterglows seen off-axis and thus too weak to be detected. Here we report the discovery of three new FXTs, XRT 170901, XRT 030511, and XRT 110919, from preliminary search over Chandra archival data. Similar to CDF-S XT2, these new FXTs had a very fast rise (less than a few ten seconds) and a plateau of X-ray flux $\sim$$1.0\times10^{-12}$ erg s$^{-1}$ cm$^{-2}$ lasting for 1-2 ks, followed by a steep decay. Their optical/IR counterparts, if present, are very weak, arguing against a stellar flare origin for these FXTs. For XRT 170901, we identified a faint host galaxy with the source at the outskirts, very similar to CDF-S XT2. Therefore, our newly discovered FXTs are also strong candidates for magnetar-powered X-ray transients resulting from binary neutron star mergers.

preprint2022arXiv

R2-D2: Roman and Rubin -- From Data to Discovery

The NASA Nancy Grace Roman Space Telescope (Roman) and the Vera C. Rubin Observatory Legacy Survey of Space and Time (Rubin), will transform our view of the wide-field sky, with similar sensitivities, but complementary in wavelength, spatial resolution, and time domain coverage. Here we present findings from the AURA Roman+Rubin Synergy Working group, charged by the STScI and NOIRLab Directors to identify frontier science questions in General Astrophysics, beyond the well-covered areas of Dark Energy and Cosmology, that can be uniquely addressed with Roman and Rubin synergies in observing strategy, data products and archiving, joint analysis, and community engagement. This analysis was conducted with input from the community in the form of brief (1-2 paragraph) "science pitches" (see Appendix), and testimony from "outside experts" (included as co-authors). We identify a rich and broad landscape of potential discoveries catalyzed by the combination of exceptional quality and quantity of Roman and Rubin data, and summarize implementation requirements that would facilitate this bounty of additional science with coordination of survey fields, joint coverage of the Galactic plane, bulge, and ecliptic, expansion of General Investigator and Target of Opportunity observing modes, co-location of Roman and Rubin data, and timely distribution of data, transient alerts, catalogs, value-added joint analysis products, and simulations to the broad astronomical community.

preprint2022arXiv

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

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

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

Deep Hubble Space Telescope Observations of GW170817: Complete Light Curves and the Properties of the Galaxy Merger of NGC 4993

We present the complete set of {\it Hubble Space Telescope} imaging of the binary neutron star merger GW170817 and its optical counterpart AT 2017gfo. Including deep template imaging in F814W, F110W, F140W, and F160W at 3.4 years post-merger, we re-analyze the full light curve of AT 2017gfo across 12 bands from 5--1273 rest-frame days after merger. We obtain four new detections of the short $γ$-ray burst (GRB) 170817A afterglow from 109--170 rest-frame days post-merger. These detections are consistent with the previously observed $β=-0.6$ spectral index in the afterglow light curve with no evidence for spectral evolution. We also analyze our limits in the case of novel late-time optical and IR emission signatures, such as a kilonova afterglow or infrared dust echo, but find our limits are not constraining in these contexts. We use the new data to construct deep optical and infrared stacks, reaching limits of $M=-6.3$ to $-4.6$ mag, to analyze the local environment around AT 2017gfo and low surface brightness features in its host galaxy NGC 4993. We rule out the presence of any globular cluster at the position of AT 2017gfo to $2.3 \times 10^{4} L_{\odot}$, including those with the reddest $V-H$ colors. Finally, we analyze the substructure of NGC 4993 in deep residual imaging, and find shell features which extend up to 71.8\arcsec\ (14.2 kpc) from the center of the galaxy. We find that the shells have a cumulative stellar mass of $6.3\times10^{8} M_{\odot}$, roughly 2% the total stellar mass of NGC 4993, and mass-weighted ages of $>$3 Gyr. We conclude that it was unlikely the GW170817 progenitor system formed in the galaxy merger, which based on dynamical signatures and the stellar population in the shells mostly likely occurred 220--685 Myr ago.

preprint2021arXiv

Impact of Massive Binary Star and Cosmic Evolution on Gravitational Wave Observations II: Double Compact Object Rates and Properties

Making the most of the rapidly increasing population of gravitational-wave detections of black hole (BH) and neutron star (NS) mergers requires comparing observations with population synthesis predictions. In this work we investigate the combined impact from the key uncertainties in population synthesis modelling of the isolated binary evolution channel: the physical processes in massive binary-star evolution and the star formation history as a function of metallicity, $Z$, and redshift $z, \mathcal{S}(Z,z)$. Considering these uncertainties we create 560 different publicly available model realizations and calculate the rate and distribution characteristics of detectable BHBH, BHNS, and NSNS mergers. We find that our stellar evolution and $\mathcal{S}(Z,z)$ variations can impact the predicted intrinsic and detectable merger rates by factors $10^2$-$10^4$. We find that BHBH rates are dominantly impacted by $\mathcal{S}(Z,z)$ variations, NSNS rates by stellar evolution variations and BHNS rates by both. We then consider the combined impact from all uncertainties considered in this work on the detectable mass distribution shapes (chirp mass, individual masses and mass ratio). We find that the BHNS mass distributions are predominantly impacted by massive binary-star evolution changes. For BHBH and NSNS we find that both uncertainties are important. We also find that the shape of the delay time and birth metallicity distributions are typically dominated by the choice of $\mathcal{S}(Z,z)$ for BHBH, BHNS and NSNS. We identify several examples of robust features in the mass distributions predicted by all 560 models, such that we expect more than 95% of BHBH detections to contain a BH $\gtrsim 8\,\rm{M}_{\odot}$ and have mass ratios $\lesssim 4$. Our work demonstrates that it is essential to consider a wide range of allowed models to study double compact object merger rates and properties.

preprint2021arXiv

Late-time Observations of Calcium-Rich Transient SN 2019ehk Reveal a Pure Radioactive Decay Power Source

We present $\textit{Hubble Space Telescope}$ imaging of the Calcium-rich supernova (SN) 2019ehk at 276 - 389 days after explosion. These observations represent the latest photometric measurements of a Calcium-rich transient to date and allows for the first opportunity to analyze the late-time evolution of an object in this observational SN class. We find that the late-time bolometric light curve of SN 2019ehk can be described predominantly through the radioactive decay of ${}^{56}\textrm{Co}$ for which we derive a mass of $M({}^{56}\textrm{Co}) = (2.8 \pm 0.1) \times 10^{-2}$$\rm{M}_\odot$. Furthermore, the rate of decline in bolometric luminosity requires the leakage of $γ$-rays on timescale $t_γ = 53.9 \pm 1.30$ days, but we find no statistical evidence for incomplete positron trapping in the SN ejecta. While our observations cannot constrain the exact masses of other radioactive isotopes synthesized in SN 2019ehk, we estimate a mass ratio limit of $M({}^{57}\textrm{Co}) / M({}^{56}\textrm{Co}) \leq 0.030$. This limit is consistent with the explosive nucleosynthesis produced in the merger of low-mass white dwarfs, which is one of the favored progenitor scenarios in early-time studies of SN 2019ehk.

preprint2021arXiv

Optical Observations and Modeling of the Superluminous Supernova 2018lfe

We present optical imaging and spectroscopy of SN\,2018lfe, which we classify as a Type I superluminous supernova (SLSN-I) at a redshift of $z = 0.3501$ with a peak absolute magnitude of $M_r\approx -22.1$ mag, one of the brightest SLSNe discovered. SN\,2018lfe was identified for follow-up using our FLEET machine learning pipeline. Both the light curve and the spectra of SN\,2018lfe are consistent with the broad population of SLSNe. We fit the light curve with a magnetar central engine model and find an ejecta mass of $M_{\rm ej}\approx 3.8$ M$_\odot$, a magnetar spin period of $P\approx 2.9$ ms and a magnetic field strength of $B_{\perp}\approx 2.8\times 10^{14}$ G. The magnetic field strength is near the top of the distribution for SLSNe, while the spin period and ejecta mass are near the median values of the distribution for SLSNe. From late-time imaging and spectroscopy we find that the host galaxy of SN\,2018lfe has an absolute magnitude of $M_r\approx -17.85$, ($L_B \approx 0.029$ $L^*$), and an inferred metallicity of $Z\approx 0.3$ Z$_\odot$, and star formation rate of $\approx 0.8$ M$_\odot$ yr$^{-1}$.

preprint2020arXiv

An extremely energetic supernova from a very massive star in a dense medium

The interaction of a supernova with a circumstellar medium (CSM) can dramatically increase the emitted luminosity by converting kinetic energy to thermal energy. In 'superluminous' supernovae (SLSNe) of Type IIn -- named for narrow hydrogen lines in their spectra -- the integrated emission can reach $\sim 10^{51}$ erg, attainable by thermalising most of the kinetic energy of a conventional SN. A few transients in the centres of active galaxies have shown similar spectra and even larger energies, but are difficult to distinguish from accretion onto the supermassive black hole. Here we present a new event, SN2016aps, offset from the centre of a low-mass galaxy, that radiated $\gtrsim 5 \times 10^{51}$ erg, necessitating a hyper-energetic supernova explosion. We find a total (SN ejecta $+$ CSM) mass likely exceeding 50-100 M$_\odot$, with energy $\gtrsim 10^{52}$ erg, consistent with some models of pair-instability supernovae (PISNe) or pulsational PISNe -- theoretically-predicted thermonuclear explosions from helium cores $>50$ M$_\odot$. Independent of the explosion mechanism, this event demonstrates the existence of extremely energetic stellar explosions, detectable at very high redshifts, and provides insight into dense CSM formation in the most massive stars.

preprint2020arXiv

Does GW190425 require an alternative formation pathway than a fast-merging channel?

The LIGO/Virgo Scientific Collaboration (LVC) recently announced the detection of a compact object binary merger, GW190425, with a total mass of $3.4^{+0.3}_{-0.1}$ M$_{\odot}$, and individual component masses in the range of about 1.1 to 2.5 $M_{\odot}$. If the constituent compact objects are neutron stars, then the total mass is five standard deviations higher than the mean of $2.66\pm 0.12$ M$_{\odot}$ for Galactic binary neutron stars (BNSs). LVC suggests such massive BNS systems are born from a fast-merging channel, and therefore, their non-detection in the Galaxy to be due to a selection effect. However, we are unable to reconcile the inferred formation efficiency from the reported merger rate, $\mathcal{R}_{\rm GW190425}=460^{+1050}_{-390}$ yr$^{-1}$ Gpc$^{-3}$, with predictions from our own study for fast-merging BNS systems. Moreover, the comparable merger rates of GW190425, and GW170817 are possibly in tension with our results for two reasons: (i) more massive systems are expected to have a lower formation rate, and (ii) fast merging channels should constitute $\lesssim 10\%$ of the total BNS systems if case BB unstable mass transfer is permitted to take place as a formation pathway. We argue that to account for the high merger rate of GW190425 as a BNS system requires: (i) revisiting our understanding of NS formation in supernova explosions, or (ii): that more massive NSs need to be preferentially born with either very weak or very high magnetic fields so that they would be undetectable in the radio surveys. Perhaps massive NSs detected in NS-white dwarf binaries are our clues to the formation path of GW190425 systems.

preprint2020arXiv

Radio Monitoring of the Tidal Disruption Event Swift J164449.3+573451. IV. Continued Fading and Non-Relativistic Expansion

We present continued radio and X-ray observations of the previously relativistic tidal disruption event (TDE) Swift J164449.3+573451 (\sw) extending to about 9.4 years post disruption, as part of ongoing campaigns with the Jansky Very Large Array (VLA) and the \textit{Chandra} X-ray observatory. We find that the X-ray emission has faded below detectable levels, with an upper limit of $\lesssim 3.5\times 10^{-15}$ erg cm$^{-2}$ s$^{-1}$ in a 100 ks observation, while the radio emission continues to be detected and steadily fade. Both are consistent with forward shock emission from a non-relativistic outflow, although we find that the radio spectral energy distribution is better fit at these late times with an electron power law index of $p\approx 3$ (as opposed to $p\approx 2.5$ at earlier times). With the revised spectral index we find $ε_B\approx 0.01$ using the radio and X-ray data, and a density of $\approx 0.04$ cm$^{3}$ at a radius of $R\approx 0.65$ pc ($R_{\rm sch}\approx 2\times 10^6$ R$_\odot$) from the black hole. The energy scale of the blastwave is $\approx 10^{52}$ erg. We also report detections of \sw\ at 3 GHz from the first two epochs of the VLA Sky Survey (VLASS), and find that $\sim 10^2$ off-axis \sw-like events to $z\sim 0.5$ may be present in the VLASS data. Finally, we find that \sw\ itself will remain detectable for decades at radio frequencies, although observations at sub-GHz frequencies will become increasingly important to characterize its dynamical evolution.

preprint2020arXiv

SN 2019ehk: A Double-Peaked Ca-rich Transient with Luminous X-ray Emission and Shock-Ionized Spectral Features

We present panchromatic observations and modeling of the Calcium-rich supernova 2019ehk in the star-forming galaxy M100 (d$\approx$16.2 Mpc) starting 10 hours after explosion and continuing for ~300 days. SN 2019ehk shows a double-peaked optical light curve peaking at $t = 3$ and $15$ days. The first peak is coincident with luminous, rapidly decaying $\textit{Swift}$-XRT discovered X-ray emission ($L_x\approx10^{41}~\rm{erg~s^{-1}}$ at 3 days; $L_x \propto t^{-3}$), and a Shane/Kast spectral detection of narrow H$α$ and He II emission lines ($v \approx 500$ km/s) originating from pre-existent circumstellar material. We attribute this phenomenology to radiation from shock interaction with extended, dense material surrounding the progenitor star at $r<10^{15}$ cm and the resulting cooling emission. We calculate a total CSM mass of $\sim$ $7\times10^{-3}$ $\rm{M_{\odot}}$ with particle density $n\approx10^{9}\,\rm{cm^{-3}}$. Radio observations indicate a significantly lower density $n < 10^{4}\,\rm{cm^{-3}}$ at larger radii. The photometric and spectroscopic properties during the second light curve peak are consistent with those of Ca-rich transients (rise-time of $t_r =13.4\pm0.210$ days and a peak B-band magnitude of $M_B =-15.1\pm0.200$ mag). We find that SN 2019ehk synthesized $(3.1\pm0.11)\times10^{-2} ~ \rm{M_{\odot}}$ of ${}^{56}\textrm{Ni}$ and ejected $M_{\rm ej} = (0.72\pm 0.040)~\rm{M_{\odot}}$ total with a kinetic energy $E_{\rm k}=(1.8\pm0.10)\times10^{50}~\rm{erg}$. Finally, deep $\textit{HST}$ pre-explosion imaging at the SN site constrains the parameter space of viable stellar progenitors to massive stars in the lowest mass bin (~10 $\rm{M_{\odot}}$) in binaries that lost most of their He envelope or white dwarfs. The explosion and environment properties of SN 2019ehk further restrict the potential WD progenitor systems to low-mass hybrid HeCO WD + CO WD binaries.

preprint2020arXiv

The Pre-Explosion Mass Distribution of Hydrogen-Poor Superluminous Supernova Progenitors and New Evidence for a Mass-Spin Correlation

Despite indications that superluminous supernovae (SLSNe) originate from massive progenitors, the lack of a uniformly analyzed statistical sample has so far prevented a detailed view of the progenitor mass distribution. Here we present and analyze the pre-explosion mass distribution of hydrogen-poor SLSN progenitors as determined from uniformly modelled light curves of 62 events. We construct the distribution by summing the ejecta mass posteriors of each event, using magnetar light curve models presented in our previous works (and using a nominal neutron star remnant mass). The resulting distribution spans $3.6-40$ M$_{\odot}$, with a sharp decline at lower masses, and is best fit by a broken power law described by ${\rm d}N/{\rm dlog}M \propto M^{-0.41 \pm 0.06}$ at $3.6-8.6$ M$_{\odot}$ and $\propto M^{-1.26 \pm 0.06}$ at $8.6-40$ M$_{\odot}$. We find that observational selection effects cannot account for the shape of the distribution. Relative to Type Ib/c SNe, the SLSN mass distribution extends to much larger masses and has a different power-law shape, likely indicating that the formation of a magnetar allows more massive stars to explode as some of the rotational energy accelerates the ejecta. Comparing the SLSN distribution with predictions from single and binary star evolution models, we find that binary models for a metallicity of $Z\lesssim 1/3$ Z$_{\odot}$ are best able to reproduce its broad shape, in agreement with the preference of SLSNe for low metallicity environments. Finally, we uncover a correlation between the pre-explosion mass and the magnetar initial spin period, where SLSNe with low masses have slower spins, a trend broadly consistent with the effects of angular momentum transport evident in models of rapidly-rotating carbon-oxygen stars.

preprint2020arXiv

The Tidal Disruption Event AT 2018hyz II: Light Curve Modeling of a Partially Disrupted Star

AT 2018hyz (=ASASSN-18zj) is a tidal disruption event (TDE) located in the nucleus of a quiescent E+A galaxy at a redshift of $z = 0.04573$, first detected by the All-Sky Automated Survey for Supernovae (ASAS-SN). We present optical+UV photometry of the transient, as well as an X-ray spectrum and radio upper limits. The bolometric light curve of AT 2018hyz is comparable to other known TDEs and declines at a rate consistent with a $t^{-5/3}$ at early times, emitting a total radiated energy of $E = 9\times10^{50}$ erg. An excess bump appears in the UV light curve about 50 days after bolometric peak, followed by a flattening beyond 250 days. The light curve shows an excess bump in the UV about 50 days after bolometric peak lasting for at least 100 days, which may be related to an outflow. We detect a constant X-ray source present for at least 86 days. The X-ray spectrum shows a total unabsorbed flux of $\sim 4\times10^{-14}$ erg cm$^{-2}$ s$^{-1}$ and is best fit by a blackbody plus power-law model with a photon index of $Γ= 0.8$. A thermal X-ray model is unable to account for photons $> 1$ keV, while the radio non-detection favors inverse-Compton scattering rather than a jet for the non-thermal component. We model the optical and UV light curves using the Modular Open-Source Fitter for Transients (MOSFiT) and find a best fit for a black hole of $5.2\times10^6$ M$_\odot$ partially disrupting a $0.1$ M$_\odot$ star (stripping a mass of $\sim 0.01$ M$_\odot$ for the inferred impact parameter, $β=0.6$). The low optical depth implied by the small debris mass may explain how we are able to see hydrogen emission with disk-like line profiles in the spectra of AT 2018hyz (see our companion paper, Short et al.~2020).

preprint2019arXiv

AT 2018cow VLBI: No Long-Lived Relativistic Outflow

We report on VLBI observations of the fast and blue optical transient (FBOT), AT 2018cow. At ~62 Mpc, AT 2018cow is the first relatively nearby FBOT. The nature of AT 2018cow is not clear, although various hypotheses from a tidal disruption event to different kinds of supernovae have been suggested. It had a very fast rise time (3.5 d) and an almost featureless blue spectrum although high photospheric velocities (40,000 km s$^{-1}$) were suggested early on. The X-ray luminosity was very high, ~$1.4 \times 10^{43}$ erg s$^{-1}$, larger than those of ordinary SNe, and more consistent with those of SNe associated with gamma-ray bursts. Variable hard X-ray emission hints at a long-lived &#34;central engine.&#34; It was also fairly radio luminous, with a peak 8.4-GHz spectral luminosity of ~$4 \times 10^{28}$ erg s$^{-1}$ Hz$^{-1}$, allowing us to make VLBI observations at ages between 22 and 287 d. We do not resolve AT 2018cow. Assuming a circularly symmetric source, our observations constrain the average apparent expansion velocity to be <0.49$c$ by t = 98 d (3$σ$ limit). We also constrain the proper motion of AT 2018cow to be <0.51$c$. Since the radio emission generally traces the fastest ejecta, our observations make the presence of a long-lived relativistic jet with a lifetime of more than one month very unlikely.

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

Formation and merging of Mass Gap Black Holes in Gravitational Wave Merger Events from Wide Hierarchical Quadruple Systems

We investigate secular evolution in hierarchical quadruple systems as a formation channel of mass-gap black holes (with masses of about $3-5$ $M_{\odot}$) in systems that will eventually lead to binary black hole mergers detectable by ground-based gravitational wave detectors (LIGO/Virgo). We show that in a 3+1 hierarchical system, two episodes of induced mergers would first cause two neutron stars to merge and form a mass-gap black hole, which will subsequently merge with another (more massive) black hole through a second induced merger. We demonstrate that such systems are stable to flybys, and their formation would predict a high mass ratio and eccentric merger of a mass-gap black hole with a more massive black hole companion. Such a formation channel may explain observed gravitational wave events such as the recently-discovered LIGO/Virgo events S190814bv and S190924h.

preprint2011arXiv

Exploring The Galaxy Mass-Metallicity Relation at 3 < z < 5

Long-duration gamma-ray bursts (GRBs) provide a premier tool for studying high-redshift star-forming galaxies thanks to their extreme brightness and association with massive stars. Here we use GRBs to study the galaxy mass-metallicity (M*-Z) relation at z~3-5, where conventional direct metallicity measurements are extremely challenging. We use the interstellar medium metallicities of long-duration GRB hosts derived from afterglow absorption pectroscopy (Z~0.01-1 solar), in conjunction with host galaxy stellar masses determined from deep Spitzer 3.6 micron observations of 20 GRB hosts. We detect about 1/4 of the hosts with I-band absolute magnitudes, M ~ -21.5 to -22.5 AB mag, and place a limit of M>-19 mag on the remaining hosts from a stacking analysis. Using a conservative range of mass-to-light ratios for simple stellar populations (with ages of 70 Myr to ~ 2 Gyr), we infer the host stellar masses and present the galaxy mass-metallicity measurements at z ~ 3-5 (<z> ~ 3.5). We find that the detected GRB hosts, with M* ~ 2e10 solar masses, display a wide range of metallicities, but that the mean metallicity at this mass scale, Z~ 0.1 solar, is lower than measurements at z < 3. Combined with stacking of the non-detected hosts (with M* < 4e9 solar masses and Z < 0.03 solar), we find evidence for the existence of an M*-Z relation at z ~ 3.5 and continued evolution of this relation to systematically lower metallicities from z ~ 2.

preprint2010arXiv

Spectroscopic Discovery of the Broad-Lined Type Ic Supernova 2010bh Associated with the Low-Redshift GRB 100316D

We present the spectroscopic discovery of a broad-lined Type Ic supernova (SN 2010bh) associated with the nearby long-duration gamma-ray burst (GRB) 100316D. At z = 0.0593, this is the third-nearest GRB-SN. Nightly optical spectra obtained with the Magellan telescopes during the first week after explosion reveal the gradual emergence of very broad spectral features superposed on a blue continuum. The supernova features are typical of broad-lined SNe Ic and are generally consistent with previous supernovae associated with low-redshift GRBs. However, the inferred velocities of SN 2010bh at 21 days after explosion are a factor of ~2 times larger than those of the prototypical SN 1998bw at similar epochs, with v ~ 26,000 km/s, indicating a larger explosion energy or a different ejecta structure. A near-infrared spectrum taken 13.8 days after explosion shows no strong evidence for He I at 1.083 microns, implying that the progenitor was largely stripped of its helium envelope. The host galaxy is of low luminosity (M_R ~ -18.5 mag) and low metallicity (Z < 0.4 Z_solar), similar to the hosts of other low-redshift GRB-SNe.

preprint2009arXiv

On The Origin Of The Highest Redshift Gamma-Ray Bursts

GRB 080913 and GRB 090423 are the most distant gamma-ray bursts (GRBs) known to-date, with spectroscopically determined redshifts of z=6.7 and z=8.1, respectively. The detection of bursts at this early epoch of the Universe significantly constrains the nature of GRBs and their progenitors. We perform population synthesis studies of the formation and evolution of early stars, and calculate the resulting formation rates of short and long-duration GRBs at high redshift. The peak of the GRB rate from Population II stars occurs at z=7 for a model with efficient/fast mixing of metals, while it is found at z=3 for an inefficient/slow metallicity evolution model. We show that in the redshift range 6<z<10 essentially all GRBs originate from Population II stars, regardless of metallicity evolution model. These stars (having small, but non-zero metallicity) are the most likely progenitors for both long GRBs (collapsars) and short GRBs (NS-NS or BH-NS mergers) at this epoch. Although the predicted intrinsic rates of long and short GRBs are similar at these high redshifts, observational selection effects lead to higher (factor of 10) observed rates for long GRBs. We conclude that the two recently observed high-z GRB events are most likely long GRBs originating from Population II collapsars.

preprint2009arXiv

Short-Term H-alpha Variability in M Dwarfs

We spectroscopically study the variability of H-alpha emission in mid- to late-M dwarfs on timescales of ~0.1-1 hr as a proxy for magnetic variability. About 80% of our sample exhibits statistically significant variability on the full range of timescales probed by the observations, and with amplitude ratios in the range of ~1.2-4. No events with an order of magnitude increase in H-alpha luminosity were detected, indicating that their rate is < 0.05 /hr (95% confidence level). We find a clear increase in variability with later spectral type, despite an overall decrease in H-alpha &#34;activity&#34; (i.e., L_{H-alpha}/L_{bol}). For the ensemble of H-alpha variability events, we find a nearly order of magnitude increase in the number of events from timescales of about 10 to 30 min, followed by a roughly uniform distribution at longer durations. The event amplitudes follow an exponential distribution with a characteristic scale of Max(EW)/Min(EW)}-1 ~ 0.7. This distribution predicts a low rate of ~ 10^{-6} /hr for events with (Max(EW)/Min(EW) > 10, but serendipitous detections of such events in the past suggests that they represent a different distribution. Finally, we find a possible decline in the amplitude of events with durations of > 0.5 hr, which may point to a typical energy release in H-alpha events for each spectral type (E_{H-alpha} ~ L_{H-alpha} x t ~ const). Longer observations of individual active objects are required to further investigate this possibility. Similarly, a larger sample may shed light on whether H-alpha variability correlates with properties such as age or rotation velocity.

preprint2008arXiv

A comprehensive study of GRB 070125, a most energetic gamma ray burst

We present a comprehensive multiwavelength analysis of the bright, long duration gamma-ray burst GRB 070125, comprised of observations in $γ$-ray, X-ray, optical, millimeter and centimeter wavebands. Simultaneous fits to the optical and X-ray light curves favor a break on day 3.78, which we interpret as the jet break from a collimated outflow. Independent fits to optical and X-ray bands give similar results in the optical bands but shift the jet break to around day 10 in the X-ray light curve. We show that for the physical parameters derived for GRB 070125, inverse Compton scattering effects are important throughout the afterglow evolution. While inverse Compton scattering does not affect radio and optical bands, it may be a promising candidate to delay the jet break in the X-ray band. Radio light curves show rapid flux variations, which are interpreted as due to interstellar scintillation, and are used to derive an upper limit of $2.4 \times 10^{17}$ cm on the radius of the fireball in the lateral expansion phase of the jet. Radio light curves and spectra suggest a high synchrotron self absorption frequency indicative of the afterglow shock wave moving in a dense medium. Our broadband modeling favors a constant density profile for the circumburst medium over a wind-like profile ($R^{-2}$). However, keeping in mind the uncertainty of the parameters, it is difficult to unambiguously distinguish between the two density profiles. Our broadband fits suggest that \event is a burst with high radiative efficiency ($> 60 %$).

preprint2007arXiv

Galaxy Clusters Associated with Short GRBs. II. Predictions for the Rate of Short GRBs in Field and Cluster Early-Type Galaxies

We determine the relative rates of short GRBs in cluster and field early-type galaxies as a function of the age probability distribution of their progenitors, P(τ) \propto τ^n. This analysis takes advantage of the difference in the growth of stellar mass in clusters and in the field, which arises from the combined effects of the galaxy stellar mass function, the early-type fraction, and the dependence of star formation history on mass and environment. This approach complements the use of the early- to late-type host galaxy ratio, with the added benefit that the star formation histories of early-type galaxies are simpler than those of late-type galaxies, and any systematic differences between progenitors in early- and late-type galaxies are removed. We find that the ratio varies from R(cluster)/R(field) ~ 0.5 for n = -2 to ~ 3 for n = 2. Current observations indicate a ratio of about 2, corresponding to n ~ 0 - 1. This is similar to the value inferred from the ratio of short GRBs in early- and late-type hosts, but it differs from the value of n ~ -1 for NS binaries in the Milky Way. We stress that this general approach can be easily modified with improved knowledge of the effects of environment and mass on the build-up of stellar mass, as well as the effect of globular clusters on the short GRB rate. It can also be used to assess the age distribution of Type Ia supernova progenitors.