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

Separating the superradiant emission from the Hawking radiation from a rotating black hole

De-Chang Dai1,2111communicating author: De-Chang Dai,
email: diedachung@gmail.com
1 Department of Physics, National Dong Hwa University, Hualien, Taiwan, Republic of China 2 CERCA, Department of Physics, Case Western Reserve University, Cleveland OH 44106-7079
   Dejan Stojkovic3 3 HEPCOS, Department of Physics, SUNY at Buffalo, Buffalo, NY 14260-1500, USA
Abstract

Emission of particles created in the background of a rotating black hole can be greatly amplified taking away rotational energy of a black hole. This amplification affects both particles created near the horizon (due to the Hawing effect), and particles created near the potential barrier far from the horizon. Only the latter effect is called the superradiance in the strict sense. We explicitly calculate the superradiant emission for scalar particles and compare it with the total scalar particle emission (Hawking radiation plus superradiance) to clarify some confusion in the literature. We clearly show that these two emissions are not the same. In particular, superradiance persists even for extremal black holes whose Hawking temperature is zero.

I Superradiance

The notion of superrdiance, or superradiant emission of particles, has been used in a wide range of situations in the literature. The term ”superradiance” was introduced by R.H. Dicke in 1954 [1], describing an effect in which disordered energy is converted into coherent energy. In classical physics, superradiance is a classical phenomenon in which an amplitude of an outgoing wave after the reflection is greater than the amplitude of the ingoing wave [2]. This phenomenon can happen in the background of a rotating black hole [3], which contains an ergosphere, i.e. the region between the infinite redshift surface and the event horizon. In such a background, an incident wave can take away some of the rotational energy of the black hole and get amplified after reflection, thus effectively yielding a reflection coefficient greater than one (i.e. negative absorption coefficient).

In the context of quantum Hawking radiation from a black hole, it was noticed in numerical studies that spontaneous emission from a black hole can also get amplified taking away rotational energy of the black hole [4]. Calculations of the black hole greybody factors indicate that this amplification is very much spin dependent, with emission of higher spin particles strongly favored. In [5] an analytic explanation of this phenomenon was given in terms of the spin-spin interaction between the spin of the rotating black hole and spin of the emitted particle.

While the spin dependent amplification of Hawking radiation is also often called superradiance in the spirit of Dicke’s definition in [1], it is different from the original superradiance in [2] or [3] which crucially rely on the negative absorption coefficient. For example the superradiance as defined in [2] or [3] is not possible for fermions [3, 6]. If the incident wave is made of fermions, then the reflected wave can not get amplified due to Pauli exclusion principle, since all the available states are already occupied. In this paper, for clarity, we will call this effect the superradiance and separate it from the Hawking effect. The crucial difference is that the Hawking effect happens in the presence of the horizon, while the superradiance does not need the horizon. Superradiant emission is simply the effect of particle creation in scattering from the potential barrier. Since the black hole contains both the horizon and the potential barrier outside the horizon, the total radiation from the rotating black hole will include both the Hawking effect and superradiance.

The effects of superradiance in the context of black holes have been extensively explored in the literature (for a comprehensive review see [7]). However, to the best of our knowledge, it appears that the superradiance has never been explicitly separated from the Hawking radiation in concrete calculations. Therefore, in this paper we will explicitly calculate the particle production due to the potential barrier away from the horizon of a rotating black hole and compare it with the total radiation. Among the other things, it will become clear that the superradiance exists even when the Hawking temperature drops to zero (i.e. for an extremal balk hole). Thus, treating a black hole as a black body emitter with a finite temperature TT, and thus intensity of radiation proportional to T4T^{4}, is a significant oversimplification. This fact can potentially have some implications even for the information loss paradox. Since suprerradiance particles are created at the barrier outside of the horizon, they should not affect the horizon and singularity physics. However, if one takes the whole process of the black hole creation and evaporation into account, potential barrier outside of the horizon still retains some information about the gravitational process that created the black hole. This implies that superradiance is relevant to black hole information[8]. This is especially important in the context of the black hole scalar hair.

We start with the metric for a rotating black hole. The geometry of a rotating black hole is described by the Kerr metric in Boyer-Lindquist coordinates

ds2\displaystyle ds^{2} =\displaystyle= (12MrΣ)dt24Mrasin2θΣdtdϕ+ΣΔdr2+Σdθ2\displaystyle-(1-\frac{2Mr}{\Sigma})dt^{2}-\frac{4Mra\sin^{2}\theta}{\Sigma}dtd\phi+\frac{\Sigma}{\Delta}dr^{2}+\Sigma d\theta^{2} (3)
+(r2+a2+2Mra2sin2θΣsin2θ)dϕ2,\displaystyle+\Big{(}r^{2}+a^{2}+\frac{2Mra^{2}\sin^{2}\theta}{\Sigma}\sin^{2}\theta\Big{)}d\phi^{2},
Δ=r22Mr+a2,\displaystyle\Delta=r^{2}-2Mr+a^{2},
Σ=r2+a2cos2θ,\displaystyle\Sigma=r^{2}+a^{2}\cos^{2}\theta,

where aa is the black hole rotation parameter, while MM is the mass of the black hole.

We now place a scalar field in this background. We will follow the Frolov’s book [11] on black hole physics. There, a complete classification of different types of bases in the black hole Penrose diagram were given (see Fig. 1). In this setup one can separate the Hawking effect from superradiance in physically and mathematically clear way. We thus decompose a scalar field ψ\psi in the spherically symmetric coordinates

ψl,m=eiωtRl,m(r,ω)Sl,m(θ,ω)eimϕ2π,\psi_{l,m}=e^{-i\omega t}R_{l,m}(r,\omega)\frac{S_{l,m}(\theta,\omega)e^{im\phi}}{\sqrt{2\pi}}, (4)

where ll and mm are angular and magnetic quantum number respectively. The equation of motion can be separated into two main equations,

ddr(ΔdRl,mdr)+(K2Δλ)Rl,m=0.\displaystyle\frac{d}{dr}\Big{(}\Delta\frac{dR_{l,m}}{dr}\Big{)}+\Big{(}\frac{K^{2}}{\Delta}-\lambda\Big{)}R_{l,m}=0. (5)
1sinθddθ(sinθdSl,mdθ)+(a2ω2cos2θm2sin2θ\displaystyle\frac{1}{\sin\theta}\frac{d}{d\theta}\Big{(}\sin\theta\frac{dS_{l,m}}{d\theta}\Big{)}+\Big{(}a^{2}\omega^{2}\cos^{2}\theta-\frac{m^{2}}{\sin^{2}\theta}
+E)Sl,m=0,\displaystyle+E\Big{)}S_{l,m}=0, (6)
K=(r2+a2)ωam,\displaystyle K=(r^{2}+a^{2})\omega-am, (7)
λ=E+a2ω22amω,\displaystyle\lambda=E+a^{2}\omega^{2}-2am\omega, (8)

here, EE is the eigenvalue of the angular equation, Eq. (6), while λ\lambda is the corresponding eigenvalue of the radial equation, Eq. (5). The radial equation can be simplified by introducing a new variable χ\chi,

χ=(r2+a2)1/2Rl,m.\chi=(r^{2}+a^{2})^{1/2}R_{l,m}. (9)

The radial part becomes

(d2dr2+K2λΔ(r2+a2)2G2dGdr)χ=0,\Big{(}\frac{d^{2}}{dr_{*}^{2}}+\frac{K^{2}-\lambda\Delta}{(r^{2}+a^{2})^{2}}-G^{2}-\frac{dG}{dr_{*}}\Big{)}\chi=0, (10)

where, rr^{*} is the tortoise coordinate

dr=r2+a2Δdr.dr_{*}=\frac{r^{2}+a^{2}}{\Delta}dr. (11)

and

G=rΔ(r2+a2)2G=\frac{r\Delta}{(r^{2}+a^{2})^{2}} (12)

The asymptotic solution near infinity is

χe±iωr,r.\chi\sim e^{\pm i\omega r_{*}},r\rightarrow\infty. (13)

The asymptotic solution near the horizon is

χe±i(ωmΩH)r,rr+,\chi\sim e^{\pm i(\omega-m\Omega_{H})r_{*}},r\rightarrow r_{+}, (14)

where, ΩH=a/(2Mr+)\Omega_{H}=a/(2Mr_{+}) and r+r_{+} is the radius of the black hole outer horizon.

In the Penrose diagram shown in fig.1, an incident wave coming from the past infinity, II^{-}, is scattered by the potential. Part of the wave is reflected to infinity, while part of it will penetrate the barrier and fall into the black hole horizon.

χ={eiωr+R~eiωr, for rT~e±i(ωmΩH)r, for rr+.\chi=\begin{cases}e^{-i\omega r_{*}}+\tilde{R}e^{i\omega r_{*}}&\text{, for }r\rightarrow\infty\\ \tilde{T}e^{\pm i(\omega-m\Omega_{H})r_{*}}&\text{, for }r\rightarrow r_{+}.\end{cases} (15)

The radial equation, Eq. (10), is a real second order ordinary differential equation which satisfies the Wronskian relation. From the Wronkian relation, one can prove that

(1mΩHω)|T~|2=1|R~|2.(1-\frac{m\Omega_{H}}{\omega})|\tilde{T}|^{2}=1-|\tilde{R}|^{2}. (16)

From here we see that if ω<mΩH\omega<m\Omega_{H}, then |R~|2>1|\tilde{R}|^{2}>1. This means that the reflected wave has a greater amplitude than the initial wave, which in turn implies that particles are created during the scattering. These new particles are not directly created by the black hole horizon, and we focus our attention to this process.

Refer to caption
Figure 1: The Penrose-Carter diagram: The space outside the horizon is presented in the right square. The upper triangle represents the space inside the future horizon. The black hole radiation due to Hawking effect involves these two regions. On the other hand, only the right square is involved in particle creation by the superradiance mechanism. The thin dashed line represents the potential barrier that induces the superradiance. Four types of bases (downdown, upup, inin and outout) are involved in the process. Five bases (including dndn) are involved in the full black hole radiation [11].

II Particle creation by the superradiance mechanism

As we already explained, the superradiance mechanism operates outside the horizon represented by the line labeled by H+ in fig.1. So all the relevant events are in the right square in fig. 1. In contrast, Hawking radiation is induced by the black hole horizon.

To study particle creation by superradiance mechanism we have to identify the basis in which we decompose the fields, and the vacuum state of the field. To describe a vacuum state of a field, we need at least two bases. For our purpose, we define four possible bases [11]. The first is the inin-coming mode. It represents a wave which goes from the past null infinity, II^{-}, to the black hole,

χJin1ωexp(iωr).\chi^{in}_{J}\sim\frac{1}{\sqrt{\omega}}\exp(-i\omega r_{*}). (17)

The second one is the outout-going mode. It represents a wave which propagates from the black hole to future null infinity, I+I^{+},

χJout1ωexp(iωr)\chi^{out}_{J}\sim\frac{1}{\sqrt{\omega}}\exp(i\omega r_{*}) (18)

The third one is the downdown mode. It represents a wave which goes into the future horizon, H+H^{+},

χJdown1|ωmΩH|exp(i(ωmΩH)r).\chi^{down}_{J}\sim\frac{1}{\sqrt{|\omega-m\Omega_{H}|}}\exp(-i(\omega-m\Omega_{H})r_{*}). (19)

The forth one is the upup mode. It represents a wave which goes away from the past horizon, HH^{-},

χHup1|ωmΩH|exp(i(ωmΩH)r).\chi^{up}_{H}\sim\frac{1}{\sqrt{|\omega-m\Omega_{H}|}}\exp(i(\omega-m\Omega_{H})r_{*}). (20)

Since a field decomposition requires two distinct bases, we can decompose χ\chi in two different ways

χ^\displaystyle\hat{\chi} =\displaystyle= Ja^JinχJin+a^Jupχ~Jup+h.c.\displaystyle\sum_{J}\hat{a}_{J}^{in}\chi^{in}_{J}+\hat{a}_{J}^{up}\tilde{\chi}^{up}_{J}+h.c. (21)
=\displaystyle= Jb^JoutχJout+b^Jdownχ~Jdown+h.c.\displaystyle\sum_{J}\hat{b}_{J}^{out}\chi^{out}_{J}+\hat{b}_{J}^{down}\tilde{\chi}^{down}_{J}+h.c. (22)

where, J={ω,l,m}J=\{\omega,l,m\}, while h.c.h.c. stands for the hermitian conjugate terms. We also have

χ~Jα\displaystyle\tilde{\chi}^{\alpha}_{J} =\displaystyle= χJα, if ωmΩH>0\displaystyle\chi^{\alpha}_{J}\text{, if }\omega-m\Omega_{H}>0 (23)
=\displaystyle= χJα, if ωmΩH<0,\displaystyle\chi^{\alpha}_{J}{}^{*}\text{, if }\omega-m\Omega_{H}<0, (24)

where α\alpha can be upup or downdown. The two types of vacuum corresponding to a^Jα\hat{a}_{J}^{\alpha} and b^Jα\hat{b}_{J}^{\alpha} are

a^Jα|in;0\displaystyle\hat{a}_{J}^{\alpha}\ket{in;0} =\displaystyle= 0\displaystyle 0 (25)
b^Jα|out;0\displaystyle\hat{b}_{J}^{\alpha}\ket{out;0} =\displaystyle= 0.\displaystyle 0. (26)

The inin mode is expressed in terms of the bases in the past, while the outout mode is expressed in terms of the bases in the future. Consider now a field which starts from vacuum in the far past, and after evolving is seen in terms of the outout bases

χJin\displaystyle\chi^{in}_{J} \displaystyle\rightarrow RJχJout+TJχJdown\displaystyle R_{J}\chi^{out}_{J}+T_{J}\chi^{down}_{J} (27)
χJup\displaystyle\chi^{up}_{J} \displaystyle\rightarrow tJχJout+rJχJdown.\displaystyle t_{J}\chi^{out}_{J}+r_{J}\chi^{down}_{J}. (28)

If we compare eqs. (13) and (14) with eqs. (17) to (20), we see that there are extra normalization factors in the latter four equations. This will lead to extra factors in the transmission coefficients. Therefore, TJT_{J} and T~\tilde{T} (eq. 15) are related by TJ=|1mΩH/ω|T~T_{J}=\sqrt{|1-m\Omega_{H}/\omega|}\tilde{T}. The creation and annihilation operators are related according to the relationship between the modes. For ωmΩH>0\omega-m\Omega_{H}>0, we have

b^Jout=RJa^Jin+tJa^Jup.\hat{b}^{out}_{J}=R_{J}\hat{a}^{in}_{J}+t_{J}\hat{a}^{up}_{J}. (29)

In this case, there is no particle creation since there is no mixing of the creation and annihilation operators. For ωmΩH<0\omega-m\Omega_{H}<0, we have

b^Jout=RJa^Jin+tJa^Jup.\hat{b}^{out}_{J}=R_{J}\hat{a}^{in}_{J}+t_{J}\hat{a}^{up}_{J}{}^{\dagger}. (30)

In this case there is particle creation because of the mixing of the creation and annihilation operators. Thus, ωmΩH<0\omega-m\Omega_{H}<0 is the necessary condition for the superradiance. From the commutation relations we can get

|tJ|2=|TJ|2.|t_{J}|^{2}=|T_{J}|^{2}. (31)

From here we can calculate the particle creation number due to the superradiance effect as

nJ=in,0|bJoutbJout|in,0=|tJ|2, if ωmΩH<0.n_{J}=\bra{in,0}b_{J}^{out}{}^{\dagger}b_{J}^{out}\ket{in,0}=|t_{J}|^{2}\text{, if }\omega-m\Omega_{H}<0. (32)

To compare the superradiance particle creation with the Hawking effect and demonstrate their difference, we calculate the total particle creation number (Hawking effect plus superradiance) characterized by the transmission coefficient |TJ|2|T_{J}|^{2}[6]

nJT=sign(1mΩω)|TJ|2exp(ωmΩHT)1,n_{J}^{T}=\frac{\text{sign}(1-\frac{m\Omega}{\omega})|T_{J}|^{2}}{\exp(\frac{\omega-m\Omega_{H}}{T})-1}, (33)

where, TT is the black hole temperature.

We can immediately see the fundamental difference between the superradiance and Hawking effect. For example, for an extremal black hole the Hawking temperature goes to zero, T0T\rightarrow 0. If ω>mΩH\omega>m\Omega_{H}, which is outside of the superradiant regime, nJT=0n_{J}^{T}=0 since both the Hawking effect and superradiance are absent. However, if ω<mΩH\omega<m\Omega_{H}, we get nJ=nJTn_{J}=n_{J}^{T}. In that case the Hawking effect is still absent (Hawking temperature is still zero), however the superradiance is present and it is the only contribution to the total radiation from a black hole.

As an illustration, we will now compare the power spectrum, dEdtdω\frac{dE}{dtd\omega}, and power, dEdt\frac{dE}{dt}, between the superradiance and total radiation. The power emitted in particles generated by the superradiance mechanism is

PsdEsdt=12πl,mnJωdω.P_{s}\equiv\frac{dE_{s}}{dt}=\frac{1}{2\pi}\int\sum_{l,m}n_{J}\omega d\omega. (34)

The total power of particles coming out of the black hole is

PTdETdt=12πl,mnJTωdω.P_{T}\equiv\frac{dE_{T}}{dt}=\frac{1}{2\pi}\int\sum_{l,m}n_{J}^{T}\omega d\omega. (35)

The power spectrum of emitted particles is defined as the emitted energy per unit time per unit frequency, i.e. dEdtdω\frac{dE}{dtd\omega}. We use the units where G=c=k=1G=c=k=1. In fig. 2, we plot the comparison between the power spectrum of particles created by the superradiance mechanism and the power spectrum of the total radiation from the black hole. For convenience, the rotation parameter aa is rescaled to a=a/Ma^{*}=a/M. We fix the value of the black hole rotation parameter to a=0.9a^{*}=0.9 where the superradiance of scalar particles is significant, but still does not dominate over the Hawking emission.

Refer to caption
Figure 2: Comparison of the scalar particle power spectra, dEdtdω\frac{dE}{dtd\omega}, for a=0.9a^{*}=0.9. The black curve represents the power spectrum for the total radiation of scalar particles from the black hole (Hawking effect plus superradiance). The dashed line represents the power spectrum for the scalar particles created by superradiance. The transmission coefficient, TJT_{J}, is taken from the BlackHawk generator [9, 10]. The spectra are clearly different. We also see that for this value of aa^{*} superradiance is significant but still does not dominate over the Hawking emission.

In fig.3, we plot the comparison of the powers as a function of the black hole rotation parameter aa^{*}. We can see that the superradiance becomes dominant for highly rotating black holes at a0.94a^{*}\approx 0.94.

Refer to caption
Figure 3: Comparison of the powers for the scalar particles, P=dEdtP=\frac{dE}{dt}, as a function of the black hole rotation parameter aa^{*}. The black solid curve is the total scalar particle radiation from the black hole (Hawking effect plus superradiance), PTP_{T}. The thin solid curve is one half of the total radiation, plotted for convenience. The dashed line is the power from the superradiance mechanism, PsP_{s}. We can see that the superradiance does not play important role for a<0.7a^{*}<0.7. However, it becomes dominant around a=0.94a^{*}=0.94.

In fig.4, we plot the ratio between the powers emitted by superradiance and total radiation. Once again we see that for the extremal black hole (i.e. a=1a^{*}=1) superradiance equals the total radiation since the Hawking effect ceases to exists at that point.

Refer to caption
Figure 4: Ratio between the powers emitted by superradiance and total radiation, Ps/PTP_{s}/P_{T}, for the scalar particles. The curve represents the ratio between the powers of scalar particles created by the superradiance mechanism and the total radiation (Hawking effect plus superradiance). We see that at a=1a^{*}=1 only superradiance contributes to the total radiation.

III Conclusion

The main point in this paper is to make a pedagogical distinction between the amplification of the Hawking emission in the background of a rotating black hole from the effect of superradiance. Both the spin dependent amplification of Hawking radiation and superradiance crucially depend on taking away rotational energy from a rotating black hole. However, superradiance is created at the potential barrier away from the horizon and should not be mixed with Hawking radiation. We explicitly calculated the superradiant emission and compared it with the total radiation. We showed that the superradiance of scalar particles is negligible for a<0.7a^{*}<0.7, but it becomes dominant around a=0.94a^{*}=0.94. The numerical values might change if one includes all types of particles since the superradiant amplification is stronger for higher spin particles. For the modes that do not satisfy the condition for superradiance, i.e. for ω>mΩH\omega>m\Omega_{H}, we have only the Hawking effect which clearly vanishes for the extremal black hole when the Hawking temperature goes to zero. However, for the superradaint modes, with ω<mΩH\omega<m\Omega_{H}, emission persists even for the extremal black hole. In this paper we considered only a scalar field, but similar calculations can be performed for the vector and graviton fields. As we mentioned in the introduction, fermions do not exhibit superradiance because of the Pauli exclusion principle.

The fact that even the extremal rotating black holes emit particles might be relevant in various cosmological scenarios that strongly depend on the existence or non-existence of such radiation (e.g. [12, 13, 14, 15]). In particular, small primordial black holes can Hawking radiate gravitons, thus contributing to the primordial stochastic gravitational wave background [12]. In our context, even the extremal rotating black holes with the zero Hawking temperature can keep contributing to this background. On the other hand, extremal primordial black holes are often mentioned as good dark matter candidates due to their lack of Hawking radiation [13, 14, 15]. Again, in our context we see that the extremal black holes will keep emitting particles, which makes them unstable and also visible.

Finally, we comment on the likelihood that small primordial black holes have large angular momentum. Density perturbations usually do not carry significant angular momentum, so the corresponding black holes created by this mechanism will not be spinning fast. However, if two black holes merge, their initial relative angular momentum gets transformed into the final black hole spin, due to the angular momentum conservation. In addition, accretion of surrounding material is very efficient in spinning the black hole up. If the last 2050%20-50\% of the black hole mass came from accretion, such a black hole would be close to extremal. Finally, black holes formed in collisions of the energetic particles in the early universe [16] are also expected to carry high spin due to the initial relative angular momentum of the colliding particles.

Acknowledgements.
D.C. Dai is supported by the National Science and Technology Council (under grant no. 111-2112-M-259-016-MY3) D.S. is partially supported by the US National Science Foundation, under Grant No. PHY-2014021.

References

  • [1] R. H. Dicke, Phys. Rev.  93, 1, 99, (1954)
  • [2] Ya. B. Zeldovich, Pisma Zh. Eksp. Teor. Fiz. 14, 270 (1971) [JETP Lett. 14, 180 (1971)]; Zh. Eksp. Teor. Fiz 62, 2076 (1972) [Sov. Phys. JETP 35, 1085 (1972)].
  • [3] A.A. Starobinsky, Sov. Phys. - JETP 37, 28 (1973); A.A. Starobinsky and S.M. Churilov, Sov. Phys. - JETP, 38, 1 (1973); W. Unruh, Phys. Rev. Lett., 31, 1265 (1973)
  • [4] D. C. Dai, G. Starkman, D. Stojkovic, C. Issever, E. Rizvi and J. Tseng, Phys. Rev. D 77, 076007 (2008) doi:10.1103/PhysRevD.77.076007 [arXiv:0711.3012 [hep-ph]].
  • [5] D. Dai and D. Stojkovic, JHEP 08, 016 (2010) doi:10.1007/JHEP08(2010)016 [arXiv:1008.4586 [gr-qc]].
  • [6] V. Frolov and I. Novikov ”Black Hole Physics: Basic Concepts and New Developments”, Springer (1998)
  • [7] R. Brito, V. Cardoso and P. Pani, Physics,” Lect. Notes Phys. 906, pp.1-237 (2015) 2020, ISBN 978-3-319-18999-4, 978-3-319-19000-6, 978-3-030-46621-3, 978-3-030-46622-0 doi:10.1007/978-3-319-19000-6 [arXiv:1501.06570 [gr-qc]].
  • [8] I. Y. Park, PTEP 2021, no.6, 6 (2021) doi:10.1093/ptep/ptab045 [arXiv:1912.07413 [hep-th]].
  • [9] A. Arbey and J. Auffinger, Eur. Phys. J. C79 (2019) 693, arXiv:1905.04268 [gr-qc]
  • [10] A. Arbey and J. Auffinger, Eur. Phys. J. C81 (2021) 910, arXiv:2108.02737 [gr-qc]
  • [11] Valeri P. Frolov , Igor D. Novikov, ”Black Hole Physics: Basic Concepts and New Developments”, by Springer (1998) . ISBN 0792351452 (ISBN13: 9780792351450)
  • [12] R. Dong, W. H. Kinney and D. Stojkovic, JCAP 10, 034 (2016) doi:10.1088/1475-7516/2016/10/034 [arXiv:1511.05642 [astro-ph.CO]].
  • [13] J. A. de Freitas Pacheco and J. Silk, Phys. Rev. D 101, no.8, 083022 (2020) doi:10.1103/PhysRevD.101.083022 [arXiv:2003.12072 [astro-ph.CO]].
  • [14] S. Chongchitnan and J. Silk, Phys. Rev. D 104, no.8, 083018 (2021) doi:10.1103/PhysRevD.104.083018 [arXiv:2109.12268 [astro-ph.CO]].
  • [15] J. A. de Freitas Pacheco, E. Kiritsis, M. Lucca and J. Silk, [arXiv:2301.13215 [astro-ph.CO]].
  • [16] A. Saini and D. Stojkovic, JCAP 05, 071 (2018) doi:10.1088/1475-7516/2018/05/071 [arXiv:1711.06732 [gr-qc]].