Search

We propose a new model of Fast Radio Bursts (FRBs) based on stellar mass black hole-massive star binaries. We argue that the inhomogeneity of the circumstellar materials or/and the time varying wind activities of the stellar companion will cause the black hole to accrete at a transient super-Eddington rate. The collision among the clumpy ejecta in the resulted jet could trigger plasma instability. As a result, the plasma in the jet will emit coherent curvature radiation. When the jet cone aims toward the observer, the apparent luminosity can be $10^{41}-10^{42}$\,erg/s. The duration of the resulted flare is $\sim$ millisecond. The high event rate of the observed non-repeating FRBs can be explained. A similar scenario in the vicinity of a supermassive black hole can be used to explain the interval distribution of the repeating source FRB121102 qualitatively., Comment: Accepted for publication in MNRAS

Among the four black hole binary merger events detected by LIGO, six progenitor black holes have masses greater than 20\,$M_\odot$. The existence of such massive BHs calls for extreme metal-poor stars as the progenitors. An alternative possibility that a pair of stellar mass black holes each with mass $\sim7\,M_\odot$ increases to $>20\,M_\odot$ via accretion from a disk surrounding a super massive black hole in an active galactic nucleus is considered. The growth of mass of the binary and the transfer of orbital angular momentum to the disk accelerates the merger. Based on the recent numerical work of Tang et al. (2017), it is found that, in the disk of a low mass AGN with mass $\sim10^6\,M_\odot$ and Eddington ratio $>0.01$, the mass of an individual BH in the binary can grow to $>20\,M_\odot$ before coalescence provided that accretion takes place at a rate more than 10 times the Eddington value. The mechanism predicts a new class of gravitational wave sources involving the merger of two extreme Kerr black holes associated with active galactic nuclei and a possible electromagnetic wave counterpart., Comment: Accepted for publication in ApJL

The $\gamma$-ray pulsar binary system PSR B1259-63 flares in GeV after each periastron. The origin of those flares is still under debate. Recently, [2017ApJ...844..114Y] proposed a mechanism that might explain the GeV flares. In that model, a transient accretion disc is expected to be formed from the matter which was gravity-captured by the neutron star from the main sequence companion's circumstellar disc. The transient accretion disc exerts a spin down torque on the neutron star (propeller effect), which might be traceable via pulsar timing observation of PSR B1259-63. In this paper, we phenomenologically consider the propeller effect with a parameter $\chi$, which describes the coupling between the disc matter and the neutron star. Comparing the expected timing residuals against the recent observation in [2014MNRAS.437.3255S], we conclude that the angular momentum transfer is very weak (with the coupling parameter $\chi\le10^{-4}$)., Comment: 5 Pages, 7 Figures,accepted for publication in MNRAS

The velocity and density distribution of $e^\pm$ in the pulsar wind are crucial distinction among magnetosphere models, and contains key parameters determining the high energy emission of pulsar binaries. In this work, a direct method is proposed, which might probe the properties of the wind from one pulsar in a double-pulsar binary. When the radio signals from the first-formed pulsar travel through the relativistic $e^\pm$ flow in the pulsar wind from the younger companion, the components of different radio frequencies will be dispersed. It will introduce an additional frequency-dependent time-of-arrival delay of pulses, which is function of the orbital phase. In this paper, we formulate the above-mentioned dispersive delay with the properties of the pulsar wind. As examples, we apply the formula to the double pulsar system PSR J0737-3039A/B and the pulsar-neutron star binary PSR B1913+16. For PSR J0737-3039A/B, the time delay in 300\,MHz is $\lesssim10\mu$s near the superior-conjunction, under the optimal pulsar wind parameters, which is $\sim$ half of the current timing accuracy. For PSR B1913+16, with the assumption that the neutron star companion has a typical spin down luminosity of $10^{33}$\,ergs/s, the time delay is as large as $10\sim20\mu$s in 300\,MHz. The best timing precision of this pulsar is $\sim5\mu$s in 1400\,MHz. Therefore, it is possible that we can find this signal in archival data. Otherwise, we can set an upper-limit on the spin down luminosity. Similar analysis can be apply to other eleven known pulsar-neutron star binaries, Comment: 6 pages, 6 figures, accepted for publication in MNRAS

The binary system HESS J0632+057 was recently detected by {Fermi} to possess orbital modulated GeV emission. In this paper, we study the possibility that the compact companion of HESS J0632+057 is a pulsar. Under such a presumption, we focus on the high energy emission mechanism of this system, which is as follows. The pulsar companion travels through the circumstellar disc of the main sequence star twice in each orbit, when some of the matter is gravity-captured. The captured matter develops an accretion disc around the pulsar, and the soft photons from which are inverse Compton scattered by the pulsar wind as the GeV emission from the system. With proper choice of parameters, SED and light curve which are in accordance with observations can be produced. We predict that the light curve of GeV emission has two peaks, the larger one is at around 0.4 after the periastron (or 0.1 after the X-ray maximum), while the smaller one is between phases 0 and 0.1, with integrated flux one forth of the larger one., Comment: 7pages, 7 figures. Accepted for publication in MNRAS

Analogy to Cherenkov radiation, when a particle moves faster than the propagation velocity of gravitational wave in matter ($v>c_{\rm{g}}$), we expect gravitational wave-Cherenkov radiation (GWCR). In the situation that a photon travels across diffuse dark matters, the GWCR condition is always satisfied, photon will thence loss its energy all the path. This effect is long been ignored in the practice of astrophysics and cosmology, without justification with serious calculation. We study this effect for the first time, and shows that this energy loss time of the photon is far longer than the Hubble time, therefore justify the practice of ignoring this effect in astrophysics context., Comment: 3 pages; We add a citation to Caves (1980) to the original journal version, after kindly reminded

PSR B1259-63/LS2883 is a binary system composed of a pulsar and a Be star. The Be star has an equatorial circumstellar disk (CD). The {\it Fermi} satellite discovered unexpected gamma-ray flares around 30 days after the last two periastron passages. The origin of the flares remain puzzling. In this work, we explore the possibility that, the GeV flares are consequences of inverse Compton-scattering of soft photons by the pulsar wind. The soft photons are from an accretion disk around the pulsar, which is composed by the matter from CD captured by the pulsar's gravity at disk-crossing before the periastron. At the other disk-crossing after the periastron, the density of the CD is not high enough so that accretion is prevented by the pulsar wind shock. This model can reproduce the observed SEDs and light curves satisfactorily., Comment: 14 pages, 8 figures, 1 table. Accepted for publication in ApJ

We explain some phenomena existing generally in the timing residuals: amplitude and sign of the second derivative of a pulsar's spin-frequency ($\ddot\nu$), some sophisticated residual patterns, which also change with the time span of data segments. The sample is taken from Hobbs et al.\,(2010), in which the pulsar's spin-frequency and its first derivative have been subtracted from the timing solution fitting. We first classify the timing residual patterns into different types based on the sign of $\ddot\nu$. Then we use the magnetic field oscillation model developed in our group \citep{zhang12a} to fit successfully the different kinds of timing residuals with the Markov Chain Monte Carlo method. Finally, we simulate the spin evolution over 20 years for a pulsar with typical parameters and analyze the data with the conventional timing solution fitting. By choosing different segments of the simulated data, we find that most of the observed residual patterns can be reproduced successfully. This is the first time that the observed residual patterns are fitted by a model and reproduced by simulations with very few parameters. From the distribution of the different residual patterns in the $P-\dot P$ diagram, we argue that (1) a single magnetic field oscillation mode exists commonly in all pulsars throughout their lifetimes; (2) there may be a transition period over the lifetimes of pulsars, in which multiple magnetic field oscillation modes exist., Comment: 21 pages, 22 figures, accepted for publication in MNRAS

XTE J1810-189 underwent an outburst in 2008, and was observed over $\sim 100$ d by RXTE. Performing a time-resolved spectral analysis on the photospheric radius expansion burst detected on 2008 May 4, we obtain the source distance in the range of 3.5--8.7 kpc for the first time. During its outburst, XTE J1810-189 did not enter into the high/soft state, and both the soft and hard colours decreased with decreasing flux. The fractional rms remained at high values ($\sim 30$ per cent). The RXTE/PCA spectra for 3-25 keV can be described by an absorbed power-law component with an additional Gaussian component, and the derived photon index $\Gamma$ increased from $1.84\pm0.01$ to $2.25\pm0.04$ when the unabsorbed X-ray luminosity in 3-25 keV dropped from $4\times10^{36}$ ergs s$^{-1}$ to $6\times10^{35}$ ergs s$^{-1}$. The relatively high flux, dense observations and broadband spectra allow us to provide strong evidence that the softening behaviour detected in the outburst of XTE J1810-189 originates from the evolution of non-thermal component rather than the thermal component (i.e. neutron star surface emission)., Comment: 7 pages, 5 figures, 3 tables, published in MNRAS

We study the possibility that the long term red timing-noise in pulsars originates from the evolution of the magnetic inclination angle $\chi$. The braking torque under consideration is a combination of the dipole radiation and the current loss. We find that the evolution of $\chi$ can give rise to extra cubic and fourth-order polynomial terms in the timing residuals. These two terms are determined by the efficiency of the dipole radiation, the relative electric-current density in the pulsar tube and $\chi$. The following observation facts can be explained with this model: a) young pulsars have positive $\ddot{\nu}$; b) old pulsars can have both positive and negative $\ddot{\nu}$; c) the absolute values of $\ddot{\nu}$ are proportional to $-\dot{\nu}$; d) the absolute values of the braking indices are proportional to the characteristic ages of pulsars. If the evolution of $\chi$ is purely due to rotation kinematics, then it can not explain the pulsars with braking index less than 3, and thus the intrinsic change of the magnetic field is needed in this case. Comparing the model with observations, we conclude that the drift direction of $\chi$ might oscillate many times during the lifetime of a pulsar. The evolution of $\chi$ is not sufficient to explain the rotation behavior of the Crab pulsar, because the observed $\chi$ and $\dot{\chi}$ are inconsistent with the values indicated from the timing residuals using this model., Comment: 5 pages, 1 figure. Accepted for publication in MNRAS

Unlike the random radial orientation distribution of field elliptical galaxies, galaxies in a cluster are expected to point preferentially towards the center of the cluster, as a result of the cluster's tidal force on its member galaxies. In this work an analytic model is formulated to simulate this effect. The deformation time scale of a galaxy in a cluster is usually much shorter than the time scale of change of the tidal force; the dynamical process of the tidal interaction within the galaxy can thus be ignored. An equilibrium shape of a galaxy is then assumed to be the surface of equipotential, which is the sum of the self-gravitational potential of the galaxy and the tidal potential of the cluster at this location. We use a Monte-Carlo method to calculate the radial orientation distribution of these galaxies, by assuming the NFW mass profile of the cluster and the initial ellipticity of field galaxies. The radial angles show a single peak distribution centered at zero. The Monte-Carlo simulations also show that a shift of the reference center from the real cluster center weakens the anisotropy of the radial angle distribution. Therefore, the expected radial alignment cannot be revealed if the distribution of spatial position angle is used instead of that of radial angle. The observed radial orientations of elliptical galaxies in cluster Abell~2744 are consistent with the simulated distribution., Comment: 8 pages, 6 figures, 2 tables. MNRAS in press