Melting of scalar mesons and black-hole quasinormal modes in a holographic QCD model

A holographic model for QCD is employed to investigate the effects of the gluon condensate on the spectrum and melting of scalar mesons. We find the evolution of the free energy density with the temperature, and the result shows that the temperature of the confinement/deconfinement transition is sensitive to the gluon-condensate parameter. The spectral functions (SPFs) are also obtained and show a series of peaks in the low-temperature regime, indicating the presence of quasiparticle states associated to the mesons, while the number of peaks decreases with the increment of the temperature, characterizing the quasiparticle melting. In the dual gravitational description, the scalar mesons are identified with the black-hole quasinormal modes (QNMs). We obtain the spectrum of QNMs and the dispersion relations corresponding to the scalar-field perturbations of the gravitational background, and find their dependence with the gluon-condensate parameter.


Introduction
It is well known in quantum chromodynamics (QCD) that the gluon condensate has a relevant role in the low-energy dynamics; in the absence of quarks it represents the vacuum of QCD [1]. It is also known that the low-energy regime of QCD cannot be studied with the usual mathematical tools employed in the regime of high energies, where the coupling constant is small and one can use perturbative techniques to obtain information about the system. Differently from the high-energy regime, the coupling constant is large for low energies and perturbative techniques are no longer applicable. There are some approaches attempting to obtain information about the low-energy dynamics. One of these approaches is the lattice field theory, where the problems are discretized and solved on a lattice, using powerful computational resources (see e.g. Ref. [2] for a review on this subject). Another technique is the operator product expansion (OPE), also called the SVZ (Shifman, Vainshtein, and Zakharov) sum rules, which tries to extend some perturbative results to the low-energy regime [3]. Both approaches have limitations and cannot provide a full real-time description of the low-energy regime of QCD.
In recent years, as an alternative tool, a variety of gravitational holographic models has been used to study the non-perturbative regime of QCD. In this context, there are two main roads to follow. The so-called top-down approach has a ten-or eleven-dimensional superstring solution as starting point. After some compactifications, it is obtained a five-dimensional effective gravitational model, which is dual to a four-dimensional conformal field theory (CFT) living at the boundary of a curved spacetime with a negative cosmological constant. This curved spacetime is known as anti-de Sitter or AdS for short (for a review on top-down models, see for instance Ref. [4]). In the second approach, the so-called bottom-up models are designed so that the dual gravitational theory reproduces some known results or properties in QCD. In these models, the gravitational dual theory does not need to follow as a classical limit of a superstring theory (for a review on bottom-up models, see Ref. [5]).
The main advantage of the gauge/gravity duality, formulated in its original form in [6][7][8], is that we can map problems in a strongly coupled field theory, living in a flat d-dimensional spacetime, to problems in a (d + 1)-dimensional classical theory of gravity. This map is implemented by associating each local operator in the quantum field theory to a classical field in the gravitational side, i.e., O(x 0 , · · · , x d ) ↔ X(x 0 , · · · , x d , x d+1 ). For example, in four-dimensional QCD the operator Tr F 2 , that characterizes the scalar sector of the gluon field F µν , is dual to a scalar field Φ, known as the dilaton. This correspondence is such that the value of the dilaton at the boundary, Φ(x 0 , · · · , x 4 , x 5 = z = 0) = φ 0 (x 0 , · · · , x 4 ), is a source for the operator Tr F 2 , which means they are coupled as dx 4 φ 0 (x)Tr F 2 (x).
The first attempt to take into account the gluon condensate in a holographic bottom-up approach for QCD was developed by Csaki and Reece [9]. By considering a dilaton that couples to the operator Tr F 2 at the AdS boundary, they showed that the asymptotic behaviour of this field close to the boundary must be of the form Φ(z) = φ 0 + G z 4 , where z is the holographic coordinate. This result is consistent with the general asymptotic solution for a massless scalar field close to the boundary: Φ(z) = φ 0 z 4−∆ + G z ∆ , where φ 0 , G and ∆ = 4 are interpreted, respectively, as the source, the vacuum expectation value (VEV) and the conformal dimension of the operator Tr F 2 (the VEV is also known as the gluon condensate). In spite of a relative success with the introduction of the gluon condensate, the authors of Ref. [9] obtained a scalar glueball spectrum which does not follow a linear Regge trajectory.
Nevertheless, there are some well-succeeded approaches to QCD, such as the soft-wall model [10] and the improved holographic model [11,12], which consider a quadratic dilaton in the IR and thus guarantee a linear behaviour for the glueball and meson spectra. Motivated by these works, in this paper we implement a model with a dilaton field which is quartic in the UV (to describe correctly the gluon condensate) and quadratic in the IR (to guarantee linear behavior of the spectrum). As in the soft-wall model, the background metric is fixed, and the conformal symmetry is explicitly broken with the introduction of an exponential dilaton-dependent term in the five-dimensional action. Such a term does not modify the gravitational background and generates a dual field-theory energy-momentum tensor with nonvanishing trace.
Among the important quantities that are obtained using our model are the spectral functions (SPFs), which are fundamental for the understanding of the hadronic properties and the vacuum structure of QCD. Additionally, the SPFs at finite temperature may shed some light about the influence of the surrounding medium on the hadronic internal structure. In Refs. [13,14] the maximum entropy method was used to construct the SPFs in lattice QCD. In holography, a dual description of a finite temperature field theory is obtained by considering a black hole in the gravitational AdS background [15]. The authors of Ref. [16] developed a prescription to calculate retarded Green functions from the five-dimensional on-shell action, allowing to get information of the finite-temperature field theory from quantities defined on the gravitational side. As an example we mention the hydrodynamic transport properties obtained in Refs. [17,18], where the poles of the low-momentum limit of retarded Green functions are identified with the low-lying quasinormal modes (QNMs). In fact, it was shown that the poles of the finite-temperature correlation functions are related to the full spectrum of QNMs, not just in the hydrodynamic regime. From the gravitational point of view, the spectrum of QNMs are produced by perturbations on the AdS black-hole spacetime that satisfy specific boundary conditions: an ingoing wave condition at the horizon and a Dirichlet condition at the boundary of the AdS spacetime. Interestingly, in recent years it has been an increasing interest on the nonhydrodynamic QNMs, because these modes dominate the dynamics of fluctuations before the system reaches the local hydrodynamic equilibrium (see Refs. [19][20][21] for a discussion). In addition, the study of the nonhydrodynamic QNMs may shed some light on the understanding of early time dynamics of some strongly coupled field theory models (for example, the QCD quark-gluon plasma). For discussion and details on quasinormal modes, see Refs. [22,23].
The present paper is organized as follows. In Sec. 2 we introduce the action and the equations of motion of the model and define the dilaton field that interpolates between the UV and IR. Then, we obtain the mass spectrum of mesons at zero temperature and explore the effects of the parameter associated with the gluon condensate on this spectrum. Sec. 3 is devoted to the study of the finite-temperature effects and, to do that, we calculate the free energy density of the thermal and black-hole states. Yet in Sec. 3 we obtain the SPFs and make explicit the dependence of the results on the gluon condensate and temperature. In Sec. 4 we present the spectrum of QNMs calculated by using three different numerical techniques: power series, Breit-Wigner, and pseudo-spectral methods. A discussion and comparison between the methods are also presented. To complement the numerical analysis, in Sec. 5 we present and discuss the results of the dispersion relations. In Sec. 6 we conclude with some remarks and the final comments.

Scalar mesons from holographic QCD
In this section we explore a holographic description of the scalar sector of the mesons. This is done by employing the action proposed in Ref. [10], whose respective approach is known as the soft-wall model. It is a five-dimensional effective action with two scalar fields, the dilaton Φ(z) and a scalar S(x µ , z), both introduced as probe fields, which means that the backreaction of the fields on the geometry is neglected. The holographic dictionary [6] establishes that the dilaton is dual to the operator Tr F 2 [9], while the scalar field S is dual to the operatorqq [24,25].

Equations of motion for the fields
We start with the metric that describes the five-dimensional anti-de Sitter spacetime. Such a metric is a solution of the Eintein equation with a negative cosmological constant, and takes the following form in Poincaré coordinates: where the five-dimensional spacetime metric has the signature (−, +, +, +, +). The warp factor is defined by A(z) = log ( /z), where is the AdS radius. From now on, in this work, we set the radius = 1 to simplify the notation. In the system of coordinates {x µ , z}, the boundary field theory lies at z = 0, which is identified as the UV fixed point, while z → ∞ is the deep IR region. This identification is possible because the warp factor and the energy scale E of the dual field theory are related by e A(z) = E [11,12].
The five-dimensional action, that describes chiral symmetry breaking in the meson sector and contains SU (2) L × SU (2) R gauge fields with a bifundamental scalar field X, can be written as [10] with A L and A R the gauge fields and X the scalar field (tachyon) responsible for the chiral symmetry breaking SU (2) L × SU (2) R → SU (2) V (for details, see Ref. [26]). Following Refs. [24,25] the scalar sector of the foregoing action can be obtained by turning off all the background fields except the fluctuation of the bifundamental field, which is decomposed in the form X = X 0 (z) + S(x, z), with S(x, z) being a perturbation on the background value X 0 (z). Hence, the action that describes the scalar sector of the mesons is given by where m 2 X = −3 is the mass of the scalar field. The dimension of the operator dual to the scalar field S(x, z) satisfies the equation ∆(∆ − 4) = m 2 X [7]. Solving this algebraic equation, we get ∆ = 3, which means the operator dual to S has dimension 3.
The equation of motion for the scalar field is obtained by varying the action (2.3) with respect to S, Since the metric and the dilaton field depend on the holographic coordinate only, it is convenient to Fourier decompose the scalar field as S(x, z) = S(z)e −ikµx µ , and then Eq. (2.4) reduces to where m 2 s = −k µ k µ is the square of the scalar-meson mass in the four-dimensional field theory. Equation (2.5) can be recast into a Schrödinger-like form by introducing the auxiliary function 2B(z) = 3A(z) − Φ(z), and using the Bogoliubov transformation S(z) = e −B(z) ψ(z). In doing so, where indicates d/dz. In terms of the dilaton field, the effective potential takes the form We now turn attention to the dilaton field Φ(z) and consider a functional dependency on z that is consistent with some aspects of QCD, namely, Φ(z) ∼ z 4 in the UV (to describe correctly the gluon condensate) [9], and Φ(z) ∼ z 2 in the IR (to guarantee confinement) [12]. To smoothly connect these regimes, we use an interpolation function in the form [27] such that the asymptotic behavior of Φ(z) close to z = 0 becomes where φ 0 is the source that couples to the dual operator Tr F 2 , G is the energy scale associated with the gluon condensate, and the ellipses indicate subleading terms. In the deep IR region, the dilaton field goes like Φ(z) = C z 2 + · · · z → ∞, (2.10) where C characterizes the confinement energy scale and the ellipses mean subleading contributions.
The parameter C can be chosen by matching the smaller eigenvalue of the vectorial sector of the action (2.2) with the experimental mass of the lightest ρ meson. This was done in Ref. [28], and the resulting value is C = 0.151(GeV ) 2 . In the forthcoming sections we are not going to fix the parameters, since the aim of this paper is to show qualitative instead of quantitative results associated to this holographic model. For numerical purposes, we rewrite Eq. (2.8) in terms of the dimensionless coordinate u = z √ C, where the information of the gluon condensate is now contained in the dimensionless parameter G = G/C 2 . After the change of coordinate z → u, the effective potential and the mass m s in the Schrödinger-like equation (2.5) are normalized by the parameter C asV = V /C andm s = m s / √ C. In the next sections, we are going to investigate if (and how) the spectrum of scalar mesons depends on the parameter G.

Analysis of the zero-temperature effective potential
Here we point out the differences between our approach and the original soft-wall model [10]. As commented above, we are considering a dilaton field that is quartic in the UV and quadratic in the IR. This difference in relation to the soft-wall model, where the dilaton is quadratic from the UV to IR, requires the introduction of a new free parameter related to the energy scale that characterizes the gluon condensate [9]. The dimensionless version of this new parameter is identified with G in Eq. (2.11).
The asymptotic behaviour of the potential in the UV depends on G. This statement can be appreciated expanding the potential (2.7) near the boundary, where the ellipses indicate subleading terms, expressed as higher order powers of u.
On the other hand, the asymptotic form of the potential (2.7) in the deep IR region can be written as V (u) = u 2 + · · · u → ∞, (2.13) where the ellipses represent subleading terms, suppressed as powers of 1/u. By comparing the above results with those of the original soft-wall model [10], we notice the influence of the quartic dilaton on the asymptotic UV behavior of the effective potential, which depends on the value of the gluon-condensate dimensionless parameter G. To see this difference quantitatively, we plot in Fig. 1 the potential obtained in the original soft-wall model [10] and the numerical results obtained from Eq. (2.7).
As expected, the asymptotic behaviors are similar and the main differences between our approach and the original soft-wall model lies in an intermediate region between the UV and IR, as it can be seen in Fig. 1.

Asymptotic solutions for the scalar field
Once the asymptotic behaviour of the effective potential is known, we can now find the asymptotic solutions of the Schrödinger-like equation (2.6). The two regimes, UV and IR, are dealt with separately. To obtain the asymptotic solutions in the UV regime, we first replace the leading terms of the potential from Eq. (2.12) into Eq. (2.6), which leads to Then, substituting ψ(u) = u α into the last equation and making some simplifications, we obtain the following algebraic equation: By taking the limit u → 0 and solving for α, we find α = −1/2 and α = 3/2, so that the general asymptotic solution for ψ(u) in the UV is given by The first term on the right hand side of (2.16) is the non-normalizable solution (see Appendix of Ref. [27] for details), while the second term is the normalizable one. As we are uniquely interested in normalizable solutions, we set c 1 = 0 and consider a wave function of the form Now one might be surprised about the fact that the parameter G does not affect the asymptotic solutions (2.16). This is because the asymptotic solution of the scalar field S(z) in Eq. (2.5) is not affected by a quartic dilaton in the UV. To see that we substitute the warp factor A(u) = log √ C/u and the dilaton field (2.9) in Eq. (2.5) and then solve the indicial problem for S(u) = u β . The resulting asymptotic solution is where a 1 is interpreted as the quark mass and a 2 as the chiral condensate [25]. Since the functions S(u) and ψ(u) are connected by a Bogoliubov transformation ψ = e B S, it follows that which reduces to the solution (2.16) as u → 0, since the contribution from the exponential factor alters only the subleading terms.
In the deep IR regime the asymptotic behaviour of the potential is given by Eq. (2.13) and the Schrödinger-like equation (2.6) reduces to Hence, the solution that guarantees convergence of the wave function in this region is given by [27] This function, which gives the asymptotic behaviour of the wave function ψ(u) in the IR region, and the solution (2.17), that represents the asymptotic behaviour of ψ(u) in the UV region, are used below in the search for a full solution of Eq. (2.6) through numerical methods.

Analysis of the mass spectrum
As pointed out above, the normalizable solutions of the differential equation (2.6) are associated to scalar-meson states. We obtain these solutions by solving numerically the Schrödinger-like equation (2.6) with a shooting method. For the numerical integration, boundary conditions need to be provided. In the present case, we use the near-boundary UV (u = 0) normalizable solution (2.17) and its derivative as "initial" conditions, and complement them by imposing regularity of the wave function in the deep IR region. These conditions are satisfied just for a discrete set of values of the mass parameterm 2 s . Notice that, alternatively, we might use the asymptotic solution (2.21) in the IR (u → ∞) and its derivative as "initial" conditions, and require a regular behaviour of the wave function in the UV boundary u = 0. Both approaches give the same solutions to the eigenvalue problem.
In Table 1 we present the first nine eigenvaluesm 2 s for two different values of the parameter G; namely, G = 1 and G = 50. These particular values are chosen arbitrarily and their separation is taken large enough to clearly display the dependence of the mass spectrum on such a parameter. These results show that the lower excited states are more sensitive to the change of the parameter G than the higher excited states. To make this fact transparent, we plot the data of Table 1 in Fig. 2, where the differences among the two spectra at low masses is clearly seen. For comparison, let us mention the previous results found in the original soft-wall model [10,29,30], whose mass spectrum has the closed formm 2 s = (4 n + 6), n = 0, 1, 2, · · · . As it can be seen, the fit for G = 50 approaches the result presented in Eq. (2.22), while the fit for G = 1 is completely different.
Although it is not considered the backreaction of the dilaton on the metric, the condensate is shown to be important for the dynamics of hadron formation in the dual field theory. It is worth mentioning that the spectrum tends to a continuum as the parameter G goes to zero. In such a limit, it is recovered the problem of a massive scalar field in an AdS background with a constant dilaton.
In the sequence of this work we explore the finite-temperature effects on the mass spectrum of scalar mesons, and in the presence of the gluon condensate (parameter G).

The finite-temperature holographic QCD model
To take into account the finite-temperature effects on the dual field theory, we need to consider a black hole on the gravity side. The standard static black-hole solutions include a horizon function, f (z), in a metric of the form (2.1). The background geometry we consider here is described by where f (z) = 1 − z 4 /z 4 h , δ ij is the Kronecker delta and the Latin indices run over the spatial transverse coordinates (i = 1, 2, 3). The domain of the holographic coordinate is 0 ≤ z ≤ z h , where z h indicates the position of the black-hole event horizon. The Hawking temperature of the black hole is given by According to the AdS/CFT dictionary, T is also the temperature of the dual thermal field theory. It is possible to observe the parameter dependence of the confinement/deconfinement temperature (or critical temperature) for gluons by analyzing the free energy of the dual thermal field theory. What is more, the free energy is obtained by calculating the Euclidean on-shell action of the background fields. In the present case, the gravitational action for both cases, the thermal AdS spacetime (2.1) and the black-hole spacetime (3.1), are given by [28] where κ is the gravitational coupling, andg mn and g mn are the corresponding metrics. To get the Euclidean on-shell action, which is related to the free energy of the dual thermal field theory, we first obtain the Einstein equations (in Ricci form) from (3.3), Then, substituting these results into Eq. (3.3) we get (for the thermal phase) where V 3 is the three-dimensional transverse volume,β corresponds to the period of the Euclidean time t ∼ t + iβ, and z 0 is the UV cutoff which lies close to the boundary. Moreover, doing the same procedure we obtain the on-shell action of the black-hole phase, where the Euclidean time belongs to the interval 0 ≤ t ≤ β, and z 0 is the UV cutoff located close to the boundary. Furthermore, the temperature T = 1/β is related to the event-horizon coordinate z h through β = πz h . The relation between the free energy F and the Euclidean on-shell action is βF = S on-shell . However, in both cases the on-shell action blows up at the UV cutoff [28]. To avoid these divergences we are going to use a prescription for which matters the difference of the black hole and thermal phases [28]. Thus, the variation of the free energy density F = F/V 3 is given by Differently from the original soft-wall model [10], here it is not possible to get an analytical expression for ∆F, so that Eq. (3.7) must be solved numerically. To guarantee the periodicity in time at the UV cutoff, we setβ = β f (z 0 ) [28]. In Fig. 3 we display the numerical results obtained for G = 1 (blue line) and G = 50 (red line). The results in this figure show the parameter dependence of the difference in the free-energy density. Furthermore, the confinement/deconfinement transition temperature is defined as the temperature for which the free-energy difference is zero. This means that we must solve the equation ∆F(z hc ) = 0, where the solution z hc is related to the critical temperature as z hc = 1/πT c . Moreover, from Fig. 3 we observe that z hc (G = 50) < z hc (G = 1), which means that T c (G = 50) > T c (G = 1) and, consequently, the critical temperature is a function of the parameter G. This result supports the idea of a temperature dependence of the gluon condensate. For a discussion in QCD see, for instance, Ref. [34] and in holography Ref. [35].
To complement this section we analyze the phase transition from the point of view of the degreees of freedom of the systems. A good thermodynamic variable to characterize this issue is the entropy. From the holographic dictionary we may determine the entropy of the dual field theory by calculating the entropy of the gravitational background by using the Bekenstein-Hawking formula, i.e., S = A/(4G 5 ), where A is the event horizon area. As discussed in Ref. [36], the degrees of freedom of the deconfined phase is proportional to N 2 (N is the number of colours), this means that S ∼ N 2 , while in the confined phase is proportional to N 0 , S ∼ N 0 . Therefore, the entropy function S(N ) has a phase transition because it is discontinuous at T c . Finally, as the dilaton field is dual to the operator TrF 2 , and the metric dual to the energy momentum tensor of the Yang-Mills SU (N ) field, consequently, in the deconfined phase the gluons are free. We have set C = 1 and κ 2 2 = 1.

Equations of motion for the scalar field
The equation of motion for the scalar field S(x, z) in the black-hole geometry is obtained by means of Eq. (2.4) with the background metric (3.1). It may be written explicitly in the form As in section 2, introducing the function 2B = 3A − Φ and the tortoise coordinate ∂ r * = −f ∂ z in Eq. (3.8), and taking the Fourier transform where q 2 = k · k is the squared modulus of the spatial part of the four-momentum vector k µ . Finally, to get a Schrödinger-like equation, we introduce the transformation S = e −B ψ and, after some simplifications, we find where the potential is given by In the present case, it is helpful to write the explicit form of the tortoise coordinate in terms of the holographic coordinate z. Such a relation is given by This coordinate ranges from r * → −∞ (for z → z h ) to r * = 0 (for z = 0). An useful form for the effective potential as a function of z is given by from which one immediately sees that the potential vanishes at the event horizon.

Asymptotic solutions for the scalar-meson field
Here we study the asymptotic solutions of the differential equation (3.11). Since the effective potential (3.14) vanishes at the horizon, because f (z h ) = 0, the solutions for ψ at z = z h are of the plane-wave form, where the first term is a purely ingoing wave, while the second term corresponds to a purely outgoing wave. However, in the neighbourhood of the event horizon the wave functions can also be expressed as power series h /4 . After substituting these series into the differential equation (3.11), we get, up to second order, where the ellipses represent higher-order contributions, and the coefficients are given by It is worth to mentioning that the leading terms in Eqs. (3.15) and (3.16) are identical. This can be shown by substituting the tortoise coordinate as a function of z into Eq. (3.15) and expanding the resulting function in power series of (1 − z/z h ). We prefer to work with Eq. (3.16) instead of Eq. (3.15) in the forthcoming sections. Now we turn attention to the solutions near the boundary z = 0. The power series expansion of the two independent solutions may be written as where the ellipses stand for higher order terms, and the coefficients are given by Notice that ψ (1) is a normalizable wave function while ψ (2) is not, since it diverges at the boundary.
In order to calculate the spectral function associated to the scalar mesons (in Sec. 3.6), we follow the same procedure applied in Refs. [37,38]. First, we write the near-horizon ingoing and outgoing solutions, ψ (−) and ψ (+) , as a linear combination of the wave functions close to the boundary, where the coefficients are functions of the wave number and frequency. Second, we may also write the near-boundary solutions as a linear combination of the near-horizon solutions, Doing so, we realized that there is a relation between Eqs. (3.22) and (3.23) which can be written in a matrix form (for details about this relation, see Refs. [37,38]). To get the spectral function associated to the mesons, we need to know some of the foregoing coefficients (A (−) and B (−) , say) and so, in the sequence of this work, we find numerically the values of these coefficients.

Analysis of the finite-temperature effective potential
In this subsection we develop a careful analysis of the effective potential (3.14) and study how it varies with the parameters of the model. First of all, let us write the potential explicitly as a function of the dilaton field Φ, the holographic coordinate z and the black-hole temperature T ,

(3.24)
Now it is helpful to consider some particular cases, beginning by the zero-temperature potential.
From Eq. (3.24), it is clear that we recover the potential (2.7) with the additional term q 2 for T = 0, Thus, by choosing the dilaton as being a quadratic function of the holographic coordinate z, we recover previous results from the literature (see, e.g, Refs. [30,39]). Now let us turn off the gluon condensate by taking G = 0. In doing so, the dilaton field (2.8) reduces to a constant, i.e., Φ(z) = φ 0 , and the potential (3.24) becomes This is the same potential as that of a pure AdS black hole metric with a massive scalar field.
To implement the numerical analysis of the effective potential, it is convenient to normalize all parameters of the model in order to make them dimensionless. The temperature is normalized by the confinement energy scale as T = πT / √ C [37,38], and the wavenumber is normalized by the temperature,q = q/πT . We introduce again the dimensionless gluon condensate G = G/C 2 . The tortoise coordinate is normalized as √ C r * and the effective potential by V /C. For a better visualization of the graphs of the potential, we separate the analysis in two regimes: one for low temperatures ( T < 1) and another for intermediate and high temperatures ( T ≥ 1). In Fig. 4 we show the results for intermediate and high temperatures, and two values of the dimensionsless condensate, G = 1 (left panel) and G = 50 (right panel). These values were chosen arbitrarily and are the same that we used to calculate the spectrum of the scalar mesons at zero temperature presented in Table 1. By comparing both figures, we observe that the effect of the dimensionless gluon condensate is more relevant for intermediate temperatures; see, for example, the results for T = 1. The results for high temperatures are less sensitive to this parameter, as can be seen by comparing the results for T = 15.
Another interesting result is the temperature for which the potential well starts to form, as is shown in the small box of Fig. 4. This temperature, that we call as T w , depends on the value of G, such that T w = 0.782 for G = 1 is greater than T w = 0.497 for G = 50. As is shown in both figures, for T > T w the presence of bound states is not expected due to the absence of a potential well. The potential well starts to form for temperatures smaller than T w . To finish this case, it is noticed the temperature effect in deforming the potential (2.7), as can be easily seen in the expression (3.24). For high temperatures, the term T 8 becomes leading and its effect is visualized in Fig. 4.  The curves for the potential at low temperatures are shown in Fig. 5. In this regime the potential presents a well and a barrier, with the height of the barrier depending on the temperature. Here we present the results for selected values of temperature such that the maximum value of the potential is equal to the mass of the zero-temperature spectrum displayed in Table 1 and plotted as horizontal lines in this figure. We did this analysis for G = 1 (left panel) and G = 50 (right panel). Our intention is always to compare with the results at zero temperature obtained in Sec. 2. From Fig. 5 we see that the potential well becomes deeper with the decreasing of the temperature, and this behaviour is the same obtained for scalar glueballs and vector mesons in Refs. [37,38]. For example, it is possible to have five bound states if the temperature is smaller than T = 0.121 (on the left panel of this figure). But since the width of the potential is finite, these bound states have a finite lifetime, and the way they decay depends on the energy they have. For a fixed value of temperature, the higher exited states decay faster that the lower exited states. A careful observation of Fig. 5 shows how the potential (2.7) is deformed in the regime of low temperatures.
As it is known from special relativity k µ k µ = −m 2 s , in our case we are neglecting the spatial part of the wave vector, this result reduces to ω 2 = m 2 s . This means that at zero temperature the frequency is equal to the mass, and the potential reduces to Fig. 1, this is an exact result. But when we add temperature the Lorentz symmetry is broken and this relation is no longer valid k µ k µ = −m 2 s (when temperature is introduced there is a preferential observer, this breaks down the Lorentz symmetry due that this is a global symmetry and the temperature is local). A consequence of this is that the frequency gain an imaginary part, that is related to the lifetime of the bound states. This result is supported by the form of the potential in Fig. 5.  From these results is more evident the conclusion that the dimensionless gluon condensate has a strong influence on the low exited states and the higher exited states are less sensitive to this parameter as can be seen by comparing the temperatures in both panels in Fig. 5.
To finish the analysis of the effective potential we write the potential close the boundary in its asymptotic form where the ellipses represent higher order contributions in the temperature and holographic coordi-nate. Or as a functions of the tortoise coordinate this result becomes V = 3 4r 2 * + 6 G 2 C r 4 * + · · · + 9 5r 2 * − 24 G 2 5C r 4 * + · · · r 4 * π 4 T 4 + · · · , (3.28) as before, the ellipses represent subleading contributions. In this result is more evident the contribution of the temperature to deform the potential. Differently from the zero temperature case, if we set G = 0 we still have bound states at finite temperature, the last term in Eq. (3.28) guarantee the existence of a potential well.

Retarded Green function
To find the real-time response in the dual field theory, we need the retarded Green function. The AdS/CFT dictionary allows us to find the correlation functions in the boundary field theory in its Euclidean version [7,8]. A prescription to find two-point correlation functions was proposed in Ref. [16], and such a prescription is equivalent to the on-shell action re-normalization strategy as developed in [40] (see also the references therein). It is important to point out that the scalar field S(x, z) satisfies specific boundary conditions: Dirichlet at the AdS boundary and incoming wave at the horizon. These two boundary conditions guarantee that the poles of the retarded Green function are precisely the black-hole QNM spectrum. Here we obtain the two-point function following Ref. [16].
Let us start by writing the action (2.3) in the form where z 0 is a point close to the boundary and z h is the position of the event horizon. Since the equation of motion is satisfied, the first term in Eq. (3.29) is zero. Hence, the on-shell action reduces to the surface term .

(3.30)
After introducing the Fourier transform (3.9) and decomposing the field as S(k, z) = S 0 (k)S k (z), where k is the four-momentum, the on-shell action (3.30) can be written as This is similar to the result obtained in Ref. [16], with the main differences being the presence of the dilaton and the fact that (3.32) is related to a massive field in the bulk. Now the asymptotic solutions for the field S k (z), which satisfies an equation of motion in Fourier space identical to Eq. (3.10), need to be found. Here we also write the generic asymptotic expansion of the solution for a massive scalar bulk field close to the boundary [40], where ∆ is the conformal dimension of the operator dual to the scalar field S, and the ellipses mean subleading terms. Then, it is possible to write the field as S k (z) = z 4−∆ f k (z), where the function We use this result to obtain the retarded Green function following the prescription of Ref. [16], In order to get the asymptotic solution for the bulk scalar field S(k, z), we use the results of the asymptotic solutions close to the boundary, Eqs. (3.19) and (3.20), together with the relations (3.37) Now it is easy to get the explicit expression for f k (z) from the relation S(k, z) where we have used ∆ = 3 and S 0 (k) = A (−) (ω, q) to guarantee the condition (3.34). Finally, by using Eqs. (3.35), (3.36) and (3.37), it is obtained the Green function where the ellipses denote power corrections in z 0 . After a renormalization process [40] we take the limit z 0 → 0 to extract the finite part of (3.39). The imaginary part of such a result is related to the spectral function (SPF), which is given by It is worth mentioning that the retarded Green function could be obtained using the prescription presented in Ref. [41]. The authors used the canonical momentum associated to the massive scalar field to get this quantity, and showed that their result is consistent with the prescription presented in Ref. [16], at least for a massless scalar field. Here we also used the prescription of the canonical momentum for a massive scalar field and obtain the same result as Eq. (3.39). In the next subsection we present the numerical results for the spectral function (3.40) and an analysis of its dependence with the temperature and the dimensionless gluon condensate.

General procedure
Here we make a brief summary of the general procedure to obtain the SPFs in holographic QCD. In the present case, the spectral function is given by Eq. (3.40). The next step is to express the coefficients B (±) (ω, q) and A (±) (ω, q) in terms of the asymptotic solutions of ψ close to the boundary, as obtained in Sec. 3.3. The idea is to write the solutions and their derivatives as (a = 1, 2), (3.41) where a = 1 (a = 2) stands for the normalizable (non-normalizable) solutions. In matrix form, Eq. (3.41) reads The aim of this procedure is to get expressions for the coefficients D (a) and C (a) as a function of the asymptotic solutions. Hence, inverting the matrix product (3.42) we get the desired result. To be more specific, we are looking for the ratio since there is a connection between the coefficients D (1) , D (2) and B (−) , A (−) , given by With the analytic expressions in hand we solve numerically the Schrödinger-like equation (3.11) from a point close to the boundary (z 0 , a sufficiently small positive number, e.g., z 0 = 0.001,as we employed in the numerical analysis), using as "initial" conditions the asymptotic solutions close the boundary, up to a point close the horizon (z nh = z h − z 0 ).

Numerical results
Here we present and discuss some numerical results obtained following the general procedure explained above. Firstly, we calculate the SPFs by setting q = 0. This means that the spatial components of the momentum are neglected, so that the four-momentum is given by k µ = (ω, 0). Hence, the bound states in the field theory do not present spatial displacement. Figure 6 shows the numerical results for the spectral function in the low temperatures regime for G = 1 (left panel) and G = 50 (right panel). Each curve in a given panel is drawn for a given  temperature, as specified in the figure. Each give curve shows several sharpen peaks, and the numerical data around each peak fit well a Breit-Wigner function [37-39, 43, 44] where ω 0 is the position of the peak, Γ is the half-width of the curve around the peak, and A and b are constant parameters that can be determined by fitting the numerical data in each case. Each peak is interpreted as a quasiparticle state. The position of each peak of the curves in Fig. 6 is related to the mass of the corresponding quasiparticle state, while the half-width of the peak is related to the inverse of the quasiparticle lifetime. For instance, in the case G = 1 the mass at zero temperature, cf. Table 1, for the fundamental state is m 2 = 4.84 C, and it is possible to get the position of the first peak at ω 2 0 = 4.78 C for the finite temperature T = 0.2. Repeating the same comparison in the case G = 50 we get m 2 = 5.99 C and ω 2 0 = 5.90 C for T = 0.2. The difference between m 2 and ω 2 0 depends on the temperature for fixed G, i.e., ∆m 2 ( T ) = m 2 − ω 2 0 . Note that this difference is positive for nonzero temperature. However, at zero temperature the spectral function reduces to delta functions with peaks located at the values of the masses presented in Table 1. In such a limit the difference is zero, i.e., ∆m 2 (0) = 0. When the temperature is introduced the peaks gain width as it can be seen in Fig.6. As the temperature increases the width of the peaks increases while the position of the peaks are shifted. We also see the number of peaks decreasing while the temperature increases, signalling the melting of the quasi-particle states.
In Ref. [43] the SPFs for scalar mesons were obtained by using the soft-wall model. The position of the first peak obtained for T = 0.22 is at ω 2 0 = 5.85 C. In order to compare to our results, we have calculated the SPFs for this temperature, but we do not show the corresponding graphics here. For G = 1 the first peak is localized at ω 2 0 = 4.75 C, while for G = 50 it is localized at ω 2 0 = 5.85 C. From these results we observe, again, that the results of the soft-wall modes are obtained when the value of the parameter G is large, just as we observed when the spectrum was calculated, cf. Sec. 2.4. If one wants to know the temperatures in MeV it is possible to fix the parameter C. This can be done by comparing the first solution of the eigenvalue problem to the ρ meson mass, see Ref. [28] for details. The value obtained in [28] is √ C = 388 MeV, and then, the temperatures in Fig. 6 are in the range 17.3 MeV ≤ T ≤ 37.1 MeV. These temperatures are smaller than the value predicted by lattice calculations T c = 192(7)(4)MeV [45], and also by the soft-wall model which is T c = 191MeV [28].
Additionally, the results shown in Fig. 6 are in agreement with the results obtained from the effective potential in Fig. 5. For instance, in the analysis of the potential it was shown the possible existence of two quasiparticle states for T = 0.2 in both cases G = 1 and G = 50. This is also supported by the SPFs where there are two peaks for the same temperature. The effect of the parameter G shows up in shifting the position of the peak and its width (compare both panels in Fig. 6 for the same temperature). Hence, the position and width now are functions of T and G, i.e., ω 0 ( T , G) and Γ( T , G).
In the left panel of Fig. 7 we present a specific comparison to observe the effects of the parameter G. It is observed that this parameter shifts ω 0 to more energetic states, while the width of the peak increases.
To finish the analysis of the SPFs we add spatial momentum to the above results, i.e., q = 0. Hence, the quasiparticle states have spatial displacement and therefore become more energetic. The results for this case are presented in Fig. 8 for G = 1 (left panel) and G = 50 (right panel). The addition of spatial momentum shifts the position of the peaks to more energetic states, and the half-widths also increase while the spatial momentum increases. This means that more energetic states, i.e., quasiparticles with higher momenta, have shorter lifetime and therefore they melt faster than low energy states.

Quasinormal modes
In this section we calculate the spectrum of the quasinormal modes (QNMs) for the scalar field, i.e., S(x, z), lying on the black hole background. The back reaction on the black hole geometry is neglected, so that the scalar field is a prove field.

Power series method
The standard procedure to use the power series method is to transform the second order eigenvalue problem into a first order eigenvalue problem [46]. To do that we introduce the following transformation into the Schrödinger-like equation (3.11) ψ(z) = e −iωr * ϕ(z). (4.1) After replacing dz/dr * = −f (z) and simplifying, the following equation shows up where V (z) is given by Eq. (3.14). As in the previous sections, we label the two independent solutions of Eq. (4.2) by ψ (+) (the outgoing wave) and ψ (−) (the ingoing wave). Now we write the solution for ψ (−) as a power expansion around the horizon z = z h , where the coefficients a n = a n (ω, q, T, C, G) are functions of the parameters of the model, and by convention we set a 0 = 1. The spectrum of the QN frequencies is obtained by imposing Dirichlet condition at the boundary, i.e., ϕ (−) (0) = 0, which is equivalent to solve the following recursive equation This method, first employed in Ref. [46], works well for large AdS black holes, i.e., for high temperatures. For low temperatures the numerical convergence gets poor, as it was previously noticed in holographyc QCD models for scalar gluons and vector mesons [37,38].

Breit-Wigner method
As commented previously, it is known that the power series method is not reliable to find QN frequencies in the regime of low temperatures. In this regime, the form of the effective potential (cf. Sec. 3.4) allow us to apply a resonance method to obtain the QN frequencies. In this subsection we use the Breit-Wigner or resonance method to calculate QN frequencies of the scalar field. This method was applied for the very first time in Ref. [47] to compute QNMs for AdS black hole. In the context of holographic QCD this numerical method was previously applied in Refs. [37,38]. Now we make a brief summary of this method, for details see for instance Refs. [37,38,47] and references therein. Firstly, we write the normalizable solution of the Schrödinger-like equation, ψ (1) close the horizon, cf. Eq. (3.15), as ψ (1) = C (1) e −iωr * + D (1) e +iωr * = α (ω) cos(ωr * ) − β (ω) sin(ωr * ) . (4.5) The ingoing wave boundary condition at the horizon requires the vanishing of the coefficient D (1) . This means that at the QN frequency we might approximate it as D (1) ∼ (ω − ω QNM ), where ω QNM = ω R − i ω I . Note also that C (1) = (D (1) ) * in Eq. (4.5). Taking this information into consideration we calculate where ω is a real parameter. Hence, the real par of ω is obtained minimizing Eq. (4.6). In turn, the imaginary part of the frequency is obtained fitting a curve, i.e., a parabola, to Eq. (4.6).
Alternatively it is possible to obtain the imaginary part of the frequency using the following formulas [37,38,47] In the numerical analysis below, we take into consideration the value of the imaginary part obtained using these two procedures.

Pseudo-spectral method
The pseudo-spectral method [48] is used here to solve the linear eigenvalue problem (4.2). In this method, the regular solution of Eq. (4.2) is expanded as a function of cardinal functions C j (z). To do that we first make the substitution ϕ(z) = (z/z h ) 3/2 g(z), such that g(z) is a regular function in the interval [0, z h ]. The differential equation for g(z), which follows form Eq. (4.2), can be written in the form λ 2 (z, q, T, C, G) g (z) + λ 1 (z, ω, q, T, C, G) g (z) + λ 0 (z, ω, q, T, C, G) g(z) = 0, (4.8) where λ 2 , λ 1 and λ 0 are regular polynomials linear in ω. Now we expand the regular function g(z) as where C j (2z/z h − 1) are the Chebyshev polynomials, and the collocation points are chosen to be the Chebyshev-Gauss-Lobatto grid [48] z j = z h 2 1 − cos j π N . (4.10) Replacing (4.9) into (4.2) gives us the matrix form of the eigenvalue problem, where A and B are (N + 1) × (N + 1) matrices. It is not difficult to implement numerically a code to solve Eq. (4.11).
It is worth mentioning that the pseudo-spectral method inevitably leads to the emergence of spurious solutions that do not have any physical meaning. To eliminate the spurious solutions we use the fact that the relevant QN frequencies do not depend on the number of Chebyshev polynomials being considered in Eq. (4.9).

Numerical results
Here we present and discuss the numerical results obtained following the methods explained above. Firstly, in Fig. 9 we display the overlap of the results obtained using the power series, pseudo-spectral and Breit-Wigner methods for different values of the parameter G setting q = 0. The results of the real and imaginary parts of the QN frequency for G = 1 are displayed in top panels. In this figure the numerical results obtained using the Breit-Wigner method are those for which T 2 1. In turn, the power series method has a poor convergence and is unable to reach this region. The contrary occurs when the temperature increases, the Breit-Wigner method has a poor convergence while the power series method has a good one [46]. The last is valid for the first mode (n = 0), where a overlap of both methods is possible. For higher modes, i.e., n ≥ 2, this is not true and the intermediate region, for example, 0.04 ≤ T 2 ≤ 0.08, has no solutions. To fill the empty regions (where both methods get poor convergence) we use the pseudo-spectral method which a good convergence from low to large temperatures. In the region of intermediate temperatures there is a superposition of the methods and is possible compare the numerical results obtained using these methods, we do this below. In Fig. 9 we observe the real part of the frequency as a function of the temperature, i.e., see left panels in this figure. We observe how the value of ω R approaches to its zero value displayed in Table 1, i.e., spectrum at zero temperature, ω 2 R = (4.05 n + 5.17) with n = 0, 1, 2, ... , for G = 1 and ω 2 R = (4.00 n + 5.99) with n = 0, 1, 2, ... , for G = 50, see Sec. 2.4. In the intermediate region there is a smooth transition between the regime of low and high temperatures. The square of the real part becomes linear for higher values of the temperature, i.e, ω R ∝ T . It is interesting to point out that in the low temperatures regime, the QN frequencies may still be interpreted as quasiparticle states.
On the other hand, the imaginary part has a different behavior as a function of the temperature, it increases monotonically with the temperature, cf. right panels in Fig. 9. For low temperatures it decreases, and at zero temperature limit the imaginary part becomes zero. This means that the lifetime of the quasiparticle state becomes infinity, due to the definition of the lifetime τ = 1/ω I . For intermediate temperatures it increases as a power law, i.e., ω I ∝ T α , where the exponent α can be obtained by fitting the numerical results, for high temperatures it grows linearly with the temperature.
Additionally, we plot on the complex plane the numerical results obtained applying the pseudospectral method in Fig. 10, where the filled circles represent the results without dilaton field, i.e., G = 0, it is worth mentioning that these results are the same found in Ref. [50]. In this figure we observe how the poles structure changes with the value of the parameter G and in Table 2 we observe how the numerical results change with it.  Table 2. Numerical results of the first six nonhydrodynamic QNMs using the pseudo-spectral method for q = 0 and selected values of the temperature setting G = 0, G = 1 and G = 50. For completeness we compare the numerical results obtained using the numerical methods above. In Table 3 we display some selected values for the QN frequencies. We observe that the Breit-Wigner and pseudo-spectral methods work well for low temperatures, where the power series method has a poor convergence. Moreover, we observe that the Breit-Wigner method does not works well when the temperature increases, in contrast, the power series has a good convergence. It is interesting to note that the pseudo-spectral method has a good convergence for all the temperatures in Table 3. Hence, the pseudo-spectral method might be used instead of Breit-Wigner and power series methods. In general, the application of each method depends in the subject the research is focused. Finally, we realized that the convergence of these methods gets poor while the overtone studied increases, i.e., n = 2, 3, · · · .

Breit-Wigner
Power Series Pseudo-spectral  Table 3. Numerical results of the first nonhydrodynamic QNM for q = 0 and selected values of the temperature setting G = 50.

Comments on the imaginary part of the QN frequencies
To finish this section we comment some interesting facts about the imaginary part of the QN frequency. It would be interesting to know a relation between the equilibration timescales, i.e., τ ∝ 1/ω I , and the temperature of the dual thermal field theory. This is the main motivation for doing this. Let us start with the relation between the imaginary part of the frequency and the radius of the event horizon, i.e., r + (r + is radial coordinate of AdS black hole space expressed in global coordinates). In Ref. [46] was suggested that the imaginary part of the frequency scales with the horizon area for small black holes, i.e., ω I ∝ A, where A is the area of the event horizon of the black hole. Moreover, in five-dimensional spacetime the area is given by A ∝ r 3 + . This behavior was also confirmed using the Breit-Wigner method in Ref. [47] (for the scalar field perturbation) and see Ref. [49] for the results extended to arbitrary dimensions. On the other hand, for large black holes the imaginary part scales with the horizon radius as ω I ∝ r + . Hence, knowing the asymptotic behaviour of the imaginary part we may construct an interpolation function which recovers the asymptotic behavior. One possibility may be where the constants a 1 , b 1 and c 1 may be determined fitting this function with numerical results [46]. This discussion is valid for AdS black hole space in global coordinates. Now, we address the same analysis for AdS black hole space in Poincaré coordinates, which is the case we are dealing with in this work. As a matter of fact, our numerical results displayed in Fig. 9 show a linear behavior in the regime of large temperatures, i.e., ω I ∝ T . In turn, this relation is no longer valid in the regime of low temperatures. However, we know (from Fig. 9) that the imaginary part vanishes in the zero temperature limit, so that the relation must be in the form ω I ∝ T α , where α must be a positive number. So that we construct an interpolation function as where the constants may be determined fitting this function with our numerical results. Hence, the best fit we have found in the interval of temperatures 0.01 ≤ T 2 ≤ 0.5 is for the values: a = 95.0 × 10 −3 , b = 1.37, c = 15.3 × 10 −4 and α = 7.60. This fit was done for the case G = 1.

Dispersion relations
Here we study the momentum dependence of the QN frequencies for selected values of the temperature. These are the dispersion relations {ω R (q), ω I (q)} for the scalar field.
The dispersion relations for the first five QN frequencies, for selected values of the temperature and setting G = 1, are displayed in Fig. 11. The results for low temperatures are displayed in the top panels, while the results for high temperatures are displayed in bottom panels of the same figure. The obtained frequencies for low temperatures are consistent with the results presented in Sec. 3.6, where the real part of the frequency, i.e., the location of the peaks of the SPFs, increases with the wave number, for the same temperature, cf. Fig. 6. Moreover, the imaginary part of the frequency, i.e., the width of the peaks in Fig. 6, also increases, but very slowly, with the wave number. This behaviour holds for low temperatures where quasiparticle states are present in the quark-gluon plasma.
The results for intermediate and high temperatures are displayed in bottom panels of Fig. 11. The real part of the frequency grows with the wave number, while the imaginary part decreases. It is worth pointing out that previous studies in the literature have calculated QN frequencies of massive scalar fields without dilaton field. The results found have a similar behavior of the imaginary part of the frequency, i.e., frequency decreasing (see for instance Ref. [50]). In general, for intermediate and large temperatures, this kind of behaviour seems to be a universal property and is shared by QNMs in several gravitational theories (see, e.g., the reviews of Refs. [22,23] and the references therein).
A second particular case is show in Fig. 12, where we plot the QN frequencies as a function of the normalized wave number for G = 50. Now a comparison between Figs. 11 and 12 is in order. The real part of the frequency calculated for G = 50 is larger than the one calculated for G = 1 (for the same low temperature). The same is true for the imaginary part of the frequency, i.e., ω I (G = 50) > ω I (G = 1). This means that the thermalization timescales defined as τ = 1/ω I   Table 4. Numerical results of the first nonhydrodynamic QNM for q = 1 and low temperatures.
becomes smaller when the value of the parameter G increases, see also Table 4. Observing carefully the top right panel of Fig. 12 we realized that for T 2 ≥ 0.40 the imaginary part decreases slightly, this is different from the corresponding result for G = 1, where the frequency increases. On the other hand, for higher values of the temperature the behavior of the imaginary part is different ω I (G = 50) < ω I (G = 1), meaning that the thermalization timescales is shorter for small G, see also Table 5.
To explicitly verify the dependence of the QN frequencies on the temperature we show in Tables 4 and 5 some selected numerical results for a fixed wave number valueq = 1, and the same two values of parameter G in Figs. 11 and 12. The dependence on the temperature is seen by comparing the two tables, while the dependence on the values of the parameter G is seen within each table.   Table 5. Numerical results of the first nonhydrodynamic QNM for q = 1 and high temperatures.

Final remarks and conclusion
In this work we have implemented the effects of the gluon condensate in the analysis of the spectrum and melting of scalar mesons in holographic QCD. Additionally, in the gravitational background, we determine the spectrum of QNMs and observe how the results depend on the parameter associated with the gluon condensate. The effects of gluon condensate were introduced by considering a quartic dilaton in the UV and quadratic in the IR (in order to guarantee confinement and linear Regge trajectories). The results obtained show that the spectrum at zero temperature is sensitive to the value of the energy scale associated with the condensate, i.e., G. For large values of this parameter we obtain the same spectrum as obtained in the original soft-wall model applied to the case of scalar mesons. On the other hand, when this parameter is zero we recover the problem of an AdS spacetime with constant dilaton field. In this case the conformal symmetry is restored and the spectrum becomes continuum.
The results are different when a black hole is embedded in the gravitational background, since the presence of the black hole introduces a temperature, which means that the dual field theory has finite temperature. In this case the conformal symmetry is broken by the dilaton and temperature. Differently from the zero temperature case, it is possible to set G = 0 and the results show a discrete spectrum, at least in the gravitational side, because the potential has a well. Our results also show that the potential well depends directly on the dilensionless parameters G and T , hence, these parameters deform the potential. This means that there are more trapped quasiparticle states at low temperatures, where the potential has wells (cf. Fig. 5), these results are supported nicely by the SPFs (cf. 6). In the low temperatures regime we observe the presence of peaks in the SPFs, these peaks disappear when the temperature increases characterizing the melting of scalar mesons. We also observe higher excited states of quasiparticles melting faster than low excited states for the same temperature when the energy increases. These effects are accelerated when spatial momentum is added.
We complemented this work by calculating the QNMs spectra, the numerical results show a finite real part of the frequency in the limit of zero temperature, while the imaginary part becomes zero (cf. Fig. 9). In this limit the QNMs become normal modes. As a complementary analysis we used three numerical techniques to calculate the spectrum, and from a comparison, we realized that they work well in a determinate regime of temperatures and they can be used for specific purposes. At the end we present the numerical results and discussion of the dispersion relations. In all the results obtained in this work we realized the relevance of the parameter G, and how it is important to take into account when the melting of scalar mesons in a holographic model for QCD is studied.