Recently, many gravitational wave events from compact binary mergers have been detected by LIGO. Determining the final mass and spin of the remnant black holes (RBHs) is a fundamental issue and is also important in astrophysics. In this paper, unified models for predicting the final mass and spin of the RBHs from compact binary mergers is proposed. The models achieve a good accuracy within the parameter range of interest. In addition, the rotational energy of the RBHs is also studied which is relevant to the electromagnetic counterparts of the mergers. It is found the distribution of the rotational energy of the RBHs from different types of mergers of compact binary has its own characteristics, which might help identify the electromagnetic counterparts associated with the mergers.

Owing to the detection of an extremely bright fast radio burst (FRB) 200428 associated with a hard X-ray counterpart from the magnetar soft gamma-ray repeater (SGR) 1935+2154, the distance of SGR 1935+2154 potentially hosted in the supernova remnant (SNR) G57.2+0.8 can be revisited. Under the assumption that the SGR and the SNR are physically related, in this Letter, by investigating the dispersion measure (DM) of the FRB contributed by the foreground medium of our Galaxy and the local environments and combining with other observational constraints, we find that the distance of SGR 1935+2154 turns out to be $9.0\pm2.5\,$kpc and the SNR radius falls into $10$ to $18\,$pc since the local DM contribution is as low as $0-18\,$pc cm$^{-3}$. These results are basically consistent with the previous studies. In addition, an estimate for the Faraday rotation measure of the SGR and SNR is also carried out.

A Kepler solver is an analytical method used to solve a two-body problem. In this paper, we propose a new correction method by slightly modifying the Kepler solver. The only change to the analytical solutions is that the obtainment of the eccentric anomaly relies on the true anomaly that is associated to a unit radial vector calculated by an integrator. This scheme rigorously conserves all integrals and orbital elements except the mean longitude. However, the Kepler energy, angular momentum vector and Laplace-Runge-Lenz vector for perturbed Kepler problems are slowly-varying quantities. However, their integral invariant relations give the quantities high-precision values that directly govern five slowly-varying orbital elements. These elements combined with the eccentric anomaly determine the desired numerical solutions. The newly proposed method can considerably reduce various errors for a post-Newtonian two-body problem compared with an uncorrected integrator, making it suitable for a dissipative two-body problem. Spurious secular changes of some elements or quasi-integrals in the outer solar system may be caused by short integration times of the fourth-order Runge-Kutta algorithm. However, they can be eliminated in a long integration time of 10^8 years by the proposed method, similar to Wisdom-Holman second-order symplectic integrator. The proposed method has an advantage over the symplectic algorithm in the accuracy but gives a larger slope to the phase error growth.

For a binary composed of a spinning black hole (BH) (with mass $\gtrsim 7M_\odot$) and a strongly magnetized neutron star (NS) (with surface magnetic field strength $B_{\rm S,NS}\gtrsim10^{12}$\,G and mass $\sim 1.4M_\odot$), the NS as a whole will possibly eventually plunge into the BH. During the inspiral phase, the spinning BH could be charged to the Wald charge quantity $Q_{\rm W}$ until merger in an electro-vacuum approximation. During the merger, if the spinning charged BH creates its own magnetosphere due to an electric field strong enough for pair cascades to spark, the charged BH would transit from electro-vacuum to force-free cases and could discharge in a time $\gtrsim1~{\rm ms}$. As the force-free magnetosphere is full of a highly conducting plasma, the magnetic flux over the NS's caps would be retained outside the BH's event horizon under the frozen-in condition. Based on this scenario, we here investigate three possible energy-dissipation mechanisms that could produce electromagnetic (EM) counterparts in a time interval of the BH's discharge post a BH-NS merger-induced gravitational wave event: (1) magnetic reconnection at the BH's poles would occur, leading to a millisecond bright EM signal, (2) a magnetic shock in the zone of closed magnetic field lines due to the detachment and reconnection of the entire BH magnetic field would probably produce a bright radio emission, e.g., a fast radio burst, and (3) the Blandford-Znajek mechanism would extract the BH's rotational energy, giving rise to a millisecond-duration luminous high-energy burst. We also calculate the luminosities due to these mechanisms as a function of BH's spin for different values of $B_{\rm S,NS}$.

Fast radio bursts (FRBs) are intense radio transients whose physical origin remains unknown. Therefore, it is of crucial importance to use a model-independent method to obtain the energy function and cosmic formation rate directly from the observational data. Based on current samples from the Parkes and ASKAP telecsopes, we determine, for the first time, the energy function and formation rate of FRBs by using the Lynden-Bell $\rm C^{-}$ method. The energy function derived from the Parkes sample is a broken power law, however it is a simple power law for the ASKAP sample. For Parkes sample, we derive the formation rate which is roughly consistent with the star formation rate up to $z\sim1.7$, with a local formation rate of $\dot \rho (0) \simeq (3.2\pm 0.3)\; \times {10^4}\;{\rm{Gp}}{{\rm{c}}^{ - 3}}{\rm{y}}{{\rm{r}}^{ - 1}}$ above a detection threshold of $2\,{\rm{Jyms}}$. For ASKAP sample, we find that the formation rate evolves much faster than the star formation rate up to $z\sim0.7$, namely $\dot \rho (z) \propto {(1 + z)^{6.9 \pm 1.9}}$, with a local formation rate of $\dot \rho (0) \simeq (4.6\pm 0.8)\; \times {10^3}\;{\rm{Gp}}{{\rm{c}}^{ - 3}}{\rm{y}}{{\rm{r}}^{ - 1}}$ above a detection threshold of $51\,{\rm{Jyms}}$. This might be a important clue for the physical origin of FRBs.

In this paper we propose the model that the coalescence of primordial black holes (PBHs) binaries with equal mass $M \sim 10^{28}$g can emit luminous gigahertz (GHz) radio transient, which may be candidate sources for the observed fast radio bursts (FRBs), if at least one black hole holds appropriate amount of net electric charge $Q$. Using a dimensionless quantity for the charge $q = Q/\sqrt{G}M$, our analyses infer that $q\sim O(10^{-4.5})$ can explain the FRBs with released energy of order $O(10^{40}) {\rm ergs}$. With the current sample of FRBs and assuming a distribution of charge $\phi(q)$ for all PBHs, we can deduce that its form is proportional to $q^{-3.0\pm0.1}$ for $q\geq 7.2\times10^{-5}$ if PBHs are sources of the observed FRBs. Furthermore, with the proposed hypothetical scenario and by estimating the local event rate of FRBs $\sim 2.6 \times 10^3 {\rm Gpc}^{-3} {\rm yr}^{-1}$, one derives a lower bound for the fraction of PBHs (at the mass of $10^{28}$g) against that of matter $f_{\rm PBH}(10^{28}{\rm g})$ $\gtrsim 10^{-5}$ needed to explain the rate. With this inspiring estimate, we expect that future observations of FRBs can help to falsify their physical origins from the PBH binaries coalescences. In the future, the gravitational waves produced by mergers of small black holes can be detected by high frequency gravitational wave detectors. We believe that this work would be a useful addition to the current literature on multimessenger astronomy and cosmology.

HI absorption studies yield information on both AGN feeding and feedback processes. This AGN activity interacts with the neutral gas in compact radio sources, which are believed to represent the young or recently re-triggered AGN population. We present the results of a survey for HI absorption in a sample of 66 compact radio sources at 0.040 < z < 0.096 with the Australia Telescope Compact Array. In total, we obtained seven detections, five of which are new, with a large range of peak optical depths (3% to 87%). Of the detections, 71% exhibit asymmetric, broad ($\Delta{v}_{\mathrm{FWHM}}$ > 100 km s$^{-1}$) features, indicative of disturbed gas kinematics. Such broad, shallow and offset features are also found within low-excitation radio galaxies which is attributed to disturbed circumnuclear gas, consistent with early-type galaxies typically devoid of a gas-rich disk. Comparing mid-infrared colours of our galaxies with HI detections indicates that narrow and deep absorption features are preferentially found in late-type and high-excitation radio galaxies in our sample. These features are attributed to gas in galactic disks. By combining XMM-Newton archival data with 21-cm data, we find support that absorbed X-ray sources may be good tracers of HI content within the host galaxy. This sample extends previous HI surveys in compact radio galaxies to lower radio luminosities and provides a basis for future work exploring the higher redshift universe.

We report our very early optical observations of GRB 110530A and investigate its jet properties together with its X-ray afterglow data. A peculiar broad onset bump followed by a plateau is observed in its early R band afterglow lightcurve. The optical data in the other bands and the X-ray data are well consistent with the temporal feature of the R band lightcurve. Our joint spectral fits of the optical and X-ray data show that they are in the same regime, with a photon index of $\sim 1.70$. The optical and X-ray afterglow lightcurves are well fitted with the standard external shock model by considering a delayed energy injection component. Based on our modeling results, we find that the radiative efficiency of the GRB jet is $\sim 1\%$ and the magnetization parameter of the afterglow jet is $<0.04$ with the derived extremely low $\epsilon_B$ (the fraction of shock energy to magnetic field) of $(1.64\pm 0.25)\times 10^{-6}$. These results indicate that the jet may be matter dominated. Discussion on delayed energy injection from accretion of late fall-back material of its pre-supernova star is also presented.

The cosmic evolution of gamma-ray burst (GRB) luminosity is essential for revealing the GRB physics and for using GRBs as cosmological probes. We investigate the luminosity evolution of long GRBs with a large sample of 258 {\em Swift}/BAT GRBs. Parameterized the peak luminosity of individual GRBs evolves as $L_{\rm p}\propto{\rm }(1+z)^{k}$, we get $k=1.49\pm0.19$ using the non-parametric $\tau$ statistics method without considering observational biases of GRB trigger and redshift measurement. By modeling these biases with the observed peak flux and characterizing the peak luminosity function of long GRBs as a smoothly broken power-law with a break that evolves as $L_{\rm b}\propto (1+z)^{k_{\rm b}}$, we obtain $k_{\rm b}=1.14^{+0.99}_{-0.47}$ through simulations based on assumption that the long GRB rate follows the star formation rate (SFR) incorporating with cosmic metallicity history. The derived $k$ and $k_b$ values are systematically smaller than that reported in previous papers. By removing the observational biases of the GRB trigger and redshift measurement based on our simulation analysis, we generate mock {\em complete} samples of 258 and 1000 GRBs to examine how these biases affects on the $\tau$ statistics method. We get $k=0.94\pm 0.14$ and $k=0.80\pm 0.09$ for the two samples, indicating that these observational biases may lead to overestimate the $k$ value. With the large uncertain of $k_b$ derived from our simulation analysis, one even cannot convincingly argue a robust evolution feature of the GRB luminosity.

The external forward shock (EFS) models have been the standard paradigm to interpret the broad-band afterglow data of gamma-ray bursts (GRBs). One prediction of the models is that some afterglow temporal breaks at different energy bands should be achromatic. Observations in the Swift era have revealed chromatic afterglow behaviors at least in some GRBs, casting doubts on the EFS origin of GRB afterglows. In this paper, we perform a systematic study to address the question: how bad/good are the external forward shock models? Our sample includes 85 GRBs well-monitored X-ray and optical lightcurves. Based on how well the data abide by the EFS models, we categorize them as: Gold sample: (Grade I and II) include 45/85 GRBs. They show evidence of, or are consistent with having, an achromatic break. The temporal/spectral behaviors in each afterglow segment are consistent with the predictions (closure relations) of the EFS models. Silver sample: (Grade III and IV) include 37/85 GRBs. They are also consistent with having an achromatic break, even though one or more afterglow segments do not comply with the closure relations. Bad sample: (Grade V), 3/85 shows direct evidence of chromatic behaviors, suggesting that the EFS models are inconsistent with the data. These are included in the Bad sample. We further perform statistical analyses of various observational properties ($\alpha$, $\beta$, $t_b$ and model parameters (energy injection index q, p, $\theta_j$, $\eta_\gamma$, etc) of the GRBs in the Gold Sample, and derive constraints on the magnetization parameter $\epsilon_B$ in the EFS. Overall, we conclude that the simplest EFS models can account for the multi-wavelength afterglow data of at least half of the GRBs. When more advanced modeling (e.g., long-lasting reverse shock, structured jets) is invoked, up to $>90 \%$ of the afterglows may be interpreted within the framework of the ESF models.

We propose a novel beam model for radio pulsars based on the scenario that the broadband and coherent emission from secondary relativistic particles, as they move along a flux tube in a dipolar magnetic field, forms a radially extended sub-beam with unique properties. The whole radio beam may consist of several sub-beams, forming a fan-shaped pattern. When only one or a few flux tubes are active, the fan beam becomes very patchy. This model differs essentially from the conal beam models in the respects of beam structure and predictions on the relationship between pulse width and impact angle $\beta$ (the angle between line of sight and magnetic pole) and the relationship between emission intensity and beam angular radius. The evidence for this model comes from the observed patchy beams of precessional binary pulsars and three statistical relationships found for a sample of 64 pulsars, of which $\beta$ were mostly constrained by fitting polarization position angle data with the Rotation Vector Model. With appropriate assumptions, the fan beam model can reproduce the relationship between 10\% peak pulse width and $|\beta|$, the anticorrelation between the emission intensity and $|\beta|$, and the upper boundary line in the scatter plot of $|\beta|$ versus pulsar distance. An extremely patchy beam model with the assumption of narrowband emission from one or a few flux tubes is studied and found unlikely to be a general model. The implications of the fan beam model to the studies on radio and gamma-ray pulsar populations and radio polarization are discussed.