This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

Constraint on the relativistic motion of fast radio bursts based on the maximal electric field

Jun-Yi Shen,1 and Yuan-Chuan Zou1
1 Department of Astronomy, School of Physics, Huazhong University of Science and Technology, Wuhan 430074, China
E-mail: zouyc@hust.edu.cn
(Accepted XXX. Received YYY; in original form ZZZ)
Abstract

Fast radio bursts (FRBs) are millisecond radio signals from cosmological distances. As they propagate, FRBs can interact with ambient photons and initiate a quantum cascade that can limit the electric field strength. This paper examines whether some observed bright and brief FRBs may challenge this limit if the source is not relativistic. The size of a static FRB source is estimated as RctR\sim ct, where tt is the time scale of the FRB and cc denotes the speed of light. But for a relativistic source moving at the Lorentz factor Γ\Gamma, the size is R2Γ2ctR\sim 2\Gamma^{2}ct. Using an FRB catalog, we plot the luminosity-duration distribution. Most FRBs fall below the limit for a static source, but two events have higher luminosity and shorter duration. This suggests these bursts may originate from relativistic sources, although more data is needed to confirm this.

keywords:
transients: fast radio bursts – plasmas
pubyear: 2023pagerange: Constraint on the relativistic motion of fast radio bursts based on the maximal electric fieldConstraint on the relativistic motion of fast radio bursts based on the maximal electric field

1 Introduction

Fast radio bursts (FRBs) were first discovered by Lorimer in 2007 as millisecond radio transient sources (Lorimer et al., 2007). They are radio signals that originate from cosmological distances. The luminosity of individual FRB event ranges from 103810^{38} erg s1\mathrm{s^{-1}}(Bochenek et al., 2020) to 104610^{46} erg s1\mathrm{s^{-1}}(Ravi et al., 2019). The origin of FRBs is still mysterious. For non-repeating events, some catastrophic events have been proposed, such as mergers of binary white dwarfs (Kashiyama et al., 2013). Observations by the Canadian Hydrogen Intensity Mapping Experiment (CHIME) collaboration suggest that at least a portion of FRBs is associated with magnetars (Bochenek et al., 2020). FRBs could be generated by magnetar’s charge starvation region and current sheet region (Zhang, 2022b). More recent review on the physics of FRBs can be seen in Zhang (2022a).

The brightness temperature of FRB is exceptionally high (Zhang, 2022a): TB1036KT_{{\rm{B}}}\simeq 10^{36}\ \rm{K}. This indicates that the FRBs must be generated by a coherent process, while the exact radiation mechanism and origin are still unknown. Several models have been proposed to address this issue, as discussed in the review by Zhang (2020). Some models propose a gamma-ray burst-like mechanism (GRB-like), that is, FRBs are produced by relativistic shocks, as suggested by Lyubarsky (2014). Some models utilize a pulsar-like model, where coherent bunches of radio signals are emitted as FRBs, as described by Melikidze et al. (2000). One difference between the GRB-like and pulsar-like models is that the GRB-like models require a relativistic jet radiation region, while this is not necessary for the pulsar-like models. It would be meaningful to investigate whether the source of FRBs is relativistic or static.

Studying extreme events characterized by high luminosity and short duration can provide insights into the source’s nature. One common estimate of the size of an astrophysical source is RctR\sim ct, with tt being the time scale and cc being the speed of light. Extreme events have high luminosity and a small source size RR, resulting in a high energy density. However, there are limits to the energy density. One of them is the strength of the corresponding electric field of the electromagnetic wave E<Esme2c3/e=1.32×1016E<E_{{\rm{s}}}\equiv m_{{\rm{e}}}^{2}c^{3}/\hbar e=1.32\times 10^{16} V cm1\mathrm{cm^{-1}}, where mem_{{\rm{e}}} is the mass of the electron and ee is the charge of the electron and \hbar is the Planck constant. Once if EE reaches this limit, the energy of electromagnetic field converts into kinetic and rest-mass energy of electron-position pairs (Schwinger, 1951). Zhang & Wu (2022) found a stronger constraint on EE based on their simulations, which is E<3×1012Vcm1E<3\times 10^{12}\ {\rm V\,cm^{-1}}. With this limit they also placed constraints on the size of the radiation site. With the known duration of the FRBs, we may also constrain the Lorentz factor of the emitting region if the source is relativistic.

This paper is organized as follows. Section 2 introduces the quantum cascade process of the interaction between FRBs and high-energy photons, demonstrating the upper limit of EE. Section 3 highlights the distinctions between a relativistic source and a non-relativistic source. In section 4, we present the selected data to examine the validity of the relativistic source hypothesis. Finally, in section 5 we draw our conclusions based on the findings presented in this paper.

2 Interaction between FRB and photons

From the observations, we can estimate the electric field strength of the FRBs to be greater than 109Vcm110^{9}\ \rm{V\,cm^{-1}} (Zhang & Wu, 2022). In such a strong field, the interaction between high-energy photons and the field has been extensively studied (Kirk et al., 2009), including processes such as the Breit-Wheeler process and nonlinear inverse Compton scattering, which are quantum-electrodynamics (QED) processes. The Breit-Wheeler process involves the interaction of two high-energy photons, resulting in the creation of an electron-positron pair. The FRB GHz photons can be upscattered by high-energy electrons through inverse Compton scattering. If this scattering process generates a significant number of high-energy photons, providing sufficient energy for subsequent pair creation, the cascade starts. The presence of high-energy photons acts as a trigger, initiating the reaction. Then the Breit-Wheeler process occurs and creates positron-electron pairs. The electrons and positrons accelerate in electromagnetic field and radiate high energy photons. These high-energy photons continue to create pairs. The reaction rate is (in weak field) (Landau & Lifshitz, 1982):

p=33/2|e|3E29/2meexp(83η),p=\frac{3^{3/2}|e|^{3}E}{2^{9/2}m_{{\rm{e}}}}\mathrm{exp(-\frac{8}{3\eta}}), (1)

with

η=2|e|γωme3c5(𝐄+𝐯×𝐁)2(𝐯𝐄)2/c2,\eta=\frac{\hbar^{2}|e|\gamma\omega}{m_{{\rm{e}}}^{3}c^{5}}\sqrt{(\mathbf{E}+\mathbf{v}\times\mathbf{B})^{2}-(\mathbf{v}\cdot\mathbf{E})^{2}/c^{2}}, (2)

where ω\omega is the angular frequency of electromagnetic field, and the γ\gamma is the Lorentz factor of electrons. There is a critical electric field EsE_{s}, which is equal to 1.32×10161.32\times 10^{16} V cm1\mathrm{cm^{-1}}. When the field strength EE of an electromagnetic wave reaches this critical value, the cascade process is expected to occur rapidly. This leads to a quick decrease in the energy of the electromagnetic field. Lu & Kumar (2019) investigated that this upper limit of EE is a constraint on the maximum luminosity of FRB. Typically, when η5×102\eta\sim 5\times 10^{-2}, the quantum cascade takes place. From Eq. (2), it can be inferred that higher field intensity EE leads to higher reaction rate, resulting in a rapid pair cascade and the formation of a dense plasma.

The cascade process causes the transfer of energy from the electromagnetic field to the plasma, resulting in an increased background plasma density and a decrease in the electromagnetic field strength. To ensure the propagation of electromagnetic waves in the medium, the condition ω>ωp\omega>\omega_{\rm{p}} must be satisfied for electromagnetic waves to propagate, where the ωp\omega_{{\rm{p}}} is the plasma frequency. We have:

ωp=4πne2me,\omega_{{\rm{p}}}=\sqrt{\frac{4\pi ne^{2}}{m_{{\rm{e}}}}}, (3)

where nn represents the plasma density, which must be sufficiently low for electromagnetic waves to propagate. However, the cascade process increases plasma density nn. Therefore, it imposes a constraint on the electric field EE of the electromagnetic wave.

Zhang & Wu (2022) used simulations to determine that the electric field EE of an FRB must be smaller than 3×10123\times 10^{12} V cm1\mathrm{cm^{-1}} (in the laboratory reference frame), ensuring that the FRBs can propagate in the medium. This new constraint is much smaller than the EsE_{s}. There should be a new constraint on the physical properties, such as the radiation radius, the Lorentz factor of the emitting region, and the maximum luminosity. Zhang & Wu (2022) has already placed a lower limit of the radiation radius. In this work, we focus on the Lorentz factor.

3 the state of source

In some models, it is predicted that FRB sources exhibit GRB-like behavior and must be in relativistic motion, with a Lorentz factor Γ102\Gamma\sim 10^{2} (Lyubarsky, 2014). These are the so-called GRB-like models, which radiate FRB from relativistic jet. In this section, we aim to explore the disparities between static and relativistic sources with a focus on their geometries. Notably, highly relativistic motion leads to the beaming effect. Considering FRBs as plane waves, we denote the energy flux density as: S0=E2c/4πS_{0}=E^{2}c/4\pi (in cgs units). Therefore, for an isotropic static source, the FRB luminosity is:

L=E2cR2.L=E^{2}cR^{2}. (4)

For the static and relativistic sources, the Eq. (4) has different forms (see the discussion below).

Assuming the source is isotropic and the radiation is simultaneous, the duration of FRB tt is the result of the time delay during propagation. The scale of the source is:

Rct.R\sim ct. (5)

Due to the differences in geometry, the scale of the source in a static configuration differs from that of a relativistic source.

In the case of relativistic sources, the radiation is concentrated within a cone of angle 1/Γ1/\Gamma, because of the relativistic beaming effect. Therefore, in areas where the angle between the line of sight and the line connecting the observer’s position and the source center is larger than 1/Γ1/\Gamma, the radiation will not propagate to observer. Consequently, not all parts of the source are observable. This geometry differs from that of static sources. In the observer’s line of sight, the source appears as a cone. In this scenario, the radius and time tt obeys:

R2Γ2ct.R\sim 2\Gamma^{2}ct. (6)

The difference between relativistic sources and non-relativistic sources lies in its geometry, as shown in Eq. (5) and Eq. (6).

4 Constraints from the FRB data

In Section 3, it is described that by imposing constraints on Γ\Gamma, it is possible to determine whether Γ\Gamma of the GRB-like model can be confirmed. Applying Eq. (5) to Eq. (4), we have

L=E2c3t2.L=E^{2}c^{3}t^{2}. (7)

Eq. (7) shows the radiation luminosity of a static source and can be plotted on a log-log scale. If the static assumption is not valid (breakthrough the Emax3×1012Vcm1E_{\max}\equiv 3\times 10^{12}\rm{V\;cm^{-1}}), the Eq. (7) should be swapped to ( isotropic equivalent luminosity) 111 Here we assume the relativistic jet opening angle is greater than 1/Γ1/\Gamma. Because of the beaming effect, the main emission is from the region inside 1/Γ1/\Gamma. Therefore, we do not need a jet angle correction for the observed the luminosity.:

L=4E2Γ4c3t2.L=4E^{2}\Gamma^{4}c^{3}t^{2}. (8)

with the same maximum EE, this relativistic source could supply higher luminous events. The difference can be shown in a diagram LtL-t, which will be shown later.

4.1 DM of FRBs

DM describes the effect of propagation for electromagnetic waves, particularly in the case of cosmological distance events such as FRBs. In the range of FRBs band, the different wavelengths propagate in the plasma of IGM with varying velocities. This can be described by the DM:

DM=ne,z1+zdl,\mathrm{DM}=\int\frac{n_{{{\rm{e}}},z}}{1+z}\mathrm{d}l, (9)

where ne,z\rm{n_{e,z}} is the particle number density at redshift zz, and dl\mathrm{d}l is an infinitesimal proper length along the line of sight. DM has three parts, which can be written as :

DMobs=DMMW+DMIGM+DMhost.\rm{DM_{obs}=DM_{MW}+DM_{IGM}+DM_{host}}. (10)

The term DMMW\rm{DM_{MW}} denotes the Milky Way contribution to DM, which can be estimated by model NE2001 (Cordes & Lazio, 2002) and YMW2016 (Yao et al., 2017). NE2001 and YMW2016 are models that describe the electron distribution in the Milky Way. Once the cardinal direction of the FRB (gl gb) has been cleared, these two models can provide an estimate of DMMW\rm{DM_{MW}} for one event. The term DMhost\rm{DM_{host}} represents the contribution to DM from the host galaxy, which is not well studied and is often assumed to follow a log-normal distribution (Macquart et al., 2020). The DMIGM\rm{DM_{IGM}} component accounts for the dispersion measure contributed by the intergalactic medium (IGM), which contains information about the IGM over long distances. DMIGM\rm{DM_{IGM}} can be used to infer the redshift (zz) of the FRB event. The relation between DMIGM\rm{DM_{IGM}} and zz is given by (Zhang, 2022a):

DMIGM=3cH0ΩbfIGM8πGmp0zχ(z)(1+z)dz[Ωm(1+z)3+ΩΛ]1/2,\mathrm{\langle DM_{IGM}\rangle}=\frac{3cH_{0}\Omega_{b}f_{\mathrm{IGM}}}{8\pi Gm_{p}}\int^{z}_{0}\frac{\chi(z)(1+z)dz}{[\Omega_{m}(1+z)^{3}+\Omega_{\Lambda}]^{1/2}}, (11)

where Ωb=0.0486,Ωm=0.3089,ΩΛ=0.6911\Omega_{b}=0.0486,\Omega_{m}=0.3089,\Omega_{\Lambda}=0.6911 are cosmology constants, chosen from the PLANCK observation (Planck Collaboration et al., 2016), and H0H_{0}, fIGMf_{\rm{IGM}}, GG, mpm_{{\rm{p}}} χ(z)\chi(z) denotes Hubble constant, the fraction of baryon mass in the IGM, gravitational constant, the proton mass, the ratio of the number of free electrons to baryons in the IGM, respectively. Due to fluctuations in the IGM, we can only obtain the expectation value of DMIGM\rm{DM_{IGM}}. The fluctuation of DMIGM\rm{DM_{IGM}} can be expressed as (Kumar & Linder, 2019):

σΔIGMDMIGM=0.2z.\mathrm{\frac{\sigma_{\Delta IGM}}{\langle DM_{IGM}\rangle}}=\frac{0.2}{\sqrt{z}}. (12)

4.2 Luminosity – duration diagram of FRBs

After obtaining DMIGM\rm{DM_{IGM}} and inferring the redshift zz using Eq. (11), we can calculate the uncertainty of DMIGM\rm{DM_{IGM}} using Eq. (12). This enables us to determine the uncertainty in the zz 222Fν,pF_{\nu,p} also will affect the uncertainty of luminosity LL. However, Fν,pF_{\nu,p} is detected precisely. The main uncertainty comes from the detection of z.. The luminosity distance DLD_{{\rm{L}}} can be inferred from zz:

L=4πDL2Fν,pν,L=4\pi D_{{{\rm{L}}}}^{2}F_{\nu,p}\nu, (13)

where DLD_{{\rm{L}}} is the luminosity distance, which can be determined from the redshift zz, and the Fν,pF_{\nu,p} is the peak flux density. Since observations cannot cover entire band of an FRB, we need to estimate the ν\nu. One common estimate is to take the center frequency νc\nu_{c} as a rough approximation. Another option is to use the telescope bandwidth νo\nu_{o} as a conservation estimate. DLD_{{\rm{L}}} can be firmed by redshift zz, and the Fν,p,νcF_{\nu,p},\nu_{c} are observational data. We use the data from the FRBSTATS Catalogue333https://www.herta-experiment.org/frbstats/catalogue, which are from CHIME, GBT, FAST, ASKAP, and others (Petroff et al., 2016; Petroff, 2020; Spanakis-Misirlis, 2021). The data provided in the catalogue include peak flux density, inferred redshift zz, bandwidth, and center frequency νc\nu_{c}. For repeating bursts, only the highest signal-to-noise event is included. There are 828 events (source) in total in this dataset. Redshift zz is induced by DM as described in Section 4.1. The data is updated to December 2022, and the latest event is FRB20221128. The scattered plot of LL (luminosity) - tt (duration) is shown in Fig. 1. We only choose those FRBs with luminosity larger than 104510^{45} erg s1\rm{s^{-1}} as the even dimmer FRBs are unlikely to break the limit line as given in Eq. (7), taking the Emax=3×1012Vcm1E_{\max}=3\times 10^{12}\,\rm{V\,cm^{-1}} into this equation. This limitation line is shown as a straight line in Fig. 1. For a certain FRB locating on the right side of this line, it indicates at least this FRB should not radiate stationarily. Since Zhang & Wu (2022) only simulated the electromagnetic wave interaction of the 1 GHz FRB, we only chose these FRBs with frequencies in the range of 0.60.6-2 GHz.

Refer to caption
Figure 1: The LL (luminosity)- tt (duration) scattered plot of FRBs. Luminosity is calculated from Eq. (13) with ν=νc\nu=\nu_{c}. The tt is full width at half maximum of the pulse. The black solid line is plotted following the Eq. (7) at E=Emax{E=E_{\max}}. Two FRBs (FRB20150807A and FRB20191108A) of them locate at the more luminous side.

Based on Fig. 1, it can be seen that the majority of the FRB events are compatible with the stationary scenario. However, there are two events FRB20150807A (duration tt=0.35 ms, redshift zz=0.27, and peak flux Fν,pF_{\nu,p}=128 Jy ms) and FRB20191108A (duration tt=0.34ms, redshift zz=0.56, and peak flux Fν,pF_{\nu,p}=27 Jy ms) (Petroff, 2020) that transcend the limit of the constraints, indicating that they may need to be relativistic. Considering that they are not clearly far from the limitation line, we cannot conclude that FRBs should radiate in relativistic velocity, and cannot constrain the Lorentz factor neither.

5 conclusion and discussion

In this work, we examined whether some bright, short duration FRBs could originate from relativistic sources, based on a quantum cascade limit on the electric field strength. The luminosity-duration distribution for 828 FRB events were analyzed. Most FRBs fell below the limit expected for a static, non-relativistic source. However, the two bursts had a higher luminosity and a shorter duration, suggesting they may originate from relativistic sources. This could support GRB-like radiation models for part of the FRBs. However, more extreme FRBs are needed to strengthen this conclusion.

If confirmed, the Lorentz factor distribution could distinguish GRB-like and pulsar-like models. The progenitor environments and emission physics may also differ between relativistic and non-relativistic FRB sources. For the GRB-like model, it predicts multiwavelength afterglow. Further multi-wavelength follow-up will shed light on the origin. The radio spectra and polarization properties could also diagnose relativistic boosting effects. Ultimately, identifying more high luminosity, short duration FRBs and modeling their radiation mechanisms will help unravel the origin of these enigmatic signals.

Acknowledgements

We thank the helpful discussions with Enping Zhou, Hao Wang, Shiyan Tian, Lin Zhou, Yuanyuan Zuo and Ruihang Dong. The English is polished with ChatGPT. This work is supported by the National Natural Science Foundation of China (Grant No. 12041306).

Data Availability

The data used are publicly available in the FRBSTATS catalog (https://www.herta-experiment.org/frbstats/catalogue).

References

  • Bochenek et al. (2020) Bochenek C. D., Ravi V., Belov K. V., Hallinan G., Kocz J., Kulkarni S. R., McKenna D. L., 2020, Nature, 587, 59
  • Cordes & Lazio (2002) Cordes J. M., Lazio T. J. W., 2002, arXiv e-prints, pp astro–ph/0207156
  • Kashiyama et al. (2013) Kashiyama K., Ioka K., Mészáros P., 2013, ApJ, 776, L39
  • Kirk et al. (2009) Kirk J. G., Bell A. R., Arka I., 2009, Plasma Physics and Controlled Fusion, 51, 085008
  • Kumar & Linder (2019) Kumar P., Linder E. V., 2019, Phys. Rev. D, 100, 083533
  • Landau & Lifshitz (1982) Landau L., Lifshitz E., 1982, Quantum Electrodynamics. Pergamon Press
  • Lorimer et al. (2007) Lorimer D. R., Bailes M., McLaughlin M. A., Narkevic D. J., Crawford F., 2007, Science, 318, 777
  • Lu & Kumar (2019) Lu W., Kumar P., 2019, MNRAS, 483, L93
  • Lyubarsky (2014) Lyubarsky Y., 2014, MNRAS, 442, L9
  • Macquart et al. (2020) Macquart J. P., et al., 2020, Nature, 581, 391
  • Melikidze et al. (2000) Melikidze G. I., Gil J. A., Pataraya A. D., 2000, ApJ, 544, 1081
  • Petroff (2020) Petroff E., 2020, Transient Name Server Fast Radio Bursts, 2470, 1
  • Petroff et al. (2016) Petroff E., et al., 2016, Publ. Astron. Soc. Australia, 33, e045
  • Planck Collaboration et al. (2016) Planck Collaboration et al., 2016, A&A, 594, A13
  • Ravi et al. (2019) Ravi V., et al., 2019, Nature, 572, 352
  • Schwinger (1951) Schwinger J., 1951, Phys. Rev., 82, 664
  • Spanakis-Misirlis (2021) Spanakis-Misirlis A., 2021, FRBSTATS: A web-based platform for visualization of fast radio burst properties, Astrophysics Source Code Library, record ascl:2106.028 (ascl:2106.028)
  • Yao et al. (2017) Yao J. M., Manchester R. N., Wang N., 2017, ApJ, 835, 29
  • Zhang (2020) Zhang B., 2020, Nature, 587, 45
  • Zhang (2022a) Zhang B., 2022a, arXiv e-prints, p. arXiv:2212.03972
  • Zhang (2022b) Zhang B., 2022b, ApJ, 925, 53
  • Zhang & Wu (2022) Zhang Y., Wu H.-C., 2022, ApJ, 929, 164