Superconducting cosmic strings as sources of cosmological fast radio bursts

TL;DR

This study models radio bursts from superconducting cosmic strings, constrains their parameters using fast radio burst data, and finds cusps likely dominate the signals, providing new insights into cosmic string properties.

Contribution

First to constrain superconducting cosmic string parameters using fast radio burst observations and analyze contributions from cusps and kinks.

Findings

Parameter space of Gμ ~ [10^{-14}, 10^{-12}] and I ~ [10^{-1}, 10^{2}] GeV fits the data.

Cusps likely dominate the radio burst signals over kinks.

Two best-fit models: one cusp-dominated at z=1.3, another kink-dominated at maximal event rate.

Abstract

In this paper we calculate the radio burst signals from three kinds of structures of superconducting cosmic strings. By taking into account the observational factors including scattering and relativistic effects, we derive the event rate of radio bursts as a function of redshift with the theoretical parameters $G\mu$ and $\mathcal{I}$ of superconducting strings. Our analyses show that cusps and kinks may have noticeable contributions to the event rate and in most cases cusps would dominate the contribution, while the kink-kink collisions tend to have secondary effects. By fitting theoretical predictions with the normalized data of fast radio bursts, we for the first time constrain the parameter space of superconducting strings and report that the parameter space of $G\mu \sim [10^{-14}, 10^{-12}]$ and $\mathcal{I} \sim [10^{-1}, 10^{2}] ~ \rm{GeV}$ fit the observation well although the statistic significance is low due to the lack of observational data. Moreover, we derive two types of best fittings, with one being dominated by cusps with a redshift $z = 1.3$, and the other dominated by kinks at the range of the maximal event rate.

Tables

Table 1. Table 1: FRB catalog as detected by Parkes Petroff:2016tcr .

Event	Redshift	$\mathrm{DM}\left[{\mathrm{cm}}^{-3}\mathrm{pc}\right]$	$W_{\mathrm{obs}}\left[\mathrm{ms}\right]$	$S_{\mathrm{peak},\mathrm{obs}}\left[\mathrm{Jy}\right]$	$F_{\mathrm{obs}}\left[\mathrm{Jy}\mathrm{ms}\right]$
FRB010125	0.57	790	${9.40}_{-0.20}^{+0.20}$	0.30	2.82
FRB010621	0.19	745	7.00	0.41	2.87
FRB010724	0.28	375	5.00	$>{30.00}_{-10.00}^{+10.00}$	$>$150.00
FRB090625	0.72	899.55	${1.92}_{-0.77}^{+0.83}$	${1.14}_{-0.21}^{+0.42}$	${2.19}_{-1.12}^{+2.10}$
FRB110220	0.76	944.38	${5.60}_{-0.10}^{+0.10}$	${1.30}_{-0.00}^{+0.00}$	${7.28}_{-0.13}^{+0.13}$
FRB110626	0.56	723.0	1.40	0.40	0.56
FRB110703	0.89	1103.6	4.30	0.50	2.15
FRB120127	0.43	553.3	1.10	0.50	0.55
FRB121002	1.3	1629.18	${5.44}_{-1.20}^{+3.50}$	${0.43}_{-0.06}^{+0.33}$	${2.34}_{-0.77}^{+4.46}$
FRB130626	0.74	952.4	${1.98}_{-0.44}^{+1.20}$	${0.74}_{-0.11}^{+0.49}$	${1.47}_{-0.50}^{+2.45}$
FRB130628	0.35	469.88	${0.64}_{-0.13}^{+0.13}$	${1.91}_{-0.23}^{+0.29}$	${1.22}_{-0.37}^{+0.47}$
FRB130729	0.69	861	${15.61}_{-6.27}^{+9.98}$	${0.22}_{-0.05}^{+0.17}$	${3.43}_{-1.81}^{+6.55}$
FRB131104	0.59	779	2.08	1.12	2.33
FRB140514	0.44	562.7	${2.80}_{-0.70}^{+0.35}$	${0.47}_{-0.08}^{+0.11}$	${1.32}_{-0.50}^{+2.34}$
FRB150418	0.49	776.2	${0.80}_{-0.30}^{+0.30}$	${2.20}_{-0.30}^{+0.60}$	${1.76}_{-0.81}^{+1.32}$

Equations

$$S=\int d^2\zeta$$

$$-\frac{1}{16\pi }\int d^{4}x\sqrt{-g}{F}_{\mu \nu }{F}^{\mu \nu },$$

$$J^a=q{ϵ}^{ab}{\varphi }_{,b}/\sqrt{-\gamma },$$

$${\dot{x}}^{\mu }{x}_{\mu }^{\prime }=0,{\dot{x}}^2+x^{\prime 2}=0,$$

$${\square }_2\varphi =-\frac{1}{2}q{ϵ}^{ab}{F}_{\mu \nu }{x}_{,a}^{\mu }{x}_{,b}^{\nu },$$

$$\displaystyle\Box_{2}x^{\mu}=\big{[}(-F^{\mu}_{~{}\sigma}x^{\sigma}_{,a}J^{a}-(\Theta^{ab}x^{\mu}_{,a})_{,b}\big{]}/\mu~{},$$

$${\square }_4{A}^{\mu }=4\pi j^{\mu }.$$

$${\square }_2={\eta }^{ab}{\partial }_a{\partial }_b,{\square }_4={\eta }^{\mu \nu }{\partial }_{\mu }{\partial }_{\nu },$$

$${\Theta }_{ab}={\varphi }_{,a}{\varphi }_{,b}-\frac{1}{2}{\gamma }_{ab}{\varphi }_{,c}{\varphi }^{,c},$$

$$\displaystyle j^{\mu}(x)=\int d\sigma d\tau x^{\mu}_{,a}\sqrt{-\gamma}~{}J^{a}~{}\delta^{(4)}\big{(}x-x(\sigma,\tau)\big{)}~{},$$

$$x^{\mu}(t,\sigma)=\frac{1}{2}\big{[}x^{\mu}_{-}(\sigma_{-})+x^{\mu}_{+}(\sigma_{+})\big{]}~{},$$

$$x_{-}^0=-{\sigma }_{-},x_{+}^0={\sigma }_{+},|{\mathbf{x}}_{+}^{\prime }|=1=|{\mathbf{x}}_{-}^{\prime }|.$$

$$j^{\mu}(t,\mathbf{x})=\mathcal{I}\int d\sigma x^{\mu}_{,\sigma}\delta^{(3)}\big{(}\mathbf{x}-\mathbf{x}(t,\sigma)\big{)}~{},$$

$$j^{\mu }(t,\mathbf{x})=\underset{\omega }{\Sigma }\int \frac{d^{3}k}{(2\pi {)}^3}{e}^{-i(\omega t-\mathbf{k}\cdot \mathbf{x})}{j}_{\omega }^{\mu }(\mathbf{k}),$$

$$j_{\omega }^{\mu }(\mathbf{k})=\frac{2\mathcal{I}}{L}{\int }_0^{L/2}dt{\int }_0^{L}d\sigma e^{i{k}_{\mu }{x}^{\mu }}{x}^{{}^{\prime }\mu }\equiv \frac{\mathcal{I}}{L}(I_{n+}^{\mu }+I_{n-}^{\mu }).$$

$$I_{n\pm }^{\mu }={\int }_0^{L/2}dt{\int }_0^{L}d{\sigma }_{\pm }{e}^{i{k}_{\mu }{x}_{\pm }^{\mu }}{x}_{\pm }^{{}^{\prime }\mu }.$$

$$P$$

$$P_{{\omega }_n}$$

$$j_{\omega }^{\mu }(\mathbf{k})$$

$$\displaystyle=\frac{2\mathcal{I}}{L}\int^{L/2}_{0}dt\int^{L}_{0}d\sigma~{}e^{i(\omega t-\mathbf{k}\cdot\mathbf{x}_{+}/2-\mathbf{k}\cdot\mathbf{x}_{-}/2)}\frac{1}{2}\big{[}x^{{}^{\prime}\mu}_{+}+x^{{}^{\prime}\mu}_{-}\big{]}~{}.$$

$$\mathbf{n}\cdot {\mathbf{x}}_{+}^{\prime }+\mathbf{n}\cdot {\mathbf{x}}_{-}^{\prime }=0,$$

$${\mathbf{x}}_{+}$$

$${\mathbf{x}}_{-}$$

$${\mathbf{x}}_{+0}{}^{\prime }=-{\mathbf{x}}_{-0}{}^{\prime },$$

$${\mathbf{x}}_{+0}{}^{\prime }\cdot {\mathbf{x}}_{+0}{}^{\prime \prime }={\mathbf{x}}_{-0}{}^{\prime }\cdot {\mathbf{x}}_{-0}{}^{\prime \prime }=0,$$

$${\mathbf{J}}_{\omega }(\mathbf{k})$$

$$\sim \frac{2\mathcal{I}}{L}|t_{\ast }|\frac{1}{L}|{\sigma }_{\ast }{|}^2.$$

$$\omega t-\mathbf{k}\cdot \mathbf{x}=$$

$$\displaystyle-\frac{1}{4}\mathbf{n}\cdot\big{[}\mathbf{x}_{+0}^{\prime\prime}(\sigma+t)^{2}+\mathbf{x}_{-0}^{\prime\prime}(\sigma-t)^{2}\big{]}$$

$$\displaystyle-\frac{1}{12}\mathbf{n}\cdot\big{[}\mathbf{x}_{+0}^{\prime\prime\prime}(\sigma+t)^{3}+\mathbf{x}_{-0}^{\prime\prime\prime}(\sigma-t)^{3}\big{]}+\cdot\cdot\cdot\Big{\}}~{}.$$

$$|t_{\ast }|,|{\sigma }_{\ast }|\sim L^{2/3}/{\omega }^{1/3}.$$

$$\omega t-\mathbf{k}\cdot \mathbf{x}\sim {\theta }^2\omega t\sim \theta \omega L^{-1}(t^2,t\sigma ,{\sigma }^2)\sim 1,$$

$$\theta \sim \frac{1}{(\omega L{)}^{1/3}}.$$

$$\Omega \sim {\theta }^2\sim \frac{1}{(\omega L{)}^{2/3}}.$$

$$P_{\omega }=\frac{{\mathcal{I}}^2{L}^{1/3}}{{\omega }^{2/3}}.$$

$$P_{\gamma }^c\sim {\mathcal{I}}^2({\omega }_{max}L{)}^{1/3}\sim \kappa \mathcal{I}\sqrt{\mu }.$$

Full text

I Introduction

Cosmic strings can be formed in phase transitions of the very early universe if the manifold of the background vacuum has nontrivial topology for the spontaneous symmetry breaking (see Refs. Kibble:1976sj ; Vilenkin:2000jqa ; Polchinski:2004ia ; Copeland:2009ga ; Ringeval:2010ca ; Copeland:2011dx ; Vachaspati:2015cma for comprehensive reviews). The potential exploration of these primordial relics could reveal important information as regards fundamental physics at extremely high energy scales, in which one class of particle physics models naturally give rise to cosmic strings that are superconducting Witten:1984eb . Superconducting cosmic strings (SCSs), as can be realized by a charged scalar field whose flux is trapped in normal cosmic strings with the electromagnetic gauge invariance broken inside the strings, can produce electromagnetic effects Vilenkin:2000jqa . These superconductive wire-like objects in the sky behave as giant antennas that can emit electromagnetic signals in a wide range of frequencies Vilenkin:1986zz ; Garfinkle:1987yw .

Due to the oscillations of cosmic strings, they are able to emit electromagnetic radiation as well as gravitational radiation and thus can give rise to various observational effects of SCSs. These effects have been extensively studied in the literature, such as the distortions and anisotropies in the cosmic microwave background (CMB) caused by electromagnetic or gravitational radiations of cosmic strings Sanchez:1988ek ; Sanchez:1990kj ; Tashiro:2012nb ; Jazayeri:2017szw , the early reionization by SCSs and the associated imprints in the CMB Tashiro:2012nv , the high energy cosmic rays due to the decay of cosmic strings Hill:1986mn ; Bhattacharjee:1989vu ; Bhattacharjee:1990js ; MacGibbon:1989kk ; Wichoski:1998kh ; Brandenberger:2009ia ; Berezinsky:2009xf ; Lunardini:2012ct , the gamma-ray bursts induced by cosmic strings Babul:1987lza ; Paczynski:1988 ; Brandenberger:1993hw ; Plaga:1993hp ; Berezinsky:2001cp ; Cheng:2010ae , gravitational waves Vachaspati:1984gt ; Damour:2000wa ; Damour:2001bk ; Damour:2004kw , as well as sparks Vachaspati:2008su and radio transients from SCSs Cai:2011bi ; Cai:2012zd ; Yu:2014gea . Accordingly, the parameter space of SCSs, which is often characterized by the string tension $\mu$ (or $G\mu$ in Planck units) and the current on the string $\mathcal{I}$, can be derived from the comparison with various observational data. For instance, the CMB analysis based on WMAP and the South Pole Telescope yields an upper bound on the string tension as $G\mu<1.7\times 10^{-7}$ Dvorkin:2011aj , and this constraint has been recently improved to be $G\mu<1.3\times 10^{-7}$ with the help of the Planck data Ade:2013xla . As cosmic strings can also lead to relic gravitational waves in the form of a stochastic background, a model-dependent constraint from the pulsar timing measurements on the string tension can be obtained in about the same order of the bound from CMB vanHaasteren:2011ni ; Pshirkov:2009vb ; Tuntsov:2010fu ; Olmez:2010bi ; Binetruy:2012ze ; Sanidas:2012ee . Further constraints on SCSs can be derived through the spectral distortions of the CMB photons Sanchez:1988ek ; Sanchez:1990kj ; Tashiro:2012nb , namely, the parameter region with $10^{-19}<G\mu<10^{-7}$ and $\mathcal{I}>10^{4}~{}{\rm GeV}$ would be ruled out Kogut:2011xw .

The signals of electromagnetic radiations from SCSs can be significantly enhanced if they were emitted from particular structures of strings, such as the cusps Vachaspati:2008su , and the kinks (Cai:2011bi, ; Cai:2012zd, ). These signals are found to be in form of radio bursts, and thus, the SCSs could serve as the possible sources of the Fast Radio Burst (FRB) events that have been discovered in astronomical observations. So far, there have been $18$ FRB events reported by various radio experiments including the Parkes Lorimer:2007qn ; Keane:2011mj ; Thornton:2013iua ; Burke-Spolaor:2014rqa ; Petroff:2014taa ; Ravi:2014mma ; Champion:2015pmj ; Keane:2016yyk ; Ravi:2016kfj , the Arecibo Pulsar ALFA Survey Spitler:2014fla , and the Green Bank Telescope Masui:2015kmb . The events of FRB are characterized as having a bandwidth of the millisecond scale, and frequency around ${\rm GHz}$ scale. All these events show extra-galactic origins and their estimated redshifts are within $z=1.4$. The flux of these events are about several ${\rm Jys}$. Theoretically there have been some attempts to provide the sources of these FRB signals, such as super-massive neutron stars Falcke:2013xpa , binary neutron stars or white dwarf mergers Totani:2013lia ; Kashiyama:2013gza , local circumnuclear magnetars Pen:2015ema , and so on.

In the present study we investigate the radio burst signals from SCSs using the method of order-of-magnitude analysis. We compare the normalized FRB event rate as a function of redshift based on theoretical models with normalized observational data, and thus derive the limits on theoretical parameters of strings. Although there are not enough data to offer sufficiently precise results, our analysis provides a theoretically sound method of constraining theoretical parameters of SCSs in the ongoing search for FRBs at the Parkes Lorimer:2007qn , ETA Patterson:2008ie , LWA Henning:2010aa , LOFAR Fender:2008zh , and FAST Nan:2011um radio telescopes. We note that an attempt on the astronomical limit of SCSs based on the specific data of a single FRB event detected by the Parkes was investigated in Miyamoto:2012ck , and later, the work of Yu:2014gea has discussed FRBs as a possible probe on model parameters of SCSs such as string length and $\mu$ by ascribing four FRB events. The model of calculating the radiations from SCSs applied in these two works only focused on the particular structure of string cusps. However, as we shall show in the present work, radiations from kinks can also lead to significant contribution under certain conditions. Moreover, we in the present study take into account the scattering effects from the interstellar medium (ISM) and intergalactic medium (IGM) as well as the statistical distributions of string loops. With the help of accumulated astronomical data, our analysis shows that the statistical bias can be greatly reduced and hence the astronomical constraints on the parameter space of SCSs from the FRB events would be statistically meaningful in near future.

The structure of the present paper is as follows. In Sec. II we describe the general characteristics of SCSs by reviewing the radiation emission from these charged wires in the sky. Then in Sec. III we perform an order-of-magnitude analysis method to solve out the equations of electromagnetic power emitted from string structures including cusps and kinks. Sec. IV is devoted to the calculation of the theoretical prediction on the event rate of radio bursts as a function of the redshift, flux and string length parameters. Afterwards, we in Sec. V perform the numerical estimation in detail by comparing the theoretical prediction with the observational data of FRBs, and then report a brand-new constraint on the parameter space of cosmic strings. We summary the results with a discussion in Sec VI. In the theoretical calculations we use natural units, $\hbar=c=1$.

II Characteristics of SCSs

To begin with, we consider the effective action of a superconducting string as follows Vilenkin:2000jqa :

[TABLE]

where $\gamma_{ab}$ is the induced metric on the worldsheet of a string with $(a,b)=0\ {\rm or}\ 1$, and $\mu$ is the string tension. The first term in the rhs of the above action is the Nambu-Goto action; the second and third terms represent for the worldsheet current and the coupling to the electromagnetic gauge potential, respectively, with $x^{\mu}$ being introduced as the string position, $\phi$ the massless real scalar field on the worldsheet, and $A_{\mu}$ the gauge field of the worldsheet current $J^{a}$; and the last term describes the dynamics of electromagnetic field in the four dimensional spacetime.

The current on the worldsheet is given by

[TABLE]

where $q$ is the charge of the current carriers. We focus on the classical production of bursts and at the moment simplify the metric of the background spacetime as the Minkowski $\eta_{\mu\nu}$. We choose the conformal gauge for the worldsheet, which can be specified by the following conditions:

[TABLE]

where “ $\dot{~{}}\equiv\partial/\partial\tau$” and “ ${}^{\prime}\equiv\partial/\partial\sigma$” has been introduced, and $(\tau,\sigma)$ correspond to the worldsheet coordinates. Under this condition, the induced worldsheet metric is then simplified as $\gamma^{ab}=\rm{diag}(1,-1)$. Consequently, we write down the equations of motion for the scalar field, the string and the gauge field respectively as,

[TABLE]

In the above equations of motion, we have introduced

[TABLE]

where $\Theta_{ab}$ is the stress energy tensor of the scalar field living on the worldsheet and $j^{\mu}$ corresponds to the current in the four dimensional spacetime. We adopt the Lorentz gauge for $A^{\mu}$: $\partial_{\mu}A^{\mu}=0$. In order to fix another degree of freedom, we also set $x^{0}=\tau=t$ with $t$ being the physical time. We use $t$ instead of $\tau$ to denote $x^{0}$ in the following discussions unless otherwise claimed.

In the center-of-mass frame of a string loop with an invariant length $L$, we can calculate the solution for Eq. (5) to be

[TABLE]

where $\sigma_{\pm}\equiv\sigma\pm t$. Since the periodicity of $\sigma$ takes the string length $L$, from the solution we can deduce that the periodicity of $t$ is $L/2$. To combine the gauge conditions in Eq. (3) with the solution, we can get

[TABLE]

As a result, for a string with length $L$, the order of magnitude of $\mathbf{x}^{\prime\prime}_{\pm}$ and $\mathbf{x}^{\prime\prime\prime}_{\pm}$ can be estimated to be ${1}/{L}$, ${1}/{L^{2}}$ respectively.

Note that the vanishing rhs of Eq. (4) implies a conservation of the current $J^{a}$ on the worldsheet. To integrate Eq. (9) over $t$ for a period, we then have

[TABLE]

where $\mathcal{I}$ is the parameter of the current on the string. After that, we perform the Fourier transform

[TABLE]

which yields

[TABLE]

In the above equation we have applied the definition of $x^{\mu}_{\pm}$ and thus there is

[TABLE]

The electromagnetic radiation power from a periodic source is then given by

[TABLE]

where $\omega_{n}$ denotes the Fourier mode with $\omega_{n}=4\pi n/L$, $n\in\mathbf{N}$. From the expression of $P_{\omega_{n}}$, we can see that the key to estimating the radiation power from SCSs is to solve out the solution of $j^{\mu}_{\omega}$. The main approach of releasing electromagnetic radiation from the string loop is to study the emission processes from certain structures including cusps, kinks and kink-kink collisions. Other than these structures, the radiation power would decrease exponentially and hence is ignored for simplicity. In the study of gravitational or electromagnetic radiations from cosmic strings, it is common to make use of the saddle point/discontinuity approximation method to find the corresponding solutions (for instance see Refs. (Damour:2004kw, ; Steer:2010jk, ; Cai:2012zd, )). In order to carry out more details of the same processes, we in the present analysis follow the method in Refs. (Vilenkin:1986zz, ; BlancoPillado:2000xy, ), which was adopted in investigating the string cusps, and, perform a much illustrative and generalized approach to solve the solutions to particular structures of cosmic strings.

III Radiations from SCSs

The radiations from string loops are mostly contributed by the areas where $j^{\mu}_{\omega}$ takes the maximal value. As we will show later, outside these areas the contributions would decrease exponentially. For convenience, one can set the position of the maximal value for $j^{\mu}_{\omega}$ to be the origin point $(t,\sigma)=(0,0)$, and then express $j^{\mu}_{\omega}$ in the following form:

[TABLE]

The wave vector $\mathbf{k}$ is also identified as $\mathbf{k}=\omega\mathbf{n}$. Note that Eq. (11) yields $j^{0}_{\omega}=0$. Therefore, in the following study we focus on the spatial part of the four current $j^{\mu}_{\omega}$: $\mathbf{J}_{\omega}$. The maximal value position requires $(\omega t-\mathbf{k}\cdot\mathbf{x})^{\prime}=0$, where $(\omega t-\mathbf{k}\cdot\mathbf{x})$ is the exponential factor in the expression of $j^{\mu}_{\omega}$. From this requirement, one derives

[TABLE]

and for $\mathbf{n}$ not perpendicular to $\mathbf{x}_{\pm}^{\prime}$, it is equivalent to that $\mathbf{x}_{+}^{\prime}+\mathbf{x}_{-}^{\prime}=0$. Such particular positions are denoted as $\it{cusps}$, if both $\mathbf{x}_{+}^{\prime}$ and $\mathbf{x}_{-}^{\prime}$ are continuous; or, $\it{kinks}$, if one of $\mathbf{x}_{\pm}^{\prime}$ is discontinuous; or even, $\it{kink-kink\ collisions}$, if both $\mathbf{x}_{\pm}^{\prime}$ are discontinuous. We analyze the electromagnetic radiation power emitted around these particular structures, respectively.

III.1 Cusps

In the neighborhood of the origin point $(t,\sigma)=(0,0)$, we can perform the Taylor expansions of $\mathbf{x}_{\pm}$ as follows:

[TABLE]

and these expressions are valid as long as the scales involved are smaller than the string length $\mathcal{O}(L)$). To substitute Eq. (11) into the above expansions, one obtains

[TABLE]

and accordingly, from Eq. (III) the form of $\mathbf{J}_{\omega}$ can be estimated to be

[TABLE]

In the above formulae, $t_{*}$ and $\sigma_{*}$ represents for the characteristic regime when the radiations dominate over. In order to estimate their magnitudes, we apply the condition $(\omega t-\mathbf{k}\cdot\mathbf{x})\simeq 1$ and then perform the Taylor expansion of $\mathbf{k}\cdot\mathbf{x}$ around $(t,\sigma)=(0,0)$ as follows:

[TABLE]

Notice that approximately $\mathbf{x}_{\pm}^{\prime\prime}\sim 1/L$ and $\mathbf{x}_{\pm}^{\prime\prime\prime}\sim 1/L^{2}$. Moreover, Eq. (22) shows that $\mathbf{x}_{+}^{\prime}$ and $\mathbf{x}_{-}^{\prime}$ are anti-parallel, and thus, from Eq. (19) one can learn that the orders of magnitude for $(t,\sigma)$ should be small quantities. In addition, if $\mathbf{n}$ is parallel to either $\mathbf{x}_{+}^{\prime}$ or $\mathbf{x}_{-}^{\prime}$, Eqs. (11), (19) and (23) then ensure $\mathcal{O}(t^{2},\sigma^{2})$ to be also negligible. As a consequence, we can reach $\omega\mathcal{O}(t^{3},\sigma^{3})/L^{2}\sim 1$, which implies that the most possible orders of magnitude for $|t_{*}|$ and $|\sigma_{*}|$ in the case of cusps can be estimated as

[TABLE]

Similar to the above process, we can also estimate the order of magnitude of the solid angle $\Omega$ as follows. We consider that $\mathbf{n}$ has a small deviation from $\mathbf{x}_{\pm}^{\prime}$ within the range of the angle $\theta$. Without loss of generality, we take $\mathbf{n}\cdot\mathbf{x}_{+}^{\prime}=cos\theta$, $\mathbf{n}\cdot\mathbf{x}_{-}^{\prime}=cos(\pi-\theta)$, and hence, from Eq. (III.1) get the following approximations:

[TABLE]

which further yields

[TABLE]

All possible deviations of the angle $\theta$ can form a cone-shape region around $\mathbf{x}_{\pm}^{\prime}$, and thus, the solid angle scope where the radiations could be observed from a cusp is limited within

[TABLE]

Moreover, the radiations from the cusp would decrease exponentially outside the solid angle $\Omega$, which implies that the electromagnetic signals emitted from a cusp are expected to be highly focused beams. On substituting Eq. (26) into Eq. (III.1), we get $\mathbf{J}_{\omega}\sim\frac{\mathcal{I}}{\omega}$ approximately. Using Eq. (17), the power emitted by string cusps at the frequency $\omega$ can be estimated as

[TABLE]

The total power emitted from this type of the string structure can be obtained by integrating $P_{\omega}$ over the frequency $\omega$. This result looks divergent. This is because that cosmic strings actually remains the inner structure characterized by its tension $\mu$. Therefore, we can provide an upper limit of $\omega$ by the theoretical constraint that, given the tension $\mu$ of a string, the maximum energy of radiations emitted from this string ought to be less than the energy of the string itself. As a result, the radiation power from cusps can then be estimated as Vilenkin:1986zz ; Cai:2012zd

[TABLE]

Numerical simulations yield a result in the same order of magnitude, with a correction of the factor $\kappa\sim 10$.

III.2 Kinks

In the case of string kinks, either ${x}^{\mu}_{+}{{}^{\prime}}$ or ${x}^{\mu}_{-}{{}^{\prime}}$ is discontinuous. Without loss of generality, we set that ${x}^{\mu}_{+}{{}^{\prime}}$ is discontinuous at $\sigma=0$. Similar to Eq. (III.1), we make the Taylor expansion of $(\omega t-\mathbf{k}\cdot\mathbf{x})$ from both sides, which is expressed as

[TABLE]

where the subscript $(0\pm)$ denotes the expansions from the positive and negative side, respectively. The coefficient $\sigma$ takes $-\frac{1}{2}(\mathbf{n}\cdot\mathbf{x}_{+}{{}^{\prime}}_{(0\pm)}+\mathbf{n}\cdot\mathbf{x}_{-}{{}^{\prime}})$, where $\mathbf{x}_{-}{{}^{\prime}}$ is parallel with one side of $\mathbf{x}_{+}{{}^{\prime}}$ by assuming that it is $\mathbf{x}_{+}{{}^{\prime}}_{(0+)}$. From Eq. (19), this coefficient must vanish, thus the direction of $\mathbf{n}$ can only be in a plane such that the angle between $\mathbf{n}$ and $\mathbf{x}_{+}{{}^{\prime}}_{(0+)}$ equals that between $\mathbf{n}$ and $\mathbf{x}_{+}{{}^{\prime}}_{(0-)}$. We can then check the $\sigma^{2}$ coefficient. Note that, for both sides of the kinks, Eqs. (19) and (23) remain applicable, but here we can only derive $\mathbf{n}~{}\cdot~{}\mathbf{x}_{+}{{}^{\prime\prime}}_{(0\pm)}~{}+~{}\mathbf{n}~{}\cdot~{}\mathbf{x}_{-}{{}^{\prime\prime}}~{}=~{}0$, and accordingly, the $\sigma^{2}$ coefficient requires $\sigma\sim L/(\omega|t_{*}|)$; however, the $\sigma^{3}$ coefficient further requires $\sigma\sim L^{2/3}/\omega^{1/3}$. To take into account both estimations, one can set the final form of $|\sigma_{*}|$ to be the minimal of these two results. Given the relations between $\mathbf{n}$ and $\mathbf{x}_{+}{{}^{\prime}}_{(0\pm)},\mathbf{x}_{-}{{}^{\prime}}$, the first-order coefficient of $t$ would not vanish. As a result, we have $|t_{*}|\sim\omega^{-1}$. Substituting $|t_{*}|$ into the $\sigma^{2}$ coefficient constraint, we find that the dominant order of magnitude for $\sigma$ comes from the $\sigma^{3}$ term in the case of kinks. To conclude, $(\omega t-\mathbf{k}\cdot\mathbf{x})\sim 1$ yields the order of $t$ and $\sigma$ for kinks as follows:

[TABLE]

Unlike cusps, for kinks, the constraint of the angle of the radiation emission, as we have discussed above, has only one dimension. By substituting Eq. (33) into (27), one gets

[TABLE]

We note that the only restriction of $\theta$ comes from the requirement of $\theta\omega L^{-1}\sigma^{2}\sim 1$. Hence the radiation from a kink is emitted in a “disk-shape” set of directions of solid angle. Then $\mathbf{J}_{\omega}$ becomes

[TABLE]

where $|\Delta\mathbf{x}{{}^{\prime}}_{+}|$ denotes the difference between $\mathbf{x}{{}^{\prime}}_{+(0+)}$ and $\mathbf{x}{{}^{\prime}}_{+(0-)}$. To apply Eq. (33), we can derive

[TABLE]

and according to Eq. (17),

[TABLE]

where $\Psi_{+}=|\Delta\mathbf{x}{{}^{\prime}}_{+}|$ denotes the sharpness of $\mathbf{x}{{}^{\prime}}_{+}$. For a loop with $N$ kinks and sharpness $\Psi_{+}$, the total radiation can be calculated to be

[TABLE]

Here $\omega_{max}\sim\sqrt{\mu}$, which is determined by the order of the inverse width of the string. The low frequency cutoff $\omega_{min}$ is set by the validity of the calculation leading to Eq. (36), which is estimated as $\omega_{min}\sim(L/N)^{-1}$ Cai:2012zd .

In the analyses of the radiation event rate in the following sections, we shall take $N\Psi_{+}\sim 1$, which implies that neither the number of kinks nor the sharpness of discontinuity would induce an order-of-magnitude change to the radiations from kinks. With this assumption, we can see that the power emitted from cusps always dominates over that of kinks in total. However, since the kinks can emit radiation in a disk-shape region other than a narrow beam for the case of cusps, it is reasonable to expect that the event rate for kinks would be larger than that of cusps by a factor of $1/\theta\sim(\omega L)^{1/3}$ in certain parameter choices. Also we point out that the calculations above hold in the same way for the case of $\mathbf{x}{{}^{\prime}}_{-}$ being discontinuous as well.

III.3 Kink-kink collision

In the situation of kink-kink collisions, two opposite-moving kinks meet each other and thus both $\mathbf{x}{{}^{\prime}}_{+}$ and $\mathbf{x}{{}^{\prime}}_{-}$ are discontinuous. In this case, both $|t_{*}|$ and $|\sigma_{*}|$ are constrained to the order of $\omega^{-1}$, the same as $|t_{*}|$ for kinks. And now $\mathbf{n}$ can be at any direction. The current for the kink-kink collisions becomes

[TABLE]

Accordingly, the power emitted from the kink-kink collisions at the frequency $\omega$ takes the form

[TABLE]

where $\Psi$ denotes $\mathcal{O}(|\Delta\mathbf{x}{{}^{\prime}}_{+}|^{2},|\Delta\mathbf{x}{{}^{\prime}}_{+}||\Delta\mathbf{x}{{}^{\prime}}_{-}|,|\Delta\mathbf{x}{{}^{\prime}}_{-}|^{2})$.

To integrate Eq. (40) over $\omega$ and to assume that there are $N$ left- and right-moving kinks with sharpness $\Psi$, the total power released from the kink-kink collisions is given by

[TABLE]

Under the assumption that $N\Psi\sim 1$, one can easily find that the radiation power from cusps dominate over the other two structures. Eventually, the total electromagnetic radiations from loops of SCSs can be calculated as

[TABLE]

IV Event rate of radio signals

As we have derived the spectrum of electromagnetic radiations emitted from various structures of cosmic strings, it is necessary to translate it to the event rate of observational signals associated with experimental instruments. To carry out this process, we need first to investigate the lifetime and number density of these strings that will be presented in the subsection below.

IV.1 String lifetime and number density

We note that cosmic strings can always emit gravitational radiation as well as electromagnetic radiation. Thus, one can estimate the lifetime of a typical string loop as follows.

The gravitational radiation power takes the form

[TABLE]

where the coefficient takes $\Gamma_{g}\approx 100$ Vachaspati:1984gt . Therefore, the string lifetime can be expressed as follows:

[TABLE]

where $P_{\gamma}\approx P^{c}_{\gamma}$ and $\Gamma={(P_{g}+P^{c}_{\gamma})}/{(G\mu^{2})}.$ As a result, for a string with the initial length $L_{i}$ at the initial moment $t_{i}$, its length would evolve as

[TABLE]

with $t\geqslant t_{i}$.

We are interested in the event rate of radio transients from string loops that exist in the matter dominated phase of the universe. Based on the scaling regime of the evolution model for cosmic strings Vilenkin:2000jqa and taking into account the evolution equation for the string length (45), for loops in the matter area, the number density of string loops can be expressed in terms of the redshift $z$ and the length $L$ as Cai:2012zd

[TABLE]

with

[TABLE]

IV.2 Burst event rate

The burst event rate from string loops of length $L$ at redshift $z$ with $N$ left- and right-moving kinks of typical sharpness $\Psi$ per unit volume can be expressed in a general form by Cai:2012zd

[TABLE]

and we have $(p=0,~{}m=-2/3)$ for cusps, $(p=1,~{}m=-1/3)$ for kinks, and $(p=1,~{}m=0)$ for kink-kink collisions, respectively. The angle $\theta_{\nu_{0}}\sim(\omega L)^{-1/3}$ is the emitted angle of radiations, where $\nu_{0}$ is the observed frequency. Due to the cosmological expansion, the emitted frequency $\nu_{e}$ is related to $\nu_{0}$ by the standard redshift relation through

[TABLE]

The physical volume element $dV(z)$ in the matter era is given by

[TABLE]

and hence, the burst event rate can be written as

[TABLE]

with the prefactor $A\sim 50$, which can be determined numerically.

IV.3 Burst flux and duration

We have denoted the burst event rate as a function of loop length $L$ and redshift $z$. In the perspective of observations, it is necessary to reformulate it in terms of the energy flux per frequency interval $S$ and duration $\Delta$ of the burst. To derive the constraint upon the model parameters, in the following we study the transformation from variables $(L,z)$ to $(S,\Delta)$ with $\mathcal{I}$ and $G\mu$ undetermined, and accordingly, the event rate of radio transients, which originally is a function of the loop length $L$ and the redshift $z$, can then be expressed as a function of the observed energy flux $S$ and the observed duration $\Delta$ with a prefactor depending on $\mathcal{I}$ and $G\mu$. In the model of SCSs, the parameter $\mathcal{I}$ depends not only on the loop length $L$ but also on the energy scale of the phase transition that was expected to occur in the very early universe such that these strings could have been generated. In this regard, it can be treated as a free parameter. Therefore, we are allowed to tune the values of $\mathcal{I}$ and $G\mu$ to accommodate with the FRB observations.

We show how $\Delta$ is related with redshift $z$ and the intrinsic frequency $\nu_{0}$, meanwhile for convenience, we change $d\dot{\mathcal{N}}(L,z)$ from variables $(L,z)$ to $(S,z)$. The observed duration consists of both the observed intrinsic duration $\Delta t$ and the time broadening induced by the cosmological medium scatterings $\Delta t_{s}$:

[TABLE]

The intrinsic duration of bursts $\Delta t_{\rm inert}$ measured at the center-of-mass frame of strings is described by Eq. (26) for cusps and Eq. (33) for kinks or kink-kink collisions. This intrinsic pulse is delayed by cosmic expansion and relativistic Doppler effects. Following Babul:1987lza , the Lorentz factor at cusps is found to be $\gamma\sim(\omega L)^{1/3}$. In the center-of-mass frame of cosmic strings, due to the relativistic motion of cusps, the pulse width $\Delta t_{\rm inert}$ is then given by $\Delta t_{\rm inert}=\gamma\Delta t_{\rm cusp}$. The observed pulse width is associated with $\Delta t_{\rm cusp}$ through $\Delta t_{\rm obs}=\gamma\Delta t_{\rm cusp}$ if the source is moving away from the observer, $\Delta t_{\rm obs}=\Delta t_{\rm cusp}/\gamma$ if the source is moving toward the observer. Meanwhile, the intrinsic frequency $\omega$ is related to the emitted frequency at the observer’s reference frame $\nu_{e}$ as follows. $\nu_{e}=\omega/\gamma$ when the source moves away from the observer, and, $\nu_{e}=\omega\gamma$ when the source moves toward the observer. In addition, the cosmological expansion can yield the redshift effect upon the frequency via $\nu_{o}=\nu_{e}/(1+z)$. By taking into account all the above effects, we can write down the observed intrinsic duration $\Delta t$ to be

[TABLE]

if the cusp moves away from the observer; and in the opposite case, we have

[TABLE]

if the cusp moves towards the observer. We note that the above analysis does not imply that the FRB production from cusps could be related to the motions of strings, but that the observed pulse width of the radiation is dependent of the Lorentz factor due to the fact that the cusps are relativistic. Moreover, the motions of strings would not affect the property that the radiation from cusps is sharply beamed along its forward direction, which can be seen from (29).

When computing the event rate, we need to take both two cases into consideration. In this regard, the ratio between the backward source and forward source can be treated as a free parameter. However, since the astronomical observations are still far from precise data fitting so far, an estimation of this ratio would be of little significance. Therefore, we would like to focus our interest only on the backward sources in the detailed estimation. Additionally, in the case of kinks, due to the extremely short intrinsic duration, even after cosmic expansion and relativistic delay, the observed intrinsic duration is still much smaller than the characteristic time of fast radio bursts, and thus, we will omit the observed intrinsic time duration for kinks.

The observed fast radio bursts are all from extra-galactic sources. Accordingly, the scatterings of these bursts are expected to be contributed by both intergalactic medium(IGM) and interstellar medium(ISM). The time broadening effect due to ISM can be related to the dispersion measures of ISM ${\rm DM}_{\rm ISM}$ using an empirical function as follows Bhat:2004xt :

[TABLE]

where ${\rm DM}_{\rm ISM}$ is often treated as a constant, which roughly equals $95~{}pc/cm^{3}$. The time broadening function due to the IGM scattering is rescaled as three orders of magnitude smaller than that due to ISM (Lorimer:2013roa, ; Caleb:2015uuk, )

[TABLE]

where ${\rm DM}_{\rm IGM}$ is a function of $z$ as follows Ioka:2003fr :

[TABLE]

As a result, the time duration broadening due to scattering can be determined by

[TABLE]

The energy flux per frequency interval $S$ of a string loop observed at a distance $r(z)$ is obtained through dividing the power radiated per unit frequency per unit solid angle by $r(z)^{2}$. The power per unit frequency emitted by cusps, kinks and kink-kink collisions can be found in Eqs. (30), (37), (40), respectively. The solid angles of cusps and kinks are provided in Eqs. (29), (34), respectively. We mention that the expressions of the power are averaged over a frequency interval $T=L/2$, and thus, we need to multiply them by $T$ and then divide them by the burst duration $\Delta$ to get the energy flux. After that, the observed energy flux $S$ for cusps, kinks and kink-kink collisions can then be derived as

[TABLE]

where $(p=0,~{}q=0)$ for cusps; $(p=1,~{}q=2/3)$ for kinks; and $(p=1,~{}q=2)$ for kink-kink collisions, respectively. The proper distance $r(z)$ in the matter dominated flat universe is given by

[TABLE]

Note that, with the help of Eq. (59), we are able to express $S$ as a function of $L$, and thus, can further translate the burst event rate from a function of $(L,z)$ to that of $(S,z)$. The previous estimation of $\Delta$ explicitly shows that $\Delta$ mainly depends on $z$. Therefore, $dLdz=|\partial L/\partial S|dSdz$. In addition, we also comment that this transformation can only work for the cases of cusps and kinks. In the situation of kink-kink collisions, however, as $S$ is also $L$-independent we need to treat this case separately as will be shown later.

For cusps and kinks, we can determine from Eq. (59) that

[TABLE]

where $L$ can be written as

[TABLE]

with $x\equiv\Delta\nu_{0}.$

Eventually we can find

[TABLE]

With this relation, the burst event rate now can be rewritten as

[TABLE]

where

[TABLE]

Now we have the general expression of the burst event rate for cusps and kinks. As we mentioned previously, however, for kink-kink collisions both $S$ and $\Delta$ depend on $z$ only. It implies that $S$ and $\Delta$ are not independent from each other any more. As a result, for kink-kink collisions, we shall estimate the burst event rate as a function of $z$ only. Correspondingly, $S$ and $\Delta$ can be derived from Eqs. (59) and (52). Note that from $S$ and $\Delta$ we cannot determine $L$, and hence we need to integrate the burst event rate $d\dot{\mathcal{N}}(L,z)$ over $L$ to have the total amount $d\dot{\mathcal{N}}(z)$ as follows:

[TABLE]

V Numerical computation and the FRB data fitting

So far we have derived the expressions of the burst event rate coming from SCSs. Accordingly, we now are able to confront the theoretical predictions of SCSs with the observational data of FRB such that model parameters of SCSs can be well constrained. Note that, from the perspective of a theoretical description as performed in previous sections, we have five parameters that characterize the application of SCSs in explaining the FRB data, which are: $\mathcal{I}$, $G\mu$, $\nu_{0}$, $\Delta$, $S$, respectively. The first two of these parameters are associated with the background theory and the latter three are connecting to observations.

Before doing the detailed data fitting of the FRB, we would like to recall that the study on the CMB background anisotropies and the stochastic gravitational wave background has imposed an upper limit on the string tension, i.e., $G\mu<1.3\times 10^{-7}$ Ade:2013xla . Moreover, for SCSs the parameter space of $\mathcal{I}>10^{4}~{}{\rm GeV}$ and $10^{-19}<G\mu<10^{-7}$ has been ruled out due to the constraint on spectral distortions of the CMB photons Kogut:2011xw . Another limit on $\mathcal{I}$ is related with $G\mu$ under the requirement that the dominant process of energy loss is due to the electromagnetic radiations rather than the gravitational radiations; otherwise, the gravitational radiations might have been detected accompanied with the electromagnetic signals observed in astronomical experiments. Accordingly, given $P^{c}_{\gamma}=P_{g}$, the critical current $\mathcal{I}_{*}$ can be determined as follows Cai:2012zd :

[TABLE]

Consequently, for a given $G\mu$, one expects that the value of $\mathcal{I}$ is larger than a critical value $\mathcal{I}_{*}$ derived above to ensure that the major channel for a superconducting string to release energy is via electromagnetic radiations. Using Eqs. (64) and (66), one can see that for each pair of $(G\mu,\mathcal{I})$ the event rate of radio bursts can be expressed as a function of redshift. In this way we are able to fit the event rate function predicted by the theory of SCSs with the normalized observational data. With this purpose, we also need to investigate the threshold flux and the threshold pulse width from the perspective of observations.

Based on the analyses in Refs. Caleb:2015uuk ; Bera:2016qfe , one finds that the signal-to-noise ratio of each event can be determined by

[TABLE]

where $G_{sys}$ is the system gain and approximately equals $0.7\ K/Jy$, $B$ is the bandwidth and approximately equals $340~{}{\rm MHz}$, $N_{pol}$ is the polarization numbers that takes the value $2$, and, $T_{sys}$ is the temperature of the system which is set as $21\ K$. All these values of the system parameters are taken from the setup of the Parkes multi-beam receiver (Caleb:2015uuk, ). Note that $S_{*}$ is the threshold flux for the detection and we choose the value of $S_{*}/N$ to be $10$. The pulse width $\Delta$ can be determined by Eq. (52) derived in the previous section. Additionally, the threshold pulse width required for the detection is set to be $0.1\ ms$.

After having settled down the system restrictions, we are able to calculate $d\dot{\mathcal{N}}(z)/dz$ case by case for cusps, kinks and kink-kink collisions. The frequency band of the Parkes multi-beam receiver is set as $\nu_{0}\sim[1.182,1.522]~{}{\rm GHz}$. We integrate Eq. (64) for $\nu_{0}$ in the Parkes bandwidth, and there is $S\sim[10^{-1},10^{2}]~{}{\rm Jy}$, which is consistent with the flux range of the detected events.

The event rate of radio bursts with $S<10^{-1}{\rm Jy}$ has little contribution to the total amount, since that their flux is lower than the threshold flux $S_{*}$. In addition, the event rate with $S>10^{2}~{}{\rm Jy}$ is also much suppressed and has almost no contribution to the total event rate. This is because the upper bound of $S$ is quite well restrained by the length of string loops and the threshold observational pulse width. Eq. (59) shows that, for both cusps and kinks, $S$ has a power-law correlation with $L$ and an inverse relation with $\Delta$. From Eq. (45), we can find an upper bound of the loop length: $L_{i}-\Gamma G\mu(t-t_{eq})$, where $L_{i}$ is the string length at some initial moment of the cosmic evolution, and $t$ is the cosmic time related to the redshift where the observed string would be located. Since $L_{i}$ is much larger than $\Gamma G\mu t_{0}$, this restriction upon the string length is often negligible. For cusps, there is another limit from Eq. (31) that requires $L<{\mu^{3/2}}/{\mathcal{I}^{3}\omega_{max}}$, and thus a radio signal from the cusp has to satisfy this inequality. For kink-kink collisions, we use Eq. (66) to calculate the event rate straightforwardly as there is no further constraint of the loop length in this case. However, we will show later that the flux from kink-kink collisions is typically much smaller than the threshold flux $S_{*}$ required by the present detection sensitivity. Therefore, we conclude that there is very little chance to observe such an occurrence by the present detection ability in radio astronomy.

Recall that the FRB parameters in association with the cosmic string model are: the derived redshift $z$, the dispersion measures DM (unit: $\rm{cm^{-3}pc}$), the pulse width $\rm{W_{obs}}$ (unit: ms), the energy flux $\rm{S_{obs}}$ or $\rm{F_{obs}}$ (unit: Jy or Jy ms), respectively. In Table. 1 we list 15 FRB events detected by the Parkes multi-beam receiver Petroff:2016tcr . Among them, the redshift of the furthest source is about $z=1.3$. For convenience of numerical analysis, we separate the regime of $z\sim[0,1.4]$ into 7 bins and count the event number $\Delta N_{obs}/\Delta z_{bin}$ in each bin. The error bars of observational data are set as ${\rm e}_{\rm obs}=\sqrt{\Delta N_{obs}/\Delta z_{bin}}$. We then normalize the data so that the sum of 7 bins satisfies

[TABLE]

where ${\rm n}_{\rm obs}$ stands for the normalized event rate per redshift bin. To compare theoretical prediction with it, for each $(\mu,\mathcal{I})$, we normalize the $d\dot{\mathcal{N}}(z)$ function so that the normalized theoretical event rate per redshift ${\rm n}_{\rm th}$ satisfies $\int{\rm n}_{\rm th}dz=\int(d\dot{\mathcal{N}}(z)/dz)dz=1$. We use the $\chi^{2}$ value to denote the quality of fitting between these two normalized data settings, where $\chi^{2}$ is defined as

[TABLE]

Afterwards, we compute the value of $\chi^{2}$ for $(\mu,\mathcal{I})$ pairs that fit the above settings, with the viable result being presented in the upper panel of Fig. 1. The black dashed line in Fig. 1 denotes the boundary where $\mathcal{I}$ equals the critical value at which the gravitational radiation is comparable with electromagnetic radiation. We are interested in the parameter space dominated by electromagnetic radiation that is below the black dashed line. In this area, the light blue contour corresponds to the regime with $\chi^{2}<10$. If we further decrease the value of $\chi^{2}$ to be $\chi^{2}<6$, we interestingly find that there exist two separate regimes allowed by the FRB data, which are shown as the dark blue contours in the upper panel of Fig. 1. Here, we would like to emphasize that, due to the very limit sample of the FRB data, the statistics of the data fitting remains very poor and thus the confidence level is quite low. According to our numerical analysis, the significance is approximately $1.3\sigma$ for the contour with $\chi^{2}<10$. From the numerical estimation, we observe that the allowed range is roughly within $G\mu\sim[10^{-14},10^{-12}]$, and, $\mathcal{I}\sim[10^{-1},10^{3}]~{}{\rm GeV}$. Given the observationally allowed parameter space as shown in the upper panel of Fig. 1, we numerically derive the corresponding distribution of the burst event rate per year in the lower panel of the same figure. Accordingly, we find that the event rate $\dot{\mathcal{N}}$ within the colored region is between $10^{5}\sim 10^{9}~{}{\rm yr}^{-1}$. This result implies that the observation of the radio bursts from SCSs is promising in various experiments of radio astronomy.

Note that best fitting contours could lie on two sides of the parameter space as observed in the upper panel of Fig 1. For each contour, one can derive a best fitting spot in the parameter space where $\chi^{2}$ is the lowest. We numerically investigate the event rate as a function of the redshift for both two “best fitting” spots as shown in Fig. 2. In this figure, the model parameters are determined as follows. For the left panel, we have

[TABLE]

with $\chi^{2}=4.72$ and the corresponding event rate is found to be $\dot{\mathcal{N}}\sim 2.3\times 10^{8}$ ${\rm yr}^{-1}$; moreover, for the right panel, the model parameters take the values

[TABLE]

with $\chi^{2}=3.88$ and the event rate takes $\dot{\mathcal{N}}\sim 1.6\times 10^{7}$ ${\rm yr}^{-1}$.

From the numerical estimations presented in Fig. 2, we can see that the event rate of the above two cases are slightly higher than the estimation in other work, namely, $\dot{\mathcal{N}}\sim 10^{6}~{}{\rm yr}^{-1}$ as was estimated in Ref. Caleb:2015uuk . For each case of our investigation, in the upper panel, the solid orange curves show the event rate emitted from the cusps and the solid green curves depict the event rate arisen from the kinks. There is no signal from the kink-kink collisions shown in the figure for the reason that either their flux is much lower than the threshold of the experimental detection or their event rate is far below detectability. In the lower panel, the blue dashed line denotes the normalized total event rate, and, the red points are normalized observational data. So far, the observed events exhibit a peak value around the redshift $z=0.5$. In the case shown in the left panel, radio bursts due to string cusps make dominant contribution throughout the whole redshift range $z\sim(0,1.6)$; while, in the right case, radio burst signals from kinks can become dominant around the redshift $z=0.5\sim 0.7$, with the corresponding amplitude being about ten times larger than that from cusps. This interesting local effect happens to overlap with the bulging profile of the observational data as shown in the lower right panel. Therefore, if the local feature of high event rate of radio bursts within the range of $z=0.5\sim 0.7$ were confirmed in the accumulated observations with high precision, the theory of SCSs could provide an interesting explanation of this possible phenomenon.

For both two cases, the burst event rate from cusps exhibits an abrupt cutoff at the high redshift regime due to the limitation of the string length. In our detailed calculation, we have constrained the frequency range to be within $1.182\sim 1.522$ ${\rm GHz}$ in order to be consistent with the experimental parameter of Parkes. As the frequency increases, we see that the cutoff redshift would also increase. For instance, $z\sim 2$ would correspond to $\nu_{0}\sim 3~{}{\rm GHz}$. In the low redshift regime, the burst event rate decreases mainly due to the fact that the pulse width is too small. Moreover, we find that the observed pulse width shows very little dependence on the redshift. However, according to the scattering model analyzed in Sec. IV.3, one expects that the pulse width broadening due to the ISM and IGM effects should be correlated with the redshift. We argue that this issue may be resolved either by a fine tuning of the background theory of SCSs or by a better version of the scattering model.

VI Conclusion and discussion

To summary, we in this paper have applied the observational data of FRB to constrain model parameters appeared in the theory of the superconducting cosmic strings. In particular, based on the method of the order-of-magnitude estimation, we have recomputed the electromagnetic emission power of three kinds of string structures, which are cusps, kinks and kink-kink collisions, respectively. After that, we have demonstrated in Fig. 1 that the SCSs models within the parameter space with $G\mu\sim[10^{-14},10^{-12}]$ and $\mathcal{I}\sim[10^{-1},10^{3}]~{}{\rm GeV}$ are consistent with the present FRB observations. By performing detailed numerical calculations, we surprisingly found two possible best fitting spots as estimated by the $\chi^{2}$ value. One best fitting spot corresponds to the string model with parameters being $G\mu=5\times 10^{-15}$ and $\mathcal{I}=0.313~{}\rm{GeV}$ and in this case the event rates of the radio bursts are mostly contributed by cusps; the other one corresponds to the model with $G\mu=3\times 10^{-13}$ and $\mathcal{I}=236.28~{}\rm{GeV}$, in which the contribution of kinks could dominate over that of cusps in the intermediate regime along the redshift. With accumulated high precision observations of radio telescopes in near future, it will be very promising to examine or to rule out the theoretical possibility of using SCSs as sources of cosmological FRB.

Recall that we have pointed out in Sec. IV.3 that radio bursts can be generated by cusps moving either backwards or towards the observers. Radio bursts from those forward sources can have larger value of the event rate than the backward sources over certain redshift by an order of magnitude. Such differences become small as $G\mu$ increases. It is then natural to consider the phenomenon that the ratio of the backward sources over forward ones would also affect the analyses performed in the present study. As future observations will provide us a more reliable relationship between the event rate and the redshift, it will also be important to take into account this ratio parameter in the follow-up study.

In addition, among all the detected FRB events, it has been reported that the event FRB 150418 and FRB 150807 are linearly polarized Keane:2016yyk ; Ravi:2016kfj , and FRB 140514 and FRB 110523 have shown the mixture of linear and circular polarizations Petroff:2014taa ; Masui:2015kmb . We note that the circular polarization could be caused by the foreground contamination, such as, the Faraday rotation effect when the radio signals pass through the magnetized region and/or the high-density medium in Milky Way as well as intergalactic environments, or even possibly the host galaxy. According to the study of Petroff:2014taa , the event FRB 140514 has shown that the linear polarization is intrinsic to the signal. In the literature there have been attempts in understanding these particular patterns of polarizations Pen:2015ema . Importantly, the study of SCSs in (Cai:2011bi, ) has already demonstrated that radiations from string cusps should be linearly polarized while no specific prediction on the polarization was made for kinks. In this regard, the precise measurement of the circular polarization can help to identify the host environment of FRB sources and the theoretical origin of SCSs. Although these environmental issues remain unclear, so far the present FRB observations yet cannot rule out the models of SCSs. Thus, it deserves to study in detail the potential co-evolution of SCSs and galaxy in the future in order to examine possible connections among the polarization of radio signals, the models of SCSs and the unknown environmental effects.

We end by noting that our present results are also consistent with the present constraints on the parameters of the SCSs models, for instance, as obtained in Miyamoto:2012ck ; Brandenberger:2017uwo . In fact, our analysis yields more stringent constraints with the starting point that SCSs could play the role of the theoretical origin of cosmological FRB. With this theoretical possibility, the current parameter of the strings could be as low as the GeV scale, and thus, might become of observational interest in high energy phenomenology as well as experiments, such as, the Large Hadron Collider, the International Linear Collider, the Circular Electron Positron Collider and so on.

Acknowledgments

We are grateful to Hong Guo, Eray Sabancilar, Daniele A. Steer, Tanmay Vachaspati and Tinggui Wang for valuable communications. We especially thank Zhi-Qiang Shen for insightful discussion on the frontier of radio astronomy. J.Y. and K.W. are supported in part by the NSFC Fund for Fostering Talents in Basic Science (No. J1310021). Y.F.C. is supported in part by the Chinese National Youth Thousand Talents Program, by the CAST Young Elite Scientists Sponsorship Program (2016QNRC001), by the NSFC (Nos. 11421303, 11653002), and by the Fundamental Research Funds for the Central Universities. Part of numerical simulations are operated on the computer cluster LINDA in the particle cosmology group at USTC.