Meson excitation at finite chemical potential

We consider a probe stable meson in the holographic quark–gluon plasma at zero temperature and chemical potential. Due to the energy injection into the plasma, the temperature and chemical potential are increased to arbitrary finite values in such a way that the plasma experiences an out-of-equilibrium process. We then observe that the meson is excited, i.e. the expectation value of Wilson loop oscillates around its static value with a specific angular frequency. By defining excitation time tex\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$t_{ex}$$\end{document} as a time at which the meson falls into the final excited state, we study the effect of various parameters of theory on the excitation time and observe that for larger values of final temperature and chemical potential the excitation time increases. Furthermore, our outcomes show that the more stable mesons, meaning that the meson with lower static potential, are excited sooner.


Introduction
Standard methods applied to calculations in quantum chromodynamics (QCD) are often based on a perturbative expansion. This makes them of limited use if applied to questions about the non-perturbative nature of the theory at low energies and where the coupling constant is large. Strong coupling techniques are in particular needed for understanding the physics of the strongly coupled quark-gluon plasma (QGP) produced at RHIC and at the LHC [1]. One such technique is lattice gauge theory that successfully describes the low energy properties of QCD especially at μ = 0 [2]. Lattice gauge theory, however, is of limited use in describing time-dependent quantities and non-equilibrium evolution. Another technique is gauge-gravity duality that we apply in this paper to explain the behavior of meson in nonequilibrium plasma in the presence of non-zero chemical potential.
The anti de-Sitter/conformal field theory (AdS/CFT) correspondence or more generally gauge-gravity duality [1,[3][4][5] is a conjectured relation between two physical theories. One of them is a strongly coupled gauge theory in d dimensional space-time and the other one is a classical gravity theory living in an extra dimension of space-time. In fact, parameters, fields and different processes in the gauge theory are translated into appropriate equivalent on the gravity side. For instance, the (thermal) vacuum state on the gauge theory side corresponds to the (black hole-AdS) pure AdS in the gravity theory. Thermalization process, which generally means evolution of a state from zero temperature to a thermal state, is dual to black hole formation in the gravity theory [6]. Moreover, as another example, the meson, quark-antiquark bound state, living in the QGP can be identified with a classical string in the gravity and by using the expectation value of the Wilson loop the static potential between a quark and antiquark has been firstly found in [7]. For more information, the interested reader is referred to [1] and references therein.
In this paper, we consider a probe stable meson in the QGP at zero temperature and chemical potential. Then the temperature and the chemical potential are simultaneously raised from zero to finite values T f and μ f , respectively. Now the questions we would like to answer are how the stable meson reacts to the energy injection into the system and what the characteristics of the new meson state are? Furthermore, it is instructive to know how much time is needed for the meson to fall into the final excited state, i.e. excitation time, and what the effect of the final values of the temperature and the chemical potential is on the excitation time? By excited state we mean a meson which has a higher energy than the ground state, i.e. the meson in the initial static gauge theory. The holographic dual of the above system is described by the dynamics of the classical string, with appropriate initial and boundary conditions, in the Reissner-Nordström-AdS Vaidya (RN-AdS-Vaidya) background, as we will review in the next section.

Review on the backgrounds
The gauge-gravity duality proposes a promising approach to investigate different properties of the strongly coupled field theory. Since we want to study the QGP as a strongly coupled system in the presence of non-zero chemical potential, we firstly review its corresponding holographic dual, i.e. RN-AdS background. We then extend our problem to the timedependent case, i.e. RN-AdS-Vaidya, which is dual to the thermalization process in the strongly coupled field theory when the temperature and the chemical potential simultaneously increase.

RN-AdS black hole background
Here we introduce charged black hole metric which is asymptotically AdS. Consider the Einstein-Maxwell anti-de Sitter action [8]: where G is Newton constant, R is Ricci scalar,F is field strength of the U (1) gauge field and R is AdS radius which we set to be one. The number of spatial directions is d and it relates to the negative cosmological constant as Equations of the motion obtaining from the above action are where μ = 0, . . . , d. The solution of the equations of motion is the RN-AdS metric 1 which can be written as [8] and the time component of the gauge field introduced in (6). M and Q are the mass and the charge of the RN-AdS black hole, respectively. In (3), z is radial coordinate and z = 0 is the AdS boundary. In addition, (t, x) are the four dimensional coordinates at the boundary. The gauge-gravity duality indicates that the Hawking temperature of the black hole corresponds to the temperature of the QGP. For d = 4, the temperature of the RN-AdS 5 black hole is Here z h is the radius of the event horizon, i.e. the smallest root of f (z) = 0. The relation between z h , M and Q is Moreover, the time component of the gauge field is given by [8] A t = − where is a constant which plays the role of the electrostatic potential between the boundary of the bulk and the event horizon. It can be defined by applying A t (z h ) = 0. Thus we obtain Gauge-gravity duality provides a correspondence between the time component of gauge field at the boundary and the chemical potential in dual boundary gauge theory [8,11], i.e.
in the AdS radius unit. Hence, it is easy to find that

RN-AdS-Vaidya background
Let's generalize the static background (3) to the timedependent case. It can be achieved by adding external source terms to the action (1). Then, the equations of motion for general d are [8] 8π The RN-AdS-Vaidya metric, in the Eddington-Finkelstein coordinates, is then given by provided that and Q(v) have been discussed in the literature and it seems that the physical results are independent of the form of the functions [12]. Therefore, we choose where I ∈ (M, Q). The transition time k is the time interval that the mass (charge) of the black hole needs to reach its final value M f (Q f ). For k 1 (k 1) the transition time is small (large) which is usually called fast (slow) quench. According to the gauge-gravity duality, study of the RN-AdS-Vaidya metric on the gravity side corresponds to the study of thermalization of strongly coupled QGP in the presence of the chemical potential in the gauge theory. Note that the relation between M f , Q f and z h is still given by (5).

Expectation value of Wilson loop
The static potential between a quark and anti-quark has been extensively studied, for example see [1,[13][14][15]. In this section, we want to review how this static potential can be obtained from the duality point of view. To do so, we use the expectation value of the Wilson loop as a gauge invariant operator. Specifically, on the one hand, we use the time-like Wilson loop in rectangular form C. One side of rectangle, l, is spatial corresponding to the distance between quarkantiquark pair and the other side is temporal, T . If we assume that T l, meaning that the world-sheet is translationally invariant along the time direction, the expectation value of the Wilson loop is [1] where m is the rest mass of quark (antiquark) and V (l) represents the static potential energy between the pair. On the other hand, the gauge-gravity duality proposes that in order to calculate the static potential in the QGP at finite chemical potential, one needs to probe the RN-AdS black hole geometry by the classical string [16,17]. In other words, the expectation value of the Wilson loop, in the saddle point approximation, is dual to the on-shell action of classical open string that its end points are located on the boundary with distance l. Therefore, where S(C) is value of the Nambu-Goto action on the rectangle C. The tension of string is proportional to α −1 ≡ l s −2 and l s is the fundamental length scale of string.
Also, τ and σ parametrize the two dimensional string worldsheet. g ab = G μν ∂ X μ ∂ξ a ∂ X ν ∂ξ b is the induced metric on the worldsheet. Here G μν and X μ (ξ a = τ, σ ) are the metric and the bulk (world-sheet) coordinates, respectively. This action, in fact, describes the dynamics of the classical string in any desired geometry. As we will see, using (14)-(16), the static potential can be found. Note that, the rest mass m is equal to , where is IR regulator in the gravity theory and according to UV/IR connection, it corresponds to the UV cut-off in the gauge theory.
In order to calculate the S(C) in the RN-AdS black hole background (3), we parametrize the two-dimensional worldsheet of the string as τ = t, σ = x 3 ≡ x. All bulk coordinates, except z and x, are chosen to be constant and therefore the shape of the string is described by z = z(x). Hence the action (16) reduces to where z = dz/dx. Since the Lagrangian does not depend explicitly on x, the associated Hamiltonian is a constant of motion. After some simple algebra, we get where z = z * at z (x) = 0. Using the new coordinate y = z * /z and the explicit form of f (z), it turns out where y h = z * /z h . The on-shell action can be then found by inserting (18) and (19) into (17). Furthermore, by applying (14) and (15) the potential energy between the quark and anti-quark will be obtained that is in agreement with results reported in [18]. Notice that, in order to get a finite value for the potential energy, we need to subtract the rest mass of the quarks from the on-shell action and we therefore have Before closing this part, two points should be noted. First, as we will see, the Eq. (19) will be utilized as an initial condition for the time evolution of the classical string in the RN-AdS-Vaidya background. Second, in the next section, we observe that the expectation value of the time-dependent Wilson loop oscillates around the static potential resulted in (20). Now our goal is to calculate the evolution of expectation value of Wilson loop in the RN-AdS-Vaidya geometry. In this geometry, the two-dimensional world-sheet of the string is no longer translationally invariant in the time direction. Thus the condition T l does not valid and therefore in this case, the expectation value of the Wilson loop (14) can be written as where, using the gauge-gravity duality, W(t) is the on-shell action of string without integrating over t-coordinate. It is important to note that although W is a function of different parameters in the gauge theory, such as temperature, chemical potential and distance l, we only show its time dependence, explicitly. As before, W(t), or equivalently the onshell string action, diverges and in order to regularize this divergence, similar to the static case, we have The subscript "R" in (22) refers to the regularized version of W(t). To calculate the string on-shell action, similar to [19][20][21], we choose the null coordinates (u, v) to parametrize the two-dimensional world-sheet of the string and therefore all the coordinates on the world-sheet depend on u and v. Then we choose the following ansatz: Note that due to rotational symmetry in spatial coordinates, there is no difference among them. Substituting this ansatz into the Nambo-Goto action (16), the equations of motion can be found. Hence, we finally have where, for an arbitrary function A(a, b), we define A ,a = ∂ A ∂a and A ,ab = ∂ 2 A ∂a∂b . Since u and v are null coordinates, we need to impose two constraints corresponding to g uu = 0 and g vv = 0. Thus, two constraint equations turn out to be In order to solve the equations of motion (24) subject to constraint equations (25), the suitable boundary and initial conditions are essential. In the Appendices B and C, we obtain the appropriate boundary and initial conditions, respectively. After solving the equations of motion (24), one can easily find the time evolution of W R (t) when the other parameters are kept fixed.

Numerical results
As we have already mentioned and shown by the Figs. 1, 2 , 3 and 4 the expectation value of the Wilson loop oscillates around the static potential. In fact, before injecting energy (t < 0), the quark-antiquark pair is in ground state. However, as the energy injection is started, the temperature and the chemical potential increases and therefore the pair is excited. When the energy injection ceases, the pair falls into a final  (26) excited state with specific frequency and amplitude of oscillation. In other words, these oscillations can be interpreted as if the energy injection puts the meson in the excited state [19][20][21]. As a consequence of conformal invariance at finite temperature, there are two independent parameters, say μ f /T f and lT f . We are now interested in investigating the effect of these parameters on the characteristic of oscillation, i.e. frequency and amplitude. Thus our results are classified into the following categories: We observe that by raising the final value of the temperature of the QGP while the distance between the quark and antiquark is kept fixed, the amplitude of the oscillation increases, also (the energy of the excited pair increases based on the relation between amplitude and energy in simple harmonic oscillator). However, the value of the frequency is independent of the final temperature in agreement with the numerical results of [22]. In Fig. 2, we have lT f = 0.0714 (left panel) and lT f = 0.1428 (right panel) for fixed μ f /T f = 3.15. In this figure with increasing the distance between the quark and antiquark while the temperature is kept fixed, the amplitude (frequency) of oscillation increases (decreases). An interesting point we would like to emphasize here is that since we are working in the probe limit, the energy of the meson does not dissipate in the plasma and therefore the oscillation remains unchanged.
Up to now, we investigate the response of the system to the time-dependent change in the temperature and the chemical potential. In fact, the response of the system is described by the behavior of the expectation value of the Wilson loop, W R (t), in terms of boundary time. We observe that the expectation value oscillates around the static potential with a specific value of frequency, ν, and amplitude of oscillation, A. The characteristics of oscillation depend on the final values of the chemical potential and the temperature as well as the distance between the pair. Now, the question we would like to answer is, at which time the expectation value of the Wilson loop starts oscillating around the static potential? And we call it excitation time. To do so, we define the following function: where A, ν and φ can be found from W R (t) at asymptotic times. We then define a time-dependent function where in (27a) T and l are kept fixed while μ changes and so on. Thus the excitation time, t ex , is defined as the time which satisfies (t ex ) < 5 × 10 −6 and (t) stays below this limit afterwards. 2 In Fig. 5 (left), blue points, the final temperature of the system decreases while the other parameters contributing to the time evolution of the expectation value are kept fixed. As it is clearly seen, at higher temperatures the meson falls into the excited state later. In other words, the higher temperature, the larger (rescaled) excitation time, (k −1 )t ex . This result can be intuitively understood: at higher temperatures the thermal fluctuations become more significant i.e. in any system with non-zero temperature, the thermal fluctuations are present and as the temperature increases, the thermal fluctuations become more highlight. In other words, the particles gain more kinetic energy and bounce more frequently off of one another, so that can lead to more deviation of the (state of) particles from their steady state. Thus this effect prevents the meson from sooner fall. In the same figure, red points, since we would like to investigate the effect of the chemical potential on the excitation time, the chemical potential is varied and the other parameters are kept fixed. Similar to the previous case, this figure indicates that the rescaled excitation time and chemical potential increase with respect to each other. As a matter of fact, for larger values of the temperature and the chemical potential, the rescaled excitation time becomes larger. In short, the effect of temperature on the rescaled excitation time is similar to the chemical potential. Apart from the fast quench, in the right panel, the rescaled excitation time is plotted as a function of μ f /T f for the case of the slow quench, i.e. k = 3. It is obviously seen that in both cases the behavior of excitation times is the same. However, k −1 t ex is larger for the case of the fast quench.
In Fig. 6, the distance l is increased while the other parameters are kept fixed. In contrast to the chemical potential and the temperature, the distance is an intrinsic characteristic of the meson in the plasma. This figure shows that for larger value of distance l, meaning that the meson is less stable, the excitation time increases. Here, by stability we mean that the value of the static potential is lower at zero temperature and chemical potential, or equivalently before the energy injection, for smaller values of distance l. Although for both quenches the rescaled excitation time behaves similarly, the value of k −1 t ex is larger for the fast quench. In other words, this figure indicates that the meson with smaller l, i.e. the more stable quark-antiquark bound states are excited sooner. It may be related to the screening of the force between color charges of the quark and antiquark due to the presence of the medium since they, quark and antiquark, can not communicate, easily. Notice that the distance between quark and antiquark can not be too large because the meson will then dissociate in the plasma.

Conclusions
Our main outcomes can be summarized as follows: • We observe that the expectation value of the timedependent Wilson loop oscillates around the static potential. It may be interpreted as the string connecting quark and antiquark (or flux tube resulted from the gluon fields between the pair) in field theory with Dirichlet boundary conditions on the endpoints. After the energy injection, the string is oscillating in one of the its normal modes or equivalently the quark-antiquark bound state has been excited. As a matter of fact, we investigate the effect of the chemical potential and the temperature on this normal mode. • Making larger each parameter under study in this paper, i.e. final chemical potential, final temperature and the quark-antiquark distance, the energy of the bound state, based on classical harmonic oscillator model, increases since the amplitude of the oscillation becomes larger. • The oscillation frequency is independent of the temperature and approximately of the chemical potential. It is intuitively comprehensible. Similar to the classical harmonic oscillator, the oscillation frequency is an intrinsic characteristic of the meson and is independent of the environmental changes. However, the larger distance between quark and antiquark, the smaller oscillation frequency. • Consider a meson in the plasma with non-zero temperature and chemical potential. For larger values of the chemical potential and higher temperatures, the excitation time of the meson increases. In other words, when the plasma is hotter or denser the meson falls into the final excited state more slowly. • Consider a plasma where its temperature and chemical potential are kept fixed. Then the more stable meson at T = μ = 0, corresponding to the smaller distance between quark and antiquark, falls into the final excited state sooner. • All of the above results are confirmed for the slow and the fast quench. By slow (fast) quench we mean the energy injection into the system is done slowly (rapidly).
We would like to emphasize that, apart from the first result in the third item, the other outcomes can not be argued by comparing to the simple harmonic oscillator. In the case of the second item, one expects that while the energy is injected into the system by an external source the energy of the bound state increases. On the other hand, since our results indicate that the amplitude of the oscillation always increases after the energy injection, this reminds us of the relation between the energy and amplitude of the oscillation in the simple harmonic oscillator and we therefore conclude that the energy of the meson increases in the plasma. However, it is important to notice that this relation is not as easy as the simple harmonic oscillator, i.e. (energy) ∝ (amplitude) 2 , as shown by our numerical results. In other words, the amplitude of oscillation, in fact, is a complicated and non-linear function of the different variables such as temperature, chemical potential and timescale of the energy injection.
In addition, studying non-local observables [8,23] such as the two-point correlation function and the expectation value of Wilson loop as probes of thermalization process, show that the larger μ/T the larger thermalization time. However, by studying the two-point correlation function, the expectation value of Wilson loop and the entanglement entropy [24], the results indicate that for a fixed small value of lT (lT 1) and small values of μ/T the thermalization time decreases by increasing μ/T , thus plasma thermalizes faster. But, for larger values of μ/T the thermalization time increases with increasing μ/T . Also, by increasing the value of lT this nonmonotonic behavior becomes less pronounced and finally disappears, so that a monotonic behavior has been observed for lT 1. More recently, equilibration of a dynamical scalar operator is considered in the charged QGP during its equilibration [25]. The numerical outcomes show that the equilibration time can be a decreasing or increasing function of μ/T and depends on the energy injection, so that one cannot report the general behavior. Furthermore, in [26] studying non-hydrodynamic quasinormal modes of a scalar field show that the equilibration time decreases with increasing the chemical potential far from the critical point, while close to the critical point it would increase. On the other hand, the thermalization time always increases with increasing the chemical potential. In summary, although the excitation time increases for larger values of the chemical potential in our case, it should be instructive to check our outcomes in other gauge theories with holographic dual to find a general behavior.
and thus, by using (8), one can find In Fig. 7 (left panel) we set μ f /T f = 3.15. However, in the right panel T f is zero, i.e. large chemical potential. We see that the static potential of the pair decreases and the amplitude of the oscillation increases. But, the frequency of oscillation does not change. In Fig. 8, we have μ f /T f = 0.1 for the left panel and large μ for the right panel. The results are in complete agreement with the Fig. 7. In addition, note that these results show that, in the large μ limit, by increasing the chemical potential, the final static potential also rises. It is intuitively correct since the matter screens the interaction between the quark and antiquark.