1 Introduction
Fast Radio Bursts (FRBs) are mysterious, bright astronomical millisecond-duration radio pulses [
1–
3], occasionally discovered in pulsar searches. Up to date, more than one hundred FRBs have been reported [
4]
†. They are expected to be of cosmological origin due to high dispersion measures (DM, 10
2–10
3 pc·cm
−3) in excess of the DM contribution from the Milky Way. And it is further confirmed by the localization of the host galaxies [
5–
8].
Among the moderately observed large sample of FRBs, most of them are one-off events, and only twenty of them sporadically show repeating bursts [
9–
11]. Moreover, FRB 180916.J0158+65 detected by CHIME [
12] exhibits a ~16 day period acitvity with unknown mechanism. Most recently, a bright FRB 200428 has been reported to be spatially coincident with the galactic Soft Gamma-ray Repeater (SGR) 1935+2154 [
13–
16], also associated with a hard X-ray burst [
14,
15,
17,
18]. The repetition of FRB 200428 is very rare, which may be different from other extragalactic repeating FRBs [
13].
The progenitors of FRBs are still under heated debates, and many theoretical models have been proposed (see [
19] and references therein
†), despite some of them have confronted serious challenges. Repeating FRBs are the best candidates to explore physical nature due to its repetition. FRB 121102 is the first observed repeating FRB and has been detected in the frequency range from 600 MHz [
20] to 8 GHz [
21]. For repeating FRBs, some models involving catastrophic processes have been ruled out, e.g., binary neutron star (NS) mergers [
22,
23], and binary white dwarf mergers [
24]. Many non-catastrophic models have been put forward, such as the flaring magnetars model [
25–
28], accretion process by a neutron star in a binary system [
29], the interaction of NS with an asteroid belt [
30], the interaction between the NS magnetosphere and cosmic winds [
31], and the giant pulses (GPs) from pulsars [
32,
33]. The central engine of FRB 121102 is questionable, which can be constrained from the persistent radio nebula of FRB 121102 for the magnetar-type model [
34–
38].
GPs were first discovered in the Crab pulsar (PSR B0531+21), and then were found in several energetic pulsars (see Table 1 of [
39] and reference therein). Different from the regular pulses, the typical duration of GPs is very short, with no clear periodicity, ranging from nanoseconds to microseconds. The flux density could be 2 or 4 orders of magnitude higher than the average pulse-integrated flux of regular energetic pulses (e.g., [
40,
41]). Among the pulsars emitting the GPs, the largest value of the magnetic field at the light cylinder
Blc are the Crab pulsar and PSR B1937+21 [
39]. The Crab pulsar is a remarkable source emitting numerous GPs and simultaneously has been detected across a broadband range of radio frequencies from 20 MHz [
42] to 46 GHz [
43]. The distinctive feature of GPs is the energy of GPs obeys a power-law distribution, whereas the regular ones follow a Gaussian or log-normal distribution (e.g., [
40,
44]), indicating that the emission mechanism of GPs may be different from that of regular pulses. Meanwhile, GPs are also crucial to understand the radiation mechanism of pulsars.
Repeating FRBs and GPs from pulsars share similarities in many observational aspects, such as brightness temperature, duration time, sub-bursts, polarization, and frequency drift. Moreover, the emission mechanism of both of them is coherent radiation, supporting that they may have a similar physical origin. Actually, in some theoretical models, it has been suggested that FRBs can be treated as GPs from young rapidly rotating pulsars [
32,
33].
The statistical properties, such as the energy distributions, of repeating FRBs and GPs might help us to understand FRBs progenitors. In this paper, we focus on the statistical properties of the repeating FRB 121102 and GPs from the Crab pulsar. The paper is organized as follows. In Section 2, we introduce the sample selection. The statistical analysis of the isotropic-equvivalent energy E is shown in Section 3. Discussion and conclusions are given in Section 4 and Section 5. We assume a flat ΛCDM universe with Ωm=0.27, and H0=70 km·s−1·Mpc−1 throughout this work.
2 Sample selection
The Crab Nebula pulsar is well known for its GPs, with detected GPs numerous enough for detailed studies (e.g., [
40,
41,
45]). Here, GPs refer to GPs generated by the main pulse (MP) and interpulse (IP) phase.
In our analysis, we adopt a complete sample of 1153 bright GPs (with fluence ≥130 Jy·ms to avoid the incomplete effect near the detection threshold) from the Crab pulsar [
45]. These GPs were detected by the 15-m telescope in Pune, India (National Centre for Radio Astrophysics) in a 260-hours observation from February to April 2019. The observation was conducted between 1280 MHz and 1380 MHz, with 65 MHz bandwidth usable. Each GP was detected with a signal-to-noise ratio (S/N) of ≥10.
Among 20 repeating FRBs have been reported, only two of repeating FRBs have host galaxy identified and redshift measured (FRB 121102 and FRB 180916.J0158+65; [
5,
12]). FRB 121102 is the first and currently the best representative of repeating bursts. More than one hundred bursts from FRB 121102 have been detected by several different surveys, including Green Bank Telescope (GBT) in C-band (4–8 GHz; [
46]) and Arecibo telescope at 1.4 GHz [
47]. The repeating FRB 121102 has a very large and variable Faraday rotation measure (RM ~10
5 rad·m
−2) which means the burst inhabits in an extreme magneto-ionic environment [
3].
To compare with the Crab GPs sample under the same condition, we use the sample of FRB 121102 detected at almost the same frequency [
47]. The observing frequency of Arecibo in the source frame is 1.4(1+
z) GHz. Due to the redshift 0.193 for FRB 121102 [
48], the frequency in the source rest frame approximatively is ~1.4 GHz. In this sample, there are 41 bursts from FRB 121102 observed by the 305 m Arecibo telescope. The data is from two observations taken on 2016 September 13/09:47:07 and 14/09:50:12, and lasting 5967 s and 5545 s. Each observation detected 18 and 23 bursts with S/N ≥10, respectively. Moreover, we also investigate the energy distribution of the repeating FRB 121102 detected by GBT telescope at C-band [
46]. There are 93 bursts in total in 5 hours with a convolutional neural network technique.
3 The cumulative energy E distributions
Energy is a crucial parameter to identify the progenitor and to explore the radiation mechanism. We perform an analysis on the energy distribution of the repeating FRB 121102 and GPs to explore the progenitor of FRB 121102, despite there are many works devoted to the statistical properties of the repeating FRBs [
46,
49–
52].
The isotropic-equvivalent energy E in the source rest frame within the observing bandwidth for the selected sample can be calculated by
where
dL is the luminosity distance of the source, with 2.0 kpc for the Crab pulsar,
is the fluence density of each burst, and
is the corresponding bandwidth of every burst. In this paper, the bandwidth for FRB 121102 is derived from the burst frequency edges (
fhigh and
flow) from the table given in Ref. [
47].
Here we adopt a simple power-law model, to fit the cumulative energy distributions of bright GPs (fluence ≥ 130 Jy·ms) from the Crab pulsar. While for the repeating FRB 121102, we fit it by the smooth broken power-law function, i.e.,
where A is the ampltitude, are slopes of the two segments, Eb is the break energy, ω is the sharpness (or smoothness) of the break of energy distribution, and it usually can be set as a constant (ω = 3).
We can now obtain the goodness-of-fit χd.o.f for the difference between the observed cumulative distribution Ncum,obs and the theoretical distribution function Ncum for the cumulative distribution (with npar = 4 parameters, E0, α1, Eb, α2)
where
Ntot is the number of bursts dectected by the same instrument at the same frequency for the same source. Note that the uncertainty of
is taken as
. The best-fitting parameters and the uncertainties can be obtained by minimizing
χd.o.f via the python package emcee (a Bayesian MCMC method by sampling the affine-invariant for Markov Chain Monte Carlo chains, [
53]).
When fitting variables via a MCMC approach, we assume the various parameters in the fitting function are uniform distribution in the prior probability distribution. Here are the range of parameters for every sample: 1 < α< 5 for GPs; 50 < A < 200, -2.5 <α1 < -0.5, 36 < lgEb <38 and 2 < α2 < 5 for Arecibo sample; 50 < A< 200, -3< α1 < -0.01, 36 < lgEb <38.5 and 1< α2 < 3 for GBT sample.
GPs are generally considered as the power-law tail of the regular pulses for the pulsar. In Fig. 1(a), the yellow-green points represent the GPs from the Crab pulsar, which is found to be well fitted by the power-law model. The rollover does not appear at the high-energy end, while it does show for the flux density distribution as stated in some previous studies [
40,
44]. The best-fit slope is 2.85±0.10 (
χd.o.f = 0.45), which is well consistent with that of the fluence distribution (2.8±0.3) reported in [
45]. In Fig. 1, we present the cumulative energy (
E) distribution for 41 Arecibo bursts from FRB 121102 (blue spots). Fitted by the broken power-law function given in Eq. (2), the best fitting results for indices of the lower and higher energy segment are
and
(
χd.o.f = 0.36), respectively, while the break energy
Eb is
erg (blue dashed line). In Ref. [
47], the threshold energy was set at
erg considering the sensitivity limit and the observed turnover. Above the
Eth, the best fit is
, which is generally consistent with our fitting result within the uncertainty. All the uncertainties quote a 68% confidence interval for the best-fitted parameters.
To investigate what causes the different power indices of the energy distribution of FRB 121102 detected by the Arecibo telescope at L-band and that by GBT telescope at C-band and the meaning of the break energy Eb in Fig. 1, we also fit the energy distribution for 93 bursts from FRB 121102 by GBT telescope with broken power-law function. As shown in Fig. 2, the best fits for the lower and higher energy are and (χd.o.f = 0.63), respectively, while the break energy is erg.
Interestingly, with the same observed frequency (~1.4 GHz), FRB 121102 and Crab GPs share the power-law form of similar indices (~2.8) (Fig. 1). At a different observed frequency (C-band, 4–8 GHz), the energy distribution of FRB 121102 detected by the GBT telescope is flatter (Fig. 2) than that of the Arecibo telescope with more energetic bursts.
There is a lack of bursts below near the break energy ~1.60×10
37 erg (the black dashed line) due to the instrument sensitivity and observational limit. So maybe the incomplete sampling below the threshold is the major reason that evokes the different power-law indices of energy distribution for FRB 121102 in Fig. 1 and Fig. 2. The power-law index of energy distribution with repeating FRB 121102 may depend on the energy band, the burst is more energetic at a higher observational frequency (see Fig. 1 and Fig. 2), which can be further tested by detecting more bursts at different energy bands. Note that the energy distribution of GPs might also depend on the energy band [
40,
44], and can have the index ~1.8 for some pulsars (such as PSR 0631+1036; [
54]).
The power-law index of energy distribution and the statistical characteristics of the other observable quantities of FRB121102 are claimed to be similar to those of soft gamma repeater bursts from magnetars [
55]. In addition, there are several statistical analyses [
50,
55–
60] dedicated to the self-organizing critical (SOC) study, especially some researches [
50,
55,
60,
61] that indicate that the repeating FRB 121102 is a SOC system [
62]. Thus, the repeating FRB 121102, is likely to be an SOC system due to some kind of instability driven above the instability threshold value. If so, the break energy in Fig. 1 should be the threshold energy, which determines the trigger mechanism of the system. It is an important quantity that affects the slope of the power-law in energy distribution, which is also discussed in Ref. [
60].
Interestingly, the power-law slope of FRB 121102 omitting bursts that fall bellow the break energy is generally consistent with that of the GPs from the Crab as shown in Fig. 1, which implies they may share a similar origin. If FRB 121102 is extragalactic GPs, which energies can power it? We will discuss this in the next section below.
4 Discussion
As mentioned in the introduction, neutron stars are the most popular FRB central engine. Especially the super GPs model involving pulsars [
32,
33,
63] and the giant flare model involving magnetars [
26–
28,
64–
66] have been widely discussed. Specifically for FRB 121102, it was found that its host is a dwarf galaxy with a high star formation rate [
5,
48], which implies that the central engine of FRB 121102 may locate in the remnant of the death of a massive star. Moreover, the polarization observations give that the Faraday rotation measure of FRB 121102 is surprisingly high ~10
5 rad·m
−2, indicating that FRB 121102 should be surrounded by a strongly magnetized environment [
3]. These observational evidences support that the central engine of FRB 121102 may be a magnetized NS. In the following, we explore the physical origin of the repeating FRB 121102 in the framework of the super GPs model.
As shown in Fig. 1(a), the typical energy of the FRB 121102 and GPs are 10
37–10
38 erg and 10
29–10
30 erg. On the other hand, the typical time scale of the FRB 121102 and GPs are ms and μs, hence the typical luminosity are 10
40–10
41 erg/s and 10
35–10
36 erg/s, respectively. Following [
67], based on the typical luminosity of FRB 121102 and that of GPs, we have
Since the GPs are powered by the spin down of the pulsar Lsd, the constraint on the magnetic field strength BFRB of the FRB source is
By scaling the power of the Crab GPs to the level of FRB 121102, i.e.,
, where
BCrab~4×10
12 G and
PCrab~33 ms are the magnetic field strength and the period of the Crab pulsar, respectively [
68]. It is seen that the millisecond pulsars may produce GPs with luminosity comparable to FRB 121102. On the other hand, the active time scale is the spin-down time scale of the pulsar, that is
which can in principle explain FRB 121102 so far. Considering the condition of t
active > 7 yrs, we constraints that the period of the initial spin
P should be smaller than 3 ms. In addition, if the radiation efficiency of FRB
ηFRB is higher than that of GPs
ηGP, for example an order of magnitude higher, then we would have a longer active time scale
tactive~500 (
ηFRB/10
ηGP). Then the constraint to the period of the initial spin is
P≤8 ms in this case. For comparison, the constraint from the persistent radio nebula gives
P≤7 ms [
37].
Moreover, we note that both repeating FRBs and GPs have high brightness temperature (≥10
35 K), inferring that their radiation mechanism should be coherent emission [
69,
70]. GPs from Crab pulsar and FRBs 121102 share similar complex pulses morphology and a similar repetition rate in a single pulsar (~1–100 per hour, see Fig. 1). Furthermore, similar to FRBs [
71,
72], GPs have strong polarization signals, which is either left-handed or right-handed polarization [
73] between linear polarization and circular polarization. Besides, spectral structures in repeating FRBs resemble those seen in Crab GPs [
74]. Therefore, the similarities between the GPs and FRB 121102 possibly link to the same origin and mechanism of these two phenomena.
5 Conclusions
In this paper, we have made a comparison of the energy distributions between the repeating FRBs and the GPs from the Crab pulsar. We find that FRB 121102 and Crab GPs share similar power-law indices in the energy distribution at the same observed frequency ~1.4 GHz†. Furthermore, we explored the physical origin of the repeating FRB 121102 in the framework of the giant pulses model. We find that the millisecond pulsar is possible to produce super GPs with luminosity comparable to FRB 121102 in a reasonable active time scale. Therefore, we argue that the super GPs model is workable and favored to explain FRB 121102 despite that the magnetar burst model is more popular.