Expulsion of neutron-rich matter following the merger of neutron star binaries is crucial to the radioactively powered electromagnetic counterparts of these events and to their relevance as sources ...of r-process nucleosynthesis. Here we explore the long-term (viscous) evolution of remnant black hole accretion discs formed in such mergers by means of two-dimensional, time-dependent hydrodynamical simulations. The evolution of the electron fraction due to charged-current weak interactions is included, and neutrino self-irradiation is modelled as a lightbulb that accounts for the disc geometry and moderate optical depth effects. Over several viscous times (∼1 s), a fraction of ∼10 per cent of the initial disc mass is ejected as a moderately neutron-rich wind (Y
e ∼ 0.2) powered by viscous heating and nuclear recombination, with neutrino self-irradiation playing a sub-dominant role. Although the properties of the outflow vary in time and direction, their mean values in the heavy-element production region are relatively robust to variations in the initial conditions of the disc and the magnitude of its viscosity. The outflow is sufficiently neutron-rich that most of the ejecta forms heavy r-process elements with mass number A 130, thus representing a new astrophysical source of r-process nucleosynthesis, distinct from that produced in the dynamical ejecta. Due to its moderately high entropy, disc outflows contain a small residual fraction ∼1 per cent of helium, which could produce a unique spectroscopic signature.
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.
The mergers of binaries containing neutron stars and stellar-mass black holes are among the most promising sources for direct detection in gravitational waves by the interferometers Advanced LIGO and ...Virgo over the next few years. The concurrent detection of electromagnetic emission from these events would greatly enhance the scientific return of these discoveries. We review the state of the art in modeling the electromagnetic signal of neutron star binary mergers across different phases of the merger and multiple wavelengths. We focus on those observables that provide the most sensitive diagnostics of the merger physics and the contribution to the synthesis of rapid neutron capture (
r
-process) elements in the Galaxy. We also outline expected future developments on the observational and theoretical sides of this rapidly evolving field.
ABSTRACT
The recurrent fast radio burst FRB 180916 was recently shown to exhibit a 16-d period (with possible aliasing) in its bursting activity. Given magnetars as widely considered FRB sources, ...this period has been attributed to precession of the magnetar spin axis or the orbit of a binary companion. Here, we make the simpler connection to a rotational period, an idea observationally motivated by the 6.7-h period of the Galactic magnetar candidate, 1E 161348–5055. We explore three physical mechanisms that could lead to the creation of ultralong period magnetars: (i) enhanced spin-down due to episodic mass-loaded charged particle winds (e.g. as may accompany giant flares), (ii) angular momentum kicks from giant flares, and (iii) fallback leading to long-lasting accretion discs. We show that particle winds and fallback accretion can potentially lead to a sub-set of the magnetar population with ultralong periods, sufficiently long to accommodate FRB 180916 or 1E 161348–5055. If confirmed, such periods implicate magnetars in relatively mature states (ages 1−10 kyr) and which possessed large internal magnetic fields at birth Bint ≳ 1016 G. In the low-twist magnetar model for FRBs, such long period magnetars may dominate FRB production for repeaters at lower isotropic-equivalent energies and broaden the energy distribution beyond that expected for a canonical population of magnetars, which terminate their magnetic activity at shorter periods P ≲ 10 s.
Rates of stellar tidal disruption events (TDEs) by supermassive black holes (SMBHs) due to two-body relaxation are calculated using a large galaxy sample (N ≈ 200) in order to explore the sensitivity ...of the TDE rates to observational uncertainties, such as the parametrization of galaxy light profiles and the stellar mass function. The largest uncertainty arises due to the poorly constrained occupation fraction of SMBHs in low-mass galaxies, which otherwise dominate the total TDE rate. The detection rate of TDE flares by optical surveys is calculated as a function of SMBH mass and other observables for several physically motivated models of TDE emission. We also quantify the fraction of galaxies that produce deeply penetrating disruption events. If the majority of the detected events are characterized by super-Eddington luminosities (such as disc winds, or synchrotron radiation from an off-axis relativistic jet), then the measured SMBH mass distribution will tightly constrain the low-end SMBH occupation fraction. If Eddington-limited emission channels dominate, however, then the occupation fraction sensitivity is much less pronounced in a flux-limited survey (although still present in a volume-complete event sample). The SMBH mass distribution of the current sample of TDEs, though highly inhomogeneous and encumbered by selection effects, already suggests that Eddington-limited emission channels dominate. Even our most conservative rate estimates appear to be in tension with much lower observationally inferred TDE rates, and we discuss several possible resolutions to this discrepancy.
Abstract
Luminous red novae (LRN) are a class of optical transients believed to originate from the mergers of binary stars, or ‘common envelope’ events. Their light curves often show secondary ...maxima, which cannot be explained in the previous models of thermal energy diffusion or hydrogen recombination without invoking multiple independent shell ejections. We propose that double-peaked light curves are a natural consequence of a collision between dynamically ejected fast shell and pre-existing equatorially focused material, which was shed from the binary over many orbits preceding the dynamical event. The fast shell expands freely in the polar directions, powering the initial optical peak through cooling envelope emission. Radiative shocks from the collision in the equatorial plane power the secondary light-curve peak on the radiative diffusion time-scale of the deeper layers, similar to luminous Type IIn supernovae and some classical novae. Using a detailed 1D analytic model, informed by complementary 3D hydrodynamical simulations, we show that shock-powered emission can explain the observed range of peak time-scales and luminosities of the secondary peaks in LRN for realistic variations in the binary parameters and fraction of the binary mass ejected. The dense shell created by the radiative shocks in the equatorial plane provides an ideal location for dust nucleation consistent with the inferred aspherical geometry of dust in LRN. For giant stars, the ejecta forms dust when the shock-powered luminosity is still high, which could explain the infrared transients recently discovered by Spitzer. Our results suggest that pre-dynamical mass-loss is common if not ubiquitous in stellar mergers, providing insight into the instabilities responsible for driving the binary merger.
The merger of a neutron star (NS) binary may result in the formation of a long-lived, or indefinitely stable, millisecond magnetar remnant surrounded by a low-mass ejecta shell. A portion of the ...magnetar's prodigious rotational energy is deposited behind the ejecta in a pulsar wind nebula, powering luminous optical/X-ray emission for hours to days following the merger. Ions in the pulsar wind may also be accelerated to ultra-high energies, providing a coincident source of high-energy cosmic rays and neutrinos. At early times, the cosmic rays experience strong synchrotron losses; however, after a day or so, pion production through photomeson interaction with thermal photons in the nebula comes to dominate, leading to efficient production of high-energy neutrinos. After roughly a week, the density of background photons decreases sufficiently for cosmic rays to escape the source without secondary production. These competing effects result in a neutrino light curve that peaks on a few day timescale near an energy of ∼1018eV. This signal may be detectable for individual mergers out to ∼10 (100) Mpc by current (next generation) neutrino telescopes, providing clear evidence for a long-lived NS remnant, the presence of which may otherwise be challenging to identify from the gravitational waves alone. Under the optimistic assumption that a sizable fraction of NS mergers produce long-lived magnetars, the cumulative cosmological neutrino background is estimated to be for an NS merger rate of , overlapping with IceCube's current sensitivity and within the reach of next-generation neutrino telescopes.
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.
Abstract
We present a toy model for the thermal optical/UV/X-ray emission from tidal disruption events (TDEs). Motivated by recent hydrodynamical simulations, we assume that the debris streams ...promptly and rapidly circularize (on the orbital period of the most tightly bound debris), generating a hot quasi-spherical pressure-supported envelope of radius
R
v
∼ 10
14
cm (photosphere radius ∼10
15
cm) surrounding the supermassive black hole (SMBH). As the envelope cools radiatively, it undergoes Kelvin–Helmholtz contraction
R
v
∝
t
−1
, its temperature rising
T
eff
∝
t
1/2
while its total luminosity remains roughly constant; the optical luminosity decays as
ν
L
ν
∝
R
v
2
T
eff
∝
t
−
3
/
2
. Despite this similarity to the mass fallback rate
M
̇
fb
∝
t
−
5
/
3
, envelope heating from fallback accretion is subdominant compared to the envelope cooling luminosity except near optical peak (where they are comparable). Envelope contraction can be delayed by energy injection from accretion from the inner envelope onto the SMBH in a regulated manner, leading to a late-time flattening of the optical/X-ray light curves, similar to those observed in some TDEs. Eventually, as the envelope contracts to near the circularization radius, the SMBH accretion rate rises to its maximum, in tandem with the decreasing optical luminosity. This cooling-induced (rather than circularization-induced) delay of up to several hundred days may account for the delayed onset of thermal X-rays, late-time radio flares, and high-energy neutrino generation, observed in some TDEs. We compare the model predictions to recent TDE light-curve correlation studies, finding both agreement and points of tension.
The neutron star (NS) merger GW170817 was followed over several days by optical-wavelength ("blue") kilonova (KN) emission likely powered by the radioactive decay of light r-process nuclei ...synthesized by ejecta with a low neutron abundance (electron fraction Ye 0.25-0.35). While the composition and high velocities of the blue KN ejecta are consistent with shock-heated dynamical material, the large quantity is in tension with the results of numerical simulations. We propose an alternative ejecta source: the neutrino-heated, magnetically accelerated wind from the strongly magnetized hypermassive NS (HMNS) remnant. A rapidly spinning HMNS with an ordered surface magnetic field of strength B (1-3) × 1014 G and lifetime trem ∼ 0.1-1 s can simultaneously explain the velocity, total mass, and electron fraction of the blue KN ejecta. The inferred HMNS lifetime is close to its Alfvén crossing time, suggesting that global magnetic torques could be responsible for bringing the HMNS into solid-body rotation and instigating its gravitational collapse. Different origins for the KN ejecta may be distinguished by their predictions for the emission in the first hours after the merger, when the luminosity is enhanced by heating from internal shocks; the latter are likely generic to any temporally extended ejecta source (e.g., magnetar or accretion disk wind) and are not unique to the emergence of a relativistic jet. The same shocks could mix and homogenize the composition to a low but nonzero lanthanide mass fraction, , as advocated by some authors, but only if the mixing occurs after neutrons are consumed in the r-process on a timescale 1 s.