The LIGO/Virgo Scientific Collaboration recently announced the detection of a compact object binary merger, GW190412, as the first asymmetric binary black hole (BBH) merger with mass ratio $q\approx0.25$. Other than the mass ratio, this BBH has shown to have a positive effective spin of around $\chi_{\rm eff}\approx0.28$. Assuming a field formation channel, associating this effective spin to either the primary or the secondary BH each has its implications: If the spin of the BBH comes form the primary BH, it has consequences for the efficiency of angular momentum transport in the formation of the BH. If, on the other hand, the spin is due to the secondary BH through tidal spin-up processes, one has to note that (i) such processes have very short delay-times, and (ii) subsequently, their local merger rate is determined by local star formation rate at assumed formation metallicity of the BBH. We show that the predicted merger rate density from this channel is $\lesssim 0.3~\rm Gpc^{-3} yr^{-1}$ and in tension with the rather high local merger rate of such systems which we estimate from this single event to be $\sim 1.7^{+2.5}_{-1.4}~\rm Gpc^{-3} yr^{-1}$ (90\% confidence interval, and assuming 50 days of observing time). Large natal kicks ($v\gtrsim 500\,{\rm km/s}$) would be required to get such BBHs with an in-plane spin component to account for the marginal detection of precession in GW190412. However, this would only exacerbate the tension as the estimated local merger rate would be further decreased. Similarly, the formation of such systems through the dynamical assembly is exceedingly rare, leaving this system a dilemma hard to account for with the currently accepted paradigms of BBH formation.

Joint multi-messenger observations with gravitational waves and electromagnetic data offer new insights into the astrophysical studies of compact objects. The third Advanced LIGO and Advanced Virgo observing run began on April 1, 2019; during the eleven months of observation, there have been 14 compact binary systems candidates for which at least one component is potentially a neutron star. Although intensive follow-up campaigns involving tens of ground and space-based observatories searched for counterparts, no electromagnetic counterpart has been detected. Following on a previous study of the first six months of the campaign, we present in this paper the next five months of the campaign from October 2019 to March 2020. We highlight two neutron star - black hole candidates (S191205ah, S200105ae), two binary neutron star candidates (S191213g and S200213t) and a binary merger with a possible neutron star and a "MassGap" component, S200115j. Assuming that the gravitational-wave candidates are of astrophysical origin and their location was covered by optical telescopes, we derive possible constraints on the matter ejected during the events based on the non-detection of counterparts. We find that the follow-up observations during the second half of the third observing run did not meet the necessary sensitivity to constrain the source properties of the potential gravitational-wave candidate. Consequently, we suggest that different strategies have to be used to allow a better usage of the available telescope time. We examine different choices for follow-up surveys to optimize sky localization coverage vs.\ observational depth to understand the likelihood of counterpart detection.

Theory and observations suggest that single-star evolution is not able to produce black holes (BHs) with masses in the range $3-5M_{\odot}$ and above $\sim 45M_{\odot}$, referred to as the lower mass gap (LMG) and the upper mas gap (UMG), respectively. However, it is possible to form BHs in these gaps through merger of compact objects in dense clusters, e.g. the LMG and the UMG can be populated through binary neutron star- and BBH mergers, respectively. This implies that if binary mergers are observed in gravitational waves (GWs) with at least one mass gap object, then either clusters are effective in assembling binary mergers, or our single-star models have to be revised. Understanding how effective clusters are at populating both mass gaps have therefore major implications for both stellar- and GW astrophysics. In this paper we present a systematic study on how efficient stellar clusters are at populating both mass gaps through in-cluster GW mergers. For this, we derive a set of closed form relations for describing the evolution of compact object binaries undergoing dynamical interactions and GW merger inside their cluster. By considering both static and time evolving populations, we find in particular that globular clusters are clearly inefficient at populating the LMG in contrast to the UMG. We further describe how these results relate to the characteristic mass, time, and length scales associated with the problem.

We study turbulent gas diffusion affects on $r$-process abundances in Milky Way stars, by a combination of an analytical approach and a Monte Carlo simulation. Higher $r$-process event rates and faster diffusion, lead to more efficient mixing corresponding to a reduced scatter of $r$-process abundances and causing $r$-process enriched stars to start appearing at lower metallicities. We use three independent observations to constrain the model parameters: (i) the scatter of radioactively stable $r$-process element abundances, (ii) the largest $r$-process enrichment values observed in any solar neighborhood stars and (iii) the isotope abundance ratios of different radioactive $r$-process elements ($^{244}$Pu/$^{238}$U and $^{247}$Cm/$^{238}$U) at the early solar system as compared to their formation. Our results indicate that the Galactic $r$-process rate and the diffusion coefficient are respectively $r<4\times 10^{-5}\mbox{ yr}^{-1}, D>0.1 \mbox{ kpc}^2\mbox{Gyr}^{-1}$ ($r<4\times 10^{-6}\mbox{ yr}^{-1}, D>0.5 \mbox{ kpc}^2\mbox{Gyr}^{-1}$ for collapsars or similarly prolific $r$-process sources) with allowed values satisfying an approximate anti-correlation such that $D\approx r^{-2/3}$, implying that the time between two $r$-process events that enrich the same location in the Galaxy, is $\tau_{\rm mix}\approx 100-200\mbox{ Myr}$. This suggests that a fraction of $\sim 0.8$ ($\sim 0.5$) of the observed $^{247}$Cm ($^{244}$Pu) abundance is dominated by one $r$-process event in the early solar system. Radioactively stable element abundances are dominated by contributions from $\sim 10$ different events in the early solar system. For metal poor stars (with [Fe/H]$\lesssim -2$), their $r$-process abundances are dominated by either a single or several events, depending on the star formation history.

We present the full panchromatic afterglow light curve data of GW170817, including new radio data as well as archival optical and X-ray data, between 0.5 and 940 days post-merger. By compiling all archival data, and reprocessing a subset of it, we have ensured that the panchromatic dataset is uniform and therefore immune to the differences in data processing or flux determination methods used by different groups. Simple power-law fits to the uniform afterglow light curve indicate a $t^{0.86\pm0.04}$ rise, a $t^{-1.90\pm0.12}$ decline, and a peak occurring at $155\pm4$ days. The afterglow is optically thin throughout its evolution, consistent with a single spectral index ($-0.569\pm0.002$) across all epochs. This gives a precise and updated estimate of the electron power-law index, $p=2.138\pm0.004$. By studying the diffuse X-ray emission from the host galaxy, we place a conservative upper limit on the hot ionized ISM density, $<$0.01 cm$^{-3}$, consistent with previous afterglow studies. Using the late-time afterglow data we rule out any long-lived neutron star remnant having magnetic field strength between 10$^{10.4}$ G and 10$^{16}$ G. Our fits to the afterglow data using an analytical model that includes VLBI proper motion from Mooley et al (2018), and a structured jet model that ignores the proper motion, indicates that the proper motion measurement needs to be considered while seeking an accurate estimate of the viewing angle.

We present low-radio-frequency follow-up observations of AT 2017gfo, the electromagnetic counterpart of GW170817, which was the first binary neutron star merger to be detected by Advanced LIGO-Virgo. These data, with a central frequency of 144 MHz, were obtained with LOFAR, the Low-Frequency Array. The maximum elevation of the target is just 13.7 degrees when observed with LOFAR, making our observations particularly challenging to calibrate and significantly limiting the achievable sensitivity. On time-scales of 130-138 and 371-374 days after the merger event, we obtain 3$\sigma$ upper limits for the afterglow component of 6.6 and 19.5 mJy beam$^{-1}$, respectively. Using our best upper limit and previously published, contemporaneous higher-frequency radio data, we place a limit on any potential steepening of the radio spectrum between 610 and 144 MHz: the two-point spectral index $\alpha^{610}_{144} \gtrsim -2.5$. We also show that LOFAR can detect the afterglows of future binary neutron star merger events occurring at more favourable elevations.

The jet opening angle and inclination of GW170817 -- the first detected binary neutron star merger -- were vital to understand its energetics, relation to short gamma-ray bursts, and refinement of the standard siren-based determination of the Hubble constant, $H_0$. These basic quantities were determined through a combination of the radio lightcurve and Very Long Baseline Interferometry (VLBI) measurements of proper motion. In this paper we discuss and quantify the prospects for the use of radio VLBI observations and observations of scintillation-induced variability to measure the source size and proper motion of merger afterglows, and thereby infer properties of the merger including inclination angle, opening angle and energetics. We show that these techniques are complementary as they probe different parts of the circum-merger density/inclination angle parameter space and different periods of the temporal evolution of the afterglow. We also find that while VLBI observations will be limited to the very closest events it will be possible to detect scintillation for a large fraction of events beyond the range of current gravitational wave detectors. Scintillation will also be detectable with next generation telescopes such as the Square Kilometre Array, 2000 antenna Deep Synoptic Array and the next generation Very Large Array, for a large fraction of events detected with third generation gravitational wave detectors. Finally, we discuss prospects for the measurement of the $H_0$ with VLBI observations of neutron star mergers and compare this technique to other standard siren methods.

The discovery of GW170817 with gravitational waves (GWs) and electromagnetic (EM) radiation is prompting new questions in strong-gravity astrophysics. Importantly, it remains unknown whether the progenitor of the merger comprised two neutron stars (NSs), or a NS and a black hole (BH). Using new numerical-relativity simulations and incorporating modeling uncertainties we produce novel GW and EM observables for NS-BH mergers with similar masses. A joint analysis of GW and EM measurements reveals that if GW170817 is a NS-BH merger, <40% of the binary parameters consistent with the GW data are compatible with EM observations.

On 2019 August 14, the Advanced LIGO and Virgo interferometers detected the high-significance gravitational wave (GW) signal S190814bv. The GW data indicated that the event resulted from a neutron star--black hole (NSBH) merger, or potentially a low-mass binary black hole merger. Due to the low false alarm rate and the precise localization (23 deg$^2$ at 90\%), S190814bv presented the community with the best opportunity yet to directly observe an optical/near-infrared counterpart to a NSBH merger. To search for potential counterparts, the GROWTH collaboration performed real-time image subtraction on 6 nights of public Dark Energy Camera (DECam) images acquired in the three weeks following the merger, covering $>$98\% of the localization probability. Using a worldwide network of follow-up facilities, we systematically undertook spectroscopy and imaging of optical counterpart candidates. Combining these data with a photometric redshift catalog, we ruled out each candidate as the counterpart to S190814bv and we placed deep, uniform limits on the optical emission associated with S190814bv. For the nearest consistent GW distance, radiative transfer simulations of NSBH mergers constrain the ejecta mass of S190814bv to be $M_\mathrm{ej} < 0.04$~$M_{\odot}$ at polar viewing angles, or $M_\mathrm{ej} < 0.03$~$M_{\odot}$ if the opacity is $\kappa < 2$~cm$^2$g$^{-1}$. Assuming a tidal deformability for the neutron star at the high end of the range compatible with GW170817 results, our limits would constrain the BH spin component aligned with the orbital momentum to be $ \chi < 0.7$ for mass ratios $Q < 6$, with weaker constraints for more compact neutron stars. We publicly release the photometry from this campaign at http://www.astro.caltech.edu/~danny/static/s190814bv.

GW170817 showed that neutron star mergers not only emit gravitational waves but also can release electromagnetic signatures in multiple wavelengths. Within the first half of the third observing run of the Advanced LIGO and Virgo detectors, there have been a number of gravitational wave candidates of compact binary systems for which at least one component is potentially a neutron star. In this article, we look at the candidates S190425z, S190426c, S190510g, S190901ap, and S190910h, predicted to have potentially a non-zero remnant mass, in more detail. All these triggers have been followed up with extensive campaigns by the astronomical community doing electromagnetic searches for their optical counterparts; however, according to the released classification, there is a high probability that some of these events might not be of extraterrestrial origin. Assuming that the triggers are caused by a compact binary coalescence and that the individual source locations have been covered during the EM follow-up campaigns, we employ three different kilonova models and apply them to derive possible constraints on the matter ejection consistent with the publicly available gravitational-wave trigger information and the lack of a kilonova detection. These upper bounds on the ejecta mass can be related to limits on the maximum mass of the binary neutron star candidate S190425z and to constraints on the mass-ratio, spin, and NS compactness for the potential black hole-neutron star candidate S190426c. Our results show that deeper electromagnetic observations for future gravitational wave events near the horizon limit of the advanced detectors are essential.

We present results from a search for a radio transient associated with the LIGO/Virgo source S190814bv, a likely neutron star-black hole (NSBH) merger, with the Australian Square Kilometre Array Pathfinder. We imaged a $30\,{\rm deg}^2$ field at $\Delta T$=2, 9 and 33 days post-merger at a frequency of 944\,MHz, comparing them to reference images from the Rapid ASKAP Continuum Survey observed 110 days prior to the event. Each epoch of our observations covers $89\%$ of the LIGO/Virgo localisation region. We conducted an untargeted search for radio transients in this field, resulting in 21 candidates. For one of these, \object[AT2019osy]{AT2019osy}, we performed multi-wavelength follow-up and ultimately ruled out the association with S190814bv. All other candidates are likely unrelated variables, but we cannot conclusively rule them out. We discuss our results in the context of model predictions for radio emission from neutron star-black hole mergers and place constrains on the circum-merger density and inclination angle of the merger. This survey is simultaneously the first large-scale radio follow-up of an NSBH merger, and the most sensitive widefield radio transients search to-date.

We discuss prospects of identifying and characterizing black hole (BH) companions to normal stars on tight but detached orbits, using photometric data from the Transiting Exoplanet Survey Satellite (TESS). We focus on the following two periodic signals from the visible stellar component: (i) in-eclipse brightening of the star due to gravitational microlensing by the BH (self-lensing), and (ii) a combination of ellipsoidal variations due to tidal distortion of the star and relativistic beaming due to its orbital motion (phase-curve variation). We evaluate the detectability of each signal in the light curves of stars in the TESS input catalog, based on a pre-launch noise model of TESS photometry as well as the actual light curves of spotted stars from the prime Kepler mission to gauge the potential impact of stellar activity arising from the tidally spun-up stellar components. We estimate that the self-lensing and phase-curve signals from BH companions, if exist, will be detectable in the light curves of effectively $\mathcal{O}(10^5)$ and $\mathcal{O}(10^6)$ low-mass stars, respectively, taking into account orbital inclination dependence of the signals. These numbers could be large enough to actually detect signals from BHs: simple population models predict some 10 and 100 detectable BHs among these "searchable" stars, although the latter may be associated with a comparable number of false-positives due to stellar variabilities and additional vetting with radial velocity measurements would be essential. Thus the TESS data could serve as a resource to study nearby BHs with stellar companions on shorter-period orbits than will potentially be probed with Gaia.

Gravitational wave (GW) sources are an excellent probe of the luminosity distance and offer a novel measure of the Hubble constant, $H_0$. This estimation of $H_0$ from standard sirens requires an accurate estimation of the cosmological redshift of the host galaxy of the GW source, after correcting for its peculiar velocity. Absence of an accurate peculiar velocity correction affects both the precision and accuracy of the measurement of $H_0$, particularly for nearby sources. We propose a framework to incorporate such a peculiar velocity correction for GW sources. The implementation of our method to the event GW170817 combined with the Very Large Baseline Interferometry (VLBI) observation leads to a revised value of $H_0= 69.3^{+ 4.5}_{-4.0}$ km s$^{-1}$Mpc$^{-1}$; more importantly, it demonstrates that our method will allow to achieve in practice unbiased, accurate measurements of $H_0$ at the precision promised by standard siren forecasts.

We study the heating rate of r-process nuclei and thermalization of decay products in neutron star merger ejecta and macronova (kilonova) light curves. Thermalization of charged decay products, i.e., electrons, $\alpha$-particles, and fission fragments is calculated according to their injection energy. The $\gamma$-ray thermalization processes are also properly calculated by taking the $\gamma$-ray spectrum of each decay into account. We show that the $\beta$-decay heating rate at later times approaches a power-law decline as $\propto t^{-2.8}$, which agrees with the result of Waxman et al. (2019). We present a new analytic model to calculate macronova light curves, in which the density structure of the ejecta is accounted for. We demonstrate that the observed bolometric light curve and temperature evolution of the macronova associated with GW170817 are reproduced well by the $\beta$-decay heating rate with the solar r-process abundance pattern. We interpret the break in the observed bolometric light curve around a week as a result of the diffusion wave crossing a significant part of the ejecta rather than a thermalization break. We also show that the time-weighted integral of the bolometric light curve (Katz integral) is useful to provide an estimate of the total r-process mass from the observed data, which is independent of the highly uncertain radiative transfer. For the macronova in GW170817, the ejecta mass is robustly estimated as $\approx 0.05M_{\odot}$ for $A_{\rm min}\leq 72$ and $85\leq A_{\rm min}\leq 130$ with the solar r-process abundance pattern. The code for computation of the heating rate and light curve for given initial nuclear abundances is publicly available.

Mergers of binary neutron stars and black hole-neutron star binaries are one of the most promising sources for the ground-based gravitational-wave (GW) detectors and also a high-energy astrophysical phenomenon as illustrated by the observations of gravitational waves and electromagnetic (EM) waves in the event of GW170817. Mergers of these neutron-star binaries are also the most promising site for r-process nucleosynthesis. Numerical simulation in full general relativity (numerical relativity) is a unique approach to the theoretical prediction of the merger process, GWs emitted, mass ejection process, and resulting EM emission. We summarize our current understanding for the processes of neutron star mergers and subsequent mass ejection based on the results of the latest numerical-relativity simulations. We emphasize that the predictions of the numerical-relativity simulations agrees broadly with the optical and infrared observations of GW170817.

We present a systematic numerical relativity study of the mass ejection and the associated electromagnetic transients and nucleosynthesis from binary neutron star (NS) mergers. We find that a few $10^{-3}\, M_\odot$ of material are ejected dynamically during the mergers. The amount and the properties of these outflow depend on binary parameters and on the NS equation of state (EOS). A small fraction of these ejecta, typically ${\sim}10^{-6}\, M_\odot$, is accelerated by shocks formed shortly after merger to velocities larger than $0.6\, {\rm c}$ and produces bright radio flares on timescales of weeks, months, or years after merger. Their observation could constrain the strength with which the NSs bounce after merger and, consequently, the EOS of matter at extreme densities. The dynamical ejecta robustly produce second and third $r$-process peak nuclei with relative isotopic abundances close to solar. The production of light $r$-process elements is instead sensitive to the binary mass ratio and the neutrino radiation treatment. Accretion disks of up to ${\sim}0.2\, M_\odot$ are formed after merger, depending on the lifetime of the remnant. In most cases, neutrino- and viscously-driven winds from these disks dominate the overall outflow. Finally, we generate synthetic kilonova light curves and find that kilonovae depend on the merger outcome and could be used to constrain the NS EOS.

We constrain the formation rate of Galactic magnetars directly from observations. Combining spin-down rates, magnetic activity, and association with supernova remnants, we put a 2$\sigma$ limit on their Galactic formation rate at $2.3-20\mbox{kyr}^{-1}$. This leads to a fraction $0.4_{-0.28}^{+0.6}$ of neutron stars being born as magnetars. We study evolutionary channels that can account for this rate as well as for the periods, period derivatives and luminosities of the observed population. We find that their typical magnetic fields at birth are $3\times 10^{14}-10^{15}$G, and that those decay on a time-scale of $\sim 10^4$years, implying a maximal magnetar period of $P_{\rm max}\approx 13$s. A sizable fraction of the magnetars' energy is released in outbursts. Giant Flares with $E\geq 10^{46}$ erg are expected to occur in the Galaxy at a rate of $\sim 5\mbox{kyr}^{-1}$. Outside our Galaxy, such flares remain observable by {\it Swift} up to a distance of $\sim 100$~Mpc, implying a detection rate of $\sim 5\mbox{ yr}^{-1}$. The specific form of magnetic energy decay is shown to be strongly tied to the total number of observable magnetars in the Galaxy. A systematic survey searching for magnetars could determine the former and inform physical models of magnetic field decay.

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 report our Spitzer Space Telescope observations and detections of the binary neutron star merger GW170817. At 4.5um, GW170817 is detected at 21.9 mag AB at +43 days and 23.9 mag AB at +74 days after merger. At 3.6um, GW170817 is not detected to a limit of 23.2 mag AB at +43 days and 23.1 mag AB at +74 days. Our detections constitute the latest and reddest constraints on the kilonova/macronova emission and composition of heavy elements. The 4.5um luminosity at this late phase cannot be explained by elements exclusively from the first abundance peak of the r-process. Moreover, the steep decline in the Spitzer band, with a power-law index of 3.4 +/- 0.2, can be explained by a few of the heaviest isotopes in the third abundance peak with half-life around 14 days dominating the luminosity (e.g. 140Ba, 143Pr, 147Nd, 156Eu, 191Os, 223Ra, 225Ra, 233Pa, 234Th) or a model with lower deposition efficiency. This data offers evidence that the heaviest elements in the second and third r-process abundance peak were indeed synthesized. Our conclusion is verified by both analytics and network simulations and robust despite intricacies and uncertainties in the nuclear physics. Future observations with Spitzer and James Webb Space Telescope will further illuminate the relative abundance of the synthesized heavy elements.

General-relativistic simulations of binary neutron star mergers with viscosity reveal a new outflow mechanism operating in unequal mass binaries on dynamical timescales and enabled by turbulent viscosity. These "viscous-dynamical" ejecta are launched during merger due to the thermalization of mass exchange streams between the secondary and the primary neutron star. They are characterized by asymptotic velocities extending up to ${\sim} 0.8\, c$, and have masses that depend on the efficiency of the viscous mechanism. Depending on the unknown strength of the effective viscosity arising from magnetohydrodynamics instabilities operating during merger, the overall mass of the dynamical ejecta could be enhanced by a factor of a few and the mass of the fast tail of the ejecta having asymptotic velocities $\geq 0.6\, c$ by up to four orders of magnitude. The radioactive decay of the expanding viscous-dynamical ejecta could produce bright kilonova transients with signatures of free neutron decay in the first hour and enhanced near infrared flux on a timescale of a few days. The synchrotron remnant produced by the interaction between the ejecta and the interstellar medium could also be significantly enhanced by viscosity. Such remnant could be detected in the case of GW170817 as a rebrightening of the radio signal in the next months to years.

The binary neutron star merger GW170817 was accompanied by radiation across the electromagnetic spectrum and localized to the galaxy NGC 4993 at a distance of 41+/-3 Mpc. The radio and X-ray afterglows of GW170817 exhibited delayed onset, a gradual rise in the emission with time as t^0.8, a peak at about 150 days post-merger, followed by a relatively rapid decline. To date, various models have been proposed to explain the afterglow emission, including a choked-jet cocoon and a successful-jet cocoon (a.k.a. structured jet). However, the observational data have remained inconclusive as to whether GW170817 launched a successful relativistic jet. Here we show, through Very Long Baseline Interferometry, that the compact radio source associated with GW170817 exhibits superluminal motion between two epochs at 75 and 230 days post-merger. This measurement breaks the degeneracy between the models and indicates that, while the early-time radio emission was powered by a wider-angle outflow (cocoon), the late-time emission was most likely dominated by an energetic and narrowly-collimated jet, with an opening angle of <5 degrees, and observed from a viewing angle of about 20 degrees. The imaging of a collimated relativistic outflow emerging from GW170817 adds substantial weight to the growing evidence linking binary neutron star mergers and short gamma-ray bursts.

We develop stellar population models to predict the number of binaries with a single luminous member is Gaia and Hipparcos. Our models yield dozens of detections of black hole - luminous companion binaries (BHLC) and hundreds to thousands of neutron star - luminous companion binaries (NSLC) with Gaia. Interestingly, our models also yield a single detection of BHLC binary with Hipparcos, and a few NSLC binaries. We also show how the statistical distribution of detected binaries with a single luminous companion can be used to constrain the formation process of neutron stars and black holes.

LIGO's detection of gravitational waves from binary black hole mergers was an unexpected surprise that immediately raised the question - what is the origin of these black hole binaries? The "simplest" scenario is evolution of field massive stellar binaries. However, other possibilities involving capture have been proposed. We explore here one of the more interesting clues on this puzzle: the relatively modest spins of the resulting black holes that imply that the progenitor black holes were not spinning rapidly. More specifically we consider the implication of observed distribution of, $\chi_{\rm eff}$, the mass weighted projected (along the orbital axis) spins on the field evolution scenario. In all cases $\chi_{\rm eff}$ is small and in two of the cases the best fit value is negative. Only in one event the spin is positive at 90\% credible. These observations are puzzling within the field binary scenario in which positive higher spins ($\chi_{\rm eff} \ge 0.5$) are expected. At first sight one may expect that this rules out the field evolutionary scenario. Indeed we show that with typical parameters a significant fraction ($\ge 25\%$) of the mergers should have high effective spin values. However, uncertainties in the outcome of the common envelope phase (the typical separation and whether the stars are rotating or not) and in the late stages of massive star evolution (the strength of the winds) make it impossible to rule out, at present, these scenarios. While observations of mergers with high effective spin will support this scenario, future observations of negative spin mergers would rule it out.

The Hubble constant ($H_0$) measures the current expansion rate of the Universe, and plays a fundamental role in cosmology. Tremendous effort has been dedicated over the past decades to measure $H_0$. Notably, Planck cosmic microwave background (CMB) and the local Cepheid-supernovae distance ladder measurements determine $H_0$ with a precision of $\sim 1\%$ and $\sim 2\%$ respectively. A $3$-$\sigma$ level of discrepancy exists between the two measurements, for reasons that have yet to be understood. Gravitational wave (GW) sources accompanied by electromagnetic (EM) counterparts offer a completely independent standard siren (the GW analogue of an astronomical standard candle) measurement of $H_0$, as demonstrated following the discovery of the neutron star merger, GW170817. This measurement does not assume a cosmological model and is independent of a cosmic distance ladder. The first joint analysis of the GW signal from GW170817 and its EM localization led to a measurement of $H_0=74^{+16}_{-8}$ km/s/Mpc (median and symmetric $68\%$ credible interval). In this analysis, the degeneracy in the GW signal between the source distance and the weakly constrained viewing angle dominated the $H_0$ measurement uncertainty. Recently, Mooley et al. (2018) obtained tight constraints on the viewing angle using high angular resolution imaging of the radio counterpart of GW170817. Here we obtain a significantly improved measurement $H_0=68.9^{+4.7}_{-4.6}$ km/s/Mpc by using these new radio observations, combined with the previous GW and EM data. We estimate that 15 more localized GW170817-like events (comparable signal-to-noise ratio, favorable orientation), having radio images and light curve data, will potentially bring resolution to the tension between the Planck and Cepheid-supernova measurements, as compared to 50-100 GW events without such data.

We present late-time radio observations of GW170817, the first binary neutron star merger discovered through gravitational waves by the advanced LIGO and Virgo detectors. Our observations, carried out with the Karl G. Jansky Very Large Array, were optimized to detect polarized radio emission, and thus to constrain the linear polarization fraction of GW170817. At an epoch of ~244 days after the merger, we rule out linearly polarized emission above a fraction of ~12% at a frequency of 2.8 GHz (99% confidence). Within the structured jet scenario (a.k.a. successful jet plus cocoon system) for GW170817, the derived upper-limit on the radio continuum linear polarization fraction strongly constrains the magnetic field configuration in the shocked ejecta. We show that our results for GW170817 are compatible with the low level of linear polarization found in afterglows of cosmological long gamma-ray bursts. Finally, we discuss our findings in the context of future expectations for the study of radio counterparts of binary neutron star mergers identified by ground-based gravitational-wave detectors.

We consider radio emission from a newborn black hole (BH), which is accompanied by a mini-disk with a mass of $\lesssim M_\odot$. Such a disk can be formed from an outer edge of the progenitor's envelope, especially for metal-poor massive stars and/or massive stars in close binaries. The disk accretion rate is typically super-Eddington and an ultrafast outflow with a velocity of $\sim 0.1\mbox{-}0.3\,c$ will be launched into the circumstellar medium. The outflow forms a collisionless shock, and electrons are accelerated and emit synchrotron emission in radio bands with a flux of $\sim 10^{26-30} \ \rm erg \ s^{-1} \ Hz^{-1}$ days to decades after the BH formation. The model predicts not only a fast UV/optical transient but also quasi-simultaneous inverse-Compton X-ray emission $\sim$ a few days after the BH formation, and the discovery of the radio counterpart with coordinated searches will enable us to identify this type of transients. The occurrence rate can be $0.1-10 \ \%$ of the core-collapse supernova rate, which makes them a promising target of dedicated radio observations such as the Jansky VLA Sky Survey.

The short Gamma-Ray Burst, GRB170817A, that followed the binary neutron star merger gravitational waves signal, GW170817, is not a usual sGRB. It is weaker by three orders of magnitude than the weakest sGRB seen before and its spectra, showing a hard early signal followed by a softer thermal spectrum, is unique. We show, first, that the $\gamma$-rays must have emerged from at least mildly relativistic outflow, implying that a relativistic jet was launched following the merger. We then show that the observations are consistent with the predictions of a mildly relativistic shock breakout: a minute $\gamma$-ray energy as compared with the total energy and a rather smooth light curve with a hard to soft evolution. We present here a novel analytic study and detailed numerical 2D and 3D relativistic hydrodynamic and radiation simulations that support the picture in which the observed $\gamma$-rays arose from a shock breakout of a cocoon from the merger's ejecta (Kasliwal 2017). The cocoon can be formed by either a choked jet which does not generate a sGRB (in any direction) or by a successful jet which generates an undetected regular sGRB along the system's axis pointing away from us. Remarkably, for the choked jet model, the macronova signal produced by the ejecta (which is partially boosted to high velocities by the cocoon's shock) and the radio that is produced by the interaction of the shocked cocoon material with the surrounding matter, agree with the observed UV/optical/IR emission and with current radio observations. Finally, we discuss the possibility that the jet propagation within the ejecta may photodissociate some of of the heavy elements and may affect the composition of a fraction of ejecta and, in turn, the opacity and the early macronova light.

We apply a newly-developed On-the-Fly mosaicing technique on the NSF's Karl G. Jansky Very Large Array (VLA) at 3 GHz in order to carry out a sensitive search for an afterglow from the Advanced LIGO binary black hole merger event GW 151226. In three epochs between 1.5 and 6 months post-merger we observed a 100 sq. deg region, with more than 80% of the survey region having a RMS sensitivity of better than 150 uJy/beam, in the northern hemisphere having a merger containment probability of 10%. The data were processed in near-real-time, and analyzed to search for transients and variables. No transients were found but we have demonstrated the ability to conduct blind searches in a time-frequency phase space where the predicted afterglow signals are strongest. If the gravitational wave event is contained within our survey region, the upper limit on any late-time radio afterglow from the merger event at an assumed mean distance of 440 Mpc is about 1e29 erg/s/Hz. Approximately 1.5% of the radio sources in the field showed variability at a level of 30%, and can be attributed to normal activity from active galactic nuclei. The low rate of false positives in the radio sky suggests that wide-field imaging searches at a few Gigahertz can be an efficient and competitive search strategy. We discuss our search method in the context of the recent afterglow detection from GW 170817 and radio follow-up in future gravitational wave observing runs.

We find, using high resolution numerical relativistic simulations, that the tail of the dynamical ejecta of neutron star mergers extends to mildly relativistic velocities faster than $0.7c$. The kinetic energy of this fast tail is $\sim 10^{47}$--$10^{49}$ erg, depending on the neutron star equation of state and on the binary masses. The synchrotron flare arising from the interaction of this fast tail with the surrounding ISM can power the observed non-thermal emission that followed GW170817, provided that the ISM density is $\sim 10^{-2}\,{\rm cm^{-3}}$, the two neutron stars had roughly equal masses and the neutron star equation of state is soft (small neutron star radii). One of the generic predictions of this scenario is that the cooling frequency crosses the X-ray band on a time scale of a few months to a year, leading to a cooling break in the X-ray light curve. If this dynamical ejecta scenario is correct, we expect that the synchrotron radio flare from the ejecta that have produced the macronova/kilonova emission will be observable on time scales of $10^3$ to $10^5$ days. Further multi-frequency observations will confirm or rule out this dynamical ejecta scenario.

Neutron star mergers have been long considered as promising sites of heavy $r$-process nucleosynthesis. We overview observational evidence supporting this scenario including: the total amount of $r$-process elements in the Galaxy, extreme metal poor stars, geological radioactive elemental abundances, dwarf galaxies, and short gamma-ray bursts (sGRBs). Recently, the advanced LIGO and Virgo observatories discovered a gravitational-wave signal of a neutron star merger, GW170817, as well as accompanying multi-wavelength electromagnetic (EM) counterparts. The ultra-violet, optical, and near infrared observations point to $r$-process elements that have been synthesized in the merger ejecta. The rate and ejected mass inferred from GW170817 and the EM counterparts are consistent with other observations. We find however that, within simple one zone chemical evolution models (based on merger rates with reasonable delay time distributions as expected from evolutionary models, or from observations of sGRBs), it is difficult to reconcile the current observations of the europium abundance history of Galactic stars for [Fe/H] $\gtrsim -1$. This implies that to account for the role of mergers in the Galactic chemical evolution, we need a Galactic model with multiple populations that have different spatial distributions and/or varying formation rates.

The recent detection of gravitational waves and electromagnetic counterparts from the double neutron star merger event GW+EM170817, supports the standard paradigm of short gamma-ray bursts (SGRBs) and kilonovae/macronovae. It is important to reveal the nature of the compact remnant left after the merger, either a black hole or neutron star, and their physical link to the origin of the long-lasting emission observed in SGRBs. The diversity of the merger remnants may also lead to different kinds of transients that can be detected in future. Here we study the high-energy emission from the long-lasting central engine left after the coalescence, under certain assumptions. In particular, we consider the X-ray emission from a remnant disk and the non-thermal nebular emission from disk-driven outflows or pulsar winds. We demonstrate that late-time X-ray and high-frequency radio emission can provide useful constraints on properties of the hidden compact remnants and their connections to long-lasting SGRB emission, and we discuss the detectability of nearby merger events through late-time observations at $\sim30-100$ d after the coalescence. We also investigate the GeV-TeV gamma-ray emission that occurs in the presence of long-lasting central engines, and show the importance of external inverse-Compton radiation due to up-scattering of X-ray photons by relativistic electrons in the jet. We also search for high-energy gamma-rays from GW170817 in the Fermi-LAT data, and report upper limits on such long-lasting emission. Finally, we consider the implications of GW+EM170817 and discuss the constraints placed by X-ray and high-frequency radio observations.

GW170817 is the first gravitational wave detection of a binary neutron star merger. It was accompanied by radiation across the electromagnetic spectrum and localized to the galaxy NGC 4993 at a distance of 40 Mpc. It has been proposed that the observed gamma-ray, X-ray and radio emission is due to an ultra-relativistic jet launched during the merger, directed away from our line of sight. The presence of such a jet is predicted from models positing neutron star mergers as the central engines driving short-hard gamma-ray bursts (SGRBs). Here we show that the radio light curve of GW170817 has no direct signature of an off-axis jet afterglow. While we cannot rule out the existence of a jet pointing elsewhere, the observed gamma-rays could not have originated from such a jet. Instead, the radio data requires a mildly relativistic wide-angle outflow moving towards us. This outflow could be the high velocity tail of the neutron-rich material dynamically ejected during the merger or a cocoon of material that breaks out when a jet transfers its energy to the dynamical ejecta. The cocoon scenario can explain the radio light curve of GW170817 as well as the gamma-rays and X-rays (possibly also ultraviolet and optical emission), and hence is the model most consistent with the observational data. Cocoons may be a ubiquitous phenomenon produced in neutron star mergers, giving rise to a heretofore unidentified population of radio, ultraviolet, X-ray and gamma-ray transients in the local universe.

Gravitational-wave observation together with a large number of electromagnetic observations shows that the source of the latest gravitational-wave event, GW170817, detected primarily by advanced LIGO, is the merger of a binary neutron star. We attempt to interpret this observational event based on our results of numerical-relativity simulations performed so far paying particular attention to the optical and infra-red observations. We finally reach a conclusion that this event is described consistently by the presence of a long-lived hypermassive or supramassive neutron star as the merger remnant, because (i) significant contamination by lanthanide elements along our line of sight to this source can be avoided by the strong neutrino irradiation from it and (ii) it could play a crucial role to produce an ejecta component of appreciable mass with fast motion in the post-merger phase. We also point out that (I) the neutron-star equation of state has to be sufficiently stiff (i.e., the maximum mass of cold spherical neutron stars, M_max, has to be appreciably higher than 2M_sun in order that a long-lived massive neutron star can be formed as the merger remnant for the binary systems of GW170817, for which the initial total mass is >~ 2.73M_sun and (II) no detection of relativistic optical counterpart suggests a not-extremely high value of M_max approximately as 2.15-2.25M_sun.

We study optical counterparts of a new-born pulsar in a double neutron star system like PSR J0737-3039A/B. This system is believed to eject a small amount of mass of $\mathcal{O}(0.1M_{\odot})$ at the second core-collapse supernova. We argue that the initial spin of the new-born pulsar can be determined by the orbital period at the time when the second supernova occurs. The spin angular momentum of the progenitor is expected to be similar to that of the He-burning core, which is tidally synchronized with the orbital motion, and then the second remnant may be born as a millisecond pulsar. If the dipole magnetic field strength of the nascent pulsar is comparable to that inferred from the current spin-down rate of PSR J0737-3039B, the initial spin-down luminosity is comparable to the luminosity of super-luminous supernovae. We consider thermal emission arising from the supernova ejecta driven by the relativistic wind from such a new-born pulsar. The resulting optical light curves have a rising time $\sim 10$ days and peak luminosity $\sim 10^{44}$ erg/s. The optical emission may last for a month to several months, due to the reprocessing of X-rays and UV photons via photoelectric absorption. These features are broadly consistent with those of the rapidly-rising optical transients. The high spin-down luminosity and small ejecta mass are favorable for the progenitor of the repeating fast radio burst, FRB 121102. We discuss a possible connection between newborn double pulsars and fast radio bursts.

Merging neutron stars offer an exquisite laboratory for simultaneously studying strong-field gravity and matter in extreme environments. We establish the physical association of an electromagnetic counterpart EM170817 to gravitational waves (GW170817) detected from merging neutron stars. By synthesizing a panchromatic dataset, we demonstrate that merging neutron stars are a long-sought production site forging heavy elements by r-process nucleosynthesis. The weak gamma-rays seen in EM170817 are dissimilar to classical short gamma-ray bursts with ultra-relativistic jets. Instead, we suggest that breakout of a wide-angle, mildly-relativistic cocoon engulfing the jet elegantly explains the low-luminosity gamma-rays, the high-luminosity ultraviolet-optical-infrared and the delayed radio/X-ray emission. We posit that all merging neutron stars may lead to a wide-angle cocoon breakout; sometimes accompanied by a successful jet and sometimes a choked jet.

Gravitational waves have been detected from a binary neutron star merger event, GW170817. The detection of electromagnetic radiation from the same source has shown that the merger occurred in the outskirts of the galaxy NGC 4993, at a distance of 40 megaparsecs from Earth. We report the detection of a counterpart radio source that appears 16 days after the event, allowing us to diagnose the energetics and environment of the merger. The observed radio emission can be explained by either a collimated ultra-relativistic jet viewed off-axis, or a cocoon of mildly relativistic ejecta. Within 100 days of the merger, the radio light curves will distinguish between these models and very long baseline interferometry will have the capability to directly measure the angular velocity and geometry of the debris.

With the first direct detection of merging black holes in 2015, the era of gravitational wave (GW) astrophysics began. A complete picture of compact object mergers, however, requires the detection of an electromagnetic (EM) counterpart. We report ultraviolet (UV) and X-ray observations by Swift and the Nuclear Spectroscopic Telescope ARray (NuSTAR) of the EM counterpart of the binary neutron star merger GW170817. The bright, rapidly fading ultraviolet emission indicates a high mass ($\approx0.03$ solar masses) wind-driven outflow with moderate electron fraction ($Y_{e}\approx0.27$). Combined with the X-ray limits, we favor an observer viewing angle of $\approx 30^{\circ}$ away from the orbital rotation axis, which avoids both obscuration from the heaviest elements in the orbital plane and a direct view of any ultra-relativistic, highly collimated ejecta (a gamma-ray burst afterglow).

One of the puzzles in the recent observations of gravitational waves from binary black hole mergers is the observed low (projected) spins of the progenitor black holes. In two of the four events, GW150914, and the recent GW170104, the observed spins are most likely negative (but consistent with zero). In the third case LVT151012 it is practically zero and only in the forth case, GW151226, the spin is positive but low. These observations are puzzling within the field binary scenario in which positive higher spins are expected. Considering the most favorable Wolfe Rayet (WR) progenitors we estimate the expected spin distribution for different evolution scenarios and compare it to the observations. With typical parameters one expects a significant fraction ($\ge 25\%$) of the mergers to have high effective spin values. However due to uncertainties in the outcome of the common envelope phase (typical separation and whether the stars are rotating or not) and in the late stages of massive star evolution (the strength of the winds) one cannot rule our scenarios in which the expected spins are low. While observations of high effective spin events will support this scenario, further observations of negative spin events would rule it out.

We explore the implications of the observed low spins aligned with the orbital axis in Advanced~LIGO~O1 run on binary black hole (BBH)~merger scenarios in which the merging BBHs have evolved from field binaries. The coalescence time determines the initial orbital separation of BBHs. This, in turn, determines whether the stars are synchronized before collapse and hence determines their projected spins. Short coalescence times imply synchronization and large spins. Among known stellar objects, Wolf-Rayet~(WR) stars seem the only progenitors consistent with the low aligned spins observed in LIGO's~O1 provided that the orbital axis maintains its direction during the collapse. We calculate the spin distribution of BBH mergers in the local Universe and its redshift evolution for WR progenitors. Assuming that the BBH formation rate peaks around a redshift of $\sim 2\,$--$\,3$, we show that BBH mergers in the local Universe are dominated by low spin events. The high spin population starts to dominate at a redshift of $\sim0.5\,$--$\,1.5$. WR stars are also progenitors of long Gamma-Ray~Bursts~(LGRBs) that take place at a comparable rate to BBH mergers. We discuss the possible connection between the two phenomena. Additionally, we show that hypothetical Population~III star progenitors are also possible. Although WR and Population~III progenitors are consistent with the current data, both models predict a non-vanishing fraction of high aligned spin BBH mergers. If those are not detected within the coming LIGO/Virgo runs, it will be unlikely that the observed BBHs formed via field binaries.

We study formation of stellar mass binary black holes (BBHs) originating from Population III (PopIII) stars, performing stellar evolution simulations for PopIII binaries with MESA. We find that a significant fraction of PopIII binaries form massive BBHs through stable mass transfer between two stars in a binary, without experiencing common envelope phases. We investigate necessary conditions required for PopIII binaries to form BBHs coalescing within the Hubble time with a semi-analytical model calibrated by the stellar evolution simulations. The formation efficiency of coalescing PopIII BBHs is estimated for two different initial conditions for PopIII binaries with large and small separations, respectively. Consequently, in both models, $\sim 10\%$ of the total PopIII binaries form BBHs only through stable mass transfer and $\sim 10\%$ of these BBHs merge due to gravitational wave emission within the Hubble time. Furthermore, the chirp mass of merging BBHs has a flat distribution over $15\lesssim M_{\rm chirp}/M_\odot \lesssim 35$. This formation pathway of PopIII BBHs is presumably robust because stable mass transfer is less uncertain than common envelope evolution, which is the main formation channel for Population II BBHs. We also test the hypothesis that the BBH mergers detected by LIGO originate from PopIII stars using our result and the total number of PopIII stars formed in the early universe as inferred from the optical depth measured by Planck. We conclude that the PopIII BBH formation scenario can explain the mass-weighted merger rate of the LIGO's O1 events with the maximal PopIII formation efficiency inferred from the Planck measurement, even without BBHs formed by unstable mass transfer or common envelope phases.

Macronovae (kilonovae) that arise in binary neutron star mergers are powered by radioactive beta decay of hundreds of $r$-process nuclides. We derive, using Fermi's theory of beta decay, an analytic estimate of the nuclear heating rate. We show that the heating rate evolves as a power law ranging between $t^{-6/5}$ to $t^{-4/3}$. The overall magnitude of the heating rate is determined by the mean values of nuclear quantities, e.g., the nuclear matrix elements of beta decay. These values are specified by using nuclear experimental data. We discuss the role of higher order beta transitions and the robustness of the power law. The robust and simple form of the heating rate suggests that observations of the late-time bolometric light curve $\propto t^{-\frac{4}{3}}$ would be a direct evidence of a $r$-process driven macronova. Such observations could also enable us to estimate the total amount of $r$-process nuclei produced in the merger.

Ultra-long gamma-ray bursts (ulGRBs) are a new population of GRBs with extreme durations of $\sim 10^{4}$ s. Leading candidates for their origin are blue supergiant collapsars, magnetars, and white dwarf tidal disruption events (WD-TDEs) caused by massive black holes (BHs). Recent observations of supernova-like (SN-like) bumps associated with ulGRBs challenged both the WD-TDE and the blue supergiant models because of the detection of SNe and the absence of hydrogen lines, respectively. We propose that WD-TDEs can accommodate the observed SN-like bumps if the fallback WD matter releases energy into the unbound WD ejecta. The observed ejecta energy, luminosity, and velocity are explained by the gravitational energy, Eddington luminosity, and escape velocity of the formed accretion disk, respectively. We also show that the observed X-rays can ionize the ejecta, eliminating lines. The SN-like light curves (SN 2011kl) for the ulGRB 111209A are consistent with all three models, although a magnetar model is unnatural because the spin-down time required to power the SN-like bump is a hundred times longer than the GRB. Our results imply that TDEs are a possible energy source for SN-like events in general and for ulGRBs in particular.

Recent observations of $r$-process material in ultra-faint dwarf galaxies (UFDs) shed light on the sources of these elements. Strong upper limits on the Eu mass in some UFDs combined with detections of much larger masses in a UFD, Reticulum II, and other dwarf galaxies imply that Eu production is dominated by rare events, and that the minimal Eu mass observed in any UFD is approximately the amount of Eu mass produced per event. This is consistent with other independent observations in the Galaxy. We estimate, using a model independent likelihood analysis, the rate and Eu (Fe) mass produced per $r$-process (Fe production) event in dwarf galaxies including classical dwarfs and UFDs. The mass and rate of the Fe production events are consistent with the normal core-collapse supernova~(ccSN) scenario. The Eu mass per event is $3\times 10^{-5}M_{\odot}<\tilde{m}_{\rm Eu}<2\times 10^{-4}M_{\odot}$, corresponding to a total $r$-process mass per event of $6\times 10^{-3}M_{\odot}<\tilde{m}_{r-process}<4\times 10^{-2}M_{\odot}$. The rate of $r$-process events is $2.5\times 10^{-4}<R_{rp/SN}<1.4\times 10^{-3}$ as compared with the ccSNe rate. These values are consistent with the total Eu mass observed in our own Galaxy, suggesting that the same mechanism is behind the production of $r$-process events in both dwarf galaxies and the Milky Way, and that it may be the dominant mechanism for production of $r$-process elements in the Universe. The results are consistent with neutron star mergers estimates, but cannot rule out other rare core collapse scenarios, provided that they produce a significant amounts of $r$-process material per event.

Mergers of binary neutron stars and black hole-neutron star binaries produce gravitational-wave (GW) emission and outflows with significant kinetic energies. These outflows result in radio emissions through synchrotron radiation. We explore the detectability of these synchrotron generated radio signals by follow-up observations of GW merger events lacking a detection of electromagnetic counterparts in other wavelengths. We model radio light curves arising from (i) sub-relativistic merger ejecta and (ii) ultra-relativistic jets. The former produces radio remnants on timescales of a few years and the latter produces $\gamma$-ray bursts in the direction of the jet and orphan-radio afterglows extending over wider angles on timescales of weeks. Based on the derived light curves, we suggest an optimized survey at $1.4$ GHz with five epochs separated by a logarithmic time interval. We estimate the detectability of the radio counterparts of simulated GW-merger events to be detected by advanced LIGO and Virgo by current and future radio facilities. The detectable distances for these GW merger events could be as high as 1 Gpc. $20$--$60\%$ of the long-lasting radio remnants will be detectable in the case of the moderate kinetic energy of $3\cdot 10^{50}$ erg and a circum-merger density of $0.1 {\rm cm^{-3}}$ or larger, while $5$--$20\%$ of the orphan radio afterglows with kinetic energy of $10^{48}$ erg will be detectable. The detection likelihood increases if one focuses on the well-localizable GW events. We discuss the background noise due to radio fluxes of host galaxies and false positives arising from extragalactic radio transients and variable Active Galactic Nuclei and we show that the quiet radio transient sky is of great advantage when searching for the radio counterparts.

Merging neutron star binaries are prime candidate sources for heavy r-process nucleosynthesis. The amount of heavy r-process material is consistent with the mass ejection and rates of mergers, and abundances of relic radioactive materials suggest that heavy r-process material is produced in rare events. Observations of possible macronovae provide further support for this model. Still, some concerns remain. One is the observation of heavy r-process elements in Ultra Faint Dwarf (UFD) galaxies. The escape velocities from UFDs are so small that the natal kicks, taking place at neutron stars birth, might eject such binaries from UFDs. Furthermore the old stellar populations of UFDs requires that r-process nucleosynthesis must have taken place very early on, while it may take several Gyr for compact binaries to merge. This last problem arises also within the Milky Way where heavy r-process materials has been observed in some low metallicity stars. We show here that since a significant fraction of neutron star binaries form with a very small proper motion those won't be ejected even from a UFD. Furthermore, a significant fraction of these binaries merge within a few hundred Myr and in some cases the merger time is even shorter than a hundred Myr. This population of "rapid mergers" explains the appearance of heavy r-process material in both UFDs and in the Milky Way.

GRB 050709 was the first short Gamma-ray Burst (sGRB) with an identified optical afterglow. In this work we report a re-analysis of the publicly available data of this event and the discovery of a Li-Paczynski macronova/kilonova that dominates the optical/IR signal at t $>$ 2.5 days. Such a signal would arise from ~0.05 M$_\odot$ r-process material launched by a compact binary merger. The implied mass ejection supports the suggestion that compact binary mergers are significant and possibly main sites of heavy r-process nucleosynthesis. We have re-analyzed all afterglow data from nearby short and hybrid GRBs. A statistical study of sGRB/macronova connection reveals that macronova may have taken place in all these GRBs though the fraction as low as 0.18 cannot be ruled out. The identification of two of the three macronova candidates in the I-band implies a more promising detection prospect for ground-based surveys.

Extracting the unique information on ultradense nuclear matter from the gravitational waves emitted by merging, neutron-star binaries requires robust theoretical models of the signal. We develop a novel effective-one-body waveform model that includes, for the first time, dynamic (instead of only adiabatic) tides of the neutron star, as well as the merger signal for neutron-star--black-hole binaries. We demonstrate the importance of the dynamic tides by comparing our model against new numerical-relativity simulations of nonspinning neutron-star--black-hole binaries spanning more than 24 gravitational-wave cycles, and to other existing numerical simulations for double neutron-star systems. Furthermore, we derive an effective description that makes explicit the dependence of matter effects on two key parameters: tidal deformability and fundamental oscillation frequency.

Compact binary mergers may have already been observed as they are the leading model for short gamma-ray bursts (sGRBs). Radioactive decay within the ejecta from these mergers is expected to produce an infra-red flare, dubbed macronova (or kilonova), on a time scale of a week. Recently two such macronova candidates were identified in followup observations of sGRBs, strengthening the possibility that those indeed arise from mergers. The same ejecta will also produce a long term (months to years) radio emission due to its interaction with the surrounding ISM. In search for this emission, we observed the two macronova candidates, GRB 130603B and GRB 060614 with the Jansky very large array (VLA) and the Australia Telescope Compact Array (ATCA). Our observations resulted in null-detections, putting strong upper limits on the kinetic energy and mass of the ejecta. A possible outcome of a merger is a highly magnetized neutron star (a magnetar), which has been suggested as the central engine for GRBs. Such a magnetar will deposit a significant fraction of its energy into the ejecta leading to a brighter radio flare. Our results, therefore, rule out magnetars in these two events.

The radioactive decay of the freshly synthesized $r$-process nuclei ejected in compact binary mergers power optical/infrared macronovae (kilonovae) that follow these events. The light curves depend critically on the energy partition among the different products of the radioactive decay and this plays an important role in estimates of the amount of ejected $r$-process elements from a given observed signal. We study the energy partition and $\gamma$-ray emission of the radioactive decay. We show that $20$-$50\%$ of the total radioactive energy is released in $\gamma$-rays on timescales from hours to a month. The number of emitted $\gamma$-rays per unit energy interval has roughly a flat spectrum between a few dozen keV and $1$ MeV so that most of this energy is carried by $\sim 1$ MeV $\gamma$-rays. However at the peak of macronova emission the optical depth of the $\gamma$-rays is $\sim 0.02$ and most of the $\gamma$-rays escape. The loss of these $\gamma$-rays reduces the heat deposition into the ejecta and hence reduces the expected macronova signals if those are lanthanides dominated. This implies that the ejected mass is larger by a factor of $2$-$3$ than what was previously estimated. Spontaneous fission heats up the ejecta and the heating rate can increase if a sufficient amount of transuranic nuclei are synthesized. Direct measurements of these escaping $\gamma$-rays may provide the ultimate proof for the macronova mechanisms and an identification of the $r$-process nucleosynthesis sites. However, the chances to detect these signals are slim with current X-ray and $\gamma$-ray missions. New detectors, more sensitive by at least a factor of ten, are needed for a realistic detection rate.

Measurements of the radioactive $^{244}$Pu abundances can break the degeneracy between high-rate/low-yield and low-rate/high-yield scenarios for the production of heavy $r$-process elements. The first corresponds to production by core collapse supernovae (cc-SNe) while the latter corresponds to production by e.g. compact binary mergers. The estimated $^{244}$Pu abundance in the current interstellar medium inferred from deep-sea measurements (Wallner et al. 2015) is significantly lower than that corresponding Early Solar System abundances (Turner et al 2007). We estimate the expected median value of the $^{244}$Pu abundances and fluctuations around this value in both models. We show that while the current and Early Solar System abundances are naturally explained within the low-rate/high-yield (e.g. merger) scenario, they are incompatible with the high-rate/low-yield (cc-SNe) model. The inferred event rate remarkably agrees with compact binary merger rates estimated from Galactic neutron star binaries and from short gamma-ray bursts. Furthermore, the ejected mass of $r$-process elements per event agrees with both theoretical and observational macronova/kilonova estimates.

Long-duration ($>2$ s) $\gamma$-ray bursts that are believed to originate from the death of massive stars are expected to be accompanied by supernovae. GRB 060614, that lasted 102 s, lacks a supernova-like emission down to very stringent limits and its physical origin is still debated. Here we report the discovery of near-infrared bump that is significantly above the regular decaying afterglow. This red bump is inconsistent with even the weakest known supernova. However, it can arise from a Li-Paczy\'{n}ski macronova $-$ the radioactive decay of debris following a compact binary merger. If this interpretation is correct GRB 060614 arose from a compact binary merger rather than from the death of a massive star and it was a site of a significant production of heavy r-process elements. The significant ejected mass favors a black hole-neutron star merger but a double neutron star merger cannot be ruled out.

We perform new longterm (15-16 orbits) simulations of coalescing binary neutron stars in numerical relativity using an updated Einstein's equation solver, employing low-eccentricity initial data, and modeling the neutron stars by a piecewise polytropic equation of state. A convergence study shows that our new results converge more rapidly than the third order and using the determined convergence order, we construct an extrapolated waveform for which the estimated total phase error should be less than 1 radian. We then compare the extrapolated waveforms with those calculated by the latest effective-one-body (EOB) formalism in which the so-called tidal deformability, higher post-Newtonian corrections, and gravitational self-force effects are taken into account. We show that for a binary of compact neutron stars with their radius 11.1 km, the waveform by the EOB formalism agrees quite well with the numerical waveform so that the total phase error is smaller than 1 radian for the total phase of $\sim 200$ radian up to the merger. By contrast, for a binary of less compact neutron stars with their radius 13.6 km, the EOB and numerical waveforms disagree with each other in the last few wave cycles, resulting in the total phase error of $\sim 3$ radian.

In addition to producing a strong gravitational signal, a short gamma-ray burst (GRB), and a compact remnant, neutron star mergers eject significant masses at significant kinetic energies. This mass ejection takes place via dynamical mass ejection and a GRB jet but other processes have also been suggested: a shock-breakout material, a cocoon resulting from the interaction of the jet with other ejecta, and viscous and neutrino driven winds from the central remnant or the accretion disk. The different components of the ejected masses include up to a few percent of a solar mass, some of which is ejected at relativistic velocities. The interaction of these ejecta with the surrounding interstellar medium will produce a long lasting radio flare, in a similar way to GRB afterglows or to radio supernovae. The relative strength of the different signals depends strongly on the viewing angle. An observer along the jet axis or close to it will detect a strong signal at a few dozen days from the radio afterglow (or the orphan radio afterglow) produced by the highly relativistic GRB jet. For a generic observer at larger viewing angles, the dynamical ejecta, whose contribution peaks a year or so after the event, will generally dominate. Depending on the observed frequency and the external density, other components may also give rise to a significant contribution. We also compare these estimates with the radio signature of the short GRB 130603B. The radio flare from the dynamical ejecta might be detectable with the EVLA and the LOFAR for the higher range of external densities $n\gtrsim 0.5{\rm cm^{-3}}$.

We perform population synthesis simulations for Population III (Pop III) coalescing binary neutron stars (NS-NSs), neutron star - black hole binaries (NS-BHs), and binary black holes (BH-BHs) which merge within the age of the universe. We found that the typical mass of Pop III BH-BHs is $\sim 30\ \rm{M}_{\odot}$ so that the inspiral chirp signal of gravitational waves can be detected up to z=0.28 by KAGRA, Adv. LIGO, Adv. Virgo and GEO network. Our simulations suggest that the detection rate of the coalescing Pop III BH-BHs is $140 (68) \cdot ({\rm SFR}_{\rm p}/10^{-2.5} \rm{M}_{\odot} {\rm yr}^{-1} {\rm Mpc}^{-3}) \cdot {\rm Err}_{\rm sys} ~{\rm events} ~{\rm yr}^{-1}$ for the flat (Salpeter) initial mass function (IMF), respectively, where $\rm SFR_p$ and $\rm Err_{sys}$ are the peak value of the Pop III star formation rate and the possible systematic errors due to the assumptions in Pop III population synthesis, respectively. $\rm Err_{sys}=1$ correspond to conventional parameters for Pop I stars. From the observation of the chirp signal of the coalescing Pop III BH-BHs, we can determine both the mass and the redshift of the binary for the cosmological parameters determined by Planck satellite. Our simulations suggest that the cumulative redshift distribution of the coalescing Pop III BH-BHs depends almost only on the cosmological parameters. We might be able to confirm the existence of Pop III massive stars of mass $\sim 30~\rm M_{\odot}$ by the detections of gravitational waves if the merger rate of the Pop III massive BH-BHs dominates that of Pop I BH-BHs.

The observations of jet breaks in the afterglows of short gamma-ray bursts (SGRBs) indicate that the jet has a small opening angle of < 10{\deg}. The collimation mechanism of the jet is a longstanding theoretical problem. We numerically analyze the jet propagation in the material ejected by double neutron star merger, and demonstrate that if the ejecta mass is > 10^{-2} M_{sun}, the jet is well confined by the cocoon and emerges from the ejecta with the required collimation angle. Our results also suggest that there are some populations of choked (failed) SGRBs or low-luminous new types of event. By constructing a model for SGRB 130603B, which is associated with the first kilonova/macronova can- didate, we infer that the equation-of-state of neutron stars would be soft enough to provide sufficient ejecta to collimate the jet, if this event was associated with a double neutron star merger.

Detection of electromagnetic counterparts of gravitational wave (GW) sources is important to unveil the nature of compact binary coalescences. We perform three-dimensional, time-dependent, multi-frequency radiative transfer simulations for radioactively powered emission from the ejecta of black hole (BH) - neutron star (NS) mergers. Depending on the BH to NS mass ratio, spin of the BH, and equations of state of dense matter, BH-NS mergers can eject more material than NS-NS mergers. In such cases, radioactively powered emission from the BH-NS merger ejecta can be more luminous than that from NS-NS mergers. We show that, in spite of the expected larger distances to BH-NS merger events, observed brightness of BH-NS mergers can be comparable to or even higher than that of NS-NS mergers. We find that, when the tidally disrupted BH-NS merger ejecta are confined to a small solid angle, the emission from BH-NS merger ejecta tends to be bluer than that from NS-NS merger ejecta for a given total luminosity. Thanks to this property, we might be able to distinguish BH-NS merger events from NS-NS merger events by multi-band observations of the radioactively powered emission. In addition to the GW observations, such electromagnetic observations can potentially provide independent information on the nature of compact binary coalescences.

An electromagnetic transient powered by the radioactive decay of r-process elements, a so-called kilonova/macronova, is one of the possible observable consequences of compact binary mergers including at east one neutron star. Recent observations strongly suggest the first discovery of the electromagnetic transient, which is associated with the short GRB 130603B. We explore a possible progenitor of this event combining the numerical-relativity simulations and radiative transfer simulations of the dynamical ejecta of binary neutron star and black hole - neutron star mergers. We show that the ejecta models within a realistic parameter range consistently reproduce the observed near-infrared excess. We also show that the soft equation of state models for binary neutron star mergers and the stiff equation of state models for black hole - neutron star mergers are favored to reproduce the observed luminosity.

Massive (hypermassive and supramassive) neutron stars are likely to be often formed after the merger of binary neutron stars. We explore the evolution process of the remnant massive neutron stars and gravitational waves emitted by them, based on numerical-relativity simulations for binary neutron star mergers employing a variety of equations of state and choosing a plausible range of the neutron-star mass of binaries. We show that the lifetime of remnant hypermassive neutron stars depends strongly on the total binary mass and also on the equations of state. Gravitational waves emitted by the remnant massive neutron stars universally have a quasiperiodic nature of an approximately constant frequency although the frequency varies with time. We also show that the frequency and time-variation feature of gravitational waves depend strongly on the equations of state. We derive a fitting formula for the quasiperiodic gravitational waveforms, which may be used for the data analysis of a gravitational-wave signal.

The merger of binary neutron stars (NSs) is among the most promising gravitational wave (GW) sources. Next-generation GW detectors are expected to detect signals from the NS merger within 200 Mpc. Detection of electromagnetic wave (EM) counterpart is crucial to understand the nature of GW sources. Among possible EM emission from the NS merger, emission powered by radioactive r-process nuclei is one of the best targets for follow-up observations. However, prediction so far does not take into account detailed r-process element abundances in the ejecta. We perform radiative transfer simulations for the NS merger ejecta including all the r-process elements from Ga to U for the first time. We show that the opacity in the NS merger ejecta is about kappa = 10 cm^2 g^{-1}, which is higher than that of Fe-rich Type Ia supernova ejecta by a factor of ~ 100. As a result, the emission is fainter and longer than previously expected. The spectra are almost featureless due to the high expansion velocity and bound-bound transitions of many different r-process elements. We demonstrate that the emission is brighter for a higher mass ratio of two NSs and a softer equation of states adopted in the merger simulations. Because of the red color of the emission, follow-up observations in red optical and near-infrared (NIR) wavelengths will be the most efficient. At 200 Mpc, expected brightness of the emission is i = 22 - 25 AB mag, z = 21 - 23 AB mag, and 21 - 24 AB mag in NIR JHK bands. Thus, observations with wide-field 4m- and 8m-class optical telescopes and wide-field NIR space telescopes are necessary. We also argue that the emission powered by radioactive energy can be detected in the afterglow of nearby short gamma-ray bursts.

We study gravitational waves emitted in the late inspiral stage of binary neutron stars by analyzing the waveform obtained in numerical-relativity simulations. For deriving the physical gravitational waveforms from the numerical results, the resolution extrapolation plays an essential role for our simulations. The extrapolated gravitational-wave phases are compared with those calculated in the post-Newtonian (PN) and effective-one-body (EOB) formalisms including corrections of tidal effects. We show that the extrapolated gravitational-wave phases in numerical relativity agree well with those by the PN and EOB calculations for most of the inspiral stage except for a tidally-dominated, final inspiral stage, in which the PN and EOB results underestimate the tidal effects. Nevertheless, the accumulated phase difference between our extrapolated results and the results by the PN/EOB calculations is at most 1--3 radian in the last 15 cycles.

Numerical-relativity simulations for the merger of binary neutron stars are performed for a variety of equations of state (EOSs) and for a plausible range of the neutron-star mass, focusing primarily on the properties of the material ejected from the system. We find that a fraction of the material is ejected as a mildly relativistic and mildly anisotropic outflow with the typical and maximum velocities $\sim 0.15$ -- $0.25c$ and $\sim 0.5$ -- $0.8c$ (where $c$ is the speed of light), respectively, and that the total ejected rest mass is in a wide range $10^{-4}$ -- $10^{-2}M_{\odot}$, which depends strongly on the EOS, the total mass, and the mass ratio. The total kinetic energy ejected is also in a wide range between $10^{49}$ and $10^{51} {\rm ergs}$. The numerical results suggest that for a binary of canonical total mass $2.7M_{\odot}$, the outflow could generate an electromagnetic signal observable by the planned telescopes through the production of heavy-element unstable nuclei via the $r$-process or through the formation of blast waves during the interaction with the interstellar matter, if the EOS and mass of the binary are favorable ones.

We perform a numerical-relativity simulation for the merger of binary neutron stars with 6 nuclear-theory-based equations of state (EOSs) described by piecewise polytropes. Our purpose is to explore the dependence of the dynamical behavior of the binary neutron star merger and resulting gravitational waveforms on the EOS of the supernuclear-density matter. The numerical results show that the merger process and the first outcome are classified into three types; (i) a black hole is promptly formed, (ii) a short-lived hypermassive neutron star (HMNS) is formed, (iii) a long-lived HMNS is formed. The type of the merger depends strongly on the EOS and on the total mass of the binaries. For the EOS with which the maximum mass is larger than 2Msun, the lifetime of the HMNS is longer than 10 ms for a total mass m_0=2.7Msun. A recent radio observation suggests that the maximum mass of spherical neutron stars is M_max \geq 1.97\pm 0.04Msun in one \sigma level. This fact and our results support the possible existence of a HMNS soon after the onset of the merger for a typical binary neutron star with m_0=2.7Msun. We also show that the torus mass surrounding the remnant black hole is correlated with the type of the merger process; the torus mass could be large, \geq 0.1Msun, in the case that a long-lived HMNS is formed. We also show that gravitational waves carry information of the merger process, the remnant, and the torus mass surrounding a black hole.