Massive, rapidly-spinning magnetar remnants produced as a result of binary neutron star (BNS) mergers may deposit a fraction of their energy into the surrounding kilonova ejecta, powering a synchrotron radio signal from the interaction of the ejecta with the circumburst medium. Here, we present 6.0 GHz Very Large Array (VLA) observations of nine, low-redshift short gamma-ray bursts (GRBs; $z<0.5$) on rest-frame timescales of $\approx 2.4-13.9$ yr following the bursts. We place $3\sigma$ limits on radio continuum emission of $F_{\nu}\lesssim 6-20\,\mu$Jy at the burst positions, or $L_{\nu} \lesssim (0.6-8.3) \times10^{28}$ erg s$^{-1}$ Hz$^{-1}$. Comparing these limits with new light curve modeling which properly incorporates relativistic effects, we obtain limits on the energy deposited into the ejecta of $E_{\rm ej}\lesssim (0.6-6.7)\times 10^{52}$ erg ($E_{\rm ej}\lesssim (1.7-14.1) \times 10^{52}$ erg) for an ejecta mass of $0.03\,M_{\odot}$ ($0.1\,M_{\odot}$). We present a uniform re-analysis of 27 short GRBs with 5.5-6.0 GHz observations, and find that $\gtrsim 50\%$ of short GRBs did not form stable magnetar remnants in their mergers. Assuming short GRBs are produced by BNS mergers drawn from the Galactic BNS population plus an additional component of high-mass GW194025-like mergers in a fraction $f_{\rm GW190425}$ of cases, we place constraints on the maximum mass of a non-rotating neutron star (NS) ($M_{\rm TOV}$), finding $M_{\rm TOV} \lesssim 2.4\,M_{\odot}$ for $f_{\rm GW190425}=0.4$. The detection (or lack thereof) of radio remnants in untargeted surveys such as the VLA Sky Survey (VLASS) could provide more stringent constraints on the fraction of mergers that produce stable remnants. If $\gtrsim 100-1000$ radio remnants are discovered in VLASS, this suggests that short GRBs are a biased population of BNS mergers in terms of the stability of the remnants they produce.

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

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.

Fast ejecta expelled in binary neutron star (NS) mergers or energetic supernovae (SNe) should produce late-time synchrotron radio emission as the ejecta shocks into the surrounding ambient medium. Models for such radio flares typically assume the ejecta expands into an unperturbed interstellar medium (ISM). However, it is also well-known that binary NS mergers and broad-lined Ic SNe can harbor relativistic jetted outflows. In this work, we show that such jets shock the ambient ISM ahead of the ejecta, thus evacuating the medium into which the ejecta subsequently collides. Using an idealized spherically-symmetric model, we illustrate that this inhibits the ejecta radio flare at early times $t < t_{\rm col} \approx 12 \, {\rm yr} \, (E_{\rm j}/10^{49} \, {\rm erg})^{1/3} (n/1 \, {\rm cm}^{-3})^{-1/3} (v_{\rm ej}/0.1c)^{-5/3}$ where $E_{\rm j}$ is the jet energy, $n$ the ISM density, and $v_{\rm ej}$ the ejecta velocity. We also show that this can produce a sharply peaked enhancement in the light-curve at $t = t_{\rm col}$. This has implications for radio observations of GW170817 and future binary NS mergers, gamma-ray burst (GRB) SNe, decade-long radio transients such as FIRST J1419, and possibly other events where a relativistic outflow precedes a slower-moving ejecta. Future numerical work will extend these analytic estimates and treat the multi-dimensional nature of the problem.

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

The properties of neutron stars are determined by the nature of the matter that they contain. These properties can be constrained by measurements of the star's size. We obtain stringent constraints on neutron-star radii by combining multimessenger observations of the binary neutron-star merger GW170817 with nuclear theory that best accounts for density-dependent uncertainties in the equation of state. We construct equations of state constrained by chiral effective field theory and marginalize over these using the gravitational-wave observations. Combining this with the electromagnetic observations of the merger remnant that imply the presence of a short-lived hyper-massive neutron star, we find that the radius of a $1.4\,\rm{M}_\odot$ neutron star is $R_{1.4\,\mathrm{M}_\odot} = 11.0^{+0.9}_{-0.6}~{\rm km}$ (90% credible interval). Using this constraint, we show that neutron stars are unlikely to be disrupted in neutron-star black-hole mergers; subsequently, such events will not produce observable electromagnetic emission.

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

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.

The electromagnetic (EM) signal of a binary neutron star (BNS) merger depends sensitively on the total binary mass, $M_{\rm tot}$, relative to various threshold masses set by the neutron star (NS) equation-of-state (EOS), parameterized through the NS maximum mass, $M_{\rm TOV}$, and characteristic radius, $R_{1.6}$. EM observations of a BNS merger detected through its gravitational wave (GW) emission, which are of sufficient quality to ascertain the identity of the merger remnant, can therefore constrain the values of $M_{\rm TOV}$ and $R_{1.6}$, given the tight connection between $M_{\rm tot}$ and the well-measured chirp mass. We elucidate the present and future landscape of EOS constraints from BNS mergers, introducing the "Multi-Messenger Matrix", a mapping between GW and EM measurables that defines the ranges of event chirp masses which provide the most leverage on constraining the EOS. By simulating a population of BNS mergers drawn from the Galactic double NS mass distribution we show that $\sim 10$ joint detections can constrain $M_{\rm TOV}$ and $R_{1.6}$ to several percent level where systematic uncertainties may become significant. Current EOS constraints imply that most mergers will produce supramassive or hypermassive remnants, a smaller minority (possibly zero) will undergo prompt-collapse, while at most only a few percent of events will form indefinitely stable NSs. In support of the envisioned program, we advocate in favor of LIGO/Virgo releasing chirp mass estimates as early as possible to the scientific community, enabling observational resources to be allocated in the most efficient way to maximize the scientific gain from multi-messenger discoveries.

We study mass ejection from accretion disks formed in the merger of a white dwarf with a neutron star or black hole. These disks are mostly radiatively-inefficient and support nuclear fusion reactions, with ensuing outflows and electromagnetic transients. Here we perform time-dependent, axisymmetric hydrodynamic simulations of these disks including a physical equation of state, viscous angular momentum transport, a coupled $19$-isotope nuclear network, and self-gravity. We find no detonations in any of the configurations studied. Our global models extend from the central object to radii much larger than the disk. We evolve these global models for several orbits, as well as alternate versions with an excised inner boundary to much longer times. We obtain robust outflows, with a broad velocity distribution in the range $10^2-10^4$ km s$^{-1}$. The outflow composition is mostly that of the initial white dwarf, with burning products mixed in at the $\lesssim 10-30\%$ level by mass, including up to $\sim 10^{-2}M_\odot$ of ${}^{56}$Ni. These heavier elements (plus ${}^{4}$He) are ejected within $\lesssim 40^\circ$ of the rotation axis, and should have higher average velocities than the lighter elements that make up the white dwarf. These results are in broad agreement with previous one- and two-dimensional studies, and point to these systems as progenitors of rapidly-rising ($\sim $ few day) transients. If accretion onto the central BH/NS powers a relativistic jet, these events could be accompanied by high energy transients with peak luminosities $\sim 10^{47}-10^{50}$ erg s$^{-1}$ and peak durations of up to several minutes, possibly accounting for events like CDF-S XT2.

The blossoming field of joint gravitational wave and electromagnetic (GW-EM) astronomy is one of the most promising in astronomy. The first, and only, joint GW-EM event GW170817 provided remarkable science returns that still continue to this day. Continued growth in this field requires increasing the sample size of joint GW-EM detections. In this white paper, we outline the case for using some percentage of LSST survey time for dedicated target-of-opportunity follow up of GW triggers in order to efficiently and rapidly identify optical counterparts. We show that the timeline for the LSST science survey is well matched to the planned improvements to ground based GW detectors in the next decade. LSST will become particularly crucial in the later half of the 2020s as more and more distant GW sources are detected. Lastly, we highlight some of the key science goals that can be addressed by a large sample of joint GW-EM detections.

Fast radio bursts (FRB) can arise from synchrotron maser emission at ultra-relativistic magnetized shocks, such as produced by flare ejecta from young magnetars. We combine PIC simulation results for the synchrotron maser with the dynamics of self-similar shock deceleration, as commonly applied to GRBs, to explore the implications for FRB emission. We assume the upstream environment into which the ultra-relativistic ejecta collides is a mildly relativistic baryon-loaded shell released following a previous flare, motivated by the high electron-ion injection rate Mdot ~ 1e19-1e21 g/s needed on larger scales to power the persistent radio nebula coincident with the repeating burster FRB 121102 and its high inferred rotation measure. The observed radio fluence peaks once the optical depth ahead of the shock to induced Compton scattering decreases to <~ few, a condition which places a GHz observer on the high frequency tail of the maser SED. Given intervals between ion shell ejection events ~1e5 s similar to the occurrence rate of the most powerful bursts from FRB 121102, we demonstrate the production of FRBs of frequency ~ 0.1-10 GHz, isotropic radiated energies ~1e37-1e40 erg and durations ~0.1-10 ms for flares of energy ~1e43-1e45 erg. Deceleration of the blast wave, and increasing transparency of the upstream medium, generates a temporal decay of the peak frequency, similar to the downward drift seen in the sub-bursts of FRB 121102 and FRB 180814.J0422+73. The delay >~ 1e5 s between major ion-injection events needed to clear sufficiently low densities around the engine for FRB emission could explain prolonged "dark" periods and clustered burst arrival times, and lead to stochastic variation in the dispersion measure. Thermal electrons heated at the shock generate a short-lived <~ 1 ms (1 s) synchrotron transient at gamma-ray (X-ray) energies, analogous to a GRB afterglow.

The last decade has seen the rapid, concurrent development of new classes of energetic astrophysical transients, including Fast Radio Bursts, Superluminous Supernovae, and ultralong gamma-ray bursts. We consider the implications of an emerging consensus on the origin of these diverse phenomena as different signatures of magnetar birth. We describe how magnetar-powered radio transients from milliseconds to gigaseconds will be used to address diverse physical questions related to the stellar evolution and death, the circum- and intergalactic medium, and fundamental physics on scales ranging from the quantum to cosmological. We argue that high-sensitivity radio interferometer concepts in development will be critical to accessing this science. There is also a growing need for scientists with expertise in new statistical tools and software techniques.

As of today, we have directly detected exactly one source in both gravitational waves (GWs) and electromagnetic (EM) radiation, the binary neutron star merger GW170817, its associated gamma-ray burst GRB170817A, and the subsequent kilonova SSS17a/AT 2017gfo. Within ten years, we will detect hundreds of events, including new classes of events such as neutron-star-black-hole mergers, core-collapse supernovae, and almost certainly something completely unexpected. As we build this sample, we will explore exotic astrophysical topics ranging from nucleosynthesis, stellar evolution, general relativity, high-energy astrophysics, nuclear matter, to cosmology. The discovery potential is extraordinary, and investments in this area will yield major scientific breakthroughs. Here we outline some of the most exciting scientific questions that can be answered by combining GW and EM observations.

The field of time-domain astrophysics has entered the era of Multi-messenger Astronomy (MMA). One key science goal for the next decade (and beyond) will be to characterize gravitational wave (GW) and neutrino sources using the next generation of Extremely Large Telescopes (ELTs). These studies will have a broad impact across astrophysics, informing our knowledge of the production and enrichment history of the heaviest chemical elements, constrain the dense matter equation of state, provide independent constraints on cosmology, increase our understanding of particle acceleration in shocks and jets, and study the lives of black holes in the universe. Future GW detectors will greatly improve their sensitivity during the coming decade, as will near-infrared telescopes capable of independently finding kilonovae from neutron star mergers. However, the electromagnetic counterparts to high-frequency (LIGO/Virgo band) GW sources will be distant and faint and thus demand ELT capabilities for characterization. ELTs will be important and necessary contributors to an advanced and complete multi-messenger network.

We present the detection of an unresolved radio source coincident with the position of the Type I superluminous supernova (SLSN) PTF10hgi ($z=0.098$) about 7.5 years post-explosion, with a luminosity of $L_\nu(6\,\,{\rm GHz})\approx 1.1\times 10^{28}$ erg s$^{-1}$ Hz$^{-1}$. This represents the first detection of radio emission coincident with a SLSN on any timescale. We investigate various scenarios for the origin of the radio emission: star formation activity, an active galactic nucleus, an off-axis jet, and a non-relativistic supernova blastwave. While any of these would be quite novel if confirmed, none appear likely when taken in context of the other properties of the host galaxy, previous radio observations of SLSNe, the sample of long gamma-ray bursts (LGRBs), and the general population of hydrogen-poor SNe. Instead, the radio emission is reminiscent of the quiescent radio source associated with the repeating FRB 121102, which has been argued to be powered by a magnetar born in a SLSN or LGRB explosion several decades ago. We show that such a central engine powered nebula is consistent with the age and luminosity of the radio source. Our directed search for FRBs from the location of PTF10hgi using 40 min of VLA phased-array data reveals no detections to a limit of $22$ mJy ($7\sigma$; 10 ms duration). We outline several follow-up observations that can conclusively establish the origin of the radio emission.

The fast radio burst FRB 121102 has repeated multiple times, enabling the identification of its host galaxy and of a spatially-coincident, compact, steady (`quiescent') radio synchrotron source. It was proposed that FRB 121102 is powered by a young flaring magnetar, embedded within a decades-old supernova remnant. Using a time-dependent one-zone model, we show that a single expanding magnetized electron-ion nebula (powered by the same outbursts likely responsible for the FRBs) can explain all the basic properties of the quiescent source (size, flux, self-absorption constraints) and the large but decreasing rotation measure (RM) of the bursts. The quiescent emission is powered by relativistic thermal electrons heated at the termination shock of the magnetar wind, while the RM originates from non-relativistic electrons injected earlier in the nebula's evolution and cooled through expansion and radiative losses. The model contains few free parameters, which are tightly constrained by observations: the total energy injected into the nebula over its history, $\sim 10^{50}-10^{51}$ erg, agrees with the magnetic energy of a millisecond magnetar; the baryon loading of the magnetar outflow (driven by intermittent flares) is close to the neutron star escape speed; the predicted source age $\sim 10-40$ years is consistent with other constraints on the nebula size. For an energy input rate $\dot{E} \propto t^{-\alpha}$ following the onset of magnetar activity, we predict secular decay of the RM and quiescent source flux, which approximately follow ${\rm RM} \propto t^{-(6+\alpha)/2}$ and $F_{\nu} \propto t^{-(\alpha^2+7\alpha-2)/4}$, respectively.

We present the first observations of a Type I superluminous supernova (SLSN) at $\gtrsim 1000$ days after maximum light. We observed SN 2015bn using the Hubble Space Telescope Advanced Camera for Surveys in the F475W, F625W and F775W filters at 721 days and 1068 days. SN 2015bn is clearly detected and resolved from its compact host, allowing reliable photometry. A galaxy template constructed from these data further enables us to isolate the SLSN flux in deep ground-based imaging. We measure a light curve decline rate at $>700$ days of $0.19 \pm 0.03$ mag (100 d)$^{-1}$, much shallower than the earlier evolution, and slower than previous SLSNe (at any phase) or the decay rate of $^{56}$Co. Neither additional radioactive isotopes nor a light echo can consistently account for the slow decline. A spectrum at 1083 days shows the same [O I] and [Ca II] lines as seen at $\sim300-400$ days, with no new features to indicate strong circumstellar interaction. Radio limits with the Very Large Array rule out an extended wind for mass-loss rates $10^{-2.7} \lesssim \dot{M}/v_{10} \lesssim 10^{-1.1}$ M$_\odot$ yr$^{-1}$ (where $v_{10}$ is the wind velocity in units of 10 km s$^{-1}$). The optical light curve is consistent with $L \propto t^{-4}$, which we show is expected for magnetar spin-down with inefficient trapping; furthermore, the evolution matches predictions from earlier magnetar model fits. The opacity to magnetar radiation is constrained at $\sim 0.01$ cm$^2$ g$^{-1}$, consistent with photon-matter pair-production over a broad $\sim$GeV-TeV range. This suggests the magnetar spectral energy distribution, and hence the 'missing energy' leaking from the ejecta, may peak in this range.

SN 2015bn is a nearby hydrogen-poor superluminous supernova (SLSN-I) that has been intensively observed in X-rays with the goal to detect the spin-down powered emission from a magnetar engine. The early-time UV/optical/infrared (UVOIR) data fit well to the magnetar model, but require leakage of energy at late times of $\lesssim 10^{43}$ erg s$^{-1}$, which is expected to be partially emitted in X-rays. Deep X-ray limits until $\sim$300 days after explosion revealed no X-ray emission. Here, we present the latest deep 0.3--10 keV X-ray limit at 805 days obtained with \textit{XMM-Newton}. We find $L_X < 10^{41}$ erg s$^{-1}$, with no direct evidence for central-engine powered emission. While the late-time optical data still follow the prediction of the magnetar model, the best-fit model to the bolometric light curve predicts that $\sim$97\% of the total input luminosity of the magnetar is escaping outside of the UVOIR bandpass at the time of observation. Our X-ray upper limit is $<$1.5\% of the input luminosity, strongly constraining the high-energy leakage, unless non-radiative losses are important. These deep X-ray observations identify a missing energy problem in SLSNe-I and we suggest future observations in hard X-rays and $\gamma$-rays for better constraints. Also, independent of the optical data, we constrain the parameter spaces of various X-ray emission scenarios, including ionization breakout by magnetar spin-down, shock interaction between the ejecta and external circumstellar medium, off-axis $\gamma$-ray burst afterglow, and black hole fallback accretion.

Young, rapidly spinning magnetars are invoked as central engines behind a diverse set of transient astrophysical phenomena, including gamma-ray bursts (GRB), super-luminous supernovae (SLSNe), fast radio bursts (FRB), and binary neutron star (NS) mergers. However, a barrier to direct confirmation of the magnetar hypothesis is the challenge of directly observing non-thermal emission from the central engine at early times (when it is most powerful and thus detectable) due to the dense surrounding ejecta. We present CLOUDY calculations of the time-dependent evolution of the temperature and ionization structure of expanding supernova or merger ejecta due to photo-ionization by a magnetar engine, in order to study the escape of X-rays (absorbed by neutral gas) and radio waves (absorbed by ionized gas), as well as to assess the evolution of the local dispersion measure due to photo-ionization. We find that ionization breakout does not occur if the engine's ionizing luminosity decays rapidly, and that X-rays typically escape the oxygen-rich ejecta of SLSNe only on $\sim 100 \, {\rm yr}$ timescales, consistent with current X-ray non-detections. We apply these results to constrain engine-driven models for the binary NS merger GW170817 and the luminous transient ASASSN-15lh. In terms of radio transparency and dispersion measure constraints, the repeating FRB 121102 is consistent with originating from a young, $\gtrsim 30-100 \, {\rm yr}$, magnetar similar to those inferred to power SLSNe. We further show that its high rotation measure can be produced within the same nebula that is proposed to power the quiescent radio source observed co-located with FRB 121102. Our results strengthen previous work suggesting that at least some FRBs may be produced by young magnetars, and motivate further study of engine powered transients.

We combine electromagnetic (EM) and gravitational wave (GW) information on the binary neutron star (NS) merger GW170817 in order to constrain the radii $R_{\rm ns}$ and maximum mass $M_{\rm max}$ of NSs. GW170817 was followed by a range of EM counterparts, including a weak gamma-ray burst (GRB), kilonova (KN) emission from the radioactive decay of the merger ejecta, and X-ray/radio emission consistent with being the synchrotron afterglow of a more powerful off-axis jet. The type of compact remnant produced in the immediate merger aftermath, and its predicted EM signal, depend sensitively on the high-density NS equation of state (EOS). For a soft EOS which supports a low $M_{\rm max}$, the merger undergoes a prompt collapse accompanied by a small quantity of shock-heated or disk wind ejecta, inconsistent with the large quantity $\gtrsim 10^{-2}M_{\odot}$ of lanthanide-free ejecta inferred from the KN. On the other hand, if $M_{\rm max}$ is sufficiently large, then the merger product is a rapidly-rotating supramassive NS (SMNS), which must spin-down before collapsing into a black hole. A fraction of the enormous rotational energy necessarily released by the SMNS during this process is transferred to the ejecta, either into the GRB jet (energy $E_{\rm GRB}$) or the KN ejecta (energy $E_{\rm ej}$), also inconsistent with observations. By combining the total binary mass of GW170817 inferred from the GW signal with conservative upper limits on $E_{\rm GRB}$ and $E_{\rm ej}$ from EM observations, we constrain the likelihood probability of a wide-range of previously-allowed EOS. These two constraints delineate an allowed region of the $M_{\rm max}-R_{\rm ns}$ parameter space, which once marginalized over NS radius places an upper limit of $M_{\rm max} \lesssim 2.17M_{\odot}$ (90\%), which is tighter or arguably less model-dependent than other current constraints.

Multiple observational lines of evidence support a connection between hydrogen-poor superluminous supernovae (SLSNe) and long duration gamma-ray bursts (GRBs). Both events require a powerful central energy source, usually attributed to a millisecond magnetar or an accreting black hole. The GRB-SLSN link raises several theoretical questions: What distinguishes the engines responsible for these different phenomena? Can a single engine power both a GRB and a luminous SN in the same event? We propose a new unifying model for magnetar thermalization and jet formation: misalignment between the rotation (${\bf \Omega}$) and magnetic dipole (${\bf \mu}$) axes thermalizes a fraction of the spindown power by reconnection in the striped equatorial wind, providing a guaranteed source of "thermal" emission to power the supernova. The remaining un-thermalized power energizes a relativistic jet. In this picture, the GRB-SLSN dichotomy is directly linked to ${\bf \Omega \cdot \mu}$. We extend earlier work to show that even weak relativistic jets of luminosity $\sim10^{46}$ erg s$^{-1}$ can escape the expanding SN ejecta hours after the explosion, implying that escaping relativistic jets may accompany many SLSNe. We calculate the observational signature of these jets. We show that they may produce transient UV cocoon emission lasting a few hours when the jet breaks out of the ejecta surface. A longer-lived optical/UV signal may originate from a mildly-relativistic wind driven from the interface between the jet and the ejecta walls. This provides a new explanation for the secondary early-time maximum observed in some SLSNe light curves, such as LSQ14bdq. This scenario also predicts a population of GRB from on-axis jets with extremely long durations, potentially similar to the population of "jetted tidal disruption events", in coincidence with a small subset of SLSNe.

Sub-arcsecond localization of the repeating fast radio burst FRB 121102 revealed its coincidence with a dwarf host galaxy and a steady (`quiescent') non-thermal radio source. We show that the properties of the host galaxy are consistent with those of long-duration gamma-ray bursts (LGRB) and hydrogen-poor superluminous supernovae (SLSNe-I). Both LGRBs and SLSNe-I were previously hypothesized to be powered by the electromagnetic spin-down of newly-formed, strongly-magnetized neutron stars with millisecond birth rotation periods (`millisecond magnetars'). This motivates considering a scenario whereby the repeated bursts from FRB 121102 originate from a young magnetar remnant embedded within a young hydrogen-poor supernova remnant. Requirements on the GHz free-free optical depth through the expanding supernova ejecta (accounting for photo-ionization by the rotationally-powered magnetar nebula), energetic constraints on the bursts, and constraints on the size of the quiescent source all point to an age of less than a few decades to a century. The quiescent radio source can be attributed to synchrotron emission from the shock interaction between the fast outer layer of the supernova ejecta with the surrounding wind of the progenitor star, or from deeper within the magnetar wind nebula. Alternatively, the radio emission could be an orphan afterglow from an initially off-axis LGRB jet, though this might require the source to be too young. The young age of the source can be tested by searching for a time derivative of the dispersion measure and predicted fading of the quiescent radio source. We propose future tests of the SLSNe-I/LGRB/FRB connection, such as searches for FRBs from nearby SLSNe-I/LGRB on timescales of decades after their explosions.

We show that the merger and tidal disruption of a C/O white dwarf (WD) by a neutron star (NS) binary companion provides a natural formation scenario for the PSR B1257+12 planetary system. Starting with initial conditions for the debris disk produced of the disrupted WD, we model its long term viscous evolution, including for the first time the effects of mass and angular momentum loss during the early radiatively inefficient accretion flow (RIAF) phase and accounting for the unusual C/O composition on the disk opacity. For plausible values of the disk viscosity $\alpha \sim 10^{-3}-10^{-2}$ and the RIAF mass loss efficiency, we find that the disk mass remaining near the planet formation radius at the time of solid condensation is sufficient to explain the pulsar planets. Rapid rocky planet formation via gravitational instability of the solid carbon-dominated disk is facilitated by the suppression of vertical shear instabilities due to the high solid-to-gas ratio. Additional evidence supporting a WD-NS merger scenario includes (1) the low observed occurrence rate of pulsar planets ($\lesssim 1\%$ of NS birth), comparable to the expected WD-NS merger rate; (2) accretion by the NS during the RIAF phase is sufficient to spin PSR B1257+12 up to its observed 6 ms period; (3) similar models of `low angular momentum' disks, such as those produced from supernova fallback, find insufficient mass reaching the planet formation radius. The unusually high space velocity of PSR B1257+12 of $\gtrsim 326\,{\rm km\,s}^{-1}$ suggests a possible connection to the Calcium-rich transients, dim supernovae which occur in the outskirts of their host galaxies and were proposed to result from mergers of WD-NS binaries receiving SN kicks. The C/O disk composition implied by our model likely results in carbon-rich planets with diamond interiors.

The distances of Fast Radio Burst (FRB) sources are currently unknown. We show that the 21-cm absorption line of hydrogen can be used to infer the redshifts of FRB sources, and determine whether they are Galactic or extragalactic. We calculate a probability of $\sim 10\%$ for the host galaxy of an FRB to exhibit a 21-cm absorption feature of equivalent width $\gtrsim 10 ~\mathrm{km}~\mathrm{s}^{-1}$. Arecibo, along with several future radio observatories, should be capable of detecting such associated 21-cm absorption signals for strong bursts of $\gtrsim {\rm several}~\mathrm{Jy}$ peak flux densities.

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

Stars with a different vertical motion relative to the galactic disk have a different average acceleration. According to Modified Newtonian Dynamics (MOND) theories they should therefore have a different average orbital velocity while revolving around the Milky Way. We show that this property can be used to constrain MOND theories by studying stars in the local neighborhood. With the {\sc Hipparcos} dataset we can only place marginal constraints. However, the forthcoming {\sc GAIA} catalogue with its significantly fainter cutoff should allow placing a stringent constraint.

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

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

The next generation gravitational waves (GW) detectors are most sensitive to GW emitted by compact (neutron star/black hole) binary mergers. If one of those is a neutron star the merger will also emit electromagnetic radiation via three possible channels: Gamma-ray bursts and their (possibly orphan) afterglows (Eichler et al. 1989), Li-Paczynski Macronovae (Li & Paczynski 1998) and radio flares (Nakar & Piran 2011). This accompanying electromagnetic radiation is vitally important in confirming the GW detections (Kochanek & Piran 1993). It could also reveal a wealth of information regarding the merger and will open a window towards multi-messenger astronomy. Identifying and characterizing these counterparts is therefore of utmost importance. In this work we explore late time radio flares emitted by the dynamically ejected outflows. We build upon previous work and consider the effect of the outflow's non-trivial geometry. Using an approximate method we estimate the radio light-curves for several ejected matter distributions obtained in numerical simulations. Our method provides an upper limit to the effect of non-sphericity. Together with the spherical estimates the resulting light curves bound the actual signal. We find that while non-spherical geometries can in principle lead to an enhanced emission, in most cases they result in an increase in the timescale compared with a corresponding spherical configuration. This would weaken somewhat these signals and might decrease the detection prospects.