Abstract
Dense circumstellar material (CSM) is thought to play an important role in observed luminous optical transients: if such CSM is shocked, e.g., by ejecta expelled from the progenitor during ...core-collapse, then radiation produced by the shock-heated CSM can power bright UV/optical emission. If the initial CSM has an “outer edge” where most of the mass is contained and at which the optical depth is large, then shock breakout—when photons are first able to escape the shocked CSM—occurs near it. The rather thin shell of shocked CSM subsequently expands, and in the ensuing cooling-envelope phase, radiative and adiabatic losses compete to expend the CSM thermal energy. Here we derive an analytic solution to the bolometric light curve produced by such shocked CSM. For the first time, we provide an analytic solution to the cooling-envelope phase that is applicable starting from shock breakout and until the expanding CSM becomes optically thin. In particular, we account for the planar CSM geometry that is relevant at early times and properly treat radiative losses within this planar phase. We show that these effects can dramatically impact the resulting light curves, particularly if the CSM optical depth is only marginally larger than
c
/
v
sh
(where
v
sh
is the shock velocity). This has important implications for interpreting observed fast optical transients, which have previously been modeled using either computationally expensive numerical simulations or more simplified models that do not properly capture the early light-curve evolution.
We combine electromagnetic (EM) and gravitational-wave (GW) information on the binary neutron star (NS) merger GW170817 in order to constrain the radii and maximum mass 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 that supports a low , the merger undergoes a prompt collapse accompanied by a small quantity of shock-heated or disk-wind ejecta, inconsistent with the large quantity of lanthanide-free ejecta inferred from the KN. On the other hand, if 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 ) or the KN ejecta (energy ), also inconsistent with observations. By combining the total binary mass of GW170817 inferred from the GW signal with conservative upper limits on and from EM observations, we constrain the likelihood probability of a wide range of previously allowed EOSs. These two constraints delineate an allowed region of the parameter space, which, once marginalized over NS radius, places an upper limit of (90%), which is tighter or arguably less model-dependent than other current constraints.
ABSTRACT
Fast radio bursts (FRBs) can arise from synchrotron maser emission at ultrarelativistic magnetized shocks, such as produced by flare ejecta from young magnetars. We combine particle-in-cell ...simulation results for the maser emission with the dynamics of self-similar shock deceleration, as commonly applied to gamma-ray bursts (GRBs), to explore the implications for FRBs. The upstream environment is a mildly relativistic baryon-loaded shell released following a previous flare, motivated by the high electron–ion injection rate $\dot{M} \sim 10^{19}\!-\!10^{21}$ g s−1 needed to power the persistent radio nebula coincident with the repeating burster FRB 121102 and its high rotation measure. The radio fluence peaks once the optical depth ahead of the shock to induced Compton scattering τc ≲ 3. Given intervals between major ion ejection events ΔT ∼ 105 s similar to the occurrence rate of the most powerful bursts from FRB 121102, we demonstrate the production of ∼0.1–10 GHz FRBs with isotropic radiated energies ∼1037–1040 erg and durations ∼0.1–10 ms for flare energies E ∼ 1043–1045 erg. Deceleration of the blast wave, and increasing transparency of the upstream medium, generates temporal decay of the peak frequency, similar to the observed downward frequency drift seen in FRB 121102 and FRB 180814.J0422+73. The delay ΔT ≳ 105 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. 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 scaled-down GRB afterglow.
The fast radio burst FRB 121102 has repeated multiple times, enabling the identification of its host galaxy and of a spatially coincident, compact, steady ("persistent") 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 (created by the same outbursts likely responsible for the fast radio bursts) can explain all of the basic properties of the persistent source (size, flux, self-absorption constraints) and the large but decreasing rotation measure (RM) of the bursts. The persistent 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, ∼1050−1051 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 ∼10-40 yr is consistent with other constraints on the nebula size. For an energy input rate following the onset of magnetar activity, we predict secular decay of the RM and persistent source flux, which approximately follow RM ∝ t−(6+ )/2 and , respectively.
Subarcsecond localization of the repeating fast radio burst FRB 121102 revealed its coincidence with a dwarf host galaxy and a steady ("quiescent") nonthermal 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 (SN) remnant. Requirements on the gigahertz free-free optical depth through the expanding SN ejecta (accounting for photoionization 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. 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 the radio source can from deeper within the magnetar wind nebula as outlined in Metzger et al. 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 the 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/LGRBs on timescales of decades after their explosions.
The electromagnetic (EM) signal of a binary neutron star (BNS) merger depends sensitively on the total binary mass, Mtot, relative to various threshold masses set by the neutron star (NS) equation of ...state (EOS), parameterized through the neutron star (NS) maximum mass, MTOV, and characteristic radius, R1.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 MTOV and R1.6, given the tight connection between Mtot 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 that 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 ∼10 joint detections can constrain MTOV and R1.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 Laser Interferometer Gravitational-Wave Observatory (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.
ABSTRACT
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 Ic can harbour 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 \lt t_{\rm col} \approx 12 \, {\rm yr} \, (E_{\rm j}/10^{49} \, {\rm erg})^{1/3} (n/1 \, {\rm cm}^{-3})^{-1/3} (\upsilon _{\rm ej}/0.1c)^{-5/3}$, where Ej is the jet energy, n the ISM density, and $\upsilon$ej the ejecta velocity. We also show that this can produce a sharply peaked enhancement in the light curve at t = tcol. 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 multidimensional nature of the problem.
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 and long gamma-ray bursts. 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 or AIC FRB channel will be characterized by distinct host galaxy and spatial offset distributions which we show are consistent with the recently reported FRB 180924, localized by the Australian Square Kilometre Array Pathfinder 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. 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.
Abstract
Numerical models of collisionless shocks robustly predict an electron distribution composed of both thermal and nonthermal electrons. Here, we explore in detail the effect of thermal ...electrons on the emergent synchrotron emission from subrelativistic shocks. We present a complete “thermal + nonthermal” synchrotron model and derive properties of the resulting spectrum and light curves. Using these results, we delineate the relative importance of thermal and nonthermal electrons for subrelativistic shock-powered synchrotron transients. We find that thermal electrons are naturally expected to contribute significantly to the peak emission if the shock velocity is ≳0.2
c
, but would be mostly undetectable in nonrelativistic shocks. This helps explain the dichotomy between typical radio supernovae and the emerging class of “AT2018cow-like” events. The signpost of thermal electron synchrotron emission is a steep optically-thin spectral index and a
ν
2
optically-thick spectrum. These spectral features are also predicted to correlate with a steep postpeak light-curve decline rate, broadly consistent with observed AT2018cow-like events. We expect that thermal electrons may be observable in other contexts where mildly relativistic shocks are present and briefly estimate this effect for gamma-ray burst afterglows and binary–neutron-star mergers. Our model can be used to fit spectra and light curves of events and accounts for both thermal and nonthermal electron populations with no additional physical degrees of freedom.
ABSTRACT 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 multimessenger astronomy. Between the plethora of modelling 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, and GRB170817 to obtain truly multimessenger constraints on the tidal deformability $\tilde{\Lambda } \in 302,860$, total binary mass M ∈ 2.722, 2.751 M⊙, 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 ≤ 1.27, all at 90 per cent confidence. Our joint novel analysis uses new phenomenological descriptions of the dynamical ejecta, debris disc mass, and remnant black hole properties, all derived from a large suite of numerical relativity simulations.