Observationally, fast radio bursts (FRBs) can be divided into repeating and apparently non-repeating (one-off) ones. It is unclear whether all FRBs repeat and whether there are genuine non-repeating FRBs. We attempt to address these questions using Monte Carlo simulations. We define a parameter $T_c$ at which the accumulated number of non-repeating sources becomes comparable to the total number of the repeating sources, which is a good proxy to denote the intrinsic repeater fraction among FRBs. Assuming that both types of sources exist, we investigate how the {\em observed} repeater fraction evolves with time for different parameters, including $T_c$, the repeating rate distribution power-law index $q$, the minimum repetition rate $r_{\rm 0,min}$, as well as the Weibull distribution index $k$ for repeating bursts. We find that unless $T_c \rightarrow \infty$ (i.e. there is no genuine non-repeating FRB source), the observed repeater fraction should increase with time first, reaching a peak, and then declines. The peak time $T_p$ and the peak fraction $F_{\rm r,obs,p}$ depend on $T_c$ and other repeating rate parameters. With the current data, one may pose a lower limit $T_c > (0.5-25)$ d for reasonable parameter values. Future continuous monitoring of FRBs using wide-field radio telescopes such as CHIME would measure or set a more stringent lower limit on $T_c$. The detection of a peak in the observed repeater fraction would disfavor the ansatz that ``all FRB sources repeat''.

The post-merger product of the first binary neutron star merger event detected in gravitational waves, GW170817, depends on neutron star equation of state (EoS) and is not well determined. We generally discuss the constraints one may pose on the maximum mass of a non-spinning neutron star, $M_{\rm TOV}$, based on the observations and some EoS-independent universal relations of rapidly-spinning neutron stars. If the merger product is a black hole after a brief hypermassive neutron star (HMNS) phase, we derive $M_{\rm TOV} < 2.09^{+0.11}_{-0.09}(^{+0.06}_{-0.04}) M_{\odot}$ at the 2$\sigma$ (1$\sigma$) level. The cases for a massive neutron star (MNS), either a supra-massive neutron star (SMNS) or even a stable neutron star (SNS), are also allowed by the data. We derive $2.09^{+0.11}_{-0.09}(^{+0.06}_{-0.04} M_{\odot}) \leq M_{\rm TOV}< 2.43^{+0.10}_{-0.08}(^{+0.06}_{-0.04}) M_{\odot}$ for the SMNS case and $M_{\rm TOV} \geq 2.43^{+0.10}_{-0.08}(^{+0.06}_{-0.04})M_{\odot}$ for the SNS case, at the $2\sigma$ ($1\sigma$) confidence level. In the MNS cases, we also discuss the constraints on the neutron star parameters (the dipolar magnetic field strength at the surface $B_p$ and the ellipticity $\epsilon$) that affect the spindown history, by considering different MNS survival times, e.g. 300 s, 1 d, and 155 d after the merger, as suggested by various observational arguments. We find that once an SMNS is formed, without violating the EM observational constraints, there always exist a set of ($B_p, \epsilon$) parameters that allow the SMNS to survive for 300s, 1 d, 155 d, or even longer.

In the literature, compact binary coalescences (CBCs) have been proposed as one of the main scenarios to explain the origin of some non-repeating fast radio bursts (FRBs). The large discrepancy between the FRB and CBC event rate densities suggest their associations, if any, should only apply at most for a small fraction of FRBs. Through a Bayesian estimation method, we show how a statistical analysis of the coincident associations of FRBs with CBC gravitational wave (GW) events may test the hypothesis of these associations. We show that during the operation period of advanced LIGO, the detection of $\sim100$ ($\sim1000$) GW-less FRBs with dispersion measure (DM) values smaller than 500 ${\rm pc~cm^{-3}}$ could reach the constraint that less than $10\%$ (or $1\%$) FRBs are related to binary black hole (BBH) mergers. The same number of FRBs with DM values smaller than 100 ${\rm pc~cm^{-3}}$ is required to reach the same constraint for binary neutron star (BNS) mergers. With the upgrade of GW detectors, the same constraints for BBH and BNS mergers can be reached with less FRBs or looser requirements for the DM values. It is also possible to pose constraints on the fraction of each type of CBCs that are able to produce observable FRBs based on the event density of FRBs and CBCs. This would further constrain the dimensionless charge of black holes in binary BH systems.

The LIGO and Virgo scientific collaboration (LVC) alerted two binary neutron star (BNS) merger candidates, S190425z and S190426c. Fermi-GBM observed 55.6\% (for S190425z) and 100\% (for S190426c) of the probability regions of both events at the respective merger times, but no gamma-ray burst (GRB) was detected in either case. The derived luminosity upper limits suggest that a short GRB similar to GRB 170817A would not be detectable for both cases due to their larger distances than GW170817. Assuming that the jet profile obtained from the GW170817/GRB 170817A is quasi-universal for all BNS-GRB associations, we derive that the viewing angles of S190425z and S190426c should be $> (0.11-0.41)$ and $> (0.09-0.39)$, respectively. Through Monte Carlo simulations, we show that with the GRB 170817A-like jet structure, all sky gamma-ray detectors, such as GBM and GECAM, are expected to detect $\sim 4.6\%$, $3.9\%$, $1.7\%$ and $6.6\%$, $5.7\%$, $2.8\%$ BNS mergers triggered by aLIGO, aLIGO A+ and ET, respectively. The joint detection fraction would be largely reduced for Swift-BAT, SVOM-ECLAIRS and Einstein Probe, whose sensitivities are better but fields of view are smaller.

With a selected sample of neutron star (NS) equation-of-states (EOSs) that are consistent with the current observations and have a range of maximum masses, we investigate the relations between NS gravitational mass $M_g$ and baryonic mass $M_b$, and the relations between the maximum NS mass supported through uniform rotation ($M_{\rm max}$) and that of nonrotating NSs ($M_{\rm TOV}$). We find that if one intends to apply an EOS-independent quadratic, universal transformation formula ($M_b=M_g+A\times M_{g}^2$) to all EOSs, the best fit $A$ value is 0.080 for non-rotating NSs only and 0.073 when different spin periods are considered. The residual error of the transformation is as large as $\sim0.1M_{\odot}$. For different EOSs, we find that the parameter $A$ for non-rotating NSs is proportional to $R_{1.4}^{-1}$ (where $R_{1.4}$ is NS radius for 1.4$M_\odot$ in unit of km). For a particular EOS, if one adopts the best-fit parameters for different spin periods, the residual error of the transformation is smaller, which is of the order of 0.01$M_\odot$ for the quadratic form and less than 0.01$M_\odot$ for the cubic form ($M_b=M_g+A_1\times M_{g}^2+A_2\times M_{g}^3$). We also find a very tight and general correlation between the normalized mass gain due to spin $\Delta m\equiv(M_{\rm max}-M_{\rm TOV})/M_{\rm TOV}$ and the spin period normalized to the Keplerian period ${\cal P}$, i.e. ${\rm log_{10}}\Delta m = (-2.74\pm0.05){\rm log_{10}}{\cal P}+{\rm log_{10}}(0.20\pm 0.01)$, which is independent of EOS models. Applications of our results to GW170817 is discussed.

The recent observations of GW170817 and its electromagnetic (EM) counterparts show that double neutron star mergers could lead to rich and bright EM emissions. Recent numerical simulations suggest that neutron star and neutron star/black hole (NS-NS/BH) mergers would leave behind a central remnant surrounded by a mildly isotropic ejecta. The central remnant could launch a collimated jet and when the jet propagating through the ejecta, a mildly relativistic cocoon would be formed and the interaction between the cocoon and the ambient medium would accelerate electrons via external shock in a wide angle. So that the merger-nova photons (i.e., thermal emission from the ejecta) would be scattered into higher frequency via inverse compton (IC) process when they propagating through the cocoon shocked region. We find that the IC scattered component peaks at X-ray band and it will reach its peak luminosity in order of days (simultaneously with the merger-nova emission). With current X-ray detectors, such a late X-ray component could be detected out to 200 Mpc, depending on the merger remnant properties. It could serve as an important electromagnetic counterpart of gravitational wave signals from NS-NS/BH mergers. Nevertheless, simultaneous detection of such a late X-ray signal and the merger-nova signal could shed light on the cocoon properties and the concrete structure of the jet.

Recently, the LIGO-Virgo collaboration reported their first detection of gravitational wave (GW) signals from a low mass compact binary merger GW170817, which is most likely due to a double neutron star (NS) merger. With the GW signals only, the chirp mass of the binary is precisely constrained to $1.188^{+0.004}_{-0.002}~\rm{M_{\odot}}$, but the mass ratio is loosely constrained in the range $0.4-1$, so that a very rough estimation of the individual NS masses ($0.86~{\rm M_{\odot}}<M_1<1.36~\rm{M_{\odot}}$ and $1.36~{\rm M_{\odot}}<M_2<2.26~\rm{M_{\odot}}$) was obtained. Here we propose that if one can constrain the dynamical ejecta mass through performing kilonova modeling of the optical/IR data, by utilizing an empirical relation between the dynamical ejecta mass and the mass ratio of NS binaries, one may place a more stringent constraint on the mass ratio of the system. For instance, considering that the red "kilonova" component is powered by the dynamical ejecta, we reach a tight constraint on the mass ratio in the range of $0.46-0.59$. Alternatively, if the blue "kilonova" component is powered by the dynamical ejecta, the mass ratio would be constrained in the range of $0.53-0.67$. Overall, such a multi-messenger approach could narrow down the mass ratio of GW170817 system to the range of $0.46-0.67$, which gives a more precise estimation of the individual NS mass than pure GW signal analysis, i.e. $0.90~{\rm M_{\odot}}<M_1<1.16~{\rm M_{\odot}}$ and $1.61~{\rm M_{\odot}}<M_2<2.11~{\rm M_{\odot}}$.

Limited by the sensitivities of the current gravitational wave (GW) detectors, the central remnant of the binary neutron star (NS) merger associated with GW170817 remains an open question. Considering the relatively large total mass, it is generally proposed that the merger of GW170817 would lead to a shortly lived hypermassive NS or directly produce a black hole (BH). There is no clear evidence to support or rule out a long-lived NS as the merger remnant. Here we utilize the GW and electromagnetic (EM) signals to comprehensively investigate the parameter space that allows a long-lived NS to survive as the merger remnant of GW170817. We find that for some stiff equations of state, the merger of GW170817 could, in principle, lead to a massive NS, which has a millisecond spin period. The post-merger GW signal could hardly constrain the ellipticity of the NS. If the ellipticity reaches $10^{-3}$, in order to be compatible with the multi-band EM observations, the dipole magnetic field of the NS ($B_p$) is constrained to the magnetar level of $\sim10^{14}$ G. If the ellipticity is smaller than $10^{-4}$, $B_p$ is constrained to the level of $\sim10^{10}-10^{12}\,$G. These conclusions weakly depend on the adoption of equations of state.