Evolution of scalar field resonances in a braneworld

In this work, we investigate the numerical evolution of massive Kaluza–Klein (KK) modes of a scalar field in a thick brane. We derive the Klein–Gordon equation in five-dimensional spacetime, and obtain the evolution equation and the Schrödinger-like equation. With the resonances of the scalar KK modes as the initial data, the scalar field is evolved with the maximally dissipative boundary condition. The results show that there are scalar KK resonant particles with long life on the brane, which indicates that these resonances might be regarded as a candidate for dark matter.


Introduction
The nature of dark matter (DM) constitutes one of the most long-standing and puzzling questions in cosmology.There is abundant evidence that nonluminous matter makes up a large fraction of all matter in our universe.Results from cosmological measurements have now determined with exquisite precision the abundance of DM [1].But the identity of DM remains a mystery.Recently, Barranco et al. proposed that the scalar dynamical resonances can form long-lived configurations around black holes [2][3][4].This ultra-light scalar resonances might be regarded as one of the candidates for DM.Then, such long-lived distribution was generalized to the Dirac field by Zhou et al. [5].After that, these massive dynamical resonances around black holes attracted a lot of interest [6][7][8][9][10][11][12][13][14][15].The study of resonances of various fields around black holes has stimulated our interest in the evolution of resonances in theories of extra dimensions and braneworld.
a Qin Tan and Yu-Peng Zhang are co-first authors of the article.b e-mail:liuyx@lzu.edu.cn,corresponding author History of extra dimensions and braneworld dates back to the last century.In the 1920s, Kaluza and Klein (KK) proposed a five-dimensional spacetime theory to unify electromagnetic and gravitational interactions [16,17].Subsequently, extra-dimensional theories remained silent for more than 70 years.Until the 1990s, to solve the huge hierarchy between the Planck and weak scales, some braneworld models were proposed.Two of them have attracted the attention of many researchers.One is the large extra dimension model proposed by Arkani-Hamed, Dimopoulos, and Dvali [18], the other is the warped extra dimension model proposed by Randall and Sundrum [19].The size of extra dimensions is finite in these braneworld models.Subsequently, Antoniadis et al. embedded the braneworld model with large extra dimensions into the string theory [20].A great development was attained in Ref. [21], which shows that even the extra dimension is infinite, the four-dimensional gravity can also be recovered on the brane.After that, extra-dimensional theories have attracted a lot of attention .
In this paper, we focus on the evolution of the scalar KK resonances on the thick brane, which is generated dynamically by a background scalar field.In a braneworld model, to recover the physics in our four-dimensional spacetime, the zero modes of various fields should be localized on the brane.But in addition to zero modes, there are massive KK modes which might propagate into extra dimensions.KK resonances are a specific class of massive KK modes in braneworld models.Usually, for a volcano-like potential, although massive KK modes cannot be localized on the brane, but KK resonances could be quasi-localized on the brane [44].In previous literatures, resonances of various fields on thick branes have been studied [44][45][46][47][48][49][50][51][52].As far as we know, evolution of KK resonances of various fields in thick brane models has not been studied.The dynamics and the final state of such long-lived modes are still unclear.In the thin brane models, scattering of KK gravitons in the Randall-Sundrum-II model has been considered [53,54].It was proved that the brane possesses a set of discrete quasi-normal modes that appear as scattering resonance, and the KK modes of graviton have a very short lifetime on the brane.Can KK modes exist in a thick brane for a long time like the long-lived resonance modes around a black hole?To answer this question, we take the scalar field as a simple example to study the evolution of scalar KK modes numerically.We will show the evolution behaviors of the scalar KK modes and obtain their half-life time on a brane.Further, we will analyze the feasibility of the KK resonances as one of the dark matter candidates.Although our research is still crude, it provides a stepping stone to further investigation of KK resonances of various fields on a thick brane.
The layout of the remaining part of this paper is as follows.In Sec. 2, we construct a thick brane solution in five-dimensional spacetime as the background of a test scalar field evolution.In Sec. 3, the scalar field is evolved with the maximally dissipative boundary condition.Both resonances and nonresonances are used as initial data and their evolution behaviors are compared.Finally, the conclusions and discussions are presented in Sec. 4.

Braneworld model in general relativity
Firstly, we consider the thick brane in five-dimensional spacetime.For the simplest case of general relativity with a canonical scalar field, the action is given by where κ 5 is the five-dimensional gravitational constant.We set κ 5 = 1 in this paper for convenience.Besides, we consider the metric of the static flat brane as Here, e 2A(y) is the wrap factor and η µν = diag(−1, 1, 1, 1) is the four-dimensional Minkowski metric.In this paper, capital Latin letters M, N, • • • = 0, 1, 2, 3, 5 label the five-dimensional indices, while Greek letters µ, ν, • • • = 0, 1, 2, 3 label the four-dimensional ones.The dynamical field equations are By substituting the metric (2) into Eqs.( 3) and (4), we can obtain the explicit equations of motion where prime denotes the derivative with respect to the extra dimensional coordinate y.Note that only two of the above equations are independent, but we need to solve three functions: A(y), φ(y), and V (φ).So we need to give one of the three functions to solve these equations.The warp factor is where the parameter b has length dimension one and the parameter k has mass dimension one.From Eqs. ( 5), ( 6), (7), and (8) we get the solution where F(y, q) is the first kind elliptic integral and Π(y, q, p) is the third kind elliptic integral.This solution was investigated in Ref. [51] in f (T ) gravity theory.Based on this braneworld background, we consider a test scalar field and study its evolution.We investigate its dynamic behavior numerically, and clarify whether their KK modes can exist on the brane for a long time.
3 Scalar field resonances and its evolution in braneworld In this section we will consider the evolution of a scalar field in the thick brane given above.Here we consider a free massless test scalar field ψ(x M ).Notice that the scalar field ψ(x M ) here is not the background scalar field φ(y) that generates the thick brane.The equation of motion for the test field is the Klein-Gordon equation With the coordinate transformation dz = e −A dy, the metric (2) becomes and Eq. ( 11) can be written as Then, we introduce the following decomposition Substituting the above decomposition ( 14) into Eq.( 13), we get the following equation where a is a constant from the separation of variables.
The effective potential U (z) has the following form: or in the y coordinate equivalently The function Φ(t, z) can be further decomposed into oscillating modes as Substituting the above decomposition in Eq. ( 15) yields where m 2 = ω 2 − a 2 is the mass of the KK modes.We shall see below that, the effective potential U (z) is a volcano-like potential with a double-well.When z → ±∞, U (±∞) → 0 + .Thus, all KK modes with m 2 > 0 are free states.Only the modes with m 2 ≤ 0 could be bound states.One can show that there is no mode with m 2 < 0. The Schrödinger-like equation ( 19) can be factorized as the supersymmetric quantum mechanics form Using this supersymmetric quantum mechanics form with the conditions that A(z) is a real function and the extra dimension z is noncompact, it can be shown that there is no tachyon mode with m 2 < 0 [55,56].
On the other hand, the solution of the zero mode with m 2 = 0 is Obviously, for the warp factor (8), the zero mode u 0 (z) is bound on the brane and it has no node along the extra dimension.According to the node theorem, the zero mode must be the ground state, which ensures that m 2 ≥ 0. Solving Eq. ( 19) we could get a series resonant modes, which can be treated as the initial data for the scalar field.The evolution is dominated by Eq. ( 15).

Scalar field resonances
In this part, we give a brief review on how to solve the KK scalar resonances.Substituting the warp factor (8) into Eq.( 17), the effective potential in the coordinate y is given by [51] U (z(y) Plots of the above effective potential are shown in Fig. 1(a).
For convenience, we define the dimensionless parameters b = kb and m = m/k.We can see that the width of the effective potential increases with b.The resonant modes can be studied by the relative probability method which was proposed in Ref. [44].The relative probability is defined as where u(z) is solved from Eq. ( 19), z b is approximately the width of the brane, and z max is a much larger width than z b , and usually set to 10z b .If the relative probability has a peak with full width at half maximum around m = m n , then there is a resonance with mass m n .In this way, the modes whose amplitudes in the quasi-well are much larger than those outside the quasi-well can be found.These modes will remain on the brane for a longer time.Note that, the wave functions can be even or odd because the potential is symmetric.Hence, the following boundary conditions can be used to solve Eq. ( 19) numerically: where u even denotes the even modes of u(z) and u odd denotes odd modes of u(z).Substituting the effective potential ( 22) into the Schrödinger-like equation ( 19), we can solve the solution of u(z) numerically for a given mass m.Then the relative probability P (m 2 ) can be obtained.The relative probability P (m 2 ) of scalar resonances for b = 5, 10, 15 are shown in Figs.1(b), 1(c), 1(d), respectively.The specific parameters of these resonances are listed in Tab. 1.It can be seen that the number of scalar resonances and their peak value increase with b.Usually, the larger peak value means smaller full width at half maximum, and longer lifetime.We will see that in the next subsection.Fig. 1 The shapes of the effective potential (22) for different parameters b.And the influence of the parameter b on the relative probability P for the odd-parity (red lines) and even-parity (blue dashed lines) massive KK modes.These figures are from Ref. [

Evolution of the scalar field resonances
Treating the scalar resonances as the initial data, we can evolve the scalar field under the evolution equation (15).Through the numerical evolution of the scalar field, we can obtain its lifetime on the brane.In this paper, we only consider the case of a 2 = 0, which means that scalar KK particles travel along the extra dimension at the speed of light at infinity.And we impose the maximally dissipative boundary condition: where n is the outward unit normal vector to the boundary.Equation ( 15) is solved numerically using fourth-order finite differences in space, and evolving in time using a method of lines with a thirdorder Runge-Kutta integrator.
In order to more intuitively display the evolution of the scalar field, we define the conserved energy of the scalar field [58] where Firstly, we consider the evolution of resonances whose parameters are given in Tab. 1.We integrate the energy density ρ E over [−z max , z max ], the resulting energy will decay due to energy loses through both left and right boundaries.We plot the evolution of the integrated scalar field energy E(t) in Fig. 2. Note that it is plotted on a logarithmic scale.From Figs. 1(b), 1(c), 1(d) we can see that there are several resonances, usually.Among these resonances, the first one will evolve the longest time.Besides, the evolution time increases with the parameter b.This is because that the larger b the more energy of the scalar field is concentrated in the potential well.In addition, we plot half-life time of the first resonance with different values of the parameter b, which can be seen from Fig 2(d).Note that here we have defined the dimensionless time t = kt.It can be clearly seen that the half-life time of the resonances increases with the parameter b.
The energy decay can be fitted as an exponential function: where s is the fitting parameter and E 0 denotes the initial energy of the KK mode.Some results of the fit are listed in Table 2.It can be seen that both the scaled mass m1 and the fitting parameter s of the first resonance decreases with b.
To better show the evolution of the scalar field over time, we analyze the result of the numerical evolution by extracting a time series for the resonance amplitude at a fixed point z ext .The results are shown in Fig. 3.It can be seen that the resonance amplitude decreases with evolutionary time, and the amplitude attenuation of the second resonance is significantly faster than that of the first resonance.This also shows that the first resonance will evolve longest.By comparing the amplitude attenuation at different positions of the same resonance, it can be seen that the attenuation rate and the overall shape are basically the same.In other words, at least for the first two resonances, there is no beating effect found in Ref. [59] for a black hole system.
As a comparison, we also consider the evolution of the nonresonances.These results are shown in Fig. 4. We find that the energy and amplitude of nonresonance decay very fast at early stage, but later they decay like those of resonances.In order to gain a better understanding of the above results, we perform a spectral analysis.We calculate the discrete Fourier transform in time of the scalar field at a fixed point z = z j .The explicit expression of the discrete Fourier transform is where A is normalization constant and t p are the discrete time values.Plots of the Fourier transform for the first resonance and the nonresonance with m2 = 0.36 for b = 15 are shown in Fig. 5.We find that, for the Fourier transform of the resonance, there is only one peak corresponding to the resonance frequency.However, for the Fourier transform of the nonresonance, there are several peaks.Thus, nonresonances can evolve into combinations of resonances, and from this point of view, resonances seem to play a similar role in the braneworld as the quasi-normal modes in black holes physics, which deserves further investigation.Finally, we consider the half-life time of the scalar resonances.If the exponential decay ( 27) is sustained throughout all evolution, then it is easy to know that t1/2 = ln (2)  s .For k = 10 −10 eV and s = 10 −10 , the halflife time t 1/2 of the first scalar resonance will reach 10 4 seconds.Admittedly, it is still short compared with the age of our universe.But note that the lifetime of resonances increases with the parameter b.Thus, for a very large b, the lifetime of resonance might reach the cosmological time scale.Here, we can simply estimate the feasibility of resonance as a candidate for dark matter.For the braneworld model in this part, the effective four-dimensional Planck scale M Pl and the fundamental five-dimensional scale M 5 have the relation:    According to the current experiment of the Large Hadron Collider, the collision energy is 13 TeV but no signal of extra dimensions is seen [60].So the fundamental fivedimensional scale M 5 should be greater than 13 TeV.Thus, combining Eq. ( 29) and the condition M 5 13 TeV, the constraint on the parameter k is given by From Tab. 2, we can know that the fitting parameter s rapidly decreases with b.In fact, when b doubles, the fitting parameter s decreases by an order of magnitude.
If we choose b = 10

Conclusion and discussion
In this paper, we investigated the evolution of a free massless scalar field in the thick brane numerically.We find that the resonances decay very slowly compared to the nonresonances and can exist on the brane for a very long time.If the lifetime of these resonances can be long enough as the cosmological time scale, they might be a candidate for dark matter.This provides a new idea for dark matter research.Firstly, we constructed a five-dimensional thick brane generated by a scalar field.Then, we considered the evolution of a test scalar field in this thick brane background.Through the coordinate transformation and the variable separation, we got the evolution equation ( 15) and the Schrödinger-like equation (19) for the extra dimension profile of the scalar field.The latter gives us the initial data of the massive KK modes, and the former evolves those initial data.Next we solved the Schrödinger-like equation (19) numerically to obtain the initial data of massive KK modes, especially the initial data of the resonances.The results were shown in Fig. 1 and Tab. 1.We found that the relative probability of the first scalar resonance increases with b, while the mass m 1 of the first scalar resonance decreases with b.Using these KK modes as initial data, we investigated their evolution.Considering their energy decay and extracting a time series for the resonance amplitude, the evolution of the scalar field was analyzed.The results were shown in Fig. 2, Fig. 3, and Tab. 2. The energy decay can be described by a decay parameter s which can be obtained by an exponential fit of E/E 0 as the function of t.On the other hand, we also considered the evolution of the nonresonances.The energies and amplitudes of nonresonances decay rapidly at early stage, but later they decay like those of the resonances.The behavior of the nonresonance evolution could be treated as a combination of resonances.Finally, we considered the half-life time of the resonance.For a very large b, the lifetime of ultra-light resonances can reach the cosmological time scale.These indicate that the scalar resonant mode could be one of the candidates for dark matter.
There is much to be improved in this paper.Firstly, the mass and interaction of the test scalar field is not taken into account, which may result in more interesting results.Secondly, other test fields such as Dirac spinor fields and gauge fields, their evolution is also worth investigating.These possibilities deserve further study.
The relative probability with b = 15

2 (d) b − t1/ 2 Fig. 2
Fig. 2 Figures.2(a), 2(b), 2(c) illustrate the energy of the scalar field vs. time for the evolution of the resonant modes for different values of the parameter b.Here o1 represents the first odd-parity resonance, e1 represents the first even-parity resonance, and so on.Figure.2(d) depicts the relation between the half-life of the first scalar resonance and the parameter b.

Fig. 3 30 Fig. 4
Fig. 3 Upper panel: Time evolution of the first odd resonance (left) and the first even resonance (right) at zext = 3 for b = 15.Lower panel: Time evolution of the first odd resonance (left) and the first even resonance (right) at zext = 30 for b = 15.Here zext = kz ext .

3 Fig. 5
Fig. 5 Left panel: Discrete Fourier transform in time vs. frequency for the evolution of the nonresonance with m2 = 0.36 for b = 15.Right panel: Discrete Fourier transform in time vs. frequency for the evolution of the first resonance with m2 1 = 0.04687 for b = 15.The blue dotted lines correspond to the frequencies of the first three odd parity resonances.

Table 1
Resonant mass spectrum m2 n , mn , and relative probability P for different values of the parameter b.

Table 2
The first resonant mass spectrum m2 1 , m1 , fitting parameter s, and half-life t 1/2 for different values of the parameter b.
8, thus, k ≥ 10 −8 eV, and s might be 10 −25 , and the half-life time t 1/2 of the first scalar resonance will reach 10 17 seconds.It is the same order of magnitude as the age of our universe (4.35 × 10 17 seconds).On the other hand, the large b means that the first resonance is very light, but they are still hard to spot in the collider.The reason is that the cross section of any process involving the interaction of the zero mode with the light continuum modes is imperceptibly low.From this perspective, these long-lived resonances with very light mass could be considered as a candidate for dark matter.