The recent discovery of electromagnetic signals in coincidence with gravitational waves from neutron-star mergers has solidified the importance of multimessenger campaigns for studying the most energetic astrophysical events. Pioneering multimessenge
r observatories, such as the LIGO/Virgo gravitational wave detectors and the IceCube neutrino observatory, record many candidate signals that fall short of the detection significance threshold. These sub-threshold event candidates are promising targets for multimessenger studies, as the information provided by these candidates may, when combined with time-coincident gamma-ray observations, lead to significant detections. In this contribution, I describe our use of sub-threshold binary neutron star merger candidates identified in Advanced LIGOs first observing run (O1) to search for transient events in very-high-energy gamma rays using archival observations from the VERITAS imaging atmospheric Cherenkov telescope array. I describe the promise of this technique for future joint sub-threshold searches.
The observed distributions of the source properties from gravitational-wave detections are biased due to the selection effects and detection criteria in the detections, analogous to the Malmquist bias. In this work, this observation bias is investiga
ted through its fundamental statistical and physical origins. An efficient semi-analytical formulation for its estimation is derived which is as accurate as the standard method of numerical simulations, with only a millionth of the computational cost. Then, the estimated bias is used for model independent inferences on the binary black hole population. These inferences show additional structures, specifically two potential mass gaps in the joint mass distribution, which were not found via modelled inferences. Example ready-to-use scripts and some produced datasets for this method are shared in an online repository.
Recent gravitational wave (GW) observations by LIGO/Virgo show evidence for hierarchical mergers, where the merging BHs are the remnants of previous BH merger events. These events may carry important clues about the astrophysical host environments of
the GW sources. In this paper, we present the distributions of the effective spin parameter ($chi_mathrm{eff}$), the precession spin parameter ($chi_mathrm{p}$), and the chirp mass ($m_mathrm{chirp}$) expected in hierarchical mergers. Under a wide range of assumptions, hierarchical mergers produce (i) a monotonic increase of the average of the typical total spin for merging binaries, which we characterize with ${bar chi}_mathrm{typ}equiv overline{(chi_mathrm{eff}^2+chi_mathrm{p}^2)^{1/2}}$, up to roughly the maximum $m_mathrm{chirp}$ among first-generation (1g) BHs, and (ii) a plateau at ${bar chi}_mathrm{typ}sim 0.6$ at higher $m_mathrm{chirp}$. We suggest that the maximum mass and typical spin magnitudes for 1g BHs can be estimated from ${bar chi}_mathrm{typ}$ as a function of $m_mathrm{chirp}$. The GW data observed in LIGO/Virgo O1--O3a prefers an increase in ${bar chi}_mathrm{typ}$ at low $m_mathrm{chirp}$, which is consistent with the growth of the BH spin magnitude by hierarchical mergers, at $sim 2 sigma$ confidence. A Bayesian analysis suggests that 1g BHs have the maximum mass of $sim 15$--$30,M_odot$ if the majority of mergers are of high-generation BHs (not among 1g-1g BHs), which is consistent with mergers in active galactic nucleus disks and/or nuclear star clusters, while if mergers mainly originate from globular clusters, 1g BHs are favored to have non-zero spin magnitudes of $sim 0.3$. We also forecast that signatures for hierarchical mergers in the ${bar chi}_mathrm{typ}$ distribution can be confidently recovered once the number of GW events increases to $gtrsim O(100)$.
We propose a novel scenario for possible electromagnetic (EM) emission by compact binary mergers in the accretion disks of active galactic nuclei (AGNs). Nuclear star clusters in AGNs are a plausible formation site of compact-stellar binaries (CSBs)
whose coalescences can be detected through gravitational waves (GWs). We investigate the accretion onto and outflows from CSBs embedded in AGN disks. We show that these outflows are likely to create outflow cavities in the AGN disks before the binaries merge, which makes EM or neutrino counterparts much less common than would otherwise be expected. We discuss the necessary conditions for detectable EM counterparts to mergers inside the outflow cavities. If the merger remnant black hole experiences a high recoil velocity and can enter the AGN disk, it can accrete gas with a super-Eddington rate, newly forming a cavity-like structure. This bubble can break out of the disk within a day to a week after the merger. Such breakout emission can be bright enough to be detectable by current soft X-ray instruments, such as Swift-XRT and Chandra.
The stellar mass binary black hole (sBBH) mergers presently detected by LIGO may originate wholly or in part from binary black hole mergers embedded in disks of gas around supermassive black holes. Determining the contribution of these active galacti
c nucleus (AGN) disks to the sBBH merger rate enables us to uniquely measure important parameters of AGN disks, including their typical density, aspect ratio, and lifetime, thereby putting unique limits on an important element of galaxy formation. For the first time, gravitational waves will allow us to reveal the properties of the hidden interior of AGN disks, while electromagnetic radiation (EM) probes the disk photosphere. The localization of sBBH merger events from LIGO is generally insufficient for association with a single EM counterpart. However, the contribution to the LIGO event rate from rare source types (such as AGNs) can be determined on a statistical basis. To determine the contribution to the sBBH rate from AGNs in the next decade requires: {em 1) a complete galaxy catalog for the LIGO search volume, 2) strategic multi-wavelength EM follow-up of LIGO events and 3) significant advances in theoretical understanding of AGN disks and the behavior of objects embedded within them.}
Fast-spinning strongly magnetized newborn neutron stars, including nascent magnetars, are popularly implemented as the engine of luminous stellar explosions. Here, we consider the scenario that they power various stripped-envelope supernovae, not onl
y super-luminous supernovae Ic but also broad-line supernovae Ibc and possibly some ordinary supernovae Ibc. This scenario is also motivated by the hypothesis that Galactic magnetars largely originate from fast-spinning neutron stars as remnants of stripped-envelope supernovae. By consistently modeling the energy injection from magnetized wind and Ni decay, we show that proto-neutron stars with >~ 10 ms rotation and B_dip >~ 5 x 10^14 G can be harbored in ordinary supernovae Ibc. On the other hand, millisecond proto-neuton stars can solely power broad-line supernovae Ibc if they are born with poloidal magnetic field of B_dip >~ 5 x 10^14 G, and superluminous supernovae Ic with B_dip >~ 10^13 G. Then, we study how multi-messenger emission can be used to discriminate such pulsar-driven supernova models from other competitive scenarios. First, high-energy x-ray and gamma-ray emission from embryonic pulsar wind nebulae is a promising smoking gun of the underlying newborn pulsar wind. Follow-up observations of stripped-envelope supernovae using NuSTAR ~ 50-100 days after the explosion is strongly encouraged for nearby objects. We also discuss possible effects of gravitational-waves on the spin-down of proto-neutron stars. If millisecond proto-neutron stars with B_dip <~ a few x 10^13 G emit gravitational waves through e.g., non-axisymmetric rotation deformed by the inner toroidal fields of B_t >~ 10^16 G, the gravitational wave signal can be detectable from ordinary supernova Ibc in the Virgo cluster by Advanced LIGO, Advanced Virgo, and KAGRA.
Kilonovae represent an important electromagnetic counterpart for compact binary mergers, which could become the most commonly detected gravitational wave (GW) source. Follow-up observations, triggered by GW events, of kilonovae are nevertheless diffi
cult due to poor localization by GW detectors and due to their faint near-infrared peak emission that has limited observational capability. We show that the Near-Infrared Camera (NIRCam) on the James Webb Space Telescope (JWST) will be able to detect kilonovae within the relevant GW-detection range of $sim$ 200 Mpc in short ($lesssim$ 12-second) exposure times for a week following the merger. Despite this sensitivity, a kilonova search fully covering a fiducial localized area of $10$ $mbox{deg}^2$ will not be viable with NIRCam due to its limited field of view. However, targeted surveys may be developed to optimize the likelihood of discovering kilonovae efficiently within limited observing time. We estimate that a survey of $10$ $mbox{deg}^2$ focused on galaxies within 200 Mpc would require about 13 hours, dominated by overhead times; a survey further focused on galaxies exhibiting high star-formation rates would require $sim$ 5 hours. The characteristic time may be reduced to as little as $sim$4 hours, without compromising the likelihood of detecting kilonovae, by surveying sky areas associated with 50%, rather than 90%, confidence regions of 3 GW events, rather than a single event. On detection and identification of a kilonova, a limited number of NIRCam follow-up observations could constrain the properties of matter ejected by the binary and the equation of state of dense nuclear matter.
The first gravitational-wave (GW) observations will greatly benefit from the detection of coincident electromagnetic counterparts. Electromagnetic follow-ups will nevertheless be challenging for GWs with poorly reconstructed directions. GW source loc
alization can be inefficient (i) if only two GW observatories are in operation; (ii) if the detectors sensitivities are highly non-uniform; (iii) for events near the detectors horizon distance. For these events, follow-up observations will need to cover 100-1000 square degrees of the sky over a limited period of time, reducing the list of suitable telescopes. We demonstrate that the Cherenkov Telescope Array will be capable of following up GW event candidates over the required large sky area with sufficient sensitivity to detect short gamma-ray bursts, which are thought to originate from compact binary mergers, out to the horizon distance of advanced LIGO/Virgo. CTA can therefore be invaluable starting with the first multimessenger detections, even with poorly reconstructed GW source directions. This scenario also provides a further scientific incentive for GW observatories to further decrease the delay of their event reconstruction.
We present the baseline multimessenger analysis method for the joint observations of gravitational waves (GW) and high-energy neutrinos (HEN), together with a detailed analysis of the expected science reach of the joint search. The analysis method co
mbines data from GW and HEN detectors, and uses the blue-luminosity-weighted distribution of galaxies. We derive expected GW+HEN source rate upper limits for a wide range of source parameters covering several emission models. Using published sensitivities of externally triggered searches, we derive joint upper limit estimates both for the ongoing analysis with the initial LIGO-Virgo GW detectors with the partial IceCube detector (22 strings) HEN detector and for projected results to advanced LIGO-Virgo detectors with the completed IceCube (86 strings). We discuss the constraints these upper limits impose on some existing GW+HEN emission models.
